Insights into the mechanism of C5aR inhibition by PMX53 via implicit solvent molecular dynamics simulations and docking

Background

The complement protein C5a acts by primarily binding and activating the G-protein coupled C5a receptor C5aR (CD88), and is implicated in many inflammatory diseases. The cyclic hexapeptide PMX53 (sequence Ace-Phe-[Orn-Pro-dCha-Trp-Arg]) is a full C5aR antagonist of nanomolar potency, and is widely used to study C5aR function in disease.

Results

We construct for the first time molecular models for the C5aR:PMX53 complex without the a priori use of experimental constraints, via a computational framework of molecular dynamics (MD) simulations, docking, conformational clustering and free energy filtering. The models agree with experimental data, and are used to propose important intermolecular interactions contributing to binding, and to develop a hypothesis for the mechanism of PMX53 antagonism.

Conclusion

This work forms the basis for the design of improved C5aR antagonists, as well as for atomic-detail mechanistic studies of complement activation and function. Our computational framework can be widely used to develop GPCR-ligand structural models in membrane environments, peptidomimetics and other chemical compounds with potential clinical use.

Apela exhibits isoform- and headgroup-dependent modulation of micelle binding, peptide conformation and dynamics

Apela (also referred to as ELABELA and toddler) is a peptide hormone that activates the apelin receptor (AR or APJ) to regulate cardiovascular system development and function. Here, we report the first biophysical characterization of three apela isoforms, apela-54, -32, and -11, alongside a monomeric C1S-apela-11 mutant, using circular dichroism (CD) spectropolarimetry and nuclear magnetic resonance (NMR) spectroscopy. The behaviour of apela-54 is consistent with a preprotein containing a hydrophobic, N-terminal signal peptide. The potential for apela-membrane binding, leading to membrane catalyzed interactions with AR, was tested comprehensively for apela-32 and -11 in the presence of membrane-mimetic dodecylphosphocholine (DPC), sodium dodecyl sulfate (SDS), and 1-palmitoyl-2-hydroxy-sn-glycero-3-[phospho-rac-(1-glycerol)] (LPPG) micelles. According to pulsed-field gradient diffusion NMR experiments, apela-32 interacts with all three micelles. Chemical shift perturbations indicate widespread interactions along apela, with DPC and LPPG micelles inducing short segments with α-helical character at distinct regions. Consistent with these data, ps-ns dynamics along the peptide backbone appear decreased in the presence of micelles. Apela-11 and C1S-apela-11, alternatively, interact preferentially with SDS and LPPG micelles, promoting β-turn character observable by CD. Distinct differences in membrane-interaction propensity are therefore apparent both as a function of apela isoform and of detergent headgroup. These results imply the potential for cell membrane involvement in apela-AR recognition and binding, with the implication that membrane catalysis has distinct functional and regulatory roles throughout the apelinergic system.

Early Events in β2AR Dimer Dynamics Mediated by Activation-Related Microswitches

G-Protein-Coupled Receptors (GPCRs) make up around 3-4% of the human genome and are the targets of one-third of FDA-approved drugs. GPCRs typically exist as monomers but also aggregate to form higher-order oligomers, including dimers. β2AR, a pharmacologically relevant GPCR, is known to be targeted for the treatment of asthma and cardiovascular diseases. The activation of β2AR at the dimer level remains under-explored. In the current study, molecular dynamics (MD) simulations have been performed to understand activation-related structural changes in β2AR at the dimer level. The transition from inactive to active and vice versa has been studied by starting the simulations in the apo, agonist-bound, and inverse agonist-bound β2AR dimers for PDB ID: 2RH1 and PDB ID: 3P0G, respectively. A cumulative total of around 21-μs simulations were performed. Residue-based distances, RMSD, and PCA calculations suggested that either of the one monomer attained activation-related features for the apo and agonist-bound β2AR dimers. The TM5 and TM6 helices within the two monomers were observed to be in significant variation in all the simulations. TM5 bulge and proximity of TM2 and TM7 helices may be contributing to one of the early events in activation. The dimeric interface between TM1 and helix 8 were observed to be well maintained in the apo and agonist-bound simulations. The presence of inverse agonists favored inactive features in both the monomers. These key features of activation known for monomers were observed to have an impact on β2AR dimers, thereby providing an insight into the oligomerization mechanism of GPCRs.