Kinetic barriers in the isomerization of substituted ureas: implications for computer-aided drug design

Discovery of Potent Multikinase Type-II Inhibitors Targeting CDK5 in the DFG-out Inactive State with Promising Potential against Glioblastoma

Kinases have proven valuable targets in successful cancer drug discovery projects, but not yet for malignant brain tumors where type-II inhibition of cyclin-dependent kinase 5 (CDK5) stabilizing the DFG-out inactive state has potential for design of selective and clinically efficient drug candidates. In the absence of crystallographic evidence for a CDK5 DFG-out inactive state protein-ligand complex, for the first time, a model was designed using metadynamics/molecular dynamics simulations. Glide docking of the ZINC15 biogenic database identified [pyrimidin-2-yl]amino-furo[3,2-b]-furyl-urea/amide hit chemical scaffolds. For four selected analogues (4, 27, 36, and 42), potent effects on glioblastoma cell viability in U87-MG, T98G, and U251-MG cell lines and patient-derived cultures were generally observed (IC50s ∼ 10-40 μM at 72 h). Selectivity profiling against 11 homologous kinases revealed multikinase inhibition (CDK2, CDK5, CDK9, and GSK-3α/β), most potent for GSK-3α in the nanomolar range (IC50s ∼ 0.23-0.98 μM). These compounds may therefore have diverse anticancer mechanisms of action and are of considerable interest for lead optimization.

Frentizole, a Nontoxic Immunosuppressive Drug, and Its Analogs Display Antitumor Activity via Tubulin Inhibition

Antimitotic agents are one of the more successful types of anticancer drugs, but they suffer from toxicity and resistance. The application of approved drugs to new indications (i.e., drug repurposing) is a promising strategy for the development of new drugs. It relies on finding pattern similarities: drug effects to other drugs or conditions, similar toxicities, or structural similarity. Here, we recursively searched a database of approved drugs for structural similarity to several antimitotic agents binding to a specific site of tubulin, with the expectation of finding structures that could fit in it. These searches repeatedly retrieved frentizole, an approved nontoxic anti-inflammatory drug, thus indicating that it might behave as an antimitotic drug devoid of the undesired toxic effects. We also show that the usual repurposing approach to searching for targets of frentizole failed in most cases to find such a relationship. We synthesized frentizole and a series of analogs to assay them as antimitotic agents and found antiproliferative activity against HeLa tumor cells, inhibition of microtubule formation within cells, and arrest at the G2/M phases of the cell cycle, phenotypes that agree with binding to tubulin as the mechanism of action. The docking studies suggest binding at the colchicine site in different modes. These results support the repurposing of frentizole for cancer treatment, especially for glioblastoma.

Conformational study into N-alkyl-N'-aryl ureas to inform drug discovery

Ureas are an important functional group in small molecule drugs as well as having wider applications in organic chemistry. Understanding of their conformation is of critical importance for rational design of urea-containing bioactive compounds. Whilst the conformational preferences of biaryl ureas have been extensively studied, very little attention has been paid to alkylated analogues. We carried out a systematic study of N-aryl (phenyl and pyridyl)-N'-cyclopentyl ureas with differing N-methylation patterns using Well Tempered Metadynamics at a semi-empirical level in implicit water (GBSA) using Well-Tempered Metadynamics to generate their conformational free-energy landscapes. Geometries and energetics of the most relevant configurations are further refined using DFT level of theory. Validation for the computation was obtained by synthesis of all 8 analogues followed by conformational studies by X-ray crystallography and NMR. These findings reveal that the methylation pattern significantly affects the conformational preference of the system. Most notably, N-phenyl-N'-cyclopentyl urea is shown to adopt both the trans-trans, and cis-trans conformations with equal energy and that the cis-trans conformation can be significantly stabilised by the presence of an internal hydrogen bond to the N'-hydrogen. This study will be of utility for the design of N-alkyl-N'-aryl ureas as drug candidates.

Novel diarylamides and diarylureas with N-substitution dependent activity against medulloblastoma

Medulloblastoma - highly aggressive and heterogeneous tumours of the cerebellum - account for 15-20% of all childhood brain tumours, and are the most common high-grade childhood embryonal tumour of the central nervous system. Herein, potent in vitro anticancer activity against two established medulloblastoma cell lines of the sonic hedgehog subgroup, namely DAOY (p53 mutant) and ONS-76 (p53 wild type), has been achieved. A number of first-generation diarylamides and diarylureas were evaluated and activity is likely to be, in-part, conformation-dependent. The most active compound from this first-generation set of compounds, 1-naphthyl derivative 4b, was selected and a second-generation of compounds were optimised and tested for activity against the medulloblastoma cell lines. This process resulted in drug-like compounds with up to sixty times the activity (sub-micromolar) of the first-generation - thus providing potent new leads for further study.

Ditopic binuclear copper(II) complexes for DNA cleavage

Herein we present the synthesis of two ligands containing two di(2-picolyl)amine (DPA) units linked by either a 1,1'-(pyridine-2,6-diyl)bis(3-ethylurea) (L1) or a 1,1'-(1,3-phenylene)bis(3-ethylurea) (L2) spacer. The corresponding binuclear Cu^II and Zn^II complexes were prepared and isolated. The X-ray structures of the L1 ligand and the [Cu₂L1Cl₂]²⁺ complex evidence an unusual cis/trans conformation of one of the urea groups stabilized by an intramolecular hydrogen bond with the nitrogen atom of the pyridyl spacer. The Cu^II complexes form rather strong ternary complexes with phosphorylated anions. The [Cu₂L1]⁴⁺ complex presents a rather high affinity for pyrophosphate (logK₁₁ = 8.19 at pH 7, 25 °C), while [Cu₂L2]⁴⁺ stands out because of its strong binding to AMP^2- (logK₁₁ = 9.3 at pH 7, 25 °C). The interaction of the Cu^II complexes with deoxyribonucleic acid from calf thymus (ct-DNA) was monitored using circular dichroism (CD) and luminescence spectroscopies. These studies revealed a quite strong interaction of the complexes with ct-DNA (K_b = (6.4 ± 0.7) × 10³ for [Cu₂L1]⁴⁺ and K_b = (6.3 ± 1.0) × 10³ for [Cu₂L2]⁴⁺). Competition experiments carried out in the presence of methyl green and BAPPA (N¹,N³-Bis(4-amidinophenyl)propane-1,3-diamine) as major and minor groove competitors, respectively, confirm that the interaction of both complexes with DNA takes place through the minor groove, in agreement with docking studies. The [Cu₂L2]⁴⁺ complex is quite efficient in promoting the cleavage of the double-stranded pUC19 plasmid DNA, by favoring the conversion of the supercoiled form to the nicked form following a hydrolytic mechanism.

A Million Crystal Structures: The Whole Is Greater than the Sum of Its Parts

The founding in 1965 of what is now called the Cambridge Structural Database (CSD) has reaped dividends in numerous and diverse areas of chemical research. Each of the million or so crystal structures in the database was solved for its own particular reason, but collected together, the structures can be reused to address a multitude of new problems. In this Review, which is focused mainly on the last 10 years, we chronicle the contribution of the CSD to research into molecular geometries, molecular interactions, and molecular assemblies and demonstrate its value in the design of biologically active molecules and the solid forms in which they are delivered. Its potential in other commercially relevant areas is described, including gas storage and delivery, thin films, and (opto)electronics. The CSD also aids the solution of new crystal structures. Because no scientific instrument is without shortcomings, the limitations of CSD research are assessed. We emphasize the importance of maintaining database quality: notwithstanding the arrival of big data and machine learning, it remains perilous to ignore the principle of garbage in, garbage out. Finally, we explain why the CSD must evolve with the world around it to ensure it remains fit for purpose in the years ahead.

A computationally designed binding mode flip leads to a novel class of potent tri-vector cyclophilin inhibitors

Cyclophilins (Cyps) are a major family of drug targets that are challenging to prosecute with small molecules because the shallow nature and high degree of conservation of the active site across human isoforms offers limited opportunities for potent and selective inhibition. Herein a computational approach based on molecular dynamics simulations and free energy calculations was combined with biophysical assays and X-ray crystallography to explore a flip in the binding mode of a reported urea-based Cyp inhibitor. This approach enabled access to a distal pocket that is poorly conserved among key Cyp isoforms, and led to the discovery of a new family of sub-micromolar cell-active inhibitors that offer unprecedented opportunities for the development of next-generation drug therapies based on Cyp inhibition. The computational approach is applicable to a broad range of organic functional groups and could prove widely enabling in molecular design.

Ligand Isomerism in Coordination Cages

Complexation reactions of palladium(II) nitrate with a set of 3-pyridyl appended nonchelating bidentate ligands possessing regioisomeric phenylene-diurea functionalities as spacers were carried out. The ligands utilized in this study are 1,1'-(1,2-phenylene)bis(3-(pyridin-3-yl)urea), L1; 1,1'-(1,3-phenylene)bis(3-(pyridin-3-yl)urea), L2; and 1,1'-(1,4-phenylene)bis(3-(pyridin-3-yl)urea), L3. The complexation reactions of the ligands (L1, L2, and L3) with palladium(II) produced single discrete isomeric cages (1, 2, and 3) of Pd₂L₄ formulation in each case and thereby illustrated ligand-isomerism in coordination cages. All 16 hydrogen atoms of eight urea moieties present in four ligand strands are delineated completely endohedrally in cage 1 and completely exohedrally in cage 3, whereas cage 2 exhibited half of the urea hydrogens in exohedral locations and the remaining half in endohedral locations. In addition to the variable number of solvent molecules, the cavities of cages 1 and 2 lodged four and two nitrate ions, respectively, using the endohedral (H)_urea atoms (i.e., NH groups) as binding sites, whereas the cavity of 3 remained anion free. The abilities of the complexes 1-3 for adsorption of CO₂ gas are demonstrated, and their behaviors are compared.

Towards the Prediction of Global Solution State Properties for Hydrogen Bonded, Self-Associating Amphiphiles

Through this extensive structure-property study we show that critical micelle concentration correlates with self-associative hydrogen bond complex formation constant, when combined with outputs from low level, widely accessible, computational models. Herein, we bring together a series of 39 structurally related molecules related by stepwise variation of a hydrogen bond donor-acceptor amphiphilic salt. The self-associative and corresponding global properties for this family of compounds have been studied in the gas, solid and solution states. Within the solution state, we have shown the type of self-associated structure present to be solvent dependent. In DMSO, this class of compound show a preference for hydrogen bonded dimer formation, however moving into aqueous solutions the same compounds are found to form larger self-associated aggregates. This observation has allowed us the unique opportunity to investigate and begin to predict self-association events at both the molecular and extended aggregate level.