A perspective on the sharing of docking data.

Computational approaches are nowadays largely applied in drug discovery projects. Among these, molecular docking is the most used for hit identification against a drug target protein. However, many scientists in the field shed light on the lack of availability and reproducibility of the data obtained from such studies to the whole community. Consequently, sustaining and developing the efforts toward a large and fully transparent sharing of those data could be beneficial for all researchers in drug discovery. The purpose of this article is first to propose guidelines and recommendations on the appropriate way to conduct virtual screening experiments and second to depict the current state of sharing molecular docking data. In conclusion, we have explored and proposed several prospects to enhance data sharing from docking experiment that could be developed in the foreseeable future.

Design and biological evaluation of substituted 5,7-dihydro-6H-indolo[2,3-c]quinolin-6-one as novel selective Haspin inhibitors.

A library of substituted indolo[2,3-c]quinolone-6-ones was developed as simplified Lamellarin isosters. Synthesis was achieved from indole after a four-step pathway sequence involving iodination, a Suzuki-Miyaura cross-coupling reaction, and a reduction/lactamization sequence. The inhibitory activity of the 22 novel derivatives was assessed on Haspin kinase. Two of them possessed an IC50 of 1 and 2 nM with selectivity towards a panel of 10 other kinases including the parent kinases DYRK1A and CLK1. The most selective compound exerted additionally a very interesting cell effect on the osteosarcoma U-2 OS cell line.

Design, Synthesis and SAR in 2,4,7-Trisubstituted Pyrido[3,2-d]Pyrimidine Series as Novel PI3K/mTOR Inhibitors.

This work describes the synthesis, enzymatic activities on PI3K and mTOR, in silico docking and cellular activities of various uncommon 2,4,7 trisubstituted pyrido[3,2-d]pyrimidines. The series synthesized offers a chemical diversity in C-7 whereas C-2 (3-hydroxyphenyl) and C-4 groups (morpholine) remain unchanged, in order to provide a better understanding of the molecular determinants of PI3K selectivity or dual activity on PI3K and mTOR. Some C-7 substituents were shown to improve the efficiency on kinases compared to the 2,4-di-substituted pyrimidopyrimidine derivatives used as references. Six novel derivatives possess IC50 values on PI3Kα between 3 and 10 nM. The compounds with the best efficiencies on PI3K and mTOR induced micromolar cytotoxicity on cancer cell lines possessing an overactivated PI3K pathway.

Comparative Assessment of Protein Kinase Inhibitors in Public Databases and in PKIDB.

Since the first approval of a protein kinase inhibitor (PKI) by the Food and Drug Administration (FDA) in 2001, 55 new PKIs have reached the market, and many inhibitors are currently being evaluated in clinical trials. This is a clear indication that protein kinases still represent major drug targets for the pharmaceutical industry. In a previous work, we have introduced PKIDB, a publicly available database, gathering PKIs that have already been approved (Phase 4), as well as those currently in clinical trials (Phases 0 to 3). This database is updated frequently, and an analysis of the new data is presented here. In addition, we compared the set of PKIs present in PKIDB with the PKIs in early preclinical studies found in ChEMBL, the largest publicly available chemical database. For each dataset, the distribution of physicochemical descriptors related to drug-likeness is presented. From these results, updated guidelines to prioritize compounds for targeting protein kinases are proposed. The results of a principal component analysis (PCA) show that the PKIDB dataset is fully encompassed within all PKIs found in the public database. This observation is reinforced by a principal moments of inertia (PMI) analysis of all molecules. Interestingly, we notice that PKIs in clinical trials tend to explore new 3D chemical space. While a great majority of PKIs is located on the area of "flatland", we find few compounds exploring the 3D structural space. Finally, a scaffold diversity analysis of the two datasets, based on frequency counts was performed. The results give insight into the chemical space of PKIs, and can guide researchers to reach out new unexplored areas. PKIDB is freely accessible from the following website: http://www.icoa.fr/pkidb.

Fr-PPIChem: An Academic Compound Library Dedicated to Protein-Protein Interactions.

Protein-protein interactions (PPIs) mediate nearly every cellular process and represent attractive targets for modulating disease states but are challenging to target with small molecules. Despite this, several PPI inhibitors (iPPIs) have entered clinical trials, and a growing number of PPIs have become validated drug targets. However, high-throughput screening efforts still endure low hit rates mainly because of the use of unsuitable screening libraries. Here, we describe the collective effort of a French consortium to build, select, and store in plates a unique chemical library dedicated to the inhibition of PPIs. Using two independent predictive models and two updated databases of experimentally confirmed PPI inhibitors developed by members of the consortium, we built models based on different training sets, molecular descriptors, and machine learning methods. Independent statistical models were used to select putative PPI inhibitors from large commercial compound collections showing great complementarity. Medicinal chemistry filters were applied to remove undesirable structures from this set (such as PAINS, frequent hitters, and toxic compounds) and to improve drug likeness. The remaining compounds were subjected to a clustering procedure to reduce the final size of the library while maintaining its chemical diversity. In practice, the library showed a 46-fold activity rate enhancement when compared to a non-iPPI-enriched diversity library in high-throughput screening against the CD47-SIRPα PPI. The Fr-PPIChem library is plated in 384-well plates and will be distributed on demand to the scientific community as a powerful tool for discovering new chemical probes and early hits for the development of potential therapeutic drugs.

Thiopurine Derivative-Induced Fpg/Nei DNA Glycosylase Inhibition: Structural, Dynamic and Functional Insights.

DNA glycosylases are emerging as relevant pharmacological targets in inflammation, cancer and neurodegenerative diseases. Consequently, the search for inhibitors of these enzymes has become a very active research field. As a continuation of previous work that showed that 2-thioxanthine (2TX) is an irreversible inhibitor of zinc finger (ZnF)-containing Fpg/Nei DNA glycosylases, we designed and synthesized a mini-library of 2TX-derivatives (TXn) and evaluated their ability to inhibit Fpg/Nei enzymes. Among forty compounds, four TXn were better inhibitors than 2TX for Fpg. Unexpectedly, but very interestingly, two dithiolated derivatives more selectively and efficiently inhibit the zincless finger (ZnLF)-containing enzymes (human and mimivirus Neil1 DNA glycosylases hNeil1 and MvNei1, respectively). By combining chemistry, biochemistry, mass spectrometry, blind and flexible docking and X-ray structure analysis, we localized new TXn binding sites on Fpg/Nei enzymes. This endeavor allowed us to decipher at the atomic level the mode of action for the best TXn inhibitors on the ZnF-containing enzymes. We discovered an original inhibition mechanism for the ZnLF-containing Fpg/Nei DNA glycosylases by disulfide cyclic trimeric forms of dithiopurines. This work paves the way for the design and synthesis of a new structural class of inhibitors for selective pharmacological targeting of hNeil1 in cancer and neurodegenerative diseases.

VSPrep: A KNIME Workflow for the Preparation of Molecular Databases for Virtual Screening.

Drug discovery is a challenging and expensive field. Hence, novel in silico tools have been developed in early discovery stage to identify and prioritize novel molecules with suitable physicochemical properties. In many in silico drug design projects, molecular databases are screened by virtual screening tools to search for potential bioactive molecules. The preparation of the molecules is therefore a key step in the success of well-established techniques such as docking, similarity or pharmacophore searching. We review here the lists of several toolkits used in different steps during the cleaning of molecular databases, integrated within a KNIME workflow. During the first step of the automatic workflow, salts are removed, and mixtures are split to get one compound per entry. Then compounds with unwanted features are filtered. Duplicated entries are then deleted while considering stereochemistry. As a compromise between exhaustiveness and computational time, most distributed tautomers at physiological pH are computed. Additionally, various flags are applied to molecules by using either classical molecular descriptors, similarity search to known libraries or substructure search rules. Moreover, stereoisomers are enumerated depending on the unassigned chiral centers. Then, three-dimensional coordinates, and optionally conformers, are generated. This workflow has been already applied to several drug design projects and can be used for molecular database preparation upon request.

Cinchona-based zwitterionic stationary phases: Exploring retention and enantioseparation mechanisms in supercritical fluid chromatography with a fragmentation approach.

Chiralpak ZWIX(+) and ZWIX(-), are brush-type bonded-silica chiral stationary phases (CSPs), based on complex diastereomeric Cinchona alkaloids derivatives bearing both a positive and a negative charge. In the present study, we aimed to improve the understanding of retention and enantioseparation mechanisms of these CSPs employed in supercritical fluid chromatography (SFC). For this purpose, 9 other stationary phases were used as comparison systems: two of them are commercially available and bear only a positive charge (Chiralpak QN-AX and QD-AX) and the 7 others were designed purposely to be structurally similar to the parent ZWIX phases, but miss some portion of the complex ligand. First, cluster analysis was employed to identify similar and dissimilar behavior among the 11 stationary phases, where ionic interactions appeared to dominate the observed differences. Secondly, the stationary phases were characterized with linear solvation energy relationships (LSER) based on the SFC analysis of 161 achiral analytes and a modified version of the solvation parameter model to include ionic interactions. This served to compare the interaction capabilities for the 11 stationary phases and showed in particular the contribution of attractive and repulsive ionic interactions. Then the ZWIX phases were characterized for their enantioseparation capabilities with a set of 58 racemic probes. Discriminant analysis was applied to explore the molecular structural features that are useful to successful enantioseparation on the ZWIX phases. In particular, it appeared that the presence of positive charges in the analyte is causing increased retention but is not necessarily a favorable feature to enantiorecognition. On the opposite, the presence of negative charges in the analyte favors early elution and enantiorecognition. Finally, a smaller set of 30 pairs of enantiomers, selected by their structural diversity and different enantioseparation values on the ZWIX phases, were analyzed on all chiral phases to observe the contribution of each structural fragment of the chiral ligand on enantioselectivity. Molecular modelling of the ligands also helped in understanding the three-dimensional arrangement of each ligand, notably the intra-molecular hydrogen bonding or the possible contribution of ionic interactions. In the end, each structural element in the ZWIX phases appeared to be a significant contributor to successful enantioresolution, whether they contribute as direct interaction groups (ion-exchange functions) or as steric constraints to orientate the interacting groups towards the analytes.

PKIDB: A Curated, Annotated and Updated Database of Protein Kinase Inhibitors in Clinical Trials.

The number of protein kinase inhibitors (PKIs) approved worldwide continues to grow steadily, with 39 drugs approved in the period between 2001 and January 2018. PKIs on the market have been the subject of many reviews, and structure-property relationships specific to this class of drugs have been inferred. However, the large number of PKIs under development is often overlooked. In this paper, we present PKIDB (Protein Kinase Inhibitor Database), a monthly-updated database gathering approved PKIs as well as PKIs currently in clinical trials. The database compiles currently 180 inhibitors ranging from phase 0 to 4 clinical trials along with annotations extracted from seven public resources. The distribution and property ranges of standard physicochemical properties are presented. They can be used as filters to better prioritize compound selection for future screening campaigns. Interestingly, more than one-third of the kinase inhibitors violate at least one Lipinski's rule. A Principal Component Analysis (PCA) reveals that Type-II inhibitors are mapped to a distinct chemical space as compared to orally administrated drugs as well as to other types of kinase inhibitors. Using a Principal Moment of Inertia (PMI) analysis, we show that PKIs under development tend to explore new shape territories as compared to approved PKIs. In order to facilitate the analysis of the protein space, the kinome tree has been annotated with all protein kinases being targeted by PKIs. Finally, we analyzed the pipeline of the pharmaceutical companies having PKIs on the market or still under development. We hope that this work will assist researchers in the kinase field in identifying and designing the next generation of kinase inhibitors for still untargeted kinases. The PKIDB database is freely accessible from a website at http://www.icoa.fr/pkidb and can be easily browsed through a user-friendly spreadsheet-like interface.

VSPrep: A General KNIME Workflow for the Preparation of Molecules for Virtual Screening.

Over the past decades, virtual screening has proved itself to be a valuable asset to identify new bioactive compounds. The vast majority of commonly used techniques can be described in three steps: pre-processing the dataset i. e. small (ligands) and eventually larger (receptors) molecules, execute the method and finally analyse the results. Hence, the preparation of ligands is a critical step for success of commonly used virtual screening approaches such as protein-ligand docking, similarity or pharmacophore search. We present here a new workflow, VSPrep, for the pre-processing of small molecules; it is based on freely accessible tools for academics and is integrated within the KNIME platform. It can be used to perform several chemoinformatics tasks such as molecular database cleaning, tautomer and stereoisomer enumeration, focused library design and conformer generation. Additionally, graphical reports of the results are provided to the user as a convenient analysis tool.

Comprehensive Network Map of ADME-Tox Databases.

In the last decade, many statistical-based approaches have been developed to improve poor pharmacokinetics (PK) and to reduce toxicity of lead compounds, which are one of the main causes of high failure rate in drug development. Predictive QSAR models are not always very efficient due to the low number of available biological data and the differences in the experimental protocols. Fortunately, the number of available databases continues to grow every year. However, it remains a challenge to determine the source and the quality of the original data. The main goal is to identify the relevant databases required to generate the most robust predictive models. In this study, an interactive network of databases was proposed to easily find online data sources related to ADME-Tox parameters data. In this map, relevant information regarding scope of application, data availability and data redundancy can be obtained for each data source. To illustrate the usage of data mining from the network, a dataset on plasma protein binding is selected based on various sources such as DrugBank, PubChem and ChEMBL databases. A total of 2,606 unique molecules with experimental values of PPB were extracted and can constitute a consistent dataset for QSAR modeling.

Novel optimization of valmerins (tetrahydropyrido[1,2-a]isoindolones) as potent dual CDK5/GSK3 inhibitors.

An efficient synthetic strategy able to modulate the structure of the tetrahydropyridine isoindolone (Valmerin) skeleton was developed. A library of more than 30 novel final structures was generated. Biological activities on CDK5 and GSK3 as well as cellular effects on cancer cell lines were measured for each novel compound. Additionally to support the SAR, a docking study was performed. A potent GSK3/CDK5 dual inhibitor (37, IC50 CDK5/GSK3 35/7 nM) was obtained. Best antiproliferative effects were obtained on lung and prostate cell lines with IC50 = 20 nM.

Synthesis of dihydropyrimidine α,γ-diketobutanoic acid derivatives targeting HIV integrase.

The synthesis and antiviral evaluation of a series of dihydropyrimidinone and thiopyrimidine derivatives bearing aryl α,γ-diketobutanoic acid moiety are described using the Biginelli multicomponent reaction as key step. The most active among 20 synthesized novel compounds were 4c, 4d and 5b, which possess nanomolar HIV-1 integrase (IN) stand transfer (ST) inhibition activities. In order to understand their mode of interactions within the IN active site, we docked all the compounds into the previously reported X-ray crystal structure of IN. We observed that compounds 4c, 4d and 5b occupied an area close to the two catalytic Mg(2+) ions surrounded by their chelating triad (E221, D128 and D185), DC16, Y212 and the ß-diketo acid moiety of 4c, 4d and 5b chelating Mg(2+). As those compounds lack anti-HIV activities in cell, their prodrugs were synthetized. The prodrug 4c' exhibited an anti-HIV activity of 0.19 µM in primary human lymphocytes with some cytotoxicity. All together, these results indicate that the new analogs potentially interact within the catalytic site with highly conserved residues important for IN catalytic activity.

Advances in tetrahydropyrido[1,2-a]isoindolone (valmerins) series: Potent glycogen synthase kinase 3 and cyclin dependent kinase 5 inhibitors.

An efficient synthetic strategy was developed to modulate the structure of the tetrahydropyridine isoindolone (Valmerin) skeleton. A library of more than 30 novel final structures was generated. Biological activities on CDK5 and GSK3 as well as cellular effects on cancer cell lines were measured for each novel compound. Additionally docking studies were performed to support medicinal chemistry efforts. A strong GSK3/CDK5 dual inhibitor (38, IC50 GSK3/CDK5 32/84 nM) was obtained. A set of highly selective GSK3 inhibitors was synthesized by fine-tuning structural modifications (29 IC50 GSK3/CDK5 32/320 nM). Antiproliferative effects on cells were correlated with the in vitro kinase activities and the best effects were obtained with lung and colon cell lines.

Synthesis and optimization of an original V-shaped collection of 4-7-disubstituted pyrido[3,2-d]pyrimidines as CDK5 and DYRK1A inhibitors.

We here report the synthesis and biological evaluation of an original collection of 4,7-disubstituted pyrido[3,2-d]pyrimidines designed as potential kinase inhibitors. The collection was generated from a single starting material, 4,7-dichloropyrido[3,2-d]pyrimidine, which afforded the final compounds after two steps: a sequential or one-pot sequence including selective cross coupling reactions in C-4, followed by the second cross-coupling in C-7. In position C-4, a Suzuki-Miyaura type reaction led to monosubstituted derivatives whereas in position C-7, synthesis was achieved via a Suzuki or a Buchwald type reaction using commercially available or undescribed boron derivatives. The biological activity of the V-shaped family was measured in protein kinase assays. The structure activity relationship (SAR) revealed that some compounds selectively inhibited DYRK1A and CDK5 without affecting GSK3. Docking studies furnished possible explanations that correlate with the SAR data. The most active compound on the two biological targets was 27 which exhibited the following IC50: 110 nM for CDK5, 24 nM for DYRK1A and only 1.2 µM for GSK3. In the C-7 amino subfamily, the best derivative was indubitably compound 48 which led to a near selective action on DYRK1A and a remarkable IC50 of 60 nM.

Design, synthesis, and biological activity of pyridopyrimidine scaffolds as novel PI3K/mTOR dual inhibitors.

The design, synthesis, and screening of dual PI3K/mTOR inhibitors that gave nanomolar enzymatic and cellular activities on both targets with an acceptable kinase selectivity profile are described. A docking study was performed to understand the binding mode of the compounds and to explain the differences in biological activity. In addition, cellular effects of the best dual inhibitors were determined on six cancer cell lines and compared to those on a healthy diploid cell line for cellular cytotoxicity. Two compounds are highly potent on cancer cells in the submicromolar range without any toxicity on healthy cells. A more detailed analysis of the cellular effect of these PI3K/mTOR dual inhibitors demonstrated that they induce G1-phase cell cycle arrest in breast cancer cells and trigger apoptosis. These compounds show an interesting kinase profile as dual PI3K/mTOR tool compounds or as a chemical series for further optimization to progress into in vivo experiments.

Synthesis and biological evaluation of enantiomerically pure cyclopropyl analogues of combretastatin A4.

To evaluate the influence of stereochemistry on biological activities of cis-cyclopropyl combretastatin A4 (CA4) analogues, we have prepared several cyclopropyl compounds in their pure enantiomeric forms. The key reactions in our synthesis are the cyclopropanation of a (Z)-alkenylboron compound bearing a chiral auxiliary, and the cross-coupling of both enantiomeric cyclopropyl trifluoroborate salts with aryl and olefinic halides. Three pairs of cis-cyclopropyl CA4 analogues were evaluated for their potential antivascular activities. The diarylcyclopropyl compounds with SR-configuration (-)-1b, (-)-2b and the cyclopropylvinyl enantiomer (+)-3a with RR-configuration were the most potent tubulin polymerization inhibitors. A correlation was noted between anti-tubulin activity and rounding up activity of endothelial cells. The cytotoxic activity on B16 melanoma cells was in the submicromolar range for most compounds, but unlike the anti-tubulin activity, there was no difference in cytotoxic activity between racemic and enantiomerically pure forms for the three series of compounds. Molecular docking studies within the colchicine binding site of tubulin were in good agreement with the tubulin polymerization inhibitory data and confirmed the importance of the configuration of the synthesized cis-cyclopropyl CA4 analogues for potential antivascular activities.

Synthesis and properties of 2'-O-neopentyl modified oligonucleotides.

2'-O-Neopentyldeoxyuridine (Un) was synthesized and incorporated into a series of oligodeoxyribonucleotides. Single and triple incorporations in various arrangements were performed. The Watson and Crick pairing properties with complementary DNA and RNA were investigated by UV melting curves, CD spectroscopy, and molecular dynamic simulations. The results were compared to those obtained with DNA-DNA and DNA-RNA duplexes involving dU at the same positions. Oligonucleotides containing Un clearly demonstrated their ability to form duplexes with both complementary DNA and RNA but with higher stabilities for the DNA-RNA duplexes similar to the one of the parent DNA-RNA duplex. Investigations into the thermodynamic properties of these 17-base-pair duplexes revealed ΔG values (37 °C) that are in line with the measured T(m) values for both the DNA-DNA and DNA-RNA duplexes. CD spectroscopic structural investigations indicated that the conformations of the DNA-DNA and DNA-RNA duplexes involving Un are similar to those of the dT-rA and dU-rA containing duplexes. Only small changes in intensities and weak blue shifts were observed when three Uns were incorporated into the duplexes. The results of the molecular dynamic simulations showed, for the six duplexes involving the modified nucleoside Un, calculated curvatures similar to those of the corresponding unmodified duplexes without base-pair disruption. The neopentyl group is able to be accommodated in the minor grooves of both the DNA-DNA and RNA-DNA duplexes. However, molecular dynamic simulations indicated that the Uns adopt a C2'-exo sugar pucker conformation close to an A-helix type without perturbing the C2'-endo sugar pucker conformations of their 2'-deoxynucleoside neighbours. These results confirm the potential of 2'-O-neopentyldeoxyuridine as a nucleoside surrogate for oligonucleotide based therapeutic strategies.

Efficient synthesis and first regioselective C-3 direct arylation of imidazo[1,2-b]pyrazoles.

Highly regioselective: An efficient synthesis of the imidazo[1,2-b]pyrazole core has been developed, and the first regioselective palladium-catalyzed direct arylation of the C-3 position is described (see scheme). Good to excellent yields were obtained for a wide range of aryl partners with electron-rich and electron-poor substituents. This methodology allows rapid access to a large variety of imidazo[1,2-b]pyrazole products and could open the way to the design of new biologically active compounds.

Novel tetrahydropyrido[1,2-a]isoindolone derivatives (valmerins): potent cyclin-dependent kinase/glycogen synthase kinase 3 inhibitors with antiproliferative activities and antitumor effects in human tumor xenografts.

The development of CDK and GSK3 inhibitors has been regarded as a potential therapeutic approach, and a substantial number of diverse structures have been reported to inhibit CDKs and GSK-3ß in recent years. Only a few molecules have gone through or are currently undergoing clinical trials as CDK and GSK inhibitors. In this paper, we prepared valmerins, a new family containing the tetrahydropyrido[1,2-a]isoindone core. The fused heterocycle was prepared with a straightforward synthesis that was functionalized by a (het)arylurea. Twelve valmerins inhibited the CDK5 and GSK3 with an IC(50) < 100 nM. A semiquantitative kinase scoring was realized, and a cellular screening was done. At the end of our study, we investigated the in vivo potency of one valmerin. Mice exhibited good tolerance to our lead, which proved its efficacy and clearly blocked tumor growth. Valmerins appear also as good candidates for further development as anticancer agents.