Potential for dramatic improvement in sequence alignment against structures of remote homologous proteins by extracting structural information from multiple structure alignment

An Improved Strategy for Task Scheduling in the Parallel Computational Alignment of Multiple Sequences

Task scheduling in parallel multiple sequence alignment (MSA) through improved dynamic programming optimization speeds up alignment processing. The increased importance of multiple matching sequences also needs the utilization of parallel processor systems. This dynamic algorithm proposes improved task scheduling in case of parallel MSA. Specifically, the alignment of several tertiary structured proteins is computationally complex than simple word-based MSA. Parallel task processing is computationally more efficient for protein-structured based superposition. The basic condition for the application of dynamic programming is also fulfilled, because the task scheduling problem has multiple possible solutions or options. Search space reduction for speedy processing of this algorithm is carried out through greedy strategy. Performance in terms of better results is ensured through computationally expensive recursive and iterative greedy approaches. Any optimal scheduling schemes show better performance in heterogeneous resources using CPU or GPU.

Parallel evolution of non-homologous isofunctional enzymes in methionine biosynthesis

Experimental validation of enzyme function is crucial for genome interpretation, but it remains challenging because it cannot be scaled up to accommodate the constant accumulation of genome sequences. We tackled this issue for the MetA and MetX enzyme families, phylogenetically unrelated families of acyl-L-homoserine transferases involved in L-methionine biosynthesis. Members of these families are prone to incorrect annotation because MetX and MetA enzymes are assumed to always use acetyl-CoA and succinyl-CoA, respectively. We determined the enzymatic activities of 100 enzymes from diverse species, and interpreted the results by structural classification of active sites based on protein structure modeling. We predict that >60% of the 10,000 sequences from these families currently present in databases are incorrectly annotated, and suggest that acetyl-CoA was originally the sole substrate of these isofunctional enzymes, which evolved to use exclusively succinyl-CoA in the most recent bacteria. We also uncovered a divergent subgroup of MetX enzymes in fungi that participate only in L-cysteine biosynthesis as O-succinyl-L-serine transferases.

Molecular function prediction for a family exhibiting evolutionary tendencies toward substrate specificity swapping: recurrence of tyrosine aminotransferase activity in the Iα subfamily

The subfamily Iα aminotransferases are typically categorized as having narrow specificity toward carboxylic amino acids (AATases), or broad specificity that includes aromatic amino acid substrates (TATases). Because of their general role in central metabolism and, more specifically, their association with liver-related diseases in humans, this subfamily is biologically interesting. The substrate specificities for only a few members of this subfamily have been reported, and the reliable prediction of substrate specificity from protein sequence has remained elusive. In this study, a diverse set of aminotransferases was chosen for characterization based on a scoring system that measures the sequence divergence of the active site. The enzymes that were experimentally characterized include both narrow-specificity AATases and broad-specificity TATases, as well as AATases with broader-specificity and TATases with narrower-specificity than the previously known family members. Molecular function and phylogenetic analyses underscored the complexity of this family's evolution as the TATase function does not follow a single evolutionary thread, but rather appears independently multiple times during the evolution of the subfamily. The additional functional characterizations described in this article, alongside a detailed sequence and phylogenetic analysis, provide some novel clues to understanding the evolutionary mechanisms at work in this family.

Dimerization and heme binding are conserved in amphibian and starfish homologues of the microRNA processing protein DGCR8

Human DiGeorge Critical Region 8 (DGCR8) is an essential microRNA (miRNA) processing factor that is activated via direct interaction with Fe(III) heme. In order for DGCR8 to bind heme, it must dimerize using a dimerization domain embedded within its heme-binding domain (HBD). We previously reported a crystal structure of the dimerization domain from human DGCR8, which demonstrated how dimerization results in the formation of a surface important for association with heme. Here, in an attempt to crystallize the HBD, we search for DGCR8 homologues and show that DGCR8 from Patiria miniata (bat star) also binds heme. The extinction coefficients (ε) of DGCR8-heme complexes are determined; these values are useful for biochemical analyses and allow us to estimate the heme occupancy of DGCR8 proteins. Additionally, we present the crystal structure of the Xenopus laevis dimerization domain. The structure is very similar to that of human DGCR8. Our results indicate that dimerization and heme binding are evolutionarily conserved properties of DGCR8 homologues not only in vertebrates, but also in at least some invertebrates.

Structural and functional analysis of a new subfamily of glycosyltransferases required for glycosylation of serine-rich streptococcal adhesins

Serine-rich repeat glycoproteins (SRRPs) are a growing family of bacterial adhesins found in many streptococci and staphylococci; they play important roles in bacterial biofilm formation and pathogenesis. Glycosylation of this family of adhesins is essential for their biogenesis. A glucosyltransferase (Gtf3) catalyzes the second step of glycosylation of a SRRP (Fap1) from an oral streptococcus, Streptococcus parasanguinis. Although Gtf3 homologs are highly conserved in SRRP-containing streptococci, they share minimal homology with functionally known glycosyltransferases. We report here the 2.3 Å crystal structure of Gtf3. The structural analysis indicates that Gtf3 forms a tetramer and shares significant structural homology with glycosyltransferases from GT4, GT5, and GT20 subfamilies. Combining crystal structural analysis with site-directed mutagenesis and in vitro glycosyltransferase assays, we identified residues that are required for UDP- or UDP-glucose binding and for oligomerization of Gtf3 and determined their contribution to the enzymatic activity of Gtf3. Further in vivo studies revealed that the critical amino acid residues identified by the structural analysis are crucial for Fap1 glycosylation in S. parasanguinis in vivo. Moreover, Gtf3 homologs from other streptococci were able to rescue the gtf3 knock-out mutant of S. parasanguinis in vivo and catalyze the sugar transfer to the modified SRRP substrate in vitro, demonstrating the importance and conservation of the Gtf3 homologs in glycosylation of SRRPs. As the Gtf3 homologs only exist in SRRP-containing streptococci, we conclude that the Gtf3 homologs represent a unique subfamily of glycosyltransferases.

Structure of the dimerization domain of DiGeorge critical region 8

Maturation of microRNAs (miRNAs, approximately 22nt) from long primary transcripts [primary miRNAs (pri-miRNAs)] is regulated during development and is altered in diseases such as cancer. The first processing step is a cleavage mediated by the Microprocessor complex containing the Drosha nuclease and the RNA-binding protein DiGeorge critical region 8 (DGCR8). We previously reported that dimeric DGCR8 binds heme and that the heme-bound DGCR8 is more active than the heme-free form. Here, we identified a conserved dimerization domain in DGCR8. Our crystal structure of this domain (residues 298-352) at 1.7 A resolution demonstrates a previously unknown use of a WW motif as a platform for extensive dimerization interactions. The dimerization domain of DGCR8 is embedded in an independently folded heme-binding domain and directly contributes to association with heme. Heme-binding-deficient DGCR8 mutants have reduced pri-miRNA processing activity in vitro. Our study provides structural and biochemical bases for understanding how dimerization and heme binding of DGCR8 may contribute to regulation of miRNA biogenesis.

Apo and ligand-bound structures of ModA from the archaeon Methanosarcina acetivorans

The trace-element oxyanion molybdate, which is required for the growth of many bacterial and archaeal species, is transported into the cell by an ATP-binding cassette (ABC) transporter superfamily uptake system called ModABC. ModABC consists of the ModA periplasmic solute-binding protein, the integral membrane-transport protein ModB and the ATP-binding and hydrolysis cassette protein ModC. In this study, X-ray crystal structures of ModA from the archaeon Methanosarcina acetivorans (MaModA) have been determined in the apoprotein conformation at 1.95 and 1.69 A resolution and in the molybdate-bound conformation at 2.25 and 2.45 A resolution. The overall domain structure of MaModA is similar to other ModA proteins in that it has a bilobal structure in which two mixed alpha/beta domains are linked by a hinge region. The apo MaModA is the first unliganded archaeal ModA structure to be determined: it exhibits a deep cleft between the two domains and confirms that upon binding ligand one domain is rotated towards the other by a hinge-bending motion, which is consistent with the 'Venus flytrap' model seen for bacterial-type periplasmic binding proteins. In contrast to the bacterial ModA structures, which have tetrahedral coordination of their metal substrates, molybdate-bound MaModA employs octahedral coordination of its substrate like other archaeal ModA proteins.

Algorithms for multiple protein structure alignment and structure-derived multiple sequence alignment

Primary amino acid content and the geometry of the folded protein 3D structure are major parameters of protein function. During the course of evolution the protein 3D structure is more preserved than its primary sequence. Thus, analysis of protein structures is expected to lead to a deep insight into protein function. Recognition of a structural core common to a set of protein structures serves as a basic tool for the studies of protein evolution and classification, analysis of similar structural motifs and functional binding sites, and for homology modeling and threading. In this chapter, we discuss several biologically related computational aspects of the multiple structure alignment and propose a method that provides solutions to these problems. Finally, we address the problem of structure-based multiple sequence alignment and propose an optimization method that unifies primary sequence and 3D structure information.

Expresso: automatic incorporation of structural information in multiple sequence alignments using 3D-Coffee

Expresso is a multiple sequence alignment server that aligns sequences using structural information. The user only needs to provide sequences. The server runs BLAST to identify close homologues of the sequences within the PDB database. These PDB structures are used as templates to guide the alignment of the original sequences using structure-based sequence alignment methods like SAP or Fugue. The final result is a multiple sequence alignment of the original sequences based on the structural information of the templates. An advanced mode makes it possible to either upload private structures or specify which PDB templates should be used to model each sequence. Providing the suitable structural information is available, Expresso delivers sequence alignments with accuracy comparable with structure-based alignments. The server is available on http://www.tcoffee.org/.

SSALN: an alignment algorithm using structure-dependent substitution matrices and gap penalties learned from structurally aligned protein pairs

In template-based modeling of protein structures, the generation of the alignment between the target and the template is a critical step that significantly affects the accuracy of the final model. This paper proposes an alignment algorithm SSALN that learns substitution matrices and position-specific gap penalties from a database of structurally aligned protein pairs. In addition to the amino acid sequence information, secondary structure and solvent accessibility information of a position are used to derive substitution scores and position-specific gap penalties. In a test set of CASP5 targets, SSALN outperforms sequence alignment methods such as a Smith-Waterman algorithm with BLOSUM50 and PSI_BLAST. SSALN also generates better alignments than PSI_BLAST in the CASP6 test set. LOOPP server prediction based on an SSALN alignment is ranked the best for target T0280_1 in CASP6. SSALN is also compared with several threading methods and sequence alignment methods on the ProSup benchmark. SSALN has the highest alignment accuracy among the methods compared. On the Fischer's benchmark, SSALN performs better than CLUSTALW and GenTHREADER, and generates more alignments with accuracy >50%, >60% or >70% than FUGUE, but fewer alignments with accuracy >80% than FUGUE. All the supplemental materials can be found at http://www.cs.cornell.edu/ approximately jianq/research.htm.

Optimization of multiple-sequence alignment based on multiple-structure alignment

Routinely used multiple-sequence alignment methods use only sequence information. Consequently, they may produce inaccurate alignments. Multiple-structure alignment methods, on the other hand, optimize structural alignment by ignoring sequence information. Here, we present an optimization method that unifies sequence and structure information. The alignment score is based on standard amino acid substitution probabilities combined with newly computed three-dimensional structure alignment probabilities. The advantage of our alignment scheme is in its ability to produce more accurate multiple alignments. We demonstrate the usefulness of the method in three applications: 1) computing more accurate multiple-sequence alignments, 2) analyzing protein conformational changes, and 3) computation of amino acid structure-sequence conservation with application to protein-protein docking prediction. The method is available at http://bioinfo3d.cs.tau.ac.il/staccato/.

Structure of the cytochrome complex SoxXA of Paracoccus pantotrophus, a heme enzyme initiating chemotrophic sulfur oxidation

The sulfur-oxidizing enzyme system (Sox) of the chemotroph Paracoccus pantotrophus is composed of several proteins, which together oxidize hydrogen sulfide, sulfur, thiosulfate or sulfite and transfers the gained electrons to the respiratory chain. The hetero-dimeric cytochrome c complex SoxXA functions as heme enzyme and links covalently the sulfur substrate to the thiol of the cysteine-138 residue of the SoxY protein of the SoxYZ complex. Here, we report the crystal structure of the c-type cytochrome complex SoxXA. The structure could be solved by molecular replacement and refined to a resolution of 1.9A identifying the axial heme-iron coordination involving an unusual Cys-251 thiolate of heme2. Distance measurements between the three heme groups provide deeper insight into the electron transport inside SoxXA and merge in a better understanding of the initial step of the aerobic sulfur oxidation process in chemotrophic bacteria.

Crystal structure of a RuBisCO-like protein from the green sulfur bacterium Chlorobium tepidum

Ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCO) catalyzes the incorporation of atmospheric CO(2) into ribulose 1,5-bisphosphate (RuBP). RuBisCOs are classified into four forms based on sequence similarity: forms I, II and III are bona fide RuBisCOs; form IV, also called the RuBisCO-like protein (RLP), lacks several of the substrate binding and catalytic residues and does not catalyze RuBP-dependent CO(2) fixation in vitro. To contribute to understanding the function of RLPs, we determined the crystal structure of the RLP from Chlorobium tepidum. The overall structure of the RLP is similar to the structures of the three other forms of RuBisCO; however, the active site is distinct from those of bona fide RuBisCOs and suggests that the RLP is possibly capable of catalyzing enolization but not carboxylation. Bioinformatic analysis of the protein functional linkages suggests that this RLP coevolved with enzymes of the bacteriochlorophyll biosynthesis pathway and may be involved in processes related to photosynthesis.

Crystal structure of the Mycobacterium tuberculosis dUTPase: insights into the catalytic mechanism

The structure of Mycobacterium tuberculosis dUTP nucleotidohydrolase (dUTPase) has been determined at 1.3 Angstrom resolution in complex with magnesium ion and the non-hydrolyzable substrate analog, alpha,beta-imido dUTP. dUTPase is an enzyme essential for depleting potentially toxic concentrations of dUTP in the cell. Given the importance of its biological role, it has been proposed that inhibiting M.tuberculosis dUTPase might be an effective means to treat tuberculosis infection in humans. The crystal structure presented here offers some insight into the potential for designing a specific inhibitor of the M.tuberculosis dUTPase enzyme. The structure also offers new insights into the mechanism of dUTP hydrolysis by providing an accurate representation of the enzyme-substrate complex in which both the metal ion and dUTP analog are included. The structure suggests that inclusion of a magnesium ion is important for stabilizing the position of the alpha-phosphorus for an in-line nucleophilic attack. In the absence of magnesium, the alpha-phosphate of dUTP can have either of the two positions which differ by 4.5 Angstrom. A transiently ordered C-terminal loop further assists catalysis by shielding the general base, Asp83, from solvent thus elevating its pK(a) so that it might in turn activate a tightly bound water molecule for nucleophilic attack. The metal ion coordinates alpha, beta, and gamma phosphate groups with tridentate geometry identical with that observed in the crystal structure of DNA polymerase beta complexed with magnesium and dNTP analog, revealing some common features in catalytic mechanism.