Conformational states dynamically populated by a kinase determine its function

The Conformational Transition Pathways and Hidden Intermediates in DFG-Flip Process of c-Met Kinase Revealed by Metadynamics Simulations

Protein kinases intrinsically translate their conformations between active and inactive states, which is key to their enzymatic activities. The conformational flipping of the three-residue conservative motif, Asp-Phe-Gly (DFG), is crucial for many kinases' biological functions. Obtaining a detailed demonstration of the DFG flipping process and its corresponding dynamical and thermodynamical features could broaden our understanding of kinases' conformation-activity relationship. In this study, we employed metadynamics simulation, a widely used enhanced sampling technique, to analyze the conformational transition pathways of the DFG flipping for the c-Met kinase. The corresponding free energy landscape suggested two distinct transition pathways between the "DFG-in" and "DFG-out" states of the DFG-flip from c-Met. On the basis of the orientation direction of the F1223 residue, we correspondingly named the two pathways the "DFG-up" path, featuring forming a commonly discovered "DFG-up" transition state, and the "DFG-down" path, a unique transition pathway with F1223 rotating along the opposite direction away from the hydrophobic cavity. The free energies along the two pathways were then calculated using the Path Collective Variable (PCV) metadynamics simulation. The simulation results showed that, though having similar free energy barriers, the free energy cuve for the DFG-down path suggested a two-step conformational transition mechanism, while that for the DFG-up path showed the one-step transition feature. The c-Met DFG flipping mechanism and the new intermediate state discovered in this work could provide a deeper understanding of the conformation-activity relationship for c-Met and, possibly, reveal a new conformational state as the drug target for c-Met and other similar kinases.

Global protein dynamics as communication sensors in peptide synthetase domains

Biological activity is governed by the timely redistribution of molecular interactions, and static structural snapshots often appear insufficient to provide the molecular determinants that choreograph communication. This conundrum applies to multidomain enzymatic systems called nonribosomal peptide synthetases (NRPSs), which assemble simple substrates into complex metabolites, where a dynamic domain organization challenges rational design to produce new pharmaceuticals. Using a nuclear magnetic resonance (NMR) atomic-level readout of biochemical transformations, we demonstrate that global structural fluctuations help promote substrate-dependent communication and allosteric responses, and impeding these global dynamics by a point-site mutation hampers allostery and molecular recognition. Our results establish global structural dynamics as sensors of molecular events that can remodel domain interactions, and they provide new perspectives on mechanisms of allostery, protein communication, and NRPS synthesis.

Protein Flexibility and Dissociation Pathway Differentiation Can Explain Onset of Resistance Mutations in Kinases

Understanding how mutations render a drug ineffective is a problem of immense relevance. Often the mechanism through which mutations cause drug resistance can be explained purely through thermodynamics. However, the more perplexing situation is when two proteins have the same drug binding affinities but different residence times. In this work, we demonstrate how all-atom molecular dynamics simulations using recent developments grounded in statistical mechanics can provide a detailed mechanistic rationale for such variances. We discover dissociation mechanisms for the anti-cancer drug Imatinib (Gleevec) against wild-type and the N368S mutant of Abl kinase. We show how this point mutation triggers far-reaching changes in the protein's flexibility and leads to a different, much faster, drug dissociation pathway. We believe that this work marks an efficient and scalable approach to obtain mechanistic insight into resistance mutations in biomolecular receptors that are hard to explain using a structural perspective.

Evaluation of residue variability in a conformation-specific context and during evolutionary sequence reconstruction narrows drug resistance selection in Abl1 tyrosine kinase

Diseases with readily available therapies may eventually prevail against the specific treatment by the acquisition of resistance. The constitutively active Abl1 tyrosine kinase known to cause chronic myeloid leukemia is an example, where patients may experience relapse after small inhibitor drug treatment. Mutations in the Abl1 tyrosine kinase domain (Abl1-KD) are a critical source of resistance and their emergence depends on the conformational states that have been observed experimentally: the inactive state, the active state, and the intermediate inactive state that resembles Src kinase. Understanding how resistant positions and amino acid identities are determined by selection pressure during drug treatment is necessary to improve future drug development or treatment decisions. We carry out in silico site-saturation mutagenesis over the Abl1-KD structure in a conformational context to evaluate the in situ and conformational stability energy upon mutation. Out of the 11 studied resistant positions, we determined that 7 of the resistant mutations favored the active conformation of Abl1-KD with respect to the inactive state. When, instead, the sequence optimization was modeled simultaneously at resistant positions, we recovered five known resistant mutations in the active conformation. These results suggested that the Abl1 resistance mechanism targeted substitutions that favored the active conformation. Further sequence variability, explored by ancestral reconstruction in Abl1-KD, showed that neutral genetic drift, with respect to amino acid variability, was specifically diminished in the resistant positions. Since resistant mutations are susceptible to chance with a certain probability of fixation, combining methodologies outlined here may narrow and limit the available sequence space for resistance to emerge, resulting in more robust therapeutic treatments over time.

Affinities and Kinetics Detection of Protein-Small Molecule Interactions with a Monolayer MoS2 -Based Optical Imaging Platform

Quantifying the binding kinetics and affinities of protein-small molecule interactions is critical for biomarker validation, drug discovery, and deep understanding of various biological processes at the molecular-scale. Novel approaches are demanded as most common label-free techniques are mass-sensitive, which are not suitable for the detection of small molecule interactions. Here, an optical imaging platform is developed to measure the binding kinetics of both protein-small molecules and protein-ions based on monolayer MoS₂ , an ultra-thin 2D material whose optical absorption is extremely sensitive to charge. A model is established to calibrate the optical response due to the charged analyte binding and it is applied to quantify the interactions between abl1 kinase and different small-molecule inhibitors. Such a presented method is capable of distinguishing different inhibitors binding to a wild or mutated kinase, which provides guidance for drug evaluation and drug mechanism exploration. The binding kinetics of calcium ions to calmodulin is also measured, further broadening the application field of the method. In addition, the imaging capability allows mapping the local binding kinetics of the molecular interactions with a high resolution, which reveals visible spatial variability and offers a promising tool for studying heterogeneous local interfacial interactions.

The structural basis of BCR-ABL recruitment of GRB2 in chronic myelogenous leukemia

BCR-ABL drives chronic myeloid leukemia (CML). BCR binding to GRB2 transduces signaling via the Ras/MAPK pathway. Despite considerable data confirming the binding, molecular-level understanding of exactly how the two proteins interact, and, especially, what are the determinants of the specificity of the SH2GRB2 domain-phosphorylated BCR (pBCR) recognition are still open questions. Yet, this is vastly important for understanding binding selectivity, and for predicting the phosphorylated receptors, or peptides, that are likely to bind. Here, we uncover these determinants and ascertain to what extent they relate to the affinity of the interaction. Toward this end, we modeled the complexes of the pBCR and SH2GRB2 and other pY/Y-peptide-SH2 complexes and compared their specificity and affinity. We observed that pBCR's 176FpYVNV180 motif is favorable and specific to SH2GRB2, similar to pEGFR, but not other complexes. SH2GRB2 contains two binding pockets: pY-binding recognition pocket triggers binding, and the specificity pocket whose interaction is governed by N179 in pBCR and W121 in SH2GRB2. Our proposed motif with optimal affinity to SH2GRB2 is E/D-pY-E/V-N-I/L. Collectively, we provide the structural basis of BCR-ABL recruitment of GRB2, outline its specificity hallmarks, and delineate a blueprint for prediction of BCR-binding scaffolds and for therapeutic peptide design.

Mapping the conformational energy landscape of Abl kinase using ClyA nanopore tweezers

Protein kinases play central roles in cellular regulation by catalyzing the phosphorylation of target proteins. Kinases have inherent structural flexibility allowing them to switch between active and inactive states. Quantitative characterization of kinase conformational dynamics is challenging. Here, we use nanopore tweezers to assess the conformational dynamics of Abl kinase domain, which is shown to interconvert between two major conformational states where one conformation comprises three sub-states. Analysis of kinase-substrate and kinase-inhibitor interactions uncovers the functional roles of relevant states and enables the elucidation of the mechanism underlying the catalytic deficiency of an inactive Abl mutant G321V. Furthermore, we obtain the energy landscape of Abl kinase by quantifying the population and transition rates of the conformational states. These results extend the view on the dynamic nature of Abl kinase and suggest nanopore tweezers can be used as an efficient tool for other members of the human kinome.

Quantitative characterization of kinase conformational dynamics remains challenging. Here, the authors show that protein nanopore tweezers allow analyzing the conformational energy landscape and ligand binding of the Abl kinase domain.

Towards autonomous analysis of chemical exchange saturation transfer experiments using deep neural networks

Macromolecules often exchange between functional states on timescales that can be accessed with NMR spectroscopy and many NMR tools have been developed to characterise the kinetics and thermodynamics of the exchange processes, as well as the structure of the conformers that are involved. However, analysis of the NMR data that report on exchanging macromolecules often hinges on complex least-squares fitting procedures as well as human experience and intuition, which, in some cases, limits the widespread use of the methods. The applications of deep neural networks (DNNs) and artificial intelligence have increased significantly in the sciences, and recently, specifically, within the field of biomolecular NMR, where DNNs are now available for tasks such as the reconstruction of sparsely sampled spectra, peak picking, and virtual decoupling. Here we present a DNN for the analysis of chemical exchange saturation transfer (CEST) data reporting on two- or three-site chemical exchange involving sparse state lifetimes of between approximately 3-60 ms, the range most frequently observed via experiment. The work presented here focuses on the ¹H CEST class of methods that are further complicated, in relation to applications to other nuclei, by anti-phase features. The developed DNNs accurately predict the chemical shifts of nuclei in the exchanging species directly from anti-phase ¹H^N CEST profiles, along with an uncertainty associated with the predictions. The performance of the DNN was quantitatively assessed using both synthetic and experimental anti-phase CEST profiles. The assessments show that the DNN accurately determines chemical shifts and their associated uncertainties. The DNNs developed here do not contain any parameters for the end-user to adjust and the method therefore allows for autonomous analysis of complex NMR data that report on conformational exchange.

Elucidating the mechanisms underlying protein conformational switching using NMR spectroscopy

How proteins switch between various ligand-free and ligand-bound structures has been a key biophysical question ever since the postulation of the Monod-Wyman-Changeux and Koshland-Nemethy-Filmer models over six decades ago. The ability of NMR spectroscopy to provide structural and kinetic information on biomolecular conformational exchange places it in a unique position as an analytical tool to interrogate the mechanisms of biological processes such as protein folding and biomolecular complex formation. In addition, recent methodological developments in the areas of saturation transfer and relaxation dispersion have expanded the scope of NMR for probing the mechanics of transitions in systems where one or more states constituting the exchange process are sparsely populated and 'invisible' in NMR spectra. In this review, we highlight some of the strategies available from NMR spectroscopy for examining the nature of multi-site conformational exchange, using five case studies that have employed NMR, either in isolation, or in conjunction with other biophysical tools.

Molecular dynamics simulations of the conformational plasticity in the active pocket of salt-inducible kinase 2 (SIK2) multi-state binding with bosutinib

The kinase domain is highly conserved among protein kinases 'in terms of both sequence and structure. Conformational rearrangements of the kinase domain are affected by the phosphorylation of residues and the binding of kinase inhibitors. Interestingly, the conformational rearrangement of the active pocket plays an important role in kinase activity and can be used to design novel kinase inhibitors. We characterized the conformational plasticity of the active pocket when bosutinib was bound to salt-inducible kinase 2 (SIK2) using homology modeling and molecular dynamics simulations. Ten different initial complex models were constructed using the Morph server, ranging from open to closed conformations of SIK2 binding with bosutinib. Our simulation showed that bosutinib binds SIK2 with up or down conformations of the P-loop and with all the conformations of the activation loop. In addition, the αC-helix conformation was induced by the conformation of the activation loop, and the salt bridge formed only with its open conformation. The binding affinity of the models was also determined using the molecular mechanics generalized Born surface area method. Bosutinib was found to form a strong binding model with SIK2 and hydrophobic interactions were the dominant factor. This discovery may help guide the design of novel SIK2 inhibitors.

Imatinib disassembles the regulatory core of Abelson kinase by binding to its ATP site and not by binding to its myristoyl pocket

It was recently reported (Xie et al., 2022) that the Abelson tyrosine kinase (Abl) ATP-site inhibitor imatinib also binds to Abl's myristoyl binding pocket, which is the target of allosteric Abl inhibitors. This was based on a crystal structure of a truncated Abl kinase domain construct in complex with imatinib bound to the allosteric site as well as further isothermal titration calorimetry (ITC), NMR, and kinase activity data. Although imatinib's affinity for the allosteric site is significantly weaker (10 µ M) than for the ATP site (10 nM), imatinib binding to the allosteric site may disassemble the regulatory core of Abl, thereby stimulating kinase activity, in particular for Abl mutants with reduced imatinib ATP-site affinity. It was argued that the previously observed imatinib-induced opening of the Abl regulatory core (Skora et al., 2013; Sonti et al., 2018) may be caused by the binding of imatinib to the allosteric site and not to the ATP site. We show here that this is not the case but that indeed imatinib binding to the ATP site induces the opening of the regulatory core at nanomolar concentrations. This agrees with findings that other type-II ATP-site inhibitors (nilotinib, ponatinib) disassemble the regulatory core despite demonstrated negligible binding to the allosteric site.

Validation of an allosteric binding site of Src kinase identified by unbiased ligand binding simulations

Allostery plays a primary role in regulating protein activity, making it an important mechanism in human disease and drug discovery. Identifying allosteric regulatory sites to explore their biological significance and therapeutic potential is invaluable to drug discovery; however, identification remains a challenge. Allosteric sites are often "cryptic" without clear geometric or chemical features. Since allosteric regulatory sites are often less conserved in protein kinases than the orthosteric ATP binding site, allosteric ligands are commonly more specific than ATP competitive inhibitors. We present a generalizable computational protocol to predict allosteric ligand binding sites based on unbiased ligand binding simulation trajectories. We demonstrate the feasibility of this protocol by revisiting our previously published ligand binding simulations using the first identified viral proto-oncogene, Src kinase, as a model system. The binding paths for kinase inhibitor PP1 uncovered three metastable intermediate states before binding the high-affinity ATP-binding pocket, revealing two previously known allosteric sites and one novel site. Herein, we validate the novel site using a combination of virtual screening and experimental assays to identify a v-type allosteric small-molecule inhibitor that targets this novel site with specificity for Src over closely related kinases. This study provides a proof-of-concept for employing unbiased ligand binding simulations to identify cryptic allosteric binding sites and is widely applicable to other protein-ligand systems.

Dynamics of protein kinases and pseudokinases by HDX-MS

Dynamics of the protein kinase fold are deeply intertwined with its structure. The past three decades of kinase biophysical studies revealed key dynamic features of the kinase domain and, more recently, how these features may endow catalytically impaired kinases-or pseudokinases-with signaling properties. Hydrogen-deuterium exchange coupled with mass spectrometry (HDX-MS) is proving to be a valuable approach for studies of kinase and pseudokinase domain dynamics. Here, we briefly discuss the methods that have provided insights into protein kinase dynamics, describe how HDX-MS is being used to answer questions in the kinase/pseudokinase field, and provide a detailed protocol for collecting an HDX-MS dataset to study the impacts of small molecule binding to a pseudokinase domain. As more small molecules are discovered that can disrupt pseudokinase conformations, HDX-MS is likely to be a powerful approach for exploring drug-induced changes in pseudokinase dynamics and structure.

The mechanism of activation of MEK1 by B-Raf and KSR1

MEK1 interactions with B-Raf and KSR1 are key steps in Ras/Raf/MEK/ERK signaling. Despite this, vital mechanistic details of how these execute signal transduction are still enigmatic. Among these is why, despite B-Raf and KSR1 kinase domains similarity, the B-Raf/MEK1 and KSR1/MEK1 complexes have distinct contributions to MEK1 activation, and broadly, what is KSR1's role. Our molecular dynamics simulations clarify these still unresolved ambiguities. Our results reveal that the proline-rich (P-rich) loop of MEK1 plays a decisive role in MEK1 activation loop (A-loop) phosphorylation. In the inactive B-Raf/MEK1 heterodimer, the collapsed A-loop of B-Raf interacts with the P-rich loop and A-loop of MEK1, minimizing MEK1 A-loop fluctuation and preventing it from phosphorylation. In the active B-Raf/MEK1 heterodimer, the P-rich loop moves in concert with the A-loop of B-Raf as it extends. This reduces the number of residues interacting with MEK1 A-loop, allowing increased A-loop fluctuation, and bringing Ser222 closer to ATP for phosphorylation. B-Raf αG-helix Arg662 promotes MEK1 activation by orienting Ser218 towards ATP. In KSR1/MEK1, the KSR1 αG-helix has Ala826 in place of B-Raf Arg662. This difference results in much fewer interactions between KSR1 αG-helix and MEK1 A-loop, thus a more flexible A-loop. We postulate that if KSR1 were to adopt an active configuration with an extended A-loop as seen in other protein kinases, then the MEK1 P-rich loop would extend in a similar manner, as seen in the active B-Raf/MEK1 heterodimer. This would result in highly flexible MEK1 A-loop, and KSR1 functioning as an active, B-Raf-like, kinase.

A topological data analytic approach for discovering biophysical signatures in protein dynamics

Identifying structural differences among proteins can be a non-trivial task. When contrasting ensembles of protein structures obtained from molecular dynamics simulations, biologically-relevant features can be easily overshadowed by spurious fluctuations. Here, we present SINATRA Pro, a computational pipeline designed to robustly identify topological differences between two sets of protein structures. Algorithmically, SINATRA Pro works by first taking in the 3D atomic coordinates for each protein snapshot and summarizing them according to their underlying topology. Statistically significant topological features are then projected back onto a user-selected representative protein structure, thus facilitating the visual identification of biophysical signatures of different protein ensembles. We assess the ability of SINATRA Pro to detect minute conformational changes in five independent protein systems of varying complexities. In all test cases, SINATRA Pro identifies known structural features that have been validated by previous experimental and computational studies, as well as novel features that are also likely to be biologically-relevant according to the literature. These results highlight SINATRA Pro as a promising method for facilitating the non-trivial task of pattern recognition in trajectories resulting from molecular dynamics simulations, with substantially increased resolution.

Structural features of proteins often serve as signatures of their biological function and molecular binding activity. Elucidating these structural features is essential for a full understanding of underlying biophysical mechanisms. While there are existing methods aimed at identifying structural differences between protein variants, such methods do not have the capability to jointly infer both geometric and dynamic changes, simultaneously. In this paper, we propose SINATRA Pro, a computational framework for extracting key structural features between two sets of proteins. SINATRA Pro robustly outperforms standard techniques in pinpointing the physical locations of both static and dynamic signatures across various types of protein ensembles, and it does so with improved resolution.

Phosphorylation-Induced Conformational Dynamics and Inhibition of Janus Kinase 1 by Suppressors of Cytokine Signaling 1

The dysfunction of the JAK/STAT (Janus kinase/signal transducers and activators of transcription) pathway results in several pathophysiological conditions, including autoimmune disorders. The negative feedback regulators of the JAK/STAT signaling pathway, suppressors of cytokine signaling (SOCS), act as a natural inhibitor of JAK and inhibit aberrant activity. SOCS1 is the most potent member of the SOCS family, whose kinase inhibitory region targets the substrate-binding groove of JAK with high affinity and blocks the phosphorylation of JAK kinases. Overall, we performed an aggregate of 13 μs molecular dynamics simulations on the activation loop's three different phosphorylation (double and single) states. Results from our simulations show that the single Tyr¹⁰³⁴ phosphorylation could stabilize the JAK1/SOCS1 complex as well as the flexible activation segment. The phosphate-binding loop (P-loop) shows conformational variability at dual and single phosphorylated states. Principal component analysis and protein structure network (PSN) analysis reveal that the different phosphorylation states and SOCS1 binding induce intermediate inactive conformations of JAK1, which could be a better target for future JAK1 selective drug design. PSN analysis suggests that the com-pY1034 system is stabilized due to higher values of network hubs than the other two complex systems. Moreover, the binding free energy calculations suggest that pTyr¹⁰³⁴ states show a higher affinity toward SOCS1 than the dual and pTyr¹⁰³⁵ states. We believe that the mechanistic understanding of JAK1/SOCS1 complexation will aid future studies related to peptide inhibitors based on SOCS1.

In-Cell Structural Biology by NMR: The Benefits of the Atomic Scale

In-cell structural biology aims at extracting structural information about proteins or nucleic acids in their native, cellular environment. This emerging field holds great promise and is already providing new facts and outlooks of interest at both fundamental and applied levels. NMR spectroscopy has important contributions on this stage: It brings information on a broad variety of nuclei at the atomic scale, which ensures its great versatility and uniqueness. Here, we detail the methods, the fundamental knowledge, and the applications in biomedical engineering related to in-cell structural biology by NMR. We finally propose a brief overview of the main other techniques in the field (EPR, smFRET, cryo-ET, etc.) to draw some advisable developments for in-cell NMR. In the era of large-scale screenings and deep learning, both accurate and qualitative experimental evidence are as essential as ever to understand the interior life of cells. In-cell structural biology by NMR spectroscopy can generate such a knowledge, and it does so at the atomic scale. This review is meant to deliver comprehensive but accessible information, with advanced technical details and reflections on the methods, the nature of the results, and the future of the field.

Assembly mechanism of early Hsp90-Cdc37-kinase complexes

Molecular chaperones have an essential role for the maintenance of a balanced protein homeostasis. Here, we investigate how protein kinases are recruited and loaded to the Hsp90-Cdc37 complex, the first step during Hsp90-mediated chaperoning that leads to enhanced client kinase stability and activation. We show that conformational dynamics of all partners is a critical feature of the underlying loading mechanism. The kinome co-chaperone Cdc37 exists primarily in a dynamic extended conformation but samples a low-populated, well-defined compact structure. Exchange between these two states is maintained in an assembled Hsp90-Cdc37 complex and is necessary for substrate loading. Breathing motions at the N-lobe of a free kinase domain partially expose the kinase segment trapped in the Hsp90 dimer downstream in the cycle. Thus, client dynamics poise for chaperone dependence. Hsp90 is not directly involved during loading, and Cdc37 is assigned the task of sensing clients by stabilizing the preexisting partially unfolded client state.

AtomNet PoseRanker: Enriching Ligand Pose Quality for Dynamic Proteins in Virtual High-Throughput Screens

Structure-based, virtual High-Throughput Screening (vHTS) methods for predicting ligand activity in drug discovery are important when there are no or relatively few known compounds that interact with a therapeutic target of interest. State-of-the-art computational vHTS necessarily relies on effective methods for pose sampling and docking and generating an accurate affinity score from the docked poses. However, proteins are dynamic; in vivo ligands bind to a conformational ensemble. In silico docking to the single conformation represented by a crystal structure can adversely affect the pose quality. Here, we introduce AtomNet PoseRanker (ANPR), a graph convolutional network trained to identify and rerank crystal-like ligand poses from a sampled ensemble of protein conformations and ligand poses. In contrast to conventional vHTS methods that incorporate receptor flexibility, a deep learning approach can internalize valid cognate and noncognate binding modes corresponding to distinct receptor conformations, thereby learning to infer and account for receptor flexibility even on single conformations. ANPR significantly enriched pose quality in docking to cognate and noncognate receptors of the PDBbind v2019 data set. Improved pose rankings that better represent experimentally observed ligand binding modes improve hit rates in vHTS campaigns and thereby advance computational drug discovery, especially for novel therapeutic targets or novel binding sites.

Imatinib can act as an Allosteric Activator of Abl Kinase

Imatinib is an ATP-competitive inhibitor of Bcr-Abl kinase and the first drug approved for chronic myelogenous leukemia (CML) treatment. Here we show that imatinib binds to a secondary, allosteric site located in the myristoyl pocket of Abl to function as an activator of the kinase activity. Abl transitions between an assembled, inhibited state and an extended, activated state. The equilibrium is regulated by the conformation of the αΙ helix, which is located nearby the allosteric pocket. Imatinib binding to the allosteric pocket elicits an αΙ helix conformation that is not compatible with the assembled state, thereby promoting the extended state and stimulating the kinase activity. Although in wild-type Abl the catalytic pocket has a much higher affinity for imatinib than the allosteric pocket does, the two binding affinities are comparable in Abl variants carrying imatinib-resistant mutations in the catalytic site. A previously isolated imatinib-resistant mutation in patients appears to be mediating its function by increasing the affinity of imatinib for the allosteric pocket, providing a hitherto unknown mechanism of drug resistance. Our results highlight the benefit of combining imatinib with allosteric inhibitors to maximize their inhibitory effect on Bcr-Abl.

SCP4-STK35/PDIK1L complex is a dual phospho-catalytic signaling dependency in acute myeloid leukemia

Acute myeloid leukemia (AML) cells rely on phospho-signaling pathways to gain unlimited proliferation potential. Here, we use domain-focused CRISPR screening and identify the nuclear phosphatase SCP4 as a dependency in AML, yet this enzyme is dispensable in normal hematopoietic progenitor cells. Using CRISPR exon scanning and gene complementation assays, we show that the catalytic function of SCP4 is essential in AML. Through mass spectrometry analysis of affinity-purified complexes, we identify the kinase paralogs STK35 and PDIK1L as binding partners and substrates of the SCP4 phosphatase domain. We show that STK35 and PDIK1L function catalytically and redundantly in the same pathway as SCP4 to maintain AML proliferation and to support amino acid biosynthesis and transport. We provide evidence that SCP4 regulates STK35/PDIK1L through two distinct mechanisms: catalytic removal of inhibitory phosphorylation and by promoting kinase stability. Our findings reveal a phosphatase-kinase signaling complex that supports the pathogenesis of AML.

Progress toward automated methyl assignments for methyl-TROSY applications

Methyl-TROSY spectroscopy has extended the reach of solution-state NMR to supra-molecular machineries over 100 kDa in size. Methyl groups are ideal probes for studying structure, dynamics, and protein-protein interactions in quasi-physiological conditions with atomic resolution. Successful implementation of the methodology requires accurate methyl chemical shift assignment, and the task still poses a significant challenge in the field. In this work, we outline the current state of technology for methyl labeling, data collection, data analysis, and nuclear Overhauser effect (NOE)-based automated methyl assignment approaches. We present MAGIC-Act and MAGIC-View, two Python extensions developed as part of the popular NMRFAM-Sparky package, and MAGIC-Net a standalone structure-based network analysis program. MAGIC-Act conducts statistically driven amino acid typing, Leu/Val pairing guided by 3D HMBC-HMQC, and NOESY cross-peak symmetry checking. MAGIC-Net provides model-based NOE statistics to aid in selection of a methyl labeling scheme. The programs provide a versatile, semi-automated framework for rapid methyl assignment.

Influence of noncovalent intramolecular and host–guest interactions on imatinib binding to MoS2 sheets: a PXRD/DFT study

The revealed pattern of imatinib drug binding to MoS2 sheets is promising for the combined exploitation of these species for therapeutic purposes.

Mechanism for the activation of the anaplastic lymphoma kinase receptor

Anaplastic lymphoma kinase (ALK) is a receptor tyrosine kinase (RTK) that regulates important functions in the central nervous system^1,2. The ALK gene is a hotspot for chromosomal translocation events that result in several fusion proteins that cause a variety of human malignancies³. Somatic and germline gain-of-function mutations in ALK were identified in paediatric neuroblastoma^4-7. ALK is composed of an extracellular region (ECR), a single transmembrane helix and an intracellular tyrosine kinase domain^8,9. ALK is activated by the binding of ALKAL1 and ALKAL2 ligands^10-14 to its ECR, but the lack of structural information for the ALK-ECR or for ALKAL ligands has limited our understanding of ALK activation. Here we used cryo-electron microscopy, nuclear magnetic resonance and X-ray crystallography to determine the atomic details of human ALK dimerization and activation by ALKAL1 and ALKAL2. Our data reveal a mechanism of RTK activation that allows dimerization by either dimeric (ALKAL2) or monomeric (ALKAL1) ligands. This mechanism is underpinned by an unusual architecture of the receptor-ligand complex. The ALK-ECR undergoes a pronounced ligand-induced rearrangement and adopts an orientation parallel to the membrane surface. This orientation is further stabilized by an interaction between the ligand and the membrane. Our findings highlight the diversity in RTK oligomerization and activation mechanisms.

Opening of a cryptic pocket in β-lactamase increases penicillinase activity

A protein is a shape-shifter, but it is currently unclear which of the many structures a protein can adopt are relevant for its function. Here, we examine conformations that contain a “cryptic” pocket (i.e., a pocket absent in ligand-free structures). Cryptic pockets have potential utility in drug discovery efforts because they provide a means to target “undruggable” proteins (i.e., proteins lacking known pockets) or enhance rather than inhibit protein function. In this study, we use a combination of thiol-labeling and kinetic assays, NMR, and molecular dynamic simulations to identify the function of the Ω-loop cryptic pocket in β-lactamase enzymes. We find that an open pocket population is beneficial for hydrolysis of the substrate benzylpenicillin.

Understanding the functional role of protein-excited states has important implications in protein design and drug discovery. However, because these states are difficult to find and study, it is still unclear if excited states simply result from thermal fluctuations and generally detract from function or if these states can actually enhance protein function. To investigate this question, we consider excited states in β-lactamases and particularly a subset of states containing a cryptic pocket which forms under the Ω-loop. Given the known importance of the Ω-loop and the presence of this pocket in at least two homologs, we hypothesized that these excited states enhance enzyme activity. Using thiol-labeling assays to probe Ω-loop pocket dynamics and kinetic assays to probe activity, we find that while this pocket is not completely conserved across β-lactamase homologs, those with the Ω-loop pocket have a higher activity against the substrate benzylpenicillin. We also find that this is true for TEM β-lactamase variants with greater open Ω-loop pocket populations. We further investigate the open population using a combination of NMR chemical exchange saturation transfer experiments and molecular dynamics simulations. To test our understanding of the Ω-loop pocket’s functional role, we designed mutations to enhance/suppress pocket opening and observed that benzylpenicillin activity is proportional to the probability of pocket opening in our designed variants. The work described here suggests that excited states containing cryptic pockets can be advantageous for function and may be favored by natural selection, increasing the potential utility of such cryptic pockets as drug targets.

Mutation in Abl kinase with altered drug-binding kinetics indicates a novel mechanism of imatinib resistance

Protein kinase inhibitors are potent anticancer therapeutics. For example, the Bcr-Abl kinase inhibitor imatinib decreases mortality for chronic myeloid leukemia by 80%, but 22 to 41% of patients acquire resistance to imatinib. About 70% of relapsed patients harbor mutations in the Bcr-Abl kinase domain, where more than a hundred different mutations have been identified. Some mutations are located near the imatinib-binding site and cause resistance through altered interactions with the drug. However, many resistance mutations are located far from the drug-binding site, and it remains unclear how these mutations confer resistance. Additionally, earlier studies on small sets of patient-derived imatinib resistance mutations indicated that some of these mutant proteins were in fact sensitive to imatinib in cellular and biochemical studies. Here, we surveyed the resistance of 94 patient-derived Abl kinase domain mutations annotated as disease relevant or resistance causing using an engagement assay in live cells. We found that only two-thirds of mutations weaken imatinib affinity by more than twofold compared to Abl wild type. Surprisingly, one-third of mutations in the Abl kinase domain still remain sensitive to imatinib and bind with similar or higher affinity than wild type. Intriguingly, we identified three clinical Abl mutations that bind imatinib with wild type-like affinity but dissociate from imatinib considerably faster. Given the relevance of residence time for drug efficacy, mutations that alter binding kinetics could cause resistance in the nonequilibrium environment of the body where drug export and clearance play critical roles.

The A39G FF domain folds on a volcano-shaped free energy surface via separate pathways

Conformational dynamics play critical roles in protein folding, misfolding, function, misfunction, and aggregation. While detecting and studying the different conformational states populated by protein molecules on their free energy surfaces (FESs) remain a challenge, NMR spectroscopy has emerged as an invaluable experimental tool to explore the FES of a protein, as conformational dynamics can be probed at atomic resolution over a wide range of timescales. Here, we use chemical exchange saturation transfer (CEST) to detect "invisible" minor states on the energy landscape of the A39G mutant FF domain that exhibited "two-state" folding kinetics in traditional experiments. Although CEST has mostly been limited to studies of processes with rates between ∼5 to 300 s^-1 involving sparse states with populations as low as ∼1%, we show that the line broadening that is often associated with minor state dips in CEST profiles can be exploited to inform on additional conformers, with lifetimes an order of magnitude shorter and populations close to 10-fold smaller than what typically is characterized. Our analysis of CEST profiles that exploits the minor state linewidths of the 71-residue A39G FF domain establishes a folding mechanism that can be described in terms of a four-state exchange process between interconverting states spanning over two orders of magnitude in timescale from ∼100 to ∼15,000 μs. A similar folding scheme is established for the wild-type domain as well. The study shows that the folding of this small domain proceeds through a pair of sparse, partially structured intermediates via two discrete pathways on a volcano-shaped FES.

Activation of Abl1 Kinase Explored Using Well-Tempered Metadynamics Simulations on an Essential Dynamics Sampled Path

Well-tempered metadynamics (wT-metaD) simulations using path collective variables (CVs) have been successfully applied in recent years to explore conformational transitions in protein kinases and other biomolecular systems. While this methodology has the advantage of describing the transitions with a limited number of predefined path CVs, it requires as an input a reference path connecting the initial and target states of the system. It is desirable to automate the path generation using approaches that do not rely on the choice of geometric CVs to describe the transition of interest. To this end, we developed an approach that couples essential dynamics sampling with wT-metaD simulations. We used this newly developed procedure to explore the activation mechanism of Abl1 kinase and compute the associated free energy barriers. Through these simulations, we identified a three-step mechanism for the activation that involved two metastable intermediates that possessed a partially open activation loop and differed primarily in the “in” or “out” conformation of the aspartate residue of the DFG motif. One of these states is similar to a conformation that was detected in previous spectroscopic studies of Abl1 kinase, albeit its mechanistic role in the activation was hitherto not well understood. The present study establishes its intermediary role in the activation and predicts a rate-determining free energy barrier of 13.8 kcal/mol that is in good agreement with previous experimental and computational estimates. Overall, our study demonstrates the usability of essential dynamics sampling as a path CV in wT-metaD to conveniently study conformational transitions and accurately calculate the associated barriers.

Structure, Function, and Thermodynamics of Lactate Dehydrogenases from Humans and the Malaria Parasite P. falciparum

Temperature adaptation is ubiquitous among all living organisms, yet the molecular basis for this process remains poorly understood. It can be assumed that for parasite-host systems, the same enzymes found in both organisms respond to the same selection factor (human body temperature) with similar structural changes. Herein, we report the existence of a reversible temperature-dependent structural transition for the glycolytic enzyme lactate dehydrogenase (LDH) from the malaria parasite Plasmodium falciparum (pfLDH) and human heart (hhLDH) occurring in the temperature range of human fever. This transition is observed for LDHs from psychrophiles, mesophiles, and moderate thermophiles in their operating temperature range. Thermodynamic analysis reveals unique thermodynamic signatures of the LDH-substrate complexes defining a specific temperature range to which human LDH is adapted and parasite LDH is not, despite their common mesophilic nature. The results of spectroscopic analysis combined with the available crystallographic data reveal the existence of an active center within pfLDH that imparts psychrophilic structural properties to the enzyme. This center consists of two pockets, one formed by the five amino acids (5AA insert) within the substrate specificity loop and the other by the active site, that mutually regulate one another in response to temperature and induce structural and functional changes in the Michaelis complex. Our findings pave the way toward a new strategy for malaria treatments and drug design using therapeutic agents that inactivate malarial LDH selectively at a specific temperature range of the cyclic malaria paroxysm.

AlphaFold 2: Why It Works and Its Implications for Understanding the Relationships of Protein Sequence, Structure, and Function

AlphaFold 2 (AF2) was the star of CASP14, the last biannual structure prediction experiment. Using novel deep learning, AF2 predicted the structures of many difficult protein targets at or near experimental resolution. Here, we present our perspective of why AF2 works and show that it is a very sophisticated fold recognition algorithm that exploits the completeness of the library of single domain PDB structures. It has also learned local side chain packing rearrangements that enable it to refine proteins to high resolution. The benefits and limitations of its ability to predict the structures of many more proteins at or close to atomic detail are discussed.

TYROSINE KINASES: COMPLEX MOLECULAR SYSTEMS CHALLENGING COMPUTATIONAL METHODOLOGIES

Classical molecular dynamics (MD) simulations based on atomic models play an increasingly important role in a wide range of applications in physics, biology, and chemistry. Nonetheless, generating genuine knowledge about biological systems using MD simulations remains challenging. Protein tyrosine kinases are important cellular signaling enzymes that regulate cell growth, proliferation, metabolism, differentiation, and migration. Due to the large conformational changes and long timescales involved in their function, these kinases present particularly challenging problems to modern computational and theoretical frameworks aimed at elucidating the dynamics of complex biomolecular systems. Markov state models have achieved limited success in tackling the broader conformational ensemble and biased methods are often employed to examine specific long timescale events. Recent advances in machine learning continue to push the limitations of current methodologies and provide notable improvements when integrated with the existing frameworks. A broad perspective is drawn from a critical review of recent studies.

Significant Loop Motions in the SsoPTP Protein Tyrosine Phosphatase Allow for Dual General Acid Functionality

Conformational dynamics are important factors in the function of enzymes, including protein tyrosine phosphatases (PTPs). Crystal structures of PTPs first revealed the motion of a protein loop bearing a conserved catalytic aspartic acid, and subsequent nuclear magnetic resonance and computational analyses have shown the presence of motions, involved in catalysis and allostery, within and beyond the active site. The tyrosine phosphatase from the thermophilic and acidophilic Sulfolobus solfataricus (SsoPTP) displays motions of its acid loop together with dynamics of its phosphoryl-binding P-loop and the Q-loop, the first instance of such motions in a PTP. All three loops share the same exchange rate, implying their motions are coupled. Further evidence of conformational flexibility comes from mutagenesis, kinetics, and isotope effect data showing that E40 can function as an alternate general acid to protonate the leaving group when the conserved acid, D69, is mutated to asparagine. SsoPTP is not the first PTP to exhibit an alternate general acid (after VHZ and TkPTP), but E40 does not correspond to the sequence or structural location of the alternate general acids in those precedents. A high-resolution X-ray structure with the transition state analogue vanadate clarifies the role of the active site arginine R102, which varied in structures of substrates bound to a catalytically inactive mutant. The coordinated motions of all three functional loops in SsoPTP, together with the function of an alternate general acid, suggest that catalytically competent conformations are present in solution that have not yet been observed in crystal structures.

Minimizing Pervasive Artifacts in 4D Covariance Maps for Protein Side Chain NMR Assignments

Nuclear magnetic resonance (NMR) is a mainstay of biophysical studies that provides atomic level readouts to formulate molecular mechanisms. Side chains are particularly important to derive mechanisms involving proteins as they carry functional groups, but NMR studies of side chains are often limited by challenges in assigning their signals. Here, we designed a novel computational method that combines spectral derivatives and matrix square-rooting to produce reliable 4D covariance maps from routinely acquired 3D spectra and facilitates side chain resonance assignments. Thus, we generate two 4D maps from 3D-HcccoNH and 3D-HCcH-TOCSY spectra that each help overcome signal overlap or sensitivity losses. These 4D maps feature HC-HSQCs of individual side chains that can be paired to assigned backbone amide resonances of individual aliphatic signals, and both are obtained from a single modified covariance calculation. Further, we present 4D maps produced using conventional triple resonance experiments to easily assign asparagine side chain amide resonances. The 4D covariance maps encapsulate the lengthy manual pattern recognition used in traditional assignment methods and distill the information as correlations that can be easily visualized. We showcase the utility of the 4D covariance maps with a 10 kDa peptidyl carrier protein and a 52 kDa cyclization domain from a nonribosomal peptide synthetase.

Mutations in PBP2 from ceftriaxone-resistant Neisseria gonorrhoeae alter the dynamics of the β3-β4 loop to favor a low-affinity drug-binding state

Resistance to the extended-spectrum cephalosporin ceftriaxone in the pathogenic bacteria Neisseria gonorrhoeae is conferred by mutations in penicillin-binding protein 2 (PBP2), the lethal target of the antibiotic, but how these mutations exert their effect at the molecular level is unclear. Using solution NMR, X-ray crystallography, and isothermal titration calorimetry, we report that WT PBP2 exchanges dynamically between a low-affinity state with an extended β3–β4 loop conformation and a high-affinity state with an inward β3–β4 loop conformation. Histidine-514, which is located at the boundary of the β4 strand, plays an important role during the exchange between these two conformational states. We also find that mutations present in PBP2 from H041, a ceftriaxone-resistant strain of N. gonorrhoeae, increase resistance to ceftriaxone by destabilizing the inward β3–β4 loop conformation or stabilizing the extended β3–β4 loop conformation to favor the low-affinity drug-binding state. These observations reveal a unique mechanism for ceftriaxone resistance, whereby mutations in PBP2 lower the proportion of target molecules in the high-affinity drug-binding state and thus reduce inhibition at lower drug concentrations.

Probing allosteric regulations with coevolution-driven molecular simulations

Protein-mediated allosteric regulations are essential in biology, but their quantitative characterization continues to posit formidable challenges for both experiments and computations. Here, we combine coevolutionary information, multiscale molecular simulations, and free-energy methods to interrogate and quantify the allosteric regulation of functional changes in protein complexes. We apply this approach to investigate the regulation of adenylyl cyclase (AC) by stimulatory and inhibitory G proteins—a prototypical allosteric system that has long escaped from in-depth molecular characterization. We reveal a surprisingly simple ON/OFF regulation of AC functional dynamics through multiple pathways of information transfer. The binding of G proteins reshapes the free-energy landscape of AC following the classical population-shift paradigm. The model agrees with structural and biochemical data and reveals previously unknown experimentally consistent intermediates. Our approach showcases a general strategy to explore uncharted functional space in complex biomolecular regulations.

Measuring radiofrequency fields in NMR spectroscopy using offset-dependent nutation profiles

The application of NMR spectroscopy for studying molecular and reaction dynamics relies crucially on the measurement of the magnitude of radiofrequency (RF) fields that are used to nutate or lock the nuclear magnetization. Here, we report a method for measuring RF field amplitudes that leverages the intrinsic modulations observed in offset-dependent NMR nutation profiles of small molecules. Such nutation profiles are exquisitely sensitive to the magnitude of the RF field, and B₁ values ranging from 1 to 2000 Hz, as well the inhomogeneity in B₁ distributions, can be determined with high accuracy and precision using this approach. In order to measure B₁ fields associated with NMR experiments carried out on protein or nucleic acids, where these modulations are obscured by the large transverse relaxation rate constants of the analyte, our approach can be used in conjunction with a suitable external small molecule standard, expanding the scope of the method for large biomolecules.

(Phenylamino)pyrimidine-1,2,3-triazole derivatives as analogs of imatinib: searching for novel compounds against chronic myeloid leukemia

The enzyme tyrosine kinase BCR-Abl-1 is the main molecular target in the treatment of chronic myeloid leukemia and can be competitively inhibited by tyrosine kinase inhibitors such as imatinib. New potential competitive inhibitors were synthesized using the (phenylamino)pyrimidine-pyridine (PAPP) group as a pharmacophoric fragment, and these compounds were biologically evaluated. The synthesis of twelve new compounds was performed in three steps and assisted by microwave irradiation in a 1,3-dipolar cycloaddition to obtain 1,2,3-triazole derivatives substituted on carbon C-4 of the triazole nucleus. All compounds were evaluated for their inhibitory activities against a chronic myeloid leukemia cell line (K562) that expresses the enzyme tyrosine kinase BCR-Abl-1 and against healthy cells (WSS-1) to observe their selectivity. Three compounds showed promising results, with IC50 values between 1.0 and 7.3 μM, and were subjected to molecular docking studies. The results suggest that such compounds can interact at the same binding site as imatinib, probably sharing a competitive inhibition mechanism. One compound showed the greatest interaction affinity for BCR-Abl-1 in the docking studies.

Dynamic equilibria in protein kinases

Structural changes involved in protein kinase activation and ligand binding have been determined from a wealth of X-ray crystallographic evidence. Recent solution studies using NMR, EPR, HX-MS, and fluorescence techniques have deepened this understanding by highlighting the underlying energetics and dynamics of multistate conformational ensembles. This new research is showing how activation mechanisms and ligand binding alter the internal motions of kinases and enable allosteric coupling between distal regulatory regions and the active site.

Synthesis of Isotopically Labeled, Spin-Isolated Tyrosine and Phenylalanine for Protein NMR Applications

Isotopically labeled amino acids are widely used to study the structure and dynamics of proteins by NMR. Herein we describe a facile, gram-scale synthesis of compounds 1b and 2b under standard laboratory conditions from the common intermediate 7. 2b is obtained via simple deprotection, while 1b is accessed through a reductive deoxygenation/deuteration sequence and deprotection. 1b and 2b provide improved signal intensity using lower amounts of labeled precursor and are alternatives to existing labeling approaches.

Making the invisible visible: Toward structural characterization of allosteric states, interaction networks, and allosteric regulatory mechanisms in protein kinases

Despite the established view of protein kinases as dynamic and versatile allosteric regulatory machines, our knowledge of allosteric functional states, allosteric interaction networks, and the intrinsic folding energy landscapes is surprisingly limited. We discuss the latest developments in structural characterization of allosteric molecular events underlying protein kinase dynamics and functions using structural, biophysical, and computational biology approaches. The recent studies highlighted progress in making the invisible aspects of protein kinase 'life' visible, including the determination of hidden allosteric states and mapping of allosteric energy landscapes, discovery of new mechanisms underlying ligand-induced modulation of allosteric activity, evolutionary adaptation of kinase allostery, and characterization of allosteric interaction networks as the intrinsic driver of kinase adaptability and signal transmission in the regulatory assemblies.

Reining in BTK: Interdomain Interactions and Their Importance in the Regulatory Control of BTK

Since Dr. Ogden Bruton's 1952 paper describing the first human primary immunodeficiency disease, the peripheral membrane binding signaling protein, aptly named Bruton's tyrosine kinase (BTK), has been the target of intense study. Dr. Bruton's description of agammaglobulinemia set the stage for ultimately understanding key signaling steps emanating from the B cell receptor. BTK is a multidomain tyrosine kinase and in the decades since Dr. Bruton's discovery it has become clear that genetic defects in the regulatory domains or the catalytic domain can lead to immunodeficiency. This finding underscores the intricate regulatory mechanisms within the BTK protein that maintain appropriate levels of signaling both in the resting B cell and during an immune challenge. In recent decades, BTK has become a target for clinical intervention in treating B cell malignancies. The survival reliance of B cell malignancies on B cell receptor signaling has allowed small molecules that target BTK to become essential tools in treating patients with hematological malignancies. The first-in-class Ibrutinib and more selective second-generation inhibitors all target the active site of the multidomain BTK protein. Therapeutic interventions targeting BTK have been successful but are plagued by resistance mutations that render drug treatment ineffective for some patients. This review will examine the molecular mechanisms that drive drug resistance, the long-range conformational effects of active site inhibitors on the BTK regulatory apparatus, and emerging opportunities to allosterically target the BTK kinase to improve therapeutic interventions using combination therapies.

Platelet-derived growth factor receptor beta activates Abl2 via direct binding and phosphorylation

Abl family kinases are nonreceptor tyrosine kinases activated by diverse cellular stimuli that regulate cytoskeleton organization, morphogenesis, and adhesion. The catalytic activity of Abl family kinases is tightly regulated in cells by a complex set of intramolecular and intermolecular interactions and post-translational modifications. For example, the platelet-derived growth factor receptor beta (PDGFRβ), important for cell proliferation and chemotaxis, is a potent activator of Abl family kinases. However, the molecular mechanism by which PDGFRβ engages and activates Abl family kinases is not known. We show here that the Abl2 Src homology 2 domain directly binds to phosphotyrosine Y771 in the PDGFRβ cytoplasmic domain. PDGFRβ directly phosphorylates multiple novel sites on the N-terminal half of Abl2, including Y116, Y139, and Y161 within the Src homology 3 domain, and Y299, Y303, and Y310 on the kinase domain. Y116, Y161, Y272, and Y310 are all located at or near the Src homology 3/Src homology 2-kinase linker interface, which helps maintain Abl family kinases in an autoinhibited conformation. We also found that PDGFRβ-mediated phosphorylation of Abl2 in vitro activates Abl2 kinase activity, but mutation of these four tyrosines (Y116, Y161, Y272, and Y310) to phenylalanine abrogated PDGFRβ-mediated activation of Abl2. These findings reveal how PDGFRβ engages and phosphorylates Abl2 leading to activation of the kinase, providing a framework to understand how growth factor receptors engage and activate Abl family kinases.

Structural analysis of the full-length human LRRK2

Mutations in leucine-rich repeat kinase 2 (LRRK2) are commonly implicated in the pathogenesis of both familial and sporadic Parkinson's disease (PD). LRRK2 regulates critical cellular processes at membranous organelles and forms microtubule-based pathogenic filaments, yet the molecular basis underlying these biological roles of LRRK2 remains largely enigmatic. Here, we determined high-resolution structures of full-length human LRRK2, revealing its architecture and key interdomain scaffolding elements for rationalizing disease-causing mutations. The kinase domain of LRRK2 is captured in an inactive state, a conformation also adopted by the most common PD-associated mutation, LRRK2^G2019S. This conformation serves as a framework for structure-guided design of conformational specific inhibitors. We further determined the structure of COR-mediated LRRK2 dimers and found that single-point mutations at the dimer interface abolished pathogenic filamentation in cells. Overall, our study provides mechanistic insights into physiological and pathological roles of LRRK2 and establishes a structural template for future therapeutic intervention in PD.

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.

Methodological advancements for characterising protein side chains by NMR spectroscopy

The surface of proteins is covered by side chains of polar amino acids that are imperative for modulating protein functionality through the formation of noncovalent intermolecular interactions. However, despite their tremendous importance, the unique structures of protein side chains require tailored approaches for investigation by nuclear magnetic resonance spectroscopy and so have traditionally been understudied compared with the protein backbone. Here, we review substantial recent methodological advancements within nuclear magnetic resonance spectroscopy to address this issue. Specifically, we consider advancements that provide new insight into methyl-bearing side chains, show the potential of using non-natural amino acids and reveal the actions of charged side chains. Combined, the new methods promise unprecedented characterisations of side chains that will further elucidate protein function.

The IDIP framework for assessing protein function and its application to the prion protein

The quest to determine the function of a protein can represent a profound challenge. Although this task is the mandate of countless research groups, a general framework for how it can be approached is conspicuously lacking. Moreover, even expectations for when the function of a protein can be considered to be 'known' are not well defined. In this review, we begin by introducing concepts pertinent to the challenge of protein function assignments. We then propose a framework for inferring a protein's function from four data categories: 'inheritance', 'distribution', 'interactions' and 'phenotypes' (IDIP). We document that the functions of proteins emerge at the intersection of inferences drawn from these data categories and emphasise the benefit of considering them in an evolutionary context. We then apply this approach to the cellular prion protein (PrP^C ), well known for its central role in prion diseases, whose function continues to be considered elusive by many investigators. We document that available data converge on the conclusion that the function of the prion protein is to control a critical post-translational modification of the neural cell adhesion molecule in the context of epithelial-to-mesenchymal transition and related plasticity programmes. Finally, we argue that this proposed function of PrP^C has already passed the test of time and is concordant with the IDIP framework in a way that other functions considered for this protein fail to achieve. We anticipate that the IDIP framework and the concepts analysed herein will aid the investigation of other proteins whose primary functional assignments have thus far been intractable.

Dasatinib-SIK2 Binding Elucidated by Homology Modeling, Molecular Docking, and Dynamics Simulations

Salt-inducible kinases (SIKs) are calcium/calmodulin-dependent protein kinase (CAMK)-like (CAMKL) family members implicated in insulin signal transduction, metabolic regulation, inflammatory response, and other processes. Here, we focused on SIK2, which is a target of the Food and Drug Administration (FDA)-approved pan inhibitor N-(2-chloro-6-methylphenyl)-2-(6-(4-(2-hydroxyethyl)piperazin-1-yl)-2-methylpyrimidin-4-ylamino)thiazole-5-carboxamide (dasatinib), and constructed four representative SIK2 structures by homology modeling. We investigated the interactions between dasatinib and SIK2 via molecular docking, molecular dynamics simulation, and binding free energy calculation and found that dasatinib showed strong binding affinity for SIK2. Binding free energy calculations suggested that the modification of various dasatinib regions may provide useful information for drug design and to guide the discovery of novel dasatinib-based SIK2 inhibitors.

Defective internal allosteric network imparts dysfunctional ATP/substrate-binding cooperativity in oncogenic chimera of protein kinase A

An aberrant fusion of the DNAJB1 and PRKACA genes generates a chimeric protein kinase (PKA-CDNAJB1) in which the J-domain of the heat shock protein 40 is fused to the catalytic α subunit of cAMP-dependent protein kinase A (PKA-C). Deceivingly, this chimeric construct appears to be fully functional, as it phosphorylates canonical substrates, forms holoenzymes, responds to cAMP activation, and recognizes the endogenous inhibitor PKI. Nonetheless, PKA-CDNAJB1 has been recognized as the primary driver of fibrolamellar hepatocellular carcinoma and is implicated in other neoplasms for which the molecular mechanisms remain elusive. Here we determined the chimera's allosteric response to nucleotide and pseudo-substrate binding. We found that the fusion of the dynamic J-domain to PKA-C disrupts the internal allosteric network, causing dramatic attenuation of the nucleotide/PKI binding cooperativity. Our findings suggest that the reduced allosteric cooperativity exhibited by PKA-CDNAJB1 alters specific recognitions and interactions between substrates and regulatory partners contributing to dysregulation.

Removal of 2H-decoupling sidebands in 13CHD2 13C-CEST profiles

A unique aspect of NMR is its capacity to provide integrated insight into both the structure and intrinsic dynamics of biomolecules. Chemical exchange phenomena that often serve as probes of dynamic processes in biological macromolecules can be quantitatively investigated with chemical exchange saturation transfer (CEST) experiments. ²H-decoupling sidebands, however, always occur in the profiles of ¹³CHD₂ ¹³C-CEST experiments when using the simple CW (continuous wave) method, which may obscure the detection of minor dips of excited states. Traditionally, these sidebands are manually eliminated from the profiles before data analysis by removing experimental points in the range of ²H-decoupling field strength ±50 Hz away from the major dips of the ground state on either side of the dips. Unfortunately, this may also eliminate potential minor dips if they overlap with the decoupling sidebands. Here, we developed methods that use pseudo-continuous waves with variable RF amplitudes distributed onto ramps for ²H decoupling. The new methods were thoroughly validated on Bruker spectrometers at a range of fields (¹H frequencies of 600, 700, and 850 MHz, and 1.1 GHz). By using these methods, we successfully removed the sidebands from the NMR profiles of ¹³CHD₂ ¹³C-CEST experiments.