Premium
Predicting the conformational states of cyclic tetrapeptides
Author(s) -
Loiseau Nicolas,
Gomis JeanMarie,
Santolini Jérôme,
Delaforge Marcel,
André François
Publication year - 2003
Publication title -
biopolymers
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.556
H-Index - 125
eISSN - 1097-0282
pISSN - 0006-3525
DOI - 10.1002/bip.10339
Subject(s) - conformational isomerism , chemistry , computational chemistry , molecular mechanics , cyclic peptide , peptide , molecule , stereochemistry , molecular dynamics , organic chemistry , biochemistry
Abstract Biologically active cyclic tetrapeptides, usually found among fungi metabolites, exhibit phytotoxic or cytostatic activities that are likely to be governed by specific conformations adopted in solution. For conformational studies and drug design, there is a strong interest in using fast and reliable methods to determine correctly the conformational population of cyclotetrapeptides. We show here that standard molecular mechanics computational approach gives satisfactory results. The method was validated step by step by experimental data either obtained after synthesis and NMR analysis, or found in the literature. The cyclo(Gly) 4 , cyclo(Ala) 4 , cyclo(Sar) 4 , and cyclo(SarGly) 2 peptides were used to evaluate the prediction of the peptide backbone conformation, and the detailed conformational analysis of tentoxin, a natural phytotoxic cyclotetrapeptide in which N‐alkylated peptide bonds alternate with regular secondary ones, was used to validate the computation of conformers proportions. From the knowledge of an initial cyclic primary structure and of the D or L configuration of the amino acids, we show that it is possible to determine the exact orientation of carbonyl groups and to predict the nature of conformers present in solution. The proportion of each conformer can be inferred from a statistical thermodynamics approach by using the potential energy values of each conformer, computed by molecular mechanics methods with the TRIPOS force field, which allowed us to account for the solvent. The solvent contribution was processed by two different methods according to the nature of the interactions: whether through the dielectric constant introduced in the electrostatic potential, when interaction with solute molecules are weak or negligible, or through the computation of free energy of solvation using the algorithm SILVERWARE© for solvents explicitly interacting with the solute. When applied to tentoxin, this conformational analysis yielded results in very good agreement with the experimental data reported by Pinet et al. ( Biopolymers , 1995, Vol. 36, pp. 135–152), on both the nature of existing conformers and their relative proportions, whatever the nature of the considered solvent. © 2003 Wiley Periodicals, Inc. Biopolymers 69:363–385, 2003