Exploring the roles of nucleobase desolvation and shape complementarity during the misreplication of O(6)-methylguanine

O⁶-methylguanine (O⁶-MeG) is a miscoding DNA lesion arising from the alkylation of guanine. This report uses the bacteriophage T4 DNA polymerase as a model to probe the roles of hydrogen-bonding interactions, shape/size, and nucleobase desolvation during the replication of this miscoding lesion. This was accomplished by using transient kinetic techniques to monitor the kinetic parameters for incorporating and extending natural and nonnatural nucleotides. In general, the efficiency of nucleotide incorporation does not depend on the hydrogen-bonding potential of the incoming nucleotide. Instead, nucleobase hydrophobicity and shape complementarity appear to be the preeminent factors controlling nucleotide incorporation. In addition, shape complementarity plays a large role in controlling the extension of various mispairs containing O⁶-MeG. This is evident as the rate constants for extension correlate with proper interglycosyl distances and symmetry between the base angles of the formed mispair. Base pairs not conforming to an acceptable geometry within the polymerase's active site are refractory to elongation and are processed via exonuclease proofreading. The collective data set encompassing nucleotide incorporation, extension, and excision is used to generate a model accounting for the mutagenic potential of O⁶-MeG observed in vivo. In addition, kinetic studies monitoring the incorporation and extension of nonnatural nucleotides identified an analog that displays high selectivity for incorporation opposite O⁶-MeG compared to unmodified purines. The unusual selectivity of this analog for replicating damaged DNA provides a novel biochemical tool to study translesion DNA synthesis.     

Cooperative hydrogen bonding in amyloid formation

Amyloid diseases, including Alzheimer's and prion diseases, are each associated with unbranched protein fibrils. Each fibril is made of a particular protein, yet they share common properties. One such property is nucleation-dependent fibril growth. Monomers of amyloid-forming proteins can remain in dissolved form for long periods, before rapidly assembly into fibrils. The lag before growth has been attributed to slow kinetics of formation of a nucleus, on which other molecules can deposit to form the fibril. We have explored the energetics of fibril formation, based on the known molecular structure of a fibril-forming peptide from the yeast prion, Sup35, using both classical and quantum (density functional theory) methods. We find that the energetics of fibril formation for the first three layers are cooperative using both methods. This cooperativity is consistent with the observation that formation of amyloid fibrils involves slow nucleation and faster growth.     

Testing electrostatic complementarity in enzyme catalysis: hydrogen bonding in the ketosteroid isomerase oxyanion hole

A longstanding proposal in enzymology is that enzymes are electrostatically and geometrically complementary to the transition states of the reactions they catalyze and that this complementarity contributes to catalysis. Experimental evaluation of this contribution, however, has been difficult. We have systematically dissected the potential contribution to catalysis from electrostatic complementarity in ketosteroid isomerase. Phenolates, analogs of the transition state and reaction intermediate, bind and accept two hydrogen bonds in an active site oxyanion hole. The binding of substituted phenolates of constant molecular shape but increasing p K _a models the charge accumulation in the oxyanion hole during the enzymatic reaction. As charge localization increases, the NMR chemical shifts of protons involved in oxyanion hole hydrogen bonds increase by 0.50–0.76 ppm/p K _a unit, suggesting a bond shortening of ˜0.02 Å/p K _a unit. Nevertheless, there is little change in binding affinity across a series of substituted phenolates (ΔΔG = −0.2 kcal/mol/p K _a unit). The small effect of increased charge localization on affinity occurs despite the shortening of the hydrogen bonds and a large favorable change in binding enthalpy (ΔΔH = −2.0 kcal/mol/p K _a unit). This shallow dependence of binding affinity suggests that electrostatic complementarity in the oxyanion hole makes at most a modest contribution to catalysis of ˜300-fold. We propose that geometrical complementarity between the oxyanion hole hydrogen-bond donors and the transition state oxyanion provides a significant catalytic contribution, and suggest that KSI, like other enzymes, achieves its catalytic prowess through a combination of modest contributions from several mechanisms rather than from a single dominant contribution.     

Association of star-shaped poly(D,L-lactide)s containing nucleobase multiple hydrogen bonding

A new family of associating polymers based on four-arm, star-shaped poly(D,L-lactide) (PDLLA) containing peripheral complementary hydrogen-bonding sites is described. Hydroxy-terminated, four-arm, star-shaped PDLLAs of controlled molar mass were functionalized with complementary DNA base pairs, adenine (A) and thymine (T), to obtain PDLLA-A and PDLLA-T, respectively. 1H NMR spectroscopy confirmed quantitative functionalization and the subsequent formation of PDLLA-A and PDLLA-T hydrogen-bonded complexes. Job's analysis revealed a 1:1 optimal stoichiometry for the hydrogen-bonded complexes, and the association constant (Ka) that was determined using the 1H NMR-based Benesi-Hildebrand treatment was 84 M(-1) for the low molar mass complementary polymers. Furthermore, the PDLLA-based hydrogen-bonded complexes exhibited higher solution viscosities compared to the corresponding non-hydrogen-bonded precursors, which further confirmed strong complementary multiple-hydrogen-bonding associations between the star-shaped polymers with terminal adenine and thymine groups. Moreover, variable-temperature 1H NMR studies demonstrated the thermoreversibility of the hydrogen-bonded PDLLA-based complexes in solution.     

A fast protein-ligand docking algorithm based on hydrogen bond matching and surface shape complementarity

With the rapid development of structural determination of target proteins for human diseases, high throughout virtual screening based drug discovery is gaining popularity gradually. In this paper, a fast docking algorithm (H-DOCK) based on hydrogen bond matching and surface shape complementarity was developed. In H-DOCK, firstly a divide-and-conquer strategy based enumeration approach is applied to rank the intermolecular modes between protein and ligand by maximizing their hydrogen bonds matching, then each docked conformation of the ligand is calculated according to the matched hydrogen bonding geometry, finally a simple but effective scoring function reflecting mainly the van der Waals interaction is used to evaluate the docked conformations of the ligand. H-DOCK is tested for rigid ligand docking and flexible one, the latter is implemented by repeating rigid docking for multiple conformations of a small molecule and ranking all together. For rigid ligands, H-DOCK was tested on a set of 271 complexes where there is at least one intermolecular hydrogen bond, and H-DOCK achieved success rate (RMSD<2.0?Å) of 91.1%. For flexible ligands, H-DOCK was tested on another set of 93 complexes, where each case was a conformation ensemble containing native ligand conformation as well as 100 decoy ones generated by AutoDock [1], and the success rate reached 81.7%. The high success rate of H-DOCK indicates that the hydrogen bonding and steric hindrance can grasp the key interaction between protein and ligand. H-DOCK is quite efficient compared with the conventional docking algorithms, and it takes only about 0.14 seconds for a rigid ligand docking and about 8.25 seconds for a flexible one on average. According to the preliminary docking results, it implies that H-DOCK can be potentially used for large scale virtual screening as a pre-filter for a more accurate but less efficient docking algorithm.     

Combination of covalent and hydrogen bonding in the formation of 3D uranyl-carboxylate networks

Four compounds containing uranyl cation [UO₂]²⁺ have been synthesized hydrothermally by reacting uranyl acetate and uranyl nitrate with various N/O donor ligands. The structure of all compounds was elucidated by single crystal X-ray diffraction study. Compound [(UO₂)(6-methylnicotinato)₃](H₃O)·4H₂O (1) is a discrete complex (0D), that gives rise in the crystal to hydrophilic channels, while [(UO₂)(OH)(μ₂-3-pyridylpropionato)]_n (2) and [(UO₂)(H₂O)(μ₃-4,4′-oxybis(benzoato)]_n (3) show the formation of 1D coordination polymers. Moreover, oxalate anions, formed in situ by using 5-methylisophthalic or 2,3-pyrazinedicarboxylic acid as reactant ligands, gave rise to a 2D coordinating network [(UO₂)₂(μ₂-oxalate)(μ₂-OH)₂(H₂O)₂]_n·nH₂O (4). All the complexes expanded their dimensionality to 3D through hydrogen bonds.     

Inhibition of the initiation of leukemic transcription by N-trifluoroacetyladriamycin-14-O-hemiadipate in vitro. Impaired formation of RNA polymerase-DNA complex

The effect of N-trifluoroacetyladriamycin-14-O-hemiadipate (AD 143), a new derivative of adriamycin, on various steps of the enzymic reaction catalyzed by chicken myeloblastosis RNA polymerase II was studied. AD 143 inhibition of RNA synthesis, which was evident at the beginning of the reaction, could not be reversed by increasing the concentrations of any one of the four nucleoside triphosphate substrates of the reaction. Furthermore, the RNA synthesis inhibition was not affected by varying the concentrations of template DNA. The AD 143-induced inhibition caused a reduction of the frequency of RNA chain initiation, whereas the average chain length of RNA synthesized at the end of the reaction remained unaltered. The susceptible step in the initiation process was found to be the formation of stable complexes between RNA polymerase and the DNA template. While AD 143 causes no inhibition of Escherichia coli RNA polymerase activity, it was found not to affect the E. coli RNA polymerase-template DNA complex formation.     

Effect of solute hydrogen bonding capacity on osmotic stability of lysosomes

The effect of solute hydrogen bonding capacity on the osmotic stability of lysosomes was examined through measurement of free enzyme activity of lysosomes after their incubation in sucrose and poly(ethylene glycol) (PEG) (1500–6000 Da molecular mass) media. Free enzyme activity of the lysosomes was less in the PEG medium than that in the sucrose medium under the same hypotonic condition. The lysosomal enzyme latency loss decreased with increasing hydrogen bonding capacity of the solute. In addition, the lysosomes lost less latency at lower incubation temperature. The results indicate that solute hydrogen bonding capacity plays an important role in the osmotic protection of an incubation medium to lysosomes.     

Contribution of hydrogen bonding to protein stability estimated from isotope effects

An unresolved issue in structural biology concerns the relative contribution of H bonds to protein stability. We use the small molecules 4-acetamidobenzoic acid and N-acetylanthranilic acid as model compounds to relate the energetic contribution from hydrogen bonds (H bonds) to the deuterium/hydrogen amide isotope effect. N-Acetylanthranilic acid models carbonyl-amide H bonds formed during protein folding; 4-acetamidobenzoic acid models the unfolded state in which the amide H bonds to water. NMR is used to measure shifts in the pK(a) of the ionizable carboxyl group when the amides of the compounds are either protonated or deuterated. From the pK(a) shift, we obtain a quantitative scale factor: SF = partial partial differential(DeltaG(HB))/partial partial differential(RT ln Phi), where DeltaG(HB) is the change in free energy of an H bond upon isotope substitution and Phi is the fractionation factor. Isotope effect data also are reported for a small globular protein, lambda repressor, using the "C(m) experiment". The protein's isotope effect, which reports on the shape of the energy well, is converted to H-bonding free energy by applying the scale factor. We estimate that amide-related H bonds (amide-carbonyl and amide-water) contribute favorably to protein stability by approximately 30-50 kcal/mol in lambda repressor, GCN4 coiled coil, and cytochrome c but unfavorably by approximately 6 kcal/mol in ubiquitin. The results indicate that H-bond strength varies from one protein to another and presumably at different sites within the same protein.     

Contribution of hydrogen bonding to the conformational stability of ribonuclease T1.

For 30 years, the prevailing view has been that the hydrophobic effect contributes considerably more than hydrogen bonding to the conformational stability of globular proteins. The results and reasoning presented here suggest that hydrogen bonding and the hydrophobic effect make comparable contributions to the conformational stability of ribonuclease T1 (RNase T1). When RNase T1 folds, 86 intramolecular hydrogen bonds with an average length of 2.95 A are formed. Twelve mutants of RNase T1 [Tyr----Phe (5), Ser----Ala (3), and Asn----Ala (4)] have been prepared that remove 17 of the hydrogen bonds with an average length of 2.93 A. On the basis of urea and thermal unfolding studies of these mutants, the average decrease in conformational stability due to hydrogen bonding is 1.3 kcal/mol per hydrogen bond. This estimate is in good agreement with results from several related systems. Thus, we estimate that hydrogen bonding contributes about 110 kcal/mol to the conformational stability of RNase T1 and that this is comparable to the contribution of the hydrophobic effect. Accepting the idea that intramolecular hydrogen bonds contribute 1.3 +/- 0.6 kcal/mol to the stability of systems in an aqueous environment makes it easier to understand the stability of the "molten globule" states of proteins, and the alpha-helical conformations of small peptides.     

Kinetic identification of a hydrogen bonding pair in the glucoamylase-maltose transition state complex.

Molecular recognition and site-directed mutagenesis are used in combination to identify kinetically, transition state interactions between glucoamylase (GA) and the substrate maltose. Earlier studies of mutant Glu180----Gln GA had indicated a role in substrate binding for Glu180 (Sierks, M.R., Ford, C., Reilly, P.J. and Svensson, B. (1990) Protein Engng, 3, 193-198). Here, changes in activation energies calculated from measured kcat/Km values for a series of deoxygenated maltose analogues indicate hydrogen bonding between the mutant enzyme and the 3-OH group of the reducing end sugar ring. Using the same substrate analogues and determining activation energies with wild-type GA an additional hydrogen bond with the 2-OH group of maltose is attributed to an interaction with the carboxylate Glu180. This novel combination of molecular recognition and site-directed mutagenesis enables an enzyme substrate transition state contact to be identified and characterized even without access to the three dimensional structure of the enzyme. Given the distant structural relationships between glucoamylases and several starch hydrolases (Svensson,B. (1988) FEBS Lett., 230, 72-76), such identified contacts may ultimately guide tailoring of the activity of these related enzymes.     

The strength of dehalogenase-substrate hydrogen bonding correlates with the rate of Meisenheimer intermediate formation

4-Chlorobenzoyl-coenzyme A (4-CBA-CoA) dehalogenase catalyzes the hydrolytic dehalogenation of 4-CBA-CoA to 4-hydroxybenzoyl-CoA by using an active site aspartate as the nucleophile. Formation of the corresponding Meisenheimer complex (EMc) is followed by chloride ion expulsion which forms the arylated intermediate (EAr). This is then hydrolyzed to the product. In this paper, we explore the relationship between active site polarizing forces acting on the benzoyl carbonyl and the rate of formation of the Meisenheimer complex. The polarizing forces at the C[double bond]O group were modulated by introducing site-selected mutations (A112V, Y65D, G113A, G113S, G113N, and F64P), near the C[double bond]O binding site. Using either the substrate, 4-CBA-CoA, or the substrate analogue, 4-methylbenzoyl-CoA (4-MBA-CoA), Raman difference spectroscopy provided the position of the C[double bond]O stretching frequency (nu(C)[double bond](O)) for a total of 10 enzyme-ligand complexes. In turn, the values of the C[double bond]O frequencies could be converted to differences in effective hydrogen bonding strengths between members of the series, based on earlier model studies [Clarkson, J., Tonge, P. J., Taylor, K. L., Dunaway-Mariano, D., and Carey, P. (1997) Biochemistry 36, 10192-10199]. Catalysis in the F64P, G113A, G113S, and G113N dehalogenase mutants was very slow with k(cat) values ranging from 8 x 10(-3) to 7.6 x 10(-6) s(-1). The EAr intermediate did not accumulate to a detectable level on these enzymes during a single turnover. Catalysis in the Y65D and A112V dehalogenase mutants were almost as efficient as catalysis in wild-type dehalogenase with k(cat) values of 0.1-0.6 s(-1). In wild-type dehalogenase, 22% of the bound substrate accumulated as the EAr intermediate during a single turnover (k(obs) for EAr formation = 24 s(-(1)); in the Y65D mutant, the level of accumulation is 17% (k(obs) for EAr formation = 3 s(-1)), and in the A112V mutant, the level is 23% (k(obs) for EAr formation = 17 s(-1)). The k(obs) for EAr formation in wild-type dehalogenase and the more active dehalogenase mutants (Y65D and A112V) was taken to be an estimate of the k for EMc formation, and the k(obs) for EP formation in a single turnover was taken to be an estimate of the k for EMc formation in the severely impaired mutants (F64P, G113A, G113S, and G113N). A plot of the log k(obs) for EMc formation versus the C[double bond]O stretching frequency of bound 4-CBA-CoA (or 4-MBA-CoA) is a straight line (R(2) = 0.9584). Throughout the series, nu(C)[double bond](O) varied by 61 cm(-1), corresponding to the change in hydrogen bonding enthalpy of 67 kJ/mol. The results show that changes in polarizing forces at the benzoyl carbonyl are transmitted to the benzoyl (4) position and correlate with the rate of aromatic nucleophilic addition five chemical bonds away. Interestingly, the relationship between effective polarizing forces and reactivity seen here for dehalogenase is similar to that reported for the addition-elimination reaction involving the hydrolysis of a series of acyl serine proteases.     

Conformation, hydrogen bonding and aggregate formation of guanosine 5'-monophosphate and guanosine in dimethylsulfoxide.

The tetrabutylammonium salt of guanosine 5'-monophosphate (5'-GMP) dissolves in DMSO-d6 forming aggregated species which exhibit some properties of reverse micelles. 1H NOESY experiments show that the 5'-GMP adopts the syn conformation about the glycosidic bond. Molecular mechanics calculations reveal a stable structure with this conformation in which the phosphate group and the amino group of the base are in close enough proximity to hydrogen bond. In contrast inosine 5'-monophosphate in DMSO-d6, which has no NH2 group for hydrogen bond stabilization of the syn conformation, is shown by NMR to have the anti structure. Guanosine in DMSO-d6 behaves differently from 5'-GMP. Guanosine adopts the anti conformation and forms a symmetric dimer via hydrogen bonding between the N3 and NH2 of the bases.     

Solid state NMR studies of hydrogen bonding in a citrate synthase inhibitor complex.

The ionization state and hydrogen bonding environment of the transition state analogue (TSA) inhibitor, carboxymethyldethia coenzyme A (CMX), bound to citrate synthase have been investigated using solid state NMR. This enzyme-inhibitor complex has been studied in connection with the postulated contribution of short hydrogen bonds to binding energies and enzyme catalysis: the X-ray crystal structure of this complex revealed an unusually short hydrogen bond between the carboxylate group of the inhibitor and an aspartic acid side chain [Usher et al. (1994) Biochemistry 33, 7753-7759]. To further investigate the nature of this short hydrogen bond, low spinning speed 13C NMR spectra of the CMX-citrate synthase complex were obtained under a variety of sample conditions. Tensor values describing the chemical shift anisotropy of the carboxyl groups of the inhibitor were obtained by simulating MAS spectra (233 +/- 4, 206 +/- 5, and 105 +/- 2 ppm vs TMS). Comparison of these values with our previously reported database and ab initio calculations of carbon shift tensor values clearly indicates that the carboxyl is deprotonated. New data from model compounds suggest that hydrogen bonds in a syn arrangement with respect to the carboxylate group have a pronounced effect upon the shift tensors for the carboxylate, while anti hydrogen bonds, regardless of their length, apparently do not perturb the shift tensors of the carboxyl group. Thus the tensor values for the enzyme-inhibitor complex could be consistent with either a very long syn hydrogen bond or an anti hydrogen bond; the latter would agree very well with previous crystallographic results. Two-dimensional 1H-13C heteronuclear correlation spectra of the enzyme-inhibitor complex were obtained. Strong cross-peaks were observed from the carboxyl carbon to proton(s) with chemical shift(s) of 22 +/- 5 ppm. Both the proton chemical shift and the intensity of the cross-peak indicate a very short hydrogen bond to the carboxyl group of the inhibitor, the C.H distance based upon the cross-peak intensity being 2.0 +/- 0.4 A. This proton resonance is assigned to Hdelta2 of Asp 375, on the basis of comparison with crystal structures and the fact that this cross-peak was absent in the heteronuclear correlation spectrum of the inhibitor-D375G mutant enzyme complex. In summary, our NMR studies support the suggestion that a very short hydrogen bond is formed between the TSA and the Asp carboxylate.     

Hydrogen bond formation between the naturally modified nucleobase and phosphate backbone

Natural RNAs, especially tRNAs, are extensively modified to tailor structure and function diversities. Uracil is the most modified nucleobase among all natural nucleobases. Interestingly, >76% of uracil modifications are located on its 5-position. We have investigated the natural 5-methoxy (5-O-CH(3)) modification of uracil in the context of A-form oligonucleotide duplex. Our X-ray crystal structure indicates first a H-bond formation between the uracil 5-O-CH(3) and its 5'-phosphate. This novel H-bond is not observed when the oxygen of 5-O-CH(3) is replaced with a larger atom (selenium or sulfur). The 5-O-CH(3) modification does not cause significant structure and stability alterations. Moreover, our computational study is consistent with the experimental observation. The investigation on the uracil 5-position demonstrates the importance of this RNA modification at the atomic level. Our finding suggests a general interaction between the nucleobase and backbone and reveals a plausible function of the tRNA 5-O-CH(3) modification, which might potentially rigidify the local conformation and facilitates translation.     

Examination of shape complementarity in docking of unbound proteins.

Here we carry out an examination of shape complementarity as a criterion in protein-protein docking and binding. Specifically, we examine the quality of shape complementarity as a critical determinant not only in the docking of 26 protein-protein "bound" complexed cases, but in particular, of 19 "unbound" protein-protein cases, where the structures have been determined separately. In all cases, entire molecular surfaces are utilized in the docking, with no consideration of the location of the active site, or of particular residues/atoms in either the receptor or the ligand that participate in the binding. To evaluate the goodness of the strictly geometry-based shape complementarity in the docking process as compared to the main favorable and unfavorable energy components, we study systematically a potential correlation between each of these components and the root mean square deviation (RMSD) of the "unbound" protein-protein cases. Specifically, we examine the non-polar buried surface area, polar buried surface area, buried surface area relating to groups bearing unsatisfied buried charges, and the number of hydrogen bonds in all docked protein-protein interfaces. For these cases, where the two proteins have been crystallized separately, and where entire molecular surfaces are considered without a predefinition of the binding site, no correlation is observed. None of these parameters appears to consistently improve on shape complementarity in the docking of unbound molecules. These findings argue that simplicity in the docking process, utilizing geometrical shape criteria may capture many of the essential features in protein-protein docking. In particular, they further reinforce the long held notion of the importance of molecular surface shape complementarity in the binding, and hence in docking. This is particularly interesting in light of the fact that the structures of the docked pairs have been determined separately, allowing side chains on the surface of the proteins to move relatively freely. This study has been enabled by our efficient, computer vision-based docking algorithms. The fast CPU matching times, on the order of minutes on a PC, allow such large-scale docking experiments of large molecules, which may not be feasible by other techniques. Proteins 1999;36:307-317.     

Significance of hydrogen bonding at the S(1)' subsite of calpain I.

alpha-Ketohydroxamates were synthesized as bioisosteres of alpha-ketoamides. The alpha-ketohydroxamates were generally more potent than the corresponding alpha-ketoamides. The potency of the compounds suggests that hydrogen bonding and steric bulk of substituents on the nitrogen atom of the ketoamide moiety influence calpain inhibition.