Prediction of protein loop structures using a local move Monte Carlo approach and a grid-based force field.
Ontology highlight
ABSTRACT:
SUBMITTER: Cui M
PROVIDER: S-EPMC2597363 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Ribonucleic acids (RNAs) play essential roles in living cells. Many of them fold into defined three-dimensional (3D) structures to perform functions. Recent advances in single-particle cryo-electron microscopy (cryo-EM) have enabled structure determinations of RNA to atomic resolutions. However, most RNA molecules are structurally flexible, limiting the resolution of their structures solved by cryo-EM. In modeling these molecules, several computational methods are limited by the requirement of massive computational resources and/or the low efficiency in exploring large-scale structural variations. Here we use hierarchical natural move Monte Carlo (HNMMC), which takes advantage of collective motions for groups of nucleic acid residues, to refine RNA structures into their cryo-EM maps, preserving atomic details in the models. After validating the method on a simulated density map of tRNA, we applied it to objectively obtain the model of the folding intermediate for the specificity domain of ribonuclease P from Bacillus subtilis and refine a flexible ribosomal RNA (rRNA) expansion segment from the Mycobacterium tuberculosis (Mtb) ribosome in different conformational states. Finally, we used HNMMC to model atomic details and flexibility for two distinct conformations of the complete genomic RNA (gRNA) inside MS2, a single-stranded RNA virus, revealing multiple pathways for its capsid assembly.
| S-EPMC7668250 | biostudies-literature
Project description:MotivationThe binding between a peptide and a major histocompatibility complex (MHC) is one of the most important processes for the induction of an adaptive immune response. Many algorithms have been developed to predict peptide/MHC (pMHC) binding. However, no approach has yet been able to give structural insight into how peptides detach from the MHC.ResultsIn this study, we used a combination of coarse graining, hierarchical natural move Monte Carlo and stochastic conformational optimization to explore the detachment processes of 32 different peptides from HLA-A*02:01. We performed 100 independent repeats of each stochastic simulation and found that the presence of experimentally known anchor amino acids affects the detachment trajectories of our peptides. Comparison with experimental binding affinity data indicates the reliability of our approach (area under the receiver operating characteristic curve 0.85). We also compared to a 1000 ns molecular dynamics simulation of a non-binding peptide (AAAKTPVIV) and HLA-A*02:01. Even in this simulation, the longest published for pMHC, the peptide does not fully detach. Our approach is orders of magnitude faster and as such allows us to explore pMHC detachment processes in a way not possible with all-atom molecular dynamics simulations.Availability and implementationThe source code is freely available for download at http://www.cs.ox.ac.uk/mosaics/.Contactbernhard.knapp@stats.ox.ac.ukSupplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC4708099 | biostudies-literature
Project description:The description of protein disordered states is important for understanding protein folding mechanisms and their functions. In this short review, we briefly describe a simulation approach to modeling protein interactions, which involve disordered peptide partners or intrinsically disordered protein regions, and unfolded states of globular proteins. It is based on the CABS coarse-grained protein model that uses a Monte Carlo (MC) sampling scheme and a knowledge-based statistical force field. We review several case studies showing that description of protein disordered states resulting from CABS simulations is consistent with experimental data. The case studies comprise investigations of protein⁻peptide binding and protein folding processes. The CABS model has been recently made available as the simulation engine of multiscale modeling tools enabling studies of protein⁻peptide docking and protein flexibility. Those tools offer customization of the modeling process, driving the conformational search using distance restraints, reconstruction of selected models to all-atom resolution, and simulation of large protein systems in a reasonable computational time. Therefore, CABS can be combined in integrative modeling pipelines incorporating experimental data and other modeling tools of various resolution.
| S-EPMC6386871 | biostudies-literature
Project description:Loop flexibility is often crucial to protein biological function in solution. We report a new Monte Carlo method for generating conformational ensembles for protein loops and cyclic peptides. The approach incorporates the triaxial loop closure method which addresses the inverse kinematic problem for generating backbone move sets that do not break the loop. Sidechains are sampled together with the backbone in a hierarchical way, making it possible to make large moves that cross energy barriers. As an initial application, we apply the method to the flexible loop in triosephosphate isomerase that caps the active site, and demonstrate that the resulting loop ensembles agree well with key observations from previous structural studies. We also demonstrate, with 3 other test cases, the ability to distinguish relatively flexible and rigid loops within the same protein.
| S-EPMC3129859 | biostudies-literature
Project description:Loops in proteins are flexible regions connecting regular secondary structures. They are often involved in protein functions through interacting with other molecules. The irregularity and flexibility of loops make their structures difficult to determine experimentally and challenging to model computationally. Conformation sampling and energy evaluation are the two key components in loop modeling. We have developed a new method for loop conformation sampling and prediction based on a chain growth sequential Monte Carlo sampling strategy, called Distance-guided Sequential chain-Growth Monte Carlo (DISGRO). With an energy function designed specifically for loops, our method can efficiently generate high quality loop conformations with low energy that are enriched with near-native loop structures. The average minimum global backbone RMSD for 1,000 conformations of 12-residue loops is 1:53 A° , with a lowest energy RMSD of 2:99 A° , and an average ensembleRMSD of 5:23 A° . A novel geometric criterion is applied to speed up calculations. The computational cost of generating 1,000 conformations for each of the x loops in a benchmark dataset is only about 10 cpu minutes for 12-residue loops, compared to ca 180 cpu minutes using the FALCm method. Test results on benchmark datasets show that DISGRO performs comparably or better than previous successful methods, while requiring far less computing time. DISGRO is especially effective in modeling longer loops (10-17 residues).
| S-EPMC3998890 | biostudies-literature
Project description:To elucidate conformational dynamics of chromatin fibers, we compared available force-spectroscopy measurements with extensive Monte Carlo simulations of nucleosome arrays under external force. Our coarse-grained model of chromatin includes phenomenological energy terms for the DNA-histone adhesion and the internucleosome stacking interactions. We found that the Monte Carlo fiber ensembles simulated with increasing degrees of DNA unwrapping and the stacking energy 8 kT can account for the intricate force-extension response observed experimentally. Our analysis shows that at low external forces (F < 3.0 picoNewtons), the DNA ends in nucleosomes breathe by ?10 bp. Importantly, under these conditions, the fiber is highly dynamic, exhibiting continuous unstacking-restacking transitions, allowing accessibility of transcription factors to DNA, while maintaining a relatively compact conformation. Of note, changing the stacking interaction by a few kT, an in silico way to mimic histone modifications, is sufficient to transform an open chromatin state into a compact fiber. The fibers are mostly two-start zigzag folds with rare occurrences of three- to five-start morphologies. The internucleosome stacking is lost during the linear response regime. At the higher forces exceeding 4 picoNewtons, the nucleosome unwrapping becomes stochastic and asymmetric, with one DNA arm opened by ?55 bp and the other arm only by ?10 bp. Importantly, this asymmetric unwrapping occurs for any kind of sequence, including the symmetric ones. Our analysis brings new, to our knowledge, insights in dynamics of chromatin modulated by histone epigenetic modifications and molecular motors such as RNA polymerase.
| S-EPMC6225046 | biostudies-literature
Project description:MHC class II molecules bind peptides derived from extracellular proteins that have been ingested by antigen-presenting cells and display them to the immune system. Peptide loading occurs within the antigen-presenting cell and is facilitated by HLA-DM. HLA-DM stabilizes the open conformation of the MHCII binding groove when no peptide is bound. While a structure of the MHCII/HLA-DM complex exists, the mechanism of stabilization is still largely unknown. Here, we applied customized Natural Move Monte Carlo to investigate this interaction. We found a possible long-range mechanism that implicates the configuration of the membrane-proximal globular domains in stabilizing the open state of the empty MHCII binding groove.
| S-EPMC7007188 | biostudies-literature
Project description:Building upon the foundational principles of the grid search algorithm and Monte Carlo numerical simulation, this article introduces an innovative epidemic monitoring and prevention plan. The plan offers the capability to accurately identify the sources of infectious diseases and predict the final scale and duration of the epidemic. The proposed plan is implemented in schools and society, utilizing computer simulation analysis. Through this analysis, the plan enables precise localization of infection sources for various demographic groups, with an error rate of less than 3%. Additionally, the plan allows for the estimation of the epidemic cycle duration, which typically spans around 14 days. Notably, higher population density enhances fault tolerance and prediction accuracy, resulting in smaller errors and more reliable simulation outcomes. Overall, this study provides highly valuable theoretical guidance for effective epidemic prevention and control efforts.
| S-EPMC10403190 | biostudies-literature
Project description:Sampling enrichment toward a target state, an analogue of the improvement of sampling efficiency (SE), is critical in both the refinement of protein structures and the generation of near-native structure ensembles for the exploration of structure-function relationships. We developed a hybrid molecular dynamics (MD)-Monte Carlo (MC) approach to enrich the sampling toward the target structures. In this approach, the higher SE is achieved by perturbing the conventional MD simulations with a MC structure-acceptance judgment, which is based on the coincidence degree of small angle x-ray scattering (SAXS) intensity profiles between the simulation structures and the target structure. We found that the hybrid simulations could significantly improve SE by making the top-ranked models much closer to the target structures both in the secondary and tertiary structures. Specifically, for the 20 mono-residue peptides, when the initial structures had the root-mean-squared deviation (RMSD) from the target structure smaller than 7 Å, the hybrid MD-MC simulations afforded, on average, 0.83 Å and 1.73 Å in RMSD closer to the target than the parallel MD simulations at 310K and 370K, respectively. Meanwhile, the average SE values are also increased by 13.2% and 15.7%. The enrichment of sampling becomes more significant when the target states are gradually detectable in the MD-MC simulations in comparison with the parallel MD simulations, and provide >200% improvement in SE. We also performed a test of the hybrid MD-MC approach in the real protein system, the results showed that the SE for 3 out of 5 real proteins are improved. Overall, this work presents an efficient way of utilizing solution SAXS to improve protein structure prediction and refinement, as well as the generation of near native structures for function annotation.
| S-EPMC4881967 | biostudies-literature
Project description:We present McGenus, an algorithm to predict RNA secondary structures with pseudoknots. The method is based on a classification of RNA structures according to their topological genus. McGenus can treat sequences of up to 1000 bases and performs an advanced stochastic search of their minimum free energy structure allowing for non-trivial pseudoknot topologies. Specifically, McGenus uses a Monte Carlo algorithm with replica exchange for minimizing a general scoring function which includes not only free energy contributions for pair stacking, loop penalties, etc. but also a phenomenological penalty for the genus of the pairing graph. The good performance of the stochastic search strategy was successfully validated against TT2NE which uses the same free energy parametrization and performs exhaustive or partially exhaustive structure search, albeit for much shorter sequences (up to 200 bases). Next, the method was applied to other RNA sets, including an extensive tmRNA database, yielding results that are competitive with existing algorithms. Finally, it is shown that McGenus highlights possible limitations in the free energy scoring function. The algorithm is available as a web server at http://ipht.cea.fr/rna/mcgenus.php.
| S-EPMC3561945 | biostudies-literature