Correlated Dynamical Crossovers of the Hydration Layer of a Single-Stranded DNA Oligomer

Solvent constraints for biopolymer folding and evolution in extraterrestrial environments

We propose that spontaneous folding and molecular evolution of biopolymers are two universal aspects that must concur for life to happen. These aspects are fundamentally related to the chemical composition of biopolymers and crucially depend on the solvent in which they are embedded. We show that molecular information theory and energy landscape theory allow us to explore the limits that solvents impose on biopolymer existence. We consider 54 solvents, including water, alcohols, hydrocarbons, halogenated solvents, aromatic solvents, and low molecular weight substances made up of elements abundant in the universe, which may potentially take part in alternative biochemistries. We find that along with water, there are many solvents for which the liquid regime is compatible with biopolymer folding and evolution. We present a ranking of the solvents in terms of biopolymer compatibility. Many of these solvents have been found in molecular clouds or may be expected to occur in extrasolar planets.

Exploring Heterogeneous Dynamical Environment around an Ensemble of Aβ42 Peptide Monomer Conformations

Exploring the conformational properties of amyloid β (Aβ) peptides and the role of solvent (water) in guiding the dynamical environment at their interfaces is crucial for microscopic understanding of Aβ misfolding, which is involved in causing the most common neurodegenerative disorder, i.e., Alzheimer's disease. While numerous studies in the past have emphasized examining the conformational states of Aβ peptides, the role of water has not received much attention. Here, we have performed all-atom molecular dynamics simulations of several full-length Aβ₄₂ peptide monomers with different initial configurations. Our efforts are directed toward probing the origin of the heterogeneous dynamics of water around various segments of the Aβ peptide, identified as the two terminal segments (N-term and C-term) and the two hydrophobic segments (hp1 and hp2), along with the central turn region interconnecting hp1 and hp2. Our results revealed that water hydrating hp1, hp2, and turn (nonterminal segments) and C-term segments exhibit nonuniformly restricted translational as well as rotational motions. The degree of such restriction has been found to be correlated with the hydrogen bond relaxation time scales at the interface. Importantly, it is revealed that the water molecules around hp1 and, to some extent, around hp2, form relatively rigid hydration layers, compared to that around the other segments. Such rigid hydration layers arise due to relatively more solid-like caging motions resulting in relatively lesser hydration entropy. As hp1 and hp2 have been demonstrated to play a central role in Aβ aggregation, we believe that distinct water dynamics in the vicinity of these two segments, as outlined in this study, can provide vital information in understanding the early stages of the onset of the aggregation process of such peptides at higher concentration that can further aid toward advances in AD therapeutics.

Dynamic Heterogeneity at the Interface of an Intrinsically Disordered Peptide

It is believed that water around an intrinsically disordered protein or peptide (IDP) in an aqueous environment plays an important role in guiding its conformational properties and aggregation behavior. However, despite its importance, only a handful of studies exploring the correlation between the conformational motions of an IDP and the microscopic properties of water at its surface are reported. Attempts have been made in this work to study the dynamic properties of water present in the vicinity of α-synuclein, an IDP associated with Parkinson's disease (PD). Room temperature molecular dynamics (MD) simulations of eight α-synuclein_1-95 peptides with a wide range of initial conformations have been carried out in aqueous media. The calculations revealed that due to solid-like caging motions, the translational and rotational mobility of water molecules near the surfaces of the peptide repeat unit segments R1 to R7 are significantly restricted. A small degree of dynamic heterogeneity in the hydration environment around the repeat units has been observed with water near the hydrophobic R6 unit exhibiting relatively more restricted diffusivity. The time scales involving the overall structural relaxations of peptide-water and water-water hydrogen bonds near the peptide have been found to be correlated with the time scale of diffusion of the interfacial water molecules. We believe that the relatively more hindered dynamic environment near R6 can help create water-mediated contacts centered around R6 between peptide monomers at a higher concentration, thereby enhancing the early stages of peptide aggregation.

Dynamical crossover of water confined within the amphiphilic nanocores of aggregated amyloid β peptides

It is believed that the self-assembly of amyloid beta (Aβ) peptides in the brain is the cause of Alzheimer's disease. Atomistic molecular dynamics simulations of aqueous solutions of Aβ protofilaments of different sizes at room temperature have been carried out to explore the dynamic properties of water confined within the core and at the exterior surface of the protofilaments. Attempts have been made to understand how the non-uniform distortion of the protofilaments associated with their structural crossover influences the diffusivity and the hydrogen bonding environment of the confined water molecules. In contrast to the homogeneous solvent dynamical environment at the exterior surface, the calculations revealed heterogeneously restricted motions of water confined within the distorted cores of the protofilaments. Importantly, it is demonstrated that the structural crossover of the aggregates observed for the decamer is associated with a dynamical transition of water confined within its core. A direct one-to-one correlation between the heterogeneously restricted core water motions and the kinetics of the breaking and formation of hydrogen bonds quantitatively demonstrated that a modified hydrogen bond arrangement within the cores of higher order Aβ protofilaments is the origin behind the crossover in core water mobility. A fraction of the water molecules forming short-lived water-water hydrogen bonds within the core of the crossover protofilament decamer are believed to diffuse away easily from the core and thus play a crucial role in further growth of the protofilament by facilitating the binding of new peptide monomers.

THz frequency spectrum of protein-solvent interaction energy using a recurrence plot-based Wiener-Khinchin method

The dynamics of a protein and the water surrounding it are coupled via nonbonded energy interactions. This coupling can exhibit a complex, nonlinear, and nonstationary nature. The THz frequency spectrum for this interaction energy characterizes both the vibration spectrum of the water hydrogen bond network, and the frequency range of large amplitude modes of proteins. We use a Recurrence Plot based Wiener-Khinchin method RPWK to calculate this spectrum, and the results are compared to those determined using the classical auto-covariance-based Wiener-Khinchin method WK. The frequency spectra for the total nonbonded interaction energy extracted from molecular dynamics simulations between the β-Lactamase Inhibitory Protein BLIP, and water molecules within a 10 Å distance from the protein surface, are calculated at 150, 200, 250, and 310 K, respectively. Similar calculations are also performed for the nonbonded interaction energy between the residues 49ASP, 53TYR, and 142PHE in BLIP, with water molecules within 10 Å from each residue respectively at 150, 200, 250, and 310 K. A comparison of the results shows that RPWK performs better than WK, and is able to detect some frequency data points that WK fails to detect. This points to the importance of using methods capable of taking the complex nature of the protein-solvent energy landscape into consideration, and not to rely on standard linear methods. In general, RPWK can be a valuable addition to the analysis tools for protein molecular dynamics simulations. Proteins 2016; 84:1549-1557. © 2016 Wiley Periodicals, Inc.

Determination of pKa values for deprotonable nucleobases in short model oligonucleotides

The deprotonation of ionizable nucleobases centrally placed in short model oligonucleotides was examined under different physical conditions, using UV absorption spectroscopy. The oligonucleotide sequences were designed so that only the central base would be ionized over the pH range examined. pKa values of 9.90±0.01 and 9.34±0.04 were determined for the guanine group in the oligomer d-ACAGCAC and 2'-deoxyguanosine, respectively, both at 25°C and 0.1M NaCl. Lengthening the oligonucleotide up to the tridecamer stage further increases the pKa of the central guanine moiety. Electrolyte concentration, temperature, and mixed water-ethanol solvents affect the acidity of the central base. Changes in the sequence surrounding the central guanine can also have a significant effect, especially in the case of strongly stacking sequences. The pKa values were also determined for the hepta(2'-O-methyl)ribonucleotide and the heptamer PNA of identical sequence, as well as for oligodeoxyribonucleotides with different deprotonable bases, viz. thymine, uracil, or hypoxanthine, in the central position. The results are interpreted in terms of the electric-field effect exerted on the departing proton by the negative electric charges located on the internucleotide phosphate groups, and calculations show this effect to approximately explain the magnitude of the pKa difference observed between the deoxyriboheptanucleotide and its electroneutral PNA analogue.

Effect of temperature on the low-frequency vibrational spectrum and relative structuring of hydration water around a single-stranded DNA

Molecular dynamics simulations of the single-stranded DNA oligomer (5'-CGCGAAT TCGCG-3') in aqueous solution have been carried out at different temperatures between 160 K and 300 K. The effects of temperature on the low-frequency vibrational spectrum and local structural arrangements of water molecules hydrating the DNA strand have been explored in detail. The low-frequency density of states distributions reveal that increasingly trapped transverse water motions play a dominant role in controlling the band corresponding to O⋯O⋯O bending or transverse oscillations of hydration water at supercooled temperatures. In addition, presence of a broad band around 260 (±20) cm(-1) under supercooled conditions indicates transformation from high density liquid-like structuring of hydration water at higher temperatures to that of a low density liquid at lower temperatures. It is found that long-range correlations between the supercooled hydration water molecules arise due to such local structural transition around the DNA oligomer.

Collective excitations in soft-sphere fluids

Despite that the thermodynamic distinction between a liquid and the corresponding gas ceases to exist at the critical point, it has been recently shown that reminiscence of gaslike and liquidlike behavior can be identified in the supercritical fluid region, encoded in the behavior of hypersonic waves dispersion. By using a combination of molecular dynamics simulations and calculations within the approach of generalized collective modes, we provide an accurate determination of the dispersion of longitudinal and transverse collective excitations in soft-sphere fluids. Specifically, we address the decreasing rigidity upon density reduction along an isothermal line, showing that the positive sound dispersion, an excess of sound velocity over the hydrodynamic limit typical for dense liquids, displays a nonmonotonic density dependence strictly correlated to that of thermal diffusivity and kinematic viscosity. This allows rationalizing recent observation parting the supercritical state based on the Widom line, i.e., the extension of the coexistence line. Remarkably, we show here that the extremals of transport properties such as thermal diffusivity and kinematic viscosity provide a robust definition for the boundary between liquidlike and gaslike regions, even in those systems without a liquid-gas binodal line. Finally, we discuss these findings in comparison with recent results for Lennard-Jones model fluid and with the notion of the "rigid-nonrigid" fluid separation lines.