The solution structure of VAT-N reveals a 'missing link' in the evolution of complex enzymes from a simple betaalphabetabeta element.

BACKGROUND: The VAT protein of the archaebacterium Thermoplasma acidophilum, like all other members of the Cdc48/p97 family of AAA ATPases, has two ATPase domains and a 185-residue amino-terminal substrate-recognition domain, VAT-N. VAT shows activity in protein folding and unfolding and thus shares the common function of these ATPases in disassembly and/or degradation of protein complexes. RESULTS: Using nuclear magnetic resonance (NMR) spectroscopy, we found that VAT-N is composed of two equally sized subdomains. The amino-terminal subdomain VAT-Nn (comprising residues Met1-Thr92) forms a double-psi beta-barrel whose pseudo-twofold symmetry is mirrored by an internal sequence repeat of 42 residues. The carboxy-terminal subdomain VAT-Nc (comprising residues Glu93-Gly185) forms a novel six-stranded beta-clam fold. Together, VAT-Nn and VAT-Nc form a kidney-shaped structure, in close agreement with results from electron microscopy. Sequence and structure analyses showed that VAT-Nn is related to numerous proteins including prokaryotic transcription factors, metabolic enzymes, the protease cofactors UFD1 and PrlF, and aspartic proteinases. These proteins map out an evolutionary path from simple homodimeric transcription factors containing a single copy of the VAT-Nn repeat to complex enzymes containing four copies. CONCLUSIONS: Our results suggest that VAT-N is a precursor of the aspartic proteinases that has acquired peptide-binding activity while remaining proteolytically incompetent. We propose that the binding site of the protein is similar to that of aspartic proteinases, in that it lies between the psi-loops of the amino-terminal beta-barrel and that it coincides with a crescent-shaped band of positive charge extending across the upper face of the molecule.

Self-consistently optimized statistical mechanical energy functions for sequence structure alignment.

A quantitative form of the principle of minimal frustration is used to obtain from a database analysis statistical mechanical energy functions and gap parameters for aligning sequences to three-dimensional structures. The analysis that partially takes into account correlations in the energy landscape improves upon the previous approximations of Goldstein et al. (1994, 1995) (Goldstein R, Luthey-Schulten Z, Wolynes P, 1994, Proceedings of the 27th Hawaii International Conference on System Sciences. Los Alamitos, California: IEEE Computer Society Press. pp 306-315; Goldstein R, Luthey-Schulten Z, Wolynes P, 1995, In: Elber R, ed. New developments in theoretical studies of proteins. Singapore: World Scientific). The energy function allows for ordering of alignments based on the compatibility of a sequence to be in a given structure (i.e., lowest energy) and therefore removes the necessity of using percent identity or similarity as scoring parameters. The alignments produced by the energy function on distant homologues with low percent identity (less than 21%) are generally better than those generated with evolutionary information. The lowest energy alignment generated with the energy function for sequences containing prosite signatures but unknown structures is a structure containing the same prosite signature, providing a check on the robustness of the algorithm. Finally, the energy function can make use of known experimental evidence as constraints within the alignment algorithm to aid in finding the correct structural alignment.

Fold recognition from sequence comparisons.

We applied a new protocol based on PSI-Blast to predict the structures of fold recognition targets during CASP4. The protocol used a back-validation step to infer biologically significant connections between sequences with PSI-Blast E-values up to 10. If connections were found to proteins of known structure, alignments were generated by using HMMer. The protocol was implemented in a fully automated version (SBauto) and in a version that allowed manual intervention (SBfold). We found that the automated version made 17 predictions for target domains, of which 8 identified the correct fold with an average alignment accuracy of 24% for alignable residues and 43% for equivalent secondary structure elements. The manual version improved predictions somewhat, with 10 of 15 predictions identifying the correct fold with alignment accuracies of 33% for alignable residues and 64% for equivalent secondary structure elements. We describe successes and failures of our approach and discuss future developments of fold recognition.

Self-consistently optimized energy functions for protein structure prediction by molecular dynamics.

The protein energy landscape theory is used to obtain optimal energy functions for protein structure prediction via simulated annealing. The analysis here takes advantage of a more complete statistical characterization of the protein energy landscape and thereby improves on previous approximations. This schema partially takes into account correlations in the energy landscape. It also incorporates the relationships between folding dynamics and characteristic energy scales that control the collapse of the proteins and modulate rigidity of short-range interactions. Simulated annealing for the optimal energy functions, which are associative memory hamiltonians using a database of folding patterns, generally leads to quantitatively correct structures. In some cases the algorithm achieves "creativity," i.e., structures result that are better than any homolog in the database.

Structural basis for thermostability and identification of potential active site residues for adenylate kinases from the archaeal genus Methanococcus.

Sequence comparisons of highly related archaeal adenylate kinases (AKs) from the mesophilic Methanococcus voltae, the moderate thermophile Methanococcus thermolithotrophicus, and two extreme thermophiles Methanococcus igneus and Methanococcus jannaschii, allow identification of interactions responsible for the large variation in temperatures for optimal catalytic activity and thermostabilities observed for these proteins. The tertiary structures of the methanococcal AKs have been predicted by using homology modeling to further investigate the potential role of specific interactions on thermal stability and activity. The alignments for the methanococcal AKs have been generated by using an energy-based sequence-structure threading procedure against high-resolution crystal structures of eukaryotic, eubacterial, and mitochondrial adenylate and uridylate (UK) kinases. From these alignments, full atomic model structures have been produced using the program MODELLER. The final structures allow identification of potential active site interactions and place a polyproline region near the active site, both of which are unique to the archaeal AKs. Based on these model structures, the additional polar residues present in the thermophiles could contribute four additional salt bridges and a higher negative surface charge. Since only one of these possible salt bridges is interior, they do not appear significantly to the thermal stability. Instead, our model structures indicate that a larger and more hydrophobic core, due to a specific increase in aliphatic amino acid content and aliphatic side chain volume, in the thermophilic AKs is responsible for increased thermal stability.

The Janus face of the archaeal Cdc48/p97 homologue VAT: protein folding versus unfolding.

Members of the AAA family of ATPases have been implicated in chaperone-like activities. We used the archaeal Cdc48/p97 homologue VAT as a model system to investigate the effect of an AAA protein on the folding and unfolding of two well-studied, heterologous substrates, cyclophilin and penicillinase. We found that, depending on the Mg2+ concentration, VAT assumes two states with maximum rates of ATP hydrolysis that differ by an order of magnitude. In the low-activity state, VAT accelerated the refolding of penicillinase, whereas in the high-activity state, it accelerated its unfolding. Both reactions were ATP-dependent. In its interaction with cyclophilin, VAT was ATP-independent and only promoted refolding. The N-terminal domain of VAT, which lacks ATPase activity, also accelerated the refolding of cyclophilin but showed no effect on penicillinase. VAT appears to be structurally equivalent over its entire length to Sec18/NSF, suggesting that these results apply more broadly to group II AAA proteins.

Fold recognition using sequence and secondary structure information.

We applied a succession of sequence search and structure prediction methods to the targets in the fold recognition part of the CASP3 experiment. For each target, we expanded an initial sequence space, obtained through PSI-BLAST, by searching for statistically significant relationships to low-scoring sequences and then by searching for conserved sequence patterns. We then divided the proteins in the sequence space into families and built an alignment hierarchically, using the multiple alignment program MACAW. If no significant similarity to a protein of known structure was apparent at this point, we submitted the alignment to the Jpred server for consensus secondary structure prediction and searched the structure space using the secondary structure mapping program MAP. Failing this, we compared the structural properties that we believed we recognized in the aligned proteins to the folds in the SCOP database, using visual inspection. If all these methods failed to uncover a plausible match, we predicted that the target would adopt a novel fold. This procedure yielded correct answers for seven of twenty-one targets and a partly correct answer for one. A retrospective analysis shows that automating the sequence search procedures would have represented a significant improvement, with at least three additional correct predictions.

Evolution of two-component signal transduction.

Two-component signal transduction (TCST) systems are the principal means for coordinating responses to environmental changes in bacteria as well as some plants, fungi, protozoa, and archaea. These systems typically consist of a receptor histidine kinase, which reacts to an extracellular signal by phosphorylating a cytoplasmic response regulator, causing a change in cellular behavior. Although several model systems, including sporulation and chemotaxis, have been extensively studied, the evolutionary relationships between specific TCST systems are not well understood, and the ancestry of the signal transduction components is unclear. Phylogenetic trees of TCST components from 14 complete and 6 partial genomes, containing 183 histidine kinases and 220 response regulators, were constructed using distance methods. The trees showed extensive congruence in the positions of 11 recognizable phylogenetic clusters. Eukaryotic sequences were found almost exclusively in one cluster, which also showed the greatest extent of domain variability in its component proteins, and archaeal sequences mainly formed species-specific clusters. Three clusters in different parts of the kinase tree contained proteins with serine-phosphorylating activity. All kinases were found to be monophyletic with respect to other members of their superfamily, such as type II topoisomerases and Hsp90. Structural analysis further revealed significant similarity to the ATP-binding domain of eukaryotic protein kinases. TCST systems are of bacterial origin and radiated into archaea and eukaryotes by lateral gene transfer. Their components show extensive coevolution, suggesting that recombination has not been a major factor in their differentiation. Although histidine kinase activity is prevalent, serine kinases have evolved multiple times independently within this family, accompanied by a loss of the cognate response regulator(s). The structural and functional similarity between TCST kinases and eukaryotic protein kinases raises the possibility of a distant evolutionary relationship.

Phylogenetic analyses do not support horizontal gene transfers from bacteria to vertebrates.

Horizontal gene transfer (HGT) has long been recognized as a principal force in the evolution of genomes. Genome sequences of Archaea and Bacteria have revealed the existence of genes whose similarity to loci in distantly related organisms is explained most parsimoniously by HGT events. In most multicellular organisms, such genetic fixation can occur only in the germ line. Therefore, it is notable that the publication of the human genome reports 113 incidents of direct HGT between bacteria and vertebrates, without any apparent occurrence in evolutionary intermediates, that is, non-vertebrate eukaryotes. Phylogenetic analysis arguably provides the most objective approach for determining the occurrence and directionality of HGT. Here we report a phylogenetic analysis of 28 proposed HGT genes, whose presence in the human genome had been confirmed by polymerase chain reaction (PCR). The results indicate that most putative HGT genes are present in more anciently derived eukaryotes (many such sequences available in non-vertebrate EST databases) and can be explained in terms of descent through common ancestry. They are, therefore, unlikely to be examples of direct HGT from bacteria to vertebrates.

The genome sequence of the thermoacidophilic scavenger Thermoplasma acidophilum.

Thermoplasma acidophilum is a thermoacidophilic archaeon that thrives at 59 degrees C and pH 2, which was isolated from self-heating coal refuse piles and solfatara fields. Species of the genus Thermoplasma do not possess a rigid cell wall, but are only delimited by a plasma membrane. Many macromolecular assemblies from Thermoplasma, primarily proteases and chaperones, have been pivotal in elucidating the structure and function of their more complex eukaryotic homologues. Our interest in protein folding and degradation led us to seek a more complete representation of the proteins involved in these pathways by determining the genome sequence of the organism. Here we have sequenced the 1,564,905-base-pair genome in just 7,855 sequencing reactions by using a new strategy. The 1,509 open reading frames identify Thermoplasma as a typical euryarchaeon with a substantial complement of bacteria-related genes; however, evidence indicates that there has been much lateral gene transfer between Thermoplasma and Sulfolobus solfataricus, a phylogenetically distant crenarchaeon inhabiting the same environment. At least 252 open reading frames, including a complete protein degradation pathway and various transport proteins, resemble Sulfolobus proteins most closely.

Systematic identification of selective essential genes in Helicobacter pylori by genome prioritization and allelic replacement mutagenesis.

A comparative genomic approach was used to identify Helicobacter pylori 26695 open reading frames (ORFs) which are conserved in H. pylori J99 but highly diverged in other eubacteria. A survey of selected pathways of central intermediary metabolism was also carried out, and genes with a potentially selective role in H. pylori were identified. Forty-five ORFs identified in these two analyses were screened using a rapid vector-free allelic replacement mutagenesis technique, and 33 were shown to be essential in vitro. Notably, 13 ORFs gave essentiality results which are unexpected in view of their known or proposed functions, and phylogenetic analysis was used to investigate the annotation of 7 such ORFs which are highly diverged. We propose that the products of a number of these H. pylori-specific essential genes may be suitable targets for novel anti-H. pylori therapies.

A genomic analysis of two-component signal transduction in Streptococcus pneumoniae.

A genomics-based approach was used to identify the entire gene complement of putative two-component signal transduction systems (TCSTSs) in Streptococcus pneumoniae. A total of 14 open reading frames (ORFs) were identified as putative response regulators, 13 of which were adjacent to genes encoding probable histidine kinases. Both the histidine kinase and response regulator proteins were categorized into subfamilies on the basis of phylogeny. Through a systematic programme of mutagenesis, the importance of each novel TCSTS was determined with respect to viability and pathogenicity. One TCSTS was identified that was essential for the growth of S. pneumoniaeThis locus was highly homologous to the yycFG gene pair encoding the essential response regulator/histidine kinase proteins identified in Bacillus subtilis and Staphylococcus aureus. Separate deletions of eight other loci led in each case to a dramatic attenuation of growth in a mouse respiratory tract infection model, suggesting that these signal transduction systems are important for the in vivo adaptation and pathogenesis of S. pneumoniae. The identification of conserved TCSTSs important for both pathogenicity and viability in a Gram-positive pathogen highlights the potential of two-component signal transduction as a multicomponent target for antibacterial drug discovery.