Interactions of 5-deazapteridine derivatives with Mycobacterium tuberculosis and with human dihydrofolate reductases

There are major differences between the structures of human dihydrofolate reductase (hDHFR) and Mycobacterium tuberculosis dihydrofolate reductase (mtDHFR). These differences may allow us to design more selective mtDHFR inhibitors. In this paper we study the reactions of six different compounds derived from 5-deazapteridine with human and bacterial enzymes. Results suggest that the addition of hydrophobic groups to the aminophenyl ring would increase mtDHFR-inhibitor affinity and selectivity.     

Interactions of 5-deazapteridine Derivatives with Mycobacterium tuberculosis and with Human Dihydrofolate Reductases

Abstract

There are major differences between the structures of human dihydrofolate reductase (hDHFR) and Mycobacterium tuberculosis dihydrofolate reductase (mtDHFR). These differences may allow us to design more selective mtDHFR inhibitors. In this paper we study the reactions of six different compounds derived from 5-deazapteridine with human and bacterial enzymes. Results suggest that the addition of hydrophobic groups to the aminophenyl ring would increase mtDHFR-inhibitor affinity and selectivity.     

Atomic structures of human dihydrofolate reductase complexed with NADPH and two lipophilic antifolates at 1.09 a and 1.05 a resolution

The crystal structures of two human dihydrofolate reductase (hDHFR) ternary complexes, each with bound NADPH cofactor and a lipophilic antifolate inhibitor, have been determined at atomic resolution. The potent inhibitors 6-([5-quinolylamino]methyl)-2,4-diamino-5-methylpyrido[2,3-d]pyrimidine (SRI-9439) and (Z)-6-(2-[2,5-dimethoxyphenyl]ethen-1-yl)-2,4-diamino-5-methylpyrido[2,3-d]pyrimidine (SRI-9662) were developed at Southern Research Institute against Toxoplasma gondii DHFR-thymidylate synthase. The 5-deazapteridine ring of each inhibitor adopts an unusual puckered conformation that enables the formation of identical contacts in the active site. Conversely, the quinoline and dimethoxybenzene moieties exhibit distinct binding characteristics that account for the differences in inhibitory activity. In both structures, a salt-bridge is formed between Arg70 in the active site and Glu44 from a symmetry-related molecule in the crystal lattice that mimics the binding of methotrexate to DHFR.     

An Innovative Strategy for Dual Inhibitor Design and Its Application in Dual Inhibition of Human Thymidylate Synthase and Dihydrofolate Reductase Enzymes

Due to the diligence of inherent redundancy and robustness in many biological networks and pathways, multitarget inhibitors present a new prospect in the pharmaceutical industry for treatment of complex diseases. Nevertheless, to design multitarget inhibitors is concurrently a great challenge for medicinal chemists. We have developed a novel computational approach by integrating the affinity predictions from structure-based virtual screening with dual ligand-based pharmacophore to discover potential dual inhibitors of human Thymidylate synthase (hTS) and human dihydrofolate reductase (hDHFR). These are the key enzymes in folate metabolic pathway that is necessary for the biosynthesis of RNA, DNA, and protein. Their inhibition has found clinical utility as antitumor, antimicrobial, and antiprotozoal agents. A druglike database was utilized to perform dual-target docking studies. Hits identified through docking experiments were mapped over a dual pharmacophore which was developed from experimentally known dual inhibitors of hTS and hDHFR. Pharmacophore mapping procedure helped us in eliminating the compounds which do not possess basic chemical features necessary for dual inhibition. Finally, three structurally diverse hit compounds that showed key interactions at both active sites, mapped well upon the dual pharmacophore, and exhibited lowest binding energies were regarded as possible dual inhibitors of hTS and hDHFR. Furthermore, optimization studies were performed for final dual hit compound and eight optimized dual hits demonstrating excellent binding features at target systems were also regarded as possible dual inhibitors of hTS and hDHFR. In general, the strategy used in the current study could be a promising computational approach and may be generally applicable to other dual target drug designs.     

Role of lysine-54 in determining cofactor specificity and binding in human dihydrofolate reductase

Lysine-54 of human dihydrofolate reductase (hDHFR) appears to be involved in the interaction with the 2'-phosphate of NADPH and is conserved as a basic residue in other species. Studies have suggested that in Lactobacillus casei dihydrofolate reductase Arg-43, the homologous residue at this position, plays an important role in the binding of NADPH and in the differentiation of Km values for NADPH and NADH. A Lys-54 to Gln-54 mutant (K54Q) of hDHFR has been constructed by oligodeoxynucleotide-directed mutagenesis in order to study the role of Lys-54 in differentiating Km and Kcat values for NADPH and NADH as well as in other functions of hDHFR. The purpose of this paper is to delineate in quantitative terms the magnitude of the effect of the Lys-54 to Gln-54 replacement on the various kinetic parameters of hDHFR. Such quantitative effects cannot be predicted solely on the basis of X-ray structures. The Km for NADPH for the K54Q mutant enzyme is 58-fold higher, while the Km for NADH for K54Q is only 3.9-fold higher than that of the wild type, indicating that the substitution of Lys-54 with Gln-54 decreases the apparent affinity of the enzyme for NADPH dramatically, but has a lesser effect on the apparent affinity for NADH.(ABSTRACT TRUNCATED AT 250 WORDS)     

Binding mode analysis of dual inhibitors of human thymidylate synthase and dihydrofolate reductase as antitumour agents via molecular docking and DFT studies

Due to the diligence of inherent redundancy and robustness in many biological networks and pathways, multitarget inhibitors present a new prospect in the pharmaceutical industry for treatment of complex diseases. Nevertheless, to design multitarget inhibitors is concurrently a great challenge for medicinal chemists. Human thymidylate synthase (hTS) and human dihydrofolate reductase (hDHFR) are the key enzymes in folate metabolic pathway that is necessary for the biosynthesis of RNA, DNA and protein. Their inhibition has found clinical utility as antitumour, antimicrobial and antiprotozoal agents. The aim of this study is to elucidate the factors which are responsible for the potent inhibition of hTS and hDHFR, respectively, through the detailed analysis of the binding modes of dual TS–DHFR inhibitors at both active sites using molecular docking study. Moreover, this study is also accompanied by the exploration of electronic features of dual inhibitors via the density functional theory approach. This study demonstrates that appropriate substitution at the sixth position of thieno[2,3-d]pyrimidines moiety in non-classical dual inhibitors of hTS and hDHFR plays a key role in the inhibition of hTS and hDHFR enzymes. In general, the outcomes of this research exertion will significantly be helpful in drug design for cancer chemotherapy.     

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. Proteins 2015; 83:2124–2136. © 2015 Wiley Periodicals, Inc.     

Effect of enzyme and ligand protonation on the binding of folates to recombinant human dihydrofolate reductase: implications for the evolution of eukaryotic enzyme efficiency.

There is marked pH dependence of the rate constant (koff) for tetrahydrofolate (H4folate) dissociation from its ternary complex with human dihydrofolate reductase (hDHFR) and NADPH. Similar pH dependence of H4folate dissociation from the ternary complex of a variant of hDHFR with the substitution Phe31----Leu (F31L hDHFR) causes this dissociation to become rate limiting in the enzyme mechanism at pH approximately 5, and this accounts for the marked decrease in kcat for this variant as the pH is decreased from 7 to 5. This decreased kcat at low pH is not seen for most DHFRs. koff for dissociation of folate, dihydrofolate (H2folate), and H4folate from their binary complexes with hDHFR is similarly pH dependent. For all the complexes examined, the pH dependence of koff in the range pH 5-7 is well described by a pKa of about 6.2 and must be due to ionization of a group on the enzyme. In the higher pH range (7-10), koff increases further as the pH is raised, and this relation is governed by a second pKa which is close to the pKa for ionization of the amide group (HN3-C4O) of the respective ligands. Thus, ionization of the ligand amide group also increases koff. Evidence is presented that the dependence of pH on koff for hDHFR accounts for the shape of the kcat versus pH curve for both hDHFR as well as its F31L variant and contributes to the higher efficiency of hDHFR compared with bacterial DHFR.     

Integration and expression of the human dihydrofolate reductase gene in the Bacillus subtilis chromosome

Integration of functionally active human dihydrofolate reductase (hDHFR) gene into the Bacillus subtilis chromosome was performed. The clones obtained contained 1 to 7 copies of hDHFR gene per chromosome equivalent and were resistant to trimethoprim. In cell lysates of such clones a protein with the molecular mass of hDHFR was detected. The hDHFR gene was stably maintained in all clones having this gene integrated into the bacterial chromosome, when grown under non-selective conditions.     

Design, synthesis, biological evaluation and X-ray crystal structure of novel classical 6,5,6-tricyclic benzo[4,5]thieno[2,3-d]pyrimidines as dual thymidylate synthase and dihydrofolate reductase inhibitors

Classical antifolates (4-7) with a tricyclic benzo[4,5]thieno[2,3-d]pyrimidine scaffold and a flexible and rigid benzoylglutamate were synthesized as dual thymidylate synthase (TS) and dihydrofolate reductase (DHFR) inhibitors. Oxidative aromatization of ethyl 2-amino-4-methyl-4,5,6,7-tetrahydro-1-benzothiophene-3-carboxylate (±)-9 to ethyl 2-amino-4-methyl-1-benzothiophene-3-carboxylate 10 with 10% Pd/C was a key synthetic step. Compounds with 2-CH? substituents inhibited human (h) TS (IC?? =0.26-0.8 μM), but not hDHFR. Substitution of the 2-CH? with a 2-NH? increases hTS inhibition by more than 10-fold and also affords excellent hDHFR inhibition (IC?? = 0.09-0.1 μM). This study shows that the tricyclic benzo[4,5]thieno[2,3-d]pyrimidine scaffold is highly conducive to single hTS or dual hTS-hDHFR inhibition depending on the 2-position substituents. The X-ray crystal structures of 6 and 7 with hDHFR reveal, for the first time, that tricyclics 6 and 7 bind with the benzo[4,5]thieno[2,3-d]pyrimidine ring in the folate binding mode with the thieno S mimicking the 4-amino of methotrexate.     

Mapping the antigenic epitopes of human dihydrofolate reductase by systematic synthesis of peptides on solid supports.

All of the 181 possible overlapping hexapeptides as well as 179 octapeptides covering the amino acid sequence of human dihydrofolate reductase (hDHFR) were synthesized on polyethylene supports. The synthetic procedure of Geysen et al. (Geysen, H. M., Rodda, S. J., Mason, T. J., Tribbick, G., and Schoofs, P. G. (1987) J. Immunol. Methods 102, 259-274) was modified to obtain up to 100 nmol of peptide on each pin. Peptides constituting antigenic epitopes on hDHFR were identified by examining the binding of antibodies raised against both native and denatured hDHFR to these peptides by enzyme-linked immunosorbent assay. The peptides bound in a similar pattern to polyclonal antibodies against both native and denatured dihydrofolate reductase (DHFR). Six major epitopes were located corresponding to residues 27-33, 45-51, 67-74, 133-139, 153-158, and 176-181 using both hexapeptides and octapeptides. An additional epitope, constituting residues 14-21, was found by the use of octapeptides. Most of the epitopes are hydrophilic and reside largely in "loop" regions at the boundaries of secondary structural elements of hDHFR. This observation is consistent with our previous results which suggested that ligand binding at the active site of the enzyme can cause a dramatic reduction in antibody binding to DHFR due to conformational constraints in flexible loop regions in various parts of the molecule. The similarity of the immunogenic profiles of native versus denatured hDHFR indicates that the two forms of the antigen share the same amino acid sequence-specific epitopes. Competitive enzyme-linked immunosorbent assay showed that the binding of anti-hDHFR antiserum to both native and denatured hDHFR was inhibited by approximately 30% by the seven antigenic peptides, indicating that a significant proportion of the antibodies elicited by this enzyme is specific for short peptides. Besides revealing the antigenic structure of DHFR our results provide a rational basis for the design of mutant DHFRs to study the importance of loop residues in the conformational dynamics of the enzyme.     

Various means of integration of the expressible human dihydrofolate reductase gene into the Bacillus subtilis genome

Integration of expressible DNA corresponding to the human dihydrofolate reductase (hDHFR) gene into the Bacillus subtilis genome has been achieved in different ways. The clones obtained contained one to seven copies of this gene per genome equivalent and were resistant to trimethoprim. Clones produce a new protein coded by the integrated hDHFR gene. In all clones, the integrated DNA was stably maintained even under nonselective growth conditions.     

Refined crystal structures of Escherichia coli and chicken liver dihydrofolate reductase containing bound trimethoprim

Refined crystal structures are reported for complexes of Escherichia coli and chicken dihydrofolate reductase containing the antibiotic trimethoprim (TMP). Structural comparison of these two complexes reveals major geometrical differences in TMP binding that may be important in understanding the stereo-chemical basis of this inhibitor's selectivity for bacterial dihydrofolate reductases. For TMP bound to chicken dihydrofolate reductase we observe an altered binding geometry in which the 2,4-diaminopyrimidine occupies a position in closer proximity (by approximately 1 A) to helix alpha B compared to the pyrimidine position for TMP or methotrexate bound to E. coli dihydrofolate reductase. One important consequence of this deeper insertion of the pyrimidine into the active site of chicken dihydrofolate reductase is the loss of a potential hydrogen bond that would otherwise form between the carbonyl oxygen of Val-115 and the inhibitor's 4-amino group. In addition, for TMP bound to E. coli dihydrofolate reductase, the inhibitor's benzyl side chain is positioned low in the active-site pocket pointing down toward the nicotinamide-binding site, whereas, in chicken dihydrofolate reductase, the benzyl group is accommodated in a side channel running upward and away from the cofactor. As a result, the torsion angles about the C5-C7 and C7-C1' bonds for TMP bound to the bacterial reductase (177 degrees, 76 degrees) differ significantly from the corresponding angles for TMP bound to chicken dihydrofolate reductase (-85 degrees, 102 degrees). Finally, when TMP binds to the chicken holoenzyme, the Tyr-31 side chain undergoes a large conformational change (average movement is 5.4 A for all atoms beyond C beta), rotating down into a new position where it hydrogen bonds via an intervening water molecule to the backbone carbonyl oxygen of Trp-24.     

Conversion of arginine to lysine at position 70 of human dihydrofolate reductase: generation of a methotrexate-insensitive mutant enzyme

Arginine-70 of human dihydrofolate reductase (hDHFR) is a highly conserved residue which X-ray crystallographic data have shown to interact with the alpha-carboxylate of the terminal L-glutamate moiety of either folic acid or methotrexate (MTX). The rationale for this study was to introduce a conservative amino acid residue change at position 70 (Arg----Lys) which might function as a titratable group and, thus, reveal possible quantitative changes in ligand binding and kinetic parameters as a function of pH. Such a mutant enzyme (R70K) has been constructed and expressed by using site-directed mutagenesis techniques. This substitution has a dramatic effect on the binding of MTX, which displays a 22,600-fold increase in the dissociation constant (KD) at pH 7.5 compared to that of the reported wild-type enzyme value. At this pH, the KD value for dihydrofolate (FAH2) for the R70K enzyme shows only a 7-fold increase over that for the wild-type hDHFR. The pH profiles of the Michaelis and dissociation constants for FAH2 and KD values for MTX for the mutant enzyme all show a 7-8-fold increase from pH 7.5 to 8.5 as compared to its wild-type counterpart. The binding of NADPH or the nonclassical inhibitor trimetrexate (TMQ) to either the wild-type or the mutant enzyme does not show such pH-dependent characteristics. Thus, since FAH2 and MTX interact with the guanidinium side chain of arginine-70 in the wild-type hDHFR, the replacement of this residue with a lysine in the R70K mutant appears to have resulted in the introduction of a titratable group with a perturbed pKa value of ca. 8.3.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

Species-specific irreversible inhibition of Neisseria gonorrhoeae dihydrofolate reductase by a substituted 2,4-diamino-5-benzylpyrimidine

Neisseria gonorrhoeae dihydrofolate reductase undergoes a time-dependent, irreversible inactivation by 2,4-diamino-5-[3,5-dimethoxy-4-(p-bromoacetamidophenoxy)benzyl] pyrimidine. The kinetics of inactivation are consistent with the reversible formation of an enzyme-inhibitor complex followed by covalent binding to the enzyme. The reversible component is competitive with dihydrofolate and has an inhibitor binding constant of 10 nM. Irreversible inactivation proceeds as a pseudo first-order process with a minimum inactivation half-time of 20 min and a Ki of 28 nM. Using radiolabeled inhibitor, it was shown that approximately 1 mol of ligand was covalently bound to the enzyme/mol of methotrexate binding site when the enzyme was completely inhibited. Radiolabeled inhibitor remained associated with the enzyme following denaturation and sodium dodecyl sulfate-polyacrylamide gel electrophoresis. Cyanogen bromide cleavage of the 14C-labeled enzyme-inhibitor complex yielded only one radioactive polypeptide, and sequence determinations showed that His-25 was modified by covalent attachment of the inhibitor. When dihydrofolate reductases from Lactobacillus casei, Streptococcus faecium, Escherichia coli, SR-1 rodent lymphoma, and chicken liver were tested with the affinity label, only the L. casei enzyme showed a time-dependent increase in inhibition. These data, along with comparisons of known amino acid sequences and x-ray crystal structures, were used to make predictions concerning the three-dimensional conformation of the gonococcal enzyme.     

Crystallographic analysis reveals a novel second binding site for trimethoprim in active site double mutants of human dihydrofolate reductase

In order to produce a more potent replacement for trimethoprim (TMP) used as a therapy for Pneumocystis pneumonia and targets dihydrofolate reductase from Pneumocystis jirovecii (pjDHFR), it is necessary to understand the determinants of potency and selectivity against DHFR from the mammalian host and fungal pathogen cells. To this end, active site residues in human (h) DHFR were replaced with those from pjDHFR. Structural data are reported for two complexes of TMP with the double mutants Gln35Ser/Asn64Phe (Q35S/N64F) and Gln35Lys/Asn64Phe (Q35K/N64F) of hDHFR that unexpectedly show evidence for the binding of two molecules of TMP: one molecule that binds in the normal folate binding site and the second molecule that binds in a novel subpocket site such that the mutated residue Phe64 is involved in van der Waals contacts to the trimethoxyphenyl ring of the second TMP molecule. Kinetic data for the binding of TMP to hDHFR and pjDHFR reveal an 84-fold selectivity of TMP against pjDHFR (K_i 49 nM) compared to hDHFR (K_i 4093 nM). Two mutants that contain one substitution from pj- and one from the closely related Pneumocystis carinii DHFR (pcDHFR) (Q35K/N64F and Q35S/N64F) show K_i values of 593 and 617 nM, respectively; these K_i values are well above both the K_i for pjDHFR and are similar to pcDHFR (Q35K/N64F and Q35S/N64F) (305 nM). These results suggest that active site residues 35 and 64 play key roles in determining selectivity for pneumocystis DHFR, but that other residues contribute to the unique binding of inhibitors to these enzymes.     

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.     

Dihydrofolate reductase. The stereochemistry of inhibitor selectivity

X-ray structural results are reported for 10 triazine and pyrimidine inhibitors of dihydrofolate reductase, each one studied as a ternary complex with NADPH and chicken dihydrofolate reductase. Analysis of these data and comparison with structural results from the preceding paper (Matthews, D.A., Bolin, J.T., Burridge, J.M., Filman, D.J., Volz, K.W., Kaufman, B. T., Beddell, C.R., Champness, J.N., Stammers, D.K., and Kraut, J. (1985) J. Biol. Chem. 260, 381-391) in which we contrasted binding of the antibiotic trimethoprim (TMP) to chicken dihydrofolate reductase on the one hand with its binding to Escherichia coli dihydrofolate reductase on the other, permit identification of differences that are important in accounting for TMP's selectivity. The crystallographic evidence strongly suggests that loss of a potential hydrogen bond between the 4-amino group of TMP and the backbone carbonyl of Val-115 when TMP binds to chicken dihydrofolate reductase but not when it binds to the E. coli reductase is the major factor responsible for this drug's more potent inhibition of bacterial dihydrofolate reductase. A key finding of the current study which is important in understanding why TMP binds differently to chicken and E. coli dihydrofolate reductases is that residues on opposite sides of the active-site cleft in chicken dihydrofolate reductase are about 1.5-2.0 A further apart than are structurally equivalent residues in the E. coli enzyme.     

Validation of a new restraint docking method for solution structure determinations of protein–ligand complexes

A new method is proposed for docking ligands into proteins in cases where an NMR-determined solution structure of a related complex is available. The method uses a set of experimentally determined values for protein–ligand, ligand–ligand, and protein–protein restraints for residues in or near to the binding site, combined with a set of protein–protein restraints involving all the other residues which is taken from the list of restraints previously used to generate the reference structure of a related complex. This approach differs from ordinary docking methods where the calculation uses fixed atomic coordinates from the reference structure rather than the restraints used to determine the reference structure. The binding site residues influenced by replacing the reference ligand by the new ligand were determined by monitoring differences in ¹H chemical shifts. The method has been validated by showing the excellent agreement between structures of L. casei dihydrofolate reductase.trimetrexate calculated by conventional methods using a full experimentally determined set of restraints and those using this new restraint docking method based on an L. casei dihydrofolate reductase.methotrexate reference structure.