1
|
Fukunishi Y, Higo J, Kasahara K. Computer simulation of molecular recognition in biomolecular system: from in silico screening to generalized ensembles. Biophys Rev 2022; 14:1423-1447. [PMID: 36465086 PMCID: PMC9703445 DOI: 10.1007/s12551-022-01015-8] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Accepted: 11/06/2022] [Indexed: 11/29/2022] Open
Abstract
Prediction of ligand-receptor complex structure is important in both the basic science and the industry such as drug discovery. We report various computation molecular docking methods: fundamental in silico (virtual) screening, ensemble docking, enhanced sampling (generalized ensemble) methods, and other methods to improve the accuracy of the complex structure. We explain not only the merits of these methods but also their limits of application and discuss some interaction terms which are not considered in the in silico methods. In silico screening and ensemble docking are useful when one focuses on obtaining the native complex structure (the most thermodynamically stable complex). Generalized ensemble method provides a free-energy landscape, which shows the distribution of the most stable complex structure and semi-stable ones in a conformational space. Also, barriers separating those stable structures are identified. A researcher should select one of the methods according to the research aim and depending on complexity of the molecular system to be studied.
Collapse
Affiliation(s)
- Yoshifumi Fukunishi
- Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 2-3-26, Aomi, Koto-Ku, Tokyo, 135-0064 Japan
| | - Junichi Higo
- Graduate School of Information Science, University of Hyogo, 7-1-28 Minatojima Minamimachi, Chuo-Ku, Kobe, Hyogo 650-0047 Japan ,Research Organization of Science and Technology, Ritsumeikan University, 1-1-1 Noji-Higashi, Kusatsu, Shiga 525-8577 Japan
| | - Kota Kasahara
- College of Life Sciences, Ritsumeikan University, 1-1-1 Noji-Higashi, Kusatsu, Shiga 525-8577 Japan
| |
Collapse
|
2
|
Dolgonosov AM. Problems of the Theory of Ion Exchange I: Describing Forces of Ion Exchange in Classical Systems. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2022. [DOI: 10.1134/s0036024422100089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
3
|
Expression for effective dielectric permittivity of polar liquid at molecular scale. Theor Chem Acc 2022. [DOI: 10.1007/s00214-022-02909-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
|
4
|
Naghsh F, Monajjemi M. Thermodynamic Study of Assembling ↔ Disassembling of Microtubules via the Monte Carlo Simulation. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2022. [DOI: 10.1134/s0036024422070111] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
|
5
|
Molecular Mean-Field Theory of Ionic Solutions: A Poisson-Nernst-Planck-Bikerman Model. ENTROPY 2020; 22:e22050550. [PMID: 33286322 PMCID: PMC7517072 DOI: 10.3390/e22050550] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/16/2020] [Revised: 05/11/2020] [Accepted: 05/12/2020] [Indexed: 12/21/2022]
Abstract
We have developed a molecular mean-field theory—fourth-order Poisson–Nernst–Planck–Bikerman theory—for modeling ionic and water flows in biological ion channels by treating ions and water molecules of any volume and shape with interstitial voids, polarization of water, and ion-ion and ion-water correlations. The theory can also be used to study thermodynamic and electrokinetic properties of electrolyte solutions in batteries, fuel cells, nanopores, porous media including cement, geothermal brines, the oceanic system, etc. The theory can compute electric and steric energies from all atoms in a protein and all ions and water molecules in a channel pore while keeping electrolyte solutions in the extra- and intracellular baths as a continuum dielectric medium with complex properties that mimic experimental data. The theory has been verified with experiments and molecular dynamics data from the gramicidin A channel, L-type calcium channel, potassium channel, and sodium/calcium exchanger with real structures from the Protein Data Bank. It was also verified with the experimental or Monte Carlo data of electric double-layer differential capacitance and ion activities in aqueous electrolyte solutions. We give an in-depth review of the literature about the most novel properties of the theory, namely Fermi distributions of water and ions as classical particles with excluded volumes and dynamic correlations that depend on salt concentration, composition, temperature, pressure, far-field boundary conditions etc. in a complex and complicated way as reported in a wide range of experiments. The dynamic correlations are self-consistent output functions from a fourth-order differential operator that describes ion-ion and ion-water correlations, the dielectric response (permittivity) of ionic solutions, and the polarization of water molecules with a single correlation length parameter.
Collapse
|
6
|
Tikhonov DB, Zhorov BS. Intersegment contacts determine geometry of the open and closed states in P-loop channels. J Biomol Struct Dyn 2019; 38:1012-1027. [DOI: 10.1080/07391102.2019.1592026] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Denis B. Tikhonov
- Sechenov Institute of Evolutionary Physiology and Biochemistry, Russian Academy of Sciences, St. Petersburg, Russian Federation
| | - Boris S. Zhorov
- Sechenov Institute of Evolutionary Physiology and Biochemistry, Russian Academy of Sciences, St. Petersburg, Russian Federation
- Department of Biochemistry and Biomedical Sciences, McMaster University, Hamilton, Ontario, Canada
| |
Collapse
|
7
|
Liu JL, Xie D, Eisenberg B. Poisson-Fermi Formulation of Nonlocal Electrostatics in Electrolyte Solutions. COMPUTATIONAL AND MATHEMATICAL BIOPHYSICS 2017. [DOI: 10.1515/mlbmb-2017-0007] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Abstract We present a nonlocal electrostatic formulation of nonuniform ions and water molecules with interstitial voids that uses a Fermi-like distribution to account for steric and correlation efects in electrolyte solutions. The formulation is based on the volume exclusion of hard spheres leading to a steric potential and Maxwell’s displacement field with Yukawa-type interactions resulting in a nonlocal electric potential. The classical Poisson-Boltzmann model fails to describe steric and correlation effects important in a variety of chemical and biological systems, especially in high field or large concentration conditions found in and near binding sites, ion channels, and electrodes. Steric effects and correlations are apparent when we compare nonlocal Poisson-Fermi results to Poisson-Boltzmann calculations in electric double layer and to experimental measurements on the selectivity of potassium channels for K+ over Na+.
Collapse
Affiliation(s)
- Jinn-Liang Liu
- 1Institute of Computational and Modeling Science, National Tsing Hua University, Hsinchu 300, Taiwan
| | - Dexuan Xie
- 2Department of Mathematical Sciences, University of Wisconsin-Milwaukee, Milwaukee, WI 53201-0413, USA
| | - Bob Eisenberg
- 3Department of Molecular Biophysics and Physiology, Rush University, Chicago, IL 60612 USA
| |
Collapse
|
8
|
Shelley MY, Selvan ME, Zhao J, Babin V, Liao C, Li J, Shelley JC. A New Mixed All-Atom/Coarse-Grained Model: Application to Melittin Aggregation in Aqueous Solution. J Chem Theory Comput 2017; 13:3881-3897. [PMID: 28636825 PMCID: PMC5551643 DOI: 10.1021/acs.jctc.7b00071] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2017] [Indexed: 11/28/2022]
Abstract
We introduce a new mixed resolution, all-atom/coarse-grained approach (AACG), for modeling peptides in aqueous solution and apply it to characterizing the aggregation of melittin. All of the atoms in peptidic components are represented, while a single site is used for each water molecule. With the full flexibility of the peptide retained, our AACG method achieves speedups by a factor of 3-4 for CPU time reduction and another factor of roughly 7 for diffusion. An Ewald treatment permits the inclusion of long-range electrostatic interactions. These characteristics fit well with the requirements for studying peptide association and aggregation, where the system sizes and time scales require considerable computational resources with all-atom models. In particular, AACG is well suited for biologics since changes in peptide shape and long-range electrostatics may play an important role. The application of AACG to melittin, a 26-residue peptide with a well-known propensity to aggregate in solution, serves as an initial demonstration of this technology for studying peptide aggregation. We observed the formation of melittin aggregates during our simulations and characterized the time-evolution of aggregate size distribution, buried surface areas, and residue contacts. Key interactions including π-cation and π-stacking involving TRP19 were also examined. Our AACG simulations demonstrated a clear salt effect and a moderate temperature effect on aggregation and support the molten globule model of melittin aggregates. As a showcase, this work illustrates the useful role for AACG in investigations of peptide aggregation and its potential to guide formulation and design of biologics.
Collapse
Affiliation(s)
- Mee Y. Shelley
- Schrödinger,
Inc., 101 SW Main Street,
Suite 1300, Portland, Oregon 97204, United States
| | - Myvizhi Esai Selvan
- Schrödinger,
Inc., 120 W. 45th Street,
17th Floor, New York, New
York 10036, United
States
| | - Jun Zhao
- Cancer
and Inflammation Program, National Cancer
Institute, Frederick, Maryland 21702, United
States
| | - Volodymyr Babin
- Schrödinger,
Inc., 101 SW Main Street,
Suite 1300, Portland, Oregon 97204, United States
| | - Chenyi Liao
- Department
of Chemistry, University of Vermont, Burlington, Vermont 05405, United States
| | - Jianing Li
- Department
of Chemistry, University of Vermont, Burlington, Vermont 05405, United States
| | - John C. Shelley
- Schrödinger,
Inc., 101 SW Main Street,
Suite 1300, Portland, Oregon 97204, United States
| |
Collapse
|
9
|
Ozaki H, Kuratani K, Sano H, Kiyobayashi T. A Monte-Carlo simulation of ionic conductivity and viscosity of highly concentrated electrolytes based on a pseudo-lattice model. J Chem Phys 2017; 147:034904. [PMID: 28734311 DOI: 10.1063/1.4993227] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
Simulating three transport phenomena-ionic conductivity, viscosity, and self-diffusion coefficient-in a common Monte-Carlo framework, we discuss their relationship to the intermolecular interactions of electrolyte solutions at high concentrations (C/mol l-1 ∼ 1). The simulation is predicated on a pseudolattice model of the solution. The ions and solvents (collectively termed "molecules") are considered dimensionless points occupying the lattice sites. The molecular transport is realized by a repetition of swapping two adjacent molecules by the stochastic Gibbs sampling process based on simple intermolecular interactions. The framework has been validated by the fact that the simulated ionic conductivity and dynamic viscosity of 1:1- and 2:1-salts qualitatively well represent the experimental data. The magnitude of the Coulombic interaction itself is not reflected in the ionic conductivity, but the extent to which the Coulombic interaction is shielded by the dielectric constant has a significant influence. On the other hand, the dielectric constant barely influences the viscosity, while the magnitude of the Coulombic interaction is directly reflected in the viscosity.
Collapse
Affiliation(s)
- Hiroyuki Ozaki
- Department of Energy and Environment, Research Institute of Electrochemical Energy, National Institute of Advanced Industrial Science and Technology (AIST), 1-8-31 Midorigaoka, Ikeda, Osaka 563-8577, Japan
| | - Kentaro Kuratani
- Department of Energy and Environment, Research Institute of Electrochemical Energy, National Institute of Advanced Industrial Science and Technology (AIST), 1-8-31 Midorigaoka, Ikeda, Osaka 563-8577, Japan
| | - Hikaru Sano
- Department of Energy and Environment, Research Institute of Electrochemical Energy, National Institute of Advanced Industrial Science and Technology (AIST), 1-8-31 Midorigaoka, Ikeda, Osaka 563-8577, Japan
| | - Tetsu Kiyobayashi
- Department of Energy and Environment, Research Institute of Electrochemical Energy, National Institute of Advanced Industrial Science and Technology (AIST), 1-8-31 Midorigaoka, Ikeda, Osaka 563-8577, Japan
| |
Collapse
|
10
|
Ceres N, Lavery R. Improving the treatment of coarse-grain electrostatics: CVCEL. J Chem Phys 2015; 143:243118. [DOI: 10.1063/1.4933434] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- N. Ceres
- Bioinformatics: Structures and Interactions, Institut de Biologie et Chimie des Protéines, BMSSI UMR CNRS 5086/Université Lyon I, 7 Passage du Vercors, Lyon 69367, France
| | - R. Lavery
- Bioinformatics: Structures and Interactions, Institut de Biologie et Chimie des Protéines, BMSSI UMR CNRS 5086/Université Lyon I, 7 Passage du Vercors, Lyon 69367, France
| |
Collapse
|
11
|
Dinpajooh M, Matyushov DV. Free energy of ion hydration: Interface susceptibility and scaling with the ion size. J Chem Phys 2015; 143:044511. [DOI: 10.1063/1.4927570] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mohammadhasan Dinpajooh
- Department of Chemistry & Biochemistry, Arizona State University, P.O. Box 871504, Tempe, Arizona 85287, USA
| | - Dmitry V. Matyushov
- Department of Physics and Department of Chemistry & Biochemistry, Arizona State University, P.O. Box 871504, Tempe, Arizona 85287, USA
| |
Collapse
|
12
|
Sahin B, Ralf B. Nonlocal and nonlinear electrostatics of a dipolar Coulomb fluid. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2014; 26:285101. [PMID: 24920153 DOI: 10.1088/0953-8984/26/28/285101] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
We study a model Coulomb fluid consisting of dipolar solvent molecules of finite extent which generalizes the point-like dipolar Poisson-Boltzmann model (DPB) previously introduced by Coalson and Duncan (1996 J. Phys. Chem. 100 2612) and Abrashkin et al (2007 Phys. Rev. Lett. 99 077801). We formulate a nonlocal Poisson-Boltzmann equation (NLPB) and study both linear and nonlinear dielectric response in this model for the case of a single plane geometry. Our results shed light on the relevance of nonlocal versus nonlinear effects in continuum models of material electrostatics.
Collapse
Affiliation(s)
- Buyukdagli Sahin
- Biological Nanosystems, Interdisciplinary Research Institute, Université des Sciences et des Technologies de Lille (USTL), USR CNRS 3078, 50 Avenue Halley, 59568 Villeneuve d'Ascq, France
| | | |
Collapse
|
13
|
Designing electrostatic interactions in biological systems via charge optimization or combinatorial approaches: insights and challenges with a continuum electrostatic framework. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1252-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
14
|
Hassan SA, Steinbach PJ. Water-exclusion and liquid-structure forces in implicit solvation. J Phys Chem B 2011; 115:14668-82. [PMID: 22007697 PMCID: PMC3415305 DOI: 10.1021/jp208184e] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A continuum model of solvation is proposed to describe (i) long-range electrostatic effects of water exclusion resulting from incomplete and anisotropic hydration in crowded environments and (ii) short-range effects of liquid-structure forces on the hydrogen-bond interactions at solute/water interfaces. The model is an extension of the phenomenological screened coulomb potential-based implicit model of solvation. The developments reported here allow a more realistic representation of highly crowded and spatially heterogeneous environments, such as those in the interior of a living cell. Only the solvent is treated as a continuum medium. It is shown that the electrostatic effects of long-range water-exclusion can strongly affect protein-protein binding energies and are then related to the thermodynamics of complex formation. Hydrogen-bond interactions modulated by the liquid structure at interfaces are calibrated based on systematic calculations of potentials of mean force in explicit water. The electrostatic component of the model is parametrized for monovalent, divalent and trivalent ions. The conceptual and practical aspects of the model are discussed based on simulations of protein complexation and peptide folding. The current implementation is ~1.5 times slower than the gas-phase force field and exhibits good parallel performance.
Collapse
Affiliation(s)
- Sergio A Hassan
- Center for Molecular Modeling, DCB/CIT, National Institutes of Health, US DHHS, Bethesda, Maryland 20892, United States
| | | |
Collapse
|
15
|
Olsson MHM, Søndergaard CR, Rostkowski M, Jensen JH. PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. J Chem Theory Comput 2011; 7:525-37. [PMID: 26596171 DOI: 10.1021/ct100578z] [Citation(s) in RCA: 2829] [Impact Index Per Article: 217.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In this study, we have revised the rules and parameters for one of the most commonly used empirical pKa predictors, PROPKA, based on better physical description of the desolvation and dielectric response for the protein. We have introduced a new and consistent approach to interpolate the description between the previously distinct classifications into internal and surface residues, which otherwise is found to give rise to an erratic and discontinuous behavior. Since the goal of this study is to lay out the framework and validate the concept, it focuses on Asp and Glu residues where the protein pKa values and structures are assumed to be more reliable. The new and improved implementation is evaluated and discussed; it is found to agree better with experiment than the previous implementation (in parentheses): rmsd = 0.79 (0.91) for Asp and Glu, 0.75 (0.97) for Tyr, 0.65 (0.72) for Lys, and 1.00 (1.37) for His residues. The most significant advance, however, is in reducing the number of outliers and removing unreasonable sensitivity to small structural changes that arise from classifying residues as either internal or surface.
Collapse
Affiliation(s)
- Mats H M Olsson
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, Copenhagen, Denmark
| | - Chresten R Søndergaard
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, Copenhagen, Denmark
| | - Michal Rostkowski
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, Copenhagen, Denmark
| | - Jan H Jensen
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, Copenhagen, Denmark
| |
Collapse
|
16
|
Paillusson F, Blossey R. Slits, plates, and Poisson-Boltzmann theory in a local formulation of nonlocal electrostatics. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:052501. [PMID: 21230529 DOI: 10.1103/physreve.82.052501] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/01/2010] [Indexed: 05/30/2023]
Abstract
Polar liquids like water carry a characteristic nanometric length scale, the correlation length of orientation polarizations. Continuum theories that can capture this feature commonly run under the name of "nonlocal" electrostatics since their dielectric response is characterized by a scale-dependent dielectric function ε(q), where q is the wave vector; the Poisson(-Boltzmann) equation then turns into an integro-differential equation. Recently, "local" formulations have been put forward for these theories and applied to water, solvated ions, and proteins. We review the local formalism and show how it can be applied to a structured liquid in slit and plate geometries, and solve the Poisson-Boltzmann theory for a charged plate in a structured solvent with counterions. Our results establish a coherent picture of the local version of nonlocal electrostatics and show its ease of use when compared to the original formulation.
Collapse
|
17
|
Popov AV, Vorob’ev YN. GUI-BioPASED: A program for molecular dynamics simulations of biopolymers with a graphical user interface. Mol Biol 2010. [DOI: 10.1134/s0026893310040217] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
18
|
Electrostatic solvation energy for two oppositely charged ions in a solvated protein system: salt bridges can stabilize proteins. Biophys J 2010; 98:470-7. [PMID: 20141761 DOI: 10.1016/j.bpj.2009.10.031] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2009] [Revised: 10/21/2009] [Accepted: 10/22/2009] [Indexed: 11/23/2022] Open
Abstract
Born-type electrostatic continuum methods have been an indispensable ingredient in a variety of implicit-solvent methods that reduce computational effort by orders of magnitude compared to explicit-solvent MD simulations and thus enable treatment using larger systems and/or longer times. An analysis of the limitations and failures of the Born approaches serves as a guide for fundamental improvements without diminishing the importance of prior works. One of the major limitations of the Born theory is the lack of a liquidlike description of the response of solvent dipoles to the electrostatic field of the solute and the changes therein, a feature contained in the continuum Langevin-Debye (LD) model applied here to investigate how Coulombic interactions depend on the location of charges relative to the protein/water boundary. This physically more realistic LD model is applied to study the stability of salt bridges. When compared head to head using the same (independently measurable) physical parameters (radii, dielectric constants, etc.), the LD model is in good agreement with observations, whereas the Born model is grossly in error. Our calculations also suggest that a salt bridge on the protein's surface can be stabilizing when the charge separation is < or =4 A.
Collapse
|
19
|
Ha-Duong T. Protein Backbone Dynamics Simulations Using Coarse-Grained Bonded Potentials and Simplified Hydrogen Bonds. J Chem Theory Comput 2010; 6:761-73. [DOI: 10.1021/ct900408s] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Tap Ha-Duong
- Laboratoire Analyse et Modélisation pour la Biologie et l’Environnement Université d’Evry-Val-d’Essonne Rue du Pere André Jarlan, 91025 Evry Cedex, France
| |
Collapse
|
20
|
Docking flexible ligands in proteins with a solvent exposure- and distance-dependent dielectric function. J Comput Aided Mol Des 2010; 24:91-105. [PMID: 20119653 DOI: 10.1007/s10822-009-9317-9] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2009] [Accepted: 12/26/2009] [Indexed: 10/19/2022]
Abstract
Physics-based force fields for ligand-protein docking usually determine electrostatic energy with distance-dependent dielectric (DDD) functions, which do not fully account for the dielectric permittivity variance between approximately 2 in the protein core and approximately 80 in bulk water. Here we propose an atom-atom solvent exposure- and distance-dependent dielectric (SEDDD) function, which accounts for both electrostatic and dehydration energy components. Docking was performed using the ZMM program, the AMBER force field, and precomputed libraries of ligand conformers. At the seeding stage, hundreds of thousands of positions and orientations of conformers from the libraries were sampled within the rigid protein. At the refinement stage, the ten lowest-energy structures from the seeding stage were Monte Carlo-minimized with the flexible ligand and flexible protein. A search was considered a success if the root mean square deviation (RMSD) of the ligand atoms in the apparent global minimum from the x-ray structure was <2 A. Calculations on an examining set of 60 ligand-protein complexes with different DDD functions and solvent-exclusion energy revealed outliers in most of which the ligand-binding site was located at the protein surface. Using a training set of 16 ligand-protein complexes, which did not overlap with the examining set, we parameterized the SEDDD function to minimize the RMSD of the apparent global minima from the x-ray structures. Recalculation of the examining set with the SEDDD function demonstrated a 20% increase in the success rate versus the best-performing DDD function.
Collapse
|
21
|
Vitalis A, Pappu RV. ABSINTH: a new continuum solvation model for simulations of polypeptides in aqueous solutions. J Comput Chem 2009; 30:673-99. [PMID: 18506808 PMCID: PMC2670230 DOI: 10.1002/jcc.21005] [Citation(s) in RCA: 264] [Impact Index Per Article: 17.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
A new implicit solvation model for use in Monte Carlo simulations of polypeptides is introduced. The model is termed ABSINTH for self-Assembly of Biomolecules Studied by an Implicit, Novel, and Tunable Hamiltonian. It is designed primarily for simulating conformational equilibria and oligomerization reactions of intrinsically disordered proteins in aqueous solutions. The paradigm for ABSINTH is conceptually similar to the EEF1 model of Lazaridis and Karplus (Proteins 1999, 35, 133). In ABSINTH, the transfer of a polypeptide solute from the gas phase into a continuum solvent is the sum of a direct mean field interaction (DMFI), and a term to model the screening of polar interactions. Polypeptide solutes are decomposed into a set of distinct solvation groups. The DMFI is a sum of contributions from each of the solvation groups, which are analogs of model compounds. Continuum-mediated screening of electrostatic interactions is achieved using a framework similar to the one used for the DMFI. Promising results are shown for a set of test cases. These include the calculation of NMR coupling constants for short peptides, the assessment of the thermal stability of two small proteins, reversible folding of both an alpha-helix and a beta-hairpin forming peptide, and the polymeric properties of intrinsically disordered polyglutamine peptides of varying lengths. The tests reveal that the computational expense for simulations with the ABSINTH implicit solvation model increase by a factor that is in the range of 2.5-5.0 with respect to gas-phase calculations.
Collapse
Affiliation(s)
- Andreas Vitalis
- Department of Biomedical Engineering, Molecular Biophysics Program, and Center for Computational Biology, Washington University in St. Louis, One Brookings Drive, Campus Box 1097, St. Louis, Missouri 63130, USA
| | | |
Collapse
|
22
|
Gong H, Freed KF. Langevin-Debye model for nonlinear electrostatic screening of solvated ions. PHYSICAL REVIEW LETTERS 2009; 102:057603. [PMID: 19257555 DOI: 10.1103/physrevlett.102.057603] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/15/2008] [Revised: 01/06/2009] [Indexed: 05/27/2023]
Abstract
Ion-ion interactions are central to numerous phenomena in geology, biology, and material science, but current understanding is based on linear theories of limited physical applicability. Using the Langevin-Debye model, we investigate how nonlinear dielectric saturation alters screening of electrostatic interactions between solvated charges and show that charge screening in liquids strongly depends on the charge magnitudes but negligibly on the ion radii. Qualitatively different universal screening functions are predicted for ions with charges of the same or opposite signs.
Collapse
Affiliation(s)
- Haipeng Gong
- The James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA
| | | |
Collapse
|
23
|
Luttmann E, Ensign DL, Vaidyanathan V, Houston M, Rimon N, Øland J, Jayachandran G, Friedrichs M, Pande VS. Accelerating molecular dynamic simulation on the cell processor and Playstation 3. J Comput Chem 2009; 30:268-74. [PMID: 18615421 DOI: 10.1002/jcc.21054] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Implementation of molecular dynamics (MD) calculations on novel architectures will vastly increase its power to calculate the physical properties of complex systems. Herein, we detail algorithmic advances developed to accelerate MD simulations on the Cell processor, a commodity processor found in PlayStation 3 (PS3). In particular, we discuss issues regarding memory access versus computation and the types of calculations which are best suited for streaming processors such as the Cell, focusing on implicit solvation models. We conclude with a comparison of improved performance on the PS3's Cell processor over more traditional processors.
Collapse
Affiliation(s)
- Edgar Luttmann
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | | | | | | | | | | | | | | | | |
Collapse
|
24
|
Vitalis A, Pappu RV. Methods for Monte Carlo simulations of biomacromolecules. ANNUAL REPORTS IN COMPUTATIONAL CHEMISTRY 2009; 5:49-76. [PMID: 20428473 PMCID: PMC2860296 DOI: 10.1016/s1574-1400(09)00503-9] [Citation(s) in RCA: 94] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The state-of-the-art for Monte Carlo (MC) simulations of biomacromolecules is reviewed. Available methodologies for sampling conformational equilibria and associations of biomacromolecules in the canonical ensemble, given a continuum description of the solvent environment, are reviewed. Detailed sections are provided dealing with the choice of degrees of freedom, the efficiencies of MC algorithms and algorithmic peculiarities, as well as the optimization of simple movesets. The issue of introducing correlations into elementary MC moves, and the applicability of such methods to simulations of biomacromolecules is discussed. A brief discussion of multicanonical methods and an overview of recent simulation work highlighting the potential of MC methods are also provided. It is argued that MC simulations, while underutilized biomacromolecular simulation community, hold promise for simulations of complex systems and phenomena that span multiple length scales, especially when used in conjunction with implicit solvation models or other coarse graining strategies.
Collapse
Affiliation(s)
- Andreas Vitalis
- Department of Biomedical Engineering, Molecular Biophysics Program, Center for Computational Biology, Washington University in St. Louis, One Brookings Drive, Campus Box 1097, St. Louis, MO 63130-4899, USA
| | | |
Collapse
|
25
|
|
26
|
Ngan SC, Hung LH, Liu T, Samudrala R. Scoring functions for de novo protein structure prediction revisited. METHODS IN MOLECULAR BIOLOGY (CLIFTON, N.J.) 2008; 413:243-81. [PMID: 18075169 DOI: 10.1007/978-1-59745-574-9_10] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/08/2023]
Abstract
De novo protein structure prediction methods attempt to predict tertiary structures from sequences based on general principles that govern protein folding energetics and/or statistical tendencies of conformational features that native structures acquire, without the use of explicit templates. A general paradigm for de novo prediction involves sampling the conformational space, guided by scoring functions and other sequence-dependent biases, such that a large set of candidate ("decoy") structures are generated, and then selecting native-like conformations from those decoys using scoring functions as well as conformer clustering. High-resolution refinement is sometimes used as a final step to fine-tune native-like structures. There are two major classes of scoring functions. Physics-based functions are based on mathematical models describing aspects of the known physics of molecular interaction. Knowledge-based functions are formed with statistical models capturing aspects of the properties of native protein conformations. We discuss the implementation and use of some of the scoring functions from these two classes for de novo structure prediction in this chapter.
Collapse
Affiliation(s)
- Shing-Chung Ngan
- Department of Microbiology, University of Washington School of Medicine, Seattle, WA, USA
| | | | | | | |
Collapse
|
27
|
Kent A, Jha AK, Fitzgerald JE, Freed KF. Benchmarking implicit solvent folding simulations of the amyloid beta(10-35) fragment. J Phys Chem B 2008; 112:6175-86. [PMID: 18348560 PMCID: PMC2719849 DOI: 10.1021/jp077099h] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A pathogenetic feature of Alzhemier disease is the aggregation of monomeric beta-amyloid proteins (Abeta) to form oligomers. Usually these oligomers of long peptides aggregate on time scales of microseconds or longer, making computational studies using atomistic molecular dynamics models prohibitively expensive and making it essential to develop computational models that are cheaper and at the same time faithful to physical features of the process. We benchmark the ability of our implicit solvent model to describe equilibrium and dynamic properties of monomeric Abeta(10-35) using all-atom Langevin dynamics (LD) simulations, since Alphabeta(10-35) is the only fragment whose monomeric properties have been measured. The accuracy of the implicit solvent model is tested by comparing its predictions with experiment and with those from a new explicit water MD simulation, (performed using CHARMM and the TIP3P water model) which is approximately 200 times slower than the implicit water simulations. The dependence on force field is investigated by running multiple trajectories for Alphabeta(10-35) using the CHARMM, OPLS-aal, and GS-AMBER94 force fields, whereas the convergence to equilibrium is tested for each force field by beginning separate trajectories from the native NMR structure, a completely stretched structure, and from unfolded initial structures. The NMR order parameter, S2, is computed for each trajectory and is compared with experimental data to assess the best choice for treating aggregates of Alphabeta. The computed order parameters vary significantly with force field. Explicit and implicit solvent simulations using the CHARMM force fields display excellent agreement with each other and once again support the accuracy of the implicit solvent model. Alphabeta(10-35) exhibits great flexibility, consistent with experiment data for the monomer in solution, while maintaining a general strand-loop-strand motif with a solvent-exposed hydrophobic patch that is believed to be important for aggregation. Finally, equilibration of the peptide structure requires an implicit solvent LD simulation as long as 30 ns.
Collapse
Affiliation(s)
- Andrew Kent
- Department of Chemistry, The University of Chicago, Chicago, IL 60637
- The James Franck Institute, The University of Chicago, Chicago, IL 60637
| | - Abhishek K Jha
- Department of Chemistry, The University of Chicago, Chicago, IL 60637
- Institute of Biophysical Dynamics, The University of Chicago, Chicago, IL 60637
- The James Franck Institute, The University of Chicago, Chicago, IL 60637
| | - James E Fitzgerald
- Institute of Biophysical Dynamics, The University of Chicago, Chicago, IL 60637
- Department of Mathematics, The University of Chicago, Chicago, IL 60637
- Department of Physics, The University of Chicago, Chicago, IL 60637
| | - Karl F Freed
- Department of Chemistry, The University of Chicago, Chicago, IL 60637
- The James Franck Institute, The University of Chicago, Chicago, IL 60637
| |
Collapse
|
28
|
Jha AK, Freed KF. Solvation effect on conformations of 1,2:dimethoxyethane: charge-dependent nonlinear response in implicit solvent models. J Chem Phys 2008; 128:034501. [PMID: 18205504 PMCID: PMC2717614 DOI: 10.1063/1.2815764] [Citation(s) in RCA: 45] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The physical content of and, in particular, the nonlinear contributions from the Langevin-Debye model are illustrated using two applications. First, we provide an improvement in the Langevin-Debye model currently used in some implicit solvent models for computer simulations of solvation free energies of small organic molecules, as well as of biomolecular folding and binding. The analysis is based on the implementation of a charge-dependent Langevin-Debye (qLD) model that is modified by subsequent corrections due to Onsager and Kirkwood. Second, the physical content of the model is elucidated by discussing the general treatment within the LD model of the self-energy of a charge submerged in a dielectric medium for three different limiting conditions and by considering the nonlinear response of the medium. The modified qLD model is used to refine an implicit solvent model (previously applied to protein dynamics). The predictions of the modified implicit solvent model are compared with those from explicit solvent molecular dynamics simulations for the equilibrium conformational populations of 1,2-dimethoxyethane (DME), which is the shortest ether molecule to reproduce the local conformational properties of polyethylene oxide, a polymer with tremendous technological importance and a wide variety of applications. Because the conformational population preferences of DME change dramatically upon solvation, DME is a good test case to validate our modified qLD model. The present analysis of the modified qLD model provides the motivation and tools for studying a wide variety of other interesting systems with heterogeneous dielectric properties and spatial anisotropy.
Collapse
Affiliation(s)
- Abhishek K Jha
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA.
| | | |
Collapse
|
29
|
Ballenegger V, Hansen JP. Dielectric permittivity profiles of confined polar fluids. J Chem Phys 2006; 122:114711. [PMID: 15836247 DOI: 10.1063/1.1845431] [Citation(s) in RCA: 147] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The dielectric response of a simple model of a polar fluid near neutral interfaces is examined by a combination of linear response theory and extensive molecular dynamics simulations. Fluctuation expressions for a local permittivity tensor epsilon(r) are derived for planar and spherical geometries, based on the assumption of a purely local relationship between polarization and electric field. While the longitudinal component of epsilon exhibits strong oscillations on the molecular scale near interfaces, the transverse component becomes ill defined and unphysical, indicating nonlocality in the dielectric response. Both components go over to the correct bulk permittivity beyond a few molecular diameters. Upon approaching interfaces from the bulk, the permittivity tends to increase, rather than decrease as commonly assumed, and this behavior is confirmed for a simple model of water near a hydrophobic surface. An unexpected finding of the present analysis is the formation of "electrostatic double layers" signaled by a dramatic overscreening of an externally applied field inside the polar fluid close to an interface. The local electric field is of opposite sign to the external field and of significantly larger amplitude within the first layer of polar molecules.
Collapse
Affiliation(s)
- V Ballenegger
- Department of Chemistry, University of Cambridge, UK
| | | |
Collapse
|
30
|
Gunner MR, Mao J, Song Y, Kim J. Factors influencing the energetics of electron and proton transfers in proteins. What can be learned from calculations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2006; 1757:942-68. [PMID: 16905113 PMCID: PMC2760439 DOI: 10.1016/j.bbabio.2006.06.005] [Citation(s) in RCA: 82] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/07/2006] [Revised: 06/07/2006] [Accepted: 06/13/2006] [Indexed: 11/15/2022]
Abstract
A protein structure should provide the information needed to understand its observed properties. Significant progress has been made in developing accurate calculations of acid/base and oxidation/reduction reactions in proteins. Current methods and their strengths and weaknesses are discussed. The distribution and calculated ionization states in a survey of proteins is described, showing that a significant minority of acidic and basic residues are buried in the protein and that most of these remain ionized. The electrochemistry of heme and quinones are considered. Proton transfers in bacteriorhodopsin and coupled electron and proton transfers in photosynthetic reaction centers, 5-coordinate heme binding proteins and cytochrome c oxidase are highlighted as systems where calculations have provided insight into the reaction mechanism.
Collapse
Affiliation(s)
- M R Gunner
- Physics Department City College of New York, New York, NY 10031, USA.
| | | | | | | |
Collapse
|
31
|
Steinbach PJ. Exploring peptide energy landscapes: a test of force fields and implicit solvent models. Proteins 2006; 57:665-77. [PMID: 15390266 DOI: 10.1002/prot.20247] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
A biased Monte Carlo-minimization/annealing conformational search was used to characterize five descriptions of the energy landscape for each of three model systems: the 20-residue "trp-cage" miniprotein, the 20-residue "BS1" peptide, and the 17-residue "U(1-17)T9D" peptide. The EEF1 and SASA energy landscapes were studied as well as those defined by using the GB/ACE implicit water model with one of three protein force fields: CHARMM19, CHARMM22, and CHARMM22/CMAP. The lowest-energy structures of the trp-cage and BS1 peptides found for the EEF1 landscape have main-chain root-mean-square deviations (rmsds) from the respective NMR structures of less than 2 A; for U(1-17)T9D, the deviation is less than 3 A using EEF1. The main-chain rmsd of the minimum-energy trp-cage conformation obtained for the GB/ACE/CHARMM22/CMAP landscape is less than 1 A. However, this energy function strongly favored helical structures for the two peptides shown by NMR to form beta-sheet structures. Brief annealing of the system following main-chain conformational changes was found to enhance the exploration of low-energy states. The thousands of simulations reported here suggest that the prediction of protein structure might be improved by the simultaneous use of a CMAP-like description of the main chain and an EEF1-like description of the solvent.
Collapse
Affiliation(s)
- Peter J Steinbach
- Center for Molecular Modeling, National Institutes of Health, DHHS, Bethesda, Maryland 20892-5624, USA.
| |
Collapse
|
32
|
Feig M, Chocholoušová J, Tanizaki S. Extending the horizon: towards the efficient modeling of large biomolecular complexes in atomic detail. Theor Chem Acc 2005. [DOI: 10.1007/s00214-005-0062-4] [Citation(s) in RCA: 48] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
|
33
|
Vizcarra CL, Mayo SL. Electrostatics in computational protein design. Curr Opin Chem Biol 2005; 9:622-6. [PMID: 16257567 DOI: 10.1016/j.cbpa.2005.10.014] [Citation(s) in RCA: 101] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2005] [Accepted: 10/11/2005] [Indexed: 11/18/2022]
Abstract
Catalytic activity and protein-protein recognition have proven to be significant challenges for computational protein design. Electrostatic interactions are crucial for these and other protein functions, and therefore accurate modeling of electrostatics is necessary for successfully advancing protein design into the realm of protein function. This review focuses on recent progress in modeling electrostatic interactions in computational protein design, with particular emphasis on continuum models.
Collapse
Affiliation(s)
- Christina L Vizcarra
- Division of Chemistry and Chemical Engineering, Division of Biology and Howard Hughes Medical Institute, California Institute of Technology, Pasadena, California 91125, USA
| | | |
Collapse
|
34
|
Priel A, Tuszynski JA, Woolf NJ. Transitions in microtubule C-termini conformations as a possible dendritic signaling phenomenon. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2005; 35:40-52. [PMID: 16184388 DOI: 10.1007/s00249-005-0003-0] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/03/2005] [Revised: 05/12/2005] [Accepted: 05/23/2005] [Indexed: 10/25/2022]
Abstract
We model the dynamical states of the C-termini of tubulin dimers that comprise neuronal microtubules. We use molecular dynamics and other computational tools to explore the time-dependent behavior of conformational states of a C-terminus of tubulin within a microtubule and assume that each C-terminus interacts via screened Coulomb forces with the surface of a tubulin dimer, with neighboring C-termini and also with any adjacent microtubule-associated protein 2 (MAP2). Each C-terminus can either bind to the tubulin surface via one of the several positively charged regions or can be allowed to explore the space available in the solution surrounding the dimer. We find that the preferential orientation of each C-terminus is away from the tubulin surface but binding to the surface may also take place, albeit at a lower probability. The results of our model suggest that perturbations generated by the C-termini interactions with counterions surrounding a MAP2 may propagate over distances greater than those between adjacent microtubules. Thus, the MAP2 structure is able to act as a kind of biological wire (or a cable) transmitting local electrostatic perturbations resulting in ionic concentration gradients from one microtubule to another. We briefly discuss the implications the current dynamic modeling may have on synaptic activation and potentiation.
Collapse
Affiliation(s)
- Avner Priel
- Department of Physics, University of Alberta Edmonton, AB, T6G 2J1, Canada
| | | | | |
Collapse
|
35
|
Affiliation(s)
- Jacopo Tomasi
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, 56126 Pisa, Italy.
| | | | | |
Collapse
|
36
|
Feig M, Brooks CL. Recent advances in the development and application of implicit solvent models in biomolecule simulations. Curr Opin Struct Biol 2005; 14:217-24. [PMID: 15093837 DOI: 10.1016/j.sbi.2004.03.009] [Citation(s) in RCA: 403] [Impact Index Per Article: 21.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
Advances have recently been made in the development of implicit solvent methodologies and their application to the modeling of biomolecules, particularly with regard to generalized Born approaches, dielectric screening function formulations and models based on solvent-accessible surface areas. Interesting new developments include more refined non-polar solvation energy estimators, and implicit methods for modeling low-dielectric and heterogeneous environments such as membrane systems. These have been successfully applied to molecular dynamics simulations, the scoring of protein conformations, and the calculation of binding affinities and folding free energy landscapes.
Collapse
Affiliation(s)
- Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824-1319, USA
| | | |
Collapse
|
37
|
Pokala N, Handel TM. Energy functions for protein design I: efficient and accurate continuum electrostatics and solvation. Protein Sci 2004; 13:925-36. [PMID: 15010542 PMCID: PMC2280065 DOI: 10.1110/ps.03486104] [Citation(s) in RCA: 65] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2003] [Revised: 01/03/2004] [Accepted: 01/09/2004] [Indexed: 10/26/2022]
Abstract
Electrostatics and solvation energies are important for defining protein stability, structural specificity, and molecular recognition. Because these energies are difficult to compute quickly and accurately, they are often ignored or modeled very crudely in computational protein design. To address this problem, we have developed a simple, fast, and accurate approximation for calculating Born radii in the context of protein design calculations. When these approximate Born radii are used with the generalized Born continuum dielectric model, energies calculated by the 10(6)-fold slower finite difference Poisson-Boltzmann model are faithfully reproduced. A similar approach can be used for estimating solvent-accessible surface areas (SASAs). As an independent test, we show that these approximations can be used to accurately predict the experimentally determined pK(a)s of >200 ionizable groups from 15 proteins.
Collapse
Affiliation(s)
- Navin Pokala
- Department of Molecular and Cell Biology, University of California, Berkeley, 237 Hilde-brand Hall, Berkeley, CA 94720-3206, USA.
| | | |
Collapse
|
38
|
Haberthür U, Majeux N, Werner P, Caflisch A. Efficient evaluation of the effective dielectric function of a macromolecule in aqueous solution. J Comput Chem 2004; 24:1936-49. [PMID: 14515376 DOI: 10.1002/jcc.10317] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
We propose an analytical approach to calculate the effective dielectric function of proteins in aqueous solution. The screening effect if quantified by a measure of enclosure which is based on the distribution of solute atomic volumes around a pair of charges in a macromolecule. For protein conformations that vary significantly in size and shape, a comparison with finite difference Poisson calculations shows that pair interaction energies, their sums and solvation energies are well reproduced. The approach rivals the speed of simple distance dependent dielectric functions and the accuracy of the generalized Born model.
Collapse
Affiliation(s)
- Urs Haberthür
- Department of Biochemistry, University of Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | | | | | | |
Collapse
|
39
|
Abstract
A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.
Collapse
Affiliation(s)
- Themis Lazaridis
- Department of Chemistry, City College of the City University of New York, New York 10031, USA.
| |
Collapse
|