Water dynamics in the large cavity of three lipid-binding proteins monitored by (17)O magnetic relaxation dispersion

Intracellular lipid-binding proteins contain a large binding cavity filled with water molecules. The role played by these water molecules in ligand binding is not well understood, but their energetic and dynamic properties must be important for protein function. Here, we use the magnetic relaxation dispersion (MRD) of the water 17O resonance to investigate the water molecules in the binding cavity of three different lipid-binding proteins: heart fatty acid-binding protein (H-FABP), ileal lipid-binding protein (I-LBP) and intestinal fatty acid-binding protein (I-FABP). Whereas about half of the crystallographically visible water molecules appear to be expelled by the ligand, we find that ligand binding actually increases the number of water molecules within the cavity. At 300 K, the water molecules in the cavity exchange positions on a time-scale of about 1ns and exchange with external water on longer time-scales (0.01-1 micros). Exchange of water molecules among hydration sites within the cavity should be strongly coupled to ligand motion. Whereas a recent MD simulation indicates that the structure of the cavity water resembles a bulk water droplet, the present MRD results show that its dynamics is more than two orders of magnitude slower than in the bulk. These findings may have significant implications for the strength, specificity and kinetics of lipid binding.     

Long-residency hydration, cation binding, and dynamics of loop E/helix IV rRNA-L25 protein complex

Molecular dynamics simulations of RNA-protein complex between Escherichia coli loop E/helix IV (LE/HeIV) rRNA and L25 protein reveal a qualitative agreement between the experimental and simulated structures. The major groove of LE is a prominent rRNA cation-binding site. Divalent cations rigidify the LE major groove geometry whereas in the absence of divalent cations LE extensively interacts with monovalent cations via inner-shell binding. The HeIV region shows bistability of its major groove explaining the observed differences between x-ray and NMR structures. In agreement with the experiments, the simulations suggest that helix-alpha1 of L25 is the least stable part of the protein. Inclusion of Mg2+ cations into the simulations causes perturbation of basepairing at the LE/HeIV junction, which does not, however, affect the protein binding. The rRNA-protein complex is mediated by a number of highly specific hydration sites with long-residing water molecules and two of them are bound throughout the entire 24-ns simulation. Long-residing water molecules are seen also outside the RNA-protein contact areas with water-binding times substantially enhanced compared to simulations of free RNA. Long-residency hydration sites thus represent important elements of the three-dimensional structure of rRNA.     

Hydration energy landscape of the active site cavity in cytochrome P450cam

Hydration of protein cavities influences protein stability, dynamics, and function. Protein active sites usually contain water molecules that, upon ligand binding, are either displaced into bulk solvent or retained to mediate protein–ligand interactions. The contribution of water molecules to ligand binding must be accounted for to compute accurate values of binding affinities. This requires estimation of the extent of hydration of the binding site. However, it is often difficult to identify the water molecules involved in the binding process when ligands bind on the surface of a protein. Cytochrome P450cam is, therefore, an ideal model system because its substrate binds in a buried active site, displacing partially disordered solvent, and the protein is well characterized experimentally. We calculated the free energy differences for having five to eight water molecules in the active site cavity of the unliganded enzyme from molecular dynamics simulations by thermodynamic integration employing a three-stage perturbation scheme. The computed free energy differences between the hydration states are small (within 12 kJ mol⁻¹) but distinct. Consistent with the crystallographic determination and studies employing hydrostatic pressure, we calculated that, although ten water molecules could in principle occupy the volume of the active site, occupation by five to six water molecules is thermodynamically most favorable. Proteins 32:381–396, 1998. © 1998 Wiley-Liss, Inc.     

Simulations of apo and holo-fatty acid binding protein: structure and dynamics of protein, ligand and internal water.

Two molecular dynamics simulations of 5 ns each have been carried out for rat intestinal fatty acid binding protein, in apo-form and with bound palmitate. The fatty acid and a number of water molecules are encapsulated in a large interior cavity of the barrel-shaped protein. The simulations are compared to experimental data and analyzed in terms of root mean square deviations, atomic B-factors, secondary structure elements, hydrogen bond patterns, and distance constraints derived from nuclear Overhauser experiments. Excellent agreement is found between simulated and experimental solution structures of holo-FABP, but a number of differences are observed for the apo-form. The ligand in holo-FABP shows considerable displacement after about 1.5 ns and displays significant configurational entropy. A novel computational approach has been employed to identify internal water and to capture exchange pathways. Orifices in the portal and gap regions of the protein, discussed in the experimental literature, have been confirmed as major openings for solvent exchange between the internal cavity and bulk water. A third opening on the opposite side of the barrel experiences significant exchange but it does not provide a pathway for further passage to the central cavity. Internal water is characterized in terms of density distributions, interaction energies, mobility, protein contact times, and water molecule coordination. A number of differences are observed between the apo and holo-forms and related to differences in the protein structure. Solvent inside apo-FABP, for example, shows characteristics of a water droplet, while solvent in holo-FABP benefits from interactions with the ligand headgroup and slightly stronger interactions with protein residues.     

Investigations of peptide hydration using NMR and molecular dynamics simulations: A study of effects of water on the conformation and dynamics of antamanide

Summary The influence of water binding on the conformational dynamics of the cyclic decapeptide antamanide dissolved in the model lipophilic environment chloroform is investigated by NMR relaxation measurements. The water-peptide complex has a lifetime of 35 s at 250 K, which is longer than typical lifetimes of water-peptide complexes reported in aqueous solution. In addition, there is a rapid intracomplex mobility that probably involves librational motions of the bound water or water molecules hopping between different binding sites. Water binding restricts the flexibility of antamanide. The experimental findings are compared with GROMOS molecular dynamics simulations of antamanide with up to eight bound water molecules. Within the simulation time of 600 ps, no water molecule leaves the complex. Additionally, the simulations show a reduced flexibility for the complex in comparison with uncomplexed antamanide. Thus, there is a qualitative agreement between the experimental NMR results and the computer simulations.     

Investigation of glucose binding sites on insulin

Possible insulin binding sites for D-glucose have been investigated theoretically by docking and molecular dynamics (MD) simulations. Two different docking programs for small molecules were used; Multiple Copy Simultaneous Search (MCSS) and Solvation Energy for Exhaustive Docking (SEED) programs. The configurations resulting from the MCSS search were evaluated with a scoring function developed to estimate the binding free energy. SEED calculations were performed using various values for the dielectric constant of the solute. It is found that scores emphasizing non-polar interactions gave a preferential binding site in agreement with that inferred from recent fluorescence and NMR NOESY experiments. The calculated binding affinity of -1.4 to -3.5 kcal/mol is within the measured range of -2.0 +/- 0.5 kcal/mol. The validity of the binding site is suggested by the dynamical stability of the bound glucose when examined with MD simulations with explicit solvent. Alternative binding sites were found in the simulations and their relative stabilities were estimated. The motions of the bound glucose during molecular dynamics simulations are correlated with the motions of the insulin side chains that are in contact with it and with larger scale insulin motions. These results raise the question of whether glucose binding to insulin could play a role in its activity. The results establish the complementarity of molecular dynamics simulations and normal mode analyses with the search for binding sites proposed with small molecule docking programs.     

The carbohydrate-binding site in galectin-3 is preorganized to recognize a sugarlike framework of oxygens: ultra-high-resolution structures and water dynamics

The recognition of carbohydrates by proteins is a fundamental aspect of communication within and between living cells. Understanding the molecular basis of carbohydrate-protein interactions is a prerequisite for the rational design of synthetic ligands. Here we report the high- to ultra-high-resolution crystal structures of the carbohydrate recognition domain of galectin-3 (Gal3C) in the ligand-free state (1.08 ? at 100 K, 1.25 ? at 298 K) and in complex with lactose (0.86 ?) or glycerol (0.9 ?). These structures reveal striking similarities in the positions of water and carbohydrate oxygen atoms in all three states, indicating that the binding site of Gal3C is preorganized to coordinate oxygen atoms in an arrangement that is nearly optimal for the recognition of β-galactosides. Deuterium nuclear magnetic resonance (NMR) relaxation dispersion experiments and molecular dynamics simulations demonstrate that all water molecules in the lactose-binding site exchange with bulk water on a time scale of nanoseconds or shorter. Nevertheless, molecular dynamics simulations identify transient water binding at sites that agree well with those observed by crystallography, indicating that the energy landscape of the binding site is maintained in solution. All heavy atoms of glycerol are positioned like the corresponding atoms of lactose in the Gal3C complexes. However, binding of glycerol to Gal3C is insignificant in solution at room temperature, as monitored by NMR spectroscopy or isothermal titration calorimetry under conditions where lactose binding is readily detected. These observations make a case for protein cryo-crystallography as a valuable screening method in fragment-based drug discovery and further suggest that identification of water sites might inform inhibitor design.     

Water in the active site of ketosteroid isomerase

Classical molecular dynamics simulations were utilized to investigate the structural and dynamical properties of water in the active site of ketosteroid isomerase (KSI) to provide insight into the role of these water molecules in the enzyme-catalyzed reaction. This reaction is thought to proceed via a dienolate intermediate that is stabilized by hydrogen bonding with residues Tyr16 and Asp103. A comparative study was performed for the wild-type (WT) KSI and the Y16F, Y16S, and Y16F/Y32F/Y57F (FFF) mutants. These systems were studied with three different bound ligands: equilenin, which is an intermediate analog, and the intermediate states of two steroid substrates. Several distinct water occupation sites were identified in the active site of KSI for the WT and mutant systems. Three additional sites were identified in the Y16S mutant that were not occupied in WT KSI or the other mutants studied. The number of water molecules directly hydrogen bonded to the ligand oxygen was approximately two in the Y16S mutant and one in the Y16F and FFF mutants, with intermittent hydrogen bonding of one water molecule in WT KSI. The molecular dynamics trajectories of the Y16F and FFF mutants reproduced the small conformational changes of residue 16 observed in the crystal structures of these two mutants. Quantum mechanical/molecular mechanical calculations of (1)H NMR chemical shifts of the protons in the active site hydrogen-bonding network suggest that the presence of water in the active site does not prevent the formation of short hydrogen bonds with far-downfield chemical shifts. The molecular dynamics simulations indicate that the active site water molecules exchange much more frequently for WT KSI and the FFF mutant than for the Y16F and Y16S mutants. This difference is most likely due to the hydrogen-bonding interaction between Tyr57 and an active site water molecule that is persistent in the Y16F and Y16S mutants but absent in the FFF mutant and significantly less probable in WT KSI.     

A nanosecond molecular dynamics study of antiparallel d(G)7 quadruplex structures: effect of the coordinated cations

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G7) quadruplex structures with different coordinated ions, namely Na+ and K+ ion, water and Na+ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of approximately 1.5 A from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K+ ions are quite immobile. MD studies at 400K indicate that K+ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K+ ion coordinated structure is more favourable than that for the Na+ ion coordinated antiparallel quadruplex structure.     

Analysis of the bacterial luciferase mobile loop by replica-exchange molecular dynamics

Bacterial luciferase contains an extended 29-residue mobile loop. Movements of this loop are governed by binding of either flavin mononucleotide (FMNH₂) or polyvalent anions. To understand this process, loop dynamics were investigated using replica-exchange molecular dynamics that yielded conformational ensembles in either the presence or absence of FMNH₂. The resulting data were analyzed using clustering and network analysis. We observed the closed conformations that are visited only in the simulations with the ligand. Yet the mobile loop is intrinsically flexible, and FMNH₂ binding modifies the relative populations of conformations. This model provides unique information regarding the function of a crystallographically disordered segment of the loop near the binding site. Structures at or near the fringe of this network were compatible with flavin binding or release. Finally, we demonstrate that the crystallographically observed conformation of the mobile loop bound to oxidized flavin was influenced by crystal packing. Thus, our study has revealed what we believe are novel conformations of the mobile loop and additional context for experimentally determined structures.     

CO migration pathways in cytochrome P450cam studied by molecular dynamics simulations

Previous laser flash photolysis investigations between 100 and 300 K have shown that the kinetics of CO rebinding with cytochrome P450(cam)(camphor) consist of up to four different processes revealing a complex internal dynamics after ligand dissociation. In the present work, molecular dynamics simulations were undertaken on the ternary complex P450(cam)(cam)(CO) to explore the CO migration pathways, monitor the internal cavities of the protein, and localize the CO docking sites. One trajectory of 1 nsec with the protein in a water box and 36 trajectories of 1 nsec in the vacuum were calculated. In each trajectory, the protein contained only one CO ligand on which no constraints were applied. The simulations were performed at 200, 300, and 320 K. The results indicate the presence of seven CO docking sites, mainly hydrophobic, located in the same moiety of the protein. Two of them coincide with xenon binding sites identified by crystallography. The protein matrix exhibits eight persistent internal cavities, four of which corresponding to the ligand docking sites. In addition, it was observed that water molecules entering the protein were mainly attracted into the polar pockets, far away from the CO docking sites. Finally, the identified CO migration pathways provide a consistent interpretation of the experimental rebinding kinetics.     

Mapping the binding sites of challenging drug targets

An increasing number of medically important proteins are challenging drug targets because their binding sites are too shallow or too polar, are cryptic and thus not detectable without a bound ligand or located in a protein–protein interface. While such proteins may not bind druglike small molecules with sufficiently high affinity, they are frequently druggable using novel therapeutic modalities. The need for such modalities can be determined by experimental or computational fragment based methods. Computational mapping by mixed solvent molecular dynamics simulations or the FTMap server can be used to determine binding hot spots. The strength and location of the hot spots provide very useful information for selecting potentially successful approaches to drug discovery.     

Crystallographic and computational studies on 4-phenyl-N-(beta-D-glucopyranosyl)-1H-1,2,3-triazole-1-acetamide, an inhibitor of glycogen phosphorylase: comparison with alpha-D-glucose, N-acetyl-beta-D-glucopyranosylamine and N-benzoyl-N'-beta-D-glucopyranosyl urea binding

4-Phenyl-N-(beta-D-glucopyranosyl)-1H-1,2,3-triazole-1-acetamide (glucosyltriazolylacetamide) has been studied in kinetic and crystallographic experiments with glycogen phosphorylase b (GPb), in an effort to utilize its potential as a lead for the design of potent antihyperglycaemic agents. Docking and molecular dynamics (MD) calculations have been used to monitor more closely the binding modes in operation and compare the results with experiment. Kinetic experiments in the direction of glycogen synthesis showed that glucosyltriazolylacetamide is a better inhibitor (K(i) = 0.18 mM) than the parent compound alpha-D-glucose (K(i) = 1.7 mM) or beta-D-glucose (K(i) = 7.4 mM) but less potent inhibitor than the lead compound N-acetyl-beta-D-glucopyranosylamine (K(i) = 32 microM). To elucidate the molecular basis underlying the inhibition of the newly identified compound, we determined the structure of GPb in complex with glucosyltriazolylacetamide at 100 K to 1.88 A resolution, and the structure of the compound in the free form. Glucosyltriazolylacetamide is accommodated in the catalytic site of the enzyme and the glucopyranose interacts in a manner similar to that observed in the GPb-alpha-D-glucose complex, while the substituent group in the beta-position of the C1 atom makes additional hydrogen bonding and van der Waals interactions to the protein. A bifurcated donor type hydrogen bonding involving O3H, N3, and N4 is seen as an important structural motif strengthening the binding of glucosyltriazolylacetamide with GP which necessitated change in the torsion about C8-N2 bond by about 62 degrees going from its free to the complex form with GPb. On binding to GP, glucosyltriazolylacetamide induces significant conformational changes in the vicinity of this site. Specifically, the 280s loop (residues 282-288) shifts 0.7 to 3.1 A (CA atoms) to accommodate glucosyltriazolylacetamide. These conformational changes do not lead to increased contacts between the inhibitor and the protein that would improve ligand binding compared with the lead compound. In the molecular modeling calculations, the GOLD docking runs with and without the crystallographic ordered cavity waters using the GoldScore scoring function, and without cavity waters using the ChemScore scoring function successfully reproduced the crystallographic binding conformation. However, the GLIDE docking calculations both with (GLIDE XP) and without (GLIDE SP and XP) the cavity water molecules were, impressively, further able to accurately reproduce the finer details of the GPb-glucosyltriazolylacetamide complex structure. The importance of cavity waters in flexible receptor MD calculations compared to "rigid" (docking) is analyzed and highlighted, while in the MD itself very little conformational flexibility of the glucosyltriazolylacetamide ligand was observed over the time scale of the simulations.     

Conformational and dynamics changes induced by bile acids binding to chicken liver bile acid binding protein

The correlation between protein motions and function is a central problem in protein science. Several studies have demonstrated that ligand binding and protein dynamics are strongly correlated in intracellular lipid binding proteins (iLBPs), in which the high degree of flexibility, principally occurring at the level of helix-II, CD, and EF loops (the so-called portal area), is significantly reduced upon ligand binding. We have recently investigated by NMR the dynamic properties of a member of the iLBP family, chicken liver bile acid binding protein (cL-BABP), in its apo and holo form, as a complex with two bile salts molecules. Binding was found to be regulated by a dynamic process and a conformational rearrangement was associated with this event. We report here the results of molecular dynamics (MD) simulations performed on apo and holo cL-BABP with the aim of further characterizing the protein regions involved in motion propagation and of evaluating the main molecular interactions stabilizing bound ligands. Upon binding, the root mean square fluctuation values substantially decrease for CD and EF loops while increase for the helix-loop-helix region, thus indicating that the portal area is the region mostly affected by complex formation. These results nicely correlate with backbone dynamics data derived from NMR experiments. Essential dynamics analysis of the MD trajectories indicates that the major concerted motions involve the three contiguous structural elements of the portal area, which however are dynamically coupled in different ways whether in the presence or in the absence of the ligands. Motions of the EF loop and of the helical region are part of the essential space of both apo and holo-BABP and sample a much wider conformational space in the apo form. Together with NMR results, these data support the view that, in the apo protein, the flexible EF loop visits many conformational states including those typical of the holo state and that the ligand acts stabilizing one of these pre-existing conformations. The present results, in agreement with data reported for other iLBPs, sharpen our knowledge on the binding mechanism for this protein family.     

Anesthetic interaction with ketosteroid isomerase: insights from molecular dynamics simulations

The nature and the sites of interactions between anesthetic halothane and homodimeric Delta5-3-ketosteroid isomerase (KSI) are characterized by flexible ligand docking and confirmed by 1H-15N NMR. The dynamics consequence of halothane interaction and the implication of the dynamic changes to KSI function are studied by multiple 5-ns molecular dynamics simulations in the presence and absence of halothane. Both docking and MD simulations show that halothane prefer the amphiphilic dimeric interface to the hydrophobic active site of KSI. Halothane occupancy at the dimer interface disrupted the intersubunit hydrogen bonding formed either directly through side chains of polar residues or indirectly through the mediation of the interfacial water molecules. Moreover, in the presence of halothane, the exchange rate of the bound waters with bulk water was increased. Halothane perturbation to the dimer interface affected the overall flexibility of the active site. This action is likely to contribute to the halothane-induced reduction of the KSI activity. The allosteric halothane modulation of the dynamics-function relationship of KSI without direct competition at the enzymatic active sites may be generalized to offer a unifying explanation of anesthetic action on a diverse range of multidomain neuronal proteins that are potentially relevant to clinical general anesthesia.     

Water and potassium dynamics inside the KcsA K(+) channel

Molecular dynamics simulations and electrostatic modeling are used to investigate structural and dynamical properties of the potassium ions and of water molecules inside the KcsA channel immersed in a membrane-mimetic environment. Two potassium ions, initially located in the selectivity filter binding sites, maintain their position during 2 ns of dynamics. A third potassium ion is very mobile in the water-filled cavity. The protein appears engineered so as to polarize water molecules inside the channel cavity. The resulting water induced dipole and the positively charged potassium ion within the cavity are the key ingredients for stabilizing the two K(+) ions in the binding sites. These two ions experience single file movements upon removal of the potassium in the cavity, confirming the role of the latter in ion transport through the channel.     

Effects of ligand binding on the dynamics of rice nonspecific lipid transfer protein 1: a model from molecular simulations

Plant nonspecific lipid transfer proteins (nsLTPs) are small, basic proteins constituted mainly of alpha-helices and stabilized by four conserved disulfide bridges. They are characterized by the presence of a tunnel-like hydrophobic cavity, capable of transferring various lipid molecules between lipid bilayers in vitro. In this study, molecular dynamics (MD) simulations were performed at room temperature to investigate the effects of lipid binding on the dynamic properties of rice nsLTP1. Rice nsLTP1, either in the free form or complexed with one or two lipids was subjected to MD simulations. The C-terminal loop was very flexible both before and after lipid binding, as revealed by calculating the root-mean-square fluctuation. After lipid binding, the flexibility of some residues that were not in direct contact with lipid molecules increased significantly, indicating an increase of entropy in the region distal from the binding site. Essential dynamics analysis revealed clear differences in motion between unliganded and liganded rice nsLTP1s. In the free form of rice nsLTP1, loop1 exhibited the largest directional motion. This specific essential motion mode diminished after binding one or two lipid molecules. To verify the origin of the essential motion observed in the free form of rice nsLTP1, we performed multiple sequence alignments to probe the intrinsic motion encoded in the primary sequence. We found that the amino acid sequence of loop1 is highly conserved among plant nsLTP1s, thus revealing its functional importance during evolution. Furthermore, the sequence of loop1 is composed mainly of amino acids with short side chains. In this study, we show that MD simulations, together with essential dynamics analysis, can be used to determine structural and dynamic differences of rice nsLTP1 upon lipid binding.     

Differences between apo and three holo forms of the intestinal fatty acid binding protein seen by molecular dynamics computer calculations

It is commonly believed that binding affinity can be estimated by consideration of local changes of ligand and protein. This paper discusses a set of molecular dynamics simulations of intestinal fatty acid binding protein addressing the protein's response to presence or absence of different ligands. A 5-ns simulation was performed of the protein without a ligand, and three simulations (one 5-ns and two 2-ns) were performed with different fatty acids bound. The results indicate that, although the basic protein structure is unchanged by the presence of the ligand, other properties are significantly affected by ligand binding. For example, zero-time covariance patterns between protein, bound waters, and ligand vary between the different simulations. Moreover, the interaction energies between ligand and specific residues indicate that different ligands are stabilized in different ways. In sum, the results suggest that binding thermodynamics within this system will need to be calculated not from a subset of nearby protein:ligand interactions, but will depend on a knowledge of the motions coupling together water, protein, and ligand.     

Protein grabs a ligand by extending anchor residues: molecular simulation for Ca2+ binding to calmodulin loop

The structural difference in proteins between unbound and bound forms directly suggests the importance of the conformational plasticity of proteins. However, pathways that connect two-end structures and how they are coupled to the binding reaction are not well understood at atomic resolution. Here, we analyzed the free-energy landscape, explicitly taking into account coupling between binding and conformational change by performing atomistic molecular dynamics simulations for Ca2+ binding to a calmodulin loop. Using the AMBER force field with explicit water solvent, we conducted umbrella sampling for the free-energy surface and steered molecular dynamics for the pathway search. We found that, at an early stage of binding, some key residue side chains extend their "arms" to catch Ca2+ and, after catching, they carry the Ca2+ to the center of the binding pocket. This grabbing motion resulted in smooth and stepwise exchange in coordination partners of Ca2+ from water oxygen to atoms in the calmodulin loop. The key residue that first caught the ion was one of the two acidic residues, which are highly conserved. In the pathway simulations, different pathways were observed between binding and dissociation reactions: The former was more diverse than the latter.