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.

Immunogenetics of invasive aspergillosis.

Invasive aspergillosis is one of the most important infections in hematopoietic stem cell transplant recipients, with an incidence rate of 5-15% and an associated mortality of 30-60%. It remains unclear why certain patients develop invasive aspergillosis while others, undergoing identical transplant regimen and similar post transplant immunosuppression, do not. Over the last decade, pattern recognition receptors such as Toll-like receptors (TLRs) and the C-type lectin receptors (CLRs) have emerged as critical components of the innate immune system. By detecting specific molecular patterns from invading microbes and initiating inflammatory and subsequent adaptive immune responses, pattern recognition receptors are strategically located at the molecular interface of hosts and pathogens. Polymorphisms in pattern recognition receptors and downstream signaling molecules have been associated with increased or decreased susceptibility to infections, suggesting that their detection may have an increasing impact on the treatment and prevention of infectious diseases in the coming years. Infectious risk stratification may be particularly relevant for patients with hematologic malignancies, because of the high prevalence and severity of infections in this population. This review summarizes the innate immune mechanisms involved in Aspergillus fumigatus detection and the role of host genetic polymorphisms in susceptibility to invasive aspergillosis.

SwissParam: a fast force field generation tool for small organic molecules.

The drug discovery process has been deeply transformed recently by the use of computational ligand-based or structure-based methods, helping the lead compounds identification and optimization, and finally the delivery of new drug candidates more quickly and at lower cost. Structure-based computational methods for drug discovery mainly involve ligand-protein docking and rapid binding free energy estimation, both of which require force field parameterization for many drug candidates. Here, we present a fast force field generation tool, called SwissParam, able to generate, for arbitrary small organic molecule, topologies, and parameters based on the Merck molecular force field, but in a functional form that is compatible with the CHARMM force field. Output files can be used with CHARMM or GROMACS. The topologies and parameters generated by SwissParam are used by the docking software EADock2 and EADock DSS to describe the small molecules to be docked, whereas the protein is described by the CHARMM force field, and allow them to reach success rates ranging from 56 to 78%. We have also developed a rapid binding free energy estimation approach, using SwissParam for ligands and CHARMM22/27 for proteins, which requires only a short minimization to reproduce the experimental binding free energy of 214 ligand-protein complexes involving 62 different proteins, with a standard error of 2.0 kcal mol(-1), and a correlation coefficient of 0.74. Together, these results demonstrate the relevance of using SwissParam topologies and parameters to describe small organic molecules in computer-aided drug design applications, together with a CHARMM22/27 description of the target protein. SwissParam is available free of charge for academic users at www.swissparam.ch.

Synthesis of a radioiodinated photoreactive MAGE-1 peptide derivative and photoaffinity labeling of cell-associated human leukocyte antigen-A1 molecules.

The synthesis of a photoreactive derivative of the human leukocyte antigen-A1 (HLA-A1)-restricted MAGE-1 peptide 161-169 (EADPTGHSY) is described. Using conventional automated solid-phase peptide synthesis, a photoreactive derivative of this peptide was synthesized by replacing histidine-167 with photo-reactive N-beta-4-azidosalicyloyl-L-2,3-diaminopropionic acid. The C-terminal tyrosine was incorporated as phosphotyrosine. This peptide derivative was radioiodinated in the presence of chloramine T. This iodination took place selectively at the photoreactive group, because the phosphate ester prevented tyrosine iodination. Following dephosphorylation with alkaline phosphatase and chromatographic purification, the radiolabeled peptide derivative was incubated with cells expressing HLA-A1 or other HLA molecules. Photoactivation resulted in efficient photoaffinity labeling of HLA-A1. Other HLA molecules or other cellular components were not detectably labeled. This labeling was inhibited by HLA-A1 but not by HLA-A2-binding peptides. This synthesis is generally applicable and can also be adapted to the synthesis of well-defined radiolabeled nonphotoreactive peptide derivatives.

Molecular characterization of HLA-C incompatibilities in HLA-ABDR-matched unrelated bone marrow donor-recipient pairs. Sequence of two new Cw alleles (Cw*02023 and Cw*0707) and recognition by cytotoxic T lymphocytes.

While the influence of HLA-AB and -DRB1 matching on the outcome of bone marrow transplantation (BMT) with unrelated donors is clear, the evaluation of HLA-C has been hampered by its poor serological definition. Because the low resolution of standard HLA-C typing could explain the significant number of positive cytotoxic T lymphocyte precursor frequency (CTLpf) tests found among HLA-AB-subtype, DRB1/B3/B5-subtype matched patient/donor pairs, we have identified by sequencing the incompatibilities recognized by CD8+ CTL clones obtained from such positive CTLpf tests. In most cases the target molecules were HLA-C antigens that had escaped detection by serology (e.g. Cw*1601, 1502 or 0702). Direct recognition of HLA-C by a CTL clone was demonstrated by lysis of the HLA class I-negative 721.221 cell line transfected with Cw*1601 cDNA. Because of the functional importance of Cw polymorphism, a PCR-SSO oligotyping procedure was set up allowing the resolution of 29 Cw alleles. Oligotyping of a panel of 382 individuals (including 101 patients and their 272 potential unrelated donors, 5 related donors and 4 platelet donors) allowed to determine HLA-C and HLA A-B-Cw-DRB1 allelic frequencies, as well as a number of A-Cw, B-Cw, and DRB1-Cw associations. Two new HLA-Cw alleles (Cw*02023 and Cw*0707) were identified by DNA sequencing of PCR-amplified exon 2-intron 2-exon 3 amplicons. Furthermore, we determined the degree of HLA-C compatibility in 287 matched pairs that could be formed from 73 patients and their 184 potential unrelated donors compatible for HLA-AB by serology and for HLA-DRB1/ B3/B5 by oligotyping. Cw mismatches were identified in 42.1% of these pairs, and AB-subtype oligotyping showed that 30% of these Cw-incompatible pairs were also mismatched for A or B-locus subtype. The degree of HLA-C incompatibility was strongly influenced by the linkage with B alleles and by the ABDR haplotypes. Cw alleles linked with B*4403, B*5101, B18, and B62 haplotypes were frequently mismatched. Apparently high resolution DNA typing for HLA-AB does not result in full matching at locus C. Since HLA-C polymorphism is recognized by alloreactive CTLs, such incompatibilities might be as relevant as AB-subtype mismatches in clinical transplantation.

Direct visualization of supercoiled DNA molecules in solution.

The shape of supercoiled DNA molecules in solution is directly visualized by cryo-electron microscopy of vitrified samples. We observe that: (i) supercoiled DNA molecules in solution adopt an interwound rather than a toroidal form, (ii) the diameter of the interwound superhelix changes from about 12 nm to 4 nm upon addition of magnesium salt to the solution and (iii) the partition of the linking deficit between twist and writhe can be quantitatively determined for individual molecules.

Reduced development of CD4-8-B220+ T cells but normal autoantibody production in lpr/lpr mice lacking major histocompatibility complex class I molecules.

The lpr gene has recently been shown to encode a functional mutation in the Fas receptor, a molecule involved in transducing apoptotic signals. Mice homozygous for the lpr gene develop an autoimmune syndrome accompanied by massive accumulation of double-negative (DN) CD4-8-B220+ T cell receptor-alpha/beta+ cells. In order to investigate the origin of these DN T cells, we derived lpr/lpr mice lacking major histocompatibility complex (MHC) class I molecules by intercrossing them with beta 2-microglobulin (beta 2m)-deficient mice. Interestingly, these lpr beta 2m-/- mice develop 13-fold fewer DNT cells in lymph nodes as compared to lpr/lpr wild-type (lprWT) mice. Analysis of anti-DNA antibodies and rheumatoid factor in serum demonstrates that lpr beta 2m-/- mice produce comparable levels of autoantibodies to lprWT mice. Collectively our data indicate that MHC class I molecules control the development of DN T cells but not autoantibody production in lpr/lpr mice and support the hypothesis that the majority of DN T cells may be derived from cells of the CD8 lineage.

Dissociation of thymic positive and negative selection in transgenic mice expressing major histocompatibility complex class I molecules exclusively on thymic cortical epithelial cells.

Thymic positive and negative selection of developing T lymphocytes confronts us with a paradox: How can a T-cell antigen receptor (TCR)-major histocompatibility complex (MHC)/peptide interaction in the former process lead to transduction of signals allowing for cell survival and in the latter induce programmed cell death or a hyporesponsive state known as anergy? One of the hypotheses put forward states that the outcome of a TCR-MHC/peptide interaction depends on the cell type presenting the selecting ligand to the developing thymocyte. Here we describe the development and lack of self-tolerance of CD8(+) T lymphocytes in transgenic mice expressing MHC class I molecules in the thymus exclusively on cortical epithelial cells. Despite the absence of MHC class I expression on professional antigen-presenting cells, normal numbers of CD8(+) cells were observed in the periphery. Upon specific activation, transgenic CD8(+) T cells efficiently lysed syngeneic MHC class I(+) targets in vitro and in vivo, indicating that thymic cortical epithelium (in contrast to medullary epithelium and antigen-presenting cells of hematopoietic origin) is incapable of tolerance induction. Thus, compartmentalization of the antigen-presenting cells involved in thymic positive selection and tolerance induction can (at least in part) explain the positive/negative selection paradox.

Phenotypic switching in Pseudomonas brassicacearum involves GacS- and GacA-dependent Rsm small RNAs.

The plant-beneficial bacterium Pseudomonas brassicacearum forms phenotypic variants in vitro as well as in planta during root colonization under natural conditions. Transcriptome analysis of typical phenotypic variants using microarrays containing coding as well as noncoding DNA fragments showed differential expression of several genes relevant to secondary metabolism and of the small RNA (sRNA) genes rsmX, rsmY, and rsmZ. Naturally occurring mutations in the gacS-gacA system accounted for phenotypic switching, which was characterized by downregulation of antifungal secondary metabolites (2,4-diacetylphloroglucinol and cyanide), indoleacetate, exoenzymes (lipase and protease), and three different N-acyl-homoserine lactone molecules. Moreover, in addition to abrogating these biocontrol traits, gacS and gacA mutations resulted in reduced expression of the type VI secretion machinery, alginate biosynthesis, and biofilm formation. In a gacA mutant, the expression of rsmX was completely abolished, unlike that of rsmY and rsmZ. Overexpression of any of the three sRNAs in the gacA mutant overruled the pleiotropic changes and restored the wild-type phenotypes, suggesting functional redundancy of these sRNAs. In conclusion, our data show that phenotypic switching in P. brassicacearum results from mutations in the gacS-gacA system.

Monitoring inflammatory bowel disease activity: clinical activity is judged to be more relevant than endoscopic severity or biomarkers.

BACKGROUND: There is increasing evidence for the clinical relevance of mucosal healing (MH) as therapeutic treatment goal in inflammatory bowel disease (IBD). We aimed to investigate by which method gastroenterologists monitor IBD activity in daily practice. METHODS: A questionnaire was sent to all board-certified gastroenterologists in Switzerland to specifically address their strategy to monitor IBD between May 2009 and April 2010. RESULTS: The response rate was 57% (153/270). Fifty-two percent of gastroenterologists worked in private practice and 48% worked in hospitals. Seventy-eight percent judged clinical activity to be the most relevant criterion for monitoring IBD activity, 15% chose endoscopic severity, and 7% chose biomarkers. Seventy percent of gastroenterologists based their therapeutic decisions on clinical activity, 24% on endoscopic severity, and 6% on biomarkers. The following biomarkers were used for IBD activity monitoring: CRP, 94%; differential blood count, 78%; fecal calprotectin (FC), 74%; iron status, 63%; blood sedimentation rate, 3%; protein electrophoresis, 0.7%; fecal neutrophils, 0.7%; and vitamin B12, 0.7%. Gastroenterologists in hospitals and those with ≤ 10 years of professional experience used FC more frequently compared with colleagues in private practice (P=0.035) and those with > 10 years of experience (P<0.001). CONCLUSIONS: Clinical activity is judged to be more relevant for monitoring IBD activity and guiding therapeutic decisions than endoscopic severity and biomarkers. As such, the accumulating scientific evidence on the clinical impact of mucosal healing does not yet seem to influence the management of IBD in daily gastroenterologic practice.

Shaping the interaction landscape of bioactive molecules.

MOTIVATION: Most bioactive molecules perform their action by interacting with proteins or other macromolecules. However, for a significant fraction of them, the primary target remains unknown. In addition, the majority of bioactive molecules have more than one target, many of which are poorly characterized. Computational predictions of bioactive molecule targets based on similarity with known ligands are powerful to narrow down the number of potential targets and to rationalize side effects of known molecules. RESULTS: Using a reference set of 224 412 molecules active on 1700 human proteins, we show that accurate target prediction can be achieved by combining different measures of chemical similarity based on both chemical structure and molecular shape. Our results indicate that the combined approach is especially efficient when no ligand with the same scaffold or from the same chemical series has yet been discovered. We also observe that different combinations of similarity measures are optimal for different molecular properties, such as the number of heavy atoms. This further highlights the importance of considering different classes of similarity measures between new molecules and known ligands to accurately predict their targets. CONTACT: olivier.michielin@unil.ch or vincent.zoete@unil.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Médecine interne générale ambulatoire: quelques enseignements de la littérature en 2005 [Primary care in general internal medicine: relevant information in 2005]

During the year 2005, the chief residents of the University Medical Outpatient Clinic of Lausanne have done a database of useful articles for daily practice, scientifically validated and with excellent didactic quality, from 10 electronic journals. They have used those selected articles in personal meetings between the chief residents on a regular basis and the possibility to access the database by the junior physicians. Six of these articles concerning different topics (depression, tuberculosis detection, anticoagulation at home, cholinesterase inhibitors, insomnia and therapies, transdermal nitroglycerin and tendinopathies) are presented.

Whose DNA is this? How relevant a question? (a note for forensic scientists).

This communication seeks to draw the attention of researchers and practitioners dealing with forensic DNA profiling analyses to the following question: is a scientist's report, offering support to a hypothesis according to which a particular individual is the source of DNA detected during the analysis of a stain, relevant from the point of view of a Court of Justice? This question relates to skeptical views previously voiced by commentators mainly in the judicial area, but is avoided by a large majority of forensic scientists. Notwithstanding, the pivotal role of this question has recently been evoked during the international conference "The hidden side of DNA profiles. Artifacts, errors and uncertain evidence" held in Rome (April 27th to 28th, 2012). Indeed, despite the fact that this conference brought together some of the world's leading forensic DNA specialists, it appeared clearly that a huge gap still exists between questions lawyers are actually interested in, and the answers that scientists deliver to Courts in written reports or during oral testimony. Participants in the justice system, namely lawyers and jurors on the one hand and forensic geneticists on the other, unfortunately talk considerably different languages. It thus is fundamental to address this issue of communication about results of forensic DNA analyses, and open a dialogue with practicing non-scientists at large who need to make meaningful use of scientific results to approach and help solve judicial cases. This paper intends to emphasize the actuality of this topic and suggest beneficial ways ahead towards a more reasoned use of forensic DNA in criminal proceedings.