Globally correlated conformational entropy underlies positive and negative cooperativity in a kinase's enzymatic cycle

Unbiased clustering of residues undergoing synchronous motions in proteins using NMR spin relaxation data

Biological macromolecules are dynamic entities that transition between various conformational states, often playing a vital role in biological functions. Their inherent flexibility spans a broad range of timescales. Motions occurring within the microsecond to millisecond range are especially important, as they are integral to processes such as enzyme catalysis, folding, ligand binding, and allostery. NMR Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion measurements are the preferred method for characterizing macromolecular motions at atomic resolution. However, it is still uncertain whether the functional motions of multiple residues in macromolecules need to be coordinated and/or synchronized within the protein matrix in order to perform the desired function. Here, we illustrate an unbiased method to analyze NMR relaxation dispersion and identify dynamic clusters of residues that fluctuate on similar timescales within proteins. The method requires relaxation dispersion data for backbone amides or side-chain methyl groups, which are globally fitted using the Bloch-McConnell equations for each pair of residues. The goodness of the pairwise fitting serves as a metric to construct two-dimensional synchronous dynamics (SyncDyn) maps, allowing us to identify residue clusters whose dynamics are influenced by ligand binding. We applied our method to the catalytic subunit of the cAMP-dependent protein kinase A (PKAC) and the T17A mutant of ribonuclease A (RNAse A). The SyncDyn maps for PKA-C showed distinct clusters of residues located in critical allosteric sites. Nucleotide binding activates the movement of residues at the interface between the two lobes and also those distal to the active site. In the case of RNAse A, the SyncDyn maps show that residues fluctuating with the same time scale are interspersed in both lobes of the enzyme. Overall, our approach eliminates arbitrary manual selection of residues for dynamic clustering and objectively identifies all possible residue pairs that fluctuate synchronously, i.e. on the same timescale.

Positive activation entropy of Bacillus circulans xylanase catalyzed ONPX2 hydrolysis: A mechanistic and engineering study

Transition state (TS) stabilization by enzymes greatly accelerates catalytic reactions. For some enzymes, the TS complex has entropy higher than enzyme substrate (ES) complex. But the origin of favorable entropy remains unclear. In this work, we studied the mechanism of Bacillus Circulans xylanase (BCX) 11 catalyzed o-nitrophenyl β-xylobioside (ONPX2) glycoside hydrolysis. The catalytic reaction exhibits a positive activation entropy, and an increase in ionic strength leads to a decrease in entropy without affecting the activation free energy, indicating that the entropy is predominantly influenced by electrostatic forces. Moreover, NMR measurements of electrostatic attractions within the active site demonstrate a positive entropy, aligning with molecular dynamics (MD) simulations showing that electrostatic interactions contribute to the entropic stabilization of the TS complex. These findings suggest that the positive entropy primarily originates from alterations in electrostatic interactions due to the formation of the oxocarbenium ion at C1 in the TS. Differences of electrostatic interactions in ES and TS modify hydrogen bonding of surrounding residues in the active site which causes their side chain dynamics and thus conformational entropy changes. Residues critical for the positive activation entropy are identified. A new BCX mutant with an increased activation entropy and catalytic activity is found.

Comprehending the Efficacy of Whitlock's Caffeine-Pincered Molecular Tweezer on β-Amyloid Aggregation

Alzheimer's disease (AD) stands as one of the most prevalent neurodegenerative conditions, leading to cognitive impairment, with no cure and preventive measures. Misfolding and aberrant aggregation of amyloid-β (Aβ) peptides are believed to be the underlying cause of AD. These amyloid aggregates culminate in the development of toxic Aβ oligomers and subsequent accumulation of β-amyloid plaques amidst neuronal cells in the brain, marking the hallmarks of AD. Drug development for the potentially curative treatment of Alzheimer's is, therefore, a tremendous challenge for the scientific community. In this study, we investigate the potency of Whitlock's caffeine-armed molecular tweezer in combating the deleterious effects of Aβ aggregation, with special emphasis on the seven residue Aβ16-22 fragment. Extensive all-atom molecular dynamics simulations are conducted to probe the various structural and conformational transitions of the peptides in an aqueous medium in both the presence and absence of tweezers. To explore the specifics of peptide-tweezer interactions, radial distribution functions, contact number calculations, binding free energies, and 2-D kernel density plots depicting the variation of distance-angle between the aromatic planes of the peptide-tweezer pair are computed. The central hydrophobic core, particularly the aromatic Phe residues, is crucial in the development of harmful amyloid oligomers. Notably, all analyses indicate reduced interpeptide interactions in the presence of the tweezer, which is attributed to the tweezer-Phe aromatic interaction. Upon increasing the tweezer concentration, the residues of the peptide are further encased in a hydrophobic environment created by the self-aggregating tweezer cluster, leading to the segregation of the peptide residues. This is further aided by the weakening of interstrand hydrogen bonding between the peptides, thereby impeding their self-aggregation and preventing the formation of neurotoxic β-amyloid. Furthermore, the study also highlights the efficacy of the molecular tweezer in destabilizing preformed amyloid fibrils as well as hindering the aggregation of the full-length Aβ1-42 peptide.

Entropy drives the ligand recognition in G-protein-coupled receptor subtypes

Achieving ligand subtype selectivity within highly homologous subtypes of G-protein-coupled receptor (GPCR) is critical yet challenging for GPCR drug discovery, primarily due to the unclear mechanism underlying ligand subtype selectivity, which hampers the rational design of subtype-selective ligands. Herein, we disclose an unusual molecular mechanism of entropy-driven ligand recognition in cannabinoid (CB) receptor subtypes, revealed through atomic-level molecular dynamics simulations, cryoelectron microscopy structure, and mutagenesis experiments. This mechanism is attributed to the distinct conformational dynamics of the receptor's orthosteric pocket, leading to variations in ligand binding entropy and consequently, differential binding affinities, which culminate in specific ligand recognition. We experimentally validated this mechanism and leveraged it to design ligands with enhanced or ablated subtype selectivity. One such ligand demonstrated favorable pharmacokinetic properties and significant efficacy in rodent inflammatory analgesic models. More importantly, it is precisely due to the high subtype selectivity obtained based on this mechanism that this ligand does not show addictive properties in animal models. Our findings elucidate the unconventional role of entropy in CB receptor subtype selectivity and suggest a strategy for rational design of ligands to achieve entropy-driven subtype selectivity for many pharmaceutically important GPCRs.

The αC-β4 loop controls the allosteric cooperativity between nucleotide and substrate in the catalytic subunit of protein kinase A

Allosteric cooperativity between ATP and substrates is a prominent characteristic of the cAMP-dependent catalytic subunit of protein kinase A (PKA-C). This long-range synergistic action is involved in substrate recognition and fidelity, and it may also regulate PKA's association with regulatory subunits and other binding partners. To date, a complete understanding of this intramolecular mechanism is still lacking. Here, we integrated NMR(Nuclear Magnetic Resonance)-restrained molecular dynamics simulations and a Markov State Model to characterize the free energy landscape and conformational transitions of PKA-C. We found that the apoenzyme populates a broad free energy basin featuring a conformational ensemble of the active state of PKA-C (ground state) and other basins with lower populations (excited states). The first excited state corresponds to a previously characterized inactive state of PKA-C with the αC helix swinging outward. The second excited state displays a disrupted hydrophobic packing around the regulatory (R) spine, with a flipped configuration of the F100 and F102 residues at the αC-β4 loop. We validated the second excited state by analyzing the F100A mutant of PKA-C, assessing its structural response to ATP and substrate binding. While PKA-CF100A preserves its catalytic efficiency with Kemptide, this mutation rearranges the αC-β4 loop conformation, interrupting the coupling of the two lobes and abolishing the allosteric binding cooperativity. The highly conserved αC-β4 loop emerges as a pivotal element to control the synergistic binding of nucleotide and substrate, explaining how mutations or insertions near or within this motif affect the function and drug sensitivity in homologous kinases.

Insights into Ligand-Mediated Activation of an Oligomeric Ring-Shaped Gene-Regulatory Protein from Solution- and Solid-State NMR

The 91 kDa oligomeric ring-shaped ligand binding protein TRAP (trp RNA binding attenuation protein) regulates the expression of a series of genes involved in tryptophan (Trp) biosynthesis in bacilli. When cellular Trp levels rise, the free amino acid binds to sites buried in the interfaces between each of the 11 (or 12, depending on the species) protomers in the ring. Crystal structures of Trp-bound TRAP show the Trp ligands are sequestered from solvent by a pair of loops from adjacent protomers that bury the bound ligand via polar contacts to several threonine residues. Binding of the Trp ligands occurs cooperatively, such that successive binding events occur with higher apparent affinity but the structural basis for this cooperativity is poorly understood. We used solution methyl-TROSY NMR relaxation experiments focused on threonine and isoleucine sidechains, as well as magic angle spinning solid-state NMR 13C-13C and 15N-13C chemical shift correlation spectra on uniformly labeled samples recorded at 800 and 1200 MHz, to characterize the structure and dynamics of the protein. Methyl 13C relaxation dispersion experiments on ligand-free apo TRAP revealed concerted exchange dynamics on the μs-ms time scale, consistent with transient sampling of conformations that could allow ligand binding. Cross-correlated relaxation experiments revealed widespread disorder on fast timescales. Chemical shifts for methyl-bearing side chains in apo- and Trp-bound TRAP revealed subtle changes in the distribution of sampled sidechain rotameric states. These observations reveal a pathway and mechanism for induced conformational changes to generate homotropic Trp-Trp binding cooperativity.

Crystal Structures Reveal Hidden Domain Mechanics in Protein Kinase A (PKA)

Cyclic-AMP-dependent protein kinase A (PKA) is a critical enzyme involved in various signaling pathways that plays a crucial role in regulating cellular processes including metabolism, gene transcription, cell proliferation, and differentiation. In this study, the mechanisms of allostery in PKA were investigated by analyzing the vast repertoire of crystal structures available in the RCSB database. From existing structures of murine and human PKA, we elucidated the conformational ensembles and protein dynamics that are altered in a ligand-dependent manner. Distance metrics to analyze conformations of the G-loop were proposed to delineate different states of PKA and were compared to existing structural metrics. Furthermore, ligand-dependent flexibility was investigated through normalized B'-factors to better understand the inherent dynamics in PKA. The presented study provides a contemporary approach to traditional methods in engaging the use of crystal structures for understanding protein dynamics. Importantly, our studies provide a deeper understanding into the conformational ensemble of PKA as the enzyme progresses through its catalytic cycle. These studies provide insights into kinase regulation that can be applied to both PKA individually and protein kinases as a class.

The αC-β4 loop controls the allosteric cooperativity between nucleotide and substrate in the catalytic subunit of protein kinase A

Allosteric cooperativity between ATP and substrates is a prominent characteristic of the cAMP-dependent catalytic (C) subunit of protein kinase A (PKA). Not only this long-range synergistic action is involved in substrate recognition and fidelity, but it is likely to regulate PKA association with regulatory subunits and other binding partners. To date, a complete understanding of the molecular determinants for this intramolecular mechanism is still lacking. Here, we used an integrated NMR-restrained molecular dynamics simulations and a Markov Model to characterize the free energy landscape and conformational transitions of the catalytic subunit of protein kinase A (PKA-C). We found that the apo-enzyme populates a broad free energy basin featuring a conformational ensemble of the active state of PKA-C (ground state) and other basins with lower populations (excited states). The first excited state corresponds to a previously characterized inactive state of PKA-C with the αC helix swinging outward. The second excited state displays a disrupted hydrophobic packing around the regulatory (R) spine, with a flipped configuration of the F100 and F102 residues at the tip of the αC-β4 loop. To experimentally validate the second excited state, we mutated F100 into alanine and used NMR spectroscopy to characterize the binding thermodynamics and structural response of ATP and a prototypical peptide substrate. While the activity of PKA-CF100A toward a prototypical peptide substrate is unaltered and the enzyme retains its affinity for ATP and substrate, this mutation rearranges the αC-β4 loop conformation interrupting the allosteric coupling between nucleotide and substrate. The highly conserved αC-β4 loop emerges as a pivotal element able to modulate the synergistic binding between nucleotide and substrate and may affect PKA signalosome. These results may explain how insertion mutations within this motif affect drug sensitivity in other homologous kinases.

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

Programming conformational cooperativity to regulate allosteric protein-oligonucleotide signal transduction

Conformational cooperativity is a universal molecular effect mechanism and plays a critical role in signaling pathways. However, it remains a challenge to develop artificial molecular networks regulated by conformational cooperativity, due to the difficulties in programming and controlling multiple structural interactions. Herein, we develop a cooperative strategy by programming multiple conformational signals, rather than chemical signals, to regulate protein-oligonucleotide signal transduction, taking advantage of the programmability of allosteric DNA constructs. We generate a cooperative regulation mechanism, by which increasing the loop lengths at two different structural modules induced the opposite effects manifesting as down- and up-regulation. We implement allosteric logic operations by using two different proteins. Further, in cell culture we demonstrate the feasibility of this strategy to cooperatively regulate gene expression of PLK1 to inhibit tumor cell proliferation, responding to orthogonal protein-signal stimulation. This programmable conformational cooperativity paradigm has potential applications in the related fields.

An NMR portrait of functional and dysfunctional allosteric cooperativity in cAMP-dependent protein kinase A

The cAMP-dependent protein kinase A (PKA) is the archetypical eukaryotic kinase. The catalytic subunit (PKA-C) structure is highly conserved among the AGC-kinase family. PKA-C is a bilobal enzyme with a dynamic N-lobe, harbouring the Adenosine-5'-triphosphate (ATP) binding site and a more rigid helical C-lobe. The substrate-binding groove resides at the interface of the two lobes. A distinct feature of PKA-C is the positive binding cooperativity between nucleotide and substrate. Several PKA-C mutations lead to the development of adenocarcinomas, myxomas, and other rare forms of liver tumours. Nuclear magnetic resonance (NMR) spectroscopy shows that these mutations disrupt the allosteric communication between the two lobes, causing a drastic decrease in binding cooperativity. The loss of cooperativity correlates with changes in substrate fidelity and reduced kinase affinity for the endogenous protein kinase inhibitor (PKI). The similarity between PKI and the inhibitory sequence of the kinase regulatory subunits suggests that the overall mechanism of regulation of the kinase may be disrupted. We surmise that a reduced or obliterated cooperativity may constitute a common trait for both orthosteric and allosteric mutations of PKA-C that may lead to dysregulation and disease.

Engineering the Active Site Lid Dynamics to Improve the Catalytic Efficiency of Yeast Cytosine Deaminase

Conformational dynamics is important for enzyme catalysis. However, engineering dynamics to achieve a higher catalytic efficiency is still challenging. In this work, we develop a new strategy to improve the activity of yeast cytosine deaminase (yCD) by engineering its conformational dynamics. Specifically, we increase the dynamics of the yCD C-terminal helix, an active site lid that controls the product release. The C-terminal is extended by a dynamical single α-helix (SAH), which improves the product release rate by up to ~8-fold, and the overall catalytic rate kcat by up to ~2-fold. It is also shown that the kcat increase is due to the favorable activation entropy change. The NMR H/D exchange data indicate that the conformational dynamics of the transition state analog complex increases as the helix is extended, elucidating the origin of the enhanced catalytic entropy. This study highlights a novel dynamics engineering strategy that can accelerate the overall catalysis through the entropy-driven mechanism.

Probing conformational landscapes and mechanisms of allosteric communication in the functional states of the ABL kinase domain using multiscale simulations and network-based mutational profiling of allosteric residue potentials

In the current study, multiscale simulation approaches and dynamic network methods are employed to examine the dynamic and energetic details of conformational landscapes and allosteric interactions in the ABL kinase domain that determine the kinase functions. Using a plethora of synergistic computational approaches, we elucidate how conformational transitions between the active and inactive ABL states can employ allosteric regulatory switches to modulate intramolecular communication networks between the ATP site, the substrate binding region, and the allosteric binding pocket. A perturbation-based network approach that implements mutational profiling of allosteric residue propensities and communications in the ABL states is proposed. Consistent with biophysical experiments, the results reveal functionally significant shifts of the allosteric interaction networks in which preferential communication paths between the ATP binding site and substrate regions in the active ABL state become suppressed in the closed inactive ABL form, which in turn features favorable allosteric coupling between the ATP site and the allosteric binding pocket. By integrating the results of atomistic simulations with dimensionality reduction methods and Markov state models, we analyze the mechanistic role of macrostates and characterize kinetic transitions between the ABL conformational states. Using network-based mutational scanning of allosteric residue propensities, this study provides a comprehensive computational analysis of long-range communications in the ABL kinase domain and identifies conserved regulatory hotspots that modulate kinase activity and allosteric crosstalk between the allosteric pocket, ATP binding site, and substrate binding regions.

Detecting the Hydrogen Bond Cooperativity in a Protein β-Sheet by H/D Exchange

The hydrogen bond (H-bond) cooperativity in the β-sheet of GB3 is investigated by a NMR hydrogen/deuterium (H/D) exchange method. It is shown that the weakening of one backbone N-H…O=C H-bond between two β-strands, β1 and β2, due to the exchange of NH to ND of the H-bond donor in β1, perturbs the chemical shift of ¹³C_α, ¹³C_β, ¹H_α, ¹H_N, and ¹⁵N of the H-bond acceptor and its following residue in β2. Quantum mechanical calculations suggest that the -H-bond chemical shift isotope effect is caused by the structural reorganization in response to the H-bond weakening. This structural reorganization perturbs four neighboring H-bonds, with three being weaker and one being stronger, indicating that three H-bonds are cooperative and one is anticooperative with the perturbed H-bond. The sign of the cooperativity depends on the relative position of the H-bonds. This H-bond cooperativity, which contributes to β-sheet stability overall, can be important for conformational coupling across the β-sheet.

M, Walker C, Kim J, Camilloni C, De Simone A, Vendruscolo M, Bernlohr DA, Taylor SS, Veglia G. ATP-competitive inhibitors modulate the substrate binding cooperativity of a kinase by altering its conformational entropy

ATP-competitive inhibitors are currently the largest class of clinically approved drugs for protein kinases. By targeting the ATP-binding pocket, these compounds block the catalytic activity, preventing substrate phosphorylation. A problem with these drugs, however, is that inhibited kinases may still recognize and bind downstream substrates, acting as scaffolds or binding hubs for signaling partners. Here, using protein kinase A as a model system, we show that chemically different ATP-competitive inhibitors modulate the substrate binding cooperativity by tuning the conformational entropy of the kinase and shifting the populations of its conformationally excited states. Since we found that binding cooperativity and conformational entropy of the enzyme are correlated, we propose a new paradigm for the discovery of ATP-competitive inhibitors, which is based on their ability to modulate the allosteric coupling between nucleotide and substrate-binding sites.

Deciphering the Role of ATP on PHF6 Aggregation

The aggregation of Tau protein, which are involved in Alzheimer's disease, are associated with the self-assembly of the hexapeptide sequence, paired helical filament 6 (PHF6) from repeat 3 of Tau. In order to treat Alzheimer's disease and other such tauopathies, one of the therapeutic strategies is to inhibit aggregation of Tau and its nucleating segments. Therefore, we have studied the effect of adenosine triphosphate (ATP) on the aggregation of PHF6. ATP has, interestingly, demonstrated its ability to inhibit and dissolve protein aggregates. Using classical molecular dynamics simulations, we observed that the hydrophobic core of PHF6 segment displays extended β-sheet conformation, which stabilizes PHF6 aggregates. However, the distribution of ATP around the vicinity of the peptides enables PHF6 to remain discrete and attain random coil conformers. The interpeptide interactions are substituted by PHF6-ATP interactions through hydrogen bonding and hydrophobic interactions (including π-π stacking). Furthermore, the adenosine moiety of ATP contributes more than the triphosphate chain toward PHF6-ATP interaction. Ultimately, this work establishes the inhibitory activity of ATP against Tau aggregation; hence, the therapeutic effect of ATP should be explored further in regard to the effective treatment of Alzheimer's disease.

Raf Kinase Inhibitory Protein regulates the cAMP-dependent protein kinase signaling pathway through a positive feedback loop

Raf Kinase Inhibitory Protein (RKIP) maintains cellular robustness and prevents the progression of diseases such as cancer and heart disease by regulating key kinase cascades including MAP kinase and protein kinase A (PKA). Phosphorylation of RKIP at S153 by Protein Kinase C (PKC) triggers a switch from inhibition of Raf to inhibition of the G protein coupled receptor kinase 2 (GRK2), enhancing signaling by the β-adrenergic receptor (β-AR) that activates PKA. Here we report that PKA-phosphorylated RKIP promotes β-AR-activated PKA signaling. Using biochemical, genetic, and biophysical approaches, we show that PKA phosphorylates RKIP at S51, increasing S153 phosphorylation by PKC and thereby triggering feedback activation of PKA. The S51V mutation blocks the ability of RKIP to activate PKA in prostate cancer cells and to induce contraction in primary cardiac myocytes in response to the β-AR activator isoproterenol, illustrating the functional importance of this positive feedback circuit. As previously shown for other kinases, phosphorylation of RKIP at S51 by PKA is enhanced upon RKIP destabilization by the P74L mutation. These results suggest that PKA phosphorylation at S51 may lead to allosteric changes associated with a higher-energy RKIP state that potentiates phosphorylation of RKIP at other key sites. This allosteric regulatory mechanism may have therapeutic potential for regulating PKA signaling in disease states.

Carbohydrate-Binding Modules of Potential Resources: Occurrence in Nature, Function, and Application in Fiber Recognition and Treatment

Great interests have recently been aroused in the independent associative domain of glycoside hydrolases that utilize insoluble polysaccharides-carbohydrate-binding module (CBM), which responds to binding while the catalytic domain reacts with the substrate. In this mini-review, we first provide a brief introduction on CBM and its subtypes including the classifications, potential sources, structures, and functions. Afterward, the applications of CBMs in substrate recognition based on different types of CBMs have been reviewed. Additionally, the progress of CBMs in paper industry as a new type of environmentally friendly auxiliary agent for fiber treatment is summarized. At last, other applications of CBMs and the future outlook have prospected. Due to the specificity in substrate recognition and diversity in structures, CBM can be a prosperous and promising 'tool' for wood and fiber processing in the future.

Allosteric pluripotency: challenges and opportunities

Allosteric pluripotency arises when the functional response of an allosteric receptor to an allosteric stimulus depends on additional allosteric modulators. Here, we discuss allosteric pluripotency as observed in the prototypical Protein Kinase A (PKA) as well as in other signaling systems, from typical multidomain signaling proteins to bacterial enzymes. We identify key drivers of pluripotent allostery and illustrate how hypothesizing allosteric pluripotency may solve apparent discrepancies currently present in the literature regarding the dual nature of known allosteric modulators. We also outline the implications of allosteric pluripotency for cellular signaling and allosteric drug design, and analyze the challenges and opportunities opened by the pluripotent nature of allostery.

Ligand binding remodels protein side-chain conformational heterogeneity

While protein conformational heterogeneity plays an important role in many aspects of biological function, including ligand binding, its impact has been difficult to quantify. Macromolecular X-ray diffraction is commonly interpreted with a static structure, but it can provide information on both the anharmonic and harmonic contributions to conformational heterogeneity. Here, through multiconformer modeling of time- and space-averaged electron density, we measure conformational heterogeneity of 743 stringently matched pairs of crystallographic datasets that reflect unbound/apo and ligand-bound/holo states. When comparing the conformational heterogeneity of side chains, we observe that when binding site residues become more rigid upon ligand binding, distant residues tend to become more flexible, especially in non-solvent-exposed regions. Among ligand properties, we observe increased protein flexibility as the number of hydrogen bonds decreases and relative hydrophobicity increases. Across a series of 13 inhibitor-bound structures of CDK2, we find that conformational heterogeneity is correlated with inhibitor features and identify how conformational changes propagate differences in conformational heterogeneity away from the binding site. Collectively, our findings agree with models emerging from nuclear magnetic resonance studies suggesting that residual side-chain entropy can modulate affinity and point to the need to integrate both static conformational changes and conformational heterogeneity in models of ligand binding.

Molecular mechanism of glycolytic flux control intrinsic to human phosphoglycerate kinase

Glycolysis plays a fundamental role in energy production and metabolic homeostasis. The intracellular [adenosine triphosphate]/[adenosine diphosphate] ([ATP]/[ADP]) ratio controls glycolytic flux; however, the regulatory mechanism underlying reactions catalyzed by individual glycolytic enzymes enabling flux adaptation remains incompletely understood. Phosphoglycerate kinase (PGK) catalyzes the reversible phosphotransfer reaction, which directly produces ATP in a near-equilibrium step of glycolysis. Despite extensive studies on the transcriptional regulation of PGK expression, the mechanism in response to changes in the [ATP]/[ADP] ratio remains obscure. Here, we report a protein-level regulation of human PGK (hPGK) by utilizing the switching ligand-binding cooperativities between adenine nucleotides and 3-phosphoglycerate (3PG). This was revealed by nuclear magnetic resonance (NMR) spectroscopy at physiological salt concentrations. MgADP and 3PG bind to hPGK with negative cooperativity, whereas MgAMPPNP (a nonhydrolyzable ATP analog) and 3PG bind to hPGK with positive cooperativity. These opposite cooperativities enable a shift between different ligand-bound states depending on the intracellular [ATP]/[ADP] ratio. Based on these findings, we present an atomic-scale description of the reaction scheme for hPGK under physiological conditions. Our results indicate that hPGK intrinsically modulates its function via ligand-binding cooperativities that are finely tuned to respond to changes in the [ATP]/[ADP] ratio. The alteration of ligand-binding cooperativities could be one of the self-regulatory mechanisms for enzymes in bidirectional pathways, which enables rapid adaptation to changes in the intracellular environment.

Enzymatic Cycle-Inspired Dynamic Biosensors Affording No False-Positive Identification

There is an urgent need for reliable biosensors to detect nucleic acid of interest in clinical samples. We propose that the accuracy of the present nucleic acid-sensing method can be advanced by avoiding false-positive identifications derived from nonspecific interactions (e.g., nonspecific binding, probe degradation). The challenge is to exploit biosensors that can distinguish false-positive from true-positive samples in nucleic acid screening. In the present study, by learning from the enzymatic cycle in nature, we raise an allostery tool displaying invertible positive/negative cooperativity for reversible or cyclic activity control of the biosensing probe. We demonstrate that the silencing and regeneration of a positive (or negative) allosteric effector can be carried out through toehold displacement or an enzymatic reaction. We, thus, have developed several dynamic biosensors that can repeatedly measure a single nucleic acid sample. The ability to distinguish a false-positive from a true-positive signal is ascribed to the nonspecific interaction presenting equivalent signal variations, while the specific target binding exhibits diverse signal variations according to repeated measurements. Given its precise identification, such consequent dynamic biosensors offer exciting opportunities in physiological and pathological diagnosis.

The N-terminus of MIF regulates the dynamic profile of residues involved in CD74 activation

Macrophage migration inhibitory factor (MIF) is an immunomodulatory protein with a pathogenic activity in various inflammatory disorders, autoimmune diseases, and cancer. The majority of MIF-triggered pathological conditions are associated with activation of the cell surface receptor CD74. In the absence of small molecule antagonists that directly target CD74, MIF variants and MIF-ligand complexes have served as modulators of CD74 activity. These molecules have been reported to have either antagonistic or agonistic effects against the receptor, although the mechanistic parameters that distinguish the two groups are largely unknown. Through molecular dynamics simulations and NMR experiments, we explored the relationship between MIF's catalytically active N-terminus and the surface residues important for the activation of CD74. We found that the two sites are connected via backbone dynamics that are propagated to the CD74 activation surface of MIF, from the β2 and β4 strands. Our results also provide mechanistic evidence that explain the functional characteristics of MIF variants, serving as CD74 agonists or antagonists. Such findings are of high importance for understanding the MIF-induced activation of CD74 as well as for the development of highly potent CD74 therapeutics.

Dynamic equilibria in protein kinases

Structural changes involved in protein kinase activation and ligand binding have been determined from a wealth of X-ray crystallographic evidence. Recent solution studies using NMR, EPR, HX-MS, and fluorescence techniques have deepened this understanding by highlighting the underlying energetics and dynamics of multistate conformational ensembles. This new research is showing how activation mechanisms and ligand binding alter the internal motions of kinases and enable allosteric coupling between distal regulatory regions and the active site.

Allosteric Type and Pathways Are Governed by the Forces of Protein-Ligand Binding

Allostery is central to many cellular processes, by up- or down-regulating target function. However, what determines the allosteric type remains elusive and currently it is impossible to predict whether the allosteric compounds would activate or inhibit target function before experimental studies. We demonstrated that the allosteric type and allosteric pathways are governed by the forces imposed by ligand binding to target protein using the anisotropic network model and developed an allosteric type prediction method (AlloType). AlloType correctly predicted 13 of the 16 allosteric systems in the data set with experimentally determined protein and complex structures as well as verified allosteric types, which was also used to identify allosteric pathways. When applied to glutathione peroxidase 4, a protein with no complex structure information, AlloType could still be able to predict the allosteric type of the recently reported allosteric activators, demonstrating its potential application in designing specific allosteric drugs and uncovering allosteric mechanisms.

CHESPA/CHESCA-SPARKY: automated NMR data analysis plugins for SPARKY to map protein allostery

MOTIVATION

Correlated Nuclear Magnetic Resonance (NMR) chemical shift changes identified through the CHEmical Shift Projection Analysis (CHESPA) and CHEmical Shift Covariance Analysis (CHESCA) reveal pathways of allosteric transitions in biological macromolecules. To address the need for an automated platform that implements CHESPA and CHESCA and integrates them with other NMR analysis software packages, we introduce here integrated plugins for NMRFAM-SPARKY that implement the seamless detection and visualization of allosteric networks.

AVAILABILITY AND IMPLEMENTATION

CHESCA-SPARKY and CHESPA-SPARKY are available in the latest version of NMRFAM-SPARKY from the National Magnetic Resonance Facility at Madison (http://pine.nmrfam.wisc.edu/download_packages.html), the NMRbox Project (https://nmrbox.org) and to subscribers to the SBGrid (https://sbgrid.org). The assigned spectra involved in this study and tutorial videos using this dataset are available at https://sites.google.com/view/chescachespa-sparky.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics Online.

Enthalpy-Entropy Compensation in Biomolecular Recognition: A Computational Perspective

This mini-review provides an overview of the enthalpy-entropy compensation phenomenon in the simulation of biomacromolecular recognition, with particular emphasis on ligand binding. We approach this complex phenomenon from the point of view of practical computational chemistry. Without providing a detailed description of the plethora of existing methodologies already reviewed in depth elsewhere, we present a series of examples to illustrate different approaches to interpret and predict compensation phenomena at an atomistic level, which is far from trivial to predict using canonical, classic textbook assumptions.

Potential of ATP toward Prevention of hIAPP Oligomerization and Destabilization of hIAPP Protofibrils: An In Silico Perspective

The aggregation of an intrinsically disordered protein, human islet amyloid polypeptide (hIAPP), leads to one of the most prevalent endocrine disorders, type II diabetes mellitus (T2DM). Hence inhibition of hIAPP aggregation provides a possible therapeutic approach for the treatment of T2DM. In this regard, a new aspect of adenosine triphosphate (ATP), which is widely known as the energy source for biological reactions, has recently been discovered, where it can inhibit the formation of protein aggregates and simultaneously dissolve preformed aggregates at a millimolar concentration scale. In this work, we investigate the effect of ATP on the aggregation of an amyloidogenic segment of hIAPP, hIAPP_22-28, and also of the full length sequence. Using all-atom classical molecular dynamics simulations, we observe that the tendency of hIAPP to oligomerize into β-sheet conformers is inhibited by ATP, due to which the peptides remain distant, loosely packed random monomers. Moreover, it can also disassemble preformed hIAPP protofibrils. ATP preferentially interacts with the hydrophobic residues of hIAPP_22-28 fragment and the terminal and turn residues of the full length peptide. The hydrogen bonding, hydrophobic, π-π, and N-H-π stacking interactions are the driving forces for the ATP induced inhibition of hIAPP aggregation. Interestingly, the hydrophobic adenosine of ATP is found to be more in contact with the peptide residues than the hydrophilic triphosphate moiety. The insight into the inhibitory mechanism of ATP on hIAPP aggregation can prove to be beneficial for the design of novel amyloid inhibitors in the future.

Activated nanoscale actin-binding domain motion in the catenin-cadherin complex revealed by neutron spin echo spectroscopy

As the core component of the adherens junction in cell-cell adhesion, the cadherin-catenin complex transduces mechanical tension between neighboring cells. Structural studies have shown that the cadherin-catenin complex exists as an ensemble of flexible conformations, with the actin-binding domain (ABD) of α-catenin adopting a variety of configurations. Here, we have determined the nanoscale protein domain dynamics of the cadherin-catenin complex using neutron spin echo spectroscopy (NSE), selective deuteration, and theoretical physics analyses. NSE reveals that, in the cadherin-catenin complex, the motion of the entire ABD becomes activated on nanosecond to submicrosecond timescales. By contrast, in the α-catenin homodimer, only the smaller disordered C-terminal tail of ABD is moving. Molecular dynamics (MD) simulations also show increased mobility of ABD in the cadherin-catenin complex, compared to the α-catenin homodimer. Biased MD simulations further reveal that the applied external forces promote the transition of ABD in the cadherin-catenin complex from an ensemble of diverse conformational states to specific states that resemble the actin-bound structure. The activated motion and an ensemble of flexible configurations of the mechanosensory ABD suggest the formation of an entropic trap in the cadherin-catenin complex, serving as negative allosteric regulation that impedes the complex from binding to actin under zero force. Mechanical tension facilitates the reduction in dynamics and narrows the conformational ensemble of ABD to specific configurations that are well suited to bind F-actin. Our results provide a protein dynamics and entropic explanation for the observed force-sensitive binding behavior of a mechanosensitive protein complex.

Defective internal allosteric network imparts dysfunctional ATP/substrate-binding cooperativity in oncogenic chimera of protein kinase A

An aberrant fusion of the DNAJB1 and PRKACA genes generates a chimeric protein kinase (PKA-CDNAJB1) in which the J-domain of the heat shock protein 40 is fused to the catalytic α subunit of cAMP-dependent protein kinase A (PKA-C). Deceivingly, this chimeric construct appears to be fully functional, as it phosphorylates canonical substrates, forms holoenzymes, responds to cAMP activation, and recognizes the endogenous inhibitor PKI. Nonetheless, PKA-CDNAJB1 has been recognized as the primary driver of fibrolamellar hepatocellular carcinoma and is implicated in other neoplasms for which the molecular mechanisms remain elusive. Here we determined the chimera's allosteric response to nucleotide and pseudo-substrate binding. We found that the fusion of the dynamic J-domain to PKA-C disrupts the internal allosteric network, causing dramatic attenuation of the nucleotide/PKI binding cooperativity. Our findings suggest that the reduced allosteric cooperativity exhibited by PKA-CDNAJB1 alters specific recognitions and interactions between substrates and regulatory partners contributing to dysregulation.

Allosteric Mechanisms of Nonadditive Substituent Contributions to Protein-Ligand Binding

Quantifying chemical substituent contributions to ligand-binding free energies is challenging due to nonadditive effects. Protein allostery is a frequent cause of nonadditivity, but the underlying allosteric mechanisms often remain elusive. Here, we propose a general NMR-based approach to elucidate such mechanisms and we apply it to the HCN4 ion channel, whose cAMP-binding domain is an archetypal conformational switch. Using NMR, we show that nonadditivity arises not only from concerted conformational transitions, but also from conformer-specific effects, such as steric frustration. Our results explain how affinity-reducing functional groups may lead to affinity gains if combined. Surprisingly, our approach also reveals that nonadditivity depends markedly on the receptor conformation. It is negligible for the inhibited state but highly significant for the active state, opening new opportunities to tune potency and agonism of allosteric effectors.

Programming Biomimetically Confined Aptamers with DNA Frameworks

Active sites of proteins are generally encapsulated within three-dimensional peptide scaffolds that provide the molecular-scale confinement microenvironment. Nevertheless, the ability to tune thermodynamic stability in biomimetic molecular confinement relies on the macromolecular crowding effect of lack of stoichiometry and reconfigurability. Here, we report a framework nucleic acid (FNA)-based strategy to increase thermodynamic stability of aptamers. We demonstrate that the molecular-scale confinement increases the thermodynamic stability of aptamers via facilitated folding kinetics, which is confirmed by the single-molecule FRET (smFRET). Unfavorable conformations of aptamers are restricted as revealed by the Monte Carlo simulation. The binding affinity of the DNA framework-confined aptamer is improved by ∼3-fold. With a similar strategy we improve the catalytic activity of hemin-binding aptamer. Our approach thus shows high potential for designing protein-mimicking DNA nanostructures with enhanced binding affinity and catalytic activity for biosensing and biomedical engineering.

Methyl TROSY spectroscopy: A versatile NMR approach to study challenging biological systems

A major goal in structural biology is to unravel how molecular machines function in detail. To that end, solution-state NMR spectroscopy is ideally suited as it is able to study biological assemblies in a near natural environment. Based on methyl TROSY methods, it is now possible to record high-quality data on complexes that are far over 100 kDa in molecular weight. In this review, we discuss the theoretical background of methyl TROSY spectroscopy, the information that can be extracted from methyl TROSY spectra and approaches that can be used to assign methyl resonances in large complexes. In addition, we touch upon insights that have been obtained for a number of challenging biological systems, including the 20S proteasome, the RNA exosome, molecular chaperones and G-protein-coupled receptors. We anticipate that methyl TROSY methods will be increasingly important in modern structural biology approaches, where information regarding static structures is complemented with insights into conformational changes and dynamic intermolecular interactions.

ATP Controls the Aggregation of Aβ16-22 Peptides

The oligomerization of Aβ_16-22 peptide, which is the hydrophobic core region of full-length Aβ_1-42, causes Alzheimer's disease (AD). This progressive neurodegenerative disease affects over 44 million people worldwide. However, very few synthesized drug molecules are available to inhibit the aggregation of Aβ. Recently, experimental studies have shown that the biological ATP molecule prevents Aβ fibrillation at the millimolar scale; however, the significance of ATP molecules on Aβ fibrillation and the mechanism behind it remain elusive. We have carried out a total of 7.5 μs extensive all-atom molecular dynamics and 8.82 μs of umbrella sampling in explicit water using AMBER14SB, AMBER99SB-ILDN, and AMBER-FB15 force fields for Aβ_16-22 peptide, to investigate the role of ATP on the disruption of Aβ_16-22 prefibrils. From various analyses, such as secondary structure analysis, residue-wise contact map, SASA, and interaction energies, we have observed that, in the presence of ATP, the aggregation of Aβ_16-22 peptide is very unfavorable. Moreover, the biological molecule ATP interacts with the Aβ_16-22 peptide via hydrogen bonding, π-π stacking, and NH-π interactions which, ultimately, prevent the aggregation of Aβ_16-22 peptide. Hence, we assume that the deficiency of ATP may cause Alzheimer's disease (AD).

Characterization of the dynamics and the conformational entropy in the binding between TAZ1 and CTAD-HIF-1α

The interaction between the C-terminal transactivation domain of HIF-1α (CTAD-HIF-1α) and the transcriptional adapter zinc binding 1 (TAZ1) domain of CREB binding protein participate in the initiation of gene transcription during hypoxia. Unbound CTAD-HIF-1α is disordered but undergoes a disorder-to-order transition upon binding to TAZ1. We have here performed NMR side chain and backbone relaxation studies on TAZ1 and side chain relaxation measurements on CTAD-HIF-1α in order to investigate the role of picosecond to nanosecond dynamics. We find that the internal motions are significantly affected upon binding, both on the side chain and the backbone level. The dynamic response corresponds to a conformational entropy change that contributes substantially to the binding thermodynamics for both binding partners. Furthermore, the conformational entropy change for the well-folded TAZ1 varies upon binding to different IDP targets. We further identify a cluster consisting of side chains in bound TAZ1 and CTAD-HIF-1α that experience extensive dynamics and are part of the binding region that involves the N-terminal end of the LPQL motif in CTAD-HIF-1α; a feature that might have an important role in the termination of the hypoxic response.

Dynamic allostery-based molecular workings of kinase:peptide complexes

A dense interplay between structure and dynamics underlies the working of proteins, especially enzymes. Protein kinases are molecular switches that are optimized for their regulation rather than catalytic turnover rates. Using long-simulations dynamic allostery analysis, this study describes an exploration of the dynamic kinase:peptide complex. We have used protein kinase A (PKA) as a model system as a generic prototype of the protein kinase superfamily of signaling enzymes. Our results explain the role of dynamic coupling of active-site residues that must work in coherence to provide for a successful activation or inhibition response from the kinase. Amino acid networks-based community analysis allows us to ponder the conformational entropy of the kinase:nucleotide:peptide ternary complex. We use a combination of 7 peptides that include 3 types of PKA-binding partners: Substrates, products, and inhibitors. The substrate peptides provide for dynamic insights into the enzyme:substrate complex, while the product phospho-peptide allows for accessing modes of enzyme:product release. Mapping of allosteric communities onto the PKA structure allows us to locate the more unvarying and flexible dynamic regions of the kinase. These distributions, when correlated with the structural elements of the kinase core, allow for a detailed exploration of key dynamics-based signatures that could affect peptide recognition and binding at the kinase active site. These studies provide a unique dynamic allostery-based perspective to kinase:peptide complexes that have previously been explored only in a structural or thermodynamic context.

GTP Hydrolysis Without an Active Site Base: A Unifying Mechanism for Ras and Related GTPases

GTP hydrolysis is a biologically crucial reaction, being involved in regulating almost all cellular processes. As a result, the enzymes that catalyze this reaction are among the most important drug targets. Despite their vital importance and decades of substantial research effort, the fundamental mechanism of enzyme-catalyzed GTP hydrolysis by GTPases remains highly controversial. Specifically, how do these regulatory proteins hydrolyze GTP without an obvious general base in the active site to activate the water molecule for nucleophilic attack? To answer this question, we perform empirical valence bond simulations of GTPase-catalyzed GTP hydrolysis, comparing solvent- and substrate-assisted pathways in three distinct GTPases, Ras, Rab, and the G_αi subunit of a heterotrimeric G-protein, both in the presence and in the absence of the corresponding GTPase activating proteins. Our results demonstrate that a general base is not needed in the active site, as the preferred mechanism for GTP hydrolysis is a conserved solvent-assisted pathway. This pathway involves the rate-limiting nucleophilic attack of a water molecule, leading to a short-lived intermediate that tautomerizes to form H₂PO₄^- and GDP as the final products. Our fundamental biochemical insight into the enzymatic regulation of GTP hydrolysis not only resolves a decades-old mechanistic controversy but also has high relevance for drug discovery efforts. That is, revisiting the role of oncogenic mutants with respect to our mechanistic findings would pave the way for a new starting point to discover drugs for (so far) "undruggable" GTPases like Ras.