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

Elucidating the Sluggish Water Dynamics at the Ice-Binding Surface of the Hyperactive Tenebrio molitor Antifreeze Protein

Quasi-ice-like hydration waters on the ice-binding surface (IBS) of an antifreeze protein (AFP) commonly exhibit sluggish dynamics especially at low temperatures. In this work, we have analyzed molecular dynamics (MD) simulation trajectories at two different temperatures for Tenebrio molitor antifreeze protein (TmAFP) to explore whether the unique quasi-ice-like structuring of hydration water has any impact on making their dynamics slower on the IBS of the protein. Our calculation reveals that, as translational dynamics is coupled with the conformational fluctuations, hydration water on the IBS exhibits sluggish translational motion due to reduced flexibility of the IBS compared to that on the non-ice-binding surface (NIBS) of the protein. Interestingly, it is noticed that rotational motion of hydration water is not coupled with the conformational fluctuations of the surfaces. In that case, structural relaxations of the protein-water (PW) and water-water (WW) hydrogen bonds compete with each other to make the rotational dynamics of hydration water around the IBS either faster or slower with respect to those around the NIBS. At low temperature, the slower structural relaxation of water-water hydrogen bonds dominates and imparts sluggish rotational motion of the hydration water on the IBS of the protein. The slower structural relaxation of water-water hydrogen bonds and hence the retarded rotational dynamics, despite the weak short-lived PW hydrogen bonds on the IBS, is clearly a manifestation of the rigid quasi-ice-like structure of the hydration shell on that surface.

Dynamical Manifestations of Supercooling in Amyloid Hydration

The effect of extreme temperature on amyloidogenic species remains sparsely explored. In a recent study (J. Phys. Chem. Lett., 2019, 10, (10)), we employed exhaustive molecular dynamics simulations to explore the cold thermal response of a putative small amyloid oligomer and to elicit the role of solvent modulation. Herein, we investigate the dynamical response of the hydration waters of the oligomer within the supercooled states. Using NMR-based formalism, we delineate the entropic response in terms of the side-chain conformational entropy that corroborates the weakening of the hydrophobic core with lowering of temperature. The translational dynamics of the protein and hydration waters reveal the coupling of protein dynamical fluctuations with solvent dynamics under supercooled conditions. Probing the translational motion as a space-time correlation indicates glassy dynamics exhibited by hydration waters in the supercooled regime. Caging of the water molecules with lowering of temperature and the resultant hopping dynamics are reflected in the longer β-relaxation timescales of translational motion. Furthermore, we utilized mode-coupling theory (MCT) and derived the ideal glass transition temperature from translational and rotational dynamics, around ∼196 and 209 K, respectively. Interestingly, rotational motion in the supercooled regime deviates from the MCT law, exhibits Arrhenius motion, and marks a fragile-to-strong crossover at 227 K. The low-frequency vibrational modes also coincide with the dynamical transition. This exposition lends dynamical insights into the hydration coupling of an amyloid aggregate under cryogenic conditions.

Hydration of a small protein under carbon nanotube confinement: Adsorbed substates induce selective separation of the dynamical response

The co-involvement of biological molecules and nanomaterials has increasingly come to the fore in modern-day applications. While the "bio-nano" (BN) interface presents physico-chemical characteristics that are manifestly different from those observed in isotropic bulk conditions, the underlying molecular reasons remain little understood; this is especially true of anomalies in interfacial hydration. In this paper, we leverage atomistic simulations to study differential adsorption characteristics of a small protein on the inner (concave) surface of a single-walled carbon nanotube whose diameter exceeds dimensions conducive to single-file water movement. Our findings indicate that the extent of adsorption is decided by the degree of foldedness of the protein conformational substate. Importantly, we find that partially folded substates, but not the natively folded one, induce reorganization of the protein hydration layer into an inner layer water closer to the nanotube axis and an outer layer water in the interstitial space near the nanotube walls. Further analyses reveal sharp dynamical differences between water molecules in the two layers as observed in the onset of increased heterogeneity in rotational relaxation and the enhanced deviation from Fickian behavior. The vibrational density of states reveals that the dynamical distinctions are correlated with differences in crucial bands in the power spectra. The current results set the stage for further systematic studies of various BN interfaces vis-à-vis control of hydration properties.

Dynamic behaviors of interfacial water on the self-assembly monolayer (SAM) heterogeneous surface

Dynamic behaviors of water molecules near the surface with mixed hydrophobic and hydrophilic areas are studied by molecular dynamics simulation. More specifically, the diffusion coefficient and hydrogen bond lifetime of interfacial water on the self-assembly monolayer composed of hydrophobic and hydrophilic groups and their dependence on the mixing ratio are studied. The diffusion dramatically slows down, and the hydrogen bond lifetime considerably increases when a few hydrophilic groups are added to the hydrophobic surface. When the percentage of hydrophilic groups increases to 25%, the behavior of interfacial water is similar to the case of the pure hydrophilic surface. The sensitivity to the hydrophilic group can be attributed to the fact that the grafted hydrophilic groups can not only retard the directly bound water molecules but also affect indirectly bound water by stabilizing hydrogen bonds among interfacial water molecules.

Heterogeneous Dynamical Environment at the Interface of a Protein-DNA Complex

Binding between protein and DNA is an essential process to regulate different biological activities. Two puzzling questions in protein-DNA recognition are (i) how the protein's binding domain identifies the DNA sequence in an aqueous solution and (ii) how the formation of the complex alters the dynamical environment around it. In this work, we present results obtained from molecular dynamics simulations of the N-terminal α-helical domain of the λ-repressor protein (in dimeric form) bound to the corresponding operator DNA. Effects of formation of the complex in modifying the microscopic dynamics of water as well as the kinetics of hydrogen bonds at the interface have been explored. Locally heterogeneous restricted water motions at the complex interface have been observed, the extent of restriction being more significant around the directly bound residues of the protein and the DNA. In particular, the calculation revealed the existence of significantly constrained motionally restricted water layer that can form either bridges around the directly bound residues of the protein and DNA or are engaged in forming water-mediated contacts between a fraction of the unbound residues. More importantly, it is observed that the restricted water motion around the complex is correlated with the hydrogen bond relaxation time scale at the interface. It is further demonstrated that the kinetics of water-water hydrogen bonds involving the bridged water are influenced more due to complex formation.

Unique Structures and Vibrational Spectra of Protic Ionic Liquids Confined in TiO2 Slits: The Role of Interfacial Hydrogen Bonds

The ionic liquid (IL)/titanium dioxide (TiO₂) interface exists in many application systems, such as nanomaterial synthesis, catalysis, and electrochemistry systems. The nanoscale interfacial properties in the above systems are a common issue. However, directly detecting the interfacial properties of nanoconfined ILs by experimental methods is still challenging. To help better learn about the interfacial issue, molecular dynamics simulations have been performed to explore the structures, vibration spectra, and hydrogen bond (HB) properties at the IL/TiO₂ interface. Ethylammonium nitrate (EAN) ILs confined in TiO₂ slit pores with different pore widths were studied. A unique vibrational spectrum appeared for EAN ILs confined in a 0.7 nm TiO₂ slit, and this phenomenon is related to interfacial hydrogen bonds (HBs). An analysis of the HB types indicated that the interfacial NH₃⁺ group of the cations was in an asymmetric HB environment in the 0.7 nm TiO₂ slit, which led to the disappearance of the symmetric N-H stretching mode. In addition, the significant increase in the HB strength between NH₃⁺ groups and the TiO₂ surface slowed down the stretching vibration of the N-H bond, resulting in one peak in the vibrational spectra at a lower frequency. For the first time, our simulation work establishes a molecular-level relationship between the vibrational spectrum and the local HB environment of nanoconfined ILs at the IL/TiO₂ interface, and this relationship is helpful for interface design in related systems.

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.

Solvation of Carbon Nanoparticles in Water/Alcohol Mixtures: Using Molecular Simulation To Probe Energetics, Structure, and Dynamics

Molecular dynamics simulations were used to examine the solvation behavior of buckminsterfullerene and single-walled carbon nanotubes (SWCNT) in a range of water/alcohol solvent compositions at 1 atm and 300 K. Results indicate that the alcohols assume the role of pseudosurfactants by shielding the nanotube from the more unfavorable interactions with polar water molecules. This is evident in both the free energies of transfer (ΔΔG_water→xOH = -68.1 kJ/mol and -86.5 kJ/mol for C₆₀ in methanol and ethanol; ΔΔG_water→xOH = -345.6 kJ/mol and -421.2 kJ/mol for the (6,5)-SWCNT in methanol and ethanol) and the composition of the solvation shell at intermediate alcohol concentrations. Additionally, we have observed the retardation of both the translational and rotational dynamics of molecules near the nanoparticle surface through use of time correlation functions. A 3-fold increase in the residence times of the alcohol molecules within the solvation shells at low concentrations further reveals their surfactant-like behavior. Such interactions are important when considering the complex molecular environment present in many schemes used for nanoparticle purification techniques.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

Characterizing Structural Stability of Amyloid Motif Fibrils Mediated by Water Molecules

In biological systems, structural confinements of amyloid fibrils can be mediated by the role of water molecules. However, the underlying effect of the dynamic behavior of water molecules on structural stabilities of amyloid fibrils is still unclear. By performing molecular dynamics simulations, we investigate the dynamic features and the effect of interior water molecules on conformations and mechanical characteristics of various amyloid fibrils. We find that a specific mechanism induced by the dynamic properties of interior water molecules can affect diffusion of water molecules inside amyloid fibrils, inducing their different structural stabilities. The conformation of amyloid fibrils induced by interior water molecules show the fibrils' different mechanical features. We elucidate the role of confined and movable interior water molecules in structural stabilities of various amyloid fibrils. Our results offer insights not only in further understanding of mechanical features of amyloids as mediated by water molecules, but also in the fine-tuning of the functional abilities of amyloid fibrils for applications.

Microscopic understanding of the conformational features of a protein–DNA complex

Protein–DNA interactions play crucial roles in different stages of genetic activities, such as replication of genome, initiation of transcription,etc.

Influence of Polarization on Carbohydrate Hydration: A Comparative Study Using Additive and Polarizable Force Fields

Carbohydrates are known to closely modulate their surrounding solvent structures and influence solvation dynamics. Spectroscopic investigations studying far-IR regions (below 1000 cm(-1)) have observed spectral shifts in the libration band (around 600 cm(-1)) of water in the presence of monosaccharides and polysaccharides. In this paper, we use molecular dynamics simulations to gain atomistic insight into carbohydrate-water interactions and to specifically highlight the differences between additive (nonpolarizable) and polarizable simulations. A total of six monosaccharide systems, α and β anomers of glucose, galactose, and mannose, were studied using additive and polarizable Chemistry at HARvard Macromolecular Mechanics (CHARMM) carbohydrate force fields. Solvents were modeled using three additive water models TIP3P, TIP4P, and TIP5P in additive simulations and polarizable water model SWM4 in polarizable simulations. The presence of carbohydrate has a significant effect on the microscopic water structure, with the effects being pronounced for proximal water molecules. Notably, disruption of the tetrahedral arrangement of proximal water molecules was observed due to the formation of strong carbohydrate-water hydrogen bonds in both additive and polarizable simulations. However, the inclusion of polarization resulted in significant water-bridge occupancies, improved ordered water structures (tetrahedral order parameter), and longer carbohydrate-water H-bond correlations as compared to those for additive simulations. Additionally, polarizable simulations also allowed the calculation of power spectra from the dipole-dipole autocorrelation function, which corresponds to the IR spectra. From the power spectra, we could identify spectral signatures differentiating the proximal and bulk water structures, which could not be captured from additive simulations.

Effect of temperature on the structure and hydration layer of TATA-box DNA: A molecular dynamics simulation study

DNA within the living cells experiences a diverse range of temperature, ranging from freezing condition to hot spring water. How the structure, the mechanical properties of DNA, and the solvation dynamics around DNA changes with the temperature is important to understand the functionality of DNA under those acute temperature conditions. In that notion, we have carried out molecular dynamics simulations of a DNA oligomer, containing TATA-box sequence for three different temperatures (250K, 300K and 350K). We observed that the structure of the DNA, in terms of backbone torsion angles, sugar pucker, base pair parameters, and base pair step parameters, did not show any unusual properties within the studied range of temperatures, but significant structural alteration was noticed between BI and BII forms at higher temperature. As expected, the flexibility of the DNA, in terms of the torsional rigidity and the bending rigidity is highly temperature dependent, confirming that flexibility increases with increase in temperature. Additionally, the groove widths of the studied DNA showed temperature sensitivity, specifically, the major groove width decreases and the minor groove width increases, respectively, with the increase in temperature. We observed that at higher temperature, water around both the major and the minor groove of the DNA is less structured. However, the water dynamics around the minor groove of the DNA is more restricted as compared to the water around the major groove throughout the studied range of temperatures, without any anomalous behavior.

Effects of protein-DNA complex formation on the intermolecular vibrational density of states of interfacial water

Single-stranded DNAs (ss-DNAs) are formed as intermediates during DNA metabolic processes. ss-DNA binding (SSB) proteins specifically bind to the single-stranded segments of the DNA and protect it from being degraded. We have performed room temperature molecular dynamics simulations of the aqueous solution of two DNA binding K homology (KH) domains (KH3 and KH4) of the far upstream element (FUSE) binding protein (FBP) complexed with two ss-DNA oligomers. Efforts have been made to explore the influence of complex formation on low-frequency vibrational density of states of the surface water molecules. It is revealed that increased back scattering of water confined around the complexed structures leads to significant blue shifts of the band corresponding to the O···O···O bending or restricted transverse motions of water, the effect being more for the bridged water molecules. Importantly, it is demonstrated that the formation of such complexed structures of a similar type may often influence the transverse and longitudinal degrees of freedom of the surrounding water molecules in a nonuniform manner.

A comprehensive model for the recognition of human telomeres by TRF1

Eukaryotic chromosomes are capped by telomeres, nucleoprotein complexes that prevent chromosome end-to-end fusions and control cell ageing. Two proteins in this complex, telomere repeat binding factors (TRF1 and TRF2), specifically recognise the double-stranded TTAGGG tandem repeat sequence. TRF1 is a homodimer with roles governing DNA architecture and negatively regulating telomere length. We explore the conformational space of this protein-DNA complex using molecular dynamics and, for the first time, generate a complete model of TRF1-DNA recognition that has not been possible on the basis of crystallographic and NMR data alone. The results reconcile previous conflicting experimental models for the sequence selectivity of the recognition process, by confirming many of the findings while identifying important new interactions and behaviour. This improved characterisation also reveals extensive indirect readout, which suggests that recognition will be affected by changes to DNA helical parameters such as bending.

Competition between Na⁺ and Rb⁺ in the minor groove of DNA

The competition between Na⁺ and Rb⁺ ions in the minor groove of a synthetic B-DNA dodecamer d (CGCGAATTCGCG) is studied using molecular dynamics simulations as the ratio of these two ions changing from 9:1 to 1:9 with the DNA merged into the solvent of water molecule at 298 K. When the concentration of Rb⁺ ions increases, from minority to majority, Na⁺ ions are gradually released from the A tract, and the binding sites in the minor groove are occupied by Rb⁺ ions, extending from the A tract to the whole minor groove. Comparing Na⁺ with Rb⁺ ions, the conformation of the minor groove is influenced strongly by Na⁺ ions.