Sapna Sarupria

Homogeneous Nucleation of Methane Hydrate in Microsecond Molecular Dynamics Simulations.

We report atomistically detailed molecular dynamics simulations of homogeneous nucleation of methane hydrate in bulk aqueous phase in the absence of any interface. Subcritical clusters of water and methane molecules are formed in the initial segment of the simulations, which then aggregate to give the critical hydrate nucleus. This occurs over time scales of several hundred nanoseconds, indicating that the formation and aggregation of subcritical clusters can contribute significantly to the overall rate of hydrate nucleation. The clusters have elements of sI hydrate structure, such as 5(12) and 5(12)6(2) cages as well as other uncommon 5(12)6(3) and 5(12)6(4) cages, but do not possess long-range order. Clusters are dynamic in nature and undergo continuous structural rearrangements.

Molecular Dynamics Study of Carbon Dioxide Hydrate Dissociation

We present results from a molecular dynamics study of the dissociation behavior of carbon dioxide (CO(2)) hydrates. We explore the effects of hydrate occupancy and temperature on the rate of hydrate dissociation. We quantify the rate of dissociation by tracking CO(2) release into the liquid water phase as well as the velocity of the hydrate-liquid water interface. Our results show that the rate of dissociation is dependent on the fractional occupancy of each cage type and cannot be described simply in terms of overall hydrate occupancy. Specifically, we find that hydrates with similar overall occupancy differ in their dissociation behavior depending on whether the small or large cages are empty. In addition, individual cages behave differently depending on their surrounding environment. For the same overall occupancy, filled small and large cages dissociate faster in the presence of empty large cages than when empty small cages are present. Therefore, hydrate dissociation is a collective phenomenon that cannot be described by focusing solely on individual cage behavior.

Enthalpy-entropy contributions to salt and osmolyte effects on molecular-scale hydrophobic hydration and interactions.

Salts and additives can significantly affect the strength of water-mediated interactions in solution. We present results from molecular dynamics simulations focused on the thermodynamics of hydrophobic hydration, association, and the folding-unfolding of a hydrophobic polymer in water and in aqueous solutions of NaCl and of an osmolyte trimethylamine oxide (TMAO). It is known that addition of NaCl makes the hydration of hydrophobic solutes unfavorable and, correspondingly, strengthens their association at the pair as well as the many-body level (Ghosh, T.; Kalra, A.; Garde, S. J. Phys. Chem. B 2005, 109, 642), whereas the osmolyte TMAO has an almost negligible effect on the hydrophobic hydration and association (Athawale, M. V.; Dordick, J. S.; Garde, S. Biophys. J. 2005, 89, 858). Whether these effects are enthalpic or entropic in origin is not fully known. Here we perform temperature-dependent simulations to resolve the free energy into entropy and enthalpy contributions. We find that in TMAO solutions, there is an almost precise entropy-enthalpy compensation leading to the negligible effect of TMAO on hydrophobic phenomena. In contrast, in NaCl solutions, changes in enthalpy dominate, making the salt-induced strengthening of hydrophobic interactions enthalpic in origin. The resolution of total enthalpy into solute-solvent and solvent-solvent terms further shows that enthalpy changes originate primarily from solvent-solvent energy terms. Our results are consistent with experimental data on the hydration of small hydrophobic solutes by Ben-Naim and Yaacobi (Ben-Naim, A.; Yaacobi, M. J. Phys. Chem. 1974, 78, 170). In combination with recent work by Zangi, Hagen, and Berne (Zangi, R.; Hagen, M.; Berne, B. J. J. Am. Chem. Soc. 2007, 129, 4678) and the experimental data on surface tensions of salt solutions by Matubayasi et al. (Matubayasi, N.; Matsuo, H.; Yamamoto, K.; Yamaguchi, S.; Matuzawa, A. J. Colloid Interface Sci. 1999, 209, 398), our results highlight interesting length scale dependences of salt effects on hydrophobic phenomena. Although NaCl strengthens hydrophobic interactions at both small and large length scales, that effect is enthalpy-dominated at small length scales and entropy-dominated for large solutes and interfaces. Our results have implications for understanding of additive effects on water-mediated interactions, as well as on biocompatibility of osmolyte molecules in aqueous solutions.

Studying pressure denaturation of a protein by molecular dynamics simulations

Many globular proteins unfold when subjected to several kilobars of hydrostatic pressure. This “unfolding‐up‐on‐squeezing” is counter‐intuitive in that one expects mechanical compression of proteins with increasing pressure. Molecular simulations have the potential to provide fundamental understanding of pressure effects on proteins. However, the slow kinetics of unfolding, especially at high pressures, eliminates the possibility of its direct observation by molecular dynamics (MD) simulations. Motivated by experimental results—that pressure denatured states are water‐swollen, and theoretical results—that water transfer into hydrophobic contacts becomes favorable with increasing pressure, we employ a water insertion method to generate unfolded states of the protein Staphylococcal Nuclease (Snase). Structural characteristics of these unfolded states—their water‐swollen nature, retention of secondary structure, and overall compactness—mimic those observed in experiments. Using conformations of folded and unfolded states, we calculate their partial molar volumes in MD simulations and estimate the pressure‐dependent free energy of unfolding. The volume of unfolding of Snase is negative (approximately −60 mL/mol at 1 bar) and is relatively insensitive to pressure, leading to its unfolding in the pressure range of 1500–2000 bars. Interestingly, once the protein is sufficiently water swollen, the partial molar volume of the protein appears to be insensitive to further conformational expansion or unfolding. Specifically, water‐swollen structures with relatively low radii of gyration have partial molar volume that are similar to that of significantly more unfolded states. We find that the compressibility change on unfolding is negligible, consistent with experiments. We also analyze hydration shell fluctuations to comment on the hydration contributions to protein compressibility. Our study demonstrates the utility of molecular simulations in estimating volumetric properties and pressure stability of proteins, and can be potentially extended for applications to protein complexes and assemblies. Proteins 2010.

Suppression of sub-surface freezing in free-standing thin films of a coarse-grained model of water

Freezing in the vicinity of water-vapor interfaces is of considerable interest to a wide range of disciplines, most notably the atmospheric sciences. In this work, we use molecular dynamics and two advanced sampling techniques, forward flux sampling and umbrella sampling, to study homogeneous nucleation of ice in free-standing thin films of supercooled water. We use a coarse-grained monoatomic model of water, known as mW, and we find that in this model a vapor-liquid interface suppresses crystallization in its vicinity. This suppression occurs in the vicinity of flat interfaces where no net Laplace pressure in induced. Our free energy calculations reveal that the pre-critical crystalline nuclei that emerge near the interface are thermodynamically less stable than those that emerge in the bulk. We investigate the origin of this instability by computing the average asphericity of nuclei that form in different regions of the film, and observe that average asphericity increases closer to the interface, which is consistent with an increase in the free energy due to increased surface-to-volume ratios.

On the Thermodynamics and Kinetics of Hydrophobic Interactions at Interfaces

We have studied how primitive hydrophobic interactions between two or more small nonpolar solutes are affected by the presence of surfaces. We show that the desolvation barriers present in the potential of mean force between the solutes in bulk water are significantly reduced near an extended hydrophobic surface. Correspondingly, the kinetics of hydrophobic contact formation and breakage are faster near a hydrophobic surface than near a hydrophilic surface or in the bulk. We propose that the reduction in the desolvation barrier is a consequence of the fact that water near extended hydrophobic surfaces is akin to that at a liquid-vapor interface and is easily displaced. We support this proposal with three independent observations. First, when small hydrophobic solutes are brought near a hydrophobic surface, they induce local dewetting, thereby facilitating the reduction of desolvation barriers. Second, our results and those of Patel et al. (Proc. Natl. Acad. Sci. U.S.A. 2011, 108, 17678-17683) show that, whereas the association of small solutes in bulk water is driven by entropy, that near hydrophobic surfaces is driven by enthalpy, suggesting that the physics of interface deformation is important. Third, moving water away from its vapor-liquid coexistence, by applying hydrostatic pressure, leads to recovery of bulklike signatures (e.g., the presence of a desolvation barrier and an entropic driving force) in the association of solutes. These observations for simple solutes also translate to end-to-end contact formation in a model peptide with hydrophobic end groups, for which lowering of the desolvation barrier and acceleration of contact formation are observed near a hydrophobic surface. Our results suggest that extended hydrophobic surfaces, such as air-water or hydrocarbon-water surfaces, could serve as excellent platforms for catalyzing hydrophobically driven assembly.

PAMAM Dendrimers and Graphene: Materials for Removing Aromatic Contaminants from Water

We present results from experiments and atomistic molecular dynamics simulations on the remediation of naphthalene by polyamidoamine (PAMAM) dendrimers and graphene oxide (GrO). Specifically, we investigate 3rd-6th generation (G3-G6) PAMAM dendrimers and GrO with different levels of oxidation. The work is motivated by the potential applications of these emerging nanomaterials in removing polycyclic aromatic hydrocarbon contaminants from water. Our experimental results indicate that GrO outperforms dendrimers in removing naphthalene from water. Molecular dynamics simulations suggest that the prominent factors driving naphthalene association to these seemingly disparate materials are similar. Interestingly, we find that cooperative interactions between the naphthalene molecules play a significant role in enhancing their association to the dendrimers and GrO. Our findings highlight that while selection of appropriate materials is important, the interactions between the contaminants themselves can also be important in governing the effectiveness of a given material. The combined use of experiments and molecular dynamics simulations allows us to comment on the possible factors resulting in better performance of GrO in removing polyaromatic contaminants from water.

Hydrate Molecular Ballet

An exceptionally long computer simulation offers a glimpse of molecular events leading to methane hydrate formation. Hydrates are crystalline solids in which guest molecules are trapped within polyhedral water cages (1). They are important in hydrocarbon processing (2) and could play a major role in sustainable energy production (3, 4). Methane hydrate occurs naturally and in vast quantities on ocean floors and in permafrost, with implications for climate change and energy recovery (2). However, the molecular mechanisms leading to hydrate formation are poorly understood; this knowledge gap affects not just the science and technology of these materials, but our comprehension of hydrophobicity (5) and of disorder-order phase transitions. On page 1095 of this issue, Walsh et al. report a computational tour de force that offers a fascinating glimpse of the molecular events leading to methane hydrate formation (6).

The surface charge distribution affects the ice nucleating efficiency of silver iodide

Heterogeneous ice nucleation is the primary pathway for ice formation. However, the detailed molecular mechanisms by which surfaces promote or hinder ice nucleation are not well understood. We present results from extensive molecular dynamics simulations of ice nucleation near modified silver iodide (AgI) surfaces. The AgI surfaces are modified to investigate the effects of the surface charge distribution on the rate of ice nucleation. We find that the surface charge distribution has significant effects on ice nucleation. Specifically, AgI surfaces with the positive charges above the negative charges in the surface promote ice nucleation, while ice nucleation is hindered for surfaces in which the negative charges are above or in-plane with the positive charges. The structure of water molecules in the interfacial region as measured by the orientations of the water molecules relative to the surface can explain the differences in the ice nucleation at the different surfaces. We suggest that the distributions of the orientations of the interfacial water molecules could be used more broadly as a measure of ice nucleating propensity.

SciFlow: A dataflow-driven model architecture for scientific computing using Hadoop

Many computational science applications utilize complex workflow patterns that generate an intricately connected set of output files for subsequent analysis. Some types of applications, such as rare event sampling, additionally require guaranteed completion of all subtasks for analysis, and place significant demands on the workflow management and execution environment. SciFlow is a user interface built over the Hadoop infrastructure that provides a framework to support the complex process and data interactions and guaranteed completion requirements of scientific workflows. It provides an efficient mechanism for building a parallel scientific application with dataflow patterns, and enables the design, deployment, and execution of data intensive, many-task computing tasks on a Hadoop platform. The design principles of this framework emphasize simplicity, scalability and fault-tolerance. A case study using the forward flux sampling rare event simulation application validates the functionality, reliability and effectiveness of the framework.