Search

Supplementary Methods and Materials from Equally Potent Inhibition of c-Src and Abl by Compounds that Recognize Inactive Kinase Conformations

Supplementary Table 1 from Equally Potent Inhibition of c-Src and Abl by Compounds that Recognize Inactive Kinase Conformations

Supplementary Figures 1-2 from Equally Potent Inhibition of c-Src and Abl by Compounds that Recognize Inactive Kinase Conformations

Imatinib is an inhibitor of the Abl tyrosine kinase domain that is effective in the treatment of chronic myelogenic leukemia. Although imatinib binds tightly to the Abl kinase domain, its affinity for the closely related kinase domain of c-Src is at least 2,000-fold lower. Imatinib recognition requires a specific inactive conformation of the kinase domain, in which a conserved Asp-Phe-Gly (DFG) motif is flipped with respect to the active conformation. The inability of c-Src to readily adopt this flipped DFG conformation was thought to underlie the selectivity of imatinib for Abl over c-Src. Here, we present a series of inhibitors (DSA compounds) that are based on the core scaffold of imatinib but which bind with equally high potency to c-Src and Abl. The DSA compounds bind to c-Src in the DFG-flipped conformation, as confirmed by crystal structures and kinetic analysis. The origin of the high affinity of these compounds for c-Src is suggested by the fact that they also inhibit clinically relevant Abl variants bearing mutations in a structural element, the P-loop, that normally interacts with the phosphate groups of ATP but is folded over a substructure of imatinib in Abl. Importantly, several of the DSA compounds block the growth of Ba/F3 cells harboring imatinib-resistant BCR-ABL mutants, including the Thr315Ile “gatekeeper” mutation, but do not suppress the growth of parental Ba/F3 cells. [Cancer Res 2009;69(6):2384–92

Effective inactivation of the HER2-HER3 tumor driver has remained elusive because of the challenging attributes of the pseudokinase HER3. We report a structure-function study of constitutive HER2-HER3 signaling to identify opportunities for targeting. The allosteric activation of the HER2 kinase domain (KD) by the HER3 KD is required for tumorigenic signaling and can potentially be targeted by allosteric inhibitors. ATP binding within the catalytically inactive HER3 KD provides structural rigidity that is important for signaling, but this is mimicked, not opposed, by small molecule ATP analogs, reported here in a bosutinib-bound crystal structure. Mutational disruption of ATP binding and molecular dynamics simulation of the apo KD of HER3 identify a conformational coupling of the ATP pocket with a hydrophobic AP-2 pocket, analogous to EGFR, that is critical for tumorigenic signaling and feasible for targeting. The value of these potential target sites is confirmed in tumor growth assays using gene replacement techniques.

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.

Upon ligand binding, bone morphogenetic protein (BMP) receptors form active tetrameric complexes, comprised of two type I and two type II receptors, which then transmit signals to SMAD proteins. The link between receptor tetramerization and the mechanism of kinase activation, however, has not been elucidated. Here, using hydrogen deuterium exchange mass spectrometry (HDX-MS), small angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations, combined with analysis of SMAD signaling, we show that the kinase domain of the type I receptor ALK2 and type II receptor BMPR2 form a heterodimeric complex via their C-terminal lobes. Formation of this dimer is essential for ligand-induced receptor signaling and is targeted by mutations in BMPR2 in patients with pulmonary arterial hypertension (PAH). We further show that the type I/type II kinase domain heterodimer serves as the scaffold for assembly of the active tetrameric receptor complexes to enable phosphorylation of the GS domain and activation of SMADs., Bone morphogenetic protein (BMP) receptors are single pass transmembrane serine/threonine kinases that form tetrameric complexes comprised of two type I and two type II BMP receptors. Here the authors characterize a structure of an active type I/type II kinase tetramer providing insight into molecular mechanism driving ligand-induced signaling.

Plexins are semaphorin receptors that play essential roles in mammalian neuronal axon guidance and in many other important mammalian biological processes. Plexin signaling depends on a semaphorin-induced dimerization mechanism and is modulated by small GTPases of the Rho family, of which RND1 serves as a plexin activator yet its close homolog RhoD an inhibitor. Using molecular dynamics (MD) simulations, we showed that RND1 reinforces the plexin dimerization interface, whereas RhoD destabilizes it due to their differential interaction with the cell membrane. Upon binding plexin at the Rho-GTPase-binding domain (RBD), RND1 and RhoD interact differently with the inner leaflet of the cell membrane and exert opposite effects on the dimerization interface via an allosteric network involving the RBD, RBD linkers, and a buttress segment adjacent to the dimerization interface. The differential membrane interaction is attributed to the fact that, unlike RND1, RhoD features a short C-terminal tail and a positively charged membrane interface.

Protein-protein interactions (PPIs) are ubiquitous biomolecular processes that are central to virtually all aspects of cellular function. Identifying small molecules that modulate specific disease-related PPIs is a strategy with enormous promise for drug discovery. The design of drugs to disrupt PPIs is challenging, however, because many potential drug-binding sites at PPI interfaces are “cryptic”: When unoccupied by a ligand, cryptic sites are often flat and featureless, and thus not readily recognizable in crystal structures, with the geometric and chemical characteristics of typical small-molecule binding sites only emerging upon ligand binding. The rational design of small molecules to inhibit specific PPIs would benefit from a better understanding of how such molecules bind at PPI interfaces. To this end, we have conducted unbiased, all-atom MD simulations of the binding of four small-molecule inhibitors (SP4206 and three SP4206 analogs) to interleukin 2 (IL2)—which performs its function by forming a PPI with its receptor—without incorporating any prior structural information about the ligands’ binding. In multiple binding events, a small molecule settled into a stable binding pose at the PPI interface of IL2, resulting in a protein–small-molecule binding site and pose virtually identical to that observed in an existing crystal structure of the IL2-SP4206 complex. Binding of the small molecule stabilized the IL2 binding groove, which when the small molecule was not bound emerged only transiently and incompletely. Moreover, free energy perturbation (FEP) calculations successfully distinguished between the native and non-native IL2–small-molecule binding poses found in the simulations, suggesting that binding simulations in combination with FEP may provide an effective tool for identifying cryptic binding sites and determining the binding poses of small molecules designed to disrupt PPI interfaces by binding to such sites.

The protein K-Ras functions as a molecular switch in signaling pathways regulating cell growth. In the human mitogen-activated protein kinase (MAPK) pathway, which is implicated in many cancers, multiple K-Ras proteins are thought to assemble at the cell membrane with Ras effector proteins from the Raf family. Here we propose an atomistic structural model for such an assembly. Our starting point was an asymmetric guanosine triphosphate-mediated K-Ras dimer model, which we generated using unbiased molecular dynamics simulations and verified with mutagenesis experiments. Adding further K-Ras monomers in a head-to-tail fashion led to a compact helical assembly, a model we validated using electron microscopy and cell-based experiments. This assembly stabilizes K-Ras in its active state and presents composite interfaces to facilitate Raf binding. Guided by existing experimental data, we then positioned C-Raf, the downstream kinase MEK1 and accessory proteins (Galectin-3 and 14-3-3σ) on and around the helical assembly. The resulting Ras-Raf signalosome model offers an explanation for a large body of data on MAPK signaling.

Plexins are semaphorin receptors that play essential roles in neuronal axon guidance and in many other important biological processes. Plexin signaling depends on a semaphorin-induced dimerization mechanism, and is modulated by small signaling GTPases of the Rho family, of which RND1 serves as a plexin activator yet its close homolog RhoD an inhibitor. Using molecular dynamics (MD) simulations we showed that RND1 reinforces plexin dimerization interface whereas RhoD destabilizes it due to their differential interaction with cell membrane. Upon binding plexin dimers at the Rho-GTPase binding (RBD) domains, RND1 and RhoD interact differently with the inner leaflet of cell membrane, and exert opposite effects on the dimerization interface via an allosteric network involving the RBD domain, RBD linkers, and a buttress segment adjacent to the dimerization interface. The differential membrane interaction is attributed to the fact that, unlike RND1, RhoD features a short C-terminal tail and a positively-charged membrane interface.

AZD9291 (Osimertinib) is highly effective in treating EGFR-mutated non-small-cell lung cancers (NSCLCs) with T790M-mediated drug resistance. Despite the remarkable success of AZD9291, its binding pose with EGFR T790M remains unclear. Here, we report unbiased, atomic-level molecular dynamics (MD) simulations in which spontaneous association of AZD9291 with EGFR kinases having WT and T790M mutant gatekeepers was observed. Simulation-generated structural models suggest that the binding pose of AZD9291 with T790M differs from its binding pose with the WT, and that AZD9291 interacts extensively with the gatekeeper residue (Met 790) in T790M but not with Thr 790 in the WT, which explains why AZD9291 binds T790M with higher affinity. The MD simulation-generated models were confirmed by experimentally determined EGFR/T790M complex crystal structures. This work may facilitate the rational design of drugs that can overcome resistance mutations to AZD9291, and more generally it suggests the potential of using unbiased MD simulation to elucidate small-molecule binding poses.

The protein K-Ras functions as a molecular switch in signaling pathways regulating cell growth. In the MAPK pathway, which is implicated in many cancers, multiple K-Ras proteins are thought to assemble at the cell membrane with Ras-effector proteins from the Raf family. Here we propose an atomistic structural model for such an assembly. Our starting point was an asymmetric, GTP-mediated K-Ras dimer model, which we generated using unbiased molecular dynamics simulations and verified with mutagenesis experiments. Adding further K-Ras monomers in a head-to-tail fashion led to a compact helical assembly, a model we validated using electron microscopy and cell-based experiments. This assembly stabilizes K-Ras in its active state and presents composite interfaces to facilitate Raf binding. Guided by existing experimental data, we then positioned C-Raf, the downstream kinase MEK1, and accessory proteins (Galectin-3 and 14-3-3σ) on the helical assembly. The resulting Ras-Raf signalosome model offers an explanation for a large body of data on MAPK signaling.

The protein JAK2 is a prototypical member of the Janus kinase family, and mediates signals from numerous cytokine receptors. The constitutively active V617F mutant of JAK2 is prevalent in many bone marrow disorders, blood cancers, and autoimmune diseases, and is an important drug target. Structures have been determined for each of the four individual domains making up JAK2, and for certain pairs of these domains, but no structure of full-length JAK2 is available, and thus the mechanisms underlying JAK2 regulation and the aberrant activity of the V617F mutant have been incompletely understood. Here we propose structural models of full-length JAK2 in both its active and inactive forms. Construction of these models was informed by long-timescale molecular dynamics simulations. Subsequent mutagenesis experiments showed that mutations at the putative interdomain interfaces modulated JAK2 activity. The models provide a structural basis for understanding JAK2 autoinhibition and activation, and suggest that the constitutive activity of the V617F mutant may arise from a dual effect of destabilizing the inactive conformation and stabilizing the active conformation.

Trans-autophosphorylation is among the most prevalent means of protein kinase activation, yet its molecular basis is poorly defined. In Toll-like receptor and interleukin-1 receptor signaling pathways, the kinase IRAK4 is recruited to the membrane-proximal adaptor MyD88 through death domain (DD) interactions, forming the oligomeric Myddosome and mediating NF-κB activation. Here we show that unphosphorylated IRAK4 dimerizes in solution with a KD of 2.5 μM and that Myddosome assembly greatly enhances IRAK4 kinase domain (KD) autophosphorylation at sub-KD concentrations. The crystal structure of the unphosphorylated IRAK4(KD) dimer captures a conformation that appears to represent the actual trans-autophosphorylation reaction, with the activation loop phosphosite of one IRAK4 monomer precisely positioned for phosphotransfer by its partner. We show that dimerization is crucial for IRAK4 autophosphorylation in vitro and ligand-dependent signaling in cells. These studies identify a mechanism for oligomerization-driven allosteric autoactivation of IRAK4 that may be general to other kinases activated by autophosphorylation.

SummaryHow the epidermal growth factor receptor (EGFR) activates is incompletely understood. The intracellular portion of the receptor is intrinsically active in solution, and to study its regulation, we measured autophosphorylation as a function of EGFR surface density in cells. Without EGF, intact EGFR escapes inhibition only at high surface densities. Although the transmembrane helix and the intracellular module together suffice for constitutive activity even at low densities, the intracellular module is inactivated when tethered on its own to the plasma membrane, and fluorescence cross-correlation shows that it fails to dimerize. NMR and functional data indicate that activation requires an N-terminal interaction between the transmembrane helices, which promotes an antiparallel interaction between juxtamembrane segments and release of inhibition by the membrane. We conclude that EGF binding removes steric constraints in the extracellular module, promoting activation through N-terminal association of the transmembrane helices.

The protein tyrosine kinase Jak2 mediates signaling through numerous cytokine receptors. Jak2 possesses a pseudokinase domain (JH2) and a tyrosine kinase domain (JH1). Through unknown mechanisms, JH2 regulates the catalytic activity of JH1, and hyperactivating mutations in the JH2 region of human Jak2 are causative for myeloproliferative neoplasms (MPNs). We showed previously that Jak2 JH2 is in fact catalytically active. Here, we present crystal structures of human Jak2 JH2, both wild-type and the most prevalent MPN mutant, V617F. The structures reveal that JH2 adopts the fold of a prototypical protein kinase but binds Mg-ATP non-canonically. The structural and biochemical data indicate that the V617F mutation rigidifies α-helix C in the N lobe of JH2, which facilitates trans-phosphorylation of JH1. The crystal structures of JH2 afford new opportunities for the design of novel Jak2 therapeutics targeting MPNs.

Following Folding Fast Many protein functions involve conformational changes that occur on time-scales between tens of microseconds and milliseconds. This has limited the usefulness of all-atom molecular dynamics simulations, which are performed over shorter time-scales. Shaw et al. (p. 341 ) now report millisecond-scale, all-atom molecular dynamics simulations in an explicitly represented solvent environment. Simulation of the folding of a WW domain showed a well-defined folding pathway and simulation of the dynamics of bovine pancreatic trypsin inhibitor showed interconversion between distinct conformational states.

Events of scientific interest in molecular dynamics (MD) simulations, including conformational changes, folding transitions, and translocations of ligands and reaction products, often correspond to high-level structural rearrangements that alter contacts between molecules or among different parts of a molecule. Due to advances in computer architecture and software, MD trajectories representing such structure-changing events have become easier to generate, but the length of these trajectories poses a challenge to scientific interpretation and analysis. In this paper, we present automated methods for the detection of potentially important structure-changing events in long MD trajectories. In contrast with traditional tools for the analysis of such trajectories, our methods provide a detailed report of broken and formed contacts that aids in the identification of specific time-dependent side-chain interactions. Our approach employs a coarse-grained representation of amino acid side chains, a contact metric based on higher order generalizations of Delaunay tetrahedralization, techniques for detecting significant shifts in the resulting contact time series, and a new kernel-based measure of contact alteration activity. The analysis methods we describe are incorporated in a newly developed package, called TimeScapes, which is freely available and compatible with trajectories generated by a variety of popular MD programs. Tests based on actual microsecond time scale simulations demonstrate that the package can be used to efficiently detect and characterize important conformational changes in realistic protein systems.

Imatinib is an inhibitor of the Abl tyrosine kinase domain that is effective in the treatment of chronic myelogenic leukemia. Although imatinib binds tightly to the Abl kinase domain, its affinity for the closely related kinase domain of c-Src is at least 2,000-fold lower. Imatinib recognition requires a specific inactive conformation of the kinase domain, in which a conserved Asp-Phe-Gly (DFG) motif is flipped with respect to the active conformation. The inability of c-Src to readily adopt this flipped DFG conformation was thought to underlie the selectivity of imatinib for Abl over c-Src. Here, we present a series of inhibitors (DSA compounds) that are based on the core scaffold of imatinib but which bind with equally high potency to c-Src and Abl. The DSA compounds bind to c-Src in the DFG-flipped conformation, as confirmed by crystal structures and kinetic analysis. The origin of the high affinity of these compounds for c-Src is suggested by the fact that they also inhibit clinically relevant Abl variants bearing mutations in a structural element, the P-loop, that normally interacts with the phosphate groups of ATP but is folded over a substructure of imatinib in Abl. Importantly, several of the DSA compounds block the growth of Ba/F3 cells harboring imatinib-resistant BCR-ABL mutants, including the Thr315Ile “gatekeeper” mutation, but do not suppress the growth of parental Ba/F3 cells. [Cancer Res 2009;69(6):2384–92

Structurally elucidating transition pathways between protein conformations gives deep mechanistic insight into protein behavior but is typically difficult. Unbiased molecular dynamics (MD) simulations provide one solution, but their computational expense is often prohibitive, motivating the development of enhanced sampling methods that accelerate conformational changes in a given direction, embodied in a collective variable. The accuracy of such methods is unclear for complex protein transitions, because obtaining unbiased MD data for comparison is difficult. Here, we use long-time scale, unbiased MD simulations of epidermal growth factor receptor kinase deactivation as a complex biological test case for two widely used methods-steered molecular dynamics (SMD) and the string method. We found that common collective variable choices, based on the root-mean-square deviation (RMSD) of the entire protein, prevented the methods from producing accurate paths, even in SMD simulations on the time scale of the unbiased transition. Using collective variables based on the RMSD of the region of the protein known to be important for the conformational change, however, enabled both methods to provide a more accurate description of the pathway in a fraction of the simulation time required to observe the unbiased transition.

Calcium-dependent protein kinases (CDPKs) comprise the major group of Ca2+-regulated kinases in plants and protists. It has long been assumed that CDPKs are activated, like other Ca2+-regulated kinases, by derepression of the kinase domain (KD). However, we found that removal of the autoinhibitory domain from Toxoplasma gondii CDPK1 is not sufficient for kinase activation. From a library of heavy chain-only antibody fragments (VHHs), we isolated an antibody (1B7) that binds TgCDPK1 in a conformation-dependent manner and potently inhibits it. We uncovered the molecular basis for this inhibition by solving the crystal structure of the complex and simulating, through molecular dynamics, the effects of 1B7-kinase interactions. In contrast to other Ca2+-regulated kinases, the regulatory domain of TgCDPK1 plays a dual role, inhibiting or activating the kinase in response to changes in Ca2+ concentrations. We propose that the regulatory domain of TgCDPK1 acts as a molecular splint to stabilize the otherwise inactive KD. This dependence on allosteric stabilization reveals a novel susceptibility in this important class of parasite enzymes.

Protein tyrosine kinases are attractive drug targets because many human diseases are associated with the deregulation of kinase activity. However, how the catalytic kinase domain integrates different signals and switches from an active to an inactive conformation remains incompletely understood. Here we identify an allosteric network of dynamically coupled amino acids in Src kinase that connects regulatory sites to the ATP- and substrate-binding sites. Surprisingly, reactants (ATP and peptide substrates) bind with negative cooperativity to Src kinase while products (ADP and phosphopeptide) bind with positive cooperativity. We confirm the molecular details of the signal relay through the allosteric network by biochemical studies. Experiments on two additional protein tyrosine kinases indicate that the allosteric network may be largely conserved among these enzymes. Our work provides new insights into the regulation of protein tyrosine kinases and establishes a potential conduit by which resistance mutations to ATP-competitive kinase inhibitors can affect their activity., Protein tyrosine kinases are subject to multiple regulatory mechanisms. Foda et al. show that reactants and products of the tyrosine kinase Src bind its catalytic domain with opposite cooperativity, and identify an allosteric network of dynamically coupled amino acids that underlie this behaviour.

Although the thermodynamic principles that control the binding of drug molecules to their protein targets are well understood, detailed experimental characterization of the process by which such binding occurs has proven challenging. We conducted relatively long, unguided molecular dynamics simulations in which a ligand (the cancer drug dasatinib or the kinase inhibitor PP1) was initially placed at a random location within a box that also contained a protein (Src kinase) to which that ligand was known to bind. In several of these simulations, the ligand correctly identified its target binding site, forming a complex virtually identical to the crystallographically determined bound structure. The simulated trajectories provide a continuous, atomic-level view of the entire binding process, revealing persistent and noteworthy intermediate conformations and shedding light on the role of water molecules. The technique we employed, which does not assume any prior knowledge of the binding site’s location, may prove particularly useful in the development of allosteric inhibitors that target previously undiscovered binding sites.

The type I insulin-like growth factor receptor (IGF1R) is involved in growth and survival of normal and neoplastic cells. A ligand-dependent conformational change is thought to regulate IGF1R activity, but the nature of this change is unclear. We point out an underappreciated dimer in the crystal structure of the related Insulin Receptor (IR) with Insulin bound that allows direct comparison with unliganded IR and suggests a mechanism by which ligand regulates IR/IGF1R activity. We test this mechanism in a series of biochemical and biophysical assays and find the IGF1R ectodomain maintains an autoinhibited state in which the TMs are held apart. Ligand binding releases this constraint, allowing TM association and unleashing an intrinsic propensity of the intracellular regions to autophosphorylate. Enzymatic studies of full-length and kinase-containing fragments show phosphorylated IGF1R is fully active independent of ligand and the extracellular-TM regions. The key step triggered by ligand binding is thus autophosphorylation. DOI: http://dx.doi.org/10.7554/eLife.03772.001, eLife digest Hormones are chemicals that are produced to carry signals around the body. Mammals, including humans, need hormones called insulin and insulin-like growth factors (or IGF for short) to grow and develop normally. These hormones bind to, and activate, specific proteins—known as receptors—that span from the outside of the cell to the inside through the cell's surface membrane. The insulin and IGF receptors are complex molecules, each composed of two identical protein subunits that are linked together. The hormones bind to the part of the proteins (known as the extracellular domains) that are on the outside surface of cells. When no hormone is bound to the receptor, these two extracellular domains form an inverted ‘V’ shape and the two regions that cross the membrane are held far apart. It was unclear, however, how this shape changes when the hormone binds, and how this change in shape activates the receptor. Kavran et al. have now compared the known three-dimensional structures of the extracellular domains of the insulin receptor, both with and without a molecule of insulin bound to it. This comparison highlighted an interaction between the two receptor subunits that was disrupted when insulin was bound to the receptor. This interaction appeared to stabilize the inverted ‘V’ structure and hold apart the parts of the proteins that span the surface—which in turn separates the regions of the proteins that are inside the cell. Kavran et al. suggest that this separation keeps the insulin receptor in an inactive state. Removing either the whole extracellular domain of the IGF receptor—or a smaller portion that keeps the regions that cross the cell membrane separate—resulted in the receptor being activated, even when insulin-like growth factor was not bound to the receptor. Kavran et al. then confirmed that when a molecule of insulin-like growth factor binds to the extracellular domain of the IGF receptor, it causes a shape change that moves the parts of the receptor that span the cell membrane closer together. This results in the regions of the receptor proteins located inside the cell adding chemical tags, called phosphate groups, to one another—which activates the receptor. When there are problems with how these receptors are activated or inactivated in humans, serious disorders including diabetes and cancer can occur. As such, the findings of Kavran et al. might help future work aimed at regulating the activation of the insulin and IGF receptors to treat these and other diseases. DOI: http://dx.doi.org/10.7554/eLife.03772.002

Janus kinase-2 (JAK2) mediates signaling by various cytokines, including erythropoietin and growth hormone. JAK2 possesses tandem pseudokinase and tyrosine kinase domains. Mutations in the pseudokinase domain are causally linked to myeloproliferative neoplasms (MPNs) in humans. The structure of the JAK2 tandem kinase domains is unknown, and therefore the molecular bases for pseudokinase-mediated autoinhibition and pathogenic activation remain obscure. Using unbiased molecular dynamics simulations of protein-protein docking, we produced a structural model for the autoinhibitory interaction between the JAK2 pseudokinase and kinase domains. A striking feature of our model, which is supported by mutagenesis experiments, is that nearly all of the disease mutations map to the domain interface. The simulations indicate that the kinase domain is stabilized in an inactive state by the pseudokinase domain, and they offer a molecular rationale for the hyperactivity of V617F, the predominant JAK2 MPN mutation.

Protein-protein interaction sites are predicted from a neural network with sequence profiles of neighboring residues and solvent exposure as input. The network was trained on 615 pairs of nonhomologous complex-forming proteins. Tested on a different set of 129 pairs of nonhomologous complex-forming proteins, 70% of the 11,004 predicted interface residues are actually located in the interfaces. These 7732 correctly predicted residues account for 65% of the 11,805 residues making up the 129 interfaces. The main strength of the network predictor lies in the fact that neighbor lists and solvent exposure are relatively insensitive to structural changes accompanying complex formation. As such, it performs equally well with bound or unbound structures of the proteins. For a set of 35 test proteins, when the input was calculated from the bound and unbound structures, the correct fractions of the predicted interface residues were 69 and 70%, respectively.

The concept of potential of mean force (PMF) is now widely used in predicting protein structures. Proteins notably differ from liquids by their inhomogeneity and chain connectivity. Does meaningful correspondence exist between PMFs in proteins and PMFs in liquids? This question was addressed in this article. We constructed “proteins” each with 90 residues selected from a system of 500 hard spheres. The residues were of two types, N and P. They interact among themselves (with energies ENN, EPP, ENP) and the 410 “solvent” spheres (with energies ENS and EPS). Out of the 500 hard spheres, we first identified all chains consisting of 90 residues that have appropriate distances between nearest neighbors. The conformation of a protein was selected as the one having the lowest total energy among the 3.7 million chains. A corresponding liquid system was constructed without imposing distance constraints among solute spheres. The PMFs obtained from the proteins and the liquid system show remarkable similarities. For...