Browse
Submit Data
Databases
API
Help

Dataset Information

22 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Utilizing Grand Canonical Monte Carlo Methods in Drug Discovery.

ABSTRACT: The concepts behind targeting waters for potency and selectivity gains have been well documented and explored, although maximizing such potential gains can prove to be challenging. This problem is exacerbated in cases where there are multiple interacting waters, wherein perturbation of one water can affect the free energy landscape of the remaining waters. Knowing the right modification a priori is challenging, and computational approaches are ideally suited to help answer the key question of which substitution is best to try. Here, we use Grand Canonical Monte Carlo and the recent Grand Canonical Alchemical Perturbation methods to both understand and predict the effect of ligand-mediated water displacement when more than one water molecule is involved, as well as to understand how exploiting water networks can help govern selectivity.

SUBMITTER: Bodnarchuk MS

PROVIDER: S-EPMC6956344 | biostudies-literature | 2020 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Utilizing Grand Canonical Monte Carlo Methods in Drug Discovery.

Bodnarchuk Michael S MS Packer Martin J MJ Haywood Alexe A

ACS medicinal chemistry letters 20191211 1

The concepts behind targeting waters for potency and selectivity gains have been well documented and explored, although maximizing such potential gains can prove to be challenging. This problem is exacerbated in cases where there are multiple interacting waters, wherein perturbation of one water can affect the free energy landscape of the remaining waters. Knowing the right modification <i>a priori</i> is challenging, and computational approaches are ideally suited to help answer the key questio ...[more]

PMID: 31938467

Similar Datasets

Reactive Grand-Canonical Monte Carlo Simulations for Modeling Hydration of MgCl2.

Project description:Thermochemical heat-storage applications, based on the reversible endo-/exothermic hydration reaction of salts, are intensively investigated to search for compact heat-storage devices. To achieve a truly valuable storage system, progressively complex salts are investigated. For these salts, the equilibrium temperature and pressure conditions are not always easy to predict. However, these conditions are crucial for the design of thermochemical heat-storage systems. A biased grand-canonical Monte Carlo (GCMC) tool is developed, enabling the study of equilibrium conditions at the molecular level. The GCMC algorithm is combined with reactive force field molecular dynamics (ReaxFF), which allows bond formation within the simulation. The Weeks-Chandler-Andersen (WCA) potential is used to scan multiple trial positions for the GCMC algorithm at a small cost. The most promising trial positions can be selected for recomputation with the more expensive ReaxFF. The developed WCA-ReaxFF-GCMC tool was used to study the hydration of MgCl2·nH2O. The simulation results show a good agreement with experimental and thermodynamic equilibriums for multiple hydration levels. The hydration shows that water, present at the surface of crystalline salt, deforms the surface layers and promotes further hydration of these deformed layers. Additionally, the WCA-ReaxFF-GCMC algorithm can be used to study other, non-TCM-related, reactive sorption processes.

| S-EPMC8655925 | biostudies-literature

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

Project description:Water molecules play a key role in many biomolecular systems, particularly when bound at protein-ligand interfaces. However, molecular simulation studies on such systems are hampered by the relatively long time scales over which water exchange between a protein and solvent takes place. Grand canonical Monte Carlo (GCMC) is a simulation technique that avoids this issue by attempting the insertion and deletion of water molecules within a given structure. The approach is constrained by low acceptance probabilities for insertions in congested systems, however. To address this issue, here, we combine GCMC with nonequilibium candidate Monte Carlo (NCMC) to yield a method that we refer to as grand canonical nonequilibrium candidate Monte Carlo (GCNCMC), in which the water insertions and deletions are carried out in a gradual, nonequilibrium fashion. We validate this new approach by comparing GCNCMC and GCMC simulations of bulk water and three protein binding sites. We find that not only is the efficiency of the water sampling improved by GCNCMC but that it also results in increased sampling of ligand conformations in a protein binding site, revealing new water-mediated ligand-binding geometries that are not observed using alternative enhanced sampling techniques.

| S-EPMC9933432 | biostudies-literature

Replica-Exchange and Standard State Binding Free Energies with Grand Canonical Monte Carlo.

Project description:The ability of grand canonical Monte Carlo (GCMC) to create and annihilate molecules in a given region greatly aids the identification of water sites and water binding free energies in protein cavities. However, acceptance rates without the application of biased moves can be low, resulting in large variations in the observed water occupancies. Here, we show that replica-exchange of the chemical potential significantly reduces the variance of the GCMC data. This improvement comes at a negligible increase in computational expense when simulations comprise of runs at different chemical potentials. Replica-exchange GCMC is also found to substantially increase the precision of water binding free energies as calculated with grand canonical integration, which has allowed us to address a missing standard state correction.

| S-EPMC5729546 | biostudies-literature

Validating the Water Flooding Approach by Comparing It to Grand Canonical Monte Carlo Simulations.

Project description:The study of the function of proteins on a quantitative level requires consideration of the water molecules in and around the protein. This requirement presents a major computational challenge due to the fact that the insertion of water molecules can have a very high activation barrier and would require a long simulation time. Recently, we developed a water flooding (WF) approach which is based on a postprocessing Monte Carlo ranking of possible water configurations. This approach appears to provide a very effective way for assessing the insertion free energies and determining the most likely configurations of the internal water molecules. Although the WF approach was used effectively in modeling challenging systems that have not been addressed reliably by other microscopic approaches, it was not validated by a comparison to the more rigorous grand canonical Monte Carlo (GCMC) method. Here we validate the WF approach by comparing its performance to that of the GCMC method. It is found that the WF approach reproduces the GCMC results in well-defined test cases but does so much faster. This established the WF approach as a useful strategy for finding correct water configurations in proteins and thus to provide a powerful way for studies of the functions of proteins.

| S-EPMC5954165 | biostudies-literature

Large-scale analysis of water stability in bromodomain binding pockets with grand canonical Monte Carlo.

Project description:Conserved water molecules are of interest in drug design, as displacement of such waters can lead to higher affinity ligands and in some cases, contribute towards selectivity. Bromodomains, small protein domains involved in the epigenetic regulation of gene transcription, display a network of four conserved water molecules in their binding pockets and have recently been the focus of intense medicinal chemistry efforts. Understanding why certain bromodomains have displaceable water molecules and others do not is extremely challenging, and it remains unclear which water molecules in a given bromodomain can be targeted for displacement. Here we estimate the stability of the conserved water molecules in 35 bromodomains via binding free energy calculations using all-atom grand canonical Monte Carlo simulations. Encouraging quantitative agreement to the available experimental evidence is found. We thus discuss the expected ease of water displacement in different bromodomains and the implications for ligand selectivity.

| S-EPMC5978690 | biostudies-literature

Determination of Ionic Hydration Free Energies with Grand Canonical Monte Carlo/Molecular Dynamics Simulations in Explicit Water.

Project description:Grand canonical Monte Carlo (GCMC) simulations of ionic solutions with explicit solvent models are known to be challenging. One challenge arises from the treatment of long-range electrostatics and finite-box size in Monte Carlo simulations when periodic boundary condition and Ewald summation methods are used. Another challenge is that constant excess chemical potential GCMC simulations for charged solutes suffer from inadequate insertion and deletion acceptance ratios. In this work, we address those problems by implementing an oscillating excess chemical potential GCMC algorithm with smooth particle mesh Ewald and finite-box-size corrections to treat the long-range electrostatics. The developed GCMC simulation program was combined with GROMACS to perform GCMC/MD simulations of ionic solutions individually containing Li+, Na+, K+, Rb+, Cs+, F-, Cl-, Br-, I-, Ca2+, and Mg2+, respectively. Our simulation results show that the combined GCMC/MD approach can approximate the ionic hydration free energies with proper treatment of long-range electrostatics. Our developed simulation approach can open up new avenues for simulating complex chemical and biomolecular systems and for drug discovery.

| S-EPMC6195813 | biostudies-literature

Screening metal-organic frameworks for adsorption-driven osmotic heat engines via grand canonical Monte Carlo simulations and machine learning.

Project description:Adsorption-driven osmotic heat engines offer an alternative way for harvesting low-grade waste heat below 80°C. In this study, we performed a high-throughput computational screening based on grand canonical Monte Carlo simulations to identify the high-performance metal-organic frameworks (MOFs) from 1322 computationally ready experimental MOF structures for adsorption-driven osmotic heat engines with LiCl-methanol as the working fluid. Structure-property relationship analysis reveals that MOFs exhibiting high energy efficiency possess large working capacity, pore size and surface area, and moderate adsorption enthalpy comparable to the evaporation enthalpy. Furthermore, machine learning is employed to accelerate the computational screening for satisfied MOFs via the structure properties. The optimal structure properties of the MOFs are further identified via the ensemble-based regression model by optimizing the energy efficiency via the genetic algorithm, which shed light on rationally designing and fabricating MOFs for desired heat-to-electricity conversion.

| S-EPMC7772570 | biostudies-literature

Influence of capillary bridge formation onto the silica nanoparticle interaction studied by grand canonical Monte Carlo simulations.

Project description:Adhesion forces between nanoparticles strongly depend on the amount of adsorbed condensed water from ambient atmosphere. Liquid water forms bridges in the cavities separating the particles, giving rise to the so-called capillary forces which in most cases dominate the van der Waals and long-range electrostatic interactions. Capillary forces promote the undesirable agglomeration of particles to large clusters, thereby hindering the flowability of dry powders in process containers. In process engineering macroscopic theories based on the Laplace pressures are used to estimate the strength of the capillary forces. However, especially for low relative humidity and when the wetting of rough or small nanoparticles is studied, those theories can fail. Molecular dynamic simulations can help to give better insight into the water-particle interface. The simulated force versus distance curve as well as adhesion forces and the adsorption isotherm for silica nanoparticles at varying relative humidity will be discussed in comparison to experiments, theories, and simulations.

| S-EPMC3793898 | biostudies-other

Ion selectivity of alpha-hemolysin with beta-cyclodextrin adapter. II. Multi-ion effects studied with grand canonical Monte Carlo/Brownian dynamics simulations.

Project description:In a previous study of ion selectivity of alpha-hemolysin (alphaHL) in complex with beta-cyclodextrin (betaCD) adapter, we calculated the potential of mean force (PMF) and characterized the self-diffusion coefficients of isolated K(+) and Cl(-) ions using molecular dynamics simulations (Y. Luo et al., "Ion Selectivity of alpha-Hemolysin with beta-Cyclodextrin Adapter: I. Single Ion Potential of Mean Force and Diffusion Coefficient"). In the present effort, these results pertaining to single isolated ions in the wide aqueous pore are extended to take into account multi-ion effects. The grand canonical Monte Carlo/Brownian dynamics (GCMC/BD) algorithm is used to simulate ion currents through the wild-type alphaHL ion channel, as well as two engineered alphaHL mutants, with and without the cyclic oligosaccaride betaCD lodged in the lumen of the pore. The GCMC/BD current-voltage curves agree well with experimental results and show that betaCD increases the anion selectivity of alphaHL. Comparisons between multi-ion PMFs from GCMC/BD simulations and single-ion PMFs demonstrate that multi-ion effects and pore shape are crucial for explaining this behavior. It is concluded that the narrow betaCD adapter increases the anion selectivity of alphaHL because it reduces the pore radius locally, which decreases the ionic screening and the dielectric shielding of the strong electrostatic field induced by a nearby ring of positively charged alphaHL side chains.

| S-EPMC2843906 | biostudies-literature

Impact of ribosomal modification on the binding of the antibiotic telithromycin using a combined grand canonical monte carlo/molecular dynamics simulation approach.

Project description:Resistance to macrolide antibiotics is conferred by mutation of A2058 to G or methylation by Erm methyltransferases of the exocyclic N6 of A2058 (E. coli numbering) that forms the macrolide binding site in the 50S subunit of the ribosome. Ketolides such as telithromycin mitigate A2058G resistance yet remain susceptible to Erm-based resistance. Molecular details associated with macrolide resistance due to the A2058G mutation and methylation at N6 of A2058 by Erm methyltransferases were investigated using empirical force field-based simulations. To address the buried nature of the macrolide binding site, the number of waters within the pocket was allowed to fluctuate via the use of a Grand Canonical Monte Carlo (GCMC) methodology. The GCMC water insertion/deletion steps were alternated with Molecular Dynamics (MD) simulations to allow for relaxation of the entire system. From this GCMC/MD approach information on the interactions between telithromycin and the 50S ribosome was obtained. In the wild-type (WT) ribosome, the 2'-OH to A2058 N1 hydrogen bond samples short distances with a higher probability, while the effectiveness of telithromycin against the A2058G mutation is explained by a rearrangement of the hydrogen bonding pattern of the 2'-OH to 2058 that maintains the overall antibiotic-ribosome interactions. In both the WT and A2058G mutation there is significant flexibility in telithromycin's imidazole-pyridine side chain (ARM), indicating that entropic effects contribute to the binding affinity. Methylated ribosomes show lower sampling of short 2'-OH to 2058 distances and also demonstrate enhanced G2057-A2058 stacking leading to disrupted A752-U2609 Watson-Crick (WC) interactions as well as hydrogen bonding between telithromycin's ARM and U2609. This information will be of utility in the rational design of novel macrolide analogs with improved activity against methylated A2058 ribosomes.

| S-EPMC3681621 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data