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.

Functions of the peroxisome proliferator-activated receptor (PPAR) alpha and beta in skin homeostasis, epithelial repair, and morphogenesis.

The three peroxisome proliferator-activated receptors (PPAR alpha, PPAR beta, and PPAR gamma) are ligand-activated transcription factors belonging to the nuclear hormone receptor superfamily. They are regarded as being sensors of physiological levels of fatty acids and fatty acid derivatives. In the adult mouse skin, they are found in hair follicle keratinocytes but not in interfollicular epidermis keratinocytes. Skin injury stimulates the expression of PPAR alpha and PPAR beta at the site of the wound. Here, we review the spatiotemporal program that triggers PPAR beta expression immediately after an injury, and then gradually represses it during epithelial repair. The opposing effects of the tumor necrosis factor-alpha and transforming growth factor-beta-1 signalling pathways on the activity of the PPAR beta promoter are the key elements of this regulation. We then compare the involvement of PPAR beta in the skin in response to an injury and during hair morphogenesis, and underscore the similarity of its action on cell survival in both situations.

The histone-interacting domain of nuclear factor I activates simian virus 40 DNA replication in vivo.

Efficient initiation of SV40 DNA replication requires transcription factors that bind auxiliary sequences flanking the minimally required origin. To evaluate the possibility that transcription factors may activate SV40 replication by acting on the chromatin structure of the origin, we used an in vivo replication system in which we targeted GAL4 fusion proteins to the minimally required origin. We found that the proline-rich transcriptional activation domain of nuclear factor I (NF-I), which has been previously shown to interact with histone H3, specifically activates replication. Evaluation of a series of deletion and point mutants of NF-I indicates that the H3-binding domain and the replication activity coincide perfectly. Assays with other transcription factors, such as Sp1, confirmed the correlation between the interaction with H3 and the activation of replication. These findings imply that transcription factors such as NF-I can activate SV40 replication via direct interaction with chromatin components, thereby contributing to the relief of nucleosomal repression at the SV40 origin.

Coregulator Control of Androgen Receptor Action by a Novel Nuclear Receptor-Binding Motif

The androgen receptor (AR) is a ligand-activated transcription factor that is essential for prostate cancer development. It is activated by androgens through its ligand-binding domain (LBD), which consists predominantly of 11 α-helices. Upon ligand binding, the last helix is reorganized to an agonist conformation termed activator function-2 (AF-2) for coactivator binding. Several coactivators bind to the AF-2 pocket through conserved LXXLL or FXXLF sequences to enhance the activity of the receptor. Recently, a small compound-binding surface adjacent to AF-2 has been identified as an allosteric modulator of the AF-2 activity and is termed binding function-3 (BF-3). However, the role of BF-3 in vivo is currently unknown, and little is understood about what proteins can bind to it. Here we demonstrate that a duplicated GARRPR motif at the N terminus of the cochaperone Bag-1L functions through the BF-3 pocket. These findings are supported by the fact that a selective BF-3 inhibitor or mutations within the BF-3 pocket abolish the interaction between the GARRPR motif(s) and the BF-3. Conversely, amino acid exchanges in the two GARRPR motifs of Bag-1L can impair the interaction between Bag-1L and AR without altering the ability of Bag-1L to bind to chromatin. Furthermore, the mutant Bag-1L increases androgen-dependent activation of a subset of AR targets in a genome-wide transcriptome analysis, demonstrating a repressive function of the GARRPR/BF-3 interaction. We have therefore identified GARRPR as a novel BF-3 regulatory sequence important for fine-tuning the activity of the AR.

Integration of Notch 1 and calcineurin/NFAT signaling pathways in keratinocyte growth and differentiation control.

The Notch and Calcineurin/NFAT pathways have both been implicated in control of keratinocyte differentiation. Induction of the p21(WAF1/Cip1) gene by Notch 1 activation in differentiating keratinocytes is associated with direct targeting of the RBP-Jkappa protein to the p21 promoter. We show here that Notch 1 activation functions also through a second Calcineurin-dependent mechanism acting on the p21 TATA box-proximal region. Increased Calcineurin/NFAT activity by Notch signaling involves downregulation of Calcipressin, an endogenous Calcineurin inhibitor, through a HES-1-dependent mechanism. Besides control of the p21 gene, Calcineurin contributes significantly to the transcriptional response of keratinocytes to Notch 1 activation, both in vitro and in vivo. In fact, deletion of the Calcineurin B1 gene in the skin results in a cyclic alopecia phenotype, associated with altered expression of Notch-responsive genes involved in hair follicle structure and/or adhesion to the surrounding mesenchyme. Thus, an important interconnection exists between Notch 1 and Calcineurin-NFAT pathways in keratinocyte growth/differentiation control.

Cytoskeletal and DNA structure abnormalities result from bypass of requirement for the cdc10 start gene in the fission yeast Schizosaccharomyces pombe.

The cdc10 gene of the fission yeast S. pombe is required for traverse of the start control in late G1 and commitment to the mitotic cell cycle. To increase our understanding of the events which occur at start, a pseudoreversion analysis was undertaken to identify genes whose products may interact with cdc10 or bypass the requirement for it. A single gene, sct1+ (suppressor of cdc ten), has been identified, mutation of which suppresses all conditional alleles and a null allele of cdc10. Bypass of the requirement for cdc10+ function by sct1-1 mutations leads to pleiotropic defects, including microtubule, microfilament and nuclear structural abnormalities. Our data suggest that sct1 encodes a protein that is dependent upon cdc10+ either for its normal function or expression, or is a component of a checkpoint that monitors execution of p85cdc10 function.

Turvallisuustoimintojen käyttö ydinvoimalaitoksilla

Tässä diplomityössä tarkastellaan turvallisuustoiminto-käsitteen määrittelyä ja käyttöä osana ydinvoimalaitosten turvallisuuden varmistamista. Työssä kuvataan paine- ja kiehutusvesilaitosten toiminnan yleispiirteet, sekä Teollisuuden Voima Oy:n (TVO) laitosten Olkiluoto 1 & 2 sekä Olkiluoto 3 tarkempi turvallisuustoiminta turvallisuusjärjestelmien ja -automaation osalta. Työssä esitellään eräs tapa määrittää turvallisuustoimintoja. Malli perustuu hierarkkiseen rakenteeseen, jossa ylimpänä ovat laitostason turvallisuustoiminnot ja alimpana turvallisuustoimintoihin osallistuvien laitteiden ja niiden osien toiminnot. Turvallisuustoimintoja on mahdollista käyttää ydinvoimalaitoksen turvallisuusluokituksen tekemiseen ja perustelemiseen. Turvallisuustoiminto kertoo suoraan luokiteltavan kohteen turvallisuusmerkityksen. Kohteen turvallisuusmerkityksen selvittäminen, eli liittäminen turvallisuustoimintoon, voi olla vaikeaa. Turvallisuustoimintoja on myös mahdollista käyttää laitoksen turvallisuusautomaation riittävän varmistamisen (mm. redundanttisuus, diversiteetti ja erotus) osoittamiseen erityisesti ohjelmoitavan automaation yhteydessä. Turvallisuustoimintoja voidaan hyödyntää laitoksen hätätilanneohjeiden kehittämisessä ja myös laitoksenmuun turvallisuusdokumentaation selkeyttämisessä. Työn tuloksena kehitettiin käyville laitoksille (OL1 & 2) aikaisempaa kattavammat turvallisuustoiminnot. Lisäksi tarkasteltiin rakenteilla olevalle laitokselle (OL3) määriteltyjä turvallisuustoimintoja.

Different concentrations of Mg++ ions affect nuclear matrix protein distribution during thermal stabilization of isolated nuclei.

The nuclear matrix, a proteinaceous network believed to be a scaffolding structure determining higher-order organization of chromatin, is usually prepared from intact nuclei by a series of extraction steps. In most cell types investigated the nuclear matrix does not spontaneously resist these treatments but must be stabilized before the application of extracting agents. Incubation of isolated nuclei at 37C or 42C in buffers containing Mg++ has been widely employed as stabilizing agent. We have previously demonstrated that heat treatment induces changes in the distribution of three nuclear scaffold proteins in nuclei prepared in the absence of Mg++ ions. We studied whether different concentrations of Mg++ (2.0-5 mM) affect the spatial distribution of nuclear matrix proteins in nuclei isolated from K562 erythroleukemia cells and stabilized by heat at either 37C or 42C. Five proteins were studied, two of which were RNA metabolism-related proteins (a 105-kD component of splicing complexes and an RNP component), one a 126-kD constituent of a class of nuclear bodies, and two were components of the inner matrix network. The localization of proteins was determined by immunofluorescent staining and confocal scanning laser microscope. Mg++ induced significant changes of antigen distribution even at the lowest concentration employed, and these modifications were enhanced in parallel with increase in the concentration of the divalent cation. The different sensitivity to heat stabilization and Mg++ of these nuclear proteins might reflect a different degree of association with the nuclear scaffold and can be closely related to their functional or structural role.

Regulation of recombination at yeast nuclear pores controls repair and triplet repeat stability.

Secondary structure-forming DNA sequences such as CAG repeats interfere with replication and repair, provoking fork stalling, chromosome fragility, and recombination. In budding yeast, we found that expanded CAG repeats are more likely than unexpanded repeats to localize to the nuclear periphery. This positioning is transient, occurs in late S phase, requires replication, and is associated with decreased subnuclear mobility of the locus. In contrast to persistent double-stranded breaks, expanded CAG repeats at the nuclear envelope associate with pores but not with the inner nuclear membrane protein Mps3. Relocation requires Nup84 and the Slx5/8 SUMO-dependent ubiquitin ligase but not Rad51, Mec1, or Tel1. Importantly, the presence of the Nup84 pore subcomplex and Slx5/8 suppresses CAG repeat fragility and instability. Repeat instability in nup84, slx5, or slx8 mutant cells arises through aberrant homologous recombination and is distinct from instability arising from the loss of ligase 4-dependent end-joining. Genetic and physical analysis of Rad52 sumoylation and binding at the CAG tract suggests that Slx5/8 targets sumoylated Rad52 for degradation at the pore to facilitate recovery from acute replication stress by promoting replication fork restart. We thereby confirmed that the relocation of damage to nuclear pores plays an important role in a naturally occurring repair process.

Nuclear receptor peroxisome proliferator activated receptor (PPAR) beta/delta in skin wound healing and cancer

We review the functions of peroxisome proliferator activated receptor (PPAR) beta/delta in skin wound healing and cancer. In particular, we highlight the roles of PPAR beta/delta in inhibiting keratinocyte apoptosis at wound edges via activation of the PI3K/PKB alpha/Akt1 pathway and its role during re-epithelialization in regulating keratinocyte adhesion and migration. In fibroblasts, PPAR beta/delta controls IL-1 signalling and thereby contributes to the homeostatic control of keratinocyte proliferation. We discuss its therapeutic potential for treating diabetic wounds and inflammatory skin diseases such as psoriasis and acne vulgaris. PPAR beta/delta is classified as a tumour growth modifier; it is activated by chronic low-grade inflammation, which promotes the production of lipids that, in turn, enhance PPAR beta/delta transcription activity. Our earlier,work unveiled a cascade of events triggered by PPAR beta/delta that involve the oncogene Src, which promotes ultraviolet-induced skin cancer in mice via enhanced EGFR/Erk1/2 signalling and the expression of epithelial-to-mesenchymal transition (EMT) markers. Interestingly, PPAR beta/delta expression is correlated with the expression of SRC and EMT markers in human skin squamous cell carcinoma. Furthermore, there is a positive interaction between PPAR beta/delta, SRC, and TGF beta 1 at the transcriptional level in various human epithelial cancers. Taken together, these observations suggest the need for evaluating PPAR beta/delta modulators that attenuate or increase its activity, depending on the therapeutic target.

Phylogeography and population genetics of north European Atlantic salmon (Salmo salar L.)

Although abundant in the number of individuals, the Atlantic salmon may be considered as a threatened species in many areas of its native distribution range. Human activities such as building of power plant dams, offshore overfishing, pollution, clearing of riverbeds for timber floating and badly designed stocking regimes have diminished the distribution of Atlantic salmon. As a result of this, many of the historical populations both in Europe and northern America have gone extinct or are severely depressed. In fact, only 1% of Atlantic salmon existing today are of natural origin, the rest being farmed salmon. All of this has lead to a vast amount of research and many restoration programmes aiming to bring Atlantic salmon back to rivers from where it has vanished. However, many of the restoration programmes conducted thus far have been unsuccessful due to inadequate scientific research or lack of its implementation, highlighting the fact that more research is needed to fully understand the biology of this complex species. The White and Barents Seas in northwest Russia are among the last regions in Europe where Atlantic salmon populations are still stable, thus forming an important source of biodiversity for the entire European region. Salmon stocks from this area are also of immense economic and social importance for the local people in the form of fishing tourism. The main aim of this thesis was to elucidate the post-glacial history and population genetic structure of north European and particularly northwest Russian Atlantic salmon, both of which are aspects of great importance for the management and conservation of the species. Throughout the whole thesis, these populations were studied by utilizing microsatellites as the main molecular tool. One of the most important discoveries of the thesis was the division of Atlantic salmon from the White and Barents Seas into four separate clusters, which has not been observed in previous studies employing nuclear markers although is supported by mtDNA studies. Populations from the western Barents Sea clustered together with the northeast Atlantic populations into a clearly distinguishable group while populations from the White Sea and eastern Barents Sea were separated into three additional groups. This has important conservation implications as this thesis clearly indicates that conservation of populations from all of the observed clusters is warranted in order to conserve as much of the genetic diversity as possible in this area. The thesis also demonstrates how differences in population life histories within a species, migratory behaviour in this case, and in their phylogeographic origin affect the genetic characteristics of populations, namely diversity and divergence levels. The anadromous populations from the Atlantic Ocean, White Sea and Barents Sea possessed higher levels of genetic diversity than the anadromous populations form the Baltic Sea basin. Among the non-anadromous populations the result was the opposite: the Baltic freshwater populations were more variable. This emphasises the importance of taking the life history of a population into consideration when developing conservation strategies: due to the limited possibilities for new genetic diversity to be generated via gene flow, it is expected that freshwater Atlantic salmon populations would be more vulnerable to extinction following a population crash and thus deserve a high conservation status. In the last chapter of this thesis immune relevant marker loci were developed and screened for signatures of natural selection along with loci linked to genes with other functions or no function at all. Also, a novel landscape genomics method, which combines environmental information with molecular data, was employed to investigate whether immune relevant markers displayed significant correlations to various environmental variables more frequently than other loci. Indications of stronger selection pressure among immune-relevant loci compared to non-immune relevant EST-linked loci was found but further studies are needed to evaluate whether it is a common phenomenon in Atlantic salmon.

Electron scattering in isotonic chains as a probe of the proton shell structure of unstable nuclei

Electron scattering on unstable nuclei is planned in future facilities of the GSI and RIKEN upgrades. Motivated by this fact, we study theoretical predictions for elastic electron scattering in the N=82, N=50, and N=14 isotonic chains from very proton-deficient to very proton-rich isotones. We compute the scattering observables by performing Dirac partial-wave calculations. The charge density of the nucleus is obtained with a covariant nuclear mean-field model that accounts for the low-energy electromagnetic structure of the nucleon. For the discussion of the dependence of scattering observables at low-momentum transfer on the gross properties of the charge density, we fit Helm model distributions to the self-consistent mean-field densities. We find that the changes shown by the electric charge form factor along each isotonic chain are strongly correlated with the underlying proton shell structure of the isotones. We conclude that elastic electron scattering experiments on isotones can provide valuable information about the filling order and occupation of the single-particle levels of protons.

Hepatic circadian clock oscillators and nuclear receptors integrate microbiome-derived signals.

The liver is a key organ of metabolic homeostasis with functions that oscillate in response to food intake. Although liver and gut microbiome crosstalk has been reported, microbiome-mediated effects on peripheral circadian clocks and their output genes are less well known. Here, we report that germ-free (GF) mice display altered daily oscillation of clock gene expression with a concomitant change in the expression of clock output regulators. Mice exposed to microbes typically exhibit characterized activities of nuclear receptors, some of which (PPARα, LXRβ) regulate specific liver gene expression networks, but these activities are profoundly changed in GF mice. These alterations in microbiome-sensitive gene expression patterns are associated with daily alterations in lipid, glucose, and xenobiotic metabolism, protein turnover, and redox balance, as revealed by hepatic metabolome analyses. Moreover, at the systemic level, daily changes in the abundance of biomarkers such as HDL cholesterol, free fatty acids, FGF21, bilirubin, and lactate depend on the microbiome. Altogether, our results indicate that the microbiome is required for integration of liver clock oscillations that tune output activators and their effectors, thereby regulating metabolic gene expression for optimal liver function.

Hepatic circadian clock oscillators and nuclear receptors integrate microbiome-derived signals.

The liver is a key organ of metabolic homeostasis with functions that oscillate in response to food intake. Although liver and gut microbiome crosstalk has been reported, microbiome-mediated effects on peripheral circadian clocks and their output genes are less well known. Here, we report that germ-free (GF) mice display altered daily oscillation of clock gene expression with a concomitant change in the expression of clock output regulators. Mice exposed to microbes typically exhibit characterized activities of nuclear receptors, some of which (PPARα, LXRβ) regulate specific liver gene expression networks, but these activities are profoundly changed in GF mice. These alterations in microbiome-sensitive gene expression patterns are associated with daily alterations in lipid, glucose, and xenobiotic metabolism, protein turnover, and redox balance, as revealed by hepatic metabolome analyses. Moreover, at the systemic level, daily changes in the abundance of biomarkers such as HDL cholesterol, free fatty acids, FGF21, bilirubin, and lactate depend on the microbiome. Altogether, our results indicate that the microbiome is required for integration of liver clock oscillations that tune output activators and their effectors, thereby regulating metabolic gene expression for optimal liver function.

Enzymes of yeast polyphosphate metabolism: structure, enzymology and biological roles.

Inorganic polyphosphate (polyP) is found in all living organisms. The known polyP functions in eukaryotes range from osmoregulation and virulence in parasitic protozoa to modulating blood coagulation, inflammation, bone mineralization and cellular signalling in mammals. However mechanisms of regulation and even the identity of involved proteins in many cases remain obscure. Most of the insights obtained so far stem from studies in the yeast Saccharomyces cerevisiae. Here, we provide a short overview of the properties and functions of known yeast polyP metabolism enzymes and discuss future directions for polyP research.