Cryo-EM structure of the human chemerin receptor 1-Gi protein complex bound to the C-terminal nonapeptide of chemerin.

Chemerin is a processed protein that acts on G protein-coupled receptors (GPCRs) for its chemotactic and adipokine activities. The biologically active chemerin (chemerin 21-157) results from proteolytic cleavage of prochemerin and uses its C-terminal peptide containing the sequence YFPGQFAFS for receptor activation. Here we report a high-resolution cryo-electron microscopy (cryo-EM) structure of human chemerin receptor 1 (CMKLR1) bound to the C-terminal nonapeptide of chemokine (C9) in complex with Gi proteins. C9 inserts its C terminus into the binding pocket and is stabilized through hydrophobic interactions involving its Y1, F2, F6, and F8, as well as polar interactions between G4, S9, and several amino acids lining the binding pocket of CMKLR1. Microsecond scale molecular dynamics simulations support a balanced force distribution across the whole ligand-receptor interface that enhances thermodynamic stability of the captured binding pose of C9. The C9 interaction with CMKLR1 is drastically different from chemokine recognition by chemokine receptors, which follow a two-site two-step model. In contrast, C9 takes an "S"-shaped pose in the binding pocket of CMKLR1 much like angiotensin II in the AT1 receptor. Our mutagenesis and functional analyses confirmed the cryo-EM structure and key residues in the binding pocket for these interactions. Our findings provide a structural basis for chemerin recognition by CMKLR1 for the established chemotactic and adipokine activities.

Enhanced-Sampling Simulations for the Estimation of Ligand Binding Kinetics: Current Status and Perspective.

The dissociation rate (k off) associated with ligand unbinding events from proteins is a parameter of fundamental importance in drug design. Here we review recent major advancements in molecular simulation methodologies for the prediction of k off. Next, we discuss the impact of the potential energy function models on the accuracy of calculated k off values. Finally, we provide a perspective from high-performance computing and machine learning which might help improve such predictions.

Structural basis for recognition of N-formyl peptides as pathogen-associated molecular patterns.

The formyl peptide receptor 1 (FPR1) is primarily responsible for detection of short peptides bearing N-formylated methionine (fMet) that are characteristic of protein synthesis in bacteria and mitochondria. As a result, FPR1 is critical to phagocyte migration and activation in bacterial infection, tissue injury and inflammation. How FPR1 distinguishes between formyl peptides and non-formyl peptides remains elusive. Here we report cryo-EM structures of human FPR1-Gi protein complex bound to S. aureus-derived peptide fMet-Ile-Phe-Leu (fMIFL) and E. coli-derived peptide fMet-Leu-Phe (fMLF). Both structures of FPR1 adopt an active conformation and exhibit a binding pocket containing the R2015.38XXXR2055.42 (RGIIR) motif for formyl group interaction and receptor activation. This motif works together with D1063.33 for hydrogen bond formation with the N-formyl group and with fMet, a model supported by MD simulation and functional assays of mutant receptors with key residues for recognition substituted by alanine. The cryo-EM model of agonist-bound FPR1 provides a structural basis for recognition of bacteria-derived chemotactic peptides with potential applications in developing FPR1-targeting agents.

Accuracy of Molecular Simulation-Based Predictions of koff Values: A Metadynamics Study.

The koff values of ligands unbinding to proteins are key parameters for drug discovery. Their predictions based on molecular simulation may under- or overestimate experiment in a system- and/or technique-dependent way. Here we use an established method-infrequent metadynamics, based on the AMBER force field-to compute the koff of the ligand iperoxo (in clinical use) targeting the muscarinic receptor M2. The ligand charges are calculated by either (i) the Amber standard procedure or (ii) B3LYP-DFT. The calculations using (i) turn out not to provide a reasonable estimation of the transition-state free energy. Those using (ii) differ from experiment by 2 orders of magnitude. On the basis of B3LYP DFT QM/MM simulations, we suggest that the observed discrepancy in (ii) arises, at least in part, from the lack of electronic polarization and/or charge transfer in biomolecular force fields. These issues might be present in other systems, such as DNA-protein complexes.

Proton Dynamics in Protein Mass Spectrometry.

Native electrospray ionization/ion mobility-mass spectrometry (ESI/IM-MS) allows an accurate determination of low-resolution structural features of proteins. Yet, the presence of proton dynamics, observed already by us for DNA in the gas phase, and its impact on protein structural determinants, have not been investigated so far. Here, we address this issue by a multistep simulation strategy on a pharmacologically relevant peptide, the N-terminal residues of amyloid-ß peptide (Aß(1-16)). Our calculations reproduce the experimental maximum charge state from ESI-MS and are also in fair agreement with collision cross section (CCS) data measured here by ESI/IM-MS. Although the main structural features are preserved, subtle conformational changes do take place in the first ∼0.1 ms of dynamics. In addition, intramolecular proton dynamics processes occur on the picosecond-time scale in the gas phase as emerging from quantum mechanics/molecular mechanics (QM/MM) simulations at the B3LYP level of theory. We conclude that proton transfer phenomena do occur frequently during fly time in ESI-MS experiments (typically on the millisecond time scale). However, the structural changes associated with the process do not significantly affect the structural determinants.