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.     

Local heterogeneous dynamics of water around lysozyme: a computer simulation study

Water present near the surface of a protein exhibits dynamic properties different from that of water in the pure bulk state. In this work, we have carried out atomistic molecular dynamics simulation of an aqueous solution of hen egg-white lysozyme. Attempts have been made to explore the correlation between the local heterogeneous mobility of water around the protein segments and the rigidity of the hydration layers with the microscopic dynamics of hydrogen bonds formed by water molecules with the protein residues. The kinetics of breaking and reformation of hydrogen bonds involving the surface water molecules have been calculated. It is found that the reformations of broken hydrogen bonds are more frequent for the hydration layers of those segments of the protein that are more rigid. The calculation of the low-frequency vibrational modes of hydration layer water molecules reveals that the protein influences the transverse and longitudinal degrees of freedom of water around it in a differential manner. These findings can be verified by appropriate experimental studies.     

Conformational fluctuations of a protein-DNA complex and the structure and ordering of water around it

Protein-DNA binding is an important process responsible for the regulation of genetic activities in living organisms. The most crucial issue in this problem is how the protein recognizes the DNA and identifies its target base sequences. Water molecules present around the protein and DNA are also expected to play an important role in mediating the recognition process and controlling the structure of the complex. We have performed atomistic molecular dynamics simulations of an aqueous solution of the protein-DNA complex formed between the DNA binding domain of human TRF1 protein and a telomeric DNA. The conformational fluctuations of the protein and DNA and the microscopic structure and ordering of water around them in the complex have been explored. In agreement with experimental studies, the calculations reveal conformational immobilization of the terminal segments of the protein on complexation. Importantly, it is discovered that both structural adaptations of the protein and DNA, and the subsequent correlation between them to bind, contribute to the net entropy loss associated with the complex formation. Further, it is found that water molecules around the DNA are more structured with significantly higher density and ordering than that around the protein in the complex.     

Conformational flexibility of a protein-carbohydrate complex and the structure and ordering of surrounding water

Protein-carbohydrate non-covalent interactions are important to understand various biological processes in living organisms. One of the important issues in protein-carbohydrate binding is how the protein identifies the target carbohydrate and recognizes its conformational features. Surrounding water molecules are expected to play a critical role not only in mediating the recognition process but also in maintaining the structure of the complex. We carried out atomistic molecular dynamics (MD) simulations of an aqueous solution of the protein-carbohydrate complex formed between the hyaluronan binding domain (HABD) of the murine Cd44 protein and the octasaccharide hyaluronan (HA(8)). The conformational flexibilities of the protein and the carbohydrate, and the microscopic structure and ordering of water molecules around them in the complexed form have been explored. It is revealed that the formation of the complex is associated with significant immobilization of the monosaccharide units of the carbohydrate moiety that are involved in binding. Further, reduction in water densities around the binding residues of the two molecules in the complex with respect to their free forms clearly demonstrated that the recognition between the protein and the carbohydrate is facilitated by removal of a fraction of water molecules from regions around the binding domains.     

Role of hydrogen bonds in protein-DNA recognition: effect of nonplanar amino groups

Amino groups are one of the various types of hydrogen bond donors, abundantly found in protein main chains, protein side chains, and DNA bases. The polar hydrogen atoms of these groups exhibit short ranged, specific, and directional hydrogen bonds, which play a decisive role in the specificity and stability of protein-DNA complexes. To date, planar amino groups are only considered for the analysis of protein-DNA interfacial hydrogen bonds. This assumption regarding hydrogen atom positions possibly failed to establish the expected role of hydrogen bonds in protein-DNA recognition. We have performed ab initio quantum chemical studies on amino acid side chains and DNA bases containing amino groups as well as on specific hydrogen bonded residue pairs selected from high-resolution protein-DNA complex crystal structures. Our results suggest that occurrences of pyramidal amino groups are more probable in comparison with the usually adopted planar geometry. This increases the quality of the existing hydrogen bonds in almost all cases. Further, detailed analysis of protein-DNA interfacial hydrogen bonds in 107 crystal structures using the in-house program "pyrHBfind" indicates that consideration of energetically more preferred nonplanar amino groups improves the geometry of hydrogen bonds and also gives rise to new contacts amounting to nearly 14.5% of the existing interactions. Large improvements have been observed specifically for the amino groups of guanine, which faces the DNA minor groove and thus helps to resolve the problem of insufficient directional contacts observed in many minor groove binding complexes. Apart from guanine, improvement observed for asparagine, glutamine, adenine, or cytosine also indicates that the consideration of nonplanar amino groups leads to a more realistic scenario of hydrogen bonds occurring between protein and DNA residues.     

Secondary structure sensitivity of hydrogen bond lifetime dynamics in the protein hydration layer

The heterogeneous nature of a protein surface plays an essential role in its biological activity and molecular recognition, and this role is mediated at least partly through the surrounding water molecules. We have performed atomistic molecular dynamics simulations of an aqueous solution of HP-36 to investigate the correlation between the dynamics of the hydration layer water molecules and the lifetimes of protein-water hydrogen bonds. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the quasi-bound states in the surface and quantification of the dynamic equilibrium between quasi-bound and free water molecules in terms of time-dependent rate of interconversion. It is noticed that, irrespective of the structural heterogeneity of different segments of the protein, namely the three alpha-helices, the positively charged amino acid residues form longer-lived hydrogen bonds with water. The overall relaxation behavior of protein-water hydrogen bonds is found to differ significantly among the three helices of the protein. Study of water number density fluctuation reveals that the hydration layer of helix-3 is much less rigid, which can be correlated with faster structural relaxation of the hydrogen bonds between its residues and water. This also agrees excellently with faster translational and rotational motions of water near helix-3, and hence the lower rigidity of its hydration layer. The lower rigidity of the helix-3 hydration layer also correlates well with the biological activity of the protein, as several of the active-site residues of HP-36 are located in helix-3.     

Hydration in drug design. 1. Multiple hydrogen-bonding features of water molecules in mediating protein-ligand interactions

Summary Water is known to play an important rôle in the recognition and stabilization of the interaction between a ligand and its site. This has important implications for drug design. Analyses of 19 high-resolution crystal structures of protein-ligand complexes reveal the multiple hydrogen-bonding feature of water molecules mediating protein-ligand interactions. Most of the water molecules (nearly 80%) involved in bridging the protein and the ligand can make three or more hydrogen bonds when distance and bond angles are used as criteria to define hydrogen-bonding interactions. Isotropic B-factors have been used to take into account the mobility of water molecules. The water molecules at binding sites bridge the protein and ligand, and interact with other water molecules to form a complex network of interconnecting hydrogen bonds. Some water molecules at the site do not directly bridge between the protein and the ligand, but may contribute indirectly to the stability of the complex by holding bridging water molecules in the right position through a network of hydrogen bonds. These water networks are probably crucial for the stability of the protein-ligand complex and are important for any site-directed drug design strategies.     

油纸复合介质中水分子扩散行为的分子动力学模拟

对不同温度下水分子在油纸复合介质中的扩散行为进行了分子动力学模拟研究. 通过分析水分子与纤维素形成的氢键发现, 油中的水分子在模拟过程中会逐渐扩散到纤维素内并与之形成氢键, 而纤维素内的水分子则与纤维素形成氢键后被束缚于纤维素中. 通过分析水分子的扩散系数发现, 由于油和纤维素的极性不同, 使得水分子在油和纤维素两种单介质中的扩散行为有较大差别, 而在复合介质中的扩散系数受水分子在油和纤维素中的比例影响较大, 两者表现出很强的相关性. 水分子和两介质的相互作用与两介质的极性也存在很大的关系, 且不同温度下水分子与两介质的相互作用能和水分子在油和纤维素中的比例也表现出了较强的相关性. 不同温度下水分子的不同分布弱化了温度对扩散系数的影响.     

Structure and dynamics of methanol in water: a quantum mechanical charge field molecular dynamics study

An ab initio quantum mechanical charge field molecular dynamics simulation was carried out for one methanol molecule in water to analyze the structure and dynamics of hydrophobic and hydrophilic groups. It is found that water molecules around the methyl group form a cage-like structure whereas the hydroxyl group acts as both hydrogen bond donor and acceptor, thus forming several hydrogen bonds with water molecules. The dynamic analyses correlate well with the structural data, evaluated by means of radial distribution functions, angular distribution functions, and coordination number distributions. The overall ligand mean residence time, τ identifies the methanol molecule as structure maker. The relative dynamics data of hydrogen bonds between hydroxyl of methanol and water molecules prove the existence of both strong and weak hydrogen bonds. The results obtained from the simulation are in excellent agreement with the experimental results for dilute solution of CH(3)OH in water. The overall hydration shell of methanol consists in average of 18 water molecules out of which three are hydrogen bonded.     

Low-frequency vibrational spectrum of water in the hydration layer of a protein: a molecular dynamics simulation study

An atomistic molecular dynamics simulation has been carried out to understand the low-frequency intermolecular vibrational spectrum of water present in the hydration layer of the protein villin headpiece subdomain or HP-36. An attempt is made to explore how the heterogeneous rigidity of the hydration layers of different segments (three alpha helices) of the protein, strength of the protein-water hydrogen bonds, and their differential relaxation behavior influence the distribution of the intermolecular vibrational density of states of water in the hydration layers. The calculations revealed that compared to bulk water these bands are nonuniformly blue-shifted for water near the helices, the extent of shifts being more pronounced for water molecules hydrogen bonded to the protein residues. It is further noticed that the larger blue shift observed for the water molecules hydrogen bonded to helix 2 residues correlates excellently with the slowest structural relaxation of these hydrogen bonds. These results can be verified by suitable experimental measurements.     

Nonspecific protein-DNA interactions: complexation of alpha-chymotrypsin with a genomic DNA

In this contribution, we report studies on nonspecific protein-DNA interactions of an enzyme protein bovine pancreatic alpha-chymotrypsin (CHT) with genomic DNA (from salmon testes) using two biologically common fluorescent probes: 1-anilinonaphthalene-8-sulfonate (ANS) and 2,6-p-toluidinonaphthalene sulfonate (TNS). TNS molecules that are nonspecifically bound to positively charged basic residues at the surface sites, not in the hydrophobic cavities of the protein, are preferentially displaced upon complexation of TNS-labeled CHT with DNA. The time-resolved fluorescence anisotropy of TNS molecules bound to hydrophobic cavities/clefts of CHT reveals that global tumbling motion of the protein is almost frozen in the protein-DNA complex. A control study on TNS-labeled human serum albumin (HSA) upon interaction with DNA clearly indicates that the ligands in the deep pockets of the protein cannot be displaced by interaction with DNA. We have also found that ANS, which binds to a specific surface site of CHT, is not displaced by DNA. The intactness of the ANS binding in CHT upon complexation with DNA offers the opportunity to measure the distance between the ANS binding site and the contact point of the ethidium bromide (EB)-labeled DNA using the F?rster resonance energy transfer (FRET) technique. Enzymatic activity studies on CHT on a substrate (Ala-Ala-Phe 7-amido-4-methyl coumarin) reveal that the active site of the enzyme remains open for the substrate even in the protein-DNA complex. Circular dichroism (CD) studies on CHT upon complexation with DNA confirm the structural integrity of CHT in the complex. Our studies have attempted to explore an application of nonspecific protein-DNA interactions in the characterization of ligand binding of a protein in solution.     

Dynamics and thermodynamics of water around EcoRI bound to a minimally mutated DNA chain

Water plays an important role in protein-DNA interactions. Here, we examine using molecular dynamics simulations the differences in the dynamic and thermodynamic properties of water in the interfacial and intercalating regions of EcoRI bound to the cognate and to a minimally mutated noncognate DNA chain. The results show that the noncognate complex is not only more hydrated than the cognate complex, but the interfacial waters in the noncognate complex exhibit a faster dynamics, which in turn reduces the hydrogen-bond lifetimes. Thus, the higher hydration, faster reorientation dynamics and faster hydrogen-bond-relaxation times of water, taken together, indicate that, even with a minimal mutation of the DNA sequence, the interfacial regions of the noncognate complex are more poised to allowing the protein to diffuse away than to promoting the formation of a stable complex. Alternatively, the results imply that the slowed water dynamics in the interfacial regions when the protein chances upon a cognate sequence allow the formation of a stable specific protein-DNA complex leading to catalytic action.     

Molecular dynamics simulation of thymine glycol-lesioned DNA reveals specific hydration at the lesion

One nanosecond molecular dynamics (MD) simulation of a thymine glycol (TG)-lesioned part of human lymphoblast AG9387 was performed to determine structural changes in DNA molecule caused by the presence of a lesion. These changes can be significant for proper recognition of lesions by a repair enzyme. Thymine glycol is the DNA oxidative lesion formed by addition of OH radicals to C5 and C6 atoms of the thymine base. This lesion is known as causing Cockayne Syndrome-inherited genetic disorder. Distribution of water molecules in a hydration shell around the DNA molecule was analyzed for its contribution to the recognition of the TG lesion by the repair enzyme. The results of MD simulation show there is a specific DNA structural configuration formed at the lesion. After 500 ps the DNA is bent in a kink at the TG site. This change dislocates the glycosyl bond at C5' to a position closer to the DNA surface, and thus its atoms are more exposed to the surrounding water shell. The increased number of water molecules that are close to the TG site indicates that the glycosyl bond may be easily contacted by the repair enzyme. In addition, the higher number of water molecules at the TG site substantiates the importance of water-mediated hydrogen bonds created between the repair enzyme and the DNA upon formation of the complex. Copyright 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1723-1731, 2001     

IR absorption spectra and association of acridinium salts with hydrogen halide molecules

It was observed that the reaction of acridine with hydrogen halides under heterogeneous conditions proceeds in two steps: An acridinium salt is formed in the first step, while a complex of the salt with hydrogen halide molecules is formed in the second step. The formation of a complex leads to substantial changes in the region of the stretching vibrations of NH⁺ bonds, and this makes it possible to propose the direct participation of these bonds in complexing. The band of the stretching vibrations of NH⁺ bonds of the acridinium ion has a complex structure that changes markedly as a function of the sort of halogen ion and replacement of hydrogen by deuterium and in the case of complexing with the hydrogen halide. It is proposed that the structure of this band is due to Fermi resonance.Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 1, pp. 99–101, January, 1981.     

Molecular dynamics simulations of the hydration of poly(vinyl methyl ether):Hydrogen bonds and quasi-hydrogen bonds

Atomistic detailed hydration structures of poly(vinyl methyl ether)(PVME) have been investigated by molecular dynamics simulations under 300 K at various concentrations. Both radial distribution functions and the distance distributions between donors and acceptors in hydrogen bonds show that the hydrogen bonds between the polymer and water are shorter by 0.005 nm than those between water molecules. The Quasi-hydrogen bonds take only 7.2% of the van der Waals interaction pairs. It was found the hydrogen bonds are not evenly distributed along the polymer chain,and there still exists a significant amount(10%) of ether oxygen atoms that are not hydrogen bonded to water at a concentration as low as 3.3%. This shows that in polymer solutions close contacts occur not only between polymer chains but also between chain segments within the polymer,which leads to inefficient contacts between ether oxygen atoms and water molecules. Variation of the quasi-hydrogen bonds with the concentration is similar to that of hydrogen bonds,but the ratio of the repeat units forming quasi-hydrogen bonds to those forming hydrogen bonds approaches 0.2. A transition was found in the demixing enthalpy at around 30% measured by dynamic testing differential scanning calorimetry(DTDSC) for aqueous solutions of a mono-dispersed low molecular weight PVME,which can be related to the transition of the fractions of hydrogen bonds and quasi-hydrogen bonds at ～27%. The transition of the fractions of hydrogen bonds and quasi-hydrogen bonds at ～27% can be used to explain the demixing enthalpy transition at 30% at a molecular scale. In addition,at the concentration of 86%,each ether oxygen atom bonded with water is assigned 1.56 water molecules on average,and 'free' water molecules emerge at the concentration of around 54%.     

Time‐dependent density functional theory study on the electronic excited‐state geometric structure,infrared spectra,and hydrogen bonding of a doubly hydrogen‐bonded complex

The geometric structures and infrared (IR) spectra in the electronically excited state of a novel doubly hydrogen‐bonded complex formed by fluorenone and alcohols, which has been observed by IR spectra in experimental study, are investigated by the time‐dependent density functional theory (TDDFT) method. The geometric structures and IR spectra in both ground state and the S₁ state of this doubly hydrogen‐bonded FN‐2MeOH complex are calculated using the DFT and TDDFT methods, respectively. Two intermolecular hydrogen bonds are formed between FN and methanol molecules in the doubly hydrogen‐bonded FN‐2MeOH complex. Moreover, the formation of the second intermolecular hydrogen bond can make the first intermolecular hydrogen bond become slightly weak. Furthermore, it is confirmed that the spectral shoulder at around 1700 cm^?1 observed in the IR spectra should be assigned as the doubly hydrogen‐bonded FN‐2MeOH complex from our calculated results. The electronic excited‐state hydrogen bonding dynamics is also studied by monitoring some vibraitonal modes related to the formation of hydrogen bonds in different electronic states. As a result, both the two intermolecular hydrogen bonds are significantly strengthened in the S₁ state of the doubly hydrogen‐bonded FN‐2MeOH complex. The hydrogen bond strengthening in the electronically excited state is similar to the previous study on the singly hydrogen‐bonded FN‐MeOH complex and play important role on the photophysics of fluorenone in solutions. © 2009 Wiley Periodicals, Inc. J Comput Chem 2009     

The effect of aqueous solutions of trimethylamine-N-oxide on pressure induced modifications of hydrophobic interactions

To understand the mechanism of protein protection by the osmolyte trimethylamine-N-oxide (TMAO) at high pressure, using molecular dynamics (MD) simulations, solvation of hydrophobic group is probed in aqueous solutions of TMAO over a wide range of pressures relevant to protein denaturation. The hydrophobic solute considered in this study is neopentane which is a considerably large molecule. The concentrations of TMAO range from 0 to 4 M and for each TMAO concentration, simulations are performed at five different pressures ranging from 1 atm to 8000 atm. Potentials of mean force are calculated and the relative stability of solvent-separated state over the associated state of hydrophobic solute are estimated. Results suggest that high pressure reduces association of hydrophobic solutes. From computations of site-site radial distribution function followed by analysis of coordination number, it is found that water molecules are tightly packed around the nonpolar particle at high pressure and the hydration number increases with increasing pressure. On the other hand, neopentane interacts preferentially with TMAO over water and although hydration of neopentane reduces in presence of this osmolyte, TMAO does not show any tendency to prevent the pressure-induced dispersion of neopentane moieties. It is also observed that TMAO molecules prefer a side-on orientation near the neopentane surface, allowing its oxygen atom to form favorable hydrogen bonds with water while maintaining some hydrophobic contacts with neopentane. Analysis of hydrogen-bond properties and solvation characteristics of TMAO reveals that TMAO can form hydrogen bonds with water and it reduces the identical nearest neighbor water molecules caused by high hydrostatic pressures. Moreover, TMAO enhances life-time of water-water hydrogen bonds and makes these hydrogen bonds more attractive. Implication of these results for counteracting effect of TMAO against protein denaturation at high pressures are discussed.     

How does the first water shell fold proteins so fast?

First shells of hydration and bulk solvent play a crucial role in the folding of proteins. Here, the role of water in the dynamics of proteins has been investigated using a theoretical protein-solvent model and a statistical physics approach. We formulate a hydration model where the hydrogen bonds between water molecules pertaining to the first shell of the protein conformation may be either mainly formed or broken. At thermal equilibrium, hydrogen bonds are formed at low temperature and are broken at high temperature. To explore the solvent effect, we follow the folding of a large sampling of protein chains, using a master-equation evolution. The dynamics shows a clear mechanism. Above the glass-transition temperature, a large ratio of chains fold very rapidly into the native structure irrespective of the temperature, following pathways of high transition rates through structures surrounded by the solvent with broken hydrogen bonds. Although these states have an infinitesimal probability, they act as strong dynamical attractors and fast folding proceeds along these routes rather than pathways with small transition rates between configurations of much higher equilibrium probabilities. At a given low temperature, a broad jump in the folding times is observed. Below this glass temperature, the pathways where hydrogen bonds are mainly formed become those of highest rates although with conformational changes of huge relaxation times. The present results reveal that folding obeys a double-funnel mechanism.     

Structure of bilirubin and its anion according to quantum chemical calculations and molecular dynamics simulation

Quantum chemical calculations of the structural characteristics of the bilirubin molecule and its anion are performed. Intramolecular hydrogen bonds are studied using NBO analysis. It is shown that hydrogen bonds in the bilirubin molecule are nonequivalent, and the bond formed by the keto oxygen of the pyrrole ring and the hydrogen of the carboxyl group belonging to the propionate residue is energetically more favorable. Structural characteristics of the molecular and ion forms of bilirubin in aqueous solution are studied by molecular dynamics simulation. It is found that intermolecular hydrogen bonds with water molecules are formed due to oxygen atoms of the carboxyl group and the keto group of bilirubin, and the probability of their formation by anions is much higher than that for molecules.     

Correlation between the dynamics of hydrogen bonds and the local density reorganization in the protein hydration layer

An atomistic molecular dynamics simulation of the protein villin headpiece subdomain or HP-36 has been carried out with explicit water to explore the microscopic inhomogeneity of local density reorganization of the hydration layers of the three alpha-helical segments of the protein. The density reorganization of the hydration layer of helix-3 is found to occur faster than that for the hydration layers of the other two helices. It is noticed that such inhomogeneous density reorganization at the surface of different secondary structures exhibits excellent correlation with the microscopic dynamics of hydrogen bonds between the protein residues and the hydration water. Further, it is observed that the reorientation of water molecules involved in the formation and breaking of protein-water or water-water hydrogen bonds plays an important role in determining the dynamics of local density of the hydration layer. The faster density reorganization of the hydration layer of helix-3 is also consistent with the functionality of HP-36, as helix-3 contains several active site residues.