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.

Molecular Magnetism: The Design, Synthesis, and Characterization of New Building Blocks for Molecule-based Magnetic Materials

Two classes of building blocks have been prepared and characterized and their coordination chemistry explored working towards the preparation of new molecule-based magnetic materials. In the first project, the amine functionality of 3,3'-diamino-2,2'- bipyridine was exploited for the preparation of a new family of ligands (H2L 1)-(H2L 4). The molecular structures of three ligands have been fully characterized by X-ray crystallography. [molecular structure diagram will not copy here, but is available in full pdf.] The coordination chemistry of these ligands with divalent first row transition metal ions was investigated. For ligand (H2L1), the molecular structures of four coordination complexes with stoichiometries [Zn2(Ll)(OAc)(MeO)]2 (I), [Cu2(L1)(OAc)2 (II), [Li(L1)]3 (III), and [Ni(L1)]3 (IV) were determined by X-ray crystallography. For ligand (H2L2), a Cu(II) complex of stoichiometry [Cu3(L2)(OAc)3MeO] (V) was determined by X-ray crystallography. The magnetic properties of complexes (II), (III), and (V) have been fully elucidated. In project two, synthetic strategies for the preparation of porphyrin molecules bearing triol substituents is presented. Following this approach, three new porphyrin derivatives have been prepared and characterized [Zn(HPTPP-CH2C(CH20H)3)] (VI), [P(TPP)(OCH2C(CH2)H)3)2]+CL- (VII), and [P(OEP)(C6H5)(OCH2C(CH2OH)3)]+Cl- (VIII). Attempts to exchange the labile methoxide bridges of a tetraironIIl single molecule magnet of stoichiometry [Fe4(OMe)6(dpm)6] (Hdpm = dipivaloylmethane) with the triol appended porphyrins will be discussed. [molecular structure diagram will not copy here, but is available in full pdf.]

What have dendritic cells ever done for adjuvant design? Cellular and molecular methods for the rational development of vaccine adjuvants

Our new molecular understanding of immune priming states that dendritic cell activation is absolutely pivotal for expansion and differentiation of naïve T lymphocytes, and it follows that understanding DC activation is essential to understand and design vaccine adjuvants. This chapter describes how dendritic cells can be used as a core tool to provide detailed quantitative and predictive immunomics information about how adjuvants function. The role of distinct antigen, costimulation, and differentiation signals from activated DC in priming is explained. Four categories of input signals which control DC activation – direct pathogen detection, sensing of injury or cell death, indirect activation via endogenous proinflammatory mediators, and feedback from activated T cells – are compared and contrasted. Practical methods for studying adjuvants using DC are summarised and the importance of DC subset choice, simulating T cell feedback, and use of knockout cells is highlighted. Finally, five case studies are examined that illustrate the benefit of DC activation analysis for understanding vaccine adjuvant function.

Insights into dietary flavonoids as molecular templates for the design of anti-platelet drugs

Flavonoids are low-molecular weight, aromatic compounds derived from fruits, vegetables, and other plant components. The consumption of these phytochemicals has been reported to be associated with reduced cardiovascular disease (CVD) risk, attributed to their anti-inflammatory, anti-proliferative, and anti-thrombotic actions. Flavonoids exert these effects by a number of mechanisms which include attenuation of kinase activity mediated at the cell-receptor level and/or within cells, and are characterized as broad-spectrum kinase inhibitors. Therefore, flavonoid therapy for CVD is potentially complex; the use of these compounds as molecular templates for the design of selective and potent small-molecule inhibitors may be a simpler approach to treat this condition. Flavonoids as templates for drug design are, however, poorly exploited despite the development of analogues based on the flavonol, isoflavonone, and isoflavanone subgroups. Further exploitation of this family of compounds is warranted due to a structural diversity that presents great scope for creating novel kinase inhibitors. The use of computational methodologies to define the flavonoid pharmacophore together with biological investigations of their effects on kinase activity, in appropriate cellular systems, is the current approach to characterize key structural features that will inform drug design. This focussed review highlights the potential of flavonoids to guide the design of clinically safer, more selective, and potent small-molecule inhibitors of cell signalling, applicable to anti-platelet therapy.

Investigating flavonoids as molecular templates for the design of small-molecule inhibitors of cell signaling

Epidemiological and clinical trials reveal compelling evidence for the ability of dietary flavonoids to lower cardiovascular disease risk. The mechanisms of action of these polyphenolic compounds are diverse, and of particular interest is their ability to function as protein and lipid kinase inhibitors. We have previously described structure-activity studies that reinforce the possibility for using flavonoid structures as templates for drug design. In the present study, we aim to begin constructing rational screening strategies for exploiting these compounds as templates for the design of clinically relevant, antiplatelet agents. We used the platelet as a model system to dissect the structural influence of flavonoids, stilbenes, anthocyanidins, and phenolic acids on inhibition of cell signaling and function. Functional groups identified as relevant for potent inhibition of platelet function included at least 2 benzene rings, a hydroxylated B ring, a planar C ring, a C ring ketone group, and a C-2 positioned B ring. Hydroxylation of the B ring with either a catechol group or a single C-4' hydroxyl may be required for efficient inhibition of collagen-stimulated tyrosine phosphorylated proteins of 125 to 130 kDa, but may not be necessary for that of phosphotyrosine proteins at approximately 29 kDa. The removal of the C ring C-3 hydroxyl together with a hydroxylated B ring (apigenin) may confer selectivity for 37 to 38 kDa phosphotyrosine proteins. We conclude that this study may form the basis for construction of maps of flavonoid inhibitory activity on kinase targets that may allow a multitargeted therapeutic approach with analogue counterparts and parent compounds.

Molecular assay optimized by Taguchi experimental design method for venous thromboembolism investigation

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Molecular hybridization: A useful tool in the design of new drug prototypes

Molecular hybridization is a new concept in drug design and development based on the combination of pharmacophoric moieties of different bioactive substances to produce a new hyrid compound with improved affinity and efficacy, when compared to the parent drugs. Additionally, this strategy can results in compounds presenting modified selectivity profile, different and/or dual modes of action and reduced undesired side effects. So, in this described several example of different strategies for drug design, discovery and pharmacomodulation focused on new innovative hybrid compounds presenting analgesic, anti-inflammatory, platelet anti-aggregating, anti-infections, anticancer, cardio- and neuroactive properties.

Nubbe natural products, source of molecular diversity for the design of new anticancer agents

Pós-graduação em Química - IQ

Anti-inflammatory drug design by molecular hybridization approach

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Molecular basis oh herpes simplex virus entry into the cell and retargeting of the viral tropism for the design of oncolytic herpesviruses

Herpes simplex virus 1 (HSV-1) infects oral epitelial cells, then spreads to the nerve endings and estabilishes latency in sensory ganglia, from where it may, or may not reactivate. Diseases caused by virus reactivation include mild diseases such as muco-cutaneous lesions, and more severe, and even life-threatening encephalitis, or systemic infections affecting diverse organs. Herpes simplex virus represents the most comprehensive example of virus receptor interaction in Herpesviridae family, and the prototype virus encoding multipartite entry genes. In fact, it encodes 11-12 glycoproteins and a number of additional membrane proteins: five of these proteins play key roles in virus entry into subsceptible cells. Thus, glycoprotein B (gB) and glycoprotein C (gC) interact with heparan sulfate proteoglycan to enable initial attachment to cell surfaces. In the next step, in the entry cascade, gD binds a specific surface receptor such as nectin1 or HVEM. The interaction of glycoprotein D with the receptor alters the conformation of gD to enable the activation of gB, glycoprotein H, and glycoprotein L, a trio of glycoproteins that execute the fusion of the viral envelope with the plasma membrane. In this thesis, I described two distinct projects: I. The retargeting of viral tropism for the design of oncolytic Herpesviruses: • capable of infecting cells through the human epitelial growth factor receptor 2 (HER2), overexpressed in highly malignant mammary and ovarian tumors and correlates with a poor prognosis; • detargeted from its natural receptors, HVEM and nectin1. To this end, we inserted a ligand to HER2 in gD. Because HER2 has no natural ligand, the selected ligand was a single chain antibody (scFv) derived from MAb4D5 (monoclonal antibody to HER2), herein designated scHER2. All recombinant viruses were targeted to HER2 receptor, but only two viruses (R-LM113 and R-LM249) were completely detargeted from HVEM and nectin1. To engineer R-LM113, we removed a large portion at the N-terminus of gD (from aa 6 to aa 38) and inserted scHER2 sequence plus 9-aa serine-glycine flexible linker at position 39. On the other hand, to engineer R-LM249, we replaced the Ig-folded core of gD (from aa 61 to aa 218) with scHER2 flanked by Ser-Gly linkers. In summary, these results provide evidence that: i. gD can tolerate an insert almost as big as gD itself; ii. the Ig-like domain of gD can be removed; iii. the large portion at the N-terminus of gD (from aa 6 to aa 38) can be removed without loss of key function; iv. R-LM113 and R-LM249 recombinants are ready to be assayed in animal models of mammary and ovary tumour. This finding and the avaibility of a large number of scFv greatly increase the collection of potential receptors to which HSV can be redirected. II. The production and purification of recombinant truncated form of the heterodimer gHgL. We cloned a stable insect cell line expressing a soluble form of gH in complex with gL under the control of a metalloprotein inducible promoter and purified the heterodimer by means of ONE-STrEP-tag system by IBA. With respect to biological function, the purified heterodimer is capable: • of reacting to antibodies that recognize conformation dependent epitopes and neutralize virion infectivity; • of binding a variety cells at cell surface. No doubt, the availability of biological active purified gHgL heterodimer, in sufficient quantities, will speed up the efforts to solve its crystal structure and makes it feasible to identify more clearly whether gHgL has a cellular partner, and what is the role of this interaction on virus entry.

Design, realization and characterization of automated millifluidic bioreactors for investigating the molecular evolution of lytic or lysogenic vector phages infecting engineered host e. Coli

In questa tesi viene presentato un bioreattore in grado di mantenere nel tempo condizioni biologiche tali che consentano di massimizzare i cicli di evoluzione molecolare di vettori di clonazione fagici: litico (T7) o lisogeno (M13). Verranno quindi introdtti concetti legati alla Teoria della Quasispecie e alla relazione tra errori di autoreplicazione e pressioni selettive naturali o artificiali su popolazioni di virus: il modello naturale del sistema evolutivo. Tuttavia, mantenere delle popolazioni di virus significa formire loro un substrato dove replicare. Per fare ciò, altri gruppi di ricerca hanno giá sviluppato complessi e costosi prototipi di macchinari per la crescita continua di popolazioni batteriche: i compartimenti dei sistemi evolutivi. Il bioreattore, oggetto di questo lavoro, fa parte del progetto europeo Evoprog: general purpose programmable machine evolution on a chip (Jaramillo’s Lab, University of Warwick) che, utilizzando tecnologie fagiche e regolazioni sintetiche esistenti, sará in grado di produrre funzionalità biocomputazionali di due ordini di grandezza più veloci rispetto alle tecniche convenzionali, riducendo allo stesso tempo i costi complessivi. Il primo prototipo consiste in uno o piú fermentatori, dove viene fatta crescere la cultura batterica in condizioni ottimizzate di coltivazione continua, e in un cellstat, un volume separato, dove avviene solo la replicazione dei virus. Entrambi i volumi sono di pochi millilitri e appropriatamente interconnessi per consentire una sorta di screening continuo delle biomolecole prodotte all’uscita. Nella parte finale verranno presentati i risultati degli esperimenti preliminari, a dimostrazione dell’affidabilità del prototipo costruito e dei protocolli seguiti per la sterilizzazione e l’assemblaggio del bioreattore. Gli esperimenti effettuati dimostrano il successo di due coltivazioni virali continue e una ricombinazione in vivo di batteriofagi litici o lisogeni ingegnerizzati. La tesi si conclude valutando i futuri sviluppi e i limiti del sistema, tenendo in considerazione, in particolare, alcune applicazioni rivolte agli studi di una terapia batteriofagica.

The Design and Synthesis of Versatile Molecular Building Blocks: Working Towards the Controlled Self-Assembly of Novel Functional Materials

Our research goals are focused on the preparation of novel molecule-based materials that possess specifically designed properties in solution or in the solid state e.g. self-assembly, magnetism, conductivity and spin crossover phenomena. Most of our systems incorporate paramagnetic transition metal ions and the search for new molecule-based magnetic materials is a prominent theme. Specific areas of research include the preparation and study of oxalate based 2D and 3D magnets, probing the versatility of octacyanometalate building blocks as precursors for new molecular magnets, and the preparation of new tetrathiafulvalene (TIF) derivatives for applications in molecular and supramolecular chemistry.

Design, synthesis and cellular and molecular pharmacology of 2$\sp\prime$-bromo-4$\sp\prime$-epidaunorubicin (WP401): A novel non-cross-resistant anthracycline antibiotic with the intrinsic ability to circumvent P -glycoprotein-mediated drug resistance

We designed and synthesized a novel daunorubicin (DNR) analogue that effectively circumvents P-glycoprotein (P-gp)-mediated drug resistance. The fully protected carbohydrate intermediate 1,2-dibromoacosamine was prepared from acosamine and effectively coupled to daunomycinone in high yield. Deprotection under alkaline conditions yielded 2$\sp\prime$-bromo-4$\sp\prime$-epidaunorubicin (WP401). The in vitro cytotoxicity and cellular and molecular pharmacology of WP401 were compared with those of DNR in a panel of wild-type cell lines (KB-3-1, P388S, and HL60S) and their multidrug-resistant (MDR) counterparts (KB-V1, P388/DOX, and HL60/DOX). Fluorescent spectrophotometry, flow cytometry, and confocal laser scanning microscopy were used to measure intracellular accumulation, retention, and subcellular distribution of these agents. All MDR cell lines exhibited reduced DNR uptake that was restored, upon incubation with either verapamil (VER) or cyclosporin A (CSA), to the level found in sensitive cell lines. In contrast, the uptake of WP401 was essentially the same in the absence or presence of VER or CSA in all tested cell lines. The in vitro cytotoxicity of WP401 was similar to that of DNR in the sensitive cell lines but significantly higher in resistant cell lines (resistance index (RI) of 2-6 for WP401 vs 75-85 for DNR). To ascertain whether drug-mediated cytotoxicity and retention were accompanied by DNA strand breaks, DNA single- and double-strand breaks were assessed by alkaline elution. High levels of such breaks were obtained using 0.1-2 $\mu$g/mL of WP401 in both sensitive and resistant cells. In contrast, DNR caused strand breaks only in sensitive cells and not much in resistant cells. We also compared drug-induced DNA fragmentation similar to that induced by DNR. However, in P-gp-positive cells, WP401 induced 2- to 5-fold more DNA fragmentation than DNR. This increased DNA strand breakage by WP401 was correlated with its increased uptake and cytotoxicity in these cell lines. Overall these results indicate that WP401 is more cytotoxic than DNR in MDR cells and that this phenomenon might be related to the reduced basicity of the amino group and increased lipophilicity of WP401. ^