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.

Nuclear translocation and DNA recognition signals colocalized within the bZIP domain of cyclic adenosine 3',5'-monophosphate response element-binding protein CREB.

CREB is a cAMP-responsive nuclear DNA-binding protein that binds to cAMP response elements and stimulates gene transcription upon activation of the cAMP signalling pathway. The protein consists of an amino-terminal transcriptional transactivation domain and a carboxyl-terminal DNA-binding domain (bZIP domain) comprised of a basic region and a leucine zipper involved in DNA recognition and dimerization, respectively. Recently, we discovered a testis-specific transcript of CREB that contains an alternatively spliced exon encoding multiple stop codons. CREB encoded by this transcript is a truncated protein lacking the bZIP domain. We postulated that the antigen detected by CREB antiserum in the cytoplasm of germinal cells is the truncated CREB that must also lack its nuclear translocation signal (NTS). To test this hypothesis we prepared multiple expression plasmids encoding carboxyl-terminal deletions of CREB and transiently expressed them in COS-1 cells. By Western immunoblot analysis as well as immunocytochemistry of transfected cells, we show that CREB proteins truncated to amino acid 286 or shorter are sequestered in the cytoplasm, whereas a CREB of 295 amino acids is translocated into the nucleus. Chimeric CREBs containing a heterologous NTS fused to the first 248 or 261 amino acids of CREB are able to drive the translocation of the protein into the nucleus. Thus, the nine amino acids in the basic region involved in DNA recognition between positions 287 and 295 (RRKKKEYVK) of CREB contain the NTS. Further, mutation of the lysine at position 290 in CREB to an asparagine diminishes nuclear translocation of the protein.(ABSTRACT TRUNCATED AT 250 WORDS)

Associations of short-term particle and noise exposures with markers of cardiovascular and respiratory health among highway maintenance workers

BACKGROUND: Highway maintenance workers are constantly and simultaneously exposed to traffic-related particle and noise emissions, and both have been linked to increased cardiovascular morbidity and mortality in population-based epidemiology studies. OBJECTIVES: We aimed to investigate short-term health effects related to particle and noise exposure. METHODS: We monitored 18 maintenance workers, during as many as five 24-hour periods from a total of 50 observation days. We measured their exposure to fine particulate matter (PM2.5), ultrafine particles, noise, and the cardiopulmonary health endpoints: blood pressure, pro-inflammatory and pro-thrombotic markers in the blood, lung function and fractional exhaled nitric oxide (FeNO) measured approximately 15 hours post-work. Heart rate variability was assessed during a sleep period approximately 10 hours post-work. RESULTS: PM2.5 exposure was significantly associated with C-reactive protein and serum amyloid A, and negatively associated with tumor necrosis factor α. None of the particle metrics were significantly associated with von Willebrand factor or tissue factor expression. PM2.5 and work noise were associated with markers of increased heart rate variability, and with increased HF and LF power. Systolic and diastolic blood pressure on the following morning were significantly associated with noise exposure after work, and non-significantly associated with PM2.5. We observed no significant associations between any of the exposures and lung function or FeNO. CONCLUSIONS: Our findings suggest that exposure to particles and noise during highway maintenance work might pose a cardiovascular health risk. Actions to reduce these exposures could lead to better health for this population of workers.

Simultaneous LC-MS/MS quantification of P-glycoprotein and cytochrome P450 probe substrates and their metabolites in DBS and plasma.

BACKGROUND: An LC-MS/MS method has been developed for the simultaneous quantification of P-glycoprotein (P-gp) and cytochrome P450 (CYP) probe substrates and their Phase I metabolites in DBS and plasma. P-gp (fexofenadine) and CYP-specific substrates (caffeine for CYP1A2, bupropion for CYP2B6, flurbiprofen for CYP2C9, omeprazole for CYP2C19, dextromethorphan for CYP2D6 and midazolam for CYP3A4) and their metabolites were extracted from DBS (10 µl) using methanol. Analytes were separated on a reversed-phase LC column followed by SRM detection within a 6 min run time. RESULTS: The method was fully validated over the expected clinical concentration range for all substances tested, in both DBS and plasma. The method has been successfully applied to a PK study where healthy male volunteers received a low dose cocktail of the here described P-gp and CYP probes. Good correlation was observed between capillary DBS and venous plasma drug concentrations. CONCLUSION: Due to its low-invasiveness, simple sample collection and minimal sample preparation, DBS represents a suitable method to simultaneously monitor in vivo activities of P-gp and CYP.

Hyperpolarizing gases via dynamic nuclear polarization and sublimation.

A high throughput method was designed to produce hyperpolarized gases by combining low-temperature dynamic nuclear polarization with a sublimation procedure. It is illustrated by applications to 129Xe nuclear magnetic resonance in xenon gas, leading to a signal enhancement of 3 to 4 orders of magnitude compared to the room-temperature thermal equilibrium signal at 7.05 T.

Effect of acute and chronic exercise on fat oxidation and hormonal and plasma metabolite kinetics in obese and severely obese men

Acute exercise increases energy expenditure (EE) during exercise and post-exercise recovery [excess post-exercise oxygen consumption (EPOC)] and therefore may be recommended as part of the multidisciplinary management of obesity. Moreover, chronic exercise (training) effectively promotes an increase in insulin sensitivity, which seems to be associated with increased fat oxidation rates (FORs). The main purpose of this thesis is to investigate 1) FORs and extra-muscular factors (hormones and plasma metabolites) that regulate fat metabolism during acute and chronic exercise; and 2) EPOC during acute post-exercise recovery in obese and severely obese men (class II and III). In the first study, we showed that obese and severely obese men present a lower exercise intensity (Fatmax) eliciting maximal fat oxidation and a lower reliance on fat oxidation at high, but not at low and moderate, exercise intensities compared to lean men. This was most likely related to an impaired muscular capacity to oxidize non-esterified fatty acids (NEFA) rather than decreased plasma NEFA availability or a change in the hormonal milieu during exercise. In the second study, we developed an accurate maximal incremental test to correctly and simultaneously evaluate aerobic fitness and fat oxidation kinetics during exercise in this population. This test may be used for the prescription of an appropriate exercise training intensity. In the third study, we demonstrated that only 2 wk of exercise training [continuous training at Fatmax and adapted high-intensity interval training (HIIT)], matched with respect to mechanical work, may be effective to improve aerobic fitness, FORs during exercise and insulin sensitivity, which suggest that FORs might be rapidly improved and that adapted HIIT is feasible in this population. The increased FORs concomitant with the lack of changes in lipolysis during exercise suggest an improvement in the mismatching between NEFA availability and oxidation, highlighting the importance of muscular (oxidative capacity) rather than extra-muscular (hormones and plasma metabolites) factors in the regulation of fat metabolism after a training program. In the fourth study, we observed a positive correlation between EE during exercise and EPOC, suggesting that a chronic increase in the volume or intensity of exercise may increase EE during exercise and during recovery. This may have an impact in weight management in obesity. In conclusion, these findings might have practical implications for exercise training prescriptions in order to improve the therapeutic approaches in obesity and severe obesity. -- L'exercice aigu augmente la dépense énergétique (DE) pendant l'exercice et la récupération post-exercice [excès de consommation d'oxygène post-exercise (EPOC)] et peut être utilisé dans la gestion multidisciplinaire de l'obésité. Quant à l'exercice chronique (entraînement), il est efficace pour augmenter la sensibilité à l'insuline, ce qui semble être associé à une amélioration du débit d'oxydation lipidique (DOL). Le but de cette thèse est d'étudier 1) le DOL et les facteurs extra-musculaires (hormones et métabolites plasmatiques) qui régulent le métabolisme lipidique pendant l'exercice aigu et chronique et 2) l'EPOC lors de la récupération aiguë post-exercice chez des hommes obèses et sévèrement obèses (classe II et III). Dans la première étude nous avons montré que les hommes obèses et sévèrement obèses présentent une plus basse intensité d'exercice (Fatmax) correspondant au débit d'oxydation lipidique maximale et un plus bas DOL à hautes, mais pas à faibles et modérées, intensités d'exercice comparé aux sujets normo-poids, ce qui est probablement lié à une incapacité musculaire à oxyder les acides gras non-estérifiés (AGNE) plutôt qu'à une diminution de leur disponibilité ou à un changement du milieu hormonal pendant l'exercice. Dans la deuxième étude nous avons développé un test maximal incrémental pour évaluer simultanément l'aptitude physique aérobie et la cinétique d'oxydation des lipides pendant l'exercice chez cette population. Dans la troisième étude nous avons montré que seulement deux semaines d'entraînement (continu à Fatmax et intermittent à haute intensité), appariés par la charge de travail, sont efficaces pour améliorer l'aptitude physique aérobie, le DOL pendant l'exercice et la sensibilité à l'insuline, ce qui suggère que le DOL peut être rapidement amélioré chez cette population. Ceci, en absence de changements de la lipolyse pendant l'exercice, suggère une amélioration de la balance entre la disponibilité et l'oxydation des AGNE, ce qui souligne l'importance des facteurs musculaires (capacité oxydative) plutôt que extra-musculaires (hormones et métabolites plasmatiques) dans la régulation du métabolisme lipidique après un entraînement. Dans la quatrième étude nous avons observé une corrélation positive entre la DE pendant l'exercice et l'EPOC, ce qui suggère qu'une augmentation chronique du volume ou de l'intensité de l'exercice pourrait augmenter la DE lors de l'exercice et lors de la récupération post-exercice. Ceci pourrait avoir un impact sur la gestion du poids chez cette population. En conclusion, ces résultats pourraient avoir des implications pratiques lors de la prescription des entraînements dans le but d'améliorer les approches thérapeutiques de l'obésité et de l'obésité sévère.

Social meets molecular: Combining phylogenetic and latent class analyses to understand HIV-1 transmission in Switzerland.

Switzerland has a complex human immunodeficiency virus (HIV) epidemic involving several populations. We examined transmission of HIV type 1 (HIV-1) in a national cohort study. Latent class analysis was used to identify socioeconomic and behavioral groups among 6,027 patients enrolled in the Swiss HIV Cohort Study between 2000 and 2011. Phylogenetic analysis of sequence data, available for 4,013 patients, was used to identify transmission clusters. Concordance between sociobehavioral groups and transmission clusters was assessed in correlation and multiple correspondence analyses. A total of 2,696 patients were infected with subtype B, 203 with subtype C, 196 with subtype A, and 733 with recombinant subtypes (mainly CRF02_AG and CRF01_AE). Latent class analysis identified 8 patient groups. Most transmission clusters of subtype B were shared between groups of gay men (groups 1-3) or between the heterosexual groups "heterosexual people of lower socioeconomic position" (group 4) and "injection drug users" (group 8). Clusters linking homosexual and heterosexual groups were associated with "older heterosexual and gay people on welfare" (group 5). "Migrant women in heterosexual partnerships" (group 6) and "heterosexual migrants on welfare" (group 7) shared non-B clusters with groups 4 and 5. Combining approaches from social and molecular epidemiology can provide insights into HIV-1 transmission and inform the design of prevention strategies.

Simultaneous determination of antidementia drugs by HPLC-MS: validation data of the method and plasma level determinations in 300 patients

Aims: A rapid and simple HPLC-MS method was developed for the simultaneousdetermination of antidementia drugs, including donepezil, galantamine, rivastigmineand its major metabolite NAP 226 - 90, and memantine, for TherapeuticDrug Monitoring (TDM). In the elderly population treated with antidementiadrugs, the presence of several comorbidities, drug interactions resulting frompolypharmacy, and variations in drug metabolism and elimination, are possiblefactors leading to the observed high interindividual variability in plasma levels.Although evidence for the benefit of TDM for antidementia drugs still remains tobe demonstrated, an individually adapted dosage through TDM might contributeto minimize the risk of adverse reactions and to increase the probability of efficienttherapeutic response. Methods: A solid-phase extraction procedure with amixed-mode cation exchange sorbent was used to isolate the drugs from 0.5 mL ofplasma. The compounds were analyzed on a reverse-phase column with a gradientelution consisting of an ammonium acetate buffer at pH 9.3 and acetonitrile anddetected by mass spectrometry in the single ion monitoring mode. Isotope-labeledinternal standards were used for quantification where possible. The validatedmethod was used to measure the plasma levels of antidementia drugs in 300patients treated with these drugs. Results: The method was validated accordingto international standards of validation, including the assessment of the trueness(-8 - 11 %), the imprecision (repeatability: 1-5%, intermediate imprecision:2 - 9 %), selectivity and matrix effects variability (less than 6 %). Furthermore,short and long-term stability of the analytes in plasma was ascertained. Themethod proved to be robust in the calibrated ranges of 1 - 300 ng/mL for rivastigmineand memantine and 2 - 300 mg/mL for donepezil, galantamine and NAP226 - 90. We recently published a full description of the method (1). We found ahigh interindividual variability in plasma levels of these drugs in a study populationof 300 patients. The plasma level measurements, with some preliminaryclinical and pharmacogenetic results, will be presented. Conclusion: A simpleLC-MS method was developed for plasma level determination of antidementiadrugs which was successfully used in a clinical study with 300 patients.

Effect of early antiretroviral therapy during primary HIV-1 infection on cell-associated HIV-1 DNA and plasma HIV-1 RNA.

Background: Early initiation of combination antiretroviral therapy (ART) during primary HIV-1 infection may prevent the establishment of large viral reservoirs, possibly resulting in improved control of plasma viraemia rebound after ART cessation.Methods: Levels of cell-associated HIV-1 DNA and plasma HIV-1 RNA were measured longitudinally in 32 acutely and recently infected patients, who started ART <= 120 days after the estimated date of infection, and interrupted ART after 18 months (median) of continuous therapy. Averages of HIV-1 DNA and RNA concentrations present in blood 30-365 days after therapy interruption (median duration 300 days, range 195-358) were compared between patients who started ART <= 60 days after the estimated date of infection (early starters), those who started between 61 and 120 days (later starters), and, for HIV-1 RNA only, with 89 untreated participants of the Swiss HIV Cohort Study with documented sero-conversion and longitudinal measurements collected 90-455 days after the first positive HIV test.Results: In early ART starters, average levels of plasma HIV-1 RNA and cell-associated HIV-1 DNA after treatment interruption were 1 log(10) (P=0.008) and 0.4 log(10) (P=0.03) lower compared with later starters. Average post-treatment plasma HIV-1 RNA levels in early starters were significantly lower, respectively, compared with untreated controls (-1.2 log(10); P<0.0004).Conclusions: Early treatment initiation within 2 months after HIV infection compared with later therapy initiation resulted in reduced levels of plasma viraemia and proviral HIV-1 DNA for >= 1 year after subsequent ART cessation. Plasma HIV-1 RNA levels in early starters were also significantly lower than in untreated controls.

Intraspecific genetic variation in arbuscular mycorrhizal fungi and its consequences for molecular biology, ecology, and development of inoculum

It has been known for some time that different arbuscular mycorrhizal fungal (AMF) taxa confer differences in plant growth. Although genetic variation within AMF species has been given less attention, it could potentially be an ecologically important source of variation. Ongoing studies on variability in AMF genes within Glomus intraradices indicate that at least for some genes, such as the BiP gene, sequence variability can be high, even in coding regions. This suggests that genetic variation within an AMF may not be selectively neutral. This clearly needs to be investigated in more detail for other coding regions of AMF genomes. Similarly, studies on AMF population genetics indicate high genetic variation in AMF populations, and a considerable amount of variation seen in phenotypes in the population can be attributed to genetic differences among the fungi. The existence of high within-species genetic variation could have important consequences for how investigations on AMF gene expression and function are conducted. Furthermore, studies of within-species genetic variability and how it affects variation in plant growth will help to identify at what level of precision ecological studies should be conducted to identify AMF in plant roots in the field. A population genetic approach to studying AMF genetic variability can also be useful for inoculum development. By knowing the amount of genetic variability in an AMF population, the maximum and minimum numbers of spores that will contain a given amount of genetic diversity can be estimated. This could be particularly useful for developing inoculum with high adaptability to different environments.