C5a receptor (C5aR) is one of the major chemoattractant receptors of the druggable proteome that binds C5a, the proinflammatory polypeptide of complement cascade, triggering inflammation and SEPSIS. Here, we report the model structures of C5aR in both inactive and peptide agonist (YSFKPMPLaR; a=D-Ala) bound meta-active state. Assembled in CYANA and evolved over molecular dynamics (MD) in POPC bilayer, the inactive C5aR demonstrates a topologically unique compact heptahelical bundle topology harboring a β-hairpin in extracellular loop 2 (ECL2), derived from the atomistic folding simulations. The peptide agonist bound meta-active C5aR deciphers the “site2” at an atomistic resolution in the extracellular surface (ECS), in contrast to the previously hypothesized inter-helical crevice. With estimated Ki≈2.75 μM, the meta-active C5aR excellently rationalizes the IC50 (0.1–13 μM) and EC50 (0.01–6 μM) values, displayed by the peptide agonist in several signaling studies. Moreover, with Ki≈5.3×105 μM, the “site2” also illustrates selectivity, by discriminating the stereochemical mutant peptide (YSFkPMPLaR; k=D-Lys), known to be inert toward C5aR, up to 1 mM concentration. Topologically juxtaposed between the structures of rhodopsin and CXCR1, the C5aR models also display excellent structural correlations with the other G-protein coupled receptors (GPCRs). The models elaborated in the current study unravel many important structural insights previously not known for regulating the agonist binding and activation mechanism of C5aR.

The interaction of hC5a with C5aR, previously hypothesized to involve a “two-site” binding, (i) recognition of the bulk of hC5a by the N-terminus (NT) of C5aR (“site1”), and (ii) recognition of C-terminus (CT) of hC5a by the extra cellular surface (ECS) of the C5aR (“site2”). However, the pharmacological landscapes of such recognition sites are yet to be illuminated at atomistic resolution. In the context, unique model complexes of C5aR, harboring pharmacophores of diverse functionality at the “site2” has recently been described. The current study provides a rational illustration of the “two-site” binding paradigm in C5aR, by recruiting the native agonist hC5a and engineered antagonist hC5a(A8). The hC5a-C5aR and hC5a(A8)-C5aR complexes studied over 250 ns of molecular dynamics (MD) each in POPC bilayer illuminate the hallmark of activation mechanism in C5aR. The intermolecular interactions in the model complexes are well supported by the molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) based binding free energy calculation, strongly correlating with the reported mutational studies. Exemplified in two unique and contrasting molecular complexes, the study provides an exceptional understanding of the pharmacological divergence observed in C5aR, which will certainly be useful for search and optimization of new generation “neutraligands” targeting the hC5a-C5aR interaction.

The C5a receptor (C5aR) is a pharmacologically important G-protein coupled receptor (GPCR) that interacts with (h)C5a, by recruiting both the "orthosteric" sites (site1 at the N-terminus and site2 at the ECS, extra cellular surface) on C5aR in a two site-binding model. However, the complex pharmacological landscape and the distinguishing chemistry operating either at the "orthosteric" site1 or at the functionally important "orthosteric" site2 of C5aR are still not clear, which greatly limits the understanding of C5aR pharmacology. One of the major bottlenecks is the lack of an experimental structure or a refined model structure of C5aR with appropriately defined active sites. The study attempts to understand the pharmacology at the "orthosteric" site2 of C5aR rationally by generating a highly refined full-blown model structure of C5aR through advanced molecular modeling techniques, and further subjecting it to automated docking and molecular dynamics (MD) studies in the POPC bilayer. The first series of structural complexes of C5aR respectively bound to a linear native peptide agonist ((h)C5a-CT), a small molecule inverse agonist (NDT) and a cyclic peptide antagonist (PMX53) are reported, apparently establishing the unique pharmacological landscape of the "orthosteric" site2, which also illustrates an energetically distinct but coherent competitive chemistry ("cation-π" vs. "π-π" interactions) involved in distinguishing the established ligands known for targeting the "orthosteric" site2 of C5aR. Over a total of 1 μs molecular dynamics (MD) simulation in the POPC bilayer, it is evidenced that while the agonist prefers a "cation-π" interaction, the inverse agonist prefers a "cogwheel/L-shaped" interaction in contrast to the "edge-to-face/T-shaped" type π-π interactions demonstrated by the antagonist by engaging the F275(7.28) of the C5aR. In the absence of a NMR or crystallographically guided model structure of C5aR, the computational model complexes not only provide valuable insights for understanding the C5aR pharmacology, but also emerge as a promising platform for the design and discovery of future potential drug candidates targeting the (h)C5a-C5aR signaling axes.

No abstract