Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

High-resolution visualisation of the states and pathways sampled in molecular dynamics simulations.

ABSTRACT: We have recently developed a scalable algorithm for ordering the instantaneous observations of a dynamical system evolving continuously in time. Here, we apply the method to long molecular dynamics trajectories. The procedure requires only a pairwise, geometrical distance as input. Suitable annotations of both structural and kinetic nature reveal the free energy basins visited by biomolecules. The profile is supplemented by a trace of the temporal evolution of the system highlighting the sequence of events. We demonstrate that the resultant SAPPHIRE (States And Pathways Projected with HIgh REsolution) plots provide a comprehensive picture of the thermodynamics and kinetics of complex, molecular systems exhibiting dynamics covering a range of time and length scales. Information on pathways connecting states and the level of recurrence are quickly inferred from the visualisation. The considerable advantages of our approach are speed and resolution: the SAPPHIRE plot is scalable to very large data sets and represents every single snapshot. This minimizes the risk of missing states because of overlap or prior coarse-graining of the data.

SUBMITTER: Blochliger N

PROVIDER: S-EPMC4151098 | biostudies-other | 2014

REPOSITORIES: biostudies-other

ACCESS DATA

Json Xml

Similar Datasets

Molecular dynamics simulations elucidate oligosaccharide recognition pathways by galectin-3 at atomic resolution.

Project description:The recognition of carbohydrates by lectins plays key roles in diverse cellular processes such as cellular adhesion, proliferation, and apoptosis, which makes it a therapeutic target of significance against cancers. One of the most functionally active lectins, galectin-3 is distinctively known for its specific binding affinity toward β-galactoside. However, despite the prevalence of high-resolution crystallographic structures, the mechanistic basis and more significantly, the dynamic process underlying carbohydrate recognition by galectin-3 are currently elusive. To this end, we employed extensive Molecular Dynamics simulations to unravel the complete binding event of human galectin-3 with its native natural ligand N-acetyllactosamine (LacNAc) at atomic precision. The simulation trajectory demonstrates that the oligosaccharide diffuses around the protein and eventually identifies and binds to the biologically designated binding site of galectin-3 in real time. The simulated bound pose correlates with the crystallographic pose with atomic-level accuracy and recapitulates the signature stabilizing galectin-3/oligosaccharide interactions. The recognition pathway also reveals a set of transient non-native ligand poses in its course to the receptor. Interestingly, kinetic analysis in combination with a residue-level picture revealed that the key to the efficacy of a more active structural variant of the LacNAc lay in the ligand's resilience against disassociation from galectin-3. By catching the ligand in the act of finding its target, our investigations elucidate the detailed recognition mechanism of the carbohydrate-binding domain of galectin-3 and underscore the importance of ligand-target binary complex residence time in understanding the structure-activity relationship of cognate ligands.

| S-EPMC8571523 | biostudies-literature

Applying molecular dynamics simulations to identify rarely sampled ligand-bound conformational states of undecaprenyl pyrophosphate synthase, an antibacterial target.

Project description:Undecaprenyl pyrophosphate synthase is a cis-prenyltransferase enzyme, which is required for cell wall biosynthesis in bacteria. Undecaprenyl pyrophosphate synthase is an attractive target for antimicrobial therapy. We performed long molecular dynamics simulations and docking studies on undecaprenyl pyrophosphate synthase to investigate its dynamic behavior and the influence of protein flexibility on the design of undecaprenyl pyrophosphate synthase inhibitors. We also describe the first X-ray crystallographic structure of Escherichia coli apo-undecaprenyl pyrophosphate synthase. The molecular dynamics simulations indicate that undecaprenyl pyrophosphate synthase is a highly flexible protein, with mobile binding pockets in the active site. By carrying out docking studies with experimentally validated undecaprenyl pyrophosphate synthase inhibitors using high- and low-populated conformational states extracted from the molecular dynamics simulations, we show that structurally dissimilar compounds can bind preferentially to different and rarely sampled conformational states. By performing structural analyses on the newly obtained apo-undecaprenyl pyrophosphate synthase and other crystal structures previously published, we show that the changes observed during the molecular dynamics simulation are very similar to those seen in the crystal structures obtained in the presence or absence of ligands. We believe that this is the first time that a rare 'expanded pocket' state, key to drug design and verified by crystallography, has been extracted from a molecular dynamics simulation.

| S-EPMC3095679 | biostudies-literature

High-resolution reversible folding of hyperstable RNA tetraloops using molecular dynamics simulations.

Project description:We report the de novo folding of three hyperstable RNA tetraloops to 1-3 Å rmsd from their experimentally determined structures using molecular dynamics simulations initialized in the unfolded state. RNA tetraloops with loop sequences UUCG, GCAA, or CUUG are hyperstable because of the formation of noncanonical loop-stabilizing interactions, and they are all faithfully reproduced to angstrom-level accuracy in replica exchange molecular dynamics simulations, including explicit solvent and ion molecules. This accuracy is accomplished using unique RNA parameters, in which biases that favor rigid, highly stacked conformations are corrected to accurately capture the inherent flexibility of ssRNA loops, accurate base stacking energetics, and purine syn-anti interconversions. In a departure from traditional quantum chemistrycentric approaches to force field optimization, our parameters are calibrated directly from thermodynamic and kinetic measurements of intra- and internucleotide structural transitions. The ability to recapitulate the signature noncanonical interactions of the three most abundant hyperstable stem loop motifs represents a significant milestone to the accurate prediction of RNA tertiary structure using unbiased all-atom molecular dynamics simulations.

| S-EPMC3801082 | biostudies-literature

Ion transfer mechanisms in Mrp-type antiporters from high resolution cryoEM and molecular dynamics simulations.

Project description:Multiple resistance and pH adaptation (Mrp) cation/proton antiporters are essential for growth of a variety of halophilic and alkaliphilic bacteria under stress conditions. Mrp-type antiporters are closely related to the membrane domain of respiratory complex I. We determined the structure of the Mrp antiporter from Bacillus pseudofirmus by electron cryo-microscopy at 2.2 Å resolution. The structure resolves more than 99% of the sidechains of the seven membrane subunits MrpA to MrpG plus 360 water molecules, including ~70 in putative ion translocation pathways. Molecular dynamics simulations based on the high-resolution structure revealed details of the antiport mechanism. We find that switching the position of a histidine residue between three hydrated pathways in the MrpA subunit is critical for proton transfer that drives gated trans-membrane sodium translocation. Several lines of evidence indicate that the same histidine-switch mechanism operates in respiratory complex I.

| S-EPMC9568556 | biostudies-literature

Capturing Membrane Phase Separation by Dual Resolution Molecular Dynamics Simulations.

Project description:Understanding the lateral organization in plasma membranes remains an open problem and is of great interest to many researchers. Model membranes consisting of coexisting domains are commonly used as simplified models of plasma membranes. The coarse-grained (CG) Martini force field has successfully captured spontaneous separation of ternary membranes into a liquid-disordered and a liquid-ordered domain. With all-atom (AA) models, however, phase separation is much harder to achieve due to the slow underlying dynamics. To remedy this problem, here, we apply the virtual site (VS) hybrid method on a ternary membrane composed of saturated lipids, unsaturated lipids, and cholesterol to investigate the phase separation. The VS scheme couples the two membrane leaflets at CG and AA resolution. We found that the rapid phase separation reached by the CG leaflet can accelerate and guide this process in the AA leaflet.

| S-EPMC8444333 | biostudies-literature

Flexible fitting of high-resolution x-ray structures into cryoelectron microscopy maps using biased molecular dynamics simulations.

Project description:A methodology for flexible fitting of all-atom high-resolution structures into low-resolution cryoelectron microscopy (cryo-EM) maps is presented. Flexibility of the modeled structure is simulated by classical molecular dynamics and an additional effective potential is introduced to enhance the fitting process. The additional potential is proportional to the correlation coefficient between the experimental cryo-EM map and a synthetic map generated for an all-atom structure being fitted to the map. The additional forces are calculated as a gradient of the correlation coefficient. During the molecular dynamics simulations under the additional forces, the molecule undergoes a conformational transition that maximizes the correlation coefficient, which results in a high-accuracy fit of all-atom structure into a cryo-EM map. Using five test proteins that exhibit structural rearrangement during their biological activity, we demonstrate performance of our method. We also test our method on the experimental cryo-EM of elongation factor G and show that the model obtained is comparable to previous studies. In addition, we show that overfitting can be avoided by assessing the quality of the fitted model in terms of correlation coefficient and secondary structure preservation.

| S-EPMC2599853 | biostudies-literature

Hybrid resolution molecular dynamics simulations of amyloid proteins interacting with membranes.

Project description:A broad range of human diseases, including Alzheimer's and Parkinson's diseases, arise from or have as key players intrinsically disordered proteins. The aggregation of these amyloid proteins into fibrillar aggregates are the key events of such diseases. Characterizing the conformation dynamics of the proteins involved is crucial for understanding the molecular mechanisms of aggregation, which in turn is important for drug development efforts against these diseases. Computational approaches have provided extensive detail about some steps of the aggregation process, however the biologically relevant elements responsible for the aggregation and or aggregation propagation have not been fully characterized. Here we describe a hybrid resolution molecular dynamics simulation method that can be employed to investigate the interaction of amyloid proteins with lipid membranes, shown to dramatically accelerate the aggregation propensity of amyloid proteins. The hybrid resolution method enables routine and accurate simulation of multi-protein and complex membrane systems, mimicking biologically relevant lipid membranes, on microsecond time scales. The hybrid resolution method was applied to computer modeling of the interactions of α -synuclein protein with a mixed lipid bilayer.

| S-EPMC8435541 | biostudies-literature

Intermediate states of the Kv1.2 voltage sensor from atomistic molecular dynamics simulations.

Project description:The response of a membrane-bound Kv1.2 ion channel to an applied transmembrane potential has been studied using molecular dynamics simulations. Channel deactivation is shown to involve three intermediate states of the voltage sensor domain (VSD), and concomitant movement of helix S4 charges 10-15 ? along the bilayer normal; the latter being enabled by zipper-like sequential pairing of S4 basic residues with neighboring VSD acidic residues and membrane-lipid head groups. During the observed sequential transitions S4 basic residues pass through the recently discovered charge transfer center with its conserved phenylalanine residue, F(233). Analysis indicates that the local electric field within the VSD is focused near the F(233) residue and that it remains essentially unaltered during the entire process. Overall, the present computations provide an atomistic description of VSD response to hyperpolarization, add support to the sliding helix model, and capture essential features inferred from a variety of recent experiments.

| S-EPMC3076833 | biostudies-other

Investigation of indazole unbinding pathways in CYP2E1 by molecular dynamics simulations.

Project description:Human microsomal cytochrome P450 2E1 (CYP2E1) can oxidize not only low molecular weight xenobiotic compounds such as ethanol, but also many endogenous fatty acids. The crystal structure of CYP2E1 in complex with indazole reveals that the active site is deeply buried into the protein center. Thus, the unbinding pathways and associated unbinding mechanisms remain elusive. In this study, random acceleration molecular dynamics simulations combined with steered molecular dynamics and potential of mean force calculations were performed to identify the possible unbinding pathways in CYP2E1. The results show that channel 2c and 2a are most likely the unbinding channels of CYP2E1. The former channel is located between helices G and I and the B-C loop, and the latter resides between the region formed by the F-G loop, the B-C loop and the ?1 sheet. Phe298 and Phe478 act as the gate keeper during indazole unbinding along channel 2c and 2a, respectively. Previous site-directed mutagenesis experiments also supported these findings.

| S-EPMC3307744 | biostudies-literature

Misfolding pathways of the prion protein probed by molecular dynamics simulations.

Project description:Although the cellular monomeric form of the benign prion protein is now well characterized, a model for the monomer of the misfolded conformation (PrP(Sc)) remains elusive. PrP(Sc) quickly aggregates into highly insoluble fibrils making experimental structural characterization very difficult. The tendency to aggregation of PrP(Sc) in aqueous solution implies that the monomer fold must be hydrophobic. Here, by using molecular dynamics simulations, we have studied the cellular mouse prion protein and its D178N pathogenic mutant immersed in a hydrophobic environment (solution of CCl4), to reveal conformational changes and/or local structural weaknesses of the prion protein fold in unfavorable structural and thermodynamic conditions. Simulations in water have been also performed. Although observing in general a rather limited conformation activity in the nanosecond timescale, we have detected a significant weakening of the antiparallel beta-sheet of the D178N mutant in CCl4 and to a less extent in water. No weakening is observed for the native prion protein. The increase of beta-structure in the monomer, recently claimed as evidence for misfolding to PrP(Sc), has been also observed in this study irrespective of the thermodynamic or structural conditions, showing that this behavior is very likely an intrinsic characteristic of the prion protein fold.

| S-EPMC1305135 | 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