Crystal structure of human dihydrofolate reductase complexed with folate

Ubiquitylation-independent cotranslational degradation of dihydrofolate reductase and ubiquitin

Nascent proteins are degraded during or immediately after synthesis, a process called cotranslational protein degradation (CTPD). Although CTPD was observed decades ago, it has never been fully explored mechanistically and functionally. We show here that dihydrofolate reductase (DHFR) and ubiquitin (Ub), two stable proteins widely used in protein degradation studies, are actually subject to CTPD. Unlike canonical posttranslational protein degradation, CTPD of DHFR and Ub does not require prior ubiquitylation. Our data also suggest that protein expression level and N-terminal folding pattern may be two critical determinants for CTPD. Thus, this study reveals that CTPD plays a role in regulating the homeostasis of long-lived proteins and provides insights into the mechanism of CTPD.

X-ray-driven chemistry and conformational heterogeneity in atomic resolution crystal structures of bacterial dihydrofolate reductases

Dihydrofolate reductase (DHFR) catalyzes the NADPH-dependent reduction of dihydrofolate to tetrahydrofolate. Bacterial DHFRs are targets of several important antibiotics as well as model enzymes for the role of protein conformational dynamics in enzyme catalysis. We collected 0.93 Å resolution X-ray diffraction data from both Bacillus subtilis (Bs) and E. coli (Ec) DHFRs bound to folate and NADP+. These oxidized ternary complexes should not be able to perform chemistry, however electron density maps suggest hydride transfer is occurring in both enzymes. Comparison of low- and high-dose EcDHFR datasets show that X-rays drive partial production of tetrahydrofolate. Hydride transfer causes the nicotinamide moiety of NADP+ to move towards the folate as well as correlated shifts in nearby residues. Higher radiation dose also changes the conformational heterogeneity of Met20 in EcDHFR, supporting a solvent gating role during catalysis. BsDHFR has a different pattern of conformational heterogeneity and an unexpected disulfide bond, illustrating important differences between bacterial DHFRs. This work demonstrates that X-rays can drive hydride transfer similar to the native DHFR reaction and that X-ray photoreduction can be used to interrogate catalytically relevant enzyme dynamics in favorable cases.

In silico screening of inhibitors against human dihydrofolate reductase to identify potential anticancer compounds

In all species, dihydrofolate reductase (DHFR) is an essential enzyme that regulates the cellular amount of tetrahydrofolate. Human DHFR (hDHFR) activity inhibition results in tetrahydrofolate depletion and cell death. This property has made hDHFR a therapeutic target for cancer. Methotrexate is a well-known hDHFR inhibitor, but its administration has shown some light to severe adverse effects. Therefore, we aimed to find new potential hDHFR inhibitors using structure-based virtual screening, ADMET prediction, molecular docking, and molecular dynamics simulations. Here, we used the PubChem database to find all compounds with at least 90% structural similarity to known natural DHFR inhibitors. To explore their interaction pattern and estimate their binding affinities, the screened compounds (2023) were subjected to structure-based molecular docking against hDHFR. The fifteen compounds that showed higher binding affinity to the hDHFR than the reference compound (methotrexate) displayed important molecular orientation and interactions with key residues in the enzyme's active site. These compounds were subjected to Lipinski and ADMET prediction. PubChem CIDs: 46886812 and 638190 were identified as putative inhibitors. In addition, molecular dynamics simulations revealed that the binding of compounds (CIDs: 46886812 and 63819) stabilized the hDHFR structure and caused minor conformational changes. Our findings suggest that two compounds (CIDs: 46886812 and 63819) could be promising potential inhibitors of hDHFR in cancer therapy.Communicated by Ramaswamy H. Sarma.

Cycloguanil and Analogues Potently Target DHFR in Cancer Cells to Elicit Anti-Cancer Activity

Dihydrofolate reductase (DHFR) is an established anti-cancer drug target whose inhibition disrupts folate metabolism and STAT3-dependent gene expression. Cycloguanil was proposed as a DHFR inhibitor in the 1950s and is the active metabolite of clinically approved plasmodium DHFR inhibitor Proguanil. The Cycloguanil scaffold was explored to generate potential cancer therapies in the 1970s. Herein, current computational and chemical biology techniques were employed to re-investigate the anti-cancer activity of Cycloguanil and related compounds. In silico modeling was employed to identify promising Cycloguanil analogues from NCI databases, which were cross-referenced with NCI-60 Human Tumor Cell Line Screening data. Using target engagement assays, it was found that these compounds engage DHFR in cells at sub-nanomolar concentrations; however, growth impairments were not observed until higher concentrations. Folinic acid treatment rescues the viability impairments induced by some, but not all, Cycloguanil analogues, suggesting these compounds may have additional targets. Cycloguanil and its most promising analogue, NSC127159, induced similar metabolite profiles compared to established DHFR inhibitors Methotrexate and Pyrimethamine while also blocking downstream signaling, including STAT3 transcriptional activity. These data confirm that Cycloguanil and its analogues are potent inhibitors of human DHFR, and their anti-cancer activity may be worth further investigation.

Towards establishment of a plant-based model to assess the novel anti-cancerous lead molecule(s): An in silico, in vivo and in vitro assessment of some potential anti-cancerous drugs on Lathyrus sativus L

The drug development process is one of the important aspects of medical biology. The classical lead identification strategy in the way of drug development based on animal cell is time-consuming, expensive and involving ethical issues. The following study aims to develop a novel plant-based screening of drugs. Study shows the efficacy of certain anti-cancerous drugs (Pemetrexed, 5-Fluorouracil, Methotrexate, Topotecan and Etoposide) on a plant-based (Lathyrus sativus L.) system. Two important characteristics of cancer cells were observed in the colchicine-treated polyploid cell and the callus, where the chromosome numbers were unusual and the division of cells were uncontrolled respectively. With increasing concentration, the drugs significantly reduced the mitotic index, ploidy level and callus growth. Increasing Pemetrexed concentration decreased the plant DHFR activity. A decrease in total RNA content was observed in 5-FU and Methotrexate with increasing concentrations of the drugs. Etoposide and Topotecan inhibited plant topoisomerase II and topoisomerase I activities, which was justified through plasmid nicking and comet assay, respectively. Molecular and biochemical study revealed similar results to the animal system. The in silico study had been done, and the structural similarity of drug binding domains of L. sativus and human beings had also been established. The binding site of the selected drugs to the domains of plant target proteins was also determined. Experimental results are significant in terms of the efficacy of known anti-cancerous drugs on the plant-based system. The proposed assay system is a cost-effective, convenient and less time-consuming process for primary screening of anti-cancerous lead molecules.

Screening of indole derivatives as the potent anticancer agents on dihydrofolate reductase: pharmaco-informatics and molecular dynamics simulation

Dihydrofolate reductase (DHFR) is a ubiquitous cellular enzyme involved in the biosynthesis of nucleotide and protein precursors, thus, the inhibition of human DHFR can be a promising strategy in cancer treatment. The design of effective anticancer drugs is an urgent need today according to the high spread of cancer. The indole molecule with diverse mechanisms of action and anticancer properties is one of the efficient pharmacophores in drug design. Hence, a virtual library of indole derivatives as a scaffold was selected for designing safer and more effective anticancer drugs against DHFR in this work. All indole derivatives utilized in the library design were selected regarding appreciable tumor growth inhibition. Structure-activity relationship (SAR), docking energy, ADMET (absorption, distribution, metabolism, excretion, and toxicity) parameters, and effective non-covalent interactions were used to identify potential anticancer with indole scaffold. Results showed a higher number of indole moieties provide a strong attachment to the DHFR binding pocket and therefore more effective anticancer activity. The indole scaffold in combination with dichlorobenzene improves DHFR inhibition whereas barbituric acid weakens inhibition activity. In the following to validate the docking results, Molecular dynamics (MD) simulation and molecular mechanics generalized-Born surface area (MM-GBSA) indicated the permanent stability of the selected ligands into the DHFR binding pocket and the key amino acids. Therefore, promising pharmacophores based on indole-DHFR interactions were discovered, and the outcome could be useful in guiding future in vitro and in vivo drug discovery in cancer medicine.Communicated by Ramaswamy H. Sarma.

Cytosolic localization and in vitro assembly of human de novo thymidylate synthesis complex

De novo thymidylate synthesis is a crucial pathway for normal and cancer cells. Deoxythymidine monophosphate (dTMP) is synthesized by the combined action of three enzymes: serine hydroxymethyltransferase (SHMT1), dihydrofolate reductase (DHFR) and thymidylate synthase (TYMS), with the latter two being targets of widely used chemotherapeutics such as antifolates and 5‐fluorouracil. These proteins translocate to the nucleus after SUMOylation and are suggested to assemble in this compartment into the thymidylate synthesis complex. We report the intracellular dynamics of the complex in cancer cells by an in situ proximity ligation assay, showing that it is also detected in the cytoplasm. This result indicates that the role of the thymidylate synthesis complex assembly may go beyond dTMP synthesis. We have successfully assembled the dTMP synthesis complex in vitro, employing tetrameric SHMT1 and a bifunctional chimeric enzyme comprising human thymidylate synthase and dihydrofolate reductase. We show that the SHMT1 tetrameric state is required for efficient complex assembly, indicating that this aggregation state is evolutionarily selected in eukaryotes to optimize protein–protein interactions. Lastly, our results regarding the activity of the complete thymidylate cycle in vitro may provide a useful tool with respect to developing drugs targeting the entire complex instead of the individual components.

Structural Insights of Three 2,4-Disubstituted Dihydropyrimidine-5-carbonitriles as Potential Dihydrofolate Reductase Inhibitors

In this report, we describe the structural characterization of three 2,4-disubstituted-dihydropyrimidine-5-carbonitrile derivatives, namely 2-{[(4-nitrophenyl)methyl]sulfanyl}-6-oxo-4-propyl-1,6-dihydropyrimidine-5-carbonitrile 1, 4-(2-methylpropyl)-2-{[(4-nitrophenyl)methyl]sulfanyl}-6-oxo-1,6-dihydropyrimidine-5-carbonitrile 2, and 2-[(2-ethoxyethyl)sulfanyl]-6-oxo-4-phenyl-1,6-dihydropyrimidine-5-carbonitrile monohydrate 3. An X-ray diffraction analysis revealed that these compounds were crystallized in the centrosymmetric space groups and adopt an L-shaped conformation. One of the compounds (3) crystallized with a water molecule. A cyclic motif (R₂²(8)) mediated by N–H···O hydrogen bond was formed in compounds 1 and 2, whereas the corresponding motif was not favorable, due to the water molecule, in compound 3. The crystal packing of these compounds was analyzed based on energy frameworks performed at the B3LYP/6-31G(d,p) level of theory. Various inter-contacts were characterized using the Hirshfeld surface and its associated 2D-fingerprint plots. Furthermore, a molecular docking simulation was carried out to assess the inhibitory potential of the title compounds against the human dihydrofolate reductase (DHFR) enzyme.

Distal Regions Regulate Dihydrofolate Reductase-Ligand Interactions

Protein motions play a fundamental role in enzyme catalysis and ligand binding. The relationship between protein motion and function has been extensively investigated in the model enzyme dihydrofolate reductase (DHFR). DHFR is an essential enzyme that catalyzes the reduction of dihydrofolate to tetrahydrofolate. Numerous experimental and computational methods have been used to probe the motions of DHFR through the catalytic cycle and to investigate the effect of distal mutations on DHFR motions and ligand binding. These experimental investigations have pushed forward the study of protein motions and their role in protein-ligand interactions. The introduction of mutations distal to the active site has been shown to have profound effects on ligand binding, hydride transfer rates and catalytic efficacy and these changes are captured by enzyme kinetics measurements. Distal mutations have been shown to exert their effects through a network of correlated amino acids and these effects have been investigated by NMR, protein dynamics, and analysis of coupled amino acids. The experimental methods and the findings that are reviewed here have broad implications for our understanding of enzyme mechanisms, ligand binding and for the future design and discovery of enzyme inhibitors.

Synthesis, Biological Activity, and Molecular Dynamics Study of Novel Series of a Trimethoprim Analogs as Multi-Targeted Compounds: Dihydrofolate Reductase (DHFR) Inhibitors and DNA-Binding Agents

Eighteen previously undescribed trimethoprim (TMP) analogs containing amide bonds (1-18) were synthesized and compared with TMP, methotrexate (MTX), and netropsin (NT). These compounds were designed as potential minor groove binding agents (MGBAs) and inhibitors of human dihydrofolate reductase (hDHFR). The all-new derivatives were obtained via solid phase synthesis using 4-nitrophenyl Wang resin. Data from the ethidium displacement test confirmed their DNA-binding capacity. Compounds 13-14 (49.89% and 43.85%) and 17-18 (41.68% and 42.99%) showed a higher binding affinity to pBR322 plasmid than NT. The possibility of binding in a minor groove as well as determination of association constants were performed using calf thymus DNA, T4 coliphage DNA, poly (dA-dT)2, and poly (dG-dC)2. With the exception of compounds 9 (IC50 = 56.05 µM) and 11 (IC50 = 55.32 µM), all of the compounds showed better inhibitory properties against hDHFR than standard, which confirms that the addition of the amide bond into the TMP structures increases affinity towards hDHFR. Derivatives 2, 6, 13, 14, and 16 were found to be the most potent hDHFR inhibitors. This molecular modelling study shows that they interact strongly with a catalytically important residue Glu-30.

In Silico Study Identified Methotrexate Analog as Potential Inhibitor of Drug Resistant Human Dihydrofolate Reductase for Cancer Therapeutics

Drug resistance is a core issue in cancer chemotherapy. A known folate antagonist, methotrexate (MTX) inhibits human dihydrofolate reductase (hDHFR), the enzyme responsible for the catalysis of 7,8-dihydrofolate reduction to 5,6,7,8-tetrahydrofolate, in biosynthesis and cell proliferation. Structural change in the DHFR enzyme is a significant cause of resistance and the subsequent loss of MTX. In the current study, wild type hDHFR and double mutant (engineered variant) F31R/Q35E (PDB ID: 3EIG) were subject to computational study. Structure-based pharmacophore modeling was carried out for wild type (WT) and mutant (MT) (variant F31R/Q35E) hDHFR structures by generating ten models for each. Two pharmacophore models, WT-pharma and MT-pharma, were selected for further computations, and showed excellent ROC curve quality. Additionally, the selected pharmacophore models were validated by the Guner-Henry decoy test method, which yielded high goodness of fit for WT-hDHFR and MT-hDHFR. Using a SMILES string of MTX in ZINC¹⁵ with the selections of 'clean', in vitro and in vivo options, 32 MTX-analogs were obtained. Eight analogs were filtered out due to their drug-like properties by applying absorption, distribution, metabolism, excretion, and toxicity (ADMET) assessment tests and Lipinski's Rule of five. WT-pharma and MT-pharma were further employed as a 3D query in virtual screening with drug-like MTX analogs. Subsequently, seven screening hits along with a reference compound (MTX) were subjected to molecular docking in the active site of WT- and MT-hDHFR. Through a clustering analysis and examination of protein-ligand interactions, one compound was found with a ChemPLP fitness score greater than that of MTX (reference compound). Finally, a simulation of molecular dynamics (MD) identified an MTX analog which exhibited strong affinity for WT- and MT-hDHFR, with stable RMSD, hydrogen bonds (H-bonds) in the binding site and the lowest MM/PBSA binding free energy. In conclusion, we report on an MTX analog which is capable of inhibiting hDHFR in wild type form, as well as in cases where the enzyme acquires resistance to drugs during chemotherapy treatment.

Detecting Protein-Protein Interaction Based on Protein Fragment Complementation Assay

Proteins are the most critical executive molecules by responding to the instructions stored in the genetic materials in any form of life. More frequently, proteins do their jobs by acting as a roleplayer that interacts with other protein(s), which is more evident when the function of a protein is examined in the real context of a cell. Identifying the interactions between (or amongst) proteins is very crucial for the biochemistry investigation of an individual protein and for the attempts aiming to draw a holo-picture for the interacting members at the scale of proteomics (or protein-protein interactions mapping). Here, we introduced the currently available reporting systems that can be used to probe the interaction between candidate protein pairs based on the fragment complementation of some particular proteins. Emphasis was put on the principles and details of experimental design. These systems are dihydrofolate reductase (DHFR), β-lactamase, tobacco etch virus (TEV) protease, luciferase, β- galactosidase, GAL4, horseradish peroxidase (HRP), focal adhesion kinase (FAK), green fluorescent protein (GFP), and ubiquitin.

Drugging the Folate Pathway in Mycobacterium tuberculosis: The Role of Multi-targeting Agents

The folate biosynthetic pathway offers many druggable targets that have yet to be exploited in tuberculosis therapy. Herein, we have identified a series of small molecules that interrupt Mycobacterium tuberculosis (Mtb) folate metabolism by dual targeting of dihydrofolate reductase (DHFR), a key enzyme in the folate pathway, and its functional analog, Rv2671. We have also compared the antifolate activity of these compounds with that of para-aminosalicylic acid (PAS). We found that the bioactive metabolite of PAS, in addition to previously reported activity against DHFR, inhibits flavin-dependent thymidylate synthase in Mtb, suggesting a multi-targeted mechanism of action for this drug. Finally, we have shown that antifolate treatment in Mtb decreases the production of mycolic acids, most likely due to perturbation of the activated methyl cycle. We conclude that multi-targeting of the folate pathway in Mtb is associated with highly potent anti-mycobacterial activity.

In Silico Study Probes Potential Inhibitors of Human Dihydrofolate Reductase for Cancer Therapeutics

Dihydrofolate reductase (DHFR) is an essential cellular enzyme and thereby catalyzes thereduction of dihydrofolate to tetrahydrofolate (THF). In cancer medication, inhibition of humanDHFR (hDHFR) remains a promising strategy, as it depletes THF and slows DNA synthesis and cellproliferation. In the current study, ligand-based pharmacophore modeling identified and evaluatedthe critical chemical features of hDHFR inhibitors. A pharmacophore model (Hypo1) was generatedfrom known inhibitors of DHFR with a correlation coefficient (0.94), root mean square (RMS)deviation (0.99), and total cost value (125.28). Hypo1 was comprised of four chemical features,including two hydrogen bond donors (HDB), one hydrogen bond acceptor (HBA), and onehydrophobic (HYP). Hypo1 was validated using Fischer's randomization, test set, and decoy setvalidations, employed as a 3D query in a virtual screening at Maybridge, Chembridge, Asinex,National Cancer Institute (NCI), and Zinc databases. Hypo1-retrieved compounds were filtered byan absorption, distribution, metabolism, excretion, and toxicity (ADMET) assessment test andLipinski's rule of five, where the drug-like hit compounds were identified. The hit compounds weredocked in the active site of hDHFR and compounds with Goldfitness score was greater than 44.67(docking score for the reference compound), clustering analysis, and hydrogen bond interactionswere identified. Furthermore, molecular dynamics (MD) simulation identified three compounds asthe best inhibitors of hDHFR with the lowest root mean square deviation (1.2 Å to 1.8 Å), hydrogenbond interactions with hDHFR, and low binding free energy (-127 kJ/mol to -178 kJ/mol). Finally,the toxicity prediction by computer (TOPKAT) affirmed the safety of the novel inhibitors of hDHFRin human body. Overall, we recommend novel hit compounds of hDHFR for cancer and rheumatoidarthritis chemotherapeutics.

Quantum mechanics implementation in drug-design workflows: does it really help?

The pharmaceutical industry is progressively operating in an era where development costs are constantly under pressure, higher percentages of drugs are demanded, and the drug-discovery process is a trial-and-error run. The profit that flows in with the discovery of new drugs has always been the motivation for the industry to keep up the pace and keep abreast with the endless demand for medicines. The process of finding a molecule that binds to the target protein using in silico tools has made computational chemistry a valuable tool in drug discovery in both academic research and pharmaceutical industry. However, the complexity of many protein-ligand interactions challenges the accuracy and efficiency of the commonly used empirical methods. The usefulness of quantum mechanics (QM) in drug-protein interaction cannot be overemphasized; however, this approach has little significance in some empirical methods. In this review, we discuss recent developments in, and application of, QM to medically relevant biomolecules. We critically discuss the different types of QM-based methods and their proposed application to incorporating them into drug-design and -discovery workflows while trying to answer a critical question: are QM-based methods of real help in drug-design and -discovery research and industry?

Pharmacoinformatic Study on the Selective Inhibition of the Protozoan Dihydrofolate Reductase Enzymes

Dihydrofolate reductase (DHFR) is an essential enzyme of the folate metabolic pathway in protozoa and it is a validated, potential drug target in many infectious diseases. Information about unique conserved residues of the DHFR enzyme is required to understand residual selectivity of the protozoan DHFR enzyme. The three dimensional crystal structures are not available for all the protozoan DHFR enzymes. Enzyme-substrate/inhibitor interaction information is required for the binding mode characterization in protozoan DHFR for selective inhibitor design. In this work, multiple sequence analysis was carried out in all the studied species. Homology models were built for protozoan DHFR enzymes, for which 3D structures are not available in PDB. The molecular docking and Prime-MMGBSA calculations of the natural substrate (dihydrofolate, DHF) and classical DHFR inhibitor (methotrexate, MTX) were performed in protozoan DHFR enzymes. Comparative sequence analysis showed that an overall sequence identity between the studied species ranging from 22.94 % (CfDHFR-BgDHFR) to 94.61 % (LdDHFR-LmDHFR). Interestingly, it was observed that most of the active site residues were conserved in all the cases and all the enzymes exhibit similar key binding interactions with DHF and MTX in molecular docking analysis, but there are a few key binding residues which differ in protozoan species that makes it suitable for target selectivity. This information can be used to design selective and potent protozoan DHFR enzyme inhibitors.

A comparative molecular docking study of curcumin and methotrexate to dihydrofolate reductase

Interaction of curcumin (CUR) with the enzyme dihydrofolate reductase (DHFR) was studied by molecular docking using AutoDock 4.2 as the docking software application. AutoDock 4.2 software serves as a valid and acceptable docking application to study the interactions of small compounds with proteins. Interactions of curcumin with DHFR were compared to those of methotrexate (MTX), a known inhibitor of the enzyme. The calculated free energy of binding (ΔG binding) shows that curcumin (ΔG = -9.02 kcal/mol; Ki = 243 nM) binds with affinity comparable to or better than MTX (ΔG = -8.78 kcal/mol; Ki = 363 nM). Binding interactions of curcumin with active site residues of the enzyme are also predicted. Curcumin appears to bind in a bent conformation making extensive VDW contacts in the active site of the enzyme. Hydrogen bonding and pi-pi interaction with key active site residues are also observed. Thus, curcumin can be considered as a good lead compound in the development of new inhibitors of DHFR, which is a potential target of anti-cancer drugs. The results of these studies can serve as a starting point for further computational and experimental studies.

Classification of ligand molecules in PDB with graph match-based structural superposition

The fast heuristic graph match algorithm for small molecules, COMPLIG, was improved by adding a structural superposition process to verify the atom-atom matching. The modified method was used to classify the small molecule ligands in the Protein Data Bank (PDB) by their three-dimensional structures, and 16,660 types of ligands in the PDB were classified into 7561 clusters. In contrast, a classification by a previous method (without structure superposition) generated 3371 clusters from the same ligand set. The characteristic feature in the current classification system is the increased number of singleton clusters, which contained only one ligand molecule in a cluster. Inspections of the singletons in the current classification system but not in the previous one implied that the major factors for the isolation were differences in chirality, cyclic conformations, separation of substructures, and bond length. Comparisons between current and previous classification systems revealed that the superposition-based classification was effective in clustering functionally related ligands, such as drugs targeted to specific biological processes, owing to the strictness of the atom-atom matching.

Arrested Hematopoiesis and Vascular Relaxation Defects in Mice with a Mutation in Dhfr

Dihydrofolate reductase (DHFR) is a critical enzyme in the folate metabolism pathway and also plays a role in regulating nitric oxide (NO) signaling in endothelial cells. Although both coding and noncoding mutations with phenotypic effects have been identified in the human DHFR gene, no mouse model is currently available to study the consequences of perturbing DHFR in vivo In order to identify genes involved in definitive hematopoiesis, we performed a forward genetic screen and produced a mouse line, here referred to as Orana, with a point mutation in the Dhfr locus leading to a Thr136Ala substitution in the DHFR protein. Homozygote Orana mice initiate definitive hematopoiesis, but expansion of progenitors in the fetal liver is compromised, and the animals die between embryonic day 13.5 (E13.5) and E14.5. Heterozygote Orana mice survive to adulthood but have tissue-specific alterations in folate abundance and distribution, perturbed stress erythropoiesis, and impaired endothelium-dependent relaxation of the aorta consistent with the role of DHFR in regulating NO signaling. Orana mice provide insight into the dual roles of DHFR and are a useful model for investigating the role of environmental and dietary factors in the context of vascular defects caused by altered NO signaling.

Structural analysis of dihydrofolate reductases enables rationalization of antifolate binding affinities and suggests repurposing possibilities

Antifolates are competitive inhibitors of dihydrofolate reductase (DHFR), a conserved enzyme that is central to metabolism and widely targeted in pathogenic diseases, cancer and autoimmune disorders. Although most clinically used antifolates are known to be target specific, some display a fair degree of cross-reactivity with DHFRs from other species. A method that enables identification of determinants of affinity and specificity in target DHFRs from different species and provides guidelines for the design of antifolates is currently lacking. To address this, we first captured the potential druggable space of a DHFR in a substructure called the 'supersite' and classified supersites of DHFRs from 56 species into 16 'site-types' based on pairwise structural similarity. Analysis of supersites across these site-types revealed that DHFRs exhibit varying extents of dissimilarity at structurally equivalent positions in and around the binding site. We were able to explain the pattern of affinities towards chemically diverse antifolates exhibited by DHFRs of different site-types based on these structural differences. We then generated an antifolate-DHFR network by mapping known high-affinity antifolates to their respective supersites and used this to identify antifolates that can be repurposed based on similarity between supersites or antifolates. Thus, we identified 177 human-specific and 458 pathogen-specific antifolates, a large number of which are supported by available experimental data. Thus, in the light of the clinical importance of DHFR, we present a novel approach to identifying differences in the druggable space of DHFRs that can be utilized for rational design of antifolates.

Substrate Activity Screening (SAS) and Related Approaches in Medicinal Chemistry

Substrate activity screening (SAS) was presented a decade ago by Ellman and co-workers as a straightforward methodology for the identification of fragment-sized building blocks for enzyme inhibitors. Ever since, SAS and variations derived from it have been successfully applied to the discovery of inhibitors of various families of enzymatically active drug targets. This review covers key achievements and challenges of SAS and related methodologies, including the modified substrate activity screening (MSAS) approach. Special attention is given to the kinetic and thermodynamic aspects of these methodologies, as a thorough understanding thereof is crucial for successfully transforming the identified fragment-sized hits into potent inhibitors.

Structural genomics for drug design against the pathogen Coxiella burnetii

Coxiella burnetii is a highly infectious bacterium and potential agent of bioterrorism. However, it has not been studied as extensively as other biological agents, and very few of its proteins have been structurally characterized. To address this situation, we undertook a study of critical metabolic enzymes in C. burnetii that have great potential as drug targets. We used high-throughput techniques to produce novel crystal structures of 48 of these proteins. We selected one protein, C. burnetii dihydrofolate reductase (CbDHFR), for additional work to demonstrate the value of these structures for structure-based drug design. This enzyme's structure reveals a feature in the substrate binding groove that is different between CbDHFR and human dihydrofolate reductase (hDHFR). We then identified a compound by in silico screening that exploits this binding groove difference, and demonstrated that this compound inhibits CbDHFR with at least 25-fold greater potency than hDHFR. Since this binding groove feature is shared by many other prokaryotes, the compound identified could form the basis of a novel antibacterial agent effective against a broad spectrum of pathogenic bacteria.

Environmental Adaptation of Dihydrofolate Reductase from Deep-Sea Bacteria

In order to elucidate the molecular adaptation mechanisms of enzymes to the high hydrostatic pressure of the deep sea, we cloned, purified, and characterized more than ten dihydrofolate reductases (DHFRs) from bacteria living in deep-sea and ambient atmospheric pressure environments. The nucleotide and amino acid sequences of these DHFRs indicate the deep-sea bacteria are adapted to their environments after the differentiation of their genus from ancestors inhabiting atmospheric pressure environments. In particular, the backbone structure of the deep-sea DHFR from Moritella profunda (mpDHFR) almost overlapped with the normal homolog from Escherichia coli (ecDHFR). Thus, those of other DHFRs would also overlap on the basis of their sequence similarities. However, the structural stability of both DHFRs was quite different: compared to ecDHFR, mpDHFR was more thermally stable but less stable against urea and pressure unfolding. The smaller volume changes due to unfolding suggest that the native structure of mpDHFR has a smaller cavity and/or enhanced hydration compared to ecDHFR. High hydrostatic pressure reduced the enzymatic activity of many DHFRs, but three deep-sea DHFRs and the D27E mutant of ecDHFR exhibited pressure-dependent activation. The inverted activation volumes from positive to negative values indicate the modification of their structural dynamics, conversion of the rate-determining step of the enzymatic reaction, and different contributions of the cavity and hydration to the transition-state structure. Since the cavity and hydration depend on amino acid side chains, DHFRs would adapt to the deep-sea environment by regulating the cavity and hydration by substituting their amino acid side chains without altering their backbone structure. The results of this study clearly indicate that the cavity and hydration play important roles in the adaptation of enzymes to the deep-sea environment.

Study on Folate Binding Domain of Dihydrofolate Reductase in Different Plant species and Human beings

Data base (NCBI and TIGR) searches are made to retrieve protein sequences of different plant species namely Medicago truncatula, Pisum sativum, Ricinus communis, Arabidopsis thaliana, Vitis vinifera, Glycine max, Daucus carota, Oryza sativa Japonica Group, Arabidopsis lyrata subsp. lyrata, Brachypodium distachyon, Oryza sativa Indica Group, Zea mays and careful alignment of derived sequences shows 95% or higher identity. Similarly, DHFR sequence of human being is also retrieved from NCBI. A phylogenetic tree is constructed from different plant and human DHFR domain using the Neighbour - Joining method in MEGA 5.05. Conservation score is performed by using PARALINE. Result suggests that folate binding domain of dihydrofolare reductase is conserved (score 8.06) and excepting some minor variations the basic structure of the domain in both plant species and human being is rather similar. Human DHFR domain contains PEKN sequence near active site, though proline is common for all the selected organisms but the other sequences are different in plants. The plant domain is always associated with TS (Thymidylate synthase). Plant based system is predicted to be an effective model for assessment of MTX (Methotrexate) and other antifolate drugs.

Structural analysis of a fungal methionine synthase with substrates and inhibitors

The cobalamin-independent methionine synthase from Candida albicans, known as Met6p, is a 90-kDa enzyme that consists of two (βα)8 barrels. The active site is located between the two domains and has binding sites for a zinc ion and substrates L-homocysteine and 5-methyl-tetrahydrofolate-glutamate3. Met6p catalyzes transfer of the methyl group of 5-methyl-tetrahydrofolate-glutamate3 to the L-homocysteine thiolate to generate methionine. Met6p is essential for fungal growth, and we currently pursue it as an antifungal drug design target. Here we report the binding of L-homocysteine, methionine, and several folate analogs. We show that binding of L-homocysteine or methionine results in conformational rearrangements at the amino acid binding pocket, moving the catalytic zinc into position to activate the thiol group. We also map the folate binding pocket and identify specific binding residues, like Asn126, whose mutation eliminates catalytic activity. We also report the development of a robust fluorescence-based activity assay suitable for high-throughput screening. We use this assay and an X-ray structure to characterize methotrexate as a weak inhibitor of fungal Met6p.

Insight into the molecular mechanism about lowered dihydrofolate binding affinity to dihydrofolate reductase-like 1 (DHFRL1)

Human dihydrofolate reductase-like 1 (DHFRL1) has been identified as a second human dihydrofolate reductase (DHFR) enzyme. Although DHFRL1 have high sequence homology with human DHFR, dihydrofolate (DHF) exhibits a lowered binding affinity to DHFRL1 and the corresponding molecular mechanism is still unknown. To address this question, we studied the binding of DHF to DHFRL1 and DHFR by using molecular dynamics simulation. Moreover, to investigate the role the 24th residue of DHFR/DHFRL1 plays in DHF binding, R24W DHFRL1 mutant was also studied. The van der Waals interaction are more crucial for the total DHF binding energies, while the difference between the DHF binding energies of human DHFR and DHFRL1 can be attributed to the electrostatic interaction and the polar desolvation free energy.More specifically, lower DHF affinity to DHFRL1 can be mainly attributed to the reduction of net electrostatic interactions of residues Arg32 and Gln35 of DHFRL1 with DHF as being affected by Arg24. The side chain of Arg24 in DHFRL1 can extend deeply into the binding sites of DHF and NADPH, and disturb the DHF binding by steric effect, which rarely happens in human DHFR and R24W DHFRL1 mutant. Additionally, the conformation of loop I in DHFRL1 was also studied in this work. Interestingly, the loop conformation resemble to normal closed state of Escherichia coli DHFR other than the closed state of human DHFR. We hope this work will be useful to understand the general characteristics of DHFRL1.

Molecular modeling study of dihydrofolate reductase inhibitors. Molecular dynamics simulations, quantum mechanical calculations, and experimental corroboration

A molecular modeling study on dihydrofolate reductase (DHFR) inhibitors was carried out. By combining molecular dynamics simulations with semiempirical (PM6), ab initio, and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of DHFR inhibitors interacting with the human enzyme is reported here, providing a clear picture of the binding interactions of these ligands from both structural and energetic viewpoints. A reduced model for the binding pocket was used. This approach allows us to perform more accurate quantum mechanical calculations as well as to obtain a detailed electronic analysis using the quantum theory of atoms in molecules (QTAIM) technique. Thus, molecular aspects of the binding interactions between inhibitors and the DHFR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental IC₅₀ values was obtained, predicting with an acceptable qualitative accuracy the potential inhibitor effect of nonsynthesized compounds. Such correlation was experimentally corroborated synthesizing and testing two new inhibitors reported in this paper.

Biological and structural evaluation of 10R- and 10S-methylthio-DDACTHF reveals a new role for sulfur in inhibition of glycinamide ribonucleotide transformylase

Glycinamide ribonucleotide transformylase (GAR Tfase) is a folate-dependent enzyme in the de novo purine biosynthesis pathway, which has long been considered a potential target for development of anti-neoplastic therapeutics. Here we report the biological and X-ray crystallographic evaluations of both independent C10 diastereomers, 10S- and 10R-methylthio-DDACTHF, bound to human GAR Tfase, including the highest-resolution apo GAR Tfase structure to date (1.52 Å). Both diastereomers are potent inhibitors (Ki = 210 nM for 10R, and Ki = 180 nM for 10S) of GAR Tfase and exhibit effective inhibition of human leukemia cell growth (IC₅₀ = 80 and 50 nM, respectively). Their inhibitory activity was surprisingly high, and these lipophilic C10-substituted analogues show distinct advantages over their hydrophilic counterparts, most strikingly in retaining potency in mutant human leukemia cell lines that lack reduced folate carrier protein activity (IC₅₀ = 70 and 60 nM, respectively). Structural characterization reveals a new binding mode for these diastereoisomers, in which the lipophilic thiomethyl groups penetrate deeper into a hydrophobic pocket within the folate-binding site. In silico docking simulations of three other sulfur-containing folate analogues also indicate that this hydrophobic cleft represents a favorable region for binding lipophilic substituents. Overall, these results suggest sulfur and its substitutions play an important role in not only the binding of anti-folates to GAR Tfase but also the selectivity and cellular activity (growth inhibition), thereby presenting new possibilities for the future design of potent and selective anti-folate drugs that target GAR Tfase.

Nonclassical antifolates, part 4. 5-(2-aminothiazol-4-yl)-4-phenyl-4H-1,2,4-triazole-3-thiols as a new class of DHFR inhibitors: synthesis, biological evaluation and molecular modeling study

A new series of compounds possessing 5-(2-aminothiazol-4-yl)-4-phenyl-4H-1,2,4-triazole-3-thiol skeleton was designed, synthesized, and evaluated for their in vitro DHFR inhibition, antimicrobial, antitumor and schistosomicidal activities. Four active compounds were allocated, the antibacterial 22 (comparable to gentamicin and ciprofloxacin), the schistosomicidal 29 (comparable to praziquantel), the DHFR inhibitor 34 (IC₅₀ 0.03 μM, 2.7 fold more active than MTX), and the antitumor 36 (comparable to doxorubicin). Molecular modeling studies concluded that recognition with key amino acid Leu4 and Val1 is essential for DHFR binding. Flexible alignment and surface mapping revealed that the obtained model could be useful for the development of new class of DHFR inhibitors.

Preferential selection of isomer binding from chiral mixtures: alternate binding modes observed for the E and Z isomers of a series of 5-substituted 2,4-diaminofuro[2,3-d]pyrimidines as ternary complexes with NADPH and human dihydrofolate reductase

The crystal structures of six human dihydrofolate reductase (hDHFR) ternary complexes with NADPH and a series of mixed E/Z isomers of 5-substituted 5-[2-(2-methoxyphenyl)-prop-1-en-1-yl]furo[2,3-d]pyrimidine-2,4-diamines substituted at the C9 position with propyl, isopropyl, cyclopropyl, butyl, isobutyl and sec-butyl (E2-E7, Z3) were determined and the results were compared with the resolved E and Z isomers of the C9-methyl parent compound. The configuration of all of the inhibitors, save one, was observed as the E isomer, in which the binding of the furopyrimidine ring is flipped such that the 4-amino group binds in the 4-oxo site of folate. The Z3 isomer of the C9-isopropyl analog has the normal 2,4-diaminopyrimidine ring binding geometry, with the furo oxygen near Glu30 and the 4-amino group interacting near the cofactor nicotinamide ring. Electron-density maps for these structures revealed the binding of only one isomer to hDHFR, despite the fact that chiral mixtures (E:Z ratios of 2:1, 3:1 and 3:2) of the inhibitors were incubated with hDHFR prior to crystallization. Superposition of the hDHFR complexes with E2 and Z3 shows that the 2'-methoxyphenyl ring of E2 is perpendicular to that of Z3. The most potent inhibitor in this series is the isopropyl analog Z3 and the least potent is the isobutyl analog E6, consistent with data that show that the Z isomer makes the most favorable interactions with the active-site residues. The isobutyl moiety of E6 is observed in two orientations and the resultant steric crowding of the E6 analog is consistent with its weaker activity. The alternative binding modes observed for the furopyrimidine ring in these E/Z isomers suggest that new templates can be designed to probe these binding regions of the DHFR active site.

Single agents with designed combination chemotherapy potential: synthesis and evaluation of substituted pyrimido[4,5-b]indoles as receptor tyrosine kinase and thymidylate synthase inhibitors and as antitumor agents

Combinations of antiangiogenic agents (AAs) with cytotoxic agents have shown significant promise in cancer treatment, and several such clinical trials are currently underway. We have designed, synthesized, and evaluated two compounds that each inhibit vascular endothelial growth factor receptor-2 (VEGFR-2) and platelet-derived growth factor receptor-beta (PDGFR-beta) for antiangiogenic effects and also inhibit human thymidylate synthase (hTS) for cytotoxic effects in single agents. The synthesis of these compounds involved the nucleophilic displacement of the common intermediate 5-chloro-9H-pyrimido[4,5-b]indole-2,4-diamine with appropriate benzenethiols. The inhibitory potency of both these single agents against VEGFR-2, PDGFR-beta, and hTS is better than or close to standards. In a COLO-205 xenograft mouse model, one of the analogs significantly decreased tumor growth (tumor growth inhibition (TGI) = 76% at 35 mg/kg), liver metastases, and tumor blood vessels compared with a standard drug and with control and thus demonstrated potent tumor growth inhibition, inhibition of metastasis, and antiangiogenic effects in vivo. These compounds afford combination chemotherapeutic potential in single agents.

Interpretable correlation descriptors for quantitative structure-activity relationships

BACKGROUND

The topological maximum cross correlation (TMACC) descriptors are alignment-independent 2D descriptors for the derivation of QSARs. TMACC descriptors are generated using atomic properties determined by molecular topology. Previous validation (J Chem Inf Model 2007, 47: 626-634) of the TMACC descriptor suggests it is competitive with the current state of the art.

RESULTS

Here, we illustrate the interpretability of the TMACC descriptors, through the analysis of the QSARs of inhibitors of angiotensin converting enzyme (ACE) and dihydrofolate reductase (DHFR). In the case of the ACE inhibitors, the TMACC interpretation shows features specific to C-domain inhibition, which have not been explicitly identified in previous QSAR studies.

CONCLUSIONS

The TMACC interpretation can provide new insight into the structure-activity relationships studied. Freely available, open source software for generating the TMACC descriptors can be downloaded from http://comp.chem.nottingham.ac.uk.

Functional role for the conformationally mobile phenylalanine 223 in the reaction of methylenetetrahydrofolate reductase from Escherichia coli

The flavoprotein methylenetetrahydrofolate reductase from Escherichia coli catalyzes the reduction of 5,10-methylenetetrahydrofolate (CH(2)-H(4)folate) by NADH via a ping-pong reaction mechanism. Structures of the reduced enzyme in complex with NADH and of the oxidized Glu28Gln enzyme in complex with CH(3)-H(4)folate [Pejchal, R., Sargeant, R., and Ludwig, M. L. (2005) Biochemistry 44, 11447-11457] have revealed Phe223 as a conformationally mobile active site residue. In the NADH complex, the NADH adopts an unusual hairpin conformation and is wedged between the isoalloxazine ring of the FAD and the side chain of Phe223. In the folate complex, Phe223 swings out from its position in the NADH complex to stack against the p-aminobenzoate ring of the folate. Although Phe223 contacts each substrate in E. coli MTHFR, this residue is not invariant; for example, a leucine occurs at this site in the human enzyme. To examine the role of Phe223 in substrate binding and catalysis, we have constructed mutants Phe223Ala and Phe223Leu. As predicted, our results indicate that Phe223 participates in the binding of both substrates. The Phe223Ala mutation impairs NADH and CH(2)-H(4)folate binding each 40-fold yet slows catalysis of both half-reactions less than 2-fold. Affinity for CH(2)-H(4)folate is unaffected by the Phe223Leu mutation, and the variant catalyzes the oxidative half-reaction 3-fold faster than the wild-type enzyme. Structures of ligand-free Phe223Leu and Phe223Leu/Glu28Gln MTHFR in complex with CH(3)-H(4)folate have been determined at 1.65 and 1.70 A resolution, respectively. The structures show that the folate is bound in a catalytically competent conformation, and Leu223 undergoes a conformational change similar to that observed for Phe223 in the Glu28Gln-CH(3)-H(4)folate structure. Taken together, our results suggest that Leu may be a suitable replacement for Phe223 in the oxidative half-reaction of E. coli MTHFR.

Multiple conformers in active site of human dihydrofolate reductase F31R/Q35E double mutant suggest structural basis for methotrexate resistance

Methotrexate is a slow, tight-binding, competitive inhibitor of human dihydrofolate reductase (hDHFR), an enzyme that provides key metabolites for nucleotide biosynthesis. In an effort to better characterize ligand binding in drug resistance, we have previously engineered hDHFR variant F31R/Q35E. This variant displays a >650-fold decrease in methotrexate affinity, while maintaining catalytic activity comparable to the native enzyme. To elucidate the molecular basis of decreased methotrexate affinity in the doubly substituted variant, we determined kinetic and inhibitory parameters for the simple variants F31R and Q35E. This demonstrated that the important decrease of methotrexate affinity in variant F31R/Q35E is a result of synergistic effects of the combined substitutions. To better understand the structural cause of this synergy, we obtained the crystal structure of hDHFR variant F31R/Q35E complexed with methotrexate at 1.7-A resolution. The mutated residue Arg-31 was observed in multiple conformers. In addition, seven native active-site residues were observed in more than one conformation, which is not characteristic of the wild-type enzyme. This suggests that increased residue disorder underlies the observed methotrexate resistance. We observe a considerable loss of van der Waals and polar contacts with the p-aminobenzoic acid and glutamate moieties. The multiple conformers of Arg-31 further suggest that the amino acid substitutions may decrease the isomerization step required for tight binding of methotrexate. Molecular docking with folate corroborates this hypothesis.

Mutational 'hot-spots' in mammalian, bacterial and protozoal dihydrofolate reductases associated with antifolate resistance: sequence and structural comparison

Human dihydrofolate reductase (DHFR) is a primary target for antifolate drugs in cancer treatment, while DHFRs from Plasmodium falciparum, Plasmodium vivax and various bacterial species are primary targets in the treatment of malaria and bacterial infections. Mutations in each of these DHFRs can result in resistance towards clinically relevant antifolates. We review the structural and functional impact of active-site mutations with respect to enzyme activity and antifolate resistance of DHFRs from mammals, protozoa and bacteria. The high structural homology between DHFRs results in a number of cross-species, active-site 'hot-spots' for broad-based antifolate resistance. In addition, we identify mutations that confer species-specific resistance, or antifolate-specific resistance. This comparative review of antifolate binding in diverse species provides new insights into the relationship between antifolate design and the development of mutational resistance. It also presents avenues for designing antifolate-resistant mammalian DHFRs as chemoprotective agents.

The sulfur atoms of the substrate CoA and the catalytic cysteine are required for a productive mode of substrate binding in bacterial biosynthetic thiolase, a thioester-dependent enzyme

Thioesters are more reactive than oxoesters, and thioester chemistry is important for the reaction mechanisms of many enzymes, including the members of the thiolase superfamily, which play roles in both degradative and biosynthetic pathways. In the reaction mechanism of the biosynthetic thiolase, the thioester moieties of acetyl-CoA and the acetylated catalytic cysteine react with each other, forming the product acetoacetyl-CoA. Although a number of studies have been carried out to elucidate the thiolase reaction mechanism at the atomic level, relatively little is known about the factors determining the affinity of thiolases towards their substrates. We have carried out crystallographic studies on the biosynthetic thiolase from Zoogloea ramigera complexed with CoA and three of its synthetic analogues to compare the binding modes of these related compounds. The results show that both the CoA terminal SH group and the side chain SH group of the catalytic Cys89 are crucial for the correct positioning of substrate in the thiolase catalytic pocket. Furthermore, calorimetric assays indicate that the mutation of Cys89 into an alanine significantly decreases the affinity of thiolase towards CoA. Thus, although the sulfur atom of the thioester moiety is important for the reaction mechanism of thioester-dependent enzymes, its specific properties can also affect the affinity and competent mode of binding of the thioester substrates to these enzymes.

Folding trajectories of human dihydrofolate reductase inside the GroEL GroES chaperonin cavity and free in solution

The chaperonin GroEL binds non-native polypeptides in an open ring via hydrophobic contacts and then, after ATP and GroES binding to the same ring as polypeptide, mediates productive folding in the now hydrophilic, encapsulated cis chamber. The nature of the folding reaction in the cis cavity remains poorly understood. In particular, it is unclear whether polypeptides take the same route to the native state in this cavity as they do when folding spontaneously free in solution. Here, we have addressed this question by using NMR measurements of the time course of acquisition of amide proton exchange protection of human dihydrofolate reductase (DHFR) during folding in the presence of methotrexate and ATP either free in solution or inside the stable cavity formed between a single ring variant of GroEL, SR1, and GroES. Recovery of DHFR refolded by the SR1/GroES-mediated reaction is 2-fold higher than in the spontaneous reaction. Nevertheless, DHFR folding was found to proceed by the same trajectories inside the cis folding chamber and free in solution. These observations are consistent with the description of the chaperonin chamber as an "Anfinsen cage" where polypeptide folding is determined solely by the amino acid sequence, as it is in solution. However, if misfolding occurs in the confinement of the chaperonin cavity, the polypeptide chain cannot undergo aggregation but rather finds its way back to a productive pathway in a manner that cannot be accomplished in solution, resulting in the observed high overall recovery.

Increasing methotrexate resistance by combination of active-site mutations in human dihydrofolate reductase

Methotrexate-resistant forms of human dihydrofolate reductase have the potential to protect healthy cells from the toxicity of methotrexate (MTX), to improve prognosis during cancer therapy. It has been shown that synergistic MTX-resistance can be obtained by combining two active-site mutations that independently confer weak MTX-resistance. In order to obtain more highly MTX-resistant human dihydrofolate reductase (hDHFR) variants for this application, we used a semi-rational approach to obtain combinatorial active-site mutants of hDHFR that are highly resistant towards MTX. We created a combinatorial mutant library encoding various amino acids at residues Phe31, Phe34 and Gln35. In vivo library selection was achieved in a bacterial system on medium containing high concentrations of MTX. We characterized ten novel MTX-resistant mutants with different amino acid combinations at residues 31, 34 and 35. Kinetic and inhibition parameters of the purified mutants revealed that higher MTX-resistance roughly correlated with a greater number of mutations, the most highly-resistant mutants containing three active site mutations (Ki(MTX)=59-180 nM; wild-type Ki(MTX)<0.03 nM). An inverse correlation was observed between resistance and catalytic efficiency, which decreased mostly as a result of increased KM toward the substrate dihydrofolate. We verified that the MTX-resistant hDHFRs can protect eukaryotic cells from MTX toxicity by transfecting the most resistant mutants into DHFR-knock-out CHO cells. The transfected variants conferred survival at concentrations of MTX between 100-fold and >4000-fold higher than the wild-type enzyme, the most resistant triple mutant offering protection beyond the maximal concentration of MTX that could be included in the medium. These highly resistant variants of hDHFR offer potential for myeloprotection during administration of MTX in cancer treatment.