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.

Exposure to bioaerosols in poultry houses at different stages of fattening: use of real-time PCR for airborne bacterial quantification

Previous studies have demonstrated that poultry-house workers are exposed to very high levels of organic dust and consequently have an increased prevalence of adverse respiratory symptoms. However, the influence of the age of broilers, on bioaerosol concentrations has not been investigated. To evaluate the evolution of bioaerosol concentration during the fattening period, bioaerosol parameters (inhalable dust, endotoxin and bacteria) were measured in 12 poultry confinement buildings in Switzerland, at 3 different stages of the birds' growth; Samples of air taken from within the breathing zones of individual poultry-house employees as they caught the chickens ready to be transported for slaughter, were also analysed. Quantitative PCR (Q-PCR) was used to assess the quantity of total airborne bacteria and total airborne Staphylococcus species. Bioaerosol levels increased significantly during the fattening period of the chickens. During the task of catching mature birds, the mean inhalable dust concentration for a worker was 31 ± 4.7 mg/m3, and endotoxin concentration was 11'080 ± 3436 UE/m3 air, more than ten-fold higher than the Swiss occupational recommended value (1000 UE/m3). The mean exposure level of bird catchers to total bacteria and Staphylococcus species measured by Q-PCR is also very high, respectively reaching values of 72 (± 11) x107 cells/m3 air and 70 (± 16) x106/m3 air. It was concluded that in the absence of wearing protective breathing apparatus, chicken catchers in Switzerland risk exposure beyond recommended limits for all measured bioaerosol parameters. Moreover, the use of Q-PCR to estimate total and specific numbers of airborne bacteria is a promising tool for evaluating any modifications intended to improve the safety of current working practices.

Veille bibliographique sur la surveillance biologique de l'exposition aux produits chimiques au poste de travail : l'expérience d'un réseau francophone multidisciplinaire

Une veille bibliographique est organisée depuis 2005 sur la surveillance biologique aux produits chimiques en milieu de travail (SBEPC MT). Elle a été mise en place par le réseau francophone multidisciplinaire, composé de l'INRS (France), l'IRSST (Québec) et l'UCL (Belgique). Cet article dresse le bilan de l'information récoltée et analysée, de 2009 à 2012, au travers de 435 articles sélectionnés. Plusieurs thèmes d'intérêt ou d'actualités font l'objet d'une analyse plus approfondie, dont notamment les pesticides, les hydrocarbures aromatiques, le benzène, le mangannèse, la variabilité biologique, les dosages cutanés et frottis de surface, les dosages dans l'air expiré ou encore la spectrométrie de masse.

Determination of 241Pu in nuclear waste slurries: a comparative study using LSC and ICP-MS.

(241)Pu was determined in slurry samples from a nuclear reactor decommissioning project at the Paul Scherrer Institute (Switzerland). To validate the results, the (241)Pu activities of five samples were determined by LSC (TriCarb and Quantulus) and ICP-MS, with each instrument at a different laboratory. In lack of certified reference materials for (241)Pu, the methods were further validated using the (241)Pu information values of two reference sediments (IAEA-300 and IAEA-384). Excellent agreement with the results was found between LSC and ICP-MS in the nuclear waste slurries and the reference sediments.

A spatially explicit life cycle inventory of the global textile chain

Life cycle analyses (LCA) approaches require adaptation to reflect the increasing delocalization of production to emerging countries. This work addresses this challenge by establishing a country-level, spatially explicit life cycle inventory (LCI). This study comprises three separate dimensions. The first dimension is spatial: processes and emissions are allocated to the country in which they take place and modeled to take into account local factors. Emerging economies China and India are the location of production, the consumption occurs in Germany, an Organisation for Economic Cooperation and Development country. The second dimension is the product level: we consider two distinct textile garments, a cotton T-shirt and a polyester jacket, in order to highlight potential differences in the production and use phases. The third dimension is the inventory composition: we track CO2, SO2, NO (x), and particulates, four major atmospheric pollutants, as well as energy use. This third dimension enriches the analysis of the spatial differentiation (first dimension) and distinct products (second dimension). We describe the textile production and use processes and define a functional unit for a garment. We then model important processes using a hierarchy of preferential data sources. We place special emphasis on the modeling of the principal local energy processes: electricity and transport in emerging countries. The spatially explicit inventory is disaggregated by country of location of the emissions and analyzed according to the dimensions of the study: location, product, and pollutant. The inventory shows striking differences between the two products considered as well as between the different pollutants considered. For the T-shirt, over 70% of the energy use and CO2 emissions occur in the consuming country, whereas for the jacket, more than 70% occur in the producing country. This reversal of proportions is due to differences in the use phase of the garments. For SO2, in contrast, over two thirds of the emissions occur in the country of production for both T-shirt and jacket. The difference in emission patterns between CO2 and SO2 is due to local electricity processes, justifying our emphasis on local energy infrastructure. The complexity of considering differences in location, product, and pollutant is rewarded by a much richer understanding of a global production-consumption chain. The inclusion of two different products in the LCI highlights the importance of the definition of a product's functional unit in the analysis and implications of results. Several use-phase scenarios demonstrate the importance of consumer behavior over equipment efficiency. The spatial emission patterns of the different pollutants allow us to understand the role of various energy infrastructure elements. The emission patterns furthermore inform the debate on the Environmental Kuznets Curve, which applies only to pollutants which can be easily filtered and does not take into account the effects of production displacement. We also discuss the appropriateness and limitations of applying the LCA methodology in a global context, especially in developing countries. Our spatial LCI method yields important insights in the quantity and pattern of emissions due to different product life cycle stages, dependent on the local technology, emphasizing the importance of consumer behavior. From a life cycle perspective, consumer education promoting air-drying and cool washing is more important than efficient appliances. Spatial LCI with country-specific data is a promising method, necessary for the challenges of globalized production-consumption chains. We recommend inventory reporting of final energy forms, such as electricity, and modular LCA databases, which would allow the easy modification of underlying energy infrastructure.

Comparison of naphthalene bioavailability determined by whole-cell biosensing and availability determined by extraction with Tenax

A rapid biological method for the determination of the bioavailability of naphthalene was developed and its value as an alternative to extraction-based chemical approaches demonstrated. Genetically engineered whole-cell biosensors are used to determine bioavailable naphthalene and their responses compared with results from Tenax extraction and chemical analysis. Results show a 1:1 correlation between biosensor results and chemical analyses for naphthalene-contaminated model materials and sediments, but the biosensor assay is much faster. This work demonstrates that biosensor technology can perform as well as standard chemical methods, though with some advantages including the inherent biological relevance of the response, rapid response time, and potential for field deployment. A survey of results from this work and the literature shows that bioavailability under non-equilibrium conditions nonetheless correlates well with K(oc) or K(d). A rationale is provided wherein chemical resistance is speculated to be operative.

La veille bibliographique. Partie 1, Un élément de surveillance

Dans cette première de quatre parties, un réseau francophone multidisciplinaire présente les principaux résultats d'une veille bibliographique sur la surveillance biologique de l'exposition aux produits chimiques en milieu de travail.

Tracing of the Rhône River, wastewater, and mixing within Lake Geneva : stable isotope compositions of water and dissolved inorganic carbon

This study presents an evaluation of the stable isotopic composition of water (hydrogen and oxygen) and dissolved inorganic carbon (DIC) of Lake Geneva, a deep, peri-alpine lake situated at the border between Switzerland and France. The research goal is to apply vertical and seasonal variations of the isotope compositions to evaluate mixing processes of pollutants, nutrients and oxygen. Depth profiles were sampled at different locations throughout Lake Geneva on a monthly and seasonal basis over the course of three years (2009-2011). The results of the oxygen isotopic composition indicate a Rhône River interflow, which can be traced for about 55 km throughout the lake during summer. The Rhône River interflow is 7 to 15 m thick and the molar fraction of Rhône water is estimated to amount up to 37 %. Calculated density of the water and measured isotopic compositions demonstrate that the interflow depth changes in conjunction with the density gradient in the water column during fall. Partial pressure of CO2 indicates that the epilimnion is taking up CO2 from the atmosphere between spring and fall. The epilimnion is most enriched in 13CDIC in September and a progressive depletion of 13CDIC can be observed in the metalimnion from spring to late fall. This stratification is dependent on the local density stratification and the results demonstrate that parameters, which are indicating photosynthesis, are not necessarily linked to δ13CDIC values. In addition, the amount of primary production shows a strong discrepancy between summer 2009 and 2010, but δ13CDIC values of the epilimnion and metalimnion do not indicate variations. In the hypolimnion of the deep lake δ13CDIC values are constant and the progressive depletion allows tracing remineralization processes. The combination of stable carbon and oxygen isotopic compositions allows furthermore tracing Rhône River water fractions, as well as wastewater, stormwater and anthropogenic induced carbon in the water column of the shallow Bay of Vidy. In combination with the results of measured micropollutants, the study underlines that concentrations of certain substances may be related to the Rhône River interflow and/or remineralization of particulate organic carbon. Water quality monitoring and research should therefore be extended to the metalimnion as well as sediment water interface.

Childhood cancer and nuclear power plants in Switzerland: a census-based cohort study.

BACKGROUND: Previous studies on childhood cancer and nuclear power plants (NPPs) produced conflicting results. We used a cohort approach to examine whether residence near NPPs was associated with leukaemia or any childhood cancer in Switzerland. METHODS: We computed person-years at risk for children aged 0-15 years born in Switzerland from 1985 to 2009, based on the Swiss censuses 1990 and 2000 and identified cancer cases from the Swiss Childhood Cancer Registry. We geo-coded place of residence at birth and calculated incidence rate ratios (IRRs) with 95% confidence intervals (CIs) comparing the risk of cancer in children born <5 km, 5-10 km and 10-15 km from the nearest NPP with children born >15 km away, using Poisson regression models. RESULTS: We included 2925 children diagnosed with cancer during 21 117 524 person-years of follow-up; 953 (32.6%) had leukaemia. Eight and 12 children diagnosed with leukaemia at ages 0-4 and 0-15 years, and 18 and 31 children diagnosed with any cancer were born <5 km from a NPP. Compared with children born >15 km away, the IRRs (95% CI) for leukaemia in 0-4 and 0-15 year olds were 1.20 (0.60-2.41) and 1.05 (0.60-1.86), respectively. For any cancer, corresponding IRRs were 0.97 (0.61-1.54) and 0.89 (0.63-1.27). There was no evidence of a dose-response relationship with distance (P > 0.30). Results were similar for residence at diagnosis and at birth, and when adjusted for potential confounders. Results from sensitivity analyses were consistent with main results. CONCLUSIONS: This nationwide cohort study found little evidence of an association between residence near NPPs and the risk of leukaemia or any childhood cancer.

Short-term increase in particulate matter blunts nocturnal blood pressure dipping and daytime urinary sodium excretion

Short-term exposure to ambient particulate matter with aerodynamic diameters<10 µm were found to be positively associated with blood pressure. Yet, little information exists regarding the association between particles and circadian rhythm of blood pressure. Hence, we analyzed the association of exposure to particulate matter with aerodynamic diameters<10 µm on the day of examination and ≤7 days before with ambulatory blood pressure and with sodium excretion in 359 adults from the general population using multiple linear regression. After controlling for potential confounders, a 10-µg/m3 increase in particulate matter with aerodynamic diameters<10 µm levels was associated with nighttime systolic blood pressure (β=1.32 mm Hg 95% CI, 0.06-2.58 mm Hg at lag 0; P=0.04), nighttime diastolic blood pressure (0.72 mm Hg 95% CI, 0.03-1.42 mm Hg at lag 2; P=0.04), nocturnal systolic blood pressure dipping (-0.96 mm Hg 95% CI, -1.89 to -0.03 mm Hg at lag 0; P=0.044), and daytime urinary sodium excretion (-0.05 log-mmol/min 95% CI, -0.10 to -0.01 log-mmol/min at lag 0; P=0.027) but not with nighttime sodium excretion. The associations with blood pressure rapidly diminished with increasing lag days, and the associations with daytime sodium excretion were maximal with particulate matter with aerodynamic diameters<10 µm in exposures 2 to 5 days before. The associations of short-term increases in particulate matter with aerodynamic diameters<10 µm with higher nighttime blood pressure and blunted systolic blood pressure dipping were preceded by associations with reduced ability of the kidney to excrete sodium during daytime. The underlying mechanism linking air pollution to increased cardiovascular risk may include disturbed circadian rhythms of renal sodium handling and blood pressure.

Bacterial variations on the methionine salvage pathway.

BACKGROUND: The thiomethyl group of S-adenosylmethionine is often recycled as methionine from methylthioadenosine. The corresponding pathway has been unravelled in Bacillus subtilis. However methylthioadenosine is subjected to alternative degradative pathways depending on the organism. RESULTS: This work uses genome in silico analysis to propose methionine salvage pathways for Klebsiella pneumoniae, Leptospira interrogans, Thermoanaerobacter tengcongensis and Xylella fastidiosa. Experiments performed with mutants of B. subtilis and Pseudomonas aeruginosa substantiate the hypotheses proposed. The enzymes that catalyze the reactions are recruited from a variety of origins. The first, ubiquitous, enzyme of the pathway, MtnA (methylthioribose-1-phosphate isomerase), belongs to a family of proteins related to eukaryotic intiation factor 2B alpha. mtnB codes for a methylthioribulose-1-phosphate dehydratase. Two reactions follow, that of an enolase and that of a phosphatase. While in B. subtilis this is performed by two distinct polypeptides, in the other organisms analyzed here an enolase-phosphatase yields 1,2-dihydroxy-3-keto-5-methylthiopentene. In the presence of dioxygen an aci-reductone dioxygenase yields the immediate precursor of methionine, ketomethylthiobutyrate. Under some conditions this enzyme produces carbon monoxide in B. subtilis, suggesting a route for a new gaseous mediator in bacteria. Ketomethylthiobutyrate is finally transaminated by an aminotransferase that exists usually as a broad specificity enzyme (often able to transaminate aromatic aminoacid keto-acid precursors or histidinol-phosphate). CONCLUSION: A functional methionine salvage pathway was experimentally demonstrated, for the first time, in P. aeruginosa. Apparently, methionine salvage pathways are frequent in Bacteria (and in Eukarya), with recruitment of different polypeptides to perform the needed reactions (an ancestor of a translation initiation factor and RuBisCO, as an enolase, in some Firmicutes). Many are highly dependent on the presence of oxygen, suggesting that the ecological niche may play an important role for the existence and/or metabolic steps of the pathway, even in phylogenetically related bacteria. Further work is needed to uncover the corresponding steps when dioxygen is scarce or absent (this is important to explore the presence of the pathway in Archaea). The thermophile T. tengcongensis, that thrives in the absence of oxygen, appears to possess the pathway. It will be an interesting link to uncover the missing reactions in anaerobic environments.

Graphitized carbon black in quartz tubes for the sampling of indoor air nicotine and analysis by microwave thermal desorption-capillary gas chromatography

Nicotine in a smoky indoor air environment can be determined using graphitized carbon black as a solid sorbent in quartz tubes. The temperature stability, high purity, and heat absorption characteristics of the sorbent, as well as the permeability of the quartz tubes to microwaves, enable the thermal desorption by means of microwaves after active sampling. Permeation and dynamic dilution procedures for the generation of nicotine in the vapor phase at low and high concentrations are used to evaluate the performances of the sampler. Tube preparation is described and the microwave desorption temperature is measured. Breakthrough volume is determined to allow sampling at 0.1-1 L/min for definite periods of time. The procedure is tested for the determination of gas and paticulate phase nicotine in sidestream smoke produced in an experimental chamber.

Miniaturized bacterial biosensor system for arsenic detection holds great promise for making integrated measurement device.

Combining bacterial bioreporters with microfluidics systems holds great promise for in-field detection of chemical or toxicity targets. Recently we showed how Escherichia coli cells engineered to produce a variant of green fluorescent protein after contact to arsenite and arsenate can be encapsulated in agarose beads and incorporated into a microfluidic chip to create a device for in-field detection of arsenic, a contaminant of well known toxicity and carcinogenicity in potable water both in industrialized and developing countries. Cell-beads stored in the microfluidics chip at -20°C retained inducibility up to one month and we were able to reproducibly discriminate concentrations of 10 and 50 μg arsenite per L (the drinking water standards for European countries and the United States, and for the developing countries, respectively) from the blank in less than 200 minutes. We discuss here the reasons for decreasing bioreporter signal development upon increased storage of cell beads but also show how this decrease can be reduced, leading to a faster detection and a longer lifetime of the device.

La veille bibliographique. Partie 2, Un élément de surveillance

Cet article est la seconde partie d'une série de quatre consacrée aux résultats d'une veille bibliographique sur la surveillance biologique de l'exposition aux produits chimiques en milieu de travail (SBEPC MT). Alors que la précédente partie présentait les objectifs et l'organisation de la veille, cette partie ainsi que la partie 3 vont donc présenter une vue d'ensemble de la base de données en fonction de l'indexation des articles analysés par différents mots clés.

A Monte Carlo approach to the inverse problem of diffuse pollution risk in agricultural catchments

The hydrological and biogeochemical processes that operate in catchments influence the ecological quality of freshwater systems through delivery of fine sediment, nutrients and organic matter. Most models that seek to characterise the delivery of diffuse pollutants from land to water are reductionist. The multitude of processes that are parameterised in such models to ensure generic applicability make them complex and difficult to test on available data. Here, we outline an alternative - data-driven - inverse approach. We apply SCIMAP, a parsimonious risk based model that has an explicit treatment of hydrological connectivity. we take a Bayesian approach to the inverse problem of determining the risk that must be assigned to different land uses in a catchment in order to explain the spatial patterns of measured in-stream nutrient concentrations. We apply the model to identify the key sources of nitrogen (N) and phosphorus (P) diffuse pollution risk in eleven UK catchments covering a range of landscapes. The model results show that: 1) some land use generates a consistently high or low risk of diffuse nutrient pollution; but 2) the risks associated with different land uses vary both between catchments and between nutrients; and 3) that the dominant sources of P and N risk in the catchment are often a function of the spatial configuration of land uses. Taken on a case-by-case basis, this type of inverse approach may be used to help prioritise the focus of interventions to reduce diffuse pollution risk for freshwater ecosystems. (C) 2012 Elsevier B.V. All rights reserved.