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.

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.

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.

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.

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.

Discovery of RGT-018: a Potent, Selective and Orally Bioavailable SOS1 Inhibitor for KRAS-driven Cancers.

KRAS is the most frequently dysregulated oncogene with high prevalence in NSCLC, colorectal cancer, and pancreatic cancer. FDA-approved sotorasib and adagrasib provide breakthrough therapies for cancer patients with KRASG12C mutation. However, there is still high unmet medical need for new agents targeting broader KRAS-driven tumors. An emerging and promising opportunity is to develop a pan KRAS inhibitor by suppressing the upstream protein SOS1. SOS1 is a key activator of KRAS and facilitates the conversion of GDP-bound KRAS state to GTP-bound KRAS state. Binding to its catalytic domain, small molecule SOS1 inhibitor has demonstrated the ability to suppress KRAS activation and cancer cell proliferation. RGT-018, a potent and selective SOS1 inhibitor, was identified with optimal drug-like properties. In vitro, RGT-018 blocked the interaction of KRAS:SOS1 with single digit nM potency and is highly selective against SOS2. RGT-018 inhibited KRAS signaling and the proliferation of a broad spectrum of KRAS-driven cancer cells as a single agent in vitro. Further enhanced anti-proliferation activity was observed when RGT-018 was combined with MEK, KRASG12C, EGFR or CDK4/6 inhibitors. Oral administration of RGT-018 inhibited tumor growth and suppressed KRAS signaling in tumor xenografts in vivo. Combination with MEK or KRASG12C inhibitors led to significant tumor regression. Furthermore, RGT-018 overcame the resistance to the approved KRASG12C inhibitors caused by clinically acquired KRAS mutations either as a single agent or in combination. RGT-018 displayed promising pharmacological properties for combination with targeted agents to treat a broader KRAS-driven patient population.

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.

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.

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.

Establishing the relationship between in vitro potency, pharmacokinetic, and pharmacodynamic parameters in a series of orally available, hydroxyethylamine-derived ß-secretase inhibitors.

Sequential proteolytic cleavage of the amyloid precursor protein (APP) by ß-site APP-cleaving enzyme 1 (BACE1) and the γ-secretase complex produces the amyloid-ß peptide (Aß), which is believed to play a critical role in the pathology of Alzheimer's disease (AD). The aspartyl protease BACE1 catalyzes the rate-limiting step in the production of Aß, and as such it is considered to be an important target for drug development in AD. The development of a BACE1 inhibitor therapeutic has proven to be difficult. The active site of BACE1 is relatively large. Consequently, to achieve sufficient potency, many BACE1 inhibitors have required unfavorable physicochemical properties such as high molecular weight and polar surface area that are detrimental to efficient passage across the blood-brain barrier. Using a rational drug design approach we have designed and developed a new series of hydroxyethylamine-based inhibitors of BACE1 capable of lowering Aß levels in the brains of rats after oral administration. Herein we describe the in vitro and in vivo characterization of two of these molecules and the overall relationship of compound properties [e.g., in vitro permeability, P-glycoprotein (P-gp) efflux, metabolic stability, and pharmacological potency] to the in vivo pharmacodynamic effect with more than 100 compounds across the chemical series. We demonstrate that high in vitro potency for BACE1 was not sufficient to provide central efficacy. A combination of potency, high permeability, low P-gp-mediated efflux, and low clearance was required for compounds to produce robust central Aß reduction after oral dosing.

Structure guided P1' modifications of HEA derived ß-secretase inhibitors for the treatment of Alzheimer's disease.

The synthesis and SAR of a series of BACE-1 hydroxyethyl amine inhibitors containing substitutions on a spirocyclobutyl moiety is described. Selectivity against cathepsin D, a related aspartyl protease with potential off target toxicity, and improved microsomal stability is exemplified.

Discovery of selective biaryl ethers as PDE10A inhibitors: improvement in potency and mitigation of Pgp-mediated efflux.

We report the discovery of a novel series of biaryl ethers as potent and selective PDE10A inhibitors. Structure-activity studies improved the potency and decreased Pgp-mediated efflux found in the initial compound 4. X-ray crystallographic studies revealed two novel binding modes to the catalytic site of the PDE10A enzyme.

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.

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.

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.

Stereoselective synthesis of anti-N-protected 3-amino-1,2-epoxides by nucleophilic addition to N-tert-butanesulfinyl imine of a glyceraldehyde synthon.

A di-O-TBS protected glyceraldehyde synthon was condensed with Ellman's reagent to form a bench-stable N-tert-butanesulfinyl imine 6, which served as a common intermediate for the stereoselective introduction of various R groups. The Ellman adducts were converted to useful multifunctional intermediates 18a-i in one pot. The alcohols 18a-i were efficiently elaborated to both known and novel anti-N-protected-3-amino-1,2-epoxides in two steps. Compound 2a is a key intermediate toward HIV protease inhibitors.

Design and synthesis of 6-oxo-1,6-dihydropyridines as CDK5 inhibitors.

Cyclin-dependent kinase 5 (CDK5) is a serine-threonine protein kinase that plays a significant role in neuronal development. In association with p25, CDK5 abnormally phosphorylates a number of cellular targets involved in neurodegenerative disorders. Using active site homology and previous structure-activity relationships, a new series of potent CDK5 inhibitors was designed. This report describes the optimization of 6-oxo-1,6-dihydropyridines as CDK5 inhibitors.

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.