Prediction of CASP6 structures using automated Robetta protocols

Computational Vaccine Design for Poxviridae Family Viruses

The computational approach to designing vaccines has several useful characteristics over traditional vaccine development, such as being highly specific, less time-consuming and less expensive. Thus, this chapter describes an immunoinformatics workflow to design a vaccine against a member of the Poxviridae family known as Monkeypox virus. The immunoinformatics approach uses several online servers to select highly antigenic and non-allergenic CTL, HTL, and B cell epitopes. Then, it links the predicted epitopes through linkers and submit them for 3D structure modeling. Afterward, the modeled vaccine is docked with TLRs to check the induction of the immune system. Finally, immune simulations are performed to check the level of several immune factors like IgG, IgM, cytokines and interleukins, among others, upon the injection of the constructed vaccine. This approach can be used to successfully design novel and effective vaccine candidates against emerging species from the Poxviridae family.

Cassava Brown Streak Viruses express second 6-kilodalton (6K2) protein with varied polarity and three dimensional (3D) structures: Basis for trait discrepancy between the virus species

Cassava Brown Streak Virus (CBSV) and Ugandan Cassava Brown Streak Virus (UCBSV) are the two among six virus species speculated to cause the most catastrophic Brown Streak Disease of Cassava (CBSD) in Africa and Asia. Cassava Brown Streak Virus (CBSV) is hard to breed resistance for compared to Ugandan Cassava Brown Streak Virus (UCBSV) species. This is exemplified by incidences of CBSV species rather than UCBSV species in elite breeding line, KBH 2006/0026 at Bagamoyo, Tanzania. It is not yet understood as to why CBSV species could breakdown CBSD-resistance in the KBH 2006/0026 unlike the UCBSV species. This marks the first in silico study conducted to understand molecular basis for the trait discrepancy between CBSV and UCBSV species from structural biology view point. Following ab initio modelling and analysis of physical-chemical properties of second 6-kilodalton (6K2) protein encoded by CBSV and UCBSV species, using ROBETTA server and Protein Parameters tool, respectively we report that; three dimensional (3D) structures and polarity of the protein differs significantly between the two virus species. (95% and 5%) and (85% and 15%) strains of 20 CBSV and 20 UCBSV species respectively, expressed the protein in homo-trimeric and homo-tetrameric forms, correspondingly. 95% and 85% of studied strain population of the two virus species expressed hydrophilic and hydrophobic 6K2, respectively. Based on findings of the curent study, we hypothesize that; (i) The hydrophilic 6K2 expressed by the CBSV species, favour its faster systemic movement via vascular tissues of cassava host and hence result into higher tissue titres than the UCBSV species encoding hydrophobic form of the protein. t and (ii) The hydrophilic 6K2 expressed byCBSV species have additional interaction advantage with Nuclear Inclusion b protease domain (NIb) and Viral genome-linked protein (VPg), components of Virus Replication Complex (VRC) and hence contributing to faster replication of viral genome than the hydrophobic 6K2 expressed by the UCBSV species. Experimental studies are needed to resolve the 3D structures of the 6K2, VPg and NIb and comprehend complex molecular interactions between them. We suggest that, the 6K2 gene should be targeted for improvement of RNA interference (RNAi)-directed transgenesis of virus-resistant cassava as a more effective way to control the CBSD besides breeding.

Annotation of Potential Vaccine Targets and Design of a Multi-Epitope Subunit Vaccine against Yersinia pestis through Reverse Vaccinology and Validation through an Agent-Based Modeling Approach

Yersinia pestis is responsible for plague and major pandemics in Asia and Europe. This bacterium has shown resistance to an array of drugs commonly used for the treatment of plague. Therefore, effective therapeutics measurements, such as designing a vaccine that can effectively and safely prevent Y. pestis infection, are of high interest. To fast-track vaccine development against Yersinia pestis, herein, proteome-wide vaccine target annotation was performed, and structural vaccinology-assisted epitopes were predicted. Among the total 3909 proteins, only 5 (rstB, YPO2385, hmuR, flaA1a, and psaB) were shortlisted as essential vaccine targets. These targets were then subjected to multi-epitope vaccine design using different linkers. EAAK, AAY, and GPGPG as linkers were used to link CTL, HTL, and B-cell epitopes, and an adjuvant (beta defensin) was also added at the N-terminal of the MEVC. Physiochemical characterization, such as determination of the instability index, theoretical pI, half-life, aliphatic index, stability profiling, antigenicity, allergenicity, and hydropathy of the ensemble, showed that the vaccine is highly stable, antigenic, and non-allergenic and produces multiple interactions with immune receptors upon docking. In addition, molecular dynamics simulation confirmed the stable binding and good dynamic properties of the vaccine-TLR complex. Furthermore, in silico and immune simulation of the developed MEVC for Y. pestis showed that the vaccine triggered strong immune response after several doses at different intervals. Neutralization of the antigen was observed at the third day of injection. Conclusively, the vaccine designed here for Y. pestis produces an immune response; however, further immunological testing is needed to unveil its real efficacy.

Identification of Immunogenic Linear B-Cell Epitopes in C. burnetii Outer Membrane Proteins Using Immunoinformatics Approaches Reveals Potential Targets of Persistent Infections

Coxiella burnetii is a global, highly infectious intracellular bacterium, able to infect a wide range of hosts and to persist for months in the environment. It is the etiological agent of Q fever—a zoonosis of global priority. Currently, there are no national surveillance data on C. burnetii’s seroprevalence for any South American country, reinforcing the necessity of developing novel and inexpensive serological tools to monitor the prevalence of infections among humans and animals—especially cattle, goats, and sheep. In this study, we used immunoinformatics and computational biology tools to predict specific linear B-cell epitopes in three C. burnetii outer membrane proteins: OMP-H (CBU_0612), Com-1 (CBU_1910), and OMP-P1 (CBU_0311). Furthermore, predicted epitopes were tested by ELISA, as synthetic peptides, against samples of patients reactive to C. burnetii in indirect immunofluorescence assay, in order to evaluate their natural immunogenicity. In this way, two linear B-cell epitopes were identified in each studied protein (OMP-H_(51–59), OMP-H_(91–106), Com-1_(57–76), Com-1_(191–206), OMP-P1_(197–209), and OMP-P1_(215–227)); all of them were confirmed as naturally immunogenic by the presence of specific antibodies in 77% of studied patients against at least one of the identified epitopes. Remarkably, a higher frequency of endocarditis cases was observed among patients who presented an intense humoral response to OMP-H and Com-1 epitopes. These data confirm that immunoinformatics applied to the identification of specific B-cell epitopes can be an effective strategy to improve and accelerate the development of surveillance tools against neglected diseases.

Neurotoxic Effects of Linalool and β-Pinene on Tribolium castaneum Herbst

Effective, ethical pest control requires the use of chemicals that are highly specific, safe, and ecofriendly. Linalool and β-pinene occur naturally as major constituents of the essential oils of many plant species distributed throughout the world, and thus meet these requirements. These monoterpenes were tested as repellents against Tribolium castaneum, using the area preference method, after four hours of exposure and the effect transcriptional of genes associated with neurotransmission. Changes in gene expression of acetylcholinesterase (Ace1), GABA-gated anion channel splice variant 3a6a (Rdl), GABA-gated ion channel (Grd), glutamate-gated chloride channel (Glucl), and histamine-gated chloride channel 2 (Hiscl2) were assessed and the interaction with proteins important for the insect using in silico methods was also studied. For linalool and β-pinene, the repellent concentration 50 (RC₅₀) values were 0.11 µL/cm² and 0.03 µL/cm², respectively. Both compounds induced overexpression of Hiscl2 gen in adult insects, and β-pinene also promoted the overexpression of Grd and the Ace1 gene. However, β-pinene and linalool had little potential to dock on computer-generated models for GABA-gated ion channel LCCH3, nicotinic acetylcholine receptor subunits alpha1 and alpha2, and putative octopamine/tyramine receptor proteins from T. castaneum as their respective binding affinities were marginal, and therefore the repellent action probably involved mechanisms other than direct interaction with these targets. Results indicated that β-pinene was more potent than linalool in inducing insect repellency, and also had a greater capacity to generate changes in the expression of genes involved in neuronal transmission.

Metal Selectivity of a Cd-, Co-, and Zn-Transporting P1B-type ATPase

The P_1B-ATPases, a family of transmembrane metal transporters important for transition metal homeostasis in all organisms, are subdivided into classes based on sequence conservation and metal specificity. The multifunctional P_1B-4-ATPase CzcP is part of the cobalt, zinc, and cadmium resistance system from the metal-tolerant, model organism Cupriavidus metallidurans. Previous work revealed the presence of an unusual soluble metal-binding domain (MBD) at the CzcP N-terminus, but the nature, extent, and selectivity of the transmembrane metal-binding site (MBS) of CzcP have not been resolved. Using homology modeling, we show that four wholly conserved amino acids from the transmembrane (TM) domain (Met²⁵⁴, Ser⁴⁷⁴, Cys⁴⁷⁶, and His⁸⁰⁷) are logical candidates for the TM MBS, which may communicate with the MBD via interactions with the first TM helix. Metal-binding analyses indicate that wild-type (WT) CzcP has three MBSs, and data on N-terminally truncated (ΔMBD) CzcP suggest the presence of a single TM MBS. Electronic absorption and electron paramagnetic resonance spectroscopic analyses of ΔMBD CzcP and variant proteins thereof provide insight into the details of Co²⁺ coordination by the TM MBS. These spectroscopic data, combined with in vitro functional studies of WT and variant CzcP proteins, show that the side chains of Met²⁵⁴, Cys⁴⁷⁶, and His⁸⁰⁷ contribute to Cd²⁺, Co²⁺, and Zn²⁺ binding and transport, whereas the side chain of Ser⁴⁷⁴ appears to play a minimal role. By comparison to other P_1B-4-ATPases, we suggest that an evolutionarily adapted flexibility in the TM region likely afforded CzcP the ability to transport Cd²⁺ and Zn²⁺ in addition to Co²⁺.

Crystal structure of the human sterol transporter ABCG5/ABCG8

ATP binding cassette (ABC) transporters play critical roles in maintaining sterol balance in higher eukaryotes. The ABCG5/ABCG8 heterodimer (G5G8) mediates excretion of neutral sterols in liver and intestines. Mutations disrupting G5G8 cause sitosterolaemia, a disorder characterized by sterol accumulation and premature atherosclerosis. Here we use crystallization in lipid bilayers to determine the X-ray structure of human G5G8 in a nucleotide-free state at 3.9 Å resolution, generating the first atomic model of an ABC sterol transporter. The structure reveals a new transmembrane fold that is present in a large and functionally diverse superfamily of ABC transporters. The transmembrane domains are coupled to the nucleotide-binding sites by networks of interactions that differ between the active and inactive ATPases, reflecting the catalytic asymmetry of the transporter. The G5G8 structure provides a mechanistic framework for understanding sterol transport and the disruptive effects of mutations causing sitosterolaemia.

Protein architecture and core residues in unwound α-helices provide insights to the transport function of plant AtCHX17

Using Arabidopsis thaliana AtCHX17 as an example, we combine structural modeling and mutagenesis to provide insights on its protein architecture and transport function which is poorly characterized. This approach is based on the observation that protein structures are significantly more conserved in evolution than linear sequences, and mechanistic similarities among diverse transporters are emerging. Two homology models of AtCHX17 were obtained that show a protein fold similar to known structures of bacterial Na(+)/H(+) antiporters, EcNhaA and TtNapA. The distinct secondary and tertiary structure models highlighted residues at positions potentially important for CHX17 activity. Mutagenesis showed that asparagine-N200 and aspartate-D201 inside transmembrane5 (TM5), and lysine-K355 inside TM10 are critical for AtCHX17 activity. We reveal previously unrecognized threonine-T170 and lysine-K383 as key residues at unwound regions in the middle of TM4 and TM11 α-helices, respectively. Mutation of glutamate-E111 located near the membrane surface inhibited AtCHX17 activity, suggesting a role in pH sensing. The long carboxylic tail of unknown purpose has an alternating β-sheet and α-helix secondary structure that is conserved in prokaryote universal stress proteins. These results support the overall architecture of AtCHX17 and identify D201, N200 and novel residues T170 and K383 at the functional core which likely participates in ion recognition, coordination and/or translocation, similar to characterized cation/H(+) exchangers. The core of AtCHX17 models according to EcNhaA and TtNapA templates faces inward and outward, respectively, which may reflect two conformational states of the alternating access transport mode for proteins belonging to the plant CHX family.

CASP10-BCL::Fold efficiently samples topologies of large proteins

During CASP10 in summer 2012, we tested BCL::Fold for prediction of free modeling (FM) and template-based modeling (TBM) targets. BCL::Fold assembles the tertiary structure of a protein from predicted secondary structure elements (SSEs) omitting more flexible loop regions early on. This approach enables the sampling of conformational space for larger proteins with more complex topologies. In preparation of CASP11, we analyzed the quality of CASP10 models throughout the prediction pipeline to understand BCL::Fold's ability to sample the native topology, identify native-like models by scoring and/or clustering approaches, and our ability to add loop regions and side chains to initial SSE-only models. The standout observation is that BCL::Fold sampled topologies with a GDT_TS score > 33% for 12 of 18 and with a topology score > 0.8 for 11 of 18 test cases de novo. Despite the sampling success of BCL::Fold, significant challenges still exist in clustering and loop generation stages of the pipeline. The clustering approach employed for model selection often failed to identify the most native-like assembly of SSEs for further refinement and submission. It was also observed that for some β-strand proteins model refinement failed as β-strands were not properly aligned to form hydrogen bonds removing otherwise accurate models from the pool. Further, BCL::Fold samples frequently non-natural topologies that require loop regions to pass through the center of the protein.

General overview on structure prediction of twilight-zone proteins

Protein structure prediction from amino acid sequence has been one of the most challenging aspects in computational structural biology despite significant progress in recent years showed by critical assessment of protein structure prediction (CASP) experiments. When experimentally determined structures are unavailable, the predictive structures may serve as starting points to study a protein. If the target protein consists of homologous region, high-resolution (typically <1.5 Å) model can be built via comparative modelling. However, when confronted with low sequence similarity of the target protein (also known as twilight-zone protein, sequence identity with available templates is less than 30%), the protein structure prediction has to be initiated from scratch. Traditionally, twilight-zone proteins can be predicted via threading or ab initio method. Based on the current trend, combination of different methods brings an improved success in the prediction of twilight-zone proteins. In this mini review, the methods, progresses and challenges for the prediction of twilight-zone proteins were discussed.

Use of a structural alphabet to find compatible folds for amino acid sequences

The structural annotation of proteins with no detectable homologs of known 3D structure identified using sequence-search methods is a major challenge today. We propose an original method that computes the conditional probabilities for the amino-acid sequence of a protein to fit to known protein 3D structures using a structural alphabet, known as "Protein Blocks" (PBs). PBs constitute a library of 16 local structural prototypes that approximate every part of protein backbone structures. It is used to encode 3D protein structures into 1D PB sequences and to capture sequence to structure relationships. Our method relies on amino acid occurrence matrices, one for each PB, to score global and local threading of query amino acid sequences to protein folds encoded into PB sequences. It does not use any information from residue contacts or sequence-search methods or explicit incorporation of hydrophobic effect. The performance of the method was assessed with independent test datasets derived from SCOP 1.75A. With a Z-score cutoff that achieved 95% specificity (i.e., less than 5% false positives), global and local threading showed sensitivity of 64.1% and 34.2%, respectively. We further tested its performance on 57 difficult CASP10 targets that had no known homologs in PDB: 38 compatible templates were identified by our approach and 66% of these hits yielded correctly predicted structures. This method scales-up well and offers promising perspectives for structural annotations at genomic level. It has been implemented in the form of a web-server that is freely available at http://www.bo-protscience.fr/forsa.

The structure and dynamics of BmR1 protein from Brugia malayi: in silico approaches

Brugia malayi is a filarial nematode, which causes lymphatic filariasis in humans. In 1995, the disease has been identified by the World Health Organization (WHO) as one of the second leading causes of permanent and long-term disability and thus it is targeted for elimination by year 2020. Therefore, accurate filariasis diagnosis is important for management and elimination programs. A recombinant antigen (BmR1) from the Bm17DIII gene product was used for antibody-based filariasis diagnosis in "Brugia Rapid". However, the structure and dynamics of BmR1 protein is yet to be elucidated. Here we study the three dimensional structure and dynamics of BmR1 protein using comparative modeling, threading and ab initio protein structure prediction. The best predicted structure obtained via an ab initio method (Rosetta) was further refined and minimized. A total of 5 ns molecular dynamics simulation were performed to investigate the packing of the protein. Here we also identified three epitopes as potential antibody binding sites from the molecular dynamics average structure. The structure and epitopes obtained from this study can be used to design a binder specific against BmR1, thus aiding future development of antigen-based filariasis diagnostics to complement the current diagnostics.

IS PROTEIN FOLDING PROBLEM REALLY A NP-COMPLETE ONE? FIRST INVESTIGATIONS

To determine the 3D conformation of proteins is a necessity to understand their functions or interactions with other molecules. It is commonly admitted that, when proteins fold from their primary linear structures to their final 3D conformations, they tend to choose the ones that minimize their free energy. To find the 3D conformation of a protein knowing its amino acid sequence, bioinformaticians use various models of different resolutions and artificial intelligence tools, as the protein folding prediction problem is a NP complete one. More precisely, to determine the backbone structure of the protein using the low resolution models (2D HP square and 3D HP cubic), by finding the conformation that minimizes free energy, is intractable exactly. Both proofs of NP-completeness and the 2D prediction consider that acceptable conformations have to satisfy a self-avoiding walk (SAW) requirement, as two different amino acids cannot occupy a same position in the lattice. It is shown in this document that the SAW requirement considered when proving NP-completeness is different from the SAW requirement used in various prediction programs, and that they are different from the real biological requirement. Indeed, the proof of NP completeness and the predictions in silico consider conformations that are not possible in practice. Consequences of this fact are investigated in this research work.

The sorting of proglucagon to secretory granules is mediated by carboxypeptidase E and intrinsic sorting signals

Proglucagon is expressed in pancreatic alpha cells, intestinal L cells and brainstem neurons. Tissue-specific processing of proglucagon yields the peptide hormones glucagon in the alpha cell and glucagon-like peptide (GLP)-1 and GLP-2 in L cells. Both glucagon and GLP-1 are secreted in response to nutritional status and are critical for regulating glycaemia. The sorting of proglucagon to the dense-core secretory granules of the regulated secretory pathway is essential for the appropriate secretion of glucagon and GLP-1. We examined the roles of carboxypeptidase E (CPE), a prohormone sorting receptor, the processing enzymes PC1/3 and PC2 and putative intrinsic sorting signals in proglucagon sorting. In Neuro 2a cells that lacked CPE, PC1/3 and PC2, proglucagon co-localised with the Golgi marker p115 as determined by quantitative immunofluorescence microscopy. Expression of CPE, but not of PC1/3 or PC2, enhanced proglucagon sorting to granules. siRNA-mediated knockdown of CPE disrupted regulated secretion of glucagon from pancreatic-derived alphaTC1-6 cells, but not of GLP-1 from intestinal cell-derived GLUTag cells. Mutation of the PC cleavage site K70R71, the dibasic R17R18 site within glucagon or the alpha-helix of glucagon, all significantly affected the sub-cellular localisation of proglucagon. Protein modelling revealed that alpha helices corresponding to glucagon, GLP-1 and GLP-2, are arranged within a disordered structure, suggesting some flexibility in the sorting mechanism. We conclude that there are multiple mechanisms for sorting proglucagon to the regulated secretory pathway, including a role for CPE in pancreatic alpha cells, initial cleavage at K70R71 and multiple sorting signals.

A dynamic model of the proteins that form the initial iron-sulfur cluster biogenesis machinery in yeast mitochondria

The assembly of iron-sulfur clusters (ISCs) in eukaryotes involves the protein Frataxin. Deficits in this protein have been associated with iron inside the mitochondria and impair ISC biogenesis as it is postulated to act as the iron donor for ISCs assembly in this organelle. A pronounced lack of Frataxin causes Friedreich's Ataxia, which is a human neurodegenerative and hereditary disease mainly affecting the equilibrium, coordination, muscles and heart. Moreover, it is the most common autosomal recessive ataxia. High similarities between the human and yeast molecular mechanisms that involve Frataxin have been suggested making yeast a good model to study that process. In yeast, the protein complex that forms the central assembly platform for the initial step of ISC biogenesis is composed by yeast frataxin homolog, Nfs1-Isd11 and Isu. In general, it is commonly accepted that protein function involves interaction with other protein partners, but in this case not enough is known about the structure of the protein complex and, therefore, how it exactly functions. The objective of this work is to model the protein complex in order to gain insight into structural details that end up with its biological function. To achieve this goal several bioinformatics tools, modeling techniques and protein docking programs have been used. As a result, the structure of the protein complex and the dynamic behavior of its components, along with that of the iron and sulfur atoms required for the ISC assembly, have been modeled. This hypothesis will help to better understand the function and molecular properties of Frataxin as well as those of its ISC assembly protein partners.

Convergence of IRBIT, phosphatidylinositol (4,5) bisphosphate, and WNK/SPAK kinases in regulation of the Na+-HCO3- cotransporters family

Fluid and HCO3(-) secretion is a vital function of secretory epithelia, involving basolateral HCO3(-) entry through the Na(+)-HCO3(-) cotransporter (NBC) NBCe1-B, and luminal HCO3(-) exit mediated by cystic fibrosis transmembrane conductance regulator (CFTR) and solute carrier family 26 (SLC26) Cl(-)/HCO3(-) exchangers. HCO3(-) secretion is highly regulated, with the WNK/SPAK kinase pathway setting the resting state and the IRBIT/PP1 pathway setting the stimulated state. However, we know little about the relationships between the WNK/SPAK and IRBIT/PP1 sites in the regulation of the transporters. The first 85 N-terminal amino acids of NBCe1-B function as an autoinhibitory domain. Here we have identified a positively charged module within NBCe1-B(37-65) that is conserved in NBCn1-A and all 20 members of the NBC superfamily except NBCe1-A. This module is required for the interaction and activation of NBCe1-B and NBCn1-A by IRBIT and their regulation by phosphatidylinositol 4,5-bisphosphate (PIP2). Activation of the transporters by IRBIT and PIP2 is nonadditive but complementary. Phosphorylation of Ser65 mediates regulation of NBCe1-B by SPAK, and phosphorylation of Thr49 is required for regulation by IRBIT and SPAK. Sequence searches using the NBCe1-B regulatory module as a template identified a homologous sequence in the CFTR R domain and Slc26a6 sulfat transporter and antisigma factor antagonist (STAS) domain. Accordingly, the R and STAS domains bind IRBIT, and the R domain is required for activation of CFTR by IRBIT. These findings reveal convergence of regulatory modalities in a conserved domain of the NBC that may be present in other HCO3(-) transporters and thus in the regulation of epithelial fluid and HCO3(-) secretion.

BCL::Fold--de novo prediction of complex and large protein topologies by assembly of secondary structure elements

Computational de novo protein structure prediction is limited to small proteins of simple topology. The present work explores an approach to extend beyond the current limitations through assembling protein topologies from idealized α-helices and β-strands. The algorithm performs a Monte Carlo Metropolis simulated annealing folding simulation. It optimizes a knowledge-based potential that analyzes radius of gyration, β-strand pairing, secondary structure element (SSE) packing, amino acid pair distance, amino acid environment, contact order, secondary structure prediction agreement and loop closure. Discontinuation of the protein chain favors sampling of non-local contacts and thereby creation of complex protein topologies. The folding simulation is accelerated through exclusion of flexible loop regions further reducing the size of the conformational search space. The algorithm is benchmarked on 66 proteins with lengths between 83 and 293 amino acids. For 61 out of these proteins, the best SSE-only models obtained have an RMSD100 below 8.0 Å and recover more than 20% of the native contacts. The algorithm assembles protein topologies with up to 215 residues and a relative contact order of 0.46. The method is tailored to be used in conjunction with low-resolution or sparse experimental data sets which often provide restraints for regions of defined secondary structure.

SIMONA 1.0: an efficient and versatile framework for stochastic simulations of molecular and nanoscale systems

Molecular simulation methods have increasingly contributed to our understanding of molecular and nanoscale systems. However, the family of Monte Carlo techniques has taken a backseat to molecular dynamics based methods, which is also reflected in the number of available simulation packages. Here, we report the development of a generic, versatile simulation package for stochastic simulations and demonstrate its application to protein conformational change, protein-protein association, small-molecule protein docking, and simulation of the growth of nanoscale clusters of organic molecules. Simulation of molecular and nanoscale systems (SIMONA) is easy to use for standard simulations via a graphical user interface and highly parallel both via MPI and the use of graphical processors. It is also extendable to many additional simulations types. Being freely available to academic users, we hope it will enable a large community of researchers in the life- and materials-sciences to use and extend SIMONA in the future. SIMONA is available for download under http://int.kit.edu/nanosim/simona.

Structure and function of ameloblastin as an extracellular matrix protein: adhesion, calcium binding, and CD63 interaction in human and mouse

The functional significance of extracellular matrix proteins in the life of vertebrates is underscored by a high level of sequence variability in tandem with a substantial degree of conservation in terms of cell-cell and cell-matrix adhesion interactions. Many extracellular matrix proteins feature multiple adhesion domains for successful attachment to substrates, such as integrin, CD63, and heparin. Here we have used homology and ab initio modeling algorithms to compare mouse ameloblastin (mAMBN) and human ameloblastin (hABMN) isoforms and to analyze their potential for cell adhesion and interaction with other matrix molecules as well as calcium binding. Sequence comparison between mAMBN and hAMBN revealed a 26-amino-acid deletion in mAMBN, corresponding to a helix-loop-helix frameshift. The human AMBN domain (174Q-201G), homologous to the mAMBN 157E-178I helix-loop-helix region, formed a helix-loop motif with an extended loop, suggesting a higher degree of flexibility of hAMBN compared with mAMBN, as confirmed by molecular dynamics simulation. Heparin-binding domains, CD63-interaction domains, and calcium-binding sites in both hAMBN and mAMBN support the concept of AMBN as an extracellular matrix protein. The high level of conservation between AMBN functional domains related to adhesion and differentiation was remarkable when compared with only 61% amino acid sequence homology.

Cross-species analyses identify the BNIP-2 and Cdc42GAP homology (BCH) domain as a distinct functional subclass of the CRAL_TRIO/Sec14 superfamily

The CRAL_TRIO protein domain, which is unique to the Sec14 protein superfamily, binds to a diverse set of small lipophilic ligands. Similar domains are found in a range of different proteins including neurofibromatosis type-1, a Ras GTPase-activating Protein (RasGAP) and Rho guanine nucleotide exchange factors (RhoGEFs). Proteins containing this structural protein domain exhibit a low sequence similarity and ligand specificity while maintaining an overall characteristic three-dimensional structure. We have previously demonstrated that the BNIP-2 and Cdc42GAP Homology (BCH) protein domain, which shares a low sequence homology with the CRAL_TRIO domain, can serve as a regulatory scaffold that binds to Rho, RhoGEFs and RhoGAPs to control various cell signalling processes. In this work, we investigate 175 BCH domain-containing proteins from a wide range of different organisms. A phylogenetic analysis with ∼100 CRAL_TRIO and similar domains from eight representative species indicates a clear distinction of BCH-containing proteins as a novel subclass within the CRAL_TRIO/Sec14 superfamily. BCH-containing proteins contain a hallmark sequence motif R(R/K)h(R/K)(R/K)NL(R/K)xhhhhHPs (‘h’ is large and hydrophobic residue and ‘s’ is small and weekly polar residue) and can be further subdivided into three unique subtypes associated with BNIP-2-N, macro- and RhoGAP-type protein domains. A previously unknown group of genes encoding ‘BCH-only’ domains is also identified in plants and arthropod species. Based on an analysis of their gene-structure and their protein domain context we hypothesize that BCH domain-containing genes evolved through gene duplication, intron insertions and domain swapping events. Furthermore, we explore the point of divergence between BCH and CRAL-TRIO proteins in relation to their ability to bind small GTPases, GAPs and GEFs and lipid ligands. Our study suggests a need for a more extensive analysis of previously uncharacterized BCH, ‘BCH-like’ and CRAL_TRIO-containing proteins and their significance in regulating signaling events involving small GTPases.

Characterization of C- and N-terminal domains of Aquifex aeolicus MutL endonuclease: N-terminal domain stimulates the endonuclease activity of C-terminal domain in a zinc-dependent manner

DNA MMR (mismatch repair) is an excision repair system that removes mismatched bases generated primarily by failure of the 3'-5' proofreading activity associated with replicative DNA polymerases. MutL proteins homologous to human PMS2 are the endonucleases that introduce the entry point of the excision reaction. Deficiency in PMS2 function is one of the major etiologies of hereditary non-polyposis colorectal cancers in humans. Although recent studies revealed that the CTD (C-terminal domain) of MutL harbours weak endonuclease activity, the regulatory mechanism of this activity remains unknown. In this paper, we characterize in detail the CTD and NTD (N-terminal domain) of aqMutL (Aquifex aeolicus MutL). On the one hand, CTD existed as a dimer in solution and showed weak DNA-binding and Mn2+-dependent endonuclease activities. On the other hand, NTD was monomeric and exhibited a relatively strong DNA-binding activity. It was also clarified that NTD promotes the endonuclease activity of CTD. NTD-mediated activation of CTD was abolished by depletion of the zinc-ion from the reaction mixture or by the substitution of the zinc-binding cysteine residue in CTD with an alanine. On the basis of these results, we propose a model for the intramolecular regulatory mechanism of MutL endonuclease activity.

Molecular cloning, sequence analysis and expression of Fein-Penaeidin from the haemocytes of Indian white shrimp Fenneropenaeus indicus

Penaeidins are members of a special family of antimicrobial peptide existing in penaeid shrimp and play an important role in the immunological defense of shrimp. Here, we report a penaeidin sequence cloned from the Indian white shrimp Fenneropenaus indicus (Fein-Penaeidin). The Fein-Penaeidin open reading frame encodes a 77 amino acid peptide including a 19 amino acid signal peptide. The deduced amino acid sequences of Fein-Penaeidin include a proline rich N-terminal domain and a carboxyl-domain that contains six cysteine residues. Structural analysis revealed an alpha-helix in its secondary structure and the predicted 3D structure indicated two-disulphide bridges in the alpha-helix. Phylogenetic analysis and sequence comparison with other known peaneidin suggest the gene shows high similarity to that of penaeidin from Peneaus monodon (95%), F. indicus (80%) and Fenneropenaeus chinensis (74%). Fein-Penaeidin was examined in normal and microbial challenged shrimp and was found to be constitutively expressed in haemocytes, Heart, gills, muscles, intestine, hepatopancreas and eyestalk. Bacterial challenge resulted in mRNA up-regulation, inducing expression at 6 h post injection indicating the penaeidin involved in the innate immunity.

Prediction of a novel RNA binding domain in crocodilepox Zimbabwe Gene 157

Background

Although the crocodilepox virus (CRV) is currently unclassified, phylogenetic analyses suggest that its closest known relatives are molluscum contagiosum virus (MCV) and the avipox viruses. The CRV genome is approximately 190 kb and contains a large number of unique genes in addition to the set of conserved Chordopoxvirus genes found in all such viruses. Upon sequencing the viral genome, others noted that this virus was also unusual because of the lack of a series of common immuno-suppressive genes. However, the genome contains multiple genes of unknown function that are likely to function in reducing the anti-viral response of the host.

Results

By using sensitive database searches for similarity, we observed that gene 157 of CRV-strain Zimbabwe (CRV-ZWE) encodes a protein with a domain that is predicted to bind dsRNA. Domain characterization supported this prediction, therefore, we tested the ability of the Robetta protein structure prediction server to model the amino acid sequence of this protein on a well-characterized RNA binding domain. The model generated by Robetta suggests that CRV-ZWE-157 does indeed contain an RNA binding domain; the model could be overlaid on the template protein structure with high confidence.

Conclusion

We hypothesize that CRV-ZWE-157 encodes a novel poxvirus RNA binding protein and suggest that as a non-core gene it may play a role in host-range determination or function to dampen host anti-viral responses. Potential targets for this CRV protein include the host interferon response and miRNA pathways.

Predicting residue-residue contacts and helix-helix interactions in transmembrane proteins using an integrative feature-based random forest approach

Integral membrane proteins constitute 25-30% of genomes and play crucial roles in many biological processes. However, less than 1% of membrane protein structures are in the Protein Data Bank. In this context, it is important to develop reliable computational methods for predicting the structures of membrane proteins. Here, we present the first application of random forest (RF) for residue-residue contact prediction in transmembrane proteins, which we term as TMhhcp. Rigorous cross-validation tests indicate that the built RF models provide a more favorable prediction performance compared with two state-of-the-art methods, i.e., TMHcon and MEMPACK. Using a strict leave-one-protein-out jackknifing procedure, they were capable of reaching the top L/5 prediction accuracies of 49.5% and 48.8% for two different residue contact definitions, respectively. The predicted residue contacts were further employed to predict interacting helical pairs and achieved the Matthew's correlation coefficients of 0.430 and 0.424, according to two different residue contact definitions, respectively. To facilitate the academic community, the TMhhcp server has been made freely accessible at http://protein.cau.edu.cn/tmhhcp.

The Proteome Folding Project: proteome-scale prediction of structure and function

The incompleteness of proteome structure and function annotation is a critical problem for biologists and, in particular, severely limits interpretation of high-throughput and next-generation experiments. We have developed a proteome annotation pipeline based on structure prediction, where function and structure annotations are generated using an integration of sequence comparison, fold recognition, and grid-computing-enabled de novo structure prediction. We predict protein domain boundaries and three-dimensional (3D) structures for protein domains from 94 genomes (including human, Arabidopsis, rice, mouse, fly, yeast, Escherichia coli, and worm). De novo structure predictions were distributed on a grid of more than 1.5 million CPUs worldwide (World Community Grid). We generated significant numbers of new confident fold annotations (9% of domains that are otherwise unannotated in these genomes). We demonstrate that predicted structures can be combined with annotations from the Gene Ontology database to predict new and more specific molecular functions.

Four small puzzles that Rosetta doesn't solve

A complete macromolecule modeling package must be able to solve the simplest structure prediction problems. Despite recent successes in high resolution structure modeling and design, the Rosetta software suite fares poorly on small protein and RNA puzzles, some as small as four residues. To illustrate these problems, this manuscript presents Rosetta results for four well-defined test cases: the 20-residue mini-protein Trp cage, an even smaller disulfide-stabilized conotoxin, the reactive loop of a serine protease inhibitor, and a UUCG RNA tetraloop. In contrast to previous Rosetta studies, several lines of evidence indicate that conformational sampling is not the major bottleneck in modeling these small systems. Instead, approximations and omissions in the Rosetta all-atom energy function currently preclude discriminating experimentally observed conformations from de novo models at atomic resolution. These molecular “puzzles” should serve as useful model systems for developers wishing to make foundational improvements to this powerful modeling suite.

Adsorption of homopolypeptides on gold investigated using atomistic molecular dynamics

We investigate the role of dynamics on adsorption of peptides to gold surfaces using all-atom molecular dynamics simulations in explicit solvent. We choose six homopolypeptides [Ala(10), Ser(10), Thr(10), Arg(10), Lys(10), and Gln(10)], for which experimental surface coverages are not correlated with amino acid level affinities for gold, with the idea that dynamic properties may also play a role. To assess dynamics we determine both conformational movement and flexibility of the peptide within a given conformation. Low conformational movement indicates stability of a given conformation and leads to less adsorption than homopolypeptides with faster conformational movement. Likewise, low flexibility within a given conformation also leads to less adsorption. Neither amino acid affinities nor dynamic considerations alone predict surface coverage; rather both quantities must be considered in peptide adsorption to gold surfaces.

Subunit organization and Rab interactions of Vps-C protein complexes that control endolysosomal membrane traffic

The Vps-C complexes, CORVET and HOPS, are key regulators of membrane traffic through late endosomes and lysosomes. In this study Vps-C intersubunit contacts, domain architecture, and interactions with Rab G proteins are systematically dissected using genetic and biochemical approaches.

Traffic through late endolysosomal compartments is regulated by sequential signaling of small G proteins of the Rab5 and Rab7 families. The Saccharomyces cerevisiae Vps-C protein complexes CORVET (class C core vacuole/endosome tethering complex) and HOPS (homotypic fusion and protein transport) interact with endolysosomal Rabs to coordinate their signaling activities. To better understand these large and intricate complexes, we performed interaction surveys to assemble domain-level interaction topologies for the eight Vps-C subunits. We identified numerous intersubunit interactions and up to six Rab-binding sites. Functional modules coordinate the major Rab interactions within CORVET and HOPS. The CORVET-specific subunits, Vps3 and Vps8, form a subcomplex and physically and genetically interact with the Rab5 orthologue Vps21. The HOPS-specific subunits, Vps39 and Vps41, also form a subcomplex. Both subunits bind the Rab7 orthologue Ypt7, but with distinct nucleotide specificities. The in vivo functions of four RING-like domains within Vps-C subunits were analyzed and shown to have distinct functions in endolysosomal transport. Finally, we show that the CORVET- and HOPS-specific subunits Vps3 and Vps39 bind the Vps-C core through a common region within the Vps11 C-terminal domain (CTD). Biochemical and genetic experiments demonstrate the importance of these regions, revealing the Vps11 CTD as a key integrator of Vps-C complex assembly, Rab signaling, and endosomal and lysosomal traffic.

Kindlin-2 regulates podocyte adhesion and fibronectin matrix deposition through interactions with phosphoinositides and integrins

Kindlin-2 is a FERM and PH domain-containing integrin-binding protein that is emerging as an important regulator of integrin activation. How kindlin-2 functions in integrin activation, however, is not known. We report here that kindlin-2 interacts with multiple phosphoinositides, preferentially with phosphatidylinositol 3,4,5-trisphosphate. Although integrin-binding is essential for focal adhesion localization of kindlin-2, phosphoinositide-binding is not required for this process. Using biologically and clinically relevant glomerular podocytes as a model system, we show that integrin activation and dependent processes are tightly regulated by kindlin-2: depletion of kindlin-2 reduced integrin activation, matrix adhesion and fibronectin matrix deposition, whereas overexpression of kindlin-2 promoted these processes. Furthermore, we provide evidence showing that kindlin-2 is involved in phosphoinositide-3-kinase-mediated regulation of podocyte-matrix adhesion and fibronectin matrix deposition. Mechanistically, kindlin-2 promotes integrin activation and integrin-dependent processes through interacting with both integrins and phosphoinositides. TGF-β1, a mediator of progressive glomerular failure, markedly increased the level of kindlin-2 and fibronectin matrix deposition, and the latter process was reversed by depletion of kindlin-2. Our results reveal important functions of kindlin-2 in the regulation of podocyte-matrix adhesion and matrix deposition and shed new light on the mechanism whereby kindlin-2 functions in these processes.

Brainstorming: weighted voting prediction of inhibitors for protein targets

The "Brainstorming" approach presented in this paper is a weighted voting method that can improve the quality of predictions generated by several machine learning (ML) methods. First, an ensemble of heterogeneous ML algorithms is trained on available experimental data, then all solutions are gathered and a consensus is built between them. The final prediction is performed using a voting procedure, whereby the vote of each method is weighted according to a quality coefficient calculated using multivariable linear regression (MLR). The MLR optimization procedure is very fast, therefore no additional computational cost is introduced by using this jury approach. Here, brainstorming is applied to selecting actives from large collections of compounds relating to five diverse biological targets of medicinal interest, namely HIV-reverse transcriptase, cyclooxygenase-2, dihydrofolate reductase, estrogen receptor, and thrombin. The MDL Drug Data Report (MDDR) database was used for selecting known inhibitors for these protein targets, and experimental data was then used to train a set of machine learning methods. The benchmark dataset (available at http://bio.icm.edu.pl/∼darman/chemoinfo/benchmark.tar.gz ) can be used for further testing of various clustering and machine learning methods when predicting the biological activity of compounds. Depending on the protein target, the overall recall value is raised by at least 20% in comparison to any single machine learning method (including ensemble methods like random forest) and unweighted simple majority voting procedures.

Oligomerization of Uukuniemi virus nucleocapsid protein

Background

Uukuniemi virus (UUKV) belongs to the Phlebovirus genus in the family Bunyaviridae. As a non-pathogenic virus for humans UUKV has served as a safe model bunyavirus in a number of studies addressing fundamental questions such as organization and regulation of viral genes, genome replication, structure and assembly. The present study is focused on the oligomerization of the UUKV nucleocapsid (N) protein, which plays an important role in several steps of virus replication. The aim was to locate the domains involved in the N protein oligomerization and study the process in detail.

Results

A set of experiments concentrating on the N- and C-termini of the protein was performed, first by completely or partially deleting putative N-N-interaction domains and then by introducing point mutations of amino acid residues. Mutagenesis strategy was based on the computer modeling of secondary and tertiary structure of the N protein. The N protein mutants were studied in chemical cross-linking, immunofluorescence, mammalian two-hybrid, minigenome, and virus-like particle-forming assays. The data showed that the oligomerization ability of UUKV-N protein depends on the presence of intact α-helices on both termini of the N protein molecule and that a specific structure in the N-terminal region plays a crucial role in the N-N interaction(s). This structure is formed by two α-helices, rich in amino acid residues with aromatic (W7, F10, W19, F27, F31) or long aliphatic (I14, I24) side chains. Furthermore, some of the N-terminal mutations (e.g. I14A, I24A, F31A) affected the N protein functionality both in mammalian two-hybrid and minigenome assays.

Conclusions

UUKV-N protein has ability to form oligomers in chemical cross-linking and mammalian two-hybrid assays. In mutational analysis, some of the introduced single-point mutations abolished the N protein functionality both in mammalian two-hybrid and minigenome assays, suggesting that especially the N-terminal region of the UUKV-N protein is essential for the N-N interaction.

In vitro effects of recombinant otoconin 90 upon calcite crystal growth. Significance of tertiary structure

Otoconia are biomineral particles of microscopic size essential for perception of gravity and maintenance of balance. Millions of older Americans are affected in their mobility, quality of life and in their health by progressive demineralization of otoconia. Currently, no effective means to prevent or counteract this process are available. Because of prohibitive anatomical and biological constraints, otoconial research is lagging far behind other systems such as bone and teeth. We have overcome these obstacles by generating otoconial matrix proteins by recombinant techniques. In the present study, we evaluated the effects of recombinant Otoconin 90 (OC90), the principal soluble matrix protein upon calcite crystal growth patterns in vitro. Our findings highlight multiple effects, including facilitation of nucleation, and inhibition of crystal growth in a concentration-dependent manner. Moreover, OC90 induces morphologic changes characteristic of native otoconia. OC90 is considerably less acidic than the prototypical invertebrate CaCO(3) -associated protein, but is nevertheless an effective modulator of calcite crystal growth. Based on homology modeling of the sPLA2-like domains of OC90, we propose that the lower density of acidic residues of the primary sequence is compensated by formation of major anionic surface clusters upon folding into tertiary conformation.

BCL::contact-low confidence fold recognition hits boost protein contact prediction and de novo structure determination

Knowledge of all residue-residue contacts within a protein allows determination of the protein fold. Accurate prediction of even a subset of long-range contacts (contacts between amino acids far apart in sequence) can be instrumental for determining tertiary structure. Here we present BCL::Contact, a novel contact prediction method that utilizes artificial neural networks (ANNs) and specializes in the prediction of medium to long-range contacts. BCL::Contact comes in two modes: sequence-based and structure-based. The sequence-based mode uses only sequence information and has individual ANNs specialized for helix-helix, helix-strand, strand-helix, strand-strand, and sheet-sheet contacts. The structure-based mode combines results from 32-fold recognition methods with sequence information to a consensus prediction. The two methods were presented in the 6(th) and 7(th) Critical Assessment of Techniques for Protein Structure Prediction (CASP) experiments. The present work focuses on elucidating the impact of fold recognition results onto contact prediction via a direct comparison of both methods on a joined benchmark set of proteins. The sequence-based mode predicted contacts with 42% accuracy (7% false positive rate), while the structure-based mode achieved 45% accuracy (2% false positive rate). Predictions by both modes of BCL::Contact were supplied as input to the protein tertiary structure prediction program Rosetta for a benchmark of 17 proteins with no close sequence homologs in the protein data bank (PDB). Rosetta created higher accuracy models, signified by an improvement of 1.3 A on average root mean square deviation (RMSD), when driven by the predicted contacts. Further, filtering Rosetta models by agreement with the predicted contacts enriches for native-like fold topologies.

RosettaBackrub--a web server for flexible backbone protein structure modeling and design

The RosettaBackrub server (http://kortemmelab.ucsf.edu/backrub) implements the Backrub method, derived from observations of alternative conformations in high-resolution protein crystal structures, for flexible backbone protein modeling. Backrub modeling is applied to three related applications using the Rosetta program for structure prediction and design: (I) modeling of structures of point mutations, (II) generating protein conformational ensembles and designing sequences consistent with these conformations and (III) predicting tolerated sequences at protein–protein interfaces. The three protocols have been validated on experimental data. Starting from a user-provided single input protein structure in PDB format, the server generates near-native conformational ensembles. The predicted conformations and sequences can be used for different applications, such as to guide mutagenesis experiments, for ensemble-docking approaches or to generate sequence libraries for protein design.

Practically useful: what the Rosetta protein modeling suite can do for you

The objective of this review is to enable researchers to use the software package Rosetta for biochemical and biomedicinal studies. We provide a brief review of the six most frequent research problems tackled with Rosetta. For each of these six tasks, we provide a tutorial that illustrates a basic Rosetta protocol. The Rosetta method was originally developed for de novo protein structure prediction and is regularly one of the best performers in the community-wide biennial Critical Assessment of Structure Prediction. Predictions for protein domains with fewer than 125 amino acids regularly have a backbone root-mean-square deviation of better than 5.0 Å. More impressively, there are several cases in which Rosetta has been used to predict structures with atomic level accuracy better than 2.5 Å. In addition to de novo structure prediction, Rosetta also has methods for molecular docking, homology modeling, determining protein structures from sparse experimental NMR or EPR data, and protein design. Rosetta has been used to accurately design a novel protein structure, predict the structure of protein−protein complexes, design altered specificity protein−protein and protein−DNA interactions, and stabilize proteins and protein complexes. Most recently, Rosetta has been used to solve the X-ray crystallographic phase problem.

Directionality in protein fold prediction

BACKGROUND

Ever since the ground-breaking work of Anfinsen et al. in which a denatured protein was found to refold to its native state, it has been frequently stated by the protein fold prediction community that all the information required for protein folding lies in the amino acid sequence. Recent in vitro experiments and in silico computational studies, however, have shown that cotranslation may affect the folding pathway of some proteins, especially those of ancient folds. In this paper aspects of cotranslational folding have been incorporated into a protein structure prediction algorithm by adapting the Rosetta program to fold proteins as the nascent chain elongates. This makes it possible to conduct a pairwise comparison of folding accuracy, by comparing folds created sequentially from each end of the protein.

RESULTS

A single main result emerged: in 94% of proteins analyzed, following the sense of translation, from N-terminus to C-terminus, produced better predictions than following the reverse sense of translation, from the C-terminus to N-terminus. Two secondary results emerged. First, this superiority of N-terminus to C-terminus folding was more marked for proteins showing stronger evidence of cotranslation and second, an algorithm following the sense of translation produced predictions comparable to, and occasionally better than, Rosetta.

CONCLUSIONS

There is a directionality effect in protein fold prediction. At present, prediction methods appear to be too noisy to take advantage of this effect; as techniques refine, it may be possible to draw benefit from a sequential approach to protein fold prediction.

Fast and accurate automatic structure prediction with HHpred

Automated protein structure prediction is becoming a mainstream tool for biological research. This has been fueled by steady improvements of publicly available automated servers over the last decade, in particular their ability to build good homology models for an increasing number of targets by reliably detecting and aligning more and more remotely homologous templates. Here, we describe the three fully automated versions of the HHpred server that participated in the community-wide blind protein structure prediction competition CASP8. What makes HHpred unique is the combination of usability, short response times (typically under 15 min) and a model accuracy that is competitive with those of the best servers in CASP8.

No Nogo66- and NgR-mediated inhibition of regenerating axons in the zebrafish optic nerve

In contrast to mammals, lesioned axons in the zebrafish (ZF) optic nerve regenerate and restore vision. This correlates with the absence of the NogoA-specific N-terminal domains from the ZF nogo/rtn-4 (reticulon-4) gene that inhibits regeneration in mammals. However, mammalian nogo/rtn-4 carries a second inhibitory C-terminal domain, Nogo-66, being 70% identical with ZF-Nogo66. The present study examines, (1) whether ZF-Nogo66 is inhibitory and effecting similar signaling pathways upon Nogo66-binding to the Nogo66 receptor NgR and its coreceptors, and (2) whether Rat-Nogo66 on fish, and ZF-Nogo66 on mouse neurons, cause inhibition via NgR. Our results from "outgrowth, collapse and contact assays" suggest, surprisingly, that ZF-Nogo66 is growth-permissive for ZF and mouse neurons, quite in contrast to its Rat-Nogo66 homolog which inhibits growth. The opposite effects of ZF- and Rat-Nogo66 are, in both fish and mouse, transmitted by GPI (glycosylphosphatidylinositol)-anchored receptors, including NgR. The high degree of sequence homology in the predicted binding site is consistent with the ability of ZF- and mammalian-Nogo66 to bind to NgRs of both species. Yet, Rat-Nogo66 elicits phosphorylation of the downstream effector cofilin whereas ZF-Nogo66 has no influence on cofilin phosphorylation--probably because of significantly different Rat- versus ZF-Nogo66 sequences outside of the receptor-binding region effecting, by speculation, recruitment of a different set of coreceptors or microdomain association of NgR. Thus, not only was the NogoA-specific domain lost in fish, but Nogo66, the second inhibitory domain in mammals, and its signaling upon binding to NgR, was modified so that ZF-Nogo/RTN-4 does not impair axon regeneration.

Phylogenesis of brain-derived neurotrophic factor (BDNF) in vertebrates

Brain-derived neurotrophic factor (BDNF) belongs to neurotrophin family, a class of molecules playing key roles in neuronal development, survival and regeneration, neurite growth and plasticity: memory processes are mainly affected, and mutations of the human BDNF gene are associated to cognitive and behavioural disturbances. All neurotrophins contain a highly conserved C-terminal domain and bind to the same receptor family. Both correct folding and post-translational processing of the entire preproprotein are pivotal for sorting to the extracellular space, dimerization and receptor binding. Evolutionary studies conducted so far demonstrate that a single ancestor gene underwent two independent duplication events at an early stage of vertebrate evolution, leading to the formation of the current neurotrophins. However, works focusing on BDNF evolution are scarce and fragmentary, mainly in lower vertebrates. In this work, we report cloning of eight DNA sequences from amphibians and teleosts, and analysis of the entire coding regions (cDNA sequences) of BDNF from 35 organisms, from teleosts to mammals. A phylogenetic tree was constructed and the analysis of non-synonymous-synonymous substitution rates performed for the different branches. Our results suggest that natural selection is acting on mammals, separating them from other classes. Since preproprotein cleavage and 3D structure of mature protein are important for functional activity of BDNF, we also propose a de novo prediction of the 3D structure of translates in at least one species for each class, in order to get hints about the functional constraints of the protein.

Borrelia burgdorferi locus BB0795 encodes a BamA orthologue required for growth and efficient localization of outer membrane proteins

The outer membrane (OM) of the pathogenic diderm spirochete, Borrelia burgdorferi, contains integral beta-barrel outer membrane proteins (OMPs) in addition to its numerous outer surface lipoproteins. Very few OMPs have been identified in B. burgdorferi, and the protein machinery required for OMP assembly and OM localization is currently unknown. Essential OM BamA proteins have recently been characterized in Gram-negative bacteria that are central components of an OM beta-barrel assembly machine and are required for proper localization and insertion of bacterial OMPs. In the present study, we characterized a putative B. burgdorferi BamA orthologue encoded by open reading frame bb0795. Structural model predictions and cellular localization data indicate that the B. burgdorferi BB0795 protein contains an N-terminal periplasmic domain and a C-terminal, surface-exposed beta-barrel domain. Additionally, assays with an IPTG-regulatable bb0795 mutant revealed that BB0795 is required for B. burgdorferi growth. Furthermore, depletion of BB0795 results in decreased amounts of detectable OMPs in the B. burgdorferi OM. Interestingly, a decrease in the levels of surface-exposed lipoproteins was also observed in the mutant OMs. Collectively, our structural, cellular localization and functional data are consistent with the characteristics of other BamA proteins, indicating that BB0795 is a B. burgdorferi BamA orthologue.

(PS)2-v2: template-based protein structure prediction server

Background

Template selection and target-template alignment are critical steps for template-based modeling (TBM) methods. To identify the template for the twilight zone of 15~25% sequence similarity between targets and templates is still difficulty for template-based protein structure prediction. This study presents the (PS)²-v2 server, based on our original server with numerous enhancements and modifications, to improve reliability and applicability.

Results

To detect homologous proteins with remote similarity, the (PS)²-v2 server utilizes the S2A2 matrix, which is a 60 × 60 substitution matrix using the secondary structure propensities of 20 amino acids, and the position-specific sequence profile (PSSM) generated by PSI-BLAST. In addition, our server uses multiple templates and multiple models to build and assess models. Our method was evaluated on the Lindahl benchmark for fold recognition and ProSup benchmark for sequence alignment. Evaluation results indicated that our method outperforms sequence-profile approaches, and had comparable performance to that of structure-based methods on these benchmarks. Finally, we tested our method using the 154 TBM targets of the CASP8 (Critical Assessment of Techniques for Protein Structure Prediction) dataset. Experimental results show that (PS)²-v2 is ranked 6^th among 72 severs and is faster than the top-rank five serves, which utilize ab initio methods.

Conclusion

Experimental results demonstrate that (PS)²-v2 with the S2A2 matrix is useful for template selections and target-template alignments by blending the amino acid and structural propensities. The multiple-template and multiple-model strategies are able to significantly improve the accuracies for target-template alignments in the twilight zone. We believe that this server is useful in structure prediction and modeling, especially in detecting homologous templates with sequence similarity in the twilight zone.

Evaluating the absolute quality of a single protein model using structural features and support vector machines

Knowing the quality of a protein structure model is important for its appropriate usage. We developed a model evaluation method to assess the absolute quality of a single protein model using only structural features with support vector machine regression. The method assigns an absolute quantitative score (i.e. GDT-TS) to a model by comparing its secondary structure, relative solvent accessibility, contact map, and beta sheet structure with their counterparts predicted from its primary sequence. We trained and tested the method on the CASP6 dataset using cross-validation. The correlation between predicted and true scores is 0.82. On the independent CASP7 dataset, the correlation averaged over 95 protein targets is 0.76; the average correlation for template-based and ab initio targets is 0.82 and 0.50, respectively. Furthermore, the predicted absolute quality scores can be used to rank models effectively. The average difference (or loss) between the scores of the top-ranked models and the best models is 5.70 on the CASP7 targets. This method performs favorably when compared with the other methods used on the same dataset. Moreover, the predicted absolute quality scores are comparable across models for different proteins. These features make the method a valuable tool for model quality assurance and ranking.

Molecular biology, genetics and biochemistry of the repulsive guidance molecule family

RGMs (repulsive guidance molecules) comprise a recently discovered family of GPI (glycosylphosphatidylinositol)-linked cell-membrane-associated proteins found in most vertebrate species. The three proteins, RGMa, RGMb and RGMc, products of distinct single-copy genes that arose early in vertebrate evolution, are approximately 40-50% identical to each other in primary amino acid sequence, and share similarities in predicted protein domains and overall structure, as inferred by ab initio molecular modelling; yet the respective proteins appear to undergo distinct biosynthetic and processing steps, whose regulation has not been characterized to date. Each RGM also displays a discrete tissue-specific pattern of gene and protein expression, and each is proposed to have unique biological functions, ranging from axonal guidance during development (RGMa) to regulation of systemic iron metabolism (RGMc). All three RGM proteins appear capable of binding selected BMPs (bone morphogenetic proteins), and interactions with BMPs mediate at least some of the biological effects of RGMc on iron metabolism, but to date no role for BMPs has been defined in the actions of RGMa or RGMb. RGMa and RGMc have been shown to bind to the transmembrane protein neogenin, which acts as a critical receptor to mediate the biological effects of RGMa on repulsive axonal guidance and on neuronal survival, but its role in the actions of RGMc remains to be elucidated. Similarly, the full spectrum of biological functions of the three RGMs has not been completely characterized yet, and will remain an active topic of ongoing investigation.