Using Markov state models to develop a mechanistic understanding of protein kinase A regulatory subunit RIα activation in response to cAMP binding

Electrophysical cardiac remodeling at the molecular level: Insights into ryanodine receptor activation and calcium-induced calcium release from a stochastic explicit-particle model

We present the first-ever, fully discrete, stochastic model of triggered cardiac Ca2+ dynamics. Using anatomically accurate subcellular cardiac myocyte geometries, we simulate the molecular players involved in Ca2+ handling using high-resolution stochastic and explicit-particle methods at the level of an individual cardiac dyadic junction. Integrating data from multiple experimental sources, the model not only replicates the findings of traditional in silico studies and complements in vitro experimental data but also reveals new insights into the molecular mechanisms driving cardiac dysfunction under stress and disease conditions. We improve upon older, nondiscrete models using the same realistic geometry by incorporating molecular mechanisms for spontaneous, as well as triggered calcium-induced calcium release (CICR). Action potentials are used to activate L-type calcium channels (LTCC), triggering CICR through ryanodine receptors (RyRs) on the surface of the sarcoplasmic reticulum. These improvements allow for the specific focus on the couplon: the structure-function relationship between LTCC and RyR. We investigate the electrophysical effects of normal and diseased action potentials on CICR and interrogate the effects of dyadic junction deformation through detubulation and orphaning of RyR. Our work demonstrates the importance of the electrophysical integrity of the calcium release unit on CICR fidelity, giving insights into the molecular basis of heart disease. Finally, we provide a unique, detailed, molecular view of the CICR process using advanced rendering techniques. This easy-to-use model comes complete with tutorials and the necessary software for use and analysis to maximize usability and reproducibility. Our work focuses on quantifying, qualifying, and visualizing the behavior of the molecular species that underlie the function and dysfunction of subcellular cardiomyocyte systems.

Fluorescent biosensors illuminate the spatial regulation of cell signaling across scales

As cell signaling research has advanced, it has become clearer that signal transduction has complex spatiotemporal regulation that goes beyond foundational linear transduction models. Several technologies have enabled these discoveries, including fluorescent biosensors designed to report live biochemical signaling events. As genetically encoded and live-cell compatible tools, fluorescent biosensors are well suited to address diverse cell signaling questions across different spatial scales of regulation. In this review, methods of examining spatial signaling regulation and the design of fluorescent biosensors are introduced. Then, recent biosensor developments that illuminate the importance of spatial regulation in cell signaling are highlighted at several scales, including membranes and organelles, molecular assemblies, and cell/tissue heterogeneity. In closing, perspectives on how fluorescent biosensors will continue enhancing cell signaling research 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.

Computational modeling approaches to cAMP/PKA signaling in cardiomyocytes

The cAMP/PKA pathway is a fundamental regulator of excitation-contraction coupling in cardiomyocytes. Activation of cAMP has a variety of downstream effects on cardiac function including enhanced contraction, accelerated relaxation, adaptive stress response, mitochondrial regulation, and gene transcription. Experimental advances have shed light on the compartmentation of cAMP and PKA, which allow for control over the varied targets of these second messengers and is disrupted in heart failure conditions. Computational modeling is an important tool for understanding the spatial and temporal complexities of this system. In this review article, we outline the advances in computational modeling that have allowed for deeper understanding of cAMP/PKA dynamics in the cardiomyocyte in health and disease, and explore new modeling frameworks that may bring us closer to a more complete understanding of this system. We outline various compartmental and spatial signaling models that have been used to understand how β-adrenergic signaling pathways function in a variety of simulation conditions. We also discuss newer subcellular models of cardiovascular function that may be used as templates for the next phase of computational study of cAMP and PKA in the heart, and outline open challenges which are important to consider in future models.

Structural analyses of the PKA RIIβ holoenzyme containing the oncogenic DnaJB1-PKAc fusion protein reveal protomer asymmetry and fusion-induced allosteric perturbations in fibrolamellar hepatocellular carcinoma

When the J-domain of the heat shock protein DnaJB1 is fused to the catalytic (C) subunit of cAMP-dependent protein kinase (PKA), replacing exon 1, this fusion protein, J-C subunit (J-C), becomes the driver of fibrolamellar hepatocellular carcinoma (FL-HCC). Here, we use cryo-electron microscopy (cryo-EM) to characterize J-C bound to RIIβ, the major PKA regulatory (R) subunit in liver, thus reporting the first cryo-EM structure of any PKA holoenzyme. We report several differences in both structure and dynamics that could not be captured by the conventional crystallography approaches used to obtain prior structures. Most striking is the asymmetry caused by the absence of the second cyclic nucleotide binding (CNB) domain and the J-domain in one of the RIIβ:J-C protomers. Using molecular dynamics (MD) simulations, we discovered that this asymmetry is already present in the wild-type (WT) RIIβ2C2 but had been masked in the previous crystal structure. This asymmetry may link to the intrinsic allosteric regulation of all PKA holoenzymes and could also explain why most disease mutations in PKA regulatory subunits are dominant negative. The cryo-EM structure, combined with small-angle X-ray scattering (SAXS), also allowed us to predict the general position of the Dimerization/Docking (D/D) domain, which is essential for localization and interacting with membrane-anchored A-Kinase-Anchoring Proteins (AKAPs). This position provides a multivalent mechanism for interaction of the RIIβ holoenzyme with membranes and would be perturbed in the oncogenic fusion protein. The J-domain also alters several biochemical properties of the RIIβ holoenzyme: It is easier to activate with cAMP, and the cooperativity is reduced. These results provide new insights into how the finely tuned allosteric PKA signaling network is disrupted by the oncogenic J-C subunit, ultimately leading to the development of FL-HCC.

Spatially compartmentalized phase regulation of a Ca2+-cAMP-PKA oscillatory circuit

Signaling networks are spatiotemporally organized to sense diverse inputs, process information, and carry out specific cellular tasks. In β cells, Ca2+, cyclic adenosine monophosphate (cAMP), and Protein Kinase A (PKA) exist in an oscillatory circuit characterized by a high degree of feedback. Here, we describe a mode of regulation within this circuit involving a spatial dependence of the relative phase between cAMP, PKA, and Ca2+. We show that in mouse MIN6 β cells, nanodomain clustering of Ca2+-sensitive adenylyl cyclases (ACs) drives oscillations of local cAMP levels to be precisely in-phase with Ca2+ oscillations, whereas Ca2+-sensitive phosphodiesterases maintain out-of-phase oscillations outside of the nanodomain. Disruption of this precise phase relationship perturbs Ca2+ oscillations, suggesting the relative phase within an oscillatory circuit can encode specific functional information. This work unveils a novel mechanism of cAMP compartmentation utilized for localized tuning of an oscillatory circuit and has broad implications for the spatiotemporal regulation of signaling networks.

Phase Separation of a PKA Regulatory Subunit Controls cAMP Compartmentation and Oncogenic Signaling

The fidelity of intracellular signaling hinges on the organization of dynamic activity architectures. Spatial compartmentation was first proposed over 30 years ago to explain how diverse G protein-coupled receptors achieve specificity despite converging on a ubiquitous messenger, cyclic adenosine monophosphate (cAMP). However, the mechanisms responsible for spatially constraining this diffusible messenger remain elusive. Here, we reveal that the type I regulatory subunit of cAMP-dependent protein kinase (PKA), RIα, undergoes liquid-liquid phase separation (LLPS) as a function of cAMP signaling to form biomolecular condensates enriched in cAMP and PKA activity, critical for effective cAMP compartmentation. We further show that a PKA fusion oncoprotein associated with an atypical liver cancer potently blocks RIα LLPS and induces aberrant cAMP signaling. Loss of RIα LLPS in normal cells increases cell proliferation and induces cell transformation. Our work reveals LLPS as a principal organizer of signaling compartments and highlights the pathological consequences of dysregulating this activity architecture.

Histamine-induced biphasic activation of RhoA allows for persistent RhoA signaling

The small GTPase RhoA is a central signaling enzyme that is involved in various cellular processes such as cytoskeletal dynamics, transcription, and cell cycle progression. Many signal transduction pathways activate RhoA—for instance, Gα_q-coupled Histamine 1 Receptor signaling via Gα_q-dependent activation of RhoGEFs such as p63. Although multiple upstream regulators of RhoA have been identified, the temporal regulation of RhoA and the coordination of different upstream components in its regulation have not been well characterized. In this study, live-cell measurement of RhoA activation revealed a biphasic increase of RhoA activity upon histamine stimulation. We showed that the first and second phase of RhoA activity are dependent on p63 and Ca²⁺/PKC, respectively, and further identified phosphorylation of serine 240 on p115 RhoGEF by PKC to be the mechanistic link between PKC and RhoA. Combined approaches of computational modeling and quantitative measurement revealed that the second phase of RhoA activation is insensitive to rapid turning off of the receptor and is required for maintaining RhoA-mediated transcription after the termination of the receptor signaling. Thus, two divergent pathways enable both rapid activation and persistent signaling in receptor-mediated RhoA signaling via intricate temporal regulation.

The small GTPase RhoA is a central signaling enzyme that is involved in various cellular processes such as cytoskeletal dynamics, transcription, and cell cycle progression. This study shows that histamine induces biphasic activation of RhoA via two divergent signaling pathways, allowing for intricate regulation of cellular processes.

Activation of PKA via asymmetric allosteric coupling of structurally conserved cyclic nucleotide binding domains

Cyclic nucleotide-binding (CNB) domains allosterically regulate the activity of proteins with diverse functions, but the mechanisms that enable the cyclic nucleotide-binding signal to regulate distant domains are not well understood. Here we use optical tweezers and molecular dynamics to dissect changes in folding energy landscape associated with cAMP-binding signals transduced between the two CNB domains of protein kinase A (PKA). We find that the response of the energy landscape upon cAMP binding is domain specific, resulting in unique but mutually coordinated tasks: one CNB domain initiates cAMP binding and cooperativity, whereas the other triggers inter-domain interactions that promote the active conformation. Inter-domain interactions occur in a stepwise manner, beginning in intermediate-liganded states between apo and cAMP-bound domains. Moreover, we identify a cAMP-responsive switch, the N3A motif, whose conformation and stability depend on cAMP occupancy. This switch serves as a signaling hub, amplifying cAMP-binding signals during PKA activation.

Two PKA RIα holoenzyme states define ATP as an isoform-specific orthosteric inhibitor that competes with the allosteric activator, cAMP

Protein kinase A (PKA) holoenzyme, comprised of a cAMP-binding regulatory (R)-subunit dimer and 2 catalytic (C)-subunits, is the master switch for cAMP-mediated signaling. Of the 4 R-subunits (RIα, RIβ, RIIα, RIIβ), RIα is most essential for regulating PKA activity in cells. Our 2 RIα2C2 holoenzyme states, which show different conformations with and without ATP, reveal how ATP/Mg2+ functions as a negative orthosteric modulator. Biochemical studies demonstrate how the removal of ATP primes the holoenzyme for cAMP-mediated activation. The opposing competition between ATP/cAMP is unique to RIα. In RIIβ, ATP serves as a substrate and facilitates cAMP-activation. The isoform-specific RI-holoenzyme dimer interface mediated by N3A-N3A' motifs defines multidomain cross-talk and an allosteric network that creates competing roles for ATP and cAMP. Comparisons to the RIIβ holoenzyme demonstrate isoform-specific holoenzyme interfaces and highlights distinct allosteric mechanisms for activation in addition to the structural diversity of the isoforms.

Markov models for the elucidation of allosteric regulation

Allosteric regulation refers to the process where the effect of binding of a ligand at one site of a protein is transmitted to another, often distant, functional site. In recent years, it has been demonstrated that allosteric mechanisms can be understood by the conformational ensembles of a protein. Molecular dynamics (MD) simulations are often used for the study of protein allostery as they provide an atomistic view of the dynamics of a protein. However, given the wealth of detailed information hidden in MD data, one has to apply a method that allows extraction of the conformational ensembles underlying allosteric regulation from these data. Markov state models are one of the most promising methods for this purpose. We provide a short introduction to the theory of Markov state models and review their application to various examples of protein allostery studied by MD simulations. We also include a discussion of studies where Markov modelling has been employed to analyse experimental data on allosteric regulation. We conclude our review by advertising the wider application of Markov state models to elucidate allosteric mechanisms, especially since in recent years it has become straightforward to construct such models thanks to software programs like PyEMMA and MSMBuilder.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

Norepinephrine stimulates glycogenolysis in astrocytes to fuel neurons with lactate

The mechanism of rapid energy supply to the brain, especially to accommodate the heightened metabolic activity of excited states, is not well-understood. We explored the role of glycogen as a fuel source for neuromodulation using the noradrenergic stimulation of glia in a computational model of the neural-glial-vasculature ensemble (NGV). The detection of norepinephrine (NE) by the astrocyte and the coupled cAMP signal are rapid and largely insensitive to the distance of the locus coeruleus projection release sites from the glia, implying a diminished impact for volume transmission in high affinity receptor transduction systems. Glucosyl-conjugated units liberated from glial glycogen by NE-elicited cAMP second messenger transduction winds sequentially through the glycolytic cascade, generating robust increases in NADH and ATP before pyruvate is finally transformed into lactate. This astrocytic lactate is rapidly exported by monocarboxylate transporters to the associated neuron, demonstrating that the astrocyte-to-neuron lactate shuttle activated by glycogenolysis is a likely fuel source for neuromodulation and enhanced neural activity. Altogether, the energy supply for both astrocytes and neurons can be supplied rapidly by glycogenolysis upon neuromodulatory stimulus.

Although efficient compared to computers, the human brain utilizes energy at 10-fold the rate of other organs by mass. How the brain is supplied with sufficient on-demand energy to support its activity in the absence of neuronal storage capacity remains unknown. Neurons are not capable of meeting their own energy requirements, instead energy supply in the brain is managed by an oligocellular cartel composed of neurons, glia and the local vasculature (NGV), wherein glia can provide the ergogenic metabolite lactate to the neuron in a process called the astrocyte-to-neuron shuttle (ANLS). The only means of energy storage in the brain is glycogen, a polymerized form of glucose that is localized largely to astrocytes, but its exact role and conditions of use are not clear. In this computational model we show that neuromodulatory stimulation by norepinephrine induces astrocytes to recover glucosyl subunits from glycogen for use in a glycolytic process that favors the production of lactate. The ATP and NADH produced support metabolism in the astrocyte while the lactate is exported to feed the neuron. Thus, rapid energy demands by both neurons and glia in a stimulated brain can be met by glycogen mobilization.

Switching of the folding-energy landscape governs the allosteric activation of protein kinase A

Protein kinases are dynamic molecular switches that sample multiple conformational states. The regulatory subunit of PKA harbors two cAMP-binding domains [cyclic nucleotide-binding (CNB) domains] that oscillate between inactive and active conformations dependent on cAMP binding. The cooperative binding of cAMP to the CNB domains activates an allosteric interaction network that enables PKA to progress from the inactive to active conformation, unleashing the activity of the catalytic subunit. Despite its importance in the regulation of many biological processes, the molecular mechanism responsible for the observed cooperativity during the activation of PKA remains unclear. Here, we use optical tweezers to probe the folding cooperativity and energetics of domain communication between the cAMP-binding domains in the apo state and bound to the catalytic subunit. Our study provides direct evidence of a switch in the folding-energy landscape of the two CNB domains from energetically independent in the apo state to highly cooperative and energetically coupled in the presence of the catalytic subunit. Moreover, we show that destabilizing mutational effects in one CNB domain efficiently propagate to the other and decrease the folding cooperativity between them. Taken together, our results provide a thermodynamic foundation for the conformational plasticity that enables protein kinases to adapt and respond to signaling molecules.

PKA-RII subunit phosphorylation precedes activation by cAMP and regulates activity termination

Activity of endogenous protein kinase A (PKA) could never be analyzed directly in the cellular environment. Isensee et al. used antibodies to quantify conformational changes leading to an open conformation of endogenous PKA-II holoenzymes, which allowed them to analyze and model its activation cycle in primary sensory neurons.

Type II isoforms of cyclic adenosine monophosphate (cAMP)–dependent protein kinase A (PKA-II) contain a phosphorylatable epitope within the inhibitory domain of RII subunits (pRII) with still unclear function. In vitro, RII phosphorylation occurs in the absence of cAMP, whereas staining of cells with pRII-specific antibodies revealed a cAMP-dependent pattern. In sensory neurons, we found that increased pRII immunoreactivity reflects increased accessibility of the already phosphorylated RII epitope during cAMP-induced opening of the tetrameric RII₂:C₂ holoenzyme. Accordingly, induction of pRII by cAMP was sensitive to novel inhibitors of dissociation, whereas blocking catalytic activity was ineffective. Also in vitro, cAMP increased the binding of pRII antibodies to RII₂:C₂ holoenzymes. Identification of an antibody specific for the glycine-rich loop of catalytic subunits facing the pRII-epitope confirmed activity-dependent binding with similar kinetics, proving that the reassociation is rapid and precisely controlled. Mechanistic modeling further supported that RII phosphorylation precedes cAMP binding and controls the inactivation by modulating the reassociation involving the coordinated action of phosphodiesterases and phosphatases.

Phosphorylation of protein kinase A (PKA) regulatory subunit RIα by protein kinase G (PKG) primes PKA for catalytic activity in cells

cAMP-dependent protein kinase (PKAc) is a pivotal signaling protein in eukaryotic cells. PKAc has two well-characterized regulatory subunit proteins, RI and RII (each having α and β isoforms), which keep the PKAc catalytic subunit in a catalytically inactive state until activation by cAMP. Previous reports showed that the RIα regulatory subunit is phosphorylated by cGMP-dependent protein kinase (PKG) in vitro, whereupon phosphorylated RIα no longer inhibits PKAc at normal (1:1) stoichiometric ratios. However, the significance of this phosphorylation as a mechanism for activating type I PKA holoenzymes has not been fully explored, especially in cellular systems. In this study, we further examined the potential of RIα phosphorylation to regulate physiologically relevant "desensitization" of PKAc activity. First, the serine 101 site of RIα was validated as a target of PKGIα phosphorylation both in vitro and in cells. Analysis of a phosphomimetic substitution in RIα (S101E) showed that modification of this site increases PKAc activity in vitro and in cells, even without cAMP stimulation. Numerous techniques were used to show that although Ser101 variants of RIα can bind PKAc, the modified linker region of the S101E mutant has a significantly reduced affinity for the PKAc active site. These findings suggest that RIα phosphorylation may be a novel mechanism to circumvent the requirement of cAMP stimulus to activate type I PKA in cells. We have thus proposed a model to explain how PKG phosphorylation of RIα creates a "sensitized intermediate" state that is in effect primed to trigger PKAc activity.

The structure of the apo cAMP-binding domain of HCN4 - a stepping stone toward understanding the cAMP-dependent modulation of the hyperpolarization-activated cyclic-nucleotide-gated ion channels

The hyperpolarization-activated cyclic-nucleotide-gated (HCN) ion channels control nerve impulse transmission and cardiac pacemaker activity. The modulation by cAMP is critical for the regulatory function of HCN in both neurons and cardiomyocytes, but the underlying mechanism is not fully understood. Here, we show how the structure of the apo cAMP-binding domain of the HCN4 isoform has contributed to a model for the cAMP-dependent modulation of the HCN ion-channel. This model recapitulates the structural and dynamical changes that occur along the thermodynamic cycle arising from the coupling of cAMP-binding and HCN self-association equilibria. The proposed model addresses some of the questions previously open about the auto-inhibition of HCN and its cAMP-induced activation, while opening new opportunities for selectively targeting HCN through allosteric ligands. A remaining challenge is the investigation of HCN dimers and their regulatory role. Overcoming this challenge will require the integration of crystallography, cryo electron microscopy, NMR, electrophysiology and simulations.

Switching Cyclic Nucleotide-Selective Activation of Cyclic Adenosine Monophosphate-Dependent Protein Kinase Holoenzyme Reveals Distinct Roles of Tandem Cyclic Nucleotide-Binding Domains

The cyclic adenosine monophosphate (cAMP)- and cyclic guanosine monophosphate (cGMP)-dependent protein kinases (PKA and PKG) are key effectors of cyclic nucleotide signaling. Both share structural features that include tandem cyclic nucleotide-binding (CNB) domains, CNB-A and CNB-B, yet their functions are separated through preferential activation by either cAMP or cGMP. Based on structural studies and modeling, key CNB contact residues have been identified for both kinases. In this study, we explored the requirements for conversion of PKA activation from cAMP-dependent to cGMP-dependent. The consequences of the residue substitutions T192R/A212T within CNB-A or G316R/A336T within CNB-B of PKA-RIα on cyclic nucleotide binding and holoenzyme activation were assessed in vitro using purified recombinant proteins, and ex vivo using RIα-deficient mouse embryonic fibroblasts genetically reconstituted with wild-type or mutant PKA-RIα. In vitro, a loss of binding and activation selectivity was observed when residues in either one of the CNB domains were mutated, while mutations in both CNB domains resulted in a complete switch of selectivity from cAMP to cGMP. The switch in selectivity was also recapitulated ex vivo, confirming their functional roles in cells. Our results highlight the importance of key cyclic nucleotide contacts within each CNB domain and suggest that these domains may have evolved from an ancestral gene product to yield two distinct cyclic nucleotide-dependent protein kinases.

cAMP-Dependent Protein Kinase and cGMP-Dependent Protein Kinase as Cyclic Nucleotide Effectors

The cAMP-dependent protein kinase (PKA) and the cGMP-dependent protein kinase (PKG) are homologous enzymes with different binding and activation specificities for cyclic nucleotides. Both enzymes harbor conserved cyclic nucleotide-binding (CNB) domains. Differences in amino acid composition of these CNB domains mediate cyclic nucleotide selectivity in PKA and PKG, respectively. Recently, the presence of the noncanonical cyclic nucleotides cCMP and cUMP in eukaryotic cells has been proven, while the existence of cellular cIMP and cXMP remains unclear. It was shown that the main effectors of cyclic nucleotide signaling, PKA and PKG, can be activated by each of these noncanonical cyclic nucleotides. With unique effector proteins still missing, such cross-activation effects might have physiological relevance. Therefore, we approach PKA and PKG as cyclic nucleotide effectors in this chapter. The focus of this chapter is the general cyclic nucleotide-binding properties of both kinases as well as the selectivity for cAMP or cGMP, respectively. Furthermore, we discuss the binding affinities and activation potencies of noncanonical cyclic nucleotides.

Molecular Simulations Reveal an Unresolved Conformation of the Type IA Protein Kinase A Regulatory Subunit and Suggest Its Role in the cAMP Regulatory Mechanism

We identify a previously unresolved, unrecognized, and highly stable conformation of the protein kinase A (PKA) regulatory subunit RIα. This conformation, which we term the "Flipback" structure, bridges conflicting characteristics in crystallographic structures and solution experiments of the PKA RIα heterotetramer. Our simulations reveal a hinge residue, G235, in the B/C helix that is conserved through all isoforms of RI. Brownian dynamics simulations suggest that the Flipback conformation plays a role in cAMP association to the A domain of the R subunit.

Tyrosine Kinase Activation and Conformational Flexibility: Lessons from Src-Family Tyrosine Kinases

Protein kinases are enzymes that catalyze the covalent transfer of the γ-phosphate of an adenosine triphosphate (ATP) molecule onto a tyrosine, serine, threonine, or histidine residue in the substrate and thus send a chemical signal to networks of downstream proteins. They are important cellular signaling enzymes that regulate cell growth, proliferation, metabolism, differentiation, and migration. Unregulated protein kinase activity is often associated with a wide range of diseases, therefore making protein kinases major therapeutic targets. A prototypical system of central interest to understand the regulation of kinase activity is provided by tyrosine kinase c-Src, which belongs to the family of Src-related non-receptor tyrosine kinases (SFKs). Although the broad picture of autoinhibition via the regulatory domains and via the phosphorylation of the C-terminal tail is well characterized from a structural point of view, a detailed mechanistic understanding at the atomic-level is lacking. Advanced computational methods based on all-atom molecular dynamics (MD) simulations are employed to advance our understanding of tyrosine kinase activation. The computational studies suggest that the isolated kinase domain (KD) is energetically most favorable in the inactive conformation when the activation loop (A-loop) of the KD is not phosphorylated. The KD makes transient visits to a catalytically competent active-like conformation. The process of bimolecular trans-autophosphorylation of the A-loop eventually locks the KD in the active state. Activating point mutations may act by slightly increasing the population of the active-like conformation, enhancing the availability of the A-loop to be phosphorylated. The Src-homology 2 (SH2) and Src-homology 3 (SH3) regulatory domains, depending upon their configuration, either promote the inactive or the active state of the kinase domain. In addition to the roles played by the SH3, SH2, and KD, the Src-homology 4-Unique domain (SH4-U) region also serves as a key moderator of substrate specificity and kinase function. Thus, a fundamental understanding of the conformational propensity of the SH4-U region and how this affects the association to the membrane surface are likely to lead to the discovery of new intermediate states and alternate strategies for inhibition of kinase activity for drug discovery. The existence of a multitude of KD conformations poses a great challenge aimed at the design of specific inhibitors. One promising computational strategy to explore the conformational flexibility of the KD is to construct Markov state models from aggregated MD data.

The quantum basis of spatiotemporality in perception and consciousness

Living systems inhabit the area of the world which is shaped by the predictable space-time of physical objects and forces that can be incorporated into their perception pattern. The process of selecting a "habitable" space-time is the internal quantum measurement in which living systems become embedded into the environment that supports their living state. This means that living organisms choose a coordinate system in which the influence of measurement is minimal. We discuss specific roles of biological macromolecules, in particular of the cytoskeleton, in shaping perception patterns formed in the internal measurement process. Operation of neuron is based on the transmission of signals via cytoskeleton where the digital output is generated that can be decoded through a reflective action of the perceiving agent. It is concluded that the principle of optimality in biology as formulated by Liberman et al. (BioSystems 22, 135-154, 1989) is related to the establishment of spatiotemporal patterns that are maximally predictable and can hold the living state for a prolonged time. This is achieved by the selection of a habitable space approximated to the conditions described by classical physics.

Structure of a PKA RIα Recurrent Acrodysostosis Mutant Explains Defective cAMP-Dependent Activation

Most disease-related mutations that impair cAMP protein kinase A (PKA) signaling are present within the regulatory (R) PKA RI alpha-subunit (RIα). Although mutations in the PRKAR1A gene are linked to Carney complex (CNC) disease and, more recently, to acrodysostosis-1 (ACRDYS1), the two diseases show contrasting phenotypes. While CNC mutations cause increased PKA activity, ACRDYS1 mutations result in decreased PKA activity and cAMP resistant holoenzymes. Mapping the ACRDYS1 disease mutations reveals their localization to the second of two tandem cAMP-binding (CNB) domains (CNB-B), and here, we characterize a recurrent deletion mutant where the last 14 residues are missing. The crystal structure of a monomeric form of this mutant (RIα92-365) bound to the catalytic (C)-subunit reveals the dysfunctional regions of the RIα subunit. Beyond the missing residues, the entire capping motif is disordered (residues 357-379) and explains the disrupted cAMP binding. Moreover, the effects of the mutation extend far beyond the CNB-B domain and include the active site and N-lobe of the C-subunit, which is in a partially open conformation with the C-tail disordered. A key residue that contributes to this crosstalk, D267, is altered in our structure, and we confirmed its functional importance by mutagenesis. In particular, the D267 interaction with Arg241, a residue shown earlier to be important for allosteric regulation, is disrupted, thereby strengthening the interaction of D267 with the C-subunit residue Arg194 at the R:C interface. We see here how the switch between active (cAMP-bound) and inactive (holoenzyme) conformations is perturbed and how the dynamically controlled crosstalk between the helical domains of the two CNB domains is necessary for the functional regulation of PKA activity.

Transition path theory analysis of c-Src kinase activation

Nonreceptor tyrosine kinases of the Src family are large multidomain allosteric proteins that are crucial to cellular signaling pathways. In a previous study, we generated a Markov state model (MSM) to simulate the activation of c-Src catalytic domain, used as a prototypical tyrosine kinase. The long-time kinetics of transition predicted by the MSM was in agreement with experimental observations. In the present study, we apply the framework of transition path theory (TPT) to the previously constructed MSM to characterize the main features of the activation pathway. The analysis indicates that the activating transition, in which the activation loop first opens up followed by an inward rotation of the αC-helix, takes place via a dense set of intermediate microstates distributed within a fairly broad "transition tube" in a multidimensional conformational subspace connecting the two end-point conformations. Multiple microstates with negligible equilibrium probabilities carry a large transition flux associated with the activating transition, which explains why extensive conformational sampling is necessary to accurately determine the kinetics of activation. Our results suggest that the combination of MSM with TPT provides an effective framework to represent conformational transitions in complex biomolecular systems.

Mapping the Free Energy Landscape of PKA Inhibition and Activation: A Double-Conformational Selection Model for the Tandem cAMP-Binding Domains of PKA RIα

Protein Kinase A (PKA) is the major receptor for the cyclic adenosine monophosphate (cAMP) secondary messenger in eukaryotes. cAMP binds to two tandem cAMP-binding domains (CBD-A and -B) within the regulatory subunit of PKA (R), unleashing the activity of the catalytic subunit (C). While CBD-A in RIα is required for PKA inhibition and activation, CBD-B functions as a “gatekeeper” domain that modulates the control exerted by CBD-A. Preliminary evidence suggests that CBD-B dynamics are critical for its gatekeeper function. To test this hypothesis, here we investigate by Nuclear Magnetic Resonance (NMR) the two-domain construct RIα (91–379) in its apo, cAMP₂, and C-bound forms. Our comparative NMR analyses lead to a double conformational selection model in which each apo CBD dynamically samples both active and inactive states independently of the adjacent CBD within a nearly degenerate free energy landscape. Such degeneracy is critical to explain the sensitivity of CBD-B to weak interactions with C and its high affinity for cAMP. Binding of cAMP eliminates this degeneracy, as it selectively stabilizes the active conformation within each CBD and inter-CBD contacts, which require both cAMP and W260. The latter is contributed by CBD-B and mediates capping of the cAMP bound to CBD-A. The inter-CBD interface is dispensable for intra-CBD conformational selection, but is indispensable for full activation of PKA as it occludes C-subunit recognition sites within CBD-A. In addition, the two structurally homologous cAMP-bound CBDs exhibit marked differences in their residual dynamics profiles, supporting the notion that conservation of structure does not necessarily imply conservation of dynamics.

Protein Kinase A (PKA) is the major receptor for the cAMP secondary messenger in eukaryotes. This study shows how PKA's regulatory subunit dynamically samples a degenerate free energy landscape that controls affinities for the catalytic subunit and cAMP; intra-domain conformational selection by cAMP controls inter-domain interactions and PKA activation.

Cyclic adenosine monophosphate (cAMP) is a messenger molecule produced within cells to control cellular metabolism in response to external stimuli. Protein Kinase A (PKA) is the major receptor for cAMP. cAMP binds to tandem cAMP-binding domains (CBD-A and -B) within the regulatory subunits of PKA (R), unleashing the activity of the catalytic subunit (C). While CBD-A is required for C-subunit inhibition and activation, in RIα CBD-B functions as a “gatekeeper” domain that modulates the control exerted by CBD-A. However, it is not currently clear how ligand binding and dynamics of CBD-B mediate its gatekeeper function. We comparatively analyzed by Nuclear Magnetic Resonance (NMR) a two-domain construct of the regulatory subunit RIα with no ligand, with cAMP₂ bound, and the C-bound form. These data show that both CBDs can exist in a system of uncorrelated conformational selection as both can independently sample activated and inactivated states (in what is known as a nearly degenerate free energy landscape). This explains why both RIα CBDs exhibit a higher cAMP-affinity than other cAMP receptors. Once cAMP has bound, the degeneracy is lost and dissociation of the kinase subunit is promoted through a combination of intra-domain conformational selection and changes in inter-CBD orientation. The proposed model—a double-conformational selection model—provides a general framework to interpret the effect of PKA mutations that have been reported in rare human disorders such as Carney complex and Acrodysostosis.

Bridging scales through multiscale modeling: a case study on protein kinase A

The goal of multiscale modeling in biology is to use structurally based physico-chemical models to integrate across temporal and spatial scales of biology and thereby improve mechanistic understanding of, for example, how a single mutation can alter organism-scale phenotypes. This approach may also inform therapeutic strategies or identify candidate drug targets that might otherwise have been overlooked. However, in many cases, it remains unclear how best to synthesize information obtained from various scales and analysis approaches, such as atomistic molecular models, Markov state models (MSM), subcellular network models, and whole cell models. In this paper, we use protein kinase A (PKA) activation as a case study to explore how computational methods that model different physical scales can complement each other and integrate into an improved multiscale representation of the biological mechanisms. Using measured crystal structures, we show how molecular dynamics (MD) simulations coupled with atomic-scale MSMs can provide conformations for Brownian dynamics (BD) simulations to feed transitional states and kinetic parameters into protein-scale MSMs. We discuss how milestoning can give reaction probabilities and forward-rate constants of cAMP association events by seamlessly integrating MD and BD simulation scales. These rate constants coupled with MSMs provide a robust representation of the free energy landscape, enabling access to kinetic, and thermodynamic parameters unavailable from current experimental data. These approaches have helped to illuminate the cooperative nature of PKA activation in response to distinct cAMP binding events. Collectively, this approach exemplifies a general strategy for multiscale model development that is applicable to a wide range of biological problems.

Allostery through the computational microscope: cAMP activation of a canonical signalling domain

Ligand-induced protein allostery plays a central role in modulating cellular signalling pathways. Here using the conserved cyclic nucleotide-binding domain of protein kinase A's (PKA) regulatory subunit as a prototype signalling unit, we combine long-timescale, all-atom molecular dynamics simulations with Markov state models to elucidate the conformational ensembles of PKA's cyclic nucleotide-binding domain A for the cAMP-free (apo) and cAMP-bound states. We find that both systems exhibit shallow free-energy landscapes that link functional states through multiple transition pathways. This observation suggests conformational selection as the general mechanism of allostery in this canonical signalling domain. Further, we expose the propagation of the allosteric signal through key structural motifs in the cyclic nucleotide-binding domain and explore the role of kinetics in its function. Our approach integrates disparate lines of experimental data into one cohesive framework to understand structure, dynamics and function in complex biological systems.

Discovery of Allostery in PKA Signaling

cAMP-dependent protein kinase (PKA) was the second protein kinase to be discovered and the PKA catalytic (C) subunit serves as a prototype for the large protein kinase superfamily that contains over 500 gene products. The protein kinases regulate much of biology in eukaryotic cells and they are now also a major therapeutic target. Although PKA was discovered nearly 50 years ago and the subsequent discovery of the regulatory subunits that bind cAMP and release the catalytic activity from the holoenzyme followed quickly. Thus in PKA we see the convergence of two major signaling mechanisms - protein phosphorylation and second messenger signaling through cAMP. Crystallography provides a foundation for understanding function, and the structure of the isolated regulatory (R) and C-subunits have been extremely informative. Yet it is the R₂C₂ holoenzyme that predominates in cells, and one can only appreciate the allosteric features of PKA signaling by seeing the full length protein. The symmetry and the quaternary constraints that one R:C hetero-dimer exerts on the other in the holoenzyme simply are not present in the isolated subunits or even in the R:C hetero-dimer.