Study of the conformational profile of the norbornane analogues of phenylalanine.

The conformational profile of the eight stereoisomeric 2-amino-3-phenylnorbornane-2-carboxylic acids (2-amino-3-phenylbicyclo[2.2.1]heptane-2-carboxylic acids) has been assessed by computational methods. These molecules constitute a series of four enantiomeric pairs that can be considered as rigid analogues of either L- or D-phenylalanine. The conformational space of their N-acetyl methylamide derivatives has been explored within the molecular mechanics framework, using the parm94 set of parameters of the AMBER force field. Local minimum energy conformations have been further investigated at the ab initio level by means of the Hartree-Fock and second order Moller-Plesset perturbation energy calculations using a 6-31G(d) basis set. The results of the present work suggest that the bulky norbornane structure induces two kinds of conformational constraints on the residues. On one hand, those of a steric nature directly imposed by the bicycle on the peptide backbone and, on the other hand, those that limit the orientations attainable by the phenyl ring which, in turn, reduces further the flexibility of the peptide backbone. A comparative analysis of the conformational profile of the phenylnorbornane amino acids with that of the norbornane amino acids devoid of the beta-phenyl substituent suggests that the norbornane system hampers the residue to adopt extended conformations in favour of C7-like structures. However, the bicycle itself does not impart a clear preference for any of the two possible C7 minima. It is the aromatic side chain, which is forced to adopt an almost eclipsed orientation, that breaks this symmetry introducing a marked preference for a single region of the (phi, psi) conformational space in each of the phenylalanine norbornane analogues investigated.     

Conformational features of the Dippu-AST 8 neuropeptide from the cockroach Diploptera punctata

Spatial structure and conformational properties of the Dippu-AST 8 (allatostatin III) neuropeptide have been investigated by the molecular mechanics method. The conformational energy and geometrical parameters corresponding to the low-energy states of the molecule are obtained. A single backbone conformation with a very restricted set of Ser 3, Phe 4 and Leu 9 amino acids positions is observed for Dippu-AST 8 neuropeptide.     

Computational Study of the Conformational Domains of Peptide T

The conformational preferences of peptide T (ASTTTNYT) were analysed by means of computational methods. A thorough exploration of the conformational space was carried out within the framework of the molecular mechanics approach, using simulated annealing as a searching strategy. Specifically, in order to obtain a subset of low-energy conformations with energies close to the global minimum as complete as possible, a simulated annealing protocol was repeated several times in a recursive fashion. The results of the search indicate that the peptide exhibits a α-helical character although most of the conformations characterized, including the global minimum, can be described as bent conformations. Conformations exhibiting β-turn motives previously proposed from NMR studies were also characterized, although they are not very predominant in the set of low-energy conformations. © 1997 European Peptide Society and John Wiley & Sons, Ltd.     

Conformation-activity relationships of cyclic dermorphin analogues

A theoretical conformational analysis (molecular mechanics study) of nine cyclic tetrapeptides, structurally related to the highly mu-receptor-selective dermorphin analogue H-Tyr-D-Orn-Phe-Asp-NH2, was performed. These compounds display considerable diversity in their mu-receptor affinity and selectivity. A systematic search and subsequent energy minimization in absence of the exocyclic Tyr1 residue and Phe3 side chain revealed the constrained nature of the 11-13-membered ring structures contained in these analogues. No more than four low-energy conformers (within 2 kcal/mol of the lowest energy conformation) were found in each case. After attachment of the Tyr1 moiety and Phe3 side chain to the "bare" low-energy ring structures, a systematic search and energy minimization of these exocyclic moieties resulted in a limited number of low-energy conformational families for all compounds. Five analogues with high mu-receptor affinity--H-Tyr-D-Orn-Phe-Asp-NH2, H-Tyr-D-Orn-Phe-D-Asp-NH2, H-Tyr-D-Asp-Phe-Orn-NH2, H-Tyr-D-Asp-Phe-A2bu-NH2 (A2 bu: alpha, gamma-diaminobutyric acid) and H-Tyr-D-Cys-Phe-Cys-NH2--all showed a tilted stacking interaction between the Tyr1 and Phe3 aromatic rings in the lowest or second lowest energy conformation found. The same kind of stacking was not possible in low-energy conformers of the four analogues with poor affinity for the mu-receptor [H-Tyr-L-Orn-Phe-Asp-NH2, H-Tyr-D-Orn-D-Phe-Asp-NH2, H-Tyr-D-Orn-Phe(NMe)-Asp-NH2 [Phe(NMe): N alpha-methylphenylalanine], and H-Tyr-D-Orn-Phg-Asp-NH2 (Phg: phenylglycine)].(ABSTRACT TRUNCATED AT 250 WORDS)     

Conformational features responsible for the binding of cyclic analogues of enkephalin to opioid receptors. III. Probable binding conformations of mu-agonists with phenylalanine in position 3.

Theoretical conformational analysis was carried out for several tetrapeptide analogues of beta-casomorphin and dermorphin containing a Phe residue in position 3. Sets of low-energy backbone structures of the mu-selective peptides [N-Me-Phe3, D-Pro4]-morphiceptin and Tyr-D-Orn-Phe-Asp-NH2 were obtained. These sets of structures were compared for geometrical similarity between themselves and with the low-energy conformations found for the delta-selective peptide Tyr-D-Cys-Phe-D-Pen-OH and nonactive peptide Tyr-Orn-Phe-Asp-NH2. Two pairs of geometrically similar conformations of mu-selective peptides, sharing no similarity with the conformations of peptides showing low affinity to the mu-receptor, were selected as two alternative models of probable mu-receptor-bound backbone conformations. Both models share geometrical similarity with the low-energy structures of the linear mu-selective peptide Tyr-D-Ala-Phe-Phe-NH2. Putative binding conformations of Tyr1 and Phe3 side chains are also discussed.     

Characterization of the conformational domains of bradykinin by computational methods

The AMBER 4.0 force field was used to perform a characterization of the conformational profile of the nonapeptide bradykinin. A thorough conformational search was carried out using molecular dynamics as sampling technique, by computing cycles of high (900 K) and low (300 K) temperature trajectories. A total of 2400 minima were generated and subsequently clustered using the root-mean-square of the backbone dihedral angles as criterium. After the use of a tolerance value of 20deg;, the conformations were clustered in 233 unique conformations with energies up to 40 kcal/mol above the lowest minimum. The analysis of the low-energy conformations indicate that the peptide exhibits a high tendency to adopt a β-turn at the C-terminus and a propensity to adopt a bent structure at the N-terminus. These results are in agreement with the experimental evidence reported in the literature and provide detailed information necessary to understand the conformational preferences of the peptide.     

Conformational studies of irreversible HIV-1 protease inhibitors containing cis-epoxide as an amide isostere.

We have carried out NMR and molecular modeling studies of peptidomimetic HIV-1 protease inhibitors, LB71116: Qc-Asn-Phepsi[(1R,2S)-cis-epoxide]Gly-NH-CH(isopropyl)2 where Qc stands for quinaldic acid and LB71148: Qc-(SMe)Pen(O)2-Phepsi[(1R,2S)-cis-epoxide]Gly-NH-CH(isoprop yl)2 where (SMe)Pen(O)2 stands for S-methyl-S-dioxo-penicillamine. Through conformational calculations and NMR data analysis, we have obtained preferred conformations of the two inhibitors in solution. To our knowledge, this work is one of the first extensive conformational studies of peptidomimetics containing cis-epoxide amide isostere. The resulting preferred conformations contain extended structures. In these conformations, the psi of Phe(cep) is maintained about 130 degrees and the phi angle of (cep)Gly prefers +/- 150 degrees [where Phe(cep) and (cep)Gly are the residues generated by the replacement of the Phe-Gly peptide bond with cis-epoxide]. Two conformations were commonly observed in the preferred conformations of each inhibitor. Through restrained molecular dynamics simulating the hydrogen bond formation between our inhibitor and a water molecule ('flap water'), one of the conformations is assumed as the conformation which can bind to the enzyme without large conformational changes. Recently, we had the opportunity to compare the selected preferred conformation with the binding conformation of LB71116 observed from the X-ray studies of the complex between LB71116 and HIV-1 protease. These two conformations are surprisingly similar to each other. Thus, we can explain high activity and selectivity of our inhibitors to the HIV-1 protease by the similarity between the preferred conformations in solution and the binding conformation.     

Conformational-functional relationships of tetragastrin analogs]

Sets of low-energy backbone conformations of the active tetragastrin analogue Boc-Trp-Leu-Asp-Phe-NH2 and two competitive antagonists Boc-Trp-Leu psi (CH2NH)-Asp-Phe-NH2 and Boc-Trp-Leu-Asp-O-CH2-CH2-C6H5 were obtained using theoretical conformational analysis methods. Groups of the conformations were selected for the three analogues, allowing a spatial matching of Trp, Asp and Phe residues responsible for the gastrin receptor binding. Three conformations possessing the lowest energies among the geometrically similar structures of these three peptides are suggested as a model for the "receptor-bound" conformations of these analogues. Backbone spatial folding resembling an alpha-helix turn is characteristic of these conformations. The correspondence of the proposed model to the available data on structure--activity relationships for tetragastrin analogues is discussed. Orientations of the putative receptor-bound conformations in a "water--lypophylic medium" two-phase system were investigated.     

Conformation-activity relationships of cyclo-constrained micro/delta opioid agonists derived from the N-terminal tetrapeptide segment of dermorphin/deltorphin.

The N-terminal tetrapeptide segments of dermorphin (Tyr-D-Ala-Phe-Gly-Tyr-Pro-Ser-NH(2)) and deltorphin (Tyr-D-Ala-Phe-Asp/Glu-Val-Val-Gly-NH(2)) are agonists at the opioid receptors micro and delta, respectively. [D-Arg(2), Lys(4)]-dermorphin-(1-4) amide (Tyr-D-Arg-Phe-Lys-NH(2), DALDA) and [Dmt(1)]DALDA (where Dmt is 2',6'-dimethyltyrosine) are among the most potent and selective micro-agonists reported to date, both in vitro (having picomolar micro receptor affinity) and in vivo. In this communication, conformation-activity studies of the following four cyclic analogs of DALDA are presented and discussed: the lead peptide S(2),S(4)-cyclo (Tyr-D-Cys-Phe-Cys-NH(2)), constrained by means of an S(4.2)--S(4.4) disulfide between Cys(2) and Cys(4); its two cis and trans C(4.2)--C(4.4)-olefinic dicarba analogs, and the product of saturation of them both. They are potent nonselective or moderately micro-selective opioid agonists in vitro.They have been synthesized and tested earlier [Berezowska I, Chung NN, Lemieux C, Wilkes BC, and Schiller PW, Acta Biochim Polon 53, 2006, 73-76]. We have studied their conformations using NMR and molecular dynamics. With major conformational constraints imposed by the 11-membered ring spanning residues 2-4, they show well defined conformations of this ring, while the exocylic Tyr(1) and Phe(3) side chains still have significant conformational freedom. The more active and selective micro versus delta disulfide and saturated dicarba agonists seem to have in common: (i) their ring structures more flexilble than those of the other two and (ii) their ring structures similar to each other and more diverse than those in the other two. Given this and the small size of the peptides having confirmed bioactivity profiles, there is a chance that their conformations determined in solution approach receptor-bound conformations. Copyright (c) 2008 European Peptide Society and John Wiley & Sons, Ltd.     

Structural and biological effects of a beta2- or beta3-amino acid insertion in a peptide.

Molecular mechanics calculations on conformers of Ac-HGly-NHMe, Ac-beta2-HAla-NHMe and Ac-beta3-HAla-NHMe indicate that low-energy conformations of the beta-amino acids backbone, corresponding to gauche rotamers around the Calpha-Cbeta bond, may overlap canonical backbone conformers observed for alpha-amino acids. Therefore, Substance P (SP) was used as a model peptide to analyse the structural and biological consequences of the substitution of Phe7 and Phe8 by (R)-beta2-HPhe and of Gly9 by HGly (R)-beta2-HAla or (S)-beta3-HAla. [(R)-beta2-HAla9]SP has pharmacological potency similar to that of SP while [HGly9]SP and [(S)-beta3-HAla9]SP show a 30- to 50-fold decrease in biological activities. The three analogues modified at position 9 are more resistant to degradation by angiotensin converting enzyme than SP and [Ala9]SP. NMR analysis of these SP analogues suggest that a beta-amino acid insertion in position 9 does not affect the overall backbone conformation. Altogether these data suggest that [HGly9]SP, [(S)-beta3-HAla9]SP and [(R)-beta2-HAla9]SP could adopt backbone conformations similar to that of SP, [Ala9]SP and [Pro9]SP. In contrast, incorporation of beta2-HPhe in position 7 and 8 of SP led to peptides that are almost devoid of biological activity. Thus, a beta-amino acid could replace an alpha-amino acid within the sequence of a bioactive peptide provided that the additional methylene group does not cause steric hindrance and does not confine orientations of the side chain to regions of space different from those permitted in the alpha-amino acid.     

Fast and accurate computation of the 13C chemical shifts for an alanine-rich peptide

The purpose of this work is, first, to present a fast and accurate technique to compute Boltzmann-averaged values of the quantum-chemical 13C chemical shifts for each amino acid in oligopeptides, demonstrated here by an application to the peptide Ac-XXAAAAAAAOO-NH2 (where X denotes diaminobutyric acid, A is alanine, and O is ornithine) [XAO] and, second, to discuss the capability of the 13Calpha and 13Cbeta chemical shifts to distinguish the PP(II) conformation from the alpha-helix and statistical-coil conformations. Use is made of a combination of approaches, summarized as follows: (1) derivation of an ensemble of conformations by using a molecular mechanics technique; (2) use of a clustering procedure to form families and build a reduced set of conformations consisting of the lowest-energy conformations of each family, and (3) computation of the 13C chemical shifts for the lowest-energy conformations of each family, using a quantum-chemical approach that treats a selected residue, or group of residues, with a 6-311+G(2d,p) locally-dense basis set, while the remaining residues in the sequence are treated with a 3-21G basis set. The whole procedure is quite accurate and speeds up the computation of the Boltzmann-averaged values of the 13C-chemical shifts by several orders of magnitude. The present application sheds some light on the conformational preference for alanine and non-alanine residues to occupy the PP(II) helical region of the Ramachandran map.     

Refinement of the thrombin-bound structure of a hirudin peptide by a restrained electrostatically driven Monte Carlo method.

Energy refinement of the structure of a linear peptide, hirudin56-65, bound to thrombin was carried out using a conformational search method in combination with restrained minimization. Five conformations originated from nmr data and distance geometry calculations having a similar global folding pattern but quite different backbone conformations were used as the starting structures. As a result of this approach, a series of low-energy conformations compatible with a set of upper and lower bounds of interproton distances determined from transferred nuclear Overhauser effects were found. A comparison among the lowest energy conformations of each run showed that the combination of energy refinement plus distance constraints led to a very well-defined structure for both the backbone and the side chains of the last 7 residues of the polypeptide. Furthermore, the low-energy conformations generated with this technique contain a segment of 3(10)-helix involving the last 5 residues at the COOH terminal end.     

Density functional conformational study of 2-O-sulfated 3,6 anhydro-α-D-galactose and of neo-κ- and ι-carrabiose molecules in gas phase and water

We examined the conformational preferences of the 2-O-sulfated-3,6-α-D-anhydrogalactose (compound I) and two 1,3 linked disaccharides constituting-κ or ι-carrageenans using density functional and ab initio methods in gas phase and aqueous solution. Systematic modifications of two torsion angles leading to 324 and 144 starting geometries for the compound I and each disaccharide were used to generate adiabatic maps using B3LYP/6-31G(d). The lower energy conformers were then fully optimized using B3LYP, B3PW91 and MP2 with several basis sets. Overall, we discuss the impact of full relaxation on the energy and structure of the dominant conformations, present the performance comparison with previous molecular mechanics calculations if available, and determine whether our results are impacted, when polarization and diffuse functions are added to the 6-31G(d) basis set, or when the MP2 level of theory is used.     

Conformational analysis of [Cpp1, Sar7, Arg8] vasopressin by 1H-NMR spectroscopy and molecular mechanics calculations.

A combined 1H-NMR and molecular mechanics study of [Cpp1, Sar7]AVP was performed in order to select the most probable conformations in DMSO solutions. The NMR constraints obtained were employed in the selection of starting conformations of the cyclic moiety of the analog. In particular, the diminished accessibility of the Asn5 NH proton to solvent and the close contact between Cpp1 and Cys6 C alpha H protons suggests a beta-turn conformation at the Phe3-Gln4 residues. Energy minimization was carried out both in the ECEPP/2 (rigid-valence geometry) and in the AMBER (flexible-valence geometry) force fields. Comparison of the experimental and calculated values of NMR characteristics has revealed that conformations containing type I, II, and III beta-turns at the Phe3-Gln4 residues are in reasonable agreement with the experimental data, with a dynamic equilibrium between the beta I (beta III) and beta II type structures of the cyclic part being the most probable. All of these conformations prefer the negative chirality of the disulfide bridge (theta 3 approximately -90 degrees). Five representative conformations were chosen for the acyclic tail: one with a beta I, one with a beta II'-turn at the Sar7-Arg8 residues, two extended-type conformations, and a conformation with a gamma-turn at Sar7. Because only high-energy extended conformations were in agreement with NMR data, it was concluded that the acyclic tail has considerable conformational flexibility in solution. The conformations obtained are discussed in terms of the structure-function relationship of the neurohypophyseal hormone analogs.     

Influence of local interactions on protein structure. I. Conformational energy studies of N-acetyl-N′-methylamides of pro-X and X-pro dipeptides

Conformational energy calculations using an Empirical Conformational Energy Program for Peptides (ECEPP) were carried out on the N-acetyl-N′-methylamides of Pro-X, where X = Ala, Asn, Asp, Gly, Leu, Phe, Ser, and Val, and of X-Pro, where X = Ala, Asn, Gly, and Pro. The conformational energy was minimized from starting conformations which included all combinations of low-energy single-residue minima and several standard bend structures. It was found that almost all resulting minima are combinations of low-energy single-residue minima, suggesting that intra residue interactions predominate in determining conformation. The calculations also indicate, however, that inter residue interactions can be important. In addition, librational entropy was found to influence the relative stabilities of some minima. Because of the existence of 10–100 low-energy minima for each dipeptide, the normalized statistical weight of an individual minimum rarely exceeds 0.3, suggesting that these dipeptides have considerable conformational flexibility and exist as statistical ensembles of low-energy structures. The propensity of each dipeptide to form bend conformations was calculated, and the results were compared with available experimental data. It was found that bends are favored in Pro-X dipeptides because ?_Pro is fixed by the pyrrolidine ring in a conformation which is frequently found in bends, but that bends are not favored in X-Pro dipeptides because interactions between the X residue and the pyrrolidine ring restrict the X residue to conformations which are not usually found in bends.     

Analysis of the conformational profile of trishomocubane amino acid dipeptide

4-Amino-(D3)-trishomocubane-4-carboxylic acid is a constrained alpha-amino acid residue that exhibits promising conformational characteristics, i.e., helical and beta-turns. As part of the development of conformational guidelines for the design of peptides and protein surrogates, the conformational energy calculations on trishomocubane using molecular mechanics and ab initio methods are presented. The C(alpha) carbon of trishomocubane forms part of the cyclic structure, and consequently a peptidic environment was simulated with an acetyl group on its N-terminus and a methylamide group on its C-terminus. Ramachandran maps computed at the molecular mechanics level using the standard AMBER (parm94) force field libraries compared reasonably well with the corresponding maps computed at the Hartree Fock level, using the 6-31G* basis set. Trishomocubane peptide (Ac-Tris-NHMe) is characterized by four low energy conformers corresponding to the C7ax, C7eq, 3(10), and alpha(L) helical structures.     

Approach to studying the structure-activity relationship of native peptides. I. Structural organization of the opioid peptide Tyr-Glu-Gly-Phe-Met-Arg-Gly-Leu and its artificial analogs]

Theoretical conformational analysis was carried out for the octapeptide Tyr1-Gly2-Gly3-Phe4-Met5-Arg6-Gly7-Leu8. Possible structure of the opioid peptide under physiological conditions may be described by a set of low-energy conformations belonging to 14 different forms of the backbone. The solution of the "reverse conformational problem" for the opioid peptide enables one to predict the modified amino acid sequences (Ala2, D-Ala2, Ala3, D-Ala3, Ala7, D-Ala7, MeMet5, MeArg6-analogues) which may assume one of the low-energy states of the native hormone. The influence of the solute was not taken into account in our calculations.     

Conformational preferences of 1,4,7-trithiacyclononane: a molecular mechanics and density functional theory study

Conformational preferences of 1,4,7-trithiacyclononane were studied using a highly efficient sampling technique based on local nonstochastic deformations and the MM2(91) force field. The results show that conformers that the molecule adopts in the crystal state were found to be low-energy conformers (LECs) within 5 kcal mol(-1) of the global minimum. A conformation with C1 symmetry was the global minimum and the C3 and C2 conformations were calculated to be 0.03 and 1.78 kcal mol(-1) higher in energy, respectively. The structures were further minimized using Density Functional Theory (DFT) calculations with two different functionals. The C2 and the C1 conformations were found to be LECs with the C3 conformation more than 4.0 kcal mol(-1) above the global minimum. The relative energies and structural ordering obtained using the BP86 functional are in agreement with the previously reported relative energies calculated using second-order Moller-Plesset (MP2) ab initio calculations. With the energy ordering being dependent on the molecular mechanics force field used, the approach of MM-->DFT (searching exhaustively the available conformational space at the MM level followed by generating the energy ordering through DFT calculations) appears to be appropriate for thiacrown ethers.     

Block-units method for conformational calculations of large nucleic acid chains. II. The two-hierarchical approach and its application to conformational arrangement of the unusual T psi C loop of rabbit tRNA(Val)

The two-level hierarchical methodology is suggested for conformational calculations of large fragments of nucleic acids. The method of the first level is intended for performing a fast screening of the conformational phase space. The high-level method may be used to refine structurally important conformations. The method of the first level is the block-units method, which has been developed specially for these purposes (see part I). It has been shown that the block-units method allows the satisfactory calculation of the structure parameters of the optimal conformations of polynucleotides. The results of the conformational rearrangement calculations of the T psi C loop of the tRNA(Phe) after modification of its sequence are represented.     

Molecular mechanics (MM3(pi)) conformational analysis of molecules containing conjugated pi-electron fragments: Leucomycin-V

The conformations of the 16-membered macrolide antibiotic leucomycin-V (1) were studied with molecular mechanics. Leucomycin-V contains a conjugated pi-electron fragment and necessitates special treatment with the MM3(pi) modeling protocol. Comparison was made with results from the standard MM3 scheme. The CONFLEX conformational search procedure was used for finding low-energy conformations. The computed data are indicative for the existence of mainly one conformation of the macro-ring of 1 and minor participation of several others. Intramolecular hydrogen bonds play important roles for the preferred geometry of the macro-ring and the conformations of the side chains. The most probable macro-ring conformation of 1 is very similar to the preferred conformation of another 16-ring macrolide antibiotic, tylosin. The same order of conformational preference for 1 was estimated with the MM3 and the MM3(pi) methods. Surprisingly, when changing the chirality of the C(9) macro-ring atom of 1, the two methods produced different order of conformational preferences for the 9-epi form (2), as well as enhanced population of several clusters of conformations.