Linkage between Fitness of Yeast Cells and Adenylate Kinase Catalysis

Protein stability is determined by single-site bias rather than pairwise covariance

The biases revealed in protein sequence alignments have been shown to provide information related to protein structure, stability, and function. For example, sequence biases at individual positions can be used to design consensus proteins that are often more stable than naturally occurring counterparts. Likewise, correlations between pairs of residue can be used to predict protein structures. Recent work using Potts models show that together, single-site biases and pair correlations lead to improved predictions of protein fitness, activity, and stability. Here we use a Potts model to design groups of protein sequences with different amounts of single-site biases and pair correlations, and determine the thermodynamic stabilities of a representative set of sequences from each group. Surprisingly, sequences excluding pair correlations maximize stability, whereas sequences that maximize pair correlations are less stable, suggesting that pair correlations contribute to another aspect of protein fitness. Consistent with this interpretation, we find that for adenylate kinase, enzyme activity is greatly increased by maximizing pair correlations. The finding that elimination of covariant residue pairs increases protein stability suggests a route to enhance stability of designed proteins; indeed, this strategy produces hyperstable homeodomain and adenylate kinase proteins that retain significant activity.

Magnesium induced structural reorganization in the active site of adenylate kinase

Phosphoryl transfer is a fundamental reaction in cellular signaling and metabolism that requires Mg2+ as an essential cofactor. While the primary function of Mg2+ is electrostatic activation of substrates, such as ATP, the full spectrum of catalytic mechanisms exerted by Mg2+ is not known. In this study, we integrate structural biology methods, molecular dynamic (MD) simulations, phylogeny, and enzymology assays to provide molecular insights into Mg2+-dependent structural reorganization in the active site of the metabolic enzyme adenylate kinase. Our results demonstrate that Mg2+ induces a conformational rearrangement of the substrates (ATP and ADP), resulting in a 30° adjustment of the angle essential for reversible phosphoryl transfer, thereby optimizing it for catalysis. MD simulations revealed transitions between conformational substates that link the fluctuation of the angle to large-scale enzyme dynamics. The findings contribute detailed insight into Mg2+ activation of enzymes and may be relevant for reversible and irreversible phosphoryl transfer reactions.

Elucidating Dynamics of Adenylate Kinase from Enzyme Opening to Ligand Release

This study explores ligand-driven conformational changes in adenylate kinase (AK), which is known for its open-to-close conformational transitions upon ligand binding and release. By utilizing string free energy simulations, we determine the free energy profiles for both enzyme opening and ligand release and compare them with profiles from the apoenzyme. Results reveal a three-step ligand release process, which initiates with the opening of the adenosine triphosphate-binding subdomain (ATP lid), followed by ligand release and concomitant opening of the adenosine monophosphate-binding subdomain (AMP lid). The ligands then transition to nonspecific positions before complete dissociation. In these processes, the first step is energetically driven by ATP lid opening, whereas the second step is driven by ATP release. In contrast, the AMP lid opening and its ligand release make minor contributions to the total free energy for enzyme opening. Regarding the ligand binding mechanism, our results suggest that AMP lid closure occurs via an induced-fit mechanism triggered by AMP binding, whereas ATP lid closure follows conformational selection. This difference in the closure mechanisms provides an explanation with implications for the debate on ligand-driven conformational changes of AK. Additionally, we determine an X-ray structure of an AK variant that exhibits significant rearrangements in the stacking of catalytic arginines, explaining its reduced catalytic activity. In the context of apoenzyme opening, the sequence of events is different. Here, the AMP lid opens first while the ATP lid remains closed, and the free energy associated with ATP lid opening varies with orientation, aligning with the reported AK opening and closing rate heterogeneity. Finally, this study, in conjunction with our previous research, provides a comprehensive view of the intricate interplay between various structural elements, ligands, and catalytic residues that collectively contribute to the robust catalytic power of the enzyme.

Mechanistic Basis for a Connection between the Catalytic Step and Slow Opening Dynamics of Adenylate Kinase

Escherichia coli adenylate kinase (AdK) is a small, monomeric enzyme that synchronizes the catalytic step with the enzyme's conformational dynamics to optimize a phosphoryl transfer reaction and the subsequent release of the product. Guided by experimental measurements of low catalytic activity in seven single-point mutation AdK variants (K13Q, R36A, R88A, R123A, R156K, R167A, and D158A), we utilized classical mechanical simulations to probe mutant dynamics linked to product release, and quantum mechanical and molecular mechanical calculations to compute a free energy barrier for the catalytic event. The goal was to establish a mechanistic connection between the two activities. Our calculations of the free energy barriers in AdK variants were in line with those from experiments, and conformational dynamics consistently demonstrated an enhanced tendency toward enzyme opening. This indicates that the catalytic residues in the wild-type AdK serve a dual role in this enzyme's function─one to lower the energy barrier for the phosphoryl transfer reaction and another to delay enzyme opening, maintaining it in a catalytically active, closed conformation for long enough to enable the subsequent chemical step. Our study also discovers that while each catalytic residue individually contributes to facilitating the catalysis, R36, R123, R156, R167, and D158 are organized in a tightly coordinated interaction network and collectively modulate AdK's conformational transitions. Unlike the existing notion of product release being rate-limiting, our results suggest a mechanistic interconnection between the chemical step and the enzyme's conformational dynamics acting as the bottleneck of the catalytic process. Our results also suggest that the enzyme's active site has evolved to optimize the chemical reaction step while slowing down the overall opening dynamics of the enzyme.

Protein dynamics: The future is bright and complicated!

Biological life depends on motion, and this manifests itself in proteins that display motion over a formidable range of time scales spanning from femtoseconds vibrations of atoms at enzymatic transition states, all the way to slow domain motions occurring on micro to milliseconds. An outstanding challenge in contemporary biophysics and structural biology is a quantitative understanding of the linkages among protein structure, dynamics, and function. These linkages are becoming increasingly explorable due to conceptual and methodological advances. In this Perspective article, we will point toward future directions of the field of protein dynamics with an emphasis on enzymes. Research questions in the field are becoming increasingly complex such as the mechanistic understanding of high-order interaction networks in allosteric signal propagation through a protein matrix, or the connection between local and collective motions. In analogy to the solution to the "protein folding problem," we argue that the way forward to understanding these and other important questions lies in the successful integration of experiment and computation, while utilizing the present rapid expansion of sequence and structure space. Looking forward, the future is bright, and we are in a period where we are on the doorstep to, at least in part, comprehend the importance of dynamics for biological function.

Insights into the evolution of enzymatic specificity and catalysis: From Asgard archaea to human adenylate kinases

Enzymatic catalysis is critically dependent on selectivity, active site architecture, and dynamics. To contribute insights into the interplay of these properties, we established an approach with NMR, crystallography, and MD simulations focused on the ubiquitous phosphotransferase adenylate kinase (AK) isolated from Odinarchaeota (OdinAK). Odinarchaeota belongs to the Asgard archaeal phylum that is believed to be the closest known ancestor to eukaryotes. We show that OdinAK is a hyperthermophilic trimer that, contrary to other AK family members, can use all NTPs for its phosphorylation reaction. Crystallographic structures of OdinAK-NTP complexes revealed a universal NTP-binding motif, while 19F NMR experiments uncovered a conserved and rate-limiting dynamic signature. As a consequence of trimerization, the active site of OdinAK was found to be lacking a critical catalytic residue and is therefore considered to be "atypical." On the basis of discovered relationships with human monomeric homologs, our findings are discussed in terms of evolution of enzymatic substrate specificity and cold adaptation.

Proteomic Shifts Reflecting Oxidative Stress and Reduced Capacity for Protein Synthesis, and Alterations to Mitochondrial Membranes in Neurospora crassa Lacking VDAC

Voltage-dependent anion-selective channels (VDAC) maintain the bidirectional flow of small metabolites across the mitochondrial outer membrane and participate in the regulation of multiple cellular processes. To understand the roles of VDAC in cellular homeostasis, preliminary proteomic analyses of S100 cytosolic and mitochondria-enriched fractions from a VDAC-less Neurospora crassa strain (ΔPor-1) were performed. In the variant cells, less abundant proteins include subunits of translation initiation factor eIF-2, enzymes in the shikimate pathway leading to precursors of aromatic amino acids, and enzymes involved in sulfate assimilation and in the synthesis of methionine, cysteine, alanine, serine, and threonine. In contrast, some of the more abundant proteins are involved in electron flow, such as the α subunit of the electron transfer flavoprotein and lactate dehydrogenase, which is involved in one pathway leading to pyruvate synthesis. Increased levels of catalase and catalase activity support predicted increased levels of oxidative stress in ΔPor-1 cells, and higher levels of protein disulfide isomerase suggest activation of the unfolded protein response in the endoplasmic reticulum. ΔPor-1 cells are cold-sensitive, which led us to investigate the impact of the absence of VDAC on several mitochondrial membrane characteristics. Mitochondrial membranes in ΔPor-1 are more fluid than those of wild-type cells, the ratio of C18:1 to C18:3n3 acyl chains is reduced, and ergosterol levels are lower. In summary, these initial results indicate that VDAC-less N. crassa cells are characterized by a lower abundance of proteins involved in amino acid and protein synthesis and by increases in some associated with pyruvate metabolism and stress responses. Membrane lipids and hyphal morphology are also impacted by the absence of VDAC.

Overexpression of Phosphoribosyl Pyrophosphate Synthase Enhances Resistance of Chlamydomonas to Ionizing Radiation

The enzyme phosphoribosyl pyrophosphate synthase (PRPS) catalyzes the conversion of ribose 5-phosphate into phosphoribosyl diphosphate; the latter is a precursor of purine and pyrimidine nucleotides. Here, we investigated the function of PRPS from the single-celled green alga Chlamydomonas reinhardtii in its response to DNA damage from gamma radiation or the alkylating agent LiCl. CrPRPS transcripts were upregulated in cells treated with these agents. We generated CrPRPS-overexpressing transgenic lines to study the function of CrPRPS. When grown in culture with LiCl or exposed to gamma radiation, the transgenic cells grew faster and had a greater survival rate than wild-type cells. CrPRPS overexpression enhanced expression of genes associated with DNA damage response, namely RAD51, RAD1, and LIG1. We observed, from transcriptome analysis, upregulation of genes that code for key enzymes in purine metabolism, namely ribonucleoside-diphosphate pyrophosphokinase subunit M1, adenylate kinase, and nucleoside-diphosphate kinase. We conclude that CrPRPS may affect DNA repair process via regulation of de novo nucleotide synthesis.

Dynamic Connection between Enzymatic Catalysis and Collective Protein Motions

Enzymes employ a wide range of protein motions to achieve efficient catalysis of chemical reactions. While the role of collective protein motions in substrate binding, product release, and regulation of enzymatic activity is generally understood, their roles in catalytic steps per se remain uncertain. Here, molecular dynamics simulations, enzyme kinetics, X-ray crystallography, and nuclear magnetic resonance spectroscopy are combined to elucidate the catalytic mechanism of adenylate kinase and to delineate the roles of catalytic residues in catalysis and the conformational change in the enzyme. This study reveals that the motions in the active site, which occur on a time scale of picoseconds to nanoseconds, link the catalytic reaction to the slow conformational dynamics of the enzyme by modulating the free energy landscapes of subdomain motions. In particular, substantial conformational rearrangement occurs in the active site following the catalytic reaction. This rearrangement not only affects the reaction barrier but also promotes a more open conformation of the enzyme after the reaction, which then results in an accelerated opening of the enzyme compared to that of the reactant state. The results illustrate a linkage between enzymatic catalysis and collective protein motions, whereby the disparate time scales between the two processes are bridged by a cascade of intermediate-scale motion of catalytic residues modulating the free energy landscapes of the catalytic and conformational change processes.

Nucleation of an Activating Conformational Change by a Cation-π Interaction

As a key molecule in biology, adenosine triphosphate (ATP) has numerous crucial functions in, for instance, energetics, post-translational modifications, nucleotide biosynthesis, and cofactor metabolism. Here, we have discovered an intricate interplay between the enzyme adenylate kinase and its substrate ATP. The side chain of an arginine residue was found to be an efficient sensor of the aromatic moiety of ATP through the formation of a strong cation-π interaction. In addition to recognition, the interaction was found to have dual functionality. First, it nucleates the activating conformational transition of the ATP binding domain and also affects the specificity in the distant AMP binding domain. In light of the functional consequences resulting from the cation-π interaction, it is possible that the mode of ATP recognition may be a useful tool in enzyme design.

Cell cycle-dependent association of polo kinase Cdc5 with CENP-A contributes to faithful chromosome segregation in budding yeast

Evolutionarily conserved polo-like kinase, Cdc5 (Plk1 in humans), associates with kinetochores during mitosis; however, the role of cell cycle-dependent centromeric ( CEN) association of Cdc5 and its substrates that exclusively localize to the kinetochore have not been characterized. Here we report that evolutionarily conserved CEN histone H3 variant, Cse4 (CENP-A in humans), is a substrate of Cdc5, and that the cell cycle-regulated association of Cse4 with Cdc5 is required for cell growth. Cdc5 contributes to Cse4 phosphorylation in vivo and interacts with Cse4 in mitotic cells. Mass spectrometry analysis of in vitro kinase assays showed that Cdc5 phosphorylates nine serine residues clustered within the N-terminus of Cse4. Strains with cse4-9SA exhibit increased errors in chromosome segregation, reduced levels of CEN-associated Mif2 and Mcd1/Scc1 when combined with a deletion of MCM21. Moreover, the loss of Cdc5 from the CEN chromatin contributes to defects in kinetochore integrity and reduction in CEN-associated Cse4. The cell cycle-regulated association of Cdc5 with Cse4 is essential for cell viability as constitutive association of Cdc5 with Cse4 at the kinetochore leads to growth defects. In summary, our results have defined a role for Cdc5-mediated Cse4 phosphorylation in faithful chromosome segregation.

Transcriptome analysis provides insights into the molecular mechanisms responsible for evisceration behavior in the sea cucumber Apostichopus japonicus

The sea cucumber Apostichopus japonicus (Selenka) is a valuable economic species in Southeast Asia. It has many fascinating behavioral characteristics, such as autolysis, aestivation, regeneration, and evisceration, thus it is a notable species for studies of special behaviors. Evisceration and autotomy are controlled by the neural network and involve a complicated physiological process. The occurrence of evisceration behavior in sea cucumbers is strongly related to their environment, and it negatively impacts their economic value. Evisceration behavior plays a pivotal role in the survival of A. japonicus, and when it is induced by dramatic changes in the coastal ecological environment and the aquaculture setting it can strongly affect the economic performance of this species. Although numerous studies have focused on intestinal regeneration of A. japonicus, less is known about evisceration behavior, especially its underlying molecular mechanisms. Thus, identification of genes that regulate evisceration in the sea cucumber likely will provide a scientific explanation for this significant specific behavior. In this study, Illumina sequencing (RNA-Seq) was performed on A. japonicus specimens in three states: normal (TCQ), eviscerating (TCZ), and 3 h after evisceration (TCH). In total, 129,905 unigenes were generated with an N50 length of 2651 base pairs, and 54,787 unigenes were annotated from seven functional databases (KEGG, KOG, GO, NR, NT, Interpro, and Swiss-Prot). Additionally, 190, 191, and 320 genes were identified as differentially expressed genes (DEGs) in the comparisons of TCQ vs. TCZ, TCZ vs. TCH, and TCQ vs. TCH, respectively. These DEGs mapped to 157, 113, and 190 signaling pathways in the KEGG database, respectively. KEGG analyses also revealed that potential DEGs enriched in the categories of "environmental information processing," "organismal system," "metabolism," and "cellular processes," and they were involved in evisceration behavior in A. japonicus. These DEGs are related to muscle contraction, hormone and neurotransmitter secretion, nerve and muscle damage, energy support, cellular stress, and apoptosis. In conclusion, through our comparative analysis of A. japonicus in different stages, we identified many candidate evisceration-related genes and signaling pathways that likely are involved in evisceration behavior. These results should help further elucidate the mechanisms underlying evisceration behavior in sea cucumbers.

Molecular mechanism of ATP versus GTP selectivity of adenylate kinase

Enzymatic substrate selectivity is critical for the precise control of metabolic pathways. In cases where chemically related substrates are present inside cells, robust mechanisms of substrate selectivity are required. Here, we report the mechanism utilized for catalytic ATP versus GTP selectivity during adenylate kinase (Adk) -mediated phosphorylation of AMP. Using NMR spectroscopy we found that while Adk adopts a catalytically competent and closed structural state in complex with ATP, the enzyme is arrested in a catalytically inhibited and open state in complex with GTP. X-ray crystallography experiments revealed that the interaction interfaces supporting ATP and GTP recognition, in part, are mediated by coinciding residues. The mechanism provides an atomic view on how the cellular GTP pool is protected from Adk turnover, which is important because GTP has many specialized cellular functions. In further support of this mechanism, a structure-function analysis enabled by synthesis of ATP analogs suggests that a hydrogen bond between the adenine moiety and the backbone of the enzyme is vital for ATP selectivity. The importance of the hydrogen bond for substrate selectivity is likely general given the conservation of its location and orientation across the family of eukaryotic protein kinases.

The Effect of Fermented Lingonberry Juice on Candida glabrata Intracellular Protein Expression

Lingonberries have a long traditional use in treating fungal infections on mucosal membranes, but very little is known about the exact antifungal mechanisms. We tested the effects of fermented lingonberry juice on Candida glabrata intracellular protein expression. A Candida glabrata clinical strain was grown in the presence of fermented lingonberry juice (FLJ). Also the effect of lowered pH was tested. Intracellular protein expression levels were analyzed by the 2D-DIGE method. Six proteins detected with ≥1.5-fold lowered expression levels from FLJ treated cells were further characterized with LC-MS/MS. Heat shock protein 9/12 and redoxin were identified with peptide coverage/scores of 68/129 and 21/26, respectively. Heat shock protein 9/12 had an oxidized methionine at position 56. We found no differences in protein expression levels at pH 3.5 compared to pH 7.6. These results demonstrate that FLJ exerts an intracellular stress response in Candida glabrata, plausibly impairing its ability to express proteins related to oxidative stress or maintaining cell wall integrity.