Hydrogen bond dynamics at vapour–water and metal–water interfaces

Discrepant Effects of Hydrated and Neat Reline on the Conformational Stability of a Knotted Protein

Although knotted proteins are rare in number, their peculiar topology has long intrigued the scientific community. In this study, we have explored the conformational stability of a trefoil-knotted protein, YbeA, in reline (choline chloride:urea in a 1:2 ratio), a well-characterized deep eutectic solvent, using classical molecular dynamics simulation. Deep eutectic solvents (DESs) are explored as a reliable alternative to conventional solvents, effectively altering a protein's structural stability and activity, either stabilizing its native state or disordering its conformation depending on the relevant interactions involved. Here, using pure and hydrated concentrations of reline, we observe the conflicting effect of the DES on the knotted protein's stability. Our studies at room temperature and elevated temperatures show that in pure reline, the protein is conformationally stable and rigid. In contrast, the protein tends to lose its structural integrity in hydrated reline. The stable knotted topology also gets untied as the protein, solvated in hydrated reline, is exposed to an elevated temperature. Using Minimum Distance Distribution Functions and Kirkwood-Buff Integrals, we analyzed the solvation pattern of the DES constituents around the protein. We expect that this study will lead to more effective strategies in developing tailored solvent systems for comprehending the conformational behavior of knotted proteins.

Probing the Degree of Restriction in Solvent Dynamics at the Interface of a Protein-RNA Complex

Protein-RNA complexation is an important step for the regulation of numerous biological processes. Water present at the interface of a protein-RNA complex plays a critical role in guiding its structure, stability, and function. Therefore, studying the microscopic properties of interfacial water is essential to gain molecular insights into the formation of such complexes. In this study, we present results obtained from molecular dynamics (MD) simulations of poly(A)-binding protein (PABP) bound with poly(A) RNA, which is an essential regulatory step to control the deadenylation process, thereby stabilizing cellular mRNAs from degradation. Efforts have been made to explore how such complexation alters the regular dynamical and hydrogen bond properties of water present at the interface. The calculations revealed restricted water dynamics at the interface, characterized by heterogeneous time scales, with the extent of restriction being more pronounced for residues directly involved in protein-RNA binding. In particular, water molecules around the protein's linker, RRM2, and the RNA strand exhibit significantly more restricted motion compared to RRM1 upon complexation. Further, longer relaxation times of hydrogen bonds at the interface due to complex formation have been found to be correlated with increasingly restricted water motions. Notably, the kinetics of hydrogen bonds around the protein's linker, RRM2, and the RNA strand are more strongly influenced by complex formation, underscoring their critical role in mediating protein-RNA interactions.

Modulating the Self-Assembly of a Camptothecin Prodrug with Paclitaxel for Anticancer Combination Therapy: A Molecular Dynamics Approach

Camptothecin (CPT) and paclitaxel (PTX), derived from natural products, are recognized for their significant efficacy in clinical cancer treatments. Despite its therapeutic advantages, CPT is challenged by issues of toxicity and solubility, necessitating its use in conjugation with other compounds for enhanced compatibility. This study delves into the coassembly mechanism of Evans blue-conjugated camptothecin (EB-CPT) with PTX, aiming to elucidate their synergistic potential in combination therapy applications, employing all-atom molecular dynamics simulations. The EB-CPT prodrug is reported to form a self-aggregated cluster. Our findings suggest that increasing the PTX concentration induces a dispersion of EB-CPT clusters, thereby disrupting their inherent self-assembly. This disruption is explained to be facilitated by the coassembly of EB-CPT and PTX. With increasing concentration of PTX, a lengthening of the coassembled structures is observed, supporting the experimental findings of tube-like coassembled structures at higher weight ratios of PTX. Hydrophobic interactions and π-π stacking are the primary forces responsible for the formation of both self- and coassembled structures. Interestingly, the structural analysis reveals that the CPT moiety of EB-CPT is less involved in assemblies due to steric hindrances. Instead, the interaction and coassembly processes are predominantly mediated by the EB derivative component of the prodrug. This research underscores the critical role of the solubilizing agent, EB derivative, in mediating the flexibility and interaction of CPT in combination therapy strategies, particularly with PTX, thus emphasizing the importance of conjugates for therapeutic developments. Furthermore, the molecular insights into the interaction sites and mechanisms facilitating coassembly between EB-CPT and PTX contribute valuable knowledge to the field, highlighting the potential of these nanomedicine combinations in advancing cancer treatment modalities.

Comparative Bindings of Glycosaminoglycans with CXCL8 Monomer and Dimer: Insights from Conformational Dynamics and Kinetics of Hydrogen Bonds

GAGs bind to both the monomeric and dimeric forms of CXCL8, helping to form a concentration gradient of the chemokine that facilitates the recruitment of neutrophils to an injury site and supports other biological functions. In this study, atomistic molecular dynamics simulations were conducted to investigate the binding behavior of two hexameric GAGs sulfated at two different positions, chondroitin sulfate (CS) and heparan sulfate (HS), with the monomer (SIL8) and dimer (DIL8) forms of the CXCL8 protein. The results support that the conformational diversity of CS and HS appeared to be more when binding with monomer SIL8 than dimer DIL8. CS gained more configurational entropy from glycosidic linkage flexibility when bound to SIL8 than DIL8, with a higher energy barrier, whereas HS exhibited a lower energy barrier for configurational entropy when bound to SIL8 and DIL8. The monomer SIL8 exhibited more favorable and preferential binding with GAGs compared to DIL8. Formation of hydrogen bonds with the basic amino acids of SIL8 and GAG was more rigid and required higher activation energy to break than the other identified hydrogen bondings. Water molecules involved in hydrogen bonding with GAGs, excluding those with basic amino acids of DIL8, showed longer lifetimes and slower relaxation compared to SIL8. This suggests that water-mediated interactions also favor binding of DIL8 with GAGs. Despite having more basic amino acids, DIL8 did not display stronger binding than SIL8, indicating the significant role of basic residues in stabilizing the GAG-protein interactions in the monomers. The reason could be that the greater number of nonbasic amino acids in dimeric CXCL8 stabilizes the complex by forming water-mediated hydrogen bonds, reducing the conformational preferences for binding with GAGs. In contrast, the monomeric form of CXCL8 exhibits a higher conformational preference for protein-GAG interactions.

Revisiting the Thickness of the Air-Water Interface from Two Extremes of Interface Hydrogen Bond Dynamics

The air-water interface plays a crucial role in many aspects of science because of its unique properties, such as a two-dimensional hydrogen bond (HB) network and completely different HB dynamics compared to bulk water. However, accurately determining the boundary of interfacial and bulk water, that is, the thickness of the air-water interface, still challenges experimentalists. Various simulation-based methods have been developed to estimate the thickness, converging on a range of approximately 3-10 (Å). In this study, we introduce a novel approach, grounded in density functional theory-based molecular dynamics and deep potential molecular dynamics simulations, to measure the air-water interface thickness, offering a different perspective based on prior research. To capture realistic HB dynamics in the air-water interface, two extreme scenarios of the interface HB dynamics are obtained: one underestimates the interface HB dynamics, while the other overestimates it. Surprisingly, our results suggest that the interface HB dynamics in both scenarios converges as the thickness of the air-water interface increases to 4 (Å). This convergence point, indicative of the realistic interface thickness, is also validated by our calculation of anisotropic decay of OH stretch and the free OH dynamics at the air-water interface.

Salt-bridge mediated conformational dynamics in the figure-of-eight knotted ketol acid reductoisomerase (KARI)

The utility of knotted proteins in biological activities has been ambiguous since their discovery. From their evolutionary significance to their functionality in stabilizing the native protein structure, a unilateral conclusion hasn't been achieved yet. While most studies have been performed to understand the stabilizing effect of the knotted fold on the protein chain, more ideas are yet to emerge regarding the interactions in stabilizing the knot. Using classical molecular dynamics (MD) simulations, we have explored the dynamics of the figure-of-eight knotted domain present in ketol acid reductoisomerase (KARI). Our main focus was on the presence of a salt bridge network evident within the knotted region and its role in shaping the conformational dynamics of the knotted chain. Through the potential of mean forces (PMFs) calculation, we have also marked the specific salt bridges that are pivotal in stabilizing the knotted structure. The correlated motions have been further monitored with the help of principal component analysis (PCA) and dynamic cross-correlation maps (DCCM). Furthermore, mutation of the specific salt bridges led to a change in their conformational stability, vindicating their importance.

Dynamics of confined water inside carbon nanotubes based on studying tetrahedral order parameters

Water dynamics inside hydrophobic confinement, such as carbon nanotubes (CNTs), has garnered significant attention, focusing on water diffusion. However, a crucial aspect remains unexplored - the influence of confinement size on water ordering and intrinsic hydrogen bond dynamics. To address this gap, we conducted extensive molecular dynamics simulations to investigate local ordering and intrinsic hydrogen bond dynamics of water molecules within CNTs of various sizes (length:20 nm, diameters: 1.0 nm to 5.0 nm) over a wide range of temperatures (260K, 280K, 300K, and 320K). A striking observation emerged: in smaller CNTs, water molecules adopt an icy structure near tube walls while maintaining liquid state towards the center. Notably, water behavior within a 2.0 nm CNT stands out as an anomaly, distinct from other CNT sizes considered in this study. This anomaly was explained through the formation of water layers inside CNTs. The hydrogen bond correlation function of water within CNTs decayed more slowly than bulk water, with an increasing rate as CNT diameter increased. In smaller CNTs, water molecules hold onto their hydrogen bond longer than larger ones. Interestingly, in larger CNTs, the innermost layer's hydrogen bond lasts a shorter time compared to the other layers, and this changes with temperature.

Second inflection point of supercooled water surface tension induced by hydrogen bonds: A molecular-dynamics study

Surface tension of supercooled water is a fundamental property in various scientific processes. In this study, we perform molecular dynamics simulations with the TIP4P-2005 model to investigate the surface tension of supercooled water down to 220 K. Our results show a second inflection point (SIP) in the surface tension at temperature TSIP ≈ 267.5 ± 2.3 K. Using an extended IAPWS-E functional fit for the water surface tension, we calculate the surface excess internal-energy and entropy terms of the excess Helmholtz free energy. Similar to prior studies [Wang et al., Phys. Chem. Chem. Phys. 21, 3360 (2019); Gorfer et al., J. Chem. Phys. 158, 054503 (2023)], our results show that the surface tension is governed by two driving forces: a surface excess entropy change above the SIP and a surface excess internal-energy change below it. We study hydrogen-bonding near the SIP because it is the main cause of water's anomalous properties. With decreasing temperature, our results show that the entropy contribution to the surface tension reaches a maximum slightly below the SIP and then decreases. This is because the number of hydrogen bonds increases more slowly below the SIP. Moreover, the strengths and lifetimes of the hydrogen bonds also rise dramatically below the SIP, causing the internal-energy term to dominate the excess surface free energy. Thus, the SIP in the surface tension of supercooled TIP4P-2005 water is associated with an increase in the strengths and lifetimes of hydrogen bonds, along with a decrease in the formation rate (#/K) of new hydrogen bonds.

Understanding Conformational Properties and Role of Hydrogen Bonds in Glycosaminoglycans-Interleukin8 Complexes in Aqueous Medium by Molecular Dynamics Simulation

Atomistic molecular dynamics simulations were performed under ambient conditions to explore the conformational features and binding affinities of hexameric glycosaminoglycans (GAGs) with chemokine Interleukin8 (IL8) in an aqueous medium. We tried to understand the role of hydrogen bonds (HBs) involving conserved water in mediating the interactions. The Luzar-Chandler model was adopted to study the kinetics of HB breaking and formation concerning different water-mediated HBs. The conformational flexibilities of bound GAGs are due to the flexible glycosidic linkages than the occasional/rare ring pucker conformation. The free energy landscape constructed with ϕ, and ψ, depicted that different conformational minima associated with the glycosidic linkage flexibility of the GAGs in bound states are separated by energy barriers. The binding affinities of IL8 towards GAGs are favored through the electrostatic and non-polar solvation interactions. 4-different types of conserved water were explored in the solvent-mediated binding of GAGs with IL8. The average lifetime of the IL8-GAG direct HB pairs was ∼ten times less than the IL8-GAG-shared water HBs. This is due to the rapid establishment of HB breaking and reformation kinetics involving water of a shared layer. We find that despite the highly negatively charged surface of GAGs, the IL8 surface populated by non-cationic amino acids could serve as a promising binding site in addition to the cationic surface of the protein.

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.

Water molecules in CNT-Si3N4 membrane: Properties and the separation effect for water-alcohol solution

Water confined in carbon nanotubes (CNTs) has been intensively studied because of its unique properties and potential for various applications and is often embedded in silicon nitride (Si₃N₄) membranes. However, the understanding of the influence of Si₃N₄ on the properties of water in CNTs lacks clarity. In this study, we performed molecular dynamics simulations to investigate the effect of the Si₃N₄ membrane on water molecules inside CNTs. The internal electric field generated in the CNTs by the point charges of the Si₃N₄ membrane changes the structure and dynamical properties of water in the nanotubes, causing it to attain a disordered structure. The Si₃N₄ membrane decreases the diffusivity of water in the CNTs; this is because the Coulomb potential energy (i.e., electrostatic interaction) of water decreases owing to the presence of Si₃N₄, whereas the Lennard-Jones potential energy (i.e., van der Waals interaction) does not change significantly. Furthermore, electrostatic interactions make the water structure more stable in the CNTs. As a result, the Si₃N₄ membrane enhances the separation effect of the water-methanol mixture with CNTs in the presence of an external electric field. Furthermore, the threshold of the external electric field strength to induce water-methanol separation with CNTs is reduced owing to the presence of a silicon nitride membrane.

Controlling the self-assembly of human calcitonin: a theoretical approach using molecular dynamics simulations

Human calcitonin (hCT) is a 32-residue amino acid poly-peptide hormone which is secreted by the C-cells (also known as parafollicular cells) of thyroid glands. It acts to inhibit osteoclast cell hormones by reducing the cell function and regulating calcium and phosphate in blood. hCT has a high tendency to assemble into protofilaments with β-sheet conformations. Amyloid fibril formation of hCT reduces its bio-activity and limits its application as a therapeutic drug. Salmon calcitonin (sCT), which also carries the same disulfide bridge at the N and C-terminus, but differs at the 16 residue position from hCT, has less propensity to aggregate than hCT. Human calcitonin has much higher bio-activity than sCT if its aggregation propensity is reduced. Substituting the key residues which are responsible for the aggregation of hCT, is one of the ways to reduce its aggregation and fibril formation. hCT analogues with less aggregation tendency can be exploited as therapeutic drugs. In this work, we study the amyloidogenic behavior of hCT and its peptide based derivatives i.e., sCT, phCT, N17H hCT, Y12L hCT and DM hCT, through classical molecular dynamics (MD) simulations. Our study reveals that sCT is the least aggregation prone derivative, and the double mutation at position 12 and 17 can reduce the aggregation propensity of this peptide. Also, we have applied these mutant variants of hCT as peptide inhibitors in the self-aggregation of hCT. This study could help in understanding and preparing peptide-based inhibitors for hCT fibrillation and their applications as therapeutic drugs.

Molecular Insight into the Effects of Enhanced Hydrophobicity on Amyloid-like Aggregation

Generally, hydrophobic amino acids provide hydrophobic interactions during peptide aggregation. However, besides the hydrophobic amino acids, some hydrophilic amino acids, such as glutamine, are also considered to be essential elements in many self-aggregating peptides. For example, huntingtin contains polyglutamine at its N-terminus and the yeast prion Sup35 protein has the GNNQQNY sequence, a peptide well-known for its ability for amyloid fibril formation. However, despite the frequent emergence of glutamine in self-assembling systems, the molecular mechanism of amyloid formation involving this unique amino acid has not been well documented. It is still not clear how this hydrophilic amino acid is responsible for the hydrophobic interaction in the self-association process. Therefore, in this study, we have carried out classical molecular dynamics simulations of the GNNQQNY peptide and its derivatives in pure water. We quantify the propensity for the formation of β-sheet conformation with an increasing glutamine number in the peptide sequence. In addition, we assess the importance of the hydrophobicity of the dimethanediyl group present in glutamine (as well as in glutamic acid) for the self-association of the peptides through nonpolar solvent medium simulations.

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.

Molecular dynamics simulation study on the inhibitory effects of choline-O-sulfate on hIAPP protofibrilation

Type 2 diabetes mellitus (T2Dm) is a neurodegenerative disease, which occurs due to the self-association of human islet amyloid polypeptide (hIAPP), also known as human amylin. It was reported experimentally that choline-O-sulfate (COS), a small organic molecule having a tertiary amino group and sulfate group, can prevent the aggregation of human amylin without providing the mechanism of the action of COS in the inhibition process. In this work, we investigate the influence of COS on the full-length hIAPP peptide by performing 500 ns classical molecular dynamics simulations. From pure water simulation (without COS), we have identified the residues 11-20 and 23-36 that mainly participate in the fibril formation, but in the presence of 1.07 M COS these residues become totally free of β-sheet conformation. Our results also show that the sulfate oxygen of COS directly interacts with the peptide backbone, which leads to the local disruption of peptide-peptide interaction. Moreover, the presence of favorable peptide-COS vdW interaction energy and high coordination number of COS molecules in the first solvation shell of the peptide indicates the hydrophobic solvation of the peptide residues by COS molecules, which also play a crucial role in the prevention of β-sheet formation. Finally, from the potential of mean force (PMFs) calculations, we observe that the free energy between two peptides is more negative in the absence of COS and with increasing concentration of COS, it becomes unfavorable significantly indicating that the peptide dimer formation is most stable in pure water, which becomes less favorable in the presence of COS. © 2019 Wiley Periodicals, Inc.

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.

Structural and Dynamic Heterogeneity of Capillary Wave Fronts at Aqueous Interfaces

Using a unique combination of slab-layering analyses and identification of truly interfacial molecules, this work examines water/vapor and water/n-hexane interfaces, specifically the structural and dynamic perturbations of the interfacial water molecules at different locations within the surface capillary waves. From both the structural and dynamic properties analyzed, it is found that these interfacial water molecules dominate the perturbations within the interfacial region, which can extend deep into the water phase relative to the Gibbs dividing surface. Of more importance is the demonstration of structural and dynamic heterogeneity of the interfacial water molecules at the capillary wave front, as indicated by the dipole orientation and the structural and dynamic behavior of hydrogen bonds and their networks.

Microscopic dynamics of charge separation at the aqueous electrochemical interface

We have used molecular simulation and methods of importance sampling to study the thermodynamics and kinetics of ionic charge separation at a liquid water-metal interface. We have considered this process using canonical examples of two different classes of ions: a simple alkali-halide pair, Na⁺I^-, or classical ions, and the products of water autoionization, H₃O⁺OH^-, or water ions. We find that for both ion classes, the microscopic mechanism of charge separation, including water's collective role in the process, is conserved between the bulk liquid and the electrode interface. However, the thermodynamic and kinetic details of the process differ between these two environments in a way that depends on ion type. In the case of the classical ion pairs, a higher free-energy barrier to charge separation and a smaller flux over that barrier at the interface result in a rate of dissociation that is 40 times slower relative to the bulk. For water ions, a slightly higher free-energy barrier is offset by a higher flux over the barrier from longer lived hydrogen-bonding patterns at the interface, resulting in a rate of association that is similar both at and away from the interface. We find that these differences in rates and stabilities of charge separation are due to the altered ability of water to solvate and reorganize in the vicinity of the metal interface.

Slow hydrogen-bond switching dynamics at the water surface revealed by theoretical two-dimensional sum-frequency spectroscopy

Using our newly developed explicit three-body (E3B) water model, we simulate the surface of liquid water. We find that the timescale for hydrogen-bond switching dynamics at the surface is about three times slower than that in the bulk. In contrast, with this model rotational dynamics are slightly faster at the surface than in the bulk. We consider vibrational two-dimensional (2D) sum-frequency generation (2DSFG) spectroscopy as a technique for observing hydrogen-bond rearrangement dynamics at the water surface. We calculate the nonlinear susceptibility for this spectroscopy for two different polarization conditions, and in each case we see the appearance of cross-peaks on the timescale of a few picoseconds, signaling hydrogen-bond rearrangement on this timescale. We thus conclude that this 2D spectroscopy will be an excellent experimental technique for observing slow hydrogen-bond switching dynamics at the water surface.

Importance of protein conformational motions and electrostatic anchoring sites on the dynamics and hydrogen bond properties of hydration water

The microscopic dynamic properties of water molecules present in the vicinity of a protein are expected to be sensitive to its local conformational motions and the presence of polar and charged groups at the surface capable of anchoring water molecules through hydrogen bonds. In this work, we attempt to understand such sensitivity by performing detailed molecular dynamics simulations of the globular protein barstar solvated in aqueous medium. Our calculations demonstrate that enhanced confinement at the protein surface on freezing its local motions leads to increasingly restricted water mobility with long residence times around the secondary structures. It is found that the inability of the surface water molecules to bind with the protein residues by hydrogen bonds in the absence of protein-water (PW) electrostatic interactions is compensated by enhanced water-water hydrogen bonds around the protein with uniform bulklike behaviors. Importantly, it is further noticed that in contrast to the PW hydrogen bond relaxation time scale, the kinetics of the breaking and formation of such bonds are not affected on freezing the protein's conformational motions.

Restricted dynamics of water around a protein-carbohydrate complex: computer simulation studies

Water-mediated protein-carbohydrate interaction is a complex phenomenon responsible for different biological processes in cellular environment. One of the unexplored but important issues in this area is the role played by water during the recognition process and also in controlling the microscopic properties of the complex. In this study, we have carried out atomistic molecular dynamics simulations of a protein-carbohydrate complex formed between the hyaluronan binding domain of the murine Cd44 protein and the octasaccharide hyaluronan in explicit water. Efforts have been made to explore the heterogeneous influence of the complex on the dynamic properties of water present in different regions around it. It is revealed from our analyses that the heterogeneous dynamics of water around the complex are coupled with differential time scales of formation and breaking of hydrogen bonds at the interface. Presence of a highly rigid thin layer of motionally restricted water molecules bridging the protein and the carbohydrate in the common region of the complex has been identified. Such water molecules are expected to play a crucial role in controlling properties of the complex. Importantly, it is demonstrated that the formation of the protein-carbohydrate complex affects the transverse and longitudinal degrees of freedom of the interfacial water molecules in a heterogeneous manner.

Influence of urea on tert-butyl alcohol aggregation in aqueous solutions

Ternary solutions consisting of urea, tert-butyl alcohol (TBA), and water are investigated employing molecular dynamics simulations. The main purpose of the present paper is to investigate the effect of urea on TBA aggregation and by extension its influence on hydrophobic interactions. The aggregation of TBA can be detected from the concentration dependence of structural properties such as first-shell TBA-water coordination numbers and TBA-TBA hydrogen-bond numbers, as well as through changes in the translational diffusion coefficients of TBA. It is found that urea acts to delay the association of TBA to concentrations greater than those required to cause TBA aggregation in binary TBA-water systems. It is shown that urea acts through a direct mechanism, whereby it preferentially binds to TBA replacing water from the first coordination shell. TBA-urea hydrogen bonds can be as strong as, or stronger than, those of TBA-water, and urea binds to both the hydrophilic and hydrophobic moieties of TBA. Our observations are qualitatively consistent with experimental results for urea-TBA-water solutions and with recent simulation studies of urea's action as a protein denaturant.

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.