Selective Aurora A-TPX2 Interaction Inhibitors Have In Vivo Efficacy as Targeted Antimitotic Agents.

Aurora A kinase, a cell division regulator, is frequently overexpressed in various cancers, provoking genome instability and resistance to antimitotic chemotherapy. Localization and enzymatic activity of Aurora A are regulated by its interaction with the spindle assembly factor TPX2. We have used fragment-based, structure-guided lead discovery to develop small molecule inhibitors of the Aurora A-TPX2 protein-protein interaction (PPI). Our lead compound, CAM2602, inhibits Aurora A:TPX2 interaction, binding Aurora A with 19 nM affinity. CAM2602 exhibits oral bioavailability, causes pharmacodynamic biomarker modulation, and arrests the growth of tumor xenografts. CAM2602 acts by a novel mechanism compared to ATP-competitive inhibitors and is highly specific to Aurora A over Aurora B. Consistent with our finding that Aurora A overexpression drives taxane resistance, these inhibitors synergize with paclitaxel to suppress the outgrowth of pancreatic cancer cells. Our results provide a blueprint for targeting the Aurora A-TPX2 PPI for cancer therapy and suggest a promising clinical utility for this mode of action.

Melodia: A Python Library for Protein Structure Analysis.

SUMMARY: Analysing protein structure similarities is an important step in protein engineering and drug discovery. Methodologies that are more advanced than simple RMSD are available but often require extensive mathematical or computational knowledge for implementation. Grouping and optimising such tools in an efficient open-source library increases accessibility and encourages the adoption of more advanced metrics. Melodia is a Python library with a complete set of components devised for describing, comparing and analysing the shape of protein structures using differential geometry of three-dimensional curves and knot theory. It can generate robust geometric descriptors for thousands of shapes in just a few minutes. Those descriptors are more sensitive to structural feature variation than RMSD deviation. Melodia also incorporates sequence structural annotation and three-dimensional visualisations. AVAILABILITY AND IMPLEMENTATION: Melodia is an open-source Python library freely available on https://github.com/rwmontalvao/Melodia_py, along with interactive Jupyter Notebook tutorials. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Evolution and host-specific adaptation of Pseudomonas aeruginosa.

The major human bacterial pathogen Pseudomonas aeruginosa causes multidrug-resistant infections in people with underlying immunodeficiencies or structural lung diseases such as cystic fibrosis (CF). We show that a few environmental isolates, driven by horizontal gene acquisition, have become dominant epidemic clones that have sequentially emerged and spread through global transmission networks over the past 200 years. These clones demonstrate varying intrinsic propensities for infecting CF or non-CF individuals (linked to specific transcriptional changes enabling survival within macrophages); have undergone multiple rounds of convergent, host-specific adaptation; and have eventually lost their ability to transmit between different patient groups. Our findings thus explain the pathogenic evolution of P. aeruginosa and highlight the importance of global surveillance and cross-infection prevention in averting the emergence of future epidemic clones.

Evaluating Representation Learning on the Protein Structure Universe.

We introduce ProteinWorkshop, a comprehensive benchmark suite for representation learning on protein structures with Geometric Graph Neural Networks. We consider large-scale pre-training and downstream tasks on both experimental and predicted structures to enable the systematic evaluation of the quality of the learned structural representation and their usefulness in capturing functional relationships for downstream tasks. We find that: (1) large-scale pretraining on AlphaFold structures and auxiliary tasks consistently improve the performance of both rotation-invariant and equivariant GNNs, and (2) more expressive equivariant GNNs benefit from pretraining to a greater extent compared to invariant models. We aim to establish a common ground for the machine learning and computational biology communities to rigorously compare and advance protein structure representation learning. Our open-source codebase reduces the barrier to entry for working with large protein structure datasets by providing: (1) storage-efficient dataloaders for large-scale structural databases including AlphaFoldDB and ESM Atlas, as well as (2) utilities for constructing new tasks from the entire PDB. ProteinWorkshop is available at: github.com/a-r-j/ProteinWorkshop.

A Fragment-Based Competitive 19F†LB-NMR Platform For Hotspot-Directed Ligand Profiling.

Ligand binding hotspots are regions of protein surfaces that form particularly favourable interactions with small molecule pharmacophores. Targeting interactions with these hotspots maximises the efficiency of ligand binding. Existing methods are capable of identifying hotspots but often lack assays to quantify ligand binding and direct elaboration at these sites. Herein, we describe a fragment-based competitive 19F Ligand Based NMR (LB-NMR) screening platform that enables routine, quantitative ligand profiling focused at ligand-binding hotspots. As a proof of concept, the method was applied to 4'-phosphopantetheine adenylyltransferase (PPAT) from Mycobacterium abscessus (Mabs). X-ray crystallographic characterisation of the hits from a 960-member fragment screen identified three ligand-binding hotspots across the PPAT active site. From the fragment hits a collection of 19F reporter candidates were designed and synthesised. By rigorous prioritisation and use of optimisation workflows, a single 19F reporter molecule was generated for each hotspot. Profiling the binding of a set of structurally characterised ligands by competitive 19F LB-NMR with this suite of 19F reporters recapitulated the binding affinity and site ID assignments made by ITC and X-ray crystallography. This quantitative mapping of ligand binding events at hotspot level resolution establishes the utility of the fragment-based competitive 19F LB-NMR screening platform for hotspot-directed ligand profiling.

Outcomes of the EMDataResource cryo-EM Ligand Modeling Challenge.

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein-nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: Escherichia coli beta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Lipoarabinomannan modification as a source of phenotypic heterogeneity in host-adapted Mycobacterium abscessus isolates.

Mycobacterium abscessus is increasingly recognized as the causative agent of chronic pulmonary infections in humans. One of the genes found to be under strong evolutionary pressure during adaptation of M. abscessus to the human lung is embC which encodes an arabinosyltransferase required for the biosynthesis of the cell envelope lipoglycan, lipoarabinomannan (LAM). To assess the impact of patient-derived embC mutations on the physiology and virulence of M. abscessus, mutations were introduced in the isogenic background of M. abscessus ATCC 19977 and the resulting strains probed for phenotypic changes in a variety of in vitro and host cell-based assays relevant to infection. We show that patient-derived mutational variations in EmbC result in an unexpectedly large number of changes in the physiology of M. abscessus, and its interactions with innate immune cells. Not only did the mutants produce previously unknown forms of LAM with a truncated arabinan domain and 3-linked oligomannoside chains, they also displayed significantly altered cording, sliding motility, and biofilm-forming capacities. The mutants further differed from wild-type M. abscessus in their ability to replicate and induce inflammatory responses in human monocyte-derived macrophages and epithelial cells. The fact that different embC mutations were associated with distinct physiologic and pathogenic outcomes indicates that structural alterations in LAM caused by nonsynonymous nucleotide polymorphisms in embC may be a rapid, one-step, way for M. abscessus to generate broad-spectrum diversity beneficial to survival within the heterogeneous and constantly evolving environment of the infected human airway.

Computational analyses of drug resistance mutations in katG and emb complexes in Mycobacterium tuberculosis.

The number of antibiotic resistant pathogens is increasing rapidly, and with this comes a substantial socioeconomic cost that threatens much of the world. To alleviate this problem, we must use antibiotics in a more responsible and informed way, further our understanding of the molecular basis of drug resistance, and design new antibiotics. Here, we focus on a key drug-resistant pathogen, Mycobacterium tuberculosis, and computationally analyze trends in drug-resistant mutations in genes of the proteins embA, embB, embC, and katG, which play essential roles in the action of the first-line drugs ethambutol and isoniazid. We use docking to predict binding modes of isoniazid to katG that agree with suggested binding sites found in our laboratory using cryo-EM. Using mutant stability predictions, we recapitulate the idea that resistance occurs when katG's heme cofactor is destabilized rather than due to a decrease in affinity to isoniazid. Conversely, we have identified resistance mutations that affect the affinity of ethambutol more drastically than the affinity of the natural substrate of embB. With this, we illustrate that we can distinguish between the two types of drug resistance-cofactor destabilization and drug affinity reduction-suggesting potential uses in the prediction of novel drug-resistant mutations.

Outcomes of the EMDataResource Cryo-EM Ligand Modeling Challenge.

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein/nucleic-acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: E. coli beta-galactosidase with inhibitor, SARS-CoV-2 RNA-dependent RNA polymerase with covalently bound nucleotide analog, and SARS-CoV-2 ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. We found that (1) the quality of submitted ligand models and surrounding atoms varied, as judged by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics, and contact scores, and (2) a composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Mutational spectra are associated with bacterial niche.

As observed in cancers, individual mutagens and defects in DNA repair create distinctive mutational signatures that combine to form context-specific spectra within cells. We reasoned that similar processes must occur in bacterial lineages, potentially allowing decomposition analysis to detect both disruption of DNA repair processes and exposure to niche-specific mutagens. Here we reconstruct mutational spectra for 84 clades from 31 diverse bacterial species and find distinct mutational patterns. We extract signatures driven by specific DNA repair defects using hypermutator lineages, and further deconvolute the spectra into multiple signatures operating within different clades. We show that these signatures are explained by both bacterial phylogeny and replication niche. By comparing mutational spectra of clades from different environmental and biological locations, we identify niche-associated mutational signatures, and then employ these signatures to infer the predominant replication niches for several clades where this was previously obscure. Our results show that mutational spectra may be associated with sites of bacterial replication when mutagen exposures differ, and can be used in these cases to infer transmission routes for established and emergent human bacterial pathogens.

Structure-guided computational insecticide discovery targeting ß-N-acetyl-D-hexosaminidase of Ostrinia furnacalis.

Ostrinia furnacalis is a species of moth in the Crambidae family that is harmful to maize and other corn crops in Southeast Asia and the Western Pacific regions. Ostrinia furnacalis causes devastating losses to economically important corn fields. The ß-N-acetyl-D-hexosaminidase is an essential enzyme in O. furnacalis and its substrate binding +1 active site is different from that of the plants and humans ß-N-acetyl-D-hexosaminidases. To develop environment-friendly insecticides against OfHex1, we conducted structure-guided computational insecticide discovery to identify potential inhibitors that can bind the active site and inhibit the substrate binding and activity of the enzyme. We adopted a three-pronged strategy to conduct virtual screening using Glide and virtual screening workflow (VSW) in Schrödinger Suite-2022-3, against crystal structures of OfHex1 (PDB Id:3NSN), its homologue in humans (PDB Id: 1NP0) and Alphafold model of ß-N-acetyl-D-hexosaminidase from Trichogramma pretiosum, an egg parasitoid that protects the crops from O. furnacalis. A library of 20,313 commercially available and "insecticide-like" compounds was extracted from published literature. LigPrep enabled 44,943 ready-to-dock conformers generation. Glide docking revealed 18 OfHex1-specific hits that were absent in human and T. pretiosum screens. Reference docking was conducted using inhibitors/natural ligands in the crystal structures and hits with better docking scores than the reference were selected for MD simulations using Desmond to understand the stability of hit-target interactions. We noted five compounds that bound to OfHex1 TMX active-site based on their docking scores, consistent binding as noted by MD simulations and their insecticide/pesticide likeliness as noted by the Comprehensive Pesticide Likeness Analysis.Communicated by Ramaswamy H. Sarma.

Structural and functional basis of inositol hexaphosphate stimulation of NHEJ through stabilization of Ku-XLF interaction.

The classical Non-Homologous End Joining (c-NHEJ) pathway is the predominant process in mammals for repairing endogenous, accidental or programmed DNA Double-Strand Breaks. c-NHEJ is regulated by several accessory factors, post-translational modifications, endogenous chemical agents and metabolites. The metabolite inositol-hexaphosphate (IP6) stimulates c-NHEJ by interacting with the Ku70-Ku80 heterodimer (Ku). We report cryo-EM structures of apo- and DNA-bound Ku in complex with IP6, at 3.5 Å and 2.74 Å resolutions respectively, and an X-ray crystallography structure of a Ku in complex with DNA and IP6 at 3.7 Å. The Ku-IP6 interaction is mediated predominantly via salt bridges at the interface of the Ku70 and Ku80 subunits. This interaction is distant from the DNA, DNA-PKcs, APLF and PAXX binding sites and in close proximity to XLF binding site. Biophysical experiments show that IP6 binding increases the thermal stability of Ku by 2°C in a DNA-dependent manner, stabilizes Ku on DNA and enhances XLF affinity for Ku. In cells, selected mutagenesis of the IP6 binding pocket reduces both Ku accrual at damaged sites and XLF enrolment in the NHEJ complex, which translate into a lower end-joining efficiency. Thus, this study defines the molecular bases of the IP6 metabolite stimulatory effect on the c-NHEJ repair activity.

Clinically relevant mutations in the PhoR sensor kinase of host-adapted Mycobacterium abscessus isolates impact response to acidic pH and virulence.

IMPORTANCE: Difficult-to-treat pulmonary infections caused by nontuberculous mycobacteria of the Mycobacterium abscessus group have been steadily increasing in the USA and globally. Owing to the relatively recent recognition of M. abscessus as a human pathogen, basic and translational research to address critical gaps in diagnosis, treatment, and prevention of diseases caused by this microorganism has been lagging behind that of the better-known mycobacterial pathogen, Mycobacterium tuberculosis. To begin unraveling the molecular mechanisms of pathogenicity of M. abscessus, we here focus on the study of a two-component regulator known as PhoPR which we found to be under strong evolutionary pressure during human lung infection. We show that PhoPR is activated at acidic pH and serves to regulate a defined set of genes involved in host adaptation. Accordingly, clinical isolates from chronically infected human lungs tend to hyperactivate this regulator enabling M. abscessus to escape macrophage killing.

Cryo-EM structure of a DNA-PK trimer: higher order oligomerisation in NHEJ.

The ability of humans to maintain the integrity of the genome is imperative for cellular survival. DNA double-strand breaks (DSBs) are considered the most critical type of DNA lesion, which can ultimately lead to diseases including cancer. Non-homologous end joining (NHEJ) is one of two core mechanisms utilized to repair DSBs. DNA-PK is a key component in this process and has recently been shown to form alternate long-range synaptic dimers. This has led to the proposal that these complexes can be formed before transitioning to a short-range synaptic complex. Here we present cryo-EM data representing an NHEJ supercomplex consisting of a trimer of DNA-PK in complex with XLF, XRCC4, and DNA Ligase IV. This trimer represents a complex of both long-range synaptic dimers. We discuss the potential role of the trimeric structure, and possible higher order oligomers, as structural intermediates in the NHEJ mechanism, or as functional DNA repair centers.

PAXX binding to the NHEJ machinery explains functional redundancy with XLF.

Nonhomologous end joining is a critical mechanism that repairs DNA double-strand breaks in human cells. In this work, we address the structural and functional role of the accessory protein PAXX [paralog of x-ray repair cross-complementing protein 4 (XRCC4) and XRCC4-like factor (XLF)] in this mechanism. Here, we report high-resolution cryo-electron microscopy (cryo-EM) and x-ray crystallography structures of the PAXX C-terminal Ku-binding motif bound to Ku70/80 and cryo-EM structures of PAXX bound to two alternate DNA-dependent protein kinase (DNA-PK) end-bridging dimers, mediated by either Ku80 or XLF. We identify residues critical for the Ku70/PAXX interaction in vitro and in cells. We demonstrate that PAXX and XLF can bind simultaneously to the Ku heterodimer and act as structural bridges in alternate forms of DNA-PK dimers. Last, we show that engagement of both proteins provides a complementary advantage for DNA end synapsis and end joining in cells.

SARS-CoV-2 Omicron subvariant spike N405 unlikely to rapidly deamidate.

The RGD motif on the SARS-CoV-2 spike protein has been suggested to interact with RGD-binding integrins αVß3 and α5ß1 to enhance viral cell entry and alter downstream signaling cascades. The D405N mutation on the Omicron subvariant spike proteins, resulting in an RGN motif, has recently been shown to inhibit binding to integrin αVß3. Deamidation of asparagines in protein ligand RGN motifs has been demonstrated to generate RGD and RGisoD motifs that permit binding to RGD-binding integrins. Two asparagines, N481 and N501, on the Wild-type spike receptor-binding domain have been previously shown to have deamidation half-lives of 16.5 and 123 days, respectively, which may occur during the viral life cycle. Deamidation of Omicron subvariant N405 may recover the ability to interact with RGD-binding integrins. Thus, herein, all-atom molecular dynamics simulations of the Wild-type and Omicron subvariant spike protein receptor-binding domains were conducted to investigate the potential for asparagines, the Omicron subvariant N405 in particular, to assume the optimized geometry for deamidation to occur. In summary, the Omicron subvariant N405 was primarily found to be stabilized in a state unfavourable for deamidation after hydrogen bonding with downstream E406. Nevertheless, a small number of RGD or RGisoD motifs on the Omicron subvariant spike proteins may restore the ability to interact with RGD-binding integrins. The simulations also provided structural clarification regarding the deamidation rates of Wild-type N481 and N501 and highlighted the utility of tertiary structure dynamics information in predicting asparagine deamidation. Further work is needed to characterize the effects of deamidation on spike-integrin interactions.

In silico analyses of isoniazid and streptomycin resistance-associated mutations in Mycobacterium tuberculosis.

Multi-drug resistant tuberculosis is categorised by the World Health Organisation (WHO) as a public health crisis. In silico techniques were used to probe the structural basis of Mycobacterium tuberculosis resistance to isoniazid and streptomycin. Isoniazid resistance-associated mutations in InhA were predicted to reduce the binding affinity of NADH to InhA, without affecting INH-NAD (competitive-inhibitor) binding. Perturbation of the mutated residues was predicted (with the AlloSigMA server) to modulate the free energy of allosteric modulation of key binding site residues F41, F149, Y158 and W222. These results suggest that allosteric modulation of the protein structure may be key to the mechanism by which isoniazid resistance-associated mutations act. Mutations in the methyltransferase glucose-inhibited division gene B (GidB) are associated with streptomycin resistance. Molecular docking was carried out to predict the structure of the GidB bound to its substrate (s-adenosyl methionine). The effects of streptomycin resistance-associated mutations in GidB on protein stability and substrate binding were predicted (using SDM and mCSM-lig). All GidB mutants were predicted to disfavour SAM binding.

Two distinct long-range synaptic complexes promote different aspects of end processing prior to repair of DNA breaks by non-homologous end joining.

Non-homologous end joining is the major double-strand break repair (DSBR) pathway in mammals. DNA-PK is the hub and organizer of multiple steps in non-homologous end joining (NHEJ). Recent high-resolution structures show how two distinct NHEJ complexes "synapse" two DNA ends. One complex includes a DNA-PK dimer mediated by XLF, whereas a distinct DNA-PK dimer forms via a domain-swap mechanism where the C terminus of Ku80 from one DNA-PK protomer interacts with another DNA-PK protomer in trans. Remarkably, the distance between the two synapsed DNA ends in both dimers is the same (∼115 Å), which matches the distance observed in the initial description of an NHEJ long-range synaptic complex. Here, a mutational strategy is used to demonstrate distinct cellular function(s) of the two dimers: one promoting fill-in end processing, while the other promotes DNA end resection. Thus, the specific DNA-PK dimer formed (which may be impacted by DNA end structure) dictates the mechanism by which ends will be made ligatable.

Chemical Validation of Mycobacterium tuberculosis Phosphopantetheine Adenylyltransferase Using Fragment Linking and CRISPR Interference.

The coenzyme A (CoA) biosynthesis pathway has attracted attention as a potential target for much-needed novel antimicrobial drugs, including for the treatment of tuberculosis (TB), the lethal disease caused by Mycobacterium tuberculosis (Mtb). Seeking to identify inhibitors of Mtb phosphopantetheine adenylyltransferase (MtbPPAT), the enzyme that catalyses the penultimate step in CoA biosynthesis, we performed a fragment screen. In doing so, we discovered three series of fragments that occupy distinct regions of the MtbPPAT active site, presenting a unique opportunity for fragment linking. Here we show how, guided by X-ray crystal structures, we could link weakly-binding fragments to produce an active site binder with a KD <20 µM and on-target anti-Mtb activity, as demonstrated using CRISPR interference. This study represents a big step toward validating MtbPPAT as a potential drug target and designing a MtbPPAT-targeting anti-TB drug.

Human DNA-dependent protein kinase activation mechanism.

DNA-dependent protein kinase (DNA-PK), a multicomponent complex including the DNA-PK catalytic subunit and Ku70/80 heterodimer together with DNA, is central to human DNA damage response and repair. Using a DNA-PK-selective inhibitor (M3814), we identified from one dataset two cryo-EM structures of the human DNA-PK complex in different states, the intermediate state and the active state. Here we show that activation of the kinase is regulated through conformational changes caused by the binding ligand and the string region (residues 802-846) of the DNA-PK catalytic subunit, particularly the helix-hairpin-helix motif (residues 816-836) that interacts with DNA. These observations demonstrate the regulatory role of the ligand and explain why DNA-PK is DNA dependent. Cooperation and coordination among binding partners, disordered flexible regions and mechanically flexible HEAT repeats modulate the activation of the kinase. Together with previous findings, these results provide a better molecular understanding of DNA-PK catalysis.