PROMALS3D: a tool for multiple protein sequence and structure alignments

GalaxyTBM: template-based modeling by building a reliable core and refining unreliable local regions

Background

Protein structures can be reliably predicted by template-based modeling (TBM) when experimental structures of homologous proteins are available. However, it is challenging to obtain structures more accurate than the single best templates by either combining information from multiple templates or by modeling regions that vary among templates or are not covered by any templates.

Results

We introduce GalaxyTBM, a new TBM method in which the more reliable core region is modeled first from multiple templates and less reliable, variable local regions, such as loops or termini, are then detected and re-modeled by an ab initio method. This TBM method is based on “Seok-server,” which was tested in CASP9 and assessed to be amongst the top TBM servers. The accuracy of the initial core modeling is enhanced by focusing on more conserved regions in the multiple-template selection and multiple sequence alignment stages. Additional improvement is achieved by ab initio modeling of up to 3 unreliable local regions in the fixed framework of the core structure. Overall, GalaxyTBM reproduced the performance of Seok-server, with GalaxyTBM and Seok-server resulting in average GDT-TS of 68.1 and 68.4, respectively, when tested on 68 single-domain CASP9 TBM targets. For application to multi-domain proteins, GalaxyTBM must be combined with domain-splitting methods.

Conclusion

Application of GalaxyTBM to CASP9 targets demonstrates that accurate protein structure prediction is possible by use of a multiple-template-based approach, and ab initio modeling of variable regions can further enhance the model quality.

Cysteine-rich domains related to Frizzled receptors and Hedgehog-interacting proteins

Frizzled and Smoothened are homologous seven-transmembrane proteins functioning in the Wnt and Hedgehog signaling pathways, respectively. They harbor an extracellular cysteine-rich domain (FZ-CRD), a mobile evolutionary unit that has been found in a number of other metazoan proteins and Frizzled-like proteins in Dictyostelium. Domains distantly related to FZ-CRDs, in Hedgehog-interacting proteins (HHIPs), folate receptors and riboflavin-binding proteins (FRBPs), and Niemann-Pick Type C1 proteins (NPC1s), referred to as HFN-CRDs, exhibit similar structures and disulfide connectivity patterns compared with FZ-CRDs. We used computational analyses to expand the homologous set of FZ-CRDs and HFN-CRDs, providing a better understanding of their evolution and classification. First, FZ-CRD-containing proteins with various domain compositions were identified in several major eukaryotic lineages including plants and Chromalveolata, revealing a wider phylogenetic distribution of FZ-CRDs than previously recognized. Second, two new and distinct groups of highly divergent FZ-CRDs were found by sensitive similarity searches. One of them is present in the calcium channel component Mid1 in fungi and the uncharacterized FAM155 proteins in metazoans. Members of the other new FZ-CRD group occur in the metazoan-specific RECK (reversion-inducing-cysteine-rich protein with Kazal motifs) proteins that are putative tumor suppressors acting as inhibitors of matrix metalloproteases. Finally, sequence and three-dimensional structural comparisons helped us uncover a divergent HFN-CRD in glypicans, which are important morphogen-binding heparan sulfate proteoglycans. Such a finding reinforces the evolutionary ties between the Wnt and Hedgehog signaling pathways and underscores the importance of gene duplications in creating essential signaling components in metazoan evolution.

Prediction of heparin binding sites in bone morphogenetic proteins (BMPs)

Heparin is a glycosaminoglycan known to bind bone morphogenetic proteins (BMPs) and the growth and differentiation factors (GDFs) and has strong and variable effects on BMP osteogenic activity. In this paper we report our predictions of the likely heparin binding sites for BMP-2 and 14. The N-terminal sequences upstream of TGF-β-type cysteine-knot domains in BMP-2, 7 and 14 contain the basic residues arginine and lysine, which are key components of the heparin/HS-binding sites, with these residues being highly non-conserved. Importantly, evolutionary conserved surfaces on the beta sheets are required for interactions with receptors and antagonists. Furthermore, BMP-2 has electropositive surfaces on two sides compared to BMP-7 and BMP-14. Molecular docking simulations suggest the presence of high and low affinity binding sites in dimeric BMP-2. Histidines were found to play a role in the interactions of BMP-2 with heparin; however, a pK(a) analysis suggests that histidines are likely not protonated. This is indicative that interactions of BMP-2 with heparin do not require acidic pH. Taken together, non-conserved amino acid residues in the N-terminus and residues protruding from the beta sheet (not overlapping with the receptor binding sites and the dimeric interface) and not C-terminal are found to be important for heparin-BMP interactions.

A two-staged model of Na+ exclusion in rice explained by 3D modeling of HKT transporters and alternative splicing

The HKT family of Na(+) and Na(+)/K(+) transporters is implicated in plant salinity tolerance. Amongst these transporters, the cereal HKT1;4 and HKT1;5 are responsible for Na(+) exclusion from photosynthetic tissues, a key mechanism for plant salinity tolerance. It has been suggested that Na(+) is retrieved from the xylem transpiration stream either in the root or the leaf sheath, protecting the leaf blades from excessive Na(+) accumulation. However, direct evidence for this scenario is scarce. Comparative modeling and evaluation of rice (Oryza sativa) HKT-transporters based on the recent crystal structure of the bacterial TrkH K(+) transporter allowed to reconcile transcriptomic and physiological data. For OsHKT1;5, both transcript abundance and protein structural features within the selectivity filter could control shoot Na(+) accumulation in a range of rice varieties. For OsHKT1;4, alternative splicing of transcript and the anatomical complexity of the sheath needed to be taken into account. Thus, Na(+) accumulation in a specific leaf blade seems to be regulated by abundance of a correctly spliced OsHKT1;4 transcript in a corresponding sheath. Overall, allelic variation of leaf blade Na(+) accumulation can be explained by a complex interplay of gene transcription, alternative splicing and protein structure.

Widespread occurrence of N-terminal acylation in animal globins and possible origin of respiratory globins from a membrane-bound ancestor

Proteins of the (hemo-)globin superfamily have been identified in many different animals but also occur in plants, fungi, and bacteria. Globins are renowned for their ability to store and to transport oxygen, but additional globin functions such as sensing, signaling, and detoxification have been proposed. Recently, we found that the zebrafish globin X protein is myristoylated and palmitoylated at its N-terminus. The addition of fatty acids results in an association with the cellular membranes, suggesting a previously unrecognized globin function. In this study, we show that N-terminal acylation likely occurs in globin proteins from a broad range of phyla. An N-terminal myristoylation site was identified in 90 nonredundant globins from Chlorophyta, Heterokontophyta, Cnidaria, Mollusca, Arthropoda, Nematoda, Echinodermata, Hemichordata, and Chordata (including Cephalochordata), of which 66 proteins carry an additional palmitoylation site. Bayesian phylogenetic analyses identified five major globin families, which may mirror the ancient globin diversity of the Metazoa. Globin X-like proteins form two related clades, which diverged before the radiation of the Eumetazoa. Vertebrate hemoglobin (Hb), myoglobin, cytoglobin, globin E, and globin Y form a strongly supported common clade, which is the sister group of a clade consisting of invertebrate Hbs and relatives. The N-terminally acylated globins do not form a single monophyletic group but are distributed to four distinct clades. This pattern may be either explained by multiple introduction of an N-terminal acylation site into distinct globin lineages or by the origin of animal respiratory globins from a membrane-bound ancestor. Similarly, respiratory globins were not monophyletic. This suggests that respiratory globins might have emerged independently several times and that the early metazoan globins might have been associated with a membrane and carried out a function that was related to lipid protection or signaling.

A position-specific distance-dependent statistical potential for protein structure and functional study

Although studied extensively, designing highly accurate protein energy potential is still challenging. A lot of knowledge-based statistical potentials are derived from the inverse of the Boltzmann law and consist of two major components: observed atomic interacting probability and reference state. These potentials mainly distinguish themselves in the reference state and use a similar simple counting method to estimate the observed probability, which is usually assumed to correlate with only atom types. This article takes a rather different view on the observed probability and parameterizes it by the protein sequence profile context of the atoms and the radius of the gyration, in addition to atom types. Experiments confirm that our position-specific statistical potential outperforms currently the popular ones in several decoy discrimination tests. Our results imply that, in addition to reference state, the observed probability also makes energy potentials different and evolutionary information greatly boost performance of energy potentials.

Alternaria alternata allergen Alt a 1: a unique β-barrel protein dimer found exclusively in fungi

BACKGROUND

Alternaria species is one of the most common molds associated with allergic diseases, and 80% of Alternaria species-sensitive patients produce IgE antibodies to a major protein allergen, Alt a 1. The structure and function of Alt a 1 is unknown.

OBJECTIVE

We sought to obtain a high-resolution structure of Alt a 1 using x-ray crystallography and to investigate structural relationships between Alt a 1 and other allergens and proteins reported in the Protein Data Bank.

METHODS

X-ray crystallography was used to determine the structure of Alt a 1 by using a custom-designed set of crystallization conditions. An initial Alt a 1 model was determined by the application of a Ta(6)Br(12)(2+) cluster and single-wavelength anomalous diffraction. Bioinformatic analyses were used to compare the Alt a 1 sequence and structure with that of other proteins.

RESULTS

Alt a 1 is a unique β-barrel comprising 11 β-strands and forms a "butterfly-like" dimer linked by a single disulfide bond with a large (1345 Å(2)) dimer interface. Intramolecular disulfide bonds are conserved among Alt a 1 homologs. Currently, the Alt a 1 structure has no equivalent in the Protein Data Bank. Bioinformatics analyses suggest that the structure is found exclusively in fungi. Four previously reported putative IgE-binding peptides have been located on the Alt a 1 structure.

CONCLUSIONS

Alt a 1 has a unique, dimeric β-barrel structure that appears to define a new protein family with unknown function found exclusively in fungi. The location of IgE antibody-binding epitopes is in agreement with the structural analysis of Alt a 1. The Alt a 1 structure will allow mechanistic structure/function studies and immunologic studies directed toward new forms of immunotherapy for Alternaria species-sensitive allergic patients.

Secreted kinase phosphorylates extracellular proteins that regulate biomineralization

Protein phosphorylation is a fundamental mechanism regulating nearly every aspect of cellular life. Several secreted proteins are phosphorylated, but the kinases responsible are unknown. We identified a family of atypical protein kinases that localize within the Golgi apparatus and are secreted. Fam20C appears to be the Golgi casein kinase that phosphorylates secretory pathway proteins within S-x-E motifs. Fam20C phosphorylates the caseins and several secreted proteins implicated in biomineralization, including the small integrin-binding ligand, N-linked glycoproteins (SIBLINGs). Consequently, mutations in Fam20C cause an osteosclerotic bone dysplasia in humans known as Raine syndrome. Fam20C is thus a protein kinase dedicated to the phosphorylation of extracellular proteins.

GalaxyWEB server for protein structure prediction and refinement

Three-dimensional protein structures provide invaluable information for understanding and regulating biological functions of proteins. The GalaxyWEB server predicts protein structure from sequence by template-based modeling and refines loop or terminus regions by ab initio modeling. This web server is based on the method tested in CASP9 (9th Critical Assessment of techniques for protein Structure Prediction) as ‘Seok-server’, which was assessed to be among top performing template-based modeling servers. The method generates reliable core structures from multiple templates and re-builds unreliable loops or termini by using an optimization-based refinement method. In addition to structure prediction, a user can also submit a refinement only job by providing a starting model structure and locations of loops or termini to refine. The web server can be freely accessed at http://galaxy.seoklab.org/.

Structural analyses of the Slm1-PH domain demonstrate ligand binding in the non-canonical site

Background

Pleckstrin homology (PH) domains are common membrane-targeting modules and their best characterized ligands are a set of important signaling lipids that include phosphatidylinositol phosphates (PtdInsPs). PH domains recognize PtdInsPs through two distinct mechanisms that use different binding pockets on opposite sides of the β-strands 1 and 2: i) a canonical binding site delimited by the β1-β2 and β3-β4loops and ii) a non-canonical binding site bordered by the β1-β2 and β5-β6loops. The PH domain-containing protein Slm1 from budding yeast Saccharomyces cerevisiae is required for actin cytoskeleton polarization and cell growth. We recently reported that this PH domain binds PtdInsPs and phosphorylated sphingolipids in a cooperative manner.

Principal Findings

To study the structural basis for the Slm1-PH domain (Slm1-PH) specificity, we co-crystallized this domain with different soluble compounds that have structures analogous to anionic lipid head groups of reported Slm1 ligands: inositol 4-phosphate, which mimics phosphatidylinositol-4-phosphate (PtdIns(4)P), and phosphoserine as a surrogate for dihydrosphingosine 1-phosphate (DHS1-P). We found electron densities for the ligands within the so-called non-canonical binding site. An additional positively charged surface that contacts a phosphate group was identified next to the canonical binding site.

Conclusions

Our results suggest that Slm1-PH utilizes a non-canonical binding site to bind PtdInsPs, similar to that described for the PH domains of β-spectrin, Tiam1 and ArhGAP9. Additionally, Slm1-PH may have retained an active canonical site. We propose that the presence of both a canonical and a non-canonical binding pocket in Slm1-PH may account for the cooperative binding to PtdInsPs and DHS-1P.

Highly dynamic exon shuffling in candidate pathogen receptors ... what if brown algae were capable of adaptive immunity?

Pathogen recognition is the first step of immune reactions. In animals and plants, direct or indirect pathogen recognition is often mediated by a wealth of fast-evolving receptors, many of which contain ligand-binding and signal transduction domains, such as leucine-rich or tetratricopeptide repeat (LRR/TPR) and NB-ARC domains, respectively. In order to identify candidates potentially involved in algal defense, we mined the genome of the brown alga Ectocarpus siliculosus for homologues of these genes and assessed the evolutionary pressures acting upon them. We thus annotated all Ectocarpus LRR-containing genes, in particular an original group of LRR-containing GTPases of the ROCO family, and 24 NB-ARC-TPR proteins. They exhibit high birth and death rates, while a diversifying selection is acting on their LRR (respectively TPR) domain, probably affecting the ligand-binding specificities. Remarkably, each repeat is encoded by an exon, and the intense exon shuffling underpins the variability of LRR and TPR domains. We conclude that the Ectocarpus ROCO and NB-ARC-TPR families are excellent candidates for being involved in recognition/transduction events linked to immunity. We further hypothesize that brown algae may generate their immune repertoire via controlled somatic recombination, so far only known from the vertebrate adaptive immune systems.

Application of DEN refinement and automated model building to a difficult case of molecular-replacement phasing: the structure of a putative succinyl-diaminopimelate desuccinylase from Corynebacterium glutamicum

Phasing by molecular replacement remains difficult for targets that are far from the search model or in situations where the crystal diffracts only weakly or to low resolution. Here, the process of determining and refining the structure of Cgl1109, a putative succinyl-diaminopimelate desuccinylase from Corynebacterium glutamicum, at ∼3 Å resolution is described using a combination of homology modeling with MODELLER, molecular-replacement phasing with Phaser, deformable elastic network (DEN) refinement and automated model building using AutoBuild in a semi-automated fashion, followed by final refinement cycles with phenix.refine and Coot. This difficult molecular-replacement case illustrates the power of including DEN restraints derived from a starting model to guide the movements of the model during refinement. The resulting improved model phases provide better starting points for automated model building and produce more significant difference peaks in anomalous difference Fourier maps to locate anomalous scatterers than does standard refinement. This example also illustrates a current limitation of automated procedures that require manual adjustment of local sequence misalignments between the homology model and the target sequence.

The structure and selectivity of the SR protein SRSF2 RRM domain with RNA

SRSF2 is a prototypical SR protein which plays important roles in the alternative splicing of pre-mRNA. It has been shown to be involved in regulatory pathways for maintaining genomic stability and play important roles in regulating key receptors in the heart. We report here the solution structure of the RNA recognition motifs (RRM) domain of free human SRSF2 (residues 9-101). Compared with other members of the SR protein family, SRSF2 structure has a longer L3 loop region. The conserved aromatic residue in the RNP2 motif is absent in SRSF2. Calorimetric titration shows that the RNA sequence 5'AGCAGAGUA3' binds SRSF2 with a K(d) of 61 ± 1 nM and a 1:1 stoichiometry. NMR and mutagenesis experiments reveal that for SFSF2, the canonical β1 and β3 interactions are themselves not sufficient for effective RNA binding; the additional loop L3 is crucial for RNA complex formation. A comparison is made between the structures of SRSF2-RNA complex with other known RNA complexes of SR proteins. We conclude that interactions involving the L3 loop, N- and C-termini of the RRM domain are collectively important for determining selectivity between the protein and RNA.

Requirements for efficient correction of ΔF508 CFTR revealed by analyses of evolved sequences

Misfolding of ΔF508 cystic fibrosis (CF) transmembrane conductance regulator (CFTR) underlies pathology in most CF patients. F508 resides in the first nucleotide-binding domain (NBD1) of CFTR near a predicted interface with the fourth intracellular loop (ICL4). Efforts to identify small molecules that restore function by correcting the folding defect have revealed an apparent efficacy ceiling. To understand the mechanistic basis of this obstacle, positions statistically coupled to 508, in evolved sequences, were identified and assessed for their impact on both NBD1 and CFTR folding. The results indicate that both NBD1 folding and interaction with ICL4 are altered by the ΔF508 mutation and that correction of either individual process is only partially effective. By contrast, combination of mutations that counteract both defects restores ΔF508 maturation and function to wild-type levels. These results provide a mechanistic rationale for the limited efficacy of extant corrector compounds and suggest approaches for identifying compounds that correct both defective steps.

Low-resolution structure of Drosophila translin

Crystals of native Drosophila melanogaster translin diffracted to 7 Å resolution. Reductive methylation of the protein improved crystal quality. The native and methylated proteins showed similar profiles in size-exclusion chromatography analyses but the methylated protein displayed reduced DNA-binding activity. Crystals of the methylated protein diffracted to 4.2 Å resolution at BM14 of the ESRF synchrotron. Crystals with 49% solvent content belonged to monoclinic space group P21 with eight protomers in the asymmetric unit. Only 2% of low-resolution structures with similar low percentage solvent content were found in the PDB. The crystal structure, solved by molecular replacement method, refined to R work (R free) of 0.24 (0.29) with excellent stereochemistry. The crystal structure clearly shows that drosophila protein exists as an octamer, and not as a decamer as expected from gel-filtration elution profiles. The similar octameric quaternary fold in translin orthologs and in translin-TRAX complexes suggests an up-down dimer as the basic structural subunit of translin-like proteins. The drosophila oligomer displays asymmetric assembly and increased radius of gyration that accounts for the observed differences between the elution profiles of human and drosophila proteins on gel-filtration columns. This study demonstrates clearly that low-resolution X-ray structure can be useful in understanding complex biological oligomers.

Membrane topology and predicted RNA-binding function of the 'early responsive to dehydration (ERD4)' plant protein

Functional annotation of uncharacterized genes is the main focus of computational methods in the post genomic era. These tools search for similarity between proteins on the premise that those sharing sequence or structural motifs usually perform related functions, and are thus particularly useful for membrane proteins. Early responsive to dehydration (ERD) genes are rapidly induced in response to dehydration stress in a variety of plant species. In the present work we characterized function of Brassica juncea ERD4 gene using computational approaches. The ERD4 protein of unknown function possesses ubiquitous DUF221 domain (residues 312-634) and is conserved in all plant species. We suggest that the protein is localized in chloroplast membrane with at least nine transmembrane helices. We detected a globular domain of 165 amino acid residues (183-347) in plant ERD4 proteins and expect this to be posited inside the chloroplast. The structural-functional annotation of the globular domain was arrived at using fold recognition methods, which suggested in its sequence presence of two tandem RNA-recognition motif (RRM) domains each folded into βαββαβ topology. The structure based sequence alignment with the known RNA-binding proteins revealed conservation of two non-canonical ribonucleoprotein sub-motifs in both the putative RNA-recognition domains of the ERD4 protein. The function of highly conserved ERD4 protein may thus be associated with its RNA-binding ability during the stress response. This is the first functional annotation of ERD4 family of proteins that can be useful in designing experiments to unravel crucial aspects of stress tolerance mechanism.

Structural characterization of Helicobacter pylori dethiobiotin synthetase reveals differences between family members

Dethiobiotin synthetase (DTBS) is involved in the biosynthesis of biotin in bacteria, fungi, and plants. As humans lack this pathway, DTBS is a promising antimicrobial drug target. We determined structures of DTBS from Helicobacter pylori (hpDTBS) bound with cofactors and a substrate analog, and described its unique characteristics relative to other DTBS proteins. Comparison with bacterial DTBS orthologs revealed considerable structural differences in nucleotide recognition. The C-terminal region of DTBS proteins, which contains two nucleotide-recognition motifs, differs greatly among DTBS proteins from different species. The structure of hpDTBS revealed that this protein is unique and does not contain a C-terminal region containing one of the motifs. The single nucleotide-binding motif in hpDTBS is similar to its counterpart in GTPases; however, isothermal titration calorimetry binding studies showed that hpDTBS has a strong preference for ATP. The structural determinants of ATP specificity were assessed with X-ray crystallographic studies of hpDTBS·ATP and hpDTBS·GTP complexes. The unique mode of nucleotide recognition in hpDTBS makes this protein a good target for H. pylori-specific inhibitors of the biotin synthesis pathway.

The intracellular Ig fold: a robust protein scaffold for the engineering of molecular recognition

Protein scaffolds that support molecular recognition have multiple applications in biotechnology. Thus, protein frames with robust structural cores but adaptable surface loops are in continued demand. Recently, notable progress has been made in the characterization of Ig domains of intracellular origin--in particular, modular components of the titin myofilament. These Ig belong to the I(intermediate)-type, are remarkably stable, highly soluble and undemanding to produce in the cytoplasm of Escherichia coli. Using the Z1 domain from titin as representative, we show that the I-Ig fold tolerates the drastic diversification of its CD loop, constituting an effective peptide display system. We examine the stability of CD-loop-grafted Z1-peptide chimeras using differential scanning fluorimetry, Fourier transform infrared spectroscopy and nuclear magnetic resonance and demonstrate that the introduction of bioreactive affinity binders in this position does not compromise the structural integrity of the domain. Further, the binding efficiency of the exogenous peptide sequences in Z1 is analyzed using pull-down assays and isothermal titration calorimetry. We show that an internally grafted, affinity FLAG tag is functional within the context of the fold, interacting with the anti-FLAG M2 antibody in solution and in affinity gel. Together, these data reveal the potential of the intracellular Ig scaffold for targeted functionalization.

SBAL: a practical tool to generate and edit structure-based amino acid sequence alignments

SUMMARY

Both alignment generation and visualization are important processes for producing biologically meaningful sequence alignments. Computational tools that combine reliable, automated and semi-automated approaches to produce secondary structure-based alignments with an appropriate visualization of the results are rare. We have developed SBAL, a tool to generate and edit secondary structure-based sequence alignments. It is easy to install and provides a user-friendly interface. Sequence alignments are displayed, with secondary structure assignments mapped to their corresponding regions in the sequence by using a simple colour scheme. The algorithm implemented for automated and semi-automated secondary structure-based alignment calculations shows a comparable performance to existing software.

AVAILABILITY AND IMPLEMENTATION

SBAL has been implemented in Java to provide cross-platform compatibility. SBAL is freely available to academic users at http://www.structuralchemistry.org/pcsb/. Users will be asked for their name, institution and email address. A manual can also be downloaded from this site. The software, manual and test sets are also available as supplementary material.

CONTACT

conan.wang@griffith.edu.au

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Structural, bioinformatic, and in vivo analyses of two Treponema pallidum lipoproteins reveal a unique TRAP transporter

Treponema pallidum, the bacterial agent of syphilis, is predicted to encode one tripartite ATP-independent periplasmic transporter (TRAP-T). TRAP-Ts typically employ a periplasmic substrate-binding protein (SBP) to deliver the cognate ligand to the transmembrane symporter. Herein, we demonstrate that the genes encoding the putative TRAP-T components from T. pallidum, tp0957 (the SBP), and tp0958 (the symporter), are in an operon with an uncharacterized third gene, tp0956. We determined the crystal structure of recombinant Tp0956; the protein is trimeric and perforated by a pore. Part of Tp0956 forms an assembly similar to those of "tetratricopeptide repeat" (TPR) motifs. The crystal structure of recombinant Tp0957 was also determined; like the SBPs of other TRAP-Ts, there are two lobes separated by a cleft. In these other SBPs, the cleft binds a negatively charged ligand. However, the cleft of Tp0957 has a strikingly hydrophobic chemical composition, indicating that its ligand may be substantially different and likely hydrophobic. Analytical ultracentrifugation of the recombinant versions of Tp0956 and Tp0957 established that these proteins associate avidly. This unprecedented interaction was confirmed for the native molecules using in vivo cross-linking experiments. Finally, bioinformatic analyses suggested that this transporter exemplifies a new subfamily of TPATs (TPR-protein-associated TRAP-Ts) that require the action of a TPR-containing accessory protein for the periplasmic transport of a potentially hydrophobic ligand(s).

Phenotype-specific adverse effects of XPD mutations on human prenatal development implicate impairment of TFIIH-mediated functions in placenta

Mutations in XPD (ERCC2), XPB (ERCC3), and TTD-A (GTF2H5), genes involved in nucleotide excision repair and transcription, can cause several disorders including trichothiodystrophy (TTD) and xeroderma pigmentosum (XP). In this study, we tested the hypothesis that mutations in the XPD gene affect placental development in a phenotype-specific manner. To test our hypothesis and decipher potential biologic mechanisms, we compared all XPD-associated TTD (n=43) and XP (n=37) cases reported in the literature with respect to frequencies of gestational complications. Our genetic epidemiologic investigations of TTD and XP revealed that the exact genetic abnormality was relevant to the mechanism leading to gestational complications such as preeclampsia. Through structural mapping, we localized the preeclampsia-associated mutations to a C-terminal motif and the helicase surfaces of XPD, most likely affecting XPD's binding to cdk-activating kinase (CAK) and p44 subunits of transcription factor (TF) IIH. Our results suggested a link between TTD- but not XP-associated XPD mutations, placental maldevelopment and risk of pregnancy complications, possibly due to impairment of TFIIH-mediated functions in placenta. Our findings highlight the importance of the fetal genotype in development of gestational complications, such as preeclampsia. Therefore, future studies of genetic associations of preeclampsia and other placental vascular complications may benefit from focusing on genetic variants within the fetal DNA.

Identification of domains and amino acids essential to the collagen galactosyltransferase activity of GLT25D1

Collagen is modified by hydroxylation and glycosylation of hydroxylysine residues. This glycosylation is initiated by the β1,O galactosyltransferases GLT25D1 and GLT25D2. The structurally similar protein cerebral endothelial cell adhesion molecule CEECAM1 was previously reported to be inactive when assayed for collagen glycosyltransferase activity. To address the cause of the absent galactosyltransferase activity, we have generated several chimeric constructs between the active human GLT25D1 and inactive human CEECAM1 proteins. The assay of these chimeric constructs pointed to a short central region and a large C-terminal region of CEECAM1 leading to the loss of collagen galactosyltransferase activity. Examination of the three DXD motifs of the active GLT25D1 by site-directed mutagenesis confirmed the importance of the first (amino acids 166–168) and second motif (amino acids 461–463) for enzymatic activity, whereas the third one was dispensable. Since the second DXD motif is incomplete in CEECAM1, we have restored the motif by introducing the substitution S461D. This change did not restore the activity of the C-terminal region, thereby showing that additional amino acids were required in this C-terminal region to confer enzymatic activity. Finally, we have introduced the substitution Q471R-V472M-N473Q-P474V in the CEECAM1-C-terminal construct, which is found in most animal GLT25D1 and GLT25D2 isoforms but not in CEECAM1. This substitution was shown to partially restore collagen galactosyltransferase activity, underlining its importance for catalytic activity in the C-terminal domain. Because multiple mutations in different regions of CEECAM1 contribute to the lack of galactosyltransferase activity, we deduced that CEECAM1 is functionally different from the related GLT25D1 protein.

MSACompro: protein multiple sequence alignment using predicted secondary structure, solvent accessibility, and residue-residue contacts

Background

Multiple Sequence Alignment (MSA) is a basic tool for bioinformatics research and analysis. It has been used essentially in almost all bioinformatics tasks such as protein structure modeling, gene and protein function prediction, DNA motif recognition, and phylogenetic analysis. Therefore, improving the accuracy of multiple sequence alignment is important for advancing many bioinformatics fields.

Results

We designed and developed a new method, MSACompro, to synergistically incorporate predicted secondary structure, relative solvent accessibility, and residue-residue contact information into the currently most accurate posterior probability-based MSA methods to improve the accuracy of multiple sequence alignments. The method is different from the multiple sequence alignment methods (e.g. 3D-Coffee) that use the tertiary structure information of some sequences since the structural information of our method is fully predicted from sequences. To the best of our knowledge, applying predicted relative solvent accessibility and contact map to multiple sequence alignment is novel. The rigorous benchmarking of our method to the standard benchmarks (i.e. BAliBASE, SABmark and OXBENCH) clearly demonstrated that incorporating predicted protein structural information improves the multiple sequence alignment accuracy over the leading multiple protein sequence alignment tools without using this information, such as MSAProbs, ProbCons, Probalign, T-coffee, MAFFT and MUSCLE. And the performance of the method is comparable to the state-of-the-art method PROMALS of using structural features and additional homologous sequences by slightly lower scores.

Conclusion

MSACompro is an efficient and reliable multiple protein sequence alignment tool that can effectively incorporate predicted protein structural information into multiple sequence alignment. The software is available at http://sysbio.rnet.missouri.edu/multicom_toolbox/.

Small heat shock proteins and α-crystallins: dynamic proteins with flexible functions

The small heat shock proteins (sHSPs) and the related α-crystallins (αCs) are virtually ubiquitous proteins that are strongly induced by a variety of stresses, but that also function constitutively in multiple cell types in many organisms. Extensive research has demonstrated that a majority of sHSPs and αCs can act as ATP-independent molecular chaperones by binding denaturing proteins and thereby protecting cells from damage due to irreversible protein aggregation. As a result of their diverse evolutionary history, their connection to inherited human diseases, and their novel protein dynamics, sHSPs and αCs are of significant interest to many areas of biology and biochemistry. However, it is increasingly clear that no single model is sufficient to describe the structure, function or mechanism of action of sHSPs and αCs. In this review, we discuss recent data that provide insight into the variety of structures of these proteins, their dynamic behavior, how they recognize substrates, and their many possible cellular roles.

Unpairing and gating: sequence-independent substrate recognition by FEN superfamily nucleases

Structure-specific 5'-nucleases form a superfamily of evolutionarily conserved phosphodiesterases that catalyse a precise incision of a diverse range of DNA and RNA substrates in a sequence-independent manner. Superfamily members, such as flap endonucleases, exonuclease 1, DNA repair protein XPG, endonuclease GEN1 and the 5'-3'-exoribonucleases, play key roles in many cellular processes such as DNA replication and repair, recombination, transcription, RNA turnover and RNA interference. In this review, we discuss recent results that highlight the conserved architectures and active sites of the structure-specific 5'-nucleases. Despite substrate diversity, a common gating mechanism for sequence-independent substrate recognition and incision emerges, whereby double nucleotide unpairing of substrates is required to access the active site.

Evolutionary and functional insights into Leishmania META1: evidence for lateral gene transfer and a role for META1 in secretion

Background

Leishmania META1 has for long been a candidate molecule for involvement in virulence: META1 transcript and protein are up-regulated in metacyclic Leishmania. Yet, how META1 contributes to virulence remains unclear. We sought insights into the possible functions of META1 by studying its evolutionary origins.

Results

Using multiple criteria including sequence similarity, nucleotide composition, phylogenetic analysis and selection pressure on gene sequence, we present evidence that META1 originated in trypanosomatids as a result of a lateral gene transfer of a bacterial heat-inducible protein, HslJ. Furthermore, within the Leishmania genome, META1 sequence is under negative selection pressure against change/substitution. Using homology modeling of Leishmania META1 based on solved NMR structure of HslJ, we show that META1 and HslJ share a similar structural fold. The best hit for other proteins with similar fold is MxiM, a protein involved in the type III secretion system in Shigella. The striking structural similarity shared by META1, HslJ and MxiM suggests a possibility of shared functions. Upon structural superposition with MxiM, we have observed a putative hydrophobic cavity in META1. Mutagenesis of select hydrophobic residues in this cavity affects the secretion of the secreted acid phosphatase (SAP), indicating META1's involvement in secretory processes in Leishmania.

Conclusions

Overall, this work uses an evolutionary biology approach, 3D-modeling and site-directed mutagenesis to arrive at new insights into functions of Leishmania META1.

The identification of antigenic proteins: 14-3-3 protein and propionyl-CoA carboxylase in Clonorchis sinensis

Clonorchis sinensis, the causative agent of clonorchiasis, is widespread in East and Southeast Asia, including China, Vietnam and the Republic of Korea. We identified antigenic proteins from adult C. sinensis liver flukes using immunoproteomic analysis. In this study, we found 23 candidate antigenic proteins with a pI in the range of 5.4-6.2 in total lysates of C. sinensis. The antigenic protein spots reacted against sera from clonorchiasis patients and were identified as cysteine proteases, glutathione transferases, gelsolin, propionyl-CoA carboxylase (PCC), prohibitin and 14-3-3 protein (14-3-3) using LC-coupled ESI-MS/MS and an EST database for C. sinensis. PCC and 14-3-3 were identified for the first time as serological antigens for the diagnosis of C. sinensis. To validate the antigenicity of PCC and 14-3-3, recombinant proteins were immunoblotted with sera from clonorchiasis patients. The structural, functional and immunological characteristics of the putative amino acid sequence were predicted by bioinformatics analysis. Our novel finding will contribute to the development of diagnostics for clonorchiasis. These results suggest that immunoproteomic approaches are valuable tools to identify antigens that could be used as targets for effective parasitic infection control strategies.

The evolution of function in strictosidine synthase-like proteins

The exponential growth of sequence data provides abundant information for the discovery of new enzyme reactions. Correctly annotating the functions of highly diverse proteins can be difficult, however, hindering use of this information. Global analysis of large superfamilies of related proteins is a powerful strategy for understanding the evolution of reactions by identifying catalytic commonalities and differences in reaction and substrate specificity, even when only a few members have been biochemically or structurally characterized. A comparison of >2500 sequences sharing the six-bladed β-propeller fold establishes sequence, structural, and functional links among the three subgroups of the functionally diverse N6P superfamily: the arylesterase-like and senescence marker protein-30/gluconolactonase/luciferin-regenerating enzyme-like (SGL) subgroups, representing enzymes that catalyze lactonase and related hydrolytic reactions, and the so-called strictosidine synthase-like (SSL) subgroup. Metal-coordinating residues were identified as broadly conserved in the active sites of all three subgroups except for a few proteins from the SSL subgroup, which have been experimentally determined to catalyze the quite different strictosidine synthase (SS) reaction, a metal-independent condensation reaction. Despite these differences, comparison of conserved catalytic features of the arylesterase-like and SGL enzymes with the SSs identified similar structural and mechanistic attributes between the hydrolytic reactions catalyzed by the former and the condensation reaction catalyzed by SS. The results also suggest that despite their annotations, the great majority of these >500 SSL sequences do not catalyze the SS reaction; rather, they likely catalyze hydrolytic reactions typical of the other two subgroups instead. This prediction was confirmed experimentally for one of these proteins.

Discovery of an orthoreovirus in the aborted fetus of a Steller sea lion (Eumetopias jubatus)

An aborted mid-gestational male Steller sea lion fetus with an attached placenta was recovered on the floor of an open floating capture trap located off Norris Rock near Denman Island, British Columbia. Viral culture of the placenta demonstrated cytopathic effect. Although no specific signal was obtained in microarray experiments using RNA obtained from viral culture, elution and sequence analysis revealed the presence of a reovirus. Complete genome pyrosequencing led to the identification of an orthoreovirus that we have tentatively named Steller sea lion reovirus (SSRV). Phylogenetic analysis revealed similarities between SSRV and orthoreoviruses of birds, bats and other mammals that suggests potential for interspecies transmission.

Aspartate 112 is the selectivity filter of the human voltage-gated proton channel

The ion selectivity of pumps and channels is central to their ability to perform a multitude of functions. Here we investigate the mechanism of the extraordinary selectivity of the human voltage gated proton channel¹, hH_V1. This selectivity is essential to its ability to regulate reactive oxygen species production by leukocytes^2–4, histamine secretion by basophils⁵, sperm capacitation⁶, and airway pH⁷. The most selective ion channel known, H_V1 shows no detectable permeability to other ions¹. Opposing classes of selectivity mechanisms postulate that (a) a titratable amino acid residue in the permeation pathway imparts proton selectivity^{1, 8–11}, or (b) water molecules “frozen” in a narrow pore conduct protons while excluding other ions¹². Here we identify Aspartate¹¹² as a crucial component of the selectivity filter of hH_V1. When a neutral amino acid replaced Asp¹¹², the mutant channel lost proton specificity and became anion selective or did not conduct. Only the glutamate mutant remained proton specific. Mutation of the nearby Asp¹⁸⁵ did not impair proton selectivity, suggesting that Asp¹¹² plays a unique role. Although histidine shuttles protons in other proteins, when histidine or lysine replaced Asp¹¹², the mutant channel was still anion permeable. Evidently, the proton specificity of hH_V1 requires an acidic group at the selectivity filter.

Voltage-gated proton channel in a dinoflagellate

Fogel and Hastings first hypothesized the existence of voltage-gated proton channels in 1972 in bioluminescent dinoflagellates, where they were thought to trigger the flash by activating luciferase. Proton channel genes were subsequently identified in human, mouse, and Ciona intestinalis, but their existence in dinoflagellates remained unconfirmed. We identified a candidate proton channel gene from a Karlodinium veneficum cDNA library based on homology with known proton channel genes. K. veneficum is a predatory, nonbioluminescent dinoflagellate that produces toxins responsible for fish kills worldwide. Patch clamp studies on the heterologously expressed gene confirm that it codes for a genuine voltage-gated proton channel, kH(V)1: it is proton-specific and activated by depolarization, its g(H)-V relationship shifts with changes in external or internal pH, and mutation of the selectivity filter (which we identify as Asp(51)) results in loss of proton-specific conduction. Indirect evidence suggests that kH(V)1 is monomeric, unlike other proton channels. Furthermore, kH(V)1 differs from all known proton channels in activating well negative to the Nernst potential for protons, E(H). This unique voltage dependence makes the dinoflagellate proton channel ideally suited to mediate the proton influx postulated to trigger bioluminescence. In contrast to vertebrate proton channels, whose main function is acid extrusion, we propose that proton channels in dinoflagellates have fundamentally different functions of signaling and excitability.

RaptorX: exploiting structure information for protein alignment by statistical inference

This work presents RaptorX, a statistical method for template-based protein modeling that improves alignment accuracy by exploiting structural information in a single or multiple templates. RaptorX consists of three major components: single-template threading, alignment quality prediction, and multiple-template threading. This work summarizes the methods used by RaptorX and presents its CASP9 result analysis, aiming to identify major bottlenecks with RaptorX and template-based modeling and hopefully directions for further study. Our results show that template structural information helps a lot with both single-template and multiple-template protein threading especially when closely-related templates are unavailable, and there is still large room for improvement in both alignment and template selection. The RaptorX web server is available at http://raptorx.uchicago.edu.

The Cryptosporidium parvum kinome

BACKGROUND

Hundreds of millions of people are infected with cryptosporidiosis annually, with immunocompromised individuals suffering debilitating symptoms and children in socioeconomically challenged regions at risk of repeated infections. There is currently no effective drug available. In order to facilitate the pursuit of anti-cryptosporidiosis targets and compounds, our study spans the classification of the Cryptosporidium parvum kinome and the structural and biochemical characterization of representatives from the CDPK family and a MAP kinase.

RESULTS

The C. parvum kinome comprises over 70 members, some of which may be promising drug targets. These C. parvum protein kinases include members in the AGC, Atypical, CaMK, CK1, CMGC, and TKL groups; however, almost 35% could only be classified as OPK (other protein kinases). In addition, about 25% of the kinases identified did not have any known orthologues outside of Cryptosporidium spp. Comparison of specific kinases with their Plasmodium falciparum and Toxoplasma gondii orthologues revealed some distinct characteristics within the C. parvum kinome, including potential targets and opportunities for drug design. Structural and biochemical analysis of 4 representatives of the CaMK group and a MAP kinase confirms features that may be exploited in inhibitor design. Indeed, screening CpCDPK1 against a library of kinase inhibitors yielded a set of the pyrazolopyrimidine derivatives (PP1-derivatives) with IC₅₀ values of < 10 nM. The binding of a PP1-derivative is further described by an inhibitor-bound crystal structure of CpCDPK1. In addition, structural analysis of CpCDPK4 identified an unprecedented Zn-finger within the CDPK kinase domain that may have implications for its regulation.

CONCLUSIONS

Identification and comparison of the C. parvum protein kinases against other parasitic kinases shows how orthologue- and family-based research can be used to facilitate characterization of promising drug targets and the search for new drugs.

Ligand discovery from a dopamine D3 receptor homology model and crystal structure

G protein-coupled receptors (GPCRs) are intensely studied as drug targets and for their role in signaling. With the determination of the first crystal structures, interest in structure-based ligand discovery increased. Unfortunately, for most GPCRs no experimental structures are available. The determination of the D(3) receptor structure and the challenge to the community to predict it enabled a fully prospective comparison of ligand discovery from a modeled structure versus that of the subsequently released crystal structure. Over 3.3 million molecules were docked against a homology model, and 26 of the highest ranking were tested for binding. Six had affinities ranging from 0.2 to 3.1 μM. Subsequently, the crystal structure was released and the docking screen repeated. Of the 25 compounds selected, five had affinities ranging from 0.3 to 3.0 μM. One of the new ligands from the homology model screen was optimized for affinity to 81 nM. The feasibility of docking screens against modeled GPCRs more generally is considered.

Global transcriptome analysis of constitutive resistance to the white pine weevil in spruce

Constitutive defense mechanisms are critical to the understanding of defense mechanisms in conifers because they constitute the first barrier to attacks by insect pests. In interior spruce, trees that are putatively resistant and susceptible to attacks by white pine weevil (Pissodes strobi) typically exhibit constitutive differences in traits such as resin duct size and number, bark thickness, and terpene content. To improve our knowledge of their genetic basis, we compared globally the constitutive expression levels of 17,825 genes between 20 putatively resistant and 20 putatively susceptible interior spruce trees from the British Columbia tree improvement program. We identified 54 upregulated and 137 downregulated genes in resistant phenotypes, relative to susceptible phenotypes, with a maximum fold change of 2.24 and 3.91, respectively. We found a puzzling increase of resistance by downregulated genes, as one would think that "procuring armaments" is the best defense. Also, although terpenes and phenolic compounds play an important role in conifer defense, we found few of these genes to be differentially expressed. We found 15 putative small heat-shock proteins (sHSP) and several other stress-related proteins to be downregulated in resistant trees. Downregulated putative sHSP belong to several sHSP classes and represented 58% of all tested putative sHSP. These proteins are well known to be involved in plant response to various kinds of abiotic stress; however, their role in constitutive resistance is not yet understood. The lack of correspondence between transcriptome profile clusters and phenotype classifications suggests that weevil resistance in spruce is a complex trait.

Molecular evolution of cyclin proteins in animals and fungi

BACKGROUND

The passage through the cell cycle is controlled by complexes of cyclins, the regulatory units, with cyclin-dependent kinases, the catalytic units. It is also known that cyclins form several families, which differ considerably in primary structure from one eukaryotic organism to another. Despite these lines of evidence, the relationship between the evolution of cyclins and their function is an open issue. Here we present the results of our study on the molecular evolution of A-, B-, D-, E-type cyclin proteins in animals and fungi.

RESULTS

We constructed phylogenetic trees for these proteins, their ancestral sequences and analyzed patterns of amino acid replacements. The analysis of infrequently fixed atypical amino acid replacements in cyclins evidenced that accelerated evolution proceeded predominantly during paralog duplication or after it in animals and fungi and that it was related to aromorphic changes in animals. It was shown also that evolutionary flexibility of cyclin function may be provided by consequential reorganization of regions on protein surface remote from CDK binding sites in animal and fungal cyclins and by functional differentiation of paralogous cyclins formed in animal evolution.

CONCLUSIONS

The results suggested that changes in the number and/or nature of cyclin-binding proteins may underlie the evolutionary role of the alterations in the molecular structure of cyclins and their involvement in diverse molecular-genetic events.

Computational analysis of DNA replicases in double-stranded DNA viruses: relationship with the genome size

Genome duplication in free-living cellular organisms is performed by DNA replicases that always include a DNA polymerase, a DNA sliding clamp and a clamp loader. What are the evolutionary solutions for DNA replicases associated with smaller genomes? Are there some general principles? To address these questions we analyzed DNA replicases of double-stranded (ds) DNA viruses. In the process we discovered highly divergent B-family DNA polymerases in phiKZ-like phages and remote sliding clamp homologs in Ascoviridae family and Ma-LMM01 phage. The analysis revealed a clear dependency between DNA replicase components and the viral genome size. As the genome size increases, viruses universally encode their own DNA polymerases and frequently have homologs of DNA sliding clamps, which sometimes are accompanied by clamp loader subunits. This pattern is highly non-random. The absence of sliding clamps in large viral genomes usually coincides with the presence of atypical polymerases. Meanwhile, sliding clamp homologs, not accompanied by clamp loaders, have an elevated positive electrostatic potential, characteristic of non-ring viral processivity factors that bind the DNA directly. Unexpectedly, we found that similar electrostatic properties are shared by the eukaryotic 9-1-1 clamp subunits, Hus1 and, to a lesser extent, Rad9, also suggesting the possibility of direct DNA binding.

Structural and functional analysis of PTPMT1, a phosphatase required for cardiolipin synthesis

PTPMT1 (PTP localized to the Mitochondrion 1) is a member of the protein tyrosine phosphatase superfamily that is localized exclusively to the mitochondrion. We recently reported that PTPMT1 dephosphorylates phosphatidylglycerol phosphate, an essential intermediate of cardiolipin biosynthesis. To gain further insights into the molecular basis of PTPMT1 function, we determined the crystal structures of the phosphatase domain of PTPMT1. PTPMT1 exhibits a canonical protein tyrosine phosphatase domain fold, resembling many dual-specificity phosphatases such as phosphatase and tensin homolog and vaccinia H1-related phosphatase. We also determined the structure of the catalytically inactive phosphatase in complex with a surrogate substrate, phosphatidylinositol 5-phosphate, which sheds light on the substrate recognition and specificity of PTPMT1. Comparison of the apo and substrate-bound structures of PTPMT1 suggests that it undergoes significant conformational change during catalysis, and we further demonstrated that an evolutionarily conserved EEYE loop is important for its activity.

Principles of activation and permeation in an anion-selective Cys-loop receptor

Fast inhibitory neurotransmission is essential for nervous system function and is mediated by binding of inhibitory neurotransmitters to receptors of the Cys-loop family embedded in the membranes of neurons. Neurotransmitter binding triggers a conformational change in the receptor, opening an intrinsic chloride channel and thereby dampening neuronal excitability. Here we present the first three-dimensional structure, to our knowledge, of an inhibitory anion-selective Cys-loop receptor, the homopentameric Caenorhabditis elegans glutamate-gated chloride channel α (GluCl), at 3.3 Å resolution. The X-ray structure of the GluCl-Fab complex was determined with the allosteric agonist ivermectin and in additional structures with the endogenous neurotransmitter L-glutamate and the open-channel blocker picrotoxin. Ivermectin, used to treat river blindness, binds in the transmembrane domain of the receptor and stabilizes an open-pore conformation. Glutamate binds in the classical agonist site at subunit interfaces, and picrotoxin directly occludes the pore near its cytosolic base. GluCl provides a framework for understanding mechanisms of fast inhibitory neurotransmission and allosteric modulation of Cys-loop receptors.

Microviridae goes temperate: microvirus-related proviruses reside in the genomes of Bacteroidetes

The Microviridae comprises icosahedral lytic viruses with circular single-stranded DNA genomes. The family is divided into two distinct groups based on genome characteristics and virion structure. Viruses infecting enterobacteria belong to the genus Microvirus, whereas those infecting obligate parasitic bacteria, such as Chlamydia, Spiroplasma and Bdellovibrio, are classified into a subfamily, the Gokushovirinae. Recent metagenomic studies suggest that members of the Microviridae might also play an important role in marine environments. In this study we present the identification and characterization of Microviridae-related prophages integrated in the genomes of species of the Bacteroidetes, a phylum not previously known to be associated with microviruses. Searches against metagenomic databases revealed the presence of highly similar sequences in the human gut. This is the first report indicating that viruses of the Microviridae lysogenize their hosts. Absence of associated integrase-coding genes and apparent recombination with dif-like sequences suggests that Bacteroidetes-associated microviruses are likely to rely on the cellular chromosome dimer resolution machinery. Phylogenetic analysis of the putative major capsid proteins places the identified proviruses into a group separate from the previously characterized microviruses and gokushoviruses, suggesting that the genetic diversity and host range of bacteriophages in the family Microviridae is wider than currently appreciated.

T-Coffee: a web server for the multiple sequence alignment of protein and RNA sequences using structural information and homology extension

This article introduces a new interface for T-Coffee, a consistency-based multiple sequence alignment program. This interface provides an easy and intuitive access to the most popular functionality of the package. These include the default T-Coffee mode for protein and nucleic acid sequences, the M-Coffee mode that allows combining the output of any other aligners, and template-based modes of T-Coffee that deliver high accuracy alignments while using structural or homology derived templates. These three available template modes are Expresso for the alignment of protein with a known 3D-Structure, R-Coffee to align RNA sequences with conserved secondary structures and PSI-Coffee to accurately align distantly related sequences using homology extension. The new server benefits from recent improvements of the T-Coffee algorithm and can align up to 150 sequences as long as 10,000 residues and is available from both http://www.tcoffee.org and its main mirror http://tcoffee.crg.cat.

A thaumarchaeal provirus testifies for an ancient association of tailed viruses with archaea

Archaeal viruses, or archaeoviruses, display a wide range of virion morphotypes. Whereas the majority of those morphotypes are unique to archaeal viruses, some are more widely distributed across different cellular domains. Tailed double-stranded DNA archaeoviruses are remarkably similar to viruses of the same morphology (order Caudovirales) that infect many bacterial hosts. They have, so far, only been found in one phylum of the archaea, the Euryarchaeota, which has led to controversial hypotheses about their origin. In the present paper, we describe the identification and analysis of a putative provirus present in the genome of a mesophilic thaumarchaeon. We show that the provirus is related to tailed bacterial and euryarchaeal viruses and encodes a full complement of proteins that are required to build a tailed virion. The recently discovered wide distribution of tailed viruses in Euryarchaeota and the identification of a related provirus in Thaumarchaeota, an archaeal phylum which might have branched off before the separation of Crenarchaeota and Euryarchaeota, suggest that an association of these viruses with Archaea might be more ancient than previously anticipated.

Conjugal plasmid transfer in Streptomyces resembles bacterial chromosome segregation by FtsK/SpoIIIE

Conjugation is a major route of horizontal gene transfer, the driving force in the evolution of bacterial genomes. Antibiotic producing soil bacteria of the genus Streptomyces transfer DNA in a unique process involving a single plasmid-encoded protein TraB and a double-stranded DNA molecule. However, the molecular function of TraB in directing DNA transfer from a donor into a recipient cell is unknown. Here, we show that TraB constitutes a novel conjugation system that is clearly distinguished from DNA transfer by a type IV secretion system. We demonstrate that TraB specifically recognizes and binds to repeated 8 bp motifs on the conjugative plasmid. The specific DNA recognition is mediated by helix α3 of the C-terminal winged-helix-turn-helix domain of TraB. We show that TraB assembles to a hexameric ring structure with a central ∼3.1 nm channel and forms pores in lipid bilayers. Structure, sequence similarity and DNA binding characteristics of TraB indicate that TraB is derived from an FtsK-like ancestor protein, suggesting that Streptomyces adapted the FtsK/SpoIIIE chromosome segregation system to transfer DNA between two distinct Streptomyces cells.

Revisiting gap locations in amino acid sequence alignments and a proposal for a method to improve them by introducing solvent accessibility

In comparative modeling, the quality of amino acid sequence alignment still constitutes a major bottleneck in the generation of high quality models of protein three-dimensional (3D) structures. Substantial efforts have been made to improve alignment quality by revising the substitution matrix, introducing multiple sequences, replacing dynamic programming with hidden Markov models, and incorporating 3D structure information. Improvements in the gap penalty have not been a major focus, however, following the development of the affine gap penalty and of the secondary structure dependent gap penalty. We revisited the correlation between protein 3D structure and gap location in a large protein 3D structure data set, and found that the frequency of gap locations approximated to an exponential function of the solvent accessibility of the inserted residues. The nonlinearity of the gap frequency as a function of accessibility corresponded well to the relationship between residue mutation pattern and residue accessibility. By introducing this relationship into the gap penalty calculation for pairwise alignment between template and target amino acid sequences, we were able to obtain a sequence alignment much closer to the structural alignment. The quality of the alignments was substantially improved on a pair of sequences with identity in the “twilight zone” between 20 and 40%. The relocation of gaps by our new method made a significant improvement in comparative modeling, exemplified here by the Bacillus subtilis yitF protein. The method was implemented in a computer program, ALAdeGAP (ALignment with Accessibility dependent GAp Penalty), which is available at http://cib.cf.ocha.ac.jp/target_protein/. Proteins 2011; © 2011 Wiley-Liss, Inc.

A multiple-template approach to protein threading

Most threading methods predict the structure of a protein using only a single template. Due to the increasing number of solved structures, a protein without solved structure is very likely to have more than one similar template structures. Therefore, a natural question to ask is if we can improve modeling accuracy using multiple templates. This article describes a new multiple-template threading method to answer this question. At the heart of this multiple-template threading method is a novel probabilistic-consistency algorithm that can accurately align a single protein sequence simultaneously to multiple templates. Experimental results indicate that our multiple-template method can improve pairwise sequence-template alignment accuracy and generate models with better quality than single-template models even if they are built from the best single templates (P-value <10(-6)) while many popular multiple sequence/structure alignment tools fail to do so. The underlying reason is that our probabilistic-consistency algorithm can generate accurate multiple sequence/template alignments. In another word, without an accurate multiple sequence/template alignment, the modeling accuracy cannot be improved by simply using multiple templates to increase alignment coverage. Blindly tested on the CASP9 targets with more than one good template structures, our method outperforms all other CASP9 servers except two (Zhang-Server and QUARK of the same group). Our probabilistic-consistency algorithm can possibly be extended to align multiple protein/RNA sequences and structures.

Study of protein complexes via homology modeling, applied to cysteine proteases and their protein inhibitors

This paper develops and evaluates large-scale calculation of 3D structures of protein complexes by homology modeling as a promising new approach for protein docking. The complexes investigated were papain-like cysteine proteases and their protein inhibitors, which play numerous roles in human and parasitic metabolisms. The structural modeling was performed in two parts. For the first part (evaluation set), nine crystal structure complexes were selected, 1325 homology models of known complexes were rebuilt by various templates including hybrids, allowing an analysis of the factors influencing the accuracy of the models. The important considerations for modeling the interface were protease coverage and inhibitor sequence identity. In the second part (study set), the findings of the evaluation set were used to select appropriate templates to model novel cysteine protease-inhibitor complexes from human and malaria parasites Plasmodium falciparum and Plasmodium vivax. The energy scores, considering the evaluation set, indicate that the models are of high accuracy.