4.7 Article

grand: A Python Module for Grand Canonical Water Sampling in OpenMM

Journal

JOURNAL OF CHEMICAL INFORMATION AND MODELING
Volume 60, Issue 10, Pages 4436-4441

Publisher

AMER CHEMICAL SOC
DOI: 10.1021/acs.jcim.0c00648

Keywords

-

Funding

  1. EPSRC [EP/L015382/1]
  2. EPSRC [EP/J010189/1] Funding Source: UKRI

Ask authors/readers for more resources

Networks of water molecules can play a critical role at the protein-ligand interface and can directly influence drug-target interactions. Grand canonical methods aid in the sampling of these water molecules, where conventional molecular dynamics equilibration times are often long, by allowing waters to be inserted and deleted from the system, according to the chemical potential. Here, we present our open source Python module, grand (https://github.com/essex-lab/grand), which allows molecular dynamics simulations to be performed in conjunction with grand canonical Monte Carlo sampling, using the OpenMM simulation engine. We demonstrate the accuracy of this module by reproducing the density of bulk water observed from constant pressure simulations. Application of this code to the bovine pancreatic trypsin inhibitor protein reproduces three buried crystallographic water sites that are poorly sampled using conventional molecular dynamics.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Chemistry, Physical

Sensitivity of Binding Free Energy Calculations to Initial Protein Crystal Structure

Miroslav Suruzhon, Michael S. Bodnarchuk, Antonella Ciancetta, Russell Viner, Ian D. Wall, Jonathan W. Essex

Summary: This study investigated the impact of initial crystal structures on binding free energy values in alchemical free energy calculations, finding that the initial structure can significantly affect values obtained at short timescales. Rare events, such as torsional ligand motions, became important factors at longer timescales, leading to higher uncertainty in the obtained values. Optimal protocols should focus on achieving convergence in the alchemical coupling parameter space, as well as longer simulations and multiple repeats.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2021)

Article Chemistry, Medicinal

ParaMol: A Package for Automatic Parameterization of Molecular Mechanics Force Fields

Joao Morado, Paul N. Mortenson, Marcel L. Verdonk, Richard A. Ward, Jonathan W. Essex, Chris-Kriton Skylaris

Summary: The quality of force field parameterization plays a crucial role in determining the accuracy of observable properties computed from molecular mechanics simulations. ParaMol is a Python package focused on parameterizing bonded and nonbonded terms of druglike molecules by fitting to ab initio data. Through case studies, ParaMol demonstrates the ability to derive near-ideal parameters within the constraints of the functional form, while also discussing best practices and weighting methods for parameterization routes.

JOURNAL OF CHEMICAL INFORMATION AND MODELING (2021)

Article Biochemistry & Molecular Biology

Rimantadine Binds to and Inhibits the Influenza A M2 Proton Channel without Enantiomeric Specificity

Jessica L. Thomaston, Marley L. Samways, Athina Konstantinidi, Chunlong Ma, Yanmei Hu, Hannah E. Bruce Macdonald, Jun Wang, Jonathan W. Essex, William F. DeGrado, Antonios Kolocouris

Summary: Studies have shown that there are slight differences in the way (R)- and (S)-rimantadine bind to the M2 WT channel, but this does not affect the drug's efficacy or binding kinetics.

BIOCHEMISTRY (2021)

Review Chemistry, Medicinal

Computational Methods and Tools in Antimicrobial Peptide Research

Pietro G. A. Aronica, Lauren M. Reid, Nirali Desai, Jianguo Li, Stephen J. Fox, Shilpa Yadahalli, Jonathan W. Essex, Chandra S. Verma

Summary: The evolution of antibiotic-resistant bacteria has led to an increase in untreatable diseases, highlighting the urgent need for alternative treatment methods. Antimicrobial peptides (AMPs) have attracted significant interest, with their design and development aided by molecular models and computational approaches can aid experimental studies in this area.

JOURNAL OF CHEMICAL INFORMATION AND MODELING (2021)

Article Chemistry, Physical

Generation of Quantum Configurational Ensembles Using Approximate Potentials

Joao Morado, Paul N. Mortenson, J. Willem M. Nissink, Marcel L. Verdonk, Richard A. Ward, Jonathan W. Essex, Chris-Kriton Skylaris

Summary: Conformational analysis is crucial in drug design, and molecular mechanics simulation methods are used to generate ensembles of structures to provide reliable structural information. Reparameterizing the force field can generate FFs that closely reproduce QM results, and the MC acceptance rate is strongly correlated with various phase space overlap measurements, serving as a robust metric to evaluate the similarity between MM and QM levels of theory.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2021)

Article Chemistry, Physical

Enhancing Sampling of Water Rehydration on Ligand Binding: A Comparison of Techniques

Yunhui Ge, David C. Wych, Marley L. Samways, Michael E. Wall, Jonathan W. Essex, David L. Mobley

Summary: This study focuses on the accuracy and efficiency of several molecular dynamics (MD simulation-based methods) in rehydrating buried water sites. The results suggest that BLUES and grand methods enhance water sampling relative to normal MD, with grand being more robust than BLUES.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2022)

Article Chemistry, Physical

Enhancing Ligand and Protein Sampling Using Sequential Monte Carlo

Miroslav Suruzhon, Michael S. Bodnarchuk, Antonella Ciancetta, Ian D. Wall, Jonathan W. Essex

Summary: The sampling problem is a widely studied topic in computational chemistry. In this work, an alchemical variation of adaptive sequential Monte Carlo (SMC) is presented, and it is applied to various test cases, showing efficient exploration of targeted degrees of freedom. Alchemical SMC is a promising tool for preparatory exploration of systems.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2022)

Article Genetics & Heredity

A novel variant in GATM causes idiopathic renal Fanconi syndrome and predicts progression to end-stage kidney disease

Eleanor G. Seaby, Steven Turner, David J. Bunyan, Fariba Seyed-Rezai, Jonathan Essex, Rodney D. Gilbert, Sarah Ennis

Summary: This study discusses a case of a 10-year-old girl who presented with leg pain and elevated creatinine. Further investigations confirmed renal Fanconi syndrome (RFS) in the girl and identified a novel pathogenic variant in the GATM gene. The findings suggest the importance of screening the GATM gene in children with RFS, in addition to adults.

CLINICAL GENETICS (2023)

Article Chemistry, Medicinal

Water Networks in Complexes between Proteins and FDA-Approved Drugs

Marley L. Samways, Hannah E. Bruce Macdonald, Richard D. Taylor, Jonathan W. Essex

Summary: Water molecules at protein-ligand interfaces play a crucial role in drug design. Predicting the location of water molecules in the absence of a crystal structure is important, and GCMC has shown promise in accurately predicting water binding locations. This study demonstrates that GCMC can correctly predict a significant proportion of nonbulk crystallographic water sites, and the number of hydrogen bonds and electron density of the water molecules are factors that influence the accuracy of prediction.

JOURNAL OF CHEMICAL INFORMATION AND MODELING (2023)

Article Chemistry, Physical

Enhanced Grand Canonical Sampling of Occluded Water Sites Using Nonequilibrium Candidate Monte Carlo

Oliver J. Melling, Marley L. Samways, Yunhui Ge, David L. Mobley, Jonathan W. Essex

Summary: Water molecules are crucial in biomolecular systems, especially at protein-ligand interfaces. However, simulating such systems is challenging due to slow water exchange between protein and solvent. To overcome this, the authors combine grand canonical Monte Carlo (GCMC) with nonequilibrium candidate Monte Carlo (NCMC) to develop grand canonical nonequilibrium candidate Monte Carlo (GCNCMC). This approach improves water sampling efficiency and enables the exploration of new ligand-binding geometries mediated by water.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2023)

Article Multidisciplinary Sciences

Reducing affinity as a strategy to boost immunomodulatory antibody agonism

Xiaojie Yu, Christian M. Orr, H. T. Claude Chan, Sonya James, Christine A. Penfold, Jinny Kim, Tatyana Inzhelevskaya, C. Ian Mockridge, Kerry L. Cox, Jonathan W. Essex, Ivo Tews, Martin J. Glennie, Mark S. Cragg

Summary: Low affinity of immunomodulatory antibodies leads to greater activity through increased clustering, resulting in higher immune cell activation, T cell expansion, and antitumor activity. This discovery reveals a new mechanism for enhancing receptor activation across diverse receptor families and sheds light on the mechanism of antibody-mediated receptor signaling. Affinity engineering offers a rational, efficient, and tunable solution for delivering antibody-mediated receptor activity for the treatment of human disease.

NATURE (2023)

Article Chemistry, Medicinal

Does a Machine-Learned Potential Perform Better Than an Optimally Tuned Traditional Force Field? A Case Study on Fluorohydrins

Joao Morado, Paul N. Mortenson, J. Willem M. Nissink, Jonathan W. Essex, Chris-Kriton Skylaris

Summary: We present a comparative study on the performance of different molecular models in simulating the stability and properties of 10 gamma-fluorohydrins. The results show that the ANI-2x model tends to predict stronger hydrogen bonding and overstabilize global minima, while conventional force fields still play an important role in condensed-phase simulations. This study provides guidelines for the future development and application of force fields and machine learning potentials.

JOURNAL OF CHEMICAL INFORMATION AND MODELING (2023)

Article Chemistry, Physical

Comparison of Grand Canonical and Conventional Molecular Dynamics Simulation Methods for Protein-Bound Water Networks

Vilhelm Ekberg, Marley L. Samways, Majda Misini Ignjatovic, Jonathan W. Essex, Ulf Ryde

Summary: Water molecules play crucial roles in biochemical processes, so it is vital to understand their structure and properties around proteins. This study compares different computational methods for modeling water molecules and their energetics. The results show that GCMC/MD simulations accelerate the sampling and equilibration of water molecules, while GIST analysis improves the precision of energy calculations.

ACS PHYSICAL CHEMISTRY AU (2022)

Article Chemistry, Physical

The automated optimisation of a coarse-grained force field using free energy data

Javier Caceres-Delpiano, Lee-Ping Wang, Jonathan W. Essex

Summary: Atomistic models provide detailed representation of molecular systems, but may not be sufficient for simulating large systems over long timescales. Coarse-grained models accelerate simulations by reducing degrees of freedom, with some loss in accuracy. New optimization processes for parameterizing these models could improve their quality and applicability range.

PHYSICAL CHEMISTRY CHEMICAL PHYSICS (2021)

Review Chemistry, Multidisciplinary

Water molecules at protein-drug interfaces: computational prediction and analysis methods

Marley L. Samways, Richard D. Taylor, Hannah E. Bruce Macdonald, Jonathan W. Essex

Summary: The review highlights the fundamental importance of water molecules at drug-protein interfaces and discusses the challenges and advantages of experimental and computational methods to analyze their role in drug binding. It provides a critical analysis of experimental data used to validate computational methods and suggests directions for future research to address the fundamental difficulties of each method.

CHEMICAL SOCIETY REVIEWS (2021)

No Data Available