Marc A. Moesser

The main protease (M pro ) of SARS-CoV-2 is central to viral maturation and is a promising drug target. In silico methods reveal structural aspects of how it binds to its 11 natural cleavage sites, the design of novel peptide inhibitors, and insights into drug design.

Abstract Endocrine resistance (EnR) in advanced prostate cancer is fatal. EnR can be mediated by androgen receptor (AR) splice variants, with AR splice variant 7 (AR-V7) arguably the most clinically important variant. In this study, we determined proteins key to generating AR-V7, validated our findings using clinical samples, and studied splicing regulatory mechanisms in prostate cancer models. Triangulation studies identified JMJD6 as a key regulator of AR-V7, as evidenced by its upregulation with in vitro EnR, its downregulation alongside AR-V7 by bromodomain inhibition, and its identification as a top hit of a targeted siRNA screen of spliceosome-related genes. JMJD6 protein levels increased (P < 0.001) with castration resistance and were associated with higher AR-V7 levels and shorter survival (P = 0.048). JMJD6 knockdown reduced prostate cancer cell growth, AR-V7 levels, and recruitment of U2AF65 to AR pre-mRNA. Mutagenesis studies suggested that JMJD6 activity is key to the generation of AR-V7, with the catalytic machinery residing within a druggable pocket. Taken together, these data highlight the relationship between JMJD6 and AR-V7 in advanced prostate cancer and support further evaluation of JMJD6 as a therapeutic target in this disease. Significance: This study identifies JMJD6 as being critical for the generation of AR-V7 in prostate cancer, where it may serve as a tractable target for therapeutic intervention.

Abstract Graph Neural Networks (GNNs) have recently gained in popularity, challenging molecular fingerprints or SMILES-based representations as the predominant way to represent molecules for binding affinity prediction. Although simple ligand-based graphs alone are already useful for affinity prediction, better performance on multi-target datasets has been achieved with models that incorporate 3D structural information. Most recent advances utilize complex GNN architectures to capture 3D protein-ligand information by incorporating ligand-interacting protein atoms as additional nodes in the graphs; or by building a second protein-based graph in parallel. This expands the graph considerably while obfuscating the shape of the underlying ligand, diminishing the advantage that GNNs have when encoding molecular structures. There is therefore a need for a simple and elegant molecular graph representation that retains the topology of the ligand while simultaneously encoding 3D protein-ligand interactions. We present Protein-Ligand Interaction Graphs (PLIGs): a simple way of representing atom-atom contacts of 3D protein-ligand complexes as node features for GNNs. PLIGs featurize an atom node in the molecular graph by describing each atom’s properties as well as all atom-atom contacts made with protein atoms within a distance threshold. The edges of the graph are therefore identical to ligand-based graphs, but the nodes encode the 3D protein-ligand contacts. Since PLIGs are applicable to any GNN architecture, we have benchmarked their performance with six different GNN architectures, and compared them to conventional ligand-based graphs and fingerprint-based multi-layer perceptron (MLP) models using the CASF-2016 benchmark set where we found PLIG-based Graph Attention Networks (GATNet) to be the best performing model ( ρ =0.84, RMSE=1.22 pK). In summary, we created a novel graph-based representation that incorporates 3D structural information into the node features of ligand-shaped molecular graphs. The PLIG representation is simple, elegant, flexible and easily customizable, opening up many possibilities of incorporating other 2D and 3D properties into the graph. Access The code and implementation for PLIGs and all models can be found at github.com/MarcMoesser/Protein-Ligand-Interaction-Graphs .

The main protease (M pro ) of SARS-CoV-2 is central to its viral lifecycle and is a promising drug target, but little is known concerning structural aspects of how it binds to its 11 natural cleavage sites. We used biophysical and crystallographic data and an array of classical molecular mechanics and quantum mechanical techniques, including automated docking, molecular dynamics (MD) simulations, linear-scaling DFT, QM/MM, and interactive MD in virtual reality, to investigate the molecular features underlying recognition of the natural M pro substrates. Analyses of the subsite interactions of modelled 11-residue cleavage site peptides, ligands from high-throughput crystallography, and designed covalently binding inhibitors were performed. Modelling studies reveal remarkable conservation of hydrogen bonding patterns of the natural M pro substrates, particularly on the N-terminal side of the scissile bond. They highlight the critical role of interactions beyond the immediate active site in recognition and catalysis, in particular at the P2/S2 sites. The binding modes of the natural substrates, together with extensive interaction analyses of inhibitor and fragment binding to M pro , reveal new opportunities for inhibition. Building on our initial M pro -substrate models, computational mutagenesis scanning was employed to design peptides with improved affinity and which inhibit M pro competitively. The combined results provide new insight useful for the development of M pro inhibitors.

Supplementary Table from JMJD6 Is a Druggable Oxygenase That Regulates AR-V7 Expression in Prostate Cancer