51
|
Kamiya N, Kayanuma M, Fujitani H, Shinoda K. A New Lipid Force Field (FUJI). J Chem Theory Comput 2020; 16:3664-3676. [PMID: 32384238 DOI: 10.1021/acs.jctc.9b01195] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
To explore inhomogeneous and anisotropic systems such as lipid bilayers, the Lennard-Jones particle mesh Ewald (LJ-PME) method has been applied without a conventional isotropic dispersion correction. As the popular AMBER and CHARMM lipid force fields were developed using a cutoff scheme, their lipid bilayers unacceptably shrink when using the LJ-PME method. In this study, a new all-atom lipid force field (FUJI) was developed on the basis of the AMBER force-field scheme including the Lipid14 van der Waals parameters. Point charges were calculated using the restrained electrostatic potentials of many lipid conformers. Further, torsion energy profiles were calculated using high-level ab initio molecular orbitals (LCCSD(T)/Aug-cc-pVTZ//LMP2/Aug-cc-pVTZ), following which the molecular mechanical dihedral parameters were derived through a fast Fourier transform. By incorporation of these parameters into a new lipid force field without fitting experimental data, the desired lipid characteristics such as the area per lipid and lateral diffusion coefficients were obtained through GROMACS molecular dynamics simulations using the LJ-PME method and virtual hydrogen sites. The calculated area per lipid and lateral diffusion coefficients showed satisfactory agreement with experimental data. Furthermore, the electron-density profiles along the membrane normal were calculated for pure lipid bilayers, and the resulting membrane thicknesses agreed well with the experimental values. As the new lipid force field is compatible with FUJI for protein and small molecules, the new FUJI force field will offer accurate modeling for complex systems consisting of various membrane proteins and lipids.
Collapse
Affiliation(s)
- Nozomu Kamiya
- Fujitsu Limited Bio-IT R&D Office, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Megumi Kayanuma
- Research Center for Computational Design of Advanced Functional Materials, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Hideaki Fujitani
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Keiko Shinoda
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| |
Collapse
|
52
|
Klein T, Lenahan FD, Kerscher M, Rausch MH, Economou IG, Koller TM, Fröba AP. Characterization of Long Linear and Branched Alkanes and Alcohols for Temperatures up to 573.15 K by Surface Light Scattering and Molecular Dynamics Simulations. J Phys Chem B 2020; 124:4146-4163. [DOI: 10.1021/acs.jpcb.0c01740] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Tobias Klein
- Institute of Advanced Optical Technologies − Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| | - Frances D. Lenahan
- Institute of Advanced Optical Technologies − Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| | - Manuel Kerscher
- Institute of Advanced Optical Technologies − Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| | - Michael H. Rausch
- Institute of Advanced Optical Technologies − Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| | - Ioannis G. Economou
- National Center for Scientific Research “Demokritos”, Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modelling of Materials Laboratory, GR-15310 Aghia Paraskevi Attikis, Greece
- Texas A&M University at Qatar, Chemical Engineering Program, Education City, PO Box 23874, Doha, Qatar
| | - Thomas M. Koller
- Institute of Advanced Optical Technologies − Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| | - Andreas P. Fröba
- Institute of Advanced Optical Technologies − Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| |
Collapse
|
53
|
Petrlova J, Petruk G, Huber RG, McBurnie EW, van der Plas MJA, Bond PJ, Puthia M, Schmidtchen A. Thrombin-derived C-terminal fragments aggregate and scavenge bacteria and their proinflammatory products. J Biol Chem 2020; 295:3417-3430. [PMID: 32034093 PMCID: PMC7076200 DOI: 10.1074/jbc.ra120.012741] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2020] [Revised: 02/04/2020] [Indexed: 12/15/2022] Open
Abstract
Thrombin-derived C-terminal peptides (TCPs), including a major 11-kDa fragment (TCP96), are produced through cleavage by human neutrophil elastase and aggregate lipopolysaccharide (LPS) and the Gram-negative bacterium Escherichia coli. However, the physiological roles of TCP96 in controlling bacterial infections and reducing LPS-induced inflammation are unclear. Here, using various biophysical methods, in silico molecular modeling, microbiological and cellular assays, and animal models, we examined the structural features and functional roles of recombinant TCP96 (rTCP96) in the aggregation of multiple bacteria and the Toll-like receptor (TLR) agonists they produce. We found that rTCP96 aggregates both Gram-negative and Gram-positive bacteria, including Staphylococcus aureus and Pseudomonas aeruginosa, and their cell-wall components LPS, lipid A, and lipoteichoic acid (LTA). The Gram-negative bacteria E. coli and P. aeruginosa were particularly sensitive to aggregation-induced bacterial permeabilization and killing. As a proof of concept, we show that rTCP96 reduces LPS-induced NF-κB activation in human monocytes, as well as in mouse models of LPS-induced subcutaneous inflammation. Moreover, in a mouse model of subcutaneous inoculation with P. aeruginosa, rTCP96 reduced bacterial levels. Together, these results link TCP-mediated aggregation of endotoxins and bacteria in vitro to attenuation of inflammation and bacterial levels in vivo.
Collapse
Affiliation(s)
- Jitka Petrlova
- Department of Clinical Sciences, Division of Dermatology and Venereology, Lund University, Lund SE-22184, Sweden.
| | - Ganna Petruk
- Department of Clinical Sciences, Division of Dermatology and Venereology, Lund University, Lund SE-22184, Sweden
| | | | - Eilish W McBurnie
- Bioinformatics Institute (A*STAR), Singapore SG-138671; Department of Chemistry, University of Southampton, Southampton UK-SO17 1BJ, United Kingdom
| | - Mariena J A van der Plas
- Department of Clinical Sciences, Division of Dermatology and Venereology, Lund University, Lund SE-22184, Sweden; Department of Pharmacy, University of Copenhagen, Copenhagen DK-2100, Denmark
| | - Peter J Bond
- Bioinformatics Institute (A*STAR), Singapore SG-138671; Department of Biological Sciences, National University of Singapore, Singapore SG-117558
| | - Manoj Puthia
- Department of Clinical Sciences, Division of Dermatology and Venereology, Lund University, Lund SE-22184, Sweden
| | - Artur Schmidtchen
- Department of Clinical Sciences, Division of Dermatology and Venereology, Lund University, Lund SE-22184, Sweden; Department of Biomedical Sciences, Copenhagen Wound Healing Center, Bispebjerg Hospital, University of Copenhagen, Copenhagen DK-2400, Denmark
| |
Collapse
|
54
|
Ding Y, Xu Y, Qian C, Chen J, Zhu J, Huang H, Shi Y, Huang J. Predicting partition coefficients of drug-like molecules in the SAMPL6 challenge with Drude polarizable force fields. J Comput Aided Mol Des 2020; 34:421-435. [PMID: 31960252 DOI: 10.1007/s10822-020-00282-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2019] [Accepted: 01/08/2020] [Indexed: 12/14/2022]
Abstract
The water-octanol partition coefficient is an important physicochemical property for small molecule drug design. Here, we report our participation in the SAMPL6 logP prediction challenge with free energy perturbation (FEP) calculations in the water phase and in the 1-octanol phase using Drude polarizable force fields. Root mean square error (RMSE) and mean absolute error (MAE) of our prediction are equal to 1.85 and 1.25 logP units. The errors are not evenly distributed. Out of eleven SAMPL6 solutes, FEP/Drude performed very badly on three molecules (deviations all larger than 2 logP units) but good on the remaining eight (deviations all less than 1 logP unit). We find while FEP converges well within one nanosecond in water, simulations in 1-octanol need much longer simulation time and possibly more independent runs for sampling. We also find out that 1-octanol, albeit being a non-polar solvent, still polarizes solute molecules and forms stable hydrogen bonds with them. At the end, we attempt to reweight FEP trajectories with QM/Drude calculations and discuss possible caveats in our simulation setup.
Collapse
Affiliation(s)
- Ye Ding
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - You Xu
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Cheng Qian
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Jinfeng Chen
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Jian Zhu
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Houhou Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Yi Shi
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Jing Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China. .,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.
| |
Collapse
|
55
|
Kubincová A, Riniker S, Hünenberger PH. Reaction-field electrostatics in molecular dynamics simulations: development of a conservative scheme compatible with an atomic cutoff. Phys Chem Chem Phys 2020; 22:26419-26437. [DOI: 10.1039/d0cp03835k] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Shifting and switching schemes are developed to enable strict energy conservation in molecular dynamics simulations relying on reaction-field electrostatic (as well as Lennard-Jones) interactions with an atom-based cutoff truncation.
Collapse
Affiliation(s)
| | - Sereina Riniker
- Laboratory of Physical Chemistry
- ETH Zurich
- 8093 Zurich
- Switzerland
| | | |
Collapse
|
56
|
Abstract
Low resolution coarse-grained (CG) models are widely adopted for investigating phenomena that cannot be effectively simulated with all-atom (AA) models. Since the development of the many-body dissipative particle dynamics method, CG models have increasingly supplemented conventional pair potentials with one-body potentials of the local density (LD) around each site. These LD potentials appear to significantly extend the transferability of CG models, while also enabling more accurate descriptions of thermodynamic properties, interfacial phenomena, and many-body correlations. In this work, we systematically examine the properties of LD potentials. We first derive and numerically demonstrate a nontrivial transformation of pair and LD potentials that leaves the total forces and equilibrium distribution invariant. Consequently, the pair and LD potentials determined via bottom-up methods are not unique. We then investigate the sensitivity of CG models for glycerol to the weighting function employed for defining the local density. We employ the multiscale coarse-graining (MS-CG) method to simultaneously parameterize both pair and LD potentials. When employing a short-ranged Lucy function that defines the local density from the first solvation shell, the MS-CG model accurately reproduces the pair structure, pressure-density equation of state, and liquid-vapor interfacial profile of the AA model. The accuracy of the model generally decreases as the range of the Lucy function increases further. The MS-CG model provides similar accuracy when a smoothed Heaviside function is employed to define the local density from the first solvation shell. However, the model performs less well when this function acts on either longer or shorter length scales.
Collapse
Affiliation(s)
- Michael R DeLyser
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| | - W G Noid
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
57
|
Liu H, Okazaki S, Shinoda W. Heteroaryldihydropyrimidines Alter Capsid Assembly By Adjusting the Binding Affinity and Pattern of the Hepatitis B Virus Core Protein. J Chem Inf Model 2019; 59:5104-5110. [DOI: 10.1021/acs.jcim.9b01010] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Huihui Liu
- Department of Materials Chemistry, Nagoya University, Furo-cho Chikusa-ku, Nagoya 464-8603, Japan
| | - Susumu Okazaki
- Department of Materials Chemistry, Nagoya University, Furo-cho Chikusa-ku, Nagoya 464-8603, Japan
| | - Wataru Shinoda
- Department of Materials Chemistry, Nagoya University, Furo-cho Chikusa-ku, Nagoya 464-8603, Japan
| |
Collapse
|
58
|
Zgorski A, Pastor RW, Lyman E. Surface Shear Viscosity and Interleaflet Friction from Nonequilibrium Simulations of Lipid Bilayers. J Chem Theory Comput 2019; 15:6471-6481. [PMID: 31476126 DOI: 10.1021/acs.jctc.9b00683] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Nonequilibrium simulation protocols based on shear deformations are applied to determine the surface viscosity and interleaflet friction of lipid bilayers. At high shear rates, a non-Newtonian shear thinning regime is observed, but lower shear rates yield a Newtonian plateau and results that are consistent with equilibrium measurements based on fluctuation-dissipation theorems. Application to all-atom bilayers modeled with the CHARMM36 parameter set yields values for the surface viscosity that are consistent with microscopic measurements based on membrane protein diffusion but are approximately 10 times lower than more macroscopic experimental measurements. The interleaflet friction is about 10 times lower than experimental measurements. Trends across different lipids, temperatures, and ternary liquid-disordered phase mixtures produce results that are consistent with experimental diffusion constants. Application of the protocol to the liquid-ordered phase fails to yield a Newtonian plateau, suggesting more complex rheology.
Collapse
Affiliation(s)
| | - Richard W Pastor
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| | | |
Collapse
|
59
|
Papavasileiou KD, Peristeras LD, Bick A, Economou IG. Molecular Dynamics Simulation of Pure n-Alkanes and Their Mixtures at Elevated Temperatures Using Atomistic and Coarse-Grained Force Fields. J Phys Chem B 2019; 123:6229-6243. [PMID: 31251061 DOI: 10.1021/acs.jpcb.9b02840] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The properties of higher n-alkanes and their mixtures is a topic of significant interest for the oil and chemical industry. However, the experimental data at high temperatures are scarce. The present study focuses on simulating n-dodecane, n-octacosane, their binary mixture at a n-dodecane mole fraction of 0.3, and a model mixture of the commercially available hydrocarbon wax SX-70 to evaluate the performance of several force fields on the reproduction of properties such as liquid densities, surface tension, and viscosities. Molecular dynamics simulations over a broad temperature range from 323.15 to 573.15 K were employed in examining a broad set of atomistic molecular models assessed for the reproduction of experimental data. The well-established united atom TraPPE (TraPPE-UA) was compared against the all atom optimized potentials for liquid simulations (OPLS) reparametrization for long n-alkanes, L-OPLS, as well as Lipid14 and MARTINI force fields. All models qualitatively reproduce the temperature dependence of the aforementioned properties, but TraPPE-UA was found to reproduce liquid densities most accurately and consistently over the entire temperature range. TraPPE-UA and MARTINI were very successful in reproducing surface tensions, and L-OPLS was found to be the most accurate in reproducing the measured viscosities as compared to the other models. Our simulations show that these widely used force fields originating from the world of biomolecular simulations are suitable candidates in the study of n-alkane properties, both in the pure and mixture states.
Collapse
Affiliation(s)
- Konstantinos D Papavasileiou
- Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modelling of Materials Laboratory , National Center for Scientific Research "Demokritos" , Aghia Paraskevi, Attikis, GR-15310 Athens , Greece.,Scienomics SARL , 16 rue de l'Arcade , 75008 , Paris , France
| | - Loukas D Peristeras
- Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modelling of Materials Laboratory , National Center for Scientific Research "Demokritos" , Aghia Paraskevi, Attikis, GR-15310 Athens , Greece
| | - Andreas Bick
- Scienomics SARL , 16 rue de l'Arcade , 75008 , Paris , France
| | - Ioannis G Economou
- Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modelling of Materials Laboratory , National Center for Scientific Research "Demokritos" , Aghia Paraskevi, Attikis, GR-15310 Athens , Greece.,Chemical Engineering Program , Texas A&M University at Qatar , Education City , P.O. Box 23874, Doha , Qatar
| |
Collapse
|
60
|
Qiu Y, Nerenberg PS, Head-Gordon T, Wang LP. Systematic Optimization of Water Models Using Liquid/Vapor Surface Tension Data. J Phys Chem B 2019; 123:7061-7073. [DOI: 10.1021/acs.jpcb.9b05455] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yudong Qiu
- Chemistry Department, University of California, Davis, Davis, California 95616, United States
| | - Paul S. Nerenberg
- Departments of Physics & Astronomy and Biological Sciences, California State University, Los Angeles, California 90032, United States
| | - Teresa Head-Gordon
- Pitzer Theory Center and Departments of Chemistry, Bioengineering and Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720, United States
| | - Lee-Ping Wang
- Chemistry Department, University of California, Davis, Davis, California 95616, United States
| |
Collapse
|
61
|
Flood E, Boiteux C, Lev B, Vorobyov I, Allen TW. Atomistic Simulations of Membrane Ion Channel Conduction, Gating, and Modulation. Chem Rev 2019; 119:7737-7832. [DOI: 10.1021/acs.chemrev.8b00630] [Citation(s) in RCA: 65] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Emelie Flood
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Céline Boiteux
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Bogdan Lev
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Igor Vorobyov
- Department of Physiology & Membrane Biology/Department of Pharmacology, University of California, Davis, 95616, United States
| | - Toby W. Allen
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| |
Collapse
|
62
|
Krämer A, Pickard FC, Huang J, Venable RM, Simmonett AC, Reith D, Kirschner KN, Pastor RW, Brooks BR. Interactions of Water and Alkanes: Modifying Additive Force Fields to Account for Polarization Effects. J Chem Theory Comput 2019; 15:3854-3867. [PMID: 31002505 DOI: 10.1021/acs.jctc.9b00016] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Atomistic biomolecular simulations predominantly utilize additive force fields (FF), where the electrostatic potential is modeled by fixed point charges. Among other consequences, the lack of polarizability in these models undermines the balance of hydrophilic/hydrophobic nonbonded interactions. Simulations of water/alkane systems using the TIP3P water model and CHARMM36 parameters reveal a 1 kcal/mol overestimate of the experimental transfer free energy of water to hexadecane; more recent optimized water models (SPC/E, TIP4P/2005, TIP4P-Ew, TIP3P-FB, TIP4P-FB, OPC, TIP4P-D) overestimate this transfer free energy by approximately 2 kcal/mol. In contrast, the polarizable SWM4-NDP and SWM6 water models reproduce experimental values to within statistical error. As an alternative to explicitly modeling polarizability, this paper develops an efficient automated workflow to optimize pair-specific Lennard-Jones parameters within an additive FF. Water/hexadecane is used as a prototype and the free energy of water transfer to hexadecane as a target. The optimized model yields quantitative agreement with the experimental transfer free energy and improves the water/hexadecane interfacial tension by 6%. Simulations of five different lipid bilayers show a strong increase of water permeabilities compared to the unmodified CHARMM36 lipid FF which consistently improves match with experiment: the order-of-magnitude underestimate for monounsaturated bilayers is rectified and the factor of 2.8-4 underestimate for saturated bilayers is turned into a factor of 1.5-3 overestimate. While agreement with experiment is decreased for the diffusion constant of water in hexadecane, alkane transfer free energies, and the bilayers' area per lipid, the method provides a permeant-specific route to achieve a wide range of heterogeneous observables via rapidly optimized pairwise parameters.
Collapse
Affiliation(s)
- Andreas Krämer
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States.,Institute of Technology, Resource and Energy-Efficient Engineering , Bonn-Rhein-Sieg University of Applied Sciences , Grantham-Allee 20 , 53757 Sankt Augustin , Germany
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| | - Jing Huang
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States.,School of Life Sciences , Westlake University , 18 Shilongshan Road , Hangzhou 310024 , Zhejiang China.,Department of Pharmaceutical Science, School of Pharmacy , University of Maryland , 20 Penn Street , Baltimore , Maryland 21201 , United States
| | - Richard M Venable
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| | - Dirk Reith
- Institute of Technology, Resource and Energy-Efficient Engineering , Bonn-Rhein-Sieg University of Applied Sciences , Grantham-Allee 20 , 53757 Sankt Augustin , Germany
| | - Karl N Kirschner
- Institute of Technology, Resource and Energy-Efficient Engineering , Bonn-Rhein-Sieg University of Applied Sciences , Grantham-Allee 20 , 53757 Sankt Augustin , Germany
| | - Richard W Pastor
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| |
Collapse
|
63
|
Abstract
This Review illustrates the evaluation of permeability of lipid membranes from molecular dynamics (MD) simulation primarily using water and oxygen as examples. Membrane entrance, translocation, and exit of these simple permeants (one hydrophilic and one hydrophobic) can be simulated by conventional MD, and permeabilities can be evaluated directly by Fick's First Law, transition rates, and a global Bayesian analysis of the inhomogeneous solubility-diffusion model. The assorted results, many of which are applicable to simulations of nonbiological membranes, highlight the limitations of the homogeneous solubility diffusion model; support the utility of inhomogeneous solubility diffusion and compartmental models; underscore the need for comparison with experiment for both simple solvent systems (such as water/hexadecane) and well-characterized membranes; and demonstrate the need for microsecond simulations for even simple permeants like water and oxygen. Undulations, subdiffusion, fractional viscosity dependence, periodic boundary conditions, and recent developments in the field are also discussed. Last, while enhanced sampling methods and increasingly sophisticated treatments of diffusion add substantially to the repertoire of simulation-based approaches, they do not address directly the critical need for force fields with polarizability and multipoles, and constant pH methods.
Collapse
Affiliation(s)
- Richard M Venable
- Laboratory of Computational Biology, National Lung, Heart, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| | - Andreas Krämer
- Laboratory of Computational Biology, National Lung, Heart, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| | - Richard W Pastor
- Laboratory of Computational Biology, National Lung, Heart, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| |
Collapse
|
64
|
Schieber NP, Shirts MR. Configurational mapping significantly increases the efficiency of solid-solid phase coexistence calculations via molecular dynamics: Determining the FCC-HCP coexistence line of Lennard-Jones particles. J Chem Phys 2019; 150:164112. [DOI: 10.1063/1.5080431] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Natalie P. Schieber
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Michael R. Shirts
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, Colorado 80309, USA
| |
Collapse
|
65
|
Singh S. Dynamics of heroin molecule inside the lipid membrane: a molecular dynamics study. J Mol Model 2019; 25:121. [PMID: 31020452 DOI: 10.1007/s00894-019-4002-y] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2018] [Accepted: 03/21/2019] [Indexed: 11/29/2022]
Abstract
Heroin, or diamorphine (C21H23NO5), is an opium product used for various pharmaceutical and euphoric purposes. In this work, the molecular dynamics simulation study of the heroin inside the two lipid bilayers, dipalmitoylphosphatidylcholine (DMPC) and dipalmitoylphosphatidylcholine (DPPC) are presented. The whole study was conducted at three different temperatures. The location of the heroin drug, the nature of the diffusion, rotational correlation function and structural variation inside both lipid bilayers is studied. Moreover, the free energy of the solvation of the drug inside both lipid bilayers is calculated. It is found that during the whole molecular dynamics study, the drug locates at the center of both lipid membranes. The effect of the temperature is not seen at the drug location. The nature of the diffusion of the heroin drug is anomalous. The radius of gyration is calculated to study the structural variations of the heroin molecule inside both lipid bilayers. It is found that the heroin molecule does not change its structure at three temperatures. From the rotational correlation function, it is seen that the drug is more hindered for rotation inside the DPPC lipid bilayer as compared to the DMPC lipid bilayer. It is applicable for all three temperatures. The rotational correlation time of the drug is decreased while the temperature of the system is increased. In the case of DMPC, there is an abrupt change in rotational correlation time while the phase is changed.
Collapse
Affiliation(s)
- Satnam Singh
- Department of Physical Sciences, Indian Institute of Science Education & Research (IISER) Mohali, Sector 81 SAS Nagar, Manauli PO, 140306, Punjab, India.
| |
Collapse
|
66
|
Leonard AN, Wang E, Monje-Galvan V, Klauda JB. Developing and Testing of Lipid Force Fields with Applications to Modeling Cellular Membranes. Chem Rev 2019; 119:6227-6269. [DOI: 10.1021/acs.chemrev.8b00384] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
67
|
Meral D, Provasi D, Filizola M. An efficient strategy to estimate thermodynamics and kinetics of G protein-coupled receptor activation using metadynamics and maximum caliber. J Chem Phys 2019; 149:224101. [PMID: 30553249 PMCID: PMC6291190 DOI: 10.1063/1.5060960] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
Computational strategies aimed at unveiling the thermodynamic and kinetic properties of G
Protein-Coupled Receptor (GPCR) activation require extensive molecular dynamics
simulations of the receptor embedded in an explicit lipid-water environment. A possible
method for efficiently sampling the conformational space of such a complex system is
metadynamics (MetaD) with path collective variables (CVs). Here, we applied well-tempered
MetaD with path CVs to one of the few GPCRs for which both inactive and fully active
experimental structures are available, the μ-opioid receptor (MOR), and assessed the
ability of this enhanced sampling method to estimate the thermodynamic properties of
receptor activation in line with those obtained by more computationally expensive adaptive
sampling protocols. While n-body information theory analysis of these simulations
confirmed that MetaD can efficiently characterize ligand-induced allosteric communication
across the receptor, standard MetaD cannot be used directly to derive kinetic rates
because transitions are accelerated by a bias potential. Applying the principle of Maximum
Caliber (MaxCal) to the free-energy landscape of morphine-bound MOR reconstructed from
MetaD, we obtained Markov state models that yield kinetic rates of MOR activation in
agreement with those obtained by adaptive sampling. Taken together, these results suggest
that the MetaD-MaxCal combination creates an efficient strategy for estimating the
thermodynamic and kinetic properties of GPCR activation at an affordable computational
cost.
Collapse
Affiliation(s)
- Derya Meral
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, One Gustave L. Levy Place, New York, New York 10029, USA
| | - Davide Provasi
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, One Gustave L. Levy Place, New York, New York 10029, USA
| | - Marta Filizola
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, One Gustave L. Levy Place, New York, New York 10029, USA
| |
Collapse
|
68
|
Abstract
All-atom, classical force fields for protein molecular dynamics (MD) simulations currently occupy a sweet spot in the universe of computational models, sufficiently detailed to be of predictive value in many cases, yet also simple enough that some biologically relevant time scales (microseconds or more) can now be sampled via specialized hardware or enhanced sampling methods. However, due to their long evolutionary history, there is now a myriad of force field branches in current use, which can make it hard for those entering the simulation field to know which would be the best set of parameters for a given application. In this chapter, I try to give an overview of the historical motivation for the different force fields available, suggestions for how to determine the most appropriate model and what to do if the results are in conflict with experimental evidence.
Collapse
Affiliation(s)
- Robert B Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD, USA.
| |
Collapse
|
69
|
Abstract
Molecular dynamics (MD) simulations have been widely applied to computer-aided drug design (CADD). While MD has been used in a variety of applications such as free energy perturbation and long-time simulations, the accuracy of the results from those methods depends strongly on the force field used. Force fields for small molecules are crucial, as they not only serve as building blocks for developing force fields for larger biomolecules but also act as model compounds that will be transferred to ligands used in CADD. Currently, a wide range of small molecule force fields based on additive or nonpolarizable models have been developed. While these nonpolarizable force fields can produce reasonable estimations of physical properties and have shown success in a variety of systems, there is still room for improvements due to inherent limitations in these models including the lack of an electronic polarization response. For this reason, incorporating polarization effects into the energy function underlying a force field is believed to be an important step forward, giving rise to the development of polarizable force fields. Recent simulations of biological systems have indicated that polarizable force fields are able to provide a better physical representation of intermolecular interactions and, in many cases, better agreement with experimental properties than nonpolarizable, additive force fields. Therefore, this chapter focuses on the development of small molecule force fields with emphasis on polarizable models. It begins with a brief introduction on the importance of small molecule force fields and their evolution from additive to polarizable force fields. Emphasis is placed on the additive CHARMM General Force Field and the polarizable force field based on the classical Drude oscillator. The theory for the Drude polarizable force field and results for small molecules are presented showing their improvements over the additive model. The potential importance of polarization for their application in a wide range of biological systems including CADD is then discussed.
Collapse
Affiliation(s)
- Fang-Yu Lin
- Department of Pharmaceutical Sciences, Computer-Aided Drug Design Center, School of Pharmacy, University of Maryland, Baltimore, MD, USA
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, Computer-Aided Drug Design Center, School of Pharmacy, University of Maryland, Baltimore, MD, USA.
| |
Collapse
|
70
|
Nhan C, Rix CJ, May BK, Hung A. Temperature-induced structural changes of apo-lactoferrin and their functional implications: a molecular dynamics simulation study. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1562187] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Carol Nhan
- School of Science, RMIT University, Bundoora Campus, Melbourne, Australia
| | - Colin J. Rix
- School of Science, RMIT University, City Campus, Melbourne, Australia
| | - Bee K. May
- School of Science, RMIT University, Bundoora Campus, Melbourne, Australia
| | - Andrew Hung
- School of Science, RMIT University, City Campus, Melbourne, Australia
| |
Collapse
|
71
|
Nhan C, Small DM, May BK, Hung A. pH-induced structural changes of apo-lactoferrin and implications for its function: a molecular dynamics simulation study. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1535178] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Carol Nhan
- School of Science, RMIT University, Melbourne, Australia
| | | | - Bee K. May
- School of Science, RMIT University, Melbourne, Australia
| | - Andrew Hung
- School of Science, RMIT University, Melbourne, Australia
| |
Collapse
|
72
|
Lawrence RL, Hughes ZE, Cendan VJ, Liu Y, Lim CK, Prasad PN, Swihart MT, Walsh TR, Knecht MR. Optical Control of Nanoparticle Catalysis Influenced by Photoswitch Positioning in Hybrid Peptide Capping Ligands. ACS APPLIED MATERIALS & INTERFACES 2018; 10:33640-33651. [PMID: 30185023 DOI: 10.1021/acsami.8b10582] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Here, we present an in-depth analysis of structural factors that modulate peptide-capped nanoparticle catalytic activity via optically driven structural reconfiguration of the biointerface present at the particle surface. Six different sets of peptide-capped Au nanoparticles were prepared, in which an azobenzene photoswitch was incorporated into one of two well-studied peptide sequences with known affinity for Au, each at one of three different positions: the N- or C-terminus or mid-sequence. Changes in the photoswitch isomerization state induce a reversible structural change in the surface-bound peptide, which modulates the catalytic activity of the material. This control of reactivity is attributed to changes in the amount of accessible metallic surface area available to drive the reaction. This research specifically focuses on the effect of the peptide sequence and photoswitch position in the biomolecule, from which potential target systems for on/off reactivity have been identified. Additionally, trends associated with photoswitch position for a peptide sequence (Pd4) have been identified. Integrating the azobenzene at the N-terminus or central region results in nanocatalysts with greater reactivity in the trans and cis conformations, respectively, however, positioning the photoswitch at the C-terminus gives rise to a unique system that is reactive in the trans conformation and partially deactivated in the cis conformation. These results provide a fundamental basis for new directions in nanoparticle catalyst development to control activity in real time, which could have significant implications in the design of catalysts for multistep reactions using a single catalyst. Additionally, such a fine level of interfacial structural control could prove to be important for applications beyond catalysis, including biosensing, photonics, and energy technologies that are highly dependent on particle surface structures.
Collapse
Affiliation(s)
- Randy L Lawrence
- Department of Chemistry , University of Miami , 1301 Memorial Drive , Coral Gables , Florida 33146 , United States
| | - Zak E Hughes
- Institute for Frontier Materials , Deakin University , Geelong , Victoria 3216 , Australia
| | - Vincent J Cendan
- Department of Chemistry , University of Miami , 1301 Memorial Drive , Coral Gables , Florida 33146 , United States
| | | | | | | | | | - Tiffany R Walsh
- Institute for Frontier Materials , Deakin University , Geelong , Victoria 3216 , Australia
| | - Marc R Knecht
- Department of Chemistry , University of Miami , 1301 Memorial Drive , Coral Gables , Florida 33146 , United States
| |
Collapse
|
73
|
Gong Z, Wu Y, Wu L, Sun H. Predicting Thermodynamic Properties of Alkanes by High-Throughput Force Field Simulation and Machine Learning. J Chem Inf Model 2018; 58:2502-2516. [DOI: 10.1021/acs.jcim.8b00407] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Zheng Gong
- School of Chemistry and Chemical Engineering, Materials Genome Initiative Center, and Key Laboratory of Scientific and Engineering Computing of Ministry of Education, Shanghai Jiao Tong University, Shanghai, China 200240
| | - Yanze Wu
- School of Chemistry and Chemical Engineering, Materials Genome Initiative Center, and Key Laboratory of Scientific and Engineering Computing of Ministry of Education, Shanghai Jiao Tong University, Shanghai, China 200240
| | - Liang Wu
- School of Chemistry and Chemical Engineering, Materials Genome Initiative Center, and Key Laboratory of Scientific and Engineering Computing of Ministry of Education, Shanghai Jiao Tong University, Shanghai, China 200240
| | - Huai Sun
- School of Chemistry and Chemical Engineering, Materials Genome Initiative Center, and Key Laboratory of Scientific and Engineering Computing of Ministry of Education, Shanghai Jiao Tong University, Shanghai, China 200240
| |
Collapse
|
74
|
Rackers JA, Liu C, Ren P, Ponder JW. A physically grounded damped dispersion model with particle mesh Ewald summation. J Chem Phys 2018; 149:084115. [PMID: 30193468 DOI: 10.1063/1.5030434] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Accurate modeling of dispersion is critical to the goal of predictive biomolecular simulations. To achieve this accuracy, a model must be able to correctly capture both the short-range and asymptotic behavior of dispersion interactions. We present here a damped dispersion model based on the overlap of charge densities that correctly captures both regimes. The overlap damped dispersion model represents a classical physical interpretation of dispersion: the interaction between the instantaneous induced dipoles of two distinct charge distributions. This model is shown to be an excellent fit with symmetry adapted perturbation theory dispersion energy calculations, yielding an RMS error on the S101x7 database of 0.5 kcal/mol. Moreover, the damping function used in this model is wholly derived and parameterized from the electrostatic dipole-dipole interaction, making it not only physically grounded but transferable as well.
Collapse
Affiliation(s)
- Joshua A Rackers
- Program in Computational and Molecular Biophysics, Washington University School of Medicine, Saint Louis, Missouri 63110, USA
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Jay W Ponder
- Program in Computational and Molecular Biophysics, Washington University School of Medicine, Saint Louis, Missouri 63110, USA
| |
Collapse
|
75
|
Comesaña A, Pérez-Rodríguez M, Fernández-Fernández AM, Piñeiro MM. A description of hydroquinone clathrates using molecular dynamics: Molecular model and crystalline structures for CH4and CO2guests. J Chem Phys 2018; 148:244502. [PMID: 29960310 DOI: 10.1063/1.5027807] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- A. Comesaña
- Universidade de Vigo, Departamento de Física Aplicada, Campus Lagoas-Marcosende, E36310 Vigo, Spain
| | - M. Pérez-Rodríguez
- Universidade de Vigo, Departamento de Física Aplicada, Campus Lagoas-Marcosende, E36310 Vigo, Spain
| | | | - M. M. Piñeiro
- Universidade de Vigo, Departamento de Física Aplicada, Campus Lagoas-Marcosende, E36310 Vigo, Spain
| |
Collapse
|
76
|
Leonard AN, Pastor RW, Klauda JB. Parameterization of the CHARMM All-Atom Force Field for Ether Lipids and Model Linear Ethers. J Phys Chem B 2018; 122:6744-6754. [PMID: 29870257 DOI: 10.1021/acs.jpcb.8b02743] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Linear ethers such as polyethylene glycol have extensive industrial and medical applications. Additionally, phospholipids containing an ether linkage between the glycerol backbone and hydrophobic tails are prevalent in human red blood cells and nerve tissue. This study uses ab initio results to revise the CHARMM additive (C36) partial-charge and dihedral parameters for linear ethers and develop parameters for the ether-linked phospholipid 1,2-di- O-hexadecyl- sn-glycero-3-phosphocholine (DHPC). The new force field, called C36e, more accurately represents the dihedral potential energy landscape and improves the densities and free energies of hydration of linear ethers. C36e allows more water to penetrate into a DHPC bilayer, increasing the surface area per lipid compared to simulations carried out with the original C36 ether parameters and improving the overall structural properties obtained from X-ray and neutron scattering. Comparison with an ester-linked DPPC bilayer (1,2-dipalmitoyl- sn-phosphatidylcholine) reveals that the ether linkage increases water organization in the headgroup region. This effect is a likely explanation for the experimentally lower water permeability of bilayers composed of ether-linked lipids.
Collapse
Affiliation(s)
| | - Richard W Pastor
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| | | |
Collapse
|
77
|
Cox SJ, Geissler PL. Interfacial ion solvation: Obtaining the thermodynamic limit from molecular simulations. J Chem Phys 2018; 148:222823. [DOI: 10.1063/1.5020563] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Stephen J. Cox
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Phillip L. Geissler
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
78
|
Huang J, Simmonett AC, Pickard FC, MacKerell AD, Brooks BR. Mapping the Drude polarizable force field onto a multipole and induced dipole model. J Chem Phys 2018; 147:161702. [PMID: 29096511 DOI: 10.1063/1.4984113] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
The induced dipole and the classical Drude oscillator represent two major approaches for the explicit inclusion of electronic polarizability into force field-based molecular modeling and simulations. In this work, we explore the equivalency of these two models by comparing condensed phase properties computed using the Drude force field and a multipole and induced dipole (MPID) model. Presented is an approach to map the electrostatic model optimized in the context of the Drude force field onto the MPID model. Condensed phase simulations on water and 15 small model compounds show that without any reparametrization, the MPID model yields properties similar to the Drude force field with both models yielding satisfactory reproduction of a range of experimental values and quantum mechanical data. Our results illustrate that the Drude oscillator model and the point induced dipole model are different representations of essentially the same physical model. However, results indicate the presence of small differences between the use of atomic multipoles and off-center charge sites. Additionally, results on the use of dispersion particle mesh Ewald further support its utility for treating long-range Lennard Jones dispersion contributions in the context of polarizable force fields. The main motivation in demonstrating the transferability of parameters between the Drude and MPID models is that the more than 15 years of development of the Drude polarizable force field can now be used with MPID formalism without the need for dual-thermostat integrators nor self-consistent iterations. This opens up a wide range of new methodological opportunities for polarizable models.
Collapse
Affiliation(s)
- Jing Huang
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, 20 Penn St., Baltimore, Maryland 21201, USA
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, 5635 Fishers Lane, Rockville, Maryland 20852, USA
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, 5635 Fishers Lane, Rockville, Maryland 20852, USA
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, 20 Penn St., Baltimore, Maryland 21201, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, 5635 Fishers Lane, Rockville, Maryland 20852, USA
| |
Collapse
|
79
|
|
80
|
Yong T, Chen S, Xie Y, Chen D, Su J, Shuai O, Jiao C, Zuo D. Cordycepin, a Characteristic Bioactive Constituent in Cordyceps militaris, Ameliorates Hyperuricemia through URAT1 in Hyperuricemic Mice. Front Microbiol 2018; 9:58. [PMID: 29422889 PMCID: PMC5788910 DOI: 10.3389/fmicb.2018.00058] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2017] [Accepted: 01/10/2018] [Indexed: 12/22/2022] Open
Abstract
Recently, we've reported the anti-hyperuricemic effects of Cordyceps militaris. As a characteristic compound of C. militaris, we hypothesized that cordycepin may play a role in preventing hyperurecimia. Remarkably, cordycepin produced important anti-hyperuricemic actions, decreasing SUA (serum uric acid) to 216, 210, and 203 μmol/L (P < 0.01) at 15, 30, and 60 mg/kg in comparison of hyperuricemic control (337 μmol/L), closing to normal control (202 μmol/L). Elisa, RT-PCR and western blot analysis demonstrated that the actions may be attributed to its downregulation of uric acid transporter 1 (URAT1) in kidney. Serum creatinine levels and blood urine nitrogen and liver, kidney, and spleen coefficients demonstrated that cordycepin may not impact liver, renal, and spleen functions. In addition, we used computational molecular simulation to investigate the binding mechanism of cordycepin. Of which, van der Waals interaction dominated the binding. Residues TRP290, ARG17, ALA408, GLY411, and MET147 contributed mainly on nonpolar energy. This provided the theoretical guidance to rationally design and synthesis novel URAT1 inhibitors.
Collapse
Affiliation(s)
- Tianqiao Yong
- State Key Laboratory of Applied Microbiology Southern China, Guangdong Provincial Key Laboratory of Microbial Culture Collection and Application and Guangdong Open Laboratory of Applied Microbiology, Guangdong Institute of Microbiology, Guangzhou, China.,Guangdong Yuewei Edible Fungi Technology Co, Guangzhou, China
| | - Shaodan Chen
- State Key Laboratory of Applied Microbiology Southern China, Guangdong Provincial Key Laboratory of Microbial Culture Collection and Application and Guangdong Open Laboratory of Applied Microbiology, Guangdong Institute of Microbiology, Guangzhou, China.,Guangdong Yuewei Edible Fungi Technology Co, Guangzhou, China
| | - Yizhen Xie
- State Key Laboratory of Applied Microbiology Southern China, Guangdong Provincial Key Laboratory of Microbial Culture Collection and Application and Guangdong Open Laboratory of Applied Microbiology, Guangdong Institute of Microbiology, Guangzhou, China.,Guangdong Yuewei Edible Fungi Technology Co, Guangzhou, China
| | - Diling Chen
- State Key Laboratory of Applied Microbiology Southern China, Guangdong Provincial Key Laboratory of Microbial Culture Collection and Application and Guangdong Open Laboratory of Applied Microbiology, Guangdong Institute of Microbiology, Guangzhou, China
| | - Jiyan Su
- State Key Laboratory of Applied Microbiology Southern China, Guangdong Provincial Key Laboratory of Microbial Culture Collection and Application and Guangdong Open Laboratory of Applied Microbiology, Guangdong Institute of Microbiology, Guangzhou, China
| | - Ou Shuai
- State Key Laboratory of Applied Microbiology Southern China, Guangdong Provincial Key Laboratory of Microbial Culture Collection and Application and Guangdong Open Laboratory of Applied Microbiology, Guangdong Institute of Microbiology, Guangzhou, China
| | - Chunwei Jiao
- State Key Laboratory of Applied Microbiology Southern China, Guangdong Provincial Key Laboratory of Microbial Culture Collection and Application and Guangdong Open Laboratory of Applied Microbiology, Guangdong Institute of Microbiology, Guangzhou, China
| | - Dan Zuo
- Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou, China
| |
Collapse
|
81
|
Leonard AN, Simmonett AC, Pickard FC, Huang J, Venable RM, Klauda JB, Brooks BR, Pastor RW. Comparison of Additive and Polarizable Models with Explicit Treatment of Long-Range Lennard-Jones Interactions Using Alkane Simulations. J Chem Theory Comput 2018; 14:948-958. [PMID: 29268012 DOI: 10.1021/acs.jctc.7b00948] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Long-range Lennard-Jones (LJ) interactions have a significant impact on the structural and thermodynamic properties of nonpolar systems. While several methods have been introduced for the treatment of long-range LJ interactions in molecular dynamics (MD) simulations, increased accuracy and extended applicability is required for anisotropic systems such as lipid bilayers. The recently refined Lennard-Jones particle-mesh Ewald (LJ-PME) method extends the particle-mesh Ewald (PME) method to long-range LJ interactions and is suitable for use with anisotropic systems. Implementation of LJ-PME with the CHARMM36 (C36) additive and CHARMM Drude polarizable force fields improves agreement with experiment for density, isothermal compressibility, surface tension, viscosity, translational diffusion, and 13C T1 relaxation times of pure alkanes. Trends in the temperature dependence of the density and isothermal compressibility of hexadecane are also improved. While the C36 additive force field with LJ-PME remains a useful model for liquid alkanes, the Drude polarizable force field with LJ-PME is more accurate for nearly all quantities considered. LJ-PME is also preferable to the isotropic long-range correction for hexadecane because the molecular order extends to nearly 20 Å, well beyond the usual 10-12 Å cutoffs used in most simulations.
Collapse
Affiliation(s)
- Alison N Leonard
- Biophysics Program, University of Maryland , College Park, Maryland 20742, United States
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States
| | - Jing Huang
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States.,Department of Pharmaceutical Science, School of Pharmacy, University of Maryland , 20 Penn Street, Baltimore, Maryland 21201, United States
| | - Richard M Venable
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States
| | - Jeffery B Klauda
- Biophysics Program, University of Maryland , College Park, Maryland 20742, United States.,Department of Chemical and Biomolecular Engineering, University of Maryland , College Park, Maryland 20742, United States
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States
| | - Richard W Pastor
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States
| |
Collapse
|
82
|
Aldeghi M, Bluck JP, Biggin PC. Absolute Alchemical Free Energy Calculations for Ligand Binding: A Beginner's Guide. Methods Mol Biol 2018; 1762:199-232. [PMID: 29594774 DOI: 10.1007/978-1-4939-7756-7_11] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Many thermodynamic quantities can be extracted from computer simulations that generate an ensemble of microstates according to the principles of statistical mechanics. Among these quantities is the free energy of binding of a small molecule to a macromolecule, such as a protein. Here, we present an introductory overview of a protocol that allows for the estimation of ligand binding free energies via molecular dynamics simulations. While we focus on the binding of organic molecules to proteins, the approach is in principle transferable to any pair of molecules.
Collapse
Affiliation(s)
- Matteo Aldeghi
- Structural Bioinformatics and Computational Biochemistry, Department of Biochemistry, University of Oxford, Oxford, UK
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Joseph P Bluck
- Structural Bioinformatics and Computational Biochemistry, Department of Biochemistry, University of Oxford, Oxford, UK
| | - Philip C Biggin
- Structural Bioinformatics and Computational Biochemistry, Department of Biochemistry, University of Oxford, Oxford, UK.
| |
Collapse
|
83
|
Klink GV, Golovin AV, Bazykin GA. Substitutions into amino acids that are pathogenic in human mitochondrial proteins are more frequent in lineages closely related to human than in distant lineages. PeerJ 2017; 5:e4143. [PMID: 29250469 PMCID: PMC5731343 DOI: 10.7717/peerj.4143] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2017] [Accepted: 11/16/2017] [Indexed: 11/23/2022] Open
Abstract
Propensities for different amino acids within a protein site change in the course of evolution, so that an amino acid deleterious in a particular species may be acceptable at the same site in a different species. Here, we study the amino acid-changing variants in human mitochondrial genes, and analyze their occurrence in non-human species. We show that substitutions giving rise to such variants tend to occur in lineages closely related to human more frequently than in more distantly related lineages, indicating that a human variant is more likely to be deleterious in more distant species. Unexpectedly, substitutions giving rise to amino acids that correspond to alleles pathogenic in humans also more frequently occur in more closely related lineages. Therefore, a pathogenic variant still tends to be more acceptable in human mitochondria than a variant that may only be fit after a substantial perturbation of the protein structure.
Collapse
Affiliation(s)
- Galya V. Klink
- Sector of Molecular Evolution, Institute for Information Transmission Problems (Kharkevich Institute) of the Russian Academy of Sciences, Moscow, Russian Federation
| | - Andrey V. Golovin
- Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, Moscow, Russian Federation
| | - Georgii A. Bazykin
- Sector of Molecular Evolution, Institute for Information Transmission Problems (Kharkevich Institute) of the Russian Academy of Sciences, Moscow, Russian Federation
- Center for Data-Intensive Biomedicine and Biotechnology, Skolkovo Institute of Science and Technology, Skolkovo, Russian Federation
| |
Collapse
|
84
|
Mohebifar M, Johnson ER, Rowley CN. Evaluating Force-Field London Dispersion Coefficients Using the Exchange-Hole Dipole Moment Model. J Chem Theory Comput 2017; 13:6146-6157. [DOI: 10.1021/acs.jctc.7b00522] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Mohamad Mohebifar
- Department
of Chemistry, Memorial University of Newfoundland, St. John’s, Newfoundland
and Labrador, Canada
| | - Erin R. Johnson
- Department
of Chemistry, Dalhousie University, Halifax, Nova Scotia, Canada
| | - Christopher N. Rowley
- Department
of Chemistry, Memorial University of Newfoundland, St. John’s, Newfoundland
and Labrador, Canada
| |
Collapse
|
85
|
Sega M, Dellago C. Long-Range Dispersion Effects on the Water/Vapor Interface Simulated Using the Most Common Models. J Phys Chem B 2017; 121:3798-3803. [DOI: 10.1021/acs.jpcb.6b12437] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Marcello Sega
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, A-1090, Vienna, Austria
| | - Christoph Dellago
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, A-1090, Vienna, Austria
| |
Collapse
|
86
|
Cerar J, Lajovic A, Jamnik A, Tomšič M. Performance of various models in structural characterization of n -butanol: Molecular dynamics and X-ray scattering studies. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2016.12.057] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
87
|
Javanainen M, Vattulainen I, Monticelli L. On Atomistic Models for Molecular Oxygen. J Phys Chem B 2017; 121:518-528. [PMID: 28004930 DOI: 10.1021/acs.jpcb.6b11183] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Molecular oxygen (O2) is key to all life on earth, as it is constantly cycled via photosynthesis and cellular respiration. Substantial scientific effort has been devoted to understanding every part of this cycle. Classical molecular dynamics (MD) simulations have been used to study some of the key processes involved in cellular respiration: O2 permeation through alveolar monolayers and cellular membranes, its binding to hemoglobin during transport in the bloodstream, as well as its transport along optimal pathways toward its reduction sites in proteins. Moreover, MD simulations can help interpret the results of several imaging techniques in which O2 is used because of its paramagnetic nature. However, despite the widespread use of computational models for the O2 molecule, their performances have never been systematically evaluated. In this paper, we assess the performances of 14 different models of O2 available in the literature by calculating four thermodynamic properties: density, heat of vaporization, free energy of hydration, and free energy of solvation in hexadecane. For each property, reliable experimental data are available. Most models perform reasonably well in predicting the correct trends, but they fail to reproduce the experimental data quantitatively. We then develop new models for O2, with and without a quadrupole moment, and compare their behavior with the behavior of previously published models. The new models show significant improvement in terms of density, heat of vaporization, and free energy of hydration. However, quantitative agreement with water-oil partitioning is not reached due to discrepancies between the calculated and measured free energies of solvation in hexadecane. We suggest that classical pairwise-additive models may be inadequate to properly describe the thermodynamics of solvation of apolar species, such as O2, in apolar solvents.
Collapse
Affiliation(s)
- Matti Javanainen
- Department of Physics, Tampere University of Technology , 33720 Tampere, Finland.,Department of Physics, University of Helsinki , 00100 Helsinki, Finland
| | - Ilpo Vattulainen
- Department of Physics, Tampere University of Technology , 33720 Tampere, Finland.,Department of Physics, University of Helsinki , 00100 Helsinki, Finland.,MEMPHYS - Centre for Biomembrane Physics, University of Southern Denmark , 5230 Odense, Denmark
| | - Luca Monticelli
- University of Lyon, CNRS, UMR 5086 MMSB , 69367 Lyon, France
| |
Collapse
|
88
|
Shirts MR, Klein C, Swails JM, Yin J, Gilson MK, Mobley DL, Case DA, Zhong ED. Lessons learned from comparing molecular dynamics engines on the SAMPL5 dataset. J Comput Aided Mol Des 2017; 31:147-161. [PMID: 27787702 PMCID: PMC5581938 DOI: 10.1007/s10822-016-9977-1] [Citation(s) in RCA: 134] [Impact Index Per Article: 19.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2016] [Accepted: 09/20/2016] [Indexed: 11/30/2022]
Abstract
We describe our efforts to prepare common starting structures and models for the SAMPL5 blind prediction challenge. We generated the starting input files and single configuration potential energies for the host-guest in the SAMPL5 blind prediction challenge for the GROMACS, AMBER, LAMMPS, DESMOND and CHARMM molecular simulation programs. All conversions were fully automated from the originally prepared AMBER input files using a combination of the ParmEd and InterMol conversion programs. We find that the energy calculations for all molecular dynamics engines for this molecular set agree to better than 0.1 % relative absolute energy for all energy components, and in most cases an order of magnitude better, when reasonable choices are made for different cutoff parameters. However, there are some surprising sources of statistically significant differences. Most importantly, different choices of Coulomb's constant between programs are one of the largest sources of discrepancies in energies. We discuss the measures required to get good agreement in the energies for equivalent starting configurations between the simulation programs, and the energy differences that occur when simulations are run with program-specific default simulation parameter values. Finally, we discuss what was required to automate this conversion and comparison.
Collapse
Affiliation(s)
- Michael R Shirts
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, CO, USA.
| | - Christoph Klein
- Department of Chemical Engineering, Vanderbilt University, Nashville, TN, USA
| | - Jason M Swails
- Department of Chemistry and Chemical Biology, Rutgers University, New Brunswick, NJ, USA
| | - Jian Yin
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, CA, USA
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, CA, USA
| | - David L Mobley
- Departments of Pharmaceutical Sciences and Chemistry, University of California, Irvine, Irvine, CA, USA
| | - David A Case
- Department of Chemistry and Chemical Biology, Rutgers University, New Brunswick, NJ, USA
| | - Ellen D Zhong
- Department of Chemical Engineering, University of Virginia, Charlottesville, VA, USA
| |
Collapse
|
89
|
Werth S, Stöbener K, Horsch M, Hasse H. Simultaneous description of bulk and interfacial properties of fluids by the Mie potential. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1206218] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Stephan Werth
- Department of Mechanical and Process Engineering, Laboratory of Engineering Thermodynamics, University of Kaiserslautern, Kaiserslautern, Germany
| | - Katrin Stöbener
- Department for Optimization, Fraunhofer Institute for Industrial Mathematics, Kaiserslautern, Germany
| | - Martin Horsch
- Department of Mechanical and Process Engineering, Laboratory of Engineering Thermodynamics, University of Kaiserslautern, Kaiserslautern, Germany
| | - Hans Hasse
- Department of Mechanical and Process Engineering, Laboratory of Engineering Thermodynamics, University of Kaiserslautern, Kaiserslautern, Germany
| |
Collapse
|
90
|
Horta BAC, Merz PT, Fuchs PFJ, Dolenc J, Riniker S, Hünenberger PH. A GROMOS-Compatible Force Field for Small Organic Molecules in the Condensed Phase: The 2016H66 Parameter Set. J Chem Theory Comput 2016; 12:3825-50. [DOI: 10.1021/acs.jctc.6b00187] [Citation(s) in RCA: 88] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Bruno A. C. Horta
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Pascal T. Merz
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Patrick F. J. Fuchs
- Institut Jacques Monod, UMR 7592 CNRS, Université Paris-Diderot, Sorbonne Paris Cité, F-75205 Paris, France
| | - Jozica Dolenc
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
- Chemistry,
Biology and Pharmacy Information Center, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | | |
Collapse
|