Application de la modélisation moléculaire à de nouvelles approches thérapeutiques du cancer [Application of molecular modeling to new therapeutic cancer approaches].

Recent progress in the experimental determination of protein structures allow to understand, at a very detailed level, the molecular recognition mechanisms that are at the basis of the living matter. This level of understanding makes it possible to design rational therapeutic approaches, in which effectors molecules are adapted or created de novo to perform a given function. An example of such an approach is drug design, were small inhibitory molecules are designed using in silico simulations and tested in vitro. In this article, we present a similar approach to rationally optimize the sequence of killer T lymphocytes receptors to make them more efficient against melanoma cells. The architecture of this translational research project is presented together with its implications both at the level of basic research as well as in the clinics.

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.

Expanding molecular modeling and design tools to non-natural sidechains.

Protein-protein interactions encode the wiring diagram of cellular signaling pathways and their deregulations underlie a variety of diseases, such as cancer. Inhibiting protein-protein interactions with peptide derivatives is a promising way to develop new biological and therapeutic tools. Here, we develop a general framework to computationally handle hundreds of non-natural amino acid sidechains and predict the effect of inserting them into peptides or proteins. We first generate all structural files (pdb and mol2), as well as parameters and topologies for standard molecular mechanics software (CHARMM and Gromacs). Accurate predictions of rotamer probabilities are provided using a novel combined knowledge and physics based strategy. Non-natural sidechains are useful to increase peptide ligand binding affinity. Our results obtained on non-natural mutants of a BCL9 peptide targeting beta-catenin show very good correlation between predicted and experimental binding free-energies, indicating that such predictions can be used to design new inhibitors. Data generated in this work, as well as PyMOL and UCSF Chimera plug-ins for user-friendly visualization of non-natural sidechains, are all available at http://www.swisssidechain.ch. Our results enable researchers to rapidly and efficiently work with hundreds of non-natural sidechains.

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.

SwissSidechain: a molecular and structural database of non-natural sidechains.

Amino acids form the building blocks of all proteins. Naturally occurring amino acids are restricted to a few tens of sidechains, even when considering post-translational modifications and rare amino acids such as selenocysteine and pyrrolysine. However, the potential chemical diversity of amino acid sidechains is nearly infinite. Exploiting this diversity by using non-natural sidechains to expand the building blocks of proteins and peptides has recently found widespread applications in biochemistry, protein engineering and drug design. Despite these applications, there is currently no unified online bioinformatics resource for non-natural sidechains. With the SwissSidechain database (http://www.swisssidechain.ch), we offer a central and curated platform about non-natural sidechains for researchers in biochemistry, medicinal chemistry, protein engineering and molecular modeling. SwissSidechain provides biophysical, structural and molecular data for hundreds of commercially available non-natural amino acid sidechains, both in l- and d-configurations. The database can be easily browsed by sidechain names, families or physico-chemical properties. We also provide plugins to seamlessly insert non-natural sidechains into peptides and proteins using molecular visualization software, as well as topologies and parameters compatible with molecular mechanics software.

The alpha 1a and alpha 1b-adrenergic receptor subtypes: molecular mechanisms of receptor activation and of drug action.

In this chapter we summarize some aspects of the structure-functional relationship of the alpha 1a and alpha 1b-adrenergic receptor subtypes related to the receptor activation process as well as the effect of different alpha-blockers on the constitutive activity of the receptor. Molecular modeling of the alpha 1a and alpha 1b-adrenergic receptor subtypes and computational simulation of receptor dynamics were useful to interpret the experimental findings derived from site directed mutagenesis studies.

Molecular basis of leukocyte rolling on PSGL-1. Predominant role of core-2 O-glycans and of tyrosine sulfate residue 51.

Interactions between the leukocyte adhesion receptor L-selectin and P-selectin glycoprotein ligand-1 play an important role in regulating the inflammatory response by mediating leukocyte tethering and rolling on adherent leukocytes. In this study, we have examined the effect of post-translational modifications of PSGL-1 including Tyr sulfation and presentation of sialylated and fucosylated O-glycans for L-selectin binding. The functional importance of these modifications was determined by analyzing soluble L-selectin binding and leukocyte rolling on CHO cells expressing various glycoforms of PSGL-1 or mutant PSGL-1 targeted at N-terminal Thr or Tyr residues. Simultaneous expression of core-2 beta1,6-N-acetylglucosaminyltransferase and fucosyltransferase VII was required for optimal L-selectin binding to PSGL-1. Substitution of Thr-57 by Ala but not of Thr-44, strongly decreased L-selectin binding and leukocyte rolling on PSGL-1. Substitution of Tyr by Phe revealed that PSGL-1 Tyr-51 plays a predominant role in mediating L-selectin binding and leukocyte rolling whereas Tyr-48 has a minor role, an observation that contrasts with the pattern seen for the interactions between PSGL-1 and P-selectin where Tyr-48 plays a key role. Molecular modeling analysis of L-selectin and P-selectin interactions with PSGL-1 further supported these observations. Additional experiments showed that core-2 O-glycans attached to Thr-57 were also of critical importance in regulating the velocity and stability of leukocyte rolling. These observations pinpoint the structural characteristics of PSGL-1 that are required for optimal interactions with L-selectin and may be responsible for the specific kinetic and mechanical bond properties of the L-selectin-PSGL-1 adhesion receptor-counterreceptor pair.

Fourier transform convolution integrals applied to generalized Born molecular volume.

Generalized Born methods are currently among the solvation models most commonly used for biological applications. We reformulate the generalized Born molecular volume method initially described by (Lee et al, 2003, J Phys Chem, 116, 10606; Lee et al, 2003, J Comp Chem, 24, 1348) using fast Fourier transform convolution integrals. Changes in the initial method are discussed and analyzed. Finally, the method is extensively checked with snapshots from common molecular modeling applications: binding free energy computations and docking. Biologically relevant test systems are chosen, including 855-36091 atoms. It is clearly demonstrated that, precision-wise, the proposed method performs as good as the original, and could better benefit from hardware accelerated boards.

Autosomal-Recessive Posterior Microphthalmos Is Caused by Mutations in PRSS56, a Gene Encoding a Trypsin-Like Serine Protease.

Posterior microphthalmos (MCOP) is a rare isolated developmental anomaly of the eye characterized by extreme hyperopia due to short axial length. The population of the Faroe Islands shows a high prevalence of an autosomal-recessive form (arMCOP) of the disease. Based on published linkage data, we refined the position of the disease locus (MCOP6) in an interval of 250 kb in chromosome 2q37.1 in two large Faroese families. We detected three different mutations in PRSS56. Patients of the Faroese families were either homozygous for c.926G>C (p.Trp309Ser) or compound heterozygous for c.926G>C and c.526C>G (p.Arg176Gly), whereas a homozygous 1 bp duplication (c.1066dupC) was identified in five patients with arMCOP from a consanguineous Tunisian family. In one patient with MCOP from the Faroe Islands and in another one from Turkey, no PRSS56 mutation was detected, suggesting nonallelic heterogeneity of the trait. Using RT-PCR, PRSS56 transcripts were detected in samples derived from the human adult retina, cornea, sclera, and optic nerve. The expression of the mouse ortholog could be first detected in the eye at E17 and was maintained into adulthood. The predicted PRSS56 protein is a 603 amino acid long secreted trypsin-like serine peptidase. The c.1066dupC is likely to result in a functional null allele, whereas the two point mutations predict the replacement of evolutionary conserved and functionally important residues. Molecular modeling of the p.Trp309Ser mutant suggests that both the affinity and reactivity of the enzyme toward in vivo protein substrates are likely to be substantially reduced.

Enzymic, phylogenetic, and structural characterization of the unusual papain-like protease domain of Plasmodium falciparum SERA5.

Serine repeat antigen 5 (SERA5) is an abundant antigen of the human malaria parasite Plasmodium falciparum and is the most strongly expressed member of the nine-gene SERA family. It appears to be essential for the maintenance of the erythrocytic cycle, unlike a number of other members of this family, and has been implicated in parasite egress and/or erythrocyte invasion. All SERA proteins possess a central domain that has homology to papain except in the case of SERA5 (and some other SERAs), where the active site cysteine has been replaced with a serine. To investigate if this domain retains catalytic activity, we expressed, purified, and refolded a recombinant form of the SERA5 enzyme domain. This protein possessed chymotrypsin-like proteolytic activity as it processed substrates downstream of aromatic residues, and its activity was reversed by the serine protease inhibitor 3,4-diisocoumarin. Although all Plasmodium SERA enzyme domain sequences share considerable homology, phylogenetic studies revealed two distinct clusters across the genus, separated according to whether they possess an active site serine or cysteine. All Plasmodia appear to have at least one member of each group. Consistent with separate biological roles for members of these two clusters, molecular modeling studies revealed that SERA5 and SERA6 enzyme domains have dramatically different surface properties, although both have a characteristic papain-like fold, catalytic cleft, and an appropriately positioned catalytic triad. This study provides impetus for the examination of SERA5 as a target for antimalarial drug design.

SwissDock, a protein-small molecule docking web service based on EADock DSS.

Most life science processes involve, at the atomic scale, recognition between two molecules. The prediction of such interactions at the molecular level, by so-called docking software, is a non-trivial task. Docking programs have a wide range of applications ranging from protein engineering to drug design. This article presents SwissDock, a web server dedicated to the docking of small molecules on target proteins. It is based on the EADock DSS engine, combined with setup scripts for curating common problems and for preparing both the target protein and the ligand input files. An efficient Ajax/HTML interface was designed and implemented so that scientists can easily submit dockings and retrieve the predicted complexes. For automated docking tasks, a programmatic SOAP interface has been set up and template programs can be downloaded in Perl, Python and PHP. The web site also provides an access to a database of manually curated complexes, based on the Ligand Protein Database. A wiki and a forum are available to the community to promote interactions between users. The SwissDock web site is available online at http://www.swissdock.ch. We believe it constitutes a step toward generalizing the use of docking tools beyond the traditional molecular modeling community.

Cathepsin X cleavage of the beta2 integrin regulates talin-binding and LFA-1 affinity in T cells.

T cell migration, essential for immune surveillance and response, is mediated by the integrin LFA-1. CatX, a cysteine carboxypeptidase, is involved in the regulation of T cell migration by interaction with LFA-1. We show that sequential cleavage of C-terminal amino acids from the β(2) cytoplasmic tail of LFA-1, by CatX, enhances binding of the adaptor protein talin to LFA-1 and triggers formation of the latter's high-affinity form. As shown by SPR analysis of peptides constituting the truncated β(2) tail, the cleavage of three C-terminal amino acids by CatX resulted in a 1.6-fold increase of talin binding. Removal of one more amino acid resulted in a 2.5-fold increase over the intact tail. CatX cleavage increased talin-binding affinity to the MD but not the MP talin-binding site on the β(2) tail. This was shown by molecular modeling of the β(2) tail/talin F3 complex to be a result of conformational changes affecting primarily the distal-binding site. Analysis of LFA-1 by conformation-specific mAb showed that CatX modulates LFA-1 affinity, promoting formation of high-affinity from intermediate-affinity LFA-1 but not the initial activation of LFA-1 from a bent to extended form. CatX post-translational modifications may thus represent a mechanism of LFA-1 fine-tuning that enables the trafficking of T cells.

Structural analysis of TCR-ligand interactions studied on H-2Kd-restricted cloned CTL specific for a photoreactive peptide derivative.

To study the interaction of the TCR with its ligand, the complex of a MHC molecule and an antigenic peptide, we modified a TCR contact residue of a H-2Kd-restricted antigenic peptide with photoreactive 4-azidobenzoic acid. The photoreactive group was a critical component of the epitope recognized by CTL clones derived from mice immunized with such a peptide derivative. The majority of these clones expressed V beta 1-encoded beta chains that were paired with J alpha TA28-encoded alpha chains. For one of these TCR, the photoaffinity labeled sites were mapped on the alpha chain as a J alpha TA28-encoded tryptophan and on the beta chain as a residue of the C' strand of V beta 1. Molecular modeling of this TCR suggested the presence of a hydrophobic pocket that harbors this tryptophan as well as a tyrosine on the C' strand of V beta 1 between which the photoreactive side chain inserts. It is concluded that this avid binding principle may account for the preferential selection of V beta 1 and J alpha TA28-encoded TCR.

Monoubiquitination and Activity of the Paracaspase MALT1 Requires Glutamate 549 in the Dimerization Interface.

The mucosa-associated lymphoid tissue protein-1 (MALT1, also known as paracaspase) is a protease whose activity is essential for the activation of lymphocytes and the growth of cells derived from human diffuse large B-cell lymphomas of the activated B-cell subtype (ABC DLBCL). Crystallographic approaches have shown that MALT1 can form dimers via its protease domain, but why dimerization is relevant for the biological activity of MALT1 remains largely unknown. Using a molecular modeling approach, we predicted Glu 549 (E549) to be localized within the MALT1 dimer interface and thus potentially relevant. Experimental mutation of this residue into alanine (E549A) led to a complete impairment of MALT1 proteolytic activity. This correlated with an impaired capacity of the mutant to form dimers of the protease domain in vitro, and a reduced capacity to promote NF-κB activation and transcription of the growth-promoting cytokine interleukin-2 in antigen receptor-stimulated lymphocytes. Moreover, this mutant could not rescue the growth of ABC DLBCL cell lines upon MALT1 silencing. Interestingly, the MALT1 mutant E549A was unable to undergo monoubiquitination, which we identified previously as a critical step in MALT1 activation. Collectively, these findings suggest a model in which E549 at the dimerization interface is required for the formation of the enzymatically active, monoubiquitinated form of MALT1.