Browse
Submit Data
Databases
API
Help

Dataset Information

26 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Hierarchical natural move Monte Carlo refines flexible RNA structures into cryo-EM densities.

ABSTRACT: 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.

SUBMITTER: Chang JY

PROVIDER: S-EPMC7668250 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Exploring peptide/MHC detachment processes using hierarchical natural move Monte Carlo.

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

ChemEM: Flexible Docking of Small Molecules in Cryo-EM Structures.

Project description:Cryo-electron microscopy (cryo-EM), through resolution advancements, has become pivotal in structure-based drug discovery. However, most cryo-EM structures are solved at 3-4 Å resolution, posing challenges for small-molecule docking and structure-based virtual screening due to issues in the precise positioning of ligands and the surrounding side chains. We present ChemEM, a software package that employs cryo-EM data for the accurate docking of one or multiple ligands in a protein-binding site. Validated against a highly curated benchmark of high- and medium-resolution cryo-EM structures and the corresponding high-resolution controls, ChemEM displayed impressive performance, accurately placing ligands in all but one case, often surpassing cryo-EM PDB-deposited solutions. Even without including the cryo-EM density, the ChemEM scoring function outperformed the well-established AutoDock Vina score. Using ChemEM, we illustrate that valuable information can be extracted from maps at medium resolution and underline the utility of cryo-EM structures for drug discovery.

| S-EPMC10788898 | biostudies-literature

Prediction of protein loop structures using a local move Monte Carlo approach and a grid-based force field.

Project description:We have developed an improved local move Monte Carlo (LMMC) loop sampling approach for loop predictions. The method generates loop conformations based on simple moves of the torsion angles of side chains and local moves of backbone of loops. To reduce the computational costs for energy evaluations, we developed a grid-based force field to represent the protein environment and solvation effect. Simulated annealing has been used to enhance the efficiency of the LMMC loop sampling and identify low-energy loop conformations. The prediction quality is evaluated on a set of protein loops with known crystal structure that has been previously used by others to test different loop prediction methods. The results show that this approach can reproduce the experimental results with the root mean square deviation within 1.8 A for all the test cases. The LMMC loop prediction approach developed here could be useful for improvement in the quality the loop regions in homology models, flexible protein-ligand and protein-protein docking studies.

| S-EPMC2597363 | biostudies-literature

Comparing cryo-EM structures.

Project description: Not available

| S-EPMC6464812 | biostudies-literature

Assessing protein loop flexibility by hierarchical Monte Carlo sampling.

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

Flexible architecture of IP3R1 by Cryo-EM.

Project description:Inositol 1,4,5-trisphosphate receptors (IP3Rs) play a fundamental role in generating Ca2+ signals that trigger many cellular processes in virtually all eukaryotic cells. Thus far, the three-dimensional (3D) structure of these channels has remained extremely controversial. Here, we report a subnanometer resolution electron cryomicroscopy (cryo-EM) structure of a fully functional type 1 IP3R from cerebellum in the closed state. The transmembrane region reveals a twisted bundle of four α helices, one from each subunit, that form a funnel shaped structure around the 4-fold symmetry axis, strikingly similar to the ion-conduction pore of K+ channels. The lumenal face of IP3R1 has prominent densities that surround the pore entrance and similar to the highly structured turrets of Kir channels. 3D statistical analysis of the cryo-EM density map identifies high variance in the cytoplasmic region. This structural variation could be attributed to genuine structural flexibility of IP3R1.

| S-EPMC3154621 | biostudies-literature

Cryo-EM structures of coagulation factors.

Project description:A State of the Art lecture titled "Cryo-EM structures of coagulation factors" was presented at the ISTH Congress in 2022. Cryogenic electron microscopy (cryo-EM) is a revolutionary technique capable of solving the structure of high molecular weight proteins and their complexes, unlike nuclear magnetic resonance (NMR), and under conditions not biased by crystal contacts, unlike X-ray crystallography. These features are particularly relevant to the analysis of coagulation factors that are too big for NMR and often recalcitrant to X-ray investigation. Using cryo-EM, we have solved the structures of coagulation factors V and Va, prothrombinase on nanodiscs, and the prothrombin-prothrombinase complex. These structures have advanced basic knowledge in the field of thrombosis and hemostasis, especially on the function of factor V and the molecular mechanism for prothrombin activation, and set the stage for exciting new lines of investigation. Finally, we summarize relevant new data on this topic presented during the 2022 ISTH Congress.

| S-EPMC9630041 | biostudies-literature

Support Vector Regression-Based Monte Carlo Simulation of Flexible Water Clusters.

Project description:Molecular simulations based on classical force fields are computationally efficient but lack accuracy due to the empirical formulation of non-bonded interactions. Quantum mechanical (QM) methods, albeit accurate, have inhibitory computational costs for large molecules and clusters. Hence, to overcome the bottleneck, machine learning (ML)-based methods have been employed in the recent years. We had earlier reported a combined scheme of many-body expansion (MBE) and ML to predict the interaction energies of rigid water clusters. In this work, we proceed toward building a flexible water model using the ML-MBE scheme. This ML-MBE scheme has an error of <1% for interaction energy prediction in comparison to the parent QM method for flexible water decamers. Machine learning-based Monte Carlo simulations (MLMC) are performed with this water model, and the structural properties of these configurations are compared with those obtained from ab initio molecular dynamics (AIMD) and the TIP3P classical force field. The radial distribution functions, tetrahedral order parameters, and number of hydrogen bonds in AIMD and MLMC have a similar qualitative and quantitative trend, whereas the classical force fields show a significant deviation.

| S-EPMC7143414 | biostudies-literature

Cryo_fit: Democratization of flexible fitting for cryo-EM.

Project description:Cryo-electron microscopy (cryo-EM) is becoming a method of choice for describing native conformations of biomolecular complexes at high resolution. The rapid growth of cryo-EM in recent years has created a high demand for automated solutions, both in hardware and software. Flexible fitting of atomic models to three-dimensional (3D) cryo-EM reconstructions by molecular dynamics (MD) simulation is a popular technique but often requires technical expertise in computer simulation. This work introduces cryo_fit, a package for the automatic flexible fitting of atomic models in cryo-EM maps using MD simulation. The package is integrated with the Phenix software suite. The module was designed to automate the multiple steps of MD simulation in a reproducible manner, as well as facilitate refinement and validation through Phenix. Through the use of cryo_fit, scientists with little experience in MD simulation can produce high quality atomic models automatically and better exploit the potential of cryo-EM.

| S-EPMC7112765 | biostudies-literature

Markov chain Monte Carlo methods for hierarchical clustering of dynamic causal models.

Project description:In this article, we address technical difficulties that arise when applying Markov chain Monte Carlo (MCMC) to hierarchical models designed to perform clustering in the space of latent parameters of subject-wise generative models. Specifically, we focus on the case where the subject-wise generative model is a dynamic causal model (DCM) for functional magnetic resonance imaging (fMRI) and clusters are defined in terms of effective brain connectivity. While an attractive approach for detecting mechanistically interpretable subgroups in heterogeneous populations, inverting such a hierarchical model represents a particularly challenging case, since DCM is often characterized by high posterior correlations between its parameters. In this context, standard MCMC schemes exhibit poor performance and extremely slow convergence. In this article, we investigate the properties of hierarchical clustering which lead to the observed failure of standard MCMC schemes and propose a solution designed to improve convergence but preserve computational complexity. Specifically, we introduce a class of proposal distributions which aims to capture the interdependencies between the parameters of the clustering and subject-wise generative models and helps to reduce random walk behaviour of the MCMC scheme. Critically, these proposal distributions only introduce a single hyperparameter that needs to be tuned to achieve good performance. For validation, we apply our proposed solution to synthetic and real-world datasets and also compare it, in terms of computational complexity and performance, to Hamiltonian Monte Carlo (HMC), a state-of-the-art Monte Carlo technique. Our results indicate that, for the specific application domain considered here, our proposed solution shows good convergence performance and superior runtime compared to HMC.

| S-EPMC8193526 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data