Anti-insulin antibody structure and conformation. II. Molecular dynamics with explicit solvent.

Molecular dynamics at 300 K was used as a conformation searching tool to analyze a knowledge-based structure prediction of an anti-insulin antibody. Solvation effects were modeled by packing water molecules around the antigen binding loops. Some loops underwent backbone and side-chain conformational changes during the 95-ps equilibration, and most of these new, lower potential energy conformations were stable during the subsequent 200-ps simulation. Alterations to the model include changes in the intraloop, main-chain hydrogen bonding network of loop H3, and adjustments of Tyr and Lys side chains of H3 induced by hydrogen bonding to water molecules. The structures observed during molecular dynamics support the conclusion of the previous paper that hydrogen bonding will play the dominant role in antibody-insulin recognition. Determination of the structure of the antibody by x-ray crystallography is currently being pursued to provide an experimental test of these results. The simulation appears to improve the model, but longer simulations at higher temperatures should be performed.     

Localization of bound water in the solution structure of the immunoglobulin binding domain of streptococcal protein G. Evidence for solvent-induced helical distortion in solution.

The presence of bound water in the solution structure of the IgG binding domain of streptococcal protein G has been investigated by nuclear magnetic resonance using three-dimensional 1H rotating frame Overhauser 1H-15N multiple quantum coherence spectroscopy. The backbone amide protons of three residues, Ala20, Gln32 and Tyr33, are found to be in close proximity to bound water. Examination of the three-dimensional structure of the IgG binding domain indicates that in the vicinity of these three residues there are no backbone groups that do not already participate in hydrogen bonding and there are no suitably placed side-chain groups available for hydrogen bonding with water. As the lifetime of the bound water detected in this nuclear magnetic resonance experiment is greater than about one nanosecond, it is likely that the two bound water molecules participate in a bifurcating hydrogen bonding network comprising a CO-NH hydrogen bonded pair, such that the water molecule accepts a hydrogen bond from the NH proton and donates one to the carbonyl oxygen with the result that the amide proton is involved in a three center hydrogen bond. On the basis of the structure, one water molecule participates in such an interaction with the Ala20(NH)-Met1(CO) hydrogen bonded pair at the beginning of an anti-parallel beta-sheet, and the other with the Tyr33(NH)-Val29(CO) hydrogen bonded pair in the single alpha-helix. The latter, which is external and solvent accessible, is associated with a distortion in the alpha-helix centered around Tyr33 which consists of a significant increase in the CO(i-4)-N(i) and CO(i-4)-NH(i) distances relative to those in the rest of the helix, as well as a significant departure in the phi, psi angles of Tyr33 relative to regular helical geometry. Such solvent induced distortions in alpha-helices have been previously noticed in crystal structures and were postulated as possible folding intermediates for helical structures. The present observation of this phenomenon in solution indicates, however, that these water molecules are tightly bound and represent an integral part of the protein framework.     

The Ribbon of Hydrogen Bonds in the Three-Dimensional Structure of Globular Proteins

We report quantum mechanical computations and experimental evidence which suggest that the backbone conformation of globular proteins depends generally on the conservation of that part of the hydrogen bond network or ribbon which is joined, in general, directly to the backbone and is largely independent of the remainder of this whole network of hydrogen bonds. The familiar hydrogen bonds of the helix and the sheet form about one-half of this ribbon of hydrogen bonds. Both water molecules and hydrogen bonding side chain groups are involved in the formation of the ribbon.This view of the three-dimensional structure of globular proteins in terms of the `molecule' allows us to deal with the non-secondary structure as well as with the familiar secondary structure. It also suggests that the ribbon contains approximately the same number of hydrogen bonds within all three structures – the helix, the sheet and the coil – and that this is the reason for the ease of interconversion of these three structures.The quantum mechanical computations on hydrogen bonding suggest that delocalised water molecules which have substantial mobility are an essential part of the ribbon. This situation arises because the hydrogen bonding groups of the protein molecule are not free to move to optimise the hydrogen bonding geometries as are the oxygen atoms in the waters and ices. Such delocalised water molecules either have high B values or are invisible in the X-ray data and yet are able to form a structure which is as strong as a normal hydrogen bond.The experimental data on the point mutations of the THRI57 residue of the T4 phage lysome provides an initial test of this model. Both the local backbone conformation and the ribbon of hydrogen bonds are conserved throughout all the mutations of residue 157,providing that the delocalised water molecules are accepted as a genuine part of the structure. These mutations include the introduction of hydrocarbon side chains at position 157 when water molecules or other side chain groups take over the formation of the hydrogen bonds.We suggest that, provided steric effects are not important, many point mutations succeed because they leave the ribbon of hydrogen bonds (and so the backbone conformation) largely unchanged.     

Peptide folding driven by Van der Waals interactions

Contrary to the widespread view that hydrogen bonding and its entropy effect play a dominant role in protein folding, folding into helical and hairpin-like structures is observed in molecular dynamics (MD) simulations without hydrogen bonding in the peptide-solvent system. In the widely used point charge model, hydrogen bonding is calculated as part of the interaction between atomic partial charges. It is removed from these simulations by setting atomic charges of the peptide and water to zero. Because of the structural difference between the peptide and water, van der Waals (VDW) interactions favor peptide intramolecular interactions and are a major contributing factor to the structural compactness. These compact structures are amino acid sequence dependent and closely resemble standard secondary structures, as a consequence of VDW interactions and covalent bonding constraints. Hydrogen bonding is a short range interaction and it locks the approximate structure into the specific secondary structure when it is included in the simulation. In contrast to standard molecular simulations where the total energy is dominated by charge-charge interactions, these simulation results will give us a new view of the folding mechanism.     

Crystal structure of N-acetylornithine transcarbamylase from Xanthomonas campestris: a novel enzyme in a new arginine biosynthetic pathway found in several eubacteria

We have identified in Xanthomonas campestris a novel N-acetylornithine transcarbamylase that replaces ornithine transcarbamylase in the canonic arginine biosynthetic pathway of several Eubacteria. The crystal structures of the protein in the presence and absence of the reaction product, N-acetylcitrulline, were determined. This new family of transcarbamylases lacks the DxxSMG motif that is characteristic of all ornithine transcarbamylases (OTCases) and contains a novel proline-rich loop that forms part of the active site. The specificity for N-acetylornithine is conferred by hydrogen bonding with residues in the proline-rich loop via water molecules and by hydrophobic interactions with residues from the adjacent 80's, 120's, and proline-rich loops. This novel protein structure provides a starting point for rational design of specific analogs that may be useful in combating human and plant pathogens that utilize acetylornithine transcarbamylase rather than ornithine transcarbamylase.     

MTH1 Substrate Recognition—An Example of Specific Promiscuity

MTH1 (NUDT1) is an oncologic target involved in the prevention of DNA damage. We investigate the way MTH1 recognises its substrates and present substrate-bound structures of MTH1 for 8-oxo-dGTP and 8-oxo-rATP as examples of novel strong and weak binding substrate motifs. Investigation of a small set of purine-like fragments using 2D NMR resulted in identification of a fragment with weak potency. The protein-ligand X-Ray structure of this fragment provides insight into the role of water molecules in substrate selectivity. Wider fragment screening by NMR resulted in three new protein structures exhibiting alternative binding configurations to the key Asp-Asp recognition element of the protein. These inhibitor binding modes demonstrate that MTH1 employs an intricate yet promiscuous mechanism of substrate anchoring through its Asp-Asp pharmacophore. The structures suggest that water-mediated interactions convey selectivity towards oxidized substrates over their non-oxidised counterparts, in particular by stabilization of a water molecule in a hydrophobic environment through hydrogen bonding. These findings may be useful in the design of inhibitors of MTH1.     

Analysis of protein main-chain solvation as a function of secondary structure

We have analysed the hydration of main-chain carbonyl and amide groups in 24 high-resolution well-refined protein structures as a function of the secondary structure in which these polar groups occur. We find that main-chain atoms in beta-sheets are as hydrated as those in alpha-helices, with most interactions involving "free" amide and carbonyl groups that do not participate in secondary structure hydrogen bonds. The distributions of water molecules around these non-bonded carbonyl groups reflect specific steric interactions due to the local secondary structure. Approximately 20% and 4%, respectively of bonded carbonyl and amide groups interact with solvent. These include interactions with carbonyl groups on the exposed faces of alpha-helices that have been correlated previously with bending of the helix. Water molecules interacting with alpha-helices occur mainly at the amino and carbonyl termini of the helices, in which case the solvent sites maintain the hydrogen bonding by bridging between residues i and i-3 or i-4 at the amino terminus and between i and i+3 or i+4 at the carbonyl terminus. We also see a number of solvent-mediated Ncap and Ccap interactions. The water molecules interacting with beta-sheets occur mainly at the edges, in which case they extend the sheet structure, or at the ends of strands, in which case they extend the beta-ladder. In summary, the solvent networks appear to extend the hydrogen-bonding structure of the secondary structures. In beta-turns, which usually occur at the surface of a protein, exposed amide and carbonyl groups are often hydrated, especially close to glycine residues. Occasionally water molecules form a bridge between residues i and i+3 in the turn and this may provide extra stabilization.     

High-resolution structures of scytalone dehydratase-inhibitor complexes crystallized at physiological pH

Scytalone dehydratase is a molecular target of inhibitor design efforts aimed at preventing the fungal disease caused by Magnaporthe grisea. A method for cocrystallization of enzyme with inhibitors at neutral pH has produced several crystal structures of enzyme-inhibitor complexes at resolutions ranging from 1.5 to 2.2 A. Four high resolution structures of different enzyme-inhibitor complexes are described. In contrast to the original X-ray structure of the enzyme, the four new structures have well-defined electron density for the loop region comprising residues 115-119 and a different conformation between residues 154 and 160. The structure of the enzyme complex with an aminoquinazoline inhibitor showed that the inhibitor is in a position to form a hydrogen bond with the amide of the Asn131 side chain and with two water molecules in a fashion similar to the salicylamide inhibitor in the original structure, thus confirming design principles. The aminoquinazoline structure also allows for a more confident assignment of donors and acceptors in the hydrogen bonding network. The structures of the enzyme complexes with two dichlorocyclopropane carboxamide inhibitors showed the two chlorine atoms nearly in plane with the amide side chain of Asn131. The positions of Phe53 and Phe158 are significantly altered in the new structures in comparison to the two structures obtained from crystals grown at acidic pH. The multiple structures help define the mobility of active site amino acids critical for catalysis and inhibitor binding.     

Do water molecules mediate protein-DNA recognition?

A comprehensive analysis of interfacial water molecules in the structures of 109 unique protein-DNA complexes is presented together with a new view on their role in protein-DNA recognition. Location of interfacial water molecules as reported in the crystal structures and as emerging from a series of molecular dynamics studies on protein-DNA complexes with explicit solvent and counterions, was analyzed based on their acceptor, donor hydrogen bond relationships with the atoms and residues of the macromolecules, electrostatic field calculations and packing density considerations. Water molecules for the purpose of this study have been categorized into four classes: viz. (I) those that contact both the protein and the DNA simultaneously and thus mediate recognition directly; (II) those that contact either the protein or the DNA exclusively via hydrogen bonds solvating each solute separately; (III) those that contact the hydrophobic groups in either the protein or the DNA; and, lastly (IV) those that contact another water molecule. Of the 17,963 crystallographic water molecules under examination, about 6% belong to class I and 76% belong to class II. About three-fourths of class I and class II water molecules are exclusively associated with hydrogen bond acceptor atoms of both protein and DNA. Noting that DNA is polyanionic, it is significant that a majority of the crystallographically observed water molecules as well as those from molecular dynamics simulations should be involved in facilitating binding by screening unfavorable electrostatics. Less than 2% of the reported water molecules occur between hydrogen bond donor atoms of protein and acceptor atoms of DNA. These represent cases where protein atoms cannot reach out to DNA to make favorable hydrogen bond interactions due to packing/structural restrictions and interfacial water molecules provide an extension to side-chains to accomplish hydrogen bonding.     

Mapping the energetics of water-protein and water-ligand interactions with the "natural" HINT forcefield: predictive tools for characterizing the roles of water in biomolecules

The energetics and hydrogen bonding pattern of water molecules bound to proteins were mapped by analyzing structural data (resolution better than 2.3A) for sets of uncomplexed and ligand-complexed proteins. Water-protein and water-ligand interactions were evaluated using hydropatic interactions (HINT), a non-Newtonian forcefield based on experimentally determined logP(octanol/water) values. Potential water hydrogen bonding ability was assessed by a new Rank algorithm. The HINT-derived binding energies and Ranks for second shell water molecules were -0.04 kcal mol(-1) and 0.0, respectively, for first shell water molecules -0.38 kcal mol(-1) and 1.6, for active site water molecules -0.45 kcal mol(-1) and 2.3, for cavity water molecules -0.55 kcal mol(-1) and 3.3, and for buried water molecules -0.56 kcal mol(-1) and 4.4. For the last four classes, similar energies indicate that internal and external water molecules interact with protein almost equally, despite different degrees of hydrogen bonding. The binding energies and Ranks for water molecules bridging ligand-protein were -1.13 kcal mol(-1) and 4.5, respectively. This energetic contribution is shared equally between protein and ligand, whereas Rank favors the protein. Lastly, by comparing the uncomplexed and complexed forms of proteins, guidelines were developed for prediction of the roles played by active site water molecules in ligand binding. A water molecule with high Rank and HINT score is unlikely to make further interactions with the ligand and is largely irrelevant to the binding process, while a water molecule with moderate Rank and high HINT score is available for ligand interaction. Water molecule displaced for steric reasons were characterized by lower Rank and HINT score. These guidelines, tested by calculating HINT score and Rank for 50 water molecules bound in the active site of four uncomplexed proteins (for which the structures of the liganded forms were also available), correctly predicted the ultimate roles (in the complex) for 76% of water molecules. Some failures were likely due to ambiguities in the structural data.     

Large-scale networks of hydration water molecules around proteins investigated by cryogenic X-ray crystallography.

The structures at protein-water interface, i.e. the hydration structure of proteins, have been investigated by cryogenic X-ray crystal structure analyses. Hydration structures appeared far clearer at cryogenic temperature than at ambient temperature, presumably because the motions of hydration water molecules were quenched by cooling. Based on the structural models obtained, the hydration structures were systematically analyzed with respect to the amount of water molecules, the interaction modes between water molecules and proteins, the local and the global distribution of them on the surface of proteins. The standard tetrahedral interaction geometry of water in bulk retained at the interface and enabled the three-dimensional chain connection of hydrogen bonds between hydration water molecules and polar protein atoms. Large-scale networks of hydrogen bonds covering the entire surface of proteins were quite flexible to accommodate to the large-scale conformational changes of proteins and seemed to have great influences on the dynamics and function of proteins. The present observation may provide a new concept for discussing the dynamics of proteins in aqueous solution.     

Water structure in [Phe4 Val6] antamanide X 12H2O crystallized from dioxane

Crystals of [Phe4 Val6] antamanide (cyclic [ValProProPhePhe]2) grown from dioxane/H2O, with space group P21212 and cell parameters a = 15.099(4), b = 22.008(5) and c = 11.024(3) A, are almost identical to crystals grown from H2O/acetone, the structure of which was determined a number of years ago. Per peptide molecule there are the equivalent of 12 water molecules occupying 16 sites in both crystals; however, in the new investigation a number of water molecules present at one-half occupancy have been found in different positions than in the earlier analysis. The interpretation of the hydrogen bonding between peptide/water and between water/water is much more satisfactory. Pentagonal water assemblies are present in the solvent channel. There is a distinct indication of the occurrence of a bifurcated bond between two water molecules, as well as the presence of three-center hydrogen bonds joining three water molecules. This may be the first experimental example of a bifurcated bond between two water molecules.     

Neutron Diffraction of Alpha,Beta and Gamma Cyclodextrins: Hydrogen Bonding Patterns

Abstract

Cyclodextrins (CD's) have proved useful as model systems for the study of hydrogen bonding. They are torus-shaped molecules composed of six(α), seven (β) or eight(γ) (1?4) linked glucoses. Because of their particular geometry, they are able to act as a “host” to form inclusion complexes with “guest” molecules very much like enzymes. Cyclodextrins have been shown to exert catalytic activity on suitable included-substrate molecules; they catalyze the hydrolysis of phenylacetates, of organic pyrophosphates and of penicillin derivatives. They also accelerate aromatic chlorinations and diazo coupling by means of their primary and/or secondary hydroxyl groups, so that the rates of hydrolysis are enhanced by up to a factor of 400. In order to understand the hydrogen bonding in these enzyme models, neutron diffraction data were collected to unambiguously determine the hydrogen atom positions, which could not be done from the x-ray diffraction data. α-CD has been shown to have two different structures with well-defined hydrogen bonds, one “tense” and the other “relaxed”. An “induced-fit”-like mechanism for α-CD complex formation has been proposed. Circular hydrogen bond networks have also been found for α-CD due to the energetically favored cooperative effect. β-CD with a disordered water structure possesses an unusual flip-flop hydrogen bonding system of the type O-H…H-O representing an equilibrium between two states: O-H…O?O…O. γ-CD with a disordered water structure similar to β-CD also possesses the flip-flop hydrogen bond. This study demonstrates that hydrogen bonds are operative in disordered systems and display dynamics even in the solid state.     

The hydration of the neurotransmitter acetylcholine in aqueous solution

Neutron diffraction augmented with hydrogen isotope substitution has been used to examine the water structure around the acetylcholine molecular ion in aqueous solution. It is shown that the nearest-neighbor water molecules in the region around the trimethylammonium headgroup are located either in a ring around the central nitrogen atom or between the carbon atoms, forming a sheath around the onium group. Moreover the water molecules in this cavity do not bond to the onium group but rather form hydrogen bonds with water molecules in the surrounding aqueous environment. Given that in the bound state the onium headgroup must be completely desolvated, the absence of bonding between the onium headgroup and the surrounding water solvent may be selectively favorable to acetylcholine-binding in the receptor site. Away from the headgroup, pronounced hydrogen-bonding of water to the carbonyl oxygen is observed, but not to the ether oxygen in the acetylcholine chain.     

Calculations on crystal packing of a flexible molecule, Leu-enkephalin

Crystal packing calculations have been carried out on a substantial number of conformations of Leu-enkephalin; namely, those obtained both from crystal structures and from energy minimizations on isolated molecules, and with and without waters of crystallization. The known crystal structures represent the most energetically stable packings found. The conformations of the enkephalin molecules in the crystal are not the most stable for an isolated molecule; i.e. intermolecular interactions force the isolated molecule to change conformation in order to achieve a small packing volume and an optimal packing energy in the crystal. It is found that the packing energy of an enkephalin molecule is a reasonably smooth function of its molecular volume in the unit cell, if structures with intermolecular hydrogen bonding are excluded, and is substantially independent of other details of the molecular conformation or of the crystal packing. Hydrogen bonding provides additional stabilization of the crystal structure, and would likely permit crystallization of the system if it is sufficiently dense. Solvent molecules further stabilize the structure when they can also provide intermolecular hydrogen bonds.     

The role of hydrogen bonding via interfacial water molecules in antigen-antibody complexation. The HyHEL-10-HEL interaction

To study the role of hydrogen bonding via interfacial water molecules in protein-protein interactions, we examined the interaction between hen egg white lysozyme (HEL) and its HyHEL-10 variable domain fragment (Fv) antibody. We constructed three antibody mutants (l-Y50F, l-S91A, and l-S93A) and investigated the interactions between the mutant Fvs and HEL. Isothermal titration calorimetry indicated that the mutations significantly decreased the negative enthalpy change (8-25 kJ mol(-1)), despite some offset by a favorable entropy change. X-ray crystallography demonstrated that the complexes had nearly identical structures, including the positions of the interfacial water molecules. Taken together, the isothermal titration calorimetric and x-ray crystallographic results indicate that hydrogen bonding via interfacial water enthalpically contributes to the Fv-HEL interaction despite the partial offset because of entropy loss, suggesting that hydrogen bonding stiffens the antigen-antibody complex.     

The hydrogen bonding in the crystal structure of raffinose pentahydrate

The crystal structure of raffinose pentahydrate, O-alpha-D-galactopyranosyl-(1----6)-O-alpha-D-glucopyranosyl-(1----2)- beta-D- fructofuranose pentahydrate, C18H32O16.5H2O, has been redetermined using low-temperature, 119 K, CuK alpha X-ray data. All hydrogen atoms were unambiguously located on difference syntheses. The final R-factor is 0.036 for 2423 observed structure amplitudes. The hydrogen bonding is composed of infinite chains, which are linked through the water molecules to form a three-dimensional network containing a chain of five linked water molecules. Three of the infinite chains extend in the directions of the crystallographic axis of the space group P2(1)2(1)2(1). Four of the water molecules accept two hydrogen bonds and one accepts one. All the hydroxyls and the ring and glycosidic oxygen atoms are involved in the hydrogen bonding. With one exception, the ring and glycosidic oxygens are hydrogen-bonded by means of the minor components of unsymmetrical three-center bonds.     

Cyclic peptide-metal salt adducts. I. crystal structure of the hexaquocopper(II) perchlorate cyclosarcosylsarcosine 1:2 adducts

The crystal structure analysis of the 2:1 adduct of cyclosarcosylsarcosine with copper(II) perchlorate shows that the independent unit is composed of six water molecules octahedrally coordinated to the Cu(II) ion, two tetrahedral perchlorate ions and four independent halves of cyclosarcosylsarcosine molecules lying on crystallographic centers of symmetry. All available hydrogens of water molecules are involved in hydrogen bonding as donors and all carbonyl oxygen atoms of the cyclic peptide molecules function as acceptors. Two other oxygen atoms for each perchlorate anion participate in the hydrogen bonding scheme, which leads to the absence of the orientational disorder usually observed in these ions. Layers of inorganic material and layers of organic material, roughly parallel to the ab plane, pack alternatively with each other. Electrostatic and ion-dipole interactions together with hydrogen bonds are responsible for the building up of the crystals.