scholarly journals Multi-level free energy simulation with a staged transformation approach

2020 ◽  
Vol 153 (4) ◽  
pp. 044115 ◽  
Author(s):  
Shingo Ito ◽  
Qiang Cui
Keyword(s):  
Free Energy ◽  
Energy Simulation ◽  
Transformation Approach ◽  
Multi Level

Determination of Base Flipping Free Energy Landscapes from Nonequilibrium Stratification

2019 ◽  
Author(s):  
Xiaohui Wang ◽  
Zhaoxi Sun
Keyword(s):  
Free Energy ◽  
Energy Landscape ◽  
Sampling Technique ◽  
Umbrella Sampling ◽  
Energy Simulation ◽  
Nonlinear Dependence ◽  
Free Energy Landscape ◽  
Convergence Criterion ◽  
Base Flipping ◽  
Convergence Behavior

<p>Correct calculation of the variation of free energy upon base flipping is crucial in understanding the dynamics of DNA systems. The free energy landscape along the flipping pathway gives the thermodynamic stability and the flexibility of base-paired states. Although numerous free energy simulations are performed in the base flipping cases, no theoretically rigorous nonequilibrium techniques are devised and employed to investigate the thermodynamics of base flipping. In the current work, we report a general nonequilibrium stratification scheme for efficient calculation of the free energy landscape of base flipping in DNA duplex. We carefully monitor the convergence behavior of the equilibrium sampling based free energy simulation and the nonequilibrium stratification and determine the empirical length of time blocks required for converged sampling. Comparison between the performances of equilibrium umbrella sampling and nonequilibrium stratification is given. The results show that nonequilibrium free energy simulation is able to give similar accuracy and efficiency compared with the equilibrium enhanced sampling technique in the base flipping cases. We further test a convergence criterion we previously proposed and it comes out that the convergence behavior determined by this criterion agrees with those given by the time-invariant behavior of PMF and the nonlinear dependence of standard deviation on the sample size. </p>

Quantum Chemical and Free Energy Simulation Analysis of Retinal Conformational Energetics

1997 ◽  
Vol 37 (6) ◽  
pp. 1018-1024 ◽  
Author(s):  
Jérôme Baudry ◽  
Serge Crouzy ◽  
Benoît Roux ◽  
Jeremy C. Smith
Keyword(s):  
Free Energy ◽  
Quantum Chemical ◽  
Simulation Analysis ◽  
Energy Simulation

scholarly journals Determination of Base Flipping Free Energy Landscapes from Nonequilibrium Stratification

2018 ◽  
Author(s):  
Xiaohui Wang ◽  
Zhaoxi Sun
Keyword(s):  
Free Energy ◽  
Energy Landscape ◽  
Sampling Technique ◽  
Umbrella Sampling ◽  
Energy Simulation ◽  
Nonlinear Dependence ◽  
Free Energy Landscape ◽  
Convergence Criterion ◽  
Base Flipping ◽  
Convergence Behavior

<p>Correct calculation of the variation of free energy upon base flipping is crucial in understanding the dynamics of DNA systems. The free energy landscape along the flipping pathway gives the thermodynamic stability and the flexibility of base-paired states. Although numerous free energy simulations are performed in the base flipping cases, no theoretically rigorous nonequilibrium techniques are devised and employed to investigate the thermodynamics of base flipping. In the current work, we report a general nonequilibrium stratification scheme for efficient calculation of the free energy landscape of base flipping in DNA duplex. We carefully monitor the convergence behavior of the equilibrium sampling based free energy simulation and the nonequilibrium stratification and determine the empirical length of time blocks required for converged sampling. Comparison between the performances of equilibrium umbrella sampling and nonequilibrium stratification is given. The results show that nonequilibrium free energy simulation is able to give similar accuracy and efficiency compared with the equilibrium enhanced sampling technique in the base flipping cases. We further test a convergence criterion we previously proposed and it comes out that the convergence behavior determined by this criterion agrees with those given by the time-invariant behavior of PMF and the nonlinear dependence of standard deviation on the sample size. </p>

scholarly journals Determination of Base Flipping Free Energy Landscapes from Nonequilibrium Stratification

2019 ◽  
Author(s):  
Xiaohui Wang ◽  
Zhaoxi Sun
Keyword(s):  
Free Energy ◽  
Energy Landscape ◽  
Sampling Technique ◽  
Umbrella Sampling ◽  
Energy Simulation ◽  
Nonlinear Dependence ◽  
Free Energy Landscape ◽  
Convergence Criterion ◽  
Base Flipping ◽  
Convergence Behavior

<p>Correct calculation of the variation of free energy upon base flipping is crucial in understanding the dynamics of DNA systems. The free energy landscape along the flipping pathway gives the thermodynamic stability and the flexibility of base-paired states. Although numerous free energy simulations are performed in the base flipping cases, no theoretically rigorous nonequilibrium techniques are devised and employed to investigate the thermodynamics of base flipping. In the current work, we report a general nonequilibrium stratification scheme for efficient calculation of the free energy landscape of base flipping in DNA duplex. We carefully monitor the convergence behavior of the equilibrium sampling based free energy simulation and the nonequilibrium stratification and determine the empirical length of time blocks required for converged sampling. Comparison between the performances of equilibrium umbrella sampling and nonequilibrium stratification is given. The results show that nonequilibrium free energy simulation is able to give similar accuracy and efficiency compared with the equilibrium enhanced sampling technique in the base flipping cases. We further test a convergence criterion we previously proposed and it comes out that the convergence behavior determined by this criterion agrees with those given by the time-invariant behavior of PMF and the nonlinear dependence of standard deviation on the sample size. </p>

scholarly journals Seeding the Multi-dimensional Nonequilibrium Pulling for Hamiltonian Variation: Indirect QM/MM Free Energy Simulations

2021 ◽  
Author(s):  
Zhaoxi Sun ◽  
Qiaole He
Keyword(s):  
Free Energy ◽  
Quantum Mechanics ◽  
Computational Cost ◽  
Selection Criterion ◽  
Energy Difference ◽  
Energy Simulation ◽  
Computationally Efficient ◽  
Free Energy Difference ◽  
Energy Simulations ◽  
Semi Empirical

<p>The combination of free energy simulations in the alchemical and configurational spaces provides a feasible route to access the thermodynamic profiles under a computationally demanding target Hamiltonian. Normally, due to the significant differences between the computational cost of ab initio quantum mechanics (QM) calculations and those of semi-empirical quantum mechanics (SQM) and molecular mechanics (MM), this indirect method is often applied to obtain the QM thermodynamics by combining the SQM or MM thermodynamics and the SQM-to-QM or MM-to-QM corrections. In our previous works, a multi-dimensional nonequilibrium pulling framework for Hamiltonian variations has been introduced based on bidirectional pulling and bidirectional reweighting. The method performs nonequilibrium free energy simulations in the configurational space to obtain the thermodynamic profile along the conformational change pathway under a selected computationally efficient Hamiltonian, and uses the nonequilibrium alchemical method to correct or perturb the thermodynamic profile to that under the target Hamiltonian. The BAR-based method is designed to achieve the best generality and transferability and thus leads to modest (~20 fold) speedup. In this work, we explore the possibility of further accelerating the nonequilibrium free energy simulation by employing unidirectional pulling and using the selection criterion to obtain the initial configurations used to initiate nonequilibrium trajectories following the idea of adaptive steered molecular dynamics (ASMD). A single initial condition is used to seed the whole multi-dimensional nonequilibrium free energy simulation and the sampling is performed fully in the nonequilibrium ensemble. The ASMD scheme estimates the free energy difference with the unidirectional exponential average (EXP), but it does not follow exactly the requirements of the EXP estimator. Another consequence of the seeding simulation is the inherently sequential or serial pulling due to the inter-segment dependency, which triggers some problems in the parallelizability of the simulation. Therefore, some tests are required to grasp some insights and guidelines for using this selection-criterion-based ASMD scheme. The ASMD method is tested thoroughly on a dihedral flipping model system and encouraging numerical results are obtained. The selection-criterion-based multi-dimensional ASMD framework follows the same perturbation framework of the BAR-based method, and thus could be used in various Hamiltonian-variation cases.</p>

scholarly journals Seeding the Multi-dimensional Nonequilibrium Pulling for Hamiltonian Variation: Indirect Nonequilibrium Free Energy Simulations at QM levels

2021 ◽  
Author(s):  
Zhaoxi Sun ◽  
Qiaole He
Keyword(s):  
Free Energy ◽  
Quantum Mechanics ◽  
Computational Cost ◽  
Selection Criterion ◽  
Energy Difference ◽  
Energy Simulation ◽  
Computationally Efficient ◽  
Free Energy Difference ◽  
Energy Simulations ◽  
Semi Empirical

The combination of free energy simulations in the alchemical and configurational spaces provides a feasible route to access the thermodynamic profiles under a computationally demanding target Hamiltonian. Normally, due to the significant differences between the computational cost of ab initio quantum mechanics (QM) calculations and those of semi-empirical quantum mechanics (SQM) and molecular mechanics (MM), this indirect method could be used to obtain the QM thermodynamics by combining the SQM or MM results and the SQM-to-QM or MM-to-QM corrections. In our previous works, a multi-dimensional nonequilibrium pulling framework for Hamiltonian variations has been introduced based on bidirectional pulling and bidirectional reweighting. The method performs nonequilibrium free energy simulations in the configurational space to obtain the thermodynamic profile along the conformational change pathway under a selected computationally efficient Hamiltonian, and uses the nonequilibrium alchemical method to correct or perturb the thermodynamic profile to that under the target Hamiltonian. The BAR-based method is designed to achieve the best generality and transferability and thus leads to modest (~20 folds) speedup. In this work, we explore the possibility of further accelerating the nonequilibrium free energy simulation by employing unidirectional pulling and using the selection criterion to obtain the initial configurations used to initiate nonequilibrium trajectories following the idea of adaptive steered molecular dynamics (ASMD). A single initial condition is used to seed the whole multi-dimensional nonequilibrium free energy simulation and the sampling is performed fully in the nonequilibrium ensemble. Introducing very short ps-length equilibrium sampling to grab more initial seeds could also be helpful. The ASMD scheme estimates the free energy difference with the unidirectional exponential average (EXP), but it does not follow exactly the requirements of the EXP estimator. Another deficiency of the seeding simulation is the inherently sequential or serial pulling due to the inter-segment dependency, which triggers some problems in the parallelizability of the simulation. Numerical tests are performed to grasp some insights and guidelines for using this selection-criterion-based ASMD scheme. The presented selection-criterion-based multi-dimensional ASMD scheme follows the same perturbation network of the BAR-based method, and thus could be used in various Hamiltonian-variation cases.

Sign in / Sign up

Export Citation Format

Share Document