Ab initio modeling and molecular dynamics simulation of the alpha 1b-adrenergic receptor activation.

This work describes the ab initio procedure employed to build an activation model for the alpha 1b-adrenergic receptor (alpha 1b-AR). The first version of the model was progressively modified and complicated by means of a many-step iterative procedure characterized by the employment of experimental validations of the model in each upgrading step. A combined simulated (molecular dynamics) and experimental mutagenesis approach was used to determine the structural and dynamic features characterizing the inactive and active states of alpha 1b-AR. The latest version of the model has been successfully challenged with respect to its ability to interpret and predict the functional properties of a large number of mutants. The iterative approach employed to describe alpha 1b-AR activation in terms of molecular structure and dynamics allows further complications of the model to allow prediction and interpretation of an ever-increasing number of experimental data.

Structural prediction of peptides bound to MHC class I.

An ab initio structure prediction approach adapted to the peptide-major histocompatibility complex (MHC) class I system is presented. Based on structure comparisons of a large set of peptide-MHC class I complexes, a molecular dynamics protocol is proposed using simulated annealing (SA) cycles to sample the conformational space of the peptide in its fixed MHC environment. A set of 14 peptide-human leukocyte antigen (HLA) A0201 and 27 peptide-non-HLA A0201 complexes for which X-ray structures are available is used to test the accuracy of the prediction method. For each complex, 1000 peptide conformers are obtained from the SA sampling. A graph theory clustering algorithm based on heavy atom root-mean-square deviation (RMSD) values is applied to the sampled conformers. The clusters are ranked using cluster size, mean effective or conformational free energies, with solvation free energies computed using Generalized Born MV 2 (GB-MV2) and Poisson-Boltzmann (PB) continuum models. The final conformation is chosen as the center of the best-ranked cluster. With conformational free energies, the overall prediction success is 83% using a 1.00 Angstroms crystal RMSD criterion for main-chain atoms, and 76% using a 1.50 Angstroms RMSD criterion for heavy atoms. The prediction success is even higher for the set of 14 peptide-HLA A0201 complexes: 100% of the peptides have main-chain RMSD values < or =1.00 Angstroms and 93% of the peptides have heavy atom RMSD values < or =1.50 Angstroms. This structure prediction method can be applied to complexes of natural or modified antigenic peptides in their MHC environment with the aim to perform rational structure-based optimizations of tumor vaccines.

Access of extracellular cations to their binding sites in Na,K-ATPase: role of the second extracellular loop of the alpha subunit.

Na,K-ATPase, the main active transport system for monovalent cations in animal cells, is responsible for maintaining Na(+) and K(+) gradients across the plasma membrane. During its transport cycle it binds three cytoplasmic Na(+) ions and releases them on the extracellular side of the membrane, and then binds two extracellular K(+) ions and releases them into the cytoplasm. The fourth, fifth, and sixth transmembrane helices of the alpha subunit of Na,K-ATPase are known to be involved in Na(+) and K(+) binding sites, but the gating mechanisms that control the access of these ions to their binding sites are not yet fully understood. We have focused on the second extracellular loop linking transmembrane segments 3 and 4 and attempted to determine its role in gating. We replaced 13 residues of this loop in the rat alpha1 subunit, from E314 to G326, by cysteine, and then studied the function of these mutants using electrophysiological techniques. We analyzed the results using a structural model obtained by homology with SERCA, and ab initio calculations for the second extracellular loop. Four mutants were markedly modified by the sulfhydryl reagent MTSET, and we investigated them in detail. The substituted cysteines were more readily accessible to MTSET in the E1 conformation for the Y315C, W317C, and I322C mutants. Mutations or derivatization of the substituted cysteines in the second extracellular loop resulted in major increases in the apparent affinity for extracellular K(+), and this was associated with a reduction in the maximum activity. The changes produced by the E314C mutation were reversed by MTSET treatment. In the W317C and I322C mutants, MTSET also induced a moderate shift of the E1/E2 equilibrium towards the E1(Na) conformation under Na/Na exchange conditions. These findings indicate that the second extracellular loop must be functionally linked to the gating mechanism that controls the access of K(+) to its binding site.

Modeling of the TCR-MHC-peptide complex.

The methodology for generating a homology model of the T1 TCR-PbCS-K(d) class I major histocompatibility complex (MHC) class I complex is presented. The resulting model provides a qualitative explanation of the effect of over 50 different mutations in the region of the complementarity determining region (CDR) loops of the T cell receptor (TCR), the peptide and the MHC's alpha(1)/alpha(2) helices. The peptide is modified by an azido benzoic acid photoreactive group, which is part of the epitope recognized by the TCR. The construction of the model makes use of closely related homologs (the A6 TCR-Tax-HLA A2 complex, the 2C TCR, the 14.3.d TCR Vbeta chain, the 1934.4 TCR Valpha chain, and the H-2 K(b)-ovalbumine peptide), ab initio sampling of CDR loops conformations and experimental data to select from the set of possibilities. The model shows a complex arrangement of the CDR3alpha, CDR1beta, CDR2beta and CDR3beta loops that leads to the highly specific recognition of the photoreactive group. The protocol can be applied systematically to a series of related sequences, permitting the analysis at the structural level of the large TCR repertoire specific for a given peptide-MHC complex.

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.

Mutational and computational analysis of the alpha(1b)-adrenergic receptor. Involvement of basic and hydrophobic residues in receptor activation and G protein coupling.

To investigate their role in receptor coupling to G(q), we mutated all basic amino acids and some conserved hydrophobic residues of the cytosolic surface of the alpha(1b)-adrenergic receptor (AR). The wild type and mutated receptors were expressed in COS-7 cells and characterized for their ligand binding properties and ability to increase inositol phosphate accumulation. The experimental results have been interpreted in the context of both an ab initio model of the alpha(1b)-AR and of a new homology model built on the recently solved crystal structure of rhodopsin. Among the twenty-three basic amino acids mutated only mutations of three, Arg(254) and Lys(258) in the third intracellular loop and Lys(291) at the cytosolic extension of helix 6, markedly impaired the receptor-mediated inositol phosphate production. Additionally, mutations of two conserved hydrophobic residues, Val(147) and Leu(151) in the second intracellular loop had significant effects on receptor function. The functional analysis of the receptor mutants in conjunction with the predictions of molecular modeling supports the hypothesis that Arg(254), Lys(258), as well as Leu(151) are directly involved in receptor-G protein interaction and/or receptor-mediated activation of the G protein. In contrast, the residues belonging to the cytosolic extensions of helices 3 and 6 play a predominant role in the activation process of the alpha(1b)-AR. These findings contribute to the delineation of the molecular determinants of the alpha(1b)-AR/G(q) interface.

TCRep 3D: an automated in silico approach to study the structural properties of TCR repertoires.

TCRep 3D is an automated systematic approach for TCR-peptide-MHC class I structure prediction, based on homology and ab initio modeling. It has been considerably generalized from former studies to be applicable to large repertoires of TCR. First, the location of the complementary determining regions of the target sequences are automatically identified by a sequence alignment strategy against a database of TCR Vα and Vβ chains. A structure-based alignment ensures automated identification of CDR3 loops. The CDR are then modeled in the environment of the complex, in an ab initio approach based on a simulated annealing protocol. During this step, dihedral restraints are applied to drive the CDR1 and CDR2 loops towards their canonical conformations, described by Al-Lazikani et. al. We developed a new automated algorithm that determines additional restraints to iteratively converge towards TCR conformations making frequent hydrogen bonds with the pMHC. We demonstrated that our approach outperforms popular scoring methods (Anolea, Dope and Modeller) in predicting relevant CDR conformations. Finally, this modeling approach has been successfully applied to experimentally determined sequences of TCR that recognize the NY-ESO-1 cancer testis antigen. This analysis revealed a mechanism of selection of TCR through the presence of a single conserved amino acid in all CDR3β sequences. The important structural modifications predicted in silico and the associated dramatic loss of experimental binding affinity upon mutation of this amino acid show the good correspondence between the predicted structures and their biological activities. To our knowledge, this is the first systematic approach that was developed for large TCR repertoire structural modeling.

Prediction of cross-recognition of peptide-HLA A2 by Melan-A-specific cytotoxic T lymphocytes using three-dimensional quantitative structure-activity relationships.

The cross-recognition of peptides by cytotoxic T lymphocytes is a key element in immunology and in particular in peptide based immunotherapy. Here we develop three-dimensional (3D) quantitative structure-activity relationships (QSARs) to predict cross-recognition by Melan-A-specific cytotoxic T lymphocytes of peptides bound to HLA A*0201 (hereafter referred to as HLA A2). First, we predict the structure of a set of self- and pathogen-derived peptides bound to HLA A2 using a previously developed ab initio structure prediction approach [Fagerberg et al., J. Mol. Biol., 521-46 (2006)]. Second, shape and electrostatic energy calculations are performed on a 3D grid to produce similarity matrices which are combined with a genetic neural network method [So et al., J. Med. Chem., 4347-59 (1997)] to generate 3D-QSAR models. The models are extensively validated using several different approaches. During the model generation, the leave-one-out cross-validated correlation coefficient (q (2)) is used as the fitness criterion and all obtained models are evaluated based on their q (2) values. Moreover, the best model obtained for a partitioned data set is evaluated by its correlation coefficient (r = 0.92 for the external test set). The physical relevance of all models is tested using a functional dependence analysis and the robustness of the models obtained for the entire data set is confirmed using y-randomization. Finally, the validated models are tested for their utility in the setting of rational peptide design: their ability to discriminate between peptides that only contain side chain substitutions in a single secondary anchor position is evaluated. In addition, the predicted cross-recognition of the mono-substituted peptides is confirmed experimentally in chromium-release assays. These results underline the utility of 3D-QSARs in peptide mimetic design and suggest that the properties of the unbound epitope are sufficient to capture most of the information to determine the cross-recognition.

Genome-wide profiling of the cardiac transcriptome after myocardial infarction identifies novel heart-specific long non-coding RNAs.

AIM: Heart disease is recognized as a consequence of dysregulation of cardiac gene regulatory networks. Previously, unappreciated components of such networks are the long non-coding RNAs (lncRNAs). Their roles in the heart remain to be elucidated. Thus, this study aimed to systematically characterize the cardiac long non-coding transcriptome post-myocardial infarction and to elucidate their potential roles in cardiac homoeostasis. METHODS AND RESULTS: We annotated the mouse transcriptome after myocardial infarction via RNA sequencing and ab initio transcript reconstruction, and integrated genome-wide approaches to associate specific lncRNAs with developmental processes and physiological parameters. Expression of specific lncRNAs strongly correlated with defined parameters of cardiac dimensions and function. Using chromatin maps to infer lncRNA function, we identified many with potential roles in cardiogenesis and pathological remodelling. The vast majority was associated with active cardiac-specific enhancers. Importantly, oligonucleotide-mediated knockdown implicated novel lncRNAs in controlling expression of key regulatory proteins involved in cardiogenesis. Finally, we identified hundreds of human orthologues and demonstrate that particular candidates were differentially modulated in human heart disease. CONCLUSION: These findings reveal hundreds of novel heart-specific lncRNAs with unique regulatory and functional characteristics relevant to maladaptive remodelling, cardiac function and possibly cardiac regeneration. This new class of molecules represents potential therapeutic targets for cardiac disease. Furthermore, their exquisite correlation with cardiac physiology renders them attractive candidate biomarkers to be used in the clinic.

The antibody response in mice to carrier-free synthetic polymers of Plasmodium falciparum circumsporozoite repetitive epitope is I-Ab-restricted: possible implications for malaria vaccines.

The immunogenicity of a novel synthetic peptide consisting of an average of 40 (Asn-Ala-Asn-Pro) repeats of the circumsporozoite protein of Plasmodium falciparum, (NANP)40, was studied in mice without using any carrier proteins. First, high titers of anti-(NANP)40 antibodies could be obtained after immunization of C57BL/6 mice. These antibodies also reacted with an extract of mosquitoes infected with P. falciparum sporozoites. C57BL/6 nu/nu mice did not produce antibodies against (NANP)40. Secondly, when 14 strains of mice with nine different H-2 haplotypes were immunized with (NANP)40 without carrier, only H-2b mice were found to produce anti-(NANP)40 antibodies, whereas all non-H-2b mice were consistently unresponsive. This response was demonstrated to be I-A-linked by using recombinant and mutant mice. I-Ab [B10.A(5R)] mice produced anti-(NANP)40 antibodies as well as H-2b inbred mice. B6CH-2bm12 I-Ab-mutant mice showed only a very low response. Third, the antibody response against (NANP)40 could be induced in nonresponder mice by immunization with the peptide coupled to a carrier protein. In view of the existence of such an exceptional H-2b restriction in the response to sporozoite synthetic peptides in mice, the triggering of peptide-specific T cell responses in humans receiving sporozoite malaria vaccines might be difficult to achieve.

Traitement du décollement de rétine du pseudophaque: étude rétrospective comparant la vitrectomie sans indentation au traitement ab-externo [Vitrectomy without scleral buckle versus ab-externo approach for pseudophakic retinal detachment: comparative retrospective study].

PURPOSE: Retinal detachment (RD) is a major complication of cataract surgery, which can be treated by either primary vitrectomy without indentation or the scleral buckling procedure. The aim of this study is to compare the results of these two techniques for the treatment of pseudophakic RD. PATIENTS AND METHODS: The charts of 40 patients (40 eyes) treated with scleral buckling for a primary pseudophakic RD were retrospectively studied and compared to the charts of 32 patients (32 eyes) treated with primary vitrectomy without scleral buckle during the same period by the same surgeons. To obtain comparable samples, patients with giant retinal tears, vitreous hemorrhage, and severe preoperative proliferative vitreoretinopathy (PVR) were not included. Minimal follow-up was 6 months. RESULTS: The primary success rate was 84% in the vitrectomy group and 82.5% in the ab-externo group. Final anatomical success was observed in 100% of cases in the vitrectomy group and in 95% of cases in the ab-externo group. Final visual acuity was 0.5 or better in 44% of cases in the vitrectomy group and 37.5% in the ab-externo group. The duration of the surgery was significantly lower in the ab-externo group, whereas the hospital stay tended to be lower in the vitrectomy group. In the vitrectomy group, postoperative PVR developed in 3 eyes and new or undetected breaks were responsible for failure of the initial procedure in 2 eyes. CONCLUSION: Primary vitrectomy appears to be as effective as scleral buckling procedures for the treatment of pseudophakic RD.

Phase-I/II radio-immunotherapy study with Iodine-131-labeled anti-CEA monoclonal antibody F6 F(ab')2 in patients with non-resectable liver metastases from colorectal cancer.

Experimental studies in nude mice with human colon-carcinoma grafts demonstrated the therapeutic efficiency of F(ab')2 fragments to carcinoembryonic antigen (CEA) labeled with a high dose of 131Iodine. A phase I/II study was designed to determine the maximum tolerated dose of 131I-labeled F(ab')2 fragments (131I-F(ab')2) from anti-CEA monoclonal antibody F6, its limiting organ toxicity and tumor uptake. Ten patients with non-resectable liver metastases from colorectal cancer (9 detected by CT scan and 1 by laparotomy) were treated with 131I-F(ab')2, doses ranging from 87 mCi to 300 mCi for the first 5 patients, with a constant 300-mCi dose for the last 5 patients. For all the patients, autologous bone marrow was harvested and stored before treatment. Circulating CEA ranged from 2 to 126 ng/ml. No severe adverse events were observed during or immediately following infusion of therapeutic doses. The 9 patients with radiologic evidence of liver metastases showed uptake of 131I-F(ab')2 in the metastases, as observed by single-photon-emission tomography. The only toxicity was hematologic, and no severe aplasia was observed when up to 250 mCi was infused. At the 300-mCi dose, 5 out of 6 patients presented grade-3 or -4 hematologic toxicity, with a nadir for neutrophils and thrombocytes ranging from 25 to 35 days after infusion. In these 5 cases, bone marrow was re-infused. No clinical complications were observed during aplasia. The tumor response could be evaluated in 9 out of 10 patients. One patient showed a partial response of one small liver metastasis (2 cm in diameter) and a stable evolution of the other metastases, 2 patients had stable disease, and 6 showed tumor progression at the time of evaluation (2 or 3 months after injection) by CT scan. This phase-I/II study demonstrated that a dose of 300 mCi of 131I-F(ab')2 from the anti-CEA Mab F6 is well tolerated with bone-marrow rescue, whereas a dose of 200 mCi can be infused without severe bone-marrow toxicity.

Adjuvant radioimmunotherapy trial with iodine-131-labeled anti-carcinoembryonic antigen monoclonal antibody F6 F(ab')2 after resection of liver metastases from colorectal cancer.

PURPOSE: To evaluate the feasibility of radioimmunotherapy (RIT) with radiolabeled anti-carcinoembryonic antigen antibodies after complete resection of liver metastases (LM) from colorectal cancer. Patients and Methods: Twenty-two patients planned for surgery of one to four LM received a preoperative diagnostic dose of a 131I-F(ab')2-labeled anti-carcinoembryonic antigen monoclonal antibody F6 (8-10 mCi/5 mg). 131I-F(ab')2 uptake was analyzed using direct radioactivity counting, and tumor-to-normal liver ratios were recorded. Ten patients with tumor-to-normal liver ratios of >5 and three others were treated with a therapeutic injection [180-200 mCi 131I/50 mg F(ab')2] 30 to 64 days after surgery. RESULTS: Median 131I-F(ab')2 immunoreactivity in patient serum remained at 91% of initial values for up to 96 hours after injection. The main and dose-limiting-toxicity was hematologic, with 92% and 85% grades 3 to 4 neutropenia and thrombocytopenia, respectively. Complete spontaneous recovery occurred in all patients. No human anti-mouse antibody response was observed after the diagnosis dose; however, 10 of the 13 treated patients developed human anti-mouse antibody approximately 3 months later. Two treated patients presented extrahepatic metastases at the time of RIT (one bone and one abdominal node) and two relapsed within 3 months of RIT (one in the lung and the other in the liver). Two patients are still alive, and one of these is disease-free at 93 months after resection. At a median follow-up of 127 months, the median disease-free survival is 12 months and the median overall survival is 50 months. CONCLUSION: RIT is feasible in an adjuvant setting after complete resection of LM from colorectal cancer and should be considered for future trials, possibly in combination with chemotherapy, because of the generally poor prognosis of these patients.