Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEquilibrium sampling by re-weighting non-equilibrium simulation trajectories
90
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
With the traditional equilibrium molecular simulations, it is usually difficult to efficiently visit the whole conformational space in complex systems, which are separated into some metastable conformational regions by high free energy barriers. The applied non-equilibrium process in simulations could enhance the transitions among these conformational regions, and the associated non-equilibrium effects can be removed by employing the Jarzynski equality (JE), then the global equilibrium distribution can be reproduced. However, the original JE requires the initial distribution of the non-equilibrium process is equilibrium, which largely limits the application of the non-equilibrium method in equilibrium sampling. By extending the previous method, the reweighted ensemble dynamics (RED), which re-weights many equilibrium simulation trajectories from arbitrary initial distribution to reproduce the global equilibrium, to non-equilibrium simulations, we present a method, named as re-weighted non-equilibrium ensemble dynamics (RNED), to generalize the JE in the non-equilibrium trajectories started from an arbitrary initial distribution, thus provide an efficient method to reproduce the equilibrium distribution based on multiple independent (short) non-equilibrium trajectories. We have illustrated the validity of the RNED in a one-dimensional toy model and in a Lennard-Jones system to detect the liquid-solid phase coexistence.
rate research
Read More
Equilibrium sampling of biomolecules remains an unmet challenge after more than 30 years of atomistic simulation. Efforts to enhance sampling capability, which are reviewed here, range from the development of new algorithms to parallelization to novel uses of hardware. Special focus is placed on classifying algorithms -- most of which are underpinned by a few key ideas -- in order to understand their fundamental strengths and limitations. Although algorithms have proliferated, progress resulting from novel hardware use appears to be more clear-cut than from algorithms alone, partly due to the lack of widely used sampling measures.
Recent thermophoretic experiments on colloidal suspensions revived an old debate, namely whether the Soret effect is properly described by thermostatics, or necessarily requires non-equilibrium thermodynamics. Based on colloidal transport theory and the entropy production of the related viscous flow, our analysis leads to the conclusion that the equilibrium approach may work for small ions, yet fails for colloidal particles and polymers. Regarding binary molecular mixtures, our results shed some doubt on the validity of thermostatic approaches that derive the Soret coefficient from equilibrium potentials.
We study the conformational dynamics within homo-polymer globules by solvent-implicit Brownian dynamics simulations. A strong dependence of the internal chain dynamics on the Lennard-Jones cohesion strength {epsilon} and the globule size NG is observed. We find two distinct dynamical regimes: a liquid- like regime (for {epsilon} < {epsilon}s) with fast internal dynamics and a solid-like regime (for {epsilon} > {epsilon}s) with slow internal dynamics. The cohesion strength {epsilon}s of this freezing transition depends on NG. Equilibrium simulations, where we investigate the diffusional chain dynamics within the globule, are compared with non-equilibrium simulations, where we unfold the globule by pulling the chain ends with prescribed velocity (encompassing low enough velocities so that the linear-response, viscous regime is reached). From both simulation protocols we derive the internal viscosity within the globule. In the liquid-like regime the internal friction increases continuously with {epsilon} and scales extensive in NG. This suggests an internal friction scenario where the entire chain (or an extensive fraction thereof) takes part in conformational reorganization of the globular structure.
We provide a non-equilibrium thermodynamic description of the life-cycle of a droplet based, chemically feasible, system of protocells. By coupling the protocells metabolic kinetics with its thermodynamics, we demonstrate how the system can be driven out of equilibrium to ensure protocell growth and replication. This coupling allows us to derive the equations of evolution and to rigorously demonstrate how growth and replication life-cycle can be understood as a non-equilibrium thermodynamic cycle. The process does not appeal to genetic information or inheritance, and is based only on non-equilibrium physics considerations. Our non-equilibrium thermodynamic description of simple, yet realistic, processes of protocell growth and replication, represents an advance in our physical understanding of a central biological phenomenon both in connection to the origin of life and for modern biology.
Colonies of bacterial cells endowed with a pili-based self-propulsion machinery represent an ideal model system for studying how active adhesion forces affect structure and dynamics of many-particle systems. As a novel computational tool, we describe here a highly parallel molecular dynamics simulation package for modeling of textit{Neisseria gonorrhoeae} colonies. Simulations are employed to investigate growth of bacterial colonies and the dependence of the colony structure on cell-cell interactions. In agreement with experimental data, active pilus retraction is found to enhance local ordering. For mixed colonies consisting of different types of cell types, the simulations show a segregation of cell types depending on the pili-mediated interactions, as seen in experiments. Using a simulated experimental setup, we study the power-spectral density of colony-shape fluctuations and the associated fluctuation-response relation. The simulations predict a strong violation of the equilibrium fluctuation-response relation across the measurable frequency range. Lastly, we illustrate the essential role of active force generation for colony dynamics by showing that pilus-mediated activity drives the spreading of colonies on surfaces and the invasion of narrow channels.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
