1
|
Macchiagodena M, Pagliai M, Procacci P. NE-RDFE: A protocol and toolkit for computing relative dissociation free energies with GROMACS between dissimilar molecules using bidirectional nonequilibrium dual topology schemes. J Comput Chem 2023; 44:1221-1230. [PMID: 36704972 DOI: 10.1002/jcc.27077] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2022] [Revised: 12/20/2022] [Accepted: 01/07/2023] [Indexed: 01/28/2023]
Abstract
We describe a step-by-step protocol and toolkit for the computation of the relative dissociation free energy (RDFE) with the GROMACS molecular dynamics package, based on a novel bidirectional nonequilibrium alchemical approach. The proposed methodology does not require any intervention on the code and allows computing with good accuracy the RDFE between small molecules with arbitrary differences in volume, charge, and chemical topology. The procedure is illustrated for the challenging SAMPL9 batch of host-guest pairs. The article is supplemented by a detailed online tutorial, available at https://procacci.github.io/vdssb_gromacs/NE-RDFE and by a public Zenodo repository available at https://zenodo.org/record/6982932.
Collapse
Affiliation(s)
- Marina Macchiagodena
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Sesto Fiorentino, Italy
| | - Marco Pagliai
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Sesto Fiorentino, Italy
| | - Piero Procacci
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Sesto Fiorentino, Italy
| |
Collapse
|
2
|
Procacci P. Relative Binding Free Energy between Chemically Distant Compounds Using a Bidirectional Nonequilibrium Approach. J Chem Theory Comput 2022; 18:4014-4026. [PMID: 35642423 PMCID: PMC9202353 DOI: 10.1021/acs.jctc.2c00295] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
In the context of
advanced hit-to-lead drug design based on atomistic
molecular dynamics simulations, we propose a dual topology alchemical
approach for calculating the relative binding free energy (RBFE) between
two chemically distant compounds. The method (termed NE-RBFE) relies
on the enhanced sampling of the end-states in bulk and in the bound
state via Hamiltonian Replica Exchange, alchemically connected by
a series of independent and fast nonequilibrium (NE) simulations.
The technique has been implemented in a bidirectional fashion, applying
the Crooks theorem to the NE work distributions for RBFE predictions.
The dissipation of the NE process, negatively affecting accuracy,
has been minimized by introducing a smooth regularization based on
shifted electrostatic and Lennard-Jones non bonded potentials. As
a challenging testbed, we have applied our method to the calculation
of the RBFEs in the recent host–guest SAMPL international contest,
featuring a macrocyclic host with guests varying in the net charge,
volume, and chemical fingerprints. Closure validation has been successfully
verified in cycles involving compounds with disparate Tanimoto coefficients,
volume, and net charge. NE-RBFE is specifically tailored for massively
parallel facilities and can be used with little or no code modification
on most of the popular software packages supporting nonequilibrium
alchemical simulations, such as Gromacs, Amber, NAMD, or OpenMM. The
proposed methodology bypasses most of the entanglements and limitations
of the standard single topology RBFE approach for strictly congeneric
series based on free-energy perturbation, such as slowly relaxing
cavity water, sampling issues along the alchemical stratification,
and the need for highly overlapping molecular fingerprints.
Collapse
Affiliation(s)
- Piero Procacci
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
| |
Collapse
|
3
|
Sun Z, He Q. Seeding the multi-dimensional nonequilibrium pulling for Hamiltonian variation: indirect nonequilibrium free energy simulations at QM levels. Phys Chem Chem Phys 2022; 24:8800-8819. [PMID: 35352744 DOI: 10.1039/d2cp00355d] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
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 work, a multi-dimensional nonequilibrium pulling framework for Hamiltonian variations was 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. 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.
Collapse
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Institute of Theoretical and Computational Chemistry, Peking University, Beijing 100871, China.
| | - Qiaole He
- AI Department of Enzymaster (Ningbo) Bio-Engineering Co., Ltd, North Century Avenue 333, 315100 Ningbo, China
| |
Collapse
|
4
|
Procacci P, Guarnieri G. SAMPL7 blind predictions using nonequilibrium alchemical approaches. J Comput Aided Mol Des 2021; 35:37-47. [PMID: 33392950 DOI: 10.1007/s10822-020-00365-3] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2020] [Accepted: 11/27/2020] [Indexed: 11/24/2022]
Abstract
In the context of the SAMPL7 challenge, we computed, employing a non-equilibrium (NE) alchemical technique, the standard binding free energy of two series of host-guest systems, involving as a host the Isaac's TrimerTrip, a Cucurbituril-like open cavitand, and the Gilson's Cyclodextrin derivatives. The adopted NE alchemy combines enhanced sampling molecular dynamics simulations with driven fast out-of-equilibrium alchemical trajectories to recover the free energy via the Jarzynski and Crooks NE theorems. The GAFF2 non-polarizable force field was used for the parametrization. Performances were acceptable and similar in accuracy to those we submitted for Gibb's Deep Cavity Cavitands in the previous SAMPL6 host-guest challenge, confirming the reliability of the computational approach and exposing, in some cases, some important deficiencies of the GAFF2 non-polarizable force field.
Collapse
Affiliation(s)
- Piero Procacci
- University of Florence, Department of Chemistry, Via Lastruccia n. 3, 50019, Sesto Fiorentino, FI, Italy.
| | - Guido Guarnieri
- ENEA, Portici Research Centre, DTE-ICT-HPC, P.le E. Fermi, 1, 80055, Portici, NA, Italy
| |
Collapse
|
5
|
Procacci P. Methodological uncertainties in drug-receptor binding free energy predictions based on classical molecular dynamics. Curr Opin Struct Biol 2020; 67:127-134. [PMID: 33220532 DOI: 10.1016/j.sbi.2020.08.001] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2020] [Revised: 08/02/2020] [Accepted: 08/02/2020] [Indexed: 12/13/2022]
Abstract
Computational approaches are becoming an essential tool in modern drug design and discovery, with fast compound triaging using a combination of machine learning and docking techniques followed by molecular dynamics binding free energies assessment using alchemical techniques. The traditional MD-based alchemical free energy perturbation (FEP) method faces severe sampling issues that may limits its reliability in automated workflows. Here we review the major sources of uncertainty in FEP protocols for drug discovery, showing how the sampling problem can be effectively tackled by switching to nonequilibrium alchemical techniques.
Collapse
Affiliation(s)
- Piero Procacci
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, dVia della Lastruccia 3, 50019 Sesto Fiorentino, Italy.
| |
Collapse
|
6
|
Macchiagodena M, Pagliai M, Karrenbrock M, Guarnieri G, Iannone F, Procacci P. Virtual Double-System Single-Box: A Nonequilibrium Alchemical Technique for Absolute Binding Free Energy Calculations: Application to Ligands of the SARS-CoV-2 Main Protease. J Chem Theory Comput 2020; 16:7160-7172. [PMID: 33090785 PMCID: PMC8015232 DOI: 10.1021/acs.jctc.0c00634] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
In the context of drug-receptor binding affinity calculations using molecular dynamics techniques, we implemented a combination of Hamiltonian replica exchange (HREM) and a novel nonequilibrium alchemical methodology, called virtual double-system single-box, with increased accuracy, precision, and efficiency with respect to the standard nonequilibrium approaches. The method has been applied for the determination of absolute binding free energies of 16 newly designed noncovalent ligands of the main protease (3CLpro) of SARS-CoV-2. The core structures of 3CLpro ligands were previously identified using a multimodal structure-based ligand design in combination with docking techniques. The calculated binding free energies for four additional ligands with known activity (either for SARS-CoV or SARS-CoV-2 main protease) are also reported. The nature of binding in the 3CLpro active site and the involved residues besides the CYS-HYS catalytic dyad have been thoroughly characterized by enhanced sampling simulations of the bound state. We have identified several noncongeneric compounds with predicted low micromolar activity for 3CLpro inhibition, which may constitute possible lead compounds for the development of antiviral agents in Covid-19 treatment.
Collapse
Affiliation(s)
- Marina Macchiagodena
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
| | - Marco Pagliai
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
| | - Maurice Karrenbrock
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
| | - Guido Guarnieri
- ENEA, Portici Research Centre, DTE-ICT-HPC P.le E. Fermi, 1, I-80055 Portici (NA), Italy
| | - Francesco Iannone
- ENEA, Portici Research Centre, DTE-ICT-HPC P.le E. Fermi, 1, I-80055 Portici (NA), Italy
| | - Piero Procacci
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
| |
Collapse
|
7
|
Macchiagodena M, Pagliai M, Andreini C, Rosato A, Procacci P. Upgraded AMBER Force Field for Zinc-Binding Residues and Ligands for Predicting Structural Properties and Binding Affinities in Zinc-Proteins. ACS OMEGA 2020; 5:15301-15310. [PMID: 32637803 PMCID: PMC7331063 DOI: 10.1021/acsomega.0c01337] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/25/2020] [Accepted: 05/13/2020] [Indexed: 05/08/2023]
Abstract
We developed a novel force field in the context of AMBER parameterization for glutamate and aspartate zinc(II)-binding residues. The interaction between the zinc ion and the coordinating atoms is represented by a spherical nonbonded parameterization. The polarization effect due to the zinc ion has been taken into account by redefining the atomic charges on the residues through accurate quantum mechanical calculations. The new zinc-binding ASP and GLU residues, along with the CYS and HIS zinc-binding residues, parameterized in a recent work [Macchiagodena M.;J. Chem. Inf. Model.2019, 59, 3803-3816], allow users to reliably simulate 96% of the Zn-proteins available in the Protein Data Bank. The upgraded force field for zinc(II)-bound residues has been tested performing molecular dynamics simulations with an explicit solvent and comparing the structural information with experimental data for five different proteins binding zinc(II) with GLU, ASP, HIS, and CYS. We further validated our approach by evaluating the binding free energy of (R)-2-benzyl-3-nitropropanoic acid to carboxypeptidase A using a recently developed nonequilibrium alchemical method. We demonstrated that in this setting it is crucial to take into account polarization effects also on the metal-bound inhibitor.
Collapse
Affiliation(s)
- Marina Macchiagodena
- Dipartimento
di Chimica “Ugo Schiff”, Università
degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
| | - Marco Pagliai
- Dipartimento
di Chimica “Ugo Schiff”, Università
degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
| | - Claudia Andreini
- Dipartimento
di Chimica “Ugo Schiff”, Università
degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
- Magnetic
Resonance Center (CERM), Università
degli Studi di Firenze, Via L. Sacconi 6, 50019 Sesto Fiorentino, Italy
| | - Antonio Rosato
- Dipartimento
di Chimica “Ugo Schiff”, Università
degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
- Magnetic
Resonance Center (CERM), Università
degli Studi di Firenze, Via L. Sacconi 6, 50019 Sesto Fiorentino, Italy
| | - Piero Procacci
- Dipartimento
di Chimica “Ugo Schiff”, Università
degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
| |
Collapse
|
8
|
Krämer A, Hudson PS, Jones MR, Brooks BR. Multi-phase Boltzmann weighting: accounting for local inhomogeneity in molecular simulations of water-octanol partition coefficients in the SAMPL6 challenge. J Comput Aided Mol Des 2020; 34:471-483. [PMID: 32060677 PMCID: PMC8750956 DOI: 10.1007/s10822-020-00285-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2019] [Accepted: 01/08/2020] [Indexed: 02/08/2023]
Abstract
Accurately computing partition coefficients is a pivotal part of drug discovery. Specifically, octanol-water partition coefficients can provide information into hydrophobicity of drug-like molecules, as well as a de facto representation of membrane permeability. However, one challenge facing the computation of partition coefficients is the need to encapsulate various microscopic environments. These include areas of largely bulk solvent (i.e., either water or octanol) or regions where octanol is saturated with water or areas of higher salt concentration. Also, tautomeric effects require consideration. Thus, we present a Boltzmann weighting approach that incorporates transfer free energies across varying microscopic media, as well as varying tautomeric state, to compute partition coefficients in the SAMPL6 challenge.
Collapse
Affiliation(s)
- Andreas Krämer
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA.
| | - Phillip S Hudson
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
- Department of Chemistry, University of South Florida, Tampa, FL, 33620, USA
| | - Michael R Jones
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| |
Collapse
|
9
|
Procacci P. A remark on the efficiency of the double-system/single-box nonequilibrium approach in the SAMPL6 SAMPLing challenge. J Comput Aided Mol Des 2020; 34:635-639. [PMID: 32277315 DOI: 10.1007/s10822-020-00312-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2019] [Accepted: 03/30/2020] [Indexed: 11/25/2022]
Abstract
The alchemical nonequilibrium switching technique was one of several methods in the top tier of performance in the recent SAMPL6 SAMPLing challenge in both accuracy and efficiency. In this paper, in the context of nonequilibrium alchemical switching, we compare the efficiency of the double-system/single-box (DSSB) approach (used in the SAMPL6 challenges) to the standard single-system/double-box method (SSDB). Exploiting the Crooks theorem in a simple but effective test case, we analytically show that the DSSB approach is almost twice as efficient as SSDB for slow near-equilibrium switching but it gives basically no gain over the conventional SSDB approach when the variance of the work distribution exceeds few [Formula: see text], with the potential of producing artifacts and entanglements if not judiciously implemented.
Collapse
Affiliation(s)
- Piero Procacci
- Department of Chemistry, University of Florence, Via Lastruccia n. 3, 50019, Sesto Fiorentino, FI, Italy.
| |
Collapse
|
10
|
Procacci P. Precision and computational efficiency of nonequilibrium alchemical methods for computing free energies of solvation. II. Unidirectional estimates. J Chem Phys 2019; 151:144115. [DOI: 10.1063/1.5120616] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Piero Procacci
- Department of Chemistry, University of Florence, Florence, Italy
| |
Collapse
|