Effect of mutation on enzyme motion in dihydrofolate reductase

Hybrid quantum-classical molecular dynamics simulations of a mutant Escherichia coli dihydrofolate reductase enzyme are presented. Although residue 121 is on the exterior of the enzyme, experimental studies have shown that the mutation of Gly-121 to valine reduces the rate of hydride transfer by a factor of 163. The simulations indicate that the decrease in the hydride transfer rate for the G121V mutant is due to an increase in the free energy barrier. The calculated free energy barrier is higher for the mutant than for the wild-type enzyme by an amount that is consistent with the experimentally observed rate reduction. The calculated transmission coefficients are comparable for the wild-type and mutant enzymes. The simulations suggest that this mutation may interrupt a network of coupled promoting motions proposed to play an important role in DHFR catalysis. This phenomenon has broad implications for protein engineering and drug design.     

Tunneling and coupled motion in the Escherichia coli dihydrofolate reductase catalysis

H-transfer was studied in the complex kinetic cascade of dihydrofolate reductase. Intrinsic kinetic isotope effects, their temperature dependence, and other temperature-dependent parameters indicated H-tunneling, but no 1 degrees to 2 degrees coupled motion. The data also suggested environmentally coupled tunneling and commitment to catalysis on pre-steady-state isotope effects.     

How dihydrofolate reductase facilitates protonation of dihydrofolate

Dihydrofolate Reductase (DHFR) catalyzes the reduction of dihydrofolate (H2F) to tetrahydrofolate. On the basis of 10-12.5 ns molecular dynamics simulations of two conformations (closed and occluded) of the ternary DHFR/NADPH/H2F complex from Escherichia coli and a free energy perturbation approach, we have calculated the pKa value for the N5 atom in H2F. Our results suggest that the N5 atom in H2F is responsible for the pH dependency of the catalyzed reaction, meaning that DHFR facilitates protonation of H2F by approximately 4 pKa units. The mechanism behind this increase is due to favorable electrostatic interactions between the Asp27 residue and a proton at the N5 atom. The electrostatic interactions are enhanced by a hydrophobic active site, which to a large extent is made hydrophobic by the M20 loop in DHFR. Moreover, we find that the conformation imposed on H2F by DHFR to some extent also favors protonation of the N5 atom. Our results add support to previous findings and suggestions by Callender and co-workers [e.g., Deng, J.; Callender, R. J. Am. Chem. Soc. 1998, 120, 7730-7737] and explain why mutation of Asp27 may lead to severely reduced activity at neutral pH.     

Conformational substates modulate hydride transfer in dihydrofolate reductase

In earlier studies of the hydride-transfer reaction catalyzed by dihydrofolate reductase (DHFR) we identified features of the protein correlated with variations in the reaction barrier. We extend the scope of those studies by carrying out potential of mean force (PMF) simulations to determine the hydride-transfer barrier in the wild-type protein as well as the G121V and G121S mutants. While our prior studies focused on the reactant state, our current work addresses the full reaction pathway and directly probes the reactive event. The free energy barriers and structural ensembles resulting from these PMF calculations exhibit the same trends reported in our previous work. Fluctuations present in these simulations also exhibit trends associated with differences in the hydride-transfer barrier height. Moreover, vibrational modes anticipated to promote hydride transfer exhibit larger amplitudes in simulations that generate lowered barriers. The results of our study indicate that discrete basins (substates) on a potential energy landscape of the enzyme give rise to distinct hydride-transfer barriers. We suggest that the long-range effects of mutations at position 121 within DHFR are mediated by differentially preorganized protein environments in the context of distinct substate distributions, with concomitant changes to the dynamic properties of the enzyme.     

Carbon-deuterium bonds as probes of dihydrofolate reductase

Much effort has been directed toward understanding the contributions of electrostatics and dynamics to protein function and especially to enzyme catalysis. Unfortunately, these studies have been limited by the absence of direct experimental probes. We have been developing the use of carbon-deuterium bonds as probes of proteins and now report the application of the technique to the enzyme dihydrofolate reductase, which catalyzes a hydride transfer and has served as a paradigm for biological catalysis. We observe that the stretching absorption frequency of (methyl- d 3) methionine carbon-deuterium bonds shows an approximately linear dependence on solvent dielectric. Solvent and computational studies support the empirical interpretation of the stretching frequency in terms of local polarity. To begin to explore the use of this technique to study enzyme function and mechanism, we report a preliminary analysis of (methyl- d 3) methionine residues within dihydrofolate reductase. Specifically, we characterize the IR absorptions at Met16 and Met20, within the catalytically important Met20 loop, and Met42, which is located within the hydrophobic core of the enzyme. The results confirm the sensitivity of the carbon-deuterium bonds to their local protein environment, demonstrate that dihydrofolate reductase is electrostatically and dynamically heterogeneous, and lay the foundation for the direct characterization protein electrostatics and dynamics and, potentially, their contribution to catalysis.     

Identification of an anti-MRSA dihydrofolate reductase inhibitor from a diversity-oriented synthesis

The screening of a diversity-oriented synthesis library followed by structure-activity relationship investigations have led to the discovery of an anti-MRSA agent which operates as an inhibitor of Staphylococcus aureus dihydrofolate reductase.     

Multiple ligand-binding modes in bacterial R67 dihydrofolate reductase

Summary R67 dihydrofolate reductase (DHFR), a bacterial plasmid-encoded enzyme associated with resistance to the drug trimethoprim, shows neither sequence nor structural homology with the chromosomal DHFR. It presents a highly symmetrical toroidal structure, where four identical monomers contribute to the unique central active-site pore. Two reactants (dihydrofolate, DHF), two cofactors (NADPH) or one of each (R67•DHF•NADPH) can be found simultaneously within the active site, the last one being the reactive ternary complex. As the positioning of the ligands has proven elusive to empirical determination, we addressed the problem from a theoretical perspective. Several potential structures of the ternary complex were generated using the docking programs AutoDock and FlexX. The variability among the final poses, many of which conformed to experimental data, prompted us to perform a comparative scoring analysis and molecular dynamics simulations to assess the stability of the complexes. Analysis of ligand–ligand and ligand–protein interactions along the 4 ns trajectories of eight different structures allowed us to identify important inter-ligand contacts and key protein residues. Our results, combined with published empirical data, clearly suggest that multipe binding modes of the ligands are possible within R67 DHFR. While the pterin ring of DHF and the nicotinamide ring of NADPH assume a stacked endo-conformation at the centre of the pore, probably assisted by V66, Q67 and I68, the tails of the molecules extend towards opposite ends of the cavity, adopting multiple configurations in a solvent rich-environment where hydrogen-bond interactions with K32 and Y69 may play important roles.     

The role of large-scale motions in catalysis by dihydrofolate reductase

Dihydrofolate reductase has long been used as a model system to study the coupling of protein motions to enzymatic hydride transfer. By studying environmental effects on hydride transfer in dihydrofolate reductase (DHFR) from the cold-adapted bacterium Moritella profunda (MpDHFR) and comparing the flexibility of this enzyme to that of DHFR from Escherichia coli (EcDHFR), we demonstrate that factors that affect large-scale (i.e., long-range, but not necessarily large amplitude) protein motions have no effect on the kinetic isotope effect on hydride transfer or its temperature dependence, although the rates of the catalyzed reaction are affected. Hydrogen/deuterium exchange studies by NMR-spectroscopy show that MpDHFR is a more flexible enzyme than EcDHFR. NMR experiments with EcDHFR in the presence of cosolvents suggest differences in the conformational ensemble of the enzyme. The fact that enzymes from different environmental niches and with different flexibilities display the same behavior of the kinetic isotope effect on hydride transfer strongly suggests that, while protein motions are important to generate the reaction ready conformation, an optimal conformation with the correct electrostatics and geometry for the reaction to occur, they do not influence the nature of the chemical step itself; large-scale motions do not couple directly to hydride transfer proper in DHFR.     

Fluorescence resonance energy transfer between unnatural amino acids in a structurally modified dihydrofolate reductase

The cleavage of a substrate protein by HIV-1 protease has been monitored in real time by the use of a dihydrofolate reductase fusion protein in which a fluorescence donor and a fluorescence acceptor were introduced into sites flanking the HIV-1 protease cleavage site. The amino acids 7-azatryptophan and dabcyl-1,2-diaminopropionic acid were introduced into specific sites of the DHFR fusion protein in an in vitro protein biosynthesizing system using two misacylated suppressor tRNAs, each of which recognized a specific, unique codon introduced into the mRNA. Excitation of the fluorescence acceptor in the initially expressed protein afforded no light production, consistent with quenching by fluorescence resonance energy transfer. Treatment of the elaborated protein with HIV-1 protease cleaved the protein between the fluorescence donor and acceptor, affording a time-dependent increase in fluorescence that was equal in magnitude to that produced by admixture of a stoichiometric amount of free 7-azatryptophan to the solution containing the intact protein.     

高效毛细管电泳法测定二氢叶酸还原酶活力的研究

通过胶束电动毛细管电泳法研究分离二氢叶酸还原酶体系中二氢叶酸、四氢叶酸、 NADP、 NADPH和酶5种组分,在含0.002%Brij-35的pH 9.18 50 mmol/L 的硼砂缓冲溶液中,5种组分在18min内得到基线分离.通过对其产物四氢叶酸峰面积的定量测定,计算出二氢叶酸还原酶的米氏常数,建立了毛细管电泳法对二氢叶酸还原酶活力的测定方法.     

Evidence for environmentally coupled hydrogen tunneling during dihydrofolate reductase catalysis

Hydride transfer during catalysis by dihydrofolate reductase from Thermotoga maritima has been studied by stopped flow spectroscopy. The reduction of dihydrofolate by NADPH showed a biphasic temperature dependence of the deuterium kinetic isotope effect. At temperatures above 25 degrees C the KIE was temperature independent, while the reaction rates were strongly temperature dependent. Below 25 degrees C the KIE becomes dependent on temperature, and the ratio of the preexponential factors is inverse, suggesting a greater role for active dynamics that modulate the tunneling distance.     

Three-dimensional quantitative structure-activity and structure-selectivity relationships of dihydrofolate reductase inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) modelling using comparative molecular similarity indices analysis (CoMSIA) was applied to a series of 406 structurally diverse dihydrofolate reductase (DHFR) inhibitors from Pneumocystis carinii (pc) and rat liver (rl). X-ray crystal structures of three inhibitors bound to pcDHFR were used for defining the alignment rule. For pcDHFR, a QSAR model containing 6 components was selected using leave-10%-out cross-validation (n= 240, q2 = 0.65), while a 4-component model was selected for rlDHFR (n= 237, q2 = 0.63); both include steric, electrostatic and hydrophobic contributions. The models were validated using a large test set, designed to maximise its diversity and to verify the predictive accuracy of models for extrapolation. The pcDHFR model has r2 = 0.60 and mean absolute error (MAE) = 0.57 for the test set after removing 4 outliers, and the rlDHFR model has r2 = 0.60 and MAE = 0.69 after removing 4 test set outliers. In addition, classification models predicting selectivity for pcDHFR over rlDHFR were developed using soft independent modelling by class analogy (SIMCA), with a selectivity ratio of 2 (IC50,rlDHFR/ IC50,pcDHFR) used for delimiting classes. A 5-component model including steric and electrostatic contributions has cross-validated and test set classification rates of 0.67 and 0.68 for selective inhibitors, and 0.85 and 0.72 for unselective inhibitors. The predictive accuracy of models, together with the identification of important contributions in QSAR and classification models, offer the possibility of designing potent selective inhibitors and estimating their activity prior to synthesis.     

Radioenzymatic determination of subpicomole quantities of folic acid and dihydrofolate based on differential reduction by dihydrofolate reductase

A sensitive method for the determination of folic acid in the presence of dihydrofolate has been described. The analysis is based on the difference in reactivity of dihydrofolate reductase toward folic acid and dihydrofolate. The tetrahydrofolate formed was measured by stoichiometric entrapment into a covalent, ternary complex with thymidylate synthase and [³H]fluorodeoxyuridylate. Differentiation between folic acid and dihydrofolate was accomplished by the use of two levels of dihydrofolate reductase. Recovery of both folic acid and dihydrofolate was quantitative in the range of 0.1 to 2.0 pmol.     

Incorporating dynamics in E. coli dihydrofolate reductase enhances structure-based drug discovery

Escherichia coli dihydrofolate reductase (DHFR) is a long-standing target for enzyme studies. The influence of protein motion on its catalytic cycle is significant, and the conformation of the M20 loop is of particular interest. We present receptor-based pharmacophore models-an equivalent of solvent-mapping of binding hotspots-based on ensembles of protein conformations from molecular dynamics simulations of DHFR.NADPH in both the closed and open conformation of the M20 loop. The optimal models identify DHFR inhibitors over druglike non-inhibitors; furthermore, high-affinity inhibitors of E. coli DHFR are preferentially identified over general DHFR inhibitors. As expected, models resulting from simulations with DHFR in the productive conformation with a closed M20 loop have better performance than those from the open-loop simulations. Model performance improves with increased dynamic sampling, indicating that including a greater degree of protein flexibility can enhance the quest for potent inhibitors. This was compared to the limited conformational sampling seen in crystal structures, which were suboptimal for this application.     

Protein flexibility and species specificity in structure-based drug discovery: dihydrofolate reductase as a test system

In structure-based drug discovery, researchers would like to identify all possible scaffolds for a given target. However, techniques that push the boundaries of chemical space could lead to many false positives or inhibitors that lack specificity for the target. Is it possible to broadly identify the appropriate chemical space for the inhibitors and yet maintain target specificity? To address this question, we have turned to dihydrofolate reductase (DHFR), a well-studied metabolic enzyme of pharmacological relevance. We have extended our multiple protein structure (MPS) method for receptor-based pharmacophore models to use multiple X-ray crystallographic structures. Models were created for DHFR from human and Pneumocystis carinii. These models incorporate a fair degree of protein flexibility and are highly selective for known DHFR inhibitors over drug-like non-inhibitors. Despite sharing a highly conserved active site, the pharmacophore models reflect subtle differences between the human and P. carinii forms, which identify species-specific, high-affinity inhibitors. We also use structures of DHFR from Candida albicans as a counter example. The available crystal structures show little flexibility, and the resulting models give poorer performance in identifying species-specific inhibitors. Therapeutic success for this system may depend on achieving species specificity between the related human host and these key fungal targets. The MPS technique is a promising advance for structure-based drug discovery for DHFR and other proteins of biomedical interest.     

Comparison of coupled motions in Escherichia coli and Bacillus subtilis dihydrofolate reductase

Hybrid quantum/classical molecular dynamics simulations are used to compare the role of protein motion in the hydride transfer reaction catalyzed by Escherichia coli and Bacillus subtilis dihydrofolate reductase (DHFR). These two enzymes have 44% sequence identity, and the experimentally determined structures and hydride transfer rates are similar. The simulations indicate that the tertiary structures of both enzymes evolve in a similar manner during the hydride transfer reaction. In both enzymes, the donor-acceptor distance decreases to approximately 2.7 Angstroms at the transition state configurations to enable hydride transfer. Zero point energy and hydrogen tunneling effects are found to be significant for both enzymes. Covariance and rank correlation analyses of motions throughout the protein and ligands illustrate that E. coli and B. subtilis DHFR exhibit both similarities and differences in the equilibrium fluctuations and the conformational changes along the collective reaction coordinate for hydride transfer. A common set of residues that play a significant role in the network of coupled motions leading to configurations conducive to hydride transfer for both E. coli and B. subtilis DHFR was identified. These results suggest a balance between conservation and flexibility in the thermal motions and conformational changes during hydride transfer. Homologous protein structures, in conjunction with conformational sampling, enable enzymes with different sequences to catalyze the same hydride transfer reaction with similar efficiency.     

Structure-guided development of efficacious antifungal agents targeting Candida glabrata dihydrofolate reductase

Candida glabrata is a lethal fungal pathogen resistant to many antifungal agents and has emerged as a critical target for drug discovery. Over the past several years, we have been developing a class of propargyl-linked antifolates as antimicrobials and hypothesized that these compounds could be effective inhibitors of dihydrofolate reductase (DHFR) from C. glabrata. We initially screened a small collection of these inhibitors and found modest levels of potency. Subsequently, we determined the crystal structure of C. glabrata DHFR bound to a representative inhibitor with data to 1.6 A resolution. Using this structure, we designed and synthesized second-generation inhibitors. These inhibitors bind the C. glabrata DHFR enzyme with subnanomolar potency, display greater than 2000-fold levels of selectivity over the human enzyme, and inhibit the growth of C. glabrata at levels observed with clinically employed therapeutics.