Syntheses, receptor bindings, in vitro and in vivo stabilities and biodistributions of DOTA-neurotensin(8-13) derivatives containing β-amino acid residues - a lesson about the importance of animal experiments

Neurotensin(8-13) (NTS(8-13)) analogs with C- and/or N-terminal β-amino acid residues and three DOTA derivatives thereof have been synthesized (i.e., 1-6). A virtual docking experiment showed almost perfect fit of one of the 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (DOTA) derivatives, 6a, into a crystallographically identified receptor NTSR1 (Fig.1). The affinities for the receptors of the NTS analogs and derivatives are low, when determined with cell-membrane homogenates, while, with NTSR1-exhibiting cancer tissues, affinities in the single-digit nanomolar range can be observed (Table 2). Most of the β-amino acid-containing NTS(8-13) analogs (Table 1 and Fig.2), including the (68) Ga complexes of the DOTA-substituted ones (6; Figs.2 and 5), are stable for ca. 1 h in human serum and plasma, and in murine plasma. The biodistributions of two (68) Ga complexes (of 6a and 6b) in HT29 tumor-bearing nude mice, in the absence and in the presence of a blocking compound, after 10, 30, and 60 min (Figs. 3 and 4) lead to the conclusion that the amount of specifically bound radioligand is rather low. This was confirmed by PET-imaging experiments with the tumor-bearing mice (Fig.6). Comparison of the in vitro plasma stability (after 1 h) with the ex vivo blood content (after 10-15 min) of the two (68) Ga complexes shows that they are rapidly cleaved in the animals (Fig.5).

Preclinical modeling of multi-drug cancer therapies

Anticancer drugs typically are administered in the clinic in the form of mixtures, sometimes called combinations. Only in rare cases, however, are mixtures approved as drugs. Rather, research on mixtures tends to occur after single drugs have been approved. The goal of this research project was to develop modeling approaches that would encourage rational preclinical mixture design. To this end, a series of models were developed. First, several QSAR classification models were constructed to predict the cytotoxicity, oral clearance, and acute systemic toxicity of drugs. The QSAR models were applied to a set of over 115,000 natural compounds in order to identify promising ones for testing in mixtures. Second, an improved method was developed to assess synergistic, antagonistic, and additive effects between drugs in a mixture. This method, dubbed the MixLow method, is similar to the Median-Effect method, the de facto standard for assessing drug interactions. The primary difference between the two is that the MixLow method uses a nonlinear mixed-effects model to estimate parameters of concentration-effect curves, rather than an ordinary least squares procedure. Parameter estimators produced by the MixLow method were more precise than those produced by the Median-Effect Method, and coverage of Loewe index confidence intervals was superior. Third, a model was developed to predict drug interactions based on scores obtained from virtual docking experiments. This represents a novel approach for modeling drug mixtures and was more useful for the data modeled here than competing approaches. The model was applied to cytotoxicity data for 45 mixtures, each composed of up to 10 selected drugs. One drug, doxorubicin, was a standard chemotherapy agent and the others were well-known natural compounds including curcumin, EGCG, quercetin, and rhein. Predictions of synergism/antagonism were made for all possible fixed-ratio mixtures, cytotoxicities of the 10 best-scoring mixtures were tested, and drug interactions were assessed. Predicted and observed responses were highly correlated (r2 = 0.83). Results suggested that some mixtures allowed up to an 11-fold reduction of doxorubicin concentrations without sacrificing efficacy. Taken together, the models developed in this project present a general approach to rational design of mixtures during preclinical drug development. ^

CHARACTERIZATION OF SELECTIVE INHIBITION OF ADENYLYL CYCLASE ACTIVITY BY SMALL MOLECULE INHIBITOR NKY80

The nine membrane-bound isoforms of adenylyl cyclase (AC), via synthesis of the signaling molecule cyclic AMP (cAMP), are involved in many isoform specific physiological functions. Decreasing AC5 activity has been shown to have potential therapeutic benefit, including reduced stress on the heart, pain relief, and attenuation of morphine dependence and withdrawal behaviors. However, AC structure is well conserved, and there are currently no isoform selective AC inhibitors in clinical use. P-site inhibitors inhibit AC directly at the catalytic site, but with an uncompetitive or noncompetitive mechanism. Due to this mechanism and nanomolar potency in cell-free systems, attempts at ligand-based drug design of novel AC inhibitors frequently use P-site inhibitors as a starting template. One small molecule inhibitor designed through this process, NKY80, is described as an AC5 selective inhibitor with low micromolar potency in vitro. P-site inhibitors reveal important ligand binding “pockets” in the AC catalytic site, but specific interactions that give NKY80 selectivity are unclear. Identifying and characterizing unique interactions between NKY80 and AC isoforms would significantly aid the development of isoform selective AC inhibitors. I hypothesized that NKY80’s selective inhibition is conferred by AC isoform specific interactions with the compound within the catalytic site. A structure-based virtual screen of the AC catalytic site was used to identify novel small molecule AC inhibitors. Identified novel inhibitors are isoform selective, supporting the catalytic site as a region capable of more potent isoform selective inhibition. Although NKY80 is touted commercially as an AC5 selective inhibitor, its characterization suggests strong inhibition of both AC5 and the closely related AC6. NKY80 was also virtually docked to AC to determine how NKY80 binds to the catalytic site. My results show a difference between NKY80 binding and the conformation of classic P-site inhibitors. The selectivity and notable differences in NKY80 binding to the AC catalytic site suggest a catalytic subregion more flexible in AC5 and AC6 that can be targeted by selective small molecule inhibitors.

Muscle-Type Nicotinic Receptor Blockade by Diethylamine, the Hydrophilic Moiety of Lidocaine

Lidocaine bears in its structure both an aromatic ring and a terminal amine, which can be protonated at physiological pH, linked by an amide group. Since lidocaine causes multiple inhibitory actions on nicotinic acetylcholine receptors (nAChRs), this work was aimed to determine the inhibitory effects of diethylamine (DEA), a small molecule resembling the hydrophilic moiety of lidocaine, on Torpedo marmorata nAChRs microtransplanted to Xenopus oocytes. Similarly to lidocaine, DEA reversibly blocked acetylcholine-elicited currents (IACh) in a dose-dependent manner (IC50 close to 70 μM), but unlike lidocaine, DEA did not affect IACh desensitization. IACh inhibition by DEA was more pronounced at negative potentials, suggesting an open-channel blockade of nAChRs, although roughly 30% inhibition persisted at positive potentials, indicating additional binding sites outside the pore. DEA block of nAChRs in the resting state (closed channel) was confirmed by the enhanced IACh inhibition when pre-applying DEA before its co-application with ACh, as compared with solely DEA and ACh co-application. Virtual docking assays provide a plausible explanation to the experimental observations in terms of the involvement of different sets of drug binding sites. So, at the nAChR transmembrane (TM) domain, DEA and lidocaine shared binding sites within the channel pore, giving support to their open-channel blockade; besides, lidocaine, but not DEA, interacted with residues at cavities among the M1, M2, M3, and M4 segments of each subunit and also at intersubunit crevices. At the extracellular (EC) domain, DEA and lidocaine binding sites were broadly distributed, which aids to explain the closed channel blockade observed. Interestingly, some DEA clusters were located at the α-γ interphase of the EC domain, in a cavity near the orthosteric binding site pocket; by contrast, lidocaine contacted with all α-subunit loops conforming the ACh binding site, both in α-γ and α-δ and interphases, likely because of its larger size. Together, these results indicate that DEA mimics some, but not all, inhibitory actions of lidocaine on nAChRs and that even this small polar molecule acts by different mechanisms on this receptor. The presented results contribute to a better understanding of the structural determinants of nAChR modulation.

Use of virtual screening, flexible docking, and molecular interaction fields to design novel HMG-CoA reductase inhibitors for the treatment of hypercholesterolemia

Dietary changes associated with drug therapy can reduce high serum cholesterol levels and dramatically decrease the risk of coronary artery disease, stroke, and overall mortality. Statins are hypolipemic drugs that are effective in the reduction of cholesterol serum levels, attenuating cholesterol synthesis in liver by competitive inhibition regarding the substrate or molecular target HMG-CoA reductase. We have herewith used computer-aided molecular design tools, i.e., flexible docking, virtual screening in large data bases, molecular interaction fields to propose novel potential HMG-CoA reductase inhibitors that are promising for the treatment of hypercholesterolemia.

Molecular dynamics, flexible docking, virtual screening, ADMET predictions, and molecular interaction field studies to design novel potential MAO-B inhibitors

Monoamine oxidase is a flavoenzyme bound to the mitochondrial outer membranes of the cells, which is responsible for the oxidative deamination of neurotransmitter and dietary amines. It has two distinct isozymic forms, designated MAO-A and MAO-B, each displaying different substrate and inhibitor specificities. They are the well-known targets for antidepressant, Parkinson`s disease, and neuroprotective drugs. Elucidation of the x-ray crystallographic structure of MAO-B has opened the way for the molecular modeling studies. In this work we have used molecular modeling, density functional theory with correlation, virtual screening, flexible docking, molecular dynamics, ADMET predictions, and molecular interaction field studies in order to design new molecules with potential higher selectivity and enzymatic inhibitory activity over MAO-B.

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.

Computational strategies to include protein flexibility in Ligand Docking and Virtual Screening

The dynamic character of proteins strongly influences biomolecular recognition mechanisms. With the development of the main models of ligand recognition (lock-and-key, induced fit, conformational selection theories), the role of protein plasticity has become increasingly relevant. In particular, major structural changes concerning large deviations of protein backbones, and slight movements such as side chain rotations are now carefully considered in drug discovery and development. It is of great interest to identify multiple protein conformations as preliminary step in a screening campaign. Protein flexibility has been widely investigated, in terms of both local and global motions, in two diverse biological systems. On one side, Replica Exchange Molecular Dynamics has been exploited as enhanced sampling method to collect multiple conformations of Lactate Dehydrogenase A (LDHA), an emerging anticancer target. The aim of this project was the development of an Ensemble-based Virtual Screening protocol, in order to find novel potent inhibitors. On the other side, a preliminary study concerning the local flexibility of Opioid Receptors has been carried out through ALiBERO approach, an iterative method based on Elastic Network-Normal Mode Analysis and Monte Carlo sampling. Comparison of the Virtual Screening performances by using single or multiple conformations confirmed that the inclusion of protein flexibility in screening protocols has a positive effect on the probability to early recognize novel or known active compounds.

Low molecular weight compounds from mushrooms as potential Bcl-2 inhibitors: docking and virtual screening studies

Mushrooms have the ability to promote apoptosis in tumor cell lines, but the mechanism of action is not quite well understood. Inhibition of the interaction between Bcl-2 and pro-apoptotic proteins could be an important step that leads to apoptosis. Therefore, the discovery of compounds with the ability to inhibit Bcl-2 is an ongoing research topic in drug discovery. In this study, we started by analyzing Bcl-2 experimental structures that are currently available in Protein Data Bank database. After analysis of the more relevant Bcl-2 structures, 4 were finally selected. An analysis of the best docking methodology was then performed using a cross-docking and re-docking approach while testing 2 docking softwares: AutoDock 4 and AutoDock Vina. Autodock4 provided the best docking results and was selected to perform a virtual screening study applied to a dataset of 40 Low Molecular Weight (LMW) compounds present in mushrooms, using the selected Bcl-2 structures as target. Results suggest that steroid are the more promising family, among the analyzed compounds, and may have the ability to interact with Bcl-2 and this way promoting tumor apoptosis. The steroids that presented lowest estimated binding energy (ΔG) were: Ganodermanondiol, Cerevisterol, Ganoderic Acid X and Lucidenic Lactone; with estimated ΔG values between -8,45 and -8,23 Kcal/mol. A detailed analysis of the docked conformation of these 4 top ranked LMW compounds was also performed and illustrates a plausible interaction between the 4 top raked steroids and Bcl-2, thus substantiating the accuracy of the predicted docked poses. Therefore, tumoral apoptosis promoted by mushroom might be related to Bcl-2 inhibition mediated by steroid family of compounds.

Homology modeling and metabolism prediction of human carboxylesterase-2 using docking analyses by GriDock: a parallelized tool based on AutoDock 4.0.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted to developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas less has been done to predict the activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES2. The study involved first a homology modeling of the hCES2 protein based on the model of hCES1 since the two proteins share a high degree of homology (congruent with 73%). A set of 40 known substrates of hCES2 was taken from the literature; the ligands were docked in both their neutral and ionized forms using GriDock, a parallel tool based on the AutoDock4.0 engine which can perform efficient and easy virtual screening analyses of large molecular databases exploiting multi-core architectures. Useful statistical models (e.g., r (2) = 0.91 for substrates in their unprotonated state) were calculated by correlating experimental pK(m) values with distance between the carbon atom of the substrate's ester group and the hydroxy function of Ser228. Additional parameters in the equations accounted for hydrophobic and electrostatic interactions between substrates and contributing residues. The negatively charged residues in the hCES2 cavity explained the preference of the enzyme for neutral substrates and, more generally, suggested that ligands which interact too strongly by ionic bonds (e.g., ACE inhibitors) cannot be good CES2 substrates because they are trapped in the cavity in unproductive modes and behave as inhibitors. The effects of protonation on substrate recognition and the contrasting behavior of substrates and products were finally investigated by MD simulations of some CES2 complexes.

Use of the FACTS solvation model for protein-ligand docking calculations. Application to EADock.

Protein-ligand docking has made important progress during the last decade and has become a powerful tool for drug development, opening the way to virtual high throughput screening and in silico structure-based ligand design. Despite the flattering picture that has been drawn, recent publications have shown that the docking problem is far from being solved, and that more developments are still needed to achieve high successful prediction rates and accuracy. Introducing an accurate description of the solvation effect upon binding is thought to be essential to achieve this goal. In particular, EADock uses the Generalized Born Molecular Volume 2 (GBMV2) solvent model, which has been shown to reproduce accurately the desolvation energies calculated by solving the Poisson equation. Here, the implementation of the Fast Analytical Continuum Treatment of Solvation (FACTS) as an implicit solvation model in small molecules docking calculations has been assessed using the EADock docking program. Our results strongly support the use of FACTS for docking. The success rates of EADock/FACTS and EADock/GBMV2 are similar, i.e. around 75% for local docking and 65% for blind docking. However, these results come at a much lower computational cost: FACTS is 10 times faster than GBMV2 in calculating the total electrostatic energy, and allows a speed up of EADock by a factor of 4. This study also supports the EADock development strategy relying on the CHARMM package for energy calculations, which enables straightforward implementation and testing of the latest developments in the field of Molecular Modeling.

Virtual screening and computational optimization for the discovery of covalent prolyl oligopeptidase inhibitors with activity in human cells.

Our docking program, Fitted, implemented in our computational platform, Forecaster, has been modified to carry out automated virtual screening of covalent inhibitors. With this modified version of the program, virtual screening and further docking-based optimization of a selected hit led to the identification of potential covalent reversible inhibitors of prolyl oligopeptidase activity. After visual inspection, a virtual hit molecule together with four analogues were selected for synthesis and made in one-five chemical steps. Biological evaluations on recombinant POP and FAPα enzymes, cell extracts, and living cells demonstrated high potency and selectivity for POP over FAPα and DPPIV. Three compounds even exhibited high nanomolar inhibitory activities in intact living human cells and acceptable metabolic stability. This small set of molecules also demonstrated that covalent binding and/or geometrical constraints to the ligand/protein complex may lead to an increase in bioactivity.

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.

rDock: A Fast, Versatile and Open Source Program for Docking Ligands to Proteins and Nucleic Acids

Identification of chemical compounds with specific biological activities is an important step in both chemical biology and drug discovery. When the structure of the intended target is available, one approach is to use molecular docking programs to assess the chemical complementarity of small molecules with the target; such calculations provide a qualitative measure of affinity that can be used in virtual screening (VS) to rank order a list of compounds according to their potential to be active. rDock is a molecular docking program developed at Vernalis for high-throughput VS (HTVS) applications. Evolved from RiboDock, the program can be used against proteins and nucleic acids, is designed to be computationally very efficient and allows the user to incorporate additional constraints and information as a bias to guide docking. This article provides an overview of the program structure and features and compares rDock to two reference programs, AutoDock Vina (open source) and Schrodinger's Glide (commercial). In terms of computational speed for VS, rDock is faster than Vina and comparable to Glide. For binding mode prediction, rDock and Vina are superior to Glide. The VS performance of rDock is significantly better than Vina, but inferior to Glide for most systems unless pharmacophore constraints are used; in that case rDock and Glide are of equal performance. The program is released under the Lesser General Public License and is freely available for download, together with the manuals, example files and the complete test sets, at http://rdock.sourceforge.net/

GRID and docking analyses reveal a molecular basis for flavonoid inhibition of src-family kinase activity

Flavonoids reduce cardiovascular disease risk through anti-inflammatory, anti-coagulant and anti-platelet actions. One key flavonoid inhibitory mechanism is blocking kinase activity that drives these processes. Flavonoids attenuate activities of kinases including phosphoinositide-3-kinase (PI3K), Fyn, Lyn, Src, Syk, PKC, PIM1/2, ERK, JNK, and PKA. X-ray crystallographic analyses of kinase-flavonoid complexes show that flavonoid ring systems and their hydroxyl substitutions are important structural features for their binding to kinases. A clearer understanding of structural interactions of flavonoids with kinases is necessary to allow construction of more potent and selective counterparts. We examined flavonoid (quercetin, apigenin and catechin) interactions with Src-family kinases (Lyn, Fyn and Hck) applying the Sybyl docking algorithm and GRID. A homology model (Lyn) was used in our analyses to demonstrate that high quality predicted kinase structures are suitable for flavonoid computational studies. Our docking results revealed potential hydrogen bond contacts between flavonoid hydroxyls and kinase catalytic site residues. Identification of plausible contacts indicated that quercetin formed the most energetically stable interactions, apigenin lacked hydroxyl groups necessary for important contacts, and the non-planar structure of catechin could not support predicted hydrogen bonding patterns. GRID analysis using a hydroxyl functional group supported docking results. Based on these findings, we predicted that quercetin would inhibit activities of Src-family kinases with greater potency than apigenin and catechin. We validated this prediction using in vitro kinase assays. We conclude that our study can be used as a basis to construct virtual flavonoid interaction libraries to guide drug discovery using these compounds as molecular templates.