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.

Regulation of human class I alcohol dehydrogenases by bile acids

Class I alcohol dehydrogenases (ADH1s) are the rate-limiting enzymes for ethanol and vitamin A (retinol) metabolism in the liver . Because previous studies have shown that human ADH1 enzymes may participate in bile acid metabolism, we investigated whether the bile acid-activated nuclear receptor farnesoid X receptor (FXR) regulates ADH1 genes. In human hepatocytes, both the endogenous FXR ligand chenodeoxycholic acid and synthetic FXR-specific agonist GW4064 increased ADH1 mRNA, protein, and activity. Moreover, overexpression of a constitutively active form of FXR induced ADH1A and ADH1B expression, whereas silencing of FXR abolished the effects of FXR agonists on ADH1 expression and activity. Transient transfection studies and electrophoretic mobility shift assays revealed functional FXR response elements in the ADH1A and ADH1B proximal promoters, thus indicating that both genes are direct targets of FXR. These findings provide the first evidence for direct connection of bile acid signaling and alcohol metabolism.

Common genetic loci influencing plasma homocysteine concentrations and their effect on risk of coronary artery disease.

BACKGROUND: The strong observational association between total homocysteine (tHcy) concentrations and risk of coronary artery disease (CAD) and the null associations in the homocysteine-lowering trials have prompted the need to identify genetic variants associated with homocysteine concentrations and risk of CAD. OBJECTIVE: We tested whether common genetic polymorphisms associated with variation in tHcy are also associated with CAD. DESIGN: We conducted a meta-analysis of genome-wide association studies (GWAS) on tHcy concentrations in 44,147 individuals of European descent. Polymorphisms associated with tHcy (P < 10(-8)) were tested for association with CAD in 31,400 cases and 92,927 controls. RESULTS: Common variants at 13 loci, explaining 5.9% of the variation in tHcy, were associated with tHcy concentrations, including 6 novel loci in or near MMACHC (2.1 Ã- 10(-9)), SLC17A3 (1.0 Ã- 10(-8)), GTPB10 (1.7 Ã- 10(-8)), CUBN (7.5 Ã- 10(-10)), HNF1A (1.2 Ã- 10(-12)), and FUT2 (6.6 Ã- 10(-9)), and variants previously reported at or near the MTHFR, MTR, CPS1, MUT, NOX4, DPEP1, and CBS genes. Individuals within the highest 10% of the genotype risk score (GRS) had 3-μmol/L higher mean tHcy concentrations than did those within the lowest 10% of the GRS (P = 1 Ã- 10(-36)). The GRS was not associated with risk of CAD (OR: 1.01; 95% CI: 0.98, 1.04; P = 0.49). CONCLUSIONS: We identified several novel loci that influence plasma tHcy concentrations. Overall, common genetic variants that influence plasma tHcy concentrations are not associated with risk of CAD in white populations, which further refutes the causal relevance of moderately elevated tHcy concentrations and tHcy-related pathways for CAD.

Transgene inheritances and genetic similarities of near isogenic lines of genetically modified common beans

The objective of the present work was to determine the inheritance and stability of transgenes of a transgenic bean line expressing the genes rep-trap-ren from Bean golden mosaic virus and the bar gene. Crosses were done between the transgenic line and four commercial bean cultivars, followed by four backcrosses to the commercial cultivars. Progenies from each cross were evaluated for the presence of the transgenes by brushing the leaves with glufosinate ammonium and by polymerase chain reaction using specific oligonucleotides. Advanced generations were rub-inoculated with an isolate of Bean common mosaic necrosis virus (BCMNV). The transgenes were inherited consistently in a Mendelian pattern in the four crosses studied. The analyzed lines recovered close to 80% of the characteristics of the recurrent parent, as determined by the random amplified DNA markers used, besides maintaining important traits such as resistance to BCMNV. The presence of the transgene did not cause any detectable undesirable effect in the evaluated progenies.

Meiotic drive favors Robertsonian metacentric chromosomes in the common shrew (Sorex araneus, Insectivora, mammalia).

Meiotic drive has attracted much interest because it concerns the robustness of Mendelian segregation and its genetic and evolutionary stability. We studied chromosomal meiotic drive in the common shrew (Sorex araneus, Insectivora, Mammalia), which exhibits one of the most remarkable chromosomal polymorphisms within mammalian species. The open question of the evolutionary success of metacentric chromosomes (Robertsonian fusions) versus acrocentrics in the common shrew prompted us to test whether a segregation distortion in favor of metacentrics is present in female and/or male meiosis. Performing crosses under controlled laboratory conditions with animals from natural populations, we found a clear trend toward a segregation distortion in favor of metacentrics during male meiosis, two chromosome combinations (gm and jl) being significantly preferred over their acrocentric homologs. Apart for one Robertsonian fusion (hi), this trend was absent in female meiosis. We propose a model based on recombination events between twin acrocentrics to explain the difference in transmission ratios of the same metacentric in different sexes and unequal drive of particular metacentrics in the same sex. Pooled data for female and male meiosis revealed a trend toward stronger segregation distortion for larger metacentrics. This is partially in agreement with the frequency of metacentrics occurring in natural populations of a chromosome race showing a high degree of chromosomal polymorphism.

ERK1/2 controls Na,K-ATPase activity and transepithelial sodium transport in the principal cell of the cortical collecting duct of the mouse kidney.

The collecting duct of normal kidney exhibits significant activity of the MEK1/2-ERK1/2 pathway as shown in vivo by immunostaining of phosphorylated active ERK1/2 (pERK1/2). The MEK1/2-ERK1/2 pathway controls many different ion transports both in proximal and distal nephron, raising the question of whether this pathway is involved in the basal and/or hormone-dependent transepithelial sodium reabsorption in the principal cell of the cortical collecting duct (CCD), a process mediated by the apical epithelial sodium channel and the basolateral sodium pump (Na,K-ATPase). To answer this question we used ex vivo microdissected CCDs from normal mouse kidney or in vitro cultured mpkCCDcl4 principal cells. Significant basal levels of pERK1/2 were observed ex vivo and in vitro. Aldosterone and vasopressin, known to up-regulate sodium reabsorption in CCDs, did not change ERK1/2 activity either ex vivo or in vitro. Basal and aldosterone- or vasopressin-stimulated sodium transport was down-regulated by the MEK1/2 inhibitor PD98059, in parallel with a decrease in pERK1/2 in vitro. The activity of Na,K-ATPase but not that of epithelial sodium channel was inhibited by MEK1/2 inhibitors in both unstimulated and aldosterone- or vasopressin-stimulated CCDs in vitro. Cell surface biotinylation showed that intrinsic activity rather than cell surface expression of Na,K-ATPase was controlled by pERK1/2. PD98059 also significantly inhibited the activity of Na,K-ATPase ex vivo. Our data demonstrate that the ERK1/2 pathway controls Na,K-ATPase activity and transepithelial sodium transport in the principal cell and indicate that basal constitutive activity of the ERK1/2 pathway is a critical component of this control.

The Epidemiology of Common Communicable Diseases, January 2016

Chart of common communicable disease, stating incubation, transmission, control measures, and public health response produced by Iowa Department of Public Health.

Phytogeographical origin of Madeiran common beans based on phaseolin patterns

The objective of this work was to determine the geographic origin of the Madeiran common bean (Phaseolus vulgaris) gene pool. Phaseolin patterns of 50 accessions representing the diversity of common bean collected in Madeira, Portugal, and conserved in the ISOPlexis Germplasm Bank, were analysed using the Experion automated electrophoresis system, based on lab-on-a-chip technology. Five common bean standard varieties with typical phaseolin patterns were used to determine the phytogeographical origin of the Madeiran common bean accessions. Ninety two percent of the accessions exhibited a phaseolin pattern consistent with the one of common bean types belonging to the Andean gene pool, while the origin of the remaining 8% of the accessions was indistinguishable. The application of a similarity coefficient of 85%, based on Pearson correlations, increases the number of accessions with uncertain pattern. The analytical approach used permitted the determination of the origin of the common bean gene pool, which is Andean in 98% of the cases, and clustering of the observed variability among the Madeiran common beans.

Overlap Between Common Genetic Polymorphisms Underpinning Kidney Traits and Cardiovascular Disease Phenotypes: The CKDGen Consortium.

BACKGROUND: Chronic kidney disease is associated with cardiovascular disease. We tested for evidence of a shared genetic basis to these traits. STUDY DESIGN: We conducted 2 targeted analyses. First, we examined whether known single-nucleotide polymorphisms (SNPs) underpinning kidney traits were associated with a series of vascular phenotypes. Additionally, we tested whether vascular SNPs were associated with markers of kidney damage. Significance was set to 1.5×10(-4) (0.05/325 tests). SETTING & PARTICIPANTS: Vascular outcomes were analyzed in participants from the AortaGen (20,634), CARDIoGRAM (86,995), CHARGE Eye (15,358), CHARGE IMT (31,181), ICBP (69,395), and NeuroCHARGE (12,385) consortia. Tests for kidney outcomes were conducted in up to 67,093 participants from the CKDGen consortium. PREDICTOR: We used 19 kidney SNPs and 64 vascular SNPs. OUTCOMES & MEASUREMENTS: Vascular outcomes tested were blood pressure, coronary artery disease, carotid intima-media thickness, pulse wave velocity, retinal venular caliber, and brain white matter lesions. Kidney outcomes were estimated glomerular filtration rate and albuminuria. RESULTS: In general, we found that kidney disease variants were not associated with vascular phenotypes (127 of 133 tests were nonsignificant). The one exception was rs653178 near SH2B3 (SH2B adaptor protein 3), which showed direction-consistent association with systolic (P = 9.3 ×10(-10)) and diastolic (P = 1.6 ×10(-14)) blood pressure and coronary artery disease (P = 2.2 ×10(-6)), all previously reported. Similarly, the 64 SNPs associated with vascular phenotypes were not associated with kidney phenotypes (187 of 192 tests were nonsignificant), with the exception of 2 high-correlated SNPs at the SH2B3 locus (P = 1.06 ×10(-07) and P = 7.05 ×10(-08)). LIMITATIONS: The combined effect size of the SNPs for kidney and vascular outcomes may be too low to detect shared genetic associations. CONCLUSIONS: Overall, although we confirmed one locus (SH2B3) as associated with both kidney and cardiovascular disease, our primary findings suggest that there is little overlap between kidney and cardiovascular disease risk variants in the overall population. The reciprocal risks of kidney and cardiovascular disease may not be genetically mediated, but rather a function of the disease milieu itself.

Quantification of the diversity among common bean accessions using Ward-MLM strategy

The present work aimed at evaluating the divergence among common bean accessions by their agronomic, morphological and molecular traits, based on the Ward-MLM procedure. A collection of 57 accessions from the gene bank of Universidade Federal do Espírito Santo was used in this study, from which: 31 were landraces belonging to the community Fortaleza, in the municipality of Muqui, ES, Brazil; 20 accessions were provided by Embrapa Trigo; and 6 were commercial cultivars. Five agronomic traits (plant cycle, number of seeds per pod, number of pods per plant, weight of 100 seeds, and grain yield), five morphological traits (growth habit, plant size, seed shape, seed color, and commercial group) and 16 microsatellite primers were evaluated. High genetic variability was detected considering morphological, agronomic and molecular traits in the 57 common bean accessions studied. The Ward-MLM procedure showed that the ideal number of groups was five, according to the pseudo F and pseudo t² criteria. The accessions from Andean origin had heavier seeds than others and formed a cluster. The Ward-MLM statistical procedure is a useful technique to detect genetic divergence and to cluster genotypes by simultaneously using morphological, agronomic and molecular data.

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.