Assessing the Catalytic Role of Native Glucagon Amyloid Fibrils.

Glucagon stands out as a pivotal peptide hormone, instrumental in controlling blood glucose levels and lipid metabolism. While the formation of glucagon amyloid fibrils has been documented, their biological functions remain enigmatic. Recently, we demonstrated experimentally that glucagon amyloid fibrils can act as catalysts in several biological reactions including esterolysis, lipid hydrolysis, and dephosphorylation. Herein, we present a multiscale quantum mechanics/molecular mechanics (QM/MM) simulation of the acylation step in the esterolysis of para-nitrophenyl acetate (p-NPA), catalyzed by native glucagon amyloid fibrils, serving as a model system to elucidate their catalytic function. This step entails a concerted mechanism, involving proton transfer from serine to histidine, followed by the nucleophilic attack of the serine oxy anion on the carbonyl carbon of p-NPA. We computed the binding energy and free-energy profiles of this reaction using the protein-dipole Langevin-dipole (PDLD) within the linear response approximation (LRA) framework (PDLD/S-LRA-2000) and the empirical valence bond (EVB) methods. This included simulations of the reaction in an aqueous environment and in the fibril, enabling us to estimate the catalytic effect of the fibril. Our EVB calculations obtained a barrier of 23.4 kcal mol-1 for the enzyme-catalyzed reaction compared to the experimental value of 21.9 kcal mol-1 (and a calculated catalytic effect of 3.2 kcal mol-1 compared to the observed effect of 4.7 kcal mol-1). This close agreement together with the barrier reduction when transitioning from the reference solution reaction to the amyloid fibril provides supporting evidence to the catalytic role of glucagon amyloid fibrils. Moreover, employing the PDLD/S-LRA-2000 approach further reinforced exclusively the enzyme's catalytic role. The results presented in this study contribute significantly to our understanding of the catalytic role of glucagon amyloid fibrils, marking, to the best of our knowledge, the first-principles mechanistic investigation of fibrils using QM/MM methods. Therefore, our findings offer fruitful insights for future research into the mechanisms of related amyloid catalysis.

Exploring the Light-Emitting Agents in Renilla Luciferases by an Effective QM/MM Approach.

Bioluminescence is a fascinating natural phenomenon, wherein organisms produce light through specific biochemical reactions. Among these organisms, Renilla luciferase (RLuc) derived from the sea pansy Renilla reniformis is notable for its blue light emission and has potential applications in bioluminescent tagging. Our study focuses on RLuc8, a variant of RLuc with eight amino acid substitutions. Recent studies have shown that the luminescent emitter coelenteramide can adopt multiple protonation states, which may be influenced by nearby residues at the enzyme's active site, demonstrating a complex interplay between protein structure and bioluminescence. Herein, using the quantum mechanical consistent force field method and the semimacroscopic protein dipole-Langevin dipole method with linear response approximation, we show that the phenolate state of coelenteramide in RLuc8 is the primary light-emitting species in agreement with experimental results. Our calculations also suggest that the proton transfer (PT) from neutral coelenteramide to Asp162 plays a crucial role in the bioluminescence process. Additionally, we reproduced the observed emission maximum for the amide anion in RLuc8-D120A and the pyrazine anion in the presence of a Na+ counterion in RLuc8-D162A, suggesting that these are the primary emitters. Furthermore, our calculations on the neutral emitter in the engineered AncFT-D160A enzyme, structurally akin to RLuc8-D162A but with a considerably blue-shifted emission peak, aligned with the observed data, possibly explaining the variance in emission peaks. Overall, this study demonstrates an effective approach to investigate chromophores' bimolecular states while incorporating the PT process in emission spectra calculations, contributing valuable insights for future studies of PT in photoproteins.

Enhancing luciferase activity and stability through generative modeling of natural enzyme sequences.

The availability of natural protein sequences synergized with generative AI provides new paradigms to engineer enzymes. Although active enzyme variants with numerous mutations have been designed using generative models, their performance often falls short of their wild type counterparts. Additionally, in practical applications, choosing fewer mutations that can rival the efficacy of extensive sequence alterations is usually more advantageous. Pinpointing beneficial single mutations continues to be a formidable task. In this study, using the generative maximum entropy model to analyze Renilla luciferase (RLuc) homologs, and in conjunction with biochemistry experiments, we demonstrated that natural evolutionary information could be used to predictively improve enzyme activity and stability by engineering the active center and protein scaffold, respectively. The success rate to improve either luciferase activity or stability of designed single mutants is ~50%. This finding highlights nature's ingenious approach to evolving proficient enzymes, wherein diverse evolutionary pressures are preferentially applied to distinct regions of the enzyme, ultimately culminating in an overall high performance. We also reveal an evolutionary preference in RLuc toward emitting blue light that holds advantages in terms of water penetration compared to other light spectra. Taken together, our approach facilitates navigation through enzyme sequence space and offers effective strategies for computer-aided rational enzyme engineering.

Enhancing Luciferase Activity and Stability through Generative Modeling of Natural Enzyme Sequences.

The availability of natural protein sequences synergized with generative artificial intelligence (AI) provides new paradigms to create enzymes. Although active enzyme variants with numerous mutations have been produced using generative models, their performance often falls short compared to their wild-type counterparts. Additionally, in practical applications, choosing fewer mutations that can rival the efficacy of extensive sequence alterations is usually more advantageous. Pinpointing beneficial single mutations continues to be a formidable task. In this study, using the generative maximum entropy model to analyze Renilla luciferase homologs, and in conjunction with biochemistry experiments, we demonstrated that natural evolutionary information could be used to predictively improve enzyme activity and stability by engineering the active center and protein scaffold, respectively. The success rate of designed single mutants is ~50% to improve either luciferase activity or stability. These finding highlights nature's ingenious approach to evolving proficient enzymes, wherein diverse evolutionary pressures are preferentially applied to distinct regions of the enzyme, ultimately culminating in an overall high performance. We also reveal an evolutionary preference in Renilla luciferase towards emitting blue light that holds advantages in terms of water penetration compared to other light spectrum. Taken together, our approach facilitates navigation through enzyme sequence space and offers effective strategies for computer-aided rational enzyme engineering.

Predicting Mutational Effects on Receptor Binding of the Spike Protein of SARS-CoV-2 Variants.

The pandemic caused by SARS-CoV-2 has cost millions of lives and tremendous social/financial loss. The virus continues to evolve and mutate. In particular, the recently emerged "UK", "South Africa", and Delta variants show higher infectivity and spreading speed. Thus, the relationship between the mutations of certain amino acids and the spreading speed of the virus is a problem of great importance. In this respect, understanding the mutational mechanism is crucial for surveillance and prediction of future mutations as well as antibody/vaccine development. In this work, we used a coarse-grained model (that was used previously in predicting the importance of mutations of N501) to calculate the free energy change of various types of single-site or combined-site mutations. This was done for the UK, South Africa, and Delta mutants. We investigated the underlying mechanisms of the binding affinity changes for mutations at different spike protein domains of SARS-CoV-2 and provided the energy basis for the resistance of the E484 mutant to the antibody m396. Other potential mutation sites were also predicted. Furthermore, the in silico predictions were assessed by functional experiments. The results establish that the faster spreading of recently observed mutants is strongly correlated with the binding-affinity enhancement between virus and human receptor as well as with the reduction of the binding to the m396 antibody. Significantly, the current approach offers a way to predict new variants and to assess the effectiveness of different antibodies toward such variants.