Cryo-EM structures reveal two allosteric inhibition modes of PI3KαH1047R involving a re-shaping of the activation loop.

PI3Kα is a lipid kinase that phosphorylates PIP2 and generates PIP3. The hyperactive PI3Kα mutation, H1047R, accounts for about 14% of breast cancer, making it a highly attractive target for drug discovery. Here, we report the cryo-EM structures of PI3KαH1047R bound to two different allosteric inhibitors QR-7909 and QR-8557 at a global resolution of 2.7 Å and 3.0 Å, respectively. The structures reveal two distinct binding pockets on the opposite sides of the activation loop. Structural and MD simulation analyses show that the allosteric binding of QR-7909 and QR-8557 inhibit PI3KαH1047R hyper-activity by reducing the fluctuation and mobility of the activation loop. Our work provides a strong rational basis for a further optimization and development of highly selective drug candidates to treat PI3KαH1047R-driven cancers.

A dual diffusion model enables 3D molecule generation and lead optimization based on target pockets.

Structure-based generative chemistry is essential in computer-aided drug discovery by exploring a vast chemical space to design ligands with high binding affinity for targets. However, traditional in silico methods are limited by computational inefficiency, while machine learning approaches face bottlenecks due to auto-regressive sampling. To address these concerns, we have developed a conditional deep generative model, PMDM, for 3D molecule generation fitting specified targets. PMDM consists of a conditional equivariant diffusion model with both local and global molecular dynamics, enabling PMDM to consider the conditioned protein information to generate molecules efficiently. The comprehensive experiments indicate that PMDM outperforms baseline models across multiple evaluation metrics. To evaluate the applications of PMDM under real drug design scenarios, we conduct lead compound optimization for SARS-CoV-2 main protease (Mpro) and Cyclin-dependent Kinase 2 (CDK2), respectively. The selected lead optimization molecules are synthesized and evaluated for their in-vitro activities against CDK2, displaying improved CDK2 activity.

Genome-wide analyses reveal the contribution of somatic variants to the immune landscape of multiple cancer types.

It has been well established that cancer cells can evade immune surveillance by mutating themselves. Understanding genetic alterations in cancer cells that contribute to immune regulation could lead to better immunotherapy patient stratification and identification of novel immune-oncology (IO) targets. In this report, we describe our effort of genome-wide association analyses across 22 TCGA cancer types to explore the associations between genetic alterations in cancer cells and 74 immune traits. Results showed that the tumor microenvironment (TME) is shaped by different gene mutations in different cancer types. Out of the key genes that drive multiple immune traits, top hit KEAP1 in lung adenocarcinoma (LUAD) was selected for validation. It was found that KEAP1 mutations can explain more than 10% of the variance for multiple immune traits in LUAD. Using public scRNA-seq data, further analysis confirmed that KEAP1 mutations activate the NRF2 pathway and promote a suppressive TME. The activation of the NRF2 pathway is negatively correlated with lower T cell infiltration and higher T cell exhaustion. Meanwhile, several immune check point genes, such as CD274 (PD-L1), are highly expressed in NRF2-activated cancer cells. By integrating multiple RNA-seq data, a NRF2 gene signature was curated, which predicts anti-PD1 therapy response better than CD274 gene alone in a mixed cohort of different subtypes of non-small cell lung cancer (NSCLC) including LUAD, highlighting the important role of KEAP1-NRF2 axis in shaping the TME in NSCLC. Finally, a list of overexpressed ligands in NRF2 pathway activated cancer cells were identified and could potentially be targeted for TME remodeling in LUAD.

Accelerated Discovery of Macrocyclic CDK2 Inhibitor QR-6401 by Generative Models and Structure-Based Drug Design.

Selective CDK2 inhibitors have the potential to provide effective therapeutics for CDK2-dependent cancers and for combating drug resistance due to high cyclin E1 (CCNE1) expression intrinsically or CCNE1 amplification induced by treatment of CDK4/6 inhibitors. Generative models that take advantage of deep learning are being increasingly integrated into early drug discovery for hit identification and lead optimization. Here we report the discovery of a highly potent and selective macrocyclic CDK2 inhibitor QR-6401 (23) accelerated by the application of generative models and structure-based drug design (SBDD). QR-6401 (23) demonstrated robust antitumor efficacy in an OVCAR3 ovarian cancer xenograft model via oral administration.

Structural insights into ligand recognition and selectivity of somatostatin receptors.

Somatostatin receptors (SSTRs) play versatile roles in inhibiting the secretion of multiple hormones such as growth hormone and thyroid-stimulating hormone, and thus are considered as targets for treating multiple tumors. Despite great progress made in therapeutic development against this diverse receptor family, drugs that target SSTRs still show limited efficacy with preferential binding affinity and conspicuous side-effects. Here, we report five structures of SSTR2 and SSTR4 in different states, including two crystal structures of SSTR2 in complex with a selective peptide antagonist and a non-peptide agonist, respectively, a cryo-electron microscopy (cryo-EM) structure of Gi1-bound SSTR2 in the presence of the endogenous ligand SST-14, as well as two cryo-EM structures of Gi1-bound SSTR4 in complex with SST-14 and a small-molecule agonist J-2156, respectively. By comparison of the SSTR structures in different states, molecular mechanisms of agonism and antagonism were illustrated. Together with computational and functional analyses, the key determinants responsible for ligand recognition and selectivity of different SSTR subtypes and multiform binding modes of peptide and non-peptide ligands were identified. Insights gained in this study will help uncover ligand selectivity of various SSTRs and accelerate the development of new molecules with better efficacy by targeting SSTRs.

Molecular insights into ago-allosteric modulation of the human glucagon-like peptide-1 receptor.

The glucagon-like peptide-1 (GLP-1) receptor is a validated drug target for metabolic disorders. Ago-allosteric modulators are capable of acting both as agonists on their own and as efficacy enhancers of orthosteric ligands. However, the molecular details of ago-allosterism remain elusive. Here, we report three cryo-electron microscopy structures of GLP-1R bound to (i) compound 2 (an ago-allosteric modulator); (ii) compound 2 and GLP-1; and (iii) compound 2 and LY3502970 (a small molecule agonist), all in complex with heterotrimeric Gs. The structures reveal that compound 2 is covalently bonded to C347 at the cytoplasmic end of TM6 and triggers its outward movement in cooperation with the ECD whose N terminus penetrates into the GLP-1 binding site. This allows compound 2 to execute positive allosteric modulation through enhancement of both agonist binding and G protein coupling. Our findings offer insights into the structural basis of ago-allosterism at GLP-1R and may aid the design of better therapeutics.

Author Correction: Cryo-EM structures of PAC1 receptor reveal ligand binding mechanism.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Cryo-EM structures of PAC1 receptor reveal ligand binding mechanism.

The pituitary adenylate cyclase-activating polypeptide type I receptor (PAC1R) belongs to the secretin receptor family and is widely distributed in the central neural system and peripheral organs. Abnormal activation of the receptor mediates trigeminovascular activation and sensitization, which is highly related to migraine, making PAC1R a potential therapeutic target. Elucidation of PAC1R activation mechanism would benefit discovery of therapeutic drugs for neuronal disorders. PAC1R activity is governed by pituitary adenylate cyclase-activating polypeptide (PACAP), known as a major vasodilator neuropeptide, and maxadilan, a native peptide from the sand fly, which is also capable of activating the receptor with similar potency. These peptide ligands have divergent sequences yet initiate convergent PAC1R activity. It is of interest to understand the mechanism of PAC1R ligand recognition and receptor activity regulation through structural biology. Here we report two near-atomic resolution cryo-EM structures of PAC1R activated by PACAP38 or maxadilan, providing structural insights into two distinct ligand binding modes. The structures illustrate flexibility of the extracellular domain (ECD) for ligands with distinct conformations, where ECD accommodates ligands in different orientations while extracellular loop 1 (ECL1) protrudes to further anchor the ligand bound in the orthosteric site. By structure-guided molecular modeling and mutagenesis, we tested residues in the ligand-binding pockets and identified clusters of residues that are critical for receptor activity. The structures reported here for the first time elucidate the mechanism of specificity and flexibility of ligand recognition and binding for PAC1R, and provide insights toward the design of therapeutic molecules targeting PAC1R.

Structure-guided discovery of a single-domain antibody agonist against human apelin receptor.

Developing antibody agonists targeting the human apelin receptor (APJ) is a promising therapeutic approach for the treatment of chronic heart failure. Here, we report the structure-guided discovery of a single-domain antibody (sdAb) agonist JN241-9, based on the cocrystal structure of APJ with an sdAb antagonist JN241, the first cocrystal structure of a class A G protein-coupled receptor (GPCR) with a functional antibody. As revealed by the structure, JN241 binds to the extracellular side of APJ, makes critical contacts with the second extracellular loop, and inserts the CDR3 into the ligand-binding pocket. We converted JN241 into a full agonist JN241-9 by inserting a tyrosine into the CDR3. Modeling and molecular dynamics simulation shed light on JN241-9-stimulated receptor activation, providing structural insights for finding agonistic antibodies against class A GPCRs.

Fragment-Based Computational Method for Designing GPCR Ligands.

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors, which is arguably the most important family of drug target. With the technology breakthroughs in X-ray crystallography and cryo-electron microscopy, more than 300 GPCR-ligand complex structures have been publicly reported since 2007, covering about 60 unique GPCRs. Such abundant structural information certainly will facilitate the structure-based drug design by targeting GPCRs. In this study, we have developed a fragment-based computational method for designing novel GPCR ligands. We first extracted the characteristic interaction patterns (CIPs) on the binding interfaces between GPCRs and their ligands. The CIPs were used as queries to search the chemical fragments derived from GPCR ligands, which were required to form similar interaction patterns with GPCR. Then, the selected chemical fragments were assembled into complete molecules by using the AutoT&T2 software. In this work, we chose ß-adrenergic receptor (ß-AR) and muscarinic acetylcholine receptor (mAChR) as the targets to validate this method. Based on the designs suggested by our method, samples of 63 compounds were purchased and tested in a cell-based functional assay. A total of 15 and 22 compounds were identified as active antagonists for ß2-AR and mAChR M1, respectively. Molecular dynamics simulations and binding free energy analysis were performed to explore the key interactions (e.g., hydrogen bonds and π-π interactions) between those active compounds and their target GPCRs. In summary, our work presents a useful approach to the de novo design of GPCR ligands based on the relevant 3D structural information.

Palladium-Promoted DNA-Compatible Heck Reaction.

Optimal conditions for palladium-promoted Heck reaction on DNA were developed with good to excellent conversions. Versatility with either DNA-conjugated styrene/acrylamide or aryl iodide and a broad substrate scope of the corresponding coupling partners were established. Furthermore, robustness of the Heck reaction conditions on single-strand DNA and feasibility for DNA-encoded library production were demonstrated.

GPCR structure and function relationship: identification of a biased apelin receptor mutant.

Biased ligands of G protein-coupled receptors (GPCRs) may have improved therapeutic benefits and safety profiles. However, the molecular mechanism of GPCR biased signaling remains largely unknown. Using apelin receptor (APJ) as a model, we systematically investigated the potential effects of amino acid residues around the orthosteric binding site on biased signaling. We discovered that a single residue mutation I109A (I1093.32) in the transmembrane domain 3 (TM3) located in the deep ligand-binding pocket was sufficient to convert a balanced APJ into a G protein signaling biased receptor. APJ I109A mutant receptor retained full capabilities in ligand binding and G protein activation, but was defective in GRK recruitment, ß-arrestin recruitment, and downstream receptor-mediated ERK activation. Based on molecular dynamics simulations, we proposed a molecular mechanism for biased signaling of I109A mutant receptor. We postulate that due to the extra space created by I109A mutation, the phenyl group of the last residue (Phe-13) of apelin rotates down and initiates a cascade of conformational changes in TM3. Phe-13 formed a new cluster of hydrophobic interactions with the sidechains of residues in TM3, including F1103.33 and M1133.36, which stabilizes the mutant receptor in a conformation favoring biased signaling. Interruption of these stabilizing interactions by double mutation F110A/I109A or M113A/I109A largely restored the ß-arrestin-mediated signaling. Taken together, we describe herein the discovery of a biased APJ mutant receptor and provide detailed molecular insights into APJ signaling selectivity, facilitating the discovery of novel therapeutics targeting APJ.

Inverse-Electron-Demand Diels-Alder Reactions for the Synthesis of Pyridazines on DNA.

The synthesis of pyridazines on DNA has been developed on the basis of inverse-electron-demand Diels-Alder (IEDDA) reactions of 1,2,4,5-tetrazines. The broad substrate scope is explored. Functionalized pyridazine products are selected for subsequent DNA-compatible Suzuki-Miyaura coupling, acylation, and SNAr substitution reactions, demonstrating the feasibility and versatility of IEDDA reactions for DNA-encoded library synthesis.

Discovery of TRPM8 Antagonist ( S)-6-(((3-Fluoro-4-(trifluoromethoxy)phenyl)(3-fluoropyridin-2-yl)methyl)carbamoyl)nicotinic Acid (AMG 333), a Clinical Candidate for the Treatment of Migraine.

Transient-receptor-potential melastatin 8 (TRPM8), the predominant mammalian cold-temperature thermosensor, is a nonselective cation channel expressed in a subpopulation of sensory neurons in the peripheral nervous system, including nerve circuitry implicated in migraine pathogenesis: the trigeminal and pterygopalatine ganglia. Genomewide association studies have identified an association between TRPM8 and reduced risk of migraine. This disclosure focuses on medicinal-chemistry efforts to improve the druglike properties of initial leads, particularly removal of CYP3A4-induction liability and improvement of pharmacokinetic properties. A novel series of biarylmethanamide TRPM8 antagonists was developed, and a subset of leads were evaluated in preclinical toxicology studies to identify a clinical candidate with an acceptable preclinical safety profile leading to clinical candidate AMG 333, a potent and highly selective antagonist of TRPM8 that was evaluated in human clinical trials.

Ruthenium-Promoted C-H Activation Reactions between DNA-Conjugated Acrylamide and Aromatic Acids.

The first example of DNA-compatible C-H activation reaction between DNA-conjugated acrylamides and aromatic acids has been developed. This new transformation enables aromatic acid, previously considered as monofunctional building block, to act like a bifunctional building block for the DNA encoded library synthesis. The general scope of aromatic acid was established for this new on-DNA C-H activation, which paved the way for its application in combinatorial library preparation.

Structural Basis for Apelin Control of the Human Apelin Receptor.

Apelin receptor (APJR) is a key regulator of human cardiovascular function and is activated by two different endogenous peptide ligands, apelin and Elabela, each with different isoforms diversified by length and amino acid sequence. Here we report the 2.6-Å resolution crystal structure of human APJR in complex with a designed 17-amino-acid apelin mimetic peptide agonist. The structure reveals that the peptide agonist adopts a lactam constrained curved two-site ligand binding mode. Combined with mutation analysis and molecular dynamics simulations with apelin-13 binding to the wild-type APJR, this structure provides a mechanistic understanding of apelin recognition and binding specificity. Comparison of this structure with that of other peptide receptors suggests that endogenous peptide ligands with a high degree of conformational flexibility may bind and modulate the receptors via a similar two-site binding mechanism.

Preclinical Evaluation of AMG 337, a Highly Selective Small Molecule MET Inhibitor, in Hepatocellular Carcinoma.

Aberrant hepatocyte growth factor (HGF)/MET signaling has been implicated in hepatocarcinogenesis, suggesting that MET may serve as an attractive therapeutic target in hepatocellular carcinoma. We sought to investigate the in vitro and in vivo antitumor activity of AMG 337, a potent and highly selective small molecule MET kinase inhibitor, in preclinical models of hepatocellular carcinoma. The antiproliferative activity of AMG 337 was evaluated across a panel of hepatocellular carcinoma cell lines in a viability assay. Daily oral administration was used to evaluate the in vivo antitumor activity of AMG 337 in two patient-derived xenograft (PDX) models of hepatocellular carcinoma (LI0612 and LI1078). AMG 337 exerted potent antiproliferative activity against 2 of 40 hepatocellular carcinoma cell lines, namely, MHCC97H (IC50, 0.015 µmol/L) and HCCLM3 (IC50, 0.025 µmol/L). Both sensitive cell lines showed MET amplification (MET/CEN-7 >2.0) assessed by FISH, and high MET expression (3+ IHC) assessed by IHC. AMG 337 potently inhibited p-MET in all cell lines with detectable levels of total MET. However, the dose-dependent inhibition of downstream effectors of HGF/MET signaling, including p-GAB1, p-AKT, and p-ERK, was limited to those cell lines sensitive to AMG 337 in a viability assay (MHCC97H and HCCLM3). AMG 337 significantly inhibited tumor growth at all doses tested in the MET-amplified and MET-high-expressing hepatocellular carcinoma PDX model LI0612 and had no effect on tumor growth in the non-MET-amplified and MET-low-expressing hepatocellular carcinoma PDX model LI1078. AMG 337 represents a promising and novel therapeutic strategy for targeting hepatocellular carcinomas with a dependence on HGF/MET signaling. Mol Cancer Ther; 15(6); 1227-37. ©2016 AACR.

Development of 2-aminooxazoline 3-azaxanthenes as orally efficacious ß-secretase inhibitors for the potential treatment of Alzheimer's disease.

The ß-site amyloid precursor protein (APP) cleaving enzyme 1 (BACE1) is one of the most hotly pursued targets for the treatment of Alzheimer's disease. We used a structure- and property-based drug design approach to identify 2-aminooxazoline 3-azaxanthenes as potent BACE1 inhibitors which significantly reduced CSF and brain Aß levels in a rat pharmacodynamic model. Compared to the initial lead 2, compound 28 exhibited reduced potential for QTc prolongation in a non-human primate cardiovascular safety model.

Hydroxyethylamine-based inhibitors of BACE1: P1-P3 macrocyclization can improve potency, selectivity, and cell activity.

We describe a systematic study of how macrocyclization in the P1-P3 region of hydroxyethylamine-based inhibitors of ß-site amyloid precursor protein (APP)-cleaving enzyme (BACE1) modulates in vitro activity. This study reveals that in a number of instances macrocyclization of bis-terminal dienes leads to improved potency toward BACE1 and selectivity against cathepsin D (CatD), as well as greater amyloid ß-peptide (Aß)-lowering activity in HEK293T cells stably expressing APPSW. However, for several closely related analogs the benefits of macrocyclization are attenuated by the effects of other structural features in different regions of the molecules. X-ray crystal structures of three of these novel macrocyclic inhibitors bound to BACE1 revealed their binding conformations and interactions with the enzyme.