Interactions of protein kinase C-α C1A and C1B domains with membranes: a combined computational and experimental study

Binding of the plant-derived toxin simplexin to bovine protein kinase C: insights from molecular dynamics

Pimelea poisoning of cattle is toxicologically linked to the activation of bovine protein kinase C (PKC) by the plant-derived toxin simplexin. To understand the affinity of PKC for simplexin, we performed molecular dynamics (MD) studies of simplexin, simplexin analogues, and several other activators of PKC. Binding enthalpy calculations indicated that simplexin had the strongest affinity for PKCα-C1B among the activators studied. Key to simplexin's affinity is its ability to form more hydrogen bonds to PKC, compared to the other activators. The C-3 carbonyl group and C-20 hydroxyl group of simplexin were identified as especially important for stabilizing the PKC binding interaction. The hydrophobic alkyl chain of simplexin induces deep membrane embedding of the PKC-simplexin complex, enhancing the protein-ligand hydrogen bonding. Our findings align with previous experiments on structure-activity relationships (SAR) for simplexin analogues, and provide insights that may guide the development of interventions or treatments for Pimelea poisoning.

PDK1:PKCα heterodimer association-dissociation dynamics in single-molecule diffusion tracks on a target membrane

Previous studies have documented the formation of a heterodimer between the two protein kinases PDK1 and PKCα on a lipid bilayer containing their target lipids. This work investigates the association-dissociation kinetics of this PDK1:PKCα heterodimer. The approach monitors the two-dimensional diffusion of single, membrane-associated PDK1 molecules for diffusivity changes as PKCα molecules bind and unbind. In the absence of PKCα, a membrane-associated PDK1 molecule exhibits high diffusivity (or large diffusion constant, D) because its membrane-contacting PH domain binds the target PIP3 lipid headgroup with little bilayer penetration, yielding minimal frictional drag against the bilayer. In contrast, membrane-associated PKCα contacts the bilayer via its C1A, C1B, and C2 domains, which each bind at least one target lipid with significant bilayer insertion, yielding a large frictional drag and low diffusivity. The present findings reveal that individual fluor-PDK1 molecules freely diffusing on the membrane surface undergo reversible switching between distinct high and low diffusivity states, corresponding to the PDK1 monomer and the PDK1:PKCα heterodimer, respectively. The observed single-molecule diffusion trajectories are converted to step length time courses, then subjected to two-state, hidden Markov modeling and dwell time analysis. The findings reveal that both the PDK1 monomer state and the PDK1:PKCα heterodimer state decay via simple exponential kinetics, yielding estimates of rate constants for state switching in both directions. Notably, the PDK1:PKCα heterodimer has been shown to competitively inhibit PDK1 phosphoactivation of AKT1, and is believed to play a tumor suppressor role by limiting excess activation of the highly oncogenic PDK1/AKT1/mTOR pathway. Thus, the present elucidation of the PDK1:PKCα association-dissociation kinetics has important biological and medical implications. More broadly, the findings illustrate the power of single-molecule diffusion measurements to reveal the kinetics of association-dissociation events in membrane signaling reactions that yield a large change in diffusive mobility.

Insertion Depth Modulates Protein Kinase C-δ-C1b Domain Interactions with Membrane Cholesterol as Revealed by MD Simulations

Protein kinase C delta (PKC-δ) is an important signaling molecule in human cells that has both proapoptotic as well as antiapoptotic functions. These conflicting activities can be modulated by two classes of ligands, phorbol esters and bryostatins. Phorbol esters are known tumor promoters, while bryostatins have anti-cancer properties. This is despite both ligands binding to the C1b domain of PKC-δ (δC1b) with a similar affinity. The molecular mechanism behind this discrepancy in cellular effects remains unknown. Here, we have used molecular dynamics simulations to investigate the structure and intermolecular interactions of these ligands bound to δC1b with heterogeneous membranes. We observed clear interactions between the δC1b-phorbol complex and membrane cholesterol, primarily through the backbone amide of L250 and through the K256 side-chain amine. In contrast, the δC1b-bryostatin complex did not exhibit interactions with cholesterol. Topological maps of the membrane insertion depth of the δC1b-ligand complexes suggest that insertion depth can modulate δC1b interactions with cholesterol. The lack of cholesterol interactions suggests that bryostatin-bound δC1b may not readily translocate to cholesterol-rich domains within the plasma membrane, which could significantly alter the substrate specificity of PKC-δ compared to δC1b-phorbol complexes.

Insights into the behavior of unsaturated diacylglycerols in mixed lipid bilayers in relation to protein kinase C activation-A molecular dynamics simulation study

The lipid second messenger diacylglycerol (DAG) is known for its involvement in many types of cellular signaling, especially as an endogenous agonist for protein kinase C (PKC). Evidence has emerged that the degree of saturation of the DAG molecules can affect PKC activation. DAG molecules with different acyl chain saturation have not only been observed to induce varying extents of PKC activation, but also to express selectivity towards different PKC isozymes. Both qualities are important for precise therapeutic activation of PKC; understanding DAG behavior at the molecular level in different environments has much potential in the development of drugs to target PKC. We used molecular dynamics simulations to study the behavior of two different unsaturated DAG species in lipid environments with varying degrees of unsaturation. We focus on phosphatidylethanolamine (PE) instead of phosphatidylcholine (PC) to more accurately model the relevant biomembranes. The effect of cholesterol (CHOL) on these systems was also explored. We found that both high level of unsaturation in the acyl chains of the DAG species and presence of CHOL in the surrounding membrane increase DAG molecule availability at the lipid-water interface. This can partially explain the previously observed differences in PKC activation strength and specificity, the complete mechanism is, however, likely to be more complex. Our simulations coupled with the current understanding of lipids highlight the need for more simulations of biologically accurate lipid environments in order to determine the correct correlations between molecular mechanisms and biological behavior when studying PKC activation.

PKCα and PKCδ: Friends and Rivals

PKC comprises a large family of serine/threonine kinases that share a requirement for allosteric activation by lipids. While PKC isoforms have significant homology, functional divergence is evident among subfamilies and between individual PKC isoforms within a subfamily. Here, we highlight these differences by comparing the regulation and function of representative PKC isoforms from the conventional (PKCα) and novel (PKCδ) subfamilies. We discuss how unique structural features of PKCα and PKCδ underlie differences in activation and highlight the similar, divergent, and even opposing biological functions of these kinases. We also consider how PKCα and PKCδ can contribute to pathophysiological conditions and discuss challenges to targeting these kinases therapeutically.

Structural anatomy of Protein Kinase C C1 domain interactions with diacylglycerol and other agonists

Diacylglycerol (DAG) is a versatile lipid whose 1,2-sn-stereoisomer serves both as second messenger in signal transduction pathways that control vital cellular processes, and as metabolic precursor for downstream signaling lipids such as phosphatidic acid. Effector proteins translocate to available DAG pools in the membranes by using conserved homology 1 (C1) domains as DAG-sensing modules. Yet, how C1 domains recognize and capture DAG in the complex environment of a biological membrane has remained unresolved for the 40 years since the discovery of Protein Kinase C (PKC) as the first member of the DAG effector cohort. Herein, we report the high-resolution crystal structures of a C1 domain (C1B from PKCδ) complexed to DAG and to each of four potent PKC agonists that produce different biological readouts and that command intense therapeutic interest. This structural information details the mechanisms of stereospecific recognition of DAG by the C1 domains, the functional properties of the lipid-binding site, and the identities of the key residues required for the recognition and capture of DAG and exogenous agonists. Moreover, the structures of the five C1 domain complexes provide the high-resolution guides for the design of agents that modulate the activities of DAG effector proteins.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Structural insights into C1-ligand interactions: Filling the gaps by in silico methods

Protein Kinase C isoenzymes (PKCs) are the key mediators of the phosphoinositide signaling pathway, which involves regulated hydrolysis of phosphatidylinositol (4,5)-bisphosphate to diacylglycerol (DAG) and inositol-1,4,5-trisphosphate. Dysregulation of PKCs is implicated in many human diseases making this class of enzymes an important therapeutic target. Specifically, the DAG-sensing cysteine-rich conserved homology-1 (C1) domains of PKCs have emerged as promising targets for pharmaceutical modulation. Despite significant progress, the rational design of the C1 modulators remains challenging due to difficulties associated with structure determination of the C1-ligand complexes. Given the dearth of experimental structural data, computationally derived models have been instrumental in providing atomistic insight into the interactions of the C1 domains with PKC agonists. In this review, we provide an overview of the in silico approaches for seven classes of C1 modulators and outline promising future directions.

Rigorous Computational Study Reveals What Docking Overlooks: Double Trouble from Membrane Association in Protein Kinase C Modulators

Increasing protein kinase C (PKC) activity is of potential therapeutic value. Its activation involves an interaction between the C1 domain and diacylglycerol (DAG) at intracellular membrane surfaces; DAG mimetics hold promise as new drugs. We previously developed the isophthalate derivative HMI-1a3, an effective but highly lipophilic (clogP = 6.46) DAG mimetic. Although a less lipophilic pyrimidine analog, PYR-1gP (clogP = 3.30), gave positive results in computational docking, it unexpectedly presented greatly diminished binding to PKC in vitro. Through more rigorous computational molecular modeling, we reveal that, unlike HMI-1a3, PYR-1gP forms an intramolecular hydrogen bond, which both obstructs binding and reorients PYR-1gP in the membrane in a fashion that prevents it from correctly accessing the PKC C1 domain. Our results highlight the great value of molecular dynamics simulations as a key component for the drug design process of ligands targeting weakly membrane-associated proteins, where simulation in the relevant membrane environment is crucial for obtaining biologically applicable results.

Anionic Lipids Impact RAS-Binding Site Accessibility and Membrane Binding Affinity of CRAF RBD-CRD

CRAF activation requires binding to membrane-anchored and active GTP-bound RAS. Whereas its RAS-binding domain (RBD) contains the main binding interface to the RAS G domain, its cysteine-rich domain (CRD) is responsible for association to anionic lipid-rich membranes. Both RAF domains are connected by a short linker, and it remains unclear if the two domains act independently or if one domain can impact the function of the other. Here, we used a combination of coarse-grained and all-atom molecular dynamics simulations of a CRAF RBD-CRD construct to investigate the dynamics of the RBD when it is tethered to CRD that is anchored to a POPC:POPS model membrane. First, we show that the RBD positioning is very dynamic with a preferential localization near the membrane surface. Next, we show that membrane-localized RBD has its RAS-binding interface mostly inaccessible because of its proximity to the membrane. Several positively charged residues in this interface were identified from simulations as important for driving RBD association to the membrane. Surface plasmon resonance (SPR) measurements confirmed that mutations of these RBD residues reduced the liposome partitioning of RBD-CRD. Last, simulations indicated that the presence of RBD near the membrane led to a local enrichment of anionic lipids that could potentially enhance the membrane affinity of the entire RBD-CRD construct. This was supported by SPR measurements that showed stronger liposome partitioning of RBD-CRD relative to CRD alone. These findings thus suggest that the RBD and CRD have synergistic effects on their membrane dynamics, with CRD bringing RBD closer to the membrane that impacts its accessibility to RAS and with RBD causing local anionic lipid enrichment that enhances the overall affinity between the membrane and RBD-CRD. These mechanisms have potential implications on the order of events of the interactions between RAS and CRAF at the membrane.

Structural and functional characterization of Solanum lycopersicum phosphatidylinositol 3-kinase C2 domain

Phosphatidylinositol 3-kinases (PI3Ks) are characterized by the presence of a C2 domain at the N-terminal end (class I, III); or at both the N-terminal and C-terminal ends (class II), sometimes including a Plextrin homology domain and/or a Ras domain. Plant PI3Ks are analogous to the class III mammalian PI3K. An N-terminal fragment (~170 aa) of the tomato PI3K regulatory domain including the C2 domain, was cloned and expressed in a bacterial system. This protein was purified to homogeneity and its physicochemical properties analyzed. The purified protein showed strong binding with monophosphorylated phosphatidylinositols, and the binding was dependent on calcium ion concentration and pH. In the overall tertiary structure of PI3K, C2 domain showed unique characteristics, having three antiparallel beta-sheets, hydrophobic regions, acidic as well as alkaline motifs, that can enable its membrane binding upon activation. To elucidate the functional significance of C2 domain, transgenic tobacco plants expressing the C2 domain of PI3K were generated. Transgenic plants showed defective pollen development and disrupted seed set. Flowers from the PI3K-C2 transgenic plants showed delayed wilting, and a decrease in ethylene production. It is likely that introduction of the PI3K-C2 segment may have interfered with the normal binding of PI3K to the membrane, delaying the onset of membrane lipid catabolism that lead to senescence.

It Takes Two to Tango: Activation of Protein Kinase D by Dimerization

The recent discovery and structure determination of a novel ubiquitin-like dimerization domain in protein kinase D (PKD) has significant implications for its activation. PKD is a serine/threonine kinase activated by the lipid second messenger diacylglycerol (DAG). It is an essential and highly conserved protein that is implicated in plasma membrane directed trafficking processes from the trans-Golgi network. However, many open questions surround its mechanism of activation, its localization, and its role in the biogenesis of cargo transport carriers. In reviewing this field, the focus is primarily on the mechanisms that control the activation of PKD at precise locations in the cell. In light of the new structural findings, the understanding of the mechanisms underlying PKD activation is critically evaluated, with particular emphasis on the role of dimerization in PKD autophosphorylation, and the provenance and recognition of the DAG that activates PKD.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Molecular recognition of RAS/RAF complex at the membrane: Role of RAF cysteine-rich domain

Activation of RAF kinase involves the association of its RAS-binding domain (RBD) and cysteine-rich domain (CRD) with membrane-anchored RAS. However, the overall architecture of the RAS/RBD/CRD ternary complex and the orientations of its constituent domains at the membrane remain unclear. Here, we have combined all-atom and coarse-grained molecular dynamics (MD) simulations with experimental data to construct and validate a model of membrane-anchored CRD, and used this as a basis to explore models of membrane-anchored RAS/RBD/CRD complex. First, simulations of the CRD revealed that it anchors to the membrane via insertion of its two hydrophobic loops, which is consistent with our NMR measurements of CRD bound to nanodiscs. Simulations of the CRD in the context of membrane-anchored RAS/RBD then show how CRD association with either RAS or RBD could play an unexpected role in guiding the membrane orientations of RAS/RBD. This finding has implications for the formation of RAS-RAS dimers, as different membrane orientations of RAS expose distinct putative dimerization interfaces.

REDOR NMR Reveals Multiple Conformers for a Protein Kinase C Ligand in a Membrane Environment

Bryostatin 1 (henceforth bryostatin) is in clinical trials for the treatment of Alzheimer's disease and for HIV/AIDS eradication. It is also a preclinical lead for cancer immunotherapy and other therapeutic indications. Yet nothing is known about the conformation of bryostatin bound to its protein kinase C (PKC) target in a membrane microenvironment. As a result, efforts to design more efficacious, better tolerated, or more synthetically accessible ligands have been limited to structures that do not include PKC or membrane effects known to influence PKC-ligand binding. This problem extends more generally to many membrane-associated proteins in the human proteome. Here, we use rotational-echo double-resonance (REDOR) solid-state NMR to determine the conformations of PKC modulators bound to the PKCδ-C1b domain in the presence of phospholipid vesicles. The conformationally limited PKC modulator phorbol diacetate (PDAc) is used as an initial test substrate. While unanticipated partitioning of PDAc between an immobilized protein-bound state and a mobile state in the phospholipid assembly was observed, a single conformation in the bound state was identified. In striking contrast, a bryostatin analogue (bryolog) was found to exist exclusively in a protein-bound state, but adopts a distribution of conformations as defined by three independent distance measurements. The detection of multiple PKCδ-C1b-bound bryolog conformers in a functionally relevant phospholipid complex reveals the inherent dynamic nature of cellular systems that is not captured with single-conformation static structures. These results indicate that binding, selectivity, and function of PKC modulators, as well as the design of new modulators, are best addressed using a dynamic multistate model, an analysis potentially applicable to other membrane-associated proteins.

Computational and theoretical approaches for studies of a lipid recognition protein on biological membranes

Many cellular functions, including cell signaling and related events, are regulated by the association of peripheral membrane proteins (PMPs) with biological membranes containing anionic lipids, e.g., phosphatidylinositol phosphate (PIP). This association is often mediated by lipid recognition modules present in many PMPs. Here, I summarize computational and theoretical approaches to investigate the molecular details of the interactions and dynamics of a lipid recognition module, the pleckstrin homology (PH) domain, on biological membranes. Multiscale molecular dynamics simulations using combinations of atomistic and coarse-grained models yielded results comparable to those of actual experiments and could be used to elucidate the molecular mechanisms of the formation of protein/lipid complexes on membrane surfaces, which are often difficult to obtain using experimental techniques. Simulations revealed some modes of membrane localization and interactions of PH domains with membranes in addition to the canonical binding mode. In the last part of this review, I address the dynamics of PH domains on the membrane surface. Local PIP clusters formed around the proteins exhibit anomalous fluctuations. This dynamic change in protein-lipid interactions cause temporally fluctuating diffusivity of proteins, i.e., the short-term diffusivity of the bound protein changes substantially with time, and may in turn contribute to the formation/dissolution of protein complexes in membranes.

Synthesis and Evaluation of Dimeric Derivatives of Diacylglycerol-Lactones as Protein Kinase C Ligands

Protein kinase C (PKC) mediates a central cellular signal transduction pathway involved in disorders such as cancer and Alzheimer's disease. PKC is regulated by binding of the second messenger sn-1,2-diacylglycerol (DAG) to its tandem C1 domains, designated C1a and C1b, leading both to PKC activation and to its translocation to the plasma membrane and to internal organelles. Depending on the isoform, there may be differences in the ligand selectivity of the C1a and C1b domains, and there is different spacing between the C1 domains of the conventional and novel PKCs. Bivalent ligands have the potential to exploit these differences between isoforms, yielding isoform selectivity. In the present study, we describe the synthesis of a series of dimeric derivatives of conformationally constrained diacylglycerol (DAG) analogs (DAG-lactones). We characterize the derivatives in vitro for their binding affinities, both to a single C1 domain (the C1b domain of PKCδ) as well as to the conventional PKCα isoform and the novel PKCδ isoform, and we measure their abilities to cause translocation of PKCδ and PKCε in intact cells. The dimeric compound with the 10-carbon linker was modestly more effective for the isolated PKCδ C1b domain than was the monomeric compound. For the intact PKCα and PKCδ, the shortest DAG-lactone dimer had similar affinity to the monomer and affinity decreased progressively up to the 16-carbon linker. The dimeric derivatives did not cause the Golgi accumulation of PKCδ. The present results provide important insights into the development of new chemical tools for biological studies on PKC.

Toggling of Diacylglycerol Affinity Correlates with Conformational Plasticity in C1 Domains

Conserved homology-1 (C1) domains are peripheral membrane domains that target their host proteins to diacylglycerol (DAG)-containing membranes. It has been previously shown that a conservative aromatic mutation of a single residue in the C1 domain has a profound effect on DAG affinity. We report that the "DAG-toggling" mutation changes the conformational dynamics of the loop region that forms the binding site for the C1 activators. Moreover, there is a correlation among the residue identity at the mutation site, DAG affinity, and loop dynamics in four C1 variants. We propose that "toggling" of DAG affinity may occur through modulation of both protein-membrane interactions and the geometry of the activator-binding cleft, with the loop dynamics being responsible for the latter.

Molecular dynamics simulations reveal ligand-controlled positioning of a peripheral protein complex in membranes

Bryostatin is in clinical trials for Alzheimer’s disease, cancer, and HIV/AIDS eradication. It binds to protein kinase C competitively with diacylglycerol, the endogenous protein kinase C regulator, and plant-derived phorbol esters, but each ligand induces different activities. Determination of the structural origin for these differing activities by X-ray analysis has not succeeded due to difficulties in co-crystallizing protein kinase C with relevant ligands. More importantly, static, crystal-lattice bound complexes do not address the influence of the membrane on the structure and dynamics of membrane-associated proteins. To address this general problem, we performed long-timescale (400–500 µs aggregate) all-atom molecular dynamics simulations of protein kinase C–ligand–membrane complexes and observed that different protein kinase C activators differentially position the complex in the membrane due in part to their differing interactions with waters at the membrane inner leaf. These new findings enable new strategies for the design of simpler, more effective protein kinase C analogs and could also prove relevant to other peripheral protein complexes.

Natural supplies of bryostatin, a compound in clinical trials for Alzheimer’s disease, cancer, and HIV, are scarce. Here, the authors perform molecular dynamics simulations to understand how bryostatin interacts with membrane-bound protein kinase C, offering insights for the design of bryostatin analogs.

Highly Coarse-Grained Representations of Transmembrane Proteins

Numerous biomolecules and biomolecular complexes, including transmembrane proteins (TMPs), are symmetric or at least have approximate symmetries. Highly coarse-grained models of such biomolecules, aiming at capturing the essential structural and dynamical properties on resolution levels coarser than the residue scale, must preserve the underlying symmetry. However, making these models obey the correct physics is in general not straightforward, especially at the highly coarse-grained resolution where multiple (∼3-30 in the current study) amino acid residues are represented by a single coarse-grained site. In this paper, we propose a simple and fast method of coarse-graining TMPs obeying this condition. The procedure involves partitioning transmembrane domains into contiguous segments of equal length along the primary sequence. For the coarsest (lowest-resolution) mappings, it turns out to be most important to satisfy the symmetry in a coarse-grained model. As the resolution is increased to capture more detail, however, it becomes gradually more important to match modular repeats in the secondary structure (such as helix-loop repeats) instead. A set of eight TMPs of various complexity, functionality, structural topology, and internal symmetry, representing different classes of TMPs (ion channels, transporters, receptors, adhesion, and invasion proteins), has been examined. The present approach can be generalized to other systems possessing exact or approximate symmetry, allowing for reliable and fast creation of multiscale, highly coarse-grained mappings of large biomolecular assemblies.

Structural Basis for the Failure of the C1 Domain of Ras Guanine Nucleotide Releasing Protein 2 (RasGRP2) to Bind Phorbol Ester with High Affinity

The C1 domain represents the recognition module for diacylglycerol and phorbol esters in protein kinase C, Ras guanine nucleotide releasing protein (RasGRP), and related proteins. RasGRP2 is exceptional in that its C1 domain has very weak binding affinity (Kd = 2890 ± 240 nm for [(3)H]phorbol 12,13-dibutyrate. We have identified four amino acid residues responsible for this lack of sensitivity. Replacing Asn(7), Ser(8), Ala(19), and Ile(21) with the corresponding residues from RasGRP1/3 (Thr(7), Tyr(8), Gly(19), and Leu(21), respectively) conferred potent binding affinity (Kd = 1.47 ± 0.03 nm) in vitro and membrane translocation in response to phorbol 12-myristate 13-acetate in LNCaP cells. Mutant C1 domains incorporating one to three of the four residues showed intermediate behavior with S8Y making the greatest contribution. Binding activity for diacylglycerol was restored in parallel. The requirement for anionic phospholipid for [(3)H]phorbol 12,13-dibutyrate binding was determined; it decreased in going from the single S8Y mutant to the quadruple mutant. The full-length RasGRP2 protein with the mutated C1 domains also showed strong phorbol ester binding, albeit modestly weaker than that of the C1 domain alone (Kd = 8.2 ± 1.1 nm for the full-length protein containing all four mutations), and displayed translocation in response to phorbol ester. RasGRP2 is a guanyl exchange factor for Rap1. Consistent with the ability of phorbol ester to induce translocation of the full-length RasGRP2 with the mutated C1 domain, phorbol ester enhanced the ability of the mutated RasGRP2 to activate Rap1. Modeling confirmed that the four mutations helped the binding cleft maintain a stable conformation.

Protein Kinase Cα (PKCα) Is Resistant to Long Term Desensitization/Down-regulation by Prolonged Diacylglycerol Stimulation

Sustained activation of PKCα is required for long term physiological responses, such as growth arrest and differentiation. However, studies with pharmacological agonists (e.g. phorbol 12-myristate 13-acetate (PMA)) indicate that prolonged stimulation leads to PKCα desensitization via dephosphorylation and/or degradation. The current study analyzed effects of chronic stimulation with the physiological agonist diacylglycerol. Repeated addition of 1,2-dioctanoyl-sn-glycerol (DiC8) resulted in sustained plasma membrane association of PKCα in a pattern comparable with that induced by PMA. However, although PMA potently down-regulated PKCα, prolonged activation by DiC8 failed to engage known desensitization mechanisms, with the enzyme remaining membrane-associated and able to support sustained downstream signaling. DiC8-activated PKCα did not undergo dephosphorylation, ubiquitination, or internalization, early events in PKCα desensitization. Although DiC8 efficiently down-regulated novel PKCs PKCδ and PKCϵ, differences in Ca(2+) sensitivity and diacylglycerol affinity were excluded as mediators of the selective resistance of PKCα. Roles for Hsp/Hsc70 and Hsp90 were also excluded. PMA, but not DiC8, targeted PKCα to detergent-resistant membranes, and disruption of these domains with cholesterol-binding agents demonstrated a role for differential membrane compartmentalization in selective agonist-induced degradation. Chronic DiC8 treatment failed to desensitize PKCα in several cell types and did not affect PKCβI; thus, conventional PKCs appear generally insensitive to desensitization by sustained diacylglycerol stimulation. Consistent with this conclusion, prolonged (several-day) membrane association/activation of PKCα is seen in self-renewing epithelium of the intestine, cervix, and skin. PKCα deficiency affects gene expression, differentiation, and tumorigenesis in these tissues, highlighting the importance of mechanisms that protect PKCα from desensitization in vivo.

Structural insights into the interactions of phorbol ester and bryostatin complexed with protein kinase C: a comparative molecular dynamics simulation study

Protein kinase C (PKC) isozymes are important regulatory enzymes that have been implicated in many diseases, including cancer, Alzheimer's disease, and in the eradication of HIV/AIDS. Given their potential clinical ramifications, PKC modulators, e.g. phorbol esters and bryostatin, are also of great interest in the drug development. However, structural details on the binding between PKC and its modulators, especially bryostatin - the highly potent and non-tumor promoting activator for PKCs, are still lacking. Here, we report the first comparative molecular dynamics study aimed at gaining structural insight into the mechanisms by which the PKC delta cys2 activator domain is used in its binding to phorbol ester and bryostatin-1. As anticipated in the phorbol ester binding, hydrogen bonds are formed through the backbone atoms of Thr242, Leu251, and Gly253 of PKC. However, the opposition of H-bond formation between Thr242 and Gly253 may cause the phorbol ester complex to become less stable when compared with the bryostatin binding. For the PKC delta-bryostatin complex, hydrogen bonds are formed between the Gly253 backbone carbonyl and the C30 carbomethoxy substituent of the ligand. Additionally, the indole Nε1 of the highly homologous Trp252 also forms an H-bond to the C20 ester group on bryostatin. Backbone fluctuations also suggest that this latter H-bond formation may abrogate the transient interaction between Trp252 and His269, thus dampening the fluctuations observed on the nearby Zn(2+)-coordinating residues. This new dynamic fluctuation dampening model can potentially benefit future design of new PKC modulators.

Dynamics and Membrane Interactions of Protein Kinase C

Protein kinase C (PKC) is a family of Ser/Thr kinases that regulate a multitude of cellular processes through participation in the phosphoinositide signaling pathway. Significant research efforts have been directed at understanding the structure, function, and regulatory modes of the enzyme since its discovery and identification as the first receptor for tumor-promoting phorbol esters. The activation of PKC involves a transition from the cytosolic autoinhibited latent form to the membrane-associated active form. The membrane recruitment step is accompanied by the conformational rearrangement of the enzyme, which relieves autoinhibitory interactions and thereby allows PKC to phosphorylate its targets. The multidomain structure and intrinsic flexibility of PKC present remarkable challenges and opportunities for the biophysical and structural biology studies of this class of enzymes and their interactions with membranes, the major focus of this Current Topic. I will highlight the recent advances in the field, outline the current challenges, and identify areas where biophysics and structural biology approaches can provide insight into the isoenzyme-specific regulation of PKC activity.

Elucidating the interaction of γ-hydroxymethyl-γ-butyrolactone substituents with model membranes and protein kinase C-C1 domains

The protein kinase C (PKC) family of proteins is an attractive drug target. Dysregulation of PKC-dependent signalling pathways is related to several human diseases like cancer, immunological and other diseases. We approached the problem of altering PKC activities by developing C1 domain-based PKC ligands. In this report γ-hydroxymethyl-γ-butyrolactone (HGL) substituents were investigated in an effort to develop small molecule-based PKC regulators with higher specificity for C1 domain than the endogenous diacylglycerols (DAGs). Extensive analysis of membrane-ligands interaction measurements revealed that the membrane-active compounds strongly interact with the lipid bilayers and the hydrophilic parts of compounds localize at the bilayer/water interface. The pharmacophores like hydroxymethyl, carbonyl groups and acyl-chain length of the compounds are crucial for their interaction with the C1 domain proteins. The potent compounds showed more than 17-fold stronger binding affinity for the C1 domains than DAG under similar experimental conditions. Nonradioactive kinase assay confirmed that these potent compounds have similar or better PKC dependent phosphorylation capabilities than DAG under similar experimental conditions. Hence, our findings reveal that these HGL analogues represent an attractive group of structurally simple C1 domain ligands that can be further structurally altered to improve their potencies.

Toward a biorelevant structure of protein kinase C bound modulators: design, synthesis, and evaluation of labeled bryostatin analogues for analysis with rotational echo double resonance NMR spectroscopy

Protein kinase C (PKC) modulators are currently of great importance in preclinical and clinical studies directed at cancer, immunotherapy, HIV eradication, and Alzheimer's disease. However, the bound conformation of PKC modulators in a membrane environment is not known. Rotational echo double resonance (REDOR) NMR spectroscopy could uniquely address this challenge. However, REDOR NMR requires strategically labeled, high affinity ligands to determine interlabel distances from which the conformation of the bound ligand in the PKC-ligand complex could be identified. Here we report the first computer-guided design and syntheses of three bryostatin analogues strategically labeled for REDOR NMR analysis. Extensive computer analyses of energetically accessible analogue conformations suggested preferred labeling sites for the identification of the PKC-bound conformers. Significantly, three labeled analogues were synthesized, and, as required for REDOR analysis, all proved highly potent with PKC affinities (∼1 nM) on par with bryostatin. These potent and strategically labeled bryostatin analogues are new structural leads and provide the necessary starting point for projected efforts to determine the PKC-bound conformation of such analogues in a membrane environment, as needed to design new PKC modulators and understand PKC-ligand-membrane structure and dynamics.