G-quadruplex formation in the sequences 5′-(TTAGGG)n and 5′(TTAGGG)nTT (n=4,8,12) was studied using circular dichroism, sedimentation velocity, differential scanning calorimetry and molecular dynamics simulations. Sequences containing 8 and 12 repeats formed higher-order structures with two and three contiguous quadruplexes, respectively. Plausible structures for these sequences were determined by molecular dynamics simulations followed by experimental testing of predicted hydrodynamic properties by sedimentation velocity. These structures featured folding of the strand into contiguous quadruplexes with mixed hybrid conformations. Thermodynamic studies showed the strands folded spontaneous to contain the maximum number contiguous quadruplexes. For the sequence 5′(TTAGGG)12TT, more than 90% of the strands contained completely folded structures with three quadruplexes. Statistical mechanical-based deconvolution of thermograms for three quadruplex structures showed that each quadruplex melted independently with unique thermodynamic parmameters. Thermodynamic analysis revealed further that quadruplexes in higher-ordered structures were destabilized relative to their monomeric counterparts, with unfavorable coupling free energies. Quadruplex stability thus depends critically on the sequence and structural context.

SUMMARY

Oxidation of unsaturated phospholipids results in the generation of aldehyde side chains that remain esterified to the phospholipid backbone. Such “core” aldehydes elicit immune responses and promote inflammation. However, the biochemical mechanisms by which phospholipid aldehydes are metabolized or detoxified are not well understood. In the studies reported here, we examined whether aldose reductase (AR), which reduces hydrophobic aldehydes, metabolizes phospholipid aldehydes. Incubation with AR led to the reduction of 5-oxovaleroyl, 7-oxo-5-heptenoyl, 5-hydroxy-6-oxo-caproyl, and 5-hydroxy-8-oxo-6-octenoyl phospholipids generated upon oxidation of 1-palmitoyl-2-arachidonoyl-sn-glycero-3-phosphocholine (PAPC). The enzyme also catalyzed the reduction of phospholipid aldehydes generated from the oxidation of 1-alkyl, and 1-alkenyl analogs of PAPC, and 1-palmitoyl-2-arachidonoyl phosphatidic acid or phosphoglycerol. Aldose reductase catalyzed the reduction of chemically synthesized 1-palmitoyl-2-(5-oxovaleroyl)-sn-glycero-3-phosphatidylcholine (POVPC) with a Km of 10 μM. Addition of POVPC to the culture medium led to incorporation and reduction of the aldehyde in COS-7 and THP-1 cells. Reduction of POVPC in these cells was prevented by the AR inhibitors sorbinil and tolrestat and was increased in COS-7 cells overexpressing AR. Together, these observations suggest that AR may be a significant participant in the metabolism of several structurally diverse phospholipid aldehydes. This metabolism may be a critical regulator of the pro-inflammatory and immunogenic effects of oxidized phospholipids.

Structural knowledge of telomeric DNA is critical for understanding telomere biology and for the utilization of telomeric DNA as a therapeutic target. Very little is known about the structure of long human DNA sequences that may form more than one quadruplex unit. Here we report a combination of molecular dynamics simulations and experimental biophysical studies to explore the structural and dynamic properties of the human telomeric sequence (TTAGGG)8TT that folds into two contiguous quadruplexes. Five higher order quadruplex models were built combining known single human telomeric quadruplex structures as unique building blocks. The biophysical properties of this sequence in K+ solution were experimentally investigated by means of analytical ultracentrifugation and UV spectroscopy. Additionally, the environments of loop adenines were probed by fluorescence studies using systematic single-substitutions of 2-aminopurine for the adenine bases. The comparison of the experimentally determined properties with the corresponding quantities predicted from the models allowed us to test the validity of each of the structural models. One model emerged whose properties are most consistent with the predictions, and which therefore is the most probable structure in solution. This structure features contiguous quadruplex units in an alternating hybrid-1-hybrid-2 conformation with a highly-ordered interface composed of loop residues from both quadruplexes.

Background

T cell activation is associated with a rapid increase in intracellular fructose-2,6-bisphosphate (F2,6BP), an allosteric activator of the glycolytic enzyme, 6-phosphofructo-1-kinase. The steady state concentration of F2,6BP in T cells is dependent on the expression of the bifunctional 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatases (PFKFB1-4) and the fructose-2,6-bisphosphatase, TIGAR. Of the PFKFB family of enzymes, PFKFB3 has the highest kinase:bisphosphatase ratio and has been demonstrated to be required for T cell proliferation. A small molecule antagonist of PFKFB3, 3-(3-pyridinyl)-1-(4-pyridinyl)-2-propen-1-one (3PO), recently has been shown to reduce F2,6BP synthesis, glucose uptake and proliferation in transformed cells. We hypothesized that the induction of PFKFB3 expression may be required for the stimulation of glycolysis in T cells and that exposure to the PFKFB3 antagonist, 3PO, would suppress T cell activation.

Methods

We examined PFKFB1-4 and TIGAR expression and F2,6BP concentration in purified CD3+ T cells stimulated with microbead-conjugated agonist antibodies specific for CD3 and the co-stimulatory receptor, CD28. We then determined the effect of 3PO on anti-CD3/anti-CD28-induced T cell activation, F2,6BP synthesis, 2-[1-14C]-deoxy-d-glucose uptake, lactate secretion, TNF-α secretion and proliferation. Finally, we examined the effect of 3PO administration on the development of delayed type hypersensitivity to methylated BSA and on imiquimod-induced psoriasis in mice.

Results

We found that purified human CD3+ T cells express PFKFB2, PFKFB3, PFKFB4 and TIGAR, and that anti-CD3/anti-CD28 conjugated microbeads stimulated a >20-fold increase in F2,6BP with a coincident increase in protein expression of the PFKFB3 family member and a decrease in TIGAR protein expression. We then found that exposure to the PFKFB3 small molecule antagonist, 3PO (1–10 μM), markedly attenuated the stimulation of F2,6BP synthesis, 2-[1-14C]-deoxy-D-glucose uptake, lactate secretion, TNF-α secretion and T cell aggregation and proliferation. We examined the in vivo effect of 3PO on the development of delayed type hypersensitivity to methylated BSA and on imiquimod-induced psoriasis in mice and found that 3PO suppressed the development of both T cell-dependent models of immunity in vivo.

Conclusions

Our data demonstrate that inhibition of the PFKFB3 kinase activity attenuates the activation of T cells in vitro and suppresses T cell dependent immunity in vivo and indicate that small molecule antagonists of PFKFB3 may prove effective as T cell immunosuppressive agents.

Computational approaches are becoming increasingly popular for the discovery of drug candidates against a target of interest. Proteins have historically been the primary targets of many virtual screening efforts. While in silico screens targeting proteins has proven successful, other classes of targets, in particular DNA, remain largely unexplored using virtual screening methods. With the realization of the functional importance of many non-cannonical DNA structures such as G-quadruplexes, increased efforts are underway to discover new small molecules that can bind selectively to DNA structures. Here, we describe efforts to build an integrated in silico and in vitro platform for discovering compounds that may bind to a chosen DNA target. Millions of compounds are initially screened in silico for selective binding to a particular structure and ranked to identify several hundred best hits. An important element of our strategy is the inclusion of an array of possible competing structures in the in silico screen. The best hundred or so hits are validated experimentally for binding to the actual target structure by a high-throughput 96-well thermal denaturation assay to yield the top ten candidates. Finally, these most promising candidates are thoroughly characterized for binding to their DNA target by rigorous biophysical methods, including isothermal titration calorimetry, differential scanning calorimetry, spectroscopy and competition dialysis.This platform was validated using quadruplex DNA as a target and a newly discovered quadruplex binding compound with possible anti-cancer activity was discovered. Some considerations when embarking on virtual screening and in silico experiments are also discussed.

Choline kinase-α expression and activity are increased in multiple human neoplasms as a result of growth factor stimulation and activation of cancer-related signaling pathways. The product of choline kinase-α, phosphocholine, serves as an essential metabolic reservoir for the production of phosphatidylcholine, the major phospholipid constituent of membranes and substrate for the production of lipid second messengers. Using in silico screening for small molecules that may interact with the choline kinase-α substrate binding domain, we identified a novel competitive inhibitor, N-(3,5-dimethylphenyl)-2-[[5-(4-ethylphenyl)-1H-1,2,4-triazol-3-yl]sulfanyl] acetamide (termed CK37) that inhibited purified recombinant human choline kinase-α activity, reduced the steady-state concentration of phosphocholine in transformed cells, and selectively suppressed the growth of neoplastic cells relative to normal epithelial cells. Choline kinase-α activity is required for the downstream production of phosphatidic acid, a promoter of several Ras signaling pathways. CK37 suppressed MAPK and PI3K/AKT signaling, disrupted actin cytoskeletal organization, and reduced plasma membrane ruffling. Finally, administration of CK37 significantly decreased tumor growth in a lung tumor xenograft mouse model, suppressed tumor phosphocholine, and diminished activating phosphorylations of ERK and AKT in vivo. Together, these results further validate choline kinase-α as a molecular target for the development of agents that interrupt Ras signaling pathways, and indicate that receptor-based computational screening should facilitate the identification of new classes of choline kinase-α inhibitors.

The in silico methods for drug discovery are becoming increasingly powerful and useful. That, in combination with increasing computer processor power, in our case using a novel distributed computing grid, has enabled us to greatly enhance our virtual screening efforts. Herein we review some of these efforts using both receptor and ligand-based virtual screening, with the goal of finding new anticancer agents. In particular, nucleic acids are a neglected set of targets, especially the different morphologies of duplex, triplex, and quadruplex DNA, many of which have increasing biological relevance. We also review examples of molecular modeling to understand receptors and using virtual screening against G-protein coupled receptor membrane proteins.

Macrophage migration inhibitory factor is a leaderless protein that is secreted from cells by a specialized, non-classical export pathway. The release of MIF nevertheless is regulated and its production in response to different inflammatory, mitogenic, and hormonal stimuli plays an important role in diverse physiologic and pathologic processes. We report herein the identification of the Golgi complex-associated protein, p115, as an intracellular binding partner for MIF. MIF interacts with p115 in the cytoplasm and the stimulated secretion of MIF results in the accumulation of both proteins in supernatants, which is consistent with MIF release from cells in conjunction with p115. The depletion of p115 from monocytes/macrophages decreases the release of MIF but not other cytokines following inflammatory stimulation or intracellular bacterial infection. Notably, the small molecule MIF inhibitor, 4-iodo-6-phenylpyrimidine, inhibits MIF secretion by targeting the interaction between MIF and p115. These data reveal p115 to be a critical intermediary component in the regulated secretion of MIF from monocytes/macrophages.

Single-stranded guanosine-rich oligodeoxyribonucleotides (GROs) have a propensity to form quadruplex structures that are stabilized by G-quartets. In addition to intense speculation about the role of G-quartet formation in vivo, there is considerable interest in the therapeutic potential of quadruplex oligonucleotides as aptamers or non-antisense antiproliferative agents. We previously have described several GROs that inhibit proliferation and induce apoptosis in cancer cell lines. The activity of these GROs was related to their ability to bind to a specific cellular protein (GRO-binding protein, which has been tentatively identified as nucleolin). In this report, we describe the physical properties and biological activity of a group of 12 quadruplex oligonucleotides whose structures have been characterized previously. This group includes the thrombin-binding aptamer, an anti-HIV oligonucleotide, and several quadruplexes derived from telomere sequences. Thermal denaturation and circular dichroism (CD) spectropolarimetry were utilized to investigate the stability, reversibility and ion dependence of G-quartet formation. The ability of each oligonucleotide to inhibit the proliferation of cancer cells and to compete for binding to the GRO-binding protein was also examined. Our results confirm that G-quartet formation is essential for biological activity of GROs and show that, in some cases, quadruplex structures formed in the presence of potassium ions are significantly more active than those formed in the presence of sodium ions. However, not all quadruplex structures exhibit antiproliferative effects, and the most accurate factor in predicting biological activity was the ability to bind to the GRO-binding protein. Our data also indicate that the CD spectra of quadruplex oligonucleotides may be more complex than previously thought.

We describe the design, synthesis and biophysical characterization of a novel DNA construct in which a folded quadruplex structure is joined to a standard double helix. Circular dichroism, gel electrophoresis, three-dimensional UV melting and differential scanning calorimetry were all used to characterize the structure. Rigorous molecular dynamics simulations were used to build a plausible atomic-level structural model of the DNA construct. This novel DNA construct provides a model for the duplex–quadruplex junction region at the end of chromosomal DNA and offers a system for the study of structure-selective ligand binding.

Human telomeric DNA consists of tandem repeats of the DNA sequence d(GGGTTA). Oligodeoxynucletotide telomere models such as d[A(GGGTTA)3GGG] (Tel22) fold in a cation-dependent manner into quadruplex structures consisting of stacked G-quartets linked by d(TTA) loops. NMR has shown that in Na+ solutions Tel22 forms a ‘basket’ topology of four antiparallel strands; in contrast, Tel22 in K+ solutions consists of a mixture of unknown topologies. Our previous studies on the mechanism of folding of Tel22 and similar telomere analogs utilized changes in UV absorption between 270 and 325 nm that report primarily on G-quartet formation and stacking showed that quadruplex formation occurs within milliseconds upon mixing with an appropriate cation. In the current study, we assessed the dynamics and equilibria of folding of specific loops by using Tel22 derivatives in which the dA residues were serially substituted with the fluorescent reporter base, 2-aminopurine (2-AP). Tel22 folding induced by Na+ or K+ assessed by changes in 2-AP fluorescence consists of at least three kinetic steps with time constants spanning a range of ms to several hundred seconds. Na+-dependent equilibrium titrations of Tel22 folding could be approximated as a cooperative two-state process. In contrast, K+-dependent folding curves were biphasic, revealing that different conformational ensembles are present in 1 mM and 30 mM K+. This conclusion was confirmed by 1H NMR. Molecular dynamics simulations revealed a K+ binding pocket in Tel22 located near dA1 that is specific for the so-called hybrid-1 conformation in which strand 1 is in a parallel arrangement. The possible presence of this topologically specific binding site suggests that K+ may play an allosteric role in regulating telomere conformation and function by modulating quadruplex tertiary structure.

Certain guanine-rich (G-rich) DNA and RNA molecules can associate intermolecularly or intramolecularly to form four stranded or “quadruplex” structures, which have unusual biophysical and biological properties. Several synthetic G-rich quadruplex-forming oligodeoxynucleotides have recently been investigated as therapeutic agents for various human diseases. We refer to these biologically active G-rich oligonucleotides as aptamers because their activities arise from binding to protein targets via shape-specific recognition (analogous to antibody-antigen binding). As therapeutic agents, the G-rich aptamers may have some advantages over monoclonal antibodies and other oligonucleotide-based approaches. For example, quadruplex oligonucleotides are non-immunogenic, heat stable and they have increased resistance to serum nucleases and enhanced cellular uptake compared to unstructured sequences. In this review, we describe the characteristics and activities of G-rich oligonucleotides. We also give a personal perspective on the discovery and development of AS1411, an antiproliferative G-rich phosphodiester oligonucleotide that is currently being tested as an anticancer agent in Phase II clinical trials. This molecule functions as an aptamer to nucleolin, a multifunctional protein that is highly expressed by cancer cells, both intracellularly and on the cell surface. Thus, the serendipitous discovery of the G-rich oligonucleotides also led to the identification of nucleolin as a new molecular target for cancer therapy.

The remarkable structural polymorphism of quadruplex-forming sequences has been a considerable impediment in the elucidation of quadruplex folds. Sequence modifications have commonly been used to perturb and purportedly select a particular form out of the ensemble of folds for nuclear magnetic resonance (NMR) or X-ray crystallographic analysis. Here we report a simple chromatographic technique that separates the individual folds without need for sequence modification. The sequence d(GGTGGTGGTGGTTGTGGTGGTGGTGG) forms a compact quadruplex according to a variety of common biophysical techniques. However, NMR and chromatography showed that this oligonucleotide produces at least eight monomeric quadruplex species that interconvert very slowly at room temperature. We have used a combination of spectroscopic, hydrodynamic and thermodynamic techniques to evaluate the physicochemical properties of the mixture and the individual species. These species have almost identical thermodynamic, hydrodynamic and electrophoretic properties, but significantly different NMR and circular dichroism (CD) spectra, as well as kinetic stability. These results demonstrate that simple standard low-resolution techniques cannot always be used for quadruplex fold determination or quality control purposes, and that simple thermodynamic analysis does not directly provide interpretable thermodynamic parameters.

SUMMARY

Severe congenital neutropenia (SCN) is characterized by a deficiency of mature neutrophils, leading to recurrent bacterial and fungal infections. Although mutations in Elastase-2, neutrophil (ELA2) predominate in human SCN, mutation of Ela2 in mice does not recapitulate SCN. The growth factor independent-1 (GFI1) transcription factor regulates ELA2. Mutations in GFI1 are associated with human SCN, and genetic deletion of Gfi1 results in murine neutropenia. We examined whether human SCN-associated GFI1N382S mutant proteins are causal in SCN and found that GFI1 functions as a rate-limiting granulopoietic molecular switch. The N382S mutation inhibited GFI1 DNA binding and resulted in a dominant-negative block to murine granulopoiesis. Moreover, Gfi1N382S selectively derepressed the monopoietic cytokine CSF1 and its receptor. Gfi1N382S–expressing Csf1−/− cells formed neutrophils. These results reveal a common transcriptional program that underlies both human and murine myelopoiesis, and that is central to the pathogenesis of SCN associated with mutations in GFI1. This shared transcriptional pathway may provide new avenues for understanding SCN caused by mutations in other genes and for clinical intervention into human neutropenias.

The structure of a higher-order G-quadruplex structure for human telomeric DNA is presented. The structure was determined by a novel integrated approach in which molecular dynamics simulations were used to produce a stable structure, from which specific experimentally accessible properties were predicted. These properties were tested by sedimentation velocity and steady-state fluorescence measurements. The structure that emerges is a dimeric structure with two quadruplex units, each with a different structure. The interface between the quadruplex units is stabilized by specific stacking interactions of loop nucleotides. The interface is a unique structure and a unique target for drug design.

The molecular docking tools Autodock and Surflex accurately reproduce the crystallographic structures of a collection of small molecule ligands that have been shown to bind nucleic acids. Docking studies were performed with the intercalators daunorubicin and ellipticine and the minor groove binders distamycin and pentamidine. Autodock and Surflex dock daunorubicin and distamycin to their nucleic acid targets within a resolution of approximately 2 Å, which is similar to the limit of the crystal structure resolution. However, for the top ranked poses, Autodock and Surflex both dock ellipticine into the correct site but in a different orientation compared to the crystal structure. This appears not only to be partly related to the symmetry of the target nucleic acid, as ellipticine is able to dock from either side of the intercalation site, but also due to the shape of the ligand and docking accuracy. Surflex docks pentamidine in a symmetrically equivalent orientation relative to the crystal structure, while Autodock was able to dock this molecule in the original orientation. In the case of the Surflex docking of pentamidine, the initial rmsd is misleading, given the symmetrical structure of pentamidine. Importantly, the ranking functions of both of these programs are able to return a top pose within approximately 2 Å rmsd for daunorubicin, distamycin, and pentamidine and approximately 3 Å rmsd for ellipticine compared to their respective crystal structures. Some docking challenges and potential pitfalls are explored, such as the importance of hydrogen treatment on ligands as well as the scoring functions of Autodock and Surflex. Overall for this set of complexes, Surflex is preferred over Autodock for virtual screening, as although the results are comparable, Surflex has significantly faster performance and ease of use under the optimal software conditions tested. These experiments show that molecular docking techniques can be successfully extended to include nucleic acid targets, a finding which has important implications for virtual screening applications and in the design of new small molecules to target therapeutically relevant morphologies of nucleic acids.

While chemokine and growth factor receptors are attractive and popular targets for cancer therapeutic intervention, structure-based targeting of the ligands themselves is generally not considered practical. New evidence indicates that a notable exception to this is macrophage migration inhibitory factor (MIF). MIF, an autocrine and paracrine acting cytokine/growth factor, plays a pivotal role in both the initiation and maintenance of neoplastic diseases. MIF possesses a non-physiologic enzymatic activity that is evolutionarily well conserved. Although small molecule antagonists of MIFs enzymatic active site have been reported to inhibit biological activities of MIF, universally high IC50s have limited their clinical appeal. Using a computational virtual screening strategy we have identified a unique small molecule inhibitor that serves as a suicide substrate for MIF resulting in the covalent modification of the catalytically active N-terminal proline. Our studies further reveal that this compound, 4-iodo-6-phenylpyrimidine (4-IPP), is ~ 5–10x times more potent in blocking MIF-dependent catalysis and lung adenocarcinoma cell migration and anchorage independent growth than the prototypical MIF inhibitor, ISO-1. Finally, using an in silico combinatorial optimization strategy, we have identified four unique congeners of 4-IPP that exhibit MIF inhibitory activity at concentrations 10–20x lower than that of parental 4-IPP.

Virtual Screening is an increasingly attractive way to discover new small molecules with potential medicinal value. We introduce a novel strategy that integrates use of the molecular docking software Surflex with experimental validation by the method of competition dialysis. This integrated approach was used to identify ligands that selectively bind to the triplex DNA poly(dA)-[poly(dT)]2. A library containing ∼2 million ligands was virtually screened to identify compounds with chemical and structural similarity to a known triplex intercalator, the napthylquinoline MHQ-12. Further molecular docking studies using compounds with high structural similarity resulted in two compounds that were then demonstrated by competition dialysis to have a superior affinity and selectivity for the triplex nucleic acid than MHQ-12. One of the compounds has a different chemical backbone than MHQ-12, which demonstrates the ability of this strategy to ‘scaffold hop’ and to identify small molecules with novel binding properties. Biophysical characterization of these compounds by circular dichroism and thermal denaturation studies confirmed their binding mode and selectivity. These studies provide a proof-of-principle for our integrated screening strategy, and suggest that this platform may be extended to discover new compounds that target therapeutically relevant nucleic acid morphologies.

A new bisintercalating anthracycline (WP762) has been designed, in which monomeric units of daunorubicin have been linked through their amino groups on the daunosamine moieties using an m-xylenyl linker. Differential scanning calorimetry and UV melting experiments were used to measure the ultratight binding of WP762 to DNA. The binding constant for the interaction of WP762 with herring sperm DNA was determined to be 7.3 (±0.2) × 1012 M−1 at 20°C. The large favorable binding free energy of −17.3 kcal mol−1 was found to result from a large negative enthalpic contribution of −33.8 kcal mol−1 and an opposing entropic term (−TΔS = +16.5 kcal mol−1). A comparative molecular modeling study rationalized the increased binding by the m-xylenyl linker of WP762 positioning in the DNA minor groove compared to the p-xylenyl linker found in WP631, the first bis-anthracycline of this type. The cytotoxicity of WP762 was compared to that of other anthracyclines in Jurkat T lymphocytes. These studies, together with an analysis of the cell-cycle traverse in the presence of WP762, suggest that in these cells the new drug is more cytotoxic than the structurally related WP631.

In this review, we give an overview of recent literature on the structure and stability of unimolecular G-rich quadruplex structures that are relevant to drug design and for in vivo function. The unifying theme in this review is energetics. The thermodynamic stability of quadruplexes has not been studied in the same detail as DNA and RNA duplexes, and there are important differences in the balance of forces between these classes of folded oligonucleotides. We provide an overview of the principles of stability and where available the experimental data that report on these principles. Significant gaps in the literature have been identified, that should be filled by a systematic study of well-defined quadruplexes not only to provide the basic understanding of stability both for design purposes, but also as it relates to in vivo occurrence of quadruplexes. Techniques that are commonly applied to the determination of the structure, stability and folding are discussed in terms of information content and limitations. Quadruplex structures fold and unfold comparatively slowly, and DNA unwinding events associated with transcription and replication may be operating far from equilibrium. The kinetics of formation and resolution of quadruplexes, and methodologies are discussed in the context of stability and their possible biological occurrence.

Arylamine N-acetyltransferase 2 (NAT2) modifies drug efficacy/toxicity and cancer risk due to its role in bioactivation and detoxification of arylamine and hydrazine drugs and carcinogens. Human NAT2 alleles possess a combination of single nucleotide polymorphisms (SNPs) associated with slow acetylation phenotypes. Clinical and molecular epidemiology studies investigating associations of NAT2 genotype with drug efficacy/toxicity and/or cancer risk are compromised by incomplete and sometimes conflicting information regarding genotype/phenotype relationships. Studies in our laboratory and others have characterized the functional effects of SNPs alone, and in combinations present in alleles or haplotypes. We extrapolate this data generated following recombinant expression in yeast and COS-1 cells to assist in the interpretation of NAT2 structure. Whereas previous structural studies used homology models based on templates of N-acetyltransferase enzyme crystal structures from various prokaryotic species, alignment scores between bacterial and mammalian N-acetyltransferase protein sequences are low (~ 30%) with important differences between the bacterial and mammalian protein structures. Recently, the crystal structure of human NAT2 was released from the Protein Data Bank under accession number 2PFR. We utilized the NAT2 crystal structure to evaluate the functional effects of SNPs resulting in the protein substitutions R64Q (G191A), R64W (C190T), I114T (T341C), D122N (G364A), L137F (A411T), Q145P (A434C), E167K (G499A), R197Q (C590A), K268R (A803G), K282T (A845C), and G286E (G857A) of NAT2. This analysis advances understanding of NAT2 structure-function relationships, important for interpreting the role of NAT2 genetic polymorphisms in bioactivation and detoxification of arylamine and hydrazine drugs and carcinogens.

Arylamine N-acetyltransferases (NATs) play an important role in the metabolism of arylamine and hydrazine drugs and many arylamine pro-carcinogens. The two human N-acetyltransferases, NAT1 and NAT2, are widely distributed in human tissues and highly polymorphic. While many xenobiotic procarcinogens and drugs are known mammalian NAT substrates, it is unclear what physiological roles these enzymes might play, what endogenous substrates they primarily act upon, or the mechanisms underlying the functional effects of specific NAT gene coding region SNPs. Analyses of mammalian NAT protein structures can greatly help to answer these questions. Homology modeling techniques can be used to approximate mammalian NAT structures using known bacterial NAT crystal structures as templates. In comparison to the bacterial template NATs used for homology modeling, mammalian NATs have a 17 residue insert of unknown structure and function. Homology modeling analyses yielded two different alignments (Modeller 8v1 or 3DCoffee algorithms) that placed this insert in two likely alternative locations. Secondary structure prediction techniques and experimental analyses of a series of human NAT2 mutants with artificial deletions/replacements of the insert region distinguished one of these alternatives as the most likely insert location and provided a better understanding of its structure and function. This study demonstrates both the utility and limitations of computational structural modeling with proteins that differ as much as the mammalian and bacterial NATs.

Human N-acetyltransferase 1 (NAT1) alleles are characterized by one or more single nucleotide polymorphisms (SNPs) associated with rapid and slow acetylation phenotypes. NAT1 both activates and deactivates arylamine drugs and carcinogens, and NAT1 polymorphisms are associated with increased frequencies of many cancers and birth defects. The recently resolved human NAT1 crystal structure was used to evaluate SNPs resulting in the protein substitutions R64W, V149I, R187Q, M205V, S214A, D251V, E261K, and I263V. The analysis enhances knowledge of NAT1 structure-function relationships, important for understanding associations of NAT1 SNPs with genetic predisposition to cancer, birth defects, and other diseases.

The structure of human telomere DNA is of intense interest because of its role in the biology of both cancer and aging. The sequence [5′-AGGG(TTAGGG)3] has been used as a model for telomere DNA in both NMR and X-ray crystallographic studies, the results of which show dramatically different structures. In Na+ solution, NMR revealed an antiparallel G-quadruplex structure that featured both diagonal and lateral TTA loops. Crystallographic studies in the presence of K+ revealed a flattened, propeller-shaped structure featuring a parallel-stranded G-quadruplex with symmetrical external TTA loops. We report the results of biophysical experiments in solution and computational studies that are inconsistent with the reported crystal structure, indicating that a different structure exists in K+ solutions. Sedimentation coefficients were determined experimentally in both Na+ and K+ solutions and were compared with values calculated using bead models for the reported NMR and crystal structures. Although the solution NMR structure accurately predicted the observed S-value in Na+ solution, the crystal structure predicted an S-value that differed dramatically from that experimentally observed in K+ solution. The environments of loop adenines were probed by quantitative fluorescence studies using strategic and systematic single-substitutions of 2-aminopurine for adenine bases. Both fluorescence intensity and quenching experiments in K+ yielded results at odds with quantitative predictions from the reported crystal structure. Circular dichroism and fluorescence quenching studies in the presence of the crowding agent polyethylene glycol showed dramatic changes in the quadruplex structure in K+ solutions, but not in Na+ solutions, suggesting that the crystal environment may have selected for a particular conformational form. Molecular dynamics simulations were performed to yield model structures for the K+ quadruplex form that are consistent with our biophysical results and with previously reported chemical modification studies. These models suggest that the biologically relevant structure of the human telomere quadruplex in K+ solution is not the one determined in the published crystalline state.

The Rps0 proteins of Saccharomyces cerevisiae are components of the 40S ribosomal subunit required for maturation of the 3′ end of 18S rRNA. Drosophila and human homologs of the Rps0 proteins physically interact with Rps21 proteins, and decreased expression of both proteins in Drosophila impairs control of cellular proliferation in hematopoietic organs during larval development. Here, we characterize the yeast RPS21A/B genes and show that strains where both genes are disrupted are not viable. Relative to the wild type, cells with disrupted RPS21A or RPS21B genes exhibit a reduction in growth rate, a decrease in free 40S subunits, an increase in the amount of free 60S subunits, and a decrease in polysome size. Ribosomal RNA processing studies reveal RPS21 and RPS0 mutants have virtually identical processing defects. The pattern of processing defects observed in RPS0 and RPS21 mutants is not a general characteristic of strains with suboptimal levels of small subunit ribosomal proteins, since disruption of the RPS18A or RPS18B genes results in related but distinct processing defects. Together, these data link the Rps0 and Rps21 proteins together functionally in promoting maturation of the 3′ end of 18S rRNA and formation of active 40S ribosomal subunits.