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

Computational Insights into Acrylamide Fragment Inhibition of SARS-CoV-2 Main Protease

The pathogen of COVID-19, SARS-CoV-2, has caused a severe global health crisis. So far, while COVID-19 has been suppressed, the continuous evolution of SARS-CoV-2 variants has reduced the effectiveness of vaccines such as mRNA-1273 and drugs such as Remdesivir. To uphold the effectiveness of vaccines and drugs prior to potential coronavirus outbreaks, it is necessary to explore the underlying mechanisms between biomolecules and nanodrugs. The experimental study reported that acrylamide fragments covalently attached to Cys145, the main protease enzyme (Mpro) of SARS-CoV-2, and occupied the substrate binding pocket, thereby disrupting protease dimerization. However, the potential mechanism linking them is unclear. The purpose of this work is to complement and validate experimental results, as well as to facilitate the study of novel antiviral drugs. Based on our experimental studies, we identified two acrylamide fragments and constructed corresponding protein-ligand complex models. Subsequently, we performed molecular dynamics (MD) simulations to unveil the crucial interaction mechanisms between these nanodrugs and SARS-CoV-2 Mpro. This approach allowed the capture of various binding conformations of the fragments on both monomeric and dimeric Mpro, revealing significant conformational dissociation between the catalytic and helix domains, which indicates the presence of allosteric targets. Notably, Compound 5 destabilizes Mpro dimerization and acts as an effective inhibitor by specifically targeting the active site, resulting in enhanced inhibitory effects. Consequently, these fragments can modulate Mpro's conformational equilibrium among extended monomeric, compact, and dimeric forms, shedding light on the potential of these small molecules as novel inhibitors against coronaviruses. Overall, this research contributes to a broader understanding of drug development and fragment-based approaches in antiviral covalent therapeutics.

Stereo-selectivity of enantiomeric inhibitors to ubiquitin-specific protease 7 (USP7) dissected by molecular docking, molecular dynamics simulations, and binding free energy calculations

The ubiquitin-specific protease 7 (USP7), as a member of deubiquitination enzymes, represents an attractive therapeutic target for various cancers, including prostate cancer and liver cancer. The change of the inhibitor stereocenter from the S to R stereochemistry (S-ALM → R-ALM34) markedly improved USP7 inhibitory activity. However, the molecular mechanism for the stereo-selectivity of enantiomeric inhibitors to USP7 is still unclear. In this work, molecular docking, molecular dynamics (MD) simulations, molecular mechanics/Generalized-Born surface area (MM/GBSA) calculations, and free energy landscapes were performed to address this mystery. MD simulations revealed that S-ALM34 showed a high degree of conformational flexibility compared to the R-ALM34 counterpart, and S-ALM34 binding led to the enhanced intradomain motions of USP7, especially the BL1 and BL2 loops and the two helices α4 and α5. MM/GBSA calculations showed that the binding strength of R-ALM34 to USP7 was stronger than that of S-ALM34 by - 4.99 kcal/mol, a similar trend observed by experimental data. MM/GBSA free energy decomposition was further performed to differentiate the ligand-residue spectrum. These analyses not only identified the hotspot residues interacting with R-ALM34, but also revealed that the hydrophobic interactions from F409, K420, H456, and Y514 play the major determinants in the binding of R-ALM34 to USP7. This result is anticipated to shed light on energetic basis and conformational dynamics information to aid in the design of more potent and selective inhibitors targeting USP7.

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.

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.

Delineating the stepwise millisecond allosteric activation mechanism of the class C GPCR dimer mGlu5

Two-thirds of signaling hormones and one-third of approved drugs exert their effects by binding and modulating the G protein-coupled receptors (GPCRs) activation. While the activation mechanism for monomeric GPCRs has been well-established, little is known about GPCRs in dimeric form. Here, by combining transition pathway generation, extensive atomistic simulation-based Markov state models, and experimental signaling assays, we reveal an asymmetric, stepwise millisecond allosteric activation mechanism for the metabotropic glutamate receptor subtype 5 receptor (mGlu5), an obligate dimeric class C GPCR. The dynamic picture is presented that agonist binding induces dimeric ectodomains compaction, amplified by the precise association of the cysteine-rich domains, ultimately loosely bringing the intracellular 7-transmembrane (7TM) domains into proximity and establishing an asymmetric TM6-TM6 interface. The active inter-domain interface enhances their intra-domain flexibility, triggering the activation of micro-switches crucial for downstream signal transduction. Furthermore, we show that the positive allosteric modulator stabilizes both the active inter-domain 7TM interface and an open, extended intra-domain ICL2 conformation. This stabilization leads to the formation of a pseudo-cavity composed of the ICL2, ICL3, TM3, and C-terminus, which facilitates G protein coordination. Our strategy may be generalizable for characterizing millisecond events in other allosteric systems.

The switch states of the GDP-bound HRAS affected by point mutations: a study from Gaussian accelerated molecular dynamics simulations and free energy landscapes

Point mutations play a vital role in the conformational transformation of HRAS. In this work, Gaussian accelerated molecular dynamics (GaMD) simulations followed by constructions of free energy landscapes (FELs) were adopted to explore the effect of mutations D33K, A59T and L120A on conformation states of the GDP-bound HRAS. The results from the post-processing analyses on GaMD trajectories suggest that mutations alter the flexibility and motion modes of the switch domains from HRAS. The analyses from FELs show that mutations induce more disordered states of the switch domains and affect interactions of GDP with HRAS, implying that mutations yield a vital effect on the binding of HRAS to effectors. The GDP-residue interaction network revealed by our current work indicates that salt bridges and hydrogen bonding interactions (HBIs) play key roles in the binding of GDP to HRAS. Furthermore, instability in the interactions of magnesium ions and GDP with the switch SI leads to the extreme disorder of the switch domains. This study is expected to provide the energetic basis and molecular mechanism for further understanding the function of HRAS.Communicated by Ramaswamy H. Sarma.

Structural and energetic insights into the selective inhibition of PKMYT1 against WEE1

Protein kinase, membrane-associated tyrosine/threonine 1 (PKMYT1), a member of the WEE family and responsible for the regulation of CDK1 phosphorylation, has been considered a promising therapeutic target for cancer therapy. However, the highly structural conservation of the ATP-binding sites of the WEE family poses a challenge to the design of selective inhibitors for PKMYT1. Here, molecular docking, multiple microsecond-length molecular dynamics (MD) simulations and end-point free energy calculations were performed to uncover the molecular mechanism of the binding selectivity of RP-6306 toward PKMYT1 over its highly homologous kinase WEE1. The binding specificity of RP-6306 reported in previous experimental bioassays was clarified by MD simulations and binding free energy calculations. Further, the binding free energy prediction indicated that the binding selectivity of RP-6306 largely derived from the difference in the protein-ligand electrostatic interactions. The per-residue free energy decomposition suggested that the non-conserved gatekeeper residue in the hinge domain of PKMYT1/WEE1, Thr187/Asn376, is the critical factor responsible for the binding selectivity of RP-6306 toward PKMYT1. In addition, a water-mediated hydrogen bond was formed between RP-6306 and Gly191 at the hinge domain in the PKMYT1/RP-6306 complex, which was absent in the WEE1/RP-6306 complex. This study is expected to offer useful information for the design of more potent and selective PKMYT1 inhibitors.Communicated by Ramaswamy H. Sarma.

The effect of phosphorylation on the conformational dynamics and allostery of the association of death-associated protein kinase with calmodulin

Protein phosphorylation plays an important role in the signal transduction and is capable of regulation of cell activity. The death-associated protein kinase 1 (DAPK1), as a Ser/Thr kinase, interacts with calmodulin (CaM) to regulate apoptotic and autophagic signaling. Autophosphorylation of DAPK1 at Ser308 located at the autoregulatory domain (ARD) blocks CaM binding and inhibits kinase catalytic activity. However, the mechanism underlying the influence of Ser308 phosphorylation (pS308) on the DAPK1 activity remains unclear. Here, we performed multiple, microsecond length molecular dynamics (MD) simulations, the molecular mechanics generalized Born/surface area (MM-GBSA) binding free energy calculations, principal component analysis, and dynamic cross-correlation analysis to unravel the conformational dynamics and allostery of the DAPK1 - CaM interaction triggered by the pS308 at the ARD. MD simulations showed that pS308 affected the conformational stability of the DAPK1 - CaM complex. Further energetic and structural exploration revealed that pS308 weakened the association of the phosphorylated DAPK1 to CaM, which lowered the susceptibility of DAPK1 to be activated by CaM. This result can provide mechanistic insights into the molecular underpinning through which the DAPK1 kinase activity is modulated by the auto-phosphorylation.Communicated by Ramaswamy H. Sarma.

Unraveling the Interplay of Extracellular Domain Conformational Changes and Parathyroid Hormone Type 1 Receptor Activation in Class B1 G Protein-Coupled Receptors: Integrating Enhanced Sampling Molecular Dynamics Simulations and Markov State Models

Parathyroid hormone (PTH) type 1 receptor (PTH1R), as a typical class B1 G protein-coupled receptor (GPCR), is responsible for regulating bone turnover and maintaining calcium homeostasis, and its dysregulation has been implicated in the development of several diseases. The extracellular domain (ECD) of PTH1R is crucial for the recognition and binding of ligands, and the receptor may exhibit an autoinhibited state with the closure of the ECD in the absence of ligands. However, the correlation between ECD conformations and PTH1R activation remains unclear. Thus, this study combines enhanced sampling molecular dynamics (MD) simulations and Markov state models (MSMs) to reveal the possible relevance between the ECD conformations and the activation of PTH1R. First, 22 intermediate structures are generated from the autoinhibited state to the active state and conducted for 10 independent 200 ns simulations each. Then, the MSM is constructed based on the cumulative 44 μs simulations with six identified microstates. Finally, the potential interplay between ECD conformational changes and PTH1R activation as well as cryptic allosteric pockets in the intermediate states during receptor activation is revealed. Overall, our findings reveal that the activation of PTH1R has a specific correlation with ECD conformational changes and provide essential insights for GPCR biology and developing novel allosteric modulators targeting cryptic sites.

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.

Conformational States of the GDP- and GTP-Bound HRAS Affected by A59E and K117R: An Exploration from Gaussian Accelerated Molecular Dynamics

The HRAS protein is considered a critical target for drug development in cancers. It is vital for effective drug development to understand the effects of mutations on the binding of GTP and GDP to HRAS. We conducted Gaussian accelerated molecular dynamics (GaMD) simulations and free energy landscape (FEL) calculations to investigate the impacts of two mutations (A59E and K117R) on GTP and GDP binding and the conformational states of the switch domain. Our findings demonstrate that these mutations not only modify the flexibility of the switch domains, but also affect the correlated motions of these domains. Furthermore, the mutations significantly disrupt the dynamic behavior of the switch domains, leading to a conformational change in HRAS. Additionally, these mutations significantly impact the switch domain's interactions, including their hydrogen bonding with ligands and electrostatic interactions with magnesium ions. Since the switch domains are crucial for the binding of HRAS to effectors, any alterations in their interactions or conformational states will undoubtedly disrupt the activity of HRAS. This research provides valuable information for the design of drugs targeting HRAS.

Roles of Accelerated Molecular Dynamics Simulations in Predictions of Binding Kinetic Parameters

Rational predictions on binding kinetics parameters of drugs to targets play significant roles in future drug designs. Full conformational samplings of targets are requisite for accurate predictions of binding kinetic parameters. In this review, we mainly focus on the applications of enhanced sampling technologies in calculations of binding kinetics parameters and residence time of drugs. The methods involved in molecular dynamics simulations are applied to not only probe conformational changes of targets but also reveal calculations of residence time that is significant for drug efficiency. For this review, special attention are paid to accelerated molecular dynamics (aMD) and Gaussian aMD (GaMD) simulations that have been adopted to predict the association or disassociation rate constant. We also expect that this review can provide useful information for future drug design.

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.

Progressive computational approaches to facilitate decryption of allosteric mechanism and drug discovery

Allostery is a ubiquitous biological phenomenon where perturbation at topologically distal areas of a protein serves as a trigger to fine-tune the orthosteric site and thus regulate protein function. The investigation of allosteric regulation greatly enhances our understanding of human diseases and broadens avenue for drug discovery. For decades, owing to the difficulty in allostery characterization through serendipitous experimental screening, researchers have developed several innovative computational approaches, which proves to accelerate the elucidation of allostery. Herein, we review the state-of-the-art advance of computational methodologies for allostery study, with particular emphasis on promising trends emerging over the past two years. We expect this review will outline the latest landscape of allostery study and inspire researchers to further facilitate this field.

Exploring conformation changes of Janus kinase 2 pseudokinase mediated by mutations through Gaussian accelerated molecular dynamics and principal component analysis

The pseudokinase domain (JH2) of the protein tyrosine kinase (Janus kinase 2, JAK2) regulates the activity of a tyrosine kinase domain (JH1) in JAK2, which is further affected by mutations in the JH2. In this work, Gaussian accelerated molecular dynamics (GaMD) simulations followed by construction of free energy landscapes (FELs) and principal component analysis (PCA) were performed to study effect of two mutations V617F and V617F/E596A on the conformations of the ATP-bound JH2. The dynamic analyses reveal that mutations affect the structural flexibility and correlated motions of the JH2, meanwhile also change the dynamics behavior of the P-loop and αC-helix of the JH2. The information from FELs unveils that mutations induce less energy states than the free JH2 and the WT one. The analyses of interaction networks uncover that mutations affect the salt bridge interactions of ATP with K581, K677 and R715 and alter hydrogen bonding interactions (HBIs) of ATP with the JH2. The changes in conformations of the JH2 and ATP-JH2 interaction networks caused by mutations in turn generate effect on the activity regulations of the JH2 on the JH1. This work is expected to provide significant theoretical helps for deeply understanding the function of the JH2 and drug design toward JAK2.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.

Deconstructing allostery by computational assessment of the binding determinants of allosteric PTP1B modulators

Fragment-based drug discovery is an established methodology for finding hit molecules that can be elaborated into lead compounds. However it is currently challenging to predict whether fragment hits that do not bind to an orthosteric site could be elaborated into allosteric modulators, as in these cases binding does not necessarily translate into a functional effect. We propose a workflow using Markov State Models (MSMs) with steered molecular dynamics (sMD) to assess the allosteric potential of known binders. sMD simulations are employed to sample protein conformational space inaccessible to routine equilibrium MD timescales. Protein conformations sampled by sMD provide starting points for seeded MD simulations, which are combined into MSMs. The methodology is demonstrated on a dataset of protein tyrosine phosphatase 1B ligands. Experimentally confirmed allosteric inhibitors are correctly classified as inhibitors, whereas the deconstructed analogues show reduced inhibitory activity. Analysis of the MSMs provide insights into preferred protein-ligand arrangements that correlate with functional outcomes. The present methodology may find applications for progressing fragments towards lead molecules in FBDD campaigns.

Negative allosteric modulation of the glucagon receptor by RAMP2

Receptor activity-modifying proteins (RAMPs) modulate the activity of many Family B GPCRs. We show that RAMP2 directly interacts with the glucagon receptor (GCGR), a Family B GPCR responsible for blood sugar homeostasis, and broadly inhibits receptor-induced downstream signaling. HDX-MS experiments demonstrate that RAMP2 enhances local flexibility in select locations in and near the receptor extracellular domain (ECD) and in the 6th transmembrane helix, whereas smFRET experiments show that this ECD disorder results in the inhibition of active and intermediate states of the intracellular surface. We determined the cryo-EM structure of the GCGR-Gs complex at 2.9 Å resolution in the presence of RAMP2. RAMP2 apparently does not interact with GCGR in an ordered manner; however, the receptor ECD is indeed largely disordered along with rearrangements of several intracellular hallmarks of activation. Our studies suggest that RAMP2 acts as a negative allosteric modulator of GCGR by enhancing conformational sampling of the ECD.

Impacts of Mutations in the P-Loop on Conformational Alterations of KRAS Investigated with Gaussian Accelerated Molecular Dynamics Simulations

G12 mutations heavily affect conformational transformation and activity of KRAS. In this study, Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the GDP-bound wild-type (WT), G12A, G12D, and G12R KRAS to probe mutation-mediated impacts on conformational alterations of KRAS. The results indicate that three G12 mutations obviously affect the structural flexibility and internal dynamics of the switch domains. The analyses of the free energy landscapes (FELs) suggest that three G12 mutations induce more conformational states of KRAS and lead to more disordered switch domains. The principal component analysis shows that three G12 mutations change concerted motions and dynamics behavior of the switch domains. The switch domains mostly overlap with the binding region of KRAS to its effectors. Thus, the high disorder states and concerted motion changes of the switch domains induced by G12 mutations affect the activity of KRAS. The analysis of interaction network of GDP with KRAS signifies that the instability in the interactions of GDP and magnesium ion with the switch domain SW1 drives the high disordered state of the switch domains. This work is expected to provide theoretical aids for understanding the function of KRAS.

The allosteric effect of the upper half of SENP1 contributes to its substrate selectivity for SUMO1 over SUMO2

SUMOylation regulates various cellular process and SENP1 (SUMO-specific protease 1) serves as a SUMO (small ubiquitin-related modifier) specific protease that participates in the SUMO cycle. Given its extensive influences on metabolic activities, SENP1 has gained more and more attentions in clinical treatments. However, there remains a question on why does the SENP1 prefer to process SUMO1 rather than SUMO2. Here, we performed molecular dynamics simulations of SENP1-SUMO1, SENP1-SUMO2, and apo SENP1 systems and observed distinct conformational dynamics in the upper half of the clamp and the three loops in the catalytic center of the SENP1. Principal component analysis revealed that the most prominent canonical variable represented the spatial distribution of the upper half of the clamp, while the openness of the cleft was closely related to the catalytic ability of SENP1. Further analysis of the SENP1-SUMO interactions revealed that the extensive and strong interactions between the SENP1 and SUMO1 were both in the interface of the upper half region and the catalytic center. Dynamic cross-correlation matrix analysis demonstrated that the inter-residue correlations in the SUMO1 system was much stronger, especially in the two essential regions belonging to the upper and lower half of cleft. Based on these observations, we proposed an allosteric propagation model and further testified it using the community analysis. These results revealed the propagation pathway of allosteric communication that contributed to the substrate discrimination of SENP1 upon SUMO1 and SUMO2.Communicated by Ramaswamy H. Sarma.

Binding modes of GDP, GTP and GNP to NRAS deciphered by using Gaussian accelerated molecular dynamics simulations

Probing binding modes of GDP, GTP and GNP to NRAS are of significance for understanding the regulation mechanism on the activity of RAS proteins. Four separate Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the apo, GDP-, GTP- and GNP-bound NRAS. Dynamics analyses suggest that binding of three ligands highly affects conformational states of the switch domains from NRAS, which disturbs binding of NRAS to its effectors. The analyses of free energy landscapes (FELs) indicate that binding of GDP, GTP and GNP induces more energetic states of NRAS compared to the apo NRAS but the presence of GNP makes the switch domains more ordered than binding of GDP and GNP. The information of interaction networks of ligands with NRAS reveals that the π-π interaction of residue F28 and the salt bridge interactions of K16 and D119 with ligands stabilize binding of GDP, GTP and GNP to NRAS. Meanwhile magnesium ion plays a bridge role in interactions of ligands with NRAS, which is favourable for associations of GDP, GTP and GNP with NRAS. This work is expected to provide useful information for deeply understanding the function and activity of NRAS.

Insights into pralsetinib resistance to the non-gatekeeper RET kinase G810C mutation through molecular dynamics simulations

OBJECTIVE

RET (rearranged during transfection) kinase, as a transmembrane receptor tyrosine kinase, is a therapeutic target for several human cancer such as non-small cell lung cancer (NSCLC) and thyroid cancer. Pralsetinib is a recently approved drug for the treatment of RET-driven NSCLC and thyroid cancers. A single point mutation G810C at the C-lobe of the RET kinase causes pralsetinib resistance to this non-gatekeeper variant. However, the detailed mechanism remains poorly understood.

METHODS

Here, multiple microsecond molecular dynamics (MD) simulations, molecular mechanics/generalized born surface area (MM/GBSA) binding free energy calculations, and community network analysis were performed to reveal the mechanism of pralsetinib resistance to the RET G810C mutant.

RESULTS

The simulations showed that the G810C mutation had a minor effect on the overall conformational dynamics of the RET kinase domain. Energetic analysis suggested that the G810C mutation reduced the binding affinity of pralsetinib to the mutant. Per-residue energy contribution and structural analyses revealed that the hydrogen bonding interactions between pralsetinib and the hinge residues Glu805 and Ala807 were disrupted in the G810C mutant, which were responsible for the decreased binding affinity of pralsetinib to the mutant.

CONCLUSIONS

The obtained results may provide understanding of the mechanism of pralsetinib resistance to the non-gatekeeper RET G810C mutant.

New insights into the structure and function of class B1 GPCRs

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors. Class B1 GPCRs constitute a subfamily of 15 receptors that characteristically contain large extracellular domains (ECDs) and respond to long polypeptide hormones. Class B1 GPCRs are critical regulators of homeostasis, and as such, many are important drug targets. While most transmembrane proteins, including GPCRs, are recalcitrant to crystallization, recent advances in electron cryo-microscopy (cryo-EM) have facilitated a rapid expansion of the structural understanding of membrane proteins. As a testament to this success, structures for all the class B1 receptors bound to G proteins have been determined by cryo-EM in the past five years. Further advances in cryo-EM have uncovered dynamics of these receptors, ligands, and signalling partners. Here, we examine the recent structural underpinnings of the class B1 GPCRs with an emphasis on structure-function relationships.

Mechanistic insights into the role of calcium in the allosteric regulation of the calmodulin-regulated death-associated protein kinase

Calcium (Ca²⁺) signaling plays an important role in the regulation of many cellular functions. Ca²⁺-binding protein calmodulin (CaM) serves as a primary effector of calcium function. Ca²⁺/CaM binds to the death-associated protein kinase 1 (DAPK1) to regulate intracellular signaling pathways. However, the mechanism underlying the influence of Ca²⁺ on the conformational dynamics of the DAPK1-CaM interactions is still unclear. Here, we performed large-scale molecular dynamics (MD) simulations of the DAPK1-CaM complex in the Ca²⁺-bound and-unbound states to reveal the importance of Ca²⁺. MD simulations revealed that removal of Ca²⁺ increased the anti-correlated inter-domain motions between DAPK1 and CaM, which weakened the DAPK1-CaM interactions. Binding free energy calculations validated the decreased DAPK1-CaM interactions in the Ca²⁺-unbound state. Structural analysis further revealed that Ca²⁺ removal caused the significant conformational changes at the DAPK1-CaM interface, especially the helices α1, α2, α4, α6, and α7 from the CaM and the basic loop and the phosphate-binding loop from the DAPK1. These results may be useful to understand the biological role of Ca²⁺ in physiological processes.

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.

Deciphering Conformational Changes of the GDP-Bound NRAS Induced by Mutations G13D, Q61R, and C118S through Gaussian Accelerated Molecular Dynamic Simulations

The conformational changes in switch domains significantly affect the activity of NRAS. Gaussian-accelerated molecular dynamics (GaMD) simulations of three separate replicas were performed to decipher the effects of G13D, Q16R, and C118S on the conformational transformation of the GDP-bound NRAS. The analyses of root-mean-square fluctuations and dynamics cross-correlation maps indicated that the structural flexibility and motion modes of the switch domains involved in the binding of NRAS to effectors are highly altered by the G13D, Q61R, and C118Smutations. The free energy landscapes (FELs) suggested that mutations induce more energetic states in NRAS than the GDP-bound WT NRAS and lead to high disorder in the switch domains. The FELs also indicated that the different numbers of sodium ions entering the GDP binding regions compensate for the changes in electrostatic environments caused by mutations, especially for G13D. The GDP-residue interactions revealed that the disorder in the switch domains was attributable to the unstable hydrogen bonds between GDP and two residues, V29 and D30. This work is expected to provide information on the energetic basis and dynamics of conformational changes in switch domains that can aid in deeply understanding the target roles of NRAS in anticancer treatment.

Markov State Models and Molecular Dynamics Simulations Reveal the Conformational Transition of the Intrinsically Disordered Hypervariable Region of K-Ras4B to the Ordered Conformation

K-Ras4B, the most frequently mutated Ras isoform in human tumors, plays a vital part in cell growth, differentiation, and survival. Its tail, the C-terminal hypervariable region (HVR), is involved in anchoring K-Ras4B at the cellular plasma membrane and in isoform-specific protein-protein interactions and signaling. In the inactive guanosine diphosphate-bound state, the intrinsically disordered HVR interacts with the catalytic domain at the effector-binding region, rendering K-Ras4B in its autoinhibited state. Activation releases the HVR from the catalytic domain, with its ensemble favoring an ordered α-helical structure. The large-scale conformational transition of the HVR from the intrinsically disordered to the ordered conformation remains poorly understood. Here, we deploy a computational scheme that integrates a transition path-generation algorithm, extensive molecular dynamics simulation, and Markov state model analysis to investigate the conformational landscape of the HVR transition pathway. Our findings reveal a stepwise pathway for the HVR transition and uncover several key conformational substates along the transition pathway. Importantly, key interactions between the HVR and the catalytic domain are unraveled, highlighting the pathogenesis of K-Ras4B mild mutations in several congenital developmental anomaly syndromes. Together, these findings provide a deeper understanding of the HVR transition mechanism and the regulation of K-Ras4B activity at an atomic level.

Understanding gilteritinib resistance to FLT3-F691L mutation through an integrated computational strategy

FMS-like tyrosine kinase 3 (FLT3) serves as an important drug target for acute myeloid leukemia (AML), and gene mutations of FLT3 have been closely associated with AML patients with an incidence rate of ~ 30%. However, the mechanism of the clinically relevant F691L gatekeeper mutation conferred resistance to the drug gilteritinib remained poorly understood. In this study, multiple microsecond molecular dynamics (MD) simulations, end-point free energy calculations, and dynamic correlated and network analyses were performed to investigate the molecular basis of gilteritinib resistance to the FLT3-F691L mutation. The simulations revealed that the resistant mutation largely induced the conformational changes of the activation loop (A-loop), the phosphate-binding loop, and the helix αC of the FLT3 protein. The binding abilities of the gilteritinib to the wild-type and the F691L mutant were different through the binding free energy prediction. The simulation results further indicated that the driving force to determine the binding affinity of gilteritinib was derived from the differences in the energy terms of electrostatic and van der Waals interactions. Moreover, the per-residue free energy decomposition suggested that the four residues (Phe803, Gly831, Leu832, and Ala833) located at the A-loop of FLT3 had a significant impact on the binding affinity of gilteritinib to the F691L mutant. This study may provide useful information for the design of novel FLT3 inhibitors specially targeting the F691L gatekeeper mutant.

Mechanistic Insights into the Mechanism of Inhibitor Selectivity toward the Dark Kinase STK17B against Its High Homology STK17A

As a member of the death-associated protein kinase (DAPK) family, STK17B plays an important role in the regulation of cellular apoptosis and has been considered as a promising drug target for hepatocellular carcinoma. However, the highly conserved ATP-binding site of protein kinases represents a challenge to design selective inhibitors for a specific DAPK isoform. In this study, molecular docking, multiple large-scale molecular dynamics (MD) simulations, and binding free energy calculations were performed to decipher the molecular mechanism of the binding selectivity of PKIS43 toward STK17B against its high homology STK17A. MD simulations revealed that STK17A underwent a significant conformational arrangement of the activation loop compared to STK17B. The binding free energy predictions suggested that the driving force to control the binding selectivity of PKIS43 was derived from the difference in the protein–ligand electrostatic interactions. Furthermore, the per-residue free energy decomposition unveiled that the energy contribution from Arg41 at the phosphate-binding loop of STK17B was the determinant factor responsible for the binding specificity of PKIS43. This study may provide useful information for the rational design of novel and potent selective inhibitors toward STK17B.

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.

Autopromotion of K-Ras4B Feedback Activation Through an SOS-Mediated Long-Range Allosteric Effect

The Ras-specific guanine nucleotide exchange factors Son of Sevenless (SOS) regulates Ras activation by converting inactive GDP-bound to active GTP-bound states. The catalytic activity of Ras is further allosterically regulated by GTP-Ras bound to a distal site through a positive feedback loop. To address the mechanism underlying the long-range allosteric activation of the catalytic K-Ras4B by an additional allosteric GTP-Ras through SOS, we employed molecular dynamics simulation of the K-Ras4BG13D•SOScat complex with and without an allosteric GTP-bound K-Ras4BG13D. We found that the binding of an allosteric GTP-K-Ras4BG13D enhanced the affinity between the catalytic K-Ras4BG13D and SOScat, forming a more stable conformational state. The peeling away of the switch I from the nucleotide binding site facilitated the dissociation of GDP, thereby contributing to the increased nucleotide exchange rate. The community networks further showed stronger edge connection upon allosteric GTP-K-Ras4BG13D binding, which represented an increased interaction between catalytic K-Ras4BG13D and SOScat. Moreover, GTP-K-Ras4BG13D binding transmitted allosteric signaling pathways though the Cdc25 domain of SOS that enhanced the allosteric regulatory from the K-Ras4BG13D allosteric site to the catalytic site. This study may provide an in-depth mechanism for abnormal activation and allosteric regulation of K-Ras4BG13D.

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.