Dynamic hydration shell restores Kauzmann's 1959 explanation of how the hydrophobic factor drives protein folding

NMR Dynamic View of the Destabilization of WW4 Domain by Chaotropic GdmCl and NaSCN

GdmCl and NaSCN are two strong chaotropic salts commonly used in protein folding and stability studies, but their microscopic mechanisms remain enigmatic. Here, by CD and NMR, we investigated their effects on conformations, stability, binding and backbone dynamics on ps-ns and µs-ms time scales of a 39-residue but well-folded WW4 domain at salt concentrations ≤200 mM. Up to 200 mM, both denaturants did not alter the tertiary packing of WW4, but GdmCl exerted more severe destabilization than NaSCN. Intriguingly, GdmCl had only weak binding to amide protons, while NaSCN showed extensive binding to both hydrophobic side chains and amide protons. Neither denaturant significantly affected the overall ps-ns backbone dynamics, but they distinctively altered µs-ms backbone dynamics. This study unveils that GdmCl and NaSCN destabilize a protein before the global unfolding occurs with differential binding properties and µs-ms backbone dynamics, implying the absence of a simple correlation between thermodynamic stability and backbone dynamics of WW4 at both ps-ns and µs-ms time scales.

Weakly Hydrated Solute of Mixed Hydrophobic-Hydrophilic Nature

Infrared (IR) spectroscopy is a commonly used and invaluable tool in studies of solvation phenomena in aqueous solutions. Concurrently, density functional theory calculations and ab initio molecular dynamics simulations deliver the solvation shell picture at the molecular detail level. The mentioned techniques allowed us to gain insights into the structure and energy of the hydrogen bonding network of water molecules around methylsulfonylmethane (MSM). In the hydration sphere of MSM, there are two types of populations of water molecules: a significant share of water molecules weakly bonded to the sulfone group and a smaller share of water molecules strongly bonded to each other around the methyl groups of MSM. The very weak hydrogen bond of water molecules with the hydrophilic group causes the extended network of water hydrogen bonds to be not "anchored" on the sulfone group, and consequently, the MSM hydration shell is labile.

Adenosine Triphosphate: The Primordial Molecule That Controls Protein Homeostasis and Shapes the Genome-Proteome Interface

Adenosine triphosphate (ATP) acts as the universal energy currency that drives various biological processes, while nucleic acids function to store and transmit genetic information for all living organisms. Liquid-liquid phase separation (LLPS) represents the common principle for the formation of membrane-less organelles (MLOs) composed of proteins rich in intrinsically disordered regions (IDRs) and nucleic acids. Currently, while IDRs are well recognized to facilitate LLPS through dynamic and multivalent interactions, the precise mechanisms by which ATP and nucleic acids affect LLPS still remain elusive. This review summarizes recent NMR results on the LLPS of human FUS, TDP-43, and the viral nucleocapsid (N) protein of SARS-CoV-2, as modulated by ATP and nucleic acids, revealing the following: (1) ATP binds to folded domains overlapping with nucleic-acid-binding interfaces; (2) ATP and nucleic acids interplay to biphasically modulate LLPS by competitively binding to overlapping pockets of folded domains and Arg/Lys within IDRs; (3) ATP energy-independently induces protein folding with the highest efficiency known so far. As ATP likely emerged in the prebiotic monomeric world, while LLPS represents a pivotal mechanism to concentrate and compartmentalize rare molecules for forming primordial cells, ATP appears to control protein homeostasis and shape genome-proteome interfaces throughout the evolutionary trajectory, from prebiotic origins to modern cells.

Hydrodynamic Control of Alzheimer Aβ Fibrillation with Glucosaminic Acid Containing Click-Cyclized β-Bodies

It is well established that the dynamic hydration shell plays a vital role in macromolecular functions such as protein-ligand, protein-protein, protein-DNA, and protein-lipid interactions. Here we investigate how the water modality affects conformational changes, solubility, and motion of fibrillar proteins. The hypothesis is that the introduction of a poly hydroxyl amino acid would increase solvation of the fibril forming peptides, preventing their misfolding and aggregation. For the amyloid β (Aβ) peptide, which is considered to be connected with nervous system diseases, including dementia and cognitive decline in Alzheimer's disease, the formation of β-sheet fibrils always occurs with a conformational change and a decrease in the dynamic hydration shell around Aβ(1-42). We present novel cyclic d-amino acid peptides that effectively inhibit fibrillation through affecting the dynamic hydration shell of Aβ(1-42) in vitro. Using de novo design within the software Molecular Operating Environment (MOE), five different peptides that recognize Alzheimer's fibrils were designed and synthesized. Three of them were cyclic all-d-amino acid peptides incorporating the same polyhydroxy building block derived from d-glucosaminic acid (GA). One peptide was the parent cyclic all d-amino acid inhibitor with no GA incorporated, and another was an all l-amino acid linear fibrillation inhibitor. The GA-containing peptides were found to show significantly improved inhibition of Aβ(1-42) aggregation. The inhibition was dramatically improved by the synergistic application of two GA peptides targeting each end of the growing fibril. The present study may facilitate future developments of intervention strategies for Alzheimer's disease and similar neurodegenerative diseases.

The temperature-dependence of host-guest binding thermodynamics: experimental and simulation studies

The thermodynamic parameters of host-guest binding can be used to describe, understand, and predict molecular recognition events in aqueous systems. However, interpreting binding thermodynamics remains challenging, even for these relatively simple molecules, as they are determined by both direct and solvent-mediated host-guest interactions. In this contribution, we focus on the contributions of water to binding by studying binding thermodynamics, both experimentally and computationally, for a series of nearly rigid, electrically neutral host-guest systems and report the temperature-dependent thermodynamic binding contributions ΔGb(T), ΔHb(T), ΔSb(T), and ΔCp,b. Combining isothermal titration calorimetry (ITC) measurements with molecular dynamics (MD) simulations, we provide insight into the binding forces at play for the macrocyclic hosts cucurbit[n]uril (CBn, n = 7-8) and β-cyclodextrin (β-CD) with a range of guest molecules. We find consistently negative changes in heat capacity on binding (ΔCp,b) for all systems studied herein - as well as for literature host-guest systems - indicating increased enthalpic driving forces for binding at higher temperatures. We ascribe these trends to solvation effects, as the solvent properties of water deteriorate as temperature rises. Unlike the entropic and enthalpic contributions to binding, with their differing signs and magnitudes for the classical and non-classical hydrophobic effect, heat capacity changes appear to be a unifying and more general feature of host-guest complex formation in water. This work has implications for understanding protein-ligand interactions and other complex systems in aqueous environments.

Quantifying Amide-Aromatic Interactions at Molecular and Atomic Levels: Experimentally Determined Enthalpic and Entropic Contributions to Interactions of Amide sp2O, N, C and sp3C Unified Atoms with Naphthalene sp2C Atoms in Water

In addition to amide hydrogen bonds and the hydrophobic effect, interactions involving π-bonded sp2 atoms of amides, aromatics, and other groups occur in protein self-assembly processes including folding, oligomerization, and condensate formation. These interactions also occur in aqueous solutions of amide and aromatic compounds, where they can be quantified. Previous analysis of thermodynamic coefficients quantifying net-favorable interactions of amide compounds with other amides and aromatics revealed that interactions of amide sp2O with amide sp2N unified atoms (presumably C═O···H-N hydrogen bonds) and amide/aromatic sp2C (lone pair π, n-π*) are particularly favorable. Sp3C-sp3C (hydrophobic), sp3C-sp2C (hydrophobic, CH-π), sp2C-sp2C (hydrophobic, π-π), and sp3C-sp2N interactions are favorable, sp2C-sp2N interactions are neutral, while sp2O-sp2O and sp2N-sp2N self-interactions and sp2O-sp3C interactions are unfavorable. Here, from determinations of favorable effects of 14 amides on naphthalene solubility at 10, 25, and 45 °C, we dissect amide-aromatic interaction free energies into enthalpic and entropic contributions and find these vary systematically with amide composition. Analysis of these results yields enthalpic and entropic contributions to intrinsic strengths of interactions of amide sp2O, sp2N, sp2C, and sp3C unified atoms with aromatic sp2C atoms. For each interaction, enthalpic and entropic contributions have the same sign and are much larger in magnitude than the interaction free energy itself. The amide sp2O-aromatic sp2C interaction is enthalpy-driven and entropically unfavorable, consistent with direct chemical interaction (e.g., lone pair-π), while amide sp3C- and sp2C-aromatic sp2C interactions are entropy-driven and enthalpically unfavorable, consistent with hydrophobic effects. These findings are relevant for interactions involving π-bonded sp2 atoms in protein processes.

Quantifying Amide-Aromatic Interactions at Molecular and Atomic Levels: Experimentally-determined Enthalpic and Entropic Contributions to Interactions of Amide sp 2 O, N, C and sp 3 C Unified Atoms with Naphthalene sp 2 C Atoms in Water

In addition to amide hydrogen bonds and the hydrophobic effect, interactions involving π-bonded sp 2 atoms of amides, aromatics and other groups occur in protein self-assembly processes including folding, oligomerization and condensate formation. These interactions also occur in aqueous solutions of amide and aromatic compounds, where they can be quantified. Previous analysis of thermodynamic coefficients quantifying net-favorable interactions of amide compounds with other amides and aromatics revealed that interactions of amide sp 2 O with amide sp 2 N unified atoms (presumably C=O···H-N hydrogen bonds) and amide/aromatic sp 2 C (lone pair-π, n-π * ) are particularly favorable. Sp 3 C-sp 3 C (hydrophobic), sp 3 C-sp 2 C (hydrophobic, CH-π), sp 2 C-sp 2 C (hydrophobic, π-π) and sp 3 C-sp 2 N interactions are favorable, sp 2 C-sp 2 N interactions are neutral, while sp 2 O-sp 2 O and sp 2 N-sp 2 N self-interactions and sp 2 O-sp 3 C interactions are unfavorable. Here, from determinations of favorable effects of fourteen amides on naphthalene solubility at 10, 25 and 45 °C, we dissect amide-aromatic interaction free energies into enthalpic and entropic contributions and find these vary systematically with amide composition. Analysis of these results yields enthalpic and entropic contributions to intrinsic strengths of interactions of amide sp 2 O, sp 2 N, sp 2 C and sp 3 C unified atoms with aromatic sp 2 C atoms. For each interaction, enthalpic and entropic contributions have the same sign and are much larger in magnitude than the interaction free energy itself. The amide sp 2 O-aromatic sp 2 C interaction is enthalpy-driven and entropically unfavorable, consistent with direct chemical interaction (e.g. lone pair-π) while amide sp 3 C- and sp 2 C-aromatic sp 2 C interactions are entropy-driven and enthalpically unfavorable, consistent with hydrophobic effects. These findings are relevant for interactions involving π-bonded sp 2 atoms in protein processes.

Table of Contents Graphic

Ordered-domain unfolding of thermophilic isolated β subunit ATP synthase

The flexibility of the ATP synthase's β subunit promotes its role in the ATP synthase rotational mechanism, but its domains stability remains unknown. A reversible thermal unfolding of the isolated β subunit (Tβ) of the ATP synthase from Bacillus thermophilus PS3, tracked through circular dichroism and molecular dynamics, indicated that Tβ shape transits from an ellipsoid to a molten globule through an ordered unfolding of its domains, preserving the β-sheet residual structure at high temperature. We determined that part of the stability origin of Tβ is due to a transversal hydrophobic array that crosses the β-barrel formed at the N-terminal domain and the Rossman fold of the nucleotide-binding domain (NBD), while the helix bundle of the C-terminal domain is the less stable due to the lack of hydrophobic residues, and thus the more flexible to trigger the rotational mechanism of the ATP synthase.

Revisiting macromolecular hydration with HullRadSAS

Hydration of biological macromolecules is important for their stability and function. Historically, attempts have been made to describe the degree of macromolecular hydration using a single parameter over a narrow range of values. Here, we describe a method to calculate two types of hydration: surface shell water and entrained water. A consideration of these two types of hydration helps to explain the "hydration problem" in hydrodynamics. The combination of these two types of hydration allows accurate calculation of hydrodynamic volume and related macromolecular properties such as sedimentation and diffusion coefficients, intrinsic viscosities, and the concentration-dependent non-ideality identified with sedimentation velocity experiments.

New insights on the catalytic center of proteins from peptidylprolyl isomerase group based on the FOD-M model

Generating the structure of the hydrophobic core is based on the orientation of hydrophobic residues towards the central part of the protein molecule with the simultaneous exposure of polar residues. Such a course of the protein folding process takes place with the active participation of the polar water environment. While the self-assembly process leading to the formation of micelles concerns freely moving bi-polar molecules, bipolar amino acids in polypeptide chain have limited mobility due to the covalent bonds. Therefore, proteins form a more or less perfect micelle-like structure. The criterion is the hydrophobicity distribution, which to a greater or lesser extent reproduces the distribution expressed by the 3D Gaussian function on the protein body. The vast majority of proteins must ensure solubility, so a certain part of it-as it is expected-should reproduce the structuring of micelles. The biological activity of proteins is encoded in the part that does not reproduce the micelle-like system. The location and quantitative assessment of the contribution of orderliness to disorder is of critical importance for the determination of biological activity. The form of maladjustment to the 3D Gauss function may be varied-hence the obtained high diversity of specific interactions with strictly defined molecules: ligands or substrates. The correctness of this interpretation was verified on the basis of the group of enzymes Peptidylprolyl isomerase-E.C.5.2.1.8. In proteins representing this class of enzymes, zones responsible for solubility-micelle-like hydrophobicity system-the location and specificity of the incompatible part in which the specific activity of the enzyme is located and coded were identified. The present study showed that the enzymes of the discussed group show two different schemes of the structure of catalytic center (taking into account the status as defined by the fuzzy oil drop model).

Light, Water, and Melatonin: The Synergistic Regulation of Phase Separation in Dementia

The swift rise in acceptance of molecular principles defining phase separation by a broad array of scientific disciplines is shadowed by increasing discoveries linking phase separation to pathological aggregations associated with numerous neurodegenerative disorders, including Alzheimer’s disease, that contribute to dementia. Phase separation is powered by multivalent macromolecular interactions. Importantly, the release of water molecules from protein hydration shells into bulk creates entropic gains that promote phase separation and the subsequent generation of insoluble cytotoxic aggregates that drive healthy brain cells into diseased states. Higher viscosity in interfacial waters and limited hydration in interiors of biomolecular condensates facilitate phase separation. Light, water, and melatonin constitute an ancient synergy that ensures adequate protein hydration to prevent aberrant phase separation. The 670 nm visible red wavelength found in sunlight and employed in photobiomodulation reduces interfacial and mitochondrial matrix viscosity to enhance ATP production via increasing ATP synthase motor efficiency. Melatonin is a potent antioxidant that lowers viscosity to increase ATP by scavenging excess reactive oxygen species and free radicals. Reduced viscosity by light and melatonin elevates the availability of free water molecules that allow melatonin to adopt favorable conformations that enhance intrinsic features, including binding interactions with adenosine that reinforces the adenosine moiety effect of ATP responsible for preventing water removal that causes hydrophobic collapse and aggregation in phase separation. Precise recalibration of interspecies melatonin dosages that account for differences in metabolic rates and bioavailability will ensure the efficacious reinstatement of the once-powerful ancient synergy between light, water, and melatonin in a modern world.

Isotope effects observed in diluted D2O/H2O mixtures identify HOD-induced low-density structures in D2O but not H2O

Normal and heavy water are solvents most commonly used to study the isotope effect. The isotope effect of a solvent significantly influences the behavior of a single molecule in a solution, especially when there are interactions between the solvent and the solute. The influence of the isotope effect becomes more significant in D₂O/H₂O since the hydrogen bond in H₂O is slightly weaker than its counterpart (deuterium bond) in D₂O. Herein, we characterize the isotope effect in a mixture of normal and heavy water on the solvation of a HOD molecule. We show that the HOD molecule affects the proximal solvent molecules, and these disturbances are much more significant in heavy water than in normal water. Moreover, in D₂O, we observe the formation of low-density structures indicative of an ordering of the solvent around the HOD molecule. The qualitative differences between HOD interaction with D₂O and H₂O were consistently confirmed with Raman spectroscopy and NMR diffusometry.

Dry Molten Globule-Like Intermediates in Protein Folding, Function, and Disease

The performance of a protein depends on its correct folding to the final functional native form. Hence, understanding the process of protein folding has remained an important field of research for the scientific community for the past five decades. Two important intermediate states, namely, wet molten globule (WMG) and dry molten globule (DMG), have emerged as critical milestones during protein folding-unfolding reactions. While much has been discussed about WMGs as a common unfolding intermediate, the evidence for DMGs has remained elusive owing to their near-native features, which makes them difficult to probe using global structural probes. This Review puts together the available literature and new evidence on DMGs to give a broader perspective on the universality of DMGs and discuss their significance in protein folding, function, and disease.

A fast-slow method to treat solute dynamics in explicit solvent

Aiming to reduce the computational cost in the current explicit solvent molecular dynamics (MD) simulation, this paper proposes a fast-slow method for the fast MD simulation of biomolecules in explicit solvent. This fast-slow method divides the entire system into two parts: a core layer (typically solute or biomolecule) and a peripheral layer (typically solvent molecules). The core layer is treated using standard MD method but the peripheral layer is treated by a slower dynamics method to reduce the computational cost. We compared four different simulation models in testing calculations for several small proteins. These include gas-phase, implicit solvent, fast-slow explicit solvent and standard explicit solvent MD simulations. Our study shows that gas-phase and implicit solvent models do not provide a realistic solvent environment and fail to correctly produce reliable dynamic structures of proteins. On the other hand, the fast-slow method can essentially reproduce the same solvent effect as the standard explicit solvent model while gaining an order of magnitude in efficiency. This fast-slow method thus provides an efficient approach for accelerating the MD simulation of biomolecules in explicit solvent.

From Protein Design to the Energy Landscape of a Cold Unfolding Protein

Understanding protein folding is crucial for protein sciences. The conformational spaces and energy landscapes of cold (unfolded) protein states, as well as the associated transitions, are hardly explored. Furthermore, it is not known how structure relates to the cooperativity of cold transitions, if cold and heat unfolded states are thermodynamically similar, and if cold states play important roles for protein function. We created the cold unfolding 4-helix bundle DCUB1 with a de novo designed bipartite hydrophilic/hydrophobic core featuring a hydrogen bond network which extends across the bundle in order to study the relative importance of hydrophobic versus hydrophilic protein-water interactions for cold unfolding. Structural and thermodynamic characterization resulted in the discovery of a complex energy landscape for cold transitions, while the heat unfolded state is a random coil. Below ∼0 °C, the core of DCUB1 disintegrates in a largely cooperative manner, while a near-native helical content is retained. The resulting cold core-unfolded state is compact and features extensive internal dynamics. Below -5 °C, two additional cold transitions are seen, that is, (i) the formation of a water-mediated, compact, and highly dynamic dimer, and (ii) the onset of cold helix unfolding decoupled from cold core unfolding. Our results suggest that cold unfolding is initiated by the intrusion of water into the hydrophilic core network and that cooperativity can be tuned by varying the number of core hydrogen bond networks. Protein design has proven to be invaluable to explore the energy landscapes of cold states and to robustly test related theories.

Analysis of water distribution in delipidated human hair by small-angle neutron scattering (SANS)

OBJECTIVE

It is known that damaged hair has a part of its internal structure damaged, and its water absorption and desorption behavior are different. In recent years, it has been reported that internal lipids play an important role in the adsorption and desorption of water to the hair. Therefore, we investigate whether the water distribution in hair and the amount of internal lipids are related.

METHODS

To investigate the effect of internal lipid on water distribution, we prepare human hair samples with and without a partial lack of internal lipids. Internal lipids have been removed using formic acid. The distribution of D₂ O in the hair is investigated using small angle neutron scattering (SANS) under the wet and dry conditions of each hair sample.

RESULTS

It is found from the obtained SANS data that formic acid-treated hairs tended to have fewer 40Å-sized water clusters that were periodically present along the fibre axis in the wet condition. On the other hand, in the dry condition, there were no differences in water distribution between samples.

CONCLUSION

These observations are believed to have been caused by the reduction of 40Å-sized water clusters existing on the constituents removed by formic acid treatment, especially the hydrophobic (lipid) constituent tissues. Consequently, it is clarified that internal lipids are deeply involved in the state of water distribution on hair in wet conditions.

Optical Kerr Effect of Liquid Acetonitrile Probed by Femtosecond Time-Resolved X-ray Liquidography

Optical Kerr effect (OKE) spectroscopy is a method that measures the time-dependent change of the birefringence induced by an optical laser pulse using another optical laser pulse and has been used often to study the ultrafast dynamics of molecular liquids. Here we demonstrate an alternative method, femtosecond time-resolved X-ray liquidography (fs-TRXL), where the microscopic structural motions related to the OKE response can be monitored using a different type of probe, i.e., X-ray solution scattering. By applying fs-TRXL to acetonitrile and a dye solution in acetonitrile, we demonstrate that different types of molecular motions around photoaligned molecules can be resolved selectively, even without any theoretical modeling, based on the anisotropy of two-dimensional scattering patterns and extra structural information contained in the q-space scattering data. Specifically, the dynamics of reorientational (libration and orientational diffusion) and translational (interaction-induced motion) motions are captured separately by anisotropic and isotropic scattering signals, respectively. Furthermore, the two different types of reorientational motions are distinguished from each other by their own characteristic scattering patterns and time scales. The measured time-resolved scattering signals are in excellent agreement with the simulated scattering signals based on a molecular dynamics simulation for plausible molecular configurations, providing the detailed structural description of the OKE response in liquid acetonitrile.

Water-mediated interactions destabilize proteins

Proteins are folded to avoid exposure of the nonpolar groups to water because water‐mediated interactions between nonpolar groups are a promising factor in the thermodynamic stabilities of proteins—which is a well‐accepted view as one of the unique effects of hydrophobic interactions. This article poses a critical question for this classical view by conducting an accurate solvation free‐energy calculation for a thermodynamic cycle of a protein folding using a liquid‐state density functional theory. Here, the solvation‐free energy for a leucine zipper formation was examined in the coiled‐coil protein GCN4‐p1, a typical model for hydrophobic interactions, which demonstrated that water‐mediated interactions were unfavorable for the association of nonpolar groups in the native state, while the dispersion forces between them were, instead, responsible for the association. Furthermore, the present analysis well predicted the isolated helical state stabilized by pressure, which was previously observed in an experiment. We reviewed the problems in the classical concept and semiempirical presumption that the energetic cost of the hydration of nonpolar groups is a driving force of folding.

Protein Surface Interactions-Theoretical and Experimental Studies

In this review, we briefly describe a theoretical discussion of protein folding, presenting the relative contribution of the hydrophobic effect versus the stabilization of proteins via direct surface forces that sometimes may be overlooked. We present NMR-based studies showing the stability of proteins lacking a hydrophobic core which in turn present hydrophobic surface clusters, such as plant defensins. Protein dynamics measurements by NMR are the key feature to understand these dynamic surface clusters. We contextualize the measurement of protein dynamics by nuclear relaxation and the information available at protein surfaces and water cavities. We also discuss the presence of hydrophobic surface clusters in multidomain proteins and their participation in transient interactions which may regulate the function of these proteins. In the end, we discuss how surface interaction regulates the reactivity of certain protein post-translational modifications, such as S-nitrosation.

Adenosine triphosphate energy-independently controls protein homeostasis with unique structure and diverse mechanisms

Proteins function in the crowded cellular environments with high salt concentrations, thus facing tremendous challenges of misfolding/aggregation which represents a pathological hallmark of aging and an increasing spectrum of human diseases. Recently, intrinsically disordered regions (IDRs) were recognized to drive liquid-liquid phase separation (LLPS), a common principle for organizing cellular membraneless organelles (MLOs). ATP, the universal energy currency for all living cells, mysteriously has concentrations of 2-12 mM, much higher than required for its previously-known functions. Only recently, ATP was decoded to behave as a biological hydrotrope to inhibit protein LLPS and aggregation at mM. We further revealed that ATP also acts as a bivalent binder, which not only biphasically modulates LLPS driven by IDRs of human and viral proteins, but also bind to the conserved nucleic-acid-binding surfaces of the folded proteins. Most unexpectedly, ATP appears to act as a hydration mediator to antagonize the crowding-induced destabilization as well as to enhance folding of proteins without significant binding. Here, this review focuses on summarizing the results of these biophysical studies and discussing their implications in an evolutionary context. By linking triphosphate with unique hydration property to adenosine, ATP appears to couple the ability for establishing hydrophobic, π-π, π-cation and electrostatic interactions to the capacity in mediating hydration of proteins, which is at the heart of folding, dynamics, stability, phase separation and aggregation. Consequently, ATP acquired a category of functions at ~mM to energy-independently control protein homeostasis with diverse mechanisms, thus implying a link between cellular ATP concentrations and protein-aggregation diseases.

Dinucleotides as simple models of the base stacking-unstacking component of DNA 'breathing' mechanisms

Regulatory protein access to the DNA duplex ‘interior’ depends on local DNA ‘breathing’ fluctuations, and the most fundamental of these are thermally-driven base stacking-unstacking interactions. The smallest DNA unit that can undergo such transitions is the dinucleotide, whose structural and dynamic properties are dominated by stacking, while the ion condensation, cooperative stacking and inter-base hydrogen-bonding present in duplex DNA are not involved. We use dApdA to study stacking-unstacking at the dinucleotide level because the fluctuations observed are likely to resemble those of larger DNA molecules, but in the absence of constraints introduced by cooperativity are likely to be more pronounced, and thus more accessible to measurement. We study these fluctuations with a combination of Molecular Dynamics simulations on the microsecond timescale and Markov State Model analyses, and validate our results by calculations of circular dichroism (CD) spectra, with results that agree well with the experimental spectra. Our analyses show that the CD spectrum of dApdA is defined by two distinct chiral conformations that correspond, respectively, to a Watson–Crick form and a hybrid form with one base in a Hoogsteen configuration. We find also that ionic structure and water orientation around dApdA play important roles in controlling its breathing fluctuations.

The role of water in the primary nucleation of protein amyloid aggregation

The understanding of the complex conformational landscape of amyloid aggregation and its modulation by relevant physicochemical and cellular factors is a prerequisite for elucidating some of the molecular basis of pathology in amyloid related diseases, and for developing and evaluating effective disease-specific therapeutics to reduce or eliminate the underlying sources of toxicity in these diseases. Interactions of proteins with solvating water have been long considered to be fundamental in mediating their function and folding; however, the relevance of water in the process of protein amyloid aggregation has been largely overlooked. Here, we provide a perspective on the role water plays in triggering primary amyloid nucleation of intrinsically disordered proteins (IDPs) based on recent experimental evidences. The initiation of amyloid aggregation likely results from the synergistic effect between both protein intermolecular interactions and the properties of the water hydration layer of the protein surface. While the self-assembly of both hydrophobic and hydrophilic IDPs would be thermodynamically favoured due to large water entropy contributions, large desolvation energy barriers are expected, particularly for the nucleation of hydrophilic IDPs. Under highly hydrating conditions, primary nucleation is slow, being facilitated by the presence of nucleation-active surfaces (heterogeneous nucleation). Under conditions of poor water activity, such as those found in the interior of protein droplets generated by liquid-liquid phase separation, however, the desolvation energy barrier is significantly reduced, and nucleation can occur very rapidly in the bulk of the solution (homogeneous nucleation), giving rise to structurally distinct amyloid polymorphs. Water, therefore, plays a key role in modulating the transition free energy of amyloid nucleation, thus governing the initiation of the process, and dictating the type of preferred primary nucleation and the type of amyloid polymorph generated, which could vary depending on the particular microenvironment that the protein molecules encounter in the cell.

Protein Denaturation, Zero Entropy Temperature, and the Structure of Water around Hydrophobic and Amphiphilic Solutes

The hydrophobic effect plays a key role in many chemical and biological processes, including protein folding. Nonetheless, a comprehensive picture of the effect of temperature on hydrophobic hydration and protein denaturation remains elusive. Here, we study the effect of temperature on the hydration of model hydrophobic and amphiphilic solutes, through molecular dynamics, aiming at getting insight on the singular behavior of water, concerning the zero-entropy temperature, T_S, and entropy convergence, T_S^*, also observed for some proteins, upon denaturation. We show that, similar to hydrocarbons, polar amphiphilic solutes exhibit a T_S, although strongly dependent on solute-water interactions, opposite to hydrocarbons. Further, the temperature dependence of the hydration entropy, normalized by the solvent accessible surface area, is shown to be nearly solute size independent for hydrophobic but not for amphiphilic solutes, for similar reasons. These results are further discussed in the light of information theory (IT) and the structure of water around hydrophobic groups. The latter shows that the tetrahedral enhancement of some water molecules around hydrophobic groups, associated with the reduction of water defects, leads to the strengthening of the weakest hydrogen bonds, relative to bulk water. In addition, a larger tetrahedrality is found in low density water populations, demonstrating that pure water has encoded structural information, similar to that associated with hydrophobic hydration. The reversal of the hydration entropy dependence on the solute size, above T_S^*, is also analyzed and shown to be associated with a greater loss of water molecules exhibiting enhanced orientational order, in the coordination sphere of large solutes. Finally, the source of the differences between Kauzmann's "hydrocarbon model" on protein denaturation and hydrophobic hydration is discussed, with relatively large amphiphilic hydrocarbons seemingly displaying a more similar behavior to some globular proteins than aliphatic hydrocarbons.

Experimentally determined strengths of favorable and unfavorable interactions of amide atoms involved in protein self-assembly in water

Folding and other protein self-assembly processes are driven by favorable interactions between O, N, and C unified atoms of the polypeptide backbone and side chains. These processes are perturbed by solutes that interact with these atoms differently than water does. Amide NH···O=C hydrogen bonding and various π-system interactions have been better characterized structurally or by simulations than experimentally in water, and unfavorable interactions are relatively uncharacterized. To address this situation, we previously quantified interactions of alkyl ureas with amide and aromatic compounds, relative to interactions with water. Analysis yielded strengths of interaction of each alkylurea with unit areas of different hybridization states of unified O, N, and C atoms of amide and aromatic compounds. Here, by osmometry, we quantify interactions of 10 pairs of amides selected to complete this dataset. An analysis yields intrinsic strengths of six favorable and four unfavorable atom-atom interactions, expressed per unit area of each atom and relative to interactions with water. The most favorable interactions are sp2O-sp2C (lone pair-π, presumably n-π*), sp2C-sp2C (π-π and/or hydrophobic), sp2O-sp2N (hydrogen bonding) and sp3C-sp2C (CH-π and/or hydrophobic). Interactions of sp3C with itself (hydrophobic) and with sp2N are modestly favorable, while sp2N interactions with sp2N and with amide/aromatic sp2C are modestly unfavorable. Amide sp2O-sp2O interactions and sp2O-sp3C interactions are more unfavorable, indicating the preference of amide sp2O to interact with water. These intrinsic interaction strengths are used to predict interactions of amides with proteins and chemical effects of amides (including urea, N-ethylpyrrolidone [NEP], and polyvinylpyrrolidone [PVP]) on protein stability.

Can the roles of polar and non-polar moieties be reversed in non-polar solvents?

Using thermodynamic integration, we study the solvation free energy of 18 amino acid side chain equivalents in solvents with different polarities, ranging from the most polar water to the most non-polar cyclohexane.

The role of water in amyloid aggregation kinetics

The role of water in protein function and aggregation is highly important and may hold some answers to understanding initiation of misfolding diseases such as Parkinson's, Alzheimer's and Huntington's where soluble intrinsically disordered proteins (IDPs) aggregate into fibrillar structures. IDPs are highly dynamic and have larger solvent exposed areas compared to globular proteins, meaning they make and break hydrogen bonds with the surrounding water more frequently. The mobility of water can be altered by presence of ions, sugars, osmolytes, proteins and membranes which differ in different cell types, cell compartments and also as cells age. A reduction in water mobility and thus protein mobility enhances the probability that IDPs can associate to form intermolecular bonds and propagate into aggregates. This poses an interesting question as to whether localised water mobility inside cells can influence the propensity of an IDP to aggregate and furthermore whether it can influence fibril polymorphism and disease outcome.

Femtosecond Hydration Map of Intrinsically Disordered α-Synuclein

Protein hydration water plays a fundamentally important role in protein folding, binding, assembly, and function. Little is known about the hydration water in intrinsically disordered proteins that challenge the conventional sequence-structure-function paradigm. Here, by combining experiments and simulations, we show the existence of dynamical heterogeneity of hydration water in an intrinsically disordered presynaptic protein, namely α-synuclein, implicated in Parkinson's disease. We took advantage of nonoccurrence of cysteine in the sequence and incorporated a number of cysteine residues at the N-terminal segment, the central amyloidogenic nonamyloid-β component (NAC) domain, and the C-terminal end of α-synuclein. We then labeled these cysteine variants using environment-sensitive thiol-active fluorophore and monitored the solvation dynamics using femtosecond time-resolved fluorescence. The site-specific femtosecond time-resolved experiments allowed us to construct the hydration map of α-synuclein. Our results show the presence of three dynamically distinct types of water: bulk, hydration, and confined water. The amyloidogenic NAC domain contains dynamically restrained water molecules that are strikingly different from the water molecules present in the other two domains. Atomistic molecular dynamics simulations revealed longer residence times for water molecules near the NAC domain and supported our experimental observations. Additionally, our simulations allowed us to decipher the molecular origin of the dynamical heterogeneity of water in α-synuclein. These simulations captured the quasi-bound water molecules within the NAC domain originating from a complex interplay between the local chain compaction and the sequence composition. Our findings from this synergistic experimental simulation approach suggest longer trapping of interfacial water molecules near the amyloidogenic hotspot that triggers the pathological conversion into amyloids via chain sequestration, chain desolvation, and entropic liberation of ordered water molecules.

Dihedral Angle Calculations To Elucidate the Folding of Peptides through Its Main Mechanical Forces

This study reports a general method to calculate dihedral angles (φ and ψ) of a given amino acid sequence, focusing on potential energy and torque moment concepts. By defining these physical measures in relation to the chemical interactions that occur on each single amino acid residue within a peptide, we analyze the folding process as the result of main mechanical forces (MMFs) exerted in the specific amino acid chain of interest. As a proof of concept, Leu-enkephalin was initially used as a model peptide to carry out the theoretical study. Our data show agreement between calculated Leu-enkephalin backbone dihedral angles and the corresponding experimentally determined X-ray values. Hence, we used calcitonin to validate our MMF-based method on a larger peptide, i.e., 32 amino acid residues forming an α-helix. Through a similar approach (although simplified with regard to electrostatic interactions), the calculations for calcitonin also demonstrate a good agreement with experimental values. This study offers new opportunities to analyze peptides' amino acid sequences and to help in the prediction of how they must fold, assisting in the development of new computational techniques in the field.

Factors Promoting the Formation of Clathrate-Like Ordering of Water in Biomolecular Structure at Ambient Temperature and Pressure

Clathrate hydrate forms when a hydrophobic molecule is entrapped inside a water cage or cavity. Although biomolecular structures also have hydrophobic patches, clathrate-like water is found in only a limited number of biomolecules. Also, while clathrate hydrates form at low temperature and moderately higher pressure, clathrate-like water is observed in biomolecular structure at ambient temperature and pressure. These indicate presence of other factors along with hydrophobic environment behind the formation of clathrate-like water in biomolecules. In the current study, we presented a systematic approach to explore the factors behind the formation of clathrate-like water in biomolecules by means of molecular dynamics simulation of a model protein, maxi, which is a naturally occurring nanopore and has clathrate-like water inside the pore. Removal of either confinement or hydrophobic environment results in the disappearance of clathrate-like water ordering, indicating a coupled role of these two factors. Apart from these two factors, clathrate-like water ordering also requires anchoring groups that can stabilize the clathrate-like water through hydrogen bonding. Our results uncover crucial factors for the stabilization of clathrate-like ordering in biomolecular structure which can be used for the development of new biomolecular structure promoting clathrate formation.

Unexpected trends in the hydrophobicity of fluorinated amino acids reflect competing changes in polarity and conformation

The hydration free energy of fluorinated amino acids is calculated with molecular simulations and explained with an analytical model.

Hydrophobic Surface Enhances Electrostatic Interaction in Water

A high dielectric constant is one of the peculiar properties of liquid water, indicating that the electrostatic interaction between charged substances is largely reduced in water. We show by molecular dynamics simulation that the dielectric constant of water is decreased near the hydrophobic surface. We further show that the decrease in the dielectric constant is due to both the decreased water density and the reduced water dipole correlation in the direction perpendicular to the surface. We finally demonstrate that electrostatic interaction in water is actually strengthened near the hydrophobic surface.

How do glycerol and dimethyl sulphoxide affect local tetrahedral structure of water around a nonpolar solute at low temperature? Importance of preferential interaction

Glycerol and dimethyl sulphoxide (DMSO) have vital roles in cryoprotection of living cells, tissues, etc. The above action has been directly linked with disruption of hydrogen (H-) bond structure and dynamics of water by these cosolvents at bulk region and around various complex units, such as peptide, amino acid, protein, and lipid membrane. However, the disruption of the local structure of the water solvent around a purely hydrophobic solute is still not studied extensively. The latter is also important in the context of stabilization of protein from cold denaturation. Through all-atom molecular dynamics simulation, we have investigated the comparative effect of glycerol and DMSO on the orientational order of water around a nonpolar solute at -5 °C. A steady reduction of the tetrahedral order of water is observed at bulk (>10 Å distance from the solute) and solute interface (<5.5 Å distance from the solute) with increasing the cosolvent concentration. Contrasting roles of glycerol and DMSO have been evidenced. While DMSO affects the H-bond structure of the interfacial water more than that of the bulk water, glycerol affects the water structure almost uniformly at all regions around the solute. Furthermore, while glycerol helps to retain water molecules at the interface, DMSO significantly reduces the water content in that region. We have put forward a plausible mechanism for these contrasting roles of these cosolvents. The solute-cosolvent hydrophobic-interaction-induced orientational alignment of an interfacial cosolvent molecule determines whether the involvement of the cosolvent molecules in H-bonding with solvent water in the interface is akin to the bulk region or not.

Protein crystal lattices are dynamic assemblies: the role of conformational entropy in the protein condensed phase

Until recently, the occurrence of conformational entropy in protein crystal contacts was considered to be a very unlikely event. A study based on the most accurately refined protein structures demonstrated that side-chain conformational entropy and static disorder might be common in protein crystal lattices. The present investigation uses structures refined using ensemble refinement to show that although paradoxical, conformational entropy is likely to be the major factor in the emergence and integrity of the protein condensed phase. This study reveals that the role of shape entropy and local entropic forces expands beyond the onset of crystallization. For the first time, the complete pattern of intermolecular interactions by protein atoms in crystal lattices is presented, which shows that van der Waals interactions dominate in crystal formation.

Dielectric screening effect of electronic polarization and intramolecular hydrogen bonding

Recent site-resolved hydrogen exchange measurements have uncovered significant discrepancies between simulations and experimental data during protein folding, including the excessive intramolecular hydrogen bonds in simulations. This finding indicates a possibility that intramolecular charge-charge interactions have not included sufficient dielectric screening effect of the electronic polarization. Scaling down peptide atomic charges according to the optical dielectric constant is tested in this study. As a result, the number of intramolecular hydrogen bonds is lower than using unscaled atomic charges while reaching the same levels of helical contents or β-hairpin backbone hydrogen bonds, because van der Waals interactions contribute substantially to peptide folding in water. Reducing intramolecular charge-charge interactions and hydrogen bonding increases conformational search efficiency. In particular, it reduces the equilibrium helical content in simulations using AMBER force field and the energy barrier in folding simulations using CHARMM force field.

Experimental Atom-by-Atom Dissection of Amide-Amide and Amide-Hydrocarbon Interactions in H2O

Quantitative information about amide interactions in water is needed to understand their contributions to protein folding and amide effects on aqueous processes and to compare with computer simulations. Here we quantify interactions of urea, alkylated ureas, and other amides by osmometry and amide-aromatic hydrocarbon interactions by solubility. Analysis of these data yields strengths of interaction of ureas and naphthalene with amide sp2O, amide sp2N, aliphatic sp3C, and amide and aromatic sp2C unified atoms in water. Interactions of amide sp2O with urea and naphthalene are favorable, while amide sp2O-alkylurea interactions are unfavorable, becoming more unfavorable with increasing alkylation. Hence, amide sp2O-amide sp2N interactions (proposed n-σ* hydrogen bond) and amide sp2O-aromatic sp2C (proposed n-π*) interactions are favorable in water, while amide sp2O-sp3C interactions are unfavorable. Interactions of all ureas with sp3C and amide sp2N are favorable and increase in strength with increasing alkylation, indicating favorable sp3C-amide sp2N and sp3C-sp3C interactions. Naphthalene results show that aromatic sp2C-amide sp2N interactions in water are unfavorable while sp2C-sp3C interactions are favorable. These results allow interactions of amide and hydrocarbon moieties and effects of urea and alkylureas on aqueous processes to be predicted or interpreted in terms of structural information. We predict strengths of favorable urea-benzene and N-methylacetamide interactions from experimental information to compare with simulations and indicate how amounts of hydrocarbon and amide surfaces buried in protein folding and other biopolymer processes and transition states can be determined from analysis of urea and diethylurea effects on equilibrium and rate constants.

Water is an active matrix of life for cell and molecular biology

Szent-Győrgi called water the "matrix of life" and claimed that there was no life without it. This statement is true, as far as we know, on our planet, but it is not clear whether it must hold throughout the cosmos. To evaluate that question requires a close consideration of the many varied and subtle roles that water plays in living cells-a consideration that must be free of both an assumed essentialism that gives water an almost mystical life-giving agency and a traditional tendency to see it as a merely passive solvent. Water is a participant in the "life of the cell," and here I describe some of the features of that active agency. Water's value for molecular biology comes from both the structural and dynamic characteristics of its status as a complex, structured liquid as well as its nature as a polar, protic, and amphoteric reagent. Any discussion of water as life's matrix must, however, begin with an acknowledgment that our understanding of it as both a liquid and a solvent is still incomplete.

Crystal structures of highly simplified BPTIs provide insights into hydration-driven increase of unfolding enthalpy

We report a thermodynamic and structural analysis of six extensively simplified bovine pancreatic trypsin inhibitor (BPTI) variants containing 19–24 alanines out of 58 residues. Differential scanning calorimetry indicated a two-state thermal unfolding, typical of a native protein with densely packed interior. Surprisingly, increasing the number of alanines induced enthalpy stabilization, which was however over-compensated by entropy destabilization. X-ray crystallography indicated that the alanine substitutions caused the recruitment of novel water molecules facilitating the formation of protein–water hydrogen bonds and improving the hydration shells around the alanine’s methyl groups, both of which presumably contributed to enthalpy stabilization. There was a strong correlation between the number of water molecules and the thermodynamic parameters. Overall, our results demonstrate that, in contrast to our initial expectation, a protein sequence in which over 40% of the residues are alanines can retain a densely packed structure and undergo thermal denaturation with a large enthalpy change, mainly contributed by hydration.

Effects of high pressure freezing (HPF) on denaturation of natural actomyosin extracted from prawn (Metapenaeus ensis)

Effects of protein denaturation caused by high pressure freezing, involving Pressure-Factors (pressure, time) and Freezing-Factors (temperature, phase transition, recrystallization, ice crystal types), are complicated. In the current study, the conformation and functional changes of natural actomyosin (NAM) under pressure assisted freezing (PAF, ^{100,150,300,400,500MPa}P_-20°C/25min), pressure shift freezing (PSF, ^200MPaP_-20°C/25min), and immersion freezing (^0.1MPaP_-20°C/5min) after pressure was released to 0.1MPa, as compared to normal immersion freezing process (IF, ^0.1MPaP_-20°C/30min). Results indicated that PSF (^200MPaP_-20°C/30min) could reduce the denaturation of frozen NAM and a pressure of 300MPa was the critical point to induce such a denaturation. During the periods of B→D in PSF or B→C→D in PAF, the generation and growth of ice crystals played an important role on changing the secondary and tertiary structure of the treated NAM.

Probing the Acid-Induced Packing Structure Changes of the Molten Globule Domains of a Protein near Equilibrium Unfolding

Using simultaneously scanning small-angle X-ray scattering (SAXS) and UV-vis absorption with integrated online size exclusion chromatography, supplemental with molecular dynamics simulations, we unveil the long-postulated global structure evolution of a model multidomain protein bovine serum albumin (BSA) during acid-induced unfolding. Our results differentiate three global packing structures of the three molten globule domains of BSA, forming three intermediates I₁, I₂, and E along the unfolding pathway. The I₁-I₂ transition, overlooked in all previous studies, involves mainly coordinated reorientations across interconnected molten globule subdomains, and the transition activates a critical pivot domain opening of the protein for entering into the E form, with an unexpectedly large unfolding free energy change of -9.5 kcal mol^-1, extracted based on the observed packing structural changes. The revealed local packing flexibility and rigidity of the molten globule domains in the E form elucidate how collective motions of the molten globule domains profoundly influence the folding-unfolding pathway of a multidomain protein.

Origin of hydrophobicity and enhanced water hydrogen bond strength near purely hydrophobic solutes

Hydrophobicity plays an important role in numerous physicochemical processes from the process of dissolution in water to protein folding, but its origin at the fundamental level is still unclear. The classical view of hydrophobic hydration is that, in the presence of a hydrophobic solute, water forms transient microscopic "icebergs" arising from strengthened water hydrogen bonding, but there is no experimental evidence for enhanced hydrogen bonding and/or icebergs in such solutions. Here, we have used the redshifts and line shapes of the isotopically decoupled IR oxygen-deuterium (O-D) stretching mode of HDO water near small purely hydrophobic solutes (methane, ethane, krypton, and xenon) to study hydrophobicity at the most fundamental level. We present unequivocal and model-free experimental proof for the presence of strengthened water hydrogen bonds near four hydrophobic solutes, matching those in ice and clathrates. The water molecules involved in the enhanced hydrogen bonds display extensive structural ordering resembling that in clathrates. The number of ice-like hydrogen bonds is 10-15 per methane molecule. Ab initio molecular dynamics simulations have confirmed that water molecules in the vicinity of methane form stronger, more numerous, and more tetrahedrally oriented hydrogen bonds than those in bulk water and that their mobility is restricted. We show the absence of intercalating water molecules that cause the electrostatic screening (shielding) of hydrogen bonds in bulk water as the critical element for the enhanced hydrogen bonding around a hydrophobic solute. Our results confirm the classical view of hydrophobic hydration.

How the hydrophobic factor drives protein folding

How hydrophobicity (HY) drives protein folding is studied. The 1971 Nozaki-Tanford method of measuring HY is modified to use gases as solutes, not crystals, and this makes the method easy to use. Alkanes are found to be much more hydrophobic than rare gases, and the two different kinds of HY are termed intrinsic (rare gases) and extrinsic (alkanes). The HY values of rare gases are proportional to solvent-accessible surface area (ASA), whereas the HY values of alkanes depend on special hydration shells. Earlier work showed that hydration shells produce the hydration energetics of alkanes. Evidence is given here that the transfer energetics of alkanes to cyclohexane [Wolfenden R, Lewis CA, Jr, Yuan Y, Carter CW, Jr (2015) Proc Natl Acad Sci USA 112(24):7484-7488] measure the release of these shells. Alkane shells are stabilized importantly by van der Waals interactions between alkane carbon and water oxygen atoms. Thus, rare gases cannot form this type of shell. The very short (approximately picoseconds) lifetime of the van der Waals interaction probably explains why NMR efforts to detect alkane hydration shells have failed. The close similarity between the sizes of the opposing energetics for forming or releasing alkane shells confirms the presence of these shells on alkanes and supports Kauzmann's 1959 mechanism of protein folding. A space-filling model is given for the hydration shells on linear alkanes. The model reproduces the n values of Jorgensen et al. [Jorgensen WL, Gao J, Ravimohan C (1985) J Phys Chem 89:3470-3473] for the number of waters in alkane hydration shells.

Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

BACKGROUND

Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics.

SCOPE OF REVIEW

This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties.

MAJOR CONCLUSIONS

In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration.

GENERAL SIGNIFICANCE

An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.