Small Molecule Allosteric Modulators of G-Protein-Coupled Receptors: Drug–Target Interactions

G protein-coupled receptors (GPCRs): advances in structures, mechanisms, and drug discovery

G protein-coupled receptors (GPCRs), the largest family of human membrane proteins and an important class of drug targets, play a role in maintaining numerous physiological processes. Agonist or antagonist, orthosteric effects or allosteric effects, and biased signaling or balanced signaling, characterize the complexity of GPCR dynamic features. In this study, we first review the structural advancements, activation mechanisms, and functional diversity of GPCRs. We then focus on GPCR drug discovery by revealing the detailed drug-target interactions and the underlying mechanisms of orthosteric drugs approved by the US Food and Drug Administration in the past five years. Particularly, an up-to-date analysis is performed on available GPCR structures complexed with synthetic small-molecule allosteric modulators to elucidate key receptor-ligand interactions and allosteric mechanisms. Finally, we highlight how the widespread GPCR-druggable allosteric sites can guide structure- or mechanism-based drug design and propose prospects of designing bitopic ligands for the future therapeutic potential of targeting this receptor family.

Enhanced mapping of small-molecule binding sites in cells

Photoaffinity probes are routinely utilized to identify proteins that interact with small molecules. However, despite this common usage, resolving the specific sites of these interactions remains a challenge. Here we developed a chemoproteomic workflow to determine precise protein binding sites of photoaffinity probes in cells. Deconvolution of features unique to probe-modified peptides, such as their tendency to produce chimeric spectra, facilitated the development of predictive models to confidently determine labeled sites. This yielded an expansive map of small-molecule binding sites on endogenous proteins and enabled the integration with multiplexed quantitation, increasing the throughput and dimensionality of experiments. Finally, using structural information, we characterized diverse binding sites across the proteome, providing direct evidence of their tractability to small molecules. Together, our findings reveal new knowledge for the analysis of photoaffinity probes and provide a robust method for high-resolution mapping of reversible small-molecule interactions en masse in native systems.

Allosteric modulation of G protein-coupled receptors as a novel therapeutic strategy in neuropathic pain

Neuropathic pain is a debilitating pathological condition that presents significant therapeutic challenges in clinical practice. Unfortunately, current pharmacological treatments for neuropathic pain lack clinical efficacy and often lead to harmful adverse reactions. As G protein-coupled receptors (GPCRs) are widely distributed throughout the body, including the pain transmission pathway and descending inhibition pathway, the development of novel neuropathic pain treatments based on GPCRs allosteric modulation theory is gaining momentum. Extensive research has shown that allosteric modulators targeting GPCRs on the pain pathway can effectively alleviate symptoms of neuropathic pain while reducing or eliminating adverse effects. This review aims to provide a comprehensive summary of the progress made in GPCRs allosteric modulators in the treatment of neuropathic pain, and discuss the potential benefits and adverse factors of this treatment. We will also concentrate on the development of biased agonists of GPCRs, and based on important examples of biased agonist development in recent years, we will describe universal strategies for designing structure-based biased agonists. It is foreseeable that, with the continuous improvement of GPCRs allosteric modulation and biased agonist theory, effective GPCRs allosteric drugs will eventually be available for the treatment of neuropathic pain with acceptable safety.

Designing drugs and chemical probes with the dualsteric approach

Traditionally, drugs are monovalent, targeting only one site on the protein surface. This includes orthosteric and allosteric drugs, which bind the protein at orthosteric and allosteric sites, respectively. Orthosteric drugs are good in potency, whereas allosteric drugs have better selectivity and are solutions to classically undruggable targets. However, it would be difficult to simultaneously reach high potency and selectivity when targeting only one site. Also, both kinds of monovalent drugs suffer from mutation-caused drug resistance. To overcome these obstacles, dualsteric modulators have been proposed in the past twenty years. Compared to orthosteric or allosteric drugs, dualsteric modulators are bivalent (or bitopic) with two pharmacophores. Each of the two pharmacophores bind the protein at the orthosteric and an allosteric site, which could bring the modulator with special properties beyond monovalent drugs. In this study, we comprehensively review the current development of dualsteric modulators. Our main effort reason and illustrate the aims to apply the dualsteric approach, including a "double win" of potency and selectivity, overcoming mutation-caused drug resistance, developments of function-biased modulators, and design of partial agonists. Moreover, the strengths of the dualsteric technique also led to its application outside pharmacy, including the design of highly sensitive fluorescent tracers and usage as molecular rulers. Besides, we also introduced drug targets, designing strategies, and validation methods of dualsteric modulators. Finally, we detail the conclusions and perspectives.

Pharmacologically targeting intracellular allosteric sites of GPCRs for drug discovery

G-protein-coupled receptors (GPCRs) are a family of cell surface proteins that can sense a variety of extracellular stimuli and mediate multiple signaling transduction pathways involved in human physiology. Recent advances in GPCR structural biology have revealed a relatively conserved intracellular allosteric site in multiple GPCRs, which can be utilized to modulate receptors from the inside. This novel intracellular site partially overlaps with the G-protein and β-arrestin coupling sites, providing a novel avenue for biological intervention. Here, we review evidence available for GPCR structures complexed with intracellular small-molecule allosteric modulators, elucidating drug-target interactions and allosteric mechanisms. Moreover, we highlight the potential of intracellular allosteric modulators in achieving biased signaling, which provides insights into biased allosteric mechanisms.

An inverse agonist of orphan receptor GPR61 acts by a G protein-competitive allosteric mechanism

GPR61 is an orphan GPCR related to biogenic amine receptors. Its association with phenotypes relating to appetite makes it of interest as a druggable target to treat disorders of metabolism and body weight, such as obesity and cachexia. To date, the lack of structural information or a known biological ligand or tool compound has hindered comprehensive efforts to study GPR61 structure and function. Here, we report a structural characterization of GPR61, in both its active-like complex with heterotrimeric G protein and in its inactive state. Moreover, we report the discovery of a potent and selective small-molecule inverse agonist against GPR61 and structural elucidation of its allosteric binding site and mode of action. These findings offer mechanistic insights into an orphan GPCR while providing both a structural framework and tool compound to support further studies of GPR61 function and modulation.

Opportunities and challenges in drug discovery targeting the orphan receptor GPR12

G-protein-coupled receptor 12 (GPR12) is a brain-specific expression orphan G-protein-coupled receptor (oGPCR) that regulates various physiological processes. It is an emerging therapeutic target for central nervous system (CNS) disorders, including Alzheimer's disease (AD), Parkinson's disease (PD), Huntington's disease (HD), attention deficit hyperactivity disorder (ADHD), and schizophrenia, as well as other human diseases, such as cancer, obesity, and metabolic disorders. GPR12 remains a less extensively investigated oGPCR, particularly in terms of its biological functions, signaling pathways, and ligand discovery. The discovery of drug-like small-molecule modulators to probe the brain functions of GPR12 or to act as a potential drug candidates, as well as the identification of reliable biomarkers, are vital to elucidate the roles of this receptor in various human diseases and develop novel target-based therapeutics.

Structure-Activity Relationship Study of the High-Affinity Neuropeptide Y4 Receptor Positive Allosteric Modulator VU0506013

Positive allosteric modulators targeting the Y₄ receptor (Y₄R), a G protein-coupled receptor (GPCR) involved in the regulation of satiety, offer great potential in anti-obesity research. In this study, we selected 603 compounds by using quantitative structure-activity relationship (QSAR) models and tested them in high-throughput screening (HTS). Here, the novel positive allosteric modulator (PAM) VU0506013 was identified, which exhibits nanomolar affinity and pronounced selectivity toward the Y₄R in engineered cell lines and mouse descending colon mucosa natively expressing the Y₄R. Based on this lead structure, we conducted a systematic SAR study in two regions of the scaffold and presented a series of 27 analogues with modifications in the N- and C-terminal heterocycles of the molecule to obtain insight into functionally relevant positions. By mutagenesis and computational docking, we present a potential binding mode of VU0506013 in the transmembrane core of the Y₄R. VU0506013 presents a promising scaffold for developing in vivo tools to move toward anti-obesity drug research focused on the Y₄R.

Advances in the allostery of angiotensin II type 1 receptor

Angiotensin II type 1 receptor (AT1R) is a promising therapeutic target for cardiovascular diseases. Compared with orthosteric ligands, allosteric modulators attract considerable attention for drug development due to their unique advantages of high selectivity and safety. However, no allosteric modulators of AT1R have been applied in clinical trials up to now. Except for the classical allosteric modulators of AT1R such as antibody, peptides and amino acids, cholesterol and biased allosteric modulators, there are non-classical allosteric modes including the ligand-independent allosteric mode, and allosteric mode of biased agonists and dimers. In addition, finding the allosteric pockets based on AT1R conformational change and interaction interface of dimers are the future of drug design. In this review, we summarize the different allosteric mode of AT1R, with a view to contribute to the development and utilization of drugs targeting AT1R allostery.

Structural insights into angiotensin receptor signaling modulation by balanced and biased agonists

The peptide hormone angiotensin II regulates blood pressure mainly through the type 1 angiotensin II receptor AT1 R and its downstream signaling proteins Gq and β-arrestin. AT1 R blockers, clinically used as antihypertensive drugs, inhibit both signaling pathways, whereas AT1 R β-arrestin-biased agonists have shown great potential for the treatment of acute heart failure. Here, we present a cryo-electron microscopy (cryo-EM) structure of the human AT1 R in complex with a balanced agonist, Sar1 -AngII, and Gq protein at 2.9 Å resolution. This structure, together with extensive functional assays and computational modeling, reveals the molecular mechanisms for AT1 R signaling modulation and suggests that a major hydrogen bond network (MHN) inside the receptor serves as a key regulator of AT1 R signal transduction from the ligand-binding pocket to both Gq and β-arrestin pathways. Specifically, we found that the MHN mutations N1113.35 A and N2947.45 A induce biased signaling to Gq and β-arrestin, respectively. These insights should facilitate AT1 R structure-based drug discovery for the treatment of cardiovascular diseases.

Application of Mixed-Solvent Molecular Dynamics Simulations for Prediction of Allosteric Sites on G Protein-Coupled Receptors

The development of small molecule allosteric modulators acting at G protein-coupled receptors (GPCRs) is becoming increasingly attractive. Such compounds have advantages over traditional drugs acting at orthosteric sites on these receptors, in particular target specificity. However, the number and locations of druggable allosteric sites within most clinically relevant GPCRs are unknown. In the present study, we describe the development and application of a mixed-solvent molecular dynamics (MixMD)-based method for the identification of allosteric sites on GPCRs. The method employs small organic probes with druglike qualities to identify druggable hotspots in multiple replicate short-timescale simulations. As proof of principle, we first applied the method retrospectively to a test set of five GPCRs (cannabinoid receptor type 1, C-C chemokine receptor type 2, M2 muscarinic receptor, P2Y purinoceptor 1, and protease-activated receptor 2) with known allosteric sites in diverse locations. This resulted in the identification of the known allosteric sites on these receptors. We then applied the method to the μ-opioid receptor. Several allosteric modulators for this receptor are known, although the binding sites for these modulators are not known. The MixMD-based method revealed several potential allosteric sites on the mu-opioid receptor. Implementation of the MixMD-based method should aid future efforts in the structure-based drug design of drugs targeting allosteric sites on GPCRs. SIGNIFICANCE STATEMENT: Allosteric modulation of G protein-coupled receptors (GPCRs) has the potential to provide more selective drugs. However, there are limited structures of GPCRs bound to allosteric modulators, and obtaining such structures is problematic. Current computational methods utilize static structures and therefore may not identify hidden or cryptic sites. Here we describe the use of small organic probes and molecular dynamics to identify druggable allosteric hotspots on GPCRs. The results reinforce the importance of protein dynamics in allosteric site identification.

Olfaction: Allosteric modulation

New research indicates that the odor-evoked responses of human olfactory receptors can be enhanced via the non-competitive binding of an allosteric modulator. This modulatory mechanism adds an additional layer of complexity to the peripheral encoding of odors.

Structural Understanding of Peptide-Bound G Protein-Coupled Receptors: Peptide-Target Interactions

The activation of G protein-coupled receptors (GPCRs) is triggered by ligand binding to their orthosteric sites, which induces ligand-specific conformational changes. Agonists and antagonists bound to GPCR orthosteric sites provide detailed information on ligand-binding modes. Among these, peptide ligands play an instrumental role in GPCR pharmacology and have attracted increased attention as therapeutic drugs. The recent breakthrough in GPCR structural biology has resulted in the remarkable availability of peptide-bound GPCR complexes. Despite the several structural similarities shared by these receptors, they exhibit distinct features in terms of peptide recognition and receptor activation. From this perspective, we have summarized the current status of peptide-bound GPCR structural complexes, largely focusing on the interactions between the receptor and its peptide ligand at the orthosteric site. In-depth structural investigations have yielded valuable insights into the molecular mechanisms underlying peptide recognition. This study would contribute to the discovery of GPCR peptide drugs with improved therapeutic effects.

Recent applications of computational methods to allosteric drug discovery

Interest in exploiting allosteric sites for the development of new therapeutics has grown considerably over the last two decades. The chief driving force behind the interest in allostery for drug discovery stems from the fact that in comparison to orthosteric sites, allosteric sites are less conserved across a protein family, thereby offering greater opportunity for selectivity and ultimately tolerability. While there is significant overlap between structure-based drug design for orthosteric and allosteric sites, allosteric sites offer additional challenges mostly involving the need to better understand protein flexibility and its relationship to protein function. Here we examine the extent to which structure-based drug design is impacting allosteric drug design by highlighting several targets across a variety of target classes.

AlphaFold2 versus experimental structures: evaluation on G protein-coupled receptors

As important drug targets, G protein-coupled receptors (GPCRs) play pivotal roles in a wide range of physiological processes. Extensive efforts of structural biology have been made on the study of GPCRs. However, a large portion of GPCR structures remain unsolved due to structural instability. Recently, AlphaFold2 has been developed to predict structure models of many functionally important proteins including all members of the GPCR family. Herein we evaluated the accuracy of GPCR structure models predicted by AlphaFold2. We revealed that AlphaFold2 could capture the overall backbone features of the receptors. However, the predicted models and experimental structures were different in many aspects including the assembly of the extracellular and transmembrane domains, the shape of the ligand-binding pockets, and the conformation of the transducer-binding interfaces. These differences impeded the use of predicted structure models in the functional study and structure-based drug design of GPCRs, which required reliable high-resolution structural information.

Discovery of small molecule agonists of the Relaxin Family Peptide Receptor 2

The relaxin/insulin-like family peptide receptor 2 (RXFP2) belongs to the family of class A G-protein coupled receptors (GPCRs) and it is the only known target for the insulin-like factor 3 peptide (INSL3). The importance of this ligand-receptor pair in the development of the gubernacular ligament during the transabdominal phase of testicular descent is well established. More recently, RXFP2 has been implicated in maintaining healthy bone formation. In this report, we describe the discovery of a small molecule series of RXFP2 agonists. These compounds are highly potent, efficacious, and selective RXFP2 allosteric agonists that induce gubernacular invagination in mouse embryos, increase mineralization activity in human osteoblasts in vitro, and improve bone trabecular parameters in adult mice. The described RXFP2 agonists are orally bioavailable and display favorable pharmacokinetic properties, which allow for future evaluation of the therapeutic benefits of modulating RXFP2 activation in disease models.

Specific small molecule RXFP2 agonists with favorable pharmacokinetic properties induce gubernacular invagination in mouse embryos, increase mineralization activity in human osteoblasts in vitro, and improve bone trabecular parameters in adult mice.

Membrane Lipids Are an Integral Part of Transmembrane Allosteric Sites in GPCRs: A Case Study of Cannabinoid CB1 Receptor Bound to a Negative Allosteric Modulator, ORG27569, and Analogs

A growing number of G-protein-coupled receptor (GPCR) structures reveal novel transmembrane lipid-exposed allosteric sites. Ligands must first partition into the surrounding membrane and take lipid paths to these sites. Remarkably, a significant part of the bound ligands appears exposed to the membrane lipids. The experimental structures do not usually account for the surrounding lipids, and their apparent contribution to ligand access and binding is often overlooked and poorly understood. Using classical and enhanced molecular dynamics simulations, we show that membrane lipids are critical in the access and binding of ORG27569 and its analogs at the transmembrane site of cannabinoid CB1 receptor. The observed differences in the binding affinity and cooperativity arise from the functional groups that interact primarily with lipids. Our results demonstrate the significance of incorporating membrane lipids as an integral component of transmembrane sites for accurate characterization, binding-affinity calculations, and lead optimization in drug discovery.

Mechanistic Understanding of the Palmitoylation of Go Protein in the Allosteric Regulation of Adhesion Receptor GPR97

Adhesion G-protein-coupled receptors (aGPCRs)—a major family of GPCRs—play critical roles in the regulation of tissue development and cancer progression. The orphan receptor GPR97, activated by glucocorticoid stress hormones, is a prototypical aGPCR. Although it has been established that the palmitoylation of the C-terminal G_o protein is essential for G_o’s efficient engagement with the active GPR97, the detailed allosteric mechanism remains to be clarified. Hence, we performed extensive large-scale molecular dynamics (MD) simulations of the GPR97−G_o complex in the presence or absence of G_o palmitoylation. The conformational landscapes analyzed by Markov state models revealed that the overall conformation of GPR97 is preferred to be fully active when interacting with palmitoylated G_o protein. Structural and energetic analyses indicated that the palmitoylation of G_o can allosterically stabilize the critical residues in the ligand-binding pocket of GPR97 and increase the affinity of the ligand for GPR97. Furthermore, the community network analysis suggests that the palmitoylation of G_o not only allosterically strengthens the internal interactions between G_αo and G_βγ, but also enhances the coupling between G_o and GPR97. Our study provides mechanistic insights into the regulation of aGPCRs via post-translational modifications of the G_o protein, and offers guidance for future drug design of aGPCRs.

Characterization of Dual-Acting A3 Adenosine Receptor Positive Allosteric Modulators That Preferentially Enhance Adenosine-Induced Gαi3 and GαoA Isoprotein Activation

The A₃ adenosine receptor (A₃AR) is a promising therapeutic target for inflammatory diseases, cancer, and chronic neuropathic pain, with agonists already in advanced clinical trials. Here we report an in-depth comparison of the pharmacological properties and structure-activity relationships of existing and expanded compound libraries of 2-substituted 1H-imidazo[4,5-c]quinolin-4-amine and 4-amino-substituted quinoline derivatives that function as A₃AR positive allosteric modulators (PAMs). We also show that our lead compound from each series enhances adenosine-induced A₃AR signaling preferentially toward activation of Gα_i3 and Gα_oA isoproteins, which are coexpressed with the A₃AR in immune cells and spinal cord neurons. Finally, utilizing an extracellular/intracellular chimeric A₃AR approach composed of sequences from a responding (human) and a nonresponding (mouse) species, we provide evidence in support of the idea that the imidazoquinolin-4-amine class of PAMs variably interacts dually with the orthosteric ligand binding site as well as with a separate allosteric site located within the inner/intracellular regions of the receptor. This study has advanced both structural and pharmacological understanding of these two classes of A₃AR PAMs, which includes leads for future pharmaceutical development.

The Role of Conformational Dynamics and Allostery in the Control of Distinct Efficacies of Agonists to the Glucocorticoid Receptor

Glucocorticoid receptor (GR) regulates various cellular functions. Given its broad influence on metabolic activities, it has been the target of drug discovery for decades. However, how drugs induce conformational changes in GR has remained elusive. Herein, we used five GR agonists (dex, AZ938, pred, cor, and dibC) with different efficacies to investigate which aspect of the ligand induced the differences in efficacy. We performed molecular dynamics simulations on the five systems (dex-, AZ938-, pred-, cor-, and dibC-bound systems) and observed a distinct discrepancy in the conformation of the cofactor TIF2. Moreover, we discovered ligand-induced differences regarding the level of conformational changes posed by the binding of cofactor TIF2 and identified a pair of essential residues D590 and T39. We further found a positive correlation between the efficacies of ligands and the interaction of the two binding pockets’ domains, where D590 and T739 were involved, implying their significance in the participation of allosteric communication. Using community network analysis, two essential communities containing D590 and T739 were identified with their connectivity correlating to the efficacy of ligands. The potential communication pathways between these two residues were revealed. These results revealed the underlying mechanism of allosteric communication between the ligand-binding and cofactor-binding pockets and identified a pair of important residues in the allosteric communication pathway, which can serve as a guide for future drug discovery.

Insights into the Allosteric Effect of SENP1 Q597A Mutation on the Hydrolytic Reaction of SUMO1 via an Integrated Computational Study

Small ubiquitin-related modifier (SUMO)-specific protease 1 (SENP1) is a cysteine protease that catalyzes the cleavage of the C-terminus of SUMO1 for the processing of SUMO precursors and deSUMOylation of target proteins. SENP1 is considered to be a promising target for the treatment of hepatocellular carcinoma (HCC) and prostate cancer. SENP1 Gln597 is located at the unstructured loop connecting the helices α4 to α5. The Q597A mutation of SENP1 allosterically disrupts the hydrolytic reaction of SUMO1 through an unknown mechanism. Here, extensive multiple replicates of microsecond molecular dynamics (MD) simulations, coupled with principal component analysis, dynamic cross-correlation analysis, community network analysis, and binding free energy calculations, were performed to elucidate the detailed mechanism. Our MD simulations showed that the Q597A mutation induced marked dynamic conformational changes in SENP1, especially in the unstructured loop connecting the helices α4 to α5 which the mutation site occupies. Moreover, the Q597A mutation caused conformational changes to catalytic Cys603 and His533 at the active site, which might impair the catalytic activity of SENP1 in processing SUMO1. Moreover, binding free energy calculations revealed that the Q597A mutation had a minor effect on the binding affinity of SUMO1 to SENP1. Together, these results may broaden our understanding of the allosteric modulation of the SENP1−SUMO1 complex.

The N Terminus of Adhesion G Protein–Coupled Receptor GPR126/ADGRG6 as Allosteric Force Integrator

The adhesion G protein–coupled receptor (aGPCR) GPR126/ADGRG6 plays an important role in several physiological functions, such as myelination or peripheral nerve repair. This renders the receptor an attractive pharmacological target. GPR126 is a mechano-sensor that translates the binding of extracellular matrix (ECM) molecules to its N terminus into a metabotropic intracellular signal. To date, the structural requirements and the character of the forces needed for this ECM-mediated receptor activation are largely unknown. In this study, we provide this information by combining classic second-messenger detection with single-cell atomic force microscopy. We established a monoclonal antibody targeting the N terminus to stimulate GPR126 and compared it to the activation through its known ECM ligands, collagen IV and laminin 211. As each ligand uses a distinct mode of action, the N terminus can be regarded as an allosteric module that can fine-tune receptor activation in a context-specific manner.

Delineating the conformational landscape and intrinsic properties of the angiotensin II type 2 receptor using a computational study

As a key regulator for the renin-angiotensin system, a class A G protein-coupled receptor (GPCR), AngII type 2 receptor (AT₂R), plays a pivotal role in the homeostasis of the cardiovascular system. Compared with other GPCRs, AT₂R has a unique antagonist-bound conformation and its mechanism is still an enigma. Here, we applied combined dynamic and evolutional approaches to investigate the conformational space and intrinsic properties of AT₂R. With molecular dynamic simulations, Markov State Models, and statistics coupled analysis, we captured the conformational landscape of AT₂R and identified its uniquity from both dynamical and evolutional viewpoints. A cryptic pocket was also discovered in the intermediate state during conformation transitions. These findings offer a deeper understanding of the AT₂R mechanism at an atomic level and provide hints for the design of novel AT₂R modulators.

Understanding the P-Loop Conformation in the Determination of Inhibitor Selectivity Toward the Hepatocellular Carcinoma-Associated Dark Kinase STK17B

As a member of the death-associated protein kinase family of serine/threonine kinases, the STK17B has been associated with diverse diseases such as hepatocellular carcinoma. However, the conformational dynamics of the phosphate-binding loop (P-loop) in the determination of inhibitor selectivity profile to the STK17B are less understood. Here, a multi-microsecond length molecular dynamics (MD) simulation of STK17B in the three different states (ligand-free, ADP-bound, and ligand-bound states) was carried out to uncover the conformational plasticity of the P-loop. Together with the analyses of principal component analysis, cross-correlation and generalized correlation motions, secondary structural analysis, and community network analysis, the conformational dynamics of the P-loop in the different states were revealed, in which the P-loop flipped into the ADP-binding site upon the inhibitor binding and interacted with the inhibitor and the C-lobe, strengthened the communication between the N- and C-lobes. These resulting interactions contributed to inhibitor selectivity profile to the STK17B. Our results may advance our understanding of kinase inhibitor selectivity and offer possible implications for the design of highly selective inhibitors for other protein kinases.

Targeting key proteins involved in transcriptional regulation for cancer therapy: Current strategies and future prospective

The key proteins involved in transcriptional regulation play convergent roles in cellular homeostasis, and their dysfunction mediates aberrant gene expressions that underline the hallmarks of tumorigenesis. As tumor progression is dependent on such abnormal regulation of transcription, it is important to discover novel chemical entities as antitumor drugs that target key tumor-associated proteins involved in transcriptional regulation. Despite most key proteins (especially transcription factors) involved in transcriptional regulation are historically recognized as undruggable targets, multiple targeting approaches at diverse levels of transcriptional regulation, such as epigenetic intervention, inhibition of DNA-binding of transcriptional factors, and inhibition of the protein-protein interactions (PPIs), have been established in preclinically or clinically studies. In addition, several new approaches have recently been described, such as targeting proteasomal degradation and eliciting synthetic lethality. This review will emphasize on accentuating these developing therapeutic approaches and provide a thorough conspectus of the drug development to target key proteins involved in transcriptional regulation and their impact on future oncotherapy.

Mechanistic Insights Into Co-Administration of Allosteric and Orthosteric Drugs to Overcome Drug-Resistance in T315I BCR-ABL1

Chronic myeloid leukemia (CML) is a myeloproliferative neoplasm, driven by the BCR-ABL1 fusion oncoprotein. The discovery of orthosteric BCR-ABL1 tyrosine kinase inhibitors (TKIs) targeting its active ATP-binding pocket, such as first-generation Imatinib and second-generation Nilotinib (NIL), has profoundly revolutionized the therapeutic landscape of CML. However, currently targeted therapeutics still face considerable challenges with the inevitable emergence of drug-resistant mutations within BCR-ABL1. One of the most common resistant mutations in BCR-ABL1 is the T315I gatekeeper mutation, which confers resistance to most current TKIs in use. To resolve such conundrum, co-administration of orthosteric TKIs and allosteric drugs offers a novel paradigm to tackle drug resistance. Remarkably, previous studies have confirmed that the dual targeting BCR-ABL1 utilizing orthosteric TKI NIL and allosteric inhibitor ABL001 resulted in eradication of the CML xenograft tumors, exhibiting promising therapeutic potential. Previous studies have demonstrated the cooperated mechanism of two drugs. However, the conformational landscapes of synergistic effects remain unclear, hampering future efforts in optimizations and improvements. Hence, extensive large-scale molecular dynamics (MD) simulations of wide type (WT), WT-NIL, T315I, T315I-NIL, T315I-ABL001 and T315I-ABL001-NIL systems were carried out in an attempt to address such question. Simulation data revealed that the dynamic landscape of NIL-bound BCR-ABL1 was significantly reshaped upon ABL001 binding, as it shifted from an active conformation towards an inactive conformation. The community network of allosteric signaling was analyzed to elucidate the atomistic overview of allosteric regulation within BCR-ABL1. Moreover, binding free energy analysis unveiled that the affinity of NIL to BCR-ABL1 increased by the induction of ABL001, which led to its favorable binding and the release of drug resistance. The findings uncovered the in-depth structural mechanisms underpinning dual-targeting towards T315I BCR-ABL1 to overcome its drug resistance and will offer guidance for the rational design of next generations of BCR-ABL1 modulators and future combinatory therapeutic regimens.

Elucidation of the conformational dynamics and assembly of Argonaute−RNA complexes by distinct yet coordinated actions of the supplementary microRNA

Argonaute (AGO) proteins, the core of RNA-induced silencing complex, are guided by microRNAs (miRNAs) to recognize target RNA for repression. The miRNA–target RNA recognition forms initially through pairing at the seed region while the additional supplementary pairing can enhance target recognition and compensate for seed mismatch. The extension of miRNA lengths can strengthen the target affinity when pairing both in the seed and supplementary regions. However, the mechanism underlying the effect of the supplementary pairing on the conformational dynamics and the assembly of AGO–RNA complex remains poorly understood. To address this, we performed large-scale molecular dynamics simulations of AGO–RNA complexes with different pairing patterns and miRNA lengths. The results reveal that the additional supplementary pairing can not only strengthen the interaction between miRNA and target RNA, but also induce the increased plasticity of the PAZ domain and enhance the domain connectivity among the PAZ, PIWI, N domains of the AGO protein. The strong community network between these domains tightens the mouth of the supplementary chamber of AGO protein, which prevents the escape of target RNA from the complex and shields it from solvent water attack. Importantly, the inner stronger matching pairs between the miRNA and target RNA can compensate for weaker mismatches at the edge of supplementary region. These findings provide guidance for the design of miRNA mimics and anti-miRNAs for both clinical and experimental use and open the way for further engineering of AGO proteins as a new tool in the field of gene regulation.

Allosteric modulation of dopamine D2L receptor in complex with Gi1 and Gi2 proteins: the effect of subtle structural and stereochemical ligand modifications

Background

Allosteric modulation of G protein-coupled receptors (GPCRs) is nowadays one of the hot topics in drug discovery. In particular, allosteric modulators of D₂ receptor have been proposed as potential modern therapeutics to treat schizophrenia and Parkinson’s disease.

Methods

To address some subtle structural and stereochemical aspects of allosteric modulation of D₂ receptor, we performed extensive in silico studies of both enantiomers of two compounds (compound 1 and compound 2), and one of them (compound 2) was synthesized as a racemate in-house and studied in vitro.

Results

Our molecular dynamics simulations confirmed literature reports that the R enantiomer of compound 1 is a positive allosteric modulator of the D_2L receptor, while its S enantiomer is a negative allosteric modulator. Moreover, based on the principal component analysis (PCA), we hypothesized that both enantiomers of compound 2 behave as silent allosteric modulators, in line with our in vitro studies. PCA calculations suggest that the most pronounced modulator-induced receptor rearrangements occur at the transmembrane helix 7 (TM7). In particular, TM7 bending at the conserved P7.50 and G7.42 was observed. The latter resides next to the Y7.43, which is a significant part of the orthosteric binding site. Moreover, the W7.40 conformation seems to be affected by the presence of the positive allosteric modulator.

Conclusions

Our work reveals that allosteric modulation of the D_2L receptor can be affected by subtle ligand modifications. A change in configuration of a chiral carbon and/or minor structural modulator modifications are solely responsible for the functional outcome of the allosteric modulator.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s43440-021-00352-x.

Delineating the activation mechanism and conformational landscape of a class B G protein-coupled receptor glucagon receptor

Class B G protein-coupled receptors (GPCRs) are important targets in the treatment of metabolic syndrome and diabetes. Although multiple structures of class B GPCRs-G protein complexes have been elucidated, the detailed activation mechanism of the receptors remains unclear. Here, we combine Gaussian accelerated molecular dynamics simulations and Markov state models (MSM) to investigate the activation mechanism of a canonical class B GPCR, human glucagon receptor-GCGR, including the negative allosteric modulator-bound inactive state, the agonist glucagon-bound active state, and both glucagon- and Gs-bound fully active state. The free-energy landscapes of GCGR show the conformational ensemble consisting of three activation-associated states: inactive, active, and fully active. The structural analysis indicates the high dynamics of GCGR upon glucagon binding with both active and inactive conformations in the ensemble. Significantly, the H8 and TM6 exhibits distinct features from the inactive to the active states. The additional simulations demonstrate the role of H8 in the recruitment of Gs. Gs binding presents a crucial function of stabilizing the glucagon binding site and MSM highlights the absolute requirement of Gs to help the GCGR reach the fully active state. Together, our results reveal the detailed activation mechanism of GCGR from the view of conformational dynamics.

Probing GPCR Dimerization Using Peptides

G protein-coupled receptors (GPCRs) are the largest class of membrane proteins and the most common and extensively studied pharmacological target. Numerous studies over the last decade have confirmed that GPCRs do not only exist and function in their monomeric form but in fact, have the ability to form dimers or higher order oligomers with other GPCRs, as well as other classes of receptors. GPCR oligomers have become increasingly attractive to investigate as they have the ability to modulate the pharmacological responses of the receptors which in turn, could have important functional roles in diseases, such as cancer and several neurological & neuropsychiatric disorders. Despite the growing evidence in the field of GPCR oligomerisation, the lack of structural information, as well as targeting the 'undruggable' protein-protein interactions (PPIs) involved in these complexes, has presented difficulties. Outside the field of GPCRs, targeting PPIs has been widely studied, with a variety of techniques being investigated; from small-molecule inhibitors to disrupting peptides. In this review, we will demonstrate several physiologically relevant GPCR dimers and discuss an array of strategies and techniques that can be employed when targeting these complexes, as well as provide ideas for future development.

CDK9 Inhibitors in Cancer Research

Cyclin dependent kinase 9 (CDK9) played an essential role in regulating transcriptional elongation. Aberrations in CDK9 activity have been observed in various cancers, which made CDK9 was an attractive therapeutic...

Deciphering the Mechanism of Gilteritinib Overcoming Lorlatinib Resistance to the Double Mutant I1171N/F1174I in Anaplastic Lymphoma Kinase

Anaplastic lymphoma kinase (ALK) is validated as a therapeutic molecular target in multiple malignancies, such as non-small cell lung cancer (NSCLC). However, the feasibility of targeted therapies exerted by ALK inhibitors is inevitably hindered owing to drug resistance. The emergence of clinically acquired drug mutations has become a major challenge to targeted therapies and personalized medicines. Thus, elucidating the mechanism of resistance to ALK inhibitors is helpful for providing new therapeutic strategies for the design of next-generation drug. Here, we used molecular docking and multiple molecular dynamics simulations combined with correlated and energetical analyses to explore the mechanism of how gilteritinib overcomes lorlatinib resistance to the double mutant ALK I1171N/F1174I. We found that the conformational dynamics of the ALK kinase domain was reduced by the double mutations I1171N/F1174I. Moreover, energetical and structural analyses implied that the double mutations largely disturbed the conserved hydrogen bonding interactions from the hinge residues Glu1197 and Met1199 in the lorlatinib-bound state, whereas they had no discernible adverse impact on the binding affinity and stability of gilteritinib-bound state. These discrepancies created the capacity of the double mutant ALK I1171N/F1174I to confer drug resistance to lorlatinib. Our result anticipates to provide a mechanistic insight into the mechanism of drug resistance induced by ALK I1171N/F1174I that are resistant to lorlatinib treatment in NSCLC.

Harnessing Reversed Allosteric Communication: A Novel Strategy for Allosteric Drug Discovery

Allostery is a fundamental and extensive mechanism of intramolecular signal transmission. Allosteric drugs possess several unique pharmacological advantages over traditional orthosteric drugs, including greater selectivity, better physicochemical properties, and lower off-target toxicity. However, owing to the complexity of allosteric regulation, experimental approaches for the development of allosteric modulators are traditionally serendipitous. Recently, the reversed allosteric communication theory has been proposed, providing a feasible tool for the unbiased detection of allosteric sites. Herein, we review the latest research on the reversed allosteric communication effect using the examples of sirtuin 6, epidermal growth factor receptor, 3-phosphoinositide-dependent protein kinase 1, and Related to A and C kinases (RAC) serine/threonine protein kinase B and recapitulate the methodologies of reversed allosteric communication strategy. The novel reversed allosteric communication strategy greatly expands the horizon of allosteric site identification and allosteric mechanism exploration and is expected to accelerate an end-to-end framework for drug discovery.

Structure-Based Virtual Screening for Ligands of G Protein-Coupled Receptors: What Can Molecular Docking Do for You?

G protein-coupled receptors (GPCRs) constitute the largest family of membrane proteins in the human genome and are important therapeutic targets. During the last decade, the number of atomic-resolution structures of GPCRs has increased rapidly, providing insights into drug binding at the molecular level. These breakthroughs have created excitement regarding the potential of using structural information in ligand design and initiated a new era of rational drug discovery for GPCRs. The molecular docking method is now widely applied to model the three-dimensional structures of GPCR-ligand complexes and screen for chemical probes in large compound libraries. In this review article, we first summarize the current structural coverage of the GPCR superfamily and the understanding of receptor-ligand interactions at atomic resolution. We then present the general workflow of structure-based virtual screening and strategies to discover GPCR ligands in chemical libraries. We assess the state of the art of this research field by summarizing prospective applications of virtual screening based on experimental structures. Strategies to identify compounds with specific efficacy and selectivity profiles are discussed, illustrating the opportunities and limitations of the molecular docking method. Our overview shows that structure-based virtual screening can discover novel leads and will be essential in pursuing the next generation of GPCR drugs. SIGNIFICANCE STATEMENT: Extraordinary advances in the structural biology of G protein-coupled receptors have revealed the molecular details of ligand recognition by this large family of therapeutic targets, providing novel avenues for rational drug design. Structure-based docking is an efficient computational approach to identify novel chemical probes from large compound libraries, which has the potential to accelerate the development of drug candidates.

Characterisation of bis(4-aminoquinoline)s as α1A adrenoceptor allosteric modulators

The development of sub-type selective α₁ adrenoceptor ligands has been hampered by the high sequence similarity of the amino acids forming the orthosteric binding pocket of the three α₁ adrenoceptor subtypes, along with other biogenic amine receptors. One possible approach to overcome this issue is to target allosteric sites on the α₁ adrenoceptors. Previous docking studies suggested that one of the quinoline moieties of a bis(4-aminoquinoline), comprising a 9-carbon methylene linker attached via the amine groups, could interact with residues outside of the orthosteric binding site while, simultaneously, the other quinoline moiety bound within the orthosteric site. We therefore hypothesized that this compound could act in a bitopic manner, displaying both orthosteric and allosteric binding properties. To test this proposition, we investigated the allosteric activity of a series of bis(4-aminoquinoline)s with linker lengths ranging from 2 to 12 methylene units (designated C2-C12). A linear trend of increasing [³H]prazosin dissociation rate with increasing linker length between C7 and C11 was observed, confirming their action as allosteric modulators. These data suggest that the optimal linker length for the bis(4-aminoquinoline)s to occupy the allosteric site of the α_1A adrenoceptor is between 7 and 11 methylene units. In addition, the ability of C9 bis(4-aminoquinoline) to modulate the activation of the α_1A adrenoceptor by norepinephrine was subsequently examined, showing that C9 acted as a non-competitive antagonist. Our findings indicate that the bis(4-aminoquinolines) are acting as allosteric modulators of orthosteric ligand binding, but not efficacy, in a bitopic manner.

Atomic-scale insights into allosteric inhibition and evolutional rescue mechanism of Streptococcus thermophilus Cas9 by the anti-CRISPR protein AcrIIA6

CRISPR-Cas systems are prokaryotic adaptive immunity against invading phages and plasmids. Phages have evolved diverse protein inhibitors of CRISPR-Cas systems, called anti-CRISPR (Acr) proteins, to neutralize this CRISPR machinery. In response, bacteria have co-evolved Cas variants to escape phage's anti-CRISPR strategies, called anti-anti-CRISPR systems. Here we explore the anti-CRISPR allosteric inhibition and anti-anti-CRISPR rescue mechanisms between Streptococcus thermophilus Cas9 (St1Cas9) and the anti-CRISPR protein AcrIIA6 at the atomic level, by generating mutants of key residues in St1Cas9. Extensive unbiased molecular dynamics simulations show that the functional motions of St1Cas9 in the presence of AcrIIA6 differ substantially from those of St1Cas9 alone. AcrIIA6 binding triggers a shift of St1Cas9 conformational ensemble towards a less catalytically competent state; this state significantly compromises protospacer adjacent motif (PAM) recognition and nuclease activity by altering interdependently conformational dynamics and allosteric signals among nuclease domains, PAM-interacting (PI) regions, and AcrIIA6 binding motifs. Via in vitro DNA cleavage assays, we further elucidate the rescue mechanism of efficiently escaping AcrIIA6 inhibition harboring St1Cas9 triple mutations (G993K/K1008M/K1010E) in the PI domain and identify the evolutionary landscape of such mutational escape within species. Our results provide mechanistic insights into Acr proteins as natural brakes for the CRISPR-Cas systems and a promising potential for the design of allosteric Acr peptidomimetics.

How does nintedanib overcome cancer drug-resistant mutation of RET protein-tyrosine kinase: insights from molecular dynamics simulations

Targeted drug therapies represent a therapeutic breakthrough in the treatment of human cancer. However, the emergence of acquired resistance inevitably compromises therapeutic drugs. Rearranged during transfection (RET) proto-oncogene, which encodes a receptor tyrosine kinase, is a target for several kinds of human cancer such as thyroid, breast, and colorectal carcinoma. A single mutation L881V at the RET kinase domain was found in familial medullary thyroid carcinoma. Nintedanib can effectively inhibit the RET L881V mutant, whereas its analog compound 1 is unable to combat this mutant. However, the underlying mechanism was still unexplored. Here, molecular dynamics (MD) simulations, binding free energy calculations, and structural analysis were performed to uncover the mechanism of overcoming the resistance of RET L881V mutant to nintedanib. Energetic analysis revealed that the L881V mutant remained sensitive to the treatment of nintedanib, whereas it was insensitive to the compound 1. Structural analysis further showed that the distribution of K758, D892, and N879 network had a detrimental effect on the binding of compound 1 to the L881V mutant. The obtained results may provide insight into the mechanism of overcoming resistance in the RET kinase.

Allosteric Modulator Leads Hiding in Plain Site: Developing Peptide and Peptidomimetics as GPCR Allosteric Modulators

Allosteric modulators (AMs) of G-protein coupled receptors (GPCRs) are desirable drug targets because they can produce fewer on-target side effects, improved selectivity, and better biological specificity (e.g., biased signaling or probe dependence) than orthosteric drugs. An underappreciated source for identifying AM leads are peptides and proteins-many of which were evolutionarily selected as AMs-derived from endogenous protein-protein interactions (e.g., transducer/accessory proteins), intramolecular receptor contacts (e.g., pepducins or extracellular domains), endogenous peptides, and exogenous libraries (e.g., nanobodies or conotoxins). Peptides offer distinct advantages over small molecules, including high affinity, good tolerability, and good bioactivity, and specific disadvantages, including relatively poor metabolic stability and bioavailability. Peptidomimetics are molecules that combine the advantages of both peptides and small molecules by mimicking the peptide's chemical features responsible for bioactivity while improving its druggability. This review 1) discusses sources and strategies to identify peptide/peptidomimetic AMs, 2) overviews strategies to convert a peptide lead into more drug-like "peptidomimetic," and 3) critically analyzes the advantages, disadvantages, and future directions of peptidomimetic AMs. While small molecules will and should play a vital role in AM drug discovery, peptidomimetics can complement and even exceed the advantages of small molecules, depending on the target, site, lead, and associated factors.

Molecular Mechanism of Ca2+ in the Allosteric Regulation of Human Parathyroid Hormone Receptor-1

Parathyroid hormone (PTH) is an endogenous ligand that activates the PTH type 1 receptor (PTH1R) signaling. Ca²⁺, a common second messenger, acts as an allosteric regulator for prolonging the activation of PTH1R. However, a clear picture of the underlying allosteric mechanism is still missing. Herein, extensive molecular dynamics (MD) simulations are performed for PTH1R-PTH complexes with and without Ca²⁺ ions, allowing us to delineate the molecular details of calcium-induced allostery. Our results indicate that acidic residues in the extracellular loop 1 (ECL1) (D251, E252, E254, and E258-E260) and PTH (E19 and E22) serve as key determinants for local Ca²⁺-coupling structures and rigidity of ECL1. Moreover, the binding of Ca²⁺ induces conformational changes of transmembrane domain 6/7 (TM6/7) that are related to PTH1R activation and strengthens the residue-residue communication within PTH and TMD allosterically. Moreover, our results demonstrate that the presence of Ca²⁺ ions potentiates the interaction between PTH and PTH1R via steered molecular dynamics (SMD) simulations, while the point mutation in the PTH (PTH^R25C) weakens the binding of PTH and PTH1R. These results support that Ca²⁺ ions might further prolong the residence time of PTH on PTH1R and facilitate the positive allostery of PTH1R. Together, the present work provides new insights into the allosteric regulation mechanism of GPCRs induced by ions and related drug design targeting the PTH1R allosteric pathway.

Activation pathway of a G protein-coupled receptor uncovers conformational intermediates as targets for allosteric drug design

G protein-coupled receptors (GPCRs) are the most common proteins targeted by approved drugs. A complete mechanistic elucidation of large-scale conformational transitions underlying the activation mechanisms of GPCRs is of critical importance for therapeutic drug development. Here, we apply a combined computational and experimental framework integrating extensive molecular dynamics simulations, Markov state models, site-directed mutagenesis, and conformational biosensors to investigate the conformational landscape of the angiotensin II (AngII) type 1 receptor (AT1 receptor) - a prototypical class A GPCR-activation. Our findings suggest a synergistic transition mechanism for AT1 receptor activation. A key intermediate state is identified in the activation pathway, which possesses a cryptic binding site within the intracellular region of the receptor. Mutation of this cryptic site prevents activation of the downstream G protein signaling and β-arrestin-mediated pathways by the endogenous AngII octapeptide agonist, suggesting an allosteric regulatory mechanism. Together, these findings provide a deeper understanding of AT1 receptor activation at an atomic level and suggest avenues for the design of allosteric AT1 receptor modulators with a broad range of applications in GPCR biology, biophysics, and medicinal chemistry.

Allosteric Mechanism of Human Mitochondrial Phenylalanyl-tRNA Synthetase: An Atomistic MD Simulation and a Mutual Information-Based Network Study

Aminoacyl-tRNA synthetases (aaRSs), a family of ubiquitous and essential enzymes, can bind target tRNAs and catalyze the aminoacylation reaction in genetic code translation. In this work, we explore the dynamic properties and allosteric communication of human mitochondrial phenylalanyl-tRNA synthetase (hmPheRS) in free and bound states to understand the mechanisms of its tRNA^Phe recognition and allostery using molecular dynamics simulations combined with the torsional mutual information-based network model. Our results reveal that hmPheRS's residue mobility and inter-residue motional coupling are significantly enhanced by tRNA^Phe binding, and there occurs a strong allosteric communication which is critical for the aminoacylation reaction, suggesting the vital role of tRNA^Phe binding in the enzyme's function. The identified signaling pathways mainly make the connections between the anticodon binding domain (ABD) and catalytic domain (CAD), as well as within the CAD composed of many functional fragments and active sites, revealing the co-regulation role of them to act coordinately and achieve hmPheRS's aminoacylation function. Besides, several key residues along the communication pathways are identified to be involved in mediating the coordinated coupling between anticodon recognition at the ABD and activation process at the CAD, showing their pivotal role in the allosteric network, which are well consistent with the experimental observation. This study sheds light on the allosteric communication mechanism in hmPheRS and can provide important information for the structure-based drug design targeting aaRSs.

Overcoming barriers in photodynamic therapy harnessing nano-formulation strategies

Photodynamic therapy (PDT) has been extensively investigated for decades for tumor treatment because of its non-invasiveness, spatiotemporal selectivity, lower side-effects, and immune activation ability. It can be a promising treatment modality in several medical fields, including oncology, immunology, urology, dermatology, ophthalmology, cardiology, pneumology, and dentistry. Nevertheless, the clinical application of PDT is largely restricted by the drawbacks of traditional photosensitizers, limited tissue penetrability of light, inefficient induction of tumor cell death, tumor resistance to the therapy, and the severe pain induced by the therapy. Recently, various photosensitizer formulations and therapy strategies have been developed to overcome these barriers. Significantly, the introduction of nanomaterials in PDT, as carriers or photosensitizers, may overcome the drawbacks of traditional photosensitizers. Based on this, nanocomposites excited by various light sources are applied in the PDT of deep-seated tumors. Modulation of cell death pathways with co-delivered reagents promotes PDT induced tumor cell death. Relief of tumor resistance to PDT with combined therapy strategies further promotes tumor inhibition. Also, the optimization of photosensitizer formulations and therapy procedures reduces pain in PDT. Here, a systematic summary of recent advances in the fabrication of photosensitizers and the design of therapy strategies to overcome barriers in PDT is presented. Several aspects important for the clinical application of PDT in cancer treatment are also discussed.

Affinity Mass Spectrometry-Based Fragment Screening Identified a New Negative Allosteric Modulator of the Adenosine A2A Receptor Targeting the Sodium Ion Pocket

Allosteric ligands provide new opportunities to modulate G protein-coupled receptor (GPCR) function and present therapeutic benefits over orthosteric molecules. Negative allosteric modulators (NAMs) can inhibit the activation of a receptor and downstream signal transduction. Screening NAMs for a GPCR target is particularly challenging because of the difficulty in distinguishing NAMs from antagonists bound to the orthosteric site as they both show inhibitory effects in receptor signaling assays. Here we report an affinity mass spectrometry (MS)-based approach tailored to screening potential NAMs of a GPCR target especially from fragment libraries. Compared to regular surface plasmon resonance or NMR-based methods for fragment screening, our approach features a reduction of the protein and compound consumption by 2-4 orders of magnitude and an increase in the data acquisition speed by 2-3 orders of magnitude. Our affinity MS-based fragment screening led to the identification of a new NAM of the adenosine A_2A receptor (A_2AAR) bearing an unprecedented azetidine moiety predicted to occupy the allosteric sodium binding site. Molecular dynamics simulations, ligand structure-activity relationship (SAR) studies, and in-solution NMR analyses further revealed the unique binding mode and antagonistic property of this compound that differs considerably from HMA (5-(N,N-hexamethylene)amiloride), a well-characterized NAM of A_2AAR. Taken together, our work would facilitate fragment-based screening of allosteric modulators, as well as guide the design of novel NAMs acting at the sodium ion pocket of class A GPCRs.

Untangling Dual-Targeting Therapeutic Mechanism of Epidermal Growth Factor Receptor (EGFR) Based on Reversed Allosteric Communication

Dual-targeting therapeutics by coadministration of allosteric and orthosteric drugs is drawing increased attention as a revolutionary strategy for overcoming the drug-resistance problems. It was further observed that the occupation of orthosteric sites by therapeutics agents has the potential to enhance allosteric ligand binding, which leads to improved potency of allosteric drugs. Epidermal growth factor receptor (EGFR), as one of the most critical anti-cancer targets belonging to the receptor tyrosine kinase family, represents a quintessential example. It was revealed that osimertinib, an ATP-competitive covalent EGFR inhibitor, remarkably enhanced the affinity of a recently developed allosteric inhibitor JBJ-04-125-02 for EGFRL858R/T790M. Here, we utilized extensive large-scale molecular dynamics simulations and the reversed allosteric communication to untangle the detailed molecular underpinning, in which occupation of osimertinib at the orthosteric site altered the overall conformational ensemble of EGFR mutant and reshaped the allosteric site via long-distance signaling. A unique intermediate state resembling the active conformation was identified, which was further stabilized by osimertinib loading. Based on the allosteric communication pathway, we predicted a novel allosteric site positioned around K867, E868, H893, and K960 within the intermediate state. Its correlation with the orthosteric site was validated by both structural and energetic analysis, and its low sequence conservation indicated the potential for selective targeting across the human kinome. Together, these findings not only provided a mechanistic basis for future clinical application of the dual-targeting therapeutics, but also explored an innovative perception of allosteric inhibition of tyrosine kinase signaling.

Mechanism of allosteric activation of SIRT6 revealed by the action of rationally designed activators

The recent discovery of activator compounds binding to an allosteric site on the NAD⁺-dependent protein lysine deacetylase, sirtuin 6 (SIRT6) has attracted interest and presents a pharmaceutical target for aging-related and cancer diseases. However, the mechanism underlying allosteric activation of SIRT6 by the activator MDL-801 remains largely elusive because no major conformational changes are observed upon activator binding. By combining molecular dynamics simulations with biochemical and kinetic analyses of wild-type SIRT6 and its variant M136A, we show that conformational rotation of 2-methyl-4-fluoro-5-bromo substituent on the right phenyl ring (R-ring) of MDL-801, which uncovers previously unseen hydrophobic interactions, contributes to increased activating deacetylation activity of SIRT6. This hypothesis is further supported by the two newly synthesized MDL-801 derivatives through the removal of the 5-Br atom on the R-ring (MDL-801-D1) or the restraint of the rotation of the R-ring (MDL-801-D2). We further propose that the 5-Br atom serves as an allosteric driver that controls the ligand allosteric efficacy. Our study highlights the effect of allosteric enzyme catalytic activity by activator binding and provides a rational approach for enhancing deacetylation activity.

Evaluating the Efficiency of the Martini Force Field to Study Protein Dimerization in Aqueous and Membrane Environments

Protein-protein complex assembly is one of the major drivers of biological response. Understanding the mechanisms of protein oligomerization/dimerization would allow one to elucidate how these complexes participate in biological activities and could ultimately lead to new approaches in designing novel therapeutic agents. However, determining the exact association pathways and structures of such complexes remains a challenge. Here, we use parallel tempering metadynamics simulations in the well-tempered ensemble to evaluate the performance of Martini 2.2P and Martini open-beta 3 (Martini 3) force fields in reproducing the structure and energetics of the dimerization process of membrane proteins and proteins in an aqueous solution in reasonable accuracy and throughput. We find that Martini 2.2P systematically overestimates the free energy of association by estimating large barriers in distinct areas, which likely leads to overaggregation when multiple monomers are present. In comparison, the less viscous Martini 3 results in a systematic underestimation of the free energy of association for proteins in solution, while it performs well in describing the association of membrane proteins. In all cases, the near-native dimer complexes are identified as minima in the free energy surface albeit not always as the lowest minima. In the case of Martini 3, we find that the spurious supramolecular protein aggregation present in Martini 2.2P multimer simulations is alleviated and thus this force field may be more suitable for the study of protein oligomerization. We propose that the use of enhanced sampling simulations with a refined coarse-grained force field and appropriately defined collective variables is a robust approach for studying the protein dimerization process, although one should be cautious of the ranking of energy minima.