Contribution of electrostatic interactions, compactness and quaternary structure to protein thermostability: lessons from structural genomics of Thermotoga maritima.

Studies of the structural basis of protein thermostability have produced a confusing picture. Small sets of proteins have been analyzed from a variety of thermophilic species, suggesting different structural features as responsible for protein thermostability. Taking advantage of the recent advances in structural genomics, we have compiled a relatively large protein structure dataset, which was constructed very carefully and selectively; that is, the dataset contains only experimentally determined structures of proteins from one specific organism, the hyperthermophilic bacterium Thermotoga maritima, and those of close homologs from mesophilic bacteria. In contrast to the conclusions of previous studies, our analyses show that oligomerization order, hydrogen bonds, and secondary structure play minor roles in adaptation to hyperthermophily in bacteria. On the other hand, the data exhibit very significant increases in the density of salt-bridges and in compactness for proteins from T.maritima. The latter effect can be measured by contact order or solvent accessibility, and network analysis shows a specific increase in highly connected residues in this thermophile. These features account for changes in 96% of the protein pairs studied. Our results provide a clear picture of protein thermostability in one species, and a framework for future studies of thermal adaptation.

Accuracy of replication in the polymerase chain reaction. Comparison between Thermotoga maritima DNA polymerase and Thermus aquaticus DNA polymerase

For certain applications of the polymerase chain reaction (PCR), it may be necessary to consider the accuracy of replication. The breakthrough that made PCR user friendly was the commercialization of Thermus aquaticus (Taq) DNA polymerase, an enzyme that would survive the high temperatures needed for DNA denaturation. The development of enzymes with an inherent 3' to 5' exonuclease proofreading activity, lacking in Taq polymerase, would be an improvement when higher fidelity is needed. We used the forward mutation assay to compare the fidelity of Taq polymerase and Thermotoga maritima (ULTMA™) DNA polymerase, an enzyme that does have proofreading activity. We did not find significant differences in the fidelity of either enzyme, even when using optimal buffer conditions, thermal cycling parameters, and number of cycles (0.2% and 0.13% error rates for ULTMA™ and Taq, respectively, after reading about 3,000 bases each). We conclude that for sequencing purposes there is no difference in using a DNA polymerase that contains an inherent 3' to 5' exonuclease activity for DNA amplification. Perhaps the specificity and fidelity of PCR are complex issues influenced by the nature of the target sequence, as well as by each PCR component.

Purification and Preliminary Characterization of the TM0727 Protein from Thermotoga maritima

The TM0727 gene of Thermotoga maritima is responsible for encoding what has been reported to be a modulator of DNA gyrase (pmbA). Although the function of pmbA is still unknown, it is believedto be involved in cell division, carbon storage regulation, and the synthesis of the antibiotic peptide microcin B17. It is suggested that it serves together with tldD, a known zinc dependent protease, tomodulate DNA gyrase. TM0727 is believed to be a zinc dependent protease that binds zinc in the central active site of the molecule, located between two equivalent monomeric units. However, thecrystal structure determined by Wilson et al. (2005) did not contain zinc. It therefore remains to be seen if TM0727 requires zinc for activity, or regulation, and if the protein is indeed a protease. To begin studying this protein, the gene was expressed in BL21(DE3) pLysS cells and the induction time was optimized. Using affinity and ion exchange chromatography, the protein has been successfully purified. The purification procedure can be replicated to obtain sufficient protein for characterization. Purification results show that the protein loses stability after 24 hours and remains stable under an imidazole-free lysis workup. Preliminary characterization of TM0727 has focused on understanding the protein’s structuralproperties through tryptophan fluorescence anisotropy measurements. The four tryptophan residues located within the TM0727 dimer fluoresce at different maximum wavelengths and with differentintensities upon excitation with 295nm light. These emission properties are highly sensitive to the environment (solvent, surrounding residues) of each tryptophan residue. The low number oftryptophans allows for a specific monitoring of the protein’s structure as it denatures. As more denaturant is added to the protein, its tryptophan environments have clearly altered. This is indicative of unfolding and increased solvent exposure of the protein. This unfolding has been confirmed with the addition of a fluorescent quencher. Additionally, fluorescence anisotropy measurements have been carried out on the protein to gain a preliminary understanding of the rotational dynamics of the tryptophan residues. These experiments excite the tryptophan residues within the sample using a polarized light source. Polarized emission is then detected, the degree of which depends on the rotational dynamics and local environment of the tryptophan residues. The protein was denatured and the changes in emission were recorded to detect these structural changes. Results have shown a large change in quaternary structure, consistent with a dimer to monomer transition, occurs at 1.5M Guandidine HCl. There has also been an examination of the crystal structure for the location of a potential active site. The inner cavity of the protein was inspected visually to locate a potential location for a catalytic triad, specifically the amino acids found in the active sites of serine, cyteine, and aspartateproteases. It was found that a potential aspartic protease active site may be located between the Asparate286 and Aspartate287 residues. Further investigation is warranted to test this remotepossibility.

Both DNA gyrase and reverse gyrase are present in the hyperthermophilic bacterium Thermotoga maritima

Like all hyperthermophiles yet tested, the bacterium Thermotoga maritima contains a reverse gyrase. Here we show that it contains also a DNA gyrase. The genes top2A and top2B encoding the two subunits of a DNA gyrase-like enzyme have been cloned and sequenced. The Top2A (type II DNA topoisomerase A protein) is more similar to GyrA (DNA gyrase A protein) than to ParC [topoisomerase IV (Topo IV) C protein]. The difference is especially striking at the C-terminal domain, which differentiates DNA gyrases from Topo IV. DNA gyrase activity was detected in T. maritima and purified to homogeneity using a novobiocin-Sepharose column. This hyperhermophilic DNA gyrase has an optimal activity around 82–86°C. In contrast to plasmids from hyperthermophilic archaea, which are from relaxed to positively supercoiled, we found that the plasmid pRQ7 from Thermotoga sp. RQ7 is negatively supercoiled. pRQ7 became positively supercoiled after addition of novobiocin to cell cultures, indicating that its negative supercoiling is due to the DNA gyrase of the host strain. The findings concerning DNA gyrase and negative supercoiling in Thermotogales put into question the role of reverse gyrase in hyperthermophiles.

Characterization of RNase P from Thermotoga maritima

The protein subunit of RNase P from a thermophilic bacterium, Thermotoga maritima, was overexpressed in and purified from Escherichia coli. The cloned protein was reconstituted with the RNA subunit transcribed in vitro. The temperature optimum of the holoenzyme is near 50°C, with no enzymatic activity at 65°C or above. This finding is in sharp contrast to the optimal growth temperature of T.maritima, which is near 80°C. However, in heterologous reconstitution experiments in vitro with RNase P subunits from other species, we found that the protein subunit from T.maritima was responsible for the comparative thermal stability of such complexes.

Thermostable variants of the recombinant xylanase A from Bacillus subtilis produced by directed evolution show reduced heat capacity changes

Directed evolution techniques have been used to improve the thermal stability of the xylanase A from Bacillus subtilis (XylA). Two generations of random mutant libraries generated by error prone PCR coupled with a single generation of DNA shuffling produced a series of mutant proteins with increasing thermostability. The most Thermostable XylA variant from the third generation contained four mutations Q7H, G13R, S22P, and S179C that showed an increase in melting temperature of 20 degrees C. The thermodynamic properties Of a representative subset of nine XylA variants showing a range of thermostabilities were measured by thermal denaturation as monitored by the change in the far ultraviolet circular dichroism signal. Analysis of the data from these thermostable variants demonstrated a correlation between the decrease in the heat capacity change (Delta C(p)) with an increase in the midpoint of the transition temperature (T(m)) on transition from the native to the unfolded state. This result could not be interpreted within the context of the changes in accessible surface area of the protein on transition from the native to unfolded states. Since all the mutations are located at the surface of the protein, these results suggest that an explanation of the decrease in Delta C(p) on should include effects arising from the prot inlsolvent interface.

On surrogate methods for detecting lateral gene transfer

Surrogate methods for detecting lateral gene transfer are those that do not require inference of phylogenetic trees. Herein I apply four such methods to identify open reading frames (ORFs) in the genome of Escherichia coli K12 that may have arisen by lateral gene transfer. Only two of these methods detect the same ORFs more frequently than expected by chance, whereas several intersections contain many fewer ORFs than expected. Each of the four methods detects a different non-random set of ORFs. The methods may detect lateral ORFs of different relative ages; testing this hypothesis will require rigorous inference of trees. (C) 2001 Federation of European Microbiological Societies. Published by Elsevier Science BN. All rights reserved.

Phylogenetic analysis of evolutionary relationships of the planctomycete division of the domain bacteria based on amino acid sequences of elongation factor Tu

Sequences from the tuf gene coding for the elongation factor EF-Tu were amplified and sequenced from the genomic DNA of Pirellula marina and Isosphaera pallida, two species of bacteria within the order Planctomycetales. A near-complete (1140-bp) sequence was obtained from Pi. marina and a partial (759-bp) sequence was obtained for I. pallida. Alignment of the deduced Pi. marina EF-Tu amino acid sequence against reference sequences demonstrated the presence of a unique Il-amino acid sequence motif not present in any other division of the domain Bacteria. Pi. marina shared the highest percentage amino acid sequence identity with I. pallida but showed only a low percentage identity with other members of the domain Bacteria. This is consistent with the concept of the planctomycetes as a unique division of the Bacteria. Neither primary sequence comparison of EF-Tu nor phylogenetic analysis supports any close relationship between planctomycetes and the chlamydiae, which has previously been postulated on the basis of 16S rRNA. Phylogenetic analysis of aligned EF-Tu amino acid sequences performed using distance, maximum-parsimony, and maximum likelihood approaches yielded contradictory results with respect to the position of planctomycetes relative to other bacteria, It is hypothesized that long-branch attraction effects due to unequal evolutionary rates and mutational saturation effects may account for some of the contradictions.

Distributional profiles of homologous open reading frames among bacterial phyla: implications for vertical and lateral transmission

If open reading frames (ORFs) have been transmitted primarily by vertical descent, the distributional profile of orthologues of each ORF should be congruent with the organismal tree or a subtree thereof. Distributional patterns not reconciled parsimoniously with tree-like descent and loss are prima facie evidence of lateral gene transfer. Herein, a rigorous criterion for recognizing ORF distributions is described and implemented; it does not require the inference of phylogenetic trees, nor does it assume any specific tree. Because lineage-specific differences in rates of sequence change can also generate unexpected distributional patterns, rate artefacts, were controlled for by requiring pairwise matches between ORFs to exceed a rigorous inclusion threshold, but absence of a match was assessed against a more-permissive exclusion threshold. Applying this dual-threshold criterion to cross-domain and cross-phylum distributional patterns for ORFs in 23 bacterial genomes, a relative abundance of ORFs was observed that find a match in exactly seven other bacterial phyla; 94-99% of these ORFs also find matches among the Archaea and/or Eukarya. In the larger (and some smaller) bacterial genomes, ORFs that find matches in exactly one other bacterial phylum are also relatively abundant, but fewer of these have non-bacterial homologues; most of their matches within the Bacteria are to the Proteobacteria and/or Firmicutes, which cannot be sister lineages to all bacteria. ORFs that are neither distributed universally among the Bacteria, nor necessarily shared with topologically adjacent lineages, are preferentially enriched in large bacterial genomes.

Dissecting structure-function-stability relationships of a thermostable GH5-CBM3 cellulase from Bacillus subtilis 168

Cellulases participate in a number of biological events, such as plant cell wall remodelling, nematode parasitism and microbial carbon uptake. Their ability to depolymerize crystalline cellulose is of great biotechnological interest for environmentally compatible production of fuels from lignocellulosic biomass. However, industrial use of cellulases is somewhat limited by both their low catalytic efficiency and stability. In the present study, we conducted a detailed functional and structural characterization of the thermostable BsCe15A (Bacillus subtilis cellulase 5A), which consists of a GH5 (glycoside hydrolase 5) catalytic domain fused to a CBM3 (family 3 carbohydrate-binding module). NMR structural analysis revealed that the Bacillus CBM3 represents a new subfamily, which lacks the classical calcium-binding motif, and variations in NMR frequencies in the presence of cellopentaose showed the importance of polar residues in the carbohydrate interaction. Together with the catalytic domain, the CBM3 forms a large planar surface for cellulose recognition, which conducts the substrate in a proper conformation to the active site and increases enzymatic efficiency. Notably, the manganese ion was demonstrated to have a hyper-stabilizing effect on BsCel5A, and by using deletion constructs and X-ray crystallography we determined that this effect maps to a negatively charged motif located at the opposite face of the catalytic site.

B12-dependent ribonucleotide reductases from deeply rooted eubacteria are structurally related to the aerobic enzyme from Escherichia coli

The ribonucleotide reductases from three ancient eubacteria, the hyperthermophilic Thermotoga maritima (TM), the radioresistant Deinococcus radiodurans (DR), and the thermophilic photosynthetic Chloroflexus aurantiacus, were found to be coenzyme-B12 (class II) enzymes, similar to the earlier described reductases from the archaebacteria Thermoplasma acidophila and Pyrococcus furiosus. Reduction of CDP by the purified TM and DR enzymes requires adenosylcobalamin and DTT. dATP is a positive allosteric effector, but stimulation of the TM enzyme only occurs close to the temperature optimum of 80–90°C. The TM and DR genes were cloned by PCR from peptide sequence information. The TM gene was sequenced completely and expressed in Escherichia coli. The deduced amino acid sequences of the two eubacterial enzymes are homologous to those of the archaebacteria. They can also be aligned to the sequence of the large protein of the aerobic E. coli ribonucleotide reductase that belongs to a different class (class I), which is not dependent on B12. Structure determinations of the E. coli reductase complexed with substrate and allosteric effectors earlier demonstrated a 10-stranded β/α-barrel in the active site. From the conservation of substrate- and effector-binding residues we propose that the B12-dependent class II enzymes contain a similar barrel.

Root of the universal tree of life based on ancient aminoacyl-tRNA synthetase gene duplications.

Universal trees based on sequences of single gene homologs cannot be rooted. Iwabe et al. [Iwabe, N., Kuma, K.-I., Hasegawa, M., Osawa, S. & Miyata, T. (1989) Proc. Natl. Acad. Sci. USA 86, 9355-9359] circumvented this problem by using ancient gene duplications that predated the last common ancestor of all living things. Their separate, reciprocally rooted gene trees for elongation factors and ATPase subunits showed Bacteria (eubacteria) as branching first from the universal tree with Archaea (archaebacteria) and Eucarya (eukaryotes) as sister groups. Given its topical importance to evolutionary biology and concerns about the appropriateness of the ATPase data set, an evaluation of the universal tree root using other ancient gene duplications is essential. In this study, we derive a rooting for the universal tree using aminoacyl-tRNA synthetase genes, an extensive multigene family whose divergence likely preceded that of prokaryotes and eukaryotes. An approximately 1600-bp conserved region was sequenced from the isoleucyl-tRNA synthetases of several species representing deep evolutionary branches of eukaryotes (Nosema locustae), Bacteria (Aquifex pyrophilus and Thermotoga maritima) and Archaea (Pyrococcus furiosus and Sulfolobus acidocaldarius). In addition, a new valyl-tRNA synthetase was characterized from the protist Trichomonas vaginalis. Different phylogenetic methods were used to generate trees of isoleucyl-tRNA synthetases rooted by valyl- and leucyl-tRNA synthetases. All isoleucyl-tRNA synthetase trees showed Archaea and Eucarya as sister groups, providing strong confirmation for the universal tree rooting reported by Iwabe et al. As well, there was strong support for the monophyly (sensu Hennig) of Archaea. The valyl-tRNA synthetase gene from Tr. vaginalis clustered with other eukaryotic ValRS genes, which may have been transferred from the mitochondrial genome to the nuclear genome, suggesting that this amitochondrial trichomonad once harbored an endosymbiotic bacterium.

Highways of gene sharing in prokaryotes

The extent to which lateral genetic transfer has shaped microbial genomes has major implications for the emergence of community structures. We have performed a rigorous phylogenetic analysis of > 220,000 proteins from genomes of 144 prokaryotes to determine the contribution of gene sharing to current prokaryotic diversity, and to identify highways of sharing between lineages. The inferred relationships suggest a pattern of inheritance that is largely vertical, but with notable exceptions among closely related taxa, and among distantly related organisms that live in similar environments.

Focusing in on structural genomics: The University of Queensland structural biology pipeline

The flood of new genomic sequence information together with technological innovations in protein structure determination have led to worldwide structural genomics (SG) initiatives. The goals of SG initiatives are to accelerate the process of protein structure determination, to fill in protein fold space and to provide information about the function of uncharacterized proteins. In the long-term, these outcomes are likely to impact on medical biotechnology and drug discovery, leading to a better understanding of disease as well as the development of new therapeutics. Here we describe the high throughput pipeline established at the University of Queensland in Australia. In this focused pipeline, the targets for structure determination are proteins that are expressed in mouse macrophage cells and that are inferred to have a role in innate immunity. The aim is to characterize the molecular structure and the biochemical and cellular function of these targets by using a parallel processing pipeline. The pipeline is designed to work with tens to hundreds of target gene products and comprises target selection, cloning, expression, purification, crystallization and structure determination. The structures from this pipeline will provide insights into the function of previously uncharacterized macrophage proteins and could lead to the validation of new drug targets for chronic obstructive pulmonary disease and arthritis. (c) 2006 Elsevier B.V. All rights reserved.

The carbohydrate-binding specificity and molecular modelling of Canavalia maritima and Dioclea grandiflora lectins

The carbohydrate-binding specificity of lectins from the seeds of Canavalia maritima and Dioclea grandiflora was studied by hapten-inhibition of haemagglutination using various sugars and sugar derivatives as inhibitors, including N-acetylneuraminic acid and N-acetylmuramic acid. Despite some discrepancies, both lectins exhibited a very similar carbohydrate-binding specificity as previously reported for other lectins from Diocleinae (tribe Phaseoleae, sub-tribe Diocleinae). Accordingly, both lectins exhibited almost identical hydropathic profiles and their three-dimensional models built up from the atomic coordinates of ConA looked very similar. However, docking experiments of glucose and mannose in their monosaccharide-binding sites, by comparison with the ConA-mannose complex used as a model, revealed conformational changes in side chains of the amino acid residues involved in the binding of monosaccharides. These results fully agree with crystallographic data showing that binding of specific ligands to ConA requires conformational chances of its monosaccharide-binding site.