Conformational properties of oxytocin in dimethyl sulfoxide solution: NMR and restrained molecular dynamics studies

Structure-Based Design of Glycosylated Oxytocin Analogues with Improved Selectivity and Antinociceptive Activity

Acute and chronic pain is often treated with opioids despite the negative side effects of constipation, physical dependence, respiratory depression, and overdose. The misuse of opioid analgesics has given rise to the opioid crisis/epidemic, and alternate nonaddictive analgesics are urgently needed. Oxytocin, a pituitary hormone, is an alternative to the small molecule treatments available and has been used as an analgesic as well as for the treatment and prevention of opioid use disorder (OUD). Clinical implementation is limited by its poor pharmacokinetic profile, a result of the labile disulfide bond between two cysteine residues in the native sequence. Stable brain penetrant oxytocin analogues have been synthesized by replacement of the disulfide bond with a stable lactam and glycosidation of the C-terminus. These analogues show exquisite selectivity for the oxytocin receptor and potent in vivo antinociception in mice following peripheral (i.v.) administration, supporting further study of their clinical potential.

Cysteines and Disulfide Bonds as Structure-Forming Units: Insights From Different Domains of Life and the Potential for Characterization by NMR

Disulfide bridges establish a fundamental element in the molecular architecture of proteins and peptides which are involved e.g., in basic biological processes or acting as toxins. NMR spectroscopy is one method to characterize the structure of bioactive compounds including cystine-containing molecules. Although the disulfide bridge itself is invisible in NMR, constraints obtained via the neighboring NMR-active nuclei allow to define the underlying conformation and thereby to resolve their functional background. In this mini-review we present shortly the impact of cysteine and disulfide bonds in the proteasome from different domains of life and give a condensed overview of recent NMR applications for the characterization of disulfide-bond containing biomolecules including advantages and limitations of the different approaches.

Systematic N-methylation of oxytocin: Impact on pharmacology and intramolecular hydrogen bonding network

Oxytocin (OT) is a peptide hormone agonist of the OT receptor (OTR) that plays an important role in social behaviors such as pair bonding, maternal bonding and trust. The pharmaceutical development of OT as an oral peptide therapeutic has been hindered historically by its unfavorable physicochemical properties, including molecular weight, polarity and number of hydrogen bond donors, which determines poor cell permeability. Here we describe the first systematic study of single and multiple N-methylations of OT and their effect on physicochemical properties as well as potency at the OT receptor. The agonist EC50 and percent effect for OTR are reported and show that most N-methylations are tolerated but with some loss in potency compared to OT. The effect of N-methylation on exposed polarity is assessed through the EPSA chromatographic method and the results validated against NMR temperature coefficient experiments and the determination of NMR solution structures. We found that backbone methylation of residues not involved in IMHB and removal of the N-terminal amine can significantly reduce the exposed polarity of peptides, and yet retain a significant OTR agonist activity. The results of this study also expose the potential challenge of using the N-methylation strategy for the OT system; while exposed polarity is reduced, in some cases backbone methylation produces a significant conformational change that compromises agonist activity. The data presented provides useful insights on the SAR of OT and suggests future design strategies that can be used to develop more permeable OTR agonists based on the OT framework.

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.

Ligand-based peptide design and combinatorial peptide libraries to target G protein-coupled receptors

G protein-coupled receptors (GPCRs) are considered to represent the most promising drug targets; it has been repeatedly said that a large fraction of the currently marketed drugs elicit their actions by binding to GPCRs (with cited numbers varying from 30-50%). Closer scrutiny, however, shows that only a modest fraction of (≈60) GPCRs are, in fact, exploited as drug targets, only ≈20 of which are peptide-binding receptors. The vast majority of receptors in the humane genome have not yet been explored as sites of action for drugs. Given the drugability of this receptor class, it appears that opportunities for drug discovery abound. In addition, GPCRs provide for binding sites other than the ligand binding sites (referred to as the "orthosteric site"). These additional sites include (i) binding sites for ligands (referred to as "allosteric ligands") that modulate the affinity and efficacy of orthosteric ligands, (ii) the interaction surface that recruits G proteins and arrestins, (iii) the interaction sites of additional proteins (GIPs, GPCR interacting proteins that regulate G protein signaling or give rise to G protein-independent signals). These sites can also be targeted by peptides. Combinatorial and natural peptide libraries are therefore likely to play a major role in identifying new GPCR ligands at each of these sites. In particular the diverse natural peptide libraries such as the venom peptides from marine cone-snails and plant cyclotides have been established as a rich source of drug leads. High-throughput screening and combinatorial chemistry approaches allow for progressing from these starting points to potential drug candidates. This will be illustrated by focusing on the ligand-based drug design of oxytocin (OT) and vasopressin (AVP) receptor ligands using natural peptide leads as starting points.

Complete NMR analysis of oxytocin in phosphate buffer

Complete NMR analysis of oxytocin (OXT) in phosphate buffer was elucidated by one-dimensional (1D)- and two-dimensional (2D)-NMR techniques, which involve the assignment of peptide amide NH protons and carbamoyl NH(2) protons. The (1)H-(15)N correlation of seven amide NH protons and three carbamoyl NH(2) protons were also shown by HSQC NMR of OXT without (15)N enrichment.

Phage-encoded combinatorial chemical libraries based on bicyclic peptides

Here we describe a phage strategy for the selection of ligands based on bicyclic or linear peptides attached covalently to an organic core. We designed peptide repertoires with three reactive cysteine residues, each spaced apart by several random amino acid residues, and we fused the repertoires to the phage gene-3-protein. Conjugation with tris-(bromomethyl)benzene via the reactive cysteines generated repertoires of peptide conjugates with two peptide loops anchored to a mesitylene core. Iterative affinity selections yielded several enzyme inhibitors; after further mutagenesis and selection, we were able to chemically synthesize a lead inhibitor (PK15; Ki =1.5 nM) specific to human plasma kallikrein that efficiently interrupted the intrinsic coagulation pathway in human plasma tested ex vivo. This approach offers a powerful means of generating and selecting bicyclic macrocycles (or if cleaved, linear derivatives thereof) as ligands poised at the interface of small-molecule drugs and biologics.

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.

Examination of structural characteristics of the potent oxytocin antagonists [dPen1,Pen6]-OT and [dPen1,Pen6, 5-tBuPro7]-OT by NMR, Raman, CD spectroscopy and molecular modeling

The synthesis and biological evaluation of penicillamine(6)-5-tert-butylproline(7)-oxytocin analogs and comparison with their proline(7)-oxytocin counterparts has led to the discovery of two potent oxytocin (OT) antagonists: [dPen(1),Pen(6)]-oxytocin (1, pA(2) = 8.22, EC(50) = 6.0 nM) and [dPen(1),Pen(6),5-tBuPro(7)]-oxytocin (2, pA(2) = 8.19, EC(50) = 6.5 nM). In an attempt to understand the conformational requirements for their biological activity, spectroscopic analyses of 1 and 2 were performed using (1)H NMR, laser Raman and CD techniques. In H(2)O, oxytocin analogs 1 and 2 exhibited cis-isomer populations of 7% and 35%, respectively. Measurement of the amide proton temperature coefficients revealed solvent shielded hydrogens for Gln(4) and Pen(6) in the major trans-conformer of 1 as well as for Gln(4) in the minor cis-conformer of 2. Few long-distance NOEs were observed, suggesting conformational averaging for analogs 1 and 2 in water; moreover, a lower barrier (16.6 +/- 0.2 kcal/mol) for isomerization of the amide N-terminal to 5-tBuPro(7) relative to OT was calculated from measuring the coalescence temperature of the Gly(9) backbone NH signals in the NMR spectra of 2. Observed bands in the Raman spectra of 1 and 2 correspond to C(beta)-S-S-C(beta) dihedral angles of +110-115 degrees and +/-90 degrees , respectively. In water, acetonitrile and methanol, the CD spectra for 1 exhibited a positive maximum around 236-239 nm; in trifluoroethanol, the spectra shifted and a negative maximum was observed at 240 nm. The CD spectra of 2 were unaffected by solvent changes and exhibited a negative maximum at 236-239 nm. The CD and Raman data both suggested that a conformation having a right-handed screw sense about the disulfide and a chi(CS-SC) dihedral angle value close to 115 degrees was favored for analog 1 in water, methanol and acetonitrile, but not trifluoroethanol, where a +/-90 degrees angle was favored. Analog 2 was more resilient to conformational change about the disulfide, and adopted a preferred disulfide geometry corresponding to a +/-90 degrees chi(CS-SC) dihedral angle. Monte Carlo conformational analysis of analogs 1 and 2 using distance restraints derived from NMR spectroscopy revealed two prominent conformational minima for analog 1 with disulfide geometries around +114 degrees and +116 degrees . Similar analysis of analog 2 revealed one conformational minimum with a disulfide geometry around +104 degrees . In sum, the conformation about the disulfide in [dPen(1),Pen(6)]-OT (1) was shown to be contingent on environment and in TFE, adopted a geometry similar to that of [dPen(1),Pen(6),5-tBuPro(7)]-OT (2) which appeared to be stabilized by hydrophobic interactions between the 5-tBuPro(7) (5R)-tert-butyl group, the Leu(8) isopropyl sidechain and the Pen(6)beta-methyl substituents. In light of the conformational rigidity of 2 about the disulfide bond, and the similar geometry adopted by 1 in TFE, a S-S dihedral angle close to +110 degrees may be a prerequisite for their binding at the receptor.

Synthesis and utilization of 13C and 15N backbone-labeled proline: NMR study of synthesized oxytocin with backbone-labeled C-terminal tripeptide amide

The 13C and 15N backbone-labeled proline was prepared using Oppolzer's method based on application of a sultam as chiral auxiliary. This isotopomer was used in the synthesis of the 13C, 15N backbone-labeled C-terminal tripeptide amide fragment of neurohypophyseal hormone oxytocin. Finally, this tripeptide amide was coupled by segment condensation with N-Boc- or N-Fmoc-tocinoic acid, followed by N-deprotection with TFA or piperidine. The labeled oxytocin exhibited biological activity identical with that of natural oxytocin. A detailed 1H, 13C and 15N NMR study confirmed the assigned oxytocin conformation containing a beta-turn in the cyclic part of the molecule, stabilized by H-bond(s) that can be perturbed by the C-terminal tripeptide amide moiety as indicated by comparison of NMR data for both the tocine ring in oxytocin and tocinoic acid.

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.

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.

Solution conformational study of Scyliorhinin I analogues with conformational constraints by two-dimensional NMR and theoretical conformational analysis

Two analogues of Scyliorhinin I (Scyl), a tachykinin with N-MeLeu in position 8 and a 1,5-disubstituted tetrazole ring between positions 7 and 8, introduced in order to generate local conformational constraints, were synthesized using the solid-phase method. Conformational studies in water and DMSO-d6 were performed on these peptides using a combination of the two-dimensional NMR technique and theoretical conformational analysis. The algorithm of conformational search consisted of the following three stages: (i) extensive global conformational analysis in order to find all low-energy conformations; (ii) calculation of the NOE effects and vicinal coupling constants for each of the low energy conformations; (iii) determining the statistical weights of these conformations by means of a nonlinear least-squares procedure, in order to obtain the best fit of the averaged simulated spectrum to the experimental one. In both solvents the three-dimensional structure of the analogues studied can be interpreted only in terms of an ensemble of multiple conformations. For [MeLeu8]Scyl, the C-terminal 6-10 fragment adopts more rigid structure than the N-terminal one. In the case of the analogue with the tetrazole ring in DMSO-d6 the three-dimensional structure is characterized by two dominant conformers with similar geometry of their backbones. They superimpose especially well (RMSD = 0.28 A) in the 6-9 fragments. All conformers calculated in both solvents superimpose in their C-terminal fragments much better than those of the first analogue. The results obtained indicate that the introduction of the tetrazole ring into the Scyl molecule rigidifies its structure significantly more than that of MeLeu.

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.

Nuclear magnetic resonance studies of the C-terminal human growth hormone fragment I179-C182-[SS]-C189-P191 and the related trisulfide peptide I179-C182-[SSS]-C189-P191

The synthetic C-terminal hGH fragment I179-C182-[SS]-C189-P191 and the related trisulfide peptide I179-C182-[SSS]-C189-P191 have been studied using homonuclear 1H-NMR methods and distance geometry calculations. The 1H-NMR spectra of both the disulfide (diS) and the trisulfide (triS) were completely assigned. Amide proton exchange rates, NOEs and the temperature dependence of the NH chemical shifts indicate a hydrogen bond in triS between Val185 and Ser188 stabilizing a turn in this region. 3JH,H coupling constants and NOEs were measured and used as input for distance geometry calculations. For triS two families of structures with averaged pairwise backbone root mean square deviations for Cys182-Cys189 of 1.3-1.5 A were found, only one of which is compatible with experimental data. For diS only one family of structures was found, but with such a low structural definition (back bone rmsd > 2 A) that no interpretation into a consensus structure is useful. The generated structures were compared to the crystal structure of the terminal loop in hGH, complexed to its binding proteins. The resemblance was low between the solution structures of the tridecapeptides and the terminal hGH loop.

NMR conformational studies on a synthetic peptide reproducing the [1-20] processing domain of the pro-ocytocin-neurophysin precursor

The combined use of several nuclear magnetic resonance and restrained molecular dynamics techniques allowed the formulation of a molecular model for the preferred solution conformation of a synthetic peptide reproducing the [1-20] processing domain of the pro-ocytocin-neurophysin precursor. In the model, the conformation of the 20-membered tocin ring, with the two Cys1 and Cys6 residues bridged by a disulphide bond, is very close to that observed for isolated ocytocin in the solid state; in addition, a type II beta-turn is postulated for the 7-10 segment of the acyclic tail containing the Lys11-Arg12 processing site, and connecting ocytocin to the neurophysin domain, while the C-terminal 13-20 segment of the molecule is believed to assume a helical structure. This particular structural organization could be important in participating as the favorable conformation for optimal substrate-enzyme active site recognition and processing by specific endoproteases.

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

Conformational analysis of the neurohypophyseal hormones oxytocin (OT) and arginine-vasopressin (AVP) has been carried out using two different computational approaches and three force fields, namely by the Electrostatically Driven Monte Carlo (EDMC) method, with the Empirical Conformational Energy Program for Peptides (ECEPP/3) force field or with the ECEPP/3 force field plus a hydration-shell model, and by simulated-annealing molecular dynamics with the Consistent Valence Force Field (CVFF). The low-energy conformations obtained for both hormones were classified using the minimal-tree clustering algorithm and characterized according to the locations of beta-turns in the cyclic moieties. Calculations with the CVFF force field located conformations with a beta-turn at residues 3 and 4 as the lowest energy ones both for OT and for AVP. In the ECEPP/3 force field the lowest energy conformation of OT contained a beta-turn at residues 2 and 3, conformations with this location of the turn being higher in energy for AVP. The latter difference can be attributed to the difference in the size of the side chain in position 3 of the sequences: the bulkier phenylalanine residue of AVP in combination with the bulky Tyr2 residue hinders the formation of a turn at residues 2 and 3. Conformations of OT and AVP with a turn at residues 3,4 were in the best agreement with the x-ray structures of deaminooxytocin and pressinoic acid (the cyclic moiety of vasopressin), respectively, and with the nmr-derived distance constraints. Generally, the low-energy conformations obtained with the hydration-shell model were in a better agreement with the experimental data than the conformations calculated in vacuo. It was found, however, that the obtained low-energy conformations do not satisfy all of the nmr-derived distance constraints and the nuclear Overhauser effect pattern observed in nmr studies can be fully explained only by assuming a dynamic equilibrium between conformations with beta-turns at residues 2,3, 3,4, and 4,5. The low-energy structures of OT with a beta-turn at residues 2,3 have the disulfide ring conformations close to the model proposed recently for a potent bicyclic antagonist of OT [M. D. Shenderovich et al. (1994) Polish Journal of Chemistry, Vol. 25, pp. 921-927], although the native hormone differs from the bicyclic analogue by the conformation of the C-terminal tripeptide. This finding confirms the hypothesis of different receptor-bound conformations of agonists and antagonists of OT.