A new angle on heat capacity changes in hydrophobic solvation

Universal cold RNA phase transitions

RNA's diversity of structures and functions impacts all life forms since primordia. We use calorimetric force spectroscopy to investigate RNA folding landscapes in previously unexplored low-temperature conditions. We find that Watson-Crick RNA hairpins, the most basic secondary structure elements, undergo a glass-like transition below [Formula: see text]C where the heat capacity abruptly changes and the RNA folds into a diversity of misfolded structures. We hypothesize that an altered RNA biochemistry, determined by sequence-independent ribose-water interactions, outweighs sequence-dependent base pairing. The ubiquitous ribose-water interactions lead to universal RNA phase transitions below TG, such as maximum stability at [Formula: see text]C where water density is maximum, and cold denaturation at [Formula: see text]C. RNA cold biochemistry may have a profound impact on RNA function and evolution.

Elastin recoil is driven by the hydrophobic effect

Elastin is an extracellular matrix material found in all vertebrates. Its reversible elasticity, robustness, and low stiffness are essential for the function of arteries, lungs, and skin. It is among the most resilient elastic materials known: During a human lifetime, arterial elastin undergoes in excess of 2 × 109 stretching/contracting cycles without replacement, and slow oxidative hardening has been identified as a limiting factor on human lifespan. For over 50 y, the mechanism of entropic recoil has been controversial. Herein, we report a combined NMR and thermomechanical study that establishes the hydrophobic effect as the primary driver of elastin function. Water ordering at the solvent:protein interface was observed as a function of stretch using double quantum 2H NMR, and the most extensive thermodynamic analysis performed to date was obtained by measuring elastin length and volume as a function of force and temperature in normal water, heavy water and with cosolvents. When stretched, elastin's heat capacity increases, water is ordered proportional to the degree of stretching, the internal energy decreases, and heat is released in excess of the work performed. These properties show that recoil in elastin under physiological conditions is primarily driven by the hydrophobic effect rather than by configurational entropy as is the case for rubber. Consistent with this conclusion are decreases in the thermodynamic signatures when cosolvents that alter the hydrophobic effect are introduced. We propose that hydrophobic effect-driven recoil, as opposed to a configurational entropy mechanism where hardening from crystallization can occur, is the origin of elastin's unusual resilience.

The Hydrophobic Effect Studied by Using Interacting Colloidal Suspensions

Interactions between nanoparticles (NPs) determine their self-organization and dynamic processes. In these systems, a quantitative description of the interparticle forces is complicated by the presence of the hydrophobic effect (HE), treatable only qualitatively, and due to the competition between the hydrophobic and hydrophilic forces. Recently, instead, a sort of crossover of HE from hydrophilic to hydrophobic has been experimentally observed on a local scale, by increasing the temperature, in pure confined water and studying the occurrence of this crossover in different water-methanol solutions. Starting from these results, we then considered the idea of studying this process in different nanoparticle solutions. By using photon correlation spectroscopy (PCS) experiments on dendrimer with OH terminal groups (dissolved in water and methanol, respectively), we show the existence of this hydrophobic-hydrophilic crossover with a well defined temperature and nanoparticle volume fraction dependence. In this frame, we have used the mode coupling theory extended model to evaluate the measured time-dependent density correlation functions (ISFs). In this context we will, therefore, show how the measured spectra are strongly dependent on the specificity of the interactions between the particles in solution. The observed transition demonstrates that just the HE, depending sensitively on the system thermodynamics, determines the hydrophobic and hydrophilic interaction properties of the studied nanostructures surface.

Exploring the role of polymer hydrophobicity in polymer-metal binding thermodynamics

Metal-chelating polymers play a key role in rare-earth element (REE) extraction and separation processes. Often, these processes occur in aqueous solution, but the interactions among water, polymer, and REE are largely under-investigated in these applications. To probe these interactions, we synthesized a series of poly(amino acid acrylamide)s with systematically varied hydrophobicity around a consistent chelating group (carboxylate). We then measured the ΔH of Eu³⁺ chelation as a function of temperature across the polymer series using isothermal titration calorimetry (ITC) to give the change in heat capacity (ΔC_P). We observed an order of magnitude variation in ΔC_P (39-471 J mol¹ K^-1) with changes in the hydrophobicity of the polymer. Atomistic simulations of the polymer-metal-water interactions revealed greater Eu³⁺ and polymer desolvation when binding to the more hydrophobic polymers. These combined experimental and computational results demonstrate that metal binding in aqueous solution can be modulated not only by directly modifying the chelating groups, but also by altering the molecular environment around the chelating site, thus suggesting a new design principle for developing increasingly effective metal-chelating materials.

Distinguishing Weak and Strong Hydrogen Bonds in Liquid Water-A Potential of Mean Force-Based Approach

The ability to form hydrogen bonds is one of the most important factors behind water's many anomalous properties. However, there is still no consensus on the hydrogen bond structure of liquid water, including the average number of hydrogen bonds in liquid water. We use molecular dynamics simulations of the polarizable iAMOEBA water model for investigating the hydrogen bond characteristics of liquid water over a wide range of temperatures and pressures. Geometric definitions of a hydrogen bond often use a rectangular region on the plane of hydrogen bond distances and angles. In this work, we find that an elliptical region is more appropriate for the identification of hydrogen bonds, based on statistically favorable molecular configurations. The two-dimensional potential of mean force (PMF) landscape along the hydrogen bond distance (O-H) and angle (O-H-O) is calculated for identifying the statistically favored molecular configurations, which is then used for defining hydrogen bond formation as well as the strength of a hydrogen bond. We further propose a new approach to characterize the hydrogen bonds as strong when the PMF is lower than -2 kT. Using this definition, a consistent explanation for the different average numbers of hydrogen bonds in water is obtained in agreement with the literature. Simulations are also performed with the rigid and nonpolarizable TIP4P/2005 water model. Both water models are qualitatively consistent in predicting the distribution of double-, single-, and non-donor configurations, in line with experimental data, while the iAMOEBA water model yields more quantitatively precise results, including a 10-15% double-donor fraction at 90 °C and 1 atm. The method is also demonstrated to be applicable to the recent, and more general, three-dimensional PMF-based definition of hydrogen bonds.

Binary structure and dynamics of the hydrogen bonds in the hydration shells of ions

Ion-specific effects of cations (Li+, Na+, K+, Mg2+, Ca2+) and anions (F-, Cl-) on the hydrogen bond structure and dynamics of the coordination waters in the hydration shells have been studied using molecular dynamics simulations. Our simulations indicate that the hydrogen bonds between the first and second hydration shell waters show binary structural and dynamic properties. The hydrogen bond with a first shell water as the donor (HD) is strengthened, while those with a first shell water as the acceptor (HA) are weakened. For a hydrated anion, this binary effect reverses, but is less significant. This ion-specific binary effect correlates with the size and the valence of the ion, and is more significant for the strong kosmotropic ions of high charge density.

Uncovering Differences in Hydration Free Energies and Structures for Model Compound Mimics of Charged Side Chains of Amino Acids

Free energies of hydration are of fundamental interest for modeling and understanding conformational and phase equilibria of macromolecular solutes in aqueous phases. Of particular relevance to systems such as intrinsically disordered proteins are the free energies of hydration and hydration structures of model compounds that mimic charged side chains of Arg, Lys, Asp, and Glu. Here, we deploy a Thermodynamic Cycle-based Proton Dissociation (TCPD) approach in conjunction with data from direct measurements to obtain estimates for the free energies of hydration for model compounds that mimic the side chains of Arg+, Lys+, Asp-, and Glu-. Irrespective of the choice made for the hydration free energy of the proton, the TCPD approach reveals clear trends regarding the free energies of hydration for Arg+, Lys+, Asp-, and Glu-. These trends include asymmetries between the hydration free energies of acidic (Asp- and Glu-) and basic (Arg+ and Lys+) residues. Further, the TCPD analysis, which relies on a combination of experimental data, shows that the free energy of hydration of Arg+ is less favorable than that of Lys+. We sought a physical explanation for the TCPD-derived trends in free energies of hydration. To this end, we performed temperature-dependent calculations of free energies of hydration and analyzed hydration structures from simulations that use the polarizable Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field and water model. At 298 K, the AMOEBA model generates estimates of free energies of hydration that are consistent with TCPD values with a free energy of hydration for the proton of ca. -259 kcal/mol. Analysis of temperature-dependent simulations leads to a structural explanation for the observed differences in free energies of hydration of ionizable residues and reveals that the heat capacity of hydration is positive for Arg+ and Lys+ and negative for Asp- and Glu-.

Hydrogen Bonding in the Liquid State of Linear Alcohols: Molecular Dynamics and Thermodynamics

Linear monohydroxy alcohols are strongly hydrogen-bonded liquids that are considered to be homologues of water. Here, we report ab initio molecular dynamics simulations of the liquid alcohols, methanol to pentanol, and from the combined radial-angular probability distribution of the intermolecular O···O distances and HO···O angles determine the geometrical parameters that define the hydrogen bonds in these systems. The key feature of hydrogen bonds in the liquid alcohols, irrespective of the size of the alkyl group, is the strong orientation dependence with the donor-acceptor HO···O angle being close to zero, similar to that observed in liquid water. Hydrogen bond formation is consequently considered to be the passage from a state where donor-acceptor pairs show no preferred orientation to one where they are almost linear. The potential of mean force, the reversible work associated with this process, is computed from the pair probability density distributions obtained from the simulations and that for a hypothetical state where donor-acceptor pairs are randomly oriented. We find that the magnitude of the free energy for hydrogen bond formation is maximum for ethanol and show that this arises from a larger electrostatic contribution to hydrogen bond formation in ethanol as compared to the other alcohols.

An ab initio molecular dynamics study of benzene in water at supercritical conditions: Structure, dynamics, and polarity of hydration shell water and the solute

Anisotropic structure and dynamics of the hydration shell of a benzene solute in supercritical water are investigated by means of ab initio molecular dynamics simulations. The polarity and structural distortion of the benzene solute in supercritical water are also investigated in this study. Calculations are done at 673 K for three different densities of the solvent. The simulations are carried out using the Becke-Lee-Yang-Parr (BLYP) and also the Becke-Lee-Yang-Parr functional including dispersion corrections of Grimme (BYLP-D). The structural anisotropy is found to exist even at supercritical conditions as elucidated by the radial distribution functions of different conical regions and also by angular and spatial distribution functions. The benzene-water πH-bond and also the water-water hydrogen bonds are found to exist even at the supercritical temperature of 673 K. However, the numbers of these hydrogen bonds are reduced substantially with a decrease in water density. The water molecules in the axial region of benzene are found to be preferably oriented with one OH vector pointing toward the benzene ring, whereas the water molecules located in the equatorial region are found to orient their dipoles mostly parallel to the ring plane. The orientational distributions, however, are found to be rather broad at the supercritical temperature due to thermal fluctuations. Although the water molecules have faster dynamics at these supercritical conditions, a slight difference is observed in the dynamics of the solvation shell and bulk molecules. The conformational flexibility of the ring is found to be enhanced which causes an increase in polarity of the benzene solute in water under supercritical conditions.

Interface-inspired formulation and molecular-level perspectives on heat conduction and energy storage of nanofluids

Aiming for the introduction of stability requirements in nanofluids processing, an interface-based three-step method is proposed in this work. It is theory-based design framework for nanofluids that aims for a minimum tension at the solid-liquid interface by adjusting the polar and dispersive components of the base fluid to meet those of disperse nanomaterial. The method was successfully tested in the preparation of aqueous nanofluids containing single-walled carbon nanotubes that resulted to be stable and to provide good thermal properties, i.e. thermal conductivity increases by 79.5% and isobaric specific heat by 8.6% for a 0.087 vol.% load of nanotubes at 70 °C. Besides, a system for these nanofluids was modelled. It was found to be thermodynamically consistent and computationally efficient, providing consistent response to changes in the state variable temperature in a classical Molecular Dynamics environment. From an analysis of the spatial components of the heat flux autocorrelation function, using the equilibrium approach, it was possible to elucidate that heat conduction through the host fluid is enhanced by phonon propagation along nanotubes longitudinal axes. From an analysis of the structural features described by radial distribution functions, it was concluded that additional heat storage arises from the hydrophobic effect.

Thermodynamic Stability of DNA Duplexes Comprising the Simplest T → dU Substitutions

Members of the uracil-DNA glycosylase (UDG) enzyme family recognize and bind uracil, sequestering it within the binding site pocket and catalyzing the cleavage of the base from the deoxyribose, leaving an abasic site. The recognition and binding are passive and rely on innate dynamic motions of DNA wherein base pairs undergo thermally induced breakage and conformational fluctuations. Once the uracil breaks from its base pair, it can be recognized and bound by the enzyme, which then alters its conformation for sequestration and catalysis. Our results suggest that the thymine to uracil substitution, which differs only by a single methyl group, causes a destabilization of the duplex thermodynamics, which would lead to an increase in the population of the extrahelical state and increase the probability of uracil being recognized and excised from DNA by UDG. This destabilization is dependent on the identity of the nearest-neighbor base-pair stacks; a G·C nearest neighbor leads to thermal and enthalpic destabilization that is weaker that that seen with two A·T neighbors. In addition, uracil substitution yields a nearest-neighbor increase in the counterion uptake of the duplexes but decreases the level of immobilization of structural water for all substituted duplexes regardless of the neighbor identity or number of substitutions.

Energetics, Ion, and Water Binding of the Unfolding of AA/UU Base Pair Stacks and UAU/UAU Base Triplet Stacks in RNA

Triplex formation occurs via interaction of a third strand with the major groove of double-stranded nucleic acid, through Hoogsteen hydrogen bonding. In this work, we use a combination of temperature-dependent UV spectroscopy and differential scanning calorimetry to determine complete thermodynamic profiles for the unfolding of polyadenylic acid (poly(rA))·polyuridylic acid (poly(rU)) (duplex) and poly(rA)·2poly(rU) (triplex). Our thermodynamic results are in good agreement with the much earlier work of Krakauer and Sturtevant using only UV melting techniques. The folding of these two helices yielded an uptake of ions, Δ n_Na⁺ = 0.15 mol Na⁺/mol base pair (duplex) and 0.30 mol Na⁺/mole base triplet (triplex), which are consistent with their polymer behavior and the higher charge density parameter of triple helices. The osmotic stress technique yielded a release of structural water, Δ n_W = 2 mol H₂O/mol base pair (duplex unfolding into single strands) and an uptake of structural water, Δ n_W = 2 mol H₂O/mole base pair (triplex unfolding into duplex and a single strand). However, an overall release of electrostricted waters is obtained for the unfolding of both complexes from pressure perturbation calorimetric experiments. In total, the Δ V values obtained for the unfolding of triplex into duplex and a single strand correspond to an immobilization of two structural waters and a release of three electrostricted waters. The Δ V values obtained for the unfolding of duplex into two single strands correspond to the release of two structural waters and the immobilization of four electrostricted water molecules.

Dynamics of water in conical solvation shells around a benzene solute under different thermodynamic conditions

Water molecules in different parts of the anisotropic hydration shell of an aromatic molecule experience different interactions. In the present study, we investigate the anisotropic dynamics of water molecules in different non-overlapping conical shells around a benzene solute at sub- and supercritical conditions by means of molecular dynamics simulations using both non-polarizable and polarizable models. In addition to the dynamical properties, the effects of polarizability on the hydration structure of benzene at varying thermodynamic conditions are also investigated in the current study. The presence of πH-bonding in the solvation shell is found to be an important factor that influences the anisotropic dynamics of the benzene hydration shell. The water molecules located axial to the benzene plane are found to be maximally influenced by the πH-bonding. The extent of πH-bonding is found to be somewhat reduced on inclusion of polarizability. The πH-bonded water molecules are found to reorient through large-amplitude angular jumps where the jump-angle amplitude increases at higher temperatures and lower densities. For both non-polarizable and polarizable models, it is found that the water molecules in the axial conical shells possess faster orientational and hydrogen bond dynamics compared to those in the equatorial plane. Water molecules in the axial conical shells are also found to diffuse at a faster rate than bulk molecules due to the relatively weaker benzene-water πH-bonding interactions in the axial region of the hydration shell. The residence dynamics of water molecules in different conical solvation shells around the solute is also investigated in the current study.

Increased Flexibility between Stems of Intramolecular Three-Way Junctions by the Insertion of Bulges

Intramolecular junctions are a ubiquitous structure within DNA and RNA; three-way junctions in particular have high strain around the junction because of the lack of flexibility, preventing the junctions from adopting conformations that would allow for optimal folding. In this work, we used a combination of calorimetric and spectroscopic techniques to study the unfolding of four intramolecular three-way junctions. The control three-way junction, 3H, has the sequence d(GAAATTGCGCT5GCGCGTGCT5GCACAATTTC), which has three arms of different sequences. We studied three other three-way junctions in which one (2HS1H), two (HS12HS1), and three (HS1HS1HS1) cytosine bulges were placed at the junction to allow the arms to adopt a wider range of conformations that may potentially relieve strain. Through calorimetric studies, it was concluded that bulges produce only minor effects on the enthalpic and thermal stability at physiological salt concentrations for 2HS1H and HS1HS1HS1. HS12HS1 displays the strongest effect, with the GTGC stem lacking a defined transition. In addition to unfolding thermodynamics, the differential binding of counterions, water, and protons was determined. It was found that with each bulge, there was a large increase in the binding of counterions; this correlated with a decrease in the immobilization of structural water molecules. The increase in counterion uptake upon folding likely displaces binding of structural water, which is measured by the osmotic stress method, in favor of electrostricted waters. The cytosine bulges do not affect the binding of protons; this finding indicates that the bulges are not forming base-triplet stacks. These results indicate that bulges in junctions do not affect the unfolding profile or the enthalpy of oligonucleotides but do affect the number and amount of molecules immobilized by the junction.

Local chemistry of the surfactant's head groups determines protein stability in reverse micelles

Protein stability in reverse micelles is determined by local chemical interactions between the surfactant molecules and the protein groups.

Mapping Hydration Water around Alcohol Chains by THz Calorimetry

THz spectroscopy was used to probe changes that occur in the dynamics of the hydrogen bond network upon solvation of alcohol chains. The THz spectra can be decomposed into the spectrum of bulk water, tetrahedral hydration water, and more disordered (or interstitial) hydration water. The tetrahedrally ordered hydration water exhibits a band at 195 cm⁻¹ and is localized around the hydrophobic moiety of the alcohol. The interstitial component yields a band at 164 cm⁻¹ which is associated with hydration water in the first hydration shell. These temperature‐dependent changes in the low‐frequency spectrum of solvated alcohol chains can be correlated with changes of heat capacity, entropy, and free energy upon solvation. Surprisingly, not the tetrahedrally ordered component but the interstitial hydration water is found to be mainly responsible for the temperature‐dependent change in ΔC_p and ΔG. The solute‐specific offset in free energy is attributed to void formation and scales linearly with the chain length.

Interfacial Water Arrangement in the Ice-Bound State of an Antifreeze Protein: A Molecular Dynamics Simulation Study

Molecular dynamics (MD) simulations have been carried out to study the heterogeneous ice nucleation on modeled peptide surfaces. Simulations show that large peptide surfaces made by TxT (threonine-x-threonine) motifs with the arrangements of threonine (Thr) residues identical to the periodic arrangements of waters on either the basal or prism plane of ice are capable of ice nucleation. Nucleated ice plane is the (0001) basal plane of hexagonal ice (Ih) or (111) plane of cubic ice (Ic). However, due to predefined simulation cell dimensions, the ice growth is only observed on the surface where the Thr residues are arranged like the water arrangement on the basal plane of ice Ih. The γ-methyl and γ-hydroxyl groups of Thr residue are necessary for such ice formation. From this ice nucleation and growth simulation, the interfacial water arrangement in the ice-bound state of Tenebrio molitor antifreeze protein (TmAFP) has been determined. The interfacial water arrangement in the ice-bound state of TmAFP is characterized by five-membered hydrogen bonded rings, where each of the hydroxyl groups of the Thr residues on the ice-binding surface (IBS) of the protein is a ring member. It is found that the water arrangement at the protein-ice interface is distorted from that in bulk ice. Our analysis further reveals that the hydroxyl groups of Thr residues on the IBS of TmAFP form maximum three hydrogen bonds each with the waters in the bound state and methyl groups of Thr residues occupy wider spaces than the normal grooves on the (111) plane of ice Ic. Methyl groups are also located above and along the 3-fold rotational axes of the chair-formed hexagonal hydrogen bonded water rings on the (111) plane.

Interaction of ice binding proteins with ice, water and ions

Ice binding proteins (IBPs) are produced by various cold-adapted organisms to protect their body tissues against freeze damage. First discovered in Antarctic fish living in shallow waters, IBPs were later found in insects, microorganisms, and plants. Despite great structural diversity, all IBPs adhere to growing ice crystals, which is essential for their extensive repertoire of biological functions. Some IBPs maintain liquid inclusions within ice or inhibit recrystallization of ice, while other types suppress freezing by blocking further ice growth. In contrast, ice nucleating proteins stimulate ice nucleation just below 0 °C. Despite huge commercial interest and major scientific breakthroughs, the precise working mechanism of IBPs has not yet been unraveled. In this review, the authors outline the state-of-the-art in experimental and theoretical IBP research and discuss future scientific challenges. The interaction of IBPs with ice, water and ions is examined, focusing in particular on ice growth inhibition mechanisms.

Characterization of Cetyltrimethylammonium Bromide/Hexanol Reverse Micelles by Experimentally Benchmarked Molecular Dynamics Simulations

Encapsulation of small molecules, proteins, and other macromolecules within the protective water core of reverse micelles is emerging as a powerful strategy for a variety of applications. The cationic surfactant cetyltrimethylammonium bromide (CTAB) in combination with hexanol as a cosurfactant is particularly useful in the context of solution NMR spectroscopy of encapsulated proteins. Small-angle X-ray and neutron scattering is employed to investigate the internal structure of the CTAB/hexanol reverse micelle particle under conditions appropriate for high-resolution NMR spectroscopy. The scattering profiles are used to benchmark extensive molecular dynamics simulations of this reverse micelle system and indicate that the parameters used in these simulations recapitulate experimental results. Scattering profiles and simulations indicate formation of homogeneous solutions of small approximately spherical reverse micelle particles at a water loading of 20 composed of ∼150 CTAB and 240 hexanol molecules. The 3000 waters comprising the reverse micelle core show a gradient of translational diffusion that reaches that of bulk water at the center. Rotational diffusion is slowed relative to bulk throughout the water core, with the greatest slowing near the CTAB headgroups. The 5 Å thick interfacial region of the micelle consists of overlapping layers of Br(-) enriched water, CTAB headgroups, and hexanol hydroxyl groups, containing about one-third of the total water. This study employs well-parametrized MD simulations, X-ray and neutron scattering, and electrostatic theory to illuminate fundamental properties of CTAB/hexanol reverse micelle size, shape, partitioning, and water behavior.

Anisotropic structure and dynamics of the solvation shell of a benzene solute in liquid water from ab initio molecular dynamics simulations

The anisotropic structure and dynamics of the hydration shell of a benzene solute in liquid water have been investigated by means of ab initio molecular dynamics simulations using the BLYP (Becke–Lee–Yang–Parr) and dispersion corrected BLYP-D functionals.

Computational Free Energy Studies of a New Ice Polymorph Which Exhibits Greater Stability than Ice Ih

The absolute free energies of several ice polymorphs were calculated using thermodynamic integration. These polymorphs are predicted by computer simulations using a variety of common water models to be stable at low pressures. A recently discovered ice polymorph that has as yet only been observed in computer simulations (Ice-i) was determined to be the stable crystalline state for all the water models investigated. Phase diagrams were generated, and phase coexistence lines were determined for all of the known low-pressure ice structures. Additionally, potential truncation was shown to play a role in the resulting shape of the free energy landscape.

The remarkable hydration of the antifreeze protein Maxi: a computational study

The long four-helix bundle antifreeze protein Maxi contains an unusual core for a globular protein. More than 400 ordered waters between the helices form a nano-pore of internal water about 150 Å long. Molecular dynamics simulations of hydrated Maxi were carried out using the CHARMM27 protein forcefield and the TIP3P water model. Solvation in the core and non-core first hydration shell was analyzed in terms of water-water H-bond distance-angle distributions. The core had an increased population of low-angle H-bonds between water pairs relative to bulk water. Enhancement of low angle H-bonds was particularly pronounced for water pairs at the interfaces between apolar and polar regions inside the protein core, characteristic of the anchored clathrate solvation structure seen previously in the ice-nuclei binding surfaces of type I, type III, and beta-helical antifreeze proteins. Anchored clathrate type solvation structure was not seen in the exterior solvation shell except around residues implicated in ice binding. Analysis of solvation dynamics using water residence times and diffusion constants showed that exterior solvation shell waters exchanged rapidly with bulk water, with no difference between ice-binding and non-binding residues. Core waters had about ten-fold slower diffusion than bulk water. Water residence times around core residues averaged about 8 ps, similar to those on the exterior surface, but they tended to exchange primarily with other core water, resulting in longer, >40 ps residence times within the core. Preferential exchange or diffusion of the water along the long axis of the water core of Maxi was not detected.

Why and How Does Pressure Unfold Proteins?

This year, 2014, marks the 100th anniversary of the first publication reporting the denaturation of proteins by high hydrostatic pressure (Bridgman 1914). Since that time a large and recently increasing number of studies of pressure effects on protein stability have been published, yet the mechanism for the action of pressure on proteins remains subject to considerable debate. This review will present an overview from this author's perspective of where this debate stands today.

Trehalose facilitates DNA melting: a single-molecule optical tweezers study

Using optical tweezers, here we show that the overstretching transition force of double-stranded DNA (dsDNA) is lowered significantly by the addition of the disaccharide trehalose as well as certain polyol osmolytes. This effect is found to depend linearly on the logarithm of the trehalose concentration. We propose an entropic driving mechanism for the experimentally observed destabilization of dsDNA that is rooted in the higher affinity of the DNA bases for trehalose than for water, which promotes base exposure and DNA melting. Molecular dynamics simulation reveals the direct interaction of trehalose with nucleobases. Experiments with other osmolytes confirm that the extent of dsDNA destabilization is governed by the ratio between polar and apolar fractions of an osmolyte.

Thermodynamic profiles at the solvated inorganic-organic interface: the case of gold-thiolate monolayers

The thermodynamic adsorption profile at a solvated organic-inorganic interface is probed by following the binding and organization of carboxylic acid-terminated alkanethiols of varying chain lengths (C2, C3, and C6) to the surface of gold nanoparticles (NPs) (5.4 ± 0.7, 9.5 ± 0.6, and 19.4 ± 1.1 nm diameter) using isothermal titration calorimetry (ITC). We discuss the effect of alkyl chain length, temperature, and Au NP size on the energetics at an organic-inorganic interface. ITC allows for the quantification of the adsorption constant, enthalpy of adsorption, entropy of adsorption, and the binding stoichiometry in a single experiment. The thermodynamic parameters support a mechanism of stepwise adsorption of thiols to the surface of Au NPs and secondary ordering of the thiols at the organic-inorganic interface. The adsorption enthalpies are chain-length dependent; enthalpy becomes more exothermic as longer chains are confined, compensating for greater decreases in entropy with increasing chain length. We observe an apparent compensation effect: the negative ΔH is compensated by a negative ΔS as the thiols self-assemble on the Au NP surface. A comparison of the thermodynamic parameters indicates thiol-Au NP association is enthalpy-driven because of the large, exothermic enthalpies accompanied by an unfavorable entropic contribution associated with confinement of alkyl chains, reduced trans-gauche interconversion, and the apparent ordering of solvent molecules around the hydrophobic organic thiols (hydrophobic effect). Understanding the thermodynamics of adsorption at NP surfaces will provide critical insight into the role of ligands in directing size and shape during NP synthesis since thiols are a common ligand choice (i.e., Brust method). The ITC technique is applicable to a larger number of structure-directing ligands and solvent combinations and therefore should become an important tool for understanding reaction mechanisms in nanostructure synthesis.

Solvation thermodynamics and heat capacity of polar and charged solutes in water

The solvation thermodynamics and in particular the solvation heat capacity of polar and charged solutes in water is studied using atomistic molecular dynamics simulations. As ionic solutes we consider a F(-) and a Na(+) ion, as an example for a polar molecule with vanishing net charge we take a SPC/E water molecule. The partial charges of all three solutes are varied in a wide range by a scaling factor. Using a recently introduced method for the accurate determination of the solvation free energy of polar solutes, we determine the free energy, entropy, enthalpy, and heat capacity of the three different solutes as a function of temperature and partial solute charge. We find that the sum of the solvation heat capacities of the Na(+) and F(-) ions is negative, in agreement with experimental observations, but our results uncover a pronounced difference in the heat capacity between positively and negatively charged groups. While the solvation heat capacity ΔC(p) stays positive and even increases slightly upon charging the Na(+) ion, it decreases upon charging the F(-) ion and becomes negative beyond an ion charge of q = -0.3e. On the other hand, the heat capacity of the overall charge-neutral polar solute derived from a SPC/E water molecule is positive for all charge scaling factors considered by us. This means that the heat capacity of a wide class of polar solutes with vanishing net charge is positive. The common ascription of negative heat capacities to polar chemical groups might arise from the neglect of non-additive interaction effects between polar and apolar groups. The reason behind this non-additivity is suggested to be related to the second solvation shell that significantly affects the solvation thermodynamics and due to its large spatial extent induces quite long-ranged interactions between solvated molecular parts and groups.