Can Simulations and Modeling Decipher NMR Data for Conformational Equilibria? Arginine–Vasopressin

Computationally Assisted Analysis of NMR Chemical Shifts as a Tool in Conformational Analysis

A key to understanding the properties of functional molecules is to determine their conformation in solution. A conformational analysis procedure that relies on quantum mechanical calculations and the widely used DP4+ probability was evaluated to decipher the structural information encoded in NMR chemical shifts. The results underscore the potential utility of using NMR chemical shifts in advancing conformational analysis studies of complex molecules in solution.

Extended Metadynamics Protocol for Binding/Unbinding Free Energies of Peptide Ligands to Class A G-Protein-Coupled Receptors

A metadynamics protocol is presented to characterize the binding and unbinding of peptide ligands to class A G-protein-coupled receptors (GPCRs). The protocol expands on the one previously presented for binding and unbinding small-molecule ligands to class A GPCRs and accounts for the more demanding nature of the peptide binding-unbinding process. It applies to almost all class A GPCRs. Exemplary simulations are described for subtypes Y1R, Y2R, and Y4R of the neuropeptide Y receptor family, vasopressin binding to the vasopressin V2 receptor (V2R), and oxytocin binding to the oxytocin receptor (OTR). Binding free energies and the positions of alternative binding sites are presented and, where possible, compared with the experiment.

Ensemble Generation for Linear and Cyclic Peptides Using a Reservoir Replica Exchange Molecular Dynamics Implementation in GROMACS

The profile of shapes presented by a cyclic peptide modulates its therapeutic efficacy and is represented by the ensemble of its sampled conformations. Although some algorithms excel at creating a diverse ensemble of cyclic peptide conformations, they seldom address the entropic contribution of flexible conformations and often have significant practical difficulty producing an ensemble with converged and reliable thermodynamic properties. In this study, an accelerated molecular dynamics (MD) method, namely, reservoir replica exchange MD (R-REMD or Res-REMD), was implemented in GROMACS ver. 4.6.7 and benchmarked on two small cyclic peptide model systems: a cyclized furin cleavage site of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (cyclo-(CGPRRARSG)) and oxytocin (disulfide-bonded CYIQNCPLG). Additionally, we also benchmarked Res-REMD on alanine dipeptide and Trpzip2 to demonstrate its validity and efficiency over REMD. For Trpzip2, Res-REMD coupled with an umbrella-sampling-derived reservoir generated similar folded fractions as regular REMD but on a much faster time scale. For cyclic peptides, Res-REMD appeared to be marginally faster than REMD in ensemble generation. Finally, Res-REMD was more effective in sampling rare events such as trans to cis peptide bond isomerization. We provide a GitHub page with the modified GROMACS source code for running Res-REMD at https://github.com/PlotkinLab/Reservoir-REMD.

Structure prediction of cyclic peptides by molecular dynamics + machine learning

Recent computational methods have made strides in discovering well-structured cyclic peptides that preferentially populate a single conformation. However, many successful cyclic-peptide therapeutics adopt multiple conformations in solution. In fact, the chameleonic properties of some cyclic peptides are likely responsible for their high cell membrane permeability. Thus, we require the ability to predict complete structural ensembles for cyclic peptides, including the majority of cyclic peptides that have broad structural ensembles, to significantly improve our ability to rationally design cyclic-peptide therapeutics. Here, we introduce the idea of using molecular dynamics simulation results to train machine learning models to enable efficient structure prediction for cyclic peptides. Using molecular dynamics simulation results for several hundred cyclic pentapeptides as the training datasets, we developed machine-learning models that can provide molecular dynamics simulation-quality predictions of structural ensembles for all the hundreds of thousands of sequences in the entire sequence space. The prediction for each individual cyclic peptide can be made using less than 1 second of computation time. Even for the most challenging classes of poorly structured cyclic peptides with broad conformational ensembles, our predictions were similar to those one would normally obtain only after running multiple days of explicit-solvent molecular dynamics simulations. The resulting method, termed StrEAMM (Structural Ensembles Achieved by Molecular Dynamics and Machine Learning), is the first technique capable of efficiently predicting complete structural ensembles of cyclic peptides without relying on additional molecular dynamics simulations, constituting a seven-order-of-magnitude improvement in speed while retaining the same accuracy as explicit-solvent simulations.

Conformational selection of vasopressin upon V1a receptor binding

The neuropeptide vasopressin (VP) and its three G protein-coupled receptors (V_1aR, V_1bR and V₂R) are of high interest in a wide array of drug discovery programs. V_1aR is of particular importance due to its cardiovascular functions and diverse roles in the central nervous system. The structure–activity relationships underpinning ligand-receptor interactions remain however largely unclear, hindering rational drug design. This is not least due to the high structural flexibility of VP in its free as well as receptor-bound states. In this work, we developed a novel approach to reveal features of conformational selectivity upon VP-V_1aR complex formation. We employed virtual screening strategies to probe VP’s conformational space for transiently adopted structures that favor binding to V_1aR. To this end, we dissected the VP conformational space into three sub-ensembles, each containing distinct structural sets for VP’s three-residue C-terminal tail. We validated the computational results with experimental nuclear magnetic resonance (NMR) data and docked each sub-ensemble to V_1aR. We observed that the conformation of VP’s three-residue tail significantly modulated the complex dissociation constants. Solvent-exposed and proline trans-configured VP tail conformations bound to the receptor with three-fold enhanced affinities compared to compacted or cis-configured conformations. The solvent-exposed and more flexible structures facilitated unique interaction patterns between VP and V_1aR transmembrane helices 3, 4, and 6 which led to high binding energies. The presented “virtual conformational space screening” approach, integrated with NMR spectroscopy, thus enabled identification and characterization of a conformational selection-type complex formation mechanism that confers novel perspectives on targeting the VP-V_1aR interactions at the level of the encounter complex – an aspect that opens novel research avenues for understanding the functionality of the evolutionary selected conformational properties of VP, as well as guidance for ligand design strategies to provide more potent and selective VP analogues.

Elucidating Solution Structures of Cyclic Peptides Using Molecular Dynamics Simulations

Protein-protein interactions are vital to biological processes, but the shape and size of their interfaces make them hard to target using small molecules. Cyclic peptides have shown promise as protein-protein interaction modulators, as they can bind protein surfaces with high affinity and specificity. Dozens of cyclic peptides are already FDA approved, and many more are in various stages of development as immunosuppressants, antibiotics, antivirals, or anticancer drugs. However, most cyclic peptide drugs so far have been natural products or derivatives thereof, with de novo design having proven challenging. A key obstacle is structural characterization: cyclic peptides frequently adopt multiple conformations in solution, which are difficult to resolve using techniques like NMR spectroscopy. The lack of solution structural information prevents a thorough understanding of cyclic peptides' sequence-structure-function relationship. Here we review recent development and application of molecular dynamics simulations with enhanced sampling to studying the solution structures of cyclic peptides. We describe novel computational methods capable of sampling cyclic peptides' conformational space and provide examples of computational studies that relate peptides' sequence and structure to biological activity. We demonstrate that molecular dynamics simulations have grown from an explanatory technique to a full-fledged tool for systematic studies at the forefront of cyclic peptide therapeutic design.

Chiral separation of diastereomers of the cyclic nonapeptides vasopressin and desmopressin by uniform field ion mobility mass spectrometry

In this study ion mobility-mass spectrometry (IM-MS) is used to distinguish chiral diastereomers of the nonapeptides desmopressin and vasopressin. The differences in gas phase cross sectional area (ca. 2%) were sufficient to directly resolve the enantiomers present in a binary mixture. Results from computational modeling indicate that chiral recognition by IM-MS for nonapeptides is possible due to their diastereomer-specific conformations adopted in the gas-phase, namely a compact ring-tail conformer specific to the l-diastereomer forms.

Improved NOE fitting for flexible molecules based on molecular mechanics data – a case study with S-adenosylmethionine

Using the important biomolecule S-adenosyl methionine as an exemplar, we provide a new, enhanced approach for fitting MD data to high-accuracy NOE data, providing improvements in structure determination.