Monte Carlo study of cooperativity in homopolypeptides

©1992 American Institute of PhysicsThe electronic version of this article is the complete one and can be found online at: http://link.aip.org/link/?JCPSA6/97/9412/1DOI:10.1063/1.463317A discretized model of globular proteins is employed in a Monte Carlo study of the helix-coil transition of polyalanine and the collapse transition of polyvaline. The present lattice realization permits real protein crystal structures to be represented at the level of 1 A resolution. Furthermore, the Monte Carlo dynamic scheme is capable of moving elements of assembled secondary and supersecondary structure. The potentials of mean force for the interactions are constructed from the statistics of a set of high resolution x-ray structures of nonhomologous proteins. The cooperativity of formation of ordered structures is found to be larger when the major contributions to the conformational energy of the low temperature states come from hydrogen bonds and short range conformational propensities. The secondary structure seen in the folded state is the result of an interplay between the short and long range interactions. Compactness itself, driven by long range, nonspecific interactions, seems to be insufficient to generate any appreciable secondary structure. A detailed examination of the dynamics of highly helical model proteins demonstrates that all elements of secondary structure are mobile in the present algorithm, and thus the folding pathways do not depend on the use of a lattice approximation. Possible applications of the present model to the prediction of protein 3D structures are briefly discussed

CABS-flex predictions of protein flexibility compared with NMR ensembles

Motivation: Identification of flexible regions of protein structures is important for understanding of their biological functions. Recently, we have developed a fast approach for predicting protein structure fluctuations from a single protein model: the CABS-flex. CABS-flex was shown to be an efficient alternative to conventional all-atom molecular dynamics (MD). In this work, we evaluate CABS-flex and MD predictions by comparison with protein structural variations within NMR ensembles. Results: Based on a benchmark set of 140 proteins, we show that the relative fluctuations of protein residues obtained from CABS-flex are well correlated to those of NMR ensembles. On average, this correlation is stronger than that between MD and NMR ensembles. In conclusion, CABS-flex is useful and complementary to MD in predicting of protein regions that undergo conformational changes and the extent of such changes

Monte Carlo dynamics of diamond-lattice multichain systems

©1986 American Institute of Physics. The electronic version of this article is the complete one and can be found online at: http://link.aip.org/link/?APCPCS/137/241/1DOI:10.1063/1.35530Presented at the 1985 La Jolla Workshop on Polymer Flow Interaction.We present preliminary results of Monte Carlo studies on the dynamics of multichain diamond-lattice systems at considerably greater densities than those done previously. Chain dynamics were simulated by a random sequence of three or four bond kink motions. The single bead autocorrelation function exhibits "slow" mode relaxation behavior with a g(t)∝ tβ. There is a smooth crossover from Rouse-like dynamics, β=1/2, at low density to smaller values of β at higher density and β=0 at the glass transition density (φG≅0.92). The simulation provides a self-diffusion coefficient D ∝ n-2, with n the number of beads, in agreement with experiment. A phenomenological model, different from the widely accepted reptation picture, is proposed

Numerical estimation of entropy loss on dimerization: improved prediction of the quaternary structure of the GCN4 leucine zipper

A lattice based model of a protein is used to study the dimerization equilibrium of the GCN4 leucine zipper. Replica exchange Monte Carlo is used to determine the free energy of both the monomeric and dimeric forms as a function of temperature. The method of coincidences is then introduced to explicitly calculate the entropy loss associated with dimerization, and from it the free energy difference between monomer and dimer, as well as the corresponding equilibrium reaction constant. We find that the entropy loss of dimerization is a strong function of energy (or temperature), and that it is much larger than previously estimated, especially for high energy states. The results confirm that it is possible to study the dimerization equilibrium of GCN4 at physiological concentrations within the reduced representation of the protein employed

Contact prediction in protein modeling: Scoring, folding and refinement of coarse-grained models

<p>Abstract</p> <p>Background</p> <p>Several different methods for contact prediction succeeded within the Sixth Critical Assessment of Techniques for Protein Structure Prediction (CASP6). The most relevant were non-local contact predictions for targets from the most difficult categories: fold recognition-analogy and new fold. Such contacts could provide valuable structural information in case a template structure cannot be found in the PDB.</p> <p>Results</p> <p>We described comprehensive tests of the effectiveness of contact data in various aspects of de novo modeling with CABS, an algorithm which was used successfully in CASP6 by the Kolinski-Bujnicki group. We used the predicted contacts in a simple scoring function for the post-simulation ranking of protein models and as a soft bias in the folding simulations and in the fold-refinement procedure. The latter approach turned out to be the most successful. The CABS force field used in the Replica Exchange Monte Carlo simulations cooperated with the true contacts and discriminated the false ones, which resulted in an improvement of the majority of Kolinski-Bujnicki's protein models. In the modeling we tested different sets of predicted contact data submitted to the CASP6 server. According to our results, the best performing were the contacts with the accuracy balanced with the coverage, obtained either from the best two predictors only or by a consensus from as many predictors as possible.</p> <p>Conclusion</p> <p>Our tests have shown that theoretically predicted contacts can be very beneficial for protein structure prediction. Depending on the protein modeling method, a contact data set applied should be prepared with differently balanced coverage and accuracy of predicted contacts. Namely, high coverage of contact data is important for the model ranking and high accuracy for the folding simulations.</p

Effect of double bonds on the dynamics of hydrocarbon chains

©1992 American Institute of PhysicsThe electronic version of this article is the complete one and can be found online at: http://link.aip.org/link/?JCPSA6/97/1240/1DOI:10.1063/1.463250Brownian dynamics simulations of isolated 18-carbon chains have been performed, both for saturated and unsaturated hydrocarbons. The effect of one or several (nonconjugated) double bonds on the properties of the chains is discussed in terms of both equilibrium and dynamic properties. The introduction of a cis double bond increases the relaxation rates of the unsaturated chain with respect to the saturated alkane. On the other hand, coupling effects in the torsional transitions around a trans double bond make the dynamics of this unsaturated chain very similar to the saturated one. Based on these results, the parameters and moves of a dynamic Monte Carlo algorithm are tuned to reproduce the observed behavior, providing an efficient method for the study of more complicated systems

Combining Coarse-Grained Protein Models with Replica-Exchange All-Atom Molecular Dynamics

We describe a combination of all-atom simulations with CABS, a well-established coarse-grained protein modeling tool, into a single multiscale protocol. The simulation method has been tested on the C-terminal beta hairpin of protein G, a model system of protein folding. After reconstructing atomistic details, conformations derived from the CABS simulation were subjected to replica-exchange molecular dynamics simulations with OPLS-AA and AMBER99sb force fields in explicit solvent. Such a combination accelerates system convergence several times in comparison with all-atom simulations starting from the extended chain conformation, demonstrated by the analysis of melting curves, the number of native-like conformations as a function of time and secondary structure propagation. The results strongly suggest that the proposed multiscale method could be an efficient and accurate tool for high-resolution studies of protein folding dynamics in larger systems.Comment: 12 pages, 4 figure

A hierarchical approach to the prediction of the quaternary structure of GCN4 and its mutants

First published in DIMACS Series in Discrete Mathematics and Theoretical Computer Science, 23 (1996) published by the American Mathematical Society.Presented at DIMACS Workshop on Global Minimization of Nonconvex Energy Functions: Molecular Conformation and Protein Folding, March 20-21, 1995.A hierarchical approach to protein folding is employed to examine the folding pathway and predict the quaternary structure of the GCN4 leucine zipper. Structures comparable in quality to experiment have been predicted. In addition, the equilibrium between dimers, trimers and tetramers of a number of GCN4 mutants has been examined. In five out of eight cases, the simulation results are in accordance with the experimental studies of Harbury, et al