Docking, virtual high throughput screening and in silico fragment-based drug design.

ABSTRACT The drug discovery process has been profoundly changed recently by the adoption of computational methods helping the design of new drug candidates more rapidly and at lower costs. In silico drug design consists of a collection of tools helping to make rational decisions at the different steps of the drug discovery process, such as the identification of a biomolecular target of therapeutical interest, the selection or the design of new lead compounds and their modification to obtain better affinities, as well as pharmacokinetic and pharmacodynamic properties. Among the different tools available, a particular emphasis is placed in this review on molecular docking, virtual high throughput screening and fragment-based ligand design.

Crystal packing modifies ligand binding affinity: the case of aldose reductase.

The relationship between the structures of protein-ligand complexes existing in the crystal and in solution, essential in the case of fragment-based screening by X-ray crystallography (FBS-X), has been often an object of controversy. To address this question, simultaneous co-crystallization and soaking of two inhibitors with different ratios, Fidarestat (FID; K(d) = 6.5 nM) and IDD594 (594; K(d) = 61 nM), which bind to h-aldose reductase (AR), have been performed. The subatomic resolution of the crystal structures allows the differentiation of both inhibitors, even when the structures are almost superposed. We have determined the occupation ratio in solution by mass spectrometry (MS) Occ(FID)/Occ(594) = 2.7 and by X-ray crystallography Occ(FID)/Occ(594) = 0.6. The occupancies in the crystal and in solution differ 4.6 times, implying that ligand binding potency is influenced by crystal contacts. A structural analysis shows that the Loop A (residues 122-130), which is exposed to the solvent, is flexible in solution, and is involved in packing contacts within the crystal. Furthermore, inhibitor 594 contacts the base of Loop A, stabilizing it, while inhibitor FID does not. This is shown by the difference in B-factors of the Loop A between the AR-594 and AR-FID complexes. A stable loop diminishes the entropic energy barrier to binding, favoring 594 versus FID. Therefore, the effect of the crystal environment should be taken into consideration in the X-ray diffraction analysis of ligand binding to proteins. This conclusion highlights the need for additional methodologies in the case of FBS-X to validate this powerful screening technique, which is widely used.

Rational design of indoleamine 2,3-dioxygenase inhibitors.

Indoleamine 2,3-dioxygenase (IDO) is an important therapeutic target for the treatment of diseases such as cancer that involve pathological immune escape. We have used the evolutionary docking algorithm EADock to design new inhibitors of this enzyme. First, we investigated the modes of binding of all known IDO inhibitors. On the basis of the observed docked conformations, we developed a pharmacophore model, which was then used to devise new compounds to be tested for IDO inhibition. We also used a fragment-based approach to design and to optimize small organic molecule inhibitors. Both approaches yielded several new low-molecular weight inhibitor scaffolds, the most active being of nanomolar potency in an enzymatic assay. Cellular assays confirmed the potential biological relevance of four different scaffolds.

EADock : design of a new molecular docking algorithm and some of its applications

3 Summary 3. 1 English The pharmaceutical industry has been facing several challenges during the last years, and the optimization of their drug discovery pipeline is believed to be the only viable solution. High-throughput techniques do participate actively to this optimization, especially when complemented by computational approaches aiming at rationalizing the enormous amount of information that they can produce. In siiico techniques, such as virtual screening or rational drug design, are now routinely used to guide drug discovery. Both heavily rely on the prediction of the molecular interaction (docking) occurring between drug-like molecules and a therapeutically relevant target. Several softwares are available to this end, but despite the very promising picture drawn in most benchmarks, they still hold several hidden weaknesses. As pointed out in several recent reviews, the docking problem is far from being solved, and there is now a need for methods able to identify binding modes with a high accuracy, which is essential to reliably compute the binding free energy of the ligand. This quantity is directly linked to its affinity and can be related to its biological activity. Accurate docking algorithms are thus critical for both the discovery and the rational optimization of new drugs. In this thesis, a new docking software aiming at this goal is presented, EADock. It uses a hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with .the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 R around the center of mass of the ligand position in the crystal structure, and conversely to other benchmarks, our algorithms was fed with optimized ligand positions up to 10 A root mean square deviation 2MSD) from the crystal structure. This validation illustrates the efficiency of our sampling heuristic, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best-ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures in this benchmark could be explained by the presence of crystal contacts in the experimental structure. EADock has been used to understand molecular interactions involved in the regulation of the Na,K ATPase, and in the activation of the nuclear hormone peroxisome proliferatoractivated receptors a (PPARa). It also helped to understand the action of common pollutants (phthalates) on PPARy, and the impact of biotransformations of the anticancer drug Imatinib (Gleevec®) on its binding mode to the Bcr-Abl tyrosine kinase. Finally, a fragment-based rational drug design approach using EADock was developed, and led to the successful design of new peptidic ligands for the a5ß1 integrin, and for the human PPARa. In both cases, the designed peptides presented activities comparable to that of well-established ligands such as the anticancer drug Cilengitide and Wy14,643, respectively. 3.2 French Les récentes difficultés de l'industrie pharmaceutique ne semblent pouvoir se résoudre que par l'optimisation de leur processus de développement de médicaments. Cette dernière implique de plus en plus. de techniques dites "haut-débit", particulièrement efficaces lorsqu'elles sont couplées aux outils informatiques permettant de gérer la masse de données produite. Désormais, les approches in silico telles que le criblage virtuel ou la conception rationnelle de nouvelles molécules sont utilisées couramment. Toutes deux reposent sur la capacité à prédire les détails de l'interaction moléculaire entre une molécule ressemblant à un principe actif (PA) et une protéine cible ayant un intérêt thérapeutique. Les comparatifs de logiciels s'attaquant à cette prédiction sont flatteurs, mais plusieurs problèmes subsistent. La littérature récente tend à remettre en cause leur fiabilité, affirmant l'émergence .d'un besoin pour des approches plus précises du mode d'interaction. Cette précision est essentielle au calcul de l'énergie libre de liaison, qui est directement liée à l'affinité du PA potentiel pour la protéine cible, et indirectement liée à son activité biologique. Une prédiction précise est d'une importance toute particulière pour la découverte et l'optimisation de nouvelles molécules actives. Cette thèse présente un nouveau logiciel, EADock, mettant en avant une telle précision. Cet algorithme évolutionnaire hybride utilise deux pressions de sélections, combinées à une gestion de la diversité sophistiquée. EADock repose sur CHARMM pour les calculs d'énergie et la gestion des coordonnées atomiques. Sa validation a été effectuée sur 37 complexes protéine-ligand cristallisés, incluant 11 protéines différentes. L'espace de recherche a été étendu à une sphère de 151 de rayon autour du centre de masse du ligand cristallisé, et contrairement aux comparatifs habituels, l'algorithme est parti de solutions optimisées présentant un RMSD jusqu'à 10 R par rapport à la structure cristalline. Cette validation a permis de mettre en évidence l'efficacité de notre heuristique de recherche car des modes d'interactions présentant un RMSD inférieur à 2 R par rapport à la structure cristalline ont été classés premier pour 68% des complexes. Lorsque les cinq meilleures solutions sont prises en compte, le taux de succès grimpe à 78%, et 92% lorsque la totalité de la dernière génération est prise en compte. La plupart des erreurs de prédiction sont imputables à la présence de contacts cristallins. Depuis, EADock a été utilisé pour comprendre les mécanismes moléculaires impliqués dans la régulation de la Na,K ATPase et dans l'activation du peroxisome proliferatoractivated receptor a (PPARa). Il a également permis de décrire l'interaction de polluants couramment rencontrés sur PPARy, ainsi que l'influence de la métabolisation de l'Imatinib (PA anticancéreux) sur la fixation à la kinase Bcr-Abl. Une approche basée sur la prédiction des interactions de fragments moléculaires avec protéine cible est également proposée. Elle a permis la découverte de nouveaux ligands peptidiques de PPARa et de l'intégrine a5ß1. Dans les deux cas, l'activité de ces nouveaux peptides est comparable à celles de ligands bien établis, comme le Wy14,643 pour le premier, et le Cilengitide (PA anticancéreux) pour la seconde.

A New Large-DNA-Fragment Delivery System Based on Integrase Activity from an Integrative and Conjugative Element.

During the past few decades, numerous plasmid vectors have been developed for cloning, gene expression analysis, and genetic engineering. Cloning procedures typically rely on PCR amplification, DNA fragment restriction digestion, recovery, and ligation, but increasingly, procedures are being developed to assemble large synthetic DNAs. In this study, we developed a new gene delivery system using the integrase activity of an integrative and conjugative element (ICE). The advantage of the integrase-based delivery is that it can stably introduce a large DNA fragment (at least 75 kb) into one or more specific sites (the gene for glycine-accepting tRNA) on a target chromosome. Integrase recombination activity in Escherichia coli is kept low by using a synthetic hybrid promoter, which, however, is unleashed in the final target host, forcing the integration of the construct. Upon integration, the system is again silenced. Two variants with different genetic features were produced, one in the form of a cloning vector in E. coli and the other as a mini-transposable element by which large DNA constructs assembled in E. coli can be tagged with the integrase gene. We confirmed that the system could successfully introduce cosmid and bacterial artificial chromosome (BAC) DNAs from E. coli into the chromosome of Pseudomonas putida in a site-specific manner. The integrase delivery system works in concert with existing vector systems and could thus be a powerful tool for synthetic constructions of new metabolic pathways in a variety of host bacteria.

Antigenicity and immunogenicity of a novel chimeric peptide antigen based on the P. vivax circumsporozoite protein.

BACKGROUND: Plasmodium vivax circumsporozoite (PvCS) protein is a major sporozoite surface antigen involved in parasite invasion of hepatocytes and is currently being considered as vaccine candidate. PvCS contains a dimorphic central repetitive fragment flanked by conserved regions that contain functional domains. METHODS: We have developed a chimeric 137-mer synthetic polypeptide (PvCS-NRC) that includes the conserved region I and region II-plus and the two natural repeat variants known as VK210 and VK247. The antigenicity of PvCS-NRC was tested using human sera from PNG and Colombia endemic areas and its immunogenicity was confirmed in mice with different genetic backgrounds, the polypeptide formulated either in Alum or GLA-SE adjuvants was assessed in inbred C3H, CB6F1 and outbred ICR mice, whereas a formulation in Montanide ISA51 was tested in C3H mice. RESULTS: Antigenicity studies indicated that the chimeric peptide is recognized by a high proportion (60-70%) of residents of malaria-endemic areas. Peptides formulated with either GLA-SE or Montanide ISA51 adjuvants induced stronger antibody responses as compared with the Alum formulation. Sera from immunized mice as well as antigen-specific affinity purified human IgG antibodies reacted with sporozoite preparations in immunofluorescence and Western blot assays, and displayed strong in vitro inhibition of sporozoite invasion (ISI) into hepatoma cells. CONCLUSIONS: The polypeptide was recognized at high prevalence when tested against naturally induced human antibodies and was able to induce significant immunogenicity in mice. Additionally, specific antibodies were able to recognize sporozoites and were able to block sporozoite invasion in vitro. Further evaluation of this chimeric protein construct in preclinical phase e.g. in Aotus monkeys in order to assess the humoral and cellular immune responses as well as protective efficacy against parasite challenge of the vaccine candidate must be conducted.

Identification of infectious agents in onychomycoses by PCR-terminal restriction fragment length polymorphism.

A fast and reliable assay for the identification of dermatophyte fungi and nondermatophyte fungi (NDF) in onychomycosis is essential, since NDF are especially difficult to cure using standard treatment. Diagnosis is usually based on both direct microscopic examination of nail scrapings and macroscopic and microscopic identification of the infectious fungus in culture assays. In the last decade, PCR assays have been developed for the direct detection of fungi in nail samples. In this study, we describe a PCR-terminal restriction fragment length polymorphism (TRFLP) assay to directly and routinely identify the infecting fungi in nails. Fungal DNA was easily extracted using a commercial kit after dissolving nail fragments in an Na(2)S solution. Trichophyton spp., as well as 12 NDF, could be unambiguously identified by the specific restriction fragment size of 5'-end-labeled amplified 28S DNA. This assay enables the distinction of different fungal infectious agents and their identification in mixed infections. Infectious agents could be identified in 74% (162/219) of cases in which the culture results were negative. The PCR-TRFLP assay described here is simple and reliable. Furthermore, it has the possibility to be automated and thus routinely applied to the rapid diagnosis of a large number of clinical specimens in dermatology laboratories.

Global Genomic Diversity of Human Papillomavirus 6 Based on 724 Isolates and 190 Complete Genome Sequences.

Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts and laryngeal papillomas and has been included in both the quadrivalent and nonavalent prophylactic HPV vaccines. This study investigated the global genomic diversity of HPV6, using 724 isolates and 190 complete genomes from six continents, and the association of HPV6 genomic variants with geographical location, anatomical site of infection/disease, and gender. Initially, a 2,800-bp E5a-E5b-L1-LCR fragment was sequenced from 492/530 (92.8%) HPV6-positive samples collected for this study. Among them, 130 exhibited at least one single nucleotide polymorphism (SNP), indel, or amino acid change in the E5a-E5b-L1-LCR fragment and were sequenced in full. A global alignment and maximum likelihood tree of 190 complete HPV6 genomes (130 fully sequenced in this study and 60 obtained from sequence repositories) revealed two variant lineages, A and B, and five B sublineages: B1, B2, B3, B4, and B5. HPV6 (sub)lineage-specific SNPs and a 960-bp representative region for whole-genome-based phylogenetic clustering within the L2 open reading frame were identified. Multivariate logistic regression analysis revealed that lineage B predominated globally. Sublineage B3 was more common in Africa and North and South America, and lineage A was more common in Asia. Sublineages B1 and B3 were associated with anogenital infections, indicating a potential lesion-specific predilection of some HPV6 sublineages. Females had higher odds for infection with sublineage B3 than males. In conclusion, a global HPV6 phylogenetic analysis revealed the existence of two variant lineages and five sublineages, showing some degree of ethnogeographic, gender, and/or disease predilection in their distribution. IMPORTANCE: This study established the largest database of globally circulating HPV6 genomic variants and contributed a total of 130 new, complete HPV6 genome sequences to available sequence repositories. Two HPV6 variant lineages and five sublineages were identified and showed some degree of association with geographical location, anatomical site of infection/disease, and/or gender. We additionally identified several HPV6 lineage- and sublineage-specific SNPs to facilitate the identification of HPV6 variants and determined a representative region within the L2 gene that is suitable for HPV6 whole-genome-based phylogenetic analysis. This study complements and significantly expands the current knowledge of HPV6 genetic diversity and forms a comprehensive basis for future epidemiological, evolutionary, functional, pathogenicity, vaccination, and molecular assay development studies.

Identification and quantification techniques in mass spectrometry-based proteomics

Résumé La protéomique basée sur la spectrométrie de masse est l'étude du proteome l'ensemble des protéines exprimées au sein d'une cellule, d'un tissu ou d'un organisme - par cette technique. Les protéines sont coupées à l'aide d'enzymes en plus petits morceaux -les peptides -, et, séparées par différentes techniques. Les différentes fractions contenant quelques centaines de peptides sont ensuite analysées dans un spectromètre de masse. La masse des peptides est enregistrée et chaque peptide est séquentiellement fragmenté pour en obtenir sa séquence. L'information de masse et séquence est ensuite comparée à une base de données de protéines afin d'identifier la protéine d'origine. Dans une première partie, la thèse décrit le développement de méthodes d'identification. Elle montre l'importance de l'enrichissement de protéines comme moyen d'accès à des protéines de moyenne à faible abondance dans le lait humain. Elle utilise des injections répétées pour augmenter la couverture en protéines et la confiance dans l'identification. L'impacte de nouvelle version de base de données sur la liste des protéines identifiées est aussi démontré. De plus, elle utilise avec succès la spectrométrie de masse comme alternative aux anticorps, pour valider la présence de 34 constructions de protéines pathogéniques du staphylocoque doré exprimées dans une souche de lactocoque. Dans une deuxième partie, la thèse décrit le développement de méthodes de quantification. Elle expose de nouvelles approches de marquage des terminus des protéines aux isotopes stables et décrit la première méthode de marquage des groupements carboxyliques au niveau protéine à l'aide de réactifs composé de carbone 13. De plus, une nouvelle méthode, appelée ANIBAL, marquant tous les groupements amines et carboxyliques au niveau de la protéine, est exposée. Summary Mass spectrometry-based proteomics is the study of the proteome -the set of all expressed proteins in a cell, tissue or organism -using mass spectrometry. Proteins are cut into smaller pieces - peptides - using proteolytic enzymes and separated using different separation techniques. The different fractions containing several hundreds of peptides are than analyzed by mass spectrometry. The mass of the peptides entering the instrument are recorded and each peptide is sequentially fragmented to obtain its amino acid sequence. Each peptide sequence with its corresponding mass is then searched against a protein database to identify the protein to which it belongs. This thesis presents new method developments in this field. In a first part, the thesis describes development of identification methods. It shows the importance of protein enrichment methods to gain access to medium-to-low abundant proteins in a human milk sample. It uses repeated injection to increase protein coverage and confidence in identification and demonstrates the impact of new database releases on protein identification lists. In addition, it successfully uses mass spectrometry as an alternative to antibody-based assays to validate the presence of 34 different recombinant constructs of Staphylococcus aureus pathogenic proteins expressed in a Lactococcus lactis strain. In a second part, development of quantification methods is described. It shows new stable isotope labeling approaches based on N- and C-terminus labeling of proteins and describes the first method of labeling of carboxylic groups at the protein level using 13C stable isotopes. In addition, a new quantitative approach called ANIBAL is explained that labels all amino and carboxylic groups at the protein level.

Molecular identification of an endemic Alpine mammal, Apodemus alpicola, using a PCR-based RFLP method

The ability of a PCR-based restriction fragment length polymorphism (RFLP) analysis of the cytochrome b (mtDNA) to distinguish Apodemus alpicola from two other Apodemus species was investigated. The partial sequencing of the cytochrome b allowed the identification of one enzyme as being potentially diagnostic. This was supported by an analysis of 131 specimens previously identified using morphometric and/or allozymic data, indicating that the PCR-based RFLP method provides a rapid and reliable tool for distinguishing A. alpicola from its two co-occurring congenerics. The method is applicable to samples taken in the field for ecological studies, and could easily be adapted to the identification of museum samples.

Practical considerations for improving the productivity of mass spectrometry-based proteomics.

Mass spectrometry (MS) is currently the most sensitive and selective analytical technique for routine peptide and protein structure analysis. Top-down proteomics is based on tandem mass spectrometry (MS/ MS) of intact proteins, where multiply charged precursor ions are fragmented in the gas phase, typically by electron transfer or electron capture dissociation, to yield sequence-specific fragment ions. This approach is primarily used for the study of protein isoforms, including localization of post-translational modifications and identification of splice variants. Bottom-up proteomics is utilized for routine high-throughput protein identification and quantitation from complex biological samples. The proteins are first enzymatically digested into small (usually less than ca. 3 kDa) peptides, these are identified by MS or MS/MS, usually employing collisional activation techniques. To overcome the limitations of these approaches while combining their benefits, middle-down proteomics has recently emerged. Here, the proteins are digested into long (3-15 kDa) peptides via restricted proteolysis followed by the MS/MS analysis of the obtained digest. With advancements of high-resolution MS and allied techniques, routine implementation of the middle-down approach has been made possible. Herein, we present the liquid chromatography (LC)-MS/MS-based experimental design of our middle-down proteomic workflow coupled with post-LC supercharging.

Evaluation of a polymerase chain reaction-restriction fragment length polymorphism assay for dermatophyte and nondermatophyte identification in onychomycosis.

BACKGROUND: Dermatophytes are the main cause of onychomycoses, but various nondermatophyte filamentous fungi are often isolated from abnormal nails. The correct identification of the aetiological agent of nail infections is necessary in order to recommend appropriate treatment. OBJECTIVE: To evaluate a rapid polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay based on 28S rDNA for fungal identification in nails on a large number of samples in comparison with cultures. METHODS: Infectious fungi were analysed using PCR-RFLP in 410 nail samples in which fungal elements were observed in situ by direct mycological examination (positive samples). The results were compared with those previously obtained by culture of fungi on Sabouraud agar from the same nail samples. RESULTS: PCR-RFLP identification of fungi in nails allowed validation of the results obtained in culture when Trichophyton spp. grew from infected samples. In addition, nondermatophyte filamentous fungi could be identified with certainty as the infectious agents in onychomycosis, and discriminated from dermatophytes as well as from transient contaminants. The specificity of the culture results relative to PCR-RFLP appeared to be 81%, 71%, 52% and 63% when Fusarium spp., Scopulariopsis brevicaulis, Aspergillus spp. and Candida spp., respectively, grew on Sabouraud agar. It was also possible to identify the infectious agent when direct nail mycological examination showed fungal elements, but negative results were obtained from fungal culture. CONCLUSIONS: Improved sensitivity for the detection of fungi in nails was obtained using the PCR-RFLP assay. Rapid and reliable molecular identification of the infectious fungus can be used routinely and presents several important advantages compared with culture in expediting the choice of appropriate antifungal therapy.

4B.05: Plasma Lasma copeptin is associated with insulin resistance in a Swiss population-based study

OBJECTIVE: Previous studies suggest that arginine vasopressin may have a role in metabolic syndrome (MetS) and diabetes by altering liver glycogenolysis, insulin, and glucagon secretion and pituitary ACTH release. We tested whether plasma copeptin, the stable C-terminal fragment of arginine vasopressin prohormone, was associated with insulin resistance and MetS in a Swiss population-based study. DESIGN AND METHOD: We analyzed data from the population-based Swiss Kidney Project on Genes in Hypertension. Copeptin was assessed by an immunoluminometric assay. Insulin resistance was derived from the HOMA model and calculated as follows: (FPI x FPG)/22.5, where FPI is fasting plasma insulin concentration (mU/L) and FPG fasting plasma glucose (mmol/L). Subjects were classified as having the MetS according to the National Cholesterol Education Program Adult Treatment Panel III criteria. Mixed multivariate linear regression models were built to explore the association of insulin resistance with copeptin. In addition, multivariate logistic regression models were built to explore the association between MetS and copeptin. In the two analyses, adjustment was done for age, gender, center, tobacco and alcohol consumption, socioeconomic status, physical activity, intake of fruits and vegetables and 24 h urine flow rate. Copeptin was log-transformed for the analyses. RESULTS: Among the 1,089 subjects included in this analysis, 47% were male. Mean (SD) age and body mass index were 47.4 (17.6) years 25.0 (4.5) kg/m2. The prevalence of MetS was 10.5%. HOMA-IR was higher in men (median 1.3, IQR 0.7-2.1) than in women (median 1.0, IQR 0.5-1.6,P < 0.0001). Plasma copeptin was higher in men (median 5.2, IQR 3.7-7.8 pmol/L) than in women (median 3.0, IQR 2.2-4.3 pmol/L), P < 0.0001. HOMA-IR was positively associated with log-copeptin after full adjustment (β (95% CI) 0.19 (0.09-0.29), P < 0.001). MetS was not associated with copeptin after full adjustment (P = 0.92). CONCLUSIONS: Insulin resistance, but not MetS, was associated with higher copeptin levels. Further studies should examine whether modifying pharmacologically the arginine vasopressin system might improve insulin resistance, thereby providing insight into the causal nature of this association.