Effectiveness and limitations of local structural entropy optimization in the thermal stabilization of mesophilic and thermophilic adenylate kinases

Evolutionary-Scale Enzymology Enables Biochemical Constant Prediction Across a Multi-Peaked Catalytic Landscape

Quantitatively mapping enzyme sequence-catalysis landscapes remains a critical challenge in understanding enzyme function, evolution, and design. Here, we expand an emerging microfluidic platform to measure catalytic constants-k cat and K M-for hundreds of diverse naturally occurring sequences and mutants of the model enzyme Adenylate Kinase (ADK). This enables us to dissect the sequence-catalysis landscape's topology, navigability, and mechanistic underpinnings, revealing distinct catalytic peaks organized by structural motifs. These results challenge long-standing hypotheses in enzyme adaptation, demonstrating that thermophilic enzymes are not slower than their mesophilic counterparts. Combining the rich representations of protein sequences provided by deep-learning models with our custom high-throughput kinetic data yields semi-supervised models that significantly outperform existing models at predicting catalytic parameters of naturally occurring ADK sequences. Our work demonstrates a promising strategy for dissecting sequence-catalysis landscapes across enzymatic evolution and building family-specific models capable of accurately predicting catalytic constants, opening new avenues for enzyme engineering and functional prediction.

Structure determination of high-energy states in a dynamic protein ensemble

Macromolecular function frequently requires that proteins change conformation into high-energy states1-4. However, methods for solving the structures of these functionally essential, lowly populated states are lacking. Here we develop a method for high-resolution structure determination of minorly populated states by coupling NMR spectroscopy-derived pseudocontact shifts5 (PCSs) with Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion6 (PCS-CPMG). Our approach additionally defines the corresponding kinetics and thermodynamics of high-energy excursions, thereby characterizing the entire free-energy landscape. Using a large set of simulated data for adenylate kinase (Adk), calmodulin and Src kinase, we find that high-energy PCSs accurately determine high-energy structures (with a root mean squared deviation of less than 3.5 angström). Applying our methodology to Adk during catalysis, we find that the high-energy excursion involves surprisingly small openings of the AMP and ATP lids. This previously unresolved high-energy structure solves a longstanding controversy about conformational interconversions that are rate-limiting for catalysis. Primed for either substrate binding or product release, the high-energy structure of Adk suggests a two-step mechanism combining conformational selection to this state, followed by an induced-fit step into a fully closed state for catalysis of the phosphoryl-transfer reaction. Unlike other methods for resolving high-energy states, such as cryo-electron microscopy and X-ray crystallography, our solution PCS-CPMG approach excels in cases involving domain rearrangements of smaller systems (less than 60 kDa) and populations as low as 0.5%, and enables the simultaneous determination of protein structure, kinetics and thermodynamics while proteins perform their function.

Rational Design of Adenylate Kinase Thermostability through Coevolution and Sequence Divergence Analysis

Protein engineering is actively pursued in industrial and laboratory settings for high thermostability. Among the many protein engineering methods, rational design by bioinformatics provides theoretical guidance without time-consuming experimental screenings. However, most rational design methods either rely on protein tertiary structure information or have limited accuracies. We proposed a primary-sequence-based algorithm for increasing the heat resistance of a protein while maintaining its functions. Using adenylate kinase (ADK) family as a model system, this method identified a series of amino acid sites closely related to thermostability. Single- and double-point mutants constructed based on this method increase the thermal denaturation temperature of the mesophilic Escherichia coli (E. coli) ADK by 5.5 and 8.3 °C, respectively, while preserving most of the catalytic function at ambient temperatures. Additionally, the constructed mutants have improved enzymatic activity at higher temperature.

Structural and mutational analyses of psychrophilic and mesophilic adenylate kinases highlight the role of hydrophobic interactions in protein thermal stability

Protein thermal stability is an important field since thermally stable proteins are desirable in many academic and industrial settings. Information on protein thermal stabilization can be obtained by comparing homologous proteins from organisms living at distinct temperatures. Here, we report structural and mutational analyses of adenylate kinases (AKs) from psychrophilic Bacillus globisporus (AKp) and mesophilic Bacillus subtilis (AKm). Sequence and structural comparison showed suboptimal hydrophobic packing around Thr26 in the CORE domain of AKp, which was replaced with an Ile residue in AKm. Mutations that improved hydrophobicity of the Thr residue increased the thermal stability of the psychrophilic AKp, and the largest stabilization was observed for a Thr-to-Ile substitution. Furthermore, a reverse Ile-to-Thr mutation in the mesophilic AKm significantly decreased thermal stability. We determined the crystal structures of mutant AKs to confirm the impact of the residue substitutions on the overall stability. Taken together, our results provide a structural basis for the stability difference between psychrophilic and mesophilic AK homologues and highlight the role of hydrophobic interactions in protein thermal stability.

Structural analyses of adenylate kinases from Antarctic and tropical fishes for understanding cold adaptation of enzymes

Psychrophiles are extremophilic organisms capable of thriving in cold environments. Proteins from these cold-adapted organisms can remain physiologically functional at low temperatures, but are structurally unstable even at moderate temperatures. Here, we report the crystal structure of adenylate kinase (AK) from the Antarctic fish Notothenia coriiceps, and identify the structural basis of cold adaptation by comparison with homologues from tropical fishes including Danio rerio. The structure of N. coriiceps AK (AKNc) revealed suboptimal hydrophobic packing around three Val residues in its central CORE domain, which are replaced with Ile residues in D. rerio AK (AKDr). The Val-to-Ile mutations that improve hydrophobic CORE packing in AKNc increased stability at high temperatures but decreased activity at low temperatures, suggesting that the suboptimal hydrophobic CORE packing is important for cold adaptation. Such linkage between stability and activity was also observed in AKDr. Ile-to-Val mutations that destabilized the tropical AK resulted in increased activity at low temperatures. Our results provide the structural basis of cold adaptation of a psychrophilic enzyme from a multicellular, eukaryotic organism, and highlight the similarities and differences in the structural adjustment of vertebrate and bacterial psychrophilic AKs during cold adaptation.