Limitations of Ab Initio Predictions of Peptide Binding to MHC Class II Molecules

Successful predictions of peptide MHC binding typically require a large set of binding data for the specific MHC molecule that is examined. Structure based prediction methods promise to circumvent this requirement by evaluating the physical contacts a peptide can make with an MHC molecule based on the highly conserved 3D structure of peptide:MHC complexes. While several such methods have been described before, most are not publicly available and have not been independently tested for their performance. We here implemented and evaluated three prediction methods for MHC class II molecules: statistical potentials derived from the analysis of known protein structures; energetic evaluation of different peptide snapshots in a molecular dynamics simulation; and direct analysis of contacts made in known 3D structures of peptide:MHC complexes. These methods are ab initio in that they require structural data of the MHC molecule examined, but no specific peptide:MHC binding data. Moreover, these methods retain the ability to make predictions in a sufficiently short time scale to be useful in a real world application, such as screening a whole proteome for candidate binding peptides. A rigorous evaluation of each methods prediction performance showed that these are significantly better than random, but still substantially lower than the best performing sequence based class II prediction methods available. While the approaches presented here were developed independently, we have chosen to present our results together in order to support the notion that generating structure based predictions of peptide:MHC binding without using binding data is unlikely to give satisfactory results

Nucleic Acids Res

Cells adapt to environmental changes by efficiently adjusting gene expression programs. Staphylococcus aureus, an opportunistic pathogenic bacterium, switches between defensive and offensive modes in response to quorum sensing signal. We identified and studied the structural characteristics and dynamic properties of the core regulatory circuit governing this switch by deterministic and stochastic computational methods, as well as experimentally. This module, termed here Double Selector Switch (DSS), comprises the RNA regulator RNAIII and the transcription factor Rot, defining a double-layered switch involving both transcriptional and post-transcriptional regulations. It coordinates the inverse expression of two sets of target genes, immuno-modulators and exotoxins, expressed during the defensive and offensive modes, respectively. Our computational and experimental analyses show that the DSS guarantees fine-tuned coordination of the inverse expression of its two gene sets, tight regulation, and filtering of noisy signals. We also identified variants of this circuit in other bacterial systems, suggesting it is used as a molecular switch in various cellular contexts and offering its use as a template for an effective switching device in synthetic biology studies

Ab initio identification of human microRNAs based on structure motifs

<p>Abstract</p> <p>Background</p> <p>MicroRNAs (miRNAs) are short, non-coding RNA molecules that are directly involved in post-transcriptional regulation of gene expression. The mature miRNA sequence binds to more or less specific target sites on the mRNA. Both their small size and sequence specificity make the detection of completely new miRNAs a challenging task. This cannot be based on sequence information alone, but requires structure information about the miRNA precursor. Unlike comparative genomics approaches, <it>ab initio </it>approaches are able to discover species-specific miRNAs without known sequence homology.</p> <p>Results</p> <p>MiRPred is a novel method for <it>ab initio </it>prediction of miRNAs by genome scanning that only relies on (predicted) secondary structure to distinguish miRNA precursors from other similar-sized segments of the human genome. We apply a machine learning technique, called linear genetic programming, to develop special classifier programs which include multiple regular expressions (motifs) matched against the secondary structure sequence. Special attention is paid to scanning issues. The classifiers are trained on fixed-length sequences as these occur when shifting a window in regular steps over a genome region. Various statistical and empirical evidence is collected to validate the correctness of and increase confidence in the predicted structures. Among other things, we propose a new criterion to select miRNA candidates with a higher stability of folding that is based on the number of matching windows around their genome location. An ensemble of 16 motif-based classifiers achieves 99.9 percent specificity with sensitivity remaining on an acceptable high level when requiring all classifiers to agree on a positive decision. A low false positive rate is considered more important than a low false negative rate, when searching larger genome regions for unknown miRNAs. 117 new miRNAs have been predicted close to known miRNAs on human chromosome 19. All candidate structures match the free energy distribution of miRNA precursors which is significantly shifted towards lower free energies. We employed a human EST library and found that around 75 percent of the candidate sequences are likely to be transcribed, with around 35 percent located in introns.</p> <p>Conclusion</p> <p>Our motif finding method is at least competitive to state-of-the-art feature-based methods for <it>ab initio </it>miRNA discovery. In doing so, it requires less previous knowledge about miRNA precursor structures while programs and motifs allow a more straightforward interpretation and extraction of the acquired knowledge.</p

Altered miRNA expression network in locus coeruleus of depressed suicide subjects

Norepinephrine (NE) is produced primarily by neurons in the locus coeruleus (LC). Retrograde and ultrastructural examinations reveal that the core of the LC and its surrounding region receives afferent projections from several brain areas which provide multiple neurochemical inputs to the LC with changes in LC neuronal firing, making it a highly coordinated event. Although NE and mediated signaling systems have been studied in relation to suicide and psychiatric disorders that increase the risk of suicide including depression, less is known about the corresponding changes in molecular network within LC. In this study, we examined miRNA networks in the LC of depressed suicide completers and healthy controls. Expression array revealed differential regulation of 13 miRNAs. Interaction between altered miRNAs and target genes showed dense interconnected molecular network. Functional clustering of predicated target genes yielded stress induced disorders that collectively showed the complex nature of suicidal behavior. In addition, 25 miRNAs were pairwise correlated specifically in the depressed suicide group, but not in the control group. Altogether, our study revealed for the first time the involvement of LC based dysregulated miRNA network in disrupting cellular pathways associated with suicidal behavior

Large-scale validation of methods for cytotoxic T-lymphocyte epitope prediction

<p>Abstract</p> <p>Background</p> <p>Reliable predictions of Cytotoxic T lymphocyte (CTL) epitopes are essential for rational vaccine design. Most importantly, they can minimize the experimental effort needed to identify epitopes. NetCTL is a web-based tool designed for predicting human CTL epitopes in any given protein. It does so by integrating predictions of proteasomal cleavage, TAP transport efficiency, and MHC class I affinity. At least four other methods have been developed recently that likewise attempt to predict CTL epitopes: EpiJen, MAPPP, MHC-pathway, and WAPP. In order to compare the performance of prediction methods, objective benchmarks and standardized performance measures are needed. Here, we develop such large-scale benchmark and corresponding performance measures and report the performance of an updated version 1.2 of NetCTL in comparison with the four other methods.</p> <p>Results</p> <p>We define a number of performance measures that can handle the different types of output data from the five methods. We use two evaluation datasets consisting of known HIV CTL epitopes and their source proteins. The source proteins are split into all possible 9 mers and except for annotated epitopes; all other 9 mers are considered non-epitopes. In the RANK measure, we compare two methods at a time and count how often each of the methods rank the epitope highest. In another measure, we find the specificity of the methods at three predefined sensitivity values. Lastly, for each method, we calculate the percentage of known epitopes that rank within the 5% peptides with the highest predicted score.</p> <p>Conclusion</p> <p>NetCTL-1.2 is demonstrated to have a higher predictive performance than EpiJen, MAPPP, MHC-pathway, and WAPP on all performance measures. The higher performance of NetCTL-1.2 as compared to EpiJen and MHC-pathway is, however, not statistically significant on all measures. In the large-scale benchmark calculation consisting of 216 known HIV epitopes covering all 12 recognized HLA supertypes, the NetCTL-1.2 method was shown to have a sensitivity among the 5% top-scoring peptides above 0.72. On this dataset, the best of the other methods achieved a sensitivity of 0.64. The NetCTL-1.2 method is available at <url>http://www.cbs.dtu.dk/services/NetCTL</url>.</p> <p>All used datasets are available at <url>http://www.cbs.dtu.dk/suppl/immunology/CTL-1.2.php</url>.</p

Analysis of Pools of Targeted Salmonella Deletion Mutants Identifies Novel Genes Affecting Fitness during Competitive Infection in Mice

Pools of mutants of minimal complexity but maximal coverage of genes of interest facilitate screening for genes under selection in a particular environment. We constructed individual deletion mutants in 1,023 Salmonella enterica serovar Typhimurium genes, including almost all genes found in Salmonella but not in related genera. All mutations were confirmed simultaneously using a novel amplification strategy to produce labeled RNA from a T7 RNA polymerase promoter, introduced during the construction of each mutant, followed by hybridization of this labeled RNA to a Typhimurium genome tiling array. To demonstrate the ability to identify fitness phenotypes using our pool of mutants, the pool was subjected to selection by intraperitoneal injection into BALB/c mice and subsequent recovery from spleens. Changes in the representation of each mutant were monitored using T7 transcripts hybridized to a novel inexpensive minimal microarray. Among the top 120 statistically significant spleen colonization phenotypes, more than 40 were mutations in genes with no previously known role in this model. Fifteen phenotypes were tested using individual mutants in competitive assays of intraperitoneal infection in mice and eleven were confirmed, including the first two examples of attenuation for sRNA mutants in Salmonella. We refer to the method as Array-based analysis of cistrons under selection (ABACUS)

Phage Lambda CIII: A Protease Inhibitor Regulating the Lysis-Lysogeny Decision

The ATP-dependent protease FtsH (HflB) complexed with HflKC participates in post-translational control of the lysis-lysogeny decision of bacteriophage lambda by rapid degradation of lambda CII. Both phage-encoded proteins, the CII transcription activator and the CIII polypeptide, are required for efficient lysogenic response. The conserved CIII is both an inhibitor and substrate of FtsH. Here we show that the protease inhibitor CIII is present as oligomeric amphipathic α helical structures and functions as a competitive inhibitor of FtsH by preventing binding of the CII substrate. We identified single alanine substitutions in CIII that abolish its activity. We characterize a dominant negative effect of a CIII mutant. Thus, we suggest that CIII oligomrization is required for its function. Real-time analysis of CII activity demonstrates that the effect of CIII is not seen in the absence of either FtsH or HflKC. When CIII is provided ectopically, CII activity increases linearly as a function of the multiplicity of infection, suggesting that CIII enhances CII stability and the lysogenic response. FtsH function is essential for cellular viability as it regulates the balance in the synthesis of phospholipids and lipopolysaccharides. Genetic experiments confirmed that the CIII bacteriostatic effects are due to inhibition of FtsH. Thus, the early presence of CIII following infection stimulates the lysogenic response, while its degradation at later times ensures the reactivation of FtsH allowing the growth of the established lysogenic cell

Three Essential Ribonucleases—RNase Y, J1, and III—Control the Abundance of a Majority of Bacillus subtilis mRNAs

Bacillus subtilis possesses three essential enzymes thought to be involved in mRNA decay to varying degrees, namely RNase Y, RNase J1, and RNase III. Using recently developed high-resolution tiling arrays, we examined the effect of depletion of each of these enzymes on RNA abundance over the whole genome. The data are consistent with a model in which the degradation of a significant number of transcripts is dependent on endonucleolytic cleavage by RNase Y, followed by degradation of the downstream fragment by the 5′–3′ exoribonuclease RNase J1. However, many full-size transcripts also accumulate under conditions of RNase J1 insufficiency, compatible with a model whereby RNase J1 degrades transcripts either directly from the 5′ end or very close to it. Although the abundance of a large number of transcripts was altered by depletion of RNase III, this appears to result primarily from indirect transcriptional effects. Lastly, RNase depletion led to the stabilization of many low-abundance potential regulatory RNAs, both in intergenic regions and in the antisense orientation to known transcripts

NetCTLpan: pan-specific MHC class I pathway epitope predictions

Reliable predictions of immunogenic peptides are essential in rational vaccine design and can minimize the experimental effort needed to identify epitopes. In this work, we describe a pan-specific major histocompatibility complex (MHC) class I epitope predictor, NetCTLpan. The method integrates predictions of proteasomal cleavage, transporter associated with antigen processing (TAP) transport efficiency, and MHC class I binding affinity into a MHC class I pathway likelihood score and is an improved and extended version of NetCTL. The NetCTLpan method performs predictions for all MHC class I molecules with known protein sequence and allows predictions for 8-, 9-, 10-, and 11-mer peptides. In order to meet the need for a low false positive rate, the method is optimized to achieve high specificity. The method was trained and validated on large datasets of experimentally identified MHC class I ligands and cytotoxic T lymphocyte (CTL) epitopes. It has been reported that MHC molecules are differentially dependent on TAP transport and proteasomal cleavage. Here, we did not find any consistent signs of such MHC dependencies, and the NetCTLpan method is implemented with fixed weights for proteasomal cleavage and TAP transport for all MHC molecules. The predictive performance of the NetCTLpan method was shown to outperform other state-of-the-art CTL epitope prediction methods. Our results further confirm the importance of using full-type human leukocyte antigen restriction information when identifying MHC class I epitopes. Using the NetCTLpan method, the experimental effort to identify 90% of new epitopes can be reduced by 15% and 40%, respectively, when compared to the NetMHCpan and NetCTL methods. The method and benchmark datasets are available at http://www.cbs.dtu.dk/services/NetCTLpan/