Unidirectional allostery in the regulatory subunit RIα facilitates efficient deactivation of protein kinase A

Impaired cAMP processivity by phosphodiesterase-protein kinase A complexes in acrodysostosis

Acrodysostosis represents a group of rare genetic disorders characterized by defective skeletal development and is often accompanied by intellectual disabilities. Mutations in the 3'5'cyclic AMP (cAMP)-dependent protein kinase (PKA) type I regulatory subunit isoform α (RIα) and phosphodiesterase (PDE) PDE4D have both been implicated in impaired PKA regulation in acrodysostosis. How mutations on PDEs and RIα interfere with the regulation of cAMP-PKA signaling is not understood. cAMP-PKA signaling can be described in two phases. In the activation phase, cAMP binding to RIα dissociates the free C-subunit (Catalytic subunit). PDEs hydrolyze cAMP bound to RIα, priming the cAMP-free RIα for reassociation with the C-subunit, thereby completing one PKA activation cycle. Signal termination is thus critical for resetting PKA to its basal state and promoting adaptation to hormonal hyperstimulation. This proceeds through formation of a transient signal termination RIα: PDE complex that facilitates cAMP channeling from the cAMP-binding domain of RIα to the catalytic site of PDE. Signal termination of cAMP-PKA proceeds in three steps: Step 1) Channeling: translocation of cAMP from the CNB of RIα to the PDE catalytic site for hydrolysis. Step 2) Processivity: binding of free cAMP from the cytosol at both CNBs of RIα. Step 3) Product (5'AMP) release from the PDE hydrolysis site through competitive displacement by a new molecule of cAMP that triggers subsequent activation cycles of PKA. We have identified the molecular basis for two acrodysostosis mutants, PDE (PDE8 T690P) and RIα (T207A), that both allosterically impair cAMP-PKA signal termination. A combination of amide hydrogen/deuterium exchange mass spectrometry (HDXMS) and fluorescence polarization (FP) reveals that PDE8 T690P and RIα T207A both blocked processive hydrolysis of cAMP by interfering with competitive displacement of product 5'AMP release from the nucleotide channel at the end of each round of cAMP hydrolysis. While T690P blocked product 5'AMP release from the PDE, T207A greatly slowed the release of the substrate from RIα. These results highlight the role of processivity in cAMP hydrolysis by RIα: PDE termination complexes for adaptation to cAMP from GPCR hyperstimulation. Impairment of the signal termination process provides an alternate molecular basis for acrodysostosis.

Conformational changes in Protein Kinase A along its activation cycle are rooted in the folding energetics of cyclic-nucleotide binding domains

Cyclic-nucleotide binding (CNB) domains are structurally and evolutionarily conserved signaling modules that regulate proteins with diverse folds and functions. Despite a wealth of structural information, the mechanisms by which CNB domains couple cyclic-nucleotide binding to conformational changes involved in signal transduction remain unknown. Here we combined single-molecule and computational approaches to investigate the conformation and folding energetics of the two CNB domains of the regulatory subunit of protein kinase A (PKA). We found that the CNB domains exhibit different conformational and folding signatures in the apo state, when bound to cAMP, or when bound to the PKA catalytic subunit, underscoring their ability to adapt to different binding partners. Moreover, we show while the two CNB domains have near-identical structures, their thermodynamic coupling signatures are divergent, leading to distinct cAMP responses and differential mutational effects. Specifically, we demonstrate the mutation W260A exerts local and allosteric effects that impact multiple steps of the PKA activation cycle. Taken together, these results highlight the complex interplay between folding energetics, conformational dynamics, and thermodynamic signatures that underlies structurally conserved signaling modules in response to ligand binding and mutational effects.

Mechanistic insight into the destabilization of p53TD tetramer by cancer-related R337H mutation: a molecular dynamics study

The p53 protein is a tumor suppressor crucial for cell cycle and genome integrity. In a very large proportion of human cancers, p53 is frequently inactivated by mutations located in its DNA-binding domain (DBD). Some experimental studies reported that the inherited R337H mutation located in the p53 tetramerization domain (p53TD) can also result in destabilization of the p53 protein, and consequently lead to an organism prone to cancer setup. However, the underlying R337H mutation-induced structural destabilization mechanism is not well understood. Herein, we investigate the structural stability and dynamic property of the wild type p53TD tetramer and its cancer-related R337H mutant by performing multiple microsecond molecular dynamics simulations. It is found that R337H mutation destroys the R337-D352 hydrogen bonds, weakens the F341-F341 π-π stacking interaction and the hydrophobic interaction between aliphatic hydrocarbons of R337 and M340, leading to more solvent exposure of all the hydrophobic cores, and thus disrupting the structural integrity of the tetramer. Importantly, our simulations show for the first time that R337H mutation results in unfolding of the α-helix starting from the N-terminal region (residues ³³⁵RER(H)FEM³⁴⁰). Consistently, community network analyses reveal that R337H mutation reduces dynamical correlation and global connectivity of p53TD tetramer, which destabilizes the structure of the p53TD tetramer. This study provides the atomistic mechanism of R337H mutation-induced destabilization of p53TD tetramer, which might be helpful for in-depth understanding of the p53 loss-of-function mechanism.

Unraveling the Allosteric Mechanism of Four Cancer-related Mutations in the Disruption of p53-DNA Interaction

The p53 protein plays active roles in the physiological regulation of cell cycle as well as in cancer developments. In more than half of human cancers, the protein is inactivated by mutations located primarily in its DNA-binding domain (DBD), and some mutations located in the β-sandwich region of DBD are reported to decrease p53-DNA binding affinities. To understand the long-range correlation between p53 β-sandwich and DNA, and the allosteric mechanism of β-sandwich mutations in the disruption of p53-DNA interactions, we first identify three regions with a strong correlation with DNA based on microsecond molecular dynamics (MD) simulations of wild-type p53-DNA complex and then perform multiple MD simulations on four cancer-related mutants L145Q, P151S, Y220C, and G266R, which are located in these three regions. Our simulations show that these mutations allosterically destabilize the structural stability of the DNA-binding groove in p53 and disrupt the p53-DNA interactions. Network analyses reveal optimal correlation paths through which the mutation-induced allosteric signal passes to DNA, and the disturbance effect of these mutations on the global connectivity and dynamical correlation of the p53-DNA complex. This work paves the way for the in-depth understanding of the mutation-induced loss in p53's DNA-recognition ability and the pathological mechanism of cancer development.

Impacts of mutations on dynamic allostery of adenylate kinase

Escherichia coli adenylate kinase (AK) is composed of CORE domain and two branch domains: LID and AMP-binding domain (AMPbd). AK exhibits considerable allostery in a reversible phosphoryl transfer reaction, which is largely attributed to the relative motion of LID and AMPbd with respect to CORE. Such an allosteric conformational change is also evident in the absence of ligands. Recent studies showed that the mutations in branch domains can adjust dynamic allostery and alter the substrate affinity and enzyme activity. In this work, we use all-atom molecular dynamics simulation to study the impacts of mutations in branch domains on AK's dynamic allostery by comparing two double mutants, i.e., LID mutant (Val135Gly, Val142Gly) and AMPbd mutant (Ala37Gly, Ala55Gly), with wild-type. Two mutants undergo considerable conformational fluctuation and exhibit deviation from the initially extended apo state to more compact structures. The LID domain in the LID mutant adjusts its relative position and firmly adheres to CORE. More strikingly, AMPbd mutations affect the relative positions of both the AMPbd domain and remote LID domain. Both domains undergo considerable movement, especially the inherent hinge swing motion of the flexible LID domain. In both mutants, the mutations can enhance the inter-domain interaction. The results about the conformation change of AK in both mutants are in line with the experiment of AK's affinity and activity. As revealed by our findings, the flexibility of branch domains and their inherent motions, especially LID domain, is highly relevant to dynamic allostery in the AK system.

Adenylate control in cAMP signaling: implications for adaptation in signalosomes

In cAMP-Protein Kinase A (PKA) signaling, A-kinase anchoring protein scaffolds assemble PKA in close proximity to phosphodiesterases (PDE), kinase-substrates to form signaling islands or 'signalosomes'. In its basal state, inactive PKA holoenzyme (R2:C2) is activated by binding of cAMP to regulatory (R)-subunits leading to dissociation of active catalytic (C)-subunits. PDEs hydrolyze cAMP-bound to the R-subunits to generate 5'-AMP for termination and resetting the cAMP signaling. Mechanistic basis for cAMP signaling has been derived primarily by focusing on the proteins in isolation. Here, we set out to simulate cAMP signaling activation-termination cycles in a signalosome-like environment with PDEs and PKA subunits in close proximity to each other. Using a combination of fluorescence polarization and amide hydrogen exchange mass spectrometry with regulatory (RIα), C-subunit (Cα) and PDE8 catalytic domain, we have tracked movement of cAMP through activation-termination cycles. cAMP signaling operates as a continuum of four phases: (1) Activation and dissociation of PKA into R- and C-subunits by cAMP and facilitated by substrate (2) PDE recruitment to R-subunits (3) Hydrolysis of cAMP to 5'-AMP (4) Reassociation of C-subunit to 5'-AMP-bound-RIα in the presence of excess ATP to reset cAMP signaling to form the inactive PKA holoenzyme. Our results demonstrate that 5'-AMP is not merely a passive hydrolysis end-product of PDE action. A 'ligand-free' state R subunit does not exist in signalosomes as previously assumed. Instead the R-subunit toggles between cAMP- or 5'-AMP bound forms. This highlights, for the first time, the importance of 5'-AMP in promoting adaptation and uncovers adenylate control in cAMP signaling.

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.

Ligand-Binding-Site Structure Shapes Allosteric Signal Transduction and the Evolution of Allostery in Protein Complexes

The structure of ligand-binding sites has been shown to profoundly influence the evolution of function in homomeric protein complexes. Complexes with multichain binding sites (MBSs) have more conserved quaternary structure, more similar binding sites and ligands between homologs, and evolve new functions slower than homomers with single-chain binding sites (SBSs). Here, using in silico analyses of protein dynamics, we investigate whether ligand-binding-site structure shapes allosteric signal transduction pathways, and whether the structural similarity of binding sites influences the evolution of allostery. Our analyses show that: 1) allostery is more frequent among MBS complexes than in SBS complexes, particularly in homomers; 2) in MBS homomers, semirigid communities and critical residues frequently connect interfaces and thus they are characterized by signal transduction pathways that cross protein-protein interfaces, whereas SBS homomers usually not; 3) ligand binding alters community structure differently in MBS and SBS homomers; and 4) except MBS homomers, allosteric proteins are more likely to have homologs with similar binding site than nonallosteric proteins, suggesting that binding site similarity is an important factor driving the evolution of allostery.

Tuning the "violin" of protein kinases: The role of dynamics-based allostery

The intricacies of allosteric regulation of protein kinases continue to engage the research community. Allostery, or control from a distance, is seen as a fundamental biomolecular mechanism for proteins. From the traditional methods of conformational selection and induced fit, the field has grown to include the role of protein motions in defining a dynamics-based allosteric approach. Harnessing of these continuous motions in the protein to exert allosteric effects can be defined by a "violin" model that focuses on distributions of protein vibrations as opposed to concerted pathways. According to this model, binding of an allosteric modifier causes global redistribution of dynamics in the protein kinase domain that leads to changes in its catalytic properties. This model is consistent with the "entropy-driven allostery" mechanism proposed by Cooper and Dryden in 1984 and does not require, but does not exclude, any major structural changes. We provide an overview of practical implementation of the violin model and how it stands amidst the other known models of protein allostery. Protein kinases have been described as the biomolecules of interest. © 2019 IUBMB Life, 71(6):685-696, 2019.

Insights into the mechanism of the PIK3CA E545K activating mutation using MD simulations

Phosphoinositide 3-kinase alpha (PI3Kα) is involved in fundamental cellular processes including cell proliferation and differentiation and is frequently mutated in human malignancies. One of the most common mutations is E545K, which results in an amino acid substitution of opposite charge. It has been recently proposed that in this oncogenic charge-reversal mutation, the interactions between the protein catalytic and regulatory subunits are abrogated, resulting in loss of regulation and constitutive PI3Kα activity, which can lead to oncogenesis. To assess the mechanism of the PI3Kα E545K activating mutation, extensive Molecular Dynamics simulations were performed to examine conformational changes differing between the wild type (WT) and mutant proteins as they occur in microsecond simulations. In the E545K mutant PI3Kα, we observe a spontaneous detachment of the nSH2 PI3Kα domain (regulatory subunit, p85α) from the helical domain (catalytic subunit, p110α) causing significant loss of communication between the regulatory and catalytic subunits. We examine the allosteric network of the two proteins and show that a cluster of residues around the mutation is important for delivering communication signals between the catalytic and regulatory subunits. Our results demonstrate the dynamical and structural effects induced by the p110α E545K mutation in atomic level detail and indicate a possible mechanism for the loss of regulation that E545K confers on PI3Kα.

Self-organization, entropy and allostery

Allostery is a fundamental regulatory mechanism in biology. Although generally accepted that it is a dynamics-driven process, the exact molecular mechanism of allosteric signal transmission is hotly debated. We argue that allostery is as a part of a bigger picture that also includes fractal-like properties of protein interior, hierarchical protein folding and entropy-driven molecular recognition. Although so far all these phenomena were studied separately, they stem from the same common root: self-organization of polypeptide chains and, thus, has to be studied collectively. This merge will allow the cross-referencing of a broad spectrum of multi-disciplinary data facilitating progress in all these fields.

Expanding the Paradigm: Intrinsically Disordered Proteins and Allosteric Regulation

Allosteric regulatory processes are implicated at all levels of biological function. Recent advances in our understanding of the diverse and functionally significant class of intrinsically disordered proteins have identified a multitude of ways in which disordered proteins function within the confines of the allosteric paradigm. Allostery within or mediated by intrinsically disordered proteins ensures robust and efficient signal integration through mechanisms that would be extremely unfavorable or even impossible for globular protein interaction partners. Here, we highlight recent examples that indicate the breadth of biological outcomes that can be achieved through allosteric regulation by intrinsically disordered proteins. Ongoing and future work in this rapidly evolving area of research will expand our appreciation of the central role of intrinsically disordered proteins in ensuring the fidelity and efficiency of cellular regulation.

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.

Free energy landscape remodeling of the cardiac pacemaker channel explains the molecular basis of familial sinus bradycardia

The hyperpolarization-activated and cyclic nucleotide-modulated ion channel (HCN) drives the pacemaker activity in the heart, and its malfunction can result in heart disorders. One such disorder, familial sinus bradycardia, is caused by the S672R mutation in HCN, whose electrophysiological phenotypes include a negative shift in the channel activation voltage and an accelerated HCN deactivation. The outcomes of these changes are abnormally low resting heart rates. However, the molecular mechanism underlying these electrophysiological changes is currently not fully understood. Crystallographic investigations indicate that the S672R mutation causes limited changes in the structure of the HCN intracellular gating tetramer, but its effects on protein dynamics are unknown. Here, we utilize comparative S672R versus WT NMR analyses to show that the S672R mutation results in extensive perturbations of the dynamics in both apo- and holo-forms of the HCN4 isoform, reflecting how S672R remodels the free energy landscape for the modulation of HCN4 by cAMP, i.e. the primary cyclic nucleotide modulator of HCN channels. We show that the S672R mutation results in a constitutive shift of the dynamic auto-inhibitory equilibrium toward inactive states of HCN4 and broadens the free-energy well of the apo-form, enhancing the millisecond to microsecond dynamics of the holo-form at sites critical for gating cAMP binding. These S672R-induced variations in dynamics provide a molecular basis for the electrophysiological phenotypes of this mutation and demonstrate that the pathogenic effects of the S672R mutation can be rationalized primarily in terms of modulations of protein dynamics.