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

Adam Liwo, Anna Tempczyk, Stanislaw Oldziej, Mark D. Shenderovich, Victor J Hruby, Sekhar Talluri, Jerzy Ciarkowski, Franciszek Kasprzykowski, Leszek Lankiewicz, Zbigniew Grzonka

Research output: Contribution to journalArticle

56 Citations (Scopus)

Abstract

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 β-turns in the cyclic moieties. Calculations with the CVFF force field located conformations with a β-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 β-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 β-turns at residues 2,3, 3,4, and 4,5. The low-energy structures of OT with a β-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.

Original languageEnglish (US)
Pages (from-to)157-175
Number of pages19
JournalBiopolymers - Peptide Science Section
Volume38
Issue number2
StatePublished - 1996

Fingerprint

Space Flight
Arginine
Arginine Vasopressin
Oxytocin
Molecular Dynamics Simulation
Conformations
Molecular dynamics
Hormones
Posterior Pituitary Hormones
Monte Carlo Method
Hydration
Vasopressins
Phenylalanine
Disulfides
Cluster Analysis
X-Rays
Peptides
Simulated annealing
Clustering algorithms
Monte Carlo methods

ASJC Scopus subject areas

  • Biochemistry, Genetics and Molecular Biology(all)
  • Biochemistry
  • Biophysics

Cite this

Exploration of the conformational space of oxytocin and arginine-vasopressin using the electrostatically driven monte carlo and molecular dynamics methods. / Liwo, Adam; Tempczyk, Anna; Oldziej, Stanislaw; Shenderovich, Mark D.; Hruby, Victor J; Talluri, Sekhar; Ciarkowski, Jerzy; Kasprzykowski, Franciszek; Lankiewicz, Leszek; Grzonka, Zbigniew.

In: Biopolymers - Peptide Science Section, Vol. 38, No. 2, 1996, p. 157-175.

Research output: Contribution to journalArticle

Liwo, A, Tempczyk, A, Oldziej, S, Shenderovich, MD, Hruby, VJ, Talluri, S, Ciarkowski, J, Kasprzykowski, F, Lankiewicz, L & Grzonka, Z 1996, 'Exploration of the conformational space of oxytocin and arginine-vasopressin using the electrostatically driven monte carlo and molecular dynamics methods', Biopolymers - Peptide Science Section, vol. 38, no. 2, pp. 157-175.
Liwo, Adam ; Tempczyk, Anna ; Oldziej, Stanislaw ; Shenderovich, Mark D. ; Hruby, Victor J ; Talluri, Sekhar ; Ciarkowski, Jerzy ; Kasprzykowski, Franciszek ; Lankiewicz, Leszek ; Grzonka, Zbigniew. / Exploration of the conformational space of oxytocin and arginine-vasopressin using the electrostatically driven monte carlo and molecular dynamics methods. In: Biopolymers - Peptide Science Section. 1996 ; Vol. 38, No. 2. pp. 157-175.
@article{847df87ce82341fbb24de5ef6513f739,
title = "Exploration of the conformational space of oxytocin and arginine-vasopressin using the electrostatically driven monte carlo and molecular dynamics methods",
abstract = "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 β-turns in the cyclic moieties. Calculations with the CVFF force field located conformations with a β-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 β-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 β-turns at residues 2,3, 3,4, and 4,5. The low-energy structures of OT with a β-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.",
author = "Adam Liwo and Anna Tempczyk and Stanislaw Oldziej and Shenderovich, {Mark D.} and Hruby, {Victor J} and Sekhar Talluri and Jerzy Ciarkowski and Franciszek Kasprzykowski and Leszek Lankiewicz and Zbigniew Grzonka",
year = "1996",
language = "English (US)",
volume = "38",
pages = "157--175",
journal = "Biopolymers",
issn = "0006-3525",
publisher = "John Wiley and Sons Inc.",
number = "2",

}

TY - JOUR

T1 - Exploration of the conformational space of oxytocin and arginine-vasopressin using the electrostatically driven monte carlo and molecular dynamics methods

AU - Liwo, Adam

AU - Tempczyk, Anna

AU - Oldziej, Stanislaw

AU - Shenderovich, Mark D.

AU - Hruby, Victor J

AU - Talluri, Sekhar

AU - Ciarkowski, Jerzy

AU - Kasprzykowski, Franciszek

AU - Lankiewicz, Leszek

AU - Grzonka, Zbigniew

PY - 1996

Y1 - 1996

N2 - 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 β-turns in the cyclic moieties. Calculations with the CVFF force field located conformations with a β-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 β-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 β-turns at residues 2,3, 3,4, and 4,5. The low-energy structures of OT with a β-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.

AB - 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 β-turns in the cyclic moieties. Calculations with the CVFF force field located conformations with a β-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 β-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 β-turns at residues 2,3, 3,4, and 4,5. The low-energy structures of OT with a β-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.

UR - http://www.scopus.com/inward/record.url?scp=0030051535&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0030051535&partnerID=8YFLogxK

M3 - Article

VL - 38

SP - 157

EP - 175

JO - Biopolymers

JF - Biopolymers

SN - 0006-3525

IS - 2

ER -