Functional role for Tyr 31 in the catalytic cycle of chicken dihydrofolate reductase

Despite much work, many key aspects of the mechanism of the dihydrofolate reductase (DHFR) catalyzed reduction of dihydrofolate remain unresolved. In bacterial forms of DHFR both substrate and water access to the active site are controlled by the conformation of the mobile M20 loop. In vertebrate DHFRs only one conformation of the residues corresponding to the M20 loop has been observed. Access to the active site was proposed to be controlled by residue 31. MD simulations of chicken DHFR complexed with substrates and cofactor revealed a closing of the side chain of Tyr 31 over the active site on binding of dihydrofolate. This conformational change was dependent on the presence of glutamate on the para-aminobenzoylamide moiety of dihydrofolate. In its absence, the conformation remained open. Although water could enter the active site and hydrogen bond to N5 of dihydrofolate, indicating the feasibility of water as the proton donor, this was not controlled by the conformation of Tyr 31. The water accessibility of the active site was low for both conformations of Tyr 31. However, when hydride was transferred from NADPH to C6 of dihydrofolate before protonation, the average time during which water was found in hydrogen bonding distance to N5 of dihydrofolate in the active site increased almost fivefold. These results indicated that water can serve as the Broensted acid for the protonation of N5 of dihydrofolate during the DHFR catalyzed reduction.     

The role of enzyme dynamics and tunnelling in catalysing hydride transfer: studies of distal mutants of dihydrofolate reductase

Residues M42 and G121 of Escherichia coli dihydrofolate reductase (ecDHFR) are on opposite sides of the catalytic centre (15 and 19 A away from it, respectively). Theoretical studies have suggested that these distal residues might be part of a dynamics network coupled to the reaction catalysed at the active site. The ecDHFR mutant G121V has been extensively studied and appeared to have a significant effect on rate, but only a mild effect on the nature of H-transfer. The present work examines the effect of M42W on the physical nature of the catalysed hydride transfer step. Intrinsic kinetic isotope effects (KIEs), their temperature dependence and activation parameters were studied. The findings presented here are in accordance with the environmentally coupled hydrogen tunnelling. In contrast to the wild-type (WT), fluctuations of the donor-acceptor distance were required, leading to a significant temperature dependence of KIEs and deflated intercepts. A comparison of M42W and G121V to the WT enzyme revealed that the reduced rates, the inflated primary KIEs and their temperature dependences resulted from an imperfect potential surface pre-arrangement relative to the WT enzyme. Apparently, the coupling of the enzyme's dynamics to the reaction coordinate was altered by the mutation, supporting the models in which dynamics of the whole protein is coupled to its catalysed chemistry.     

Conformational selection and induced changes along the catalytic cycle of Escherichia coli dihydrofolate reductase

Protein function often involves changes between different conformations. Central questions are how these conformational changes are coupled to the binding or catalytic processes during which they occur, and how they affect the catalytic rates of enzymes. An important model system is the enzyme dihydrofolate reductase (DHFR) from Escherichia coli, which exhibits characteristic conformational changes of the active‐site loop during the catalytic step and during unbinding of the product. In this article, we present a general kinetic framework that can be used (1) to identify the ordering of events in the coupling of conformational changes, binding, and catalysis and (2) to determine the rates of the substeps of coupled processes from a combined analysis of nuclear magnetic resonance R₂ relaxation dispersion experiments and traditional enzyme kinetics measurements. We apply this framework to E. coli DHFR and find that the conformational change during product unbinding follows a conformational‐selection mechanism, that is, the conformational change occurs predominantly prior to unbinding. The conformational change during the catalytic step, in contrast, is an induced change, that is, the change occurs after the chemical reaction. We propose that the reason for these conformational changes, which are absent in human and other vertebrate DHFRs, is robustness of the catalytic rate against large pH variations and changes to substrate/product concentrations in E. coli. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

Cooperative effects of potassium, magnesium, and magnesium-ADP on the release of Escherichia coli dihydrofolate reductase from the chaperonin GroEL.

Previous investigation has shown that at 22 degrees C and in the presence of the chaperonin GroEL, the slowest step in the refolding of Escherichia coli dihydrofolate reductase (EcDHFR) reflects release of a late folding intermediate from the cavity of GroEL (Clark AC, Frieden C, 1997, J Mol Biol 268:512-525). In this paper, we investigate the effects of potassium, magnesium, and MgADP on the release of the EcDHFR late folding intermediate from GroEL. The data demonstrate that GroEL consists of at least two conformational states, with apparent rate constants for EcDHFR release that differ by four- to fivefold. In the absence of potassium, magnesium, and ADP, approximately 80-90% of GroEL resides in the form with the faster rate of release. Magnesium and potassium both shift the distribution of GroEL forms toward the form with the slower release rate, though cooperativity for the magnesium-induced transition is observed only in the presence of potassium. MgADP at low concentrations (0-50 microM) shifts the distribution of GroEL forms toward the form with the faster release rate, and this effect is also potassium dependent. Nearly identical results were obtained with a GroEL mutant that forms only a single ring, demonstrating that these effects occur within a single toroid of GroEL. In the presence of saturating magnesium, potassium, and MgADP, the apparent rate constant for the release of EcDHFR from wild-type GroEL at 22 degrees C reaches a limiting value of 0.014 s(-1). For the single ring mutant of GroEL, the rate of EcDHFR release under the same conditions reaches a limiting value of 0.024 s(-1), suggesting that inter-ring negative cooperativity exists for MgADP-induced substrate release. The data suggest that MgADP preferentially binds to one conformation of GroEL, that with the faster apparent rate constant for EcDHFR release, and induces a conformational change leading to more rapid release of substrate protein.     

The EF-hand domain: A globally cooperative structural unit

EF-hand Ca(2+)-binding proteins participate in both modulation of Ca(2+) signals and direct transduction of the ionic signal into downstream biochemical events. The range of biochemical functions of these proteins is correlated with differences in the way in which they respond to the binding of Ca(2+). The EF-hand domains of calbindin D(9k) and calmodulin are homologous, yet they respond to the binding of calcium ions in a drastically different manner. A series of comparative analyses of their structures enabled the development of hypotheses about which residues in these proteins control the calcium-induced changes in conformation. To test our understanding of the relationship between protein sequence and structure, we specifically designed the F36G mutation of the EF-hand protein calbindin D(9k) to alter the packing of helices I and II in the apoprotein. The three-dimensional structure of apo F36G was determined in solution by nuclear magnetic resonance spectroscopy and showed that the design was successful. Surprisingly, significant structural perturbations also were found to extend far from the site of mutation. The observation of such long-range effects provides clear evidence that four-helix EF-hand domains should be treated as a single globally cooperative unit. A hypothetical mechanism for how the long-range effects are transmitted is described. Our results support the concept of energetic and structural coupling of the key residues that are crucial for a protein's fold and function.     

Fluorescence tools to investigate riboswitch structural dynamics

Riboswitches are novel regulatory elements that respond to cellular metabolites to control gene expression. They are constituted of highly conserved domains that have evolved to recognize specific metabolites. Such domains, so-called aptamers, are folded into intricate structures to enable metabolite recognition. Over the years, the development of ensemble and single-molecule fluorescence techniques has allowed to probe most of the mechanistic aspects of aptamer folding and ligand binding. In this review, we summarize the current fluorescence toolkit available to study riboswitch structural dynamics. We fist describe those methods based on fluorescent nucleotide analogues, mostly 2-aminopurine (2AP), to investigate short-range conformational changes, including some key steady-state and time-resolved examples that exemplify the versatility of fluorescent analogues as structural probes. The study of long-range structural changes by Förster resonance energy transfer (FRET) is mostly discussed in the context of single-molecule studies, including some recent developments based on the combination of single-molecule FRET techniques with controlled chemical denaturation methods. This article is part of a Special Issue entitled: Riboswitches.     

Conformational change of the methionine 20 loop of Escherichia coli dihydrofolate reductase modulates pKa of the bound dihydrofolate

We evaluate the pK(a) of dihydrofolate (H(2)F) at the N(5) position in three ternary complexes with Escherichia coli dihydrofolate reductase (ecDHFR), namely ecDHFR(NADP(+):H(2)F) in the closed form (1), and the Michaelis complexes ecDHFR(NADPH:H(2)F) in the closed (2) and occluded (3) forms, by performing free energy perturbation with molecular dynamics simulations (FEP/MD). Our simulations suggest that in the Michaelis complex the pK(a) is modulated by the Met20 loop fluctuations, providing the largest pK(a) shift in substates with a "tightly closed" loop conformation; in the "partially closed/open" substates, the pK(a) is similar to that in the occluded complex. Conducive to the protonation, tightly closing the Met20 loop enhances the interactions of the cofactor and the substrate with the Met20 side chain and aligns the nicotinamide ring of the cofactor coplanar with the pterin ring of the substrate. Overall, the present study favors the hypothesis that N(5) is protonated directly from solution and provides further insights into the mechanism of the substrate protonation.     

Structure and hydride transfer mechanism of a moderate thermophilic dihydrofolate reductase from Bacillus stearothermophilus and comparison to its mesophilic and hyperthermophilic homologues

Dihydrofolate reductase (DHFR) from a moderate thermophilic organism, Bacillus stearothermophilus, has been cloned and expressed. Physical characterization of the protein (BsDHFR) indicates that it is a monomeric protein with a molecular mass of 18,694.6 Da (0.8), coincident with the mass of 18 694.67 Da calculated from the primary sequence. Determination of the X-ray structure of BsDHFR provides the first structure for a monomeric DHFR from a thermophilic organism, indicating a high degree of conservation of structure in relation to all chromosomal DHFRs. Structurally based sequence alignment of DHFRs indicates the following levels of sequence identity and similarity for BsDHFR: 38 and 58% with Escherichia coli, 35 and 56% with Lactobacillus casei, and 23 and 40% with Thermotoga maritima, respectively. Steady state kinetic isotope effect studies indicate an ordered kinetic mechanism at elevated temperatures, with NADPH binding first to the enzyme. This converts to a more random mechanism at reduced temperatures, reflected in a greatly reduced K(m) for dihydrofolate at 20 degrees C in relation to that at 60 degrees C. A reduction in either temperature or pH reduces the degree to which the hydride transfer step is rate-determining for the second-order reaction of DHF with the enzyme-NADPH binary complex. Transient state kinetics have been used to study the temperature dependence of the isotope effect on hydride transfer at pH 9 between 10 and 50 degrees C. The data support rate-limiting hydride transfer with a moderate enthalpy of activation (E(a) = 5.5 kcal/mol) and a somewhat greater temperature dependence for the kinetic isotope effect than predicted from classical behavior [A(H)/A(D) = 0.57 (0.15)]. Comparison of kinetic parameters for BsDHFR to published data for DHFR from E. coli and T. maritima shows a decreasing trend in efficiency of hydride transfer with increasing thermophilicity of the protein. These results are discussed in the context of the capacity of each enzyme to optimize H-tunneling from donor (NADPH) to acceptor (DHF) substrates.     

Binding of fexofenadine hydrochloride to bovine serum albumin: structural considerations by spectroscopic techniques and molecular docking

The binding interaction of peripheral H₁ receptor antagonist drug, fexofenadine hydrochloride to bovine serum albumin (BSA) is investigated by fluorescence spectroscopy in combination with UV-absorption spectroscopy under physiological conditions. The Stern–Volmer plots at different temperatures and the steady state fluorescence suggested a static type of interaction between fexofenadine and BSA. Binding constants were determined to provide a measure of the binding affinity between fexofenadine and BSA. It was found that BSA has one binding site for fexofenadine. On the basis of the competitive site marker experiments and thermodynamic results, it was considered that fexofenadine was primarily bound to the site I of BSA mainly by hydrogen bond and van der Waals force. Utilising Förster resonance energy transfer the distance, r between the donor, BSA and acceptor fexofenadine was obtained. Furthermore, the results of circular dichroism and synchronous fluorescence spectrum indicated that the secondary structure of BSA was changed in the presence of fexofenadine. Molecular docking was applied to further define the interaction of fexofenadine with BSA.     

Crystal structure of reduced thioredoxin reductase from Escherichia coli: structural flexibility in the isoalloxazine ring of the flavin adenine dinucleotide cofactor

Catalysis by thioredoxin reductase (TrxR) from Escherichia coli requires alternation between two domain arrangements. One of these conformations has been observed by X-ray crystallography (Waksman G, Krishna TSR, Williams CH Jr, Kuriyan J, 1994, J Mol Biol 236:800-816). This form of TrxR, denoted FO, permits the reaction of enzyme-bound reduced FAD with a redox-active disulfide on TrxR. As part of an investigation of conformational changes and intermediates in catalysis by TrxR, an X-ray structure of the FO form of TrxR with both the FAD and active site disulfide reduced has been determined. Reduction after crystallization resulted in significant local conformation changes. The isoalloxazine ring of the FAD cofactor, which is essentially planar in the oxidized enzyme, assumes a 34 degree "butterfly" bend about the N(5)-N(10) axis in reduced TrxR. Theoretical calculations reported by others predict ring bending of 15-28 degrees for reduced isoalloxazines protonated at N(1). The large bending in reduced TrxR is attributed in part to steric interactions between the isoalloxazine ring and the sulfur of Cys138, formed by reduction of the active site disulfide, and is accompanied by changes in the positions and interactions of several of the ribityl side-chain atoms of FAD. The bending angle in reduced TrxR is larger than that for any flavoprotein in the Protein Data Bank. Distributions of bending angles in published oxidized and reduced flavoenzyme structures are different from those found in studies of free flavins, indicating that the protein environment has a significant effect on bending.     

Protein conformational transitions coupled to binding in molecular recognition of unstructured proteins: hierarchy of structural loss from all-atom Monte Carlo simulations of p27Kip1 unfolding-unbinding and structural determinants of the binding mechanism

Conformational transitions coupled to binding are studied for the p27(Kip1) protein which undergoes a functional disorder-to-order folding transition during tertiary complex formation with the phosphorylated cyclin A-cyclin-dependent kinase 2 (Cdk2) binary complex. Temperature-induced Monte Carlo simulations of p27(Kip1) unfolding-unbinding carried out from the crystal structure of the tertiary complex have revealed a systematic trend in the hierarchy of structural loss for p27(Kip1) and a considerable difference in mobility of p27(Kip1) secondary structure elements. The most persistent interactions of p27(Kip1) at the intermolecular interface during unfolding-unbinding simulations are formed by beta-hairpin and beta-strand that on average maintain their structural integrity considerably longer than other p27(Kip1) elements. We have found that the ensemble of unfolded p27(Kip1) conformations is characterized by transitions between mostly unbound, collapsed conformations and entropically favorable p27(Kip1) conformations, which are weakly bound to the cyclin A side of the binary complex. The results of this study are consistent with the experimental evidence pointing to this region of the intermolecular interface as a potential initiation docking site during binding reaction and may reconcile conflicting experimental hypotheses on the recognition of substrate recruitment motifs.     

The cleavable prepiece of an imported mitochondrial protein is sufficient to direct cytosolic dihydrofolate reductase into the mitochondrial matrix

The cleavable prepiece of the precursor to yeast cytochrome c oxidase subunit IV (an imported mitochondrial protein) was attached to the amino-terminus of mouse dihydrofolate reductase (a cytosolic protein) by gene fusion. The resulting fusion protein was imported into the matrix of isolated, energized yeast mitochondria and cleaved to a polypeptide whose size was similar to that of authentic dihydrofolate reductase.     

Analysis of hydride transfer and cofactor fluorescence decay in mutants of dihydrofolate reductase: possible evidence for participation of enzyme molecular motions in catalysis.

A remarkable correlation has been discovered between fluorescence lifetimes of bound NADPH and rates of hydride transfer among mutants of dihydrofolate reductase (DHFR) from Escherichia coli. Rates of hydride transfer from NADPH to dihydrofolate change by a factor of 1,000 for the series of mutant enzymes. Since binding constants for the initial complex between coenzyme and DHFR change by only a factor of 10, the major portion of the change in hydride transfer must be attributed to losses in transition-state stabilization. The time course of fluorescence decay for NADPH bound to DHFR is biphasic. Lifetimes ranging from 0.3 to 0.5 ns are attributed to a solvent-exposed dihydronicotinamide conformation of bound coenzyme which is presumably not active in catalysis, while decay times (tau 2) in the range of 1.3 to 2.3 ns are assigned to a more tightly bound species of NADPH in which dihydronicotinamide is sequestered from solvent. It is this slower component that is of interest. Ternary complexes with three different inhibitors, methotrexate, 5-deazafolate, and trimethoprim, were investigated, along with the holoenzyme complex; 3-acetylNADPH was also investigated. Fluorescence polarization decay, excitation polarization spectra, the temperature variation of fluorescence lifetimes, fluorescence amplitudes, and wavelength of absorbance maxima were measured. We suggest that dynamic quenching or internal conversion promotes decay of the excited state in NADPH-DHFR. When rates of hydride transfer are plotted against the fluorescence lifetime (tau 2) of tightly bound NADPH, an unusual correlation is observed. The fluorescence lifetime becomes longer as the rate of catalysis decreases for most mutants studied. However, the fluorescence lifetime is unchanged for those mutations that principally alter the binding of dihydrofolate while leaving most dihydronicotinamide interactions relatively undisturbed. The data are interpreted in terms of possible dynamic motions of a flexible loop region in DHFR which closes over both substrate and coenzyme binding sites. These motions could lead to faster rates of fluorescence decay in holoenzyme complexes and, when correlated over time, may be involved in other motions which give rise to enhanced rates of catalysis in DHFR.     

The noncatalytic triad of alpha-amylases: a novel structural motif involved in conformational stability

Chloride-activated alpha-amylases contain a noncatalytic triad, independent of the glycosidic active site, perfectly mimicking the catalytic triad of serine-proteases and of other active serine hydrolytic enzymes. Mutagenesis of Glu, His, and Ser residues in various alpha-amylases shows that this pattern is a structural determinant of the enzyme conformation that cannot be altered without losing the intrinsic stability of the protein. (1)H-(15)N NMR spectra of a bacterial alpha-amylase reveal proton signals that are identical with the NMR signature of catalytic triads and especially a deshielded proton involving a protonated histidine and displaying properties similar to that of a low barrier hydrogen bond. It is proposed that the H-bond between His and Glu of the noncatalytic triad is an unusually strong interaction, responsible for the observed NMR signal and for the weak stability of the triad mutants. Furthermore, a stringent template-based search of the Protein Data Bank demonstrated that this motif is not restricted to alpha-amylases, but is also found in 80 structures from 33 different proteins, amongst which SH2 domain-containing proteins are the best representatives.     

A hydrophobic cluster forms early in the folding of dihydrofolate reductase

The rapid kinetic phase that leads from unfolded species to transient folding intermediates in dihydrofolate reductase from Escherichia coli was examined by site-directed mutagenesis and by physicochemical means. The absence of this fluorescence-detected phase in the refolding of the Trp-74Phe mutant protein strongly implies that this early phase in refolding can be assigned to just one of the five Trp residues in the protein, Trp-74. In addition, water-soluble fluorescence quenching agents, iodide and cesium, have a much less significant effect on this early step in refolding than on the slower phases that lead to native and native-like conformers. These and other data imply that an important early event in the folding of dihydrofolate reductase is the formation of a hydrophobic cluster which protects Trp-74 from solvent.     

The influence of glycerol on the binding of methotrexate to Mycobacterium tuberculosis dihydrofolate reductase: a molecular modelling study

The enzyme, dihydrofolate reductase (DHFR), from Mycobacterium tuberculosis (mt-DHFR) is believed to be a potential drug target for the treatment of tuberculosis. The co-crystal structure of mt-DHFR bound with glycerol (GOL), NAPDH and methotrexate (MTX) reveals a GOL binding site on the enzyme. This GOL binding site could be very important for the design of novel, selective mt-DHFR inhibitors, because this binding site is absent on human DHFR (h-DHFR). We have performed molecular dynamic simulations and free energy calculations to evaluate the binding affinity of GOL and its free energy contribution to the binding of MTX to mt-DHFR. The results showed that GOL does not bind tightly to mt-DHFR. Although GOL itself contributed free energy on MTX binding to mt-DHFR, GOL also increased the flexibilities of MTX, so that MTX cannot maintain strong electronic interactions with ARG32 and ARG60, which caused the total binding free energy to decrease. These data suggest that GOL binding is weak and it could be expelled from the binding site, to allow inhibitors containing appropriate side chains to bind. This observation can be used to inform future drug design studies, especially those aimed at improving drug selectivity against h-DHFR.     

A flexible loop in tyrosine hydroxylase controls coupling of amino acid hydroxylation to tetrahydropterin oxidation

The role of a polypeptide loop in tyrosine hydroxylase (TyrH) whose homolog in phenylalanine hydroxylase (PheH) takes on a different conformation when substrates are bound has been studied using site-directed mutagenesis. The loop spans positions 177 to 191; alanine was introduced into those positions, introducing one alanine substitution per TyrH variant. Mutagenesis of residues in the center of the loop resulted in alterations in the KM values for substrates, the Vmax value for dihydroxyphenylalanine (DOPA) synthesis, and the coupling of tetrahydropterin oxidation to tyrosine hydroxylation. The variant with the most altered KM value for 6-methyltetrahydropterin was TyrH F184A. The variants with the most affected K(tyr) values were those with substitutions in the center of the loop, TyrH K183A, F184A, D185A, P186A and D187A. These five variants also had the most reduced Vmax values for DOPA synthesis. Alanine substitution in positions 182-186 resulted in lowered ratios of tyrosine hydroxylation to tetrahydropterin oxidation. TyrH F184Y and PheH Y138F, variants with the residue at the center of the loop substituted with the residue present at the homologous position in the other hydroxylase, were also studied. The V/K(tyr) to V/K(phe) ratios for these variants were altered significantly, but the results did not suggest that F184 of TyrH or Y138 of PheH plays a dominant role in determining amino acid substrate specificity.     

Effects of the difference in the unfolded-state ensemble on the folding of Escherichia coli dihydrofolate reductase

The unfolded state of a protein is an ensemble of a large number of conformations ranging from fully extended to compact structures. To investigate the effects of the difference in the unfolded-state ensemble on protein folding, we have studied the structure, stability, and folding of "circular" dihydrofolate reductase (DHFR) from Escherichia coli in which the N and C-terminal regions are cross-linked by a disulfide bond, and compared the results with those of disulfide-reduced "linear" DHFR. Equilibrium studies by circular dichroism, difference absorption spectra, solution X-ray scattering, and size-exclusion chromatography show that whereas the native structures of both proteins are essentially the same, the unfolded state of circular DHFR adopts more compact conformations than the unfolded state of the linear form, even with the absence of secondary structure. Circular DHFR is more stable than linear DHFR, which may be due to the decrease in the conformational entropy of the unfolded state as a result of circularization. Kinetic refolding measurements by stopped-flow circular dichroism and fluorescence show that under the native conditions both proteins accumulate a burst-phase intermediate having the same structures and both fold by the same complex folding mechanism with the same folding rates. Thus, the effects of the difference in the unfolded state of circular and linear DHFRs on the refolding reaction are not observed after the formation of the intermediate. This suggests that for the proteins with close termini in the native structure, early compaction of a protein molecule to form a specific folding intermediate with the N and C-terminal regions in close proximity is a crucial event in folding. If there is an enhancement in the folding reflecting the reduction in the breadth of the unfolded-state ensemble for circular DHFR, this acceleration must occur in the sub-millisecond time-range.     

Quantifying structural changes and stoichiometry of protein interactions using size and density profiling

Summary The use of Dual Polarisation Interferometry, and emerging analytical biophysical technique, is described for the determination of the optogeometrical properties (thickness and density) at high resolution of adsorbed protein layers at the solid-liquid interface. The technique has been used to quantify, in real time and at subatomic resolution, the structural changes occurring in two well-characterised protein interaction systems, an antibody-antigen interaction and the biotin-streptavidin interaction. The realtime data obtained on structural changes during the interactions is in excellent agreement with previously reported X-ray crystallography and neutron reflection data. The precision of the measurements taken was of the order of 0.01 nm with respect to protein size. The dual-parameter approach also allowed the stoichiometry of both of these interactions to be calculated, giving values that confirm the current understanding of the interactions. This approach provides detailed insights into the inherent and subtle link between structural change and function in proteins, to a degree not previously possible through mass change measurements alone. The technique is expected to find utility in the increasingly important study of protein structure and function.