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

Towards the high-resolution protein structure prediction. Fast refinement of reduced models with all-atom force field

<p>Abstract</p> <p>Background</p> <p>Although experimental methods for determining protein structure are providing high resolution structures, they cannot keep the pace at which amino acid sequences are resolved on the scale of entire genomes. For a considerable fraction of proteins whose structures will not be determined experimentally, computational methods can provide valuable information. The value of structural models in biological research depends critically on their quality. Development of high-accuracy computational methods that reliably generate near-experimental quality structural models is an important, unsolved problem in the protein structure modeling.</p> <p>Results</p> <p>Large sets of structural decoys have been generated using reduced conformational space protein modeling tool CABS. Subsequently, the reduced models were subject to all-atom reconstruction. Then, the resulting detailed models were energy-minimized using state-of-the-art all-atom force field, assuming fixed positions of the alpha carbons. It has been shown that a very short minimization leads to the proper ranking of the quality of the models (distance from the native structure), when the all-atom energy is used as the ranking criterion. Additionally, we performed test on medium and low accuracy decoys built via classical methods of comparative modeling. The test placed our model evaluation procedure among the state-of-the-art protein model assessment methods.</p> <p>Conclusion</p> <p>These test computations show that a large scale high resolution protein structure prediction is possible, not only for small but also for large protein domains, and that it should be based on a hierarchical approach to the modeling protocol. We employed Molecular Mechanics with fixed alpha carbons to rank-order the all-atom models built on the scaffolds of the reduced models. Our tests show that a physic-based approach, usually considered computationally too demanding for large-scale applications, can be effectively used in such studies.</p

Type II restriction endonuclease R.Eco29kI is a member of the GIY-YIG nuclease superfamily

<p>Abstract</p> <p>Background</p> <p>The majority of experimentally determined crystal structures of Type II restriction endonucleases (REases) exhibit a common PD-(D/E)XK fold. Crystal structures have been also determined for single representatives of two other folds: PLD (R.BfiI) and half-pipe (R.PabI), and bioinformatics analyses supported by mutagenesis suggested that some REases belong to the HNH fold. Our previous bioinformatic analysis suggested that REase R.Eco29kI shares sequence similarities with one more unrelated nuclease superfamily, GIY-YIG, however so far no experimental data were available to support this prediction. The determination of a crystal structure of the GIY-YIG domain of homing endonuclease I-TevI provided a template for modeling of R.Eco29kI and prompted us to validate the model experimentally.</p> <p>Results</p> <p>Using protein fold-recognition methods we generated a new alignment between R.Eco29kI and I-TevI, which suggested a reassignment of one of the putative catalytic residues. A theoretical model of R.Eco29kI was constructed to illustrate its predicted three-dimensional fold and organization of the active site, comprising amino acid residues Y49, Y76, R104, H108, E142, and N154. A series of mutants was constructed to generate amino acid substitutions of selected residues (Y49A, R104A, H108F, E142A and N154L) and the mutant proteins were examined for their ability to bind the DNA containing the Eco29kI site 5'-CCGCGG-3' and to catalyze the cleavage reaction. Experimental data reveal that residues Y49, R104, E142, H108, and N154 are important for the nuclease activity of R.Eco29kI, while H108 and N154 are also important for specific DNA binding by this enzyme.</p> <p>Conclusion</p> <p>Substitutions of residues Y49, R104, H108, E142 and N154 predicted by the model to be a part of the active site lead to mutant proteins with strong defects in the REase activity. These results are in very good agreement with the structural model presented in this work and with our prediction that R.Eco29kI belongs to the GIY-YIG superfamily of nucleases. Our study provides the first experimental evidence for a Type IIP REase that does not belong to the PD-(D/E)XK or HNH superfamilies of nucleases, and is instead a member of the unrelated GIY-YIG superfamily.</p

Comparison of α-Helix and β-Sheet Structure Adaptation to a Quantum Dot Geometry: Toward the Identification of an Optimal Motif for a Protein Nanoparticle Cover

While quantum dots (QDs) are useful as fluorescent labels, their application in biosciences is limited due to the stability and hydrophobicity of their surface. In this study, we tested two types of proteins for use as a cover for spherical QDs, composed of cadmium selenide. Pumilio homology domain (Puf), which is mostly α-helical, and leucine-rich repeat (LRR) domain, which is rich in β-sheets, were selected to determine if there is a preference for one of these secondary structure types for nanoparticle covers. The protein sequences were optimized to improve their interaction with the surface of QDs. The solubilization of the apoproteins and their assembly with nanoparticles required the application of a detergent, which was removed in subsequent steps. Finally, only the Puf-based cover was successful enough as a QD hydrophilic cover. We showed that a single polypeptide dimer of Puf, PufPuf, can form a cover. We characterized the size and fluorescent properties of the obtained QD:protein assemblies. We showed that the secondary structure of the Puf proteins was not destroyed upon contact with the QDs. We demonstrated that these assemblies do not promote the formation of reactive oxygen species during illumination of the nanoparticles. The data represent advances in the effort to obtain a stable biocompatible cover for QDs

Comparison of three Monte Carlo conformational search strategies for a proteinlike homopolymer model: Folding thermodynamics and identification of low-energy structures

©2000 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/113/5065/1DOI:10.1063/1.1289533Entropy sampling Monte Carlo, the replica method, and the classical Metropolis scheme were applied in numerical studies of the collapse transition in a simple face-centered cubic lattice polymer. The force field of the model consists of pairwise, contact-type, long-range interactions and a short-range potential based on the β -sheet definition assumed in the model. The ability to find the lowest energy conformation by various Monte Carlo methods and the computational cost associated with each was examined. It is shown that all of the methods generally provide the same picture of the collapse transition. However, the most complete thermodynamic description of the transition derives from the results of entropy sampling Monte Carlo simulations, but this is the most time-consuming method. The replica method is shown to be the most effective and efficient in searching for the lowest energy conformation. The possible consequences of these findings for the development of simulation strategies for the folding of model proteins are discussed briefly

Protein structure prediction by tempering spatial constraints

The probability to predict correctly a protein structure can be enhanced through introduction of spatial constraints - Either from NMR experiments or from homologous structures. However, the additional constraints lead often to new local energy minima and worse sampling efficiency in simulations. In this work, we present a new parallel tempering variant that alleviates the energy barriers resulting from spatial constraints and therefore yields to an enhanced sampling in structure prediction simulations. © Springer Science+Business Media, Inc. 2005

Exploring protein energy landscapes with hierarchical clustering

In this work we present a new method for investigating local energy minima on a protein energy landscape. The Cα(CA), Cβ(B), and the center of mass of the side chain (CABS) method was employed for generating protein models, but any other method could be used instead. Cα traces from an ensemble of models are hierarchical clustered with the hierarchical clustering of protein models (HCPM) method. The efficiency of this method for sampling and analyzing energy landscapes is shown

Towards the high-resolution protein structure prediction. Fast refinement of reduced models with all-atom force field-2

<p><b>Copyright information:</b></p><p>Taken from "Towards the high-resolution protein structure prediction. Fast refinement of reduced models with all-atom force field"</p><p>http://www.biomedcentral.com/1472-6807/7/43</p><p>BMC Structural Biology 2007;7():43-43.</p><p>Published online 29 Jun 2007</p><p>PMCID:PMC1933428.</p><p></p>h abnormal high energy values resulted from structural inaccuracies (right panels). On the left panels, energies of the native structures are denoted by asterisks. The native structures were subjected to the same rebuilding procedure from the Cα-traces as that applied to the decoys. Proteins are ordered in respect to their chain lengths (Table 1) – from the smallest on top (2GR8A) to the largest (2CJPA) on the bottom