Lattice‐sum methods for calculating electrostatic interactions in molecular simulations

1995 ◽  
Vol 103 (8) ◽  
pp. 3014-3021 ◽  
Author(s):  
Brock A. Luty ◽  
Ilario G. Tironi ◽  
Wilfred F. van Gunsteren
Keyword(s):  
Electrostatic Interactions ◽  
Molecular Simulations ◽  
Lattice Sum

scholarly journals Long-ranged contributions to solvation free energies from theory and short-ranged models

2016 ◽  
Vol 113 (11) ◽  
pp. 2819-2826 ◽  
Author(s):  
Richard C. Remsing ◽  
Shule Liu ◽  
John D. Weeks
Keyword(s):  
Electrostatic Interactions ◽  
Large Contribution ◽  
Thermodynamic Quantities ◽  
Coulomb Interactions ◽  
Free Energies ◽  
Lattice Sum ◽  
Solvation Thermodynamics ◽  
Simulation Cell ◽  
Solvation Free Energies ◽  
Ionic Solutes

Long-standing problems associated with long-ranged electrostatic interactions have plagued theory and simulation alike. Traditional lattice sum (Ewald-like) treatments of Coulomb interactions add significant overhead to computer simulations and can produce artifacts from spurious interactions between simulation cell images. These subtle issues become particularly apparent when estimating thermodynamic quantities, such as free energies of solvation in charged and polar systems, to which long-ranged Coulomb interactions typically make a large contribution. In this paper, we develop a framework for determining very accurate solvation free energies of systems with long-ranged interactions from models that interact with purely short-ranged potentials. Our approach is generally applicable and can be combined with existing computational and theoretical techniques for estimating solvation thermodynamics. We demonstrate the utility of our approach by examining the hydration thermodynamics of hydrophobic and ionic solutes and the solvation of a large, highly charged colloid that exhibits overcharging, a complex nonlinear electrostatic phenomenon whereby counterions from the solvent effectively overscreen and locally invert the integrated charge of the solvated object.

scholarly journals Enhancing Sidechain Rotamer Sampling Using Non-Equilibrium Candidate Monte Carlo

2019 ◽  
Author(s):  
Kalistyn H. Burley ◽  
Samuel C. Gill ◽  
Nathan M. Lim ◽  
David Mobley
Keyword(s):  
Molecular Dynamics ◽  
Monte Carlo ◽  
Ligand Binding ◽  
Electrostatic Interactions ◽  
Molecular Simulations ◽  
Adequate Sampling ◽  
Binding Pockets ◽  
Sidechain Rotamers ◽  
Dynamics Simulations ◽  
Non Equilibrium

<div>Molecular simulations are a valuable tool for studying biomolecular motions and thermodynamics. However, such motions can be slow compared to simulation timescales, yet critical. Specifically, adequate sampling of sidechain motions in protein binding pockets proves crucial for obtaining accurate estimates of ligand binding free energies from molecular simulations. The timescale of sidechain rotamer flips can range from a few ps to several hundred ns or longer, particularly in crowded environments like the interior of proteins. Here, we apply a mixed non-equilibrium candidate Monte Carlo (NCMC)/molecular dynamics (MD) method to enhance sampling of sidechain rotamers. The NCMC portion of our method applies a switching protocol wherein the steric and electrostatic interactions between target sidechain atoms and the surrounding environment are cycled off and then back on during the course of a move proposal. Between NCMC move proposals, simulation of the system continues via traditional molecular dynamics. Here, we first validate this approach on a simple, solvated valine-alanine dipeptide system and then apply it to a well-studied model ligand binding site in T4 lysozyme L99A. We compute the rate of rotamer transitions for a valine sidechain using our approach and compare it to that of traditional molecular dynamics simulations. Here, we show that our NCMC/MD method substantially enhances sidechain sampling, especially in systems where the torsional barrier to rotation is high (>10 kcal/mol). These barriers can be intrinsic torsional barriers or steric barriers imposed by the environment.</div><div>Overall, this may provide a promising strategy to selectively improve sidechain sampling in molecular simulations.</div>

scholarly journals Enhancing Sidechain Rotamer Sampling Using Non-Equilibrium Candidate Monte Carlo

2018 ◽  
Author(s):  
Kalistyn H. Burley ◽  
Samuel C. Gill ◽  
Nathan M. Lim ◽  
David Mobley
Keyword(s):  
Molecular Dynamics ◽  
Monte Carlo ◽  
Ligand Binding ◽  
Electrostatic Interactions ◽  
Molecular Simulations ◽  
Adequate Sampling ◽  
Binding Pockets ◽  
Sidechain Rotamers ◽  
Dynamics Simulations ◽  
Non Equilibrium

<div>Molecular simulations are a valuable tool for studying biomolecular motions and thermodynamics. However, such motions can be slow compared to simulation timescales, yet critical. Specifically, adequate sampling of sidechain motions in protein binding pockets proves crucial for obtaining accurate estimates of ligand binding free energies from molecular simulations. The timescale of sidechain rotamer flips can range from a few ps to several hundred ns or longer, particularly in crowded environments like the interior of proteins. Here, we apply a mixed non-equilibrium candidate Monte Carlo (NCMC)/molecular dynamics (MD) method to enhance sampling of sidechain rotamers. The NCMC portion of our method applies a switching protocol wherein the steric and electrostatic interactions between target sidechain atoms and the surrounding environment are cycled off and then back on during the course of a move proposal. Between NCMC move proposals, simulation of the system continues via traditional molecular dynamics. Here, we first validate this approach on a simple, solvated valine-alanine dipeptide system and then apply it to a well-studied model ligand binding site in T4 lysozyme L99A. We compute the rate of rotamer transitions for a valine sidechain using our approach and compare it to that of traditional molecular dynamics simulations. Here, we show that our NCMC/MD method substantially enhances sidechain sampling, especially in systems where the torsional barrier to rotation is high (>10 kcal/mol). These barriers can be intrinsic torsional barriers or steric barriers imposed by the environment.</div><div>Overall, this may provide a promising strategy to selectively improve sidechain sampling in molecular simulations.</div>

scholarly journals Enhancing Sidechain Rotamer Sampling Using Non-Equilibrium Candidate Monte Carlo

2019 ◽  
Author(s):  
Kalistyn H. Burley ◽  
Samuel C. Gill ◽  
Nathan M. Lim ◽  
David Mobley
Keyword(s):  
Molecular Dynamics ◽  
Monte Carlo ◽  
Ligand Binding ◽  
Electrostatic Interactions ◽  
Molecular Simulations ◽  
Adequate Sampling ◽  
Binding Pockets ◽  
Sidechain Rotamers ◽  
Dynamics Simulations ◽  
Non Equilibrium

<div>Molecular simulations are a valuable tool for studying biomolecular motions and thermodynamics. However, such motions can be slow compared to simulation timescales, yet critical. Specifically, adequate sampling of sidechain motions in protein binding pockets proves crucial for obtaining accurate estimates of ligand binding free energies from molecular simulations. The timescale of sidechain rotamer flips can range from a few ps to several hundred ns or longer, particularly in crowded environments like the interior of proteins. Here, we apply a mixed non-equilibrium candidate Monte Carlo (NCMC)/molecular dynamics (MD) method to enhance sampling of sidechain rotamers. The NCMC portion of our method applies a switching protocol wherein the steric and electrostatic interactions between target sidechain atoms and the surrounding environment are cycled off and then back on during the course of a move proposal. Between NCMC move proposals, simulation of the system continues via traditional molecular dynamics. Here, we first validate this approach on a simple, solvated valine-alanine dipeptide system and then apply it to a well-studied model ligand binding site in T4 lysozyme L99A. We compute the rate of rotamer transitions for a valine sidechain using our approach and compare it to that of traditional molecular dynamics simulations. Here, we show that our NCMC/MD method substantially enhances sidechain sampling, especially in systems where the torsional barrier to rotation is high (>10 kcal/mol). These barriers can be intrinsic torsional barriers or steric barriers imposed by the environment.</div><div>Overall, this may provide a promising strategy to selectively improve sidechain sampling in molecular simulations.</div>

Sign in / Sign up

Export Citation Format

Share Document