1
|
Fass J, York F, Wittmann M, Kaus J, Zhao Y. Local Resampling Trick for Focused Molecular Dynamics. J Chem Theory Comput 2023; 19:6139-6150. [PMID: 37706456 DOI: 10.1021/acs.jctc.3c00539] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/15/2023]
Abstract
We describe a method that focuses sampling effort on a user-defined selection of a large system, which can lead to substantial decreases in computational effort by speeding up the calculation of nonbonded interactions. A naive approach can lead to incorrect sampling if the selection depends on the configuration in a way that is not accounted for. We avoid this pitfall by introducing appropriate auxiliary variables. This results in an implementation that is closely related to "configurational freezing" and "elastic barrier dynamical freezing." We implement the method and validate that it can be used to supplement conventional molecular dynamics in free energy calculations (absolute hydration and relative binding).
Collapse
Affiliation(s)
- Joshua Fass
- Computation, Relay Therapeutics, Cambridge, Massachusetts 02139, United States
| | - Forrest York
- Computation, Relay Therapeutics, Cambridge, Massachusetts 02139, United States
| | - Matthew Wittmann
- Computation, Relay Therapeutics, Cambridge, Massachusetts 02139, United States
| | - Joseph Kaus
- Computation, Relay Therapeutics, Cambridge, Massachusetts 02139, United States
| | - Yutong Zhao
- Computation, Relay Therapeutics, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
2
|
Boubeta FM, Contestín García RM, Lorenzo EN, Boechi L, Estrin D, Sued M, Arrar M. Lessons learned about steered molecular dynamics simulations and free energy calculations. Chem Biol Drug Des 2019; 93:1129-1138. [PMID: 30793836 DOI: 10.1111/cbdd.13485] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2018] [Revised: 12/18/2018] [Accepted: 12/22/2018] [Indexed: 01/30/2023]
Abstract
The calculation of free energy profiles is central in understanding differential enzymatic activity, for instance, involving chemical reactions that require QM-MM tools, ligand migration, and conformational rearrangements that can be modeled using classical potentials. The use of steered molecular dynamics (sMD) together with the Jarzynski equality is a popular approach in calculating free energy profiles. Here, we first briefly review the application of the Jarzynski equality to sMD simulations, then revisit the so-called stiff-spring approximation and the consequent expectation of Gaussian work distributions and, finally, reiterate the practical utility of the second-order cumulant expansion, as it coincides with the parametric maximum-likelihood estimator in this scenario. We illustrate this procedure using simulations of CO, both in aqueous solution and in a carbon nanotube as a model system for biologically relevant nanoheterogeneous environments. We conclude the use of the second-order cumulant expansion permits the use of faster pulling velocities in sMD simulations, without introducing bias due to large dispersion in the non-equilibrium work distribution.
Collapse
Affiliation(s)
- Fernando Martín Boubeta
- CONICET-Facultad de Ciencias Exactas y Naturales, Instituto de Química-Física de los Materiales, Medio Ambiente y Energía, Universidad de Buenos Aires, Buenos Aires, Argentina.,Departamento de Química Inorgánica, Analítica y Química Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Rocío María Contestín García
- CONICET-Facultad de Ciencias Exactas y Naturales, Instituto de Química-Física de los Materiales, Medio Ambiente y Energía, Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Ezequiel Norberto Lorenzo
- CONICET-Facultad de Ciencias Exactas y Naturales, Instituto de Química-Física de los Materiales, Medio Ambiente y Energía, Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Leonardo Boechi
- CONICET-Facultad de Ciencias Exactas y Naturales, Instituto de Cálculo, Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Dario Estrin
- CONICET-Facultad de Ciencias Exactas y Naturales, Instituto de Química-Física de los Materiales, Medio Ambiente y Energía, Universidad de Buenos Aires, Buenos Aires, Argentina.,Departamento de Química Inorgánica, Analítica y Química Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Mariela Sued
- CONICET-Facultad de Ciencias Exactas y Naturales, Instituto de Cálculo, Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Mehrnoosh Arrar
- CONICET-Facultad de Ciencias Exactas y Naturales, Instituto de Química-Física de los Materiales, Medio Ambiente y Energía, Universidad de Buenos Aires, Buenos Aires, Argentina.,Departamento de Química Inorgánica, Analítica y Química Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina
| |
Collapse
|
3
|
Arrar M, Boubeta FM, Szretter ME, Sued M, Boechi L, Rodriguez D. On the accurate estimation of free energies using the jarzynski equality. J Comput Chem 2018; 40:688-696. [DOI: 10.1002/jcc.25754] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2018] [Revised: 06/18/2018] [Accepted: 09/03/2018] [Indexed: 11/09/2022]
Affiliation(s)
- Mehrnoosh Arrar
- Instituto de Química-Física de los Materiales, Medio Ambiente y Energía, CONICET-Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina, Departamento de Química Inorgánica, Analítica y Química Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires; Buenos Aires Argentina
| | - Fernando Martín Boubeta
- Instituto de Química-Física de los Materiales, Medio Ambiente y Energía, CONICET-Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina, Departamento de Química Inorgánica, Analítica y Química Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires; Buenos Aires Argentina
| | - Maria Eugenia Szretter
- Departamento de Matemática, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina, Instituto de Cálculo, CONICET-Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires; Buenos Aires Argentina
| | - Mariela Sued
- Instituto de Cálculo, CONICET-Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires; Buenos Aires Argentina
| | - Leonardo Boechi
- Instituto de Cálculo, CONICET-Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires; Buenos Aires Argentina
| | - Daniela Rodriguez
- Instituto de Cálculo, CONICET-Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires; Buenos Aires Argentina
| |
Collapse
|
4
|
Giovannelli E, Gellini C, Pietraperzia G, Cardini G, Procacci P, Pagliai M, Volkov V, Chelli R. Correspondence between light-absorption spectrum and nonequilibrium work distribution as a mean to access free energy differences between electronic states. J Chem Phys 2018; 149:084101. [PMID: 30193485 DOI: 10.1063/1.5028116] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The problem of recovering the free energy difference between two electronic states has been investigated by Frezzato [Chem. Phys. Lett. 533, 106 (2012)], exploring the equivalence between light-absorption spectra and work distribution, hence opening to the application of a spectroscopic version of the Jarzynski equality (JE) [Phys. Rev. Lett. 78, 2690 (1997)]. Here, assuming the validity of the time-dependent perturbation theory, we demonstrate that such equivalence does not lead to the known form of the JE. This is ascribed to the fact that light-absorption processes cannot be described as stochastic processes. To emphasize such an aspect, we devise a stochastic model for the UV-vis (ultraviolet and visible) absorption, suitable for determining the free energy difference between two generic quantum manifolds in a JE-like fashion. However, the model would require explicit knowledge of the transition dipole moments, which are in general not available. Nonetheless, we derive a spectroscopic version of the JE that allows us to recover the free energy difference between the ground and an excited electronic state when the latter state is the only one observed in the spectrum.
Collapse
Affiliation(s)
- Edoardo Giovannelli
- Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Cristina Gellini
- Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Giangaetano Pietraperzia
- Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Gianni Cardini
- Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Piero Procacci
- Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Marco Pagliai
- Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Victor Volkov
- Interdisciplinary Biomedical Research Center, School of Science and Technology, Nottingham Trent University, Clifton Lane, Nottingham NG11 8NS, United Kingdom
| | - Riccardo Chelli
- Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
5
|
Kurut A, Fonseca R, Boomsma W. Driving Structural Transitions in Molecular Simulations Using the Nonequilibrium Candidate Monte Carlo. J Phys Chem B 2018; 122:1195-1204. [PMID: 29260565 DOI: 10.1021/acs.jpcb.7b11426] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Hybrid simulation procedures which combine molecular dynamics with Monte Carlo are attracting increasing attention as tools for improving the sampling efficiency in molecular simulations. In particular, encouraging results have been reported for nonequilibrium candidate protocols, in which a Monte Carlo move is applied gradually, and interleaved with a process that equilibrates the remaining degrees of freedom. Although initial studies have uncovered a substantial potential of the method, its practical applicability for sampling structural transitions in macromolecules remains incompletely understood. Here, we address this issue by systematically investigating the efficiency of the nonequilibrium candidate Monte Carlo on the sampling of rotameric distributions of two peptide systems at atomistic resolution both in vacuum and explicit solvent. The studied systems allow us to directly probe the efficiency with which a single or a few slow degrees of freedom can be driven between well-separated free-energy minima and to explore the sensitivity of the method toward the involved free parameters. In line with results on other systems, our study suggests that order-of-magnitude gains can be obtained in certain scenarios but also identifies challenges that arise when applying the procedure in explicit solvent.
Collapse
Affiliation(s)
- Anıl Kurut
- Department of Computer Science, University of Copenhagen , 2100 Copenhagen Ø, Denmark
| | - Rasmus Fonseca
- Department of Molecular and Cellular Physiology, Stanford University , Stanford, California 94305, United States
| | - Wouter Boomsma
- Department of Computer Science, University of Copenhagen , 2100 Copenhagen Ø, Denmark
| |
Collapse
|
6
|
Giovannelli E, Procacci P, Cardini G, Pagliai M, Volkov V, Chelli R. Binding Free Energies of Host–Guest Systems by Nonequilibrium Alchemical Simulations with Constrained Dynamics: Theoretical Framework. J Chem Theory Comput 2017; 13:5874-5886. [DOI: 10.1021/acs.jctc.7b00594] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Edoardo Giovannelli
- Dipartimento
di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Piero Procacci
- Dipartimento
di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Gianni Cardini
- Dipartimento
di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Marco Pagliai
- Dipartimento
di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Victor Volkov
- Interdisciplinary
Biomedical Research Center, School of Science and Technology, Nottingham Trent University, Clifton Lane, Nottingham NG11 8NS, U.K
| | - Riccardo Chelli
- Dipartimento
di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
7
|
Giovannelli E, Cioni M, Procacci P, Cardini G, Pagliai M, Volkov V, Chelli R. Binding Free Energies of Host–Guest Systems by Nonequilibrium Alchemical Simulations with Constrained Dynamics: Illustrative Calculations and Numerical Validation. J Chem Theory Comput 2017; 13:5887-5899. [DOI: 10.1021/acs.jctc.7b00595] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Edoardo Giovannelli
- Dipartimento
di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Matteo Cioni
- Dipartimento
di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Piero Procacci
- Dipartimento
di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Gianni Cardini
- Dipartimento
di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Marco Pagliai
- Dipartimento
di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Victor Volkov
- Dipartimento
di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Riccardo Chelli
- Interdisciplinary
Biomedical Research Center, School of Science and Technology, Nottingham Trent University, Clifton Lane, Nottingham NG11 8NS, U.K
| |
Collapse
|
8
|
Zhang Z, Santos AP, Zhou Q, Liang L, Wang Q, Wu T, Franzen S. Steered molecular dynamics study of inhibitor binding in the internal binding site in dehaloperoxidase-hemoglobin. Biophys Chem 2016; 211:28-38. [DOI: 10.1016/j.bpc.2016.01.003] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2015] [Revised: 12/24/2015] [Accepted: 01/12/2016] [Indexed: 10/22/2022]
|
9
|
Giovannelli E, Cardini G, Chelli R. Elastic Barrier Dynamical Freezing in Free Energy Calculations: A Way To Speed Up Nonequilibrium Molecular Dynamics Simulations by Orders of Magnitude. J Chem Theory Comput 2016; 12:1029-39. [PMID: 26771534 DOI: 10.1021/acs.jctc.5b01117] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
An important issue concerning computer simulations addressed to free energy estimates via nonequilibrium work theorems, such as the Jarzynski equality [Phys. Rev. Lett. 1997, 78, 2690], is the computational effort required to achieve results with acceptable accuracy. In this respect, the dynamical freezing approach [Phys. Rev. E 2009, 80, 041124] has been shown to improve the efficiency of this kind of simulations, by blocking the dynamics of particles located outside an established mobility region. In this report, we show that dynamical freezing produces a systematic spurious decrease of the particle density inside the mobility region. As a consequence, the requirements to apply nonequilibrium work theorems are only approximately met. Starting from these considerations, we have developed a simulation scheme, called "elastic barrier dynamical freezing", according to which a stiff potential-energy barrier is enforced at the boundaries of the mobility region, preventing the particles from leaving this region of space during the nonequilibrium trajectories. The method, tested on the calculation of the distance-dependent free energy of a dimer immersed into a Lennard-Jones fluid, provides an accuracy comparable to the conventional steered molecular dynamics, with a computational speedup exceeding a few orders of magnitude.
Collapse
Affiliation(s)
- Edoardo Giovannelli
- Dipartimento di Chimica, Università di Firenze , Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Gianni Cardini
- Dipartimento di Chimica, Università di Firenze , Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Riccardo Chelli
- Dipartimento di Chimica, Università di Firenze , Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
10
|
Giovannelli E, Gellini C, Pietraperzia G, Cardini G, Chelli R. Nonequilibrium Candidate Monte Carlo Simulations with Configurational Freezing Schemes. J Chem Theory Comput 2015; 10:4273-83. [PMID: 26588124 DOI: 10.1021/ct500340b] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Nonequilibrium Candidate Monte Carlo simulation [Nilmeier et al., Proc. Natl. Acad. Sci. U.S.A. 2011, 108, E1009-E1018] is a tool devised to design Monte Carlo moves with high acceptance probabilities that connect uncorrelated configurations. Such moves are generated through nonequilibrium driven dynamics, producing candidate configurations accepted with a Monte Carlo-like criterion that preserves the equilibrium distribution. The probability of accepting a candidate configuration as the next sample in the Markov chain basically depends on the work performed on the system during the nonequilibrium trajectory and increases with decreasing such a work. It is thus strategically relevant to find ways of producing nonequilibrium moves with low work, namely moves where dissipation is as low as possible. This is the goal of our methodology, in which we combine Nonequilibrium Candidate Monte Carlo with Configurational Freezing schemes developed by Nicolini et al. (J. Chem. Theory Comput. 2011, 7, 582-593). The idea is to limit the configurational sampling to particles of a well-established region of the simulation sample, namely the region where dissipation occurs, while leaving fixed the other particles. This allows to make the system relaxation faster around the region perturbed by the finite-time switching move and hence to reduce the dissipated work, eventually enhancing the probability of accepting the generated move. Our combined approach enhances significantly configurational sampling, as shown by the case of a bistable dimer immersed in a dense fluid.
Collapse
Affiliation(s)
- Edoardo Giovannelli
- Dipartimento di Chimica, Università di Firenze , Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Cristina Gellini
- Dipartimento di Chimica, Università di Firenze , Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy.,European Laboratory for Nonlinear Spectroscopy (LENS), Via Nello Carrara 1, I-50019 Sesto Fiorentino, Italy
| | - Giangaetano Pietraperzia
- Dipartimento di Chimica, Università di Firenze , Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy.,European Laboratory for Nonlinear Spectroscopy (LENS), Via Nello Carrara 1, I-50019 Sesto Fiorentino, Italy
| | - Gianni Cardini
- Dipartimento di Chimica, Università di Firenze , Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy.,European Laboratory for Nonlinear Spectroscopy (LENS), Via Nello Carrara 1, I-50019 Sesto Fiorentino, Italy
| | - Riccardo Chelli
- Dipartimento di Chimica, Università di Firenze , Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy.,European Laboratory for Nonlinear Spectroscopy (LENS), Via Nello Carrara 1, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
11
|
Siders PD. Conformational free energy of alkylsilanes by nonequilibrium-pulling Monte Carlo simulation. MOLECULAR SIMULATION 2015. [DOI: 10.1080/08927022.2015.1083101] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
12
|
Giovannelli E, Cardini G, Gellini C, Pietraperzia G, Chelli R. Annealed importance sampling with constant cooling rate. J Chem Phys 2015; 142:074102. [PMID: 25701997 DOI: 10.1063/1.4907883] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Annealed importance sampling is a simulation method devised by Neal [Stat. Comput. 11, 125 (2001)] to assign weights to configurations generated by simulated annealing trajectories. In particular, the equilibrium average of a generic physical quantity can be computed by a weighted average exploiting weights and estimates of this quantity associated to the final configurations of the annealed trajectories. Here, we review annealed importance sampling from the perspective of nonequilibrium path-ensemble averages [G. E. Crooks, Phys. Rev. E 61, 2361 (2000)]. The equivalence of Neal's and Crooks' treatments highlights the generality of the method, which goes beyond the mere thermal-based protocols. Furthermore, we show that a temperature schedule based on a constant cooling rate outperforms stepwise cooling schedules and that, for a given elapsed computer time, performances of annealed importance sampling are, in general, improved by increasing the number of intermediate temperatures.
Collapse
Affiliation(s)
- Edoardo Giovannelli
- Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Gianni Cardini
- Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Cristina Gellini
- Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Giangaetano Pietraperzia
- Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| | - Riccardo Chelli
- Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
13
|
Zerbetto M, Frezzato D. Towards bulk thermodynamics via non-equilibrium methods: gaseous methane as a case study. Phys Chem Chem Phys 2015; 17:1966-79. [DOI: 10.1039/c4cp03815k] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
The equation of state of bulk materials is achieved via thermodynamic derivatives of the free energy yielded by nonequilibrium transformations and Jarzynski equality.
Collapse
Affiliation(s)
- Mirco Zerbetto
- Dipartimento di Scienze Chimiche
- Università degli Studi di Padova
- I-35131 Padova
- Italy
| | - Diego Frezzato
- Dipartimento di Scienze Chimiche
- Università degli Studi di Padova
- I-35131 Padova
- Italy
| |
Collapse
|
14
|
Zerbetto M, Piserchia A, Frezzato D. Looking for some free energy? Call JEFREE (…). J Comput Chem 2014; 35:1865-81. [DOI: 10.1002/jcc.23701] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2013] [Revised: 04/09/2014] [Accepted: 07/13/2014] [Indexed: 11/09/2022]
Affiliation(s)
- Mirco Zerbetto
- Dipartimento di Scienze Chimiche; Università degli Studi di Padova; via Marzolo 1 I-35131 Padova Italy
| | - Andrea Piserchia
- Dipartimento di Scienze Chimiche; Università degli Studi di Padova; via Marzolo 1 I-35131 Padova Italy
| | - Diego Frezzato
- Dipartimento di Scienze Chimiche; Università degli Studi di Padova; via Marzolo 1 I-35131 Padova Italy
| |
Collapse
|
15
|
Hartmann AK. High-precision work distributions for extreme nonequilibrium processes in large systems. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:052103. [PMID: 25353736 DOI: 10.1103/physreve.89.052103] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/06/2013] [Indexed: 06/04/2023]
Abstract
The distributions of work for strongly nonequilibrium processes are studied using a very general form of a large-deviation approach, which allows one to study distributions down to extremely small probabilities of almost arbitrary quantities of interest for equilibrium, nonequilibrium stationary, and even nonstationary processes. The method is applied to quickly vary the external field in a wide range B = 3 ↔ 0 for a critical (T = 2.269) two-dimensional Ising system of size L × L = 128 × 128. To obtain free-energy differences from the work distributions, they must be studied in ranges where the probabilities are as small as 10^{-240}, which is not possible using direct simulation approaches. By comparison with the exact free energies, which are available for this model for the zero-field case, one sees that the present approach allows one to obtain the free energy with a very high relative precision of 10^{-4}. This works well also for a nonzero field, i.e., for a case where standard umbrella-sampling methods are not efficient to calculate free energies. Furthermore, for the present case it is verified that the resulting distributions of work for forward and backward processes fulfill Crooks theorem with high precision. Finally, the free energy for the Ising magnet as a function of the field strength is obtained.
Collapse
|
16
|
Giovannelli E, Gellini C, Pietraperzia G, Cardini G, Chelli R. Combining path-breaking with bidirectional nonequilibrium simulations to improve efficiency in free energy calculations. J Chem Phys 2014; 140:064104. [DOI: 10.1063/1.4863999] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
|
17
|
Zhang Z, Wu T, Wang Q, Pan H, Tang R. Impact of interfacial high-density water layer on accurate estimation of adsorption free energy by Jarzynski's equality. J Chem Phys 2014; 140:034706. [DOI: 10.1063/1.4858428] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
|
18
|
Chelli R, Gellini C, Pietraperzia G, Giovannelli E, Cardini G. Path-breaking schemes for nonequilibrium free energy calculations. J Chem Phys 2014; 138:214109. [PMID: 23758360 DOI: 10.1063/1.4808037] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We propose a path-breaking route to the enhancement of unidirectional nonequilibrium simulations for the calculation of free energy differences via Jarzynski's equality [C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997)]. One of the most important limitations of unidirectional nonequilibrium simulations is the amount of realizations necessary to reach suitable convergence of the work exponential average featuring the Jarzynski's relationship. In this respect, a significant improvement of the performances could be obtained by finding a way of stopping trajectories with negligible contribution to the work exponential average, before their normal end. This is achieved using path-breaking schemes which are essentially based on periodic checks of the work dissipated during the pulling trajectories. Such schemes can be based either on breaking trajectories whose dissipated work exceeds a given threshold or on breaking trajectories with a probability increasing with the dissipated work. In both cases, the computer time needed to carry out a series of nonequilibrium trajectories is reduced up to a factor ranging from 2 to more than 10, at least for the processes under consideration in the present study. The efficiency depends on several aspects, such as the type of process, the number of check-points along the pathway and the pulling rate as well. The method is illustrated through radically different processes, i.e., the helix-coil transition of deca-alanine and the pulling of the distance between two methane molecules in water solution.
Collapse
Affiliation(s)
- Riccardo Chelli
- Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy.
| | | | | | | | | |
Collapse
|
19
|
Computing Equilibrium Free Energies Using Non-Equilibrium Molecular Dynamics. ENTROPY 2013. [DOI: 10.3390/e16010041] [Citation(s) in RCA: 50] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
20
|
Nicolini P, Frezzato D, Gellini C, Bizzarri M, Chelli R. Toward quantitative estimates of binding affinities for protein-ligand systems involving large inhibitor compounds: a steered molecular dynamics simulation route. J Comput Chem 2013; 34:1561-76. [PMID: 23620471 DOI: 10.1002/jcc.23286] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2012] [Revised: 03/01/2013] [Accepted: 03/05/2013] [Indexed: 11/09/2022]
Abstract
Understanding binding mechanisms between enzymes and potential inhibitors and quantifying protein-ligand affinities in terms of binding free energy is of primary importance in drug design studies. In this respect, several approaches based on molecular dynamics simulations, often combined with docking techniques, have been exploited to investigate the physicochemical properties of complexes of pharmaceutical interest. Even if the geometric properties of a modeled protein-ligand complex can be well predicted by computational methods, it is still challenging to rank with chemical accuracy a series of ligand analogues in a consistent way. In this article, we face this issue calculating relative binding free energies of a focal adhesion kinase, an important target for the development of anticancer drugs, with pyrrolopyrimidine-based ligands having different inhibitory power. To this aim, we employ steered molecular dynamics simulations combined with nonequilibrium work theorems for free energy calculations. This technique proves very powerful when a series of ligand analogues is considered, allowing one to tackle estimation of protein-ligand relative binding free energies in a reasonable time. In our cases, the calculated binding affinities are comparable with those recovered from experiments by exploiting the Michaelis-Menten mechanism with a competitive inhibitor.
Collapse
Affiliation(s)
- Paolo Nicolini
- Departament de Física i Enginyeria Nuclear, Universitat Politècnica de Catalunya, Campus Nord B4-B5, E-08034 Barcelona, Spain
| | | | | | | | | |
Collapse
|