Computer module design of the parametric structure of water

The algorithms of module design and the results of constructing the parametric structures of water were considered. The tetrahedral electron structure of the oxygen atom of water, which is a double symmetric donor and acceptor of protons, was taken as the main postulate. As opposed to the crystal lattices of diamond and ice, the hexacycle twist-bath was considered as the basic element of the original structure. A great variety of possible structures were obtained, which involve quasi-unidimensional (helices, rods), dendrite-like, ring-shaped, planar, and fractal structures, as well as combinatoric structures composed of these forms. The functions of distribution of valence angles in fractal structures optimized with respect to energy have two main maxima at 104 and 112 degrees, as differentiated from the starting ideal tetrahedral angle of 109.5 degrees. Changes in, and the accumulation of elastic energy of distortions of rod structure of different symmetry were analyzed. This energy can be utilized during the conformational changes of biopolymers. The substantial difference in the solubility of two anomeric forms of glucose is explained by a different degree of conformity of the glucose molecule structure to the structure of bound water.     

Phase diagram of deep rough mutant lipopolysaccharide from Salmonella minnesota R595.

The structural polymorphism of deep rough mutant lipopolysaccharide--in many biological systems the most active endotoxin--from Salmonella minnesota strain R595 was investigated as function of temperature, water content, and Mg2+ concentration. Differential scanning calorimetry was used to determine the amount of bound water and the enthalpy change at the beta<==>alpha gel to liquid crystalline acyl chain melting. The onset, midtemperature Tc, and completion of the beta<==>alpha phase transition were studied with Fourier-transform infrared spectroscopy. Synchrotron radiation X-ray diffraction was used to characterize the supramolecular three-dimensional structures in each phase state. The results indicate an extremely complex dependence of the structural behavior of LPS on ambient conditions. The beta<==>alpha acyl chain melting temperature Tc lying at 30 degrees C at high water content (95%) increases with decreasing water content reaching a value of 50 degrees C at 30% water content. Concomitantly, a broadening of the transition range takes place. At still lower water content, no distinct phase transition can be observed. This behavior is even more clearly expressed in the presence of Mg2+. In the lower water concentration range (< 50%) at temperatures below 70 degrees C, only lamellar structures can be observed independent of the Mg2+ concentration. This correlates with the absence of free water. Above 50% water concentration, the supramolecular structure below 70 degrees C strongly depends on the [LPS]:[Mg2+] ratio. For large [LPS]:[Mg2+] ratios, the predominant structure is nonlamellar, for smaller [LPS]:[Mg2+] ratios there is a superposition of lamellar and nonlamellar structures. At an equimolar ratio of LPS and Mg2+ a multibilayered organization is observed. The nonlamellar structures can be assigned in various cases to structures with cubic symmetry with periodicities between 12 and 16 nm. Under all investigated conditions, a transition into the hexagonal II structure takes place between 70 and 80 degrees C. These observations are discussed in relation to the biological importance of LPS as constituent of the outer membrane of gram-negative bacteria and as potent inducer of biological effects in mammals.     

Protonation status of metal-bound ligands can be determined by quantum refinement

The protonation status of key residues and bound ligands are often important for the function of a protein. Unfortunately, protons are not discerned in normal protein crystal structures, so their positions have to be determined by more indirect methods. We show that the recently developed quantum refinement method can be used to determine the position of protons in crystal structures. By replacing the molecular-mechanics potential, normally used in crystallographic refinement, by more accurate quantum chemical calculations, we get information about the ideal structure of a certain protonation state. By comparing the refined structures of different protonation states, the one that fits the crystallographic raw data best can be decided using four criteria: the R factors, electron density maps, strain energy, and divergence from the unrestrained quantum chemical structure. We test this method on alcohol dehydrogenase, for which the pK(a) of the zinc-bound solvent molecule is experimentally known. We show that we can predict the correct protonation state for both a deprotonated alcohol and a neutral water molecule.     

On the contribution of water-mediated interactions to protein-complex stability

Protein-water interactions have long been recognized as a major determinant of chain folding, conformational stability, binding specificity and catalysis. However, the detailed effects of water on stabilizing protein-protein interactions remain elusive. A way to test experimentally the contribution of water-mediated interactions is by applying double mutant cycle analysis on pairs of residues that do not form direct interactions, but are bridged by water. Seven such interactions within the interface between TEM1 and BLIP proteins were evaluated. No significant interaction free energy was found between either of them. Water can bridge interactions, but also stabilize the structure of the monomer. To distinguish between these, we performed a bioinformatic analysis using AQUAPROT (http://bioinfo.weizmann.ac.il/aquaprot) to determine the degree of water conservation between the bound and unbound states. 29 structures of twelve complexes and 20 related monomers were analyzed. Of the 262 water molecules located within the interfaces, 145 were conserved between the unbound and bound structures. Strikingly, all 50 buried or partially buried waters in the monomer structures were conserved at the same location in the bound structures. Thus, buried waters have an important role in stabilizing the monomer fold rather than contributing to protein-protein binding, and are not replaced by residues from the incoming protein. Taking together the experimental and bioinformatics evidence suggests that exposed waters within the interface may be good sites for protein engineering, while buried or mostly buried waters should be left unchanged.     

How similar are enzyme active site geometries derived from quantum mechanical theozymes to crystal structures of enzyme‐inhibitor complexes? Implications for enzyme design

Quantum mechanical optimizations of theoretical enzymes (theozymes), which are predicted catalytic arrays of biological functionalities stabilizing a transition state, have been carried out for a set of nine diverse enzyme active sites. For each enzyme, the theozyme for the rate-determining transition state plus the catalytic groups modeled by side-chain mimics was optimized using B3LYP/6-31G(d) or, in one case, HF/3-21G(d) quantum mechanical calculations. To determine if the theozyme can reproduce the natural evolutionary catalytic geometry, the positions of optimized catalytic atoms, i.e., covalent, partial covalent, or stabilizing interactions with transition state atoms, are compared to the positions of the atoms in the X-ray crystal structure with a bound inhibitor. These structure comparisons are contrasted to computed substrate-active site structures surrounded by the same theozyme residues. The theozyme/transition structure is shown to predict geometries of active sites with an average RMSD of 0.64 A from the crystal structure, while the RMSD for the bound intermediate complexes are significantly higher at 1.42 A. The implications for computational enzyme design are discussed.     

Enzyme substrate and inhibitor interactions.

An enzyme is designed to bind most tightly to a substrate when it is in the transition state of the reaction which the enzyme catalyses. The consequent reduction of the activation energy of the reaction constitutes the catalytic mechanism. The energetic contributions of different features of the interaction can only be crudely assessed, but they are dominated by entropically driven effects. The binding site of trypsin orients the substrate so that the reacting groups are correctly placed for reaction to occur. Apart from two side chains which take part in chemical steps of the reaction, the enzyme behaves almost as a rigid body. The full binding interactions are only developed when the substrate is in an intermediate stage of the reaction. The tightly bound complexes of trypsin with protein trypsin inhibitors have proved amenable to structural analysis. Enzyme inhibitor interactions, which account for almost 80 kJ mol-1 of interaction energy, are known fairly accurately. The similarity of the two known trypsin inhibitor structures, close to the primary binding site, indicates a high specificity, even for this simple interaction. In cases where no large conformational changes occur the specificity of an enzyme should be predictable from accurate knowledge of its tertiary structure.     

Protein-inhibitor flexible docking by a multicanonical sampling: native complex structure with the lowest free energy and a free-energy barrier distinguishing the native complex from the others

Flexible docking between a protein (lysozyme) and an inhibitor (tri-N-acetyl-D-glucosamine, tri-NAG) was carried out by an enhanced conformational sampling method, multicanonical molecular dynamics simulation. We used a flexible all-atom model to express lysozyme, tri-NAG, and water molecules surrounding the two bio-molecules. The advantages of this sampling method are as follows: the conformation of system is widely sampled without trapping at energy minima, a thermally equilibrated conformational ensemble at an arbitrary temperature can be reconstructed from the simulation trajectory, and the thermodynamic weight can be assigned to each sampled conformation. During the simulation, exchanges between the binding and free (i.e., unbinding) states of the protein and the inhibitor were repeatedly observed. The conformational ensemble reconstructed at 300 K involved various conformational clusters. The main outcome of the current study is that the most populated conformational cluster (i.e., the cluster of the lowest free energy) was assigned to the native complex structure (i.e., the X-ray complex structure). The simulation also produced non-native complex structures, where the protein and the inhibitor bound with different modes from that of the native complex structure, as well as the unbinding structures. A free-energy barrier (i.e., activation free energy) was clearly detected between the native complex structures and the other structures. The thermal fluctuations of tri-NAG in the lowest free-energy complex correlated well with the X-ray B-factors of tri-NAG in the X-ray complex structure. The existence of the free-energy barrier ensures that the lowest free-energy structure can be discriminated naturally from the other structures. In other words, the multicanonical molecular dynamics simulation can predict the native complex structure without any empirical objective function. The current study also manifested that the flexible all-atom model and the physico-chemically defined atomic-level force field can reproduce the native complex structure. A drawback of the current method is that it requires a time consuming computation due to the exhaustive conformational sampling. We discussed a possibility for combining the current method with conventional docking methods.     

Stereochemistry of binding of the tetrapeptide acetyl-Pro-Ala-Pro-Tyr-NH2 to porcine pancreatic elastase. Combined use of two-dimensional transferred nuclear Overhauser enhancement measurements, restrained molecular dynamics, X-ray crystallography and molecular modelling

A nuclear magnetic resonance study of the conformation of the tetrapeptide acetyl-Pro-Ala-Pro-Tyr-NH2 bound to porcine pancreatic elastase is presented. From two-dimensional transferred nuclear Overhauser enhancement measurements, a set of 23 approximate distance restraints between pairs of bound ligand protons, indicative of an extended type structure, is derived. The structure of the bound tetrapeptide is then refined from two different starting structures (an extended beta-strand and a polyproline helix) by restrained molecular dynamics, in which the interproton distances are incorporated into the total energy of the system in the form of effective potentials. Convergence to essentially the same average restrained dynamics structures is achieved. The refined structures are then modelled into the active site of elastase by interactive molecular graphics. The determination of the anchor point of the bound tetrapeptide on the enzyme was aided by a simultaneous crystallographic study which, despite the fact that only electron density for a Pro-X dipeptide fragment was visible, enabled both the approximate position and orientation of binding to be determined. It is found that the tetrapeptide is bound in the S' binding site in the reverse orientation found in other serine protease-inhibitor complexes and is stabilized both by hydrogen-bonding and by van der Waals' interactions.     

Water structures of differing order and mobility in aqueous solutions of schizophyllan, a triple-helical polysaccharide as revealed by dielectric dispersion measurements

Dielectric dispersion measurements were made on aqueous solutions of a triple-helical polysaccharide schizophyllan over a wide concentration range 10-50 wt % at -45 to +30 degrees C. In the solution state, three different water structures with the different relaxation times tau were found, namely, bound water (taul), structured water (taus), and loosely structured water (tauls) in addition to free water (tauP). Structured water is less mobile and loosely structured water is nearly as mobile as free water, but bound water with taul is much less mobile, thus taul > taus > tauls greater, similar tauP. The order-disorder transition accompanies the conversion between structured water and loosely structured water. However, the species with taus remains even in the disordered state and constitutes part of bound water in the entire temperature range. In the frozen state, in addition to bulk water formed by partial melting, two mobile species existed, which were assigned to liquidlike bound water and found to be a continuation of bound water in the solution state. These relaxation time data are discussed in connection with the entropy levels of the four structures deduced from heat capacity data (cf. Yoshiba, K.; et al. Biomacromolecules 2003, 4, 1348-1356).     

Structural studies show energy transfer within stabilized phycobilisomes independent of the mode of rod–core assembly

The major light harvesting complex in cyanobacteria and red algae is the phycobilisome (PBS), comprised of hundreds of seemingly similar chromophores, which are protein bound and assembled in a fashion that enables highly efficient uni-directional energy transfer to reaction centers. The PBS is comprised of a core containing 2–5 cylinders surrounded by 6–8 rods, and a number of models have been proposed describing the PBS structure. One of the most critical steps in the functionality of the PBS is energy transfer from the rod substructures to the core substructure. In this study we compare the structural and functional characteristics of high-phosphate stabilized PBS (the standard fashion of stabilization of isolated complexes) with cross-linked PBS in low ionic strength buffer from two cyanobacterial species, Thermosynechococcus vulcanus and Acaryochloris marina. We show that chemical cross-linking preserves efficient energy transfer from the phycocyanin containing rods to the allophycocyanin containing cores with fluorescent emission from the terminal emitters. However, this energy transfer is shown to exist in PBS complexes of different structures as characterized by determination of a 2.4 Å structure by X-ray crystallography, single crystal confocal microscopy, mass spectrometry and transmission electron microscopy of negatively stained and cryogenically preserved complexes. We conclude that the PBS has intrinsic structural properties that enable efficient energy transfer from rod substructures to the core substructures without requiring a single unique structure. We discuss the significance of our observations on the functionality of the PBS in vivo.     

High-resolution structures of adenylate kinase from yeast ligated with inhibitor Ap5A, showing the pathway of phosphoryl transfer.

The structure of adenylate kinase from yeast ligated with the two-substrate-mimicking inhibitor Ap5A and Mg2+ has been refined to 1.96 A resolution. In addition, the refined structure of the same complex with a bound imidazole molecule replacing Mg2+ has been determined at 1.63 A. These structures indicate that replacing Mg2+ by imidazole disturbs the water structure and thus the complex. A comparison with the G-proteins shows that Mg2+ is exactly at the same position with respect to the phosphates. However, although the Mg2+ ligand sphere of the G-proteins is a regular octahedron containing peptide ligands, the reported adenylate kinase has no such ligands and an open octahedron leaving space for the Mg2+ to accompany the transferred phosphoryl group. A superposition of the known crystalline and therefore perturbed phosphoryl transfer geometries in the adenylate kinases demonstrates that all of them are close to the start of the forward reaction with bound ATP and AMP. Averaging all observed perturbed structures gives rise to a close approximation of the transition state, indicating in general how to establish an elusive transition state geometry. The average shows that the in-line phosphoryl transfer is associative, because there is no space for a dissociative metaphosphate intermediate. As a side result, the secondary dipole interaction in the alpha-helices of both protein structures has been quantified.     

Modulation of an IDP binding mechanism and rates by helix propensity and non-native interactions: association of HIF1α with CBP

Intrinsically disordered proteins that acquire their three dimensional structures only upon binding to their targets are very important in cellular signal regulation. While experimental studies have been made on the structures of both bound (structured) and unbound (disordered) states, less is known about the actual folding-binding transition. Coarse grained simulations using native-centric (i.e. Gō) potentials have been particularly useful in addressing this problem, given the large search space for IDP binding, but have well-known deficiencies in reproducing the unfolded state structure and dynamics. Here, we investigate the interaction of HIF1α with CBP using a hierarchy of coarse-grained models, in each case matching the binding affinity at 300 K to the experimental value. Starting from a pure Gō-like model based on the native structure of the complex we go on to consider a more realistic model of helix propensity in the HIF1α, and finally the effect of non-native interactions between binding partners. We find structural disorder (i.e."fuzziness") in the bound state of HIF1α in all models which is supported by the results of atomistic simulations. Correcting the over-stabilized helices in the unbound state gives rise to a more cooperative folding-binding transition (destabilizing partially bound intermediates). Adding non-native contacts lowers the free energy barrier for binding to an almost barrierless scenario, leading to higher binding/unbinding rates relative to the other models, in better agreement with the near diffusion-limited binding rates measured experimentally. Transition state structures for the three models are highly disordered, supporting a fly-casting mechanism for binding.     

Ion-induced folding of the hammerhead ribozyme: a fluorescence resonance energy transfer study.

The ion-induced folding transitions of the hammerhead ribozyme have been analysed by fluorescence resonance energy transfer. The hammerhead ribozyme may be regarded as a special example of a three-way RNA junction, the global structure of which has been studied by comparing the distances (as energy transfer efficiencies) between the ends of pairs of labelled arms for the three possible end-to-end vectors as a function of magnesium ion concentration. The data support two sequential ion-dependent transitions, which can be interpreted in the light of the crystal structures of the hammerhead ribozyme. The first transition corresponds to the formation of a coaxial stacking between helices II and III; the data can be fully explained by a model in which the transition is induced by a single magnesium ion which binds with an apparent association constant of 8000-10 000 M-1. The second structural transition corresponds to the formation of the catalytic domain of the ribozyme, induced by a single magnesium ion with an apparent association constant of approximately 1100 M-1. The hammerhead ribozyme provides a well-defined example of ion-dependent folding in RNA.     

Properties of Double-Stranded DNA as a Polyelectrolyte

The stability of the structure of double-stranded DNA in the salt-free solution is discussed on the basis of the polyelectrolyte theory. Assuming that DNA is an infinitely long rod, and the formation of double strands is divided into combining process and folding process, the free energy changes required in these processes are calculated by the use of the exact solutions of two-dimensional Poisson-Boltzmann equation for the one rod and the two rod systems.

By strong depression of electrostatic interaction due to counter-ion condensation phenomena, the free energy change is remarkably decreased so that the double-stranded structure of DNA can be stabilized by energy of hydrogen bonds between base pairs. The increase of the activity coefficient of a counterion upon heat denaturation of DNA is also explained.

     

Effect of sugars on the hydration of serum albumin at low temperatures

The influence of sugars (glucose, fructose, and sucrose) on the hydration characteristics of serum albumin was studied using 1H NMR spectroscopy in combination with layer-by-layer freezing-out of bulk and interfacial water. It was found that the presence of sugars in protein solution leads to a considerable decrease in the concentration of bound water at T < 273 K; i. e., sugars cause the dehydration of protein molecules, which may be caused by those alterations in albumin structure that are associated with the formation of more compact globular structures. The most considerable effect was recorded in case of sucrose, which causes a decrease in the dehydration of albumin by at least one order of magnitude. The interfacial energy values for the protein/water system were calculated.     

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.     

Factors affecting the structure of substituted tris(pyrazolyl)borate rhodium dicarbonyl complexes

Predictions by density functional calculations of the structure and relative energy of various isomers of the hydridotris(pyrazol-1-yl)borate ligand in Tp^3R,5R rhodium(I) dicarbonyl complexes (R=H, Me) and their IR and ¹¹B NMR spectra are compared to experimental observations. The lowest energy structure of Tp^3,5-Me-, Tp^3-Me-, and TpRh(CO)₂ is a non-classical square pyramidal (SPy) structure with a long metal apical ligand distance in rapid exchange with an equivalent SPy structure through a low energy trigonal bipyramid (TBP) transition state (a ‘reverse’ Berry pseudorotation). A second higher energy minimum, a pseudo square-planar complex with the third uncoordinated pyrazolyl arm rotated approximately parallel with the metal ligand pseudo-plane (SP1), is accessed through a second low energy transition state. Another pseudo square-planar minimum structure (SP2) is produced by a transition state, which lengthens the rhodium-apical nitrogen (of the third pyrazolyl arm) bond distance. The relative stability of SP2 depends on the degree of tris(pyrazolyl)borate (Tp) substitution, where 5-substituents larger than hydrogen disfavor SP2 because of steric interactions. The previously reported empirical correlation between ¹¹B NMR chemical shifts, ν_BH stretching frequencies, and the crystallographic Tp ligand denticity is reproduced by our calculations. The variety of structures observed by experiment can be explained by the calculated relative energies of the structures, the bulk dielectric of the solvents when in solution, specific interaction by certain solvents, and conditions of crystallization when in the solid-state.     

Binding of the Bacteriophage P22 N-Peptide to the boxB RNA Motif Studied by Molecular Dynamics Simulations

Protein-RNA interactions are important for many cellular processes. The Nut-utilization site (N)-protein of bacteriophages contains an N-terminal arginine-rich motif that undergoes a folding transition upon binding to the boxB RNA hairpin loop target structure. Molecular dynamics simulations were used to investigate the dynamics of the P22 N-peptide-boxB complex and to elucidate the energetic contributions to binding. In addition, the free-energy changes of RNA and peptide conformational adaptation to the bound forms, as well as the role of strongly bound water molecules at the peptide-RNA interface, were studied. The influence of peptide amino acid substitutions and the salt dependence of interaction were investigated and showed good agreement with available experimental results. Several tightly bound water molecules were found at the RNA-binding interface in both the presence and absence of N-peptide. Explicit consideration of the waters resulted in shifts of calculated contributions during the energetic analysis, but overall similar binding energy contributions were found. Of interest, it was found that the electrostatic field of the RNA has a favorable influence on the coil-to-α-helix transition of the N-peptide already outside of the peptide-binding site. This result may have important implications for understanding peptide-RNA complex formation, which often involves coupled folding and association processes. It indicates that electrostatic interactions near RNA molecules can lead to a shift in the equilibrium toward the bound form of an interacting partner before it enters the binding pocket.     

Structure of mandelate racemase with bound intermediate analogues benzohydroxamate and cupferron

Mandelate racemase (MR, EC 5.1.2.2) from Pseudomonas putida catalyzes the Mg(2+)-dependent interconversion of the enantiomers of mandelate, stabilizing the altered substrate in the transition state by 26 kcal/mol relative to the substrate in the ground state. To understand the origins of this binding discrimination, we determined the X-ray crystal structures of wild-type MR complexed with two analogues of the putative aci-carboxylate intermediate, benzohydroxamate and Cupferron, to 2.2-? resolution. Benzohydroxamate is shown to be a reasonable mimic of the transition state and/or intermediate because its binding affinity for 21 MR variants correlates well with changes in the free energy of transition state stabilization afforded by these variants. Both benzohydroxamate and Cupferron chelate the active site divalent metal ion and are bound in a conformation with the phenyl ring coplanar with the hydroxamate and diazeniumdiolate moieties, respectively. Structural overlays of MR complexed with benzohydroxamate, Cupferron, and the ground state analogue (S)-atrolactate reveal that the para carbon of the substrate phenyl ring moves by 0.8-1.2 ? between the ground state and intermediate state, consistent with the proposal that the phenyl ring moves during MR catalysis while the polar groups remain relatively fixed. Although the overall protein structure of MR with bound intermediate analogues is very similar to that of MR with bound (S)-atrolactate, the intermediate-Mg(2+) distance becomes shorter, suggesting a tighter complex with the catalytic Mg(2+). In addition, Tyr 54 moves closer to the phenyl ring of the bound intermediate analogues, contributing to an overall constriction of the active site cavity. However, site-directed mutagenesis experiments revealed that the role of Tyr 54 in MR catalysis is relatively minor, suggesting that alterations in enzyme structure that contribute to discrimination between the altered substrate in the transition state and the ground state by this proficient enzyme are extremely subtle.