Exploration of the conformational space of oxytocin and arginine-vasopressin using the electrostatically driven Monte Carlo and molecular dynamics methods

Design, synthesis, conformational analysis, and biological activity of Cα1-to-Cα6 1,4- and 4,1-disubstituted 1H-[1,2,3]triazol-1-yl-bridged oxytocin analogues

Oxytocin (OT) is a neurohypophyseal peptide hormone containing a disulphide-bridged pseudocyclic conformation. The biomedical use of OT peptides is limited amongst others by disadvantageous pharmacokinetic parameters. To increase the stability of OT by replacing the disulphide bridge with the stable and more rigid [1,2,3]triazol-1-yl moiety, we employed the Cu2+-catalysed side chain-to-side chain azide-alkyne 1,3-cycloaddition. Here we report the design, synthesis, conformational analysis, and in vitro pharmacological activity of a homologous series of Cα1-to-Cα6 side chain-to-side chain [1,2,3]triazol-1-yl-containing OT analogues differing in the length of the bridge, location, and orientation of the linking moiety. Exploiting this macrocyclisation approach, it was possible to generate a systematic series of compounds providing interesting insight into the structure-conformation-function relationship of OT. Most analogues were able to adopt similar conformation to endogenous OT in water, namely, a type I β-turn. This approach may in the future generate stabilised pharmacological peptide tools to advance understanding of OT physiology.

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.

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.

In Silico design of AVP (4-5) peptide and synthesis, characterization and in vitro activity of chitosan nanoparticles

BACKGROUND

Arginine-vasopressin (AVP) is a neuropeptide and provides learning and memory modulation. The AVP (4-5) dipeptide corresponds to the N-terminal fragment of the major vasopressin metabolite AVP (4-9), has a neuroprotective effect and used in the treatment of Alzheimer's and Parkinson's disease.

METHODS

The main objective of the present study is to evaluate the molecular mechanism of AVP (4-5) dipeptide and to develop and synthesize chitosan nanoparticle formulation using modified version of ionic gelation method, to increase drug effectiveness. For peptide loaded chitosan nanoparticles, the synthesized experiment medium was simulated for the first time by molecular dynamics method and used to determine the stability of the peptide, and the binding mechanism to protein (HSP70) was also investigated by molecular docking calculations. A potential pharmacologically features of the peptide was also characterized by ADME (Absorption, Distribution, Metabolism and Excretion) analysis. The characterization, in vitro release study, encapsulation efficiency and loading capacity of the peptide loaded chitosan nanoparticles (CS NPs) were performed by Dynamic Light Scattering (DLS), UV-vis absorption (UV), Scanning Electron Microscopy (SEM), Fourier transform infrared (FT-IR) spectroscopy techniques. Additionally, in vitro cytotoxicity of the peptide on human neuroblastoma cells (SH-SY5Y) was examined with XTT assay and the statistical analysis was evaluated.

RESULTS

The results showed that; hydrodynamic size, zeta potential and polydispersity index (PdI) of the peptide-loaded CS NPs were 167.6 nm, +13.2 mV, and 0.211, respectively. In vitro release study of the peptide-loaded CS NPs showed that 17.23% of the AVP (4-5)-NH2 peptide was released in the first day, while 61.13% of AVP (4-5)-NH2 peptide was released in the end of the 10th day. The encapsulation efficiency and loading capacity were 99% and 10%, respectively. According to the obtained results from XTT assay, toxicity on SHSY-5Y cells in the concentration from 0.01 μg/μL to 30 μg/μL were evaluated and no toxicity was observed. Also, neuroprotective effect was showed against H₂O₂ treatment.

CONCLUSION

The experimental medium of peptide-loaded chitosan nanoparticles was created for the first time with in silico system and the stability of the peptide in this medium was carried out by molecular dynamics studies. The binding sites of the peptide with the HSP70 protein were determined by molecular docking analysis. The size and morphology of the prepared NPs capable of crossing the blood-brain barrier (BBB) were monitored using DLS and SEM analyses, and the encapsulation efficiency and loading capacity were successfully performed with UV Analysis. In vitro release studies and in vitro cytotoxicity analysis on SHSY-5Y cell lines of the peptide were conducted for the first time. Grapical abstract.

Potent antidiuretic agonists, deamino-vasopressin and desmopressin, and their inverso analogs: NMR structure and interactions with micellar and liposomic models of cell membrane

Deamination of vasopressin (AVP) enhances its antidiuretic activity. Moreover, introduction of D-Arg8 instead of its L enantiomer in deamino-vasopressin (dAVP) results in an extremely potent and selective antidiuretic agonist - desmopressin (dDAVP). In this study we describe the synthesis, pharmacological properties and structures of these two potent antidiuretic agonists, and their inverso analogs. The structures of the peptides are studied in micellar and liposomic models of cell membrane using CD spectroscopy. Additionally, three-dimensional structures in mixed anionic-zwitterionic micelles are obtained using NMR spectroscopy supported by molecular dynamics simulations. Our conformational studies have shown that desmopressin in a membrane mimicking environment adopts one of the characteristic for vasopressin-like peptides β-turn - in position 3,4. Furthermore, dDAVP shows the tendency to create a β-turn in the Cys6-Gly9 C-tail, considered to be important for the antidiuretic activity, and also some tendency to adopt a 5,6 β-turn. In desmopressin, in contrast to the native vasopressin, deamino-vasopressin and [D-Arg8]-vasopressin (DAVP), the Arg8 side chain, crucial for the pressor and antidiuretic activities, is very well exposed for interaction with the receptor, whereas Gly9, crucial for the pressor and uterotonic activities, is situated together with the C-terminal amide group very close to the tocin ring. The arrangements of the Gln4 and Asn5 side chains, being crucial for OT activity, also differ in desmopressin as compared to those of AVP, dAVP and DAVP. These differences in arrangement of the important for activities side chains are likely to explain extremely potent and selective antidiuretic activities of desmopressin. © 2016 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 106: 245-259, 2016.

Arginine-, D-arginine-vasopressin, and their inverso analogues in micellar and liposomic models of cell membrane: CD, NMR, and molecular dynamics studies

We describe the synthesis, pharmacological properties, and structures of antidiuretic agonists, arginine vasopressin (AVP) and [d-Arg⁸]-vasopressin (DAVP), and their inverso analogues. The structures of the peptides are studied based on micellar and liposomic models of cell membranes using CD spectroscopy. Additionally, three-dimensional structures in mixed anionic–zwitterionic micelles are obtained using NMR spectroscopy and molecular dynamics simulations. NMR data have shown that AVP and DAVP tend to adopt typical of vasopressin-like peptides β-turns: in the 2–5 and 3–6 fragments. The inverso-analogues also adopt β-turns in the 3–6 fragments. For this reason, their inactivity seems to be due to the difference in side chains orientations of Tyr², Phe³, and Arg⁸, important for interactions with the receptors. Again, the potent antidiuretic activity of DAVP can be explained by CD data suggesting differences in mutual arrangement of the aromatic side chains of Tyr² and Phe³ in this peptide in liposomes rather than of native AVP. In the presence of liposomes, the smallest conformational changes of the peptides are noticed with DPPC and the largest with DPPG liposomes. This suggests that electrostatic interactions are crucial for the peptide–membrane interactions. We obtained similar, probably active, conformations of the antidiuretic agonists in the mixed DPC/SDS micelles (5:1) and in the mixed DPPC/DPPG (7:3) liposomes. Thus it can be speculated that the anionic–zwitterionic liposomes as well as the anionic–zwitterionic micelles, mimicking the eukaryotic cell membrane environment, partially restrict conformational freedom of the peptides and probably induce conformations resembling those of biologically relevant ones.

Conformations of a model cyclic hexapeptide, CYIQNC: (1)H-NMR and molecular dynamics studies

Solution conformation of the cyclic hexapeptide sequence, [cyclo-S-Cys-Tyr-Ile-Gln-Asn-Cys-S] (CYIQNC) - a disulfide-linked fragment of a neurohypophyseal peptide hormone oxytocin (OT) - has been investigated by high-field one-dimensional (1D) and two-dimensional (2D) NMR spectroscopic methods and compared with the results obtained from computer simulation studies. (1)H-NMR results based on temperature dependence of amide proton chemical shifts and nuclear Overhauser effect indicate that peptide in solution populates different conformations, characterized by two fused β-turns. The segment Ile(3)-Gln(4)-Asn(5)-Cys(6) yields a preferred type-III β-turn at residues 4, 5 (HB, 3HN → 6CO), while the segment Cys(6), Cys(1)-Tyr(2)-Ile(3) exhibits inherently weaker, flexible β-turn either of type I/II'/III/half-turn at residues 1, 2 (HB, 6HN → 3CO). The computer simulation studies using a mixed protocol of distance geometry-simulated annealing followed by constrained minimization, restrained molecular dynamics, and energy minimization showed the possibility of existence of additional conformations with the hydrogen bonds, (a) 5HN → 3CO and (b) 2HN → 6CO. These results, therefore, indicate that the additional conformations obtained from both NMR and simulation studies can also be possible to the peptide. These additional conformations might have very small population in the solution and did not show their signatures in these conditions. These findings will be helpful in designing more analogs with modifications in the cyclic moiety of OT.

Constraining cyclic peptides to mimic protein structure motifs

Many proteins exert their biological activities through small exposed surface regions called epitopes that are folded peptides of well-defined three-dimensional structures. Short synthetic peptide sequences corresponding to these bioactive protein surfaces do not form thermodynamically stable protein-like structures in water. However, short peptides can be induced to fold into protein-like bioactive conformations (strands, helices, turns) by cyclization, in conjunction with the use of other molecular constraints, that helps to fine-tune three-dimensional structure. Such constrained cyclic peptides can have protein-like biological activities and potencies, enabling their uses as biological probes and leads to therapeutics, diagnostics and vaccines. This Review highlights examples of cyclic peptides that mimic three-dimensional structures of strand, turn or helical segments of peptides and proteins, and identifies some additional restraints incorporated into natural product cyclic peptides and synthetic macrocyclic peptidomimetics that refine peptide structure and confer biological properties.

Disordered structural ensembles of vasopressin and oxytocin and their mutants

Vasopressin and oxytocin are intrinsically disordered cyclic nonapeptides belonging to a family of neurohypophysial hormones. Although unique in their functions, these peptides differ only by two residues and both feature a tocin ring formed by the disulfide bridge between first and sixth cysteine residues. This sequence and structural similarity are experimentally linked to oxytocin agonism at vasopressin receptors and vasopressin antagonism at oxytocin receptors. Yet single- or double-residue mutations in both peptides have been shown to have drastic impacts on their activities at either receptor, and possibly the ability to bind to their neurophysin carrier protein. In this study we perform molecular dynamics simulations of the unbound native and mutant sequences of the oxytocin and vasopressin hormones to characterize their structural ensembles. We classify the subpopulations of these structural ensembles on the basis of the distributions of radius of gyration and secondary structure and hydrogen-bonding features of the canonical tocin ring and disordered tail region. We then relate the structural changes observed in the unbound form of the different hormone sequences to experimental information about peptide receptor binding, and more indirectly, carrier protein binding affinity, receptor activity, and protease degradation. This study supports the hypothesis that the structural characteristics of the unbound form of an IDP can be used to predict structural or functional preferences of its functional bound form.

Thermodynamics of the interaction between oxytocin and its myometrial receptor in sheep: a stepwise binding mechanism

Entropy (ΔS), enthalpy (ΔH) and heat capacity (ΔCp) changes attending the oxytocin interaction with its two binding sites on myometrial cell membranes in sheep were derived from the temperature dependence of Kd values. The high affinity oxytocin site (Kd on the order of 10(-9)mol l(-1), 25 °C), ascribed to the oxytocin receptor (OXTR), is entropy-driven in the temperature range 0-37 °C. Enthalpy component prevails as a driving force in the binding to the low affinity site (Kd ≈ 10(-7)) within the higher temperature range. ΔCp values in both cases do not differ significantly from zero but become highly relevant in the presence of a GTP analog (10(-4)M GTP-γS). Under these conditions, ΔCp in the low site interaction becomes negative and ΔS is shifted toward negative values (enthalpy drift); ΔCp of the high affinity site rises to a high positive value and the interaction is even more strongly entropy driven. Atosiban, a competitive antagonist of oxytocin at OXTR displays a single significant binding site on myometrial cells (Kd about 10(-7)mol l(-1)). Thermodynamic profiles of atosiban and the low affinity oxytocin site show conspicuous similarities, indicating that the inhibitor is bound to the low affinity site, and not, with a lower affinity, to the putative receptor protein. It is suggested that the interaction of oxytocin with its responding system on myometrial membranes follows in two distinct steps that are likely to be associated with several independent binding domains in the GPCR receptor.

NMR studies of new arginine vasopressin analogs modified with alpha-2-indanylglycine enantiomers at position 2 bound to sodium dodecyl sulfate micelles

In this paper, we use NMR spectroscopy and molecular modeling to examine four new vasopressin analogs modified with alpha-2-indanylglycine (Igl) at position 2, [L-Igl(2)]AVP (I), [D-Igl(2)]AVP (II), [Mpa(1),L-Igl(2)]AVP (III) and [Mpa(1),D-Igl(2)]AVP (IV), embedded in a sodium dodecyl sulfate (SDS) micelle. All the analogs display antiuterotonic activity. In addition, the analogs with D-Igl reveal antipressor properties. Each analog exhibits the tendency to adopt beta-turns at positions 2, 3 and/or 3, 4, which is characteristic of oxytocin-like peptides. Mutual arrangement of aromatic residues at positions 2 and 3 has been found to be crucial for binding antagonists with the OT and V(1a) receptors. The orientation of the Gln(4) side chain seems to be important for the V(1a) receptor affinity. In each of the peptides studied, the Gln(4) side chain is folded back over the ring moiety. However, it lies on the opposite face of the tocin moiety in analogs with L and D enantiomers of Igl.

Structure-antifungal activity relationship of His-Phe-Arg-Trp-Gly-Lys-Pro-Val-NH2 and analogues

The synthesis, in vitro evaluation and conformational study of His-Phe-Arg-Trp-Gly-Lys-Pro-Val-NH(2) and analogues acting as antifungal agents are reported. Among them, His-Phe-Lys-Trp-Gly-Arg-Phe-Val-NH(2) exhibited a moderate but significant antifungal activity against Cryptococcus neoformans, Candida albicans and Candida tropicalis. A theoretical study allows us to propose a biologically relevant conformation for these octapeptides acting as antifungal agents. In addition, these theoretical calculations allow us to determine the minimal structural requirements to produce the antifungal response and can provide a guide for the design of compounds with this biological activity.

A vasopressin/oxytocin-related conopeptide with gamma-carboxyglutamate at position 8

Vasopressins and oxytocins are homologous, ubiquitous and multifunctional peptides present in animals. Conopressins are vasopressin/oxytocin-related peptides that have been found in the venom of cone snails, a genus of marine predatory molluscs that envenom their prey with a complex mixture of neuroactive peptides. In the present paper, we report the purification and characterization of a unique conopressin isolated from the venom of Conus villepinii, a vermivorous cone snail species from the western Atlantic Ocean. This novel peptide, designated gamma-conopressin-vil, has the sequence CLIQDCPgammaG* (gamma is gamma-carboxyglutamate and * is C-terminal amidation). The unique feature of this vasopressin/oxytocin-like peptide is that the eighth residue is gamma-carboxyglutamate instead of a neutral or basic residue; therefore it could not be directly classified into either the vasopressin or the oxytocin peptide families. Nano-NMR spectroscopy of the peptide isolated directly from the cone snails revealed that the native gamma-conopressin-vil undergoes structural changes in the presence of calcium. This suggests that the peptide binds calcium, and the calcium-binding process is mediated by the gamma-carboxyglutamate residue. However, the negatively charged residues in the sequence of gamma-conopressin-vil may mediate calcium binding by a novel mechanism not observed in other peptides of this family.

Ureido group containing cyclic dermorphin(1–7) analogues: synthesis, biology and conformation

Six cyclic peptides related to dermorphin(1-7) have been synthesized. The synthesis of linear peptides containing diamino acid residues in positions 2 and 4 was carried out on a 4-methylbenzhydrylamine resin, and cyclization was achieved by treatment with bis-(4-nitrophenyl)carbonate to form a urea unit. The peptides were tested in the guinea-pig ileum (GPI) and mouse vas deferens (MVD) assays. Diverse opioid agonist activities were observed, depending on the size of the ring. The results were compared with those obtained earlier for 1-4 dermorphin analogues. The conformations of all six dermorphin analogues were studied. The conformational space of the peptides was examined using the electrostatically driven Monte Carlo method. On the basis of NMR data, an ensemble of conformations was obtained for each peptide. The opioid activity profiles of the compounds are discussed in the light of the structural data.

Synthesis and conformational analysis of His-Phe-Arg-Trp-NH2 and analogues with antifungal properties

The synthesis, in vitro evaluation, and conformational study of His-Phe-Arg-Trp-NH2 and related derivatives acting as antifungal agents are reported. Among them, His-Phe-Arg-Trp-NH2 and His-Tyr-Arg-Trp-NH2 exhibited antifungal activity against Cryptococcus neoformans. Antifungal activity of these compounds appears to be closely related to the alpha-MSH effect. A conformational and electronic study allows us to propose a biologically relevant conformation for these tetrapeptides acting as antifungal agents. In addition, these theoretical calculations permit us to determine the minimal structural requirements to produce the antifungal response and may provide a guide for the design of compounds with this biological activity.

Raman and surface-enhanced Raman spectroscopy investigation of vasopressin analogues containing 1-aminocyclohexane-1-carboxylic acid residue

In this work, Raman spectroscopy (RS) was employed to characterize molecular structures of [Arg8]vasopressin (AVP) and its [Acc2,D-Arg8]AVP, [Acc3]AVP, and [Cpa1, Acc3]AVP analogues. The RS band assignments have been proposed. To determine the mechanism of adsorption of the above-mentioned compounds adsorbed on a colloidal silver surface, surface-enhanced Raman spectra (SERS) were measured. The SERS spectra were used to determine relative proximity of the adsorbed functional groups of [corrected] investigated peptides and their orientation on the silver surface. The AVP and [Acc3]AVP SERS spectra (Acc: 1-aminocyclohexane-1-carboxylic acid) show that the L-tyrosine (Tyr) lies far from the metal surface, whereas the [Cpa1,Acc3]AVP spectrum (Cpa: 1-mercaptocyclohexaneacetic acid) provides evidence that Tyr interacts with the silver surface. These results suggest that [corrected] the binding of the Tyr-ionized phenolic group might be responsible for the selectivity of the analogues. We show that the aromatic ring of L-phenylalanine (Phe) of AVP and [Acc2,D-Arg8]AVP interacts with the silver surface. The strength of this interaction is considerably weaker for [Acc2,D-Arg8]AVP than for AVP. This might be due either to a longer distance between the Phe ring and the silver surface, or to the almost perpendicular orientation of the Phe ring towards the surface. The carbonyl group of the L-glutamine [corrected] (Gln) or L-asparagine [corrected](Asn) of AVP, [Acc2,D-Arg8]AVP, and [Acc3]AVP is strongly bound to the silver surface. We have also found that all peptides adsorb on the silver surface via sulfur atoms of the disulfide bridge, adopting a "GGG" conformation, except [Cpa1,Acc3]AVP, which accepts a "TGG" geometry.

Investigation ofcis/trans ratios of peptide bonds in AVP analogues containingN-methylphenylalanine enantiomers

The solution conformation of vasopressin analogues, modified at positions 2 and 3 with N-methylphenylalanine or its enantiomer, [D-MePhe2,MePhe3]AVP and [MePhe2,D-MePhe3]AVP, were studied by 2D NMR spectroscopy in H2O/D2O and theoretical calculations (EDMC/ANALYZE). In the case of [MePhe2,D-MePhe3]AVP, the synthesis afforded two products, A and B, which had identical molecular ions and similar retention times on HPLC. This finding was explained by racemization of Cys1, which gave an additional analogue, [D-Cys1,MePhe2,D-MePhe3]AVP (B). The possibility is not excluded of racemization of Cys1 in the remaining analogues of this series. However, only in the case of [MePhe2,D-MePhe3]AVP was this process so distinct that two strong peaks in the HPLC chromatogram were observed. The NMR spectra of all the analogues showed several distinct sets of residual proton resonances. This suggests that the peptides adopt more than two groups of conformations in H2O/D2O. This fact is due to cis/trans isomerization. Two more populated isomers arise from the cis/trans isomerization across the 2-3 peptide bond in [D-MePhe2,MePhe3]AVP and [MePhe2,D-MePhe3]AVP (A) and across the 1-2 peptide bond in [D-Cys1,MePhe2,D-MePhe3]AVP (B).

Cyclic enkephalin and dermorphin analogues containing a carbonyl bridge

Four cyclic enkephalin analogues and four cyclic dermorphin analogues have been synthesized. Cyclization of linear peptides containing basic amino acid residues of various side chain length in position 2 and 5 (enkephalin analogues) or 2 and 4 (dermorphin analogues) was achieved by treatment with bis-(4-nitrophenyl) carbonate to form a urea unit. The peptides were tested in the guinea-pig ileum (GPI) and mouse vas deferens (MVD) assays. Diverse activity was observed, depending on the size of the ring and the location of the urea unit. The conformation of two dermorphin analogues has been studied: one of high activity (IC(50) = 4.15 nM in the GPI assay) and a second of low activity (IC(50) = 6700 nM in the GPI assay). The conformational space of these peptides was examined using the EDMC method. Using data from the NMR spectra, each peptide was described as an ensemble of conformers. Biological activity was discussed in light of the structural data.

Potent side-chain to side-chain cyclized dermorphin analogues containing a carbonyl bridge

A new family of cyclic opioid peptide analogues related to the 1-4 sequence of dermorphin/deltorphin (Tyr-D-Aaa2-Phe-Aaa4-NH2) has been synthesized. The synthesis of the linear precursor peptides was accomplished by the solid-phase method and ring formation was achieved via a ureido group incorporating the side chain amino functions of D-Aaa2 (D-Lys, D-Orn) and Aaa4 (Lys, Orn, Dab, Dap). The peptides were tested in the guinea-pig ileum (GPI) and mouse vas deferens (MVD) assays. Most showed very high agonist potency in the GPI assay. The peptide containing D-Lys in position 2 and Dab in position 4 was 210 times more active than enkephalin, and that containing Orn and Dab, respectively, was 150 times more active than enkephalin. The latter peptide was also very active in the MVD assay, and showed an IC50 MVD/GPI ratio of 0.816. NMR spectra of selected peptides were recorded, and structural parameters were determined. The conformational space of the peptides was examined using the electrostatically driven Monte Carlo method. With the help of the NMR spectra each peptide was described as an ensemble of conformations. The conformations have been interpreted with regard to the opioid activities, and comparisons have been made with a model proposed earlier for enkephalin analogues.

Molecular dynamics study of a gelsolin-derived peptide binding to a lipid bilayer containing phosphatidylinositol 4,5-bisphosphate

Gelsolin is an actin-severing protein whose action is initiated by Ca(2+) and inhibited by binding to phosphorylated inositol lipid or phosphoinositides. The regions of gelsolin responsible for phosphoinositide binding are comprised of residues 150-169 (G150-169) and 135-142 (G135-142). The corresponding peptides possess similar binding potency as native gelsolin. Their common feature is the presence of arginine and lysine residues that can bind to negatively charged phosphate groups of phosphoinositides. In this work the binding of the G150-169 peptide to a phosphatidylinositol 4,5-bisphosphate (PIP2) cluster in a lipid membrane model was investigated by molecular dynamics calculations (MD) with the AMBER 4.1 force field, taking into account explicit solvent molecules. Initially the structure of G150-169 was simulated by using the electrostatically driven Monte Carlo (EDMC) and MD methods, and the resulting structure agreed within 3.7 A backbone-atom root mean square deviation with the corresponding experimentally derived structure (PDB code: 1SOL). Using this model for the peptide, a subsequent MD simulation of G150-169 in a periodic box containing a model of dimyristoyl-phosphatidylcholine (DMPC) lipids with a cluster of four PIP2 molecules was carried out. During the simulation G150-169 interacted strongly with PIP2 molecules, initially by formation of salt bridges between its N-terminal basic groups and the phosphate groups of PIP2, followed by formation of hydrophobic bonds between the hydrophobic side chains of the peptide and the fatty acid tail of the lipid. As a result of the formation of hydrophobic bonds, the PIP2 molecules were pulled out from the lipid bilayer. This mode of binding differs from those of other PIP2-binding protein motifs such as PH domains that interact solely with the hydrophilic head group of PIP2. These results suggest that dissociation of gelsolin from actin by PIP2 lipids may involve entering of the PIP2 molecules to the gelsolin-actin interface, thereby weakening the interactions between these proteins.

Conformational solution studies of (Sar7)desamino- and (MeAla7)desamino-vasopressin analogues using NMR spectroscopy

Solution structures of two analogues of vasopressin with an amino acid sequence of c[Mpa1-Tyr2-Phe3-Gln4-Asn5-Cys6]-Xaa7-Arg8-Gly9-NH2 (Xaa = Sar [I] or MeAla [II]) were established using ROE and the 3J(HNH alpha) couplings. Each of the peptides was found to exist in two stable isomers, pertaining to the cis or trans status of the Cys6-Xaa7 peptide bond, thus giving rise to four study objects. Two methods for studies of the conformational properties of the structures were compared. In the first, the algorithm consisted of three steps: (i) An Electrostatically Driven Monte-Carlo (EDMC) search for low-energy conformations. (ii) Simulations of the NOESY spectra and the vicinal couplings for these conformations. (iii) Determination of the statistical weights of the conformations with the ANALYZE package, so as to meet the best fit of the averaged NOE intensities and couplings to the experimental data. In the second method, the distance constraints and the torsion angles were used as the usual constraints in the Distance Geometry and Simulated Annealing algorithms. The flexibility of the pressin ring and the C-terminus was characterized by a large number of families of conformations. The presence of the beta-turn at position 4,5 was confirmed for all low energy structures found. The use of the EDMC method for the elaboration of the NMR data for small flexible peptides yielded an adequate description of their conformational diversity and is the method of choice for the analysis of their solution structures.

Conformational studies of cyclic enkephalin analogues with L- or D-proline in position 3

The conformation of a series of cyclic enkephalin analogues of a general formula X(1)-cyclo[Y(2)-Z(3)-Nal(4)-Leu(5)] (Nal: beta-(2-naphthyl)alanine), where X = Tyr, Phe, or Phe(NO(2)), Y = D-Dab or L-Dab (Dab: 2,4-diaminobutyric acid), and Z = D-Pro or L-Pro, was studied by means of NMR spectroscopy and theoretical conformational analysis with the Empirical Conformational Energy Program for Peptides and Proteins force field plus solvation. The NMR measurements were performed in dimethyl sulfoxide solution. The nuclear Overhauser effect intensities and coupling constants were used to compute the statistical weights of the conformations of the ensemble generated in global conformational searches. The purpose of this study was to determine whether introducing the D- or L-proline residue in position 3 can produce peptides with both rigid backbone and significant separation of the pharmacophore groups in position 1 and 4 (as required for high affinity for the mu-type opioid receptors). It was found that the analogues with D-Dab in position 2 and D-Pro in position 3 possess a stable type II' beta-turn at positions 3 and 4, which rigidifies the cyclic backbone; this finding was confirmed by independent measurements of the temperature coefficients of the amide protons, which indicated very significant screening of the Leu(5) amide proton from the solvent. However, these analogues were found to possess a short interchromophore distance. The analogues containing both Dab and Pro in the L-configuration are characterized by a larger interchromophore distance; however, they do not possess a stable beta-turn and have therefore a higher conformational flexibility. The modifications proposed in this work are therefore not likely to lead to enkephalin analogues with a high affinity for the mu-receptors.

Evolution of physics-based methodology for exploring the conformational energy landscape of proteins

The evolution of our physics-based computational methods for determining protein conformation without the introduction of secondary-structure predictions, homology modeling, threading, or fragment coupling is described. Initial use of a hard-sphere potential captured much of the structural properties of polypeptide chains, and subsequent more refined force fields, together with efficient methods of global optimization provide indications that progress is being made toward an understanding of the interresidue interactions that underlie protein folding.

Determination of conformational equilibrium of peptides in solution by NMR spectroscopy and theoretical conformational analysis: application to the calibration of mean-field solvation models

Peptides occur in solution as ensembles of conformations rather than in a fixed conformation. The existing energy functions are usually inadequate to predict the conformational equilibrium in solution, because of failure to account properly for solvation, if the solvent is not considered explicitly (which is usually prohibitively expensive). NMR data are therefore widely incorporated into theoretical conformational analysis. Because of conformational flexibility, restrained molecular dynamics (with restraints derived from NMR data), which is usually applied to determine protein conformation is of limited use in the case of peptides. Instead, (a) the restraints are averaged within predefined time windows during molecular dynamics (MD) simulations (time averaging), (b) multiple-copy MD simulations are carried out and the restraints are averaged over the copies (ensemble averaging), or (c) a representative ensemble of sterically feasible conformations is generated and the weights of the conformations are then fitted so that the computed average observables match the experimental data (weight fitting). All these approaches are briefly discussed in this article. If an adequate force field is used, conformations with large statistical weights obtained from the weight-fitting procedure should also have low energies, which can be implemented in force field calibration. Such a procedure is particularly attractive regarding the parameterization of the solvation energy in nonaqueous solvents, e.g., dimethyl sulfoxide, for which thermodynamic solvation data are scarce. A method for calibration of solvation parameters in dimethyl sulfoxide, which is based on this principle was recently proposed by C. Baysal and H. Meirovitch (Journal of the American Chemical Society, 1998, Vol. 120, pp. 800--812), in which the energy gap between the conformations compatible with NMR data and the alternative conformations is maximized. In this work we propose an alternative method based on the principle that the best-fitting statistical weights of conformations should match the Boltzmann weights computed with the force field applied. Preliminary results obtained using three test peptides of varying conformational mobility: H-Ser(1)-Pro(2)-Lys(3)-Leu(4)-OH, Ac-Tyr(1)-D-Phe(2)-Ser(3)-Pro(4)-Lys(5)-Leu(6)-NH(2), and cyclo(Tyr(1)-D-Phe(2)-Ser(3)-Pro(4)-Lys(5)-Leu(6)) are presented.

Influence of solvents and leucine configuration at position 5 on tryptophan fluorescence in cyclic enkephalin analogues

The fluorescence decay of tryptophan is a sensitive indicator of its local environment within a peptide or protein. In this study we carried out fluorescence measurements of the tryptophan residue of cyclic enkephalin analogues of a general formula X-c[D-Dab(2)-Gly(3)-Trp(4)-Y(5)] where X = Cbz or H and Y = D- or L-Leu, in four solvents [water, methanol, acetonitrile, and dimethyl sulfoxide (DMSO)]. An analysis of the tryptophan fluorescence decays using a discrete-exponential model indicates that tryptophan fluorescence decay can be described by a double exponential function in all solvents studied. Lifetime distribution analysis yields a bimodal distribution in protic solvents (water and methanol), whereas an asymmetric, unimodal distribution in an aprotic solvent (DMSO) and uni- or bimodal distributions in acetonitrile solution, depending on leucine configuration. The data are interpreted in terms of the rotamer model, in which the modality and the relative proportions of the lifetime components are related to the population distribution of tryptophan chi(1) rotamers about the C(alpha)--C(beta) bond. The chirality of the Leu(5) residue and solvent properties affect the local environment of the tryptophan residue and therefore influence the distribution of side-chain rotamers. These results are consistent with the results of theoretical conformational calculations.

Highly potent side-chain to side-chain cyclized enkephalin analogues containing a carbonyl bridge: synthesis, biology and conformation

Six novel cyclic enkephalin analogues have been synthesized. Cyclization of the linear peptides containing basic amino acid residues in position 2 and 5 was achieved by treatment with bis(4-nitrophenyl)carbonate. It was found that some of the compounds exibit unusually high mu-opioid activity in the guinea pig ileum (GPI) assay. The 18-membered analogue cyclo(N(epsilon),N(beta)-carbonyl-D-Lys2,Dap5)-enkephalinamide turned out to be one of the most potent mu-agonists reported so far. NMR spectra of the peptides were recorded and structural parameters were determined. The conformational space was exhaustively examined for each of them using the electrostatically driven Monte Carlo method. Each peptide was finally described as an ensemble of conformations. A model of the bioactive conformation of this class of opioid peptides was proposed.

Towards molecularly imprinted polymers selective to peptides and proteins. The epitope approach

In this paper, we describe the epitope approach to molecular imprinting. The applicability of molecular imprinting, a method that allows the preparation of biomimetic compounds (artificial receptors and antibodies), is extended by this approach. Our approach makes it possible to obtain imprinted polymers selective to peptides and proteins whereas, to date, molecular imprinting has been used primarily for the preparation of polymers that selectively bind to relatively low molecular weight substances. The epitope approach is based on using (as a template) a short peptide that represents only part of a larger peptide or protein (as an epitope represents an antigen), which in turn can be recognized by the synthesized polymer. It is demonstrated that although other parts of peptides can influence the process of molecular recognition, the polymers imprinted with a short peptide efficiently recognize both the template and larger peptides (for example, oxytocin) that possess the same C-terminal part of the structure.

Theoretical conformational analysis of six arginine vasopressin analogs with the L-naphthylalanine in position 3

Introduction of the naphthylalanine residue into either position 3 of arginine vasopressin (AVP), or its analogs results in peptides with interesting pharmacological properties. The single substituted analog of AVP with L-2-Nal in position 3 causes moderate antiduretic activity, whereas [Mpa1, (L-1-Nal)3, (D-Arg)8] VP and [Mpa1, (L-2-Nal)3, (D-Arg)8] VP are potent and selective V2 agonists. Moreover [(L-2-Nal)3, (D-Arg)8] VP is among the most potent and selective antagonists of V1a receptors. In this study we carried out conformational calculations on [(L-1-Nal)3] AVP, [(L-2-Nal)3] AVP, [(L-1-Nal)3, (D-Arg)8] VP, [(L-2-Nal)3, (D-Arg)8] VP, [Mpa1, (L-1-Nal)3, (D-Arg)8] VP, [Mpa1, (L-2-Nal)3, (D-Arg)8] VP, using the ECEPP/3 force field with and without including hydration to simulate aqueous and nonpolar environments. It was found that in all six compound studied, the low-energy conformations have common geometry and relative energies. Therefore, the modifications of the Phe in position 3 influence the binding to the receptor by changing the size of the third residue, rather than by changing the conformational space. The lowest-energy conformations in the presence and absence of water had beta-turns at residues Phe3-Gln4 and Gln4-Asn5 and Gln4-Asn5, respectively. The conformation at the Gln4-Asn5 turn was most similar to the crystal structure of the pressinoic acid (the cyclic moiety of vasopressin).

Interaction of oxytocin with Ca2+: II. Proton magnetic resonance and molecular modeling studies of conformations of the hormone and its Ca2+ complex

Drastic changes in the CD and fluorescence spectra of oxytocin [cyclo(Cys1-Tyr2-Ile3-Gln4-Asn5-Cys6)-Pro7-Leu8-Gly 9-NH2] occur on binding Ca2+ in trifluoroethanol (Ananthanarayanan and Brimble, preceding paper). To further characterize the conformation of the Ca(2+)-bound hormone, we carried out 1H-nmr measurements in deuterated trifluorethanol of oxytocin and its 1:1 Ca2+ complex. The one-dimensional nmr data identified residues involved in Ca2+ binding and the extent of their perturbation on Ca2+ addition. The 3JNH-CH coupling constants and two-dimensional nuclear Overhauser effect (NOE) spectral cross peaks confirmed the helical nature of the Ca2+ complex deduced from CD data. Interproton distances in the free hormone and its Ca2+ complex were estimated from the respective NOE data. Apparent global minimum-energy conformations of free and Ca2+ bound oxytocin were computed using the Monte Carlo with energy minimization protocol, with and without incorporating the NOE-derived distance constraints. Taken together, our results show Ca2+ binding to oxytocin to be a two-step process. The binding of the first Ca2+ brings the otherwise extended tail segment of oxytocin closer to the ring moiety so that it wraps around the cation. This causes the maximal extent of change in all the spectral parameters. The subsequent formation of the 2:1 Ca-oxytocin complex results in the tail detaching itself away from the ring so as to bind the second Ca2+ ion. This leads to further spectral changes in the hormone molecule. The tail segment plays a major role in both steps. These observations may be useful in understanding the structural basis of oxytocin action.

Fluorescence study of neurohypophyseal hormones and their analogues. Distance distributions in a series of arginine-vasopressin analogues

Analogues of arginine-vasopressin (AVP) in which substitution of the proline residue in position 7 (by either sarcosine or N-methylalanine) combined with replacement of the cysteine residue in position 1 were the subject of a fluorescence and molecular mechanics study. We obtained two groups of analogues: selective antidiuretic agonists (cysteine or beta-mercaptopropionic acid in position 1) and pressor and uterotonic antagonists (deaminopenicillamine or beta-mercapto-beta, beta-cyclopentamethylenepropionic acid in position 1). Using frequency-domain measurements of fluorescence resonance energy transfer (FRET) we estimated the distance distribution between the phenolic ring of Tyr2 and the disulphide bridge Cys1-Cys6. We also analyzed acrylamide quenching of tyrosyl fluorescence to determine the exposure of the tyrosyl ring to the solvent. Results from fluorescence experiments were compared with those from Monte Carlo simulation (ECEPP/3 force-field).