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 solution studies of neuropeptide gamma using CD and NMR spectroscopy.

Neuropeptide gamma is one of the largest members of the tachykinin family of peptides, exhibiting strong agonistic activity towards the NK-2 tachykinin receptor. This peptide was synthesized by the solid-phase method using the Fmoc chemistry. Circular-dichroism spectroscopy (CD) investigations of this peptide were performed in phosphate buffer, in the presence of sodium dodecylsulphate (SDS) micelles and trifluoroethanol (TFE) solutions and in DMSO-d6 using the 2D NMR technique in conjunction with two different theoretical approaches. The first assumes multiconformational equilibrium of the peptide studied characterized by the values of statistical weights of low-energy conformations. These calculations were performed using three different force fields ECEPP/3, AMBER4.1 and CHARMM (implemented in the X-PLOR program). The second method incorporates interproton distance and dihedral angle constraints into the starting conformation using the Simulated Annealing algorithm (X-PLOR program). The CD experiments revealed that although the peptide studied is flexible in polar solvents, a tendency to adopt a helical structure was observed in the hydrophobic environment. The NMR data (NOE effects) indicate a helical or reverse structure in the Ile7-His12 fragment of the peptide studied in DMSO-d6 solution. The results obtained cannot be interpreted in terms of a single conformation. Most of the conformations obtained with the ECEPP/3 force field possess a high content of a helical structure. None of the conformers, obtained with the AMBER4.1 and CHARMM force fields, can be considered as the dominant one. In all conformations several beta-turns were detected and in some cases gamma-turns were also found. But in fact, it is rather difficult to select the position of the secondary element(s) present in the structure of NPgamma in solution. All conformers calculated with the X-PLOR program (with using NMR derived distance and torsion angle constraints) are stabilized by several beta-turns. Common structural motives are a type IV beta-turn in the Gln6-His12 fragment. All conformations obtained using two approaches adopt very similar turn shapes in the middle region of molecule and a random structure on the N- and C-terminal fragments.     

The Ac-RGD-NH2 peptide as a probe of slow conformational exchange of short linear peptides in DMSO

According to general belief, the conformational information on short linear peptides in solution derived at ambient temperature from NMR spectrometry represents a population-weighted average over all members of an ensemble of rapidly interconverting conformations. Usually the search for discrete conformations is concentrated at low temperatures especially when sharp NMR resonances are detected at room temperature. Using the peptide Ac-RGD-NH(2) (Ac-Arg-Gly-Asp-NH(2), Ac: acetyl) as a model system and following a new approach, we have been able to demonstrate that short linear peptides can adopt discrete conformational states in DMSO-d(6) (DMSO: dimethylsulfoxide) which vary in a way critically dependent on the reconstitution conditions used before their dissolution in DMSO-d(6). The conformers are stabilized by intramolecular hydrogen bonds, which persist at high temperatures and undergo a very slow exchange with their extended structures in the NMR chemical shift time scale. The reported findings provide clear evidence for the occurrence of solvent-induced conformational exchange and point to DMSO as a valuable medium for folding studies of short linear peptides.     

Structural studies by 1H NMR and molecular modeling of peptide substrates of thermolysin in relation with its proteasic activity in water and glycerol.

In an attempt to explain the relationship between conformations of peptide substrates of thermolysin in natural form and the experimental enzymatic cleavages, five peptides of various length were studied in two solvents H2O and glycerol, which may mimic the catalytic environmental conditions. As NMR failed to define sufficiently rough constraints to ensure a convergence of a refinement process for such short and flexible peptides, the conformational space was first searched using the MCMM method. The generated structures were then clustered in families using a 0.3A rmsd criterion and the derived structural characteristics were compared to the experimental NMR parameters. In a first approach, the NMR consistent conformations were compared with the structure of a thermolysin bound peptidic inhibitor ZG(P)LL to characterize the free-ligand predisposition to be cleaved. Further molecular dynamic calculations were performed at 300 K on the conformations corresponding to families in agreement with the ZG(P)LL structure in order to obtain information on their stability and on the trajectories of the torsion angles involved in the active site recognition. In conclusion, for four studied peptides, some conformations were found to be in agreement with 5 of the 8 cleavages experimentally observed.     

Conformational transition states of a beta-hairpin peptide between the ordered and disordered conformations in explicit water

The conformational transition states of a beta-hairpin peptide in explicit water were identified from the free energy landscapes obtained from the multicanonical ensemble, using an enhanced conformational sampling calculation. The beta-hairpin conformations were significant at 300 K in the landscape, and the typical nuclear Overhauser effect signals were reproduced, consistent with the previously reported experiment. In contrast, the disordered conformations were predominant at higher temperatures. Among the stable conformations at 300 K, there were several free energy barriers, which were not visible in the landscapes formed with the conventional parameters. We identified the transition states around the saddle points along the putative folding and unfolding paths between the beta-hairpin and the disordered conformations in the landscape. The characteristic features of these transition states are the predominant hydrophobic contacts and the several hydrogen bonds among the side-chains, as well as some of the backbone hydrogen bonds. The unfolding simulations at high temperatures, 400 K and 500 K, and their principal component analyses also provided estimates for the transition state conformations, which agreed well with those at 400 K and 500 K deduced from the current free energy landscapes at 400 K and 500 K, respectively. However, the transition states at high temperatures were much more widely distributed on the landscape than those at 300 K, and their conformations were different.     

Conformational analysis of the 20-residue membrane-bound portion of melittin by conformational space annealing

The conformational space of the 20-residue membrane-bound portion of melittin has been investigated extensively with the conformational space annealing (CSA) method and the ECEPP/3 (Empirical Conformational Energy Program for Peptides) algorithm. Starting from random conformations, the CSA method finds that there are at least five different classes of conformations, within 4 kcal/mol, which have distinct backbone structures. We find that the lowest energy conformation of this peptide from previous investigations is not the global minimum-energy conformation (GMEC); but it belongs to the second lowest energy class of the five classes found here. In four independent runs, one conformation is found repeatedly as the lowest energy conformation of the peptide (two of the four lowest energy conformations are identical; the other two have essentially identical backbone conformations but slightly different side-chain conformations). We propose this conformation, whose energy is lower than that found previously by 1.9 kcal/mol, as the GMEC of the ECEPP/3 force field. The structure of the proposed GMEC is less helical and more compact than the previous one. It appears that the CSA method can find several classes of conformations of a 20-residue peptide starting from random conformations utilizing only its amino acid sequence information. The proposed GMEC has also been found with a modified electrostatically driven Monte Carlo method [D. R. Ripoll, A. Liwo, and H.A. Scheraga (1998) “New Developments of the Electrostatically Driven Monte Carlo Method: Test on the Membrane-Bound Portion of Melittin,” Biopolymers, Vol. 46, pp. 117–126]. © 1998 John Wiley & Sons, Inc. Biopoly 46: 103–115, 1998     

Constructing ensembles of flexible fragments in native proteins by iterative stochastic elimination is relevant to protein-protein interfaces

We investigate the extent to which ensembles of flexible fragments (FF), generated by our loop conformational search method, include conformations that are near experimental and reflect conformational changes that these FFs undergo when binary protein-protein complexes are formed. Twenty-eight FFs, which are located in protein-protein interfaces and have different conformations in the bound structure (BS) and unbound structure (UbS) were extracted. The conformational space of these fragments in the BS and UbS was explored with our method which is based on the iterative stochastic elimination (ISE) algorithm. Conformational search of BSs generated bound ensembles and conformational search of UbSs produced unbound ensembles. ISE samples conformations near experimental (less than 1.05 A root mean square deviation, RMSD) for 51 out of the 56 examined fragments in the bound and unbound ensembles. In 14 out of the 28 unbound fragments, it also samples conformations within 1.05 A from the BS in the unbound ensemble. Sampling the bound conformation in the unbound ensemble demonstrates the potential biological relevance of the predicted ensemble. The 10 lowest energy conformations are the best choice for docking experiments, compared with any other 10 conformations of the ensembles. We conclude that generating conformational ensembles for FFs with ISE is relevant to FF conformations in the UbS and BS. Forming ensembles of the isolated proteins with our method prior to docking represents more comprehensively their inherent flexibility and is expected to improve docking experiments compared with results obtained by docking only UbSs.     

三肽和四肽构象空间的可视化方法

研究蛋白质寡肽构象在构象空间中的分布情况,对提取寡肽模式并构建短肽库具有重要意义。通过构建一个保距映射,将以主链原子均方根距离(root mean square distance,RMSD)为距离测度的三肽构象空间变换为一维直线上的欧氏距离空间,从而直观地展现三肽构象的聚集情况,表明三肽主链构象可以用单一变量编码。应用该特性对四肽的构象空间加以分析,将四肽构象映射到三维空间中,从而以可视的方式描述四肽构象空间的聚集情况。对短肽构象空间的初步分析表明,短肽的聚集性和二级结构有着密切的联系。在四肽构象空间中存在有自然边界的离散区域(与螺旋等结构相关),也有一些区域(与折叠等结构有关)难以进一步划分。这种方法也为以可视方式分析高维空间中肽段的聚集性给出了一种可能的方案。     

Modeling the anti-CEA antibody combining site by homology and conformational search.

A model for an antibody specific for the carcinoembryonic antigen (CEA) has been constructed using a method which combines the concept of canonical structures with conformational search. A conformational search technique is introduced which couples random generation of backbone loop conformations to a simulated annealing method for assigning side chain conformations. This technique was used both to verify conformations selected from the set of known canonical structures and to explore conformations available to the H3 loop in CEA ab initio. Canonical structures are not available for H3 due to its variability in length, sequence, and observed conformation in known antibody structures. Analysis of the results of conformational search resulted in three equally probable conformations for H3 loop in CEA. Force field energies, solvation free energies, exposure of charged residues and burial of hydrophobic residues, and packing of hydrophobic residues at the base of the loop were used as selection criteria. The existence of three equally plausible structures may reflect the high degree of flexibility expected for an exposed loop of this length. The nature of the combining site and features which could be important to interaction with antigen are discussed.     

Structure and energy landscape of a photoswitchable peptide: a replica exchange molecular dynamics study

A replica exchange molecular dynamics (REMD) simulation of a bicyclic azobenzene peptide in explicit dimethyl sulfoxide solution is presented in order to characterize the conformational structures and energy landscape of a photoswitchable peptide. It is shown that an enhanced-sampling technique such as the REMD method is essential to obtain a converged conformational sampling of the peptide at room temperature. This is because conventional MD simulations of less than approximately 100-ns length are either trapped in local minima (at 295 K) or-if run at high temperature-do not resemble the room-temperature REMD results. Calculating various nuclear Overhauser effects (NOEs) and (3)J-couplings, a good overall agreement between the REMD simulations and the NMR experiments of Renner et al. (Biopolymers 2000;54:501-514) is found. In particular, the REMD study confirms the general picture drawn by Renner et al. that the trans-isomer of the azobenzene peptide exhibits a well-defined structure, while the cis-isomer is a conformational heterogeneous system; that is, the trans-isomer occurs in 2 well-defined conformers, while the cis-isomer represents an energetically frustrated system that leads to an ensemble of conformational structures. Employing a principal component analysis of the REMD data, the free energy landscape of the systems is studied at various temperatures. The implications for the folding and unfolding pathways of the system are discussed.     

Energy landscape of a peptide consisting of alpha-helix, 3(10)-helix, beta-turn, beta-hairpin, and other disordered conformations

The energy landscape of a peptide [Ace-Lys-Gln-Cys-Arg-Glu-Arg-Ala-Nme] in explicit water was studied with a multicanonical molecular dynamics simulation, and the AMBER parm96 force field was used for the energy calculation. The peptide was taken from the recognition helix of the DNA-binding protein, c-MYB: A rugged energy landscape was obtained, in which the random-coil conformations were dominant at room temperature. The CD spectra of the synthesized peptide revealed that it is in the random state at room temperature. However, the 300 K canonical ensemble, Q(300K), contained alpha-helix, 3(10)-helix, beta-turn, and beta-hairpin structures with small but notable probabilities of existence. The complete alpha-helix, imperfect alpha-helix, and random-coil conformations were separated from one another in the conformational space. This means that the peptide must overcome energy barriers to form the alpha-helix. The overcoming process may correspond to the hydrogen-bond rearrangements from peptide-water to peptide-peptide interactions. The beta-turn, imperfect 3(10)-helix, and beta-hairpin structures, among which there are no energy barriers at 300 K, were embedded in the ensemble of the random-coil conformations. Two types of beta-hairpin with different beta-turn regions were observed in Q(300K). The two beta-hairpin structures may have different mechanisms for the beta-hairpin formation. The current study proposes a scheme that the random state of this peptide consists of both ordered and disordered conformations. In contrast, the energy landscape obtained from the parm94 force field was funnel like, in which the peptide formed the helical conformation at room temperature and random coil at high temperature.     

Structural model of the anti-snake-toxin antibody,Mα2,3

We present results of structural modeling of the variable fragment of Mα2,3, an antibody capable of neutralizing all short snake toxins. Three different methods were used to model the hypervariable loops: the conformational search algorithm CONGEN (Bruccoleri and Karplus, Biopolymers 26:137–168, 1987), high-temperature molecular dynamics (Bruccoleri and Karplus, Biopolymers 29:1847–1862, 1990), and a combined knowledge-based and energy-based algorithm (Martin et al., Proc. Natl. Acad. Sci. USA 86:9268–9272, 1989). Ninety plausible conformations were generated and were clustered into 13 classes. The clustering results indicate that there was little overlap of the conformational space explored by the different methods. Canonical loop structures were found by all methods for two of the loops, in agreement with previously established empirical modeling criteria. Nine of the 13 classes of structure were rejected on the ground of their lacking common features of antibody combining-site structure. The remaining four models were refined using restrained molecular dynamics. It was found that interconversion between the four resulting structures is possible with no significant energy barriers, suggesting that they are in thermodynamic equilibrium at 300 K. Features of the combining-site structure likely to be particularly important for antigen binding are discussed. © 1996 Wiley-Liss, Inc.     

Structure of the 21-30 fragment of amyloid beta-protein

Folding and self-assembly of the 42-residue amyloid beta-protein (Abeta) are linked to Alzheimer's disease (AD). The 21-30 region of Abeta, Abeta(21-30), is resistant to proteolysis and is believed to nucleate the folding of full-length Abeta. The conformational space accessible to the Abeta(21-30) peptide is investigated by using replica exchange molecular dynamics simulations in explicit solvent. Conformations belonging to the global free energy minimum (the "native" state) from simulation are in good agreement with reported NMR structures. These conformations possess a bend motif spanning the central residues V24-K28. This bend is stabilized by a network of hydrogen bonds involving the side chain of residue D23 and the amide hydrogens of adjacent residues G25, S26, N27, and K28, as well as by a salt bridge formed between side chains of K28 and E22. The non-native states of this peptide are compact and retain a native-like bend topology. The persistence of structure in the denatured state may account for the resistance of this peptide to protease degradation and aggregation, even at elevated temperatures.     

Solution conformational study of nociceptin an(d its 1-13 and 1-11 fragments using circular dichroism and two-dimensional NMR in conjunction with theoretical conformational analysis.

Conformational studies of nociceptin (NC-NH2), its fully active fragment, NC(1-13)-NH2, and two significantly less potent fragments, NC(1-13)-OH and NC(1-11)-OH, were conducted in water and TFE solutions by the employment of circular dichroism, and in DMSO-d6 by 2DNMR spectroscopy in conjunction with theoretical conformational analysis. The conformations of all thepeptides studied were calculated taking two approaches. The first assumes multiconformational equilibrium of the peptide studied, which is characterized by a set of conformations (and their statistical weight values)obtained from a global conformational analysis using three methods: the electrostatically driven Monte-Carlo (EDMC) with the ECEPP/3 force field, the simulated annealing (SA) protocols in the AMBER and CHARMM force fields. The second approach incorporates the interproton distance and dihedral angle constraints into the starting conformation. Calculations were performed using the distance geometry and SA protocol in the CHARMM force field implemented in the X-PLOR program. The CD experiments indicated that for the active peptides, hydrophobic solvents induced a significantly higher (compared with those remaining)content order, probably a helical structure. Unfortunately, as a result of the conformational flexibility of thepeptides, the analysis of conformations obtained with both approaches and different force fields did not alllow the selection of any structural elements of the NC peptides that might be connected with their bioactivity. The only common element found in most conformations of the active peptides was a helical character of fragment 8-13, which allowed the side chains of basic amino acid residues to be exposed to the outside of the molecule and probably to interact with the ORL1 receptor.     

The relationship between the flexibility of proteins and their conformational states on forming protein-protein complexes with an application to protein-protein docking

We investigate the extent to which the conformational fluctuations of proteins in solution reflect the conformational changes that they undergo when they form binary protein-protein complexes. To do this, we study a set of 41 proteins that form such complexes and whose three-dimensional structures are known, both bound in the complex and unbound. We carry out molecular dynamics simulations of each protein, starting from the unbound structure, and analyze the resulting conformational fluctuations in trajectories of 5 ns in length, comparing with the structure in the complex. It is found that fluctuations take some parts of the molecules into regions of conformational space close to the bound state (or give information about it), but at no point in the simulation does each protein as whole sample the complete bound state. Subsequent use of conformations from a clustered MD ensemble in rigid-body docking is nevertheless partially successful when compared to docking the unbound conformations, as long as the unbound conformations are themselves included with the MD conformations and the whole globally rescored. For one key example where sub-domain motion is present, a ribonuclease inhibitor, principal components analysis of the MD was applied and was also able to produce conformations for docking that gave enhanced results compared to the unbound. The most significant finding is that core interface residues show a tendency to be less mobile (by size of fluctuation or entropy) than the rest of the surface even when the other binding partner is absent, and conversely the peripheral interface residues are more mobile. This surprising result, consistent across up to 40 of the 41 proteins, suggests different roles for these regions in protein recognition and binding, and suggests ways that docking algorithms could be improved by treating these regions differently in the docking process.     

Folding landscapes of the Alzheimer amyloid-beta(12-28) peptide

The energy landscape for folding of the 12-28 fragment of the Alzheimer amyloid beta (Abeta) peptide is characterized using replica-exchange molecular dynamics simulations with an all-atom peptide model and explicit solvent. At physiological temperatures, the peptide exists mostly as a collapsed random coil, populating a small fraction (less than 10%) of hairpins with a beta-turn at position V18F19, with another 10% of hairpin-like conformations possessing a bend rather than a turn in the central VFFA positions. A small fraction of the populated states, approximately 14%, adopt polyproline II (PPII) conformations. Folding of the structured hairpin states proceeds through the assembly of two locally stable segments, VFFAE and EDVGS. The interactions stabilizing these locally folded structural motifs are in conflict with those stabilizing the global fold of A12-28, a signature of underlying residual frustration in this peptide. At increased temperature, the population of both beta-strand and PPII conformations diminishes in favor of beta-turn and random-coil states. On the basis of the conformational preferences of Abeta 12-28 monomers, two models for the molecular structure of amyloid fibrils formed by this peptide are proposed.     

Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing-potential

During replica exchange molecular dynamics (RexMD) simulations, several replicas of a system are simulated at different temperatures in parallel allowing for exchange between replicas at frequent intervals. This technique allows significantly improved sampling of conformational space and is increasingly being used for structure prediction of peptides and proteins. A drawback of the standard temperature RexMD is the rapid increase of the replica number with increasing system size to cover a desired temperature range. In an effort to limit the number of replicas, a new Hamiltonian-RexMD method has been developed that is specifically designed to enhance the sampling of peptide and protein conformations by applying various levels of a backbone biasing potential for each replica run. The biasing potential lowers the barrier for backbone dihedral transitions and promotes enhanced peptide backbone transitions along the replica coordinate. The application on several peptide cases including in all cases explicit solvent indicates significantly improved conformational sampling when compared with standard MD simulations. This was achieved with a very modest number of 5-7 replicas for each simulation system making it ideally suited for peptide and protein folding simulations as well as refinement of protein model structures in the presence of explicit solvent.     

The human prion protein alpha2 helix: a thermodynamic study of its conformational preferences

We have synthesized both free and terminally-blocked peptide corresponding to the second helical region of the globular domain of normal human prion protein, which has recently gained the attention of structural biologists because of a possible role in the nucleation process and fibrillization of prion protein. The profile of the circular dichroism spectrum of the free peptide was that typical of alpha-helix, but was converted to that of beta-structure in about 16 h. Instead, below 2.1 x 10(-5) M, the spectrum of the blocked peptide exhibited a single band centered at 200 nm, unequivocally associated to random conformations, which did not evolve even after 24 h. Conformational preferences of this last peptide have been investigated as a function of temperature, using trifluoroethanol or low-concentration sodium dodecyl sulfate as alpha- or beta-structure inducers, respectively. Extrapolation of free energy data to zero concentration of structuring agent highlighted that the peptide prefers alpha-helical to beta-type organization, in spite of results from prediction algorithms. However, the free energy difference between the two forms, as obtained by a thermodynamic cycle, is subtle (roughly 5-8 kJ mol(-1) at any temperature from 280 K to 350 K), suggesting conformational ambivalence. This result supports the view that, in the prion protein, the structural behavior of the peptide is governed by the cellular microenvironment.     

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.