2-tier bacterial and in vitro selection of active and methotrexate-resistant variants of human dihydrofolate reductase

Preparation, characterization and in vitro antitumor activity evaluation of hyaluronic acid-alendronate-methotrexate nanoparticles

As an anti-metabolic drug, methotrexate (MTX) plays an important role in the treatment of various malignant tumors. However, several side effects such as low selectivity and high toxic of MTX limited its further applications. With aims to increase its accumulation in the tumor sites and reduce the toxicity of normal tissue nonspecific uptake, a self-assembled hyaluronic acid-alendronate-methotrexate nanoparticle (HA-ALN-MTX NPs) with a dual-tumor-targeted drug loaded system was designed and synthesized with an average particle size of 265.6 ± 13.3 nm. The advantage of this nanosystem is that the anticancer drug MTX can be used as a tumor-targeted ligand for folate acid receptors (FA), and hyaluronic acid (HA) can be used as another tumor targeted ligand for CD44 receptors. In vitro experiments confirmed that HA-ALN-MTX NPs has lower toxic effect on normal tissue cells HUVECs and has relatively high proliferation inhibition effect on tumor cells A549. Moreover, the inhibition effect could be adjusted by altering the dose of given drugs. All these results revealed that the prepared HA-ALN-MTX NPs could be selectively taken up by tumor cells by FA and CD44 receptor-mediated endocytosis. Therefore, self-assembled HA-ALN-MTX NPs targeted by these FA/CD44 receptors for anticancer drugs could act as effective antitumor drugs.

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.

Investigation of Classical Organic and Ionic Liquid Cosolvents for Early-Stage Screening in Fragment-Based Inhibitor Design with Unrelated Bacterial and Human Dihydrofolate Reductases

Drug design by methods such as fragment screening requires effective solubilization of millimolar concentrations of small organic compounds while maintaining the properties of the biological target. We investigate four organic solvents and three 1-butyl-3-methylimidazolium (BMIm)-based ionic liquids (ILs) as cosolvents to establish conditions for screening two structurally unrelated dihydrofolate reductases (DHFRs) that are prime drug targets. Moderate concentrations (10%-15%) of cosolvents had little effect on inhibition of the microbial type II R67 DHFR and of human DHFR (hDHFR), while higher concentrations of organic cosolvents generally decreased activity of both DHFRs. In contrast, a specific IL conserved the activity of one DHFR, while severely reducing the activity of the other, and vice versa, illustrating the differing effect of ILs on distinct protein folds. Most of the cosolvents investigated preserved the fold of R67 DHFR and had little effect on binding of the cofactor NADPH, but reduced the productive affinity for its substrate. In contrast, cosolvents resulted in modest structural destabilization of hDHFR with little effect on productive affinity. We conclude that the organic cosolvents, methanol, dimethylformamide, and dimethylsulfoxide, offer the most balanced conditions for early-stage compound screening as they maintain sufficient biological activity of both DHFRs while allowing for compound dissolution in the millimolar range. However, IL cosolvents showed poor capacity to solubilize organic compounds at millimolar concentrations, mitigating their utility in early-stage screening. Nonetheless, ILs could provide an alternative to classical organic cosolvents when low concentrations of inhibitors are used, as when characterizing higher affinity inhibitors.

Reverse Chemical Genetics: Comprehensive Fitness Profiling Reveals the Spectrum of Drug Target Interactions

The emergence and prevalence of drug resistance demands streamlined strategies to identify drug resistant variants in a fast, systematic and cost-effective way. Methods commonly used to understand and predict drug resistance rely on limited clinical studies from patients who are refractory to drugs or on laborious evolution experiments with poor coverage of the gene variants. Here, we report an integrative functional variomics methodology combining deep sequencing and a Bayesian statistical model to provide a comprehensive list of drug resistance alleles from complex variant populations. Dihydrofolate reductase, the target of methotrexate chemotherapy drug, was used as a model to identify functional mutant alleles correlated with methotrexate resistance. This systematic approach identified previously reported resistance mutations, as well as novel point mutations that were validated in vivo. Use of this systematic strategy as a routine diagnostics tool widens the scope of successful drug research and development.

One of the most profound outcomes of fast, reliable genome sequencing is the ability to tailor drug therapy to an individual’s genotype. This ‘personalized’ or ‘precision medicine’ is the realization of a decades-long effort to maximize drug effect and limit unwanted side effects. An undesirable consequence of such targeted therapies, however, is the emergence of drug resistance. This outcome is the result of an evolutionary process where mutations in the drug target render the drug perturbation allow such mutant cells to proliferate. Because of the unbiased, and stochastic nature of the emergence of drug resistance, it is impossible to predict. We developed a test where hundreds of thousands of mutant cells are exposed to a drug simultaneously and those cells that modulate resistance survive. This method is innovative because it partners a high-throughput experimental protocol with a tailored statistical model to identify all mutations that modulate resistance. Finally, we used synthetic biology to re-create these mutations and demonstrate that they were, in fact, bona fide drug-resistant variants. These mutations were further extended and confirmed to also be resistant in the human orthologue. This combined biological-computational approach allows one to identify drug’s degree of resistance to both guide treatments and future drug discovery.

Engineering human T cells for resistance to methotrexate and mycophenolate mofetil as an in vivo cell selection strategy

Gene transfer and drug selection systems that enforce ongoing transgene expression in vitro and in vivo which are compatible with human pharmaceutical drugs are currently underdeveloped. Here, we report on the utility of incorporating human enzyme muteins that confer resistance to the lymphotoxic/immunosuppressive drugs methotrexate (MTX) and mycophenolate mofetil (MMF) in a multicistronic lentiviral vector for in vivo T lymphocyte selection. We found that co-expression of human dihydrofolate reductase (DHFR(FS); L22F, F31S) and inosine monophosphate dehydrogenase II (IMPDH2(IY); T333I, S351Y) conferred T cell resistance to the cytocidal and anti-proliferative effects of these drugs at concentrations that can be achieved clinically (up to 0.1 µM MTX and 1.0 µM MPA). Furthermore, using a immunodeficient mouse model that supports the engraftment of central memory derived human T cells, in vivo selection studies demonstrate that huEGFRt(+)DHFR(FS+)IMPDH2(IY+) T cells could be enriched following adoptive transfer either by systemic administration of MTX alone (4.4 -fold), MMF alone (2.9-fold), or combined MTX and MMF (4.9-fold). These findings demonstrate the utility of both DHFR(FS)/MTX and IMPDH2(IY)/MMF for in vivo selection of lentivirally transduced human T cells. Vectors incorporating these muteins in combination with other therapeutic transgenes may facilitate the selective engraftment of therapeutically active cells in recipients.

Efficient selection of genetically modified human T cells using methotrexate-resistant human dihydrofolate reductase

Genetic modification of human T cells to express transgene-encoded polypeptides, such as tumor targeting chimeric antigen receptors, is an emerging therapeutic modality showing promise in clinical trials. The development of simple and efficient techniques for purifying transgene⁺ T cells is needed to facilitate the derivation of cell products with uniform potency and purity. Unlike selection platforms that utilize physical methods (immunomagnetic or sorting) that are technically cumbersome and limited by the expense and availability of clinical-grade components, we focused on designing a selection system based on the pharmaceutical drug methotrexate (MTX), a potent allosteric inhibitor of human dihydrofolate reductase (DHFR). Here, we describe the development of SIN lentiviral vectors that direct the coordinated expression of a CD19-specific CAR, the human EGFRt tracking/suicide construct, and a methotrexate-resistant human DHFR mutein (huDHFR^FS; L22F, F31S). Our results demonstrate that huDHFR^FS co-expression renders lentivirally transduced primary human CD45RO⁺CD62L⁺ central memory T cells resistant to lymphotoxic concentrations of MTX up to 0.1 µM. Our modular cDNA design insures that selected MTX-resistant T cells co-express functionally relevant levels of the CD19-specific CAR and EGFRt. This selection system based on huDHFR^FS and MTX has considerable potential utility in the manufacturing of clinical-grade T cell products.

Selectively weakened binding of methotrexate by human dihydrofolate reductase allows rapid ex vivo selection of mammalian cells

Ex vivo selection of transduced hematopoietic stem cells (HSC) with drug-resistance genes offers the possibility to enrich transduced cells prior to engraftment, toward increased reconstitution in transplant recipients. We evaluated the potential of highly methotrexate (MTX)-resistant variants of human dihydrofolate reductase (hDHFR) for this application. Two subsets of hDHFR variants with reduced affinity for MTX that had been previously identified in a bacterial system were considered: those with substitutions at positions 31, 34, and/or 35, and those with substitutions at position 115. The variants were characterized for their resistance to pemetrexed (PMTX), an antifolate that is related to MTX. We observed a strong correlation between decreased binding to both antifolates, although the identity of specific sequence variations modulated the correlation. We chose a subset of hDHFR variants for tests of ex vivo MTX resistance, taking into consideration their residual specific activity and their decrease in affinity for the related antifolates. Murine myeloid progenitors and other differentiated hematopoietic cells were transduced and exposed to MTX in a nucleotide-free medium. Bone marrow (BM) cells including 15% cells infected with F31R/Q35E were enriched to 98% transduced cells within 6 days of ex vivo selection. hDHFR variant F31R/Q35E allowed a strong ex vivo enrichment upon a short exposure to MTX relative to a less resistant variant of hDHFR, L22Y. We have thus demonstrated that bacterial selection of highly antifolate-resistant hDHFR variants can provide selectable markers for rapid ex vivo enrichment of hematopoietic cells.

Correlations of inhibitor kinetics for Pneumocystis jirovecii and human dihydrofolate reductase with structural data for human active site mutant enzyme complexes

To understand the role of specific active site residues in conferring selective dihydrofolate reductase (DHFR) inhibition from pathogenic organisms such as Pneumocystis carinii (pc) or Pneumocystis jirovecii (pj), the causative agent in AIDS pneumonia, it is necessary to evaluate the role of these residues in the human enzyme. We report the first kinetic parameters for DHFR from pjDHFR and pcDHFR with methotrexate (MTX), trimethoprim (TMP), and its potent analogue, PY957. We also report the mutagenesis and kinetic analysis of active site mutant proteins at positions 35 and 64 of human (h) DHFR and the crystal structure determinations of hDHFR ternary complexes of NADPH and PY957 with the wild-type DHFR enzyme, the single mutant protein, Gln35Lys, and two double mutant proteins, Gln35Ser/Asn64Ser and Gln35Ser/Asn64Phe. These substitutions place into human DHFR amino acids found at those sites in the opportunistic pathogens pcDHFR (Q35K/N64F) and pjDHFR (Q35S/N64S). The K(i) inhibition constant for PY957 showed greatest potency of the compound for the N64F single mutant protein (5.2 nM), followed by wild-type pcDHFR (K(i) 22 nM) and then wild-type hDHFR enzyme (K(i) 230 nM). Structural data reveal significant conformational changes in the binding interactions of PY957 in the hDHFR Q35S/N64F mutant protein complex compared to the other hDHFR mutant protein complexes and the pcDHFR ternary complex. The conformation of PY957 in the wild-type DHFR is similar to that observed for the single mutant protein. These data support the hypothesis that the enhanced selectivity of PY957 for pcDHFR is in part due to the contributions at positions 37 and 69 (pcDHFR numbering). This insight will help in the design of more selective inhibitors that target these opportunistic pathogens.

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.