Nucleotide sequences of genes encoding heat-stable and heat-labile glyceraldehyde-3-phosphate dehydrogenases; amino acid sequence and protein thermostability

Ornithine carbamoyltransferase from the extreme thermophile Thermus thermophilus--analysis of the gene and characterisation of the protein

The ornithine carbamoyltransferase (OTC) gene from Thermus thermophilus was cloned from a lambda-ZAP genomic library. An ORF of 903 bp was found coding for a protein of Mr 33,200. The coding region has a very high overall G+C content of 68.0%. T. thermophilus OTC displays 38-48% amino acid identity with other OTC, the most closely related proteins being OTC from the archaeon Pyrococcus furiosus and from Bacillus subtilis. The enzyme was expressed in Escherichia coli and purified to homogeneity using a thermoshock followed by affinity chromatography on delta-N-phosphonoacetyl-L-ornithine-Sepharose. The native enzyme has an Mr of about 110,000, suggesting a trimeric structure, as for most anabolic OTC from various organisms. T. thermophilus OTC exhibits Michaelis-Menten kinetics for carbamoyl phosphate and ornithine with a Km(app) of 0.10 mM for both substrates. The pH optimum was dependent on ornithine concentration with an optimum at pH 8 for ornithine concentrations around Km values. Higher concentrations shift the optimum towards lower pH. The optimal temperature was above 65 degrees C and the activation energy 39.1 kJ/mol. The enzyme is highly thermostable. In the presence of its substrates the half-life time was several hours at 85 degrees C. Ionic and hydrophobic interactions contribute to the stability. The expression of T. thermophilus OTC was negatively regulated by arginine.

Genomic organization of a mouse glyceraldehyde 3-phosphate dehydrogenase gene (Gapd-s) expressed in post-meiotic spermatogenic cells

The Gapd-s gene encodes an isoform of the glyceraldehyde 3-phosphate dehydrogenase enzyme expressed only in post-meiotic spermatogenic cells. Two clones containing the Gapd-s gene were isolated from a mouse genomic library. Sequencing and restriction enzyme analysis demonstrated that this single-copy gene contains 11 exons and spans 9596 base pairs. The locations of Gapd-s exons and introns are conserved when compared to the corresponding portions of the chicken and human somatic Gapd genes. The promoter region contains no TATA box, although there is a potential SP1 recognition site within exon 1. Like other TATA-less genes, primer extension analysis reveals some heterogeneity in the site of transcription initiation with Gapd-s transcripts initiating from three discrete sites. Northern analysis demonstrated that a 1.5-kb Gapd-s mRNA is expressed in the testis in at least three mammalian orders, indicating that the Gapd-s gene appeared early in mammalian evolution. Using GAPD-deficient bacteria, mouse GAPD-S was shown to be capable of functioning as a glycolytic enzyme. Since GAPD has been proposed to be a key enzyme regulating glycolysis in spermatogenic cells, GAPD-S may represent a potential target for toxicological or contraceptive agents affecting fertility by interfering with glycolysis.

Analysis of a cis-active sequence mediating catabolite repression in gram-positive bacteria

One form of catabolite repression (CR) in the Gram-positive genus, Bacillus, is mediated by a cis-acting element (CRE). We use here a consensus sequence to identify such elements in sequenced genes of Gram-positive bacteria. These are analysed with respect to position and type of gene in which they occur. CRE sequences near the promoter region are mainly identified in genes encoding carbon catabolic enzymes, which are thus likely to be subject to CR by a global mechanism. Functional aspects of CREs are evaluated.

Relevance of sequence statistics for the properties of extremophilic proteins

The amino acid composition of proteins from mesophilic and extremophilic organisms is commonly assumed to reflect the mechanisms of molecular adaptation to extremes of physical conditions. In this context, halophilic behaviour has been attributed to significantly increased numbers of aspartic and glutamic acid residues. However, extending the analysis to a statistically relevant set of related proteins, dihydrofolate reductase from Halobacterium volcanii, as an example, shows that the increase in negative charge is found to be less significant than other exchanges of amino acids (e.g., Ala, Asn, Arg, Lys, Phe, Ser). Thus, the high water binding capacity of negatively charged residues cannot be unambiguously correlated with the anomalous stability of halophilic proteins. A similar caveat holds for generalizations regarding the thermal stability of proteins. In this case, D-glyceraldehyde-3-phosphate dehydrogenase from the hyperthermophilic bacterium Thermotoga maritima was compared with a number of mesophilic and moderately thermophilic homologs. Again, 'traffic rules of stabilization', in terms of amino acid changes in going from mesophilic to thermophilic proteins, cannot be given.

A glyceraldehyde-3-phosphate dehydrogenase with eubacterial features in the amitochondriate eukaryote, Trichomonas vaginalis

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), localized in the cytosol of Trichomonas vaginalis, was partially purified. The enzyme is specific for NAD+ and is similar in most of its catalytic properties to glycolytic GAPDHs from other organisms. Its sensitivity to koningic acid is similar to levels observed in GAPDHs from eubacteria and two orders of magnitude lower than those observed for eukaryotic GAPDHs. The complete amino acid sequence of T. vaginalis GAPDH was derived from the N-terminal sequence of the purified protein and the deduced sequence of a cDNA clone. It showed great similarity to other eubacterial and eukaryotic GAPDH sequences. The sequence of the S-loop displayed a eubacterial signature. The overall sequence was more similar to eubacterial sequences than to cytosolic and glycosomal eukaryotic sequences. In phylogenetic trees obtained with distance matrix and parsimony methods T. vaginalis GAPDH clustered with its eubacterial homologs. GAPDHs of other amitochondriate protists, belonging to early branches of the eukaryotic lineage (Giardia lamblia and Entamoeba histolytica--Smith M.W. and Doolittle R.F., unpublished data in GenBank), showed typical eukaryotic signatures and clustered with other eukaryotic sequences, indicating that T. vaginalis GAPDH occupies an anomalous position, possibly due to horizontal gene transfer from a eubacterium.

Prediction and analysis of structure, stability and unfolding of thermolysin-like proteases

Bacillus neutral proteases (NPs) form a group of well-characterized homologous enzymes, that exhibit large differences in thermostability. The three-dimensional (3D) structures of several of these enzymes have been modelled on the basis of the crystal structures of the NPs of B. thermoproteolyticus (thermolysin) and B. cereus. Several new techniques have been developed to improve the model-building procedures. Also a 'model-building by mutagenesis' strategy was used, in which mutants were designed just to shed light on parts of the structures that were particularly hard to model. The NP models have been used for the prediction of site-directed mutations aimed at improving the thermostability of the enzymes. Predictions were made using several novel computational techniques, such as position-specific rotamer searching, packing quality analysis and property-profile database searches. Many stabilizing mutations were predicted and produced: improvement of hydrogen bonding, exclusion of buried water molecules, capping helices, improvement of hydrophobic interactions and entropic stabilization have been applied successfully. At elevated temperatures NPs are irreversibly inactivated as a result of autolysis. It has been shown that this denaturation process is independent of the protease activity and concentration and that the inactivation follows first-order kinetics. From this it has been conjectured that local unfolding of (surface) loops, which renders the protein susceptible to autolysis, is the rate-limiting step. Despite the particular nature of the thermal denaturation process, normal rules for protein stability can be applied to NPs. However, rather than stabilizing the whole protein against global unfolding, only a small region has to be protected against local unfolding. In contrast to proteins in general, mutational effects in proteases are not additive and their magnitude is strongly dependent on the location of the mutation. Mutations that alter the stability of the NP by a large amount are located in a relatively weak region (or more precisely, they affect a local unfolding pathway with a relatively low free energy of activation). One weak region, that is supposedly important in the early steps of NP unfolding, has been determined in the NP of B. stearothermophilus. After eliminating this weakest link a drastic increase in thermostability was observed and the search for the second-weakest link, or the second-lowest energy local unfolding pathway is now in progress. Hopefully, this approach can be used to unravel the entire early phase of unfolding.

Nucleotide sequence and molecular evolution of the gene coding for glyceraldehyde-3-phosphate dehydrogenase in the thermoacidophilic archaebacterium Sulfolobus solfataricus

A Sulfolobus solfataricus genomic library cloned in the EMBL3 phage was screened using as probes synthetic oligonucleotides designed from the known amino acid sequence of a peptide obtained from the purified glyceraldehyde-3-phosphate dehydrogenase (aGAPD) protein. The screening led to the isolation of six recombinant phages (lambda G1-lambda G6) and one of them (lambda G4) contained the entire GAPD gene. The deduced amino acid sequence accounts for a protein made of 341 amino acids and the initial methionine is encoded by a GTG triplet. Alignment of the S. solfataricus aGAPD sequence versus GAPD from archaea, eukarya, and bacteria showed that aGAPD is very similar to other archaebacterial but not to eukaryotic or eubacterial GAPD. For known archaebacterial GAPD sequences, the rate of nucleotide substitutions per site per year showed that these sequences are homologous not only at the amino acid but also at the nucleotide level. The evolutionary rates are nearly similar to those reported for other eukaryotic genes.

Structural consequences of sequence patterns in the fingerprint region of the nucleotide binding fold. Implications for nucleotide specificity

The dinucleotide binding beta alpha beta motif in the crystal structures of seven different enzymes has been analysed in terms of their three-dimensional structures and primary sequences. We have identified that the hydrogen bonding of the adenine ribose to the glycine-rich turn containing the fingerprint sequence GXGXXG/A occurs via a direct or indirect mechanism, depending on the nature of the fingerprint sequence but independent of coenzyme specificity. The major determinant of the type of interaction is the nature of the residue occupying the last position of the above fingerprint. In the NAD(+)-linked dehydrogenases, an acidic residue is commonly used to form important hydrogen bonds to the adenine ribose hydroxyls and, hitherto, this residue has been thought to be an indicator of NAD+ specificity. However, on the basis of the three-dimensional structure of the NAD(+)-linked glutamate dehydrogenase (GDH) from Clostridium symbiosum we have demonstrated that this residue is not a universal requirement for the construction of an NAD+ binding site. Furthermore, considerations of sequence homology unambiguously identify an equivalent acidic residue in both NADP+ and dual specificity glutamate dehydrogenases. The conservation of this residue in these enzymes, coupled to its close proximity to the 2' phosphate implied by the necessary similarity in three-dimensional structure to C. symbiosum GDH, implicates this residue in the recognition of the 2' phosphate either via water-mediated or direct hydrogen-bonding schemes. Analysis of the latter has led us to suggest that two patterns of recognition for the 2' phosphate group of NADP(+)-binding enzymes may exist, which are distinguished by the ionization state of the 2' phosphate.

Thermal analysis of bacteria by differential scanning calorimetry: relationship of protein denaturation in situ to maximum growth temperature

Differential scanning calorimetry (DSC) was used to analyze thermal transitions in two strains of the thermophile Bacillus stearothermophilus (ATCC 12016 and WAT), the mesophile Bacillus megaterium and the psychrotroph Bacillus psychrophilus. The observed transitions, representing lipid melting and DNA and protein unfolding, are compared to the maximum growth temperature (Tmax) in each species as a means of identifying critical, thermolabile targets responsible for heat-induced inhibition of growth. A low temperature, lipid transition was detected in B. stearothermophilus and B. megaterium which varied slightly with Tmax but whose high temperature end is always 22-33 degrees C below Tmax. The transition temperature (Tm) of the main melting of DNA varies from 88 to 92 degrees C, 23-32 degrees C above Tmax. The main part of the profile representing irreversible transitions is resolvable into at least three distinct peaks and is identified primarily with protein denaturation. The onset temperature for denaturation (Tl), i.e., minimum temperature of detectable denaturation, is somewhat dependent on growth temperature (Tg). Tmax for B. stearothermophilus ATCC and WAT is 69 and 56 degrees C, respectively. For cells grown between 4 and 20 degrees C below Tmax, Tl is 2-4 degrees C lower than Tmax, demonstrating that some denaturation can be tolerated before complete inhibition of growth and suggesting that inhibition of growth is due to the denaturation of a critical protein with a Tm a few degrees above Tl or to the accumulation of denatured protein to a critical level. A similar pattern holds for B. megaterium and B. psychrophilus, except that Tmax is 48 and 32.5 degrees C (Tl = 45-46 degrees C and 30 degrees C), respectively. Thus, there is an excellent correlation between the onset of protein denaturation and maximum growth temperature for these three species of the same genus. This study also demonstrates the applicability of DSC for resolving transitions in intact cells on the basis of thermostability of cellular constituents and for obtaining an overall view of macromolecular stability.