Design, synthesis and biological evaluation of novel pyrrolo[2,3-d]pyrimidine as tumor-targeting agents with selectivity for tumor uptake by high affinity folate receptors over the reduced folate carrier.

Tumor-targeted 6-substituted pyrrolo[2,3-d]pyrimidine benzoyl compounds based on 2 were isosterically modified at the 4-carbon bridge by replacing the vicinal (C11) carbon by heteroatoms N (4), O (5) or S (6), or with an N-substituted formyl (7), trifluoroacetyl (8) or acetyl (9). Replacement with sulfur (6) afforded the most potent KB tumor cell inhibitor, ~6-fold better than the parent 2. In addition, 6 retained tumor transport selectivity via folate receptor (FR) α and -ß over the ubiquitous reduced folate carrier (RFC). FRα-mediated cell inhibition for 6 was generally equivalent to 2, while the FRß-mediated activity was improved by 16-fold over 2. N (4) and O (5) substitutions afforded similar tumor cell inhibitions as 2, with selectivity for FRα and -ß over RFC. The N-substituted analogs 7-9 also preserved transport selectivity for FRα and -ß over RFC. For FRα-expressing CHO cells, potencies were in the order of 8 > 7 > 9. Whereas 8 and 9 showed similar results with FRß-expressing CHO cells, 7 was ~16-fold more active than 2. By nucleoside rescue experiments, all the compounds inhibited de novo purine biosynthesis, likely at the step catalyzed by glycinamide ribonucleotide formyltransferase. Thus, heteroatom replacements of the CH2 in the bridge of 2 afford analogs with increased tumor cell inhibition that could provide advantages over 2, as well as tumor transport selectivity over clinically used antifolates including methotrexate and pemetrexed.

Structural and Enzymatic Analysis of Tumor-Targeted Antifolates That Inhibit Glycinamide Ribonucleotide Formyltransferase.

Pemetrexed and methotrexate are antifolates used for cancer chemotherapy and inflammatory diseases. These agents have toxic side effects resulting, in part, from nonspecific cellular transport by the reduced folate carrier (RFC), a ubiquitously expressed facilitative transporter. We previously described 2-amino-4-oxo-6-substituted pyrrolo[2,3-d]pyrimidine antifolates with modifications of the side chain linker and aromatic ring that are poor substrates for RFC but are efficiently transported via folate receptors (FRs) and the proton-coupled folate transporter (PCFT). These targeted antifolates are cytotoxic in vitro toward FR- and PCFT-expressing tumor cells and in vivo with human tumor xenografts in immune-compromised mice, reflecting selective cellular uptake. Antitumor efficacy is due to inhibition of glycinamide ribonucleotide (GAR) formyltransferase (GARFTase) activity in de novo synthesis of purine nucleotides. This study used purified human GARFTase (formyltransferase domain) to assess in vitro inhibition by eight novel thieno- and pyrrolo[2,3-d]pyrimidine antifolates. Seven analogues (AGF23, AGF71, AGF94, AGF117, AGF118, AGF145, and AGF147) inhibited GARFTase with Ki values in the low- to mid-nanomolar concentration range, whereas AGF50 inhibited GARFTase with micromolar potency similar to that of PMX. On the basis of crystal structures of ternary complexes with GARFTase, ß-GAR, and the monoglutamyl antifolates, differences in inhibitory potencies correlated well with antifolate binding and the positions of the terminal carboxylates. Our data provide a mechanistic basis for differences in inhibitory potencies between these novel antifolates and a framework for future structure-based drug design. These analogues could be more efficacious than clinically used antifolates, reflecting their selective cellular uptake by FRs and PCFT and potent GARFTase inhibition.

6-Substituted Pyrrolo[2,3-d]pyrimidine Thienoyl Regioisomers as Targeted Antifolates for Folate Receptor α and the Proton-Coupled Folate Transporter in Human Tumors.

2-Amino-4-oxo-6-substituted-pyrrolo[2,3-d]pyrimidine antifolate thiophene regioisomers of AGF94 (4) with a thienoyl side chain and three-carbon bridge lengths [AGF150 (5) and AGF154 (7)] were synthesized as potential antitumor agents. These analogues inhibited proliferation of Chinese hamster ovary (CHO) sublines expressing folate receptors (FRs) α or ß (IC50s < 1 nM) or the proton-coupled folate transporter (PCFT) (IC50 < 7 nM). Compounds 5 and 7 inhibited KB, IGROV1, and SKOV3 human tumor cells at subnanomolar concentrations, reflecting both FRα and PCFT uptake. AGF152 (6) and AGF163 (8), 2,4-diamino-5-substituted-furo[2,3-d]pyrimidine thiophene regioisomers, also inhibited growth of FR-expressing CHO and KB cells. All four analogues inhibited glycinamide ribonucleotide formyltransferase (GARFTase). Crystal structures of human GARFTase complexed with 5 and 7 were reported. In severe combined immunodeficient mice bearing SKOV3 tumors, 7 was efficacious. The selectivity of these compounds for PCFT and for FRα and ß over the ubiquitously expressed reduced folate carrier is a paradigm for selective tumor targeting.

Structures and reaction mechanisms of the two related enzymes, PurN and PurU.

The crystal structures of glycinamide ribonucleotide transformylases (PurNs) from Aquifex aeolicus (Aa), Geobacillus kaustophilus (Gk) and Symbiobacterium toebii (St), and of formyltetrahydrofolate hydrolase (PurU) from Thermus thermophilus (Tt) were determined. The monomer structures of the determined PurN and PurU were very similar to the known structure of PurN, but oligomeric states were different; AaPurN and StPurN formed dimers, GkPurN formed monomer and PurU formed tetramer in the crystals. PurU had a regulatory ACT domain in its N-terminal side. So far several structures of PurUs have been determined, yet, the mechanisms of the catalysis and the regulation of PurU have not been elucidated. We, therefore, modelled ligand-bound structures of PurN and PurU, and performed molecular dynamics simulations to elucidate the reaction mechanisms. The evolutionary relationship of the two enzymes is discussed based on the comparisons of the structures and the catalytic mechanisms.

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 (IC50 = 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 (IC50 = 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.

Predicting flexible loop regions that interact with ligands: the challenge of accurate scoring.

Flexible loop regions play a critical role in the biological function of many proteins and have been shown to be involved in ligand binding. In the context of structure-based drug design, using or predicting an incorrect loop configuration can be detrimental to the study if the loop is capable of interacting with the ligand. Three protein systems, each with at least one flexible loop region in close proximity to the known binding site, were selected for loop prediction using the CorLps program; a six residue loop region from phosphoribosylglycinamide formyltransferase (GART), two nine residue loop regions from cytochrome P450 (CYP) 119, and an 11 residue loop region from enolase were selected for loop prediction. The results of this study indicate that the statistically based DFIRE scoring function implemented in the CorLps program did not accurately rank native-like predicted loop configurations in any protein system. In an attempt to improve the ranking of the native-like predicted loop configurations, the MM/GBSA and the optimized MM/GBSA-dsr scoring functions were used to re-rank the predicted loops with and without bound ligand. In general, single snapshot MM/GBSA scoring provided the best ranking of native-like loop configurations. Based on the scoring function analyses presented, the optimal ranking of native-like loop configurations is still a difficult challenge and the choice of the "best" scoring function appears to be system dependent.

Purine biosynthesis in archaea: variations on a theme.

BACKGROUND: The ability to perform de novo biosynthesis of purines is present in organisms in all three domains of life, reflecting the essentiality of these molecules to life. Although the pathway is quite similar in eukaryotes and bacteria, the archaeal pathway is more variable. A careful manual curation of genes in this pathway demonstrates the value of manual curation in archaea, even in pathways that have been well-studied in other domains. RESULTS: We searched the Integrated Microbial Genome system (IMG) for the 17 distinct genes involved in the 11 steps of de novo purine biosynthesis in 65 sequenced archaea, finding 738 predicted proteins with sequence similarity to known purine biosynthesis enzymes. Each sequence was manually inspected for the presence of active site residues and other residues known or suspected to be required for function.Many apparently purine-biosynthesizing archaea lack evidence for a single enzyme, either glycinamide ribonucleotide formyltransferase or inosine monophosphate cyclohydrolase, suggesting that there are at least two more gene variants in the purine biosynthetic pathway to discover. Variations in domain arrangement of formylglycinamidine ribonucleotide synthetase and substantial problems in aminoimidazole carboxamide ribonucleotide formyltransferase and inosine monophosphate cyclohydrolase assignments were also identified.Manual curation revealed some overly specific annotations in the IMG gene product name, with predicted proteins without essential active site residues assigned product names implying enzymatic activity (21 proteins, 2.8% of proteins inspected) or Enzyme Commission (E. C.) numbers (57 proteins, 7.7%). There were also 57 proteins (7.7%) assigned overly generic names and 78 proteins (10.6%) without E.C. numbers as part of the assigned name when a specific enzyme name and E. C. number were well-justified. CONCLUSIONS: The patchy distribution of purine biosynthetic genes in archaea is consistent with a pathway that has been shaped by horizontal gene transfer, duplication, and gene loss. Our results indicate that manual curation can improve upon automated annotation for a small number of automatically-annotated proteins and can reveal a need to identify further pathway components even in well-studied pathways.

Expression, crystallization and preliminary X-ray analysis of the phosphoribosylglycinamide formyltransferase from Streptococcus mutans.

Phosphoribosylglycinamide formyltransferase (PurN) from Streptococcus mutans was recombinantly expressed in Escherichia coli. An effective purification protocol was established. The purified protein, which had a purity of >95%, was identified by SDS-PAGE and MALDI-TOF MS. The protein was crystallized using the vapour-diffusion method in hanging-drop mode with PEG 3350 as the primary precipitant. X-ray diffraction data were collected to 2.1 Šresolution. Preliminary X-ray analysis indicated that the crystal belonged to space group P2(1)2(1)2(1), with unit-cell parameters a = 52.25, b = 63.29, c = 131.81 Å.

Structural studies of tri-functional human GART.

Human purine de novo synthesis pathway contains several multi-functional enzymes, one of which, tri-functional GART, contains three enzymatic activities in a single polypeptide chain. We have solved structures of two domains bearing separate catalytic functions: glycinamide ribonucleotide synthetase and aminoimidazole ribonucleotide synthetase. Structures are compared with those of homologous enzymes from prokaryotes and analyzed in terms of the catalytic mechanism. We also report small angle X-ray scattering models for the full-length protein. These models are consistent with the enzyme forming a dimer through the middle domain. The protein has an approximate seesaw geometry where terminal enzyme units display high mobility owing to flexible linker segments. This resilient seesaw shape may facilitate internal substrate/product transfer or forwarding to other enzymes in the pathway.

Protein preparation, crystallization and preliminary crystallographic studies of Bacillus subtilis glycinamide ribonucleotide transformylase.

Glycinamide ribonucleotide transformylase (GART) catalyzes the transfer of a formyl group from formyl tetrahydrofolate (FTHF) to glycinamide ribonucleotide (GAR), which is an essential step in the de novo synthesis pathway of purines. In Bacillus subtilis, GART is encoded by the gene purN. In order to study the structure and function of B. subtilis GART, the purN gene was amplified, cloned into an expression vector and expressed in soluble form in Escherichia coli. The protein was purified to homogeneity and crystals suitable for X-ray data collection were obtained. These crystals diffracted to 2.5 A resolution and belonged to space group P3(1)21, with unit-cell parameters a = b = 95.5, c = 64.0 A.

Structures of glycinamide ribonucleotide transformylase (PurN) from Mycobacterium tuberculosis reveal a novel dimer with relevance to drug discovery.

Enzymes from the de novo purine biosynthetic pathway have been exploited for the development of anti-cancer drugs, and represent novel targets for anti-bacterial drug development. In Mycobacterium tuberculosis, the cause of tuberculosis, this pathway has been identified as essential for growth and survival. The structure of M. tuberculosis PurN (MtPurN) has been determined in complex with magnesium and iodide at 1.30 A resolution, and with cofactor analogue, 5-methyltetrahydrofolate (5MTHF) at 2.2 A resolution. The structure shows a Rossmann-type fold that is very similar to the known structures of the human and E. coli PurN proteins. In contrast, MtPurN forms a dimer that is quite different from that formed by the Escherichia coli PurN, and which suggests a mechanism whereby communication could take place between the two active sites. Differences are seen in two active site loops and in the binding mode of the 5MTHF cofactor analogue between the two MtPurN molecules of the dimer. A binding site for halide ions is found in the dimer interface, and bound magnesium and iodide ions in the active site suggest sites that might be exploited in potential drug discovery strategies.

Information theory-based scoring function for the structure-based prediction of protein-ligand binding affinity.

The development and validation of a new knowledge based scoring function (SIScoreJE) to predict binding energy between proteins and ligands is presented. SIScoreJE efficiently predicts the binding energy between a small molecule and its protein receptor. Protein-ligand atomic contact information was derived from a Non-Redundant Data set (NRD) of over 3000 X-ray crystal structures of protein-ligand complexes. This information was classified for individual "atom contact pairs" (ACP) which is used to calculate the atomic contact preferences. In addition to the two schemes generated in this study we have assessed a number of other common atom-type classification schemes. The preferences were calculated using an information theoretic relationship of joint entropy. Among 18 different atom-type classification schemes "ScoreJE Atom Type set2" (SATs2) was found to be the most suitable for our approach. To test the sensitivity of the method to the inclusion of solvent, Single-body Solvation Potentials (SSP) were also derived from the atomic contacts between the protein atom types and water molecules modeled using AQUARIUS2. Validation was carried out using an evaluation data set of 100 protein-ligand complexes with known binding energies to test the ability of the scoring functions to reproduce known binding affinities. In summary, it was found that a combined SSP/ScoreJE (SIScoreJE) performed significantly better than ScoreJE alone, and SIScoreJE and ScoreJE performed better than GOLD::GoldScore, GOLD::ChemScore, and XScore.

Human glycinamide ribonucleotide transformylase: active site mutants as mechanistic probes.

Human glycinamide ribonucleotide transformylase (GART) (EC 2.1.2.2) is a validated target for cancer chemotherapy, but mechanistic studies of this therapeutically important enzyme are limited. Site-directed mutagenesis, initial velocity studies, pH-rate studies, and substrate binding studies have been employed to probe the role of the strictly conserved active site residues, N106, H108, and D144, and the semiconserved K170 in substrate binding and catalysis. Only two conservative substitutions, N106Q and K170R, resulted in catalytically active enzymes, and these active mutant enzymes gave pH-rate profiles and a steady-state kinetic mechanism essentially identical to those of the native enzyme. All inactive mutants were able to bind both substrates, ruling out disrupted formation of the ternary complex as the source of inactivity. Differences between human and Escherichia coli GART, previously used as a model for the human enzyme, were evident.

The study of helical distortions due to environmental changes: choice of parameters.

Parameters like interhelical angles, helical parameters, levels of distortions, etc., have been analysed to test their sensitivity to environmental changes using a method developed in this laboratory. This analysis was done on protein structures solved under different environmental conditions like temperature and pH, and ligand binding. The study reveals that the helical parameters are not sensitive enough to study the effect of environmental changes on protein helices. On the other hand the helical distortions as well as changes in the interhelical angles are more sensitive to these changes. The study also provides with additional information like the origin of distortions in a helix when a ligand binds to a protein, bending in helical axis, identification and extent of domain movements, etc.

Discovery of a potent, nonpolyglutamatable inhibitor of glycinamide ribonucleotide transformylase.

Glycinamide ribonucleotide transformylase (GAR Tfase) catalyzes the first of two formyl transfer steps in the de novo purine biosynthetic pathway that require folate cofactors. Herein we report the discovery of a potent, nonpolyglutamatable, and selective inhibitor of GAR Tfase. Compound 12, which possesses a tetrazole in place of the gamma-carboxylic acid in the l-glutamate subunit of the potent GAR Tfase inhibitor 1, was active in cellular-based functional assays exhibiting purine-sensitive cytotoxic activity (IC(50) = 40 nM, CCRF-CEM) and was selective for inhibition of rhGAR Tfase (K(i) = 130 nM). Notably, 12 was only 2.5-fold less potent than 1 in cellular assays and 4-fold less potent against rhGAR Tfase. Like 1, this functional activity of 12 in the cell-based assay benefits from and requires transport into the cell by the reduced folate carrier but, unlike 1, is independent of folyl polyglutamate synthase (FPGS) expression levels and polyglutamation.