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.

A hybrid algorithm combining metaheuristic with Monte Carlo simulation for solving the Stochastic Flow Shop problem

In this paper, a hybrid simulation-based algorithm is proposed for the StochasticFlow Shop Problem. The main idea of the methodology is to transform the stochastic problem into a deterministic problem and then apply simulation to the latter. In order to achieve this goal, we rely on Monte Carlo Simulation and an adapted version of a deterministic heuristic. This approach aims to provide flexibility and simplicity due to the fact that it is not constrained by any previous assumption and relies in well-tested heuristics.

Selection strategies of segregant soybean populations for resistance to Asian rust

The objective of this work was to identify the best selection strategies for the more promising parental combinations to obtain lines with good resistance to soybean Asian rust (Phakopsora pachyrhizi). Two experiments were carried out in the field during the 2006/2007 and 2007/2008 growing seasons, to determine the percentage of infected leaf area of individual plants of five parents and their segregant F2 and F3 populations. The data obtained indicates that additive genetic variance predominates in the control of soybean resistance to Asian rust, and that the year and time of assessment do not significantly influence the estimates of the genetic parameters obtained. The narrow-sense heritability (h²r) ranged from 23.12 to 55.83%, and indicates the possibility of successful selection of resistant individuals in the early generations of the breeding program. All the procedures used to select the most promising populations to generate superior inbred lines for resistance to P. pachyrhizi presented similar results and identified the BR01-18437 x BRS 232 population as the best for inbred line selection.

Leaf plasticity in successive selection cycles of 'Saracura' maize in response to periodic soil flooding

The objective of this work was to assess the effect of successive selection cycles on leaf plasticity of 'Saracura' maize BRS-4154 under periodical flooding in field conditions. Soil flooding started at the six-leaf stage with the application of a 20-cm depth water layer three times a week. At flowering, samples of leaves were collected and fixed. Paradermic and transverse sections were observed under photonic microscope. Several changes were observed throughout the selection cycles, such as modifications in the number and size of the stomata, higher amount of vascular bundles and the resulting decrease of the distance between them, smaller diameter of the metaxylem, decrease of cuticle and epidermis thickness, decrease of number and size of bulliform cells, increase of phloem thickness, smaller sclerenchyma area. Therefore, the successive selection cycles of 'Saracura' maize resulted in changes in the leaf anatomy, which might be favorable to the plant's tolerance to the intermittent flooding of the soil.

Habitat fragmentation, ecology and sexual selection in forest bird species in Monteverde, Costa Rica

Abstract Forest fragmentation is often associated with local extinction and changes in species abundance patterns. The main topic of this thesis is the effect of forest fragmentation on birds in Monteverde, Costa Rica. This thesis also studies aspects of sexual selection and ecology of Long-tailed Manakins, Chiroxiphia Linearis. Chapter 1 investigates bird species assemblages in two degrees of forest fragmentation. It is shown that the distribution, abundance and diversity of forest bird species are strongly influenced by the amount of forest in the landscape matrix. Presence of cattle within the forest influences the presence of some bird species. The prevalence and intensity of ticks and blood parasites on birds in relation to fragmentation is described in Chapter 2. Overall tick prevalence is 3%. Understory birds are significantly more infested with ticks than species at intermediate heights. Tick prevalence on birds does not differ significantly between two degrees of forest fragmentation and individual tick loads tend to be higher in High- than in Low-fragmentation sites. Infestations by the blood parasites Haemoproteus sp. was low except in white-eared ground sparrow, Melozone leucotis, that is 28% and is significantly higher in High- than in Low-fragmentation sites. In chapter 3 results on the ecology and habitat movements of the Bare-necked Umbrellabird, Cephalopterus glabricollis, are presented. The abundance of umbrellabirds at high elevations during the breeding season coincides with the highest peak of fruit abundance. Birds leave the protected area during the non-breeding season moving to unprotected forest fragments. In chapter 4 ontogenetic changes in feather morphology through sexual maturity in Long-tailed Manakins are described. In adult males, rectrices length is positively correlated to testis volume. Changes in male morphology during ontogeny in the long-tailed manakin appear to be associated with their specific-display behaviours. Significant interpopulation differences in the morphology of Long-tailed Manakins are shown in chapter 5. These differences are more accentuated in morphological traits related to flight displays. A field experiment demonstrates that long rectrices impose flying costs for males and females. A reduction in flying ability was found to be strongest in males from a population presenting the highest degree of sexual dimorphism. Résumé La fragmentation des forêts est souvent associée avec des modifications dans l'abondance des espèces et des extinctions locales. Le thème principale de cette thèse est l'étude de l'effet de la fragmentation des forêts sur les oiseaux de Monteverde, Costa Rica. Elle décrit par ailleurs certains aspects de la sélection sexuelle et l'écologie du manakin à longue queue, Chiroxiphia linearis. Dans le Chapitre 1 je montre que la distribution, l'abondance et la diversité des assemblages d'oiseaux vivant dans la forêt sont fortement influencées pas le degré de fragmentation de celle ci. Par ailleurs, la présence ou l'absence de bétail dans les forêts influence la présence de certaines espèces d'oiseaux. Dans le chapitre 2 j'ai étudié la prévalence et l'intensité d'infestation par des tiques ainsi que la présence de parasites sanguins chez les oiseaux en relation avec la fragmentation des forêts. La prévalence globale de tiques est de 3 %, les oiseaux vivant au niveau du sol étaient plus souvent infectés par des tiques que les espèces se déplaçant à un niveau plus élevé. La prévalence de tiques sur les oiseaux n'était pas significativement différente entre les paysages avec différentes fragmentations. Les parasites sanguins du genre Haemoproteus sp. étaient présents à très basse fréquence à l'exception chez Melozone leucotis ou la prevalence était de 28% et significativement plus élevée chez les oiseaux vivant dans les forêts à forte fragmentation. Dans le Chapitre 3 je présente des résultats sur l'écologie et les mouvements entre habitats chez le "Bare-necked umbrellabird", Cephalopterus glabricollis. Cette espèce endémique du Costa Rica niche à haute altitude durant la période d'abondance des fruits et réalise une migration altitudinale vers des zones basses durant la saison de non reproduction. Dans le chapitre 4 je présente les changements ontogénétiques dans la morphologie du plumage des manakins à longue queue. Chez les mâles, les changements de morphologie semblent être associés avec leurs comportements de parade spécifiques. Dans le chapitre 5 je présente des différences morphologiques significative entre deux populations chez le manakin à longue queue et je montre que la capacité de vols chez les mâles est plus fortement influencée dans la population avec le degré de dimorphisme sexuel le plus prononcé.

Parent selection for cocoa resistance to witches'-broom

The objective of this work was to identify genotypes with high general combining ability for resistance to witches'-broom (Moniliophthora perniciosa) in populations formed from a first cycle of recurrent selection. Highly productive and resistant clones from different origins were interbred using the North Carolina II design. The clones SCA 6, CSUL 7, RB 39, CEPEC 89, OC 67, BE 4, EEG 29 and ICS 98 were used as paternal parents, while the maternal ones were NA 33, CCN 10, IMC 67, P 4B, CCN 51, CEPEC 86, SGU 54 and ICS 9. Twenty days after germination, 56 seedlings of each cross (four replicates of 14 seedlings) received the inoculation of a 1-mL suspension with 7.5x10(4 ) basidiospores mL-1. Symptoms were evaluated 60 days after inoculation. Significant differences were observed among paternal and among maternal parents, for resistance to witches'-broom assessed according to the proportion of progeny seedlings with the disease symptoms. Differences were also observed between groups of mothers or fathers previously defined as resistant, and groups previously defined as susceptible. It is possible to obtain a combination of genes that can increase the level of resistance to witches'-broom directly from the first cycle of recurrent selection.

Transcriptome of a mouse kidney cortical collecting duct cell line: effects of aldosterone and vasopressin.

Aldosterone and vasopressin are responsible for the final adjustment of sodium and water reabsorption in the kidney. In principal cells of the kidney cortical collecting duct (CCD), the integral response to aldosterone and the long-term functional effects of vasopressin depend on transcription. In this study, we analyzed the transcriptome of a highly differentiated mouse clonal CCD principal cell line (mpkCCD(cl4)) and the changes in the transcriptome induced by aldosterone and vasopressin. Serial analysis of gene expression (SAGE) was performed on untreated cells and on cells treated with either aldosterone or vasopressin for 4 h. The transcriptomes in these three experimental conditions were determined by sequencing 169,721 transcript tags from the corresponding SAGE libraries. Limiting the analysis to tags that occurred twice or more in the data set, 14,654 different transcripts were identified, 3,642 of which do not match known mouse sequences. Statistical comparison (at P < 0.05 level) of the three SAGE libraries revealed 34 AITs (aldosterone-induced transcripts), 29 ARTs (aldosterone-repressed transcripts), 48 VITs (vasopressin-induced transcripts) and 11 VRTs (vasopressin-repressed transcripts). A selection of the differentially-expressed, hormone-specific transcripts (5 VITs, 2 AITs and 1 ART) has been validated in the mpkCCD(cl4) cell line either by Northern blot hybridization or reverse transcription-PCR. The hepatocyte nuclear transcription factor HNF-3-alpha (VIT39), the receptor activity modifying protein RAMP3 (VIT48), and the glucocorticoid-induced leucine zipper protein (GILZ) (AIT28) are candidate proteins playing a role in physiological responses of this cell line to vasopressin and aldosterone.

Selection of shrimp breeders free of white spot syndrome and infectious hypodermal and hematopoietic necrosis

The objective of this work was to select surviving breeders of Litopenaeus vannamei from white spot syndrome virus (WSSV) outbreak, adapted to local climatic conditions and negatively diagnosed for WSSV and infectious hypodermal and hematopoietic necrosis virus (IHHNV), and to evaluate if this strategy is a viable alternative for production in Santa Catarina, Brazil. A total of 800 males and 800 females were phenotypically selected in a farm pond. Nested-PCR analyses of 487 sexually mature females and 231 sexually mature males showed that 63% of the females and 55% of the males were infected with IHHNV. Animals free of IHHNV were tested for WSSV, and those considered double negative were used for breeding. The post-larvae produced were stocked in nine nursery tanks for analysis. From the 45 samples, with 50 post-larvae each, only two were positive for IHHNV and none for WSSV. Batches of larvae diagnosed free of virus by nested-PCR were sent to six farms. A comparative analysis was carried out in growth ponds, between local post-larvae and post-larvae from Northeast Brazil. Crabs (Chasmagnathus granulata), blue crabs (Callinectes sapidus), and sea hares (Aplysia brasiliana), which are possible vectors of these viruses, were also evaluated. The mean survival was 55% for local post-larvae against 23.4% for post-larvae from the Northeast. Sea hares showed prevalence of 50% and crabs of 67% of WSSV.

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.

Detection and mapping of a lethal locus in a eucalyptus hybrid population

The objective of this work was to verify the existence of a lethal locus in a eucalyptus hybrid population, and to quantify the segregation distortion in the linkage group 3 of the Eucalyptus genome. A E. grandis x E. urophylla hybrid population, which segregates for rust resistance, was genotyped with 19 microsatellite markers belonging to linkage group 3 of the Eucalyptus genome. To quantify the segregation distortion, maximum likelihood (ML) models, specific to outbreeding populations, were used. These models consider the observed marker genotypes and the lethal locus viability as parameters. The ML solutions were obtained using the expectation‑maximization algorithm. A lethal locus in the linkage group 3 was verified and mapped, with high confidence, between the microssatellites EMBRA 189 e EMBRA 122. This lethal locus causes an intense gametic selection from the male side. Its map position is 25 cM from the locus which controls the rust resistance in this population.

Altered thymic selection by overexpressing cellular FLICE inhibitory protein in T cells causes lupus-like syndrome in a BALB/c but not C57BL/6 strain.

The cellular FLICE inhibitory protein (c-FLIP) is an endogenous inhibitor of the caspase-8 proapoptotic signaling pathway downstream of death receptors. Recent evidence indicates that the long form of c-FLIP (c-FLIP(L)) is required for proliferation and effector T-cell development. However, the role of c-FLIP(L) in triggering autoimmunity has not been carefully analyzed. We now report that c-FLIP(L) transgenic (Tg) mice develop splenomegaly, lymphadenopathy, multiorgan infiltration, high titers of auto-antibodies, and proliferative glomerulonephritis with immune complex deposition in a strain-dependent manner. The development of autoimmunity requires CD4(+) T cells and may result from impaired thymic selection. At the molecular level, c-FLIP(L) overexpression inhibits the zeta chain-associated protein tyrosine kinase of 70 kDa (ZAP-70) activation, thus impairing the signaling pathway derived from ZAP-70 required for thymic selection. Therefore, we have identified c-FLIP(L) as a susceptibility factor under the influence of epistatic modifiers for the development of autoimmunity.

Optimization strategies for fast detection of positive selection on phylogenetic trees.

MOTIVATION: The detection of positive selection is widely used to study gene and genome evolution, but its application remains limited by the high computational cost of existing implementations. We present a series of computational optimizations for more efficient estimation of the likelihood function on large-scale phylogenetic problems. We illustrate our approach using the branch-site model of codon evolution. RESULTS: We introduce novel optimization techniques that substantially outperform both CodeML from the PAML package and our previously optimized sequential version SlimCodeML. These techniques can also be applied to other likelihood-based phylogeny software. Our implementation scales well for large numbers of codons and/or species. It can therefore analyse substantially larger datasets than CodeML. We evaluated FastCodeML on different platforms and measured average sequential speedups of FastCodeML (single-threaded) versus CodeML of up to 5.8, average speedups of FastCodeML (multi-threaded) versus CodeML on a single node (shared memory) of up to 36.9 for 12 CPU cores, and average speedups of the distributed FastCodeML versus CodeML of up to 170.9 on eight nodes (96 CPU cores in total).Availability and implementation: ftp://ftp.vital-it.ch/tools/FastCodeML/. CONTACT: selectome@unil.ch or nicolas.salamin@unil.ch.

The contribution of a pollinating seed predator to selection on Silene latifolia females.

Interactions, antagonistic or mutualistic, can exert selection on plant traits. We explored the role of Hadena bicruris, a pollinating seed predator, as a selective agent on its host, the dioecious plant Silene latifolia. We exposed females from artificial-selection lines (many, small flowers (SF) vs. few, large flowers (LF)) to this moth. Infestation did not differ significantly between lines, but the odds of attacked fruits aborting were higher in SF females. We partitioned selection between that caused by moth attack and that resulting from all other factors. In both lines, selection via moth attack for fewer, smaller flowers contrasted with selection via other factors for more flowers. In LF females, selection via the two components was strongest and selection via moth attack also favoured increased fruit abortion. This suggests that the moths act as more of a selective force on flower size and number via their predating than their pollinating role.