Evaluation of 3D-Jury on CASP7 models

MQAPRank: improved global protein model quality assessment by learning-to-rank

BACKGROUND

Protein structure prediction has achieved a lot of progress during the last few decades and a greater number of models for a certain sequence can be predicted. Consequently, assessing the qualities of predicted protein models in perspective is one of the key components of successful protein structure prediction. Over the past years, a number of methods have been developed to address this issue, which could be roughly divided into three categories: single methods, quasi-single methods and clustering (or consensus) methods. Although these methods achieve much success at different levels, accurate protein model quality assessment is still an open problem.

RESULTS

Here, we present the MQAPRank, a global protein model quality assessment program based on learning-to-rank. The MQAPRank first sorts the decoy models by using single method based on learning-to-rank algorithm to indicate their relative qualities for the target protein. And then it takes the first five models as references to predict the qualities of other models by using average GDT_TS scores between reference models and other models. Benchmarked on CASP11 and 3DRobot datasets, the MQAPRank achieved better performances than other leading protein model quality assessment methods. Recently, the MQAPRank participated in the CASP12 under the group name FDUBio and achieved the state-of-the-art performances.

CONCLUSIONS

The MQAPRank provides a convenient and powerful tool for protein model quality assessment with the state-of-the-art performances, it is useful for protein structure prediction and model quality assessment usages.

Sorting protein decoys by machine-learning-to-rank

Much progress has been made in Protein structure prediction during the last few decades. As the predicted models can span a broad range of accuracy spectrum, the accuracy of quality estimation becomes one of the key elements of successful protein structure prediction. Over the past years, a number of methods have been developed to address this issue, and these methods could be roughly divided into three categories: the single-model methods, clustering-based methods and quasi single-model methods. In this study, we develop a single-model method MQAPRank based on the learning-to-rank algorithm firstly, and then implement a quasi single-model method Quasi-MQAPRank. The proposed methods are benchmarked on the 3DRobot and CASP11 dataset. The five-fold cross-validation on the 3DRobot dataset shows the proposed single model method outperforms other methods whose outputs are taken as features of the proposed method, and the quasi single-model method can further enhance the performance. On the CASP11 dataset, the proposed methods also perform well compared with other leading methods in corresponding categories. In particular, the Quasi-MQAPRank method achieves a considerable performance on the CASP11 Best150 dataset.

Structural models of mammalian mitochondrial transcription factor B2

Mammalian mitochondrial DNA (mtDNA) encodes 13 core proteins of oxidative phosphorylation, 12S and 16S ribosomal RNAs, and 22 transfer RNAs. Mutations and deletions of mtDNA and/or nuclear genes encoding mitochondrial proteins have been implicated in a wide range of diseases. Thus, cell survival and health of the organism require some steady-state level of the mitochondrial genome and its expression. In mammalian systems, the mitochondrial transcription factor B2 (mtTFB2 or TFB2M) is indispensable for transcription initiation. TFB2M along with two other proteins, mitochondrial RNA polymerase (mtRNAP or POLRMT) and mitochondrial transcription factor A (mtTFA or TFAM), are key components of the core mitochondrial transcription apparatus. Structural information for POLRMT and TFAM from humans is available; however, there is no available structure for TFB2M. In the present study, three-dimensional structure of TFB2M from humans was modeled using a combination of homology modeling and small-angle X-ray scattering (SAXS). The TFB2M structural model adds substantively to our understanding of TFB2M function. An explanation for the low or absent RNA methyltransferase activity is provided. A putative nucleic acid-binding site is revealed. The amino and carboxy termini, while likely lacking defined secondary structure, appear to adopt compact, globular conformations, thus "capping" the ends of the protein. Finally, sites of interaction of TFB2M with other factors, protein and/or nucleic acid, are suggested by the identification of species-specific clusters on the surface of the protein.

Designing and evaluating the MULTICOM protein local and global model quality prediction methods in the CASP10 experiment

BACKGROUND

Protein model quality assessment is an essential component of generating and using protein structural models. During the Tenth Critical Assessment of Techniques for Protein Structure Prediction (CASP10), we developed and tested four automated methods (MULTICOM-REFINE, MULTICOM-CLUSTER, MULTICOM-NOVEL, and MULTICOM-CONSTRUCT) that predicted both local and global quality of protein structural models.

RESULTS

MULTICOM-REFINE was a clustering approach that used the average pairwise structural similarity between models to measure the global quality and the average Euclidean distance between a model and several top ranked models to measure the local quality. MULTICOM-CLUSTER and MULTICOM-NOVEL were two new support vector machine-based methods of predicting both the local and global quality of a single protein model. MULTICOM-CONSTRUCT was a new weighted pairwise model comparison (clustering) method that used the weighted average similarity between models in a pool to measure the global model quality. Our experiments showed that the pairwise model assessment methods worked better when a large portion of models in the pool were of good quality, whereas single-model quality assessment methods performed better on some hard targets when only a small portion of models in the pool were of reasonable quality.

CONCLUSIONS

Since digging out a few good models from a large pool of low-quality models is a major challenge in protein structure prediction, single model quality assessment methods appear to be poised to make important contributions to protein structure modeling. The other interesting finding was that single-model quality assessment scores could be used to weight the models by the consensus pairwise model comparison method to improve its accuracy.

Dynactin helps target Polo-like kinase 1 to kinetochores via its left-handed beta-helical p27 subunit

Dynactin is a protein complex required for the in vivo function of cytoplasmic dynein, a microtubule (MT)-based motor. Dynactin binds both dynein and MTs via its p150(Glued) subunit, but little is known about the 'pointed-end complex' that includes the protein subunits Arp11, p62 and the p27/p25 heterodimer. Here, we show that the p27/p25 heterodimer undergoes mitotic phosphorylation by cyclin-dependent kinase 1 (Cdk1) at a single site, p27 Thr186, to generate an anchoring site for polo-like kinase 1 (Plk1) at kinetochores. Removal of p27/p25 from dynactin results in reduced levels of Plk1 and its phosphorylated substrates at kinetochores in prometaphase, which correlates with aberrant kinetochore-MT interactions, improper chromosome alignment and abbreviated mitosis. To investigate the structural implications of p27 phosphorylation, we determined the structure of human p27. This revealed an unusual left-handed β-helix domain, with the phosphorylation site located within a disordered, C-terminal segment. We conclude that dynactin plays a previously undescribed regulatory role in the spindle assembly checkpoint by recruiting Plk1 to kinetochores and facilitating phosphorylation of important downstream targets.

An assembly model of rift valley Fever virus

Rift Valley fever virus (RVFV) is a bunyavirus endemic to Africa and the Arabian Peninsula that infects humans and livestock. The virus encodes two glycoproteins, Gn and Gc, which represent the major structural antigens and are responsible for host cell receptor binding and fusion. Both glycoproteins are organized on the virus surface as cylindrical hollow spikes that cluster into distinct capsomers with the overall assembly exhibiting an icosahedral symmetry. Currently, no experimental three-dimensional structure for any entire bunyavirus glycoprotein is available. Using fold recognition, we generated molecular models for both RVFV glycoproteins and found significant structural matches between the RVFV Gn protein and the influenza virus hemagglutinin protein and a separate match between RVFV Gc protein and Sindbis virus envelope protein E1. Using these models, the potential interaction and arrangement of both glycoproteins in the RVFV particle was analyzed, by modeling their placement within the cryo-electron microscopy density map of RVFV. We identified four possible arrangements of the glycoproteins in the virion envelope. Each assembly model proposes that the ectodomain of Gn forms the majority of the protruding capsomer and that Gc is involved in formation of the capsomer base. Furthermore, Gc is suggested to facilitate intercapsomer connections. The proposed arrangement of the two glycoproteins on the RVFV surface is similar to that described for the alphavirus E1-E2 proteins. Our models will provide guidance to better understand the assembly process of phleboviruses and such structural studies can also contribute to the design of targeted antivirals.

Tri-domain bifunctional inhibitor of metallocarboxypeptidases A and serine proteases isolated from marine annelid Sabellastarte magnifica

This study describes a novel bifunctional metallocarboxypeptidase and serine protease inhibitor (SmCI) isolated from the tentacle crown of the annelid Sabellastarte magnifica. SmCI is a 165-residue glycoprotein with a molecular mass of 19.69 kDa (mass spectrometry) and 18 cysteine residues forming nine disulfide bonds. Its cDNA was cloned and sequenced by RT-PCR and nested PCR using degenerated oligonucleotides. Employing this information along with data derived from automatic Edman degradation of peptide fragments, the SmCI sequence was fully characterized, indicating the presence of three bovine pancreatic trypsin inhibitor/Kunitz domains and its high homology with other Kunitz serine protease inhibitors. Enzyme kinetics and structural analyses revealed SmCI to be an inhibitor of human and bovine pancreatic metallocarboxypeptidases of the A-type (but not B-type), with nanomolar K(i) values. SmCI is also capable of inhibiting bovine pancreatic trypsin, chymotrypsin, and porcine pancreatic elastase in varying measures. When the inhibitor and its nonglycosylated form (SmCI N23A mutant) were overproduced recombinantly in a Pichia pastoris system, they displayed the dual inhibitory properties of the natural form. Similarly, two bi-domain forms of the inhibitor (recombinant rSmCI D1-D2 and rSmCI D2-D3) as well as its C-terminal domain (rSmCI-D3) were also overproduced. Of these fragments, only the rSmCI D1-D2 bi-domain retained inhibition of metallocarboxypeptidase A but only partially, indicating that the whole tri-domain structure is required for such capability in full. SmCI is the first proteinaceous inhibitor of metallocarboxypeptidases able to act as well on another mechanistic class of proteases (serine-type) and is the first of this kind identified in nature.

Pleiotropic roles of a ribosomal protein in Dictyostelium discoideum

The cell cycle phase at starvation influences post-starvation differentiation and morphogenesis in Dictyostelium discoideum. We found that when expressed in Saccharomyces cerevisiae, a D. discoideum cDNA that encodes the ribosomal protein S4 (DdS4) rescues mutations in the cell cycle genes cdc24, cdc42 and bem1. The products of these genes affect morphogenesis in yeast via a coordinated moulding of the cytoskeleton during bud site selection. D. discoideum cells that over- or under-expressed DdS4 did not show detectable changes in protein synthesis but displayed similar developmental aberrations whose intensity was graded with the extent of over- or under-expression. This suggested that DdS4 might influence morphogenesis via a stoichiometric effect – specifically, by taking part in a multimeric complex similar to the one involving Cdc24p, Cdc42p and Bem1p in yeast. In support of the hypothesis, the S. cerevisiae proteins Cdc24p, Cdc42p and Bem1p as well as their D. discoideum cognates could be co-precipitated with antibodies to DdS4. Computational analysis and mutational studies explained these findings: a C-terminal domain of DdS4 is the functional equivalent of an SH3 domain in the yeast scaffold protein Bem1p that is central to constructing the bud site selection complex. Thus in addition to being part of the ribosome, DdS4 has a second function, also as part of a multi-protein complex. We speculate that the existence of the second role can act as a safeguard against perturbations to ribosome function caused by spontaneous variations in DdS4 levels.

Structural features and chaperone activity of the NudC protein family

The NudC family consists of four conserved proteins with representatives in all eukaryotes. The archetypal nudC gene from Aspergillus nidulans is a member of the nud gene family that is involved in the maintenance of nuclear migration. This family also includes nudF, whose human orthologue, Lis1, codes for a protein essential for brain cortex development. Three paralogues of NudC are known in vertebrates: NudC, NudC-like (NudCL), and NudC-like 2 (NudCL2). The fourth distantly related member of the family, CML66, contains a NudC-like domain. The three principal NudC proteins have no catalytic activity but appear to play as yet poorly defined roles in proliferating and dividing cells. We present crystallographic and NMR studies of the human NudC protein and discuss the results in the context of structures recently deposited by structural genomics centers (i.e., NudCL and mouse NudCL2). All proteins share the same core CS domain characteristic of proteins acting either as cochaperones of Hsp90 or as independent small heat shock proteins. However, while NudC and NudCL dimerize via an N-terminally located coiled coil, the smaller NudCL2 lacks this motif and instead dimerizes as a result of unique domain swapping. We show that NudC and NudCL, but not NudCL2, inhibit the aggregation of several target proteins, consistent with an Hsp90-independent heat shock protein function. Importantly, and in contrast to several previous reports, none of the three proteins is able to form binary complexes with Lis1. The availability of structural information will be of help in further studies on the cellular functions of the NudC family.

La-related protein 4 binds poly(A), interacts with the poly(A)-binding protein MLLE domain via a variant PAM2w motif, and can promote mRNA stability

The conserved RNA binding protein La recognizes UUU-3'OH on its small nuclear RNA ligands and stabilizes them against 3'-end-mediated decay. We report that newly described La-related protein 4 (LARP4) is a factor that can bind poly(A) RNA and interact with poly(A) binding protein (PABP). Yeast two-hybrid analysis and reciprocal immunoprecipitations (IPs) from HeLa cells revealed that LARP4 interacts with RACK1, a 40S ribosome- and mRNA-associated protein. LARP4 cosediments with 40S ribosome subunits and polyribosomes, and its knockdown decreases translation. Mutagenesis of the RNA binding or PABP interaction motifs decrease LARP4 association with polysomes. Several translation and mRNA metabolism-related proteins use a PAM2 sequence containing a critical invariant phenylalanine to make direct contact with the MLLE domain of PABP, and their competition for the MLLE is thought to regulate mRNA homeostasis. Unlike all ∼150 previously analyzed PAM2 sequences, LARP4 contains a variant PAM2 (PAM2w) with tryptophan in place of the phenylalanine. Binding and nuclear magnetic resonance (NMR) studies have shown that a peptide representing LARP4 PAM2w interacts with the MLLE of PABP within the affinity range measured for other PAM2 motif peptides. A cocrystal of PABC bound to LARP4 PAM2w shows tryptophan in the pocket in PABC-MLLE otherwise occupied by phenylalanine. We present evidence that LARP4 expression stimulates luciferase reporter activity by promoting mRNA stability, as shown by mRNA decay analysis of luciferase and cellular mRNAs. We propose that LARP4 activity is integrated with other PAM2 protein activities by PABP as part of mRNA homeostasis.

Fast geometric consensus approach for protein model quality assessment

Model quality assessment (MQA) is an integral part of protein structure prediction methods that typically generate multiple candidate models. The challenge lies in ranking and selecting the best models using a variety of physical, knowledge-based, and geometric consensus (GC)-based scoring functions. In particular, 3D-Jury and related GC methods assume that well-predicted (sub-)structures are more likely to occur frequently in a population of candidate models, compared to incorrectly folded fragments. While this approach is very successful in the context of diversified sets of models, identifying similar substructures is computationally expensive since all pairs of models need to be superimposed using MaxSub or related heuristics for structure-to-structure alignment. Here, we consider a fast alternative, in which structural similarity is assessed using 1D profiles, e.g., consisting of relative solvent accessibilities and secondary structures of equivalent amino acid residues in the respective models. We show that the new approach, dubbed 1D-Jury, allows to implicitly compare and rank N models in O(N) time, as opposed to quadratic complexity of 3D-Jury and related clustering-based methods. In addition, 1D-Jury avoids computationally expensive 3D superposition of pairs of models. At the same time, structural similarity scores based on 1D profiles are shown to correlate strongly with those obtained using MaxSub. In terms of the ability to select the best models as top candidates 1D-Jury performs on par with other GC methods. Other potential applications of the new approach, including fast clustering of large numbers of intermediate structures generated by folding simulations, are discussed as well.

Similarity of molecular phenotype between known epilepsy gene LGI1 and disease candidate gene LGI2

Background

The LGI2 (leucine-rich, glioma inactivated 2) gene, a prime candidate for partial epilepsy with pericentral spikes, belongs to a family encoding secreted, beta-propeller domain proteins with EPTP/EAR epilepsy-associated repeats. In another family member, LGI1 (leucine-rich, glioma inactivated 1) mutations are responsible for autosomal dominant lateral temporal epilepsy (ADLTE). Because a few LGI1 disease mutations described in the literature cause secretion failure, we experimentally analyzed the secretion efficiency and subcellular localization of several LGI1 and LGI2 mutant proteins corresponding to observed non-synonymous single nucleotide polymorphisms (nsSNPs) affecting the signal peptide, the leucine-rich repeats and the EAR propeller.

Results

Mapping of disease-causing mutations in the EAR domain region onto a 3D-structure model shows that many of these mutations co-localize at an evolutionary conserved surface region of the propeller. We find that wild-type LGI2 is secreted to the extracellular medium in glycosylated form similarly to LGI1, whereas several mutant proteins tested in this study are secretion-deficient and accumulate in the endoplasmic reticulum. Interestingly, mutations at structurally homologous positions in the EAR domain have the same effect on secretion in LGI1 and LGI2.

Conclusions

This similarity of experimental mislocalization phenotypes for mutations at homologous positions of LGI2 and the established epilepsy gene LGI1 suggests that both genes share a potentially common molecular pathogenesis mechanism that might be the reason for genotypically distinct but phenotypically related forms of epilepsy.

A software tool to accelerate design of protein constructs for recombinant expression

Structural and biochemical analysis of proteins requires access to purified protein material. Modern molecular biology technologies facilitate straightforward molecular cloning and expression analysis of multiple protein constructs in parallel, and such approaches have proven very efficient to identify samples suitable for further analysis. A variety of information can be used to support rational design of protein constructs. This includes, e.g. prediction of secondary structure elements, protein domain predictions, and structure prediction methods such as threading. To fully access the available information, collation of data extracted from several different sources is required. This can be cumbersome and sometimes also confusing due to for example different implementation of amino acid residue numbering schemes. The SGC Domain Boundary Analyser tool provides a graphical interface that simplifies and accelerates rational design of protein expression constructs.

Specific coiled-coil interactions contribute to a global model of the structure of the spindle pole body

As the microtubule-organizing center of yeast, the spindle pole body (SPB) is essential for cell viability. Structural studies of the SPB are limited by its low copy number in the cell, its large size and heterogeneous composition, and its association with the nuclear membrane. However, low-resolution or indirect structural information about the SPB may be deciphered through a variety of techniques. Interestingly, a large proportion of SPB proteins are predicted to contain one or more coiled coils, a common protein interaction motif. The high frequency of coiled coils suggests that this structure is important for establishing the overall architecture of the complex. Support for this hypothesis was reported previously for coiled coils from some SPB proteins. Here, we extend this approach of isolating and characterizing additional SPB coiled coils to improve our understanding of SPB structure and organization. Self-associating coiled coils from Bbp1, Mps2, and Nbp1 were observed to form stable parallel homodimers in solution. Coiled-coil peptides from Bbp1 and Mps2 were also observed to hetero-associate. Experimental coiled-coil interaction data from this work and previous studies, as well as predicted and experimental structures for other SPB protein fragments and domains, were integrated to generate a model of the SPB structure.

Innate immune detection of the type III secretion apparatus through the NLRC4 inflammasome

The mammalian innate immune system uses Toll-like receptors (TLRs) and Nod-LRRs (NLRs) to detect microbial components during infection. Often these molecules work in concert; for example, the TLRs can stimulate the production of the proforms of the cytokines IL-1beta and IL-18, whereas certain NLRs trigger their subsequent proteolytic processing via caspase 1. Gram-negative bacteria use type III secretion systems (T3SS) to deliver virulence factors to the cytosol of host cells, where they modulate cell physiology to favor the pathogen. We show here that NLRC4/Ipaf detects the basal body rod component of the T3SS apparatus (rod protein) from S. typhimurium (PrgJ), Burkholderia pseudomallei (BsaK), Escherichia coli (EprJ and EscI), Shigella flexneri (MxiI), and Pseudomonas aeruginosa (PscI). These rod proteins share a sequence motif that is essential for detection by NLRC4; a similar motif is found in flagellin that is also detected by NLRC4. S. typhimurium has two T3SS: Salmonella pathogenicity island-1 (SPI1), which encodes the rod protein PrgJ, and SPI2, which encodes the rod protein SsaI. Although PrgJ is detected by NLRC4, SsaI is not, and this evasion is required for virulence in mice. The detection of a conserved component of the T3SS apparatus enables innate immune responses to virulent bacteria through a single pathway, a strategy that is divergent from that used by plants in which multiple NB-LRR proteins are used to detect T3SS effectors or their effects on cells. Furthermore, the specific detection of the virulence machinery permits the discrimination between pathogenic and nonpathogenic bacteria.

An integrated methodology for mining promiscuous proteins: a case study of an integrative bioinformatics approach for hepatitis C virus non-structural 5A protein

A methodology for elucidation of structural, functional, and mechanistic knowledge on promiscuous proteins is proposed that constitutes a workflow of integrated bioinformatics analysis. Sequence alignments with closely related homologues can reveal conserved regions which are functionally important. Scanning protein motif databases, along with secondary and surface accessibility predictions integrated with post-translational modification sites (PTMs) prediction reveal functional and protein-binding motifs. Integrating this information about the protein with the GO, SCOP, and CATH annotations of the templates can help to formulate a 3D model with reasonable accuracy even in the case of distant sequence homology. A novel integrative model of the non-structural protein 5A of Hepatitis C virus: a hub promiscuous protein with roles in virus replication and host interactions is proposed. The 3D structure for domain II was predicted based on, the Homo sapiens Replication factor-A protein-1 (RPA1), as a template using consensus meta-servers results. Domain III is an intrinsically unstructured domain with a fold from the retroviral matrix protein, which conducts diverse protein interactions and is involved in viral replication and protein interactions. It also has a single-stranded DNA-binding protein motif (SSDP) signature for pyrimidine binding during viral replication. Two protein-binding motifs with high sequence conservation and disordered regions are proposed; the first corresponds to an Interleukin-8B receptor signature (IL-8R-B), while the second has a lymphotoxin beta receptor (LTβR) high local similarity. A mechanism is proposed to their contribution to NS5A Interferon signaling pathway interception. Lastly, the overlapping between LTβR and SSDP is considered as a sign for NS5A date hubs.

Natural genetic engineering of hepatitis C virus NS5A for immune system counterattack

The Hepatitis C virus nonstructural 5A (NS5A) protein is a hydrophilic phosphoprotein with diverse functions. The domain assignment of NS5A had been refined using a systematic in silico bioinformatics approach using DOMAC, the protein is divided into three domains and domain III is subdivided into two subdomains using ProDom and SSEP servers. The fold structure for domains II and III were predicted using the meta-server 3D-Jury. Scanning motif databases (SMART, BLOCKS, and PROSITE) gave new motifs. Two important motifs, the interleukins 1 and 8 interaction motifs, relating to NS5A function in inducing the interleukin 8 promoter, were discovered from the BLOCKS scan. Protein-protein interaction motifs were predicted as hot loops and disordered regions, corresponding to binding regions with the ds-protein kinase R, viral polymerase, and Src homology 3 signaling proteins binding motif. Other hot loops were predicted in the V3 region and in the single-stranded DNA-binding protein motif. The different mechanisms by which the NS5A protein leads to immune system signaling dysfunction points to the natural genetic engineering of this protein.

PreDisorder: ab initio sequence-based prediction of protein disordered regions

Background

Disordered regions are segments of the protein chain which do not adopt stable structures. Such segments are often of interest because they have a close relationship with protein expression and functionality. As such, protein disorder prediction is important for protein structure prediction, structure determination and function annotation.

Results

This paper presents our protein disorder prediction server, PreDisorder. It is based on our ab initio prediction method (MULTICOM-CMFR) which, along with our meta (or consensus) prediction method (MULTICOM), was recently ranked among the top disorder predictors in the eighth edition of the Critical Assessment of Techniques for Protein Structure Prediction (CASP8). We systematically benchmarked PreDisorder along with 26 other protein disorder predictors on the CASP8 data set and assessed its accuracy using a number of measures. The results show that it compared favourably with other ab initio methods and its performance is comparable to that of the best meta and clustering methods.

Conclusion

PreDisorder is a fast and reliable server which can be used to predict protein disordered regions on genomic scale. It is available at http://casp.rnet.missouri.edu/predisorder.html.

Rapid model quality assessment for protein structure predictions using the comparison of multiple models without structural alignments

MOTIVATION

The accurate prediction of the quality of 3D models is a key component of successful protein tertiary structure prediction methods. Currently, clustering- or consensus-based Model Quality Assessment Programs (MQAPs) are the most accurate methods for predicting 3D model quality; however, they are often CPU intensive as they carry out multiple structural alignments in order to compare numerous models. In this study, we describe ModFOLDclustQ--a novel MQAP that compares 3D models of proteins without the need for CPU intensive structural alignments by utilizing the Q measure for model comparisons. The ModFOLDclustQ method is benchmarked against the top established methods in terms of both accuracy and speed. In addition, the ModFOLDclustQ scores are combined with those from our older ModFOLDclust method to form a new method, ModFOLDclust2, that aims to provide increased prediction accuracy with negligible computational overhead.

RESULTS

The ModFOLDclustQ method is competitive with leading clustering-based MQAPs for the prediction of global model quality, yet it is up to 150 times faster than the previous version of the ModFOLDclust method at comparing models of small proteins (<60 residues) and over five times faster at comparing models of large proteins (>800 residues). Furthermore, a significant improvement in accuracy can be gained over the previous clustering-based MQAPs by combining the scores from ModFOLDclustQ and ModFOLDclust to form the new ModFOLDclust2 method, with little impact on the overall time taken for each prediction.

AVAILABILITY

The ModFOLDclustQ and ModFOLDclust2 methods are available to download from http://www.reading.ac.uk/bioinf/downloads/.

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.

Fido, a novel AMPylation domain common to fic, doc, and AvrB

BACKGROUND

The Vibrio parahaemolyticus type III secreted effector VopS contains a fic domain that covalently modifies Rho GTPase threonine with AMP to inhibit downstream signaling events in host cells. The VopS fic domain includes a conserved sequence motif (HPFx[D/E]GN[G/K]R) that contributes to AMPylation. Fic domains are found in a variety of species, including bacteria, a few archaea, and metazoan eukaryotes.

METHODOLOGY/PRINCIPAL FINDINGS

We show that the AMPylation activity extends to a eukaryotic fic domain in Drosophila melanogaster CG9523, and use sequence and structure based computational methods to identify related domains in doc toxins and the type III effector AvrB. The conserved sequence motif that contributes to AMPylation unites fic with doc. Although AvrB lacks this motif, its structure reveals a similar topology to the fic and doc folds. AvrB binds to a peptide fragment of its host virulence target in a similar manner as fic binds peptide substrate. AvrB also orients a phosphate group from a bound ADP ligand near the peptide-binding site and in a similar position as a bound fic phosphate.

CONCLUSIONS/SIGNIFICANCE

The demonstrated eukaryotic fic domain AMPylation activity suggests that the VopS effector has exploited a novel host posttranslational modification. Fic domain-related structures give insight to the AMPylation active site and to the VopS fic domain interaction with its host GTPase target. These results suggest that fic, doc, and AvrB stem from a common ancestor that has evolved to AMPylate protein substrates.

Chaperone therapy for neuronopathic lysosomal diseases: competitive inhibitors as chemical chaperones for enhancement of mutant enzyme activities

Chaperone therapy is a newly developed molecular approach to lysosomal diseases, a group of human genetic diseases causing severe brain damage. We found two valienamine derivatives, N-octyl-4-epi-β-valienamine (NOEV) and N-octyl-β-valienamine (NOV), as promising therapeutic agents for human β-galactosidase deficiency disorders (mainly G_M1-gangliosidosis) and β-glucosidase deficiency disorders (Gaucher disease), respectively. We briefly reviewed the historical background of research in carbasugar glycosidase inhibitors. Originally NOEV and NOV had been discovered as competitive inhibitors, and then their paradoxical bioactivities as chaperones were confirmed in cultured fibroblasts from patients with these disorders. Subsequently G_M1-gangliosidosis model mice were developed and useful for experimental studies. Orally administered NOEV entered the brain through the blood-brain barrier, enhanced β-galactosidase activity, reduced substrate storage, and improved neurological deterioration clinically. Furthermore, we executed computational analysis for prediction of molecular interactions between β-galactosidase and NOEV. Some preliminary results of computational analysis of molecular interaction mechanism are presented in this article. NOV also showed the chaperone effect toward several β-glucosidase gene mutations in Gaucher disease. We hope chaperone therapy will become available for some patients with G_M1-gangliosidosis, Gaucher disease, and potentially other lysosomal storage diseases with central nervous system involvement.

Synthesizing non-natural parts from natural genomic template

Background

The current knowledge of genes and proteins comes from 'naturally designed' coding and non-coding regions. It would be interesting to move beyond natural boundaries and make user-defined parts. To explore this possibility we made six non-natural proteins in E. coli. We also studied their potential tertiary structure and phenotypic outcomes.

Results

The chosen intergenic sequences were amplified and expressed using pBAD 202/D-TOPO vector. All six proteins showed significantly low similarity to the known proteins in the NCBI protein database. The protein expression was confirmed through Western blot. The endogenous expression of one of the proteins resulted in the cell growth inhibition. The growth inhibition was completely rescued by culturing cells in the inducer-free medium. Computational structure prediction suggests globular tertiary structure for two of the six non-natural proteins synthesized.

Conclusion

To our best knowledge, this is the first study that demonstrates artificial synthesis of non-natural proteins from existing genomic template, their potential tertiary structure and phenotypic outcome. The work presented in this paper opens up a new avenue of investigating fundamental biology. Our approach can also be used to synthesize large numbers of non-natural RNA and protein parts for useful applications.