In silico prediction of peptide-MHC binding affinity using SVRMHC

The binding between peptide epitopes and major histocompatibility complex (MHC) proteins is a major event in the cellular immune response. Accurate prediction of the binding between short peptides and class I or class II MHC molecules is an important task in immunoinformatics. SVRMHC which is a novel method to model peptide-MHC binding affinities based on support rector machine regression (SVR) is described in this chapter. SVRMHC is among a small handful of quantitative modeling methods that make predictions about precise binding affinities between a peptide and an MHC molecule. As a kernel-based learning method, SVRMHC has rendered models with demonstrated appealing performance in the practice of modeling peptide-MHC binding.

Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity

Background - MHC Class I molecules present antigenic peptides to cytotoxic T cells, which forms an integral part of the adaptive immune response. Peptides are bound within a groove formed by the MHC heavy chain. Previous approaches to MHC Class I-peptide binding prediction have largely concentrated on the peptide anchor residues located at the P2 and C-terminus positions. Results - A large dataset comprising MHC-peptide structural complexes was created by re-modelling pre-determined x-ray crystallographic structures. Static energetic analysis, following energy minimisation, was performed on the dataset in order to characterise interactions between bound peptides and the MHC Class I molecule, partitioning the interactions within the groove into van der Waals, electrostatic and total non-bonded energy contributions. Conclusion - The QSAR techniques of Genetic Function Approximation (GFA) and Genetic Partial Least Squares (G/PLS) algorithms were used to identify key interactions between the two molecules by comparing the calculated energy values with experimentally-determined BL50 data. Although the peptide termini binding interactions help ensure the stability of the MHC Class I-peptide complex, the central region of the peptide is also important in defining the specificity of the interaction. As thermodynamic studies indicate that peptide association and dissociation may be driven entropically, it may be necessary to incorporate entropic contributions into future calculations.

Large-scale molecular dynamics simulations of HLA-A*0201 complexed with a tumor-specific antigenic peptide:can the α3 and β2m domains be neglected?

Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histo-compatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVY-DGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide–MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent-accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the α1–α2 domain, the α3 and β2m domains cannot be neglected. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1803–1813, 2004

Visualisation of bioinformatics datasets

Analysing the molecular polymorphism and interactions of DNA, RNA and proteins is of fundamental importance in biology. Predicting functions of polymorphic molecules is important in order to design more effective medicines. Analysing major histocompatibility complex (MHC) polymorphism is important for mate choice, epitope-based vaccine design and transplantation rejection etc. Most of the existing exploratory approaches cannot analyse these datasets because of the large number of molecules with a high number of descriptors per molecule. This thesis develops novel methods for data projection in order to explore high dimensional biological dataset by visualising them in a low-dimensional space. With increasing dimensionality, some existing data visualisation methods such as generative topographic mapping (GTM) become computationally intractable. We propose variants of these methods, where we use log-transformations at certain steps of expectation maximisation (EM) based parameter learning process, to make them tractable for high-dimensional datasets. We demonstrate these proposed variants both for synthetic and electrostatic potential dataset of MHC class-I. We also propose to extend a latent trait model (LTM), suitable for visualising high dimensional discrete data, to simultaneously estimate feature saliency as an integrated part of the parameter learning process of a visualisation model. This LTM variant not only gives better visualisation by modifying the project map based on feature relevance, but also helps users to assess the significance of each feature. Another problem which is not addressed much in the literature is the visualisation of mixed-type data. We propose to combine GTM and LTM in a principled way where appropriate noise models are used for each type of data in order to visualise mixed-type data in a single plot. We call this model a generalised GTM (GGTM). We also propose to extend GGTM model to estimate feature saliencies while training a visualisation model and this is called GGTM with feature saliency (GGTM-FS). We demonstrate effectiveness of these proposed models both for synthetic and real datasets. We evaluate visualisation quality using quality metrics such as distance distortion measure and rank based measures: trustworthiness, continuity, mean relative rank errors with respect to data space and latent space. In cases where the labels are known we also use quality metrics of KL divergence and nearest neighbour classifications error in order to determine the separation between classes. We demonstrate the efficacy of these proposed models both for synthetic and real biological datasets with a main focus on the MHC class-I dataset.

Estudio en aves silvestres de canto ("Passeriformes") de las moléculas MHC y reclasificación molecular de la familia "Carduelini"

The Major Histocompatibility Complex (MHC) comprises the most polymorphic loci in animals. MHC plays an important role during the first steps of the immune response in vertebrates. In humans, MHC molecules (also named human leukocyte antigens, HLA) were initially regarded as class I or class II molecules. Each of them, presents to different T cells subsets. MHC class I molecules, are heterodimers in which the heavy chain (alpha) has three extracellular domains, two of which (alpha 1 and alpha 2) are polymorphic and conform the antigen recognition sites (ARS). The ARS is thought to be subjected to balancing selection for variability, which is the cause of the very high polymorphism of the MHC molecules. Different pathogenic epitopes would be the evolutionary force causing balancing selection. MHC class I genes have been completely sequenced (α1 and α2 protein domains) and thoroughly studied in Gallus gallus (chicken) as well as in mammals. In fact, the MHC locus was first defined in chicken, specifically in the highly consanguineous variety „Leghorn‟. It has been found that, in the case of chickens the MHC genetic region is considerably smaller than it is in mammals (remarkably shorter introns were found in chickens), and is organized quite differently. The noteworthy presence of short introns in chickens; supported the hypothesis that chicken‟s MHC represented a „minimal essential MHC‟. Until now, it has been assumed that chicken (order Galliformes) MHC was similar to all species included in the whole class Aves...

Doença enxerto-versus hospedeiro: determinação da expressão de genes e alterações de vias associadas ao sistema imunológico usando metodos da biofísica

O transplante de medula óssea (TMO) é um procedimento terapêutico importante em casos relacionados à pacientes com leucemia ou linfoma. Em decorrência desse processo, uma reação conhecida como doença enxerto-versus-hospedeiro (GVHD) pode ocorrer em pacientes susceptíveis como conseqüência da presença de células imunocompetentes do doador. Entretanto, não existe um modelo para descrever completamente as ações relacionadas ao mecanismo imunológico da GVHD desde a fase que inicializa a doença até a fase efetora. O Objetivo geral deste estudo é a investigação da resposta imunológica considerando-se o sistema HLA (antígenos leucocitários humano) em pacientes que desenvolveram a GVHD em decorrência do TMO. O National Cancer Institute (NCI) – Pathway interaction Database e Reactome foram usados como bases de dados com o objetivo de se estudar a expressão de genes e vias relacionados às Classes I e II do sistema HLA (antígenos leucocitários humano). O estudo considerou a mudança de expressão de genes relacionados às 17 vias do sistema imunológico com potencialidade para se expressar em pacientes que desenvolveram a GVHD associada à TMO. Dados referentes aos transcriptomas foram obtidos utilizando-se a plataforma GPL570 Affymetrix Genoma Humano U133 Plus. A atividade relativa foi usada para determinar as alterações das vias em amostras de GVHD em relação ao controle. As análises foram realizadas utilizando-se o software Via Complex e Bioconductor. Observou-se aumento significativo da expressão de genes ralacionados às vias do sistema imune adaptativo, antígenos associados às Classe I e II do HLA, fosforilação de CD3 e CD247, sinalização dos receptores de células T em CD4+ nativas e ativação de NF-kapa β nas células B. Também observou-se alterações significativas na mudança de expressão dos genes associados às vias relacionadas à super família de moléculas B7:CD28\CTLA-4 quando comparadas ao controle. Isso pode indicar a necessidade de geração de um segundo sinal co-estimulador em GVHD, acionado pelas moléculas dessa super família. O aumento da expressão do gene CD69 nas amostras experimentais caracteriza a ativação celular e, portanto, a sinalização de estímulos em GVHD. Os achados obtidos neste estudo contribuem para melhor elucidar o mecanismo imunopatogênico associado à GVHD. P

Scheduling of track inspection and maintenance activities in railroad networks

The U.S. railroad companies spend billions of dollars every year on railroad track maintenance in order to ensure safety and operational efficiency of their railroad networks. Besides maintenance costs, other costs such as train accident costs, train and shipment delay costs and rolling stock maintenance costs are also closely related to track maintenance activities. Optimizing the track maintenance process on the extensive railroad networks is a very complex problem with major cost implications. Currently, the decision making process for track maintenance planning is largely manual and primarily relies on the knowledge and judgment of experts. There is considerable potential to improve the process by using operations research techniques to develop solutions to the optimization problems on track maintenance. In this dissertation study, we propose a range of mathematical models and solution algorithms for three network-level scheduling problems on track maintenance: track inspection scheduling problem (TISP), production team scheduling problem (PTSP) and job-to-project clustering problem (JTPCP). TISP involves a set of inspection teams which travel over the railroad network to identify track defects. It is a large-scale routing and scheduling problem where thousands of tasks are to be scheduled subject to many difficult side constraints such as periodicity constraints and discrete working time constraints. A vehicle routing problem formulation was proposed for TISP, and a customized heuristic algorithm was developed to solve the model. The algorithm iteratively applies a constructive heuristic and a local search algorithm in an incremental scheduling horizon framework. The proposed model and algorithm have been adopted by a Class I railroad in its decision making process. Real-world case studies show the proposed approach outperforms the manual approach in short-term scheduling and can be used to conduct long-term what-if analyses to yield managerial insights. PTSP schedules capital track maintenance projects, which are the largest track maintenance activities and account for the majority of railroad capital spending. A time-space network model was proposed to formulate PTSP. More than ten types of side constraints were considered in the model, including very complex constraints such as mutual exclusion constraints and consecution constraints. A multiple neighborhood search algorithm, including a decomposition and restriction search and a block-interchange search, was developed to solve the model. Various performance enhancement techniques, such as data reduction, augmented cost function and subproblem prioritization, were developed to improve the algorithm. The proposed approach has been adopted by a Class I railroad for two years. Our numerical results show the model solutions are able to satisfy all hard constraints and most soft constraints. Compared with the existing manual procedure, the proposed approach is able to bring significant cost savings and operational efficiency improvement. JTPCP is an intermediate problem between TISP and PTSP. It focuses on clustering thousands of capital track maintenance jobs (based on the defects identified in track inspection) into projects so that the projects can be scheduled in PTSP. A vehicle routing problem based model and a multiple-step heuristic algorithm were developed to solve this problem. Various side constraints such as mutual exclusion constraints and rounding constraints were considered. The proposed approach has been applied in practice and has shown good performance in both solution quality and efficiency.

Evaluación de la posición hioidea en radiografías de pacientes clase I, II y III esqueletales entre 9 y 18 años de edad registrados en el período 2012 y 2015, de un Centro Radiológico de la Ciudad de Cuenca

Objetivo: El objetivo del presente estudio descriptivo, fue evaluar la posición del hueso hioides en los diferentes patrones esqueletales de Clase I, II y III mediante el trazado cefalométrico del triángulo hioideo propuesto por Bibby y Preston, estableciendo diferencias entre cada clase esqueletal. Materiales y métodos: La muestra consistió en 161 radiografías cefálicas laterales digitales, correspondientes a individuos de ambos sexos (75 hombres y 86 mujeres), entre edades de 9 y 18 años, las mismas que fueron divididas en tres subgrupos (Clase I, clase II y clase III) de acuerdo a los ángulos ANB y APDI. Se determinó la posición anteroposterior, vertical y angular del hueso hioides mediante el trazado cefalométrico del triángulo hioideo siendo el mentón, la tercera vértebra cervical y el hueso hioides las estructuras anatómicas utilizadas para el trazado del mismo. Se obtuvieron medidas estándar para cada clase esqueletal. Resultados: Se observaron diferencias estadísticamente significativas en la medida de H-Rgn entre clase I y II y entre clase II y III (p<0,005). El valor del ángulo del plano hioidal presentó diferencias estadísticamente significativas entre clase I y III y entre clase II y III (p<0,005). Se evidenciaron diferencias estadísticamente significativas entre hombres y mujeres con clase I esqueletal en la medida H-Rgn (p<0,005). Conclusiones: La posición del hueso hioides varía en los diferentes patrones esqueletales. Sin embargo, su posición en relación a la columna cervical presenta menos variabilidad que su relación con la mandíbula

Relación cráneo cervical en pacientes clase I, II, y III esqueletal entre 9 y 18 años de edad, mediante radiografías cefálicas laterales de un centro radiológico de la ciudad de cuenca.

El objetivo de este estudio fue determinar la relación cráneo cervical en pacientes clases I, II y III esqueletal entre 9 y 18 años de edad, mediante radiografías cefálicas laterales de un centro radiológico de la ciudad de Cuenca, utilizando el análisis cráneo cervical propuesto por Rocabado. Materiales y métodos: Fueron analizadas 161 radiografías cefálicas laterales digitales, de ambos sexos, con edad promedio de 12.3 años (DE± 2.4). Se incluyeron radiografías de individuos con dentición mixta y permanente, sin tratamiento ortodóncico y en donde se observe hasta la sexta vértebra cervical. Fueron excluidas las radiografías de pacientes con mordida abierta, traumatismos maxilofaciales y radiografías de mala calidad. Las telerradiografías fueron analizadas mediante el programa cefalométrico Nemoceph NX, donde se determinó el patrón esqueletal mediante los ángulos SNA, SNB, ANB y APDI. La evaluación de la postura cervical, se realizó mediante el análisis cráneo cervical propuesto por Rocabado. Se obtuvo el índice de concordancia (ICC=0.94). Mediante estadística descriptiva se analizaron las relaciones entre variables usando la prueba de Chi cuadrado y T de Student. Resultados: Se encontró mayor rotación posterior de cráneo en clase I y II esqueletal, encontrándose diferencias estadísticamente significativas respecto al ángulo cráneo vertebral entre hombres y mujeres en individuos clase II esqueletal. Las mujeres presentaron mayor rotación posterior de cráneo a diferencia de los hombres. (p=0.004). En clase III se encontró una relación normal. El espacio suboccipital en las tres clases esqueletales se presento con normalidad. No se encontró diferencia significativa respecto a la edad. Conclusiones: La relación cráneo cervical se presenta con una tendencia a la rotación posterior de cráneo, influida fuertemente por el sexo del individuo. El espacio suboccipital es normal en clase I y II esqueletal y con tendencia al aumento en clase III.

Homotheties and topology of tangent sphere bundles

We prove a Theorem on homotheties between two given tangent sphere bundles SrM of a Riemannian manifold (M,g) of dim ≥ 3, assuming different variable radius functions r and weighted Sasaki metrics induced by the conformal class of g. New examples are shown of manifolds with constant positive or with constant negative scalar curvature which are not Einstein. Recalling results on the associated almost complex structure I^G and symplectic structure ω^G on the manifold TM , generalizing the well-known structure of Sasaki by admitting weights and connections with torsion, we compute the Chern and the Stiefel-Whitney characteristic classes of the manifolds TM and SrM.

Clinical Consequences Of Tityus Bahiensis And Tityus Serrulatus Scorpion Stings In The Region Of Campinas, Southeastern Brazil.

Scorpion stings account for most envenomations by venomous animals in Brazil. A retrospective study (1994-2011) of the clinical consequences of Tityus scorpion stings in 1327 patients treated at a university hospital in Campinas, southeastern Brazil, is reported. The clinical classification, based on outcome, was: dry sting (no envenoming), class I (only local manifestations), class II (systemic manifestations), class III (life-threatening manifestations, such as shock and/or cardiac failure requiring inotropic/vasopressor agents, and/or respiratory failure), and fatal. The median patient age was 27 years (interquartile interval = 15-42 years). Scorpions were brought for identification in 47.2% of cases (Tityus bahiensis 27.7%; Tityus serrulatus 19.5%). Sting severity was classified and each accounted for the following percentage of cases: dry stings - 3.4%, class I - 79.6%, class II - 15.1%, class III - 1.8% and fatal - 0.1%. Pain was the primary local manifestation (95.5%). Systemic manifestations such as vomiting, agitation, sweating, dyspnea, bradycardia, tachycardia, tachypnea, somnolence/lethargy, cutaneous paleness, hypothermia and hypotension were detected in class II or class III + fatal groups, but were significantly more frequent in the latter group. Class III and fatal cases occurred only in children <15 years old, with scorpions being identified in 13/25 cases (T. serrulatus, n = 12; T. bahiensis, n = 1). Laboratory blood abnormalities (hyperglycemia, hypokalemia, leukocytosis, elevations in serum total CK, CK-MB and troponin T, bicarbonate consumption and an increase in base deficit and blood lactate), electrocardiographic changes (ST segment) and echocardiographic alterations (ventricular ejected fraction <54%) were frequently detected in class III patients. Seventeen patients developed pulmonary edema, 16 had cardiac failure and seven had cardiogenic shock. These results indicate that most scorpion stings involved only local manifestations, mainly pain; the greatest severity was associated with stings by T. serrulatus and in children <15 years old.

The Inadequacy Of The Y-axis Of Growth (sngn) For The Vertical Pattern Assessment In Patients With Sagittal Discrepancies.

The aim of this cephalometric study was to evaluate the influence of the sagittal skeletal pattern on the 'Y-axis of growth' measurement in patients with different malocclusions. Lateral head films from 59 patients (mean age 16y 7m, ranging from 11 to 25 years) were selected after a subjective analysis of 1630 cases. Sample was grouped as follows: Group 1 - class I facial pattern; group 2 - class II facial pattern; and Group 3 - class III facial pattern. Two angular measurements, SNGoGn and SNGn, were taken in order to determine skeletal vertical facial pattern. A logistic regression with errors distributed according to a binomial distribution was used to test the influence of the sagittal relationship (Class I, II, III facial patterns) on vertical diagnostic measurement congruence (SNGoGn and SNGn). RESULTS show that the probability of congruence between the patterns SNGn and SNGoGn was relatively high (70%) for group 1, but for groups II (46%) and III (37%) this congruence was relatively low. The use of SNGn appears to be inappropriate to determine the vertical facial skeletal pattern of patients, due to Gn point shifting throughout sagittal discrepancies. Clinical Significance: Facial pattern determined by SNGn must be considered carefully, especially when severe sagittal discrepancies are present.

A influência do protocolo de extração de dois pré-molares superiores na estabilidade oclusal do tratamento da Classe II

OBJETIVO: com o propósito de avaliar a influência da extração de dois pré-molares superiores na estabilidade oclusal do tratamento da má oclusão de Classe II completa, foi realizada uma comparação com o protocolo de tratamento sem extrações. MÉTODOS: selecionou-se, a partir das documentações do arquivo da Disciplina de Ortodontia da Faculdade de Odontologia de Bauru, uma amostra composta pelas documentações de 59 pacientes com má oclusão de Classe II completa. Em seguida, dividiu-se essa amostra em dois grupos, apresentando as seguintes características: Grupo 1, constituído por 29 pacientes, tratados sem extrações; e Grupo 2, composto por 30 pacientes, tratados com extrações de dois pré-molares superiores. Os modelos ao início do tratamento, ao final do tratamento e em um período mínimo de 2,4 anos após o tratamento foram medidos e avaliados por meio dos índices oclusais IPT e PAR. As condições oclusais ao final do tratamento e no estágio pós-tratamento, o percentual de recidiva e as alterações oclusais pós-tratamento foram comparados por meio do teste t. RESULTADOS: os resultados demonstraram que os protocolos de tratamento sem extração e com extrações de dois pré-molares superiores não apresentaram, em nenhuma das variáveis avaliadas, diferenças estatisticamente significativas em relação à estabilidade oclusal do tratamento da má oclusão de Classe II completa. CONCLUSÃO: a extração de dois pré-molares superiores no tratamento da má oclusão de Classe II completa não influenciou a estabilidade dos resultados oclusais alcançados ao final da correção ortodôntica. Portanto, terminar o tratamento com uma relação molar em Classe II ou em Classe I proporciona estabilidade semelhante.

Comparative study of dental arch width in plaster models, photocopies and digitized images

The aim of this study was to comparatively assess dental arch width, in the canine and molar regions, by means of direct measurements from plaster models, photocopies and digitized images of the models. The sample consisted of 130 pairs of plaster models, photocopies and digitized images of the models of white patients (n = 65), both genders, with Class I and Class II Division 1 malocclusions, treated by standard Edgewise mechanics and extraction of the four first premolars. Maxillary and mandibular intercanine and intermolar widths were measured by a calibrated examiner, prior to and after orthodontic treatment, using the three modes of reproduction of the dental arches. Dispersion of the data relative to pre- and posttreatment intra-arch linear measurements (mm) was represented as box plots. The three measuring methods were compared by one-way ANOVA for repeated measurements (α = 0.05). Initial / final mean values varied as follows: 33.94 to 34.29 mm / 34.49 to 34.66 mm (maxillary intercanine width); 26.23 to 26.26 mm / 26.77 to 26.84 mm (mandibular intercanine width); 49.55 to 49.66 mm / 47.28 to 47.45 mm (maxillary intermolar width) and 43.28 to 43.41 mm / 40.29 to 40.46 mm (mandibular intermolar width). There were no statistically significant differences between mean dental arch widths estimated by the three studied methods, prior to and after orthodontic treatment. It may be concluded that photocopies and digitized images of the plaster models provided reliable reproductions of the dental arches for obtaining transversal intra-arch measurements.

Confiabilidade da eficiência mastigatória com beads e correlação com a atividade muscular

TEMA: a avaliação da eficiência mastigatória pela análise colorimétrica com beads, pode ser um método promissor, mas não há relatos sobre a sua confiabilidade. OBJETIVO: investigar a confiabiabilidade das beads para teste de eficiência mastigatória e a correlação com a atividade eletromiográfica dos músculos masseter e temporal anterior. MÉTODO: participaram dezenove sujeitos adultos jovens, nove do gênero masculino e dez do feminino com idades entre dezoito e vinte-oito anos, com dentição completa, sem histórico de desordem temporomandibular, trauma, cirurgia na região de cabeça e pescoço, tratamento ortodôntico ou fonoaudiológico. O teste de eficiência mastigatória foi realizado com beads nas condições: mastigação habitual, mastigação unilateral direita e esquerda, com duração de 20 segundos. Simultaneamente, foi realizada a eletromiografia. A atividade em máxima intercuspidação habitual dos dentes também foi registrada. A quantidade de fucsina liberada após a mastigação foi medida usando o espectrofotômetro Beckman DU-7 UV-Visible (Beckman Inc., Palo Alto, CA, USA). RESULTADOS: houve alta confiabilidade do teste de eficiência mastigatória (r = 0,86, p < 0,01) e correlação significante com a atividade eletromiográfica (r = 0,76, p < 0,01). Também houve correlações positivas quando as provas foram analisadas separadamente. CONCLUSÃO: o teste de eficiência mastigatória realizado com beads mostrou-se um método confiável e correlacionado positivamente à atividade eletromiográfica dos músculos temporal anterior e músculos masseter.