Structure-based drug design studies of the interactions ofent-kaurane diterpenes derived from Wedelia paludosa with the Plasmodium falciparumsarco/endoplasmic reticulum Ca2+-ATPase PfATP6

Malaria is responsible for more deaths around the world than any other parasitic disease. Due to the emergence of strains that are resistant to the current chemotherapeutic antimalarial arsenal, the search for new antimalarial drugs remains urgent though hampered by a lack of knowledge regarding the molecular mechanisms of artemisinin resistance. Semisynthetic compounds derived from diterpenes from the medicinal plant Wedelia paludosawere tested in silico against the Plasmodium falciparumCa2+-ATPase, PfATP6. This protein was constructed by comparative modelling using the three-dimensional structure of a homologous protein, 1IWO, as a scaffold. Compound 21 showed the best docking scores, indicating a better interaction with PfATP6 than that of thapsigargin, the natural inhibitor. Inhibition of PfATP6 by diterpene compounds could promote a change in calcium homeostasis, leading to parasite death. These data suggest PfATP6 as a potential target for the antimalarial ent-kaurane diterpenes.

Computational drug design strategies applied to the modelling of human immunodeficiency virus-1 reverse transcriptase inhibitors

Reverse transcriptase (RT) is a multifunctional enzyme in the human immunodeficiency virus (HIV)-1 life cycle and represents a primary target for drug discovery efforts against HIV-1 infection. Two classes of RT inhibitors, the nucleoside RT inhibitors (NRTIs) and the nonnucleoside transcriptase inhibitors are prominently used in the highly active antiretroviral therapy in combination with other anti-HIV drugs. However, the rapid emergence of drug-resistant viral strains has limited the successful rate of the anti-HIV agents. Computational methods are a significant part of the drug design process and indispensable to study drug resistance. In this review, recent advances in computer-aided drug design for the rational design of new compounds against HIV-1 RT using methods such as molecular docking, molecular dynamics, free energy calculations, quantitative structure-activity relationships, pharmacophore modelling and absorption, distribution, metabolism, excretion and toxicity prediction are discussed. Successful applications of these methodologies are also highlighted.

Comparative Study of a Real-Time PCR Assay Targeting senX3-regX3 versus Other Molecular Strategies Commonly Used in the Diagnosis of Tuberculosis.

BACKGROUND Nucleic acid amplification tests are increasingly used for the rapid diagnosis of tuberculosis. We undertook a comparative study of the efficiency and diagnostic yield of a real-time PCR senX3-regX3 based assay versus the classical IS6110 target and the new commercial methods. METHODS This single-blind prospective comparative study included 145 consecutive samples: 76 from patients with culture-confirmed tuberculosis (86.8% pulmonary and 13.2% extrapulmonary tuberculosis: 48.7% smear-positive and 51.3% smear-negative) and 69 control samples (24 from patients diagnosed with non-tuberculous mycobacteria infections and 45 from patients with suspected tuberculosis which was eventually ruled out). All samples were tested by two CE-marked assays (Xpert®MTB/RIF and AnyplexTM plus MTB/NTM) and two in-house assays targeting senX3-regX3 and the IS6110 gene. RESULTS The detection limit ranged from 1.00E+01 fg for Anyplex, senX3-regX3 and IS6110 to 1.00E+04 fg for Xpert. All three Xpert, senX3-regX3 and IS6110 assays detected all 37 smear-positive cases. Conversely, Anyplex was positive in 34 (91.9%) smear-positive cases. In patients with smear-negative tuberculosis, differences were observed between the assays; Xpert detected 22 (56.41%) of the 39 smear-negative samples, Anyplex 24 (61.53%), senX3-regX3 28 (71.79%) and IS6110 35 (89.74%). Xpert and senX3-regX3 were negative in all control samples; however, the false positive rate was 8.7% and 13% for Anyplex and IS6110, respectively. The overall sensitivity was 77.6%, 85.7%, 77.3% and 94.7% and the specificity was 100%, 100%, 90.8% and 87.0% for the Xpert, senX3-regX3, Anyplex and IS6110 assays, respectively. CONCLUSION Real-time PCR assays targeting IS6110 lack the desired specificity. The Xpert MTB/RIF and in-house senX3-regX3 assays are both sensitive and specific for the detection of MTBC in both pulmonary and extrapulmonary samples. Therefore, the real time PCR senX3-regX3 based assay could be a useful and complementary tool in the diagnosis of tuberculosis.

Design of potent inhibitors of human RAD51 recombinase based on BRC motifs of BRCA2 protein: modeling and experimental validation of a chimera peptide.

We have previously shown that a 28-amino acid peptide derived from the BRC4 motif of BRCA2 tumor suppressor inhibits selectively human RAD51 recombinase (HsRad51). With the aim of designing better inhibitors for cancer treatment, we combined an in silico docking approach with in vitro biochemical testing to construct a highly efficient chimera peptide from eight existing human BRC motifs. We built a molecular model of all BRC motifs complexed with HsRad51 based on the crystal structure of the BRC4 motif-HsRad51 complex, computed the interaction energy of each residue in each BRC motif, and selected the best amino acid residue at each binding position. This analysis enabled us to propose four amino acid substitutions in the BRC4 motif. Three of these increased the inhibitory effect in vitro, and this effect was found to be additive. We thus obtained a peptide that is about 10 times more efficient in inhibiting HsRad51-ssDNA complex formation than the original peptide.

Social meets molecular: Combining phylogenetic and latent class analyses to understand HIV-1 transmission in Switzerland.

Switzerland has a complex human immunodeficiency virus (HIV) epidemic involving several populations. We examined transmission of HIV type 1 (HIV-1) in a national cohort study. Latent class analysis was used to identify socioeconomic and behavioral groups among 6,027 patients enrolled in the Swiss HIV Cohort Study between 2000 and 2011. Phylogenetic analysis of sequence data, available for 4,013 patients, was used to identify transmission clusters. Concordance between sociobehavioral groups and transmission clusters was assessed in correlation and multiple correspondence analyses. A total of 2,696 patients were infected with subtype B, 203 with subtype C, 196 with subtype A, and 733 with recombinant subtypes (mainly CRF02_AG and CRF01_AE). Latent class analysis identified 8 patient groups. Most transmission clusters of subtype B were shared between groups of gay men (groups 1-3) or between the heterosexual groups "heterosexual people of lower socioeconomic position" (group 4) and "injection drug users" (group 8). Clusters linking homosexual and heterosexual groups were associated with "older heterosexual and gay people on welfare" (group 5). "Migrant women in heterosexual partnerships" (group 6) and "heterosexual migrants on welfare" (group 7) shared non-B clusters with groups 4 and 5. Combining approaches from social and molecular epidemiology can provide insights into HIV-1 transmission and inform the design of prevention strategies.

Multifunctional ligands for application in Alzheimer’s Disease: molecular modelling and biological evaluation

Projecte de recerca elaborat a partir d’una estada a la University of British Columbia, Canadà, entre 2010 i 2012 La malaltia d'Alzheimer (MA) representa avui la forma més comuna de demència en la població envellida. Malgrat fa 100 anys que va ser descoberta, encara avui no existeix cap tractament preventiu i/o curatiu ni cap agent de diagnòstic que permeti valorar quantitativament l'evolució d'aquesta malaltia. L'objectiu en el que s'emmarca aquest treball és contribuir a aportar solucions al problema de la manca d'agents terapèutics i de diagnosi, unívocs i rigorosos, per a la MA. Des del camp de la química bioinorgànica és fàcil fixar-se en l'excessiva concentració d'ions Zn(II) i Cu(II) en els cervells de malalts de MA, plantejar-se la seva utilització com a dianes terapèutica i, en conseqüència, cercar agents quelants que evitin la formació de plaques senils o contribueixin a la seva dissolució. Si bé aquest va ser el punt de partida d’aquest projecte, els múltiples factors implicats en la patogènesi de la MA fan que el clàssic paradigma d’ ¨una molècula, una diana¨ limiti la capacitat de la molècula de combatre aquesta malaltia tan complexa. Per tant, un esforç considerable s’ha dedicat al disseny d’agentsmultifuncionals que combatin els múltiples factors que caracteritzen el desenvolupament de la MA. En el present treball s’han dissenyat agents multifuncionals inspirats en dos esquelets moleculars ben establers i coneguts en el camp de la química medicinal: la tioflavina-T (ThT) i la deferiprona (DFP). La utilització de tècniques in silico que inclouen càlculs farmacocinètics i modelatge molecular ha estat un procés cabdal per a l’avaluació dels millors candidats en base als següents requeriments: (a) compliment de determinades propietats farmacocinètiques que estableixin el seu possible ús com a fàrmac (b) hidrofobicitat adequada per travessar la BBB i (c) interacció amb el pèptid Aen solució.

Molecular modeling study of the binding and signaling properties of the TCRpMHC complex

The activation of the specific immune response against tumor cells is based on the recognition by the CD8+ Cytotoxic Τ Lymphocytes (CTL), of antigenic peptides (p) presented at the surface of the cell by the class I major histocompatibility complex (MHC). The ability of the so-called T-Cell Receptors (TCR) to discriminate between self and non-self peptides constitutes the most important specific control mechanism against infected cells. The TCR/pMHC interaction has been the subject of much attention in cancer therapy since the design of the adoptive transfer approach, in which Τ lymphocytes presenting an interesting response against tumor cells are extracted from the patient, expanded in vitro, and reinfused after immunodepletion, possibly leading to cancer regression. In the last decade, major progress has been achieved by the introduction of engineered lypmhocytes. In the meantime, the understanding of the molecular aspects of the TCRpMHC interaction has become essential to guide in vitro and in vivo studies. In 1996, the determination of the first structure of a TCRpMHC complex by X-ray crystallography revealed the molecular basis of the interaction. Since then, molecular modeling techniques have taken advantage of crystal structures to study the conformational space of the complex, and understand the specificity of the recognition of the pMHC by the TCR. In the meantime, experimental techniques used to determine the sequences of TCR that bind to a pMHC complex have been used intensively, leading to the collection of large repertoires of TCR sequences that are specific for a given pMHC. There is a growing need for computational approaches capable of predicting the molecular interactions that occur upon TCR/pMHC binding without relying on the time consuming resolution of a crystal structure. This work presents new approaches to analyze the molecular principles that govern the recognition of the pMHC by the TCR and the subsequent activation of the T-cell. We first introduce TCRep 3D, a new method to model and study the structural properties of TCR repertoires, based on homology and ab initio modeling. We discuss the methodology in details, and demonstrate that it outperforms state of the art modeling methods in predicting relevant TCR conformations. Two successful applications of TCRep 3D that supported experimental studies on TCR repertoires are presented. Second, we present a rigid body study of TCRpMHC complexes that gives a fair insight on the TCR approach towards pMHC. We show that the binding mode of the TCR is correctly described by long-distance interactions. Finally, the last section is dedicated to a detailed analysis of an experimental hydrogen exchange study, which suggests that some regions of the constant domain of the TCR are subject to conformational changes upon binding to the pMHC. We propose a hypothesis of the structural signaling of TCR molecules leading to the activation of the T-cell. It is based on the analysis of correlated motions in the TCRpMHC structure. - L'activation de la réponse immunitaire spécifique dirigée contre les cellules tumorales est basée sur la reconnaissance par les Lymphocytes Τ Cytotoxiques (CTL), d'un peptide antigénique (p) présenté à la suface de la cellule par le complexe majeur d'histocompatibilité de classe I (MHC). La capacité des récepteurs des lymphocytes (TCR) à distinguer les peptides endogènes des peptides étrangers constitue le mécanisme de contrôle le plus important dirigé contre les cellules infectées. L'interaction entre le TCR et le pMHC est le sujet de beaucoup d'attention dans la thérapie du cancer, depuis la conception de la méthode de transfer adoptif: les lymphocytes capables d'une réponse importante contre les cellules tumorales sont extraits du patient, amplifiés in vitro, et réintroduits après immunosuppression. Il peut en résulter une régression du cancer. Ces dix dernières années, d'importants progrès ont été réalisés grâce à l'introduction de lymphocytes modifiés par génie génétique. En parallèle, la compréhension du TCRpMHC au niveau moléculaire est donc devenue essentielle pour soutenir les études in vitro et in vivo. En 1996, l'obtention de la première structure du complexe TCRpMHC à l'aide de la cristallographie par rayons X a révélé les bases moléculaires de l'interaction. Depuis lors, les techniques de modélisation moléculaire ont exploité les structures expérimentales pour comprendre la spécificité de la reconnaissance du pMHC par le TCR. Dans le même temps, de nouvelles techniques expérimentales permettant de déterminer la séquence de TCR spécifiques envers un pMHC donné, ont été largement exploitées. Ainsi, d'importants répertoires de TCR sont devenus disponibles, et il est plus que jamais nécessaire de développer des approches informatiques capables de prédire les interactions moléculaires qui ont lieu lors de la liaison du TCR au pMHC, et ce sans dépendre systématiquement de la résolution d'une structure cristalline. Ce mémoire présente une nouvelle approche pour analyser les principes moléculaires régissant la reconnaissance du pMHC par le TCR, et l'activation du lymphocyte qui en résulte. Dans un premier temps, nous présentons TCRep 3D, une nouvelle méthode basée sur les modélisations par homologie et ab initio, pour l'étude de propriétés structurales des répertoires de TCR. Le procédé est discuté en détails et comparé à des approches standard. Nous démontrons ainsi que TCRep 3D est le plus performant pour prédire des conformations pertinentes du TCR. Deux applications à des études expérimentales des répertoires TCR sont ensuite présentées. Dans la seconde partie de ce travail nous présentons une étude de complexes TCRpMHC qui donne un aperçu intéressant du mécanisme d'approche du pMHC par le TCR. Finalement, la dernière section se concentre sur l'analyse détaillée d'une étude expérimentale basée sur les échanges deuterium/hydrogène, dont les résultats révèlent que certaines régions clés du domaine constant du TCR sont sujettes à un changement conformationnel lors de la liaison au pMHC. Nous proposons une hypothèse pour la signalisation structurelle des TCR, menant à l'activation du lymphocyte. Celle-ci est basée sur l'analyse des mouvements corrélés observés dans la structure du TCRpMHC.

iLOGP: A Simple, Robust, and Efficient Description of n-Octanol/Water Partition Coefficient for Drug Design Using the GB/SA Approach.

The n-octanol/water partition coefficient (log Po/w) is a key physicochemical parameter for drug discovery, design, and development. Here, we present a physics-based approach that shows a strong linear correlation between the computed solvation free energy in implicit solvents and the experimental log Po/w on a cleansed data set of more than 17,500 molecules. After internal validation by five-fold cross-validation and data randomization, the predictive power of the most interesting multiple linear model, based on two GB/SA parameters solely, was tested on two different external sets of molecules. On the Martel druglike test set, the predictive power of the best model (N = 706, r = 0.64, MAE = 1.18, and RMSE = 1.40) is similar to six well-established empirical methods. On the 17-drug test set, our model outperformed all compared empirical methodologies (N = 17, r = 0.94, MAE = 0.38, and RMSE = 0.52). The physical basis of our original GB/SA approach together with its predictive capacity, computational efficiency (1 to 2 s per molecule), and tridimensional molecular graphics capability lay the foundations for a promising predictor, the implicit log P method (iLOGP), to complement the portfolio of drug design tools developed and provided by the SIB Swiss Institute of Bioinformatics.

Clinicopathologic and molecular analysis of a choroidal pigmented schwannoma in the context of a PTEN hamartoma tumor syndrome.

PURPOSE: To report the first case of choroidal schwannoma in a patient affected by PTEN hamartoma tumor syndrome (PHTS) and investigate the molecular involvement of the phosphatase and tensin homolog (PTEN) and neurofibromin 2 (NF2) genes in this rare intraocular tumor. DESIGN: Observational case report. PARTICIPANT: A 10-year-old girl diagnosed with PHTS. METHODS: The enucleated specimen underwent histologic, immunohistochemical, and transmission electronic microscopy. The expression of PTEN and NF2 and their protein products were evaluated by reverse transcription-polymerase chain reaction and immunohistochemistry. Somatic mutations of PTEN and NF2, as well as allelic loss, were investigated by direct sequencing of DNA extracted from the tumor. PTEN epigenetic silencing was investigated by pyrosequencing. MAIN OUTCOME MEASURES: Histopathologic and molecular characterization of a choroidal pigmented schwannoma. RESULTS: Histopathologic, immunohistochemical, and electron microscopic analysis demonstrated features consistent with a pigmented cellular schwannoma of the choroid. We found no loss of heterozygosity at the genomic level for the PTEN germline mutation and no promoter hypermethylation or other somatic intragenic mutations. However, we observed an approximate 40% reduction of PTEN expression at both the mRNA and the protein level, indicating that the tumor was nonetheless functionally deficient for PTEN. Although DNA sequencing of NF2 failed to identify any pathologic variants, its expression was abolished within the tumor. CONCLUSIONS: We report the first description of a pigmented choroidal schwannoma in the context of a PHTS. This rare tumor showed a unique combination of reduction of PTEN and absence of NF2 expression.

Design of new promoters and of a dual-bioreporter based on cross-activation by the two regulatory proteins XylR and HbpR.

The HbpR protein is the sigma54-dependent transcription activator for 2-hydroxybiphenyl degradation in Pseudomonas azelaica. The ability of HbpR and XylR, which share 35% amino acid sequence identity, to cross-activate the PhbpC and Pu promoters was investigated by determining HbpR- or XylR-mediated luciferase expression and by DNA binding assays. XylR measurably activated the PhbpC promoter in the presence of the effector m-xylene, both in Escherichia coli and Pseudomonas putida. HbpR weakly stimulated the Pu promoter in E. coli but not in P. azelaica. Poor HbpR-dependent activation from Pu was caused by a weak binding to the operator region. To create promoters efficiently activated by both regulators, the HbpR binding sites on PhbpC were gradually changed into the XylR binding sites of Pu by site-directed mutagenesis. Inducible luciferase expression from mutated promoters was tested in E. coli on a two plasmid system, and from mono copy gene fusions in P. azelaica and P. putida. Some mutants were efficiently activated by both HbpR and XylR, showing that promoters can be created which are permissive for both regulators. Others achieved a higher XylR-dependent transcription than from Pu itself. Mutants were also obtained which displayed a tenfold lower uninduced expression level by HbpR than the wild-type PhbpC, while keeping the same maximal induction level. On the basis of these results, a dual-responsive bioreporter strain of P. azelaica was created, containing both XylR and HbpR, and activating luciferase expression from the same single promoter independently with m-xylene and 2-hydroxybiphenyl.

Homeostatic and circadian contribution to EEG and molecular state variables of sleep regulation.

STUDY OBJECTIVES: Besides their well-established role in circadian rhythms, our findings that the forebrain expression of the clock-genes Per2 and Dbp increases and decreases, respectively, in relation to time spent awake suggest they also play a role in the homeostatic aspect of sleep regulation. Here, we determined whether time of day modulates the effects of elevated sleep pressure on clock-gene expression. Time of day effects were assessed also for recognized electrophysiological (EEG delta power) and molecular (Homer1a) markers of sleep homeostasis. DESIGN: EEG and qPCR data were obtained for baseline and recovery from 6-h sleep deprivation starting at ZT0, -6, -12, or -18. SETTING: Mouse sleep laboratory. PARTICIPANTS: Male mice. INTERVENTIONS: Sleep deprivation. RESULTS: The sleep-deprivation induced changes in Per2 and Dbp expression importantly varied with time of day, such that Per2 could even decrease during sleep deprivations occurring at the decreasing phase in baseline. Dbp showed similar, albeit opposite dynamics. These unexpected results could be reliably predicted assuming that these transcripts behave according to a driven damped harmonic oscillator. As expected, the sleep-wake distribution accounted for a large degree of the changes in EEG delta power and Homer1a. Nevertheless, the sleep deprivation-induced increase in delta power varied also with time of day with higher than expected levels when recovery sleep started at dark onset. CONCLUSIONS: Per2 and delta power are widely used as exclusive state variables of the circadian and homeostatic process, respectively. Our findings demonstrate a considerable cross-talk between these two processes. As Per2 in the brain responds to both sleep loss and time of day, this molecule is well positioned to keep track of and to anticipate homeostatic sleep need. CITATION: Curie T; Mongrain V; Dorsaz S; Mang GM; Emmenegger Y; Franken P. Homeostatic and circadian contribution to EEG and molecular state variables of sleep regulation. SLEEP 2013;36(3):311-323.

The Coordinated Radio and Infrared Survey for High-Mass Star Formation (The CORNISH Survey). I. Survey Design

We describe the motivation, design, and implementation of the CORNISH survey, an arcsecondresolution radio continuum survey of the inner galactic plane at 5 GHz using the Very Large Array (VLA). It is a blind survey coordinated with the northern SpitzerGLIMPSE I region covering 10° detail the strategy that we employed to control the shape of the synthesised beam across this survey, which covers a wide range of fairly low declinations. Two snapshots separated by 4 h kept the beam elongation to less that 1.5 over 75% of the survey area and less than 2 over 98% of the survey. The prime scientific motivation is to provide an unbiased survey for ultra-compact H II regions to study this key phase in massive star formation. A sensitivity around 2 mJy will allow the automatic distinction between radio-loud and radio-quiet midIR sources found in the Spitzersurveys. This survey has many legacy applications beyond star formation, including evolved stars, active stars and binaries, and extragalactic sources. The CORNISH survey for compact ionized sources complements other Galactic plane surveys that target diffuse and nonthermal sources, as well as atomic and molecular phases to build up a complete picture of the interstellar medium in the Galaxy.

Influence of ionization state on the activation of temocapril by hCES1: a molecular-dynamics study.

Temocapril is a prodrug whose hydrolysis by carboxylesterase 1 (CES1) yields the active ACE inhibitor temocaprilat. This molecular-dynamics (MD) study uses a resolved structure of the human CES1 (hCES1) to investigate some mechanistic details of temocapril hydrolysis. The ionization constants of temocapril (pK1 and pK3) and temocaprilat (pK1, pK2, and pK3) were determined experimentally and computationally using commercial algorithms. The constants so obtained were in good agreement and revealed that temocapril exists mainly in three ionic forms (a cation, a zwitterion, and an anion), whereas temocaprilat exists in four major ionic forms (a cation, a zwitterion, an anion, and a dianion). All these ionic forms were used as ligands in 5-ns MS simulations. While the cationic and zwitterionic forms of temocapril were involved in an ion-pair bond with Glu255 suggestive of an inhibitor behavior, the anionic form remained in a productive interaction with the catalytic center. As for temocaprilat, its cation appeared trapped by Glu255, while its zwitterion and anion made a slow departure from the catalytic site and a partial egress from the protein. Only its dianion was effectively removed from the catalytic site and attracted to the protein surface by Lys residues. A detailed mechanism of product egress emerges from the simulations.

The effects of dominance, regular inbreeding and sampling design on Q(ST), an estimator of population differentiation for quantitative traits.

To test whether quantitative traits are under directional or homogenizing selection, it is common practice to compare population differentiation estimates at molecular markers (F(ST)) and quantitative traits (Q(ST)). If the trait is neutral and its determinism is additive, then theory predicts that Q(ST) = F(ST), while Q(ST) > F(ST) is predicted under directional selection for different local optima, and Q(ST) < F(ST) is predicted under homogenizing selection. However, nonadditive effects can alter these predictions. Here, we investigate the influence of dominance on the relation between Q(ST) and F(ST) for neutral traits. Using analytical results and computer simulations, we show that dominance generally deflates Q(ST) relative to F(ST). Under inbreeding, the effect of dominance vanishes, and we show that for selfing species, a better estimate of Q(ST) is obtained from selfed families than from half-sib families. We also compare several sampling designs and find that it is always best to sample many populations (>20) with few families (five) rather than few populations with many families. Provided that estimates of Q(ST) are derived from individuals originating from many populations, we conclude that the pattern Q(ST) > F(ST), and hence the inference of directional selection for different local optima, is robust to the effect of nonadditive gene actions.

Molecular modeling of peptide-MHC class I complexes : applications to peptide based immunotherapy

Summary The specific CD8+ T cell immune response against tumors relies on the recognition by the T cell receptor (TCR) on cytotoxic T lymphocytes (CTL) of antigenic peptides bound to the class I major histocompatibility complex (MHC) molecule. Such tumor associated antigenic peptides are the focus of tumor immunotherapy with peptide vaccines. The strategy for obtaining an improved immune response often involves the design of modified tumor associated antigenic peptides. Such modifications aim at creating higher affinity and/or degradation resistant peptides and require precise structures of the peptide-MHC class I complex. In addition, the modified peptide must be cross-recognized by CTLs specific for the parental peptide, i.e. preserve the structure of the epitope. Detailed structural information on the modified peptide in complex with MHC is necessary for such predictions. In this thesis, the main focus is the development of theoretical in silico methods for prediction of both structure and cross-reactivity of peptide-MHC class I complexes. Applications of these methods in the context of immunotherapy are also presented. First, a theoretical method for structure prediction of peptide-MHC class I complexes is developed and validated. The approach is based on a molecular dynamics protocol to sample the conformational space of the peptide in its MHC environment. The sampled conformers are evaluated using conformational free energy calculations. The method, which is evaluated for its ability to reproduce 41 X-ray crystallographic structures of different peptide-MHC class I complexes, shows an overall prediction success of 83%. Importantly, in the clinically highly relevant subset of peptide-HLAA*0201 complexes, the prediction success is 100%. Based on these structure predictions, a theoretical approach for prediction of cross-reactivity is developed and validated. This method involves the generation of quantitative structure-activity relationships using three-dimensional molecular descriptors and a genetic neural network. The generated relationships are highly predictive as proved by high cross-validated correlation coefficients (0.78-0.79). Together, the here developed theoretical methods open the door for efficient rational design of improved peptides to be used in immunotherapy. Résumé La réponse immunitaire spécifique contre des tumeurs dépend de la reconnaissance par les récepteurs des cellules T CD8+ de peptides antigéniques présentés par les complexes majeurs d'histocompatibilité (CMH) de classe I. Ces peptides sont utilisés comme cible dans l'immunothérapie par vaccins peptidiques. Afin d'augmenter la réponse immunitaire, les peptides sont modifiés de façon à améliorer l'affinité et/ou la résistance à la dégradation. Ceci nécessite de connaître la structure tridimensionnelle des complexes peptide-CMH. De plus, les peptides modifiés doivent être reconnus par des cellules T spécifiques du peptide natif. La structure de l'épitope doit donc être préservée et des structures détaillées des complexes peptide-CMH sont nécessaires. Dans cette thèse, le thème central est le développement des méthodes computationnelles de prédiction des structures des complexes peptide-CMH classe I et de la reconnaissance croisée. Des applications de ces méthodes de prédiction à l'immunothérapie sont également présentées. Premièrement, une méthode théorique de prédiction des structures des complexes peptide-CMH classe I est développée et validée. Cette méthode est basée sur un échantillonnage de l'espace conformationnel du peptide dans le contexte du récepteur CMH classe I par dynamique moléculaire. Les conformations sont évaluées par leurs énergies libres conformationnelles. La méthode est validée par sa capacité à reproduire 41 structures des complexes peptide-CMH classe I obtenues par cristallographie aux rayons X. Le succès prédictif général est de 83%. Pour le sous-groupe HLA-A*0201 de complexes de grande importance pour l'immunothérapie, ce succès est de 100%. Deuxièmement, à partir de ces structures prédites in silico, une méthode théorique de prédiction de la reconnaissance croisée est développée et validée. Celle-ci consiste à générer des relations structure-activité quantitatives en utilisant des descripteurs moléculaires tridimensionnels et un réseau de neurones couplé à un algorithme génétique. Les relations générées montrent une capacité de prédiction remarquable avec des valeurs de coefficients de corrélation de validation croisée élevées (0.78-0.79). Les méthodes théoriques développées dans le cadre de cette thèse ouvrent la voie du design de vaccins peptidiques améliorés.