A vibrational structure of 7,8-dihydrobiopterin bound to dihydroneopterin aldolase

Bioenergetic mapping of 'healthy microbiomes' via compound processing potential imprinted in gut and soil metagenomes

Despite mounting evidence of their importance in human health and ecosystem functioning, the definition and measurement of 'healthy microbiomes' remain unclear. More advanced knowledge exists on health associations for compounds used or produced by microbes. Environmental microbiome exposures (especially via soils) also help shape, and may supplement, the functional capacity of human microbiomes. Given the synchronous interaction between microbes, their feedstocks, and micro-environments, with functional genes facilitating chemical transformations, our objective was to examine microbiomes in terms of their capacity to process compounds relevant to human health. Here we integrate functional genomics and biochemistry frameworks to derive new quantitative measures of in silico potential for human gut and environmental soil metagenomes to process a panel of major compound classes (e.g., lipids, carbohydrates) and selected biomolecules (e.g., vitamins, short-chain fatty acids) linked to human health. Metagenome functional potential profile data were translated into a universal compound mapping 'landscape' based on bioenergetic van Krevelen mapping of function-level meta-compounds and corresponding functional relative abundances, reflecting imprinted genetic capacity of microbiomes to metabolize an array of different compounds. We show that measures of 'compound processing potential' associated with human health and disease (examining atherosclerotic cardiovascular disease, colorectal cancer, type 2 diabetes and anxious-depressive behavior case studies), and displayed seemingly predictable shifts along gradients of ecological disturbance in plant-soil ecosystems (three case studies). Ecosystem quality explained 60-92 % of variation in soil metagenome compound processing potential measures in a post-mining restoration case study dataset. With growing knowledge of the varying proficiency of environmental microbiota to process human health associated compounds, we might design environmental interventions or nature prescriptions to modulate our exposures, thereby advancing microbiota-oriented approaches to human health. Compound processing potential offers a simplified, integrative approach for applying metagenomics in ongoing efforts to understand and quantify the role of microbiota in environmental- and human-health.

Solution structures of purine base analogues 9-deazaguanine and 9-deazahypoxanthine

Deaza analogues of nucleobases are potential drugs against infectious diseases caused by parasites. A caveat is that apart from binding their target parasite enzymes, they also bind and inhibit enzymes of the host. In order to design derivatives of deaza analogues which specifically bind target enzymes, knowledge of their molecular structure, protonation state, and predominant tautomers at physiological conditions is essential. We have employed resonance Raman spectroscopy at an excitation wavelength of 260 nm, to decipher solution structure of 9-deazaguanine (9DAG) and 9-deazahypoxanthine (9DAH). These are analogues of guanine and hypoxanthine, respectively, and have been exploited to study static complexes of nucleobase binding enzymes. Such enzymes are known to perturb pKa of their ligands, and thus, we also determined solution structures of these analogues at two, acidic and alkaline, pH. Structure of each possible protonation state and tautomer was computed using density functional theoretical calculations. Species at various pHs were identified based on isotopic shifts in experimental wavenumbers and by comparing these shifts with corresponding computed isotopic shifts. Our results show that at physiological pH, N1 of pyrimidine ring in 9DAG and 9DAH bears a proton. At lower pH, N3 is place of protonation, and at higher pH, deprotonation occurs at N1 position. The proton at N7 of purine ring remains intact even at pH 12.5. We have further compared these results with naturally occurring nucleotides. Our results identify key vibrational modes which can report on hydrogen bonding interactions, protonation and deprotonation in purine rings upon binding to the active site of enzymes.

Biochemical and structural studies of 6-carboxy-5,6,7,8-tetrahydropterin synthase reveal the molecular basis of catalytic promiscuity within the tunnel-fold superfamily

6-Pyruvoyltetrahydropterin synthase (PTPS) homologs in both mammals and bacteria catalyze distinct reactions using the same 7,8-dihydroneopterin triphosphate substrate. The mammalian enzyme converts 7,8-dihydroneopterin triphosphate to 6-pyruvoyltetrahydropterin, whereas the bacterial enzyme catalyzes the formation of 6-carboxy-5,6,7,8-tetrahydropterin. To understand the basis for the differential activities we determined the crystal structure of a bacterial PTPS homolog in the presence and absence of various ligands. Comparison to mammalian structures revealed that although the active sites are nearly structurally identical, the bacterial enzyme houses a His/Asp dyad that is absent from the mammalian protein. Steady state and time-resolved kinetic analysis of the reaction catalyzed by the bacterial homolog revealed that these residues are responsible for the catalytic divergence. This study demonstrates how small variations in the active site can lead to the emergence of new functions in existing protein folds.

Vibrational spectroscopic analyses of unique yellow feather pigments (spheniscins) in penguins

Many animals extract, synthesize and refine chemicals for colour display, where a range of compounds and structures can produce a diverse colour palette. Feather colours, for example, span the visible spectrum and mostly result from pigments in five chemical classes (carotenoids, melanins, porphyrins, psittacofulvins and metal oxides). However, the pigment that generates the yellow colour of penguin feathers appears to represent a sixth, poorly characterized class of feather pigments. This pigment class, here termed 'spheniscin', is displayed by half of the living penguin genera; the larger and richer colour displays of the pigment are highly attractive. Using Raman and mid-infrared spectroscopies, we analysed yellow feathers from two penguin species (king penguin, Aptenodytes patagonicus; macaroni penguin, Eudyptes chrysolophus) to further characterize spheniscin pigments. The Raman spectrum of spheniscin is distinct from spectra of other feather pigments and exhibits 17 distinctive spectral bands between 300 and 1700 cm(-1). Spectral bands from the yellow pigment are assigned to aromatically bound carbon atoms, and to skeletal modes in an aromatic, heterocyclic ring. It has been suggested that the penguin pigment is a pterin compound; Raman spectra from yellow penguin feathers are broadly consistent with previously reported pterin spectra, although we have not matched it to any known compound. Raman spectroscopy can provide a rapid and non-destructive method for surveying the distribution of different classes of feather pigments in the avian family tree, and for correlating the chemistry of spheniscin with compounds analysed elsewhere. We suggest that the sixth class of feather pigments may have evolved in a stem-lineage penguin and endowed modern penguins with a costly plumage trait that appears to be chemically unique among birds.

Mechanism of dihydroneopterin aldolase

Dihydroneopterin aldolase (DHNA) catalyzes both the cleavage of 7,8-dihydro-D-neopterin (DHNP) to form 6-hydroxymethyl-7,8-dihydropterin (HP) and glycolaldehyde and the epimerization of DHNP to form 7,8-dihydro-L-monapterin (DHMP). Whether the epimerization reaction uses the same reaction intermediate as the aldol reaction or the deprotonation and reprotonation of C2' of DHNP has been investigated by NMR analysis of the reaction products in a D2O solvent. No deuteration of C2' was observed for the newly formed DHMP. This result strongly suggests that the epimerization reaction uses the same reaction intermediate as the aldol reaction. In contrast with an earlier observation, the DHNA-catalyzed reaction is reversible, which also supports a nonstereospecific retroaldol/aldol mechanism for the epimerization reaction. The binding and catalytic properties of DHNAs from both Staphylococcus aureus (SaDHNA) and Escherichia coli (EcDHNA) were determined by equilibrium binding and transient kinetic studies. A complete set of kinetic constants for both the aldol and epimerization reactions according to a unified kinetic mechanism was determined for both SaDHNA and EcDHNA. The results show that the two enzymes have significantly different binding and catalytic properties, in accordance with the significant sequence differences between them.

Mechanism of dihydroneopterin aldolase: a molecular dynamics study of the apo enzyme and its product complex

Dihydroneopterin aldolase (DHNA), an enzyme in the pathway that generates folic acid in bacteria, is investigated by a series of molecular dynamics simulations in its free form and complexed with its product, 6-hydroxymethyl-7,8-dihydropterin (HP). The active sites in DHNA are formed at the interface between pairs of protomers in this octameric protein. On the basis of root-mean-square deviation and root-mean-square fluctuation analyses of the trajectories, which take advantage of the presence of eight active sites, flexible regions of the apo protein surrounding the active site are identified and, upon binding HP, show that the active site is rigidified. Specific residues, associated with binding and the catalytic mechanism of DHNA, are associated with these flexible regions, and their interactions with HP account for most of the binding energy. A Principal Component Analysis shows rigidification of DHNA upon HP binding and that only a few modes of motion capture most of the atomic fluctuations in both apo and HP-bound forms. HP is pushed out of the active site in a series of simulations with different restrained positions between HP and DHNA to obtain a view of the exit pathway and energetic barrier to product release. The chosen pathway leads to a minimal disturbance of the system and provides a barrier consistent with the experimentally determined rate of product release. An analysis of the various components that contribute to the exit path energy and entropy provides insight into the energy-entropy compensation for product release.

Biosynthesis of Tetrahydrofolate in Plants: Crystal Structure of 7,8-Dihydroneopterin Aldolase from Arabidopsis thaliana Reveals a Novel Adolase Class

Dihydroneopterin aldolase (DHNA) catalyses a retroaldol reaction yielding 6-hydroxymethyl-7,8-dihydropterin, a biosynthetic precursor of the vitamin, tetrahydrofolate. The enzyme is a potential target for antimicrobial and anti-parasite chemotherapy. A gene specifying a dihydroneopterin aldolase from Arabidopsis thaliana was expressed in a recombinant Escherichia coli strain. The recombinant protein was purified to apparent homogeneity and crystallised using polyethylenglycol as the precipitating agent. The crystal structure was solved by X-ray diffraction analysis at 2.2A resolution. The enzyme forms a D(4)-symmetric homooctamer. Each polypeptide chain is folded into a single domain comprising an antiparallel four-stranded beta-sheet and two long alpha-helices. Four monomers are arranged in a tetrameric ring, and two of these rings form a hollow cylinder. Well defined purine derivatives are found at all eight topologically equivalent active sites. The subunit fold of the enzyme is related to substructures of dihydroneopterin triphosphate epimerase, GTP cyclohydrolase I, and pyruvoyltetrahydropterin synthase, which are all involved in the biosynthesis of pteridine type cofactors, and to urate oxidase, although some members of that superfamily have no detectable sequence similarity. Due to structural and mechanistical differences of DHNA in comparison with class I and class II aldolases, a new aldolase class is proposed.

The folic acid biosynthesis pathway in bacteria: evaluation of potential for antibacterial drug discovery

The potential of the folic acid biosynthesis pathway as a target for the development of antibiotics has been acknowledged for many years and validated by the clinical use of several drugs. Recently, the crystal structures of all but one of the enzymes in the pathway from GTP to dihydrofolate have been determined. Given that structure-based drug design strategies are now widely employed, these recent developments have prompted a re-evaluation of the potential of each of the enzymes in the pathway as a target for development of specific inhibitors. Here, we review the current knowledge of the structure and mechanism of each enzyme in the bacterial folic acid biosynthesis pathway from GTP to dihydrofolate and draw conclusions regarding the potential of each enzyme as a target for therapeutic intervention.

Vibrational structure of dihydrofolate bound to R67 dihydrofolate reductase

R67 is a Type II dihydrofolate reductase (DHFR) that catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate by facilitating the addition of a proton to N5 of DHF and the transfer of a hydride ion from NADPH to C6. Because this enzyme is a plasmid-encoded DHFR from trimethoprim-resistant bacteria, extensive studies on R67 with various methods have been performed to elucidate its reaction mechanism. Here, Raman difference measurements, conducted on the ternary complex of R67.NADP(+).DHF believed to be an accurate mimic of the productive DHFR.NADPH.DHF complex, show that the pK(a) of N5 in the complex is less than 4. This is in clear contrast to the behavior observed in Escherichia coli DHFR, a substantially more efficient enzyme, where the pK(a) of bound DHF at N5 is increased to 6.5 compared with its solution value of 2.6. A comparison of the ternary complexes in R67 and E. coli DHFRs suggests that enzymic raising of the pK(a) at N5 can significantly increase the catalytic efficiency of the hydride transfer step. However, R67 shows that even without such a strategy an effective DHFR can still be designed.