Heterogeneous preferential solvation of water and trifluoroethanol in homologous lysozymes

Nanoclustering in non-ideal ethanol/heptane solutions alters solvation dynamics

Alcohol/alkane solutions widely used in chemical synthesis and as transportation fuels are highly non-ideal due to the nanoscale clustering of the amphiphilic alcohol molecules within the nonpolar alkanes. Besides impacting reactivity, such as combustion, non-ideal solutions are likely to exhibit unusual solvation dynamics on ultrafast time scales arising from the structurally heterogeneous nature of molecular-scale association. Using a convenient transition metal carbonyl vibrational probe [(C5H5)Mn(CO)3, CMT], linear absorption and nonlinear two-dimensional infrared (2D-IR) spectroscopy reveal composition-dependent solvation dynamics as reported by the frequency fluctuation correlation function in a series of ethanol/heptane solutions. Slow spectral diffusion with dilute ethanol indicates preferential solvation of the polar solute by the alcohol with a mechanism largely dominated by solvent exchange. Comparison with an ethanol/acetonitrile solution series yields no substantial preferential solvation or solvent exchange signatures in the linear or 2D-IR spectra. In ethanol/heptane solutions, increasing the ethanol concentration speeds up the solvation dynamics, which is largely consistent with a model that includes solvent exchange and single-solvent spectral diffusion. Detailed analysis of the deviation from the experimental time constants from the model's optimal parameters yields a remarkable resemblance of the concentration-weighted Kirkwood-Buff integrals for ethanol/heptane solutions. This trend indicates that solution non-ideality alters the spectral diffusion dynamics of the probe solute. Given that nanoscale clustering drives the non-ideality, these experiments reveal a dynamical consequence of nanoscale heterogeneity on the ultrafast dynamics of the solution. Refined understanding of the structural and dynamical aspects of mixed solvents will be necessary for predictive solution strategies in chemistry.

Alcohol-Induced Denaturation of Hen Egg White Lysozyme Studied by Infrared, Circular Dichroism, and Small-Angle Neutron Scattering

In aqueous binary solvents with fluorinated alcohols, 2,2,2-trifluoroethanol (TFE) and 1,1,1,3,3,3-hexafluoroisopropanol (HFIP), and aliphatic alcohols, ethanol (EtOH) and 2-propanol (2-PrOH), the denaturation of hen egg white lysozyme (HEWL) with increasing alcohol mole fraction xA has been investigated in a wide view from the molecular vibration to the secondary and ternary structures. Circular dichroism (CD) measurement showed that the secondary structure of α-helix content of HEWL increases on adding a small amount of the fluorinated alcohol to the aqueous solution, while the β-sheet content decreases. On the contrary, the secondary structure does not significantly change by the addition of the aliphatic alcohols. Correspondingly, the infrared (IR) spectroscopic measurements revealed that the amide I band red-shifts on the addition of the fluorinated alcohol. However, the band remains unchanged in the aliphatic alcohol systems with increasing alcohol content. To observe the ternary structure of HEWL, small-angle neutron scattering (SANS) experiments with H/D substitution technique have been applied to the HEWL solutions. The SANS experiments were successful in revealing the details of how the geometry of the HEWL changes as a function of xA. The SANS profiles indicated the spherical structure of HEWL in all of the alcohol systems in the xA range examined. The mean radius of HEWL in the two fluorinated alcohol systems increases from ∼16 to ∼18 Å during the change in the secondary structure against the increase in the fluorinated alcohol content. On contrast, the radius does not significantly change in both aliphatic alcohol systems below xA = 0.3 but expands to ∼19 Å as the alcohol content is close to the limitation of the HEWL solubility. According to the present results, together with our knowledge of the alcohol cluster formation and the interaction of the trifluoromethyl (CF3) groups with the hydrophobic moieties of biomolecules, the effects of alcohols on the denaturation of the protein have been discussed on a molecular scale.

Examination of Solvent Interactions with Trp-Cage in 1,1,1,3,3,3-Hexafluoro-2-propanol-water at 298 K through MD Simulations and Intermolecular Nuclear Overhauser Effects

MD simulations of the peptide Trp-cage dissolved in 28% hexafluoroisopropanol (HFIP)-water have been carried out at 298 K with the goal of exploring peptide hydrogen-solvent fluorine nuclear spin cross-relaxation. The work was motivated by the observation that most experimental fluoroalcohol-peptide cross-relaxation terms at 298 K are small, both positive and negative, and not always well predicted from simulations. The cross-relaxation terms for hydrogens of the caged tryptophan residue of Trp-cage are substantially negative, a result consistent with simulations. It was concluded that hexafluoroisopropanol interactions near this part of the peptide are particularly long-lived. While both HFIP and water are present in all regions of the simulation box, the composition of the solvent mixture is not homogeneous throughout the system. HFIP generally accumulates near the peptide surface, while water molecules are preferentially found in regions that are more than 1.5 nm from the surface of the peptide. However, some water remains in higher-than-expected amounts in the solvent layer surrounding 6Trp, 9Asp, Ser13, and Ser14 residues in the helical region of Trp-cage. As observed in simulations of this system at 278 K, HFIP molecules aggregate into clusters that continually form and re-form. Translational diffusion of both HFIP and water appears to be slowed near the surface of the peptide with reduction in diffusion near the 6Trp residue 2- to 3-fold larger than calculated for solvent interactions with other regions of Trp-cage.

Molecular Picture of the Effect of Cosolvent Crowding on Ligand Binding and Dispersed Solvation Dynamics in G-Quadruplex DNA

Understanding molecular interactions and dynamics of proteins and DNA in a cell-like crowded environment is crucial for predicting their functions within the cell. Noncanonical G-quadruplex DNA (GqDNA) structures adopt various topologies that were shown to be strongly affected by molecular crowding. However, it is unknown how such crowding affects the solvation dynamics in GqDNA. Here, we study the effect of cosolvent (acetonitrile) crowding on ligand (DAPI) solvation dynamics within human telomeric antiparallel GqDNA through direct comparison of time-resolved fluorescence Stokes shift (TRFSS) experiments and molecular dynamics (MD) simulations results. We show that ligand binding affinity to GqDNA is drastically affected by acetonitrile (ACN). Solvation dynamics probed by DAPI in GqDNA groove show dispersed dynamics from ∼100 fs to 10 ns in the absence and presence of 20% and 40% (v/v) ACN. The nature of dynamics remain similar in buffer and 20% ACN, although in 40% ACN, distinct dynamics is observed in <100 ps. MD simulations performed on GqDNA/DAPI complex reveal preferential solvation of ligand by ACN, particularly in 40% ACN. Simulated solvation time-correlation functions calculated from MD trajectories compare very well to the overall solvation dynamics of DAPI in GqDNA, observed in experiments. Linear response decomposition of simulated solvation correlation functions unfolds the origin of dispersed dynamics, showing that the slower dynamics is dominated by DNA-motion in the presence of ACN (and also by the ACN dynamics at higher concentration). However, water-DNA coupled motion controls the slow dynamics in the absence of ACN. Our data, thus, unravel a detailed molecular picture showing that though ACN crowding affect ligand binding affinity to GqDNA significantly, the overall dispersed solvation dynamics in GqDNA remain similar in the absence and the presence of 20% ACN, albeit with a small effect on the dynamics in the presence of 40% ACN due to preferential solvation of ligand by ACN.

Two different regimes in alcohol-induced coil-helix transition: effects of 2,2,2-trifluoroethanol on proteins being either independent of or enhanced by solvent structural fluctuations

Inhomogeneous distribution of constituent molecules in a mixed solvent has been known to give remarkable effects on the solute, e.g., conformational changes of biomolecules in an alcohol-water mixture. We investigated the general effects of 2,2,2-trifluoroethanol (TFE) on proteins/peptides in a mixture of water and TFE using melittin as a model protein. Fluctuations and Kirkwood-Buff integrals (KBIs) in the TFE-H2O mixture, quantitative descriptions of inhomogeneity, were determined by small-angle X-ray scattering investigation and compared with those in the aqueous solutions of other alcohols. The concentration fluctuation for the mixtures ranks as methanol < ethanol ≪ TFE < tert-butanol < 1-propanol, indicating that the inhomogeneity of molecular distribution in the TFE-H2O mixture is unexpectedly comparable to those in the series of mono-ols. On the basis of the concentration dependence of KBIs between the TFE molecules, it was found that a strong attraction between the TFE molecules is not necessarily important to induce helix conformation, which is inconsistent with the previously proposed mechanism. To address this issue, by combining the KBIs and the helix contents reported by the experimental spectroscopic studies, we quantitatively evaluated the change in the preferential binding parameter of TFE to melittin attributed to the coil-helix transition. As a result, we found two different regimes on TFE-induced helix formation. In the dilute concentration region of TFE below ∼2 M, where the TFE molecules are not aggregated among themselves, the excess preferential binding of TFE to the helix occurs due to the direct interaction between them, namely independent of the solvent fluctuation. In the higher concentration region above ∼2 M, in addition to the former effect, the excess preferential binding is significantly enhanced by the solvent fluctuation. This scheme should be held as general cosolvent effects of TFE on proteins/peptides.

Protein Preferential Solvation in Water:Glycerol Mixtures

For proteins in solvent mixtures, the relative abundances of each solvent in their solvation shell have a critical impact on their properties. Preferential solvation of a series of proteins in water-glycerol mixtures is studied here over a broad range of solvent compositions via classical molecular dynamics simulations. Our simulation results reveal that the differences between shell and bulk compositions exhibit dramatic changes with solvent composition, temperature, and protein nature. In contrast with the simple and widely used picture where glycerol is completely excluded from the protein interface, we show that for aqueous solutions with less than 50% glycerol in volume, protein solvation shells have approximately the same composition as the bulk solvent and proteins are in direct contact with glycerol. We further demonstrate that at high glycerol concentration, glycerol depletion from the solvation shell is due to an entropic factor arising from the reduced accessibility of bulky glycerol molecules in protein cavities. The resulting molecular picture is important to understand protein activity and cryopreservation in mixed aqueous solvents.

Examination of Trifluoroethanol Interactions with Trp-Cage in Trifluoroethanol-Water at 298 K through Molecular Dynamics Simulations and Intermolecular Nuclear Overhauser Effects

Molecular dynamics simulations of the protein model Trp-cage in 42% trifluoroethanol (TFE)-water at 298 K have been carried out with the goal of exploring peptide hydrogen-solvent fluorine nuclear spin cross-relaxation. The TFE5 model of TFE developed in a previous work was used with the TIP5P-Ew model of water. System densities and component translational diffusion coefficients predicted by the simulations were within 20% of the experimental values. Consideration of the calculated relative amounts of TFE and water surrounding the hydrogens of Trp-cage indicated that the composition of the solvent mixture beyond ∼1.5 nm from the van der Waals surface of the peptide is close to the composition of the bulk solvent, but as observed by others, TFE accumulates preferentially near the peptide surface. In the simulations, both TFE and water molecules make contacts with the peptide surface; water molecules predominate in contacts with the peptide backbone atoms and TFE molecules generally preferentially interact with side chains. Translational diffusion of solvent molecules appears to be slowed near the surface of the peptide. Depending on the location in the structure, TFE molecules form complexes with the peptide that may persist for up to ∼7 ns. Many of the peptide spin-solvent fluorine cross-relaxation parameters (Σ_HF) for which experimental values are available are reasonably well-predicted from the simulations. However, the calculated Σ_HF values were too small for some hydrogens of the 6Trp indole ring and the amino acid hydrogens near this residue in the native structure, whereas Σ_HF values for hydrogens on the side chains of 1Asn, 4Ile, and 7Leu are too large. In 42% TFE-water, persistent conformations of Trp-cage are found, which differ from the conformation found in water by the orientation of the 3Tyr ring.

Heterogeneity of water structure and dynamics at the protein-water interface

In this molecular dynamics simulation study, we analyze the local structural and dynamic properties of water hydrating the protein ubiquitin on a spatial grid with 1 Å resolution. This allows for insights into the spatial distribution of water number densities, molecular orientations, translations, and rotations as a function of distance from the protein surface. Water molecule orientations follow a heterogeneous distribution with preferred local orientations of water dipoles and O-H bond vectors up to 10-15 Å distances from the protein, while local variations of the water number density converge to homogeneous bulk-like values within less than 8 Å. Interestingly, we find that the long-ranged orientational structure of water does not impact either the translational or rotational dynamics of water. Instead, heterogeneous distributions of local dynamical parameters and averaged dynamical retardation factors are only found close to the protein surface and follow a distance dependence comparable to heterogeneities in the local water number density. This study shows that the formation of nanodomains of preferred water orientations far from the protein does not significantly impact dynamical processes probed as a non-local average in most experiments.

Spatially Heterogeneous Surface Water Diffusivity around Structured Protein Surfaces at Equilibrium

Hydration water on the surface of a protein is thought to mediate the thermodynamics of protein-ligand interactions. For hydration water to play a role beyond modulating global protein solubility or stability, the thermodynamic properties of hydration water must reflect on the properties of the heterogeneous protein surface and thus spatially vary over the protein surface. A potent read-out of local variations in thermodynamic properties of hydration water is its equilibrium dynamics spanning picosecond to nanosecond time scales. In this study, we employ Overhauser dynamic nuclear polarization (ODNP) to probe the equilibrium hydration water dynamics at select sites on the surface of Chemotaxis Y (CheY) in dilute solution. ODNP reports on site-specific hydration water dynamics within 5-10 Å of a label tethered to the biomolecular surface on two separate time scales of motion, corresponding to diffusive water (DW) and protein-water coupled motions, referred to as bound water (BW). We find DW dynamics to be highly heterogeneous across the surface of CheY. We identify a significant correlation between DW dynamics and the local hydropathy of the CheY protein surface, empirically determined by molecular dynamics (MD) simulations, and find the more hydrophobic sites to be hydrated with slower diffusing water. Furthermore, we compare the hydration water dynamics on different polypeptides and liposome surfaces and find the DW dynamics on globular proteins to be significantly more heterogeneous than on intrinsically disordered proteins (IDPs), peptides, and liposomes. The heterogeneity in the hydration water dynamics suggests that structured proteins have the capacity to encode information into the surrounding hydration shell.

Parallelization and improvements of the generalized born model with a simple sWitching function for modern graphics processors

Two fundamental challenges of simulating biologically relevant systems are the rapid calculation of the energy of solvation and the trajectory length of a given simulation. The Generalized Born model with a Simple sWitching function (GBSW) addresses these issues by using an efficient approximation of Poisson-Boltzmann (PB) theory to calculate each solute atom's free energy of solvation, the gradient of this potential, and the subsequent forces of solvation without the need for explicit solvent molecules. This study presents a parallel refactoring of the original GBSW algorithm and its implementation on newly available, low cost graphics chips with thousands of processing cores. Depending on the system size and nonbonded force cutoffs, the new GBSW algorithm offers speed increases of between one and two orders of magnitude over previous implementations while maintaining similar levels of accuracy. We find that much of the algorithm scales linearly with an increase of system size, which makes this water model cost effective for solvating large systems. Additionally, we utilize our GPU-accelerated GBSW model to fold the model system chignolin, and in doing so we demonstrate that these speed enhancements now make accessible folding studies of peptides and potentially small proteins.