Expression and purification of imidazole glycerol phosphate synthase from Saccharomyces cerevisiae

MDiGest: A Python package for describing allostery from molecular dynamics simulations

Many biological processes are regulated by allosteric mechanisms that communicate with distant sites in the protein responsible for functionality. The binding of a small molecule at an allosteric site typically induces conformational changes that propagate through the protein along allosteric pathways regulating enzymatic activity. Elucidating those communication pathways from allosteric sites to orthosteric sites is, therefore, essential to gain insights into biochemical processes. Targeting the allosteric pathways by mutagenesis can allow the engineering of proteins with desired functions. Furthermore, binding small molecule modulators along the allosteric pathways is a viable approach to target reactions using allosteric inhibitors/activators with temporal and spatial selectivity. Methods based on network theory can elucidate protein communication networks through the analysis of pairwise correlations observed in molecular dynamics (MD) simulations using molecular descriptors that serve as proxies for allosteric information. Typically, single atomic descriptors such as α-carbon displacements are used as proxies for allosteric information. Therefore, allosteric networks are based on correlations revealed by that descriptor. Here, we introduce a Python software package that provides a comprehensive toolkit for studying allostery from MD simulations of biochemical systems. MDiGest offers the ability to describe protein dynamics by combining different approaches, such as correlations of atomic displacements or dihedral angles, as well as a novel approach based on the correlation of Kabsch-Sander electrostatic couplings. MDiGest allows for comparisons of networks and community structures that capture physical information relevant to allostery. Multiple complementary tools for studying essential dynamics include principal component analysis, root mean square fluctuation, as well as secondary structure-based analyses.

Influence of tetraconazole on the proteome profile of Saccharomyces cerevisiae Lalvin T73™ strain

This work aimed to evaluate the modifications on the proteome profile of Saccharomyces cerevisiae T73™ strain as a consequence of its adaptive response to the presence of tetraconazole molecules in the fermentation medium. Pasteurised grape juices were separately supplemented with tetraconazole or a commercial formulation containing 12.5% w/v of tetraconazole at two concentration levels. In addition, experiments without fungicides were developed for comparative purposes. Proteome profiles of yeasts cultured in the presence or absence of fungicide molecules were different. Independently of the fungicide treatment applied, the highest variations concerning the control sample were observed for those proteins involved in metabolic processes, especially in the metabolism of nitrogen compounds. Tetraconazole molecules altered the abundance of several enzymes involved in the biosynthesis of amino acids, purines, and ergosterol. Moreover, differences in the abundance of several enzymes of the TCA cycle were found. Changes observed were different between the active substance and the commercial formulation. SIGNIFICANCE: The presence of fungicide residues in grape juice has direct implications on the development of the aromatic profile of the wine. These alterations could be related to changes in the secondary metabolism of yeasts. However, the molecular mechanisms involved in the response of yeasts to fungicide residues remains quite unexplored. Through this exhaustive proteomic study, alterations in the amino acids biosynthesis pathways due to the presence of the tetraconazole molecules were observed. Amino acids are precursors of some important higher alcohols and ethyl acetates (such as methionol, 2-phenylethanol, isoamyl alcohol or 2-phenylacetate). Besides, the effect of tetraconazole on the ergosterol biosynthesis pathway could be related to a higher production of medium-chain fatty acids and their corresponding ethyl acetates.

Combining ancestral sequence reconstruction with protein design to identify an interface hotspot in a key metabolic enzyme complex

It is important to identify hotspot residues that determine protein-protein interactions in interfaces of macromolecular complexes. We have applied a combination of ancestral sequence reconstruction and protein design to identify hotspots within imidazole glycerol phosphate synthase (ImGPS). ImGPS is a key metabolic enzyme complex, which links histidine and de novo purine biosynthesis and consists of the cyclase subunit HisF and the glutaminase subunit HisH. Initial fluorescence titration experiments showed that HisH from Zymomonas mobilis (zmHisH) binds with high affinity to the reconstructed HisF from the last universal common ancestor (LUCA-HisF) but not to HisF from Pyrobaculum arsenaticum (paHisF), which differ by 103 residues. Subsequent titration experiments with a reconstructed evolutionary intermediate linking LUCA-HisF and paHisF and inspection of the subunit interface of a contemporary ImGPS allowed us to narrow down the differences crucial for zmHisH binding to nine amino acids of HisF. Homology modeling and in silico mutagenesis studies suggested that at most two of these nine HisF residues are crucial for zmHisH binding. These computational results were verified by experimental site-directed mutagenesis, which finally enabled us to pinpoint a single amino acid residue in HisF that is decisive for high-affinity binding of zmHisH. Our work shows that the identification of protein interface hotspots can be very efficient when reconstructed proteins with different binding properties are included in the analysis. Proteins 2017; 85:312-321. © 2016 Wiley Periodicals, Inc.

Allosteric Communication Disrupted by a Small Molecule Binding to the Imidazole Glycerol Phosphate Synthase Protein-Protein Interface

Allosteric enzymes regulate a wide range of catalytic transformations, including biosynthetic mechanisms of important human pathogens, upon binding of substrate molecules to an orthosteric (or active) site and effector ligands at distant (allosteric) sites. We find that enzymatic activity can be impaired by small molecules that bind along the allosteric pathway connecting the orthosteric and allosteric sites, without competing with endogenous ligands. Noncompetitive allosteric inhibitors disrupted allostery in the imidazole glycerol phosphate synthase (IGPS) enzyme from Thermotoga maritima as evidenced by nuclear magnetic resonance, microsecond time-scale molecular dynamics simulations, isothermal titration calorimetry, and kinetic assays. The findings are particularly relevant for the development of allosteric antibiotics, herbicides, and antifungal compounds because IGPS is absent in mammals but provides an entry point to fundamental biosynthetic pathways in plants, fungi, and bacteria.

Millisecond dynamics in the allosteric enzyme imidazole glycerol phosphate synthase (IGPS) from Thermotoga maritima

IGPS is a 51 kDa heterodimeric enzyme comprised of two proteins, HisH and HisF, that catalyze the hydrolysis of glutamine to produce NH(3) in the HisH active site and the cyclization of ammonia with N'-[(5'-phosphoribulosyl)formimino]-5-aminoimidazole-4-carboxamide-ribonucleotide (PRFAR) in HisF to produce imidazole glycerol phosphate (IGP) and 5-aminoimidazole-4-carboxamide ribotide (AICAR). Binding of PRFAR and IGP stimulates glutaminase activity in the HisH enzyme over 5,000 and 100-fold, respectively, despite the active sites being >25 A apart. The details of this long-range protein communication process were investigated by solution NMR spectroscopy and CPMG relaxation dispersion experiments. Formation of the heterodimer enzyme results in a reduction in millisecond motions in HisF that extend throughout the protein. Binding of lGP results in an increase in protein-wide millisecond dynamics evidenced as severe NMR line broadening and elevated R (ex) values. Together, these data demonstrate a grouping of flexible residues that link the HisF active site with the protein interface to which HisH binds and provide a model for the path of communication between the IGPS active sites.

Two independent routes of de novo vitamin B6 biosynthesis: not that different after all

Vitamin B6 is well known in its biochemically active form as pyridoxal 5'-phosphate, an essential cofactor of numerous metabolic enzymes. The vitamin is also implicated in numerous human body functions ranging from modulation of hormone function to its recent discovery as a potent antioxidant. Its de novo biosynthesis occurs only in bacteria, fungi and plants, making it an essential nutrient in the human diet. Despite its paramount importance, its biosynthesis was predominantly investigated in Escherichia coli, where it is synthesized from the condensation of deoxyxylulose 5-phosphate and 4-phosphohydroxy-L-threonine catalysed by the concerted action of PdxA and PdxJ. However, it has now become clear that the majority of organisms capable of producing this vitamin do so via a different route, involving precursors from glycolysis and the pentose phosphate pathway. This alternative pathway is characterized by the presence of two genes, Pdx1 and Pdx2. Their discovery has sparked renewed interest in vitamin B6, and numerous studies have been conducted over the last few years to characterize the new biosynthesis pathway. Indeed, enormous progress has been made in defining the nature of the enzymes involved in both pathways, and important insights have been provided into their mechanisms of action. In the present review, we summarize the recent advances in our knowledge of the biosynthesis of this versatile molecule and compare the two independent routes to the biosynthesis of vitamin B6. Surprisingly, this comparison reveals that the key biosynthetic enzymes of both pathways are, in fact, very similar both structurally and mechanistically.

Insight into gene fusion from molecular dynamics simulation of fused and un-fused IGPS (Imidazole Glycerol Phosphate Synthetase)

Gene fusion produces proteins with novel structural architectures during evolution. Recent comparative genome analysis shows several cases of fusion/fission across distant phylogeny. However, the selection forces driving gene fusion are not fully understood due to the lack of structural, dynamics and kinetics data. Available structural data at PDB (protein databank) contains limited cases of structural pairs describing fused and un-fused structures. Nonetheless, we identified a pair of IGPS (imidazole glycerol phosphate synthetase) structures (comprising of HisF - glutaminase unit and HisH - cyclase unit) from S. cerevisiae (SC) and T. thermophilus (TT). The HisF-HisH structural units are domains in SC and subunits in TT. Hence, they are fused in SC and un-fused in TT. Subsequently, a domain-domain interface is formed in SC and a subunit-subunit interface in TT between HisF and HisH. Our interest is to document the structure and dynamics differences between fused and un-fused IGPS. Therefore, we probed into the structures of fused IGPS in SC and un-fused IGPS in TT using molecular dynamics simulation for 5ns. Simulation shows that fused IGPS in SC has larger interface area between HisF-HisH and greater radius of gyration compared to un-fused IGPS in TT. These structural features for the first time demonstrate the evolutionary advantage in generating proteins with novel structural architecture through gene fusion.

Vitamin B6 biosynthesis by the malaria parasite Plasmodium falciparum: biochemical and structural insights

Vitamin B6 is one of nature's most versatile cofactors. Most organisms synthesize vitamin B6 via a recently discovered pathway employing the proteins Pdx1 and Pdx2. Here we present an in-depth characterization of the respective orthologs from the malaria parasite, Plasmodium falciparum. Expression profiling of Pdx1 and -2 shows that blood-stage parasites indeed possess a functional vitamin B6 de novo biosynthesis. Recombinant Pdx1 and Pdx2 form a complex that functions as a glutamine amidotransferase with Pdx2 as the glutaminase and Pdx1 as pyridoxal-5 '-phosphate synthase domain. Complex formation is required for catalytic activity of either domain. Pdx1 forms a chimeric bi-enzyme with the bacterial YaaE, a Pdx2 ortholog, both in vivo and in vitro, although this chimera does not attain full catalytic activity, emphasizing that species-specific structural features govern the interaction between the protein partners of the PLP synthase complexes in different organisms. To gain insight into the activation mechanism of the parasite bi-enzyme complex, the three-dimensional structure of Pdx2 was determined at 1.62 A. The obstruction of the oxyanion hole indicates that Pdx2 is in a resting state and that activation occurs upon Pdx1-Pdx2 complex formation.

Structural elements in IGP synthase exclude water to optimize ammonia transfer

In the complex pathway of histidine biosynthesis, a key branch point linking amino acid and purine biosynthesis is catalyzed by the bifunctional enzyme imidazole glycerol phosphate (IGP) synthase. The first domain of IGP synthase, a triad glutamine amidotransferase, hydrolyzes glutamine to form glutamate and ammonia. Its activity is tightly regulated by the binding of the substrate PRFAR to its partner synthase domain. Recent crystal structures and molecular dynamics simulations strongly suggest that the synthase domain, a (beta/alpha)(8) barrel protein, mediates the insertion of ammonia and ring formation in IGP by channeling ammonia from one remote active site to the other. Here, we combine both mutagenesis experiments and computational investigations to gain insight into the transfer of ammonia and the mechanism of conduction. We discover an alternate route for the entrance of ammonia into the (beta/alpha)(8) barrel and argue that water acts as both agonist and antagonist to the enzymatic function. Our results indicate that the architecture of the two subdomains, most notably the strict conservation of key residues at the interface and within the (beta/alpha)(8) barrel, has been optimized to allow the efficient passage of ammonia, and not water, between the two remote active sites.

Structural evidence for ammonia tunneling across the (beta alpha)(8) barrel of the imidazole glycerol phosphate synthase bienzyme complex

Since reactive ammonia is not available under physiological conditions, glutamine is used as a source for the incorporation of nitrogen in a number of metabolic pathway intermediates. The heterodimeric ImGP synthase that links histidine and purine biosynthesis belongs to the family of glutamine amidotransferases in which the glutaminase activity is coupled with a subsequent synthase activity specific for each member of the enzyme family. Its X-ray structure from the hyperthermophile Thermotoga maritima shows that the glutaminase subunit is associated with the N-terminal face of the (beta alpha)(8) barrel cyclase subunit. The complex reveals a putative tunnel for the transfer of ammonia over a distance of 25 A. Although ammonia tunneling has been reported for glutamine amidotransferases, the ImGP synthase has evolved a novel mechanism, which extends the known functional properties of the versatile (beta alpha)(8) barrel fold.

Imidazole glycerol phosphate synthase from Thermotoga maritima. Quaternary structure, steady-state kinetics, and reaction mechanism of the bienzyme complex

Imidazole glycerol phosphate synthase, which links histidine and de novo purine biosynthesis, is a member of the glutamine amidotransferase family. In bacteria, imidazole glycerol phosphate synthase constitutes a bienzyme complex of the glutaminase subunit HisH and the synthase subunit HisF. Nascent ammonia produced by HisH reacts at the active site of HisF with N'-((5'-phosphoribulosyl)formimino)-5-aminoimidazole-4-carboxamide-ribonucleotide to yield the products imidazole glycerol phosphate and 5-aminoimidazole-4-carboxamide ribotide. In order to elucidate the interactions between HisH and HisF and the catalytic mechanism of the HisF reaction, the enzymes tHisH and tHisF from Thermotoga maritima were produced in Escherichia coli, purified, and characterized. Isolated tHisH showed no detectable glutaminase activity but was stimulated by complex formation with tHisF to which either the product imidazole glycerol phosphate or a substrate analogue were bound. Eight conserved amino acids at the putative active site of tHisF were exchanged by site-directed mutagenesis, and the purified variants were investigated by steady-state kinetics. Aspartate 11 appeared to be essential for the synthase activity both in vitro and in vivo, and aspartate 130 could be partially replaced only by glutamate. The carboxylate groups of these residues could provide general acid/base catalysis in the proposed catalytic mechanism of the synthase reaction.

Subunit interactions and glutamine utilization by Escherichia coli imidazole glycerol phosphate synthase

A selection strategy has been developed to identify amino acid residues involved in subunit interactions that coordinate the two half-reactions catalyzed by glutamine amidotransferases. The protein structures known for this class of enzymes have revealed that ammonia is shuttled over long distances and that each amidotransferase evolved different molecular tunnels for this purpose. The heterodimeric Escherichia coli imidazole glycerol phosphate (IGP) synthase was probed to assess if residues in the substrate amination subunit (HisF) are critical for the glutaminase activity in the HisH subunit. The activity of the HisH subunit is dependent upon binding of the nucleotide substrate at the HisF active site. This regulatory function has been exploited as a biochemical selection of mutant HisF subunits that retain full activity with ammonia as a substrate but, when constituted as a holoenzyme with wild-type HisH, impair the glutamine-dependent activity of IGP synthase. The steady-state kinetic constants for these IGP synthases with HisF alleles showed three distinct effects depending upon the site of mutation. For example, mutation of the R5 residue has similar effects on the glutamine-dependent amidotransfer reaction; however, k(cat)/K(m) for the glutaminase half-reaction was increased 10-fold over that for the wild-type enzyme with nucleotide substrate. This site appears essential for coupling of the glutamine hydrolysis and ammonia transfer steps and is the first example of a site remote to the catalytic triad that modulates the process. The results are discussed in the context of recent X-ray crystal structures of glutamine amidotransferases that relate the glutamine binding and acceptor binding sites.