The binding site distance test score: a robust method for the assessment of predicted protein binding sites

We propose a novel method for scoring the accuracy of protein binding site predictions – the Binding-site Distance Test (BDT) score. Recently, the Matthews Correlation Coefficient (MCC) has been used to evaluate binding site predictions, both by developers of new methods and by the assessors for the community wide prediction experiment – CASP8. Whilst being a rigorous scoring method, the MCC does not take into account the actual 3D location of the predicted residues from the observed binding site. Thus, an incorrectly predicted site that is nevertheless close to the observed binding site will obtain an identical score to the same number of nonbinding residues predicted at random. The MCC is somewhat affected by the subjectivity of determining observed binding residues and the ambiguity of choosing distance cutoffs. By contrast the BDT method produces continuous scores ranging between 0 and 1, relating to the distance between the predicted and observed residues. Residues predicted close to the binding site will score higher than those more distant, providing a better reflection of the true accuracy of predictions. The CASP8 function predictions were evaluated using both the MCC and BDT methods and the scores were compared. The BDT was found to strongly correlate with the MCC scores whilst also being less susceptible to the subjectivity of defining binding residues. We therefore suggest that this new simple score is a potentially more robust method for future evaluations of protein-ligand binding site predictions.

S-nitrosylation of proteins at the leading edge of migrating trophoblasts by inducible nitric oxide synthase promotes trophoblast invasion

Nitric oxide regulates many important cellular processes including motility and invasion. Many of its effects are mediated through the modification of specific cysteine residues in target proteins, a process called S-nitrosylation. Here we show that S-nitrosylation of proteins occurs at the leading edge of migrating trophoblasts and can be attributed to the specific enrichment of inducible nitric oxide synthase (iNOS/NOS2) in this region. Localisation of iNOS to the leading edge is co-incidental with a site of extensive actin polymerisation and is only observed in actively migrating cells. In contrast endothelial nitric oxide synthase (eNOS/NOS3) shows distribution that is distinct and non-colocalised with iNOS, suggesting that the protein S-nitrosylation observed at the leading edge is caused only by iNOS and not eNOS. We have identified MMP-9 as a potential target for S-nitrosylation in these cells and demonstrate that it co-localises with iNOS at the leading edge of migrating cells. We further demonstrate that iNOS plays an important role in promoting trophoblast invasion, which is an essential process in the establishment of a successful pregnancy.

Fas ligand-induced apoptosis is regulated by nitric oxide through the inhibition of fas receptor clustering and the nitrosylation of protein kinase Cepsilon

Apoptosis induced by the death-inducing ligand FasL (CD95L) is a major mechanism of cell death. Trophoblast cells express the Fas receptor yet survive in an environment that is rich in the ligand. We report that basal nitric oxide (NO) production is responsible for the resistance of trophoblasts to FasL-induced apoptosis. In this study we demonstrate that basal NO production resulted in the inhibition of receptor clustering following ligand binding. In addition NO also protected cells through the selective nitrosylation, and inhibition, of protein kinase Cepsilon (PKCepsilon) but not PKCalpha. In the absence of NO production PKCepsilon interacted with, and phosphorylated, the anti-apoptotic protein cFLIP. The interaction is predominantly with the short form of cFLIP and its phosphorylation reduces its recruitment to the death-inducing signaling complex (DISC) that is formed following binding of a death-inducing ligand to its receptor. Inhibition of cFLIP recruitment to the DISC leads to increased activation of caspase 8 and subsequently to apoptosis. Inhibition of PKCepsilon using siRNA significantly reversed the sensitivity to apoptosis induced by inhibition of NO synthesis suggesting that NO-mediated inhibition of PKCepsilon plays an important role in the regulation of Fas-induced apoptosis.

Teachers' networks in and out of school

A 'mapping task' was used to explore the networks available to head teachers, school coordinators and local authority staff. Beginning from an ego-centred perspective on networks, we illustrate a number of key analytic categories, including brokerage, formality, and strength and weakness of links with reference to a single UK primary school. We describe how teachers differentiate between the strength of network links and their value, which is characteristically related to their potential impact on classroom practice.

Proteogenomics and in silico structural and functional annotation of the barley powdery mildew Blumeria graminis f. sp. hordei

Blumeria graminis is an economically important obligate plant-pathogenic fungus, whose entire genome was recently sequenced and manually annotated using ab initio in silico predictions [7]. Employing large scale proteogenomic analysis we are now able to verify independently the existence of proteins predicted by 24% of open reading frame models. We compared the haustoria and sporulating hyphae proteomes and identified 71 proteins exclusively in haustoria, the feeding and effector-delivery organs of the pathogen. These proteins are ‘significantly smaller than the rest of the protein pool and predicted to be secreted. Most do not share any similarities with Swiss–Prot or Trembl entries nor possess any identifiable Pfam domains. We used a novel automated prediction pipeline to model the 3D structures of the proteins, identify putative ligand binding sites and predict regions of intrinsic disorder. This revealed that the protein set found exclusively in haustoria is significantly less disordered than the rest of the identified Blumeria proteins or random (and representative) protein sets generated from the yeast proteome. For most of the haustorial proteins with unknown functions no good templates could be found, from which to generate high quality models. Thus, these unknown proteins present potentially new protein folds that can be specific to the interaction of the pathogen with its host.

Automated tertiary structure prediction with accurate local model quality assessment using the IntFOLD-TS method

The IntFOLD-TS method was developed according to the guiding principle that the model quality assessment would be the most critical stage for our template based modelling pipeline. Thus, the IntFOLD-TS method firstly generates numerous alternative models, using in-house versions of several different sequence-structure alignment methods, which are then ranked in terms of global quality using our top performing quality assessment method – ModFOLDclust2. In addition to the predicted global quality scores, the predictions of local errors are also provided in the resulting coordinate files, using scores that represent the predicted deviation of each residue in the model from the equivalent residue in the native structure. The IntFOLD-TS method was found to generate high quality 3D models for many of the CASP9 targets, whilst also providing highly accurate predictions of their per-residue errors. This important information may help to make the 3D models that are produced by the IntFOLD-TS method more useful for guiding future experimental work

A novel, non-canonical mechanism of regulation of MST3 (mammalian Sterile20-related kinase 3)

The canonical pathway of regulation of the germinal centre kinase (GCK) III subgroup member, mammalian Sterile20-related kinase 3 (MST3), involves a caspase-mediated cleavage between N-terminal catalytic and C-terminal regulatory domains with possible concurrent autophosphorylation of the activation loop MST3(Thr178-), induction of Ser-/Thr-protein kinase activity and nuclear localisation. We identified an alternative ‘non-canonical’ pathway of MST3 activation (regulated primarily through dephosphorylation) which may also be applicable to other GCKIII (and GCKVI) subgroup members. In the basal state, inactive MST3 co-immunoprecipitated with the Golgi protein, GOLGA2/gm130. Activation of MST3 by calyculin A (a protein Ser-/Thr- phosphatase 1/2A inhibitor) stimulated (auto)phosphorylation of MST3(Thr178-) in the catalytic domain with essentially simultaneous cis-autophosphorylation of MST3(Thr328-) in the regulatory domain, an event also requiring the MST3(341-376) sequence which acts as a putative docking domain. MST3(Thr178-) phosphorylation increased MST3 kinase activity but this activity was independent of MST3(Thr328-) phosphorylation. Interestingly, MST3(Thr328-) lies immediately C-terminal to a STRAD pseudokinase-like site recently identified as being involved in binding of GCKIII/GCKVI members to MO25 scaffolding proteins. MST3(Thr178- /Thr328-) phosphorylation was concurrent with dissociation of MST3 from GOLGA2/gm130 and association of MST3 with MO25, and MST3(Thr328-) phosphorylation was necessary for formation of the activated MST3-MO25 holocomplex.

Improving sequence-based fold recognition by using 3D model quality assessment

Motivation: The ability of a simple method (MODCHECK) to determine the sequence–structure compatibility of a set of structural models generated by fold recognition is tested in a thorough benchmark analysis. Four Model Quality Assessment Programs (MQAPs) were tested on 188 targets from the latest LiveBench-9 automated structure evaluation experiment. We systematically test and evaluate whether the MQAP methods can successfully detect native-likemodels. Results: We show that compared with the other three methods tested MODCHECK is the most reliable method for consistently performing the best top model selection and for ranking the models. In addition, we show that the choice of model similarity score used to assess a model's similarity to the experimental structure can influence the overall performance of these tools. Although these MQAP methods fail to improve the model selection performance for methods that already incorporate protein three dimension (3D) structural information, an improvement is observed for methods that are purely sequence-based, including the best profile–profile methods. This suggests that even the best sequence-based fold recognition methods can still be improved by taking into account the 3D structural information.

Prediction of novel and analogous folds using fragment assembly and fold recognition

A number of new and newly improved methods for predicting protein structure developed by the Jones–University College London group were used to make predictions for the CASP6 experiment. Structures were predicted with a combination of fold recognition methods (mGenTHREADER, nFOLD, and THREADER) and a substantially enhanced version of FRAGFOLD, our fragment assembly method. Attempts at automatic domain parsing were made using DomPred and DomSSEA, which are based on a secondary structure parsing algorithm and additionally for DomPred, a simple local sequence alignment scoring function. Disorder prediction was carried out using a new SVM-based version of DISOPRED. Attempts were also made at domain docking and “microdomain” folding in order to build complete chain models for some targets.

Protein structure prediction servers at University College London

A number of state-of-the-art protein structure prediction servers have been developed by researchers working in the Bioinformatics Unit at University College London. The popular PSIPRED server allows users to perform secondary structure prediction, transmembrane topology prediction and protein fold recognition. More recent servers include DISOPRED for the prediction of protein dynamic disorder and DomPred for domain boundary prediction.