1
|
Hosseni A, Ashbaugh HS. Osmotic Force Balance Evaluation of Aqueous Electrolyte Osmotic Pressures and Chemical Potentials. J Chem Theory Comput 2023; 19:8826-8838. [PMID: 37978934 PMCID: PMC10720338 DOI: 10.1021/acs.jctc.3c00982] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2023] [Revised: 10/26/2023] [Accepted: 10/27/2023] [Indexed: 11/19/2023]
Abstract
Concentrated aqueous salt solutions are ubiquitous in problems of biological and environmental relevance. The development of accurate force fields that capture the interactions between dissolved species in solution is crucial to simulating these systems to gain molecular insights into the underlying processes under saline conditions. The osmotic pressure is a relatively simple thermodynamic property connecting the experimental and simulation measurements of the associative properties of the ions in solution. Milner [C. Gillespie and S. T. Milner, Soft Matter, 16, 9816 (2020)] proposed a simulation approach to evaluate the osmotic pressures of salts in solution by applying a restraint potential to the ions alone in solution and determining the resulting pressure required to balance that potential, referred to here as the osmotic force balance. Here, we expand Milner's approach, demonstrating that the chemical potentials of the salts in solution as a function of concentration can be fitted to the concentration profiles determined from simulation, additionally providing an analytical expression for the osmotic pressure. This approach is used to determine the osmotic pressures of 15 alkali halide salts in water from simulations. The cross interactions between cations and anions in solution are subsequently optimized to capture their experimental osmotic pressures.
Collapse
Affiliation(s)
- Alireza Hosseni
- Department of Chemical and
Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, United States
| | - Henry S. Ashbaugh
- Department of Chemical and
Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, United States
| |
Collapse
|
2
|
Lüking M, van der Spoel D, Elf J, Tribello GA. Can molecular dynamics be used to simulate biomolecular recognition? J Chem Phys 2023; 158:2889489. [PMID: 37158325 DOI: 10.1063/5.0146899] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2023] [Accepted: 04/19/2023] [Indexed: 05/10/2023] Open
Abstract
There are many problems in biochemistry that are difficult to study experimentally. Simulation methods are appealing due to direct availability of atomic coordinates as a function of time. However, direct molecular simulations are challenged by the size of systems and the time scales needed to describe relevant motions. In theory, enhanced sampling algorithms can help to overcome some of the limitations of molecular simulations. Here, we discuss a problem in biochemistry that offers a significant challenge for enhanced sampling methods and that could, therefore, serve as a benchmark for comparing approaches that use machine learning to find suitable collective variables. In particular, we study the transitions LacI undergoes upon moving between being non-specifically and specifically bound to DNA. Many degrees of freedom change during this transition and that the transition does not occur reversibly in simulations if only a subset of these degrees of freedom are biased. We also explain why this problem is so important to biologists and the transformative impact that a simulation of it would have on the understanding of DNA regulation.
Collapse
Affiliation(s)
- Malin Lüking
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, SE-75124 Uppsala, Sweden
| | - David van der Spoel
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, SE-75124 Uppsala, Sweden
| | - Johan Elf
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, SE-75124 Uppsala, Sweden
| | - Gareth A Tribello
- Centre for Quantum Materials and Technologies, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, United Kingdom
| |
Collapse
|
3
|
Tillotson MJ, Diamantonis NI, Buda C, Bolton LW, Müller EA. Molecular modelling of the thermophysical properties of fluids: expectations, limitations, gaps and opportunities. Phys Chem Chem Phys 2023; 25:12607-12628. [PMID: 37114325 DOI: 10.1039/d2cp05423j] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
This manuscript provides an overview of the current state of the art in terms of the molecular modelling of the thermophysical properties of fluids. It is intended to manage the expectations and serve as guidance to practising physical chemists, chemical physicists and engineers in terms of the scope and accuracy of the more commonly available intermolecular potentials along with the peculiarities of the software and methods employed in molecular simulations while providing insights on the gaps and opportunities available in this field. The discussion is focused around case studies which showcase both the precision and the limitations of frequently used workflows.
Collapse
Affiliation(s)
- Marcus J Tillotson
- Department of Chemical Engineering, Imperial College London, London, UK.
| | | | | | | | - Erich A Müller
- Department of Chemical Engineering, Imperial College London, London, UK.
| |
Collapse
|
4
|
Blazquez S, Conde MM, Vega C. Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water. J Chem Phys 2023; 158:054505. [PMID: 36754806 DOI: 10.1063/5.0136498] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na+ and Cl- using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties.
Collapse
Affiliation(s)
- S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - C Vega
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
5
|
Paloncýová M, Pykal M, Kührová P, Banáš P, Šponer J, Otyepka M. Computer Aided Development of Nucleic Acid Applications in Nanotechnologies. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2022; 18:e2204408. [PMID: 36216589 DOI: 10.1002/smll.202204408] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/23/2022] [Revised: 09/12/2022] [Indexed: 06/16/2023]
Abstract
Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.
Collapse
Affiliation(s)
- Markéta Paloncýová
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Martin Pykal
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Petra Kührová
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Pavel Banáš
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Jiří Šponer
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
- Institute of Biophysics of the Czech Academy of Sciences, v. v. i., Královopolská 135, Brno, 612 65, Czech Republic
| | - Michal Otyepka
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
- IT4Innovations, VŠB - Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba, 708 00, Czech Republic
| |
Collapse
|
6
|
Abstract
![]()
We extend the modular AMBER lipid
force field to include anionic
lipids, polyunsaturated fatty acid (PUFA) lipids, and sphingomyelin,
allowing the simulation of realistic cell membrane lipid compositions,
including raft-like domains. Head group torsion parameters are revised,
resulting in improved agreement with NMR order parameters, and hydrocarbon
chain parameters are updated, providing a better match with phase
transition temperature. Extensive validation runs (0.9 μs per
lipid type) show good agreement with experimental measurements. Furthermore,
the simulation of raft-like bilayers demonstrates the perturbing effect
of increasing PUFA concentrations on cholesterol molecules. The force
field derivation is consistent with the AMBER philosophy, meaning
it can be easily mixed with protein, small molecule, nucleic acid,
and carbohydrate force fields.
Collapse
Affiliation(s)
- Callum J Dickson
- Computer-Aided Drug Discovery, Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Ross C Walker
- GlaxoSmithKline PLC, 1250 S. Collegeville Road, Collegeville, Pennsylvania 19426, United States.,Department of Chemistry and Biochemistry, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093, United States
| | - Ian R Gould
- Department of Chemistry, Imperial College London, London, SW7 2AZ, U.K
| |
Collapse
|
7
|
Blazquez S, Conde MM, Abascal JLF, Vega C. The Madrid-2019 force field for electrolytes in water using TIP4P/2005 and scaled charges: Extension to the ions F−, Br−, I−, Rb+, and Cs+. J Chem Phys 2022; 156:044505. [DOI: 10.1063/5.0077716] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- S. Blazquez
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - M. M. Conde
- Departamento de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - J. L. F. Abascal
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Vega
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
8
|
Silva CV, Santos JR, Melia Rodrigo M, Ribeiro AC, Valente AJ, Abreu PE, Marques JM, Esteso MA. On the transport and dynamics of disaccharides: H-bonding effect in sucrose and sucralose. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.117855] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
9
|
Elts E, Luxenburger F, Briesen H. Influence of Monovalent Salts on α-Glycine Crystal Growth from Aqueous Solution: Molecular Dynamics Simulations at Constant Supersaturation Conditions. J Phys Chem B 2021; 125:11732-11741. [PMID: 34643406 DOI: 10.1021/acs.jpcb.1c07168] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The growth of α-glycine crystals from aqueous solution is investigated at constant supersaturations by utilizing the constant chemical potential molecular dynamics method. The study considers two faces (010) and (011) that predominantly determine the α-glycine crystal morphology. The general Amber force field (GAFF) with two different charge sets derived from semi-empirical calculations using the complete neglect of differential overlap method (CNDO) and from density functional calculations using the double-numerical plus d- and p-polarization basis set (DNP) is applied to describe α-glycine. The extended simple point charge model is used to simulate water. It is observed that the GAFF/DNP set leads to a much slower integration of glycine molecules into the crystal structure than the GAFF/CNDO set. The GAFF/CNDO set, however, causes the growth even at concentrations well below the experimental solubility. For the GAFF/DNP set, the influence of potassium chloride (KCl) and sodium chloride (NaCl) on the face growth rates is investigated. The parameters recently proposed by Yagasaki et al. [J. Chem. Theory Comput. 2020, 16, 2460-2473] are used to describe salt ions, as standard GAFF parameters lead to the unexpected formation of salt clusters at a concentration lower than the experimental solubility value. According to our simulation results, both salts suppress the growth of the (011) and (010) faces. The inhibiting effect of NaCl is much stronger than that of KCl for the (011) face, while both salts have a similar inhibiting effect on the (010) face. The results are in line with the experimental observations of the impact of salt ions on the α-glycine growth rates for the (011) face reported in literature.
Collapse
Affiliation(s)
- Ekaterina Elts
- Chair of Process Systems Engineering, Technical University of Munich, 85354 Freising, Germany
| | - Frederik Luxenburger
- Chair of Process Systems Engineering, Technical University of Munich, 85354 Freising, Germany
| | - Heiko Briesen
- Chair of Process Systems Engineering, Technical University of Munich, 85354 Freising, Germany
| |
Collapse
|
10
|
Vodopivec AA, Chen Y, Russo PS, Hung FR. Molecular Dynamics Simulations of Nanostructures Formed by Hydrophobins and Oil in Seawater. J Phys Chem B 2021; 125:7886-7899. [PMID: 34236182 DOI: 10.1021/acs.jpcb.1c02040] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Classical molecular dynamics simulations using the Martini coarse-grained force field were performed to study oil nanodroplets surrounded by fungal hydrophobin (HP) proteins in seawater. The class I EAS and the class II HFBII HPs were studied along with two model oils, namely, benzene and n-decane. Both HPs exhibit free energy minima at the oil-seawater interface, which is deeper in benzene compared to the n-decane systems. Larger constraint forces are required to keep both HPs within the n-decane phase compared to inside benzene, with HFBII being more affine to benzene compared to EAS. Smaller surface tensions are observed at benzene-seawater interfaces coated with HPs compared to their n-decane counterparts. In the latter the surface tension remains unchanged upon increases in the concentration of HPs, whereas in benzene systems adding more HPs lead to decreases in surface tension. EAS has a larger tendency to cluster together in the interface compared to HFBII, with both HPs having larger coordination numbers when surrounding benzene droplets compared to when they are around n-decane nanoblobs. The HP-oil nanostructures in seawater examined have radii of gyration ranging between 2 and 12 nm, where the n-decane structures are larger and have more irregular shapes compared to the benzene systems. The n-decane molecules within the nanostructures form a compact spherical core, with the HPs partially covering its surface and clustering together, conferring irregular shapes to the nanostructures. The EAS with n-decane structures are larger and have more irregular shapes compared to their HFBII counterparts. In contrast, in the HP-benzene structures both HPs tend to penetrate the oil part of the droplet. The HFBII-benzene structures having the larger oil/HP ratios examined tend to be more compact and spherical compared to their EAS counterparts; however, some of the HFBII-benzene systems that have smaller oil/HP ratios have a more elongated structure compared to their EAS counterparts. This simulation study provides insights into HP-oil nanostructures that are smaller than the oil droplets and gas bubbles recently studied in experiments and, thus, might be challenging to examine with experimental techniques.
Collapse
Affiliation(s)
- Andrés A Vodopivec
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Yuwu Chen
- Center for Computation and Technology, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Paul S Russo
- School of Materials Science and Engineering and School of Chemistry & Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Francisco R Hung
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| |
Collapse
|
11
|
Spoel D, Zhang J, Zhang H. Quantitative predictions from molecular simulations using explicit or implicit interactions. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1560] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- David Spoel
- Uppsala Center for Computational Chemistry, Science for Life Laboratory, Department of Cell and Molecular Biology Uppsala University Uppsala Sweden
| | - Jin Zhang
- Department of Chemistry Southern University of Science and Technology Shenzhen China
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering University of Science and Technology Beijing Beijing China
| |
Collapse
|
12
|
Self-assembly of the surfactant mixtures on graphene in the presence of electrolyte: a molecular simulation study. ADSORPTION 2020. [DOI: 10.1007/s10450-020-00264-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
13
|
Kolafa J. Pressure in Molecular Simulations with Scaled Charges. 1. Ionic Systems. J Phys Chem B 2020; 124:7379-7390. [PMID: 32790401 DOI: 10.1021/acs.jpcb.0c02641] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Charge scaling, rationalized as MDEC (molecular dynamics in electronic continuum) or ECC (electronic continuum correction), has become a widely used simple approach to how to avoid self-consistent induced dipoles yet approximately take into account the effects of electronic polarizability. It has been assumed that the continuum permittivity does not depend on density; in turn, pressure is calculated by standard formulas. In this work, we elaborate a complementary approximation of density-independent molecular polarizability and derive formulas for pressure corrections within the MDEC framework; real behavior lies between these two extremes. The pressure corrections for test ionic systems are huge and negative, leading to sizable densities in constant-pressure MDEC simulations. A comparison of MDEC results with equivalent polarizable systems gives a good pressure match for a crystal but very low MDEC pressures for ionic liquids. These results witness about the importance of a correct density dependence not only of continuum permittivity in MDEC simulations but also of polarizability in polarizable simulations.
Collapse
Affiliation(s)
- Jiří Kolafa
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28 Praha 6, Czech Republic
| |
Collapse
|
14
|
Saric D, Kohns M, Vrabec J. Dielectric constant and density of aqueous alkali halide solutions by molecular dynamics: A force field assessment. J Chem Phys 2020; 152:164502. [PMID: 32357782 DOI: 10.1063/1.5144991] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The concentration dependence of the dielectric constant and the density of 11 aqueous alkali halide solutions (LiCl, NaCl, KCl, RbCl, CsCl, LiI, NaI, KI, CsI, KF, and CsF) is investigated by molecular simulation. Predictions using eight non-polarizable ion force fields combined with the TIP4P/ε water model are compared to experimental data. The influence of the water model and the temperature on the results for the NaCl brine are also addressed. The TIP4P/ε water model improves the accuracy of dielectric constant predictions compared to the SPC/E water model. The solution density is predicted well by most ion models. Almost all ion force fields qualitatively capture the decline of the dielectric constant with the increase of concentration for all solutions and with the increase of temperature for NaCl brine. However, the sampled dielectric constant is mostly in poor quantitative agreement with experimental data. These results are related to the microscopic solution structure, ion pairing, and ultimately the force field parameters. Ion force fields with excessive contact ion pairing and precipitation below the experimental solubility limit generally yield higher dielectric constant values. An adequate reproduction of the experimental solubility limit should therefore be a prerequisite for further investigations of the dielectric constant of aqueous electrolyte solutions by molecular simulation.
Collapse
Affiliation(s)
- Denis Saric
- Thermodynamics and Energy Technology, University of Paderborn, 33098 Paderborn, Germany
| | - Maximilian Kohns
- Laboratory of Engineering Thermodynamics, Technische Universität Kaiserslautern, 67633 Kaiserslautern, Germany
| | - Jadran Vrabec
- Thermodynamics and Process Engineering, Technical University Berlin, 10587 Berlin, Germany
| |
Collapse
|
15
|
Yagasaki T, Matsumoto M, Tanaka H. Lennard-Jones Parameters Determined to Reproduce the Solubility of NaCl and KCl in SPC/E, TIP3P, and TIP4P/2005 Water. J Chem Theory Comput 2020; 16:2460-2473. [PMID: 32207974 DOI: 10.1021/acs.jctc.9b00941] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Most classical nonpolarizable ion potential models underestimate the solubility values of NaCl and KCl in water significantly. We determine Lennard-Jones parameters of Na+, K+, and Cl- that reproduce the solubility as well as the hydration free energy in dilute aqueous solutions for three water potential models, SPC/E, TIP3P, and TIP4P/2005. The ion-oxygen distance in the solution and the cation-anion distance in salt are also considered in the parametrization. In addition to the target properties, the hydration enthalpy, hydration entropy, self-diffusion coefficient, coordination number, lattice energy, enthalpy of solution, density, viscosity, and number of contact ion pairs are calculated for comparison with 17 frequently used or recently developed ion potential models. The overall performance of each ion model is represented by a global score using a scheme that was originally developed for comparison of water potential models. The global score is better for our models than for the other 17 models not only because of the quite good prediction for the solubility but also because of the relatively small deviation from the experimental value for many of the other properties.
Collapse
Affiliation(s)
- Takuma Yagasaki
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan.,Department of Chemistry, Faculty of Science, Okayama University, Okayama 700-8530, Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan.,Department of Chemistry, Faculty of Science, Okayama University, Okayama 700-8530, Japan
| | - Hideki Tanaka
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan.,Department of Chemistry, Faculty of Science, Okayama University, Okayama 700-8530, Japan
| |
Collapse
|
16
|
Oliveira Bortot L, Bashardanesh Z, van der Spoel D. Making Soup: Preparing and Validating Models of the Bacterial Cytoplasm for Molecular Simulation. J Chem Inf Model 2019; 60:322-331. [DOI: 10.1021/acs.jcim.9b00971] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Leandro Oliveira Bortot
- Laboratory of Biological Physics, School of Pharmaceutical Sciences of Ribeirão Preto, University of São Paulo, Avenida do Café s/n, 14040-903 Ribeirão Preto-SP, Brazil
| | - Zahedeh Bashardanesh
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124 Uppsala, Sweden
| | - David van der Spoel
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124 Uppsala, Sweden
| |
Collapse
|
17
|
Tolmachev DA, Boyko OS, Lukasheva NV, Martinez-Seara H, Karttunen M. Overbinding and Qualitative and Quantitative Changes Caused by Simple Na+ and K+ Ions in Polyelectrolyte Simulations: Comparison of Force Fields with and without NBFIX and ECC Corrections. J Chem Theory Comput 2019; 16:677-687. [DOI: 10.1021/acs.jctc.9b00813] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
- D. A. Tolmachev
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, St. Petersburg 199004, Russia
| | - O. S. Boyko
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, St. Petersburg 199004, Russia
| | - N. V. Lukasheva
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, St. Petersburg 199004, Russia
| | - H. Martinez-Seara
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 542/2, Prague 6 CZ166 10, Czech Republic
| | - Mikko Karttunen
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, St. Petersburg 199004, Russia
| |
Collapse
|
18
|
Zeron IM, Abascal JLF, Vega C. A force field of Li +, Na +, K +, Mg 2+, Ca 2+, Cl -, and SO 4 2- in aqueous solution based on the TIP4P/2005 water model and scaled charges for the ions. J Chem Phys 2019; 151:134504. [PMID: 31594349 DOI: 10.1063/1.5121392] [Citation(s) in RCA: 131] [Impact Index Per Article: 26.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
In this work, a force field for several ions in water is proposed. In particular, we consider the cations Li+, Na+, K+, Mg2+, and Ca2+ and the anions Cl- and SO4 2-. These ions were selected as they appear in the composition of seawater, and they are also found in biological systems. The force field proposed (denoted as Madrid-2019) is nonpolarizable, and both water molecules and sulfate anions are rigid. For water, we use the TIP4P/2005 model. The main idea behind this work is to further explore the possibility of using scaled charges for describing ionic solutions. Monovalent and divalent ions are modeled using charges of 0.85 and 1.7, respectively (in electron units). The model allows a very accurate description of the densities of the solutions up to high concentrations. It also gives good predictions of viscosities up to 3 m concentrations. Calculated structural properties are also in reasonable agreement with the experiment. We have checked that no crystallization occurred in the simulations at concentrations similar to the solubility limit. A test for ternary mixtures shows that the force field provides excellent performance at an affordable computer cost. In summary, the use of scaled charges, which could be regarded as an effective and simple way of accounting for polarization (at least to a certain extend), improves the overall description of ionic systems in water. However, for purely ionic systems, scaled charges will not adequately describe neither the solid nor the melt.
Collapse
Affiliation(s)
- I M Zeron
- Depto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J L F Abascal
- Depto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Vega
- Depto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
19
|
Bachler J, Handle PH, Giovambattista N, Loerting T. Glass polymorphism and liquid-liquid phase transition in aqueous solutions: experiments and computer simulations. Phys Chem Chem Phys 2019; 21:23238-23268. [PMID: 31556899 DOI: 10.1039/c9cp02953b] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
One of the most intriguing anomalies of water is its ability to exist as distinct amorphous ice forms (glass polymorphism or polyamorphism). This resonates well with the possible first-order liquid-liquid phase transition (LLPT) in the supercooled state, where ice is the stable phase. In this Perspective, we review experiments and computer simulations that search for LLPT and polyamorphism in aqueous solutions containing salts and alcohols. Most studies on ionic solutes are devoted to NaCl and LiCl; studies on alcohols have mainly focused on glycerol. Less attention has been paid to protein solutions and hydrophobic solutes, even though they reveal promising avenues. While all solutions show polyamorphism and an LLPT only in dilute, sub-eutectic mixtures, there are differences regarding the nature of the transition. Isocompositional transitions for varying mole fractions are observed in alcohol but not in ionic solutions. This is because water can surround alcohol molecules either in a low- or high-density configuration whereas for ionic solutes, the water ion hydration shell is forced into high-density structures. Consequently, the polyamorphic transition and the LLPT are prevented near the ions, but take place in patches of water within the solutions. We highlight discrepancies and different interpretations within the experimental community as well as the key challenges that need consideration when comparing experiments and simulations. We point out where reinterpretation of past studies helps to draw a unified, consistent picture. In addition to the literature review, we provide original experimental results. A list of eleven open questions that need further consideration is identified.
Collapse
Affiliation(s)
- Johannes Bachler
- Institute of Physical Chemistry, University of Innsbruck, A-6020 Innsbruck, Austria.
| | | | | | | |
Collapse
|
20
|
Nikitin A, Del Frate G. Development of Nonbonded Models for Metal Cations Using the Electronic Continuum Correction. J Comput Chem 2019; 40:2464-2472. [PMID: 31301182 DOI: 10.1002/jcc.26021] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2019] [Revised: 05/14/2019] [Accepted: 06/06/2019] [Indexed: 12/24/2022]
Abstract
The parametrization of classical nonbonded models of metal ions has been widely addressed in the recent years. Despite the continuous development of novel and more physically inspired functional forms, the 12-6 Lennard-Jones plus Coulomb potential is still the most adopted force field in molecular dynamics (MD) codes, owing to its simple form and easy implementation. However, due to the integer formal charge, unpolarizable force fields of ions may suffer from overestimated interatomic electrostatic interactions, leading to nonphysical clustering or repulsion between such full charges. The electronic continuum correction (ECC) can fix this problem through a simple inclusion of solvent polarization effects via ionic charge rescaling. In this work, the development of novel nonbonded models for mono, divalent, and highly charged metal ions is presented. For each metal species, the ionic charge has been scaled, according to the ECC. Lennard-Jones parameters have been optimized using experimental structural and thermodynamic properties as target quantities. Performances of the proposed models are discussed and compared with the literature data, while transferability attitudes among different and well-known water models are evaluated. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Alexei Nikitin
- Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, Moscow, 119991, Russian Federation.,Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| | - Gianluca Del Frate
- IMT School for Advanced Studies Lucca, Piazza S. Francesco 19, I-55100, Lucca, Italy
| |
Collapse
|
21
|
Huraskin D, Horn AHC. Alkali ion influence on structure and stability of fibrillar amyloid-β oligomers. J Mol Model 2019; 25:37. [PMID: 30637529 DOI: 10.1007/s00894-018-3920-4] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2018] [Accepted: 12/26/2018] [Indexed: 12/25/2022]
Abstract
Alzheimer's disease is characterized by the aggregation of Amyloid-β (Aβ) peptide into oligomers, fibrils and plaques. Many factors influencing this process as well as the stability of the various Aβ aggregates are known to date, and include the concentration and type of metal ions. Most experimental and theoretical studies have concentrated on heavy metal ions, like Fe2+, Zn2+, or Cu2+, while the smaller alkali ions Li+, Na+, and K+ have not gained much attention notwithstanding their role and ubiquity in physiological environments. In this work, we applied atomistic molecular dynamics simulations to investigate the potential role of these alkali ions in stabilizing fibrillar Aβ oligomers of different size and topology, i.e., single and double filament systems comprising 3-24 peptide chains per filament. We find a pronounced difference on the molecular level in the interaction behavior with free carboxylate groups of the Aβ oligomer: Li+ forms stable bridged interactions, whereas K+ interacts more transiently and lacks bridging. The behavior of Na+ is in between, so that this ion-protein interaction obeys the renowned Hofmeister series. These differences are also reflected in the ability of the alkali ions to stabilize the oligomer secondary structure. The stabilizing effect is most pronounced for the smaller fibrillar oligomers, suggesting that the type of alkali ion critically affects the initial stages of fibril formation. Our findings thus offer a molecular explanation for the observation that the polymorphisms of Aβ fibril structures are caused by differences in the surrounding ionic environment. Graphical abstract Influence of alkali ions on the structure and stability of fibrillar amyloid-β oligomers.
Collapse
Affiliation(s)
- Danyil Huraskin
- Bioinformatik Institut für Biochemie Emil-Fischer-Zentrum, Friedrich-Alexander-Universität Erlangen-Nürnberg, Fahrstr. 17, 91054, Erlangen, Germany
| | - Anselm H C Horn
- Bioinformatik Institut für Biochemie Emil-Fischer-Zentrum, Friedrich-Alexander-Universität Erlangen-Nürnberg, Fahrstr. 17, 91054, Erlangen, Germany.
| |
Collapse
|
22
|
Walz MM, Ghahremanpour MM, van Maaren PJ, van der Spoel D. Phase-Transferable Force Field for Alkali Halides. J Chem Theory Comput 2018; 14:5933-5948. [DOI: 10.1021/acs.jctc.8b00507] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Marie-Madeleine Walz
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - Mohammad M. Ghahremanpour
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - Paul J. van Maaren
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - David van der Spoel
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| |
Collapse
|
23
|
Leonarski F, Jasiński M, Trylska J. Thermodynamics of the fourU RNA thermal switch derived from molecular dynamics simulations and spectroscopic techniques. Biochimie 2018; 156:22-32. [PMID: 30244136 DOI: 10.1016/j.biochi.2018.09.005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2018] [Accepted: 09/10/2018] [Indexed: 10/28/2022]
Abstract
Bacterial 5' untranslated regions of mRNA, termed thermal switches or thermometers, change their structure in response to temperature change. This structural change provides for the regulation of gene expression. One of such thermal switches, called fourU, is present in the Salmonella species. Mutations of fourU were found to abrogate its regulatory properties. We investigated the thermodynamics of the fourU fragment responsible for its structural changes. All-atom molecular dynamics simulations at various temperatures and spectroscopic experiments in solution were performed for the wild-type fourU and its mutants. We found that the U11C and A8C mutations stabilize the fourU structure in comparison to the wild-type fourU, and the double-point G14A/C25U mutant has the most destabilizing effect on the fourU hairpin 2 responsible for temperature sensing. The G14A/C25U mutant is also the easiest to strand-invade by a complementary oligonucleotide as indicated by fluorescence spectroscopy experiments.
Collapse
Affiliation(s)
- Filip Leonarski
- Centre of New Technologies, University of Warsaw, Banacha 2c, 02-097, Warsaw, Poland
| | - Maciej Jasiński
- Centre of New Technologies, University of Warsaw, Banacha 2c, 02-097, Warsaw, Poland
| | - Joanna Trylska
- Centre of New Technologies, University of Warsaw, Banacha 2c, 02-097, Warsaw, Poland.
| |
Collapse
|
24
|
Kulkarni M, Yang C, Pak Y. Refined Alkali Metal Ion Parameters for the OPC Water Model. B KOREAN CHEM SOC 2018. [DOI: 10.1002/bkcs.11527] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Mandar Kulkarni
- Department of Chemistry and Institute of Functional Materials; Pusan National University; Busan 609-735 South Korea
| | - Changwon Yang
- Department of Chemistry and Institute of Functional Materials; Pusan National University; Busan 609-735 South Korea
| | - Youngshang Pak
- Department of Chemistry and Institute of Functional Materials; Pusan National University; Busan 609-735 South Korea
| |
Collapse
|
25
|
Fischer NM, Polêto MD, Steuer J, van der Spoel D. Influence of Na+ and Mg2+ ions on RNA structures studied with molecular dynamics simulations. Nucleic Acids Res 2018; 46:4872-4882. [PMID: 29718375 PMCID: PMC6007214 DOI: 10.1093/nar/gky221] [Citation(s) in RCA: 62] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2017] [Revised: 02/16/2018] [Accepted: 04/23/2018] [Indexed: 01/11/2023] Open
Abstract
The structure of ribonucleic acid (RNA) polymers is strongly dependent on the presence of, in particular Mg2+ cations to stabilize structural features. Only in high-resolution X-ray crystallography structures can ions be identified reliably. Here, we perform molecular dynamics simulations of 24 RNA structures with varying ion concentrations. Twelve of the structures were helical and the others complex folded. The aim of the study is to predict ion positions but also to evaluate the impact of different types of ions (Na+ or Mg2+) and the ionic strength on structural stability and variations of RNA. As a general conclusion Mg2+ is found to conserve the experimental structure better than Na+ and, where experimental ion positions are available, they can be reproduced with reasonable accuracy. If a large surplus of ions is present the added electrostatic screening makes prediction of binding-sites less reproducible. Distinct differences in ion-binding between helical and complex folded structures are found. The strength of binding (ΔG‡ for breaking RNA atom-ion interactions) is found to differ between roughly 10 and 26 kJ/mol for the different RNA atoms. Differences in stability between helical and complex folded structures and of the influence of metal ions on either are discussed.
Collapse
Affiliation(s)
- Nina M Fischer
- Uppsala Centre for Computational Chemistry, Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - Marcelo D Polêto
- Uppsala Centre for Computational Chemistry, Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
- Center of Biotechnology, Universidade Federal do Rio Grande do Sul, Bento Gonçalves 9500, BR-91500-970 Porto Alegre, Brazil
| | - Jakob Steuer
- Uppsala Centre for Computational Chemistry, Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
- Department of Chemistry, University of Konstanz, Universitätstraße 10, D-78457 Konstanz, Germany
| | - David van der Spoel
- Uppsala Centre for Computational Chemistry, Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| |
Collapse
|
26
|
Šponer J, Bussi G, Krepl M, Banáš P, Bottaro S, Cunha RA, Gil-Ley A, Pinamonti G, Poblete S, Jurečka P, Walter NG, Otyepka M. RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview. Chem Rev 2018; 118:4177-4338. [PMID: 29297679 PMCID: PMC5920944 DOI: 10.1021/acs.chemrev.7b00427] [Citation(s) in RCA: 336] [Impact Index Per Article: 56.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Indexed: 12/14/2022]
Abstract
With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology , University of Copenhagen , Copenhagen 2200 , Denmark
| | - Richard A Cunha
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Nils G Walter
- Single Molecule Analysis Group and Center for RNA Biomedicine, Department of Chemistry , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| |
Collapse
|
27
|
Young JM, Panagiotopoulos AZ. System-Size Dependence of Electrolyte Activity Coefficients in Molecular Simulations. J Phys Chem B 2018; 122:3330-3338. [DOI: 10.1021/acs.jpcb.7b09861] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jeffrey M. Young
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | | |
Collapse
|
28
|
Sinha S, Bhattacharyya PK. Understanding the influence of external perturbation on aziridinium ion formation. Mol Phys 2018. [DOI: 10.1080/00268976.2017.1363922] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Sourab Sinha
- Department of Chemistry, Arya Vidyapeeth College, Guwahati, India
| | | |
Collapse
|
29
|
Benavides AL, Portillo MA, Chamorro VC, Espinosa JR, Abascal JLF, Vega C. A potential model for sodium chloride solutions based on the TIP4P/2005 water model. J Chem Phys 2017; 147:104501. [DOI: 10.1063/1.5001190] [Citation(s) in RCA: 60] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- A. L. Benavides
- Departamento de Ingeniería Física, División de Ciencias e Ingenierías, Universidad de Guanajuato, Loma del Bosque 103, Col. Lomas del Campestre, CP 37150 León, Mexico
| | - M. A. Portillo
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - V. C. Chamorro
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J. R. Espinosa
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J. L. F. Abascal
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Vega
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
30
|
Abstract
The 97 amino acid bitopic membrane protein M2 of influenza A forms a tetrameric bundle in which two of the monomers are covalently linked via a cysteine bridge. In its tetrameric assembly the protein conducts protons across the viral envelope and within intracellular compartments during the infectivity cycle of the virus. A key residue in the translocation of the protons is His-37 which forms a planar tetrad in the configuration of the bundle accepting and translocating the incoming protons from the N terminal side, exterior of the virus, to the C terminal side, inside the virus. With experimentally available data from NMR spectroscopy of the transmembrane domains of the tetrameric M2 bundle classical MD simulations are conducted with the protein bundle in different protonation stages in respect to His-37. A full correlation analysis (FCA) of the data sets with the His-37 tetrad either in a fully four times unprotonated or protonated state, assumed to mimic high and low pH in vivo, respectively, in both cases reveal asymmetric backbone dynamics. His-37 side chain rotation dynamics is increased at full protonation of the tetrad compared to the dynamics in the fully unprotonated state. The data suggest that proton translocation can be achieved by decoupled side chain or backbone dynamics. Graphical abstract Visualization of the tetrameric bundle of the transmembrane domains of M2 of influenza A after 200 ns of MD simulations (upper left). The four histidine residues 37 are either not protonated as in M20 or fully protonated is in M24+. The asymmetric dynamics of the backbones are shown after a full correlation analysis (FCA) in blue (lower left). The rotational dynamics of the χ2 dihedral angles of the histidines in M20 (upper right) are less than those in M24+ (lower right).
Collapse
Affiliation(s)
- Monoj Mon Kalita
- Institute of Biophotonics and Biophotonics & Molecular Imaging Research Center (BMIRC), School of Biomedical Science and Engineering, National Yang-Ming University, 155, Li-Nong St., Sec. 2, Taipei, 112, Taiwan
| | - Wolfgang B Fischer
- Institute of Biophotonics and Biophotonics & Molecular Imaging Research Center (BMIRC), School of Biomedical Science and Engineering, National Yang-Ming University, 155, Li-Nong St., Sec. 2, Taipei, 112, Taiwan.
| |
Collapse
|
31
|
Abstract
Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.
Collapse
Affiliation(s)
| | - Kenneth M. Merz
- Department of Chemistry, Department of Biochemistry and Molecular Biology, and Institute of Cyber-Enabled Research, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
32
|
Giri AK, Spohr E. Cluster formation of NaCl in bulk solutions: Arithmetic vs. geometric combination rules. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2016.09.089] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
33
|
Kolafa J. Solubility of NaCl in water and its melting point by molecular dynamics in the slab geometry and a new BK3-compatible force field. J Chem Phys 2016; 145:204509. [DOI: 10.1063/1.4968045] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
|
34
|
Choi JH, Cho M. Ion aggregation in high salt solutions. VI. Spectral graph analysis of chaotropic ion aggregates. J Chem Phys 2016; 145:174501. [DOI: 10.1063/1.4966246] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jun-Ho Choi
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 02841, South Korea
- Department of Chemistry, Korea University, Seoul 02841, South Korea
| | - Minhaeng Cho
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 02841, South Korea
- Department of Chemistry, Korea University, Seoul 02841, South Korea
| |
Collapse
|
35
|
Benavides AL, Aragones JL, Vega C. Consensus on the solubility of NaCl in water from computer simulations using the chemical potential route. J Chem Phys 2016; 144:124504. [PMID: 27036458 DOI: 10.1063/1.4943780] [Citation(s) in RCA: 83] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
The solubility of NaCl in water is evaluated by using three force field models: Joung-Cheatham for NaCl dissolved in two different water models (SPC/E and TIP4P/2005) and Smith Dang NaCl model in SPC/E water. The methodology based on free-energy calculations [E. Sanz and C. Vega, J. Chem. Phys. 126, 014507 (2007)] and [J. L. Aragones et al., J. Chem. Phys. 136, 244508 (2012)] has been used, except, that all calculations for the NaCl in solution were obtained by using molecular dynamics simulations with the GROMACS package instead of homemade MC programs. We have explored new lower molalities and made longer runs to improve the accuracy of the calculations. Exploring the low molality region allowed us to obtain an analytical expression for the chemical potential of the ions in solution as a function of molality valid for a wider range of molalities, including the infinite dilute case. These new results are in better agreement with recent estimations of the solubility obtained with other methodologies. Besides, two empirical simple rules have been obtained to have a rough estimate of the solubility of a certain model, by analyzing the ionic pairs formation as a function of molality and/or by calculating the difference between the NaCl solid chemical potential and the standard chemical potential of the salt in solution.
Collapse
Affiliation(s)
- A L Benavides
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J L Aragones
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - C Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
36
|
Rebič M, Laaksonen A, Šponer J, Uličný J, Mocci F. Molecular Dynamics Simulation Study of Parallel Telomeric DNA Quadruplexes at Different Ionic Strengths: Evaluation of Water and Ion Models. J Phys Chem B 2016; 120:7380-91. [PMID: 27379924 DOI: 10.1021/acs.jpcb.6b06485] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Most molecular dynamics (MD) simulations of DNA quadruplexes have been performed under minimal salt conditions using the Åqvist potential parameters for the cation with the TIP3P water model. Recently, this combination of parameters has been reported to be problematic for the stability of quadruplex DNA, especially caused by the ion interactions inside or near the quadruplex channel. Here, we verify how the choice of ion parameters and water model can affect the quadruplex structural stability and the interactions with the ions outside the channel. We have performed a series of MD simulations of the human full-parallel telomeric quadruplex by neutralizing its negative charge with K(+) ions. Three combinations of different cation potential parameters and water models have been used: (a) Åqvist ion parameters, TIP3P water model; (b) Joung and Cheatham ion parameters, TIP3P water model; and (c) Joung and Cheatham ion parameters, TIP4Pew water model. For the combinations (b) and (c), the effect of the ionic strength has been evaluated by adding increasing amounts of KCl salt (50, 100, and 200 mM). Two independent simulations using the Åqvist parameters with the TIP3P model show that this combination is clearly less suited for the studied quadruplex with K(+) as counterions. In both simulations, one ion escapes from the channel, followed by significant deformation of the structure, leading to deviating conformation compared to that in the reference crystallographic data. For the other combinations of ion and water potentials, no tendency is observed for the channel ions to escape from the quadruplex channel. In addition, the internal mobility of the three loops, torsion angles, and counterion affinity have been investigated at varied salt concentrations. In summary, the selection of ion and water models is crucial as it can affect both the structure and dynamics as well as the interactions of the quadruplex with its counterions. The results obtained with the TIP4Pew model are found to be closest to the experimental data at all of the studied ion concentrations.
Collapse
Affiliation(s)
- Matúš Rebič
- Department of Materials and Environmental Chemistry, Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University , 10691 Stockholm, Sweden.,Science for Life Laboratory (SciLifelab) , 17121 Solna, Sweden
| | - Aatto Laaksonen
- Department of Materials and Environmental Chemistry, Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University , 10691 Stockholm, Sweden.,Science for Life Laboratory (SciLifelab) , 17121 Solna, Sweden.,Stellenbosch Institute of Advanced Study (STIAS), Wallenberg Research Centre at Stellenbosch University , 7600 Stellenbosch, South Africa.,Department of Chemical and Geological Sciences, University of Cagliari , I-09042 Monserrato, Italy
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic (AVČR) , Kralovopolska 135, 612 65 Brno, Czech Republic.,Central European Institute of Technology (CEITEC) , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | | | - Francesca Mocci
- Department of Materials and Environmental Chemistry, Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University , 10691 Stockholm, Sweden.,Science for Life Laboratory (SciLifelab) , 17121 Solna, Sweden.,Department of Chemical and Geological Sciences, University of Cagliari , I-09042 Monserrato, Italy
| |
Collapse
|
37
|
Choi JH, Cho M. Ion aggregation in high salt solutions. V. Graph entropy analyses of ion aggregate structure and water hydrogen bonding network. J Chem Phys 2016; 144:204126. [DOI: 10.1063/1.4952648] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Jun-Ho Choi
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 02841, South Korea and Department of Chemistry, Korea University, Seoul 02841, South Korea
| | - Minhaeng Cho
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 02841, South Korea and Department of Chemistry, Korea University, Seoul 02841, South Korea
| |
Collapse
|
38
|
Affiliation(s)
- Roman Shevchuk
- Department
for Mathematics
and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
- Institute
for Microbiology
and Genetics, Georg-August-University Göttingen, Justus-von-Liebig-Weg 11, 37077 Göttingen, Germany
| |
Collapse
|
39
|
Kalita MM, Fischer WB. Asymmetric dynamics of ion channel forming proteins - Hepatitis C virus (HCV) p7 bundles. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:1462-70. [PMID: 27079148 DOI: 10.1016/j.bbamem.2016.04.004] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/17/2015] [Revised: 04/07/2016] [Accepted: 04/08/2016] [Indexed: 01/01/2023]
Abstract
Protein p7 of hepatitis C virus (HCV) is a short 63 amino acid membrane protein which homo-oligomerises in the lipid membrane to form ion and proton conducting bundles. Two different genotypes (GTs) of p7, 1a and 5a, are used to simulate hexameric bundles of the protein embedded in a fully hydrated lipid bilayer during 400 ns molecular dynamics (MD) simulations. Whilst the bundle of GT 1a is based on a fully computational derived structure, the bundle of GT 5a is based on NMR spectroscopic data. Results of a full correlation analysis (FCA) reveal that albeit structural differences both bundles screen local minima during the simulation. The collective motion of the protein domains is asymmetric. No 'breathing-mode'-like dynamics is observed. The presence of divalent ions, such as Ca-ions affects the dynamics of especially solvent exposed parts of the protein, but leaves the asymmetric domain motion unaffected.
Collapse
Affiliation(s)
- Monoj Mon Kalita
- Institute of Biophotonics and Biophotonics and Molecular Imaging Research Center (BMIRC), School of Biomedical Science and Engineering, National Yang-Ming University, Taipei, Taiwan, ROC
| | - Wolfgang B Fischer
- Institute of Biophotonics and Biophotonics and Molecular Imaging Research Center (BMIRC), School of Biomedical Science and Engineering, National Yang-Ming University, Taipei, Taiwan, ROC.
| |
Collapse
|
40
|
Nezbeda I, Moučka F, Smith WR. Recent progress in molecular simulation of aqueous electrolytes: force fields, chemical potentials and solubility. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1165296] [Citation(s) in RCA: 75] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- Ivo Nezbeda
- Faculty of Science, J.E. Purkinje University, Ústí nad Labem, Czech Republic
- Institute of Chemical Process Fundamentals, Academy of Sciences, Prague 6, Czech Republic
| | - Filip Moučka
- Faculty of Science, J.E. Purkinje University, Ústí nad Labem, Czech Republic
| | - William R. Smith
- Department of Mathematics and Statistics, University of Guelph, Guelph, ON, Canada
- Faculty of Science, University of Ontario Institute of Technology, Oshawa, ON, Canada
| |
Collapse
|
41
|
Daily MD, Baer MD, Mundy CJ. Divalent Ion Parameterization Strongly Affects Conformation and Interactions of an Anionic Biomimetic Polymer. J Phys Chem B 2016; 120:2198-208. [DOI: 10.1021/acs.jpcb.5b12277] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Michael D. Daily
- Physical
Science Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
- Department
of Biochemistry and Molecular Biology, University of Texas Medical Branch, Galveston, Texas 77555, United States
| | - Marcel D. Baer
- Physical
Science Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Christopher J. Mundy
- Physical
Science Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| |
Collapse
|
42
|
Abstract
Metal ions are essential cofactors for the structure and functions of nucleic acids. Yet, the early discovery in the 70s of the crucial role of Mg(2+) in stabilizing tRNA structures has occulted for a long time the importance of monovalent cations. Renewed interest in these ions was brought in the late 90s by the discovery of specific potassium metal ions in the core of a group I intron. Their importance in nucleic acid folding and catalytic activity is now well established. However, detection of K(+) and Na(+) ions is notoriously problematic and the question about their specificity is recurrent. Here we review the different methods that can be used to detect K(+) and Na(+) ions in nucleic acid structures such as X-ray crystallography, nuclear magnetic resonance or molecular dynamics simulations. We also discuss specific versus non-specific binding to different structures through various examples.
Collapse
Affiliation(s)
- Pascal Auffinger
- Architecture et Réactivité de l'ARN, Université de Strasbourg, IBMC, CNRS, 15 rue René Descartes, F-67084, Strasbourg, France.
| | - Luigi D'Ascenzo
- Architecture et Réactivité de l'ARN, Université de Strasbourg, IBMC, CNRS, 15 rue René Descartes, F-67084, Strasbourg, France.
| | - Eric Ennifar
- Architecture et Réactivité de l'ARN, Université de Strasbourg, IBMC, CNRS, 15 rue René Descartes, F-67084, Strasbourg, France.
| |
Collapse
|
43
|
Do TN, Choy WY, Karttunen M. Binding of Disordered Peptides to Kelch: Insights from Enhanced Sampling Simulations. J Chem Theory Comput 2015; 12:395-404. [PMID: 26636721 DOI: 10.1021/acs.jctc.5b00868] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Keap1 protein plays an essential role in regulating cellular oxidative stress response and is a crucial binding hub for multiple proteins, several of which are intrinsically disordered proteins (IDP). Among Kelch's IDP binding partners, NRF2 and PTMA are the two most interesting cases. They share a highly similar binding motif; however, NRF2 binds to Kelch with a binding affinity of approximately 100-fold higher than that of PTMA. In this study, we perform an exhaustive sampling composed of 6 μs well-tempered metadynamics and 2 μs unbiased molecular dynamics (MD) simulations aiming at characterizing the binding mechanisms and structural properties of these two peptides. Our results agree with previous experimental observations that PTMA is remarkably more disordered than NRF2 in both the free and bound states. This explains PTMA's lower binding affinity. Our extensive sampling also provides valuable insights into the vast conformational ensembles of both NRF2 and PTMA, supports the hypothesis of coupled folding-binding, and confirms the essential role of linear motifs in IDP binding.
Collapse
Affiliation(s)
- Trang Nhu Do
- Department of Chemistry and Waterloo Institute for Nanotechnology, University of Waterloo , 200 University Avenue West, Waterloo, ON, Canada N2L 3G1
| | - Wing-Yiu Choy
- Department of Biochemistry, University of Western Ontario , 1151 Richmond Street, London, ON, Canada N6A 3K7
| | - Mikko Karttunen
- Department of Mathematics and Computer Science & Institute for Complex Molecular Systems, Eindhoven University of Technology , P.O. Box 513, MetaForum, 5600 MB, Eindhoven, The Netherlands
| |
Collapse
|
44
|
Do TN, Choy WY, Karttunen M. Accelerating the Conformational Sampling of Intrinsically Disordered Proteins. J Chem Theory Comput 2015; 10:5081-94. [PMID: 26584388 DOI: 10.1021/ct5004803] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Intrinsically disordered proteins (IDPs) are a class of proteins lacking a well-defined secondary structure. Instead, they are able to attain multiple conformations, bind to multiple targets, and respond to changes in their surroundings. Functionally, IDPs have been associated with molecular recognition, cell regulation, and signal transduction. The dynamic conformational ensemble of IDPs is highly environmental and binding partner dependent, rendering the characterization of IDPs extremely challenging. Here, we compare the sampling efficiencies of conventional molecular dynamics (MD), well-tempered metadynamics (WT-META), and bias-exchange metadynamics (BE-META). The total simulation time was over 10 μs, and a 20-mer peptide derived from the Neh2 domain of the Nuclear factor erythroid 2-related factor 2 (Nrf2) protein was simulated. BE-META, with a neutral replica and seven biased replicas employing a set of seven relevant collective variables (CVs), provided the most reliable and efficient sampling. Finally, we propose a free-energy reconstruction method based on the probability distribution of the secondary structure contents. This postprocessing analysis confirms the presence of not only the β-hairpin conformation of the free Neh2 peptide but also its rare bound-state-like conformation, both of that have been experimentally observed. In addition, our simulations also predict other possible conformations to be verified with future experiments.
Collapse
Affiliation(s)
- Trang Nhu Do
- Department of Chemistry and Waterloo Institute for Nanotechnology, University of Waterloo , 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| | - Wing-Yiu Choy
- Department of Biochemistry, University of Western Ontario , 1151 Richmond Street, London, Ontario N6A 3K7, Canada
| | - Mikko Karttunen
- Department of Chemistry and Waterloo Institute for Nanotechnology, University of Waterloo , 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
45
|
Cordomí A, Edholm O, Perez JJ. Effect of Force Field Parameters on Sodium and Potassium Ion Binding to Dipalmitoyl Phosphatidylcholine Bilayers. J Chem Theory Comput 2015; 5:2125-34. [PMID: 26613152 DOI: 10.1021/ct9000763] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
The behavior of electrolytes in molecular dynamics simulations of zwitterionic phospholipid bilayers is very sensitive to the force field parameters used. Here, several 200 ns molecular dynamics of simulations of dipalmitoyl phosphotidylcholine (PC) bilayers in 0.2 M sodium or potassium chloride using various common force field parameters for the cations are presented. All employed parameter sets give a larger number of Na(+) ions than K(+) ions that bind to the lipid heads, but depending on the parameter choice quite different results are seen. A wide range of coordination numbers for the Na(+) and K(+) ions is also observed. These findings have been analyzed and compared to published experimental data. Some simulations produce aggregates of potassium chloride, indicating (in accordance with published simulations) that these force fields do not reproduce the delicate balance between salt and solvated ions. The differences between the force fields can be characterized by one single parameter, the electrostatic radius of the ion, which is correlated to σMO (M represents Na(+)/K(+)), the Lennard-Jones radius. When this parameter exceeds a certain threshold, binding to the lipid heads is no longer observed. One would, however, need more accurate experimental data to judge or rank the different force fields precisely. Still, reasons for the poor performance of some of the parameter sets are clearly demonstrated, and a quality control procedure is provided.
Collapse
Affiliation(s)
- Arnau Cordomí
- Department d'Enginyeria Química, Technical University of Catalonia (UPC), Avenue Diagonal 647, 08028 Barcelona, Spain, and, Theoretical Biological Physics, Royal Institute of Technology (KTH), SE-10691 Stockholm, Sweden
| | - Olle Edholm
- Department d'Enginyeria Química, Technical University of Catalonia (UPC), Avenue Diagonal 647, 08028 Barcelona, Spain, and, Theoretical Biological Physics, Royal Institute of Technology (KTH), SE-10691 Stockholm, Sweden
| | - Juan J Perez
- Department d'Enginyeria Química, Technical University of Catalonia (UPC), Avenue Diagonal 647, 08028 Barcelona, Spain, and, Theoretical Biological Physics, Royal Institute of Technology (KTH), SE-10691 Stockholm, Sweden
| |
Collapse
|
46
|
Giambaşu GM, Gebala MK, Panteva MT, Luchko T, Case DA, York DM. Competitive interaction of monovalent cations with DNA from 3D-RISM. Nucleic Acids Res 2015; 43:8405-15. [PMID: 26304542 PMCID: PMC4787805 DOI: 10.1093/nar/gkv830] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2015] [Revised: 08/05/2015] [Accepted: 08/06/2015] [Indexed: 12/15/2022] Open
Abstract
The composition of the ion atmosphere surrounding nucleic acids affects their folding, condensation and binding to other molecules. It is thus of fundamental importance to gain predictive insight into the formation of the ion atmosphere and thermodynamic consequences when varying ionic conditions. An early step toward this goal is to benchmark computational models against quantitative experimental measurements. Herein, we test the ability of the three dimensional reference interaction site model (3D-RISM) to reproduce preferential interaction parameters determined from ion counting (IC) experiments for mixed alkali chlorides and dsDNA. Calculations agree well with experiment with slight deviations for salt concentrations >200 mM and capture the observed trend where the extent of cation accumulation around the DNA varies inversely with its ionic size. Ion distributions indicate that the smaller, more competitive cations accumulate to a greater extent near the phosphoryl groups, penetrating deeper into the grooves. In accord with experiment, calculated IC profiles do not vary with sequence, although the predicted ion distributions in the grooves are sequence and ion size dependent. Calculations on other nucleic acid conformations predict that the variation in linear charge density has a minor effect on the extent of cation competition.
Collapse
Affiliation(s)
- George M Giambaşu
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University 174 Frelinghuysen Road, Piscataway, NJ 08854, USA
| | - Magdalena K Gebala
- Department of Biochemistry, Stanford University, Stanford, CA 94305, USA
| | - Maria T Panteva
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University 174 Frelinghuysen Road, Piscataway, NJ 08854, USA
| | - Tyler Luchko
- Department of Physics & Astronomy, California State University, Northridge, CA 91330, USA
| | - David A Case
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University 174 Frelinghuysen Road, Piscataway, NJ 08854, USA
| | - Darrin M York
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University 174 Frelinghuysen Road, Piscataway, NJ 08854, USA
| |
Collapse
|
47
|
Antila HS, Härkönen M, Sammalkorpi M. Chemistry specificity of DNA-polycation complex salt response: a simulation study of DNA, polylysine and polyethyleneimine. Phys Chem Chem Phys 2015; 17:5279-89. [PMID: 25607687 DOI: 10.1039/c4cp04967e] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
In this work, the chemistry specific stability determining factors of DNA-polycation complexes are examined by performing all-atom molecular dynamics simulations. To this end, we conduct a systematic variation of polycation line charge through polyethyleneimine (PEI) protonation and polycation chemistry via comparison with poly-l-lysine (PLL). Our simulations show that increasing line charge of the polycation alone does not lead to more salt tolerant complexes. Instead, the effective charge compensation by the polycation correlates with the increased stability of the complex against additional salt. The salt stability of PEI-DNA complexes also links to the proton sponge property of weak polycations, commonly assumed to be behind the effectivity of PEI as a gene delivery vector. Examination of the complexes reveals the mechanism behind this behaviour; more Cl(-) ions are attracted by the protonated complexes but, in contrast to the common depiction of the proton sponge behaviour, the ion influx does not cause swelling of the complex structure itself. However, PEI protonation leads to release of PEI while DNA remains tightly bound to the complex. Jointly, these findings shed light on the stability determining factors of DNA-polycation complexes, raise charge distribution as an important stability determining contributor, and indicate that the effectivity of PEI in gene delivery is likely to result from the freed PEI facilitating gene transfection.
Collapse
Affiliation(s)
- Hanne S Antila
- Aalto University School of Chemical Technology, Department of Chemistry, P.O. Box 16100, FI-00076, Aalto, Finland.
| | | | | |
Collapse
|
48
|
Choi JH, Cho M. Ion aggregation in high salt solutions. IV. Graph-theoretical analyses of ion aggregate structure and water hydrogen bonding network. J Chem Phys 2015; 143:104110. [DOI: 10.1063/1.4930608] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jun-Ho Choi
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
| | - Minhaeng Cho
- Department of Chemistry, Korea University, Seoul 136-713, South Korea
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Korea University, Seoul 136-713, South Korea
| |
Collapse
|
49
|
Bochicchio A, Rossetti G, Tabarrini O, Krauβ S, Carloni P. Molecular view of ligands specificity for CAG repeats in anti-Huntington therapy. J Chem Theory Comput 2015; 11:4911-22. [PMID: 26574279 DOI: 10.1021/acs.jctc.5b00208] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Huntington's disease is a fatal and devastating neurodegenerative genetic disorder for which there is currently no cure. It is characterized by Huntingtin protein's mRNA transcripts with 36 or more CAG repeats. Inhibiting the formation of pathological complexes between these expanded transcripts and target proteins may be a valuable strategy against the disease. Yet, the rational design of molecules specifically targeting the expanded CAG repeats is limited by the lack of structural information. Here, we use well-tempered metadynamics-based free energy calculations to investigate pose and affinity of two ligands targeting CAG repeats for which affinities have been previously measured. The first consists of two 4-guanidinophenyl rings linked by an ester group. It is the most potent ligand identified so far, with Kd = 60(30) nM. The second consists of a 4-phenyl dihydroimidazole and 4-1H-indole dihydroimidazole connected by a C-C bond (Kd = 700(80) nM). Our calculations reproduce the experimental affinities and uncover the recognition pattern between ligands' and their RNA target. They also provide a molecular basis for the markedly different affinity of the two ligands for CAG repeats as observed experimentally. These findings may pave the way for a structure-based hit-to-lead optimization to further improve ligand selectivity toward CAG repeat-containing mRNAs.
Collapse
Affiliation(s)
| | - Giulia Rossetti
- Department of Oncology, Hematology and Stem Cell Transplantation, RWTH Aachen University , D-52074 Aachen, North Rhine-Westphalia, Germany
| | - Oriana Tabarrini
- Department of Pharmaceutical Sciences, Università di Perugia , Via del Liceo 1, I-06123 Perugia, Perugia, Italy
| | - Sybille Krauβ
- German Center for Neurodegenerative Diseases (DZNE), Sigmund-Freud-Strasse 25, D-53127 Bonn, North Rhine-Westphalia, Germany
| | | |
Collapse
|
50
|
Henriksen NM, Fenley A, Gilson MK. Computational Calorimetry: High-Precision Calculation of Host-Guest Binding Thermodynamics. J Chem Theory Comput 2015; 11:4377-94. [PMID: 26523125 PMCID: PMC4614838 DOI: 10.1021/acs.jctc.5b00405] [Citation(s) in RCA: 76] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2015] [Indexed: 11/29/2022]
Abstract
We present a strategy for carrying out high-precision calculations of binding free energy and binding enthalpy values from molecular dynamics simulations with explicit solvent. The approach is used to calculate the thermodynamic profiles for binding of nine small molecule guests to either the cucurbit[7]uril (CB7) or β-cyclodextrin (βCD) host. For these systems, calculations using commodity hardware can yield binding free energy and binding enthalpy values with a precision of ∼0.5 kcal/mol (95% CI) in a matter of days. Crucially, the self-consistency of the approach is established by calculating the binding enthalpy directly, via end point potential energy calculations, and indirectly, via the temperature dependence of the binding free energy, i.e., by the van't Hoff equation. Excellent agreement between the direct and van't Hoff methods is demonstrated for both host-guest systems and an ion-pair model system for which particularly well-converged results are attainable. Additionally, we find that hydrogen mass repartitioning allows marked acceleration of the calculations with no discernible cost in precision or accuracy. Finally, we provide guidance for accurately assessing numerical uncertainty of the results in settings where complex correlations in the time series can pose challenges to statistical analysis. The routine nature and high precision of these binding calculations opens the possibility of including measured binding thermodynamics as target data in force field optimization so that simulations may be used to reliably interpret experimental data and guide molecular design.
Collapse
Affiliation(s)
- Niel M. Henriksen
- Skaggs School of Pharmacy
and Pharmaceutical Sciences, University
of California San Diego, La Jolla, California 92093-0736, United States
| | - Andrew
T. Fenley
- Skaggs School of Pharmacy
and Pharmaceutical Sciences, University
of California San Diego, La Jolla, California 92093-0736, United States
| | - Michael K. Gilson
- Skaggs School of Pharmacy
and Pharmaceutical Sciences, University
of California San Diego, La Jolla, California 92093-0736, United States
| |
Collapse
|