Temperature Dependence and Energetics of Single Ions at the Aqueous Liquid–Vapor Interface

We investigate temperature-dependence of free energetics with two single halide anions, I^– and Cl^–, crossing the aqueous liquid–vapor interface through molecular dynamics simulations. The result shows that I^– has a modest surface stability of 0.5 kcal/mol at 300 K and the stability decreases as the temperature increases, indicating the surface adsorption process for the anion is entropically disfavored. In contrast, Cl^– shows no such surface state at all temperatures. Decomposition of free energetics reveals that water–water interactions provide a favorable enthalpic contribution, while the desolvation of ion induces an increase in free energy. Calculations of surface fluctuations demonstrate that I^– generates significantly greater interfacial fluctuations compared to Cl^–. The fluctuation is attributed to the malleability of the solvation shells, which allows for more long-ranged perturbations and solvent density redistribution induced by I^– as the anion approaches the liquid–vapor interface. The increase in temperature of the solvent enhances the inherent thermally excited fluctuations and consequently reduces the relative contribution from anion to surface fluctuations, which is consistent with the decrease in surface stability of I^–. Our results indicate a strong correlation with induced interfacial fluctuations and anion surface stability; moreover, resulting temperature dependent behavior of induced fluctuations suggests the possibility of a critical level of induced fluctuations associated with surface stability

Free Energetics of Arginine Permeation into Model DMPC Lipid Bilayers: Coupling of Effective Counterion Concentration and Lateral Bilayer Dimensions

Mechanisms and underlying thermodynamic determinants of translocation of charged cationic peptides such as cell-penetrating peptides across the cellular membrane continue to receive much attention. Two widely held views include endocytotic and non-endocytotic (diffusive) processes of permeant transfer across the bilayer. Considering a purely diffusive process, we consider the free energetics of translocation of a monoarginine peptide mimic across a model DMPC bilayer. We compute potentials of mean force for the transfer of a charged monoarginine peptide unit from water to the center of a 1,2-dimyristoyl-<i>sn</i>-glycero-3-phosphocholine (DMPC) model lipid bilayer. We use fully atomistic molecular dynamics simulations coupled with the adaptive biasing force (ABF) method for free energy estimation. The estimated potential of mean force difference from bulk to bilayer center is 6.94 ± 0.28 kcal/mol. The order of magnitude of this prediction is consistent with past experimental estimates of arginine partitioning into physiological bilayers in the context of translocon-based experiments, though the correlation between the bench and computer experiments is not unambiguous. Moreover, the present value is roughly one-half of previous estimates based on all-atom molecular dynamics free energy calculations. We trace the differences between the present and earlier calculations to system sizes used in the simulations and the dependence of the contributions to the free energy from various system components (water, lipids, ions, peptide) on overall system size. By varying the bilayer lateral dimensions in simulations using only sufficient numbers of counterions to maintain overall system charge neutrality, we find the possibility of an inherent convergent transfer free energy value

Liquid–Vapor Interfacial Properties of Aqueous Solutions of Guanidinium and Methyl Guanidinium Chloride: Influence of Molecular Orientation on Interface Fluctuations

The guanidinium cation (C­(NH₂)₃⁺) is a highly stable cation in aqueous solution due to its efficient solvation by water molecules and resonance stabilization of the charge. Its salts increase the solubility of nonpolar molecules (“salting-in”) and decrease the ordering of water. It is one of the strongest denaturants used in biophysical studies of protein folding. We investigate the behavior of guanidinium and its derivative, methyl guanidinium (an amino acid analogue) at the air–water surface, using atomistic molecular dynamics (MD) simulations and calculation of potentials of mean force. Methyl guanidinium cation is less excluded from the air–water surface than guanidinium cation, but both cations show orientational dependence of surface affinity. Parallel orientations of the guanidinium ring (relative to the Gibbs dividing surface) show pronounced free energy minima in the interfacial region, while ring orientations perpendicular to the GDS exhibit no discernible surface stability. Calculations of surface fluctuations demonstrate that, near the air–water surface, the parallel-oriented cations generate significantly greater interfacial fluctuations compared to other orientations, which induces more long-ranged perturbations and solvent density redistribution. Our results suggest a strong correlation with induced interfacial fluctuations and ion surface stability. These results have implications for interpreting molecular-level, mechanistic action of this osmolyte’s interaction with hydrophobic interfaces as they impact protein denaturation (solubilization)

Free Energetics of Carbon Nanotube Association in Pure and Aqueous Ionic Solutions

Carbon nanotubes are a promising platform across a broad spectrum of applications ranging from separations technology, drug delivery, to bio­(electronic) sensors. Proper dispersion of carbon nanotube materials is important to retaining the electronic properties of nanotubes. Experimentally it has been shown that salts can regulate the dispersing properties of CNTs in aqueous system with surfactants (Niyogi, S.; Densmore, C. G.; Doorn, S. K. <i>J. Am. Chem. Soc.</i> <b>2009</b>, <i>131</i>, 1144–1153); details of the physicochemical mechanisms underlying such effects continue to be explored. We address the effects of inorganic monovalent salts (NaCl and NaI) on dispersion stability of carbon nanotubes.We perform all-atom molecular dynamics simulations using nonpolarizable interaction models to compute the potential of mean force between two (10,10) single-walled carbon nanotubes (SWNTs) in the presence of NaCl/NaI and compare to the potential of mean force between SWNTs in pure water. Addition of salts enhances stability of the contact state between two SWNT’s on the order of 4 kcal/mol. The ion-specific spatial distribution of different halide anions gives rise to starkly different contributions to the free energy stability of nanotubes in the contact state. Iodide anion directly stabilizes the contact state to a much greater extent than chloride anion. The enhanced stability arises from the locally repulsive forces imposed on nanotubes by the surface-segregated iodide anion. Within the time scale of our simulations, both NaI and NaCl solutions stabilize the contact state by equivalent amounts. The marginally higher stability for contact state in salt solutions recapitulates results for small hydrophobic solutes in NaCl solutions (Athawale, M. V.; Sarupria, S.; Garde, S. <i>J. Phys. Chem. B</i> <b>2008</b>, <i>112</i>, 5661–5670) as well as single-walled carbon nanotubes in NaCl and CaCl₂ aqueous solutions

Free Energetics and the Role of Water in the Permeation of Methyl Guanidinium across the Bilayer–Water Interface: Insights from Molecular Dynamics Simulations Using Charge Equilibration Potentials

Combining umbrella sampling molecular dynamics (MD) simulations, the weighted histogram analysis method (WHAM) for unbiasing probabilities, and polarizable charge equilibration force fields, we compute the potential of mean force for the reversible transfer of methyl guanidinium from bulk solution to the center of a model DPPC bilayer. A 5 kcal/mol minimum in the potential of mean force profile for membrane permeation suggests that the analogue will preferentially reside in the headgroup region of the lipid, qualitatively in agreement with previously published results. We find the potential of mean force for permeation to be approximately 28 kcal/mol (relative to the minimum in the headgroups), within the range of values reported for similar types of simulations using fixed-charge force fields. From analysis of the lipid structure, we find that the lipid deformation leads to a substantial destabilizing contribution to the free energy of the methyl guanidinium as it resides in the bilayer center, though this deformation allows more efficient stabilization by water defects and transient pores. Water in the bilayer core stabilizes the charged residue. The role of water in stabilizing or destabilizing the solute as it crosses the bilayer depends on bulk electrolyte concentration. In 1 M KCl solution, the water contribution to the potential of mean force is stabilizing over the entire range of the permeation coordinate, with the sole destabilizing force originating from the anionic species in solution. Conversely, methyl guanidinium experiences net destabilization from water in the absence of electrolyte. The difference in solvent contributions to permeation free energy is traced to a local effect arising from differences in water density in the bilayer–water solution interface, thus leading to starkly opposite net forces on the permeant. The origin of the local water density differential rests with the penetration of hydrated chloride anions into the solution–bilayer interface. Finally, water permeation into the bilayer is required for the deformation of individual lipid molecules and permeation of ions into the membrane. From simulations where water is first excluded from the bilayer center where methyl guanidinium is restrained and then, after equilibration, allowed to enter the bilayer, we find that in the absence of any water defects/permeation into the bilayer, the lipid headgroups do not follow the methyl guanidinium. Only when water enters the bilayer do we see deformation of individual lipid molecules to associate with the amino acid analogue at bilayer center