Enhanced Tetrahedral Ordering of Water Molecules in Minor Grooves of DNA: Relative Role of DNA Rigidity, Nanoconfinement, and Surface Specific Interactions

Hierarchical Hydration Dynamics of RNA with Nano-Water-Pool at Its Core

Many functional RNAs fold into a compact, roughly globular shape by minimizing the electrostatic repulsion between their negatively charged phosphodiester backbone. The fold of such close, compact RNA architecture is often so designed that its outer surface and complex core both are predominately populated by phosphate groups loosely sequestering bases in the intermediate layers. A number of helical junctions maintain the RNA core and its nano-water-pool. While the folding of RNA is manifested by its counterion environment composed of mixed mono- and divalent salts, the concerted role of ion and water in maintaining an RNA fold is yet to be explored. In this work, detailed atomistic simulations of SAM-I and Add Adenine riboswitch aptamers, and subgenomic flavivirus RNA (sfRNA) have been performed in a physiological mixed mono- and divalent salt environment. All three RNA systems have compact folds with a core diameter of range 1-1.7 nm. The spatiotemporal heterogeneity of RNA hydration was probed in a layer-wise manner by distinguishing the core, the intermediate, and the outer layers. The layer-wise decomposition of hydrogen bonds and collective single-particle reorientational dynamics reveal a nonmonotonic relaxation pattern with the slowest relaxation observed at the intermediate layers that involves functionally important tertiary motifs. The slowness of this intermediate layer is attributed to two types of long-resident water molecules: (i) water from ion-hydration layers and (ii) structurally trapped water (distant from ions). The relaxation kinetics of the core and the surface water essentially exposed to the phosphate groups show well-separated time scales from the intermediate layers. In the slow intermediate layers, site-specific ions and water control the functional dynamics of important RNA motifs like kink-turn, observed in different structure-probing experiments. Most interestingly, we find that as the size of the RNA core increases (SAM1 core < sfRNAcore < Add adenine core), its hydration tends to show faster relaxation. The hierarchical hydration and the layer-wise base-phosphate composition uniquely portray the globular RNA to act like a soft vesicle with a quasi-dynamic nano-water-pool at its core.

From structure and dynamics to biomolecular functions: The ubiquitous role of solvent in biology

Biological activity requires a solvent that can provide a suitable environment, which satisfies the twin need for stability and the ability to change. Among all the solvents water plays the most important role. We review, analyze, and comment on recent works on the structure and dynamics of water around biomolecules and their role in specific biological functions. While studies in the past have focused on understanding the biomolecule-water interactions through a hydration layer; recently the attention has shifted towards understanding functions at a molecular level. Such a microscopic understanding clearly requires elucidation of detailed dynamical processes where solvent molecules play an important role. Finally, we comment on the advances made in understanding the role of water inside a biological cell.

Molecular dynamics simulations of structural and dynamical aspects of DNA hydration water

Water is a remarkable liquid, both because of it is intriguing but also because of its importance. Water plays a key role on the structure and function of biological molecules, but on the other hand also the structure and dynamics of water are deeply influenced by its interactions with biological molecules, specially at low temperatures, where water's anomalies are enhanced. Here we present extensive molecular dynamics simulations of water hydrating a oligonucleotide down to very low temperatures (supercooled water), comparing four water models and analyzing the water structure and dynamics in different domains: water in the minor groove, water in the major groove and bulk water. We found that the water in the grooves is slowed down by the interactions with the nucleic acid and a hints of a dynamic transition regarding translational and orientational dynamics were found, specially for the water models TIP4P/2005 and TIP4P-Ew, which also showed the closest agreement with available experimental data. The behavior of water in such extreme conditions is relevant for the study of cryopreservation of biological tissues.

Influence of the extent of hydrophobicity on water organization and dynamics on 2D graphene oxide surfaces

Graphene oxide (GO) nanomaterials are being extensively explored for a wide spectrum of applications, ranging from water desalination to fuel cell applications due to their tunable mechanical, thermal, and electrical...

Microscopic structural features of water in aqueous-reline mixtures of varying compositions

Reline, a mixture of urea and choline chloride in a 2 : 1 molar ratio, is one of the most frequently used deep eutectic solvents. Pure reline and its aqueous solution have large scale industrial use. Owing to the presence of active hydrogen bond formation sites, urea and choline cations can disrupt the hydrogen-bonded network in water. However, a quantitative understanding of the microscopic structural features of water in the presence of reline is still lacking. We carry out extensive all-atom molecular dynamics simulations to elucidate the effect of the gradual addition of co-solvents on the microscopic arrangements of water molecules. We consider four aqueous solutions of reline, between 26.3 and 91.4 wt%. A disruption of the local hydrogen-bonded structure in water is observed upon inclusion of urea and choline chloride. The extent of deviation of the water structure from tetrahedrality is quantified using the tetrahedral order parameter (q_tet). Our analyses show a monotonic increase in the structural disorder as the co-solvents are added. Increase in the q_tet values are observed when highly electro-negative hetero-atoms like nitrogen, oxygen of urea and choline cations are counted as partners of the central water molecules. Further insights are drawn from the characterization of the hydrogen-bonded network in water and we observe the gradual rupturing of water-water hydrogen bonds and their subsequent replacement by the water-urea hydrogen bonds. A negligible contribution from the hydrogen bonds between water and bulky choline cations has also been found. Considering all the constituents as the hydrogen bond partners we calculate the possibility of a successful hydrogen bond formation with a central water molecule. This gives a clear picture of the underlying mechanism of water replacement by urea.

Implications of relaxation dynamics of collapsed conjugated polymeric nanoparticles for light-harvesting applications

The mechanism of the formation of nanoparticles (collapsed state) from the extended state of polymers and their ultrafast excited state relaxation dynamics are illustrated.

Water in biodegradable glucose–water–urea deep eutectic solvent: modifications of structure and dynamics in a crowded environment

Molecular dynamics simulations have been performed on a highly viscous (η ∼ 255 cP) naturally abundant deep eutectic solvent (NADES) composed of glucose, urea and water in a weight ratio of 6 : 4 : 1 at 328 K.

Origin of heat capacity increment in DNA folding: The hydration effect

BACKGROUND

Understanding DNA folding thermodynamics is crucial for prediction of DNA thermal stability. It is now well established that DNA folding is accompanied by a decrease of the heat capacity ∆c_{p, F}, however its molecular origin is not understood. In analogy to protein folding it has been assumed that this is due to dehydration of DNA constituents, however no evidence exists to support this conclusion.

METHODS

Here we analyze partial molar heat capacity of nucleic bases and nucleosides in aqueous solutions obtained from calorimetric experiments and calculate the hydration heat capacity contribution ∆c_p^hyd.

RESULTS

We present hydration heat capacity contributions of DNA constituents and show that they correlate with the solvent accessible surface area. The average contribution for nucleic base dehydration is +0.56 J mol^-1 K^-1 Å^-2 and can be used to estimate the ∆c_{p, F} contribution for DNA folding.

CONCLUSIONS

We show that dehydration is one of the major sources contributing to the observed ∆c_{p, F} increment in DNA folding. Other possible sources contributing to the overall ∆c_{p, F} should be significant but appear to compensate each other to high degree. The calculated ∆c_p^hyd for duplexes and noncanonical DNA structures agree excellently with the overall experimental ∆c_{p, F} values. By contrast, empirical parametrizations developed for proteins result in poor ∆c_p^hyd predictions and should not be applied to DNA folding.

GENERAL SIGNIFICANCE

Heat capacity is one of the main thermodynamic quantities that strongly affects thermal stability of macromolecules. At the molecular level the heat capacity in DNA folding stems from removal of water from nucleobases.

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.

Molecular Insights into the Unusual Structure of an Antifreeze Protein with a Hydrated Core

The primary driving force for protein folding is the formation of a well-packed, anhydrous core. However, recently, the crystal structure of an antifreeze protein, maxi, has been resolved where the core of the protein is filled with water, which apparently contradicts the existing notion of protein folding. Here, we have performed standard molecular dynamics (MD) simulation, replica exchange MD (REMD) simulation, and umbrella sampling using TIP4P water at various temperatures (300, 260, and 240 K) to explore the origin of this unusual structural feature. It is evident from standard MD and REMD simulations that the protein is found to be stable at 240 K in its unusual state. The core of protein has two layers of semi-clathrate water separating the methyl groups of alanine residues from different helical strands. However, with increasing temperature (260 and 300 K), the stability decreases as the core becomes dehydrated, and methyl groups of alanine are tightly packed driven by hydrophobic interactions. Calculation of the potential of mean force by an umbrella sampling technique between a pair of model hydrophobes resembling maxi protein at 240 K shows the stabilization of second solvent-separated minima (SSM), which provides a thermodynamic rationale of the unusual structural feature in terms of weakening of the hydrophobic interaction. Because the stabilization of SSMs is implicated for cold denaturation, it suggests that the maxi protein is so designed by nature where the cold denatured-like state becomes the biologically active form as it works near or below the freezing point of water.

Optimum Number of Anchored Clathrate Water and Its Instantaneous Fluctuations Dictate Ice Plane Recognition Specificities of Insect Antifreeze Protein

Ice recognition by antifreeze proteins (AFPs) is a subject of topical interest. Among several classes of AFPs, insect AFPs are hyperactive presumably due to their ability to adsorb on basal plane. However, the origin of the basal plane binding specificity is not clearly known. Present work aims to provide atomistic insight into the origin of basal plane recognition by an insect antifreeze protein. Free energy calculations reveal that the order of binding affinity of the AFP toward different ice planes is basal plane > prism plane > pyramidal plane. Critical insight reveals that the observed plane specificity is strongly correlated with the number and their instantaneous fluctuations of clathrate water forming hydrogen bonds with both ice binding surface (IBS) of AFP and ice surface, thus anchoring AFP to the ice surface. On basal plane, anchored clathrate water array is highly stable due to exact match in the periodicity of oxygen atom repeat distances of the ice surface and the threonine repeat distances at the IBS. The stability of anchored clathrate water array progressively decreases upon prism and pyramidal plane adsorption due to mismatch between the threonine ladder and oxygen atom repeat distance. Further analysis reveals that hydration around the methyl side-chains of threonine residues becomes highly significant at low temperature which stabilizes the anchored clathrate water array and dual hydrogen-bonding is a consequence of this stability. Structural insight gained from this study paves the way for rational designing of highly potent antifreeze-mimetic with potential industrial applications.

Anomalous water dynamics at surfaces and interfaces: synergistic effects of confinement and surface interactions

In nature, water is often found in contact with surfaces that are extended on the scale of molecule size but small on a macroscopic scale. Examples include lipid bilayers and reverse micelles as well as biomolecules like proteins, DNA and zeolites, to name a few. While the presence of surfaces and interfaces interrupts the continuous hydrogen bond network of liquid water, confinement on a mesoscopic scale introduces new features. Even when extended on a molecular scale, natural and biological surfaces often have features (like charge, hydrophobicity) that vary on the scale of the molecular diameter of water. As a result, many new and exotic features, which are not seen in the bulk, appear in the dynamics of water close to the surface. These different behaviors bear the signature of both water-surface interactions and of confinement. In other words, the altered properties are the result of the synergistic effects of surface-water interactions and confinement. Ultrafast spectroscopy, theoretical modeling and computer simulations together form powerful synergistic approaches towards an understanding of the properties of confined water in such systems as nanocavities, reverse micelles (RMs), water inside and outside biomolecules like proteins and DNA, and also between two hydrophobic walls. We shall review the experimental results and place them in the context of theory and simulations. For water confined within RMs, we discuss the possible interference effects propagating from opposite surfaces. Similar interference is found to give rise to an effective attractive force between two hydrophobic surfaces immersed and kept fixed at a separation of d, with the force showing an exponential dependence on this distance. For protein and DNA hydration, we shall examine a multitude of timescales that arise from frustration effects due to the inherent heterogeneity of these surfaces. We pay particular attention to the role of orientational correlations and modification of the same due to interaction with the surfaces.

[KCl] Dependence of B-DNA Groove Bending Anisotropy

The energetics of B-DNA bending toward the major and minor grooves were quantified by free energy simulations at four different KCl concentrations. Increased [KCl] led to more flexible DNA, with persistence lengths that agreed well with experimental values. At all salt concentrations, major groove bending was preferred, although preferences for major and minor groove bending were similar for the A-tract containing sequence. Since the phosphate repulsions and DNA internal energy favored minor groove bending, the preference for major groove bending was thought to originate from differences in solvation. Water in the minor groove was tighter bound than water in the major groove, and harder to displace than major groove water, which favored the compression of the major groove upon bending. Higher [KCl] decreased the persistence length for both major and minor groove bending but did not greatly affect the free energy spacing between the minor and major groove bending curves. For sequences without A-tracts, salt affected major and minor bending to nearly the same degree, and did not change the preference for major groove bending. For the A-tract containing sequence, an increase in salt concentration decreased the already small energetic difference between major and minor groove bending. Since salts did not significantly affect the relative differences in bending energetics and hydration, it is likely that the increased bending flexibilities upon salt increase are simply due to screening.

Conformational and hydration properties modulate ice recognition by type I antifreeze protein and its mutants

Mutation of wfAFP changes the intrinsic dynamics in such a way that it significantly influences water mediated AFP adsorption on ice.

Anisotropy of B-DNA Groove Bending

DNA bending is critical for DNA packaging, recognition, and repair, and occurs toward either the major or the minor groove. The anisotropy of B-DNA groove bending was quantified for eight DNA sequences by free energy simulations employing a novel reaction coordinate. The simulations show that bending toward the major groove is preferred for non-A-tracts while the A-tract has a high tendency of bending toward the minor groove. Persistence lengths were generally larger for bending toward the minor groove, which is thought to originate from differences in groove hydration. While this difference in stiffness is one of the factors determining the overall preference of bending direction, the dominant contribution is shown to be a free energy offset between major and minor groove bending. The data suggests that, for the A-tract, this offset is largely determined by inherent structural properties, while differences in groove hydration play a large role for non-A-tracts. By quantifying the energetics of DNA groove bending and rationalizing the origins of the anisotropy, the calculations provide important new insights into a key biological process.

Mechanistic insights into temperature-dependent regulation of the simple cyanobacterial hsp17 RNA thermometer at base-pair resolution

The cyanobacterial hsp17 ribonucleicacid thermometer (RNAT) is one of the smallest naturally occurring RNAT. It forms a single hairpin with an internal 1×3-bulge separating the start codon in stem I from the ribosome binding site (RBS) in stem II. We investigated the temperature-dependent regulation of hsp17 by mapping individual base-pair stabilities from solvent exchange nuclear magnetic resonance (NMR) spectroscopy. The wild-type RNAT was found to be stabilized by two critical CG base pairs (C14-G27 and C13-G28). Replacing the internal 1×3 bulge by a stable CG base pair in hsp17^rep significantly increased the global stability and unfolding cooperativity as evidenced by circular dichroism spectroscopy. From the NMR analysis, remote stabilization and non-nearest neighbour effects exist at the base-pair level, in particular for nucleotide G28 (five nucleotides apart from the side of mutation). Individual base-pair stabilities are coupled to the stability of the entire thermometer within both the natural and the stabilized RNATs by enthalpy–entropy compensation presumably mediated by the hydration shell. At the melting point the Gibbs energies of the individual nucleobases are equalized suggesting a consecutive zipper-type unfolding mechanism of the RBS leading to a dimmer-like function of hsp17 and switch-like regulation behaviour of hsp17^rep. The data show how minor changes in the nucleotide sequence not only offset the melting temperature but also alter the mode of temperature sensing. The cyanobacterial thermosensor demonstrates the remarkable adjustment of natural RNATs to execute precise temperature control.

Hydration properties of natural and synthetic DNA sequences with methylated adenine or cytosine bases in the R.DpnI target and BDNF promoter studied by molecular dynamics simulations

Adenine and cytosine methylation are two important epigenetic modifications of DNA sequences at the levels of the genome and transcriptome. To characterize the differential roles of methylating adenine or cytosine with respect to their hydration properties, we performed conventional MD simulations and free energy perturbation calculations for two particular DNA sequences, namely the brain-derived neurotrophic factor (BDNF) promoter and the R.DpnI-bound DNA that are known to undergo methylation of C5-methyl cytosine and N6-methyl adenine, respectively. We found that a single methylated cytosine has a clearly favorable hydration free energy over cytosine since the attached methyl group has a slightly polar character. In contrast, capping the strongly polar N6 of adenine with a methyl group gives a slightly unfavorable contribution to its free energy of solvation. Performing the same demethylation in the context of a DNA double-strand gave quite similar results for the more solvent-accessible cytosine but much more unfavorable results for the rather buried adenine. Interestingly, the same demethylation reactions are far more unfavorable when performed in the context of the opposite (BDNF or R.DpnI target) sequence. This suggests a natural preference for methylation in a specific sequence context. In addition, free energy calculations for demethylating adenine or cytosine in the context of B-DNA vs. Z-DNA suggest that the conformational B-Z transition of DNA transition is rather a property of cytosine methylated sequences but is not preferable for the adenine-methylated sequences investigated here.

Small-molecule binding to the DNA minor groove is mediated by a conserved water cluster

High-resolution crystal structures of the DNA duplex sequence d(CGCGAATTCGCG)(2) complexed with three minor-groove ligands are reported. A highly conserved cluster of 11 linked water molecules has been found in the native and all 3 ligand-bound structures, positioned at the boundary of the A/T and G/C regions where the minor groove widens. This cluster appears to play a key structural role in stabilizing noncovalently binding small molecules in the AT region of the B-DNA minor groove. The cluster extends from the backbone phosphate groups along the mouth of the groove and links to DNA and ligands by a network of hydrogen bonds that help to maintain the ligands in position. This arrangement of water molecules is distinct from, but linked by, hydrogen bonding to the well-established spine of hydration, which is displaced by bound ligands. Features of the water cluster and observed differences in binding modes help to explain the measured binding affinities and thermodynamic characteristics of these ligands on binding to AT sites in DNA.

Molecular dynamics study of the role of the spine of hydration in DNA A-tracts in determining nucleosome occupancy

A-tracts in DNA are generally associated with reduced nucleosome occupancy relative to other sequences, such that the longer the A-tract, the less likely that nucleosomes are found. In this paper, we use molecular dynamics methods to study the structural properties of A-tracts, and in particular the role that the spine of hydration in A-tracts plays in allowing DNA to distort to the highly bent structure needed to form nucleosomes. This study includes a careful assessment of the ability of the Amber (parmbsc0), CHARMM27, and BMS force fields to describe these structural waters for the AAATTT sequence (here capped with CGC and GCG), including comparisons with X-ray results. All three force fields show a spine of hydration, but BMS and Amber show better correlation with measured properties, such as in narrowing of the minor groove width associated with the A-tract. We have used Amber to study the spine properties for several 6 and 14 base-pair A-tracts (all capped with CGC and GCG). These calculations show that the structural waters are tightly bound for "pure" A-tracts that allow for A-water-T links, and for AT steps that allow for a T-water-T link, but other sequences disfavor structural water, especially those that lead to A-water-A, G-water-G, and C-water-A structures. In addition, we show that pure A-tracts favor roll values close to the Watson-Crick value for linear DNA, while A-tract sequences containing embedded T's, C's, or G's that are less favorable to structural water are more flexible. This implies the essential role of the spine of hydration in disfavoring nucleosome formation.

Hydrophobicity of proteins and interfaces: insights from density fluctuations

Macroscopic characterizations of hydrophobicity (e.g., contact angle measurements) do not extend to the surfaces of proteins and nanoparticles. Molecular measures of hydrophobicity of such surfaces need to account for the behavior of hydration water. Theory and state-of-the-art simulations suggest that water density fluctuations provide such a measure; fluctuations are enhanced near hydrophobic surfaces and quenched with increasing surface hydrophilicity. Fluctuations affect conformational equilibria and dynamics of molecules at interfaces. Enhanced fluctuations are reflected in enhanced cavity formation, more favorable binding of hydrophobic solutes, increased compressibility of hydration water, and enhanced water-water correlations at hydrophobic surfaces. These density fluctuation-based measures can be used to develop practical methods to map the hydrophobicity/philicity of heterogeneous surfaces including those of proteins. They highlight that the hydrophobicity of a group is context dependent and is significantly affected by its environment (e.g., chemistry and topography) and especially by confinement. The ability to include information about hydration water in mapping hydrophobicity is expected to significantly impact our understanding of protein-protein interactions as well as improve drug design and discovery methods and bioseparation processes.

Dynamic properties of water around a protein-DNA complex from molecular dynamics simulations

Formation of protein-DNA complex is an important step in regulation of genes in living organisms. One important issue in this problem is the role played by water in mediating the protein-DNA interactions. In this work, we have carried out atomistic molecular dynamics simulations to explore the heterogeneous dynamics of water molecules present in different regions around a complex formed between the DNA binding domain of human TRF1 protein and a telomeric DNA. It is demonstrated that such heterogeneous water motions around the complex are correlated with the relaxation time scales of hydrogen bonds formed by those water molecules with the protein and DNA. The calculations reveal the existence of a fraction of extraordinarily restricted water molecules forming a highly rigid thin layer in between the binding motifs of the protein and DNA. It is further proved that higher rigidity of water layers around the complex originates from more frequent reformations of broken water-water hydrogen bonds. Importantly, it is found that the formation of the complex affects the transverse and longitudinal degrees of freedom of surrounding water molecules in a nonuniform manner.

Effects of an unnatural base pair replacement on the structure and dynamics of DNA and neighboring water and ions

Incorporating small molecule probes into biomolecular systems to report on local structure and dynamics is a powerful strategy that underlies a wide variety of experimental techniques, including fluorescence, electron paramagnetic resonance (EPR), and Forster resonance energy transfer (FRET) measurements. When an unnatural probe is inserted into a protein or DNA, the degree to which the presence of the probe has perturbed the local structure and dynamics it was intended to study is always an important concern. Here, molecular dynamics (MD) simulations are used to systematically study the effect of replacing a DNA base pair with a fluorescent probe, coumarin 102 deoxyriboside, at six unique sites along an A-tract DNA dodecamer. While the overall structure of the DNA oligonucleotide remains intact, replacement of A*T base pairs leads to widespread structural and dynamic perturbations up to four base pairs away from the probe site, including widening of the minor groove and increased DNA flexibility. New DNA conformations, not observed in the native sequence, are sometimes found in the vicinity of the probe and its partner abasic site analog. Strong correlations are demonstrated between DNA surface topology and water mobility.