EADock : design of a new molecular docking algorithm and some of its applications

3 Summary 3. 1 English The pharmaceutical industry has been facing several challenges during the last years, and the optimization of their drug discovery pipeline is believed to be the only viable solution. High-throughput techniques do participate actively to this optimization, especially when complemented by computational approaches aiming at rationalizing the enormous amount of information that they can produce. In siiico techniques, such as virtual screening or rational drug design, are now routinely used to guide drug discovery. Both heavily rely on the prediction of the molecular interaction (docking) occurring between drug-like molecules and a therapeutically relevant target. Several softwares are available to this end, but despite the very promising picture drawn in most benchmarks, they still hold several hidden weaknesses. As pointed out in several recent reviews, the docking problem is far from being solved, and there is now a need for methods able to identify binding modes with a high accuracy, which is essential to reliably compute the binding free energy of the ligand. This quantity is directly linked to its affinity and can be related to its biological activity. Accurate docking algorithms are thus critical for both the discovery and the rational optimization of new drugs. In this thesis, a new docking software aiming at this goal is presented, EADock. It uses a hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with .the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 R around the center of mass of the ligand position in the crystal structure, and conversely to other benchmarks, our algorithms was fed with optimized ligand positions up to 10 A root mean square deviation 2MSD) from the crystal structure. This validation illustrates the efficiency of our sampling heuristic, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best-ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures in this benchmark could be explained by the presence of crystal contacts in the experimental structure. EADock has been used to understand molecular interactions involved in the regulation of the Na,K ATPase, and in the activation of the nuclear hormone peroxisome proliferatoractivated receptors a (PPARa). It also helped to understand the action of common pollutants (phthalates) on PPARy, and the impact of biotransformations of the anticancer drug Imatinib (Gleevec®) on its binding mode to the Bcr-Abl tyrosine kinase. Finally, a fragment-based rational drug design approach using EADock was developed, and led to the successful design of new peptidic ligands for the a5ß1 integrin, and for the human PPARa. In both cases, the designed peptides presented activities comparable to that of well-established ligands such as the anticancer drug Cilengitide and Wy14,643, respectively. 3.2 French Les récentes difficultés de l'industrie pharmaceutique ne semblent pouvoir se résoudre que par l'optimisation de leur processus de développement de médicaments. Cette dernière implique de plus en plus. de techniques dites "haut-débit", particulièrement efficaces lorsqu'elles sont couplées aux outils informatiques permettant de gérer la masse de données produite. Désormais, les approches in silico telles que le criblage virtuel ou la conception rationnelle de nouvelles molécules sont utilisées couramment. Toutes deux reposent sur la capacité à prédire les détails de l'interaction moléculaire entre une molécule ressemblant à un principe actif (PA) et une protéine cible ayant un intérêt thérapeutique. Les comparatifs de logiciels s'attaquant à cette prédiction sont flatteurs, mais plusieurs problèmes subsistent. La littérature récente tend à remettre en cause leur fiabilité, affirmant l'émergence .d'un besoin pour des approches plus précises du mode d'interaction. Cette précision est essentielle au calcul de l'énergie libre de liaison, qui est directement liée à l'affinité du PA potentiel pour la protéine cible, et indirectement liée à son activité biologique. Une prédiction précise est d'une importance toute particulière pour la découverte et l'optimisation de nouvelles molécules actives. Cette thèse présente un nouveau logiciel, EADock, mettant en avant une telle précision. Cet algorithme évolutionnaire hybride utilise deux pressions de sélections, combinées à une gestion de la diversité sophistiquée. EADock repose sur CHARMM pour les calculs d'énergie et la gestion des coordonnées atomiques. Sa validation a été effectuée sur 37 complexes protéine-ligand cristallisés, incluant 11 protéines différentes. L'espace de recherche a été étendu à une sphère de 151 de rayon autour du centre de masse du ligand cristallisé, et contrairement aux comparatifs habituels, l'algorithme est parti de solutions optimisées présentant un RMSD jusqu'à 10 R par rapport à la structure cristalline. Cette validation a permis de mettre en évidence l'efficacité de notre heuristique de recherche car des modes d'interactions présentant un RMSD inférieur à 2 R par rapport à la structure cristalline ont été classés premier pour 68% des complexes. Lorsque les cinq meilleures solutions sont prises en compte, le taux de succès grimpe à 78%, et 92% lorsque la totalité de la dernière génération est prise en compte. La plupart des erreurs de prédiction sont imputables à la présence de contacts cristallins. Depuis, EADock a été utilisé pour comprendre les mécanismes moléculaires impliqués dans la régulation de la Na,K ATPase et dans l'activation du peroxisome proliferatoractivated receptor a (PPARa). Il a également permis de décrire l'interaction de polluants couramment rencontrés sur PPARy, ainsi que l'influence de la métabolisation de l'Imatinib (PA anticancéreux) sur la fixation à la kinase Bcr-Abl. Une approche basée sur la prédiction des interactions de fragments moléculaires avec protéine cible est également proposée. Elle a permis la découverte de nouveaux ligands peptidiques de PPARa et de l'intégrine a5ß1. Dans les deux cas, l'activité de ces nouveaux peptides est comparable à celles de ligands bien établis, comme le Wy14,643 pour le premier, et le Cilengitide (PA anticancéreux) pour la seconde.

Awakening dormant haematopoietic stem cells.

Haematopoietic stem cells (HSCs) in mouse bone marrow are located in specialized niches as single cells. During homeostasis, signals from this environment keep some HSCs dormant, which preserves long-term self-renewal potential, while other HSCs actively self renew to maintain haematopoiesis. In response to haematopoietic stress, dormant HSCs become activated and rapidly replenish the haematopoietic system. Interestingly, three factors - granulocyte colony-stimulating factor, interferon-alpha and arsenic trioxide - have been shown to efficiently activate dormant stem cells and thereby could break their resistance to anti-proliferative chemotherapeutics. Thus, we propose that two-step strategies could target resistant leukaemic stem cells by priming tumours with activators of dormancy followed by chemotherapy or targeted therapies.

Gene expression profiling reveals consistent differences between clinical samples of human leukaemias and their model cell lines.

Microarray gene expression profiles of fresh clinical samples of chronic myeloid leukaemia in chronic phase, acute promyelocytic leukaemia and acute monocytic leukaemia were compared with profiles from cell lines representing the corresponding types of leukaemia (K562, NB4, HL60). In a hierarchical clustering analysis, all clinical samples clustered separately from the cell lines, regardless of leukaemic subtype. Gene ontology analysis showed that cell lines chiefly overexpressed genes related to macromolecular metabolism, whereas in clinical samples genes related to the immune response were abundantly expressed. These findings must be taken into consideration when conclusions from cell line-based studies are extrapolated to patients.

Telomere length dynamics in normal individuals and in patients with hematopoietic stem cell-associated disorders.

The telomere length in nucleated peripheral blood (PB) cells indirectly reflects the mitotic history of their precursors: the hematopoietic stem cells (HSCs). The average length of telomeres in PB leukocytes can be measured using fluorescence in situ hybridization and flow cytometry (flow FISH). We previously used flow FISH to characterize the age-related turnover of HSCs in healthy individuals. In this review, we describe results of recent flow FISH studies in patients with selected hematopoietic stem cell-associated disorders: chronic myelogenous leukemia (CML) and several bone marrow failure syndromes. CML is characterized by a marked expansion of myeloid Philadelphia chromosome positive (Ph+) cells. Nevertheless, nonmalignant (Ph-) HSCs typically coexist in the bone marrow of CML patients. We analyzed the telomere length in > 150 peripheral blood leukocytes (PBLs) and bone marrow samples of patients with CML as well as samples of Ph- T-lymphocytes. Compared to normal controls, the overall telomere fluorescence in PBLs of patients with CML was significantly reduced. However, no telomere shortening was observed in Ph- T-lymphocytes. Patients in late chronic phase (CP) had significantly shorter telomeres than those assessed earlier in CP. Our data suggest that progressive telomere shortening is correlated with disease progression in CML. Within the group of patients with bone marrow failure syndromes, we only found significantly shortened telomeres (compared to age-adjusted controls) in granulocytes from patients with aplastic anemia (AA). Strikingly, the telomere length in granulocytes from AA patients who had recovered after immunosuppressive therapy (recAA) did not differ significantly from controls, whereas untreated patients and nonresponders with persistent severe pancytopenia (sAANR) showed marked and significant telomere shortening compared to healthy donors and patients with recAA. Furthermore, an inverse correlation between age-adjusted telomere length and peripheral blood counts was found in support of a model in which the degree of cytopenia and the amount of telomere shortening are correlated. These results support the concept of extensive proliferation of HSCs in subgroups of AA patients and suggest a potential use of telomere-length measurements as a prognostic tool in this group of disorders as well.

BCR and TLR signalling pathways are recurrently targeted by genetic changes in splenic marginal zone lymphomas.

The genetics and pathogenesis of splenic marginal zone lymphoma are poorly understood. The lymphoma lacks chromosome translocation, and ~30% of cases are featured by 7q deletion, but the gene targeted by the deletion is unknown. A recent study showed inactivation of A20, a 'global' NF-kB negative regulator, in 1 of 12 splenic marginal zone lymphoma. To investigate further whether deregulation of the NF-kB pathway plays a role in the pathogenesis of splenic marginal zone lymphoma, we screened several NF-kB regulators for genetic changes by PCR and sequencing. Somatic mutations were found in A20 (6/46=13%), MYD88 (6/46=13%), CARD11 (3/34=8.8%), but not in CD79A, CD79B and ABIN1. Interestingly, these genetic changes are largely mutually exclusive from each other and MYD88 mutation was also mutually exclusive from 7q deletion. These results strongly suggest that deregulation of the TLR (toll like receptor) and BCR (B-cell receptor) signalling pathway may play an important role in the pathogenesis of splenic marginal zone lymphoma.

Expression of the proto-oncogene c-fos in three-dimensional fetal brain cell cultures and the lack of correlation with maturation-inducing stimuli.

Previous work has shown that aggregating fetal brain cell cultures are able to attain a highly differentiated state, and that their development is greatly enhanced by growth and/or differentiation factors such as epidermal growth factor (EGF), basic fibroblast growth factor (bFGF), and the protein kinase C-activating tumor promoter mezerein. The present study shows that in these 3-dimensional cultures the peptide growth factors EGF and bFGF as well as mezerein are able to induce the expression of the proto-oncogene c-fos. This induction was rapid and transient, in good agreement with observations reported from a wide variety of cell types in vitro. The maximal levels of c-fos mRNA found after stimulation were low in immature cultures and increased greatly as maturation progressed. Of the three factors tested, mezerein was the most potent inducer of c-fos. In contrast to the peptide growth factors EGF and bFGF which were found to induce c-fos only in glial cells, mezerein was stimulatory in glial cells as well as in neurons. A similar cell type specificity has been observed previously for the maturation-enhancing response in immature aggregate cultures. However, in the present study no correlation was found between the degree of c-fos induction and the extent of the maturation-enhancing stimulation. Immature cultures known to be most sensitive and responsive to these maturation-enhancing agents required relatively high doses of peptide growth factors for the induction of c-fos, and the maximal levels of c-fos mRNA elicited were much lower than those in differentiated cultures which did not show any long-term response to these stimuli.(ABSTRACT TRUNCATED AT 250 WORDS)

Calcineurin inhibitors activate the proto-oncogene Ras and promote protumorigenic signals in renal cancer cells.

The development of cancer is a major problem in immunosuppressed patients, particularly after solid organ transplantation. We have recently shown that calcineurin inhibitors (CNI) used to treat transplant patients may play a critical role in the rapid progression of renal cancer. To examine the intracellular signaling events for CNI-mediated direct tumorigenic pathway(s), we studied the effect of CNI on the activation of proto-oncogenic Ras in human normal renal epithelial cells (REC) and renal cancer cells (786-0 and Caki-1). We found that CNI treatment significantly increased the level of activated GTP-bound form of Ras in these cells. In addition, CNI induced the association of Ras with one of its effector molecules, Raf, but not with Rho and phosphatidylinositol 3-kinase; CNI treatment also promoted the phosphorylation of the Raf kinase inhibitory protein and the downregulation of carabin, all of which may lead to the activation of the Ras-Raf pathway. Blockade of this pathway through either pharmacologic inhibitors or gene-specific small interfering RNA significantly inhibited CNI-mediated augmented proliferation of renal cancer cells. Finally, it was observed that CNI treatment increased the growth of human renal tumors in vivo, and the Ras-Raf pathway is significantly activated in the tumor tissues of CNI-treated mice. Together, targeting the Ras-Raf pathway may prevent the development/progression of renal cancer in CNI-treated patients.

Post-transcriptional down-regulation of expression of transcription factor NF1 by Ha-ras oncogene.

NF1 is a family of polypeptides that binds to discrete DNA motifs and plays varying roles in the regulation of gene expression. These polypeptides are also thought to mediate the expression of differentiation-specific markers such as adipocyte and mammary cell type-specific genes. The expression of a number of cellular differentiation-specific markers is down-regulated during neoplastic transformation. We therefore investigated whether oncogenic transformation interferes with the action of NF1. Stable transfection of activated Ha-ras into a number of murine cells correlated with a down-regulation of the expression of the NF1 genes NF1/CTF and NF1/X. The down-regulation was not at the transcriptional level but at the level of stability of the NF1 mRNAs. The level of the DNA binding activity of the NF1 proteins was also reduced in Ha-v-ras-transformed cells, and the expression of a gene that depends on this family of transcription factors was specifically repressed. These results demonstrate that an activated Ha-ras-induced pathway destabilizes the half-life of mRNAs encoding specific members in the NF1 family of transcription factors, which leads to a decrease in NF1-dependent gene expression.

Induction of human papillomavirus oncogene-specific CD8 T-cell effector responses in the genital mucosa of vaccinated mice.

Cervical cancer, the second leading cause of cancer mortality in women worldwide, results from infection with a subset of human papillomaviruses (HPV), HPV-16 being the most prevalent type. The available prophylactic vaccines are an effective strategy to prevent this cancer in the long term. However, they only target 70-80% of all cervical cancers and cannot control existing HPV infections and associated lesions. Therapeutic vaccines are thus necessary for women who cannot benefit from prophylactic vaccination. Induction of protective immune responses in the genital mucosa (GM) may be crucial for efficacy of HPV therapeutic vaccines. We report here that mice that received a single subcutaneous (s.c.) vaccination of an adjuvanted long synthetic HPV16 E7(1-98) polypeptide showed induction of 100% tumor protection against s.c. TC-1 tumors and that tumor regression was mainly provided by CD8 T cells. In vivo cytotoxic assay revealed high E7-specific cytolytic T lymphocytes activity in spleen and in genital draining lymph nodes (LN), and E7-specific CD8 T cells could be detected in GM by tetramer staining. More importantly, high-avidity E7-specific INF-gamma secreting CD8 T cells were induced not only in blood, spleen and LN but also in GM of vaccinated mice, thus providing evidence that a parenteral vaccination may be sufficient to provide regression of genital tumors. In addition, there was no correlation between the responses measured in blood with those measured in GM, highlighting the necessity and relevance to determine the immune responses in the mucosa where HPV-tumors reside.

The Oncogene ATF3 Is Potentiated by Cyclosporine A and Ultraviolet Light A.

Cutaneous squamous cell carcinoma (SCC) represents the most important cutaneous complication following organ transplantation. It develops mostly on sun-exposed areas. A recent study showed the role of activating transcription factor 3 (ATF3) in SCC development following treatment with calcineurin inhibitors. It has been reported that ATF3, which may act as an oncogene, is under negative calcineurin/nuclear factor of activated T cells (NFAT) control and is upregulated by calcineurin inhibitors. Still, these findings do not fully explain the preferential appearance of SCC on chronically sun-damaged skin. We analyzed the influence of UV radiation on ATF3 expression and its potential role in SCC development. We found that ATF3 is a specifically induced AP1 member in SCC of transplanted patients. Its expression was strongly potentiated by combination of cyclosporine A and UVA treatment. UVA induced ATF3 expression through reactive oxygen species-mediated nuclear factor erythroid 2-related factor 2 (NRF2) activation independently of calcineurin/NFAT inhibition. Activated NRF2 directly binds to ATF3 promoter, thus inducing its expression. These results demonstrate two mechanisms that independently induce and, when combined together, potentiate the expression of ATF3, which may then force SCC development. Taking into account the previously defined role of ATF3 in the SCC development, these findings may provide an explanation and a mechanism for the frequently observed burden on SCCs on sun-exposed areas of the skin in organ transplant recipients treated by calcineurin inhibitors.

The Oncogene Metadherin Modulates the Apoptotic Pathway Based on the Tumor Necrosis Factor Superfamily Member TRAIL (Tumor Necrosis Factor-related Apoptosis-inducing Ligand) in Breast Cancer.

Metadherin (MTDH), the newly discovered gene, is overexpressed in more than 40% of breast cancers. Recent studies have revealed that MTDH favors an oncogenic course and chemoresistance. With a number of breast cancer cell lines and breast tumor samples, we found that the relative expression of MTDH correlated with tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) sensitivity in breast cancer. In this study, we found that knockdown of endogenous MTDH cells sensitized the MDA-MB-231 cells to TRAIL-induced apoptosis both in vitro and in vivo. Conversely, stable overexpression of MTDH in MCF-7 cells enhanced cell survival with TRAIL treatment. Mechanically, MTDH down-regulated caspase-8, decreased caspase-8 recruitment into the TRAIL death-inducing signaling complex, decreased caspase-3 and poly(ADP-ribose) polymerase-2 processing, increased Bcl-2 expression, and stimulated TRAIL-induced Akt phosphorylation, without altering death receptor status. In MDA-MB-231 breast cancer cells, sensitization to TRAIL upon MTDH down-regulation was inhibited by the caspase inhibitor Z-VAD-fmk (benzyloxycarbonyl-VAD-fluoromethyl ketone), suggesting that MTDH depletion stimulates activation of caspases. In MCF-7 breast cancer cells, resistance to TRAIL upon MTDH overexpression was abrogated by depletion of Bcl-2, suggesting that MTDH-induced Bcl-2 expression contributes to TRAIL resistance. We further confirmed that MTDH may control Bcl-2 expression partly by suppressing miR-16. Collectively, our results point to a protective function of MTDH against TRAIL-induced death, whereby it inhibits the intrinsic apoptosis pathway through miR-16-mediated Bcl-2 up-regulation and the extrinsic apoptosis pathway through caspase-8 down-regulation.

Combined loss of the BH3-only proteins Bim and Bmf restores B-cell development and function in TACI-Ig transgenic mice.

Terminal differentiation of B cells depends on two interconnected survival pathways, elicited by the B-cell receptor (BCR) and the BAFF receptor (BAFF-R), respectively. Loss of either signaling pathway arrests B-cell development. Although BCR-dependent survival depends mainly on the activation of the v-AKT murine thymoma viral oncogene homolog 1 (AKT)/PI3-kinase network, BAFF/BAFF-R-mediated survival engages non-canonical NF-κB signaling as well as MAPK/extracellular-signal regulated kinase and AKT/PI3-kinase modules to allow proper B-cell development. Plasma cell survival, however, is independent of BAFF-R and regulated by APRIL that signals NF-κB activation via alternative receptors, that is, transmembrane activator and CAML interactor (TACI) or B-cell maturation (BCMA). All these complex signaling events are believed to secure survival by increased expression of anti-apoptotic B-cell lymphoma 2 (Bcl2) family proteins in developing and mature B cells. Curiously, how lack of BAFF- or APRIL-mediated signaling triggers B-cell apoptosis remains largely unexplored. Here, we show that two pro-apoptotic members of the 'Bcl2 homology domain 3-only' subgroup of the Bcl2 family, Bcl2 interacting mediator of cell death (Bim) and Bcl2 modifying factor (Bmf), mediate apoptosis in the context of TACI-Ig overexpression that effectively neutralizes BAFF as well as APRIL. Surprisingly, although Bcl2 overexpression triggers B-cell hyperplasia exceeding the one observed in Bim(-/-)Bmf(-/-) mice, Bcl2 transgenic B cells remain susceptible to the effects of TACI-Ig expression in vivo, leading to ameliorated pathology in Vav-Bcl2 transgenic mice. Together, our findings shed new light on the molecular machinery restricting B-cell survival during development, normal homeostasis and under pathological conditions. Our data further suggest that Bcl2 antagonists might improve the potency of BAFF/APRIL-depletion strategies in B-cell-driven pathologies.

Redirection of T cells by delivering a transgenic mouse-derived MDM2 tumor antigen-specific TCR and its humanized derivative is governed by the CD8 coreceptor and affects natural human TCR expression.

Retroviral transfer of T cell antigen receptor (TCR) genes selected by circumventing tolerance to broad tumor- and leukemia-associated antigens in human leukocyte antigen (HLA)-A*0201 (A2.1) transgenic (Tg) mice allows the therapeutic reprogramming of human T lymphocytes. Using a human CD8 x A2.1/Kb mouse derived TCR specific for natural peptide-A2.1 (pA2.1) complexes comprising residues 81-88 of the human homolog of the murine double-minute 2 oncoprotein, MDM2(81-88), we found that the heterodimeric CD8 alpha beta coreceptor, but not normally expressed homodimeric CD8 alpha alpha, is required for tetramer binding and functional redirection of TCR- transduced human T cells. CD8+T cells that received a humanized derivative of the MDM2 TCR bound pA2.1 tetramers only in the presence of an anti-human-CD8 anti-body and required more peptide than wild-type (WT) MDM2 TCR+T cells to mount equivalent cytotoxicity. They were, however, sufficiently effective in recognizing malignant targets including fresh leukemia cells. Most efficient expression of transduced TCR in human T lymphocytes was governed by mouse as compared to human constant (C) alphabeta domains, as demonstrated with partially humanized and murinized TCR of primary mouse and human origin, respectively. We further observed a reciprocal relationship between the level of Tg WT mouse relative to natural human TCR expression, resulting in T cells with decreased normal human cell surface TCR. In contrast, natural human TCR display remained unaffected after delivery of the humanized MDM2 TCR. These results provide important insights into the molecular basis of TCR gene therapy of malignant disease.

Remyelination in vitro following protein kinase C activator-induced demyelination.

In previous work we found that mezerein, a C kinase activator, as well as basic fibroblast growth factor (FGF-2) induce demyelination and partial oligodendrocyte dedifferentiation in highly differentiated aggregating brain cell cultures. Here we show that following protein kinase C activator-induced demyelination, effective remyelination occurs. We found that mezerein or FGF-2 caused a transient increase in DNA synthesis following a pronounced decrease of the myelin markers myelin basic protein and 2',3'-cyclic nucleotide 3'-phosphohydrolase. Both oligodendrocytes and astrocytes were involved in this mitogenic response. Within 17 days after demyelination, myelin was restored to the level of the untreated controls. Transient mitotic activity was indispensable for remyelination. The present results suggest that myelinating oligodendrocytes retain the capacity to reenter the cell cycle, and that this plasticity is important for the regeneration of the oligodendrocyte lineage and remyelination. Although it cannot be excluded that a quiescent population of oligodendrocyte precursor cells was present in the aggregates and able to proliferate, differentiate and remyelinate, we could not find evidence supporting this view.