Crystallographic fragment screening in academic cancer drug discovery.

Fragment-based drug discovery (FBDD) has brought several drugs to the clinic, notably to target proteins once considered to be challenging, or even undruggable. Screening in FBDD relies upon observing and/or measuring weak (millimolar-scale) binding events using biophysical techniques or crystallographic fragment screening. This latter structural approach provides no information about binding affinity but can reveal binding mode and atomic detail on protein-fragment interactions to accelerate hit-to-lead development. In recent years, high-throughput platforms have been developed at synchrotron facilities to screen thousands of fragment-soaked crystals. However, using accessible manual techniques it is possible to run informative, smaller-scale screens within an academic lab setting. This chapter describes general protocols for home laboratory-scale fragment screening, from fragment soaking through to structure solution and, where appropriate, signposts to background, protocols or alternatives elsewhere.

An open-source molecular builder and free energy preparation workflow.

Automated free energy calculations for the prediction of binding free energies of congeneric series of ligands to a protein target are growing in popularity, but building reliable initial binding poses for the ligands is challenging. Here, we introduce the open-source FEgrow workflow for building user-defined congeneric series of ligands in protein binding pockets for input to free energy calculations. For a given ligand core and receptor structure, FEgrow enumerates and optimises the bioactive conformations of the grown functional group(s), making use of hybrid machine learning/molecular mechanics potential energy functions where possible. Low energy structures are optionally scored using the gnina convolutional neural network scoring function, and output for more rigorous protein-ligand binding free energy predictions. We illustrate use of the workflow by building and scoring binding poses for ten congeneric series of ligands bound to targets from a standard, high quality dataset of protein-ligand complexes. Furthermore, we build a set of 13 inhibitors of the SARS-CoV-2 main protease from the literature, and use free energy calculations to retrospectively compute their relative binding free energies. FEgrow is freely available at https://github.com/cole-group/FEgrow, along with a tutorial.

Targeting Mre11 overcomes platinum resistance and induces synthetic lethality in XRCC1 deficient epithelial ovarian cancers.

Platinum resistance is a clinical challenge in ovarian cancer. Platinating agents induce DNA damage which activate Mre11 nuclease directed DNA damage signalling and response (DDR). Upregulation of DDR may promote chemotherapy resistance. Here we have comprehensively evaluated Mre11 in epithelial ovarian cancers. In clinical cohort that received platinum- based chemotherapy (n = 331), Mre11 protein overexpression was associated with aggressive phenotype and poor progression free survival (PFS) (p = 0.002). In the ovarian cancer genome atlas (TCGA) cohort (n = 498), Mre11 gene amplification was observed in a subset of serous tumours (5%) which correlated highly with Mre11 mRNA levels (p < 0.0001). Altered Mre11 levels was linked with genome wide alterations that can influence platinum sensitivity. At the transcriptomic level (n = 1259), Mre11 overexpression was associated with poor PFS (p = 0.003). ROC analysis showed an area under the curve (AUC) of 0.642 for response to platinum-based chemotherapy. Pre-clinically, Mre11 depletion by gene knock down or blockade by small molecule inhibitor (Mirin) reversed platinum resistance in ovarian cancer cells and in 3D spheroid models. Importantly, Mre11 inhibition was synthetically lethal in platinum sensitive XRCC1 deficient ovarian cancer cells and 3D-spheroids. Selective cytotoxicity was associated with DNA double strand break (DSB) accumulation, S-phase cell cycle arrest and increased apoptosis. We conclude that pharmaceutical development of Mre11 inhibitors is a viable clinical strategy for platinum sensitization and synthetic lethality in ovarian cancer.

FEN1 Blockade for Platinum Chemo-Sensitization and Synthetic Lethality in Epithelial Ovarian Cancers.

FEN1 plays critical roles in long patch base excision repair (LP-BER), Okazaki fragment maturation, and rescue of stalled replication forks. In a clinical cohort, FEN1 overexpression is associated with aggressive phenotype and poor progression-free survival after platinum chemotherapy. Pre-clinically, FEN1 is induced upon cisplatin treatment, and nuclear translocation of FEN1 is dependent on physical interaction with importin ß. FEN1 depletion, gene inactivation, or inhibition re-sensitizes platinum-resistant ovarian cancer cells to cisplatin. BRCA2 deficient cells exhibited synthetic lethality upon treatment with a FEN1 inhibitor. FEN1 inhibitor-resistant PEO1R cells were generated, and these reactivated BRCA2 and overexpressed the key repair proteins, POLß and XRCC1. FEN1i treatment was selectively toxic to POLß deficient but not XRCC1 deficient ovarian cancer cells. High throughput screening of 391,275 compounds identified several FEN1 inhibitor hits that are suitable for further drug development. We conclude that FEN1 is a valid target for ovarian cancer therapy.

Discriminative SKP2 Interactions with CDK-Cyclin Complexes Support a Cyclin A-Specific Role in p27KIP1 Degradation.

The SCFSKP2 ubiquitin ligase relieves G1 checkpoint control of CDK-cyclin complexes by promoting p27KIP1 degradation. We describe reconstitution of stable complexes containing SKP1-SKP2 and CDK1-cyclin B or CDK2-cyclin A/E, mediated by the CDK regulatory subunit CKS1. We further show that a direct interaction between a SKP2 N-terminal motif and cyclin A can stabilize SKP1-SKP2-CDK2-cyclin A complexes in the absence of CKS1. We identify the SKP2 binding site on cyclin A and demonstrate the site is not present in cyclin B or cyclin E. This site is distinct from but overlapping with features that mediate binding of p27KIP1 and other G1 cyclin regulators to cyclin A. We propose that the capacity of SKP2 to engage with CDK2-cyclin A by more than one structural mechanism provides a way to fine tune the degradation of p27KIP1 and distinguishes cyclin A from other G1 cyclins to ensure orderly cell cycle progression.

Chatterboxes: the structural and functional diversity of cyclins.

Proteins of the cyclin family have divergent sequences and execute diverse roles within the cell while sharing a common fold: the cyclin box domain. Structural studies of cyclins have played a key role in our characterization and understanding of cellular processes that they control, though to date only ten of the 29 CDK-activating cyclins have been structurally characterized by X-ray crystallography or cryo-electron microscopy with or without their cognate kinases. In this review, we survey the available structures of human cyclins, highlighting their molecular features in the context of their cellular roles. We pay particular attention to how cyclin activity is regulated through fine control of degradation motif recognition and ubiquitination. Finally, we discuss the emergent roles of cyclins independent of their roles as cyclin-dependent protein kinase activators, demonstrating the cyclin box domain to be a versatile and generalized scaffolding domain for protein-protein interactions across the cellular machinery.

Paradoxical activation of the protein kinase-transcription factor ERK5 by ERK5 kinase inhibitors.

The dual protein kinase-transcription factor, ERK5, is an emerging drug target in cancer and inflammation, and small-molecule ERK5 kinase inhibitors have been developed. However, selective ERK5 kinase inhibitors fail to recapitulate ERK5 genetic ablation phenotypes, suggesting kinase-independent functions for ERK5. Here we show that ERK5 kinase inhibitors cause paradoxical activation of ERK5 transcriptional activity mediated through its unique C-terminal transcriptional activation domain (TAD). Using the ERK5 kinase inhibitor, Compound 26 (ERK5-IN-1), as a paradigm, we have developed kinase-active, drug-resistant mutants of ERK5. With these mutants, we show that induction of ERK5 transcriptional activity requires direct binding of the inhibitor to the kinase domain. This in turn promotes conformational changes in the kinase domain that result in nuclear translocation of ERK5 and stimulation of gene transcription. This shows that both the ERK5 kinase and TAD must be considered when assessing the role of ERK5 and the effectiveness of anti-ERK5 therapeutics.

The diverse structural landscape of quadruplexes.

G-quadruplexes are secondary structures formed in G-rich sequences in DNA and RNA. Considerable research over the past three decades has led to in-depth insight into these unusual structures in DNA. Since the more recent exploration into RNA G-quadruplexes, such structures have demonstrated their in cellulo existence, function and roles in pathology. In comparison to Watson-Crick-based secondary structures, most G-quadruplexes display highly redundant structural characteristics. However, numerous reports of G-quadruplex motifs/structures with unique features (e.g. bulges, long loops, vacancy) have recently surfaced, expanding the repertoire of G-quadruplex scaffolds. This review addresses G-quadruplex formation and structure, including recent reports of non-canonical G-quadruplex structures. Improved methods of detection will likely further expand this collection of novel structures and ultimately change the face of quadruplex-RNA targeting as a therapeutic strategy.

Relative Binding Energies Predict Crystallographic Binding Modes of Ethionamide Booster Lead Compounds.

Transcriptional repressor EthR from Mycobacterium tuberculosis is a valuable target for antibiotic booster drugs. We previously reported a virtual screening campaign to identify EthR inhibitors for development. Two ligand binding orientations were often proposed, though only the top scoring pose was utilized for filtering of the large data set. We obtained biophysically validated hits, some of which yielded complex crystal structures. In some cases, the crystallized binding mode and top scoring mode agree, while for others an alternate ligand binding orientation was found. In this contribution, we combine rigid docking, molecular dynamics simulations, and the linear interaction energy method to calculate binding free energies and derive relative binding energies for a number of EthR inhibitors in both modes. This strategy allowed us to correctly predict the most favorable orientation. Therefore, this widely applicable approach will be suitable to triage multiple binding modes within EthR and other potential drug targets with similar characteristics.

Differences in the Conformational Energy Landscape of CDK1 and CDK2 Suggest a Mechanism for Achieving Selective CDK Inhibition.

Dysregulation of the cell cycle characterizes many cancer subtypes, providing a rationale for developing cyclin-dependent kinase (CDK) inhibitors. Potent CDK2 inhibitors might target certain cancers in which CCNE1 is amplified. However, current CDK2 inhibitors also inhibit CDK1, generating a toxicity liability. We have used biophysical measurements and X-ray crystallography to investigate the ATP-competitive inhibitor binding properties of cyclin-free and cyclin-bound CDK1 and CDK2. We show that these kinases can readily be distinguished by such inhibitors when cyclin-free, but not when cyclin-bound. The basis for this discrimination is unclear from either inspection or molecular dynamics simulation of ligand-bound CDKs, but is reflected in the contacts made between the kinase N- and C-lobes. We conclude that there is a subtle but profound difference between the conformational energy landscapes of cyclin-free CDK1 and CDK2. The unusual properties of CDK1 might be exploited to differentiate CDK1 from other CDKs in future cancer therapeutic design.

Probing biological activity through structural modelling of ligand-receptor interactions of 2,4-disubstituted thiazole retinoids.

Retinoids, such as all-trans-retinoic acid (ATRA), regulate cellular differentiation and signalling pathways in chordates by binding to nuclear retinoic acid receptors (RARα/ß/γ). Polar interactions between receptor and ligand are important for binding and facilitating the non-polar interactions and conformational changes necessary for RAR-mediated transcriptional regulation. The constraints on activity and RAR-type specificity with respect to the structural link between the polar and non-polar functions of synthetic retinoids are poorly understood. To address this, predictions from in silico ligand-RAR docking calculations and molecular dynamics simulations for a small library of stable, synthetic retinoids (designated GZ series) containing a central thiazole linker structure and different hydrophobic region substituents, were tested using a ligand binding assay and a range of cellular biological assays. The docking analysis showed that these thiazole-containing retinoids were well suited to the binding pocket of RARα, particularly via a favorable hydrogen bonding interaction between the thiazole and Ser232 of RARα. A bulky hydrophobic region (i.e., present in compounds GZ23 and GZ25) was important for interaction with the RAR binding pockets. Ligand binding assays generally reflected the findings from in silico docking, and showed that GZ25 was a particularly strongly binding ligand for RARα/ß. GZ25 also exhibited higher activity as an inducer of neuronal differentiation than ATRA and other GZ derivatives. These data demonstrate that GZ25 is a stable synthetic retinoid with improved activity which efficiently regulates neuronal differentiation and help to define the key structural requirements for retinoid activity enabling the design and development of the next generation of more active, selective synthetic retinoids as potential therapeutic regulators of neurogenesis.

New active leads for tuberculosis booster drugs by structure-based drug discovery.

The transcriptional repressor EthR from Mycobacterium tuberculosis, a member of the TetR family of prokaryotic homo-dimeric transcription factors, controls the expression of the mycobacterial mono-oxygenase EthA. EthA is responsible for the bio-activation of the second-line tuberculosis pro-drug ethionamide, and consequently EthR inhibitors boost drug efficacy. Here, we present a comprehensive in silico structure-based screening protocol that led to the identification of a number of novel scaffolds of EthR inhibitors in subsequent biophysical screening by thermal shift assay. Growth inhibition assays demonstrated that five of the twenty biophysical hits were capable of boosting ethionamide activity in vitro, with the best novel scaffold displaying an EC50 of 34 µM. In addition, the co-crystal structures of EthR with four new ligands at resolution ranging from 2.1 to 1.4 Å confirm the binding and inactivation mode, and will enable future lead development.

Structural and docking studies of potent ethionamide boosters.

Tuberculosis remains the second only to HIV as the leading cause of death by infectious disease worldwide, and was responsible for 1.4 million deaths globally in 2011. One of the essential drugs of the second-line antitubercular regimen is the prodrug ethionamide, introduced in the 1960s. Ethionamide is primarily used in cases of multi-drug resistant (MDR) and extensively drug resistant (XDR) TB due to severe adverse side effects. As a prodrug, ethionamide is bioactivated by EthA, a mono-oxygenase whose activity is repressed by EthR, a member of the TetR family of regulators. Previous studies have established that inhibition of EthR improves ethionamide potency. We report here the crystal structures of three EthR inhibitors at 0.8 Šresolution (3-oxo-3-{4-[3-(thiophen-2-yl)-1,2,4-oxadiazol-5-yl]piperidin-1-yl}propanenitrile (BDM31343), 4,4,4-trifluoro-1-{4-[3-(6-methoxy-1,3-benzothiazol-2-yl)-1,2,4-oxadiazol-5-yl]piperidin-1-yl}butanone (BDM41325) and 5,5,5-trifluoro-1-{4-[3-(4-methanesulfonylphenyl)-1,2,4-oxadiazol-5-yl]piperidin-1-yl}pentanone (BDM41907)), and the docking studies undertaken to investigate possible binding modes. The results revealed two distinct orientations of the three compounds in the binding channel, a direct consequence of the promiscuous nature of the largely lipophilic binding site.