Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

AWE-WQ: fast-forwarding molecular dynamics using the accelerated weighted ensemble.

ABSTRACT: A limitation of traditional molecular dynamics (MD) is that reaction rates are difficult to compute. This is due to the rarity of observing transitions between metastable states since high energy barriers trap the system in these states. Recently the weighted ensemble (WE) family of methods have emerged which can flexibly and efficiently sample conformational space without being trapped and allow calculation of unbiased rates. However, while WE can sample correctly and efficiently, a scalable implementation applicable to interesting biomolecular systems is not available. We provide here a GPLv2 implementation called AWE-WQ of a WE algorithm using the master/worker distributed computing WorkQueue (WQ) framework. AWE-WQ is scalable to thousands of nodes and supports dynamic allocation of computer resources, heterogeneous resource usage (such as central processing units (CPU) and graphical processing units (GPUs) concurrently), seamless heterogeneous cluster usage (i.e., campus grids and cloud providers), and support for arbitrary MD codes such as GROMACS, while ensuring that all statistics are unbiased. We applied AWE-WQ to a 34 residue protein which simulated 1.5 ms over 8 months with peak aggregate performance of 1000 ns/h. Comparison was done with a 200 μs simulation collected on a GPU over a similar timespan. The folding and unfolded rates were of comparable accuracy.

SUBMITTER: Abdul-Wahid B

PROVIDER: S-EPMC4210180 | biostudies-other | 2014 Oct

REPOSITORIES: biostudies-other

ACCESS DATA

Json Xml

Similar Datasets

Gaussian-Accelerated Molecular Dynamics with the Weighted Ensemble Method: A Hybrid Method Improves Thermodynamic and Kinetic Sampling.

Project description:Gaussian-accelerated molecular dynamics (GaMD) is a well-established enhanced sampling method for molecular dynamics simulations that effectively samples the potential energy landscape of the system by adding a boost potential, which smoothens the surface and lowers the energy barriers between states. GaMD is unable to give time-dependent properties such as kinetics directly. On the other hand, the weighted ensemble (WE) method can efficiently sample transitions between states with its many weighted trajectories, which directly yield rates and pathways. However, convergence to equilibrium conditions remains a challenge for the WE method. Hence, we have developed a hybrid method that combines the two methods, wherein GaMD is first used to sample the potential energy landscape of the system and WE is subsequently used to further sample the potential energy landscape and kinetic properties of interest. We show that the hybrid method can sample both thermodynamic and kinetic properties more accurately and quickly compared to using either method alone.

| S-EPMC8983023 | biostudies-literature

Bimolecular reaction simulation using Weighted Ensemble Brownian dynamics and the University of Houston Brownian Dynamics program.

Project description:We discuss here the implementation of the Weighted Ensemble Brownian (WEB) dynamics algorithm of Huber and Kim in the University of Houston Brownian Dynamics (UHBD) suite of programs and its application to bimolecular association problems. WEB dynamics is a biased Brownian dynamics (BD) algorithm that is more efficient than the standard Northrup-Allison-McCammon (NAM) method in cases where reaction events are infrequent because of intervening free energy barriers. Test cases reported here include the Smoluchowski rate for association of spheres, the association of the enzyme copper-zinc superoxide dismutase with superoxide anion, and the binding of the superpotent sweetener N-(p-cyanophenyl)-N'-(diphenylmethyl)-guanidinium acetic acid to a monoclonal antibody fragment, NC6.8. Our results show that the WEB dynamics algorithm is a superior simulation method for enzyme-substrate reaction encounters with large free energy barriers.

| S-EPMC1300969 | biostudies-other

Exploring CCRL2 chemerin binding using accelerated molecular dynamics.

Project description:Chemokine (C-C motif) receptor-like 2 (CCRL2), is a seven transmembrane receptor closely related to the chemokine receptors CCR1, CCR2, CCR3, and CCR5. Nevertheless, CCRL2 is unable to activate conventional G-protein dependent signaling and to induce cell directional migration. The only commonly accepted CCRL2 ligand is the nonchemokine chemotactic protein chemerin (RARRES2). The chemerin binding to CCLR2 does induce leukocyte chemotaxis, yet, genetic targeting of CCRL2 was shown to modulate the inflammatory response in different experimental models. This mechanism was shown to be crucial for lung dendritic cell migration, neutrophil recruitment, and Natural Killer cell-dependent immune surveillance in lung cancer. To gain more insight in the interactions involved in the CCRL2-chemerin, the binding complexes were generated by protein-protein docking, then submitted to accelerated molecular dynamics. The obtained trajectories were inspected by principal component analyses followed by kernel density estimation to identify the ligand-receptor regions most frequently involved in the binding. To conclude, the reported analyses led to the identification of the putative hot-spot residues involved in CCRL2-chemerin binding.

| S-EPMC9543397 | biostudies-literature

Enhancing protein adsorption simulations by using accelerated molecular dynamics.

Project description:The atomistic modeling of protein adsorption on surfaces is hampered by the different time scales of the simulation ([Formula: see text][Formula: see text]s) and experiment (up to hours), and the accordingly different 'final' adsorption conformations. We provide evidence that the method of accelerated molecular dynamics is an efficient tool to obtain equilibrated adsorption states. As a model system we study the adsorption of the protein BMP-2 on graphite in an explicit salt water environment. We demonstrate that due to the considerably improved sampling of conformational space, accelerated molecular dynamics allows to observe the complete unfolding and spreading of the protein on the hydrophobic graphite surface. This result is in agreement with the general finding of protein denaturation upon contact with hydrophobic surfaces.

| S-EPMC3670854 | biostudies-literature

Binding Analysis Using Accelerated Molecular Dynamics Simulations and Future Perspectives.

Project description:Biomolecular recognition such as binding of small molecules, nucleic acids, peptides and proteins to their target receptors plays key roles in cellular function and has been targeted for therapeutic drug design. Molecular dynamics (MD) is a computational approach to analyze these binding processes at an atomistic level, which provides valuable understandings of the mechanisms of biomolecular recognition. However, the rather slow biomolecular binding events often present challenges for conventional MD (cMD), due to limited simulation timescales (typically over hundreds of nanoseconds to tens of microseconds). In this regard, enhanced sampling methods, particularly accelerated MD (aMD), have proven useful to bridge the gap and enable all-atom simulations of biomolecular binding events. Here, we will review the recent method developments of Gaussian aMD (GaMD), ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD), which have greatly expanded our capabilities to simulate biomolecular binding processes. Spontaneous binding of various biomolecules to their receptors has been successfully simulated by GaMD. Microsecond LiGaMD and Pep-GaMD simulations have captured repetitive binding and dissociation of small-molecule ligands and highly flexible peptides, and thus enabled ligand/peptide binding thermodynamics and kinetics calculations. We will also present relevant application studies in simulations of important drug targets and future perspectives for rational computer-aided drug design.

| S-EPMC8747661 | biostudies-literature

Accelerated Estimation of Long-Timescale Kinetics from Weighted Ensemble Simulation via Non-Markovian "Microbin" Analysis.

Project description:The weighted ensemble (WE) simulation strategy provides unbiased sampling of nonequilibrium processes, such as molecular folding or binding, but the extraction of rate constants relies on characterizing steady-state behavior. Unfortunately, WE simulations of sufficiently complex systems will not relax to steady state on observed simulation times. Here, we show that a postsimulation clustering of molecular configurations into "microbins" using methods developed in the Markov State Model (MSM) community can yield unbiased kinetics from WE data before steady-state convergence of the WE simulation itself. Because WE trajectories are directional and not equilibrium distributed, the history-augmented MSM (haMSM) formulation can be used, which yields the mean first-passage time (MFPT) without bias for arbitrarily small lag times. Accurate kinetics can be obtained while bypassing the often prohibitive convergence requirements of the nonequilibrium weighted ensemble. We validate the method in a simple diffusive process on a two-dimensional (2D) random energy landscape and then analyze atomistic protein folding simulations using WE molecular dynamics. We report significant progress toward the unbiased estimation of protein folding times and pathways, though key challenges remain.

| S-EPMC8045600 | biostudies-literature

Large-scale conformational sampling of proteins using temperature-accelerated molecular dynamics.

Project description:We show how to apply the method of temperature-accelerated molecular dynamics (TAMD) in collective variables [Maragliano L, Vanden-Eijnden E (2006) Chem Phys Lett 426:168-175] to sample the conformational space of multidomain proteins in all-atom, explicitly solvated molecular dynamics simulations. The method allows the system to hyperthermally explore the free-energy surface in a set of collective variables computed at the physical temperature. As collective variables, we pick Cartesian coordinates of centers of contiguous subdomains. The method is applied to the GroEL subunit, a 55-kDa, three-domain protein, and HIV-1 gp120. For GroEL, the method induces in about 40 ns conformational changes that recapitulate the t --> r('') transition and are not observed in unaccelerated molecular dynamics: The apical domain is displaced by 30 A, with a twist of 90 degrees relative to the equatorial domain, and the root-mean-squared deviation relative to the r('') conformer is reduced from 13 to 5 A, representing fairly high predictive capability. For gp120, the method predicts both counterrotation of inner and outer domains and disruption of the so-called bridging sheet. In particular, TAMD on gp120 initially in the CD4-bound conformation visits conformations that deviate by 3.6 A from the gp120 conformer in complex with antibody F105, again reflecting good predictive capability. TAMD generates plausible all-atom models of the so-far structurally uncharacterized unliganded conformation of HIV-1 gp120, which may prove useful in the development of inhibitors and immunogens. The fictitious temperature employed also gives a rough estimate of 10 kcal/mol for the free-energy barrier between conformers in both cases.

| S-EPMC2841907 | biostudies-other

Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics.

Project description:An extension of the constant pH method originally implemented by Mongan et al. (J. Comput. Chem.2004, 25, 2038-2048) is proposed in this study. This adapted version of the method couples the constant pH methodology with the enhanced sampling technique of accelerated molecular dynamics, in an attempt to overcome the sampling issues encountered with current standard constant pH molecular dynamics methods. Although good results were reported by Mongan et al. on application of the standard method to the hen egg-white lysozyme (HEWL) system, residues which possess strong interactions with neighboring groups tend to converge slowly, resulting in the reported inconsistencies for predicted pK(a) values, as highlighted by the authors. The application of the coupled method described in this study to the HEWL system displays improvements over the standard version of the method, with the improved sampling leading to faster convergence and producing pK(a) values in closer agreement to those obtained experimentally for the more slowly converging residues.

| S-EPMC2817915 | biostudies-literature

Analysis of mutations of defensin protein using accelerated molecular dynamics simulations.

Project description:Plant defensins possess diverse biological functions that include antifungal and antibacterial activities and α-amylase and trypsin inhibitory properties. Two mutations, G9R and V39R, were confirmed to increase the antifungal activity of Raphanus sativus antifungal protein 2 (RsAFP2). Accelerated Molecular Dynamics (aMD) were carried out to examine the conformational changes present in these RsAFP2 mutants, and its two closest homologs compared to the wild-type protein. Specifically, the root mean square fluctuation values for the eight cysteine amino acids involved in the four disulfide bonds were low in the V39R mutant compared to the wild-type. Additionally, analysis of the free energy change revealed that G9R and V39R mutations exert a neutral and stabilizing effect on RsAFP2 conformation, and this is supported by the observed lower total energy of mutants compared to the wild-type, suggesting that enhanced stability of the mutants. However, MD simulations to a longer time scale would aid in capturing more conformational state of the wild-type and mutants defensin protein. Furthermore, the aMD simulations on fungal mimic membranes with RsAFP2 and its mutants and homologs showed that the mutant proteins caused higher deformation and water diffusion than the native RsAFP2, especially the V39R mutant. The mutant variants seem to interact by specifically targeting the POPC and POPI lipids amongst others. This work highlights the stabilizing effect of mutations at the 9th and 39th positions of RsAFP2 and their increased membrane deformation activity.

| S-EPMC7703945 | biostudies-literature

Exploring ligand binding pathways on proteins using hypersound-accelerated molecular dynamics.

Project description:Capturing the dynamic processes of biomolecular systems in atomistic detail remains difficult despite recent experimental advances. Although molecular dynamics (MD) techniques enable atomic-level observations, simulations of "slow" biomolecular processes (with timescales longer than submilliseconds) are challenging because of current computer speed limitations. Therefore, we developed a method to accelerate MD simulations by high-frequency ultrasound perturbation. The binding events between the protein CDK2 and its small-molecule inhibitors were nearly undetectable in 100-ns conventional MD, but the method successfully accelerated their slow binding rates by up to 10-20 times. Hypersound-accelerated MD simulations revealed a variety of microscopic kinetic features of the inhibitors on the protein surface, such as the existence of different binding pathways to the active site. Moreover, the simulations allowed the estimation of the corresponding kinetic parameters and exploring other druggable pockets. This method can thus provide deeper insight into the microscopic interactions controlling biomolecular processes.

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