Toward the prediction of class I and II mouse major histocompatibility complex-peptide-binding affinity:in silico bioinformatic step-by-step guide using quantitative structure-activity relationships

Quantitative structure-activity relationship (QSAR) analysis is a cornerstone of modern informatics. Predictive computational models of peptide-major histocompatibility complex (MHC)-binding affinity based on QSAR technology have now become important components of modern computational immunovaccinology. Historically, such approaches have been built around semiqualitative, classification methods, but these are now giving way to quantitative regression methods. We review three methods--a 2D-QSAR additive-partial least squares (PLS) and a 3D-QSAR comparative molecular similarity index analysis (CoMSIA) method--which can identify the sequence dependence of peptide-binding specificity for various class I MHC alleles from the reported binding affinities (IC50) of peptide sets. The third method is an iterative self-consistent (ISC) PLS-based additive method, which is a recently developed extension to the additive method for the affinity prediction of class II peptides. The QSAR methods presented here have established themselves as immunoinformatic techniques complementary to existing methodology, useful in the quantitative prediction of binding affinity: current methods for the in silico identification of T-cell epitopes (which form the basis of many vaccines, diagnostics, and reagents) rely on the accurate computational prediction of peptide-MHC affinity. We have reviewed various human and mouse class I and class II allele models. Studied alleles comprise HLA-A*0101, HLA-A*0201, HLA-A*0202, HLA-A*0203, HLA-A*0206, HLA-A*0301, HLA-A*1101, HLA-A*3101, HLA-A*6801, HLA-A*6802, HLA-B*3501, H2-K(k), H2-K(b), H2-D(b) HLA-DRB1*0101, HLA-DRB1*0401, HLA-DRB1*0701, I-A(b), I-A(d), I-A(k), I-A(S), I-E(d), and I-E(k). In this chapter we show a step-by-step guide into predicting the reliability and the resulting models to represent an advance on existing methods. The peptides used in this study are available from the AntiJen database (http://www.jenner.ac.uk/AntiJen). The PLS method is available commercially in the SYBYL molecular modeling software package. The resulting models, which can be used for accurate T-cell epitope prediction, will be made are freely available online at the URL http://www.jenner.ac.uk/MHCPred.

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.

New horizons in mouse immunoinformatics:reliable in silico prediction of mouse class I histocompatibility major complex peptide binding affinity

Quantitative structure–activity relationship (QSAR) analysis is a main cornerstone of modern informatic disciplines. Predictive computational models, based on QSAR technology, of peptide-major histocompatibility complex (MHC) binding affinity have now become a vital component of modern day computational immunovaccinology. Historically, such approaches have been built around semi-qualitative, classification methods, but these are now giving way to quantitative regression methods. The additive method, an established immunoinformatics technique for the quantitative prediction of peptide–protein affinity, was used here to identify the sequence dependence of peptide binding specificity for three mouse class I MHC alleles: H2–Db, H2–Kb and H2–Kk. As we show, in terms of reliability the resulting models represent a significant advance on existing methods. They can be used for the accurate prediction of T-cell epitopes and are freely available online (http://www.jenner.ac.uk/MHCPred).

Life without complex I: proteome analyses of an Arabidopsis mutant lacking the mitochondrial NADH dehydrogenase complex

The mitochondrial NADH dehydrogenase complex (complex I) is of particular importance for the respiratory chain in mitochondria. It is the major electron entry site for the mitochondrial electron transport chain (mETC) and therefore of great significance for mitochondrial ATP generation. We recently described an Arabidopsis thaliana double-mutant lacking the genes encoding the carbonic anhydrases CA1 and CA2, which both form part of a plant-specific 'carbonic anhydrase domain' of mitochondrial complex I. The mutant lacks complex I completely. Here we report extended analyses for systematically characterizing the proteome of the ca1ca2 mutant. Using various proteomic tools, we show that lack of complex I causes reorganization of the cellular respiration system. Reduced electron entry into the respiratory chain at the first segment of the mETC leads to induction of complexes II and IV as well as alternative oxidase. Increased electron entry at later segments of the mETC requires an increase in oxidation of organic substrates. This is reflected by higher abundance of proteins involved in glycolysis, the tricarboxylic acid cycle and branched-chain amino acid catabolism. Proteins involved in the light reaction of photosynthesis, the Calvin cycle, tetrapyrrole biosynthesis, and photorespiration are clearly reduced, contributing to the significant delay in growth and development of the double-mutant. Finally, enzymes involved in defense against reactive oxygen species and stress symptoms are much induced. These together with previously reported insights into the function of plant complex I, which were obtained by analysing other complex I mutants, are integrated in order to comprehensively describe 'life without complex I'.

Mutation in NDUFA13/GRIM19 leads to early onset hypotonia, dyskinesia and sensorial deficiencies, and mitochondrial complex I instability

International audience

Rapamycin reverses age-related increases in mitochondrial ROS production at complex I, oxidative stress, accumulation of mtDNA fragments inside nuclear DNA, and lipofuscin level, and increases autophagy, in the liver of middle-aged mice

Rapamycin consistently increases longevity in mice although the mechanism of action of this drug is unknown. In the present investigation we studied the effect of rapamycin on mitochondrial oxidative stress at the same dose that is known to increase longevity in mice (14 mg of rapamycin/kg of diet). Middle aged mice (16 months old) showed significant age-related increases in mitochondrial ROS production at complex I, accumulation of mtDNA fragments inside nuclear DNA, mitochondrial protein lipoxidation, and lipofuscin accumulation compared to young animals (4 months old) in the liver. After 7 weeks of dietary treatment all those increases were totally or partially (lipofuscin) abolished by rapamycin, middle aged rapamycin-treated animals showing similar levels in those parameters to young animals. The decrease in mitochondrial ROS production was due to qualitative instead of quantitative changes in complex I. The decrease in mitochondrial protein lipoxidation was not due to decreases in the amount of highly oxidizable unsaturated fatty acids. Rapamycin also decreased the amount of RAPTOR (of mTOR complex) and increased the amounts of the PGC1-α and ATG13 proteins. The results are consistent with the possibility that rapamycin increases longevity in mice at least in part by lowering mitochondrial ROS production and increasing autophagy, decreasing the derived final forms of damage accumulated with age which are responsible for increased longevity. The decrease in lipofuscin accumulation induced by rapamycin adds to previous information suggesting that the increase in longevity induced by this drug can be due to a decrease in the rate of aging. © 2016 Elsevier Inc.

Má oclusão Classe I de Angle, com mordida aberta anterior, tratada com extração de dentes permanentes

A mordida aberta é uma anomalia com características distintas que, além da complexidade dos múltiplos fatores etiológicos, traz consequências estéticas e funcionais. Muitas alternativas têm sido utilizadas em seu tratamento, entre elas a grade palatina, forças ortopédicas, ajuste oclusal, camuflagem com ou sem exodontias, mini-implantes ou miniplacas e cirurgia ortognática. O diagnóstico preciso e a determinação da etiologia permitem estabelecer os objetivos e o plano de tratamento ideal para essa má oclusão. O presente relato descreve o tratamento de uma má oclusão Classe I de Angle, com padrão esquelético de Classe II e mordida aberta anterior, realizado em duas fases e que foi apresentado à diretoria do Board Brasileiro de Ortodontia e Ortopedia Facial (BBO), representando a categoria 2, como parte dos requisitos para a obtenção do título de Diplomado pelo BBO.

Oxidative DNA damage induced by S(IV) in the presence of Cu(II) and Cu(I) complexes

The DNA damage induced by S(IV) in the presence of some Cu(II) complexes in air saturated solution was investigated. The addition of S(IV) to an air saturated solution containing CuII GGA (GGA = glycylglycyl-L-alanine), CuII G3 (G3 = triglycine) or CuII G4 (G4 = tetraglycine) and Ni(II) traces, causes rapid formation of the respective Cu(III) complex, with simultaneous O2 uptake and S(IV) oxidation. SO3•- and HO• were detected by EPR-spin trapping experiments. The DNA strand breaks were attributed to the oxysulfur radicals formed. In the reduction of Cu(II)/BCA (BCA = 4,4' dicarboxy-2-2'-biquinoline) by S(IV), with CuI BCA complex formation, there is the possible formation of carbon centered radical of BCA or peroxyl radical (ROO•) capable of oxidizing DNA bases. The intensity of DNA damage in the presence of these Cu(II) complexes and S(IV) (10-300 µmol L-1) followed the order: CuII BCA ∼ CuII G4 ∼ Cu(II) (added as Cu(NO3)2) > CuII G3 ∼ CuII GGA. Specifically for CuII BCA the damage occurred even at lower S(IV) concentration (0.1 µmol L-1). For the Cu(II) complexes with glycylglycylhistidine, glycylhistidylglycine, glycylhistidyllysine and glycylglycyltyrosylarginine the Cu(III) formation and the DNA damage was not observed.

Reactivity, photolability, and computational studies of the ruthenium nitrosyl complex with a substituted cyclam fac-[Ru(NO)Cl(2)(kappa(3)N(4),N(8),N(11)(1-carboxypropyl)cyclam)]Cl center dot H(2)O

Chemical reactivity, photolability, and computational studies of the ruthenium nitrosyl complex with a substituted cyclam, fac-[Ru(NO)Cl(2)(kappa(3)N(4),N(8),N(11)(1-carboxypropyl)cyclam)]Cl center dot H(2)O ((1-carboxypropyl) cyclam = 3-(1,4,8,11-tetraazacyclotetradecan-1-yl) propionic acid)), (I) are described. Chloride ligands do not undergo aquation reactions (at 25 degrees C, pH 3). The rate of nitric oxide (NO) dissociation (k(obs-NO)) upon reduction of I is 2.8 s(-1) at 25 +/- 1 degrees C (in 0.5 mol L(-1) HCl), which is close to the highest value found for related complexes. The uncoordinated carboxyl of I has a pK(a) of similar to 3.3, which is close to that of the carboxyl of the non coordinated (1-carboxypropyl) cyclam (pK(a) = 3.4). Two additional pK(a) values were found for I at similar to 8.0 and similar to 11.5. Upon electrochemical reduction or under irradiation with light (lambda(irr) = 350 or 520 nm; pH 7.4), I releases NO in aqueous solution. The cyclam ring N bound to the carboxypropyl group is not coordinated, resulting in a fac configuration that affects the properties and chemical reactivities of I, especially as NO donor, compared with analogous trans complexes. Among the computational models tested, the B3LYP/ECP28MDF, cc-pVDZ resulted in smaller errors for the geometry of I. The computational data helped clarify the experimental acid-base equilibria and indicated the most favourable site for the second deprotonation, which follows that of the carboxyl group. Furthermore, it showed that by changing the pH it is possible to modulate the electron density of I with deprotonation. The calculated NO bond length and the Ru/NO charge ratio indicated that the predominant canonical structure is [Ru(III)NO], but the Ru-NO bond angles and bond index (b.i.) values were less clear; the angles suggested that [Ru(II)NO(+)] could contribute to the electronic structure of I and b.i. values indicated a contribution from [Ru(IV)NO(-)]. Considering that some experimental data are consistent with a [Ru(II)NO(+)] description, while others are in agreement with [Ru(III)NO], the best description for I would be a linear combination of the three canonical forms, with a higher weight for [Ru(II)NO(+)] and [Ru(III)NO].

An empirical evaluation of imputation accuracy for association statistics reveals increased type-I error rates in genome-wide associations

Background: Genome wide association studies (GWAS) are becoming the approach of choice to identify genetic determinants of complex phenotypes and common diseases. The astonishing amount of generated data and the use of distinct genotyping platforms with variable genomic coverage are still analytical challenges. Imputation algorithms combine directly genotyped markers information with haplotypic structure for the population of interest for the inference of a badly genotyped or missing marker and are considered a near zero cost approach to allow the comparison and combination of data generated in different studies. Several reports stated that imputed markers have an overall acceptable accuracy but no published report has performed a pair wise comparison of imputed and empiric association statistics of a complete set of GWAS markers. Results: In this report we identified a total of 73 imputed markers that yielded a nominally statistically significant association at P < 10(-5) for type 2 Diabetes Mellitus and compared them with results obtained based on empirical allelic frequencies. Interestingly, despite their overall high correlation, association statistics based on imputed frequencies were discordant in 35 of the 73 (47%) associated markers, considerably inflating the type I error rate of imputed markers. We comprehensively tested several quality thresholds, the haplotypic structure underlying imputed markers and the use of flanking markers as predictors of inaccurate association statistics derived from imputed markers. Conclusions: Our results suggest that association statistics from imputed markers showing specific MAF (Minor Allele Frequencies) range, located in weak linkage disequilibrium blocks or strongly deviating from local patterns of association are prone to have inflated false positive association signals. The present study highlights the potential of imputation procedures and proposes simple procedures for selecting the best imputed markers for follow-up genotyping studies.

Group I introns and associated homing endonuclease genes reveals a clinal structure for Porphyra spiralis var. amplifolia (Bangiales, Rhodophyta) along the Eastern coast of South America

Background: Group I introns are found in the nuclear small subunit ribosomal RNA gene (SSU rDNA) of some species of the genus Porphyra (Bangiales, Rhodophyta). Size polymorphisms in group I introns has been interpreted as the result of the degeneration of homing endonuclease genes (HEG) inserted in peripheral loops of intron paired elements. In this study, intron size polymorphisms were characterized for different Porphyra spiralis var. amplifolia (PSA) populations on the Southern Brazilian coast, and were used to infer genetic relationships and genetic structure of these PSA populations, in addition to cox2-3 and rbcL-S regions. Introns of different sizes were tested qualitatively for in vitro self-splicing. Results: Five intron size polymorphisms within 17 haplotypes were obtained from 80 individuals representing eight localities along the distribution of PSA in the Eastern coast of South America. In order to infer genetic structure and genetic relationships of PSA, these polymorphisms and haplotypes were used as markers for pairwise Fst analyses, Mantel's test and median joining network. The five cox2-3 haplotypes and the unique rbcL-S haplotype were used as markers for summary statistics, neutrality tests Tajima's D and Fu's Fs and for median joining network analyses. An event of demographic expansion from a population with low effective number, followed by a pattern of isolation by distance was obtained for PSA populations with the three analyses. In vitro experiments have shown that introns of different lengths were able to self-splice from pre-RNA transcripts. Conclusion: The findings indicated that degenerated HEGs are reminiscent of the presence of a full-length and functional HEG, once fixed for PSA populations. The cline of HEG degeneration determined the pattern of isolation by distance. Analyses with the other markers indicated an event of demographic expansion from a population with low effective number. The different degrees of degeneration of the HEG do not refrain intron self-splicing. To our knowledge, this was the first study to address intraspecific evolutionary history of a nuclear group I intron; to use nuclear, mitochondrial and chloroplast DNA for population level analyses of Porphyra; and intron size polymorphism as a marker for population genetics.

Characterization of subgraph relationships and distribution in complex networks

A network can be analyzed at different topological scales, ranging from single nodes to motifs, communities, up to the complete structure. We propose a novel approach which extends from single nodes to the whole network level by considering non-overlapping subgraphs (i.e. connected components) and their interrelationships and distribution through the network. Though such subgraphs can be completely general, our methodology focuses on the cases in which the nodes of these subgraphs share some special feature, such as being critical for the proper operation of the network. The methodology of subgraph characterization involves two main aspects: (i) the generation of histograms of subgraph sizes and distances between subgraphs and (ii) a merging algorithm, developed to assess the relevance of nodes outside subgraphs by progressively merging subgraphs until the whole network is covered. The latter procedure complements the histograms by taking into account the nodes lying between subgraphs, as well as the relevance of these nodes to the overall subgraph interconnectivity. Experiments were carried out using four types of network models and five instances of real-world networks, in order to illustrate how subgraph characterization can help complementing complex network-based studies.

Fluoxetine induces preventive and complex effects against colon cancer development in epithelial and stromal areas in rats

Fluoxetine (FIX) is a drug commonly used as antidepressant. However, its effects on tumorigenesis remain controversial. Aiming to evaluate the effects of FIX treatment on early malignant changes, we analyzed serotonin (5-HT) metabolism and recognition, aberrant crypt foci (ACF), proliferative process, microvessels, vascular endothelial growth factor (VEGF), and cyclooxygenase-2 (COX-2) expression in colon tissue. Male Wistar rats received a daily FLX-gavage (30 mg kg(-1)) and, a single dose of 1.2 dimethylhydrazine (DMH; i.p., 125 mg kg(-1)). After 6 weeks of FIX-treatment, our results revealed that FIX and nor-fluoxetine (N-FIX) are present in colon tissue, which was related to significant increase in serotonin (5-HT) levels (P < 0.05) possibly through a blockade in SERT mRNA (serotonin reuptake transporter; P < 0.05) resulting in lower 5-hydroxyindoleacetic acid (5-HIAA) levels (P < 0.01) and, 5-HT2C receptor mRNA expressions. FIX-treatment decreased dysplastic ACF development (P < 0.01) and proliferative process (P < 0.001) in epithelia. We observed a significant decrease in the development of malignant microvessels (P < 0.05), VEGF (P < 0.001), and COX-2 expression (P < 0.01). These findings suggest that FIX may have oncostatic effects on carcinogenic colon tissue, probably due to its modulatory activity on 5-HT metabolism and/or its ability to reduce colonic malignant events. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

The resolution of complex spectral patterns by cochlear implant and normal-hearing listeners

The differences in spectral shape resolution abilities among cochlear implant ~CI! listeners, and between CI and normal-hearing ~NH! listeners, when listening with the same number of channels ~12!, was investigated. In addition, the effect of the number of channels on spectral shape resolution was examined. The stimuli were rippled noise signals with various ripple frequency-spacings. An adaptive 4IFC procedure was used to determine the threshold for resolvable ripple spacing, which was the spacing at which an interchange in peak and valley positions could be discriminated. The results showed poorer spectral shape resolution in CI compared to NH listeners ~average thresholds of approximately 3000 and 400 Hz, respectively!, and wide variability among CI listeners ~range of approximately 800 to 8000 Hz!. There was a significant relationship between spectral shape resolution and vowel recognition. The spectral shape resolution thresholds of NH listeners increased as the number of channels increased from 1 to 16, while the CI listeners showed a performance plateau at 4–6 channels, which is consistent with previous results using speech recognition measures. These results indicate that this test may provide a measure of CI performance which is time efficient and non-linguistic, and therefore, if verified, may provide a useful contribution to the prediction of speech perception in adults and children who use CIs.

Papillomavirus virus-like particles can deliver defined CTL epitopes to the MHC class I pathway

To evaluate an antigen delivery system in which exogenous antigen can target the major histocompatibility complex (MHC) class I pathway, a single human papillomavirus (HPV) 16 E7 cytotoxic T lymphocyte (CTL) epitope and a single HIV gp160 CTL epitope were separately fused to the C-terminus or bovine papillomavirus 1 (BPV1) L1 sequence to form hybrid BPV1L1 VLPs. Mice immunized with these hybrid VLPs mounted strong CTL responses against the relevant target cells in the absence of any adjuvants. In addition, the CTL responses induced by immunization with BPV1L1/HPV16E7CTL VLPs protected mice against challenge with E7-transformed tumor cells. Furthermore, a high titer-specific antibody response against BPV1L1 VLPs was also induced, and this antiserum could inhibit papillomavirus-induced agglutination of mouse erythrocytes, suggesting that the antibody may recognize conformational determinates relevant to virus neutralization. These data demonstrate that hybrid BPV1L1 VLPs can be used as carriers to target antigenic epitopes to both the MHC class I and class II pathways, providing a promising strategy for the design of vaccines to prevent virus infection, with the potential to elicit therapeutic virus-specific CTL responses. (C) 1998 Academic Press.