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

Decoding the selective mechanism behind a monobody inhibitor to the phosphatase domain of SHP2: insights from molecular dynamics simulations

The Src-homology 2 (SH2) domain-containing phosphatase 2 (SHP2), encoded by PTPN11, is a critical tyrosine phosphatase that regulates key cellular processes, including cell proliferation, survival, and migration. The catalytic activity of its protein tyrosine phosphatase (PTP) domain plays a pivotal role in cancer progression by activating oncogenic signaling pathways. In contrast, SHP1, another SH2 phosphatase encoded by PTPN6, generally functions as a tumor suppressor. Given their structural similarity yet distinct biological functions, developing selective SHP2-PTP inhibitors is crucial for targeted cancer therapy. Recently, a monobody, Mb (SHP2PTP_13) (Mb13), has been designed to bind to the SHP2-PTP structure specifically. However, the detailed mechanism involved in selective inhibition remains to be clarified. To achieve this objective, we conducted extensive molecular dynamics simulations of the Mb13-SHP2-PTP and Mb13-SHP1-PTP systems, together with multiple analyses, including cluster analysis, principal component analysis, free energy landscape evaluation, a cross-correlation matrix and binding free energy calculation. Our results demonstrated that Mb13 bound more stably to SHP2-PTP compared to SHP1-PTP. The SHP2 complex exhibited conformational stability and reduced flexibility, indicating a more substantial interaction. Detailed analysis revealed that key residues within SHP2-PTP formed more robust interactions with Mb13, enhancing the complex's overall stability. These findings suggested that the selective binding mechanism was primarily driven by specific stabilizing interactions at the molecular level. Overall, the enhanced understanding of SHP2-PTP's binding dynamics and stability offers valuable guidance for advancing drug design strategies targeting SHP2-mediated pathways.

Computational Elucidation of a Monobody Targeting the Phosphatase Domain of SHP2

Src homology 2 (SH2) domain-containing phosphatase 2 (SHP2) is a key regulator in cellular signaling pathways because its dysregulation has been implicated in various pathological conditions, including cancers and developmental disorders. Despite its importance, the molecular basis of SHP2's regulatory mechanism remains poorly understood, hindering the development of effective targeted therapies. In this study, we utilized the high-specificity monobody Mb11 to investigate its interaction with the SHP2 phosphatase domain (PTP) using multiple replica molecular dynamics simulations. Our analyses elucidate the precise mechanisms through which Mb11 achieves selective recognition and stabilization of the SHP2-PTP domain, identifying key residues and interaction networks essential for its high binding specificity and regulatory dynamics. Furthermore, the study highlights the pivotal role of residue C459 in preserving the structural integrity and functional coherence of the complex, acting as a central node within the interaction network and underpinning its stability and efficiency. These findings have significantly advanced the understanding of the mechanisms underlying SHP2's involvement in disease-related signaling and pathology while simultaneously paving the way for the rational design of targeted inhibitors, offering significant implications for therapeutic strategies in SHP2-associated diseases and contributing to the broader scope of precision medicine.

Targeting cryptic allosteric sites of G protein-coupled receptors as a novel strategy for biased drug discovery

G protein-coupled receptors (GPCRs) represent the largest family of membrane receptors and are highly effective targets for therapeutic drugs. GPCRs couple different downstream effectors, including G proteins (such as Gi/o, Gs, G12, and Gq) and β-arrestins (such as β-arrestin 1 and β-arrestin 2) to mediate diverse cellular and physiological responses. Biased signaling allows for the specific activation of certain pathways from the full range of receptors' signaling capabilities. Targeting more variable allosteric sites, which are spatially different from the highly conserved orthosteric sites, represents a novel approach in biased GPCR drug discovery, leading to innovative strategies for targeting GPCRs. Notably, the emergence of cryptic allosteric sites on GPCRs has expanded the repertoire of available drug targets and improved receptor subtype selectivity. Here, we conduct a summary of recent progress in the structural determination of cryptic allosteric sites on GPCRs and elucidate the biased signaling mechanisms induced by allosteric modulators. Additionally, we discuss means to identify cryptic allosteric sites and design biased allosteric modulators based on cryptic allosteric sites through structure-based drug design, which is an advanced pharmacotherapeutic approach for treating GPCR-associated diseases.

Exploring the constitutive activation mechanism of the class A orphan GPR20

GPR20, an orphan G protein-coupled receptor (GPCR), shows significant expression in intestinal tissue and represents a potential therapeutic target to treat gastrointestinal stromal tumors. GPR20 performs high constitutive activity when coupling with Gi. Despite the pharmacological importance of GPCR constitutive activation, determining the mechanism has long remained unclear. In this study, we explored the constitutive activation mechanism of GPR20 through large-scale unbiased molecular dynamics simulations. Our results unveil the allosteric nature of constitutively activated GPCR signal transduction involving extracellular and intracellular domains. Moreover, the constitutively active state of the GPR20 requires both the N-terminal cap and Gi protein. The N-terminal cap of GPR20 functions like an agonist and mediates long-range activated conformational shift. Together with the previous study, this study enhances our knowledge of the self-activation mechanism of the orphan receptor, facilitates the drug discovery efforts that target GPR20.

Targeting adhesion G protein-coupled receptors. Current status and future perspectives

G protein-coupled receptors (GPCRs) orchestrate many physiological functions and are a crucial target in drug discovery. Adhesion GPCRs (aGPCRs), the second largest family within this superfamily, are promising yet underexplored targets for treating various diseases, including obesity, psychiatric disorders, and cancer. However, the receptors' unique and complex structure and miscellaneous interactions complicate comprehensive pharmacological studies. Despite recent progress in determining structures and elucidation of the activation mechanism, the function of many receptors remains to be determined. This review consolidates current knowledge on aGPCR ligands, focusing on small molecule orthosteric ligands and allosteric modulators identified for the ADGRGs subfamily (subfamily VIII), (GPR56/ADGRG1, GPR64/ADGRG2, GPR97/ADGRG3, GPR114/ADGRG5, GPR126/ADGRG6, and GPR128/ADGRG7). We discuss challenges in hit identification, target validation, and drug discovery, highlighting molecular compositions and recent structural breakthroughs. ADGRG ligands can offer new insights into aGPCR modulation and have significant potential for novel therapeutic interventions targeting various diseases.

Exploring Protein S-Palmitoylation: Mechanisms, Detection, and Strategies for Inhibitor Discovery

S-palmitoylation is a reversible and dynamic process that involves the addition of long-chain fatty acids to proteins. This protein modification regulates various aspects of protein function, including subcellular localization, stability, conformation, and biomolecular interactions. The zinc finger DHHC (ZDHHC) domain-containing protein family is the main group of enzymes responsible for catalyzing protein S-palmitoylation, and 23 members have been identified in mammalian cells. Many proteins that undergo S-palmitoylation have been linked to disease pathogenesis and progression, suggesting that the development of effective inhibitors is a promising therapeutic strategy. Reducing the protein S-palmitoylation level can target either the PATs directly or their substrates. However, there are rare clinically effective S-palmitoylation inhibitors. This review aims to provide an overview of the S-palmitoylation field, including the catalytic mechanism of ZDHHC, S-palmitoylation detection methods, and the functional impact of protein S-palmitoylation. Additionally, this review focuses on current strategies for expanding the chemical toolbox to develop novel and effective inhibitors that can reduce the level of S-palmitoylation of the target protein.

Mechanistic insights into the allosteric inactivation mechanism of ZAP-70 induced by the hot spot W165C mutation

Zeta chain-associated protein kinase 70 (ZAP-70) is a non-receptor tyrosine kinase that interacts with the activated T-cell receptor to transduce downstream signals, and thus plays an important role in the adaptive immune system. The biphosphorylated immunotyrosine-based activation motifs (ITAM-Y2P) binds to the N-SH2 and C-SH2 domains of ZAP-70 to promote the activation of ZAP-70. The present study explores molecular mechanisms of allosteric inactivation of ZAP-70 induced by the hot spot W165C mutation through atomically detailed molecular dynamics simulation approaches. We report microsecond-length simulations of two states of the tandem SH2 domains of ZAP-70 in complex with the ITAM-Y2P motif, including the wild-type and W165C mutant. Extensive analysis of local flexibility and dynamical correlated motions show that W165C mutation changes coupled motions of protein domains and community networks. The binding affinities of the ITAM-Y2P motif to the wild-type and W165C mutant of ZAP-70 are predicted using binding free energy calculations. The results suggest that the driving force to decrease the binding affinity in the W165C mutant derives from the difference in the protein-protein electrostatic interactions. Moreover, the per-residue free energy decomposition unravels that the contributions from residues in the phosphorylated Tyr315 (pY315) binding site, in particular pY315 of ITAM-Y2P, and Arg43, Tyr240 of ZAP-70, are the key determinants for the loss of binding affinity. This study may insights into our understanding of the pathological mechanism of ZAP-70.Communicated by Ramaswamy H. Sarma.

The structure of the Caenorhabditis elegans TMC-2 complex suggests roles of lipid-mediated subunit contacts in mechanosensory transduction

Mechanotransduction is the process by which a mechanical force, such as touch, is converted into an electrical signal. Transmembrane channel-like (TMC) proteins are an evolutionarily conserved family of membrane proteins whose function has been linked to a variety of mechanosensory processes, including hearing and balance sensation in vertebrates and locomotion in Drosophila. TMC1 and TMC2 are components of ion channel complexes, but the molecular features that tune these complexes to diverse mechanical stimuli are unknown. Caenorhabditis elegans express two TMC homologs, TMC-1 and TMC-2, both of which are the likely pore-forming subunits of mechanosensitive ion channels but differ in their expression pattern and functional role in the worm. Here, we present the single-particle cryo-electron microscopy structure of the native TMC-2 complex isolated from C. elegans. The complex is composed of two copies of the pore-forming TMC-2 subunit, the calcium and integrin binding protein CALM-1 and the transmembrane inner ear protein TMIE. Comparison of the TMC-2 complex to the recently published cryo-EM structure of the C. elegans TMC-1 complex highlights conserved protein-lipid interactions, as well as a π-helical structural motif in the pore-forming helices, that together suggest a mechanism for TMC-mediated mechanosensory transduction.

Decoupling the dynamic mechanism revealed by FGFR2 mutation-induced population shift

The fibroblast growth factor receptor 2 (FGFR2) is a key component in cellular signaling networks, and its dysfunctional activation has been implicated in various diseases including cancer and developmental disorders. Mutations at the activation loop (A-loop) have been suggested to trigger an increased basal kinase activity. However, the molecular mechanism underlying this highly dynamic process has not been fully understood due to the limitation of static structural information. Here, we conducted multiple, large-scale Gaussian accelerated molecular dynamics simulations of five (K659E, K659N, K659M, K659Q, and K659T) FGFR2 mutants at the A-loop, and comprehensively analyzed the dynamic molecular basis of FGFR2 activation. The results quantified the population shift of each system, revealing that all mutants had a higher proportion of active-like states. Using Markov state models, we extracted the representative structure of different conformational states and identified key residues related to the increased kinase activity. Furthermore, community network analysis showed enhanced information connections in the mutants, highlighting the long-range allosteric communication between the A-loop and the hinge region. Our findings may provide insights into the dynamic mechanism for FGFR2 dysfunctional activation and allosteric drug discovery.Communicated by Ramaswamy H. Sarma.

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.

Decoding the deactivation mechanism of R192W mutation of ZAP-70 using molecular dynamics simulations and binding free energy calculations

CONTEXT

ZAP-70 (zeta-chain-associated protein of 70 kDa), serving as a critical regulator for T cell antigen receptor signaling, represents an attractive therapeutic target for autoimmunity disease. How the mechanistical mechanism of ZAP-70 to a human autoimmune syndrome-associated R192W mutation remains unclear. The results indicated that the R192W mutation of ZAP-70 clearly affected the conformational flexibility of the N-terminal ITAM-Y2P. Structural analysis unveiled that the R192W mutation of ZAP-70 caused the exposure of the N-terminal ITAM-Y2P to the solvent. MM-GBSA binding free energy calculations exhibited that the R192W mutation decreased the binding affinity of ITAM-Y2P to the ZAP-70 mutant. Residue-based free energy decomposition further revealed that the protein-peptide interaction networks involving electrostatic interactions provide significant contributions for complex formation. The energy unfavorable residues include Arg43, Arg192, Tyr240, and Lys244 from ZAP-70 and Asn301, Leu303, pY304, and pY315 from ITAM-Y2P in the R192W mutant. Our obtained results may help the understanding of the deactivation mechanism of ZAP-70 induced by the R192W mutation.

METHODS

In the work, multiple replica molecular dynamics simulations and molecular mechanics-generalized Born surface area (MM-GBSA) method were performed to reveal the doubly phosphorylated ITAMs (ITAM-Y2P)-mediated deactivation mechanism of ZAP-70 induced by the R192W mutation.

Decoding molecular mechanism of species-selective targeting of fungal versus human HSP90 using multiple replica molecular dynamics simulations and binding free energy calculations

As a highly evolutionarily conserved molecular chaperone, heat shock protein (HSP90), plays an important role in virulence traits, representing a therapeutic target for the treatment of fungal infections. The close evolutionary relationship between fungi and their human hosts poses a key challenge for the development of selective antifungal agents. In this work, molecular docking, multiple replica microsecond-based molecular dynamics (MD) simulations, and binding free energy calculations were performed to decode molecular mechanism of species-selective targeting of fungal versus human HSP90 triggered by the compound A11. MD simulations reveal that binding of compound A11 to human HSP90 nucleotide-binding domain (NBD) leads to obvious conformational changes relative to fungal HSP90 NBD. Binding free energy calculations show that the binding of compound A11 to fungal HSP90 NBD is stronger than that to human HSP90 NBD. Per residue-based free energy decomposition analysis was used to evaluate the inhibitor - residue interaction profile. The results efficiently identify the hot spot residues that play vital roles in favorable binding of compound A11 to fungal HSP90 NBD. This study is expected to provide a useful guidance for the development of selective inhibitors toward fungal HSP90.Communicated by Ramaswamy H. Sarma.

Mechanistic Elucidation of Activation/Deactivation Signal Transduction within Neurotensin Receptor 1 Triggered by 'Driver Chemical Groups' of Modulators: A Comparative Molecular Dynamics Simulation

Small-molecule modulators of neurotensin receptor 1 (NTSR1), a class A G-protein-coupled receptor (GPCR), has emerged as promising therapeutic agent for psychiatric disorders and cancer. Interestingly, a chemical group substitution in NTSR1 modulators can launch different types of downstream regulation, highlighting the significance of deciphering the internal fine-tuning mechanism. Here, we conducted a synergistic application of a Gaussian accelerated molecular dynamics simulation, a conventional molecular dynamics simulation, and Markov state models (MSM) to investigate the underlying mechanism of 'driver chemical groups' of modulators triggering inverse signaling. The results indicated that the flexibility of the leucine moiety in NTSR1 agonists contributes to the inward displacement of TM7 through a loosely coupled allosteric pathway, while the rigidity of the adamantane moiety in NTSR1 antagonists leads to unfavorable downward transduction of agonistic signaling. Furthermore, we found that R322^6.54, Y319^6.51, F353^7.42, R148^3.32, S356^7.45, and S357^7.46 may play a key role in inducing the activation of NTSR1. Together, our findings not only highlight the ingenious signal transduction within class A GPCRs but also lay a foundation for the development of targeted drugs harboring different regulatory functions of NTSR1.

Decoding the Conformational Selective Mechanism of FGFR Isoforms: A Comparative Molecular Dynamics Simulation

Fibroblast growth factor receptors (FGFRs) play critical roles in the regulation of cell growth, differentiation, and proliferation. Specifically, FGFR2 gene amplification has been implicated in gastric and breast cancer. Pan-FGFR inhibitors often cause large toxic side effects, and the highly conserved ATP-binding pocket in the FGFR1/2/3 isoforms poses an immense challenge in designing selective FGFR2 inhibitors. Recently, an indazole-based inhibitor has been discovered that can selectively target FGFR2. However, the detailed mechanism involved in selective inhibition remains to be clarified. To this end, we performed extensive molecular dynamics simulations of the apo and inhibitor-bound systems along with multiple analyses, including Markov state models, principal component analysis, a cross-correlation matrix, binding free energy calculation, and community network analysis. Our results indicated that inhibitor binding induced the phosphate-binding loop (P-loop) of FGFR2 to switch from the open to the closed conformation. This effect enhanced extensive hydrophobic FGFR2-inhibitor contacts, contributing to inhibitor selectivity. Moreover, the key conformational intermediate states, dynamics, and driving forces of this transformation were uncovered. Overall, these findings not only provided a structural basis for understanding the closed P-loop conformation for therapeutic potential but also shed light on the design of selective inhibitors for treating specific types of cancer.