Predicting class I major histocompatibility complex (MHC) binders using multivariate statistics:comparison of discriminant analysis and multiple linear regression

The accurate in silico identification of T-cell epitopes is a critical step in the development of peptide-based vaccines, reagents, and diagnostics. It has a direct impact on the success of subsequent experimental work. Epitopes arise as a consequence of complex proteolytic processing within the cell. Prior to being recognized by T cells, an epitope is presented on the cell surface as a complex with a major histocompatibility complex (MHC) protein. A prerequisite therefore for T-cell recognition is that an epitope is also a good MHC binder. Thus, T-cell epitope prediction overlaps strongly with the prediction of MHC binding. In the present study, we compare discriminant analysis and multiple linear regression as algorithmic engines for the definition of quantitative matrices for binding affinity prediction. We apply these methods to peptides which bind the well-studied human MHC allele HLA-A*0201. A matrix which results from combining results of the two methods proved powerfully predictive under cross-validation. The new matrix was also tested on an external set of 160 binders to HLA-A*0201; it was able to recognize 135 (84%) of them.

EpiJen:a server for multistep T cell epitope prediction

Background - The main processing pathway for MHC class I ligands involves degradation of proteins by the proteasome, followed by transport of products by the transporter associated with antigen processing (TAP) to the endoplasmic reticulum (ER), where peptides are bound by MHC class I molecules, and then presented on the cell surface by MHCs. The whole process is modeled here using an integrated approach, which we call EpiJen. EpiJen is based on quantitative matrices, derived by the additive method, and applied successively to select epitopes. EpiJen is available free online. Results - To identify epitopes, a source protein is passed through four steps: proteasome cleavage, TAP transport, MHC binding and epitope selection. At each stage, different proportions of non-epitopes are eliminated. The final set of peptides represents no more than 5% of the whole protein sequence and will contain 85% of the true epitopes, as indicated by external validation. Compared to other integrated methods (NetCTL, WAPP and SMM), EpiJen performs best, predicting 61 of the 99 HIV epitopes used in this study. Conclusion - EpiJen is a reliable multi-step algorithm for T cell epitope prediction, which belongs to the next generation of in silico T cell epitope identification methods. These methods aim to reduce subsequent experimental work by improving the success rate of epitope prediction.

Quantitative structure-activity relationships and the prediction of MHC supermotifs

The underlying assumption in quantitative structure–activity relationship (QSAR) methodology is that related chemical structures exhibit related biological activities. We review here two QSAR methods in terms of their applicability for human MHC supermotif definition. Supermotifs are motifs that characterise binding to more than one allele. Supermotif definition is the initial in silico step of epitope-based vaccine design. The first QSAR method we review here—the additive method—is based on the assumption that the binding affinity of a peptide depends on contributions from both amino acids and the interactions between them. The second method is a 3D-QSAR method: comparative molecular similarity indices analysis (CoMSIA). Both methods were applied to 771 peptides binding to 9 HLA alleles. Five of the alleles (A*0201, A* 0202, A*0203, A*0206 and A*6802) belong to the HLA-A2 superfamily and the other four (A*0301, A*1101, A*3101 and A*6801) to the HLA-A3 superfamily. For each superfamily, supermotifs defined by the two QSAR methods agree closely and are supported by many experimental data.

Transporter associated with antigen processing preselection of peptides binding to the MHC:a bioinformatic evaluation

TAP is responsible for the transit of peptides from the cytosol to the lumen of the endoplasmic reticulum. In an immunological context, this event is followed by the binding of peptides to MHC molecules before export to the cell surface and recognition by T cells. Because TAP transport precedes MHC binding, TAP preferences may make a significant contribution to epitope selection. To assess the impact of this preselection, we have developed a scoring function for TAP affinity prediction using the additive method, have used it to analyze and extend the TAP binding motif, and have evaluated how well this model acts as a preselection step in predicting MHC binding peptides. To distinguish between MHC alleles that are exclusively dependent on TAP and those exhibiting only a partial dependence on TAP, two sets of MHC binding peptides were examined: HLA-A*0201 was selected as a representative of partially TAP-dependent HLA alleles, and HLA-A*0301 represented fully TAP-dependent HLA alleles. TAP preselection has a greater impact on TAP-dependent alleles than on TAP-independent alleles. The reduction in the number of nonbinders varied from 10% (TAP-independent) to 33% (TAP-dependent), suggesting that TAP preselection is an important component in the successful in silico prediction of T cell epitopes.

Designing immunogenic peptides

Peptides fulfill many roles in immunology, yet none are more important than their role as immunogenic epitopes driving the adaptive immune response, our ultimate bulwark against infectious disease. Peptide epitopes are mediated primarily by their interaction with major histocompatibility complexes (T-cell epitopes) and antibodies (B-cell epitopes). As pathogen genomes continue to be revealed, both experimental and computational epitope mapping are becoming crucial tools in vaccine discovery1,2. Immunoinformatics offers many tools, techniques and approaches for in silico epitope characterization, which is capable of greatly accelerating epitope design. © 2013 Nature America, Inc. All rights reserved.

The handedness-associated PCSK6 locus spans an intronic promoter regulating novel transcripts

We recently reported the association of the PCSK6 gene with handedness through a quantitative genome-wide association study (GWAS; P < 0.5 × 10(-8)) for a relative hand skill measure in individuals with dyslexia. PCSK6 activates Nodal, a morphogen involved in regulating left-right body axis determination. Therefore, the GWAS data suggest that the biology underlying the patterning of structural asymmetries may also contribute to behavioural laterality, e.g. handedness. The association is further supported by an independent study reporting a variable number tandem repeat (VNTR) within the same PCSK6 locus to be associated with degree of handedness in a general population cohort. Here, we have conducted a functional analysis of the PCSK6 locus combining further genetic analysis, in silico predictions and molecular assays. We have shown that the previous GWAS signal was not tagging a VNTR effect, suggesting that the two markers have independent effects. We demonstrated experimentally that one of the top GWAS-associated markers, rs11855145, directly alters the binding site for a nuclear factor. Furthermore, we have shown that the predicted regulatory region adjacent to rs11855415 acts as a bidirectional promoter controlling the expression of novel RNA transcripts. These include both an antisense long non-coding RNA (lncRNA) and a short PCSK6 isoform predicted to be coding. This is the first molecular characterization of a handedness-associated locus that supports the role of common variants in non-coding sequences in influencing complex phenotypes through gene expression regulation.

Better primer design for metagenomics applications by increasing taxonomic distinguishability

Current methods of understanding microbiome composition and structure rely on accurately estimating the number of distinct species and their relative abundance. Most of these methods require an efficient PCR whose forward and reverse primers bind well to the same, large number of identifiable species, and produce amplicons that are unique. It is therefore not surprising that currently used universal primers designed many years ago are not as efficient and fail to bind to recently cataloged species. We propose an automated general method of designing PCR primer pairs that abide by primer design rules and uses current sequence database as input. Since the method is automated, primers can be designed for targeted microbial species or updated as species are added or deleted from the database. In silico experiments and laboratory experiments confirm the efficacy of the newly designed primers for metagenomics applications.

Mechanisms of Arsenic Detoxification and Resistance

Arsenic is a ubiquitous environmental toxic substance. As a consequence of continual exposure to arsenic, nearly every organism, from Escherichia coli to humans have evolved arsenic detoxification pathways. One of the pathways is extrusion of arsenic from inside the cells, thereby conferring resistance. The R773 arsRDABC operon in E. coli encodes an ArsAB efflux pump that confers resistance to arsenite. ArsA is the catalytic subunit of the pump, while ArsB forms the oxyanion conducting pathway. ArsD is an arsenite metallochaperone that binds arsenite and transfers it to ArsA. The interaction of ArsA and ArsD allows for resistance to As(III) at environmental concentrations. The interaction between ArsA ATPase and ArsD metallochaperone was examined. A quadruple mutant in the arsD gene encoding a K2A/K37A/K62A/K104A ArsD is unable to interact with ArsA. An error-prone mutagenesis approach was used to generate random mutations in the arsA gene that restored interaction with the quadruple arsD mutant in yeast two-hybrid assays. Three such mutants encoding Q56R, F120I and D137V ArsA were able to restore interaction with the quadruple ArsD mutant. Structural models generated by in silico docking suggest that an electrostatic interface favors reversible interaction between ArsA and ArsD. Mutations in ArsA that propagate changes in hydrogen bonding and salt bridges to the ArsA-ArsD interface also affect their interactions. The second objective was to examine the mechanism of arsenite resistance through methylation and subsequent volatilization. Microbial ArsM (As(III) S-adenosylmethyltransferase) catalyzes the formation of trimethylarsine as the volatile end product. The net result is loss of arsenic from cells. The gene for CrArsM from the eukaryotic green alga Chlamydomonas reinhardtii was chemically synthesized and expressed in E. coli. The purified protein catalyzed the methylation of arsenite into methyl-, dimethyl- and trimethyl products. Synthetic purified CrArsM was crystallized in an unliganded form. Biochemical and biophysical studies conducted on CrArsM sheds new light on the pathways of biomethylation. While in microbes ArsM detoxifies arsenic, the human homolog, hAS3MT, converts inorganic arsenic into more toxic and carcinogenic forms. An understanding of the enzymatic mechanism of ArsM will be critical in deciphering its parallel roles in arsenic detoxification and carcinogenesis.

Filogenia molecular das enzimas isocitrato liase e malato sintase e sua evolução em Viridiplanta

The plant metabolism consists of a complex network of physical and chemical events resulting in photosynthesis, respiration, synthesis and degradation of organic compounds. This is only possible due to the different kinds of responses to many environmental variations that a plant could be subject through evolution, leading also to conquering new surroundings. The glyoxylate cycle is a metabolic pathway found in glyoxysomes plant, which has unique role in the seedling establishment. Considered as a variation of the citric acid cycle, it uses an acetyl coenzyme A molecule, derived from lipids beta-oxidation to synthesize compounds which are used in carbohydrate synthesis. The Malate synthase (MLS) and Isocitrate lyase (ICL) enzyme of this cycle are unique and essential in regulating the biosynthesis of carbohydrates. Because of the absence of decarboxylation steps as rate-limiting steps, detailed studies of molecular phylogeny and evolution of these proteins enables the elucidation of the effects of this route presence in the evolutionary processes involved in their distribution across the genome from different plant species. Therefore, the aim of this study was to establish a relationship between the molecular evolution of the characteristics of enzymes from the glyoxylate cycle (isocitrate lyase and malate synthase) and their molecular phylogeny, among green plants (Viridiplantae). For this, amino acid and nucleotide sequences were used, from online repositories as UniProt and Genbank. Sequences were aligned and then subjected to an analysis of the best-fit substitution models. The phylogeny was rebuilt by distance methods (neighbor-joining) and discrete methods (maximum likelihood, maximum parsimony and Bayesian analysis). The identification of structural patterns in the evolution of the enzymes was made through homology modeling and structure prediction from protein sequences. Based on comparative analyzes of in silico models and from the results of phylogenetic inferences, both enzymes show significant structure conservation and their topologies in agreement with two processes of selection and specialization of the genes. Thus, confirming the relevance of new studies to elucidate the plant metabolism from an evolutionary perspective

Utilização de peptídeos sintéticos derivados de moléculas imunodominantes de Toxoplasma gondii para diagnóstico sorológico da toxoplasmose em ovinos e galinhas caipiras

Chapter I - The obligate intracellular parasite Toxoplasma gondii is the causative agent of toxoplasmosis, a disease that affects humans and generates economic losses in farm animals. When prevention fails disease refers to the diagnosis and subsequent treatment if the individual is diagnosed as positive. Therefore, the development of new accurate diagnostic tools for detecting T. gondii infection is a need in particular to determine the environmental source of infection which can result in more appropriate public health policies against different routes of infection and prevent potential damage that toxoplasmosis can cause when animals are infected. Chapter II - The domestic chicken (Gallus gallus domesticus) are considered epidemiological sentinels, still representing a major source of recombinant strains when predated by cats, it is common to find them with multiple infections. We evaluate the diagnostic potential of six synthetic peptides (SAG2Y, MIC1, M2AP, GRA10, ROP2 and ROP7) predicted in silico from tachyzoites immunodominant markers of T. gondii in samples from naturally infected chickens, comparing synthetic peptides with antigen total soluble (STAg). In general, our results demonstrated that reactivity rates and positivity for these peptides are similar to the STAg, and the ROP7 peptide and the combination of peptides MIC1+ROP2 have significant sensitivity, confirming them as potential diagnostic tools for the diagnosis of toxoplasmosis in chickens. Chapter III - Sheep (Ovis aries) are commonly infected with Toxoplasma gondii due to his eating habits. Infection in pregnant sheep can have serious consequences such as embryonic death, fetal resorption, mummification, and neonatal death. One concern regarding the infection in these animals is that the meat can be a source of contamination to humans and other carnivores. Therefore perform accurate diagnosis in these animals is of fundamental importance. In the present study we evaluated the potential of new synthetic peptides as a diagnostic tool. Synthetic peptides (SAG2Y, SRS52A, MIC14, GRA4, GRA10 and GRA15) were predicted in silico from immunodominant proteins of T. gondii. We determine the levels of IgG antibodies using sera obtained from two farms in the city of Uberlândia. Analyzing the results together, the peptide combination GRA10+GRA15 (Accuracy = 0,82) showed better characteristics compared with the other mixtures. This preparation could be better analyzed with an antigenic mixture potential use in the diagnosis of toxoplasmosis in sheep and other species.

Développement d'approches quantitatives de type "structure-activité" pour la modélisation pharmacocinétique

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Identification de gènes ciblès par ETV6-AML1, un facteur de transcription chimérique retrouvé dans la leucémie de l'enfant

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

An investigation of the molecular interactions between statins and bacterial pathogens, and their combined impact on the human immune system

Statins are a class of drug that inhibits cholesterol biosynthesis, and are used to treat patients with high serum cholesterol levels. They exert this function by competitively binding to the enzyme 3-hydroxy-3-methylglutaryl-CoenzymeA reductase (HMGR), which catalyses the formation of mevalonate, a rate-limiting step in cholesterol biosynthesis. In addition, statins have what are called “pleiotropic effects”, which include the reduction of inflammation, immunomodulation, and antimicrobial effects. Statins can also improve survival of patients with sepsis and pneumonia. Cystic fibrosis (CF) is the most common recessive inherited disease in the Caucasian population, which is characterised by factors including, but not limited to, excessive lung inflammation and increased susceptibility to infection. Therefore, the overall objective of this study was to examine the effects of statins on CFassociated bacterial pathogens and the host response. In this work, the prevalence of HMGR was examined in respiratory pathogens, and several CF-associated pathogens were found to possess homologues of this enzyme. HMGR homology was analysed in Staphylococcus aureus, Burkholderia cenocepacia and Streptococcus pneumoniae, and the HMGR of B. cenocepacia was found to have significant conservation to that of Pseudomonas mevalonii, which is the most widely-characterised bacterial HMGR. However, in silico analysis revealed that, unlike S. aureus and S. pneumoniae, B. cenocepacia did not possess homologues of other mevalonate pathway proteins, and that the HMGR of B. cenocepacia appeared to be involved in an alternative metabolic pathway. The effect of simvastatin was subsequently tested on the growth and virulence of S. aureus, B. cenocepacia and S. pneumoniae. Simvastatin inhibited the growth of all 3 species in a dose-dependent manner. In addition, statin treatment also attenuated biofilm formation of all 3 species, and reduced in vitro motility of S. aureus. Interestingly, simvastatin also increased the potency of the aminoglycoside antibiotic gentamicin against B. cenocepacia. The impact of statins was subsequently tested on the predominant CF-associated pathogen Pseudomonas aeruginosa, which does not possess a HMGR homologue. Mevastatin, lovastatin and simvastatin did not influence the growth of this species. However, sub-inhibitory statin concentrations reduced the swarming motility and biofilm formation of P. aeruginosa. The influence of statins was also examined on Type 3 toxin secretion, quorum sensing and chemotaxis, and no statin effect was observed on any of these phenotypes. Statins did not appear to have a characteristic effect on the P. aeruginosa transcriptome. However, a mutant library screen revealed that the effect of statins on P. aeruginosa biofilm was mediated through the PvrR regulator and the Cup fimbrial biosynthesis genes. Furthermore, proteomic analysis demonstrated that 6 proteins were reproducibly induced by simvastatin in the P. aeruginosa swarming cells. The effect of statins on the regulation of the host-P. aeruginosa immune response was also investigated. Statin treatment increased expression of the pro-inflammatory cytokines IL-8 and CCL20 in lung epithelial cells, but did not attenuate P. aeruginosa-mediated inflammatory gene induction. In fact, simvastatin and P. aeruginosa caused a synergistic effect on CCL20 expression. The expression of the transcriptional regulators KLF2 and KLF6 was also increased by statins and P. aeruginosa, with the induction of KLF6 by simvastatin proving to be a novel effect. Interestingly, both statins and P. aeruginosa were capable of inducing alternative splicing of KLF6. P. aeruginosa was found to induce KLF6 alternative splicing by way of the type 3 secreted toxin ExoS. In addition, a mechanistic role was elucidated for KLF6 in the lung, as it was determined that statin-mediated induction of this protein was responsible for the induction of the host response genes CCL20 and iNOS. Moreover, statin treatment caused a slight increase in infection-related cytotoxicity, and increased bacterial adhesion to cells. Taken together, these data demonstrate that statins can reduce the virulence of CFassociated bacterial pathogens and alter host response effectors. Furthermore, novel statin effectors were identified in both bacterial and host cells.

Neural Network Evidence for the Coupling of Presaccadic Visual Remapping to Predictive Eye Position Updating.

As we look around a scene, we perceive it as continuous and stable even though each saccadic eye movement changes the visual input to the retinas. How the brain achieves this perceptual stabilization is unknown, but a major hypothesis is that it relies on presaccadic remapping, a process in which neurons shift their visual sensitivity to a new location in the scene just before each saccade. This hypothesis is difficult to test in vivo because complete, selective inactivation of remapping is currently intractable. We tested it in silico with a hierarchical, sheet-based neural network model of the visual and oculomotor system. The model generated saccadic commands to move a video camera abruptly. Visual input from the camera and internal copies of the saccadic movement commands, or corollary discharge, converged at a map-level simulation of the frontal eye field (FEF), a primate brain area known to receive such inputs. FEF output was combined with eye position signals to yield a suitable coordinate frame for guiding arm movements of a robot. Our operational definition of perceptual stability was "useful stability," quantified as continuously accurate pointing to a visual object despite camera saccades. During training, the emergence of useful stability was correlated tightly with the emergence of presaccadic remapping in the FEF. Remapping depended on corollary discharge but its timing was synchronized to the updating of eye position. When coupled to predictive eye position signals, remapping served to stabilize the target representation for continuously accurate pointing. Graded inactivations of pathways in the model replicated, and helped to interpret, previous in vivo experiments. The results support the hypothesis that visual stability requires presaccadic remapping, provide explanations for the function and timing of remapping, and offer testable hypotheses for in vivo studies. We conclude that remapping allows for seamless coordinate frame transformations and quick actions despite visual afferent lags. With visual remapping in place for behavior, it may be exploited for perceptual continuity.

Interaction between Candida albicans and Pseudomonas aeruginosa

Fungal pathogen Candida albicans causes serious nosocomial infections in patients, in part, due to formation of drug-resistant biofilms. Protein kinases (PK) and transcription factors (TF) mediate signal transduction and transcription of proteins involved in biofilm development. To discover biofilm-related PKs, a collection of 63 C. albicans PK mutants was screened twice independently with microtiter plate-based biofilm assay (XTT). Thirty-eight (60%) mutants showed different degrees of biofilm impairment with the poor biofilm formers additionally possessing filamentation defects. Most of these genes were already known to encode proteins associated with Candida morphology and biofilms but VPS15, PKH3, PGA43, IME2 and CEX1, were firstly associated with both processes in this study. Previous studies of Holcombe et al. (2010) had shown that bacterial pathogen, Pseudomonas aeruginosa can impair C. albicans filamentation and biofilm development. To investigate their interaction, the good biofilm former PK mutants of C. albicans were assessed for their response to P. aeruginosa supernatants derived from two strains, wildtype PAO1 and homoserine lactone (HSL)-free mutant ΔQS, without finding any nonresponsive mutants. This suggested that none of the PKs in this study was implicated in Candida-Pseudomonas signaling. To screen promoter sequences for overrepresented TFs across C. albicans gene sets significantly up/downregulated in presence of bacterial supernatants from Holcombe et al. (2010) study, TFbsST database was created online. The TFbsST database integrates experimentally verified TFs of Candida to analyse promoter sequences for TF binding sites. In silico studies predicted that Efg1p was overrepresented in C. albicans and C. parapsilosis RBT family genes.