Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Langevin dynamics simulation of polymer-assisted virus-like assembly.

ABSTRACT: Starting from a coarse grained representation of the building units of the minute virus of mice and a flexible polyelectrolyte molecule, we have explored the mechanism of assembly into icosahedral structures with the help of Langevin dynamics simulations and the parallel tempering technique. Regular icosahedra with appropriate symmetry form only in a narrow range of temperature and polymer length. Within this region of parameters where successful assembly would proceed, we have systematically investigated the growth kinetics. The assembly of icosahedra is found to follow the classical nucleation and growth mechanism in the absence of the polymer, with the three regimes of nucleation, linear growth, and slowing down in the later stage. The calculated average nucleation time obeys the laws expected from the classical nucleation theory. The linear growth rate is found to obey the laws of secondary nucleation as in the case of lamellar growth in polymer crystallization. The same mechanism is seen in the simulations of the assembly of icosahedra in the presence of the polymer as well. The polymer reduces the nucleation barrier significantly by enhancing the local concentration of subunits via adsorbing them on their backbone. The details of growth in the presence of the polymer are also found to be consistent with the classical nucleation theory, despite the smallness of the assembled structures.

SUBMITTER: Mahalik JP

PROVIDER: S-EPMC3331873 | biostudies-other | 2012 Apr

REPOSITORIES: biostudies-other

ACCESS DATA

Json Xml

Similar Datasets

Conformational Transitions of Amyloid-β: A Langevin and Generalized Langevin Dynamics Simulation Study.

Project description:The dynamics of conformational transitions of the disordered protein, amyloid-β, is studied via Langevin and generalized Langevin dynamics simulations. The transmission coefficient for the unfold-misfold transition of amyloid-β is calculated from multiple independent trajectories that originate at the transition state with different initial velocities and are directly correlated to Kramers and Grote-Hynes theories. For lower values of the frictional coefficient, a well-defined rate constant is obtained, whereas, for higher values, the transmission coefficient decays with time, indicating a breakdown of the Kramers and Grote-Hynes theories and the emergence of a dynamic disorder, which demonstrates the presence of multiple local minima in the misfolding potential energy surface. The calculated free energy profile describes a two-state transition of amyloid-β in the energy landscape. The transition path time distribution computed from these simulations is compared with the related experimental and theoretical results for the unfold-misfold transition of amyloid-β. The high free energy barrier for this transition confirms the misfolding of amyloid-β. These findings offer an insight into the dynamics of the unfold-misfold transition of this protein.

| S-EPMC8173568 | biostudies-literature

Toward canonical ensemble distribution from self-guided Langevin dynamics simulation.

Project description:This work derives a quantitative description of the conformational distribution in self-guided Langevin dynamics (SGLD) simulations. SGLD simulations employ guiding forces calculated from local average momentums to enhance low-frequency motion. This enhancement in low-frequency motion dramatically accelerates conformational search efficiency, but also induces certain perturbations in conformational distribution. Through the local averaging, we separate properties of molecular systems into low-frequency and high-frequency portions. The guiding force effect on the conformational distribution is quantitatively described using these low-frequency and high-frequency properties. This quantitative relation provides a way to convert between a canonical ensemble and a self-guided ensemble. Using example systems, we demonstrated how to utilize the relation to obtain canonical ensemble properties and conformational distributions from SGLD simulations. This development makes SGLD not only an efficient approach for conformational searching, but also an accurate means for conformational sampling.

| S-EPMC3087419 | biostudies-other

Variance Reduction in Stochastic Gradient Langevin Dynamics.

Project description:Stochastic gradient-based Monte Carlo methods such as stochastic gradient Langevin dynamics are useful tools for posterior inference on large scale datasets in many machine learning applications. These methods scale to large datasets by using noisy gradients calculated using a mini-batch or subset of the dataset. However, the high variance inherent in these noisy gradients degrades performance and leads to slower mixing. In this paper, we present techniques for reducing variance in stochastic gradient Langevin dynamics, yielding novel stochastic Monte Carlo methods that improve performance by reducing the variance in the stochastic gradient. We show that our proposed method has better theoretical guarantees on convergence rate than stochastic Langevin dynamics. This is complemented by impressive empirical results obtained on a variety of real world datasets, and on four different machine learning tasks (regression, classification, independent component analysis and mixture modeling). These theoretical and empirical contributions combine to make a compelling case for using variance reduction in stochastic Monte Carlo methods.

| S-EPMC5508544 | biostudies-literature

Inferring temporal dynamics from cross-sectional data using Langevin dynamics.

Project description:Cross-sectional studies are widely prevalent since they are more feasible to conduct compared with longitudinal studies. However, cross-sectional data lack the temporal information required to study the evolution of the underlying dynamics. This temporal information is essential to develop predictive computational models, which is the first step towards causal modelling. We propose a method for inferring computational models from cross-sectional data using Langevin dynamics. This method can be applied to any system where the data-points are influenced by equal forces and are in (local) equilibrium. The inferred model will be valid for the time span during which this set of forces remains unchanged. The result is a set of stochastic differential equations that capture the temporal dynamics, by assuming that groups of data-points are subject to the same free energy landscape and amount of noise. This is a 'baseline' method that initiates the development of computational models and can be iteratively enhanced through the inclusion of domain expert knowledge as demonstrated in our results. Our method shows significant predictive power when compared against two population-based longitudinal datasets. The proposed method can facilitate the use of cross-sectional datasets to obtain an initial estimate of the underlying dynamics of the respective systems.

| S-EPMC8580443 | biostudies-literature

Normal mode partitioning of Langevin dynamics for biomolecules.

Project description:We propose a novel normal mode multiple time stepping Langevin dynamics integrator called NML. The aim is to approximate the kinetics or thermodynamics of a biomolecule by a reduced model based on a normal mode decomposition of the dynamical space. Our basis set uses the eigenvectors of a mass reweighted Hessian matrix calculated with a biomolecular force field. This particular choice has the advantage of an ordering according to the eigenvalues, which have a physical meaning of being the square of the mode frequency. Low frequency eigenvalues correspond to more collective motions, whereas the highest frequency eigenvalues are the limiting factor for the stability of the integrator. In NML, the higher frequency modes are overdamped and relaxed near their energy minimum while respecting the subspace of low frequency dynamical modes. Our numerical results confirm that both sampling and rates are conserved for an implicitly solvated alanine dipeptide model, with only 30% of the modes propagated, when compared to the full model. For implicitly solvated systems, NML gives a twofold improvement in efficiency over plain Langevin dynamics for sampling a small 22 atom (alanine dipeptide) model and in excess of an order of magnitude for sampling an 882 atom (bovine pancreatic trypsin inhibitor) model, with good scaling with system size subject to the number of modes propagated. NML has been implemented in the open source software PROTOMOL.

| S-EPMC2673189 | biostudies-literature

Multiscale dynamics of macromolecules using normal mode Langevin.

Project description:Proteins and other macromolecules have coupled dynamics over multiple time scales (from femtosecond to millisecond and beyond) that make resolving molecular dynamics challenging. We present an approach based on periodically decomposing the dynamics of a macromolecule into slow and fast modes based on a scalable coarse-grained normal mode analysis. A Langevin equation is used to propagate the slowest degrees of freedom while minimizing the nearly instantaneous degrees of freedom. We present numerical results showing that time steps of up to 1000 fs can be used, with real speedups of up to 200 times over plain molecular dynamics. We present results of successfully folding the Fip35 mutant of WW domain.

| S-EPMC4308582 | biostudies-literature

Vulnerability in Popular Molecular Dynamics Packages Concerning Langevin and Andersen Dynamics.

Project description:We report a serious problem associated with a number of current implementations of Andersen and Langevin dynamics algorithms. When long simulations are run in many segments, it is sometimes possible to have a repeating sequence of pseudorandom numbers enter the calcuation. We show that, if the sequence repeats rapidly, the resulting artifacts can quickly denature biomolecules and are then easily detectable. However, if the sequence repeats less frequently, the artifacts become subtle and easily overlooked. We derive a formula for the underlying cause of artifacts in the case of the Langevin thermostat, and find it vanishes slowly as the inverse square root of the number of time steps per simulation segment. Numerous examples of simulation artifacts are presented, including dissociation of a tetrameric protein after 110 ns of dynamics, reductions in atomic fluctuations for a small protein in implicit solvent, altered thermodynamic properties of a box of water molecules, and changes in the transition free energies between dihedral angle conformations. Finally, in the case of strong thermocoupling, we link the observed artifacts to previous work in nonlinear dynamics and show that it is possible to drive a 20-residue, implicitly solvated protein into periodic trajectories if the thermostat is not used properly. Our findings should help other investigators re-evaluate simulations that may have been corrupted and obtain more accurate results.

| S-EPMC2632580 | biostudies-literature

Efficient and flexible implementation of Langevin simulation for gene burst production.

Project description:Gene expression involves bursts of production of both mRNA and protein, and the fluctuations in their number are increased due to such bursts. The Langevin equation is an efficient and versatile means to simulate such number fluctuation. However, how to include these mRNA and protein bursts in the Langevin equation is not intuitively clear. In this work, we estimated the variance in burst production from a general gene expression model and introduced such variation in the Langevin equation. Our approach offers different Langevin expressions for either or both transcriptional and translational bursts considered and saves computer time by including many production events at once in a short burst time. The errors can be controlled to be rather precise (<2%) for the mean and <10% for the standard deviation of the steady-state distribution. Our scheme allows for high-quality stochastic simulations with the Langevin equation for gene expression, which is useful in analysis of biological networks.

| S-EPMC5715166 | biostudies-literature

Molecular dynamics simulation in virus research.

Project description:Virus replication in the host proceeds by chains of interactions between viral and host proteins. The interactions are deeply influenced by host immune molecules and anti-viral compounds, as well as by mutations in viral proteins. To understand how these interactions proceed mechanically and how they are influenced by mutations, one needs to know the structures and dynamics of the proteins. Molecular dynamics (MD) simulation is a powerful computational method for delineating motions of proteins at an atomic-scale via theoretical and empirical principles in physical chemistry. Recent advances in the hardware and software for biomolecular simulation have rapidly improved the precision and performance of this technique. Consequently, MD simulation is quickly extending the range of applications in biology, helping to reveal unique features of protein structures that would be hard to obtain by experimental methods alone. In this review, we summarize the recent advances in MD simulations in the study of virus-host interactions and evolution, and present future perspectives on this technique.

| S-EPMC3400276 | biostudies-literature

Template-Assisted Self-Assembly of Conductive Polymer Electrodes for Ionic Electroactive Polymers.

Project description:Ionic electroactive polymers (ionic EAPs) can greatly aid in biomedical applications where micro-sized actuators are required for delicate procedures. Since these types of actuators generally require platinum or gold metallic electrodes, they tend to be expensive and susceptible to delamination. Previous research has solved this problem by replacing the metallic electrodes with conductive polymers (CP) and forming an interpenetrating polymer network (IPN) between the conductive polymer (CP) and the solid polymer electrolyte (SPE). Since these actuators contain toxic ionic liquids, they are unsuitable for biological applications. In this study, we present a novel and facile method of fabricating a biocompatible and ionic liquid-free actuator that uses semi-IPN to hold the CP and Nafion-based SPE layers together. Surface activated fabrication treatment (SAFT) is applied to the precursor-Nafion membrane in order to convert the sulfonyl fluoride groups on the surface to sulfonate. Through template-assisted self-assembly, the CP electrodes from either polyaniline (PANI) or poly(3,4-ethylenedioxythiophene) (PEDOT) interlock with the surface treated precursor-Nafion membrane so that no delamination can occur. The electrodes growth pattern, interfacial layer's thickness, and shape can be controlled by adjusting the SAFT concentration and duration.

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