Alignment of distantly related protein structures: algorithm, bound and implications to homology modeling

Comprehensive analysis of Novel mutations in CCM1/KRIT1 and CCM2/MGC4607 and their clinical implications in Cerebral Cavernous malformations

BACKGROUND

Cerebral Cavernous Malformations (CCM) is a genetic disease characterized by vascular abnormalities in the brain and spinal cord, affecting 0.4-0.5 % of the population. We identified two novel pathogenic mutations, CCM1/KRIT1 c.811delT (p.Trp271GlyfsTer5) and CCM2/MGC4607 c.613_614insGG p.Glu205GlyfsTer31), which disrupt crucial protein domains and potentially alter disease progression.

OBJECTIVE

The study aims to comprehensively analyze a Brazilian cohort of CCM patients, integrating genetic, clinical, and structural aspects. Specifically, we sought to identify novel mutations within the CCM complex, and explore their potential impact on disease progression.

METHODS

We conducted a detailed examination of neuroradiological and clinical features in both symptomatic and asymptomatic CCM patients, performing genetic analyses through sequencing of the CCM1/KRIT1, CCM2/MGC4607, and CCM3/PDCD10 genes In silico structural predictions were carried out using PolyPhen-2, SIFT, and Human Genomics Community tools. Protein-protein interactions and docking analyses were explored using the STRING database.

RESULTS

Genetic analysis identifies 6 pathogenic mutations, 4 likely pathogenic, 1 variants of uncertain significance, and 7 unclassified mutations, including the novel mutations in CCM1 c.811delT and CCM2 c.613_614insGG. In silico structural analysis revealed significant alterations in protein structure, supporting their pathogenicity. Protein-protein interaction analysis indicated nuanced impacts on cellular processes. Clinically, we observed a broad spectrum of symptoms, including seizures and focal neurological deficits. However, no statistically significant differences were found in lesion burden, age of first symptom onset, or sex between the identified CCM1/KRIT1 and CCM2/MGC4607 mutations among all patients studied.

CONCLUSION

This study enhances the understanding of CCM by linking clinical variability, genetic mutations, and structural effects. The identification of these novel mutations opens new avenues for research and potential therapeutic strategies.

AARS Online: A collaborative database on the structure, function, and evolution of the aminoacyl-tRNA synthetases

The aminoacyl-tRNA synthetases (aaRS) are a large group of enzymes that implement the genetic code in all known biological systems. They attach amino acids to their cognate tRNAs, moonlight in various translational and non-translational activities beyond aminoacylation, and are linked to many genetic disorders. The aaRS have a subtle ontology characterized by structural and functional idiosyncrasies that vary from organism to organism, and protein to protein. Across the tree of life, the 22 coded amino acids are handled by 16 evolutionary families of Class I aaRS and 21 families of Class II aaRS. We introduce AARS Online, an interactive Wikipedia-like tool curated by an international consortium of field experts. This platform systematizes existing knowledge about the aaRS by showcasing a taxonomically diverse selection of aaRS sequences and structures. Through its graphical user interface, AARS Online facilitates a seamless exploration between protein sequence and structure, providing a friendly introduction to the material for non-experts and a useful resource for experts. Curated multiple sequence alignments can be extracted for downstream analyses. Accessible at www.aars.online, AARS Online is a free resource to delve into the world of the aaRS.

Meta-omics analysis reveals the marine arsenic cycle driven by bacteria

Arsenic is a toxic element widely distributed in the Earth's crust and ranked as a class I human carcinogen. Microbial metabolism makes significant contributions to arsenic detoxification, migration and transformation. Nowadays, research on arsenic is primarily in areas affected by arsenic pollution associated with human health activities. However, the biogeochemical traits of arsenic in the global marine ecosystem remain to be explicated. In this study, we revealed that seawater environments were primarily governed by the process of arsenate reduction to arsenite, while arsenite methylation was predominant in marine sediments which may serve as significant sources of arsenic emission into the atmosphere. Significant disparities existed in the distribution patterns of the arsenic cycle between surface and deep seawaters at middle and low latitudes, whereas these situations tend to be similar in the Arctic and Antarctic oceans. Significant variations were also observed in the taxonomic diversity and core microbial community of arsenic cycling across different marine environments. Specifically, γ-proteobacteria played a pivotal role in the arsenic cycle in the whole marine environment. Temperature, dissolved oxygen and phosphate were the crucial factors that related to these differentiations in seawater environments. Overall, our study contributes to a deeper understanding of the marine arsenic cycle.

High-Throughput Screening Method Using Escherichia coli Keio Mutants for Assessing Primary Damage Mechanism of Antimicrobials

The Escherichia coli Keio mutant collection has been a tool for assessing the role of specific genes and determining their role in E. coli physiology and uncovering novel functions. In this work, specific mutants in the DNA repair pathways and oxidative stress response were evaluated to identify the primary targets of silver nanoparticles (NPs) and their mechanism of action. The results presented in this work suggest that NPs mainly target DNA via double-strand breaks and base modifications since the recA, uvrC, mutL, and nfo mutants rendered the most susceptible phenotype, rather than involving the oxidative stress response. Concomitantly, during the establishment of the control conditions for each mutant, the katG and sodA mutants showed a hypersensitive phenotype to mitomycin C, an alkylating agent. Thus, we propose that KatG catalase plays a key role as a cellular chaperone, as reported previously for the filamentous fungus Neurospora crassa, a large subunit catalase. The Keio collection mutants may also be a key tool for assessing the resistance mechanism to metallic NPs by using their potential to identify novel pathways involved in the resistance to NPs.

Comprehensive CCM3 Mutational Analysis in Two Patients with Syndromic Cerebral Cavernous Malformation

Cerebral cavernous malformation (CCM) is a vascular disease that affects the central nervous system, which familial form is due to autosomal dominant mutations in the genes KRIT1(CCM1), MGC4607(CCM2), and PDCD10(CCM3). Patients affected by the PDCD10 mutations usually have the onset of symptoms at an early age and a more aggressive phenotype. The aim of this study is to investigate the molecular mechanism involved with CCM3 disease pathogenesis. Herein, we report two typical cases of CCM3 phenotype and compare the clinical and neuroradiological findings with five patients with a familial form of KRIT1 or CCM2 mutations and six patients with a sporadic form. In addition, we evaluated the PDCD10 gene expression by qPCR and developed a bioinformatic pipeline to understand the structural changes of mutations. The two CCM3 patients had an early onset of symptoms and a high lesion burden. Furthermore, the sequencing showed that Patient 1 had a frameshift mutation in c.222delT; p.(Asn75Thrfs*14) that leads to lacking the last 124 C-terminal amino acids on its primary structure and Patient 2 had a variant on the splicing site region c.475-2A > G. The mRNA expression was fourfold lower in both patients with PDCD10 mutation. Using in silico analysis, we identify that the frameshift mutation transcript lacks the C-terminal FAT-homology domain compared to the wild-type PDCD10 and preserves the N-terminal dimerization domain. The two patients studied here allow estimating the potential impact of mutations in clinical interpretation as well as support to better understand the mechanism and pathogenesis of CCM3.

Enzymic recognition of amino acids drove the evolution of primordial genetic codes

How genetic information gained its exquisite control over chemical processes needed to build living cells remains an enigma. Today, the aminoacyl-tRNA synthetases (AARS) execute the genetic codes in all living systems. But how did the AARS that emerged over three billion years ago as low-specificity, protozymic forms then spawn the full range of highly-specific enzymes that distinguish between 22 diverse amino acids? A phylogenetic reconstruction of extant AARS genes, enhanced by analysing modular acquisitions, reveals six AARS with distinct bacterial, archaeal, eukaryotic, or organellar clades, resulting in a total of 36 families of AARS catalytic domains. Small structural modules that differentiate one AARS family from another played pivotal roles in discriminating between amino acid side chains, thereby expanding the genetic code and refining its precision. The resulting model shows a tendency for less elaborate enzymes, with simpler catalytic domains, to activate amino acids that were not synthesised until later in the evolution of the code. The most probable evolutionary route for an emergent amino acid type to establish a place in the code was by recruiting older, less specific AARS, rather than adapting contemporary lineages. This process, retrofunctionalisation, differs from previously described mechanisms through which amino acids would enter the code.

Genomic database furnishes a spontaneous example of a functional Class II glycyl-tRNA synthetase urzyme

The chief barrier to studies of how genetic coding emerged is the lack of experimental models for ancestral aminoacyl-tRNA synthetases (AARS). We hypothesized that conserved core catalytic sites could represent such ancestors. That hypothesis enabled engineering functional "urzymes" from TrpRS, LeuRS, and HisRS. We describe here a fourth urzyme, GlyCA, detected in an open reading frame from the genomic record of the arctic fox, Vulpes lagopus. GlyCA is homologous to a bacterial heterotetrameric Class II GlyRS-B. Alphafold2 predicted that the N-terminal 81 amino acids would adopt a 3D structure nearly identical to the HisRS urzyme (HisCA1). We expressed and purified that N-terminal segment. Enzymatic characterization revealed a robust single-turnover burst size and a catalytic rate for ATP consumption well in excess of that previously published for HisCA1. Time-dependent aminoacylation of tRNAGly proceeds at a rate consistent with that observed for amino acid activation. In fact, GlyCA is actually 35 times more active in glycine activation by ATP than the full-length GlyRS-B α-subunit dimer. ATP-dependent activation of the 20 canonical amino acids favors Class II amino acids that complement those favored by HisCA and LeuAC. These properties reinforce the notion that urzymes represent the requisite ancestral catalytic activities to implement a reduced genetic coding alphabet.

Dual function of the O-antigen WaaL ligase of Aggregatibacter actinomycetemcomitans

Protein glycosylation is critical to the quaternary structure and collagen-binding activity of the extracellular matrix protein adhesin A (EmaA) associated with Aggregatibacter actinomycetemcomitans. The glycosylation of this large, trimeric autotransporter adhesin is postulated to be mediated by WaaL, an enzyme with the canonical function to ligate the O-polysaccharide (O-PS) antigen with a terminal sugar of the lipid A-core oligosaccharide of lipopolysaccharide (LPS). In this study, we have determined that the Escherichia coli waaL ortholog (rflA) does not restore collagen binding of a waaL mutant strain of A. actinomycetemcomitans but does restore O-PS ligase activity following transformation of a plasmid expressing waaL. Therefore, a heterologous E. coli expression system was developed constituted of two independently replicating plasmids expressing either waaL or emaA of A. actinomycetemcomitans to directly demonstrate the necessity of ligase activity for EmaA collagen binding. Proper expression of the protein encoded by each plasmid was characterized, and the individually transformed strains did not promote collagen binding. However, coexpression of the two plasmids resulted in a strain with a significant increase in collagen binding activity and a change in the biochemical properties of the protein. These results provide additional data supporting the novel hypothesis that the WaaL ligase of A. actinomycetemcomitans shares a dual role as a ligase in LPS biosynthesis and is required for collagen binding activity of EmaA.

Structural similarities between the metacyclic and bloodstream form variant surface glycoproteins of the African trypanosome

During infection of mammalian hosts, African trypanosomes thwart immunity using antigenic variation of the dense Variant Surface Glycoprotein (VSG) coat, accessing a large repertoire of several thousand genes and pseudogenes, and switching to antigenically distinct copies. The parasite is transferred to mammalian hosts by the tsetse fly. In the salivary glands of the fly, the pathogen adopts the metacyclic form and expresses a limited repertoire of VSG genes specific to that developmental stage. It has remained unknown whether the metacyclic VSGs possess distinct properties associated with this particular and discrete phase of the parasite life cycle. We present here three novel metacyclic form VSG N-terminal domain crystal structures (mVSG397, mVSG531, and mVSG1954) and show that they mirror closely in architecture, oligomerization, and surface diversity the known classes of bloodstream form VSGs. These data suggest that the mVSGs are unlikely to be a specialized subclass of VSG proteins, and thus could be poor candidates as the major components of prophylactic vaccines against trypanosomiasis.

On the stability and layered organization of protein-DNA condensates

Multi-component phase separation is emerging as a key mechanism for the formation of biological condensates that play essential roles in signal sensing and transcriptional regulation. The molecular factors that dictate these condensates' stability and spatial organization are not fully understood, and it remains challenging to predict their microstructures. Using a near-atomistic, chemically accurate force field, we studied the phase behavior of chromatin regulators that are crucial for heterochromatin organization and their interactions with DNA. Our computed phase diagrams recapitulated previous experimental findings on different proteins. They revealed a strong dependence of condensate stability on the protein-DNA mixing ratio as a result of balancing protein-protein interactions and charge neutralization. Notably, a layered organization was observed in condensates formed by mixing HP1, histone H1, and DNA. This layered organization may be of biological relevance, as it enables cooperative DNA packaging between the two chromatin regulators: histone H1 softens the DNA to facilitate the compaction induced by HP1 droplets. Our study supports near-atomistic models as a valuable tool for characterizing the structure and stability of biological condensates.

The Structure and Function of Modular Escherichia coli O157:H7 Bacteriophage FTBEc1 endolysin, LysT84: Defining a New Endolysin Catalytic Subfamily

Bacteriophage endolysins degrade peptidoglycan and have been identified as antibacterial candidates to combat antimicrobial resistance. Considering the catalytic and structural diversity of endolysins, there is a paucity of structural data to inform how these enzymes work at the molecular level-key data that is needed to realize the potential of endolysin-based antibacterial agents. Here, we determine the atomic structure and define the enzymatic function of Escherichia coli O157:H7 phage FTEBc1 endolysin, LysT84. Bioinformatic analysis reveals that LysT84 is a modular endolysin, which is unusual for Gram-negative endolysins, comprising a peptidoglycan binding domain and an enzymatic domain. The crystal structure of LysT84 (2.99 Å) revealed a mostly α-helical protein with two domains connected by a linker region but packed together. LysT84 was determined to be a monomer in solution using analytical ultracentrifugation. Small-angle X-ray scattering data revealed that LysT84 is a flexible protein but does not have the expected bimodal P(r) function of a multidomain protein, suggesting that the domains of LysT84 pack closely creating a globular protein as seen in the crystal structure. Structural analysis reveals two key glutamate residues positioned on either side of the active site cavity; mutagenesis demonstrating these residues are critical for peptidoglycan degradation. Molecular dynamic simulations suggest that the enzymatically active domain is dynamic, allowing the appropriate positioning of these catalytic residues for hydrolysis of the β(1-4) bond. Overall, our study defines the structural basis for peptidoglycan degradation by LysT84 which supports rational engineering of related endolysins into effective antibacterial agents.

Catalytic Fields as a Tool to Analyze Enzyme Reaction Mechanism Variants and Reaction Steps

Catalytic fields representing the topology of the optimal molecular environment charge distribution that reduces the activation barrier have been used to examine alternative reaction variants and to determine the role of conserved catalytic residues for two consecutive reactions catalyzed by the same enzyme. Until now, most experimental and conventional top-down theoretical studies employing QM/MM or ONIOM methods have focused on the role of enzyme electric fields acting on broken bonds of reactants. In contrast, our bottom-up approach dealing with a small reactant and transition-state model allows the analysis of the opposite effects: how the catalytic field resulting from the charge redistribution during the enzyme reaction acts on conserved amino acid residues and contributes to the reduction of the activation barrier. This approach has been applied to the family of histidyl tRNA synthetases involved in the translation of the genetic code into the protein amino acid sequence. Activation energy changes related to conserved charged amino acid residues for 12 histidyl tRNA synthetases from different biological species allowed to compare on equal footing the catalytic residues involved in ATP aminoacylation and tRNA charging reactions and to analyze different reaction mechanisms proposed in the literature. A scan of the library of atomic multipoles for amino acid side-chain rotamers within the catalytic field pointed out the change in the Glu83 conformation as the critical catalytic effect, providing, at low computational cost, insight into the electrostatic preorganization of the enzyme catalytic site at a level of detail that has not yet been accessible in conventional experimental or theoretical methods. This opens the way for rational reverse biocatalyst design at a very limited computational cost without resorting to empirical methods.

DAMA: a method for computing multiple alignments of protein structures using local structure descriptors

MOTIVATION

The well-known fact that protein structures are more conserved than their sequences forms the basis of several areas of computational structural biology. Methods based on the structure analysis provide more complete information on residue conservation in evolutionary processes. This is crucial for the determination of evolutionary relationships between proteins and for the identification of recurrent structural patterns present in biomolecules involved in similar functions. However, algorithmic structural alignment is much more difficult than multiple sequence alignment. This study is devoted to the development and applications of DAMA-a novel effective environment capable to compute and analyze multiple structure alignments.

RESULTS

DAMA is based on local structural similarities, using local 3D structure descriptors and thus accounts for nearest-neighbor molecular environments of aligned residues. It is constrained neither by protein topology nor by its global structure. DAMA is an extension of our previous study (DEDAL) which demonstrated the applicability of local descriptors to pairwise alignment problems. Since the multiple alignment problem is NP-complete, an effective heuristic approach has been developed without imposing any artificial constraints. The alignment algorithm searches for the largest, consistent ensemble of similar descriptors. The new method is capable to capture most of the biologically significant similarities present in canonical test sets and is discriminatory enough to prevent the emergence of larger, but meaningless, solutions. Tests performed on the test sets, including protein kinases, demonstrate DAMA's capability of identifying equivalent residues, which should be very useful in discovering the biological nature of proteins similarity. Performance profiles show the advantage of DAMA over other methods, in particular when using a strict similarity measure QC, which is the ratio of correctly aligned columns, and when applying the methods to more difficult cases.

AVAILABILITY AND IMPLEMENTATION

DAMA is available online at http://dworkowa.imdik.pan.pl/EP/DAMA. Linux binaries of the software are available upon request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Molecular Evolutionary Analysis of Nematode Zona Pellucida (ZP) Modules Reveals Disulfide-Bond Reshuffling and Standalone ZP-C Domains

Zona pellucida (ZP) modules mediate extracellular protein-protein interactions and contribute to important biological processes including syngamy and cellular morphogenesis. Although some biomedically relevant ZP modules are well studied, little is known about the protein family's broad-scale diversity and evolution. The increasing availability of sequenced genomes from "nonmodel" systems provides a valuable opportunity to address this issue and to use comparative approaches to gain new insights into ZP module biology. Here, through phylogenetic and structural exploration of ZP module diversity across the nematode phylum, I report evidence that speaks to two important aspects of ZP module biology. First, I show that ZP-C domains-which in some modules act as regulators of ZP-N domain-mediated polymerization activity, and which have never before been found in isolation-can indeed be found as standalone domains. These standalone ZP-C domain proteins originated in independent (paralogous) lineages prior to the diversification of extant nematodes, after which they evolved under strong stabilizing selection, suggesting the presence of ZP-N domain-independent functionality. Second, I provide a much-needed phylogenetic perspective on disulfide bond variability, uncovering evidence for both convergent evolution and disulfide-bond reshuffling. This result has implications for our evolutionary understanding and classification of ZP module structural diversity and highlights the usefulness of phylogenetics and diverse sampling for protein structural biology. All told, these findings set the stage for broad-scale (cross-phyla) evolutionary analysis of ZP modules and position Caenorhabditis elegans and other nematodes as important experimental systems for exploring the evolution of ZP modules and their constituent domains.

LJ, Franco Castrillón J, Rojas-Montoya M, Castaño Osorio JC. New Insect Host Defense Peptides (HDP) From Dung Beetle (Coleoptera: Scarabaeidae) Transcriptomes

The Coleoptera Scarabaeidae family is one of the most diverse groups of insects on the planet, which live in complex microbiological environments. Their immune systems have evolved diverse families of Host Defense Peptides (HDP) with strong antimicrobial and immunomodulatory activities. However, there are several peptide sequences that await discovery in this group of organisms. This would pave the way to identify molecules with promising therapeutic potential. This work retrieved two sources of information: 1) De-novo transcriptomic data from two species of neotropical Scarabaeidae (Dichotomius satanas and Ontophagus curvicornis); 2) Sequence data deposited in available databases. A Blast-based search was conducted against the transcriptomes with a subset of sequences representative of the HDP. This work reports 155 novel HDP sequences identified in nine transcriptomes from seven species of Coleoptera: D. satanas (n = 76; 49.03%), O. curvicornis (n = 23; 14.83%), (Trypoxylus dichotomus) (n = 18; 11.61%), (Onthophagus nigriventris) (n = 10; 6.45%), (Heterochelus sp) (n = 6; 3.87%), (Oxysternon conspicillatum) (n = 18; 11.61%), and (Popillia japonica) (n = 4; 2.58%). These sequences were identified based on similarity to known HDP insect families. New members of defensins (n = 58; 37.42%), cecropins (n = 18; 11.61%), attancins (n = 41; 26.45%), and coleoptericins (n = 38; 24.52%) were described based on their physicochemical and structural characteristics, as well as their sequence relationship to other insect HDPs. Therefore, the Scarabaeidae family is a complex and rich group of insects with a great diversity of antimicrobial peptides with potential antimicrobial activity.

The Molecular Basis for Escherichia coli O157:H7 Phage FAHEc1 Endolysin Function and Protein Engineering to Increase Thermal Stability

Bacteriophage-encoded endolysins have been identified as antibacterial candidates. However, the development of endolysins as mainstream antibacterial agents first requires a comprehensive biochemical understanding. This study defines the atomic structure and enzymatic function of Escherichia coli O157:H7 phage FAHEc1 endolysin, LysF1. Bioinformatic analysis suggests this endolysin belongs to the T4 Lysozyme (T4L)-like family of proteins and contains a highly conserved catalytic triad. We then solved the structure of LysF1 with x-ray crystallography to 1.71 Å. LysF1 was confirmed to exist as a monomer in solution by sedimentation velocity experiments. The protein architecture of LysF1 is conserved between T4L and related endolysins. Comparative analysis with related endolysins shows that the spatial orientation of the catalytic triad is conserved, suggesting the catalytic mechanism of peptidoglycan degradation is the same as that of T4L. Differences in the sequence illustrate the role coevolution may have in the evolution of this fold. We also demonstrate that by mutating a single residue within the hydrophobic core, the thermal stability of LysF1 can be increased by 9.4 °C without compromising enzymatic activity. Overall, the characterization of LysF1 provides further insight into the T4L-like class of endolysins. Our study will help advance the development of related endolysins as antibacterial agents, as rational engineering will rely on understanding mutable positions within this protein fold.

DZIP1 regulates mammalian cardiac valve development through a Cby1-β-catenin mechanism

Background

Mitral valve prolapse (MVP) is a common and progressive cardiovascular disease with developmental origins. How developmental errors contribute to disease pathogenesis are not well understood.

Results

A multimeric complex was identified that consists of the MVP gene Dzip1, Cby1, and β‐catenin. Co‐expression during valve development revealed overlap at the basal body of the primary cilia. Biochemical studies revealed a DZIP1 peptide required for stabilization of the complex and suppression of β‐catenin activities. Decoy peptides generated against this interaction motif altered nuclear vs cytosolic levels of β‐catenin with effects on transcriptional activity. A mutation within this domain was identified in a family with inherited non‐syndromic MVP. This novel mutation and our previously identified DZIP1^S24R variant resulted in reduced DZIP1 and CBY1 stability and increased β‐catenin activities. The β‐catenin target gene, MMP2 was up‐regulated in the Dzip1^S14R/+ valves and correlated with loss of collagenous ECM matrix and myxomatous phenotype.

Conclusion

Dzip1 functions to restrain β‐catenin signaling through a CBY1 linker during cardiac development. Loss of these interactions results in increased nuclear β‐catenin/Lef1 and excess MMP2 production, which correlates with developmental and postnatal changes in ECM and generation of a myxomatous phenotype.

pmTM-align: scalable pairwise and multiple structure alignment with Apache Spark and OpenMP

BACKGROUND

Structure comparison can provide useful information to identify functional and evolutionary relationship between proteins. With the dramatic increase of protein structure data in the Protein Data Bank, computation time quickly becomes the bottleneck for large scale structure comparisons. To more efficiently deal with informative multiple structure alignment tasks, we propose pmTM-align, a parallel protein structure alignment approach based on mTM-align/TM-align. pmTM-align contains two stages to handle pairwise structure alignments with Spark and the phylogenetic tree-based multiple structure alignment task on a single computer with OpenMP.

RESULTS

Experiments with the SABmark dataset showed that parallelization along with data structure optimization provided considerable speedup for mTM-align. The Spark-based structure alignments achieved near ideal scalability with large datasets, and the OpenMP-based construction of the phylogenetic tree accelerated the incremental alignment of multiple structures and metrics computation by a factor of about 2-5.

CONCLUSIONS

pmTM-align enables scalable pairwise and multiple structure alignment computing and offers more timely responses for medium to large-sized input data than existing alignment tools such as mTM-align.

Topology-independent and global protein structure alignment through an FFT-based algorithm

MOTIVATION

Protein structure alignment is one of the fundamental problems in computational structure biology. A variety of algorithms have been developed to address this important issue in the past decade. However, due to their heuristic nature, current structure alignment methods may suffer from suboptimal alignment and/or over-fragmentation and thus lead to a biologically wrong alignment in some cases. To overcome these limitations, we have developed an accurate topology-independent and global structure alignment method through an FFT-based exhaustive search algorithm, which is referred to as FTAlign.

RESULTS

Our FTAlign algorithm was extensively tested on six commonly used datasets and compared with seven state-of-the-art structure alignment approaches, TMalign, DeepAlign, Kpax, 3DCOMB, MICAN, SPalignNS and CLICK. It was shown that FTAlign outperformed the other methods in reproducing manually curated alignments and obtained a high success rate of 96.7 and 90.0% on two gold-standard benchmarks, MALIDUP and MALISAM, respectively. Moreover, FTAlign also achieved the overall best performance in terms of biologically meaningful structure overlap (SO) and TMscore on both the sequential alignment test sets including MALIDUP, MALISAM and 64 difficult cases from HOMSTRAD, and the non-sequential sets including MALIDUP-NS, MALISAM-NS, 199 topology-different cases, where FTAlign especially showed more advantage for non-sequential alignment. Despite its global search feature, FTAlign is also computationally efficient and can normally complete a pairwise alignment within one second.

AVAILABILITY AND IMPLEMENTATION

http://huanglab.phys.hust.edu.cn/ftalign/.

Benchmarking Methods of Protein Structure Alignment

The function of a protein is primarily determined by its structure and amino acid sequence. Many biological questions of interest rely on being able to accurately determine the group of structures to which domains of a protein belong; this can be done through alignment and comparison of protein structures. Dozens of different methods for Protein Structure Alignment (PSA) have been proposed that use a wide range of techniques. The aim of this study is to determine the ability of PSA methods to identify pairs of protein domains known to share differing levels of structural similarity, and to assess their utility for clustering domains from several different folds into known groups. We present the results of a comprehensive investigation into eighteen PSA methods, to our knowledge the largest piece of independent research on this topic. Overall, SP-AlignNS (non-sequential) was found to be the best method for classification, and among the best performing methods for clustering. Methods (where possible) were split into the algorithm used to find the optimal alignment and the score used to assess similarity. This allowed us to largely separate the algorithm from the score it maximizes and thus, to assess their effectiveness independently of each other. Surprisingly, we found that some hybrids of mismatched scores and algorithms performed better than either of the native methods at classification and, in some cases, clustering as well. It is hoped that this investigation and the accompanying discussion will be useful for researchers selecting or designing methods to align protein structures.

Evidence of nuclear transport mechanisms in the protozoan parasite Giardia lamblia

Nuclear-cytoplasmic trafficking of proteins is a highly regulated process that modulates multiple biological processes in eukaryotic cells. In Giardia lamblia, shuttling has been described from the cytoplasm to nuclei of proteins during the biological cell cycle of the parasite. This suggests that a mechanism of nucleocytoplasmic transport is present and functional in G. lamblia. By means of computational biology analyses, we found that there are only two genes for nuclear transport in this parasite, named Importin α and Importin β. When these transporters were overexpressed, both localized close to the nuclear envelope, and no change was observed in trophozoite growth rate. However, during the encystation process, both transporters induced an increase in the number of cysts produced. Importazole and Ivermectin, two known specific inhibitors of importins, separately influenced the encysting process by inducing an arrest in the trophozoite stage that prevents the production of cysts. This effect was more noticeable when Ivermectin, an anti-parasitic drug, was used. Finally, we tested whether the enzyme arginine deiminase, which shuttles from the cytoplasm to the nuclei during encystation, was influenced by these transporters. We found that treatment with each of the inhibitors abrogates arginine deiminase nuclear translocation and favors perinuclear localization. This suggests that Importin α and Importin β are key transporters during the encystation process and are involved, at least, in the transport of arginine deiminase into the nuclei. Considering the effect produced by Ivermectin during growth and encystation, we postulate that this drug could be used to treat giardiasis.

mTM-align: an algorithm for fast and accurate multiple protein structure alignment

Motivation

As protein structure is more conserved than sequence during evolution, multiple structure alignment can be more informative than multiple sequence alignment, especially for distantly related proteins. With the rapid increase of the number of protein structures in the Protein Data Bank, it becomes urgent to develop efficient algorithms for multiple structure alignment.

Results

A new multiple structure alignment algorithm (mTM-align) was proposed, which is an extension of the highly efficient pairwise structure alignment program TM-align. The algorithm was benchmarked on four widely used datasets, HOMSTRAD, SABmark_sup, SABmark_twi and SISY-multiple, showing that mTM-align consistently outperforms other algorithms. In addition, the comparison with the manually curated alignments in the HOMSTRAD database shows that the automated alignments built by mTM-align are in general more accurate. Therefore, mTM-align may be used as a reliable complement to construct multiple structure alignments for real-world applications.

Availability and implementation

http://yanglab.nankai.edu.cn/mTM-align.

Contact

zhng@umich.edu or yangjy@nankai.edu.cn.

Supplementary information

Supplementary data are available at Bioinformatics online.

Tertiary Structural Motif Sequence Statistics Enable Facile Prediction and Design of Peptides that Bind Anti-apoptotic Bfl-1 and Mcl-1

Understanding the relationship between protein sequence and structure well enough to design new proteins with desired functions is a longstanding goal in protein science. Here, we show that recurring tertiary structural motifs (TERMs) in the PDB provide rich information for protein-peptide interaction prediction and design. TERM statistics can be used to predict peptide binding energies for Bcl-2 family proteins as accurately as widely used structure-based tools. Furthermore, design using TERM energies (dTERMen) rapidly and reliably generates high-affinity peptide binders of anti-apoptotic proteins Bfl-1 and Mcl-1 with just 15%-38% sequence identity to any known native Bcl-2 family protein ligand. High-resolution structures of four designed peptides bound to their targets provide opportunities to analyze the strengths and limitations of the computational design method. Our results support dTERMen as a powerful approach that can complement existing tools for protein engineering.

Mechanism of action of the moonlighting protein EfTu as a Substance P sensor in Bacillus cereus

The striking feature of the ubiquitous protein EfTu (Thermo unstable ribosomal Elongation factor) is its moonlighting (multifunctional) activity. Beyond its function at the ribosomal level it should be exported to the bacterial surface and act as an environmental sensor. In Bacillus cereus, and other cutaneous bacteria, it serves as a Substance P (SP) receptor and is essential for bacterial adaptation to the host. However, the modus operandi of EfTu as a bacterial sensor remains to be investigated. Studies realized by confocal and transmission electron microscopy revealed that, in the absence of an exogenous signal, EfTu is not exposed on the bacterial surface but is recruited under the effect of SP. In addition, SP acts as a transcriptional regulator of the tuf gene encoding for EfTu. As observed using gadolinium chloride, an inhibitor of membrane mechanosensitive channels (Msc), Msc control EfTu export and subsequently the bacterial response to SP both in terms of cytotoxicity and biofilm formation activity. Microscale thermophoresis revealed that in response to SP, EfTu can form homopolymers. This event should occur after EfTu export and, as shown by proteo-liposome reconstruction studies, SP appears to promote EfTu polymers association to the membrane, leading subsequently to the bacterial response. Molecular modeling suggests that this mechanism should involve EfTu unfolding and insertion into the bacterial cytoplasmic membrane, presumably through formation of homopolymers. This study is unraveling the original mechanism action of EfTu as a bacterial sensor but also reveals that this protein should have a broader role, including in eukaryotes.

HipA-Mediated Phosphorylation of SeqA Does not Affect Replication Initiation in Escherichia coli

The SeqA protein of Escherichia coli is required to prevent immediate re-initiation of chromosome replication from oriC. The SeqA protein is phosphorylated at the serine-36 (Ser36) residue by the HipA kinase. The role of phosphorylation was addressed by mutating the Ser36 residue to alanine, which cannot be phosphorylated and to aspartic acid, which mimics a phosphorylated serine residue. Both mutant strains were similar to wild-type with respect to origin concentration and initiation synchrony. The minimal time between successive initiations was also unchanged. We therefore suggest that SeqA phosphorylation at the Ser36 residue is silent, at least with respect to SeqA's role in replication initiation.

Characterization of two related Erwinia myoviruses that are distant relatives of the PhiKZ-like Jumbo phages

Bacteriophages are a major force in the evolution of bacteria due to their sheer abundance as well as their ability to infect and kill their hosts and to transfer genetic material. Bacteriophages that infect the Enterobacteriaceae family are of particular interest because this bacterial family contains dangerous animal and plant pathogens. Herein we report the isolation and characterization of two jumbo myovirus Erwinia phages, RisingSun and Joad, collected from apple trees. These two genomes are nearly identical with Joad harboring two additional putative gene products. Despite mass spectrometry data that support the putative annotation, 43% of their gene products have no significant BLASTP hit. These phages are also more closely related to Pseudomonas and Vibrio phages than to published Enterobacteriaceae phages. Of the 140 gene products with a BLASTP hit, 81% and 63% of the closest hits correspond to gene products from Pseudomonas and Vibrio phages, respectively. This relatedness may reflect their ecological niche, rather than the evolutionary history of their host. Despite the presence of over 800 Enterobacteriaceae phages on NCBI, the uniqueness of these two phages highlights the diversity of Enterobacteriaceae phages still to be discovered.

Searching whole genome sequences for biochemical identification features of emerging and reemerging pathogenic Corynebacterium species

Biochemical tests are traditionally used for bacterial identification at the species level in clinical microbiology laboratories. While biochemical profiles are generally efficient for the identification of the most important corynebacterial pathogen Corynebacterium diphtheriae, their ability to differentiate between biovars of this bacterium is still controversial. Besides, the unambiguous identification of emerging human pathogenic species of the genus Corynebacterium may be hampered by highly variable biochemical profiles commonly reported for these species, including Corynebacterium striatum, Corynebacterium amycolatum, Corynebacterium minutissimum, and Corynebacterium xerosis. In order to identify the genomic basis contributing for the biochemical variabilities observed in phenotypic identification methods of these bacteria, we combined a comprehensive literature review with a bioinformatics approach based on reconstruction of six specific biochemical reactions/pathways in 33 recently released whole genome sequences. We used data retrieved from curated databases (MetaCyc, PathoSystems Resource Integration Center (PATRIC), The SEED, TransportDB, UniProtKB) associated with homology searches by BLAST and profile Hidden Markov Models (HMMs) to detect enzymes participating in the various pathways and performed ab initio protein structure modeling and molecular docking to confirm specific results. We found a differential distribution among the various strains of genes that code for some important enzymes, such as beta-phosphoglucomutase and fructokinase, and also for individual components of carbohydrate transport systems, including the fructose-specific phosphoenolpyruvate-dependent sugar phosphotransferase (PTS) and the ribose-specific ATP-binging cassette (ABC) transporter. Horizontal gene transfer plays a role in the biochemical variability of the isolates, as some genes needed for sucrose fermentation were seen to be present in genomic islands. Noteworthy, using profile HMMs, we identified an enzyme with putative alpha-1,6-glycosidase activity only in some specific strains of C. diphtheriae and this may aid to understanding of the differential abilities to utilize glycogen and starch between the biovars.

PixelDB: Protein-peptide complexes annotated with structural conservation of the peptide binding mode

PixelDB, the Peptide Exosite Location Database, compiles 1966 non-redundant, high-resolution structures of protein-peptide complexes filtered to minimize the impact of crystal packing on peptide conformation. The database is organized to facilitate study of structurally conserved versus non-conserved elements of protein-peptide engagement. PixelDB clusters complexes based on the structural similarity of the peptide-binding protein, and by comparing complexes within a cluster highlights examples of domains that engage peptides using more than one binding mode. PixelDB also identifies conserved peptide core structural motifs characteristic of each binding mode. Peptide regions that flank core motifs often make non-structurally conserved interactions with the protein surface in regions we call exosites. Many examples establish that exosite contacts can be important for enhancing protein binding and interaction specificity. PixelDB provides a resource for computational and structural biologists to study, model, and predict core-motif and exosite-contacting peptide interactions. PixelDB is available to the community without restriction in a convenient flat-file format with accompanying visualization tools.

The evolution of function within the Nudix homology clan

The Nudix homology clan encompasses over 80,000 protein domains from all three domains of life, defined by homology to each other. Proteins with a domain from this clan fall into four general functional classes: pyrophosphohydrolases, isopentenyl diphosphate isomerases (IDIs), adenine/guanine mismatch-specific adenine glycosylases (A/G-specific adenine glycosylases), and nonenzymatic activities such as protein/protein interaction and transcriptional regulation. The largest group, pyrophosphohydrolases, encompasses more than 100 distinct hydrolase specificities. To understand the evolution of this vast number of activities, we assembled and analyzed experimental and structural data for 205 Nudix proteins collected from the literature. We corrected erroneous functions or provided more appropriate descriptions for 53 annotations described in the Gene Ontology Annotation database in this family, and propose 275 new experimentally-based annotations. We manually constructed a structure-guided sequence alignment of 78 Nudix proteins. Using the structural alignment as a seed, we then made an alignment of 347 "select" Nudix homology domains, curated from structurally determined, functionally characterized, or phylogenetically important Nudix domains. Based on our review of Nudix pyrophosphohydrolase structures and specificities, we further analyzed a loop region downstream of the Nudix hydrolase motif previously shown to contact the substrate molecule and possess known functional motifs. This loop region provides a potential structural basis for the functional radiation and evolution of substrate specificity within the hydrolase family. Finally, phylogenetic analyses of the 347 select protein domains and of the complete Nudix homology clan revealed general monophyly with regard to function and a few instances of probable homoplasy. Proteins 2017; 85:775-811. © 2016 Wiley Periodicals, Inc.

Structural divergence of essential triad ribbon synapse proteins among placental mammals - Implications for preclinical trials in photoreceptor transplantation therapy

As photoreceptor transplantation rapidly moves closer to the clinic, verifying graft efficacy in animal models may have unforeseen xenogeneic barriers. Although photoreceptor transplants have most convincingly exhibited functional synaptogenesis in conspecific studies, such evidence (while ruling out false-positives due to: viral graft labeling, fusion/cytosolic transfer, or neuroprotection) has not yet been shown for discordant xenografts. From this, a fundamental question should be raised: is useful xenosynaptogenesis likely between human photoreceptors and mouse retina? The triad ribbon synapse (TRS) that would normally form is unique and contains trans-synaptic proteins essential to its formation and function. Thus, could interspecific structural divergence be present that may inhibit this trans-synaptic bridge in discordant xenografts? In an effort to address this question computationally, we compared eight recently confirmed (including subcellular location) TRS specific (or predominantly expressed at the TRS) proteins among placental mammals (1-to-1 orthologs) using HyPhy selection analysis (a predictive measure of structural divergence) and by using Phyre2 tertiary structural modeling. Here, selection analysis revealed strong positive (diversifying) selection acting on a particularly important TRS protein: pikachurin. This positive selection was localized to its second Laminin-G (LG)-like domain and on its N-terminal domain - a putative region of trans-synaptic interaction. Localization of structural divergence to the N-terminus of each putative post-translational cleavage (PTC) product may suggest neofunctionalization from ancestral uncleaved pikachurin - this would be consistent with a recent counter-paradigm report of pikachurin cleavage predominating at the TRS. From this, we suggest a dual role after cleavage where the N-terminal fragment can still mediate the trans-synaptic bridge, while the C-terminal fragment may act as a diffusible trophic or "homing" factor for bipolar cell dendrite migration. Tertiary structural models mirrored the conformational divergence predicted by selection analysis. With human and mouse pikachurin (as well as other TRS proteins) likely to diverge considerably in structure among placental mammals - alongside known inter-mammalian variation in TRS phenotype and protein repertoire, high levels of diversifying selection acting on genes involving sensation, considerable timespans allowing for genetic drift that can create xenogeneic epistasis, and uncertainty surrounding the extent of xenosynaptogenesis in PPC transplant studies to date - use of distantly related hosts to test human photoreceptor graft therapeutic efficacy should be considered with caution.

Annotating Protein Functional Residues by Coupling High-Throughput Fitness Profile and Homologous-Structure Analysis

Identification and annotation of functional residues are fundamental questions in protein sequence analysis. Sequence and structure conservation provides valuable information to tackle these questions. It is, however, limited by the incomplete sampling of sequence space in natural evolution. Moreover, proteins often have multiple functions, with overlapping sequences that present challenges to accurate annotation of the exact functions of individual residues by conservation-based methods. Using the influenza A virus PB1 protein as an example, we developed a method to systematically identify and annotate functional residues. We used saturation mutagenesis and high-throughput sequencing to measure the replication capacity of single nucleotide mutations across the entire PB1 protein. After predicting protein stability upon mutations, we identified functional PB1 residues that are essential for viral replication. To further annotate the functional residues important to the canonical or noncanonical functions of viral RNA-dependent RNA polymerase (vRdRp), we performed a homologous-structure analysis with 16 different vRdRp structures. We achieved high sensitivity in annotating the known canonical polymerase functional residues. Moreover, we identified a cluster of noncanonical functional residues located in the loop region of the PB1 β-ribbon. We further demonstrated that these residues were important for PB1 protein nuclear import through the interaction with Ran-binding protein 5. In summary, we developed a systematic and sensitive method to identify and annotate functional residues that are not restrained by sequence conservation. Importantly, this method is generally applicable to other proteins about which homologous-structure information is available.

IMPORTANCE

To fully comprehend the diverse functions of a protein, it is essential to understand the functionality of individual residues. Current methods are highly dependent on evolutionary sequence conservation, which is usually limited by sampling size. Sequence conservation-based methods are further confounded by structural constraints and multifunctionality of proteins. Here we present a method that can systematically identify and annotate functional residues of a given protein. We used a high-throughput functional profiling platform to identify essential residues. Coupling it with homologous-structure comparison, we were able to annotate multiple functions of proteins. We demonstrated the method with the PB1 protein of influenza A virus and identified novel functional residues in addition to its canonical function as an RNA-dependent RNA polymerase. Not limited to virology, this method is generally applicable to other proteins that can be functionally selected and about which homologous-structure information is available.

Crystal Structures of SgcE6 and SgcC, the Two-Component Monooxygenase That Catalyzes Hydroxylation of a Carrier Protein-Tethered Substrate during the Biosynthesis of the Enediyne Antitumor Antibiotic C-1027 in Streptomyces globisporus

C-1027 is a chromoprotein enediyne antitumor antibiotic produced by Streptomyces globisporus. In the last step of biosynthesis of the (S)-3-chloro-5-hydroxy-β-tyrosine moiety of the C-1027 enediyne chromophore, SgcE6 and SgcC compose a two-component monooxygenase that hydroxylates the C-5 position of (S)-3-chloro-β-tyrosine. This two-component monooxygenase is remarkable for two reasons. (i) SgcE6 specifically reacts with FAD and NADH, and (ii) SgcC is active with only the peptidyl carrier protein (PCP)-tethered substrate. To address the molecular details of substrate specificity, we determined the crystal structures of SgcE6 and SgcC at 1.66 and 2.63 Å resolution, respectively. SgcE6 shares a similar β-barrel fold with the class I HpaC-like flavin reductases. A flexible loop near the active site of SgcE6 plays a role in FAD binding, likely by providing sufficient space to accommodate the AMP moiety of FAD, when compared to that of FMN-utilizing homologues. SgcC shows structural similarity to a few other known FADH₂-dependent monooxygenases and sheds light on some biochemically but not structurally characterized homologues. The crystal structures reported here provide insights into substrate specificity, and comparison with homologues provides a catalytic mechanism of the two-component, FADH₂-dependent monooxygenase (SgcE6 and SgcC) that catalyzes the hydroxylation of a PCP-tethered substrate.

Flavoprotein-Mediated Tellurite Reduction: Structural Basis and Applications to the Synthesis of Tellurium-Containing Nanostructures

The tellurium oxyanion tellurite (TeO₃^2-) is extremely harmful for most organisms. It has been suggested that a potential bacterial tellurite resistance mechanism would consist of an enzymatic, NAD(P)H-dependent, reduction to the less toxic form elemental tellurium (Te⁰). To date, a number of enzymes such as catalase, type II NADH dehydrogenase and terminal oxidases from the electron transport chain, nitrate reductases, and dihydrolipoamide dehydrogenase (E3), among others, have been shown to display tellurite-reducing activity. This activity is generically referred to as tellurite reductase (TR). Bioinformatic data resting on some of the abovementioned enzymes enabled the identification of common structures involved in tellurite reduction including vicinal catalytic cysteine residues and the FAD/NAD(P)⁺-binding domain, which is characteristic of some flavoproteins. Along this line, thioredoxin reductase (TrxB), alkyl hydroperoxide reductase (AhpF), glutathione reductase (GorA), mercuric reductase (MerA), NADH: flavorubredoxin reductase (NorW), dihydrolipoamide dehydrogenase, and the putative oxidoreductase YkgC from Escherichia coli or environmental bacteria were purified and assessed for TR activity. All of them displayed in vitro TR activity at the expense of NADH or NADPH oxidation. In general, optimal reducing conditions occurred around pH 9–10 and 37°C. Enzymes exhibiting strong TR activity produced Te-containing nanostructures (TeNS). While GorA and AhpF generated TeNS of 75 nm average diameter, E3 and YkgC produced larger structures (>100 nm). Electron-dense structures were observed in cells over-expressing genes encoding TrxB, GorA, and YkgC.

Calculating and scoring high quality multiple flexible protein structure alignments

MOTIVATION

Calculating multiple protein structure alignments (MSAs) is important for understanding functional and evolutionary relationships between protein families, and for modeling protein structures by homology. While incorporating backbone flexibility promises to circumvent many of the limitations of rigid MSA algorithms, very few flexible MSA algorithms exist today. This article describes several novel improvements to the Kpax algorithm which allow high quality flexible MSAs to be calculated. This article also introduces a new Gaussian-based MSA quality measure called 'M-score', which circumvents the pitfalls of RMSD-based quality measures.

RESULTS

As well as calculating flexible MSAs, the new version of Kpax can also score MSAs from other aligners and from previously aligned reference datasets. Results are presented for a large-scale evaluation of the Homstrad, SABmark and SISY benchmark sets using Kpax and Matt as examples of state-of-the-art flexible aligners and 3DCOMB as an example of a state-of-the-art rigid aligner. These results demonstrate the utility of the M-score as a measure of MSA quality and show that high quality MSAs may be achieved when structural flexibility is properly taken into account.

AVAILABILITY AND IMPLEMENTATION

Kpax 5.0 may be downloaded for academic use at http://kpax.loria.fr/

CONTACT

dave.ritchie@inria.fr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

TAPO: A combined method for the identification of tandem repeats in protein structures

In recent years, there has been an emergence of new 3D structures of proteins containing tandem repeats (TRs), as a result of improved expression and crystallization strategies. Databases focused on structure classifications (PDB, SCOP, CATH) do not provide an easy solution for selection of these structures from PDB. Several approaches have been developed, but no best approach exists to identify the whole range of 3D TRs. Here we describe the TAndem PrOtein detector (TAPO) that uses periodicities of atomic coordinates and other types of structural representation, including strings generated by conformational alphabets, residue contact maps, and arrangements of vectors of secondary structure elements. The benchmarking shows the superior performance of TAPO over the existing programs. In accordance with our analysis of PDB using TAPO, 19% of proteins contain 3D TRs. This analysis allowed us to identify new families of 3D TRs, suggesting that TAPO can be used to regularly update the collection and classification of existing repetitive structures.

"PP2C7s", Genes Most Highly Elaborated in Photosynthetic Organisms, Reveal the Bacterial Origin and Stepwise Evolution of PPM/PP2C Protein Phosphatases

Mg⁺²/Mn⁺²-dependent type 2C protein phosphatases (PP2Cs) are ubiquitous in eukaryotes, mediating diverse cellular signaling processes through metal ion catalyzed dephosphorylation of target proteins. We have identified a distinct PP2C sequence class (“PP2C7s”) which is nearly universally distributed in Eukaryotes, and therefore apparently ancient. PP2C7s are by far most prominent and diverse in plants and green algae. Combining phylogenetic analysis, subcellular localization predictions, and a distillation of publically available gene expression data, we have traced the evolutionary trajectory of this gene family in photosynthetic eukaryotes, demonstrating two major sequence assemblages featuring a succession of increasingly derived sub-clades. These display predominant expression moving from an ancestral pattern in photosynthetic tissues toward non-photosynthetic, specialized and reproductive structures. Gene co-expression network composition strongly suggests a shifting pattern of PP2C7 gene functions, including possible regulation of starch metabolism for one homologue set in Arabidopsis and rice. Distinct plant PP2C7 sub-clades demonstrate novel amino terminal protein sequences upon motif analysis, consistent with a shifting pattern of regulation of protein function. More broadly, neither the major events in PP2C sequence evolution, nor the origin of the diversity of metal binding characteristics currently observed in different PP2C lineages, are clearly understood. Identification of the PP2C7 sequence clade has allowed us to provide a better understanding of both of these issues. Phylogenetic analysis and sequence comparisons using Hidden Markov Models strongly suggest that PP2Cs originated in Bacteria (Group II PP2C sequences), entered Eukaryotes through the ancestral mitochondrial endosymbiosis, elaborated in Eukaryotes, then re-entered Bacteria through an inter-domain gene transfer, ultimately producing bacterial Group I PP2C sequences. A key evolutionary event, occurring first in ancient Eukaryotes, was the acquisition of a conserved aspartate in classic Motif 5. This has been inherited subsequently by PP2C7s, eukaryotic PP2Cs and bacterial Group I PP2Cs, where it is crucial to the formation of a third metal binding pocket, and catalysis.

AcconPred: Predicting Solvent Accessibility and Contact Number Simultaneously by a Multitask Learning Framework under the Conditional Neural Fields Model

Motivation. The solvent accessibility of protein residues is one of the driving forces of protein folding, while the contact number of protein residues limits the possibilities of protein conformations. The de novo prediction of these properties from protein sequence is important for the study of protein structure and function. Although these two properties are certainly related with each other, it is challenging to exploit this dependency for the prediction. Method. We present a method AcconPred for predicting solvent accessibility and contact number simultaneously, which is based on a shared weight multitask learning framework under the CNF (conditional neural fields) model. The multitask learning framework on a collection of related tasks provides more accurate prediction than the framework trained only on a single task. The CNF method not only models the complex relationship between the input features and the predicted labels, but also exploits the interdependency among adjacent labels. Results. Trained on 5729 monomeric soluble globular protein datasets, AcconPred could reach 0.68 three-state accuracy for solvent accessibility and 0.75 correlation for contact number. Tested on the 105 CASP11 domain datasets for solvent accessibility, AcconPred could reach 0.64 accuracy, which outperforms existing methods.

DeepCNF-D: Predicting Protein Order/Disorder Regions by Weighted Deep Convolutional Neural Fields

Intrinsically disordered proteins or protein regions are involved in key biological processes including regulation of transcription, signal transduction, and alternative splicing. Accurately predicting order/disorder regions ab initio from the protein sequence is a prerequisite step for further analysis of functions and mechanisms for these disordered regions. This work presents a learning method, weighted DeepCNF (Deep Convolutional Neural Fields), to improve the accuracy of order/disorder prediction by exploiting the long-range sequential information and the interdependency between adjacent order/disorder labels and by assigning different weights for each label during training and prediction to solve the label imbalance issue. Evaluated by the CASP9 and CASP10 targets, our method obtains 0.855 and 0.898 AUC values, which are higher than the state-of-the-art single ab initio predictors.

Defining the cytosolic pathway of glutathione degradation in Arabidopsis thaliana: role of the ChaC/GCG family of γ-glutamyl cyclotransferases as glutathione-degrading enzymes and AtLAP1 as the Cys-Gly peptidase

Glutathione homoeostasis is critical to plant life and its adaptation to stress. The γ-glutamyl cycle of glutathione biosynthesis and degradation plays a pre-eminent role in glutathione homoeostasis. The genes encoding two enzymatic steps of glutathione degradation, the γ-glutamyl cyclotransferase (GGCT; acting on γ-glutamyl amino acids) and the Cys-Gly dipeptidase, have, however, lacked identification. We have investigated the family of GGCTs in Arabidopsis thaliana. We show through in vivo functional assays in yeast that all three members of the ChaC/GCG subfamily show significant activity towards glutathione but no detectable activity towards γ-glutamyl methionine. Biochemical characterization of the purified recombinant enzymes GGCT2;2 and GGCT2;3 further confirmed that they act specifically to degrade glutathione to yield 5-oxoproline and Cys-Gly peptide and show no significant activity towards γ-glutamyl cysteine. The Km for glutathione was 1.7 and 4.96 mM for GGCT2;2 and GGCT2;3 respectively and was physiologically relevant. Evaluation of representative members of other subfamilies indicates the absence of GGCTs from plants showing significant activity towards γ-glutamyl-amino acids as envisaged in the classical γ-glutamyl cycle. To identify the Cys-Gly peptidase, we evaluated leucine aminopeptidases (LAPs) as candidate enzymes. The cytosolic AtLAP1 (A. thaliana leucine aminopeptidase 1) and the putative chloroplastic AtLAP3 displayed activity towards Cys-Gly peptide through in vivo functional assays in yeast. Biochemical characterization of the in vitro purified hexameric AtLAP1 enzyme revealed a Km for Cys-Gly of 1.3 mM that was physiologically relevant and indicated that AtLAP1 represents a cytosolic Cys-Gly peptidase activity of A. thaliana. The studies provide new insights into the functioning of the γ-glutamyl cycle in plants.

POSA: a user-driven, interactive multiple protein structure alignment server

POSA (Partial Order Structure Alignment), available at http://posa.godziklab.org, is a server for multiple protein structure alignment introduced in 2005 (Ye,Y. and Godzik,A. (2005) Multiple flexible structure alignment using partial order graphs. Bioinformatics, 21, 2362–2369). It is free and open to all users, and there is no login requirement, albeit there is an option to register and store results in individual, password-protected directories. In the updated POSA server described here, we introduce two significant improvements. First is an interface allowing the user to provide additional information by defining segments that anchor the alignment in one or more input structures. This interface allows users to take advantage of their intuition and biological insights to improve the alignment and guide it toward a biologically relevant solution. The second improvement is an interactive visualization with options that allow the user to view all superposed structures in one window (a typical solution for visualizing results of multiple structure alignments) or view them individually in a series of synchronized windows with extensive, user-controlled visualization options. The user can rotate structure(s) in any of the windows and study similarities or differences between structures clearly visible in individual windows.

A local average distance descriptor for flexible protein structure comparison

BACKGROUND

Protein structures are flexible and often show conformational changes upon binding to other molecules to exert biological functions. As protein structures correlate with characteristic functions, structure comparison allows classification and prediction of proteins of undefined functions. However, most comparison methods treat proteins as rigid bodies and cannot retrieve similarities of proteins with large conformational changes effectively.

RESULTS

In this paper, we propose a novel descriptor, local average distance (LAD), based on either the geodesic distances (GDs) or Euclidean distances (EDs) for pairwise flexible protein structure comparison. The proposed method was compared with 7 structural alignment methods and 7 shape descriptors on two datasets comprising hinge bending motions from the MolMovDB, and the results have shown that our method outperformed all other methods regarding retrieving similar structures in terms of precision-recall curve, retrieval success rate, R-precision, mean average precision and F1-measure.

CONCLUSIONS

Both ED- and GD-based LAD descriptors are effective to search deformed structures and overcome the problems of self-connection caused by a large bending motion. We have also demonstrated that the ED-based LAD is more robust than the GD-based descriptor. The proposed algorithm provides an alternative approach for blasting structure database, discovering previously unknown conformational relationships, and reorganizing protein structure classification.

Protein structure alignment beyond spatial proximity

Protein structure alignment is a fundamental problem in computational structure biology. Many programs have been developed for automatic protein structure alignment, but most of them align two protein structures purely based upon geometric similarity without considering evolutionary and functional relationship. As such, these programs may generate structure alignments which are not very biologically meaningful from the evolutionary perspective. This paper presents a novel method DeepAlign for automatic pairwise protein structure alignment. DeepAlign aligns two protein structures using not only spatial proximity of equivalent residues (after rigid-body superposition), but also evolutionary relationship and hydrogen-bonding similarity. Experimental results show that DeepAlign can generate structure alignments much more consistent with manually-curated alignments than other automatic tools especially when proteins under consideration are remote homologs. These results imply that in addition to geometric similarity, evolutionary information and hydrogen-bonding similarity are essential to aligning two protein structures.

mulPBA: an efficient multiple protein structure alignment method based on a structural alphabet

The increasing number of available protein structures requires efficient tools for multiple structure comparison. Indeed, multiple structural alignments are essential for the analysis of function, evolution and architecture of protein structures. For this purpose, we proposed a new web server called multiple Protein Block Alignment (mulPBA). This server implements a method based on a structural alphabet to describe the backbone conformation of a protein chain in terms of dihedral angles. This 'sequence-like' representation enables the use of powerful sequence alignment methods for primary structure comparison, followed by an iterative refinement of the structural superposition. This approach yields alignments superior to most of the rigid-body alignment methods and highly comparable with the flexible structure comparison approaches. We implement this method in a web server designed to do multiple structure superimpositions from a set of structures given by the user. Outputs are given as both sequence alignment and superposed 3D structures visualized directly by static images generated by PyMol or through a Jmol applet allowing dynamic interaction. Multiple global quality measures are given. Relatedness between structures is indicated by a distance dendogram. Superimposed structures in PDB format can be also downloaded, and the results are quickly obtained. mulPBA server can be accessed at www.dsimb.inserm.fr/dsimb_tools/mulpba/ .

Fast protein binding site comparison via an index-based screening technology

We present TrixP, a new index-based method for fast protein binding site comparison and function prediction. TrixP determines binding site similarities based on the comparison of descriptors that encode pharmacophoric and spatial features. Therefore, it adopts the efficient core components of TrixX, a structure-based virtual screening technology for large compound libraries. TrixP expands this technology by new components in order to allow a screening of protein libraries. TrixP accounts for the inherent flexibility of proteins employing a partial shape matching routine. After the identification of structures with matching pharmacophoric features and geometric shape, TrixP superimposes the binding sites and, finally, assesses their similarity according to the fit of pharmacophoric properties. TrixP is able to find analogies between closely and distantly related binding sites. Recovery rates of 81.8% for similar binding site pairs, assisted by rejecting rates of 99.5% for dissimilar pairs on a test data set containing 1331 pairs, confirm this ability. TrixP exclusively identifies members of the same protein family on top ranking positions out of a library consisting of 9802 binding sites. Furthermore, 30 predicted kinase binding sites can almost perfectly be classified into their known subfamilies.