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.

Mobilization of hemopoietic stem cells with high-dose methotrexate plus granulocyte-colony-stimulating factor in patients with primary central nervous system lymphoma.

BACKGROUND: High-dose therapy with autologous stem cell support after standard dose induction is a promising approach for therapy of primary central nervous system lymphoma (PCNSL). High-dose methotrexate (HD-MTX) is a standard drug for induction of PCNSL; however, data about the capacity of HD-MTX plus granulocyte-colony-stimulating factor (G-CSF) to mobilize hemopoietic progenitors are lacking. STUDY DESIGN AND METHODS: This investigation describes the data from stem cell mobilization and apheresis procedures after one or two cycles of HD-MTX for induction of PCNSL within the East German Study Group for Haematology and Oncology 053 trial. Eligible patients proceeded to high-dose busulfan/thiotepa after induction therapy and mobilization. RESULTS: Data were available from nine patients with a median age of 58 years. The maximal CD34+ cell count per microL of blood after the first course of HD-MTX was 13.89 (median). Determination was repeated in six patients after the second course with a significantly higher median CD34+ cell count of 33.69 per microL. Five patients required two apheresis procedures and in four patients a single procedure was sufficient. The total yield of CD34+ cells per kg of body weight harvested by one or two leukapheresis procedures was 6.60 x 10(6) (median; range, 2.68 x 10(6)-15.80 x 10(6)). The yield of CD34+ cells exceeded the commonly accepted lower threshold of 3 x 10(6) cells per kg of body weight in eight of nine cases. Even in the ninth, hemopoietic recovery after stem cell reinfusion was rapid and safe. CONCLUSION: HD-MTX plus G-CSF is a powerful combination for stem cell mobilization in patients with PCNSL and permits safe conduction of time-condensed and dose-intense protocols with high-dose therapy followed by stem cell reinfusion after HD-MTX induction.

Molecular and social regulation of worker division of labour in fire ants.

Reproductive and worker division of labour (DOL) is a hallmark of social insect societies. Despite a long-standing interest in worker DOL, the molecular mechanisms regulating this process have only been investigated in detail in honey bees, and little is known about the regulatory mechanisms operating in other social insects. In the fire ant Solenopsis invicta, one of the most studied ant species, workers are permanently sterile and the tasks performed are modulated by the worker's internal state (age and size) and the outside environment (social environment), which potentially includes the effect of the queen presence through chemical communication via pheromones. However, the molecular mechanisms underpinning these processes are unknown. Using a whole-genome microarray platform, we characterized the molecular basis for worker DOL and we explored how a drastic change in the social environment (i.e. the sudden loss of the queen) affects global gene expression patterns of worker ants. We identified numerous genes differentially expressed between foraging and nonforaging workers in queenright colonies. With a few exceptions, these genes appear to be distinct from those involved in DOL in bees and wasps. Interestingly, after the queen was removed, foraging workers were no longer distinct from nonforaging workers at the transcriptomic level. Furthermore, few expression differences were detected between queenright and queenless workers when we did not consider the task performed. Thus, the social condition of the colony (queenless vs. queenright) appears to impact the molecular pathways underlying worker task performance, providing strong evidence for social regulation of DOL in S. invicta.

Management of Pharyngeal Fistulas After Anterior Cervical Spine Surgery: A Treatment Algorithm for Severe Complications.

STUDY DESIGN:: Retrospective database- query to identify all anterior spinal approaches. OBJECTIVES:: To assess all patients with pharyngo-cutaneous fistulas after anterior cervical spine surgery. SUMMARY OF BACKGROUND DATA:: Patients treated in University of Heidelberg Spine Medical Center, Spinal Cord Injury Unit and Department of Otolaryngology (Germany), between 2005 and 2011 with the diagnosis of pharyngo-cutaneous fistulas. METHODS:: We conducted a retrospective study on 5 patients between 2005 and 2011 with PCF after ACSS, their therapy management and outcome according to radiologic data and patient charts. RESULTS:: Upon presentation 4 patients were paraplegic. 2 had PCF arising from one piriform sinus, two patients from the posterior pharyngeal wall and piriform sinus combined and one patient only from the posterior pharyngeal wall. 2 had previous unsuccessful surgical repair elsewhere and 1 had prior radiation therapy. In 3 patients speech and swallowing could be completely restored, 2 patients died. Both were paraplegic. The patients needed an average of 2-3 procedures for complete functional recovery consisting of primary closure with various vascularised regional flaps and refining laser procedures supplemented with negative pressure wound therapy where needed. CONCLUSION:: Based on our experience we are able to provide a treatment algorithm that indicates that chronic as opposed to acute fistulas require a primary surgical closure combined with a vascularised flap that should be accompanied by the immediate application of a negative pressure wound therapy. We also conclude that particularly in paraplegic patients suffering this complication the risk for a fatal outcome is substantial.

Mate availability and male dispersal in the Argentine ant Linepithema humile (Mayr) (=Iridomyrmex humilis)

In the Argentine ant Linepithema humile (=Iridomyrmex humilis) only males disperse whereas female sexuals (unmated winged queens) stay in their mother nest where they mate. This study investigated (1) whether dispersing males are accepted into foreign colonies, (2) whether they can mate with resident female sexuals, and (3) whether the propensity of males to disperse is affected by the expectation of mating in their mother nest. Field experiments demonstrated that males were accepted into foreign colonies only when these colonies contained female sexuals or queen pupae. Before and after the time of (sic) female sexuals, workers attacked and killed most of the foreign males. Laboratory experiments snowed that males that successfully enter foreign colonies can mate with resident female sexuals. The propensity of males to disperse was significantly influenced by the presence of female sexuals in their nest. Males were more likely to fly out from colonies containing no female sexuals than from those with them. These results are consistent with males preferentially dispersing when there is little or no opportunity to mate in their mother nest. Thus there are two mating strategies available for males: staying in their mother nest when an opportunity to mate arises or dispersing and attempting to mate in a foreign nest when there are no female sexuals in their mother nest. This latter behaviour could mediate gene flow between colonies and account for the lack of significant inbreeding previously documented in this species.

Unicoloniality, recognition and genetic differentiation in a native Formica ant.

Some ants have an extraordinary form of social organization, called unicoloniality, whereby individuals mix freely among physically separated nests. This mode of social organization has been primarily studied in introduced and invasive ant species, so that the recognition ability and genetic structure of ants forming unicolonial populations in their native range remain poorly known. We investigated the pattern of aggression and the genetic structure of six unicolonial populations of the ant Formica paralugubris at four hierarchical levels: within nests, among nests within the same population, among nests of populations within the Alps or Jura Mountains and among nests of the two mountain ranges. Ants within populations showed no aggressive behaviour, but recognized nonnestmates as shown by longer antennation bouts. Overall, the level of aggression increased with geographic and genetic distance but was always considerably lower than between species. No distinct behavioural supercolony boundaries were found. Our study provides evidence that unicoloniality can be maintained in noninvasive ants despite significant genetic differentiation and the ability to discriminate between nestmates and nonnestmates.

Odorant binding proteins of the red imported fire ant, Solenopsis invicta: an example of the problems facing the analysis of widely divergent proteins.

We describe the odorant binding proteins (OBPs) of the red imported fire ant, Solenopsis invicta, obtained from analyses of an EST library and separate 454 sequencing runs of two normalized cDNA libraries. We identified a total of 18 putative functional OBPs in this ant. A third of the fire ant OBPs are orthologs to honey bee OBPs. Another third of the OBPs belong to a lineage-specific expansion, which is a common feature of insect OBP evolution. Like other OBPs, the different fire ant OBPs share little sequence similarity (∼ 20%), rendering evolutionary analyses difficult. We discuss the resulting problems with sequence alignment, phylogenetic analysis, and tests of selection. As previously suggested, our results underscore the importance for careful exploration of the sensitivity to the effects of alignment methods for data comprising widely divergent sequences.

Safety of a new algorithm for the management of childhood illness (Almanach) to improve quality of care and rational use of drugs

Introduction New evidence from randomized controlled and etiology of fever studies, the availability of reliable RDT for malaria, and novel technologies call for revision of the IMCI strategy. We developed a new algorithm based on (i) a systematic review of published studies assessing the safety and appropriateness of RDT and antibiotic prescription, (ii) results from a clinical and microbiological investigation of febrile children aged <5 years, (iii) international expert IMCI opinions. The aim of this study was to assess the safety of the new algorithm among patients in urban and rural areas of Tanzania.Materials and Methods The design was a controlled noninferiority study. Enrolled children aged 2-59 months with any illness were managed either by a study clinician using the new Almanach algorithm (two intervention health facilities), or clinicians using standard practice, including RDT (two control HF). At day 7 and day 14, all patients were reassessed. Patients who were ill in between or not cured at day 14 were followed until recovery or death. Primary outcome was rate of complications, secondary outcome rate of antibiotic prescriptions.Results 1062 children were recruited. Main diagnoses were URTI 26%, pneumonia 19% and gastroenteritis (9.4%). 98% (531/541) were cured at D14 in the Almanach arm and 99.6% (519/521) in controls. Rate of secondary hospitalization was 0.2% in each. One death occurred in controls. None of the complications was due to withdrawal of antibiotics or antimalarials at day 0. Rate of antibiotic use was 19% in the Almanach arm and 84% in controls.Conclusion Evidence suggests that the new algorithm, primarily aimed at the rational use of drugs, is as safe as standard practice and leads to a drastic reduction of antibiotic use. The Almanach is currently being tested for clinician adherence to proposed procedures when used on paper or a mobile phone

Granulocyte macrophage colony stimulating factor improves mucosal repair in a mouse model of acute colitis (P278)

Background and Aims: Granulocyte-macrophage colonystimulating factor (GM-CSF), a cytokine modulating the number and function of innate immune cells, has been shown to provide symptomatic benefit in some patients with Crohn's disease (CD). Since, it becomes widely appreciated that a timely and spatially regulated action of innate immune cells is critical for tissue regeneration, we tested whether GM-CSF therapy may favours intestinal mucosal repair in the acute mouse model of dextran sulfate sodium (DSS)-induced colitis. Methods: Mice treated with GM-CSF or saline were exposed for 7 days to DSS to induce colitis. On day 5, 7 and 10, mice were subjected to colonoscopy or sacrificed for evaluation of inflammatory reaction and mucosal healing. Results: GM-CSF therapy prevented body weight loss, diarrhea, dampened inflammatory reactions and ameliorated mucosal damages. Mucosal repair improvement in GM-CSF-treated mice was observed from day 7 on both by colonoscopy (ulceration score 1.2}0.3 (GM-CSF-treated) vs 3.1}0.5 (untreated), p = 0.01) and histological analysis (percentage of reepithelialized ulcers 55%}4% (GM-CSF-treated) vs 18%}13% (untreated), p = 0.01). GM-CSF therapy can still improve the colitis when hematopoietic, but not non-hematopoietic cells, are responsive to GM-CSF, as shown in WT→GM-CSFRKO chimeras. Lastly, we observed that GM-CSF-induced promotion of wound healing is associated with a modification of the cellular composition of DSS-induced colonic inflammatory infiltrate, characterized by the reduction of neutrophil numbers and early accumulation of CD11b+Gr1lo myeloid cells. Conclusion: Our study shows that GM-CSF therapy accelerates the complex program leading to tissue repair during acute colitis and suggests that GM-CSF promotion of mucosal repair might contribute to the symptomatic benefits of GM-CSF therapy observed in some CD patients.

Population genomics of eusocial insects: the costs of a vertebrate-like effective population size

The evolution of reproductive division of labour and social life in social insects has lead to the emergence of several life-history traits and adaptations typical of larger organisms: social insect colonies can reach masses of several kilograms, they start reproducing only when they are several years old, and can live for decades. These features and the monopolization of reproduction by only one or few individuals in a colony should affect molecular evolution by reducing the effective population size. We tested this prediction by analysing genome-wide patterns of coding sequence polymorphism and divergence in eusocial vs. noneusocial insects based on newly generated RNA-seq data. We report very low amounts of genetic polymorphism and an elevated ratio of nonsynonymous to synonymous changes - a marker of the effective population size - in four distinct species of eusocial insects, which were more similar to vertebrates than to solitary insects regarding molecular evolutionary processes. Moreover, the ratio of nonsynonymous to synonymous substitutions was positively correlated with the level of social complexity across ant species. These results are fully consistent with the hypothesis of a reduced effective population size and an increased genetic load in eusocial insects, indicating that the evolution of social life has important consequences at both the genomic and population levels.

Division of labor in insect societies : genetic components and physiological regulation

The shift from solitary to social organisms constitutes one of the major transitions in evolution. The highest level of sociality is found in social insects (ants, termites and some species of bees and wasps). Division of labor is central to the organization of insect societies and is thought to be at the root of their ecological success. There are two main levels of division of labor in social insect colonies. The first relates to reproduction and involves the coexistence of queen and worker castes: while reproduction is usually monopolized by one or several queens, functionally sterile workers perform all the tasks to maintain the colony, such as nest building, foraging or brood care. The second level of division of labor, relating to such non-reproductive duties, is characterized by the performance of different tasks or roles by different groups of workers. This PhD aims to better understand the mechanisms underlying division of labor in insect societies, by investigating how genes and physiology influence caste determination and worker behavior in ants. In the first axis of this PhD, we studied the nature of genetic effects on division of labor. We used the Argentine ant Linepithema humile to conduct controlled crosses in the laboratory, which revealed the existence of non-additive genetic effects, such as parent-of-origin and genetic compatibility effects, on caste determination and worker behavior. In the second axis, we focused on the physiological regulation of division of labor. Using Pogonomyrmex seed- harvester ants, we performed experimental manipulation of hibernation, hormonal treatments, gene expression analyses and protein quantification to identify the physiological pathways regulating maternal effects on caste determination. Finally, comparing gene expression between nurses and foragers allowed us to reveal the association between vitellogenin and worker behavior in Pogonomyrmex ants. This PhD provides important insights into the role of genes and physiology in the regulation of division of labor in social insect colonies, helping to better understand the organization, evolution and ecological success of insect societies. - L'une des principales transitions évolutives est le passage de la vie solitaire à la vie sociale. La socialité atteint son paroxysme chez les insectes sociaux que sont les fourmis, les termites et certaines espèces d'abeilles et de guêpes. La division du travail est la clé de voûte de l'organisation de ces sociétés d'insectes et la raison principale de leur succès écologique. La division du travail s'effectue à deux niveaux dans les colonies d'insectes sociaux. Le premier niveau concerne la reproduction et implique la coexistence de deux castes : les reines et les ouvrières. Tandis que la reproduction est le plus souvent monopolisée par une ou plusieurs reines, les ouvrières stériles effectuent les tâches nécessaires au bon fonctionnement de la colonie, telles que la construction du nid, la recherche de nourriture ou le soin au couvain. Le second niveau de division du travail, qui concerne les tâches autres que la reproduction, implique la réalisation de différents travaux par différents groupes d'ouvrières. Le but de ce doctorat est de mieux comprendre les mécanismes sous-jacents de la division du travail dans les sociétés d'insectes en étudiant comment les gènes et la physiologie influencent la détermination de la caste et le comportement des ouvrières chez les fourmis. Dans le premier axe de ce doctorat, nous avons étudié la nature des influences génétiques sur la division du travail. Nous avons utilisé la fourmi d'Argentine, Linepithema humile, pour effectuer des croisements contrôlés en laboratoire. Cette méthode nous a permis de révéler l'existence d'influences génétiques non additives, telles que des influences dépendantes de l'origine parentale ou des effets de compatibilité génétique, sur la détermination de la caste et le comportement des ouvrières. Dans le second axe, nous nous sommes intéressés à la régulation physiologique de la division du travail. Nous avons utilisé des fourmis moissonneuses du genre Pogonomyrmex pour effectuer des hibernations artificieHes, des traitements hormonaux, des analyses d'expression de gènes et des mesures de vitellogénine, ce qui nous a permis d'identifier les mécanismes physiologiques régulant les effets maternels sur la détermination de la caste. Enfin, la comparaison d'expression de gènes entre nourrices et fourrageuses suggère un rôle de la vitellogénine dans la régulation du comportement des ouvrières chez les fourmis moissonneuses. En détaillant les influences des gènes et de la physiologie dans la régulation de la division du travail dans les colonies d'insectes sociaux, ce doctorat fournit d'importantes informations permettant de mieux comprendre l'organisation, l'évolution et le succès écologique des sociétés d'insectes.

Conflict and co-operation in ant societies

Cette thèse se propose d?étudier les forces évolutives qui influencent la bal- ance entre coopération et conflit. Deux exemples sont analysés en détail, les sociétés d?insectes et l?association entre la fourmi Formica exsecta et son sym- bionte intracellulaire Wolbachia. Chacune de ces deux associations animales est caractérisée par la coexistence de coopération et de conflit. Les colonies d?insectes sociaux sont un exemple typique de coopération. Néanmoins, la détermination haplodiploïde du sexe crée des apparentements asymétriques entre les membres de la colonie qui se traduisent par des con- flits entre reines et ouvrières. Un des conflits les plus importants porte sur le sexe-ratio, l?investissement optimal en sexués males et femelles étant équilibré pour la reine, mais biaisé en faveur des femelles pour les ouvrières. Jusqu?ici, les modèles étudiant ce conflit se basaient uniquement sur les asymétries d?ap- parentement. La première partie de cette thèse (Chapitres 1 et 2) propose d?introduire dans des modèles de sexe-ratio les mécanismes proximaux utilisés par les deux parties pour influencer le sexe-ratio. Ces mécanismes, par leurs ef- fets négatifs sur la productivité de la colonie, peuvent entraîner des déviations importantes du sexe-ratio par rapport à la valeur prédite par des modèles clas- siques. De plus, la prise en compte de ces mécanismes nous permet d?étudier dans quelle mesure les parties opposées sont sélectionnées pour influencer le sexe-ratio. Enfin, les modèles génèrent des prédictions quant à l?issue du con- flit, c?est-à-dire un compromis où reines et ouvrières partagent le contrôle de l?investissement. La seconde partie du travail porte sur l?infection de la fourmi Formica ex- secta par le symbionte cytoplasmique Wolbachia. Ce symbionte, présent chez un grand nombre d?espèces d?Arthropodes, est transmis maternellement. Afin de se répandre dans la population d?hôtes, Wolbachia manipule la physiolo- gie reproductive de l?organisme qui l?abrite. Le Chapitre 3 décrit une étude consacrée aux variations du sexe-ratio produit dans des colonies de fourmis? en fonction de la présence ou de l?absence de Wolbachia. Le Chapitre 4 étudie l?association entre les lignées du symbionte et les haplotypes mitochondriaux des hôtes. L?étude de l?association entre les deux éléments maternellement transmis est d?un intérêt particulier chez Formica exsecta car le taux de migra- tion femelle très faible engendre une forte strucure génétique mitochondriale.