Christopher L. McClendon

Turning a protein kinase on or off from a single allosteric site via disulfide trapping.

There is significant interest in identifying and characterizing allosteric sites in enzymes such as protein kinases both for understanding allosteric mechanisms as well as for drug discovery. Here, we apply a site-directed technology, disulfide trapping, to interrogate structurally and functionally how an allosteric site on the Ser/Thr kinase, 3-phosphoinositide-dependent kinase 1 (PDK1)—the PDK1-interacting-fragment (PIF) pocket—is engaged by an activating peptide motif on downstream substrate kinases (PIFtides) and by small molecule fragments. By monitoring pairwise disulfide conjugation between PIFtide and PDK1 cysteine mutants, we defined the PIFtide binding orientation in the PIF pocket of PDK1 and assessed subtle relationships between PIFtide positioning and kinase activation. We also discovered a variety of small molecule fragment disulfides (< 300 Da) that could either activate or inhibit PDK1 by conjugation to the PIF pocket, thus displaying greater functional diversity than is displayed by PIFtides conjugated to the same sites. Biochemical data and three crystal structures provided insight into the mechanism of action of the best fragment activators and inhibitors. These studies show that disulfide trapping is useful for characterizing allosteric sites on kinases and that a single allosteric site on a protein kinase can be exploited for both activation and inhibition by small molecules.

Dynamic architecture of a protein kinase.

Significance Protein kinases represent a critically important family of regulatory enzymes. Their activity can be altered by mutations and binding events distant from the active site. To understand the nature of these long-distance effects, we used microsecond-timescale molecular dynamic simulation to subdivide a prototypical kinase, protein kinase A, into contiguous communities that exhibit internally correlated motions. Surprisingly, most of these unconventional structural entities were centered around known protein kinase functions. We thus propose a new framework for analysis of protein kinase structure and function that differs from traditional representations based simply on sequence motifs and secondary structure elements. These results extend our view on the dynamic nature of protein kinases and open a door to understanding of allosteric signaling in these enzymes. Protein kinases are dynamically regulated signaling proteins that act as switches in the cell by phosphorylating target proteins. To establish a framework for analyzing linkages between structure, function, dynamics, and allostery in protein kinases, we carried out multiple microsecond-scale molecular-dynamics simulations of protein kinase A (PKA), an exemplar active kinase. We identified residue–residue correlated motions based on the concept of mutual information and used the Girvan–Newman method to partition PKA into structurally contiguous “communities.” Most of these communities included 40–60 residues and were associated with a particular protein kinase function or a regulatory mechanism, and well-known motifs based on sequence and secondary structure were often split into different communities. The observed community maps were sensitive to the presence of different ligands and provide a new framework for interpreting long-distance allosteric coupling. Communication between different communities was also in agreement with the previously defined architecture of the protein kinase core based on the “hydrophobic spine” network. This finding gives us confidence in suggesting that community analyses can be used for other protein kinases and will provide an efficient tool for structural biologists. The communities also allow us to think about allosteric consequences of mutations that are linked to disease.

Modeling Conformational Ensembles of Slow Functional Motions in Pin1-WW

Protein-protein interactions are often mediated by flexible loops that experience conformational dynamics on the microsecond to millisecond time scales. NMR relaxation studies can map these dynamics. However, defining the network of inter-converting conformers that underlie the relaxation data remains generally challenging. Here, we combine NMR relaxation experiments with simulation to visualize networks of inter-converting conformers. We demonstrate our approach with the apo Pin1-WW domain, for which NMR has revealed conformational dynamics of a flexible loop in the millisecond range. We sample and cluster the free energy landscape using Markov State Models (MSM) with major and minor exchange states with high correlation with the NMR relaxation data and low NOE violations. These MSM are hierarchical ensembles of slowly interconverting, metastable macrostates and rapidly interconverting microstates. We found a low population state that consists primarily of holo-like conformations and is a “hub” visited by most pathways between macrostates. These results suggest that conformational equilibria between holo-like and alternative conformers pre-exist in the intrinsic dynamics of apo Pin1-WW. Analysis using MutInf, a mutual information method for quantifying correlated motions, reveals that WW dynamics not only play a role in substrate recognition, but also may help couple the substrate binding site on the WW domain to the one on the catalytic domain. Our work represents an important step towards building networks of inter-converting conformational states and is generally applicable.

Synchronous Opening and Closing Motions Are Essential for cAMP-Dependent Protein Kinase A Signaling

Conformational fluctuations play a central role in enzymatic catalysis. However, it is not clear how the rates and the coordination of the motions affect the different catalytic steps. Here, we used NMR spectroscopy to analyze the conformational fluctuations of the catalytic subunit of the cAMP-dependent protein kinase (PKA-C), a ubiquitous enzyme involved in a myriad of cell signaling events. We found that the wild-type enzyme undergoes synchronous motions involving several structural elements located in the small lobe of the kinase, which is responsible for nucleotide binding and release. In contrast, a mutation (Y204A) located far from the active site desynchronizes the opening and closing of the active cleft without changing the enzymes structure, rendering it catalytically inefficient. Since the opening and closing motions govern the rate-determining product release, we conclude that optimal and coherent conformational fluctuations are necessary for efficient turnover of protein kinases.

A New Coarse-Grained Model for E. coli Cytoplasm: Accurate Calculation of the Diffusion Coefficient of Proteins and Observation of Anomalous Diffusion

A new coarse-grained model of the E. coli cytoplasm is developed by describing the proteins of the cytoplasm as flexible units consisting of one or more spheres that follow Brownian dynamics (BD), with hydrodynamic interactions (HI) accounted for by a mean-field approach. Extensive BD simulations were performed to calculate the diffusion coefficients of three different proteins in the cellular environment. The results are in close agreement with experimental or previously simulated values, where available. Control simulations without HI showed that use of HI is essential to obtain accurate diffusion coefficients. Anomalous diffusion inside the crowded cellular medium was investigated with Fractional Brownian motion analysis, and found to be present in this model. By running a series of control simulations in which various forces were removed systematically, it was found that repulsive interactions (volume exclusion) are the main cause for anomalous diffusion, with a secondary contribution from HI.

Conformational equilibrium of N-myristoylated cAMP-dependent protein kinase A by molecular dynamics simulations.

The catalytic subunit of protein kinase A (PKA-C) is subject to several post- or cotranslational modifications that regulate its activity both spatially and temporally. Among those, N-myristoylation increases the kinase affinity for membranes and might also be implicated in substrate recognition and allosteric regulation. Here, we investigated the effects of N-myristoylation on the structure, dynamics, and conformational equilibrium of PKA-C using atomistic molecular dynamics simulations. We found that the myristoyl group inserts into the hydrophobic pocket and leads to a tighter packing of the A-helix against the core of the enzyme. As a result, the conformational dynamics of the A-helix are reduced and its motions are more coupled with the active site. Our simulations suggest that cation-π interactions among W30, R190, and R93 are responsible for coupling these motions. Two major conformations of the myristoylated N-terminus are the most populated: a long loop (LL conformation), similar to Protein Data Bank (PDB) entry 1CMK , and a helix-turn-helix structure (HTH conformation), similar to PDB entry 4DFX , which shows stronger coupling between the conformational dynamics observed at the A-helix and active site. The HTH conformation is stabilized by S10 phosphorylation of the kinase via ionic interactions between the protonated amine of K7 and the phosphate group on S10, further enhancing the dynamic coupling to the active site. These results support a role of N-myristoylation in the allosteric regulation of PKA-C.

Substrate and Inhibitor-induced Dimerization and Cooperativity in Caspase-1 but Not Caspase-3

Background: The inflammatory caspase-1 shows positive cooperativity not seen for the apoptotic caspase-3. Results: Substrate binding increases the dimerization affinity and activity of caspase-1 but not for caspase-3. Conclusion: Caspase-1 is regulated by concentration with substrates, whereas caspase-3 is not. Significance: Subcellular co-localization of caspase-1 with substrates in inflammosomes may explain its more restricted family of substrates observed. Caspases are intracellular cysteine-class proteases with aspartate specificity that is critical for driving processes as diverse as the innate immune response and apoptosis, exemplified by caspase-1 and caspase-3, respectively. Interestingly, caspase-1 cleaves far fewer cellular substrates than caspase-3 and also shows strong positive cooperativity between the two active sites of the homodimer, unlike caspase-3. Biophysical and kinetic studies here present a molecular basis for this difference. Analytical ultracentrifugation experiments show that mature caspase-1 exists predominantly as a monomer under physiological concentrations that undergoes dimerization in the presence of substrate; specifically, substrate binding shifts the KD for dimerization by 20-fold. We have created a hemi-active site-labeled dimer of caspase-1, where one site is blocked with the covalent active site inhibitor, benzyloxycarbonyl-Val-Ala-Asp-fluoromethylketone. This hemi-labeled enzyme is about 9-fold more active than the apo-dimer of caspase-1. These studies suggest that substrate not only drives dimerization but also, once bound to one site in the dimer, promotes an active conformation in the other monomer. Steady-state kinetic analysis and modeling independently support this model, where binding of one substrate molecule not only increases substrate binding in preformed dimers but also drives the formation of heterodimers. Thus, the cooperativity in caspase-1 is driven both by substrate-induced dimerization as well as substrate-induced activation. Substrate-induced dimerization and activation seen in caspase-1 and not in caspase-3 may reflect their biological roles. Whereas caspase-1 cleaves a dramatically smaller number of cellular substrates that need to be concentrated near inflammasomes, caspase-3 is a constitutively active dimer that cleaves many more substrates located diffusely throughout the cell.

A dynamic hydrophobic core orchestrates allostery in protein kinases

The synchronized motions of inner core residues allosterically modulate the activity of the protein kinase A catalytic subunit. Eukaryotic protein kinases (EPKs) constitute a class of allosteric switches that mediate a myriad of signaling events. It has been postulated that EPKs’ active and inactive states depend on the structural architecture of their hydrophobic cores, organized around two highly conserved spines: C-spine and R-spine. How the spines orchestrate the transition of the enzyme between catalytically uncommitted and committed states remains elusive. Using relaxation dispersion nuclear magnetic resonance spectroscopy, we found that the hydrophobic core of the catalytic subunit of protein kinase A, a prototypical and ubiquitous EPK, moves synchronously to poise the C subunit for catalysis in response to binding adenosine 5′-triphosphate. In addition to completing the C-spine, the adenine ring fuses the β structures of the N-lobe and the C-lobe. Additional residues that bridge the two spines (I150 and V104) are revealed as part of the correlated hydrophobic network; their importance was validated by mutagenesis, which led to inactivation. Because the hydrophobic architecture of the catalytic core is conserved throughout the EPK superfamily, the present study suggests a universal mechanism for dynamically driven allosteric activation of kinases mediated by coordinated signal transmission through ordered motifs in their hydrophobic cores.

Hydrogen Bond Strengths in Phosphorylated and Sulfated Amino Acid Residues

Post-translational modification by the addition of an oxoanion functional group, usually a phosphate group and less commonly a sulfate group, leads to diverse structural and functional consequences in protein systems. Building upon previous studies of the phosphoserine residue (pSer), we address the distinct nature of hydrogen bonding interactions in phosphotyrosine (pTyr) and sulfotyrosine (sTyr) residues. We derive partial charges for these modified residues and then study them in the context of molecular dynamics simulation of model tripeptides and sulfated protein complexes, potentials of mean force for interacting residue pairs, and a survey of the interactions of modified residues among experimental protein structures. Overall, our findings show that for pTyr, bidentate interactions with Arg are particularly dominant, as has been previously demonstrated for pSer. sTyr interactions with Arg are significantly weaker, even as compared to the same interactions made by the Glu residue. Our work sheds light on the distinct nature of these modified tyrosine residues, and provides a physical-chemical foundation for future studies with the goal of understanding their roles in systems of biological interest.

The role of tyrosine sulfation in the dimerization of the CXCR4:SDF-1 complex

Oligomerization of G protein‐coupled receptors is a recognized mode of regulation of receptor activities, with alternate oligomeric states resulting in different signaling functions. The CXCR4 chemokine receptor is a G protein‐coupled receptor that is post‐translationally modified by tyrosine sulfation at three sites on its N‐terminus (Y7, Y12, Y21), leading to enhanced affinity for its ligand, stromal cell derived factor (SDF‐1, also called CXCL12). The complex has been implicated in cancer metastasis and is a therapeutic target in cancer treatment. Using molecular dynamics simulation of NMR‐derived structures of the CXCR4 N‐terminus in complex with SDF‐1, and calculations of electrostatic binding energies for these complexes, we address the role of tyrosine sulfation in this complex. Our results show that sulfation stabilizes the dimeric state of the CXCR4:SDF‐1 complex through hydrogen bonding across the dimer interface, conformational changes in residues at the dimer interface, and an enhancement in electrostatic binding energies associated with dimerization. These findings suggest a mechanism through which post‐translational modifications such as tyrosine sulfation might regulate downstream function through modulation of the oligomeric state of the modified system.