The Influence of Starting Coordinates in Free Energy Simulations of Ligand Binding to Dihydrofolate Reductase

Abstract

Molecular dynamics (MD) simulation combined with free energy perturbation (FEP) methods have been used to study the key structural differences and relative free energies for the binding of 6-methyl-N5-deazapterin (N8 protonated) and the 8-substituted compound, 6,8-dimethyl-N5-deazapterin (N3 protonated), to dihydrofolate reductase (DHFR). The free energy changes have been calculated using a variety of initial X-ray coordinates derived from bacterial and vertebrate (including human) DHFRs, and both with and without the reduced cofactor nicotinamide adenine dinucleotide (NADPH) bound. Given a sufficiently long simulation time for the FEP calculations (ca. 200 ps), all structures obtained after mutating 6,8-methyl-N5-deazapterin to 6-methyl-N5-deazapterin exhibited hydrogen bond formation between a backbone carbonyl group of DHFR and H(N8) of 6-methyl-N5-deazapterin, analogous to that found in the X-ray crystal structure of N5-deazafolate(N8 protonated) bound to human DHFR. However, both simulation and experiment suggest this additional H-bonding does not greatly enhance thermodynamic stability, with experiment indicating at most a factor of 2 difference in the relative affinities of the two ligand cations for vertebrate DHFR. Moreover, a binding differential of 10 in favour of the protonated 8-substituted compound is found experimentally for bacterial DHFR. The MD/FEP calculations suggest that the relative cost of ligand desolvation may largely cancel the lowering of free energy obtained in the active site, resulting in predicted binding differences within the range indicated by the vertebrate and bacterial DHFR experiments. However, the theoretical free energy changes could not be obtained with the accuracy required for the rationalization of the observed species dependence. While sampling difficulties are known to be inherent in MD simulation methodologies, these studies with several initial coordinate sets have demonstrated the contribution of coordinate choice to this problem. The results indicate that for demanding protein-ligand binding problems such as this one, the accuracy of the method may be no better than ± 2 kcal/mol.     

Analysis of the binding of phenylalanine to phenylalanyl-tRNA synthetase

Using the complete rate equation for the PP_i-ATP exchange reaction at equilibrium, the dissociation constants of phenylalanine (10^?5m), phenylalanine butyl ester (8 × 10^?5m), benzyl alcohol (6 × 10^?4m), phenylalaninol (2 × 10^?4m), hydrocinnamic acid (3 × 10^?3m) and glycine (>1 m) with the phenylalanyl-tRNA synthetase (Escherichia coli K12) were determined. Taking the model of Koshland (1962) for the estimation of the configurational free energy change due to proximity and orientation, and decomposing the process of binding into several thermodynamic steps, the contribution to binding of the benzyl group, glycine unit, protonated amino group, carboxylate group and joint interactions were estimated. The results are: (1) the standard free energy contributions for binding phenylalanine are benzyl group (?8.2 kcal/mol), glycine unit (?2.5 kcal/mol), protonated amino group (?0.8 kcal/mol) and carboxylate group (1 kcal/mol). (2) The standard free energy change due to the change in the interaction between the protonated amino group and carboxylate group when they are transferred from the aqueous environment to the enzyme environment is ?2.7 kcal/mol. (3) A dissociation constant for glycine of 7.5 m is calculated without the hypothesis that a conformational change occurs in the enzyme when the benzyl unit of phenylalanine binds, permitting an interaction of the enzyme with the protonated amino and/or carboxylate groups.The detection of E·AA² and E·ATP shows that a sequential addition of substrates is not necessary for binding. A comparison of the dissociation constants of E·AA (10^?5m), E·ATP (1.5 × 10^?3m), E·PP (5.5 × 10^?4m), E·I (8 × 10^?5m) and the mixed complexes E·I·ATP (6 × 10^?8m²), E·I·PP (5 × 10^?8m²) and E·AA·PP (7 × 10^?9m²), with phenylalanine butyl ester as the inhibitor, indicates no strong interaction between the binding of ATP or PP_i with the binding of phenylalanine.     

Effects of pressure on enzyme function of Escherichia coli dihydrofolate reductase

To elucidate the effects of pressure on the function of Escherichia coli dihydrofolate reductase (DHFR), the enzyme activity and the dissociation constants of substrates and cofactors were measured at pressures up to 250 MPa at 25 degrees C and pH 7.0. The enzyme activity decreased with increasing pressure, accompanying the activation volume of 7.8 ml mol(-1). The values of the Michaelis constant (K(m)) for dihydrofolate and NADPH were slightly higher at 200 MPa than at atmospheric pressure. The hydride-transfer step was insensitive to pressure, as monitored by the effects of the deuterium isotope of NADPH on the reaction velocity. The dissociation constants of substrates and cofactors increased with pressure, producing volume reductions from 6.5 ml mol(-1) (tetrahydrofolate) to 33.5 ml mol(-1) (NADPH). However, the changes in Gibbs free energy with dissociation of many ligands showed different pressure dependences below and above 50 MPa, suggesting conformational changes of the enzyme at high pressure. The enzyme function at high pressure is discussed based on the volume levels of the intermediates and the candidates for the rate-limiting process.     

Reciprocal cooperative effects of multiple ligand binding to pyruvate kinase

The formation of multiple ligand complexes with muscle pyruvate kinase was measured in terms of dissociation constants and the standard free energies of formation were calculated. The binding of Mn2+ to the enzyme (KA = 55 +/- 5 X 10(-6) M; deltaF degrees = -5.75 +/- 0.05 kcal/mol) and to the enzyme saturated with phosphoenolpyruvate (conditional free energy) KA' = 0.8 +/- 0.4 X 10(-6) M; deltaF degrees = -8.22 +/- 0.34 kcal/mol) has been measured under identical conditions giving a free energy of coupling, delta(deltaF degrees) = -2.47 +/- 0.34 kcal/mol. Such a large negative free energy of coupling is diagnostic of a strong positively cooperative effect in ligand binding. The binding of the substrate phosphoenolpyruvate to free enzyme and the enzyme-Mn2+ complex was, by necessity, measured by different methods. The free energy of phosphoenolpyruvate binding to free enzyme (KS = 1.58 +/- 0.10 X 10(-4)M; deltaF degrees = -5.13 +/- 0.04 kcal/mol) and to the enzyme-Mn2+ complex (K3 = 0.75 +/- 0.10 X 10(-6)M; deltaF degrees = -8.26 +/- 0.07 kcal/mol) also gives a large negative free energy of coupling, delta(deltaF degrees) = -3.16 +/- 0.08 kcal/mol. Such a large negative value confirms reciprocal binding effects between the divalent cation and the substrate phosphoenolpyruvate. The binding of Mn2+ to the enzyme-ADP complex was also investigated and a free energy of coupling, delta(deltaF degrees) = -0.08 +/- 0.08 kcal/mol, was measured, indicative of little or no cooperativity in binding. The free energy of coupling with Mn2+ and pyruvate was measured as -1.52 +/- 0.14 kcal/mol, showing a significant amount of cooperativity in ligand binding but a substantially smaller effect than that observed for phosphoenolpyruvate binding. The magnitude of the coupling free energy may be related to the role of the divalent cation in the formation of the enzyme-substrate complexes. In the absence of the activating monovalent cation, the coupling free energies for phosphoenolpyruvate and pyruvate binding decrease by 40-60% and 25%, respectively, substantiating a role for the monovalent cation in the formation of enzyme-substrate complexes with phosphoenolpyruvate and with pyruvate.     

Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in six 22-base pair deoxyoligonucleotides.

Effects of different end sequences on stability, circular dichroism spectra (CD), and enzyme binding properties were investigated for six 22-base pair, non-self-complementary duplex DNA oligomers. The center sequences of these deoxyoligonucleotides have 8-14 base pairs in common and are flanked on both sides by sequences differing in context and A-T content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming a two-state melting transition. Melting free energies (20 degrees C) of the six DNAs evaluated from DSC experiments ranged from -18.7 to -32.7 kcal/mol. van't Hoff estimates of the free energies ranged from -18.5 to -48.0 kcal/mol. With either method, the trends in free energy as a function of sequence were identical. Equilibrium binding by BamHI restriction endonuclease to the 22-base pair DNAs was also investigated. The central eight base pairs of all six molecules, 5'-A-GGATCC-A-3', contained a BamHI recognition sequence bounded by A-T base pairs. Magnesium free binding assays were performed by titering BamHI against a constant concentration of each of the deoxyoligonucleotide substrates and analyzing reaction products by gel retardation. Binding isotherms of the total amount of bound DNA versus protein concentration were constructed which provided semiquantitative estimates of the equilibrium dissociation constants for dissociation of BamHI from the six DNA oligomers. Dissociation constants ranged from 0.5 x 10(-)(9) to 12.0 x 10(-)(9) M with corresponding binding free energies of -12.5 to -10.6 (+/-0. 1) kcal/mol. An inverse relationship is found when binding and stability are compared.     

Role of aspartate 27 of dihydrofolate reductase from Escherichia coli in interconversion of active and inactive enzyme conformers and binding of NADPH

The apoenzyme of wild-type (WT) dihydrofolate reductase (DHRF) from Escherichia coli exists in two conformational states, Et and Ew, which differ in affinity for NADPH and in kinetic competence. Dissociation constants for the binary complex of NADPH with the two conformers differ by over 100-fold (KDt = 0.17 microM, KDw = 22 microM). Rate constants governing the interconversion of conformers are small (t1/2 for Ew----Et = 71 s), and since Ew is not catalytically competent, this conversion is accompanied by an increase in catalytic velocity. The equilibrium proportion of Et in the absence of ligands is 63%, but binding of NADPH greatly increases this proportion, and t1/2 for conversion of Ew.NADPH to Et.NADPH is 30 s. This conformational equilibrium has also been examined in mutant enzyme in which aspartate 27 is replaced by asparagine (D27N E. coli DHFR). Although ASp27 is an active site residue, it does not interact directly with bound NADPH, and in the mutant the rate constant for NADPH binding to Et is unchanged as are the dissociation constants for NADPH complexes with Et or Ew. However, for mutant apoenzyme, the proportion of Et is decreased to 18% in the absence of ligands so that the overall KD for NADPH is increased (0.15 microM for WT E. coli DHFR, 0.68 microM for D27N E. coli DHFR). The lower proportion of Et is due to a decreased rate for Ew----Et (t1/2 = 221 s) and an increased rate for Et----Ew (t1/2 = 50 s versus 120 s for WT E. coli DHFR).     

Reaction of Escherichia coli and yeast photolyases with homogeneous short-chain oligonucleotide substrates

Similar rates have been observed for dimer repair with Escherichia coli photolyase and the heterogeneous mixtures generated by UV irradiation of oligothymidylates [UV-oligo(dT)n, n greater than or equal to 4] or DNA. Comparable stability was observed for ES complexes formed with UV-oligo(dT)n, (n greater than or equal to 9) or dimer-containing DNA. In this paper, binding studies with E. coli photolyase and a series of homogeneous oligonucleotide substrates (TpT, TpTp, pTpT, TpTpT, TpTpT, TpTpTpT, TpTpTpT, TpTpTpT, TpTpTpT) show that about 80% of the binding energy observed with DNA as substrate (delta G approximately 10 kcal/mol) can be attributed to the interaction of the enzyme with a dimer-containing region that spans only four nucleotides in length. This major binding determinant (TpTpTpT) coincides with the major conformational impact region of the dimer and reflects contributions from the dimer itself (TpT, delta G = 4.6 kcal/mol), adjacent phosphates (5'p, 0.8 kcal/mol; 3'p, 1.1 kcal/mol), and adjacent thymine residues (5'T, 0.8 kcal/mol; 3'T, 1.3 kcal/mol). Similar turnover rates (average kcat = 6.7 min-1) are observed with short-chain oligonucleotide substrates and UV-oligo(dT)18, despite a 25,000-fold variation in binding constants (Kd). In contrast, the ratio Km/Kd decreases as binding affinity decreases and appears to plateau at a value near 1. Turnover with oligonucleotide substrates occurs at a rate similar to that estimated for the photochemical step (5.1 min-1), suggesting that this step is rate determining. Under these conditions, Km will approach Kd when the rate of ES complex dissociation exceeds kcat.(ABSTRACT TRUNCATED AT 250 WORDS)     

Kinetics of the inhibition of bovine liver dihydrofolate reductase by tea catechins: origin of slow-binding inhibition and pH studies

Dihydrofolate reductase (DHFR) is the subject of intensive investigation since it appears to be the primary target enzyme for "antifolate" drugs, such as methotrexate and trimethoprim. Fluorescence quenching and stopped-flow fluorimetry show that the ester bond-containing tea polyphenols (-)-epigallocatechin gallate (EGCG) and (-)-epicatechin gallate (ECG) are potent and specific inhibitors of DHFR with inhibition constants (K(I)) of 120 and 82 nM, respectively. Both tea compounds showed the characteristics of slow-binding inhibitors of bovine liver DHFR. In this work, we have determined a complete kinetic scheme to explain the slow-binding inhibition and the pH effects observed during the inhibition of bovine liver DHFR by these tea polyphenols. Experimental data, based on fluorimetric titrations, and transient phase and steady-state kinetic studies confirm that EGCG and ECG are competitive inhibitors with respect to 7,8-dihydrofolate, which bind preferentially to the free form of the enzyme. The origin of their slow-binding inhibition is proposed to be the formation of a slow dissociation ternary complex by the reaction of NADPH with the enzyme-inhibitor complex. The pH controls both the ionization of critical catalytic residues of the enzyme and the protonation state of the inhibitors. At acidic pH, EGCG and ECG are mainly present as protonated species, whereas near neutrality, they evolve toward deprotonated species due to ionization of the ester-bonded gallate moiety (pK = 7.8). Although DHFR exhibits different affinities for the protonated and deprotonated forms of EGCG and ECG, it appears that the ionization state of Glu-30 in DHFR is critical for its inhibition. The physiological implications of these pH dependencies are also discussed.     

Energy of binding of Aspergillus oryzae beta-glucosidase with the substrate, and the mechanism of its enzymic action

Several beta-D-glucopyranosides (p-nitrophenyl, phenyl, and ethyl), 1-thio-beta-D-glucopyranosides, and phenyl 2-deoxy, 3-deoxy, 4-deoxy, and 6-deoxy beta-D-glucopyranosides were synthesized and used to study the mechanism of the enzymatic action of Taka-beta-glucosidase [EC 3.2.1.21 Aspergillus oryzae]. Kinetic constants of the enzyme for these glycosides were determined from S/V-S or 1/V-1/S plots, and the hydrolysis rates of these compounds with the enzyme, acid (3 N HCl) and alkali (3 N NaOH) were compared. Inhibition of the enzyme by 1,5-anhydroglucitol, glucal, dihydroglucal, and 1,6-anhydroglucopyranose was also examined. Glucal and 1,5-anhydroglucitol showed strong competitive inhibition. Free energy of binding of each hydroxyl group of glucosidic glucose with the enzyme was estimated from Kms of phenyl beta-glucoside and its deoxy analogues, and also Ki values of some inhibitors. The free energies of binding of 2-OH, 3-OH, 4-OH, and 6-OH were calculated to be 1.1, 2.4, 0.7, and 1.8 kcal/mol, respectively. The free energy of binding of phenoxide at C-1 (0.3 kcal/mol) was calculated from the Km of Ph-beta-Glc and Ki of 1,5-anhydroglucitol. The energy of binding of 5-CH2OH (2.3 kcal/mol) was obtained from the Km of Ph-beta-Glc and that of Ph-beta-Xyl. The sum (6.8 kcal/mol) of each partial binding free energy was close to the value of binding free energy of Ph-beta-Glc (7.0 kcal/mol) calculated by the equation; -delta Gbind = -RT ln Km-T delta Smix, showing that the methods of estimation of each binding energy used in the present study seemed reasonable. Glucal, having a pyranose form distorted slightly, showed strong competitive inhibition and the Ki of this inhibitor was smaller than the Km of Ph-beta-Glc, suggesting that the sugar ring bound to the active site was distored to a half chair form which is labile to acid hydrolysis.     

Human platelet factor 4 subunit association/dissociation thermodynamics and kinetics

Platelet factor 4 (PF4) monomers (7800 daltons) form dimers and tetramers in varying molar ratios under certain solution conditions [Mayo, K. H., & Chen, M. J. (1989) Biochemistry 28, 9469]. The presence of a simplified aromatic region (one Tyr and two His) and resolved monomer, dimer, and tetramer Y60 3,5 ring proton resonances makes study of PF4 aggregate association/dissociation thermodynamics and kinetics possible. PF4 protein subunit association/dissociation equilibrium thermodynamic parameters have been derived by 1H NMR (500MHz) resonance line-fitting analysis of steady-state Y60 3,5 ring proton resonance monomer-dimer-tetramer populations as a function of temperature from 10 to 40 degrees C. Below 10 degrees C and above 40 degrees C, resonance broadening and overlap severely impaired analysis. Enthalpic and entropic contributions to dimer association Gibb's free energy [-5.1 kcal/mol (30 degrees C)] are +2.5 +/- 1 kcal/mol and +26 +/- 7 eu, respectively, and for tetramer association Gibb's free energy [-5.7 kcal/mol (30 degrees C)], they are -7.5 +/- 1 kcal/mol and -7 +/- 3 eu, respectively. These thermodynamic parameters are consistent with low dielectric medium electrostatic/hydrophobic interactions governing dimer formation and hydrogen bonding governing tetramer formation. Association/dissociation kinetic parameters, i.e., steady-state jump rates, have been derived from exchange-induced line-width increases and from 1H NMR (500 MHz) saturation-transfer and spin-lattice (Tl) relaxation experiments. From dissociation jump rates and equilibrium constants, association rate constants were estimated. For dimer and tetramer equilibria at 30 degrees C, unimolecular dissociation rate constants are 35 +/- 10 s-1 for dimer dissociation and 6 +/- 2 s-1 for tetramer dissociation.(ABSTRACT TRUNCATED AT 250 WORDS)     

How fast monoamine oxidases decompose adrenaline? Kinetics of isoenzymes A and B evaluated by empirical valence bond simulation

This work scrutinizes kinetics of decomposition of adrenaline catalyzed by monoamine oxidase (MAO) A and B enzymes, a process controlling the levels of adrenaline in the central nervous system and other tissues. Experimental kinetic data for MAO A and B catalyzed decomposition of adrenaline are reported only in the form of the maximum reaction rate. Therefore, we estimated the experimental free energy barriers form the kinetic data of closely related systems using regression method, as was done in our previous study. By using multiscale simulation on the Empirical Valence Bond (EVB) level, we studied the chemical reactivity of the MAO A catalyzed decomposition of adrenaline and we obtained a value of activation free energy of 17.3 ± 0.4 kcal/mol. The corresponding value for MAO B is 15.7 ± 0.7 kcal/mol. Both values are in good agreement with the estimated experimental barriers of 16.6 and 16.0 kcal/mol for MAO A and MAO B, respectively. The fact that we reproduced the kinetic data and preferential catalytic effect of MAO B over MAO A gives additional support to the validity of the proposed hydride transfer mechanism. Furthermore, we demonstrate that adrenaline is preferably involved in the reaction in a neutral rather than in a protonated form due to considerably higher barriers computed for the protonated adrenaline substrate. The results are discussed in the context of chemical mechanism of MAO enzymes and possible applications of multiscale simulation to rationalize the effects of MAO activity on adrenaline level.     

Temperature dependence of anion transport inhibitor binding to human red cell membranes

The binding characteristics of the inhibitor of anion transport in human red cells, 4,4'-dibenzamido-2,2'-disulfonic stilbene (DBDS), to the anion transport protein of red cell ghost membranes in buffer containing 150 mM NaCl have been measured over the temperature range 0-30 degrees C by equilibrium and stopped-flow fluorescence methods. The equilibrium dissociation constant Keq, increased with temperature. No evidence of a 'break' in the ln(Keq) vs. 1/T plot was found. The standard dissociation enthalpy and entropy changes calculated from the temperature dependence are 9.1 +/- 0.9 kcal/mol and 3.2 +/- 0.3 e.u., respectively. Stopped-flow kinetic studies resolve the overall binding into two steps: a bimolecular association of DBDS with the anion transport protein, followed by a unimolecular rearrangement of the DBDS-protein complex. The rate constants for the individual steps in the binding mechanism can be determined from an analysis of the concentration dependence of the binding time course. Arrhenius plots of the rate constants showed no evidence of a break. Activation energies for the individual steps in the binding mechanism are 11.6 +/- 0.9 kcal/mol (bimolecular, forward step), 17 +/- 2 kcal/mol (bimolecular, reverse step), 6.4 +/- 2.3 kcal/mol (unimolecular, forward step), and 10.6 +/- 1.9 kcal/mol (unimolecular, reverse step). Our results indicate that there is an appreciable enthalpic energy barrier for the bimolecular association of DBDS with the transport protein, and appreciable enthalpic and entropic barriers for the unimolecular rearrangement of the DBDS-protein complex.     

Theoretical study of the ligand-CYP2B4 complexes: effect of structure on binding free energies and heme spin state

The molecular origins of temperature-dependent ligand-binding affinities and ligand-induced heme spin state conversion have been investigated using free energy analysis and DFT calculations for substrates and inhibitors of cytochrome P450 2B4 (CYP2B4), employing models of CYP2B4 based on CYP2C5(3LVdH)/CYP2C9 crystal structures, and the results compared with experiment. DFT calculations indicate that large heme-ligand interactions (ca. -15 kcal/mol) are required for inducing a high to low spin heme transition, which is correlated with large molecular electrostatic potentials (approximately -45 kcal/mol) at the ligand heteroatom. While type II ligands often contain oxygen and nitrogen heteroatoms that ligate heme iron, DFT results indicate that BP and MF heme complexes, with weak substrate-heme interactions (ca. -2 kcal/mol), and modest MEPS minima (>-35 kcal/mol) are high spin. In contrast, heme complexes of the CYP2B4 inhibitor, 4PI, the product of benzphetamine metabolism, DMBP, and water are low spin, have substantial heme-ligand interaction energies (<-15 kcal/mol) and deep MEPS minima (<-45 kcal/mol) near their heteroatoms. MMPBSA analysis of MD trajectories were made to estimate binding free energies of these ligands at the heme binding site of CYP2B4. In order to initially assess the realism of this approach, the binding free energy of 4PI inhibitor was computed and found to be a reasonable agreement with experiment: -7.7 kcal/mol [-7.2 kcal/mol (experiment)]. BP was determined to be a good substrate [-6.3 kcal/mol (with heme-ligand water), -7.3 kcal/mol (without ligand water)/-5.8 kcal/mol (experiment)], whereas the binding of MF was negligible, with only marginal binding binding free energy of -1.7 kcal/mol with 2-MF bound [-3.8 kcal/mol (experiment)], both with and without retained heme-ligand water. Analysis of the free energy components reveal that hydrophobic/nonpolar contributions account for approximately 90% of the total binding free energy of these substrates and are the source of their differential and temperature-dependent CYP2B4 binding. The results indicate the underlying origins of the experimentally observed differential binding affinities of BP and MF, and indicate the plausibility of the use of models derived from moderate sequence identity templates in conjunction with approximate free energy methods in the estimation of ligand-P450 binding affinities.     

Kinetics and thermodynamics of diacetyl reduction with NADPH by alpha-dicarbonyl reductase from pigeon liver

alpha-dicarbonyl reductase from pigeon liver catalyzes diacetyl reduction with NADPH via an ordered Bi-Bi mechanism in which the coenzyme is the leading substrate, as deduced from the inhibition pattern by products and by acetone. The activation energy of the reaction has been calculated as 16.6 kcal/mol, delta H and delta F as 15.6 and 15.3 kcal/mol, respectively, and delta S as 1 cal/mol per k. Kinetic constants obtained for substrates (KmNADPH = 15 microM, KsNADPH = 10 microM; Kmdiacetyl = 0.5 mM, Ksdiacetyl = 0.35 mM) and products (KiNADP 50 microM; Kiacetoin = 100 mM) are about 10 times lower than those reported for this enzyme in the reduction of diacetyl with NADH. This confirms that NADPH is its physiological coenzyme.     

Effects of temperature and viscosity on R67 dihydrofolate reductase catalysis

R67 dihydrofolate reductase (DHFR) is a novel homotetrameric protein that possesses 222 symmetry and a single, voluminous active site pore. This symmetry poses numerous limitations on catalysis; for example, two dihydrofolate (DHF) molecules or two NADPH molecules, or one substrate plus one cofactor can bind. Only the latter combination leads to catalysis. To garner additional information on how this enzyme facilitates transition-state formation, the temperature dependence of binding and catalysis was monitored. The binding of NADPH and DHF is enthalpy-driven. Previous primary isotope effect studies indicate hydride transfer is at least partially rate-determining. Accordingly, the activation energy associated with transition-state formation was measured and is found to be 6.9 kcal/mol (DeltaH(++)(25) = 6.3 kcal/mol). A large entropic component is also found associated with catalysis, TDeltaS(++)(25) = -11.3 kcal/mol. The poor substrate, dihydropteroate, binds more weakly than dihydrofolate (DeltaDeltaG = 1.4 kcal/mol) and displays a large loss in the binding enthalpy value (DeltaDeltaH = 3.8 kcal/mol). The k(cat) value for dihydropteroate reduction is decreased 1600-fold compared to DHF usage. This effect appears to derive mostly from the DeltaDeltaH difference in binding, demonstrating that the glutamate tail is important for catalysis. This result is surprising, as the para-aminobenzoyl-glutamate tail of DHF has been previously shown to be disordered by both NMR and crystallography studies. Viscosity studies were also performed and confirmed that the hydride transfer rate is not sensitive to sucrose addition. Surprisingly, binding of DHF, by both K(m) and K(d) determination, was found to be sensitive to added viscogens, suggesting a role for water in DHF binding.     

NMR structures of apo L. casei dihydrofolate reductase and its complexes with trimethoprim and NADPH: contributions to positive cooperative binding from ligand-induced refolding, conformational changes, and interligand hydrophobic interactions

In order to examine the origins of the large positive cooperativity (ΔG(0)(coop) = -2.9 kcal mol(-1)) of trimethoprim (TMP) binding to a bacterial dihydrofolate reductase (DHFR) in the presence of NADPH, we have determined and compared NMR solution structures of L. casei apo DHFR and its binary and ternary complexes with TMP and NADPH and made complementary thermodynamic measurements. The DHFR structures are generally very similar except for the A-B loop region and part of helix B (residues 15-31) which could not be directly detected for L. casei apo DHFR because of line broadening from exchange between folded and unfolded forms. Thermodynamic and NMR measurements suggested that a significant contribution to the cooperativity comes from refolding of apo DHFR on binding the first ligand (up to -0.95 kcals mol(-1) if 80% of A-B loop requires refolding). Comparisons of Cα-Cα distance differences and domain rotation angles between apo DHFR and its complexes indicated that generally similar conformational changes involving domain movements accompany formation of the binary complexes with either TMP or NADPH and that the binary structures are approaching that of the ternary complex as would be expected for positive cooperativity. These favorable ligand-induced structural changes upon binding the first ligand will also contribute significantly to the cooperative binding. A further substantial contribution to cooperative binding results from the proximity of the bound ligands in the ternary complex: this reduces the solvent accessible area of the ligand and provides a favorable entropic hydrophobic contribution (up to -1.4 kcal mol(-1)).     

Pterin deaminase from Bacillus megaterium. Purification and properties.

A pterin deaminase catalyzing the hydrolytic deamination of various pteridines was found in the bacterium, Bacillus megaterium, and partially purified from bacterial extract. The specific activity was raised 90-fold over that of the crude extract. The pH optimum is around 7.3, and the Km value for 6-carboxypterin is 1.3 mM. The molecular weight of the enzyme was estimated by gel filtration to be about 110,000. The enzyme deaminated pterin, 6-carboxypterin, biopterin, 6-methylpterin, 7-methylpterin, xanthopterin, 6-hydroxymethylpterin, sepiapterin, isosepiapterin, folic acid, and 6,7-dimethylpterin to their corresponding lumazines, whereas guanine, 7-carboxypterin, leucopterin, isoxanthopterin, and 6-methylisoxanthopterin did not serve as substrates. The enzyme was inhibited by PCMB and 8-azaguanine.     

Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.     

Nuclear magnetic resonance study of the state of protonation of inhibitors bound to mutant dihydrofolate reductase lacking the active-site carboxyl

13C nuclear magnetic resonance spectra have been obtained for complexes of [2-13C]methotrexate and [2-13C]trimethoprim with wild-type dihydrofolate reductase (DHFR) from Escherichia coli and with two mutant enzymes in which aspartic acid-27 is replaced by asparagine and by serine, respectively. In both the wild-type and mutated enzymes, exchange between the free inhibitor and the enzyme-complexed inhibitor is slow on the NMR time scale; hence, despite the considerably increased dissociation constants for binary complexes with the enzymes, the dissociation rate remains small relative to the frequency separation of the resonances. In all cases but one, the pKa of an inhibitor that is complexed to enzyme differs greatly from that of the free inhibitor. However, while the pKa of both inhibitors in complexes with the wild-type enzyme is elevated to above 10, the pKa of the inhibitors complexed with the Asn-27 and Ser-27 enzymes is lowered to a value below 4. Exact determinations of bound pKa values are limited by the solubility of the enzyme and the dissociation constants of the complexes. The single exception to these general conclusions is the ternary complex of the Ser-27 DHFR with trimethoprim and NADPH. In this complex, both free and enzyme-complexed trimethoprim exhibit similar pKa values (approximately equal to 7.6). However, both the exchange between free and enzyme-complexed inhibitor and the protonation of the enzyme-complexed inhibitor are slow in the NMR time scale, so that the spectra reveal three resonances corresponding to free inhibitor, to protonated enzyme-complexed inhibitor, and to unprotonated enzyme-complexed inhibitor.(ABSTRACT TRUNCATED AT 250 WORDS)     

Dissociation of minor groove binders from DNA: insights from metadynamics simulations

We have used metadynamics to investigate the mechanism of noncovalent dissociation from DNA by two representatives of alkylating and noncovalent minor groove (MG) binders. The compounds are anthramycin in its anhydrous form (IMI) and distamycin A (DST), which differ in mode of binding, size, flexibility and net charge. This choice enables to evaluate the influence of such factors on the mechanism of dissociation. Dissociation of IMI requires an activation free energy of approximately 12 kcal/mol and occurs via local widening of the MG and loss of contacts between the drug and one DNA strand, along with the insertion of waters in between. The detachment of DST occurs at a larger free energy cost, approximately 16.5 or approximately 18 kcal/mol depending on the binding mode. These values compare well with that of 16.6 kcal/mol extracted from stopped-flow experiments. In contrast to IMI, an intermediate is found in which the ligand is anchored to the DNA through its amidinium tail. From this conformation, binding and unbinding occur almost at the same rate. Comparison between DST and with kinetic models for the dissociation of Hoechst 33258 from DNA uncovers common characteristics across different classes of noncovalent MG ligands.