Maximum entropy approach to the determination of solution conformation of flexible polypeptides by global conformational analysis and NMR spectroscopy – Application to DNS1-c-[d-A2bu2, Trp4,Leu5]- enkephalin and DNS1-c-[d-A2bu2, Trp4, d-Leu5]enkephalin

A method is proposed to determine the conformational equilibrium of flexible polypeptides in solution, using the data provided by NMR spectroscopy and theoretical conformational calculations. The algorithm consists of the following three steps: (i) search of the conformational space in order to find conformations with reasonably low energy; (ii) simulation of the NOE spectrum and vicinal coupling constants for each of the low energy conformations; and (iii) determining the statistical weights of the conformations, by means of the maximum-entropy method, in order to obtain the best fit of the averaged NOE intensities and coupling constants to the experimental quantities. The method has been applied to two cyclic enkephalin analogs: DNS1-c-[d-A2bu2,Trp4,Leu5]enkephalin (ENKL) and DNS1-c-[d-A2bu2,Trp4,d-Leu5]enkephalin (ENKD). NMR measurements were carried out in deuterated dimethyl sulfoxide. Two techniques were used in conformational search: the electrostatically driven Monte Carlo method (EDMC), which results in extensive search of the conformational space, but gives only energy minima, and the molecular dynamics method (MD), which results in a more accurate, but also more confined search. In the case of EDMC calculations, conformational energy was evaluated using the ECEPP/3 force field augmented with the SRFOPT solvation-shell model, while in the case of MD the AMBER force field was used with explicit solvent molecules. Both searches and subsequent fitting of conformational weights to NMR data resulted in similar conformations of the cyclic part of the peptides studied. For both ENKL and ENKD a common feature of the low-energy solution conformations is the presence of a type II′ or type IV β-turn at residues 3 and 4; the ECEPP/3 force field also gives a remarkable content of type III β-turn. These β-turns are tighter in the case of ENKL, which is reflected in different distributions of the d-A2bu(NγH)...d-A2bu(CO) and d-A2bu(NγH)...Gly3(CO) hydrogen-bonding distances, indicating that the d-A2bu(NγH) amide proton is more shielded from the solvent than in the case of ENKD. This finding conforms with the results of temperature coefficient data of the d-A2bu(NγH) proton. It has also been found that direct (MD) or Boltzmann (EDMC) averages of the observables do not exactly conform with the measured values, even when explicit solvent molecules are included. This suggests that improving force-field parameters might be necessary in order to obtain reliable conformational ensembles in computer simulations, without the aid of experimental data.

[1]  D. Gorenstein,et al.  On the evaluation of interproton distances for three-dimensional structure determination by NMR using a relaxation rate matrix analysis , 1990 .

[2]  R. Lieu,et al.  A MAXIMUM ENTROPY METHOD WHEN PRIOR INFORMATION CONSISTS OF INEXACT CONSTRAINTS , 1994 .

[3]  Hagai Meirovitch,et al.  DETERMINATION OF THE STABLE MICROSTATES OF A PEPTIDE FROM NOE DISTANCE CONSTRAINTS AND OPTIMIZATION OF ATOMIC SOLVATION PARAMETERS , 1998 .

[4]  K. Wüthrich NMR of proteins and nucleic acids , 1988 .

[5]  Data analysis in frequency-domain fluorometry by the maximum entropy method — recovery of fluorescence lifetime distributions , 1990 .

[6]  J. Stewart Optimization of parameters for semiempirical methods I. Method , 1989 .

[7]  W. F. Gunsteren,et al.  Time-dependent distance restraints in molecular dynamics simulations , 1989 .

[8]  K. Neubert,et al.  Cyclization studies with tetra- and pentapeptide sequences corresponding to beta-casomorphins. , 2009, International journal of peptide and protein research.

[9]  G. Chipens,et al.  Determination of cyclopeptide conformations in solution using NMR data and conformational energy calculations. , 1988, Biophysical chemistry.

[10]  H. Scheraga,et al.  On the multiple‐minima problem in the conformational analysis of polypeptides. II. An electrostatically driven Monte Carlo method—tests on poly(L‐alanine) , 1988, Biopolymers.

[11]  Richard G. Lyon,et al.  A Maximum Entropy Method with a Priori Maximum Likelihood Constraints , 1997 .

[12]  R. Bhaskaran,et al.  Conformational properties of oxytocin in dimethyl sulfoxide solution: NMR and restrained molecular dynamics studies , 1992, Biopolymers.

[13]  T. Yamazaki,et al.  Conformational studies of stereoisomeric 14‐membered cyclic enkephalin analogues containing 1‐naphthyialanine at the fourth position: Chirality effect of leucine at the fifth position on biological activity and receptor selectivity , 1991, Biopolymers.

[14]  Professor Dr. Miklos Bodanszky,et al.  The Practice of Peptide Synthesis , 1994, Springer Lab Manual.

[15]  Combined use of NMR, distance geometry, and restrained energy minimization for the conformational analysis of 8-lysine–vasopressin , 1992 .

[16]  V. Hruby,et al.  Conformations of the dermenkephalin backbone in DMSO solution by a new approach to the solution conformations of flexible peptides , 1993 .

[17]  M. Billeter,et al.  Calibration of the angular dependence of the amide proton-Cα proton coupling constants, 3JHNα, in a globular protein: Use of 3JHNα for identification of helical secondary structure , 1984 .

[18]  P. Schiller,et al.  Conformational studies of diastereomeric cyclic enkephalins by 1H‐NMR and computer simulations , 1987, Biopolymers.

[19]  Olga Kennard,et al.  Systematic analysis of structural data as a research technique in organic chemistry , 1983 .

[20]  P. Schiller,et al.  Enkephalin analogues containing β‐naphthylalanine at the fourth position , 1990 .

[21]  E. Meirovitch,et al.  New theoretical methodology for elucidating the solution structure of peptides from NMR data. II. Free energy of dominant microstates of Leu-enkephalin and population-weighted average nuclear Overhauser effects intensities. , 1996, Biopolymers.

[22]  H. Scheraga,et al.  Chain reversals in proteins. , 1973, Biochimica et biophysica acta.

[23]  E. Meirovitch,et al.  New theoretical methodology for elucidating the solution structure of peptides from NMR data. 1. The relative contribution of low-energy microstates to the partition function , 1995 .

[24]  H. Scheraga,et al.  Empirical solvation models can be used to differentiate native from near‐native conformations of bovine pancreatic trypsin inhibitor , 1991, Proteins.

[25]  G. Fields,et al.  Solid phase peptide synthesis utilizing 9-fluorenylmethoxycarbonyl amino acids. , 2009, International journal of peptide and protein research.

[26]  A. Liwo,et al.  Conformational analysis of [Cpp1, Sar7, Arg8] vasopressin by 1H-NMR spectroscopy and molecular mechanics calculations. , 2009, International journal of peptide and protein research.

[27]  W. V. van Gunsteren,et al.  Time-averaged nuclear Overhauser effect distance restraints applied to tendamistat. , 1990, Journal of molecular biology.

[28]  D. Ripoll,et al.  On the multiple‐minima problem in the conformational analysis of polypeptides. IV. Application of the electrostatically driven monte carlo method to the 20‐residue membrane‐bound portion of melittin , 1990, Biopolymers.

[29]  Arnold Weissberger,et al.  Organic solvents;: Physical properties and methods of purification , 1970 .

[30]  P E Wright,et al.  Distributed torsion angle grid search in high dimensions: A systematic approach to NMR structure determination , 1998, Journal of biomolecular NMR.

[31]  T. Creighton,et al.  Local structure due to an aromatic-amide interaction observed by 1H-nuclear magnetic resonance spectroscopy in peptides related to the N terminus of bovine pancreatic trypsin inhibitor. , 1993, Journal of molecular biology.

[32]  Wilfred F. van Gunsteren,et al.  A Force Field for Liquid Dimethyl Sulfoxide and Physical Properties of Liquid Dimethyl Sulfoxide Calculated Using Molecular Dynamics Simulation , 1995 .

[33]  H A Scheraga,et al.  New developments of the electrostatically driven Monte Carlo method: test on the membrane-bound portion of melittin. , 1998, Biopolymers.

[34]  M. Billeter,et al.  MOLMOL: a program for display and analysis of macromolecular structures. , 1996, Journal of molecular graphics.

[35]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[36]  H. Scheraga,et al.  Energy parameters in polypeptides. 10. Improved geometrical parameters and nonbonded interactions for use in the ECEPP/3 algorithm, with application to proline-containing peptides , 1994 .

[37]  J. Brochon,et al.  Analyzing the distribution of decay constants in pulse-fluorimetry using the maximum entropy method. , 1987, Biophysical journal.

[38]  V J Hruby,et al.  Exploration of the conformational space of oxytocin and arginine-vasopressin using the electrostatically driven Monte Carlo and molecular dynamics methods. , 1998, Biopolymers.

[39]  R Brüschweiler,et al.  Conformational backbone dynamics of the cyclic decapeptide antamanide. Application of a new multiconformational search algorithm based on NMR data. , 1993, Biochemistry.

[40]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[41]  D. Ruppert,et al.  Transformation and Weighting in Regression , 1988 .