Divya Nayar

Water and water-like liquids: relationships between structure, entropy and mobility.

Liquids with very diverse underlying interactions share the thermodynamic and transport anomalies of water, including metalloids, ionic melts and mesoscopic fluids. The generic feature that characterises such water-like liquids is a density-driven shift in the nature of local order in the condensed phases. The key semiquantitative relationships between structural order, thermodynamics and transport that are necessary in order to map out the consequences of this common qualitative feature for liquid-state properties and phase transformations of such systems are reviewed here. The application of these ideas to understand and model tetrahedral liquids, especially water, is discussed and possible extensions to other complex fluids are considered.

Comparison of Tetrahedral Order, Liquid State Anomalies, and Hydration Behavior of mTIP3P and TIP4P Water Models.

The relationship between local tetrahedral order, tagged particle potential energy, and coordination number is studied for mTIP3P and TIP4P models of water in the bulk as well as in the neighborhood of a small peptide. The tendency of water molecules with different binding or tagged particle potential energies to occupy environments with different degrees of disorder can be effectively illustrated by constructing tetrahedral order distributions and corresponding entropy metrics conditional on restricted ranges of local binding energy. At the state point corresponding to the onset of the density anomaly, the correlation between tetrahedral entropy versus tagged potential energy is strong and virtually identical for mTIP3P and TIP4P. In TIP4P, this correlation is retained up to temperatures as high as 300 K, while it is lost by 250 K in mTIP3P. In the 250-300 K regime that is important for biomolecular simulations, mTIP3P behaves essentially as a simple liquid while TIP4P shows the density and related anomalies characteristic of water. We also study the number of water molecules, the tetrahedral order, and the tagged molecule potential energies for water molecules as a function of the distance from the peptide for the 16-residue β-hairpin fragment of 2GB1 in mTIP3P and TIP4P solvents. The hydration shell coordination profiles (n(r)) of the number of water molecules are almost identical in the two solvents, but the radial variation in the local energies and local order show significant differences. The residue-wise variation in the tagged potential energy of water molecules within the first hydration shell is qualitatively similar in the two models. A comparison of the tetrahedral order distributions of water molecules lying at different distances from the biomolecular solute shows that the perturbation in the local tetrahedral order distributions of the bulk solvent due to the presence of the solute is marginal. Thus, in the 250-300 K regime, the mTIP3P and TIP4P water models show qualitatively different behavior in terms of the relationship between tetrahedral order and local energy, but as solvents in the neighborhood of a biomolecular solute, the differences between the two models are only quantitative and not qualitative.

Comparison of hydration behavior and conformational preferences of the Trp-cage mini-protein in different rigid-body water models

The secondary structure conformational properties and hydration shell metrics of the Trp-cage mini-protein are examined in the folded and unfolded ensembles in mTIP3P, TIP4P, and TIP4P-Ew water models with the CHARMM22 force-field using molecular dynamics simulations at 250 K. Upon changing the water model, the conformational order metrics of the peptide show significant differences in the unfolded rather than in the folded ensemble. The unfolding temperatures for Trp-cage are observed to be around 460, 470, and 430 K in mTIP3P, TIP4P, and TIP4P-Ew, respectively. Upon comparing the results with a previous study on a 16-residue β-hairpin fragment of the 2GB1 protein, the same set of conformational order metrics are found to be insufficient in describing the free energy landscape of peptides having a distinct native secondary structure. However, the hydration shell properties of the peptide have been found to be independent of the sequence of the peptide and it changes in conformation upon unfolding. Our calculations reveal that for a particular water model, the secondary structure preferences in the unfolded ensembles of the two peptides are qualitatively different. The unfolded structures of Trp-cage prefer extended and compact structures in TIP4P-Ew and mTIP3P water, respectively, whereas the β-hairpin peptide prefers extended unfolded structures in mTIP3P. The conformational preferences of the unfolded peptide in a given water model have been found to depend on the peptide sequence, where the binding energies of the water molecules around the polar residues in the unfolded conformations show sensitivity to the multipole moments of the water models. The significance of an accurate description of peptide-solvent interactions in the parametrization of biomolecular force-fields, to obtain an accurate description of conformational preferences, in particular in the unfolded ensembles of proteins, is highlighted.

The Hydrophobic Effect and the Role of Cosolvents

Cosolvents modulate aqueous solubility, hydrophobic interactions, and the stability and function of most proteins in the living cell. Our molecular-level understanding of cosolvent effects is incomplete, not only at the level of complex systems such as proteins, but also at the level of very fundamental interactions that underlie the hydrophobic effect. This Feature Article discusses cosolvent effects on the aqueous solubility of nonpolar solutes, hydrophobic interactions, and hydrophobic self-assembly/collapse of aqueous polymers, recently studied with molecular dynamics simulations. It is shown that direct interactions of cosolvents with nonpolar solutes and aqueous polymers can strengthen hydrophobic interactions and can contribute to stabilizing collapsed globular structures. The molecular-level explanation of these observations requires a better understanding of the entropy associated with fluctuations of attractive solute-solvent interactions and of length-scale dependencies of this quantity.

Convergence of Kirkwood–Buff Integrals of Ideal and Nonideal Aqueous Solutions Using Molecular Dynamics Simulations

The computation of Kirkwood-Buff integrals (KBIs) using molecular simulations of closed systems is challenging due to finite system-size effects. One of the problems involves the incorrect asymptotic behavior of the radial distribution function. Corrections to rectify such effects have been proposed in the literature. This study reports a systematic comparison of the proposed corrections (as given by Ganguly et al. J. Chem. Theory Comput. 2013, 9, 1347-1355 and Krüger et al. J. Phys. Chem. Lett. 2013, 4, 4-7) to assess the asymptotic behavior of the RDFs, the KBIs, as well as the estimation of thermodynamic quantities for ideal urea-water and nonideal modified-urea-water mixtures using molecular dynamics simulations. The results show that applying the KBI correction suggested by Krüger et al. on the RDF corrected with the Ganguly et al. correction (denoted as B-KBI) yields improved KBI convergence for the ideal and nonideal aqueous mixtures. Different averaging regions in the running KBIs (correlated or long-range) are assessed, and averaging over the correlated region for large system sizes is found to be robust toward the change in the degree of solvent nonideality and concentration, providing good estimates of thermodynamic quantities. The study provides new insights into improving the KBI convergence, the suitability of different averaging regions in KBIs to estimate thermodynamic properties, as well as the applicability of correction methods to achieve KBI convergence for nonideal aqueous binary mixtures.

Free Energy Landscapes of Alanine Oligopeptides in Rigid-Body and Hybrid Water Models

Replica exchange molecular dynamics is used to study the effect of different rigid-body (mTIP3P, TIP4P, SPC/E) and hybrid (H1.56, H3.00) water models on the conformational free energy landscape of the alanine oligopeptides (acAnme and acA5nme), in conjunction with the CHARMM22 force field. The free energy landscape is mapped out as a function of the Ramachandran angles. In addition, various secondary structure metrics, solvation shell properties, and the number of peptide-solvent hydrogen bonds are monitored. Alanine dipeptide is found to have similar free energy landscapes in different solvent models, an insensitivity which may be due to the absence of possibilities for forming i-(i + 4) or i-(i + 3) intrapeptide hydrogen bonds. The pentapeptide, acA5nme, where there are three intrapeptide backbone hydrogen bonds, shows a conformational free energy landscape with a much greater degree of sensitivity to the choice of solvent model, though the three rigid-body water models differ only quantitatively. The pentapeptide prefers nonhelical, non-native PPII and β-sheet populations as the solvent is changed from SPC/E to the less tetrahedral liquid (H1.56) to an LJ-like liquid (H3.00). The pentapeptide conformational order metrics indicate a preference for open, solvent-exposed, non-native structures in hybrid solvent models at all temperatures of study. The possible correlations between the properties of solvent models and secondary structure preferences of alanine oligopeptides are discussed, and the competition between intrapeptide, peptide-solvent, and solvent-solvent hydrogen bonding is shown to be crucial in the relative free energies of different conformers.

Relating structure, entropy, and energy of solvation of nanoscale solutes: application to gold nanoparticle dispersions.

Structural estimators for the entropy are combined with an analysis of the different contributions to the energy of solvation to understand the molecular basis of the thermodynamics of solvation of passivated nanoparticles. Molecular dynamics simulations of thiolated gold clusters in ethane are performed over a wide range of densities close to the critical isotherm. The entropic changes associated with solvent reorganization around the passivated nanoparticle are estimated from the nanoparticle-solvent pair correlation function, while the entropy of the ligand shell is estimated from the covariance in the positional fluctuations of the ligand atoms. The ligand-shell entropy (S(L)) is shown to be fairly insensitive to variations in solvent density ranging from vacuum to twice the critical density (ρ(c)). In contrast, the entropy change due to solvent reorganization (ΔS(ns)(ord)) shows a minimum around the critical point where the solvent excess shows a maximum. Combining the entropic estimates with the nanoparticle-solvent interaction energies, the free energy of solvation is shown to decrease with density once the critical point is crossed in a manner qualitatively consistent with available experimental data. The results suggest that such an approach to obtain structural insights into the thermodynamics of solvation of passivated nanoparticles could be useful in understanding the stability of nanoparticle dispersions of widely varying chemistries. This study also demonstrates that the theoretical analysis of solvation and self-assembly developed in the context of biomolecular hydration can be very usefully extended to understand the behavior of inorganic nanoparticle dispersions.

Cosolvent Effects on Polymer Hydration Drive Hydrophobic Collapse

Water-mediated hydrophobic interactions play an important role in self-assembly processes, aqueous polymer solubility, and protein folding, to name a few. Cosolvents affect these interactions; however, the implications for hydrophobic polymer collapse and protein folding equilibria are not well-understood. This study examines cosolvent effects on the hydrophobic collapse equilibrium of a generic 32-mer hydrophobic polymer in urea, trimethylamine- N-oxide (TMAO), and acetone aqueous solutions using molecular dynamics simulations. Our results unveil a remarkable cosolvent-concentration-dependent behavior. Urea, TMAO, and acetone all shift the equilibrium toward collapsed structures below 2 M cosolvent concentration and, in turn, to unfolded structures at higher cosolvent concentrations, irrespective of the differences in cosolvent chemistry and the nature of cosolvent-water interactions. We find that weakly attractive polymer-water van der Waals interactions oppose polymer collapse in pure water, corroborating related observations reviewed by Ben-Amotz ( Annu. Rev. Phys. Chem. 2016, 67, 617-638). The cosolvents studied in the present work adsorb at the polymer/water interface and expel water molecules into the bulk, thereby effectively removing the dehydration energy penalty that opposes polymer collapse in pure water. At low cosolvent concentrations, this leads to cosolvent-induced stabilization of collapsed polymer structures. Only at sufficiently high cosolvent concentrations, polymer-cosolvent interactions favor polymer unfolding.

Intrinsic Conformational Preferences and Interactions in α-Synuclein Fibrils: Insights from Molecular Dynamics Simulations

Amyloid formation by the intrinsically disordered α-synuclein protein is the hallmark of Parkinsons disease. We present atomistic Molecular Dynamics simulations of the core of α-synuclein using enhanced sampling techniques to describe the conformational and binding free energy landscapes of fragments implicated in fibril stabilization. The theoretical framework is derived to combine the free energy profiles of the fragments into the reaction free energy of a protein binding to a fibril. Our study shows that individual fragments in solution have a propensity toward attaining non-β conformations, indicating that in a fibril β-strands are stabilized by interactions with other strands. We show that most dimers of hydrogen-bonded fragments are unstable in solution, while hydrogen bonding stabilizes the collective binding of five fragments to the end of a fibril. Hydrophobic effects make further contributions to the stability of fibrils. This study is the first of its kind where structural and binding preferences of the five major fragments of the hydrophobic core of α-synuclein have been investigated. This approach improves sampling of intrinsically disordered proteins, provides information on the binding mechanism between the core sequences of α-synuclein, and enables the parametrization of coarse grained models.