Procedure for selecting starting conformations for energy minimization of nucleosides and nucleotides

The purpose of this study was to carry out a thorough search of the conformational space of various adenine-containing nucleotides, applying a previously published searching procedure, known as the representative method. This method, which reduces the number of starting conformations required to explore all the important regions of conformational space, appears to be successful in finding all (or nearly all) the putative low-energy conformations of each molecule.     

Molecular Mechanics Evaluation of the Proposed Mechanisms for the Degradation of Urea by Urease

Abstract

A thorough conformational search of all the conformations available to oxygen-bound urea within wild-type urease was carried out. Identical low energy urea conformations were obtained by a Ramachandran type plot for the N_His272-Ni1-O-C_urea and Ni1-O-Curea-N_urea dihedral angles. Ramachandran plots, with active sites and protonation states modified to model the different urease mechanisms, were used to evaluate the different mechanisms. Based upon the low energy conformations available to urea in the active site of wild-type urease one can conclude that the traditional “His320 acts as a base” mechanism is unlikely, while the N,O urea bridged and the reverse protonation mechanisms cannot be ruled out. A consensus hydrogen-bonding network that does not favor any of the mechanisms has been reconfirmed by the extensive conformational search.     

On the Bioactive Conformation of a Small Peptide and its Set of Thermodynamically Accessible Conformations

Abstract

In order to investigate the relationship between the bioactive conformation of a peptide and its set of thermodynamically accessible structures in solution, the conformational profile of the tetrapeptide Ac-Pro-Ala-Pro-Tyr-OH was characterized by computational methods. Search of the conformational space was performed within the molecular mechanics framework using the AMBER4.0 force field with an effective dielectric constant of 80. Unique structures of the peptide were compared with its bioactive conformation for the protein Streptomyces griseus Protease A, as taken from the crystal structure of the enzyme-peptide complex. The results show that the bound conformation is close to one of the unique conformations characterized in the conformational search of the isolated peptide. Moreover, the lowest energy minimum characterized in the conformational search exhibits large deviations when compared to the bound conformation of the crystal structure.     

Conformational ensemble of an intrinsically flexible loop in mitochondrial import protein Tim21 studied by modeling and molecular dynamics simulations

BackgroundTim21, a subunit of a highly dynamic translocase of the inner mitochondrial membrane (TIM23) complex, translocates proteins by interacting with subunits in the translocase of the outer membrane (TOM) complex and Tim23 channel in the TIM23 complex. A loop segment in Tim21, which is in close proximity of the binding site of Tim23, has different conformations in X-ray, NMR and new crystal contact-free space (CCFS) structures. MD simulations can provide information on the structure and dynamics of the loop in solution.MethodsThe conformational ensemble of the loop was characterized using loop modeling and molecular dynamics (MD) simulations.ResultsMD simulations confirmed mobility of the loop. Multidimensional scaling and clustering were used to characterize the dynamic conformational ensemble of the loop. Free energy landscape showed that the CCFS crystal structure occupied a low energy region as compared to the conventional X-ray crystal structure. Analysis of crystal packing indicates that the CCFS provides larger conformational space for the motions of the loop.ConclusionsOur work reported the conformational ensemble of the loop in solution, which is in agreement with the structure obtained from CCFS approach. The combination of the experimental techniques and computational methods is beneficial for studying highly flexible regions of proteins.General significanceComputational methods, such as loop modeling and MD simulations, have proved to be useful for studying conformational flexibility of proteins. These methods in integration with experimental techniques such as CCFS has the potential to transform the studies on flexible regions of proteins.     

Conformational Space Comparison of GnRH and lGnRH-III using Molecular Dynamics,Cluster Analysis and Monte Carlo Thermodynamic Integration

Abstract

The conformational space available to GnRH and 1GnRH-III was compared using 5.2 ns constant temperature and pressure molecular dynamics simulations with explicit TIP3P solvation and the AMBER v. 5.0 force field. Cluster analysis of both trajectories resulted in two groups of conformations. Results of free energy calculations, in agreement with previous experimental data, indicate that a conformation with a turn from residues 5 through 8 is preferred for GnRH in an aqueous environment. By contrast, a conformation with a helix from residues 2 through 7 with a bend from residues 6 through 10 is preferred for 1GnRH-III in an aqueous environment. The side chains of His² and Trp³ in 1GnRH-III occupy different regions of phase space and participate in weakly polar interactions different from those in GnRH. The unique conformational properties of 1GnRH-III may account for its specific anti cancer activity.     

Differences in Conformational Behavior of ATA and TAT Sequences in Single Strand DNA Trimer

Abstract

The conformational behavior of single strand (ss) TAT and ATA trimers of DNA have been studied by computational chemistry tools including CICADA software interfaced with AMBER molecular mechanics and dynamics. The Single-Coordinate-Driving (SCD) method has been used in conjunction with molecular dynamics simulated annealing. It has been revealed that the conformational flexibility of each sequence differs substantially from the other one. Four common conformational families have been found for both trimers. These are: helical, reverse-stacked (base 3), half-stacked (base 3), reverse-stacked (base 1). However, the energies of conformers representing the families are different for both the studied systems. An additional conformational family, bulged, has been found for ss(ATA), while ss(TAT) has been found also in half-stacked (base 1) conformation. In general, ss(TAT) exhibits a higher number of low energy conformations while ss(ATA) shows one interesting low energy conformational interconversion between reverse-stacked (A3) family and half-stacked (A3) family. The high conformational variability of the trimers has been confirmed by flexibility analysis and by molecular dynamics simulations, which have also shown the conformational stability of single conformational families. It has been concluded that the methodology used is able to provide a very detailed picture of the conformational space of these molecules.     

Prediction of the native conformation of a polypeptide by a statistical-mechanical procedure. III. Probable and average conformations of enkephalin

The program SMAPPS (Statistical-Mechanical Algorithm for Predicting Protein Structure) was originally designed to determine the probable and average backbone (?, ψ) conformations of a polypeptide by the application of equilibrium statistical mechanics in conjunction with an adaptive importance sampling Monte Carlo procedure. In the present paper, the algorithm has been extended to include the variation of all side-chain (χ) and peptide-bond (ω) dihedral angles of a polypeptide during the Monte Carlo search of the conformational space. To test the effectiveness of the generalized algorithm, SMAPPS was used to calculate the probable and average conformations of Met-enkephalin for which all dihedral angles of the pentapeptide were allowed to vary. The total conformational energy for each randomly generated structure of Met-enkephalin was obtained by summing over the interaction energies of all pairs of nonbonded atoms of the whole molecule. The interaction energies were computed by the program ECEPP /2 (Empirical Conformational Energy Program for Peptides). Solvent effects were not included in the computation. The results of the Monte Carlo calculation of the structure of Met-enkephalin indicate that the thermodynamically preferred conformation of the pentapeptide contains a γ-turn involving the three residues Gly²-Gly³-Phe⁴. The γ-turn conformation, however, does not correspond to the structure of lowest conformational energy. Rather, the global minimum-energy conformation, recently determined by a new optimization technique developed in this laboratory, contains a type II′ β-bend that is formed by the interaction of the four residues Gly²-Gly³-Phe⁴-Met⁵. A similar minimum-energy conformation is found by the SMAPPS procedure. The thermodynamically preferred γ-turn structure has a conformational energy of 4.93 kcal/mole higher than the β-bend structure of lowest energy but, because of the inclusion of entropy in the SMAPPS procedure, it is estimated to be ～ 9 kcal/mole lower in free energy. The calculation of the average conformation of Met-enkephalin was repeated until a total of ten independent average conformations were established. As far as the phenylalanine residue of the pentapeptide is concerned, the results of the ten independent average conformations were all found to lie in the region of conformational space corresponding to the γ-turn. These results further support the conclusion that the γturn conformation is thermodynamically favored.     

Conformational Studies on Nucleosides with Furanose Ring Modifications. 1.

Abstract

Conformational energy calculations have been presented on adenine and thymine nucleosides in which the furanose ring is replaced by 2′,3′-dideoxy-2′,3′-didehydrofuran using molecular mechanics and conformational analysis. Conformational energies have been evaluated using the MM2 and AMBER94 force field parameters at two different dielectric constants. The results are presented in terms of isoenergy contours in the conformational space of the glycosidic (χ) and C4′-C5′ (γ) bonds torsions. In general, the χ-γ interrelationships exhibit similarities with the corresponding plots for unmodified nucleosides and nucleotides, reported previously. Consistency of the calculated preferred conformations with the X-ray data is sensitive to the force field employed.     

Structure and Dynamics of Double Helical DNA in Torsion Angle Hyperspace: A Molecular Mechanics Approach

Abstract

Analysis of the conformational space populated by the torsion angles and the correlation between the conformational energy and the sequence of DNA are important for fully understanding DNA structure and function. Presence of seven variable torsion angles about single covalent bonds in DNA main chain puts a big challenge for such analysis. We have carried out restrained energy minimization studies for four representative dinucleosides, namely d(ApA):d(TpT), d(CpG):d(CpG), d(GpC):d(GpC) and d(CpA):d(TpG) to determine the energy hyperspace of DNA in context to the values of the torsion angles and the structural properties of the DNA conformations populating the favorable regions of this energy hyperspace. The torsion angles were manipulated by constraining their values at the reference points and then performing energy minimization. The energy minima obtained on the potential energy contour plots mostly correspond to the conformations populated in crystal structures of DNA. Some novel favorable conformations that are not present in crystal structure data are also found. The plots also suggest few low energy routes for conformational transitions or the associated energy barrier heights. Analyses of base pairing and stacking possibility reveal structural changes accompanying these transitions as well as the flexibility of different base steps towards variations in different torsion angles.     

New developments of the electrostatically driven monte carlo method: Test on the membrane-bound portion of melittin

The electrostatically driven Monte Carlo (EDMC) method has been greatly improved by adding a series of new features, including a procedure for cluster analysis of the accepted conformations. This information is used to guide the search for the global energy minimum. Alternative procedures for generating perturbed conformations to sample the conformational space were also included. These procedures enhance the efficiency of the method by generating a larger number of low-energy conformations. The improved EDMC method has been used to explore the conformational space of a 20-residue polypeptide chain whose sequence corresponds to the membrane-bound portion of melittin. The ECEPP/3 (Empirical Conformational Energy Program for Peptides) algorithm was used to describe the conformational energy of the chain. After an exhaustive search involving 14 independent runs, the lowest energy conformation (LEC) (−91.0 kcal/mol) of the entire study was encountered in four of the runs, while conformations higher in energy by no more than 1.8 kcal/mol were found in the remaining runs with the exception of one of them (run 8). The LEC is identical to the conformation found recently by J. Lee, H.A. Scheraga, and S. Rackovsky [(1998) “Conformational Analysis of the 20-Residue Membrane-Bound Portion of Melittin by Conformational Space Annealing,” Biopolymers, Vol. 46, pp. 103–115] as the lowest energy conformation obtained in their study using the conformational space annealing method. These results suggest that this conformation corresponds to the global energy minimum of the ECEPP/3 potential function for this specific sequence; it also appears to be the conformation of lowest free energy. © 1998 John Wiley & Sons, Inc. Biopoly 46: 117–126, 1998     

Statistical Mechanical Approach for Predicting the Transition to Non-B DNA Structures in Supercoiled DNA

Abstract

Supercoiling causes global twist of DNA structure and the supercoiled state has wide influence on conformational transition. A statistical mechanical approach was made for prediction of the transition probability to non-B DNA structures under torsional stress. A conditional partition function was defined as the sum over all possible states of the DNA sequence with basepair 1 and basepair n being in B-form helix and a recurrence formula was developed which expressed the partition function for basepair n with those for less number of pairs. This new definition permits a quick enumeration of every configuration of secondary structures. Energetic parameters of all conformations concerned, involving B-form, interior loop, cruciform and Z-form, were included in the equation. The probability of transition to each non-B conformation could be derived from these conditional partition functions. For treatment of effects of superhelicity, supercoiling energy was considered, and a twist of each conformation was determined to minimize the supercoiling energy. As the twist itself affects the transition probability, the whole scheme of equations was solved by renormalization technique. The present method permits a simultaneous treatment of serveral types of conformations under a common torsional stress.

A set of energetic parameters of DNA secondary structures has been chosen for calculation. Some DNA sequences were submitted to the calculation, and all the sequences that we submitted gave stable convergence. Some of them have been investigated the critical supercoil density for the transition to non-B DNA structures. Even though the reliability of the set of parameters was not enough, the prediction of secondary structure transition showed good agreement with reported observation. Hence, the present algorithm can estimate the probability of local conformational change of DNA under a given supercoil density, and also be employed to predict some specific sequences in which conformational change is sensitive to superhelicity.     

Beyond rotamers: a generative,probabilistic model of side chains in proteins

Background  

Accurately covering the conformational space of amino acid side chains is essential for important applications such as protein design, docking and high resolution structure prediction. Today, the most common way to capture this conformational space is through rotamer libraries - discrete collections of side chain conformations derived from experimentally determined protein structures. The discretization can be exploited to efficiently search the conformational space. However, discretizing this naturally continuous space comes at the cost of losing detailed information that is crucial for certain applications. For example, rigorously combining rotamers with physical force fields is associated with numerous problems.     

Data bank of three-dimensional structures of disaccharides,a tool to build 3-D structures of oligosaccharides

This work presents the first part of a database of conformations for all the disaccharide fragments that are found inN-glycans. The conformational study of the five disaccharides found in the oligo-mannose type are presented here. For each disaccharide, several possible conformations are described. A method is presented to obtain realistic models of oligosaccharides using molecular mechanic methods. Analysis of some possible conformations of the oligo-mannose type glycan Man₆-GlcNAc₂ is given as an illustration of the possibilities.On post-doctoral leave from CERMAV-CNRS, Grenoble, France     

Molecular Simulation of Solution Conformations of the Amyloid β-Peptide Aβ(1–42) by a Backpropagation Neural Network Model

Abstract

The predictive power of solution-dependent conformational states of the Aβ(1–42) peptide of Alzheimer's disease by an optimized backpropagation neural network was tested. It was found that the neural network simulates well the solution-dependent conformations. The model was also examined by using geometry-optimized conformations (hybrid approach of Gasteiger charges plus MM+ molecular-mechanics) where the initial coordinates were obtained by NMR solution spectroscopy.     

Complete basis set,hybrid-DFT study and NBO interpretation of conformational analysis of 2-methoxytetrahydropyran and its thiopyran and selenopyran analogues in relation to the anomeric effect

2-Methoxytetrahydropyran (1), -thiopyran (2) and -selenopyran (3) have been chosen as model compounds to investigate the origin of the anomeric effect (AE). The impacts of the hyperconjugation, electrostatic and steric interactions on the conformational preferences of compounds 1–3 have been analysed by means of complete basis set-4, hybrid-density functional theory (B3LYP/6-311+G**) based methods and natural bond orbital (NBO) interpretation. Both levels of theory showed that the axial conformations of compounds 1–3 are more stable than their equatorial conformations. The Gibbs free energy difference (G _eq–G _ax) values (i.e. ΔG _eq–ax) between the axial and equatorial conformations increase from compound 1 to compound 2 but decrease from compound 2 to compound 3. Based on the NBO results obtained, the AE associated with the electron delocalisation [i.e. Σ(endo-AE_eq + exo-AE_eq) ? Σ(endo-AE_ax + exo-AE_ax)] increase slightly from compound 1 to compound 2 but decrease from compound 2 to compound 3. Similar trend is also observed for the differences between the calculated total steric exchange energy values [i.e. Δ(TSEE)_eq–ax]. On the other hand, the calculated differences between the dipole moment values of the axial and equatorial conformations [i.e. Δ(μ_eq–μ_ax)] decrease from compound 1 to compound 3. These findings led to the proposal that the AE associated with the electron delocalisation (the hyperconjugation effect) is more significant for the explanation of the conformational preferences of compounds 1–3 than the electrostatic model. The correlations between the AE associated with the electron delocalisation, bond orders, TSEE, ΔG _eq–ax, dipole–dipole interactions, structural parameters and conformational behaviours of compounds 1–3 have been investigated.     

CRANKITE: A fast polypeptide backbone conformation sampler

Background

CRANKITE is a suite of programs for simulating backbone conformations of polypeptides and proteins. The core of the suite is an efficient Metropolis Monte Carlo sampler of backbone conformations in continuous three-dimensional space in atomic details.

Methods

In contrast to other programs relying on local Metropolis moves in the space of dihedral angles, our sampler utilizes local crankshaft rotations of rigid peptide bonds in Cartesian space.

Results

The sampler allows fast simulation and analysis of secondary structure formation and conformational changes for proteins of average length.     

Conformational Analysis of Receptor Selective Tachykinin Analogs: Senktide and Septide

Abstract

The conformational behavior in solution of two receptor selective tachykinin agonists, senktide (succiny1-D-F-MeF-G-L-M-NH₂) and septide (pQ-F-F-P-L-M-NH₂) is described. Two dimensional cross relaxation NMR spectroscopy is used together with coupling constant data to obtain interproton distance constraints. These results are used in conjunction with semi-empirical energy computations to indicate favorable conformations. Senktide is found to have a high degree of conformational order which is attributed to rotational restriction associated with the N-methylation of phenylalanine. The lowest energy conformation in accord with the experimental interproton distances contains a β-turn. Interproton distances indicate that septide exists as a random coil or in an extended chain conformation. Energy computations suggest that septide is primarily an extended chain with internal reorientation restricted by the proline residue. These results may be related to the selectivity of these peptides for different receptors, in that the analogs, with conformations more stable than tachykinins, are more receptor selective.     

Global and Local Structural Properties of the Principal Neutralizing Determinant of the HIV-1 Envelope Protein gp120

Abstract

The model of spatial structure for the principal neutralizing determinant (PND) of the HIV-1 envelope protein gpl20 is proposed in terms of two-dimensional nuclear Overhauser effect (NOE) spectroscopy data. To build the model, the NMR-based theoretical conformational analysis of synthetic PND peptides of length 40, 24, and 12 residues is carried out. The modeling of the molecular spatial structures is performed by a new approach to research of conformationally mobile peptides using the algorithms of the restrained molecular mechanics method developed earlier. The following major conclusions are made based on the analysis of the simulated peptide conformations: i) there is not unique PND structure in solution, ii) there are seven different PND structures each of which agrees with the experimental data and stereochemical criteria used in computing its spatial model, iii) the PND is characterized by irregular conformation containing a number of reverse turns, iv) all of the selected conformations are conserved in the Gly-Pro-Gly-Arg-Ala-Phe stretch, the most provable viral immunodominant epitope. These data allow to suppose that binding properties of this site are determined by the structural motif which forms the conformation of a double β-turn and appears common for all hexapeptide structures.     

CDR-H3 loop ensemble in solution – conformational selection upon antibody binding

ABSTRACT

We analyzed pairs of protein-binding, peptide-binding and hapten-binding antibodies crystallized as complex and in the absence of the antigen with and without conformational differences upon binding in the complementarity-determining region (CDR)-H3 loop. Here, we introduce a molecular dynamics-based approach to capture a diverse conformational ensemble of the CDR-H3 loop in solution. The results clearly indicate that the inherently flexible CDR-H3 loop indeed needs to be characterized as a conformational ensemble. The conformational changes of the CDR-H3 loop in all antibodies investigated follow the paradigm of conformation selection, because we observe the experimentally determined binding competent conformation without the presence of the antigen within the ensemble of pre-existing conformational states in solution before binding. We also demonstrate for several examples that the conformation observed in the antibody crystal structure without antigen present is actually selected to bind the carboxyterminal tail region of the antigen-binding fragment (Fab). Thus, special care must be taken when characterizing antibody CDR-H3 loops by Fab X-ray structures, and the possibility that pre-existing conformations are present should always be considered.     

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

The conformational space of the membrane-bound portion of melittin has been searched using the electrostatically driven Monte Carlo (EDMC) method with the ECEPP/2 (empirical conformational energy program for peptides) algorithm. The former methodology assumes that a polypeptide or protein molecule is driven toward the native structure by the combined action of electrostatic interactions and stochastic conformational changes associated with thermal movements. The algorithm produces a Monte Carlo search in the conformational hyperspace of the polypeptide using electrostatic predictions and a random sampling technique, combined with local minimization of the energy function, to locate low-energy conformations. As a result of 8 test calculations on the 20-residue membrane-bound portion of melittin, starting from six arbitrary and two completely random conformations, the method was able to locate a very low-energy region of the potential with a well-defined structure for the backbone. In all of the cases under study, the method found a cluster of similar low-energy conformations that agree well with the structure deduced from x-ray diffraction experiments and with one computed earlier by the build-up procedure.