101
|
Fyta M, Kalcher I, Dzubiella J, Vrbka L, Netz RR. Ionic force field optimization based on single-ion and ion-pair solvation properties. J Chem Phys 2010; 132:024911. [DOI: 10.1063/1.3292575] [Citation(s) in RCA: 66] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
|
102
|
Hughes ZE, Gale JD. A computational investigation of the properties of a reverse osmosis membrane. ACTA ACUST UNITED AC 2010. [DOI: 10.1039/c0jm01545h] [Citation(s) in RCA: 64] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
103
|
Maiti PK, Bagchi B. Diffusion of flexible, charged, nanoscopic molecules in solution: Size and pH dependence for PAMAM dendrimer. J Chem Phys 2009; 131:214901. [DOI: 10.1063/1.3266512] [Citation(s) in RCA: 59] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
104
|
Chen C, Li WZ, Song YC, Yang J. A molecular dynamics study of cryoprotective agent – Water–sodium chloride ternary solutions. ACTA ACUST UNITED AC 2009. [DOI: 10.1016/j.theochem.2009.09.007] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
105
|
Dzubiella J. Salt-Specific Stability of Short and Charged Alanine-Based α-Helices. J Phys Chem B 2009; 113:16689-94. [DOI: 10.1021/jp9077932] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Joachim Dzubiella
- Physics Department T37, Technical University Munich, 85748 Garching, Germany
| |
Collapse
|
106
|
Kalcher I, Horinek D, Netz RR, Dzubiella J. Ion specific correlations in bulk and at biointerfaces. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2009; 21:424108. [PMID: 21715843 DOI: 10.1088/0953-8984/21/42/424108] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
Ion specific effects are ubiquitous in any complex colloidal or biological fluid in bulk or at interfaces. The molecular origins of these 'Hofmeister effects' are not well understood and their theoretical description poses a formidable challenge to the modeling and simulation community. On the basis of the combination of atomistically resolved molecular dynamics (MD) computer simulations and statistical mechanics approaches, we present a few selected examples of specific electrolyte effects in bulk, at simple neutral and charged interfaces, and on a short α-helical peptide. The structural complexity in these strongly Coulomb-correlated systems is highlighted and analyzed in the light of available experimental data. While in general the comparison of MD simulations to experiments often lacks quantitative agreement, mostly because molecular force fields and coarse-graining procedures remain to be optimized, the consensus as regards trends provides important insights into microscopic hydration and binding mechanisms.
Collapse
Affiliation(s)
- I Kalcher
- Physics Department T37, Technical University Munich, 85748 Garching, Germany
| | | | | | | |
Collapse
|
107
|
Miettinen MS, Gurtovenko AA, Vattulainen I, Karttunen M. Ion dynamics in cationic lipid bilayer systems in saline solutions. J Phys Chem B 2009; 113:9226-34. [PMID: 19534449 DOI: 10.1021/jp810233q] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Positively charged lipid bilayer systems are a promising class of nonviral vectors for safe and efficient gene and drug delivery. Detailed understanding of these systems is therefore not only of fundamental but also of practical biomedical interest. Here, we study bilayers comprising a binary mixture of cationic dimyristoyltrimethylammoniumpropane (DMTAP) and zwitterionic (neutral) dimyristoylphosphatidylcholine (DMPC) lipids. Using atomistic molecular dynamics simulations, we address the effects of bilayer composition (cationic to zwitterionic lipid fraction) and of NaCl electrolyte concentration on the dynamical properties of these cationic lipid bilayer systems. We find that, despite the fact that DMPCs form complexes via Na(+) ions that bind to the lipid carbonyl oxygens, NaCl concentration has a rather minute effect on lipid diffusion. We also find the dynamics of Cl(-) and Na(+) ions at the water-membrane interface to differ qualitatively. Cl(-) ions have well-defined characteristic residence times of nanosecond scale. In contrast, the binding of Na(+) ions to the carbonyl region appears to lack a characteristic time scale, as the residence time distributions displayed power-law features. As to lateral dynamics, the diffusion of Na(+) ions within the water-membrane interface consists of two qualitatively different modes of motion: very slow diffusion when ions are bound to DMPC, punctuated by fast rapid jumps when detached from the lipids. Overall, the prolonged dynamics of the Na(+) ions are concluded to be interesting for the physics of the whole membrane, especially considering its interaction dynamics with charged macromolecular surfaces.
Collapse
Affiliation(s)
- Markus S Miettinen
- Department of Applied Physics, Helsinki University of Technology, Finland
| | | | | | | |
Collapse
|
108
|
Kaszuba K, Rog T, St Pierre JF, Mannisto PT, Karttunen M, Bunker A. Molecular dynamics study of prolyl oligopeptidase with inhibitor in binding cavity. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2009; 20:595-609. [PMID: 20024801 DOI: 10.1080/10629360903438198] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
We used the crystal structure of prolyl oligopeptidase (POP) with bound Z-pro-prolinal (ZPP) inhibitor (Protein Data Bank (PDB) structure 1QFS) to perform an intensive molecular dynamics study of the POP-ZPP complex. We performed 100 ns of simulation with the hemiacetal bond, through which the ZPP is bound to the POP, removed in order to better investigate the binding cavity environment. From basic analysis, measuring the radius of gyration, root mean square deviation, solvent accessible surface area and definition of the secondary structure of protein, we determined that the protein structure is highly stable and maintains its structure over the entire simulation time. This demonstrates that such long time simulations can be performed without the protein structure losing stability. We found that water bridges and hydrogen bonds play a negligible role in binding the ZPP thus indicating the importance of the hemiacetal bond. The two domains of the protein are bound by a set of approximately 12 hydrogen bonds, specific to the particular POP protein.
Collapse
Affiliation(s)
- K Kaszuba
- Department of Physics, Tampere University of Technology, FI-33101 Tampere, Finland
| | | | | | | | | | | |
Collapse
|
109
|
Kerisit S, Rosso KM. Transition path sampling of water exchange rates and mechanisms around aqueous ions. J Chem Phys 2009; 131:114512. [DOI: 10.1063/1.3224737] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
|
110
|
Fujita T, Fukuzawa K, Mochizuki Y, Nakano T, Tanaka S. Accuracy of fragmentation in ab initio calculations of hydrated sodium cation. Chem Phys Lett 2009. [DOI: 10.1016/j.cplett.2009.07.060] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
111
|
Hénin J, Shinoda W, Klein ML. Models for phosphatidylglycerol lipids put to a structural test. J Phys Chem B 2009; 113:6958-63. [PMID: 19371035 DOI: 10.1021/jp900645z] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Three atomistic empirical models for phosphatidylglycerol (PG) lipids are tested against structural data in the crystal and liquid crystal states. Simulations of the anhydrous crystal of dimyristoyl-phosphatidylglycerol (DMPG) show that only the CHARMM force field describes the conformation and interactions of PG head groups accurately. The other two models do not reproduce the native network of hydrogen bonds, suggesting the presence of biases in their conformational and nonbonded interaction properties. The CHARMM model is further validated in the biologically relevant liquid crystal phase by comparing experimental small-angle X-ray scattering spectra from DMPG unilamellar vesicles with data calculated from fluid bilayer simulations. The good agreement found in this model-free comparison implies that liquid crystal PG bilayers as described by CHARMM exhibit realistic bilayer thickness and lateral packing. Last, this model is used to simulate a fluid bilayer of palmitoyl-oleoyl-phosphatidylglycerol (POPG). The resulting view of the POPG bilayer structure is at variance with that proposed previously based on simulations, in particular, with respect to lateral packing of head groups and the role of counterions.
Collapse
Affiliation(s)
- Jérôme Hénin
- Center for Molecular Modeling, Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104-6323, USA.
| | | | | |
Collapse
|
112
|
Khandelia H, Mouritsen OG. Lipid gymnastics: evidence of complete acyl chain reversal in oxidized phospholipids from molecular simulations. Biophys J 2009; 96:2734-43. [PMID: 19348756 DOI: 10.1016/j.bpj.2009.01.007] [Citation(s) in RCA: 104] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2008] [Revised: 12/31/2008] [Accepted: 01/02/2009] [Indexed: 01/17/2023] Open
Abstract
In oxidative environments, biomembranes contain oxidized lipids with short, polar acyl chains. Two stable lipid oxidation products are PoxnoPC and PazePC. PoxnoPC has a carbonyl group, and PazePC has an anionic carboxyl group pendant at the end of the short, oxidized acyl chain. We have used MD simulations to explore the possibility of complete chain reversal in OXPLs in POPC-OXPL mixtures. The polar AZ chain of PazePC undergoes chain reversal without compromising the lipid bilayer integrity at concentrations up to 25% OXPL, and the carboxyl group points into the aqueous phase. Counterintuitively, the perturbation of overall membrane structural and dynamic properties is stronger for PoxnoPC than for PazePC. This is because of the overall condensing and ordering effect of sodium ions bound strongly to the lipids in the PazePC simulations. The reorientation of AZ chain is similar for two different lipid force fields. This work provides the first molecular evidence of the "extended lipid conformation" in phospholipid membranes. The chain reversal of PazePC lipids decorates the membrane interface with reactive, negatively charged functional groups. Such chain reversal is likely to exert a profound influence on the structure and dynamics of biological membranes, and on membrane-associated biological processes.
Collapse
Affiliation(s)
- Himanshu Khandelia
- MEMPHYS-Center for Biomembrane Physics, Department of Physics and Chemistry, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark.
| | | |
Collapse
|
113
|
Fennell CJ, Bizjak A, Vlachy V, Dill KA. Ion pairing in molecular simulations of aqueous alkali halide solutions. J Phys Chem B 2009; 113:6782-91. [PMID: 19206510 PMCID: PMC2698044 DOI: 10.1021/jp809782z] [Citation(s) in RCA: 220] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Using classical molecular dynamics simulations, we study ion-ion interactions in water. We study the potentials of mean force (PMF) for the full set of alkali halide ion pairs, and in each case, we test different parameter sets for modeling both the water and the ions. Altogether, we compared 300 different PMFs. We also calculate association equilibrium constants (KA) and compare them to two types of experiments. Of additional interest here was the proposition of Collins called the "law of matching water affinities", where the relative affinity of ions in solution depends on the matching of cation and anion sizes. From observations on the relative depths of the free energies of the contact ion pair (CIP) and the solvent-shared ion pair (SIP), along with related solvent structure analyses, we find a good correlation with this proposition: small-small and large-large should associate in water, and small-large should be more dissociated.
Collapse
Affiliation(s)
- Christopher J Fennell
- Department of Pharmaceutical Chemistry, University of California at San Francisco, San Francisco, CA 94143, USA.
| | | | | | | |
Collapse
|
114
|
Chialvo AA, Gruszkiewicz MS, Simonson JM, Palmer DA, Cole DR. Ion Pair Association in Extreme Aqueous Environments: Molecular-based and Electrical Conductance Approaches. J SOLUTION CHEM 2009. [DOI: 10.1007/s10953-009-9414-3] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
115
|
Alejandre J, Chapela GA, Bresme F, Hansen JP. The short range anion-H interaction is the driving force for crystal formation of ions in water. J Chem Phys 2009; 130:174505. [DOI: 10.1063/1.3124184] [Citation(s) in RCA: 60] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
116
|
Corradini D, Gallo P, Rovere M. Effect of concentration on the thermodynamics of sodium chloride aqueous solutions in the supercooled regime. J Chem Phys 2009; 130:154511. [DOI: 10.1063/1.3119634] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
117
|
Chen C, Li W, Song Y, Yang J. Molecular dynamics simulation studies of cryoprotective agent solutions: the relation between melting temperature and the ratio of hydrogen bonding acceptor to donor number. Mol Phys 2009. [DOI: 10.1080/00268970902852632] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
118
|
Kalcher I, Dzubiella J. Structure-thermodynamics relation of electrolyte solutions. J Chem Phys 2009; 130:134507. [DOI: 10.1063/1.3097530] [Citation(s) in RCA: 93] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
|
119
|
Sammalkorpi M, Karttunen M, Haataja M. Ionic Surfactant Aggregates in Saline Solutions: Sodium Dodecyl Sulfate (SDS) in the Presence of Excess Sodium Chloride (NaCl) or Calcium Chloride (CaCl2). J Phys Chem B 2009; 113:5863-70. [DOI: 10.1021/jp901228v] [Citation(s) in RCA: 174] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Maria Sammalkorpi
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544, and Department of Applied Mathematics, The University of Western Ontario, London, Ontario, Canada N6A 5B7
| | - Mikko Karttunen
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544, and Department of Applied Mathematics, The University of Western Ontario, London, Ontario, Canada N6A 5B7
| | - Mikko Haataja
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544, and Department of Applied Mathematics, The University of Western Ontario, London, Ontario, Canada N6A 5B7
| |
Collapse
|
120
|
Song Y, Gunner M. Using Multiconformation Continuum Electrostatics to Compare Chloride Binding Motifs in α-Amylase, Human Serum Albumin, and Omp32. J Mol Biol 2009; 387:840-56. [DOI: 10.1016/j.jmb.2009.01.038] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
121
|
Horinek D, Mamatkulov SI, Netz RR. Rational design of ion force fields based on thermodynamic solvation properties. J Chem Phys 2009; 130:124507. [PMID: 19334851 DOI: 10.1063/1.3081142] [Citation(s) in RCA: 188] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Dominik Horinek
- Physik Department, Technische Universität München, 85748 Garching, Germany.
| | | | | |
Collapse
|
122
|
Yamaguchi T, Matsuoka T, Koda S. A theoretical study on the frequency-dependent electric conductivity of electrolyte solutions. II. Effect of hydrodynamic interaction. J Chem Phys 2009; 130:094506. [DOI: 10.1063/1.3085717] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
|
123
|
A short guide for molecular dynamics simulations of RNA systems. Methods 2009; 47:187-97. [DOI: 10.1016/j.ymeth.2008.09.020] [Citation(s) in RCA: 64] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2008] [Accepted: 09/19/2008] [Indexed: 01/13/2023] Open
|
124
|
Wu Y, Joseph S, Aluru NR. Effect of Cross-Linking on the Diffusion of Water, Ions, and Small Molecules in Hydrogels. J Phys Chem B 2009; 113:3512-20. [DOI: 10.1021/jp808145x] [Citation(s) in RCA: 141] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Yanbin Wu
- Department of Mechanical Science and Engineering, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana−Champaign, Urbana, IL-61801
| | - Sony Joseph
- Department of Mechanical Science and Engineering, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana−Champaign, Urbana, IL-61801
| | - N. R. Aluru
- Department of Mechanical Science and Engineering, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana−Champaign, Urbana, IL-61801
| |
Collapse
|
125
|
Cruz-Chu ER, Aksimentiev A, Schulten K. Ionic Current Rectification Through Silica Nanopores. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2009; 113:1850. [PMID: 20126282 PMCID: PMC2658614 DOI: 10.1021/jp804724p] [Citation(s) in RCA: 60] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
Nanopores immersed in electrolytic solution and under the influence of an electric field can produce ionic current rectification, where ionic currents are higher for one voltage polarity than for the opposite polarity, resulting in an asymmetric current-voltage (I-V) curve. This behavior has been observed in polymer and silicon-based nanopores as well as in theoretically studied continuum models. By means of atomic level molecular dynamics (MD) simulations, we have performed a systematic investigation of KCl conductance in silica nanopores with a total simulation time of 680 ns. We found that ion-binding spots at the silica surfaces, such as dangling atoms, have effects on the ion concentration and electrostatic potential inside the nanopore, producing asymmetric I-V curves. Conversely, silica surfaces without ion-binding spots produce symmetric I-V curves.
Collapse
|
126
|
Dong H, Liu W, Doren DJ, Wood RH. Structure of an Accurate ab Initio Model of the Aqueous Na+ Ion at High Temperatures. J Phys Chem B 2008; 112:13552-60. [DOI: 10.1021/jp8057405] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Haitao Dong
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716
| | - Wenbin Liu
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716
| | - Douglas J. Doren
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716
| | - Robert H. Wood
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716
| |
Collapse
|
127
|
Dzubiella J. Salt-Specific Stability and Denaturation of a Short Salt-Bridge-Forming α-Helix. J Am Chem Soc 2008; 130:14000-7. [DOI: 10.1021/ja805562g] [Citation(s) in RCA: 79] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Joachim Dzubiella
- Physics Department T37, Technical University Munich, 85748 Garching, Germany
| |
Collapse
|
128
|
Hassan SA. Computer simulation of ion cluster speciation in concentrated aqueous solutions at ambient conditions. J Phys Chem B 2008; 112:10573-84. [PMID: 18680338 PMCID: PMC2561909 DOI: 10.1021/jp801147t] [Citation(s) in RCA: 56] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Dynamic simulations are used to investigate ion cluster formation in unsaturated aqueous NaCl at 25 degrees C. Statistical, structural, and dynamic properties are reported. An effort is made to identify general behaviors that are expected to hold beyond the limitations of the force field. Above approximately 1 M, clusters with more than ten ions begin to form after approximately 10-20 ns of simulation time, but no evidence of irreversible ion aggregation is observed. Cluster survival times are estimated, showing that the kinetics become increasingly complex as salt is added, leading to multiple decay rates. Cluster dipole moment distributions show characteristic peaks that reflect the preferred conformations of clusters in solution. These are modulated by electrostatic and liquid-structure forces and are described in detail for clusters of up to five ions. For a given size and charge, the cluster morphology is independent of salt concentration. Below approximately 2 M, clusters affect the structure of water in their first hydration shells, so dipole moments parallel to the cluster macrodipoles are induced. These effects show a weak dependence with concentration below approximately 2 M, but vanish in the 2-3 M range. A possible connection with the structural transition recently suggested by NMR data in concentrated electrolytes is discussed. The effects of electrostatics on cluster speciation and morphology are discussed based on results from a set of simulations carried out with the ionic charges removed.
Collapse
Affiliation(s)
- Sergio A Hassan
- Center for Molecular Modeling, Division of Computational Bioscience, CIT National Institutes of Health, U.S. DHHS, Bethesda, MD 20892, USA
| |
Collapse
|
129
|
Dahlberg M, Maliniak A. Molecular dynamics simulations of cardiolipin bilayers. J Phys Chem B 2008; 112:11655-63. [PMID: 18712912 DOI: 10.1021/jp803414g] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Cardiolipin is a key lipid component in the inner mitochondrial membrane, where the lipid is involved in energy production, cristae structure, and mechanisms in the apoptotic pathway. In this article we used molecular dynamics computer simulations to investigate cardiolipin and its effect on the structure of lipid bilayers. Three cardiolipin/POPC bilayers with different lipid compositions were simulated: 100, 9.2, and 0% cardiolipin. We found strong association of sodium counterions to the carbonyl groups of both lipid types, leaving in the case of 9.2% cardiolipin virtually no ions in the aqueous compartment. Although binding occurred primarily at the carbonyl position, there was a preference to bind to the carbonyl groups of cardiolipin. Ion binding and the small headgroup of cardiolipin gave a strong ordering of the hydrocarbon chains. We found significant effects in the water dipole orientation and water dipole potential which can compensate for the electrostatic repulsion that otherwise should force charged lipids apart. Several parameters relevant for the molecular structure of cardiolipin were calculated and compared with results from analyses of coarse-grained simulations and available X-ray structural data.
Collapse
Affiliation(s)
- Martin Dahlberg
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, Sweden.
| | | |
Collapse
|
130
|
Joung IS, Cheatham, TE. Determination of alkali and halide monovalent ion parameters for use in explicitly solvated biomolecular simulations. J Phys Chem B 2008; 112:9020-41. [PMID: 18593145 PMCID: PMC2652252 DOI: 10.1021/jp8001614] [Citation(s) in RCA: 2345] [Impact Index Per Article: 146.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2008] [Revised: 04/30/2008] [Indexed: 12/11/2022]
Abstract
Alkali (Li(+), Na(+), K(+), Rb(+), and Cs(+)) and halide (F(-), Cl(-), Br(-), and I(-)) ions play an important role in many biological phenomena, roles that range from stabilization of biomolecular structure, to influence on biomolecular dynamics, to key physiological influence on homeostasis and signaling. To properly model ionic interaction and stability in atomistic simulations of biomolecular structure, dynamics, folding, catalysis, and function, an accurate model or representation of the monovalent ions is critically necessary. A good model needs to simultaneously reproduce many properties of ions, including their structure, dynamics, solvation, and moreover both the interactions of these ions with each other in the crystal and in solution and the interactions of ions with other molecules. At present, the best force fields for biomolecules employ a simple additive, nonpolarizable, and pairwise potential for atomic interaction. In this work, we describe our efforts to build better models of the monovalent ions within the pairwise Coulombic and 6-12 Lennard-Jones framework, where the models are tuned to balance crystal and solution properties in Ewald simulations with specific choices of well-known water models. Although it has been clearly demonstrated that truly accurate treatments of ions will require inclusion of nonadditivity and polarizability (particularly with the anions) and ultimately even a quantum mechanical treatment, our goal was to simply push the limits of the additive treatments to see if a balanced model could be created. The applied methodology is general and can be extended to other ions and to polarizable force-field models. Our starting point centered on observations from long simulations of biomolecules in salt solution with the AMBER force fields where salt crystals formed well below their solubility limit. The likely cause of the artifact in the AMBER parameters relates to the naive mixing of the Smith and Dang chloride parameters with AMBER-adapted Aqvist cation parameters. To provide a more appropriate balance, we reoptimized the parameters of the Lennard-Jones potential for the ions and specific choices of water models. To validate and optimize the parameters, we calculated hydration free energies of the solvated ions and also lattice energies (LE) and lattice constants (LC) of alkali halide salt crystals. This is the first effort that systematically scans across the Lennard-Jones space (well depth and radius) while balancing ion properties like LE and LC across all pair combinations of the alkali ions and halide ions. The optimization across the entire monovalent series avoids systematic deviations. The ion parameters developed, optimized, and characterized were targeted for use with some of the most commonly used rigid and nonpolarizable water models, specifically TIP3P, TIP4P EW, and SPC/E. In addition to well reproducing the solution and crystal properties, the new ion parameters well reproduce binding energies of the ions to water and the radii of the first hydration shells.
Collapse
Affiliation(s)
| | - Thomas E. Cheatham,
- To whom correspondence should be addressed: Phone: (801) 587-9652. Fax: (801) 585-9119. E-mail:
| |
Collapse
|
131
|
Savelyev A, Papoian GA. Polyionic Charge Density Plays a Key Role in Differential Recognition of Mobile Ions by Biopolymers. J Phys Chem B 2008; 112:9135-45. [DOI: 10.1021/jp801448s] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Alexey Savelyev
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599-3290
| | - Garegin A. Papoian
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599-3290
| |
Collapse
|
132
|
Corradini D, Gallo P, Rovere M. Thermodynamic behavior and structural properties of an aqueous sodium chloride solution upon supercooling. J Chem Phys 2008; 128:244508. [DOI: 10.1063/1.2939118] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
133
|
Eriksson M, Lindhorst TK, Hartke B. Differential effects of oligosaccharides on the hydration of simple cations. J Chem Phys 2008; 128:105105. [PMID: 18345929 DOI: 10.1063/1.2873147] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Changed ion hydration properties near surfaces, proteins, and deoxyribose nucleic acid have been reported before in the literature. In the present work, we extend this work to carbohydrates: We have performed classical-mechanical molecular dynamics simulations to study solvation properties of simple cations of biological relevance (Na(+),K(+),Mg(2+),Ca(2+)) in explicit water, near single and multiple oligosaccharides as glycocalyx models. We find that our oligosaccharides prefer direct contact with K(+) over Na(+), but that the Na(+) contacts are longer lived. These interactions also lead to strong but short-lived changes in oligosaccharide conformations, with oligosaccharides wrapping around K(+) with multiple contacts. These findings may have implications for current hypotheses on glycocalyx functions.
Collapse
Affiliation(s)
- Mats Eriksson
- Institut für Physikalische Chemie, Christian-Albrechts-Universität, Olshausenstrasse 40, Kiel, Germany
| | | | | |
Collapse
|
134
|
|
135
|
Wang M, Liu J, Chen S. Electric potential distribution in nanoscale electroosmosis: from molecules to continuum. MOLECULAR SIMULATION 2008. [DOI: 10.1080/08927020701663321] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
136
|
Gurtovenko AA, Vattulainen I. Effect of NaCl and KCl on phosphatidylcholine and phosphatidylethanolamine lipid membranes: insight from atomic-scale simulations for understanding salt-induced effects in the plasma membrane. J Phys Chem B 2008; 112:1953-62. [PMID: 18225878 DOI: 10.1021/jp0750708] [Citation(s) in RCA: 207] [Impact Index Per Article: 12.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
To gain a better understanding of how monovalent salt under physiological conditions affects plasma membranes, we have performed 200 ns atomic-scale molecular dynamics simulations of phosphatidylcholine (PC) and phosphatidylethanolamine (PE) lipid bilayers. These two systems provide representative models for the outer and inner leaflets of the plasma membrane, respectively. The implications of cation-lipid interactions in these lipid systems have been considered in two different aqueous salt solutions, namely NaCl and KCl, and the sensitivity of the results on the details of interactions used for ions is determined by repeating the simulations with two distinctly different force fields. We demonstrate that the main effect of monovalent salt on a phospholipid membrane is determined by cations binding to the carbonyl region of a membrane, while chloride anions mostly stay in the water phase. It turns out that the strength and character of the cation-lipid interactions are quite different for different types of lipids and cations. PC membranes and Na+ ions demonstrate strongest interactions, leading to notable membrane compression. This finding was confirmed by both force fields (Gromacs and Charmm) employed for the ions. The binding of potassium ions to PC membranes (and the overall effect of KCl), in turn, was found to be much weaker mainly due to the larger size of a K+ ion compared to Na+. Furthermore, the effect of KCl on PC membranes was found to be force-field sensitive: The binding of a potassium ion was not observed at all in simulations performed with the Gromacs force-field, which seems to exaggerate the size of a K+ ion. As far as PE lipid bilayers are concerned, they are found to be influenced by monovalent salt to a significantly lesser extent compared to PC bilayers, which is a direct consequence of the ability of PE lipids to form both intra- and intermolecular hydrogen bonds and hence to adopt a more densely packed bilayer structure. Whereas for NaCl we observed weak binding of Na+ cations to the PE lipid-water interface, in the case of KCl we witnessed almost complete lack of cation binding. Overall, our findings indicate that monovalent salt ions affect lipids in the inner and outer leaflets of plasma cell membranes in substantially different ways.
Collapse
Affiliation(s)
- Andrey A Gurtovenko
- Computational Biophysics Laboratory, Institute of Pharmaceutical Innovation, University of Bradford, West Yorkshire BD7 1DP, UK.
| | | |
Collapse
|
137
|
Cordomí A, Edholm O, Perez JJ. Effect of ions on a dipalmitoyl phosphatidylcholine bilayer. a molecular dynamics simulation study. J Phys Chem B 2008; 112:1397-408. [PMID: 18193856 DOI: 10.1021/jp073897w] [Citation(s) in RCA: 116] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The effect of physiological concentrations of different chlorides on the structure of a dipalmitoyl phosphatidylcholine (DPPC) bilayer has been investigated through atomistic molecular dynamics simulations. These calculations provide support to the concept that Li+, Na+, Ca2+, Mg2+, Sr2+, Ba2+, and Ac3+, but not K+, bind to the lipid-head oxygens. Ion binding exhibits an influence on lipid order, area per lipid, orientation of the lipid head dipole, the charge distribution in the system, and therefore the electrostatic potential across the head-group region of the bilayer. These structural effects are sensitive to the specific characteristics of each cation, i.e., radius, charge, and coordination properties. These results provide evidence aimed at shedding some light into the apparent contradictions among different studies reported recently regarding the ordering effect of ions on zwitterionic phosphatidylcholine lipid bilayers.
Collapse
Affiliation(s)
- Arnau Cordomí
- Department of Chemical Engineering, Technical University of Catalonia (UPC), Av. Diagonal 647, 08028 Barcelona, Spain
| | | | | |
Collapse
|
138
|
Karttunen M, Rottler J, Vattulainen I, Sagui C. Chapter 2 Electrostatics in Biomolecular Simulations: Where Are We Now and Where Are We Heading? CURRENT TOPICS IN MEMBRANES 2008. [DOI: 10.1016/s1063-5823(08)00002-1] [Citation(s) in RCA: 55] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
139
|
Structure of spheroidal HDL particles revealed by combined atomistic and coarse-grained simulations. Biophys J 2007; 94:2306-19. [PMID: 18065479 DOI: 10.1529/biophysj.107.115857] [Citation(s) in RCA: 70] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Spheroidal high-density lipoprotein (HDL) particles circulating in the blood are formed through an enzymatic process activated by apoA-I, leading to the esterification of cholesterol, which creates a hydrophobic core of cholesteryl ester molecules in the middle of the discoidal phospholipid bilayer. In this study, we investigated the conformation of apoA-I in model spheroidal HDL (ms-HDL) particles using both atomistic and coarse-grained molecular dynamics simulations, which are found to provide consistent results for all HDL properties we studied. The observed small contribution of cholesteryl oleate molecules to the solvent-accessible surface area of the entire ms-HDL particle indicates that palmitoyloleoylphosphatidylcholines and apoA-I molecules cover the hydrophobic core comprised of cholesteryl esters particularly well. The ms-HDL particles are found to form a prolate ellipsoidal shape, with sizes consistent with experimental results. Large rigid domains and low mobility of the protein are seen in all the simulations. Additionally, the average number of contacts of cholesteryl ester molecules with apoA-I residues indicates that cholesteryl esters interact with protein residues mainly through their cholesterol moiety. We propose that the interaction of annular cholesteryl oleate molecules contributes to apoA-I rigidity stabilizing and regulating the structure and function of the ms-HDL particle.
Collapse
|
140
|
Thomas AS, Elcock AH. Molecular Dynamics Simulations of Hydrophobic Associations in Aqueous Salt Solutions Indicate a Connection between Water Hydrogen Bonding and the Hofmeister Effect. J Am Chem Soc 2007; 129:14887-98. [DOI: 10.1021/ja073097z] [Citation(s) in RCA: 135] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Andrew S. Thomas
- Contribution from the Department of Biochemistry, University of Iowa, Iowa City, Iowa 52242
| | - Adrian H. Elcock
- Contribution from the Department of Biochemistry, University of Iowa, Iowa City, Iowa 52242
| |
Collapse
|
141
|
Stimson L, Dong L, Karttunen M, Wisniewska A, Dutka M, Róg T. Stearic Acid Spin Labels in Lipid Bilayers: Insight through Atomistic Simulations. J Phys Chem B 2007; 111:12447-53. [DOI: 10.1021/jp0746796] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Lorna Stimson
- Laboratory of Physics, Helsinki University of Technology, Finland, Department of Applied Mathematics, The University of Western Ontario, London (ON), Canada, Department of Biophysics, Faculty of Biochemistry, Biophysics and Biotechnology, Jagiellonian University, Kraków, Poland, and Biophysics and Statistical Mechanics Group, Department of Electrical and Communication Engineering, Helsinki University of Technology, Finland
| | - Lei Dong
- Laboratory of Physics, Helsinki University of Technology, Finland, Department of Applied Mathematics, The University of Western Ontario, London (ON), Canada, Department of Biophysics, Faculty of Biochemistry, Biophysics and Biotechnology, Jagiellonian University, Kraków, Poland, and Biophysics and Statistical Mechanics Group, Department of Electrical and Communication Engineering, Helsinki University of Technology, Finland
| | - Mikko Karttunen
- Laboratory of Physics, Helsinki University of Technology, Finland, Department of Applied Mathematics, The University of Western Ontario, London (ON), Canada, Department of Biophysics, Faculty of Biochemistry, Biophysics and Biotechnology, Jagiellonian University, Kraków, Poland, and Biophysics and Statistical Mechanics Group, Department of Electrical and Communication Engineering, Helsinki University of Technology, Finland
| | - Anna Wisniewska
- Laboratory of Physics, Helsinki University of Technology, Finland, Department of Applied Mathematics, The University of Western Ontario, London (ON), Canada, Department of Biophysics, Faculty of Biochemistry, Biophysics and Biotechnology, Jagiellonian University, Kraków, Poland, and Biophysics and Statistical Mechanics Group, Department of Electrical and Communication Engineering, Helsinki University of Technology, Finland
| | - Małgorzata Dutka
- Laboratory of Physics, Helsinki University of Technology, Finland, Department of Applied Mathematics, The University of Western Ontario, London (ON), Canada, Department of Biophysics, Faculty of Biochemistry, Biophysics and Biotechnology, Jagiellonian University, Kraków, Poland, and Biophysics and Statistical Mechanics Group, Department of Electrical and Communication Engineering, Helsinki University of Technology, Finland
| | - Tomasz Róg
- Laboratory of Physics, Helsinki University of Technology, Finland, Department of Applied Mathematics, The University of Western Ontario, London (ON), Canada, Department of Biophysics, Faculty of Biochemistry, Biophysics and Biotechnology, Jagiellonian University, Kraków, Poland, and Biophysics and Statistical Mechanics Group, Department of Electrical and Communication Engineering, Helsinki University of Technology, Finland
| |
Collapse
|
142
|
Sammalkorpi M, Karttunen M, Haataja M. Structural Properties of Ionic Detergent Aggregates: A Large-Scale Molecular Dynamics Study of Sodium Dodecyl Sulfate. J Phys Chem B 2007; 111:11722-33. [PMID: 17877384 DOI: 10.1021/jp072587a] [Citation(s) in RCA: 161] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The properties of sodium dodecyl sulfate (SDS) aggregates were studied through extensive molecular dynamics simulations with explicit solvent. First, we provide a parametrization of the model within Gromacs. Then, we probe the kinetics of aggregation by starting from a random solution of SDS molecules and letting the system explore its kinetic pathway during the aggregation of multiple units. We observe a structural transition for the surfactant aggregates brought upon by a change in temperature. Specifically, at low temperatures, the surfactants form crystalline aggregates, whereas at elevated temperatures, they form micelles. We also investigate the dependence of aggregation kinetics on surfactant concentration and report on the molecular level structural changes involved in the transition.
Collapse
Affiliation(s)
- Maria Sammalkorpi
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, NJ, USA
| | | | | |
Collapse
|
143
|
Dong H, Liu W, Doren D, Wood R. Structure of an accurate ab initio model of the aqueous Cl- ion at high temperatures. J Phys Chem B 2007; 110:18504-14. [PMID: 16970478 DOI: 10.1021/jp0628333] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The structure of an accurate ab initio model of aqueous chloride ion was calculated at two high-temperature state points (573 K, 0.725 g/cm(3) and 723 K, 0.0098 g/cm(3)) by a two-step procedure. First, the structure of an approximate model was calculated from a molecular dynamics simulation of the model. Then the difference between the structure of the ab initio model and the approximate model was calculated by non-Boltzmann weighting of a sample of configurations taken from the approximate model simulation. Radial distribution functions, average coordination numbers, the distribution of coordination numbers, an analysis of orientations of water in the first coordination shell, and the free energy of hydration of the chloride ion are reported for both state points. The most common water structure has one hydrogen close to the chloride ion and one pointing away (46% at 573 K and 57% at 723 K). Waters in the first coordination shell that are not strongly bound to the chloride ions are common. Several variations of the method were tested. Models in which the water-water interaction is calculated with ab initio methods predict only a slightly different structure than models in which water-water interactions are determined from the approximate models. Similarly, using the approximate model for solute-water interactions when the water is far from the chloride ion did not affect the results. Uncertainties due to the limited sample of configurations are estimated and found to be small. The results are in qualitative agreement with X-ray and neutron diffraction experiments and with simulations of approximate models.
Collapse
Affiliation(s)
- Haitao Dong
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, USA
| | | | | | | |
Collapse
|
144
|
Baker CM, Grant GH. The effect of solvation on biomolecular conformation: 2-amino-1-phenylethanol. J Phys Chem B 2007; 111:9940-54. [PMID: 17672488 DOI: 10.1021/jp071059w] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Small molecule neurotransmitters form one the most important classes of pharmaceutical molecules. While the behavior of these molecules in their neutral forms in the gas phase is well understood, their behavior in more biologically relevant scenarios (protonated and in aqueous solution) has received comparatively little attention. Here we address this problem by using molecular mechanics simulations to build up a detailed picture of the conformational behavior of 2-amino-1-phenylethanol, a noradrenaline analogue, in aqueous solution in both its neutral and protonated forms. For the sake of comparison, equivalent simulations are also performed on the gas-phase molecules and gas-phase hydrated clusters. These calculations reveal the important role that water has to play in determining the conformational preferences and dynamic behavior of the molecules. Water molecules are found to bridge between the various functional groups within the molecule, significantly affecting their relative stabilities in comparison to the gas-phase values. The reorganization of these solvation structures also provides a mechanism for conformational interconversion. The role of the solvent in mediating interactions between the various functional groups within the molecule suggests that in noradrenaline the catechol groups will be able to interact, albeit indirectly, with the other functional groups, thereby influencing the behavior of the molecule.
Collapse
Affiliation(s)
- Christopher M Baker
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, The University of Oxford, South Parks Road, Oxford, United Kingdom OX1 3QZ
| | | |
Collapse
|
145
|
Chen AA, Pappu RV. Quantitative characterization of ion pairing and cluster formation in strong 1:1 electrolytes. J Phys Chem B 2007; 111:6469-78. [PMID: 17518490 DOI: 10.1021/jp0708547] [Citation(s) in RCA: 97] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Aqueous solutions of 1:1 strong electrolytes are considered to be the prototype for complete ionic dissociation. Nonetheless, clustering of strong 1:1 electrolytes has been widely reported in all atom molecular dynamics simulations, and their presence is indirectly implicated in a diverse range of experimental results. Is there a physical basis for nonidealities such as ion pairing and cluster formation in aqueous solutions of strong 1:1 electrolytes? We attempt to answer this question by direct comparison of results from detailed molecular dynamics simulations to experimentally observed properties of 1:1 electrolytes. We report the analysis of a series of lengthy molecular dynamics simulations of alkali-halide solutions carried out over a wide range of physiologically relevant concentrations using explicit representations of water molecules. We find evidence for pronounced nonideal behavior of ions at all concentrations in the form of ion pairs and clusters which are in rapid equilibrium with dissociated ions. The phenomenology for ion pairing seen in these simulations is congruent with the multistep scheme proposed by Eigen and Tamm based on data from ultrasonic absorption experiments. For a given electrolyte, we show that the dependence of cluster populations on concentration can be described through a single set of equilibrium constants. We assess the accuracy of calculated ion pairing constants by favorable comparison to estimates obtained by Fuoss and co-workers and based on conductometric experiments. Ion pairs and clusters form on length scales where the size of individual water molecules is as important as the hard core radius of ions. Ion pairing results as a balance between the favorable Coulomb interactions and the unfavorable partial desolvation of ions needed to form a pair.
Collapse
Affiliation(s)
- Alan A Chen
- Molecular Biophysics Program & Center for Computational Biology, Washington University School of Medicine, St. Louis, Missouri 63110, USA
| | | |
Collapse
|
146
|
Savelyev A, Papoian GA. Inter-DNA Electrostatics from Explicit Solvent Molecular Dynamics Simulations. J Am Chem Soc 2007; 129:6060-1. [PMID: 17455935 DOI: 10.1021/ja070207t] [Citation(s) in RCA: 59] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Alexey Savelyev
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599-3290, USA
| | | |
Collapse
|
147
|
Lenart PJ, Jusufi A, Panagiotopoulos AZ. Effective potentials for 1:1 electrolyte solutions incorporating dielectric saturation and repulsive hydration. J Chem Phys 2007; 126:044509. [PMID: 17286489 DOI: 10.1063/1.2431169] [Citation(s) in RCA: 55] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Implicit water potentials are developed for the study of thermodynamic and structural properties of solutions of NaCl, LiCl, and KCl. The interaction potential between cations and anions is parametrized from the ionic crystal potential. Two short-range corrections were added to the system to account for the water solvent. The first is due to dielectric saturation which reduces the dielectric permittivity in the vicinity of an ion. The second is a repulsive Gaussian potential which represents the first hydration shell around the ions. Grand canonical Monte Carlo simulations were performed to calculate the mean ionic activity coefficients. Molecular dynamics simulations were performed to calculate the radial distribution functions of 1.0 molal solutions at 298 K which were used to compare the structure of the explicit and implicit water simulations. The implementation of dielectric saturation and a repulsive hydration potential results in an excellent description of the mean activity coefficient and is able to capture structural features of contact ion pairs and solvent separated ions.
Collapse
Affiliation(s)
- Philip J Lenart
- Department of Chemical Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | | | | |
Collapse
|
148
|
Yang Y, Meng S. Atomistic nature of NaCl nucleation at the solid-liquid interface. J Chem Phys 2007; 126:044708. [PMID: 17286500 DOI: 10.1063/1.2431363] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
The early stage of heterogeneous nucleation of NaCl from supersaturated NaCl aqueous solution at the water-NaCl (001) interface has been investigated by molecular dynamics simulations. The critical size of the nuclei for spontaneous growth was found to be as small as two atoms (a Na(+)-Cl(-) ion pair) at high supersaturation. Due to the presence of a relatively stable water network and the effect of the hydration force at the interface, the stable nuclei formed on the NaCl (001) are found to contain more Na(+) ions than Cl(-) ions. The different deposition characteristics of the Na(+) and Cl(-) solutes lead to a positively charged substrate and thus may introduce another driving force for nucleation besides the level of solution supersaturation. The role of water was further confirmed by comparison with NaCl epitaxy growth in the vacuum.
Collapse
Affiliation(s)
- Yong Yang
- Institute of Physics, Chinese Academy of Sciences, Box 603, Beijing 100080, China
| | | |
Collapse
|
149
|
Gurtovenko AA, Vattulainen I. Ion leakage through transient water pores in protein-free lipid membranes driven by transmembrane ionic charge imbalance. Biophys J 2007; 92:1878-90. [PMID: 17208976 PMCID: PMC1861780 DOI: 10.1529/biophysj.106.094797] [Citation(s) in RCA: 100] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2006] [Accepted: 12/04/2006] [Indexed: 11/18/2022] Open
Abstract
We have employed atomic-scale molecular dynamics simulations to address ion leakage through transient water pores in protein-free phospholipid membranes. Our results for phospholipid membranes in aqueous solution with NaCl and KCl salts show that the formation of transient water pores and the consequent ion leakage can be induced and be driven by a transmembrane ionic charge imbalance, an inherent feature in living cells. These processes take place if the gradient is large enough to develop a sufficiently significant potential difference across the membrane. The transport of cations and anions through the water pores is then seen; it discharges the transmembrane potential, considerably reduces the size of a water pore, and makes the water pore metastable, leading eventually to its sealing. The ion transport is found to be sensitive to the type of ions. It turns out that Na(+) and Cl(-) ions leak through a membrane at approximately the same ratio despite the fact that Na(+) ions are expected to experience a lower potential barrier for the permeation through the pore. This is because of strong interactions of sodium ions with the carbonyl region of a phospholipid membrane as well as with lipid headgroups forming pore "walls," considerably slowing down the permeation of sodium ions. In contrast, we observed a pronounced selectivity of a phospholipid membrane to the permeation of potassium ions as compared to chloride ions: Potassium ions, being larger than sodium ions, interact only weakly with phospholipid headgroups, so that these interactions are not able to compensate for a large difference in free-energy barriers for permeation of K(+) and Cl(-) ions. These findings are found to be robust to a choice of force-field parameters for ions (tested by Gromacs and Charmm force-fields for ions). What is more, a potassium ion is found to be able to permeate a membrane along an alternate, "water-defect-mediated" pathway without actual formation of a pore. The "water-defect-mediated" leakage involves formation of a single water defect only and is found to be at least one order of magnitude faster than the pore-mediated ion leakage.
Collapse
Affiliation(s)
- Andrey A Gurtovenko
- Computational Laboratory, Institute of Pharmaceutical Innovation, University of Bradford, Bradford, West Yorkshire, BD7 1DP, United Kingdom.
| | | |
Collapse
|
150
|
Free energy calculations of counterion partitioning between DNA and chloride solutions. MENDELEEV COMMUNICATIONS 2007. [DOI: 10.1016/j.mencom.2007.03.015] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|