Prediction of a 12-residue loop in bovine pancreatic trypsin inhibitor: effects of buried water

The x-ray conformations of 5-, 7-, 9-, and 12-residue loops in bovine pancreatic trypsin inhibitor (BPTI) were predicted by the use of multiple independent Monte Carlo simulating annealing (MCSA) runs starting from random conformations. Four buried water molecules interacted with a 12-residue loop that started at residue 8 and ended at residue 19, and that included the binding region. The final conformation at the end of an MCSA run was characterized. Solvation free energy based on the solvent accessible surface area was included in the energy function at low simulated annealing temperatures. Conformational states were interactively separated by a recently developed algorithm. Computed loops were characterized in terms of total energy, and backbone and side chain root mean square deviations (RMSDs) between computed native loop conformations and the x-ray conformation. The 12-residue loop was computed with and without buried water [called WL12(8-19) and L12(8-19), respectively]. The backbone was reliably and reproducibly computed to within 1.1 A in L12(8-19) and 0.9 A in WL12(8-19). L12(8-19) required significantly more MCSA runs to achieve the same level of reproducibility as WL12(8-19). Based on the size of the cluster of low energy native loop conformations, and the computational effort, WL12(8-19) had greater entropy. In calculations of 7-, 9-, and 12-residue loops without buried water, the effects of buried water became obvious in the 12-residue loop calculation, which interacted with all four buried water molecules. Nearly all conformations of the native loop conformer had a hydrogen bond between the Lys 15 side chain and the backbone of Gly 12, Pro 13, and Cys 14, which may have implications in the rate of exchange of buried water with bulk solvent and in protein folding. The present version of MCSA program was more efficient than earlier versions.     

A general Monte Carlo/simulated annealing algorithm for resonance assignment in NMR of uniformly labeled biopolymers

We describe a general computational approach to site-specific resonance assignments in multidimensional NMR studies of uniformly ¹⁵N,¹³C-labeled biopolymers, based on a simple Monte Carlo/simulated annealing (MCSA) algorithm contained in the program MCASSIGN2. Input to MCASSIGN2 includes lists of multidimensional signals in the NMR spectra with their possible residue-type assignments (which need not be unique), the biopolymer sequence, and a table that describes the connections that relate one signal list to another. As output, MCASSIGN2 produces a high-scoring sequential assignment of the multidimensional signals, using a score function that rewards good connections (i.e., agreement between relevant sets of chemical shifts in different signal lists) and penalizes bad connections, unassigned signals, and assignment gaps. Examination of a set of high-scoring assignments from a large number of independent runs allows one to determine whether a unique assignment exists for the entire sequence or parts thereof. We demonstrate the MCSA algorithm using two-dimensional (2D) and three-dimensional (3D) solid state NMR spectra of several model protein samples (α-spectrin SH3 domain and protein G/B1 microcrystals, HET-s_218–289 fibrils), obtained with magic-angle spinning and standard polarization transfer techniques. The MCSA algorithm and MCASSIGN2 program can accommodate arbitrary combinations of NMR spectra with arbitrary dimensionality, and can therefore be applied in many areas of solid state and solution NMR.     

Prediction of conformation of rat galanin in the presence and absence of water with the use of monte Carlo methods and the ECEPP/3 force field

The conformation of the 29-residue rat galanin neuropeptide was studied using the Monte Carlo with energy minimization (MCM) and electrostatically driven Monte Carlo (EDMC) methods. According to a previously elaborated procedure, the polypeptide chain was first treated in a united-residue approximation, in order to enable extensive exploration of the conformational space to be carried out (with the use of MCM), Then the low-energy united-residue conformations were converted to the all-atom representations, and EDMC simulations were carried out for the all-atom polypeptide chains, using the ECEPP/3 force field with hydration included. In order to estimate the effect of environment on galanin conformation, the low-energy conformations obtained as a result of these simulations were taken as starting structures for further EDMC runs that did not include hydration. The lowest-energy conformation obtained in aqueous solution calculations had a nonhelical N-terminal part packed against the nonpolar face of a residual helix that extended from Pro¹³ toward the C-terminus. One next lowest-energy structure was a nearly-all-helical conformation, but with a markedly higher energy. In contrast, all of the low-energy conformations in the absence of water were all-helical differing only by the extent to which the helix was kinked around Pro¹³. These results are in qualitative agreement with the available NMR and CD data of galanin in aqueous and nonaqueous solvents.     

Folding simulations of small proteins

Understanding how a protein folds is a long-standing challenge in modern science. We have used an optimized atomistic model (united-residue force field) to simulate folding of small proteins of various structures: HP-36 (alpha protein), protein A (beta), 1fsd (alpha+beta), and betanova (beta). Extensive Monte Carlo folding simulations (ten independent runs with 10(9) Monte Carlo steps at a temperature) starting from non-native conformations are carried out for each protein. In all cases, proteins fold into their native-like conformations at appropriate temperatures, and glassy transitions occur at low temperatures. To investigate early folding trajectories, 200 independent runs with 10(6) Monte Carlo steps are also performed at a fixed temperature for a protein. There are a variety of possible pathways during non-equilibrium early processes (fast process, approximately 10(4) Monte Carlo steps). Finally, these pathways converge to the point unique for each protein. The convergence point of the early folding pathways can be determined only by direct folding simulations. The free energy surface, an equilibrium thermodynamic property, dictates the rest of the folding (slow process, approximately 10(8) Monte Carlo steps).     

Conformational analysis of an alpha-galactosyl trisaccharide epitope involved in hyperacute rejection upon xenotransplantation.

alpha-Galactosyl epitopes are carbohydrate structures bearing an alpha-Gal-(1-->3)-Gal terminus (alpha-Gal epitopes). The interaction of these epitopes on the surface of animal cells with anti alpha-Gal antibodies in human serum is believed to be the main cause in antibody-mediated hyperacute rejection in xenotransplantation. In this paper, conformational analysis of an N-linked alpha-D-Galp-(1-->3)-beta-D-Galp-(1-->4)-beta-D-Glcp trisaccharide epitope was conducted in terms of each monosaccharide residue conformation, primary hydroxymethyl group configuration, and interglycosidic conformations. Selective 2D J-delta INEPT experiments have been carried out at three different temperatures to evaluate three-bond, long-range 13C-1H coupling constants for the crucial alpha-(1-->3) linkage. The NMR experimental data were complemented by theoretical calculations. The flexibility and dynamics of the trisaccharide have been studied by Metropolis Monte Carlo simulations. Ensemble-averaged three-bond, long-range 13C-1H coupling constants and nuclear Overhauser effects were in good agreement with the experimental data. The alpha-(1-->3) glycosidic linkage has shown a restricted flexibility as indicated by NMR spectroscopy and molecular modeling.     

Simulations on the Primitive Electrolyte Environment of a High Charge-density Polyelectrolyte. A Sampling Problem and its Solution

Abstract

We show that the classical Metropolis Monte Carlo (MMC) algorithm converges very slowly when applied to the primitive electrolyte environment for a high charge-density polyelectrolyte. This slowness of convergence, which is due to the large density inhomogeneity around the polyelectrolyte, produces noticeable errors in the ion distribution functions for MMC runs of 1.3 × 10⁶ trial steps started from nonequilibrium distributions. We report that an algorithm which we call DSMC (for density-scaled Monte Carlo) overcomes this problem and provides relatively rapid convergence in this application. We suggest that DSMC should be well-suited for other Monte Carlo simulations on physical systems where large density inhomogeneities occur.     

The loop problem in proteins: A monte carlo simulated annealing approach

A Monte Carlo simulated annealing (MCSA) algorithm was used to generate the conformations of local regions in bovine pancreatic trypsin inhibitor (BPTI) starting from random initial conformations. In the approach explored, only the conformation of the segment is computed; the rest of the protein is fixed in the known native conformation. Rather than follow a single simulation exhaustively, computer time is better used by performing multiple independent MCSA simulations in which different starting temperatures are employed and the number of conformations sampled is varied. The best computed conformation is chosen on the basis of lowest total energy and refined further. The total energy used in the annealing is the sum of the intrasegment energy, the interaction energy of the segment with the local surrounding region, and a distance constraint to generate a smooth connection of the initially randomized segment with the rest of the protein. The rms deviations between the main-chain conformations of the computed segments in BPTI and those of the native x-ray structure are 0.94 Å for a 5-residue α-helical segment, 1.11 Å for a 5-residue β-strand segment, and 1.03, 1.61, and 1.87 Ã for 5-, 7-, and 9-residue loop segments. Side-chain deviations are comparable to the main-chain deviations for those side chains that interact strongly with the fixed part of the protein. A detailed view of the deviations at an atom-resolved level is obtained by comparing the predicted segments with their known conformations in the crystal structure of BPTI. These results emphasize the value of predetermined fixed structure against which the computed segment can nest. © 1993 John Wiley & Sons, Inc.     

Synthesis and conformational analysis of a pentasaccharide corresponding to the cell-wall polysaccharide of the Group A Streptococcus

The synthesis and conformational analysis of a pentasaccharide corresponding to a fragment of the cell-wall polysaccharide (CWPS) of the bacteria Streptococcus Group A are described. The polysaccharide consists of alternating alpha-(1 --> 2)- and alpha-(1 --> 3)-linked L-rhamnopyranose (Rhap) residues with branching 2-acetamido-2-deoxy-D-glucopyranose (GlcpNAc) residues linked beta-(1 --> 3) to alternate rhamnose rings. The pentasaccharide is of interest as a possible terminal unit on the CWPS, for use in a vaccine. The syntheses employed a trichloroacetimidate glycosyl donor. Molecular dynamics (MD) calculations of the pentasaccharide with the force fields CVFF and PARM22, both in gas phase and with explicit water present, gave different predictions for the flexibility and preferred conformational space. Metropolis Monte Carlo (MMC) calculations with the HSEA force field were also performed. Experimental data were obtained from 1D transient NOE measurements. Complete build-up curves were compared to those obtained by full relaxation matrix calculations in order to derive a model of the conformation. Overall, the best fit between experimental and calculated data was obtained with MMC simulations using the HSEA force field. Molecular dynamics and MMC simulations of a tetrasaccharide corresponding to the Group A-variant polysaccharide, which differs in structure from Group A in lacking the GlcpNAc residues, were also performed for purposes of comparison.     

Determination of the conformation of folding initiation sites in proteins by computer simulation

Experimental evidence and theoretical models both suggest that protein folding begins by specific short regions of the polypeptide chain intermittently assuming conformations close to their final ones. The independent folding properties and small size of these folding initiation sites make them suitable subjects for computational methods aimed at deriving structure from sequence. We have used a torsion space Monte Carlo procedure together with an all-atom free energy function to investigate the folding of a set of such sites. The free energy function is derived by a potential of mean force analysis of experimental protein structures. The most important contributions to the total free energy are the local main chain electrostatics, main chain hydrogen bonds, and the burial of nonpolar area. Six proposed independent folding units and four control peptides 11–14 residues long have been investigated. Thirty Monte Carlo simulations were performed on each peptide, starting from different random conformations. Five of the six folding units adopted conformations close to the experimental ones in some of the runs. None of the controls did so, as expected. The generated conformations which are close to the experimental ones have among the lowest free energies encountered, although some less native like low free energy conformations were also found. The effectiveness of the method on these peptides, which have a wide variety of experimental conformations, is encouraging in two ways: First, it provides independent evidence that these regions of the sequences are able to adopt native like conformations early in folding, and therefore are most probably key components of the folding pathways. Second, it demonstrates that available simulation methods and free energy functions are able to produce reasonably accurate structures. Extensions of the methods to the folding of larger portions of proteins are suggested. © 1995 Wiley-Liss, Inc.     

The hydration shell of myoglobin

The space in the unit cell of a metmyoglobin crystal not occupied by myoglobin atoms was filled with water using Monte Carlo calculations. Independent calculations with different amounts of water have been performed. Structure factors were calculated using the water coordinates thus obtained and the known coordinates of the myoglobin atoms. A comparison with experimental structure factors showed that both the low and the high resolution regime could be well reproduced with 814 Monte Carlo water molecules per unit cell with a B-value of 50 Ų. The Monte Carlo water molecules yield a smaller standard R-value (0.166) than using a homogeneous electron density for the simulation of the crystal water (R = 0.212). A reciprocal space refinement of the water and the protein coordinates has been performed. Monte Carlo calculations can be used to obtain information for crystallographically invisible parts of the unit cell and yield better coordinates for the visible part in the refinement. Correspondence to: F. Parak     

Conformational analysis of a 12-residue analogue of mastoparan and of mastoparan X.

We have investigated the conformational properties of a truncated analogue of mastoparan and of mastoparan X, both peptides from wasp venom. The electrostatically driven Monte Carlo method was used to explore the conformational space of these short peptides. The initial conformations used in this study, mainly random ones, led to alpha-helical conformations. The alpha-helical conformations thus found exhibit an amphipathic character. These results are in accord with experimental data from NMR and CD spectroscopy.     

Ureido group containing cyclic dermorphin(1-7) analogues: synthesis, biology and conformation.

Six cyclic peptides related to dermorphin(1-7) have been synthesized. The synthesis of linear peptides containing diamino acid residues in positions 2 and 4 was carried out on a 4-methylbenzhydrylamine resin, and cyclization was achieved by treatment with bis-(4-nitrophenyl)carbonate to form a urea unit. The peptides were tested in the guinea-pig ileum (GPI) and mouse vas deferens (MVD) assays. Diverse opioid agonist activities were observed, depending on the size of the ring. The results were compared with those obtained earlier for 1-4 dermorphin analogues. The conformations of all six dermorphin analogues were studied. The conformational space of the peptides was examined using the electrostatically driven Monte Carlo method. On the basis of NMR data, an ensemble of conformations was obtained for each peptide. The opioid activity profiles of the compounds are discussed in the light of the structural data.     

Study of bradykinin conformation in the presence of model membrane by Nuclear Magnetic Resonance and molecular modelling

The conformation of bradykinin (BK), Arg1-Pro2-Pro3-Gly4-Phe5-Ser6-Pro7-Phe8-Arg9, was investigated by Nuclear Magnetic Resonance (NMR) spectroscopy and Monte Carlo simulation in two different media, i.e. in pure aqueous solution and in the presence of phospholipid vesicles. Monolamellar liposomes are a good model for biological membranes and mimic the environment experienced by bradykinin when interacting with G-protein coupled receptors (GPCRs). The NMR spectra showed that lipid bilayers induced a secondary structure in the otherwise inherently flexible peptide. The results of ensemble calculations revealed conformational changes occurring rapidly on the NMR time scale and allowed for the identification of different families of conformations that were averaged to reproduce the NMR observables. These structural results supported the hypothesis of the central role played by the peptide C-terminal domain in biological environments, and provided an explanation for the different biological behaviours observed for bradykinin.     

Monte Carlo loop refinement and virtual screening of the thyroid-stimulating hormone receptor transmembrane domain

Metropolis Monte Carlo (MMC) loop refinement has been performed on the three extracellular loops (ECLs) of rhodopsin and opsin-based homology models of the thyroid-stimulating hormone receptor transmembrane domain, a class A type G protein-coupled receptor. The Monte Carlo sampling technique, employing torsion angles of amino acid side chains and local moves for the six consecutive backbone torsion angles, has previously reproduced the conformation of several loops with known crystal structures with accuracy consistently less than 2?Å. A grid-based potential map, which includes van der Waals, electrostatics, hydrophobic as well as hydrogen-bond potentials for bulk protein environment and the solvation effect, has been used to significantly reduce the computational cost of energy evaluation. A modified sigmoidal distance-dependent dielectric function has been implemented in conjunction with the desolvation and hydrogen-bonding terms. A long high-temperature simulation with 2?kcal/mol repulsion potential resulted in extensive sampling of the conformational space. The slow annealing leading to the low-energy structures predicted secondary structure by the MMC technique. Molecular docking with the reported agonist reproduced the binding site within 1.5?Å. Virtual screening performed on the three lowest structures showed that the ligand-binding mode in the inter-helical region is dependent on the ECL conformations.     

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     

Solution Structures of Rat Amylin Peptide: Simulation, Theory, and Experiment

Amyloid deposits of amylin in the pancreas are an important characteristic feature found in patients with Type-2 diabetes. The aggregate has been considered important in the disease pathology and has been studied extensively. However, the secondary structures of the individual peptide have not been clearly identified. In this work, we present detailed solution structures of rat amylin using a combination of Monte Carlo and molecular dynamics simulations. A new Monte Carlo method is presented to determine the free energy of distinct biomolecular conformations. Both folded and random-coil conformations of rat amylin are observed in water and their relative stability is examined in detail. The former contains an α-helical segment comprised of residues 7-17. We find that at room temperature the folded structure is more stable, whereas at higher temperatures the random-coil structure predominates. From the configurations and weights we calculate the α-carbon NMR chemical shifts, with results that are in reasonable agreement with experiments of others. We also calculate the infrared spectrum in the amide I stretch regime, and the results are in fair agreement with the experimental line shape presented herein.     

Conformational analysis of perezone and dihydroperezone using vibrational circular dichroism

The vibrational circular dichroism (VCD) spectra of perezone and dihydroperezone measured from CDCl₃ solutions were quite similar, suggesting analogous conformations for both molecules. Their absolute configurations were confirmed by comparison of the experimental VCD spectrum of each compound with curves generated from theoretical calculations using density functional theory (DFT) at the B3LYP/DGDZVP level of theory taking into account their conformational mobility. Conformational analysis of the 8-(R) enantiomer showed 19 low energy conformers in a 2.4 kcal/mol energy range, while for 8-(R), with the saturated side alkyl chain, 34 conformers were considered in the first 2 kcal/mol. Initial analyses were carried out using a Monte Carlo searching with the MMFF94 molecular mechanics force field, all MMFF94 conformers were geometrically optimized using DFT at the B3LYP/6-31G(d) level of theory, followed by reoptimization and calculations of their vibrational frequencies at the B3LYP/DGDZVP level. Good agreement between the theoretical 8-(R) enantiomers and experimental VCD curves were observed for both.     

Computer simulations of de novo designed helical proteins.

In the context of reduced protein models, Monte Carlo simulations of three de novo designed helical proteins (four-member helical bundle) were performed. At low temperatures, for all proteins under consideration, protein-like folds having different topologies were obtained from random starting conformations. These simulations are consistent with experimental evidence indicating that these de novo designed proteins have the features of a molten globule state. The results of Monte Carlo simulations suggest that these molecules adopt four-helix bundle topologies. They also give insight into the possible mechanism of folding and association, which occurs in these simulations by on-site assembly of the helices. The low-temperature conformations of all three sequences have the features of a molten globule state.     

Structure and electronic properties of azadirachtin

We performed a combined DFT and Monte Carlo ¹³C NMR chemical-shift study of azadirachtin A, a triterpenoid that acts as a natural insect antifeedant. A conformational search using a Monte Carlo technique based on the RM1 semiempirical method was carried out in order to establish its preferred structure. The B3LYP/6-311++G(d,p), wB97XD/6-311++G(d,p), M06/6-311++G(d,p), M06-2X/6-311++G(d,p), and CAM-B3LYP/6-311++G(d,p) levels of theory were used to predict NMR chemical shifts. A Monte Carlo population-weighted average spectrum was produced based on the predicted Boltzmann contributions. In general, good agreement between experimental and theoretical data was obtained using both methods, and the ¹³C NMR chemical shifts were predicted highly accurately. The geometry was optimized at the semiempirical level and used to calculate the NMR chemical shifts at the DFT level, and these shifts showed only minor deviations from those obtained following structural optimization at the DFT level, and incurred a much lower computational cost. The theoretical ultraviolet spectrum showed a maximum absorption peak that was mainly contributed by the tiglate group.     

Metropolis Monte Carlo calculations of DNA structure using internal coordinates and NMR distance restraints: An alternative method for generating a high-resolution solution structure

Summary A new method, a restrained Monte Carlo (rMC) calculation, is demonstrated for generating high-resolution structures of DNA oligonucleotides in solution from interproton distance restraints and bounds derived from complete relaxation matrix analysis of two-dimensional nuclear Overhauser effect (NOE) spectral peak intensities. As in the case of restrained molecular dynamics (rMD) refinement of structures, the experimental distance restraints and bounds are incorporated as a pseudo-energy term (or penalty function) into the mathematical expression for the molecular energy. However, the use of generalized helical parameters, rather than Cartesian coordinates, to define DNA conformation increases efficiency by decreasing by an order of magnitude the number of parameters needed to describe a conformation and by simplifying the potential energy profile. The Metropolis Monte Carlo method is employed to simulate an annealing process. The rMC method was applied to experimental 2D NOE data from the octamer duplex d(GTA-TAATG)·d(CATTATAC). Using starting structures from different locations in conformational space (e.g. A-DNA and B-DNA), the rMC calculations readily converged, with a root-mean-square deviation (RMSD) of <0.3 Å between structures generated using different protocols and starting structures. Theoretical 2D NOE peak intensities were calculated for the rMC-generated structures using the complete relaxation matrix program CORMA, enabling a comparison with experimental intensities via residual indices. Simulation of the vicinal proton coupling constants was carried out for the structures generated, enabling a comparison with the experimental deoxyribose ring coupling constants, which were not utilized in the structure determination in the case of the rMC simulations. Agreement with experimental 2D NOE and scalar coupling data was good in all cases. The rMC structures are quite similar to that refined by a traditional restrained MD approach (RMSD<0.5 Å) despite the different force fields used and despite the fact that MD refinement was conducted with additional restraints imposed on the endocyclic torsion angles of deoxyriboses. The computational time required for the rMC and rMD calculations is about the same. A comparison of structural parameters is made and some limitations of both methods are discussed with regard to the average nature of the experimental restraints used in the refinement.Abbreviations MC Monte Carlo - rMC restrained Monte Carlo - MD molecular dynamics - rMD restrained molecular dynamics - DG distance geometry - EM energy minimization - 2D NOE two-dimensional nuclear Overhauser effect - DQF-COSY double-quantum-filtered correlation spectroscopy - RMSD root-mean-square deviation To whom correspondence should be addressed.