Conformation and Dynamics of Human Urotensin II and Urotensin Related Peptide in Aqueous 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.

Urotensin II Receptor Modulation with 1,3,4-Benzotriazepin-2-one Tetrapeptide Mimics

Urotensin II receptor (UT) modulators that differentiate the effects of the endogenous cyclic peptide ligands urotensin II (UII) and urotensin II-related peptide (URP) offer potential for dissecting their respective biological roles in disease etiology. Selective modulators of hUII and URP activities were obtained using 1,3,4-benzotriazepin-2-one mimics of a purported bioactive γ-turn conformation about the Bip-Lys-Tyr tripeptide sequence of urocontrin ([Bip4]URP). Considering an active β-turn conformer about the shared Phe-Trp-Lys-Tyr sequence of UII and URP, 8-substituted 1,3,4-benzotriazepin-2-ones were designed to mimic the Phe-Bip-Lys-Tyr tetrapeptide sequence of urocontrin, synthesized, and examined for biological activity. Subtle 5- and 8-position modifications resulted in biased signaling and selective modulation of hUII- or URP-induced vasoconstriction. For example, p-hydroxyphenethyl analogs 17b-d were strong Gα13 and βarr1 activators devoid of Gαq-mediated signaling. Tertiary amides 15d and 17d negatively modulated hUII-induced vasoconstriction without affecting URP-mediated responses. Benzotriazepinone carboxamides proved to be exceptional tools for elucidating the pharmacological complexity of UT.

1,3,5,8-Tetrasubstituted 1,3,4-Benzotriazepin-2-one Scaffolds for β-Turn Mimicry without Stereogenic Carbon Centers: Synthesis and Conformational Analysis

Topological mimicry of peptide β-turn secondary structures has been investigated using a 1,3,5,8-tetrasubstituted 1,3,4-benzotriazepin-2-one scaffold. Approaches were conceived for the synthesis of tetrasubstituted benzotriazepinones from 4-acetyl-3-aminobenzoate based on aza-amino acid chemistry and different orthogonal protection strategies. Installation of an 8-position carboxylate on the aromatic ring enabled a diverse array of substituents to be introduced for mimicry of the i-position residue. Benzotriazepin-2-one crystallization and X-ray analysis demonstrated that in spite the absence of a stereogenic carbon center, the scaffold could serve as type I and I' β-turn mimics, because pyramidalization of the N³-nitrogen in the benzotriazepin-2-one provides potential for adoptive chirality. 1,3,5,8-Tetrasubstituted 1,3,4-benzotriazepin-2-one scaffolds offer interesting potential for the cost-effective synthesis of nonpeptide β-turn surrogates for peptide mimicry in various recognition events.

Higher Affinity Antibodies Bind With Lower Hydration and Flexibility in Large Scale Simulations

We have carried out a long-timescale simulation study on crystal structures of nine antibody-antigen pairs, in antigen-bound and antibody-only forms, using molecular dynamics with enhanced sampling and an explicit water model to explore interface conformation and hydration. By combining atomic level simulation and replica exchange to enable full protein flexibility, we find significant numbers of bridging water molecules at the antibody-antigen interface. Additionally, a higher proportion of interactions excluding bulk waters and a lower degree of antigen bound CDR conformational sampling are correlated with higher antibody affinity. The CDR sampling supports enthalpically driven antibody binding, as opposed to entropically driven, in that the difference between antigen bound and unbound conformations do not correlate with affinity. We thus propose that interactions with waters and CDR sampling are aspects of the interface that may moderate antibody-antigen binding, and that explicit hydration and CDR flexibility should be considered to improve antibody affinity prediction and computational design workflows.

How well does molecular simulation reproduce environment-specific conformations of the intrinsically disordered peptides PLP, TP2 and ONEG?

Understanding the conformational ensembles of intrinsically disordered proteins and peptides (IDPs) in their various biological environments is essential for understanding their mechanisms and functional roles in the proteome, leading to a greater knowledge of, and potential treatments for, a broad range of diseases. To determine whether molecular simulation is able to generate accurate conformational ensembles of IDPs, we explore the structural landscape of the PLP peptide (an intrinsically disordered region of the proteolipid membrane protein) in aqueous and membrane-mimicking solvents, using replica exchange with solute scaling (REST2), and examine the ability of four force fields (ff14SB, ff14IDPSFF, CHARMM36 and CHARMM36m) to reproduce literature circular dichroism (CD) data. Results from variable temperature (VT) 1H and Rotating frame Overhauser Effect SpectroscopY (ROESY) nuclear magnetic resonance (NMR) experiments are also presented and are consistent with the structural observations obtained from the simulations and CD. We also apply the optimum simulation protocol to TP2 and ONEG (a cell-penetrating peptide (CPP) and a negative control peptide, respectively) to gain insight into the structural differences that may account for the observed difference in their membrane-penetrating abilities. Of the tested force fields, we find that CHARMM36 and CHARMM36m are best suited to the study of IDPs, and accurately predict a disordered to helical conformational transition of the PLP peptide accompanying the change from aqueous to membrane-mimicking solvents. We also identify an α-helical structure of TP2 in the membrane-mimicking solvents and provide a discussion of the mechanistic implications of this observation with reference to the previous literature on the peptide. From these results, we recommend the use of CHARMM36m with the REST2 protocol for the study of environment-specific IDP conformations. We believe that the simulation protocol will allow the study of a broad range of IDPs that undergo conformational transitions in different biological environments.

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.

Insights into the Molecular Determinants Involved in Urocontrin and Urocontrin A Action

In the past few years, we have identified two allosteric modulators of the urotensinergic system with probe-dependent action, termed Urocontrin (UC) and Urocontrin A (UCA). Such action is atypical and important since it will allow us to understand the specific function of the functionally selective cognate ligands of this system, namely urotensin II and urotensin II-related peptide. Delineating the molecular determinants involved in this particular behavior would represent an important step toward designing small molecules suitable for pharmacologic and/or therapeutic intervention. Hence, we undertook an exploratory research by replacing the Trp⁴ residue of URP with several para-substituted phenylalanine amino acids in order to get a grasp on the required nature, distance, and orientation of the side chain of this residue for allosteric modulatory action. We found that the position of the second aromatic group is crucial, and we identified two new allosteric modulators: [Trip⁴]URP and [Phe(pPy-4)⁴]URP with probe-dependent action.

Total flavones of Rhododendron simsii Planch flower protect isolated rat heart from ischaemia-reperfusion injury and its mechanism of UTR-RhoA-ROCK pathway inhibition

Objectives

Total flavones of Rhododendron simsii Planch flower (TFR) are an effective part extracted from the flower. The present study was designed to investigate the protective effect of TFR in isolated rat heart following global ischaemia-reperfusion and the possible underlying mechanisms.

Methods

Langendorff perfusion apparatus was used to perfuse isolated rat heart which was subjected to global ischaemia-reperfusion. The hemodynamic parameters were continuously monitored. Coronary flow as well as lactate dehydrogenase (LDH), creatine phosphokinase-MB (CK-MB) and cardiac troponin I (cTnI) in coronary effluents was measured. RhoA activity and urotensin receptor (UTR) and Rho-related coiled-coil-forming protein kinase (ROCK) protein expressions in rat myocardium were examined, respectively. Cardiac dysfunction was indicated by the alterations of hemodynamic parameters and the reduced coronary flow.

Key findings

Total flavones of Rhododendron simsii Planch flower significantly improved ischaemia-reperfusion–induced cardiac dysfunction and leakages of LDH, CK-MB and cTnI, and inhibited myocardial ischaemia-reperfusion–increased RhoA activity and UTR, ROCK1 and ROCK2 protein expressions. The improvement of TFR in the cardiac dysfunction and the leakage of LDH, CK-MB and cTnI were markedly attenuated under the UTR blockade and ROCK inhibition. TFR-inhibited RhoA activity was decreased under the UTR blockade.

Conclusions

Total flavones of Rhododendron simsii Planch flower had a protective effect on ischaemia-reperfusion injury in isolated rat heart, which may be attributed to the blocking of UTR and subsequent inhibition of the RhoA-ROCK pathway.