Principles of Protein Stability and Their Application in Computational Design

Engineering the next-generation synthetic cell factory driven by protein engineering

Synthetic cell factory offers substantial advantages in economically efficient production of biofuels, chemicals, and pharmaceutical compounds. However, to create a high-performance synthetic cell factory, precise regulation of cellular material and energy flux is essential. In this context, protein components including enzymes, transcription factor-based biosensors and transporters play pivotal roles. Protein engineering aims to create novel protein variants with desired properties by modifying or designing protein sequences. This review focuses on summarizing the latest advancements of protein engineering in optimizing various aspects of synthetic cell factory, including: enhancing enzyme activity to eliminate production bottlenecks, altering enzyme selectivity to steer metabolic pathways towards desired products, modifying enzyme promiscuity to explore innovative routes, and improving the efficiency of transporters. Furthermore, the utilization of protein engineering to modify protein-based biosensors accelerates evolutionary process and optimizes the regulation of metabolic pathways. The remaining challenges and future opportunities in this field are also discussed.

Enhanced Thermostability of Nattokinase by Computation-Based Rational Redesign of Flexible Regions

Nattokinase is a nutrient in healthy food natto that has the function of preventing and treating blood thrombus. However, its low thermostability and fibrinolytic activity limit its application in food and pharmaceuticals. In this study, we used bioinformatics analysis to identify two loops (loop10 and loop12) in the flexible region of nattokinase rAprY. Using this basis, we screened the G131S-S161T variant, which showed a 2.38-fold increase in half-life at 55 °C, and the M3 variant, which showed a 2.01-fold increase in activity, by using a thermostability prediction algorithm. Bioinformatics analysis revealed that the enhanced thermostability of the G131S-S161T variant was due to the increased rigidity and structural shrinkage of the overall structure. Additionally, the increased rigidity of the local region surrounding the active center and its mutated sites helps maintain its normal conformation in high-temperature environments. The increased catalytic activity of the M3 variant may be due to its more efficient substrate binding mechanism. We investigated strategies to improve the thermostability and fibrinolytic activity of nattokinase, and the resulting variants show promise for industrial production and application.

Application of an alchemical free energy method for the prediction of thermostable DuraPETase variants

Non-equilibrium (NEQ) alchemical free energy calculations are an emerging tool for accurately predicting changes in protein folding free energy resulting from amino acid mutations. In this study, this method in combination with the Rosetta ddg monomer tool was applied to predict more thermostable variants of the polyethylene terephthalate (PET) degrading enzyme DuraPETase. The Rosetta ddg monomer tool efficiently enriched promising mutations prior to more accurate prediction by NEQ alchemical free energy calculations. The relative change in folding free energy of 96 single amino acid mutations was calculated by NEQ alchemical free energy calculation. Experimental validation of ten of the highest scoring variants identified two mutations (DuraPETaseS61M and DuraPETaseS223Y) that increased the melting temperature (Tm) of the enzyme by up to 1 °C. The calculated relative change in folding free energy showed an excellent correlation with experimentally determined Tm resulting in a Pearson's correlation coefficient of r =  - 0.84. Limitations in the prediction of strongly stabilizing mutations were, however, encountered and are discussed. Despite these challenges, this study demonstrates the practical applicability of NEQ alchemical free energy calculations in prospective enzyme engineering projects. KEY POINTS: • Rosetta ddg monomer enriches stabilizing mutations in a library of DuraPETase variants • NEQ free energy calculations accurately predict changes in Tm of DuraPETase • The DuraPETase variants S223Y, S42M, and S61M have increased Tm.

Combining Rosetta Sequence Design with Protein Language Model Predictions Using Evolutionary Scale Modeling (ESM) as Restraint

Computational protein sequence design has the ambitious goal of modifying existing or creating new proteins; however, designing stable and functional proteins is challenging without predictability of protein dynamics and allostery. Informing protein design methods with evolutionary information limits the mutational space to more native-like sequences and results in increased stability while maintaining functions. Recently, language models, trained on millions of protein sequences, have shown impressive performance in predicting the effects of mutations. Assessing Rosetta-designed sequences with a language model showed scores that were worse than those of their original sequence. To inform Rosetta design protocols with language model predictions, we added a new metric to restrain the energy function during design using the Evolutionary Scale Modeling (ESM) model. The resulting sequences have better language model scores and similar sequence recovery, with only a minor decrease in the fitness as assessed by Rosetta energy. In conclusion, our work combines the strength of recent machine learning approaches with the Rosetta protein design toolbox.

Opportunities and challenges in design and optimization of protein function

The field of protein design has made remarkable progress over the past decade. Historically, the low reliability of purely structure-based design methods limited their application, but recent strategies that combine structure-based and sequence-based calculations, as well as machine learning tools, have dramatically improved protein engineering and design. In this Review, we discuss how these methods have enabled the design of increasingly complex structures and therapeutically relevant activities. Additionally, protein optimization methods have improved the stability and activity of complex eukaryotic proteins. Thanks to their increased reliability, computational design methods have been applied to improve therapeutics and enzymes for green chemistry and have generated vaccine antigens, antivirals and drug-delivery nano-vehicles. Moreover, the high success of design methods reflects an increased understanding of basic rules that govern the relationships among protein sequence, structure and function. However, de novo design is still limited mostly to α-helix bundles, restricting its potential to generate sophisticated enzymes and diverse protein and small-molecule binders. Designing complex protein structures is a challenging but necessary next step if we are to realize our objective of generating new-to-nature activities.

A Novel Trypsin Kunitz-Type Inhibitor from Cajanus cajan Leaves and Its Inhibitory Activity on New Cancer Serine Proteases and Its Effect on Tumor Cell Growth

A novel trypsin inhibitor from Cajanus cajan (TIC) fresh leaves was partially purified by affinity chromatography. SDS-PAGE revealed one band with about 15 kDa with expressive trypsin inhibitor activity by zymography. TIC showed high affinity for trypsin (Ki = 1.617 μM) and was a competitive inhibitor for this serine protease. TIC activity was maintained after 24 h of treatment at 70 °C, after 1 h treatments with different pH values, and β-mercaptoethanol increasing concentrations, and demonstrated expressive structural stability. However, the activity of TIC was affected in the presence of oxidizing agents. In order to study the effect of TIC on secreted serine proteases, as well as on the cell culture growth curve, SK-MEL-28 metastatic human melanoma cell line and CaCo-2 colon adenocarcinoma was grown in supplemented DMEM, and the extracellular fractions were submitted salting out and affinity chromatography to obtain new secreted serine proteases. TIC inhibited almost completely, 96 to 89%, the activity of these serine proteases and reduced the melanoma and colon adenocarcinoma cells growth of 48 and 77% respectively. Besides, it is the first time that a trypsin inhibitor was isolated and characterized from C. cajan leaves and cancer serine proteases were isolated and partial characterized from SK-MEL-28 and CaCo-2 cancer cell lines. Furthermore, TIC shown to be potent inhibitor of tumor protease affecting cell growth, and can be one potential drug candidate to be employed in chemotherapy of melanoma and colon adenocarcinoma.

CyclicPepedia: a knowledge base of natural and synthetic cyclic peptides

Cyclic peptides offer a range of notable advantages, including potent antibacterial properties, high binding affinity and specificity to target molecules, and minimal toxicity, making them highly promising candidates for drug development. However, a comprehensive database that consolidates both synthetically derived and naturally occurring cyclic peptides is conspicuously absent. To address this void, we introduce CyclicPepedia (https://www.biosino.org/iMAC/cyclicpepedia/), a pioneering database that encompasses 8744 known cyclic peptides. This repository, structured as a composite knowledge network, offers a wealth of information encompassing various aspects of cyclic peptides, such as cyclic peptides' sources, categorizations, structural characteristics, pharmacokinetic profiles, physicochemical properties, patented drug applications, and a collection of crucial publications. Supported by a user-friendly knowledge retrieval system and calculation tools specifically designed for cyclic peptides, CyclicPepedia will be able to facilitate advancements in cyclic peptide drug development.

Stepwise introduction of stabilizing mutations reveals nonlinear additive effects in de novo TIM barrels

Over the past decades, the TIM-barrel fold has served as a model system for the exploration of how changes in protein sequences affect their structural, stability, and functional characteristics, and moreover, how this information can be leveraged to design proteins from the ground up. After numerous attempts to design de novo proteins with this specific fold, sTIM11 was the first validated de novo design of an idealized four-fold symmetric TIM barrel. Subsequent efforts to enhance the stability of this initial design resulted in the development of DeNovoTIMs, a family of de novo TIM barrels with various stabilizing mutations. In this study, we present an investigation into the biophysical and thermodynamic effects upon introducing a varying number of stabilizing mutations per quarter along the sequence of a four-fold symmetric TIM barrel. We compared the base design DeNovoTIM0 without any stabilizing mutations with variants containing mutations in one, two, three, and all four quarters-designated TIM1q, TIM2q, TIM3q, and DeNovoTIM6, respectively. This analysis revealed a stepwise and nonlinear change in the thermodynamic properties that correlated with the number of mutated quarters, suggesting positive nonadditive effects. To shed light on the significance of the location of stabilized quarters, we engineered two variants of TIM2q which contain the same number of mutations but positioned in different quarter locations. Characterization of these TIM2q variants revealed that the mutations exhibit varying effects on the overall protein stability, contingent upon the specific region in which they are introduced. These findings emphasize that the amount and location of stabilized interfaces among the four quarters play a crucial role in shaping the conformational stability of these four-fold symmetric TIM barrels. Analysis of de novo proteins, as described in this study, enhances our understanding of how sequence variations can finely modulate stability in both naturally occurring and computationally designed proteins.

Understanding activity-stability tradeoffs in biocatalysts by enzyme proximity sequencing

Understanding the complex relationships between enzyme sequence, folding stability and catalytic activity is crucial for applications in industry and biomedicine. However, current enzyme assay technologies are limited by an inability to simultaneously resolve both stability and activity phenotypes and to couple these to gene sequences at large scale. Here we present the development of enzyme proximity sequencing, a deep mutational scanning method that leverages peroxidase-mediated radical labeling with single cell fidelity to dissect the effects of thousands of mutations on stability and catalytic activity of oxidoreductase enzymes in a single experiment. We use enzyme proximity sequencing to analyze how 6399 missense mutations influence folding stability and catalytic activity in a D-amino acid oxidase from Rhodotorula gracilis. The resulting datasets demonstrate activity-based constraints that limit folding stability during natural evolution, and identify hotspots distant from the active site as candidates for mutations that improve catalytic activity without sacrificing stability. Enzyme proximity sequencing can be extended to other enzyme classes and provides valuable insights into biophysical principles governing enzyme structure and function.

Computational Enzyme Redesign Enhances Tolerance to Denaturants for Peptide C-Terminal Amidation

The escalating demand for biocatalysts in pharmaceutical and biochemical applications underscores the critical imperative to enhance enzyme activity and durability under high denaturant concentrations. Nevertheless, the development of a practical computational redesign protocol for improving enzyme tolerance to denaturants is challenging due to the limitations of relying solely on model-driven approaches to adequately capture denaturant-enzyme interactions. In this study, we introduce an enzyme redesign strategy termed GRAPE_DA, which integrates multiple data-driven and model-driven computational methods to mitigate the sampling biases inherent in a single approach and comprehensively predict beneficial mutations on both the protein surface and backbone. To illustrate the methodology's effectiveness, we applied it to engineer a peptidylamidoglycolate lyase, resulting in a variant exhibiting up to a 24-fold increase in peptide C-terminal amidation activity under 2.5 M guanidine hydrochloride. We anticipate that this integrated engineering strategy will facilitate the development of enzymatic peptide synthesis and functionalization under denaturing conditions and highlight the role of engineering surface residues in governing protein stability.

In silico design of multipoint mutants for enhanced performance of Thermomyces lanuginosus lipase for efficient biodiesel production

BACKGROUND

Biodiesel, an emerging sustainable and renewable clean energy, has garnered considerable attention as an alternative to fossil fuels. Although lipases are promising catalysts for biodiesel production, their efficiency in industrial-scale application still requires improvement.

RESULTS

In this study, a novel strategy for multi-site mutagenesis in the binding pocket was developed via FuncLib (for mutant enzyme design) and Rosetta Cartesian_ddg (for free energy calculation) to improve the reaction rate and yield of lipase-catalyzed biodiesel production. Thermomyces lanuginosus lipase (TLL) with high activity and thermostability was obtained using the Pichia pastoris expression system. The specific activities of the mutants M11 and M21 (each with 5 and 4 mutations) were 1.50- and 3.10-fold higher, respectively, than those of the wild-type (wt-TLL). Their corresponding melting temperature profiles increased by 10.53 and 6.01 °C, [Formula: see text] (the temperature at which the activity is reduced to 50% after 15 min incubation) increased from 60.88 to 68.46 °C and 66.30 °C, and the optimum temperatures shifted from 45 to 50 °C. After incubation in 60% methanol for 1 h, the mutants M11 and M21 retained more than 60% activity, and 45% higher activity than that of wt-TLL. Molecular dynamics simulations indicated that the increase in thermostability could be explained by reduced atomic fluctuation, and the improved catalytic properties were attributed to a reduced binding free energy and newly formed hydrophobic interaction. Yields of biodiesel production catalyzed by mutants M11 and M21 for 48 h at an elevated temperature (50 °C) were 94.03% and 98.56%, respectively, markedly higher than that of the wt-TLL (88.56%) at its optimal temperature (45 °C) by transesterification of soybean oil.

CONCLUSIONS

An integrating strategy was first adopted to realize the co-evolution of catalytic efficiency and thermostability of lipase. Two promising mutants M11 and M21 with excellent properties exhibited great potential for practical applications for in biodiesel production.

Leveraging Artificial Intelligence to Expedite Antibody Design and Enhance Antibody-Antigen Interactions

This perspective sheds light on the transformative impact of recent computational advancements in the field of protein therapeutics, with a particular focus on the design and development of antibodies. Cutting-edge computational methods have revolutionized our understanding of protein-protein interactions (PPIs), enhancing the efficacy of protein therapeutics in preclinical and clinical settings. Central to these advancements is the application of machine learning and deep learning, which offers unprecedented insights into the intricate mechanisms of PPIs and facilitates precise control over protein functions. Despite these advancements, the complex structural nuances of antibodies pose ongoing challenges in their design and optimization. Our review provides a comprehensive exploration of the latest deep learning approaches, including language models and diffusion techniques, and their role in surmounting these challenges. We also present a critical analysis of these methods, offering insights to drive further progress in this rapidly evolving field. The paper includes practical recommendations for the application of these computational techniques, supplemented with independent benchmark studies. These studies focus on key performance metrics such as accuracy and the ease of program execution, providing a valuable resource for researchers engaged in antibody design and development. Through this detailed perspective, we aim to contribute to the advancement of antibody design, equipping researchers with the tools and knowledge to navigate the complexities of this field.

Transfer learning to leverage larger datasets for improved prediction of protein stability changes

Amino acid mutations that lower a protein's thermodynamic stability are implicated in numerous diseases, and engineered proteins with enhanced stability can be important in research and medicine. Computational methods for predicting how mutations perturb protein stability are, therefore, of great interest. Despite recent advancements in protein design using deep learning, in silico prediction of stability changes has remained challenging, in part due to a lack of large, high-quality training datasets for model development. Here, we describe ThermoMPNN, a deep neural network trained to predict stability changes for protein point mutations given an initial structure. In doing so, we demonstrate the utility of a recently released megascale stability dataset for training a robust stability model. We also employ transfer learning to leverage a second, larger dataset by using learned features extracted from ProteinMPNN, a deep neural network trained to predict a protein's amino acid sequence given its three-dimensional structure. We show that our method achieves state-of-the-art performance on established benchmark datasets using a lightweight model architecture that allows for rapid, scalable predictions. Finally, we make ThermoMPNN readily available as a tool for stability prediction and design.

Revisiting structural organization of proteins at high temperature from a network perspective

Interactions between distantly placed amino acids in the primary chain (long-range) play a very crucial role in the formation and stabilization of the tertiary structure of a protein, while interactions between closely placed amino acids in the primary chain (short-range) mostly stabilize the secondary structures. Every protein needs to maintain marginal stability in order to perform its physiological functions in its native environment. The requirements for this stability in mesophilic and thermophilic proteins are different. Thermophilic proteins need to form more interactions as well as more stable interactions to survive in the extreme environment, they live in. Here, we aim to find out how the interacting amino acids in three-dimensional space are positioned in the primary chains in thermophilic and mesophilic. How does this arrangement help thermophiles to maintain their structural integrity at high temperatures? Working on a dataset of 1560 orthologous pairs we perceive that thermophiles are not only enriched with long-range interactions, they feature bigger connected clusters and higher network densities compared to their mesophilic orthologs, at higher interaction strengths between the amino acids. Moreover, we have observed the enrichment of different types of interactions at different secondary structural regions.

ProSTAGE: Predicting Effects of Mutations on Protein Stability by Using Protein Embeddings and Graph Convolutional Networks

Protein thermodynamic stability is essential to clarify the relationships among structure, function, and interaction. Therefore, developing a faster and more accurate method to predict the impact of the mutations on protein stability is helpful for protein design and understanding the phenotypic variation. Recent studies have shown that protein embedding will be particularly powerful at modeling sequence information with context dependence, such as subcellular localization, variant effect, and secondary structure prediction. Herein, we introduce a novel method, ProSTAGE, which is a deep learning method that fuses structure and sequence embedding to predict protein stability changes upon single point mutations. Our model combines graph-based techniques and language models to predict stability changes. Moreover, ProSTAGE is trained on a larger data set, which is almost twice as large as the most used S2648 data set. It consistently outperforms all existing state-of-the-art methods on mutation-affected problems as benchmarked on several independent data sets. The protein embedding as the prediction input achieves better results than the previous results, which shows the potential of protein language models in predicting the effect of mutations on proteins. ProSTAGE is implemented as a user-friendly web server.

Statistical Analysis of the Role of Cavity Flexibility in Thermostability of Proteins

Conventional statistical investigations have primarily focused on the comparison of the simple one-dimensional characteristics of protein cavities, such as number, surface area, and volume. These studies have failed to discern the crucial distinctions in cavity properties between thermophilic and mesophilic proteins that contribute to protein thermostability. In this study, the significance of cavity properties, i.e., flexibility and location, in protein thermostability was investigated by comparing structural differences between homologous thermophilic and mesophilic proteins. Three dimensions of protein structure were categorized into three regions (core, boundary, and surface) and a comparative analysis of cavity properties using this structural index was conducted. The statistical analysis revealed that cavity flexibility is closely related to protein thermostability. The core cavities of thermophilic proteins were less flexible than those of mesophilic proteins (averaged B' factor values, -0.6484 and -0.5111), which might be less deleterious to protein thermostability. Thermophilic proteins exhibited fewer cavities in the boundary and surface regions. Notably, cavities in mesophilic proteins, across all regions, exhibited greater flexibility than those in thermophilic proteins (>95% probability). The increased flexibility of cavities in the boundary and surface regions of mesophilic proteins, as opposed to thermophilic proteins, may compromise stability. Recent protein engineering investigations involving mesophilic xylanase and protease showed results consistent with the findings of this study, suggesting that the manipulation of flexible cavities in the surface region can enhance thermostability. Consequently, our findings suggest that a rational or computational approach to the design of flexible cavities in surface or boundary regions could serve as an effective strategy to enhance the thermostability of mesophilic proteins.

Enzyme engineering for functional lipids synthesis: recent advance and perspective

Functional lipids, primarily derived through the modification of natural lipids by various processes, are widely acknowledged for their potential to impart health benefits. In contrast to chemical methods for lipid modification, enzymatic catalysis offers distinct advantages, including high selectivity, mild operating conditions, and reduced byproduct formation. Nevertheless, enzymes face challenges in industrial applications, such as low activity, stability, and undesired selectivity. To address these challenges, protein engineering techniques have been implemented to enhance enzyme performance in functional lipid synthesis. This article aims to review recent advances in protein engineering, encompassing approaches from directed evolution to rational design, with the goal of improving the properties of lipid-modifying enzymes. Furthermore, the article explores the future prospects and challenges associated with enzyme-catalyzed functional lipid synthesis.

Assessing computational tools for predicting protein stability changes upon missense mutations using a new dataset

Insight into how mutations affect protein stability is crucial for protein engineering, understanding genetic diseases, and exploring protein evolution. Numerous computational methods have been developed to predict the impact of amino acid substitutions on protein stability. Nevertheless, comparing these methods poses challenges due to variations in their training data. Moreover, it is observed that they tend to perform better at predicting destabilizing mutations than stabilizing ones. Here, we meticulously compiled a new dataset from three recently published databases: ThermoMutDB, FireProtDB, and ProThermDB. This dataset, which does not overlap with the well-established S2648 dataset, consists of 4038 single-point mutations, including over 1000 stabilizing mutations. We assessed these mutations using 27 computational methods, including the latest ones utilizing mega-scale stability datasets and transfer learning. We excluded entries with overlap or similarity to training datasets to ensure fairness. Pearson correlation coefficients for the tested tools ranged from 0.20 to 0.53 on unseen data, and none of the methods could accurately predict stabilizing mutations, even those performing well in anti-symmetric property analysis. While most methods present consistent trends for predicting destabilizing mutations across various properties such as solvent exposure and secondary conformation, stabilizing mutations do not exhibit a clear pattern. Our study also suggests that solely addressing training dataset bias may not significantly enhance accuracy of predicting stabilizing mutations. These findings emphasize the importance of developing precise predictive methods for stabilizing mutations.

Evaluating the industrial potential of naturally occurring proteases: A focus on kinetic and thermodynamic parameters

Thermodynamic and kinetic parameters, such as enthalpy, entropy, and free energy, are crucial in evaluating enzyme stability and activity. These parameters, including the free energy of activation (ΔG#) and the Gibbs free energy of inactivation (ΔG*), are important for predicting energy requirements and reaction rates. However, relying solely on these parameters is insufficient in selecting an enzyme for industrial processes. Numerous studies have explored the measurement of thermodynamic parameters for proteases. Unfortunately, some of the definitions and calculations of key parameters such as ΔG#, ΔG*, and substrate-binding free energy have contained significant errors. In this study, these mistakes have been addressed and corrected. Additionally, a new parameter called δ, defined as the difference between ΔG* and ΔG#, has been introduced for the first time. It is argued that δ provides a more reliable measure for predicting the potential industrial application of enzymes. The highest calculated value for δ was found to be 39.6 kJ·mol-1 at 55 °C. Furthermore, this study also presents a comprehensive collection and determination of all thermodynamic and kinetic parameters for proteases, providing researchers and professionals in the field with a valuable resource to compare and understand the relationships between these parameters and the industrial potential of enzymes.

Calculation of Protein Folding Thermodynamics Using Molecular Dynamics Simulations

Despite advances in artificial intelligence methods, protein folding remains in many ways an enigma to be solved. Accurate computation of protein folding energetics could help drive fields such as protein and drug design and genetic interpretation. However, the challenge of calculating the state functions governing protein folding from first-principles remains unaddressed. We present here a simple approach that allows us to accurately calculate the energetics of protein folding. It is based on computing the energy of the folded and unfolded states at different temperatures using molecular dynamics simulations. From this, two essential quantities (ΔH and ΔCp) are obtained and used to calculate the conformational stability of the protein (ΔG). With this approach, we have successfully calculated the energetics of two- and three-state proteins, representatives of the major structural classes, as well as small stability differences (ΔΔG) due to changes in solution conditions or variations in an amino acid residue.

Natural compound targeting BDNF V66M variant: insights from in silico docking and molecular analysis

Brain-Derived Neurotrophic Factor (BDNF) is a neurotrophin gene family gene that encodes proteins vital for the growth, maintenance, and survival of neurons in the nervous system. The study aimed to screen natural compounds against BDNF variant (V66M), which affects memory, cognition, and mood regulation. BDNF variant (V66M) as a target structure was selected, and Vitamin D, Curcumin, Vitamin C, and Quercetin as ligands structures were taken from PubChem database. Multiple tools like AUTODOCK VINA, BIOVIA discovery studio, PyMOL, CB-dock, IMOD server, Swiss ADEMT, and Swiss predict ligands target were used to analyze binding energy, interaction, stability, toxicity, and visualize BDNF-ligand complexes. Compounds Vitamin D3, Curcumin, Vitamin C, and Quercetin with binding energies values of - 5.5, - 6.1, - 4.5, and - 6.7 kj/mol, respectively, were selected. The ligands bind to the active sites of the BDNF variant (V66M) via hydrophobic bonds, hydrogen bonds, and electrostatic interactions. Furthermore, ADMET analysis of the ligands revealed they exhibited sound pharmacokinetic and toxicity profiles. In addition, an MD simulation study showed that the most active ligand bound favorably and dynamically to the target protein, and protein-ligand complex stability was determined. The finding of this research could provide an excellent platform for discovering and rationalizing novel drugs against stress related to BDNF (V66M). Docking, preclinical drug testing and MD simulation results suggest Quercetin as a more potent BDNF variant (V66M) inhibitor and forming a more structurally stable complex.

Computational design of highly efficient thermostable MHET hydrolases and dual enzyme system for PET recycling

Recently developed enzymes for the depolymerization of polyethylene terephthalate (PET) such as FAST-PETase and LCC-ICCG are inhibited by the intermediate PET product mono(2-hydroxyethyl) terephthalate (MHET). Consequently, the conversion of PET enzymatically into its constituent monomers terephthalic acid (TPA) and ethylene glycol (EG) is inefficient. In this study, a protein scaffold (1TQH) corresponding to a thermophilic carboxylesterase (Est30) was selected from the structural database and redesigned in silico. Among designs, a double variant KL-MHETase (I171K/G130L) with a similar protein melting temperature (67.58 °C) to that of the PET hydrolase FAST-PETase (67.80 °C) exhibited a 67-fold higher activity for MHET hydrolysis than FAST-PETase. A fused dual enzyme system comprising KL-MHETase and FAST-PETase exhibited a 2.6-fold faster PET depolymerization rate than FAST-PETase alone. Synergy increased the yield of TPA by 1.64 fold, and its purity in the released aromatic products reached 99.5%. In large reaction systems with 100 g/L substrate concentrations, the dual enzyme system KL36F achieved over 90% PET depolymerization into monomers, demonstrating its potential applicability in the industrial recycling of PET plastics. Therefore, a dual enzyme system can greatly reduce the reaction and separation cost for sustainable enzymatic PET recycling.

Cooperative Assembly of Self-Adjusting α-Helical Coiled Coils along the Length of an mRNA Chain to Form a Thermodynamically Stable Nanotube Carrier

Although mRNA delivery technology is very promising, problems in safety and transport arise due to the intrinsically low thermodynamic stability of the current mRNA carriers. Considering that mRNAs are filamentous and a nanotube is one of the most thermodynamically stable shapes among nanoassemblies, a nanotube is one of the most stable supramolecular structures that can be assembled with mRNA. Here, we develop a nanotube-shaped filamentous mRNA delivery platform that shows exceptionally high thermodynamic stability. The key to the development of the mRNA nanotube is the design of self-adjusting supramolecular building blocks (SABs) that have two disparate properties, i.e., dynamic property and stiffness, in a single molecule. The counterbalance of the dynamic property and stiffness in SABs enables the coating of mRNA by winding its way through the flexible and irregular mRNA chain via cooperative interactions. SAB nanotubes with targeting ligands installed show a high uptake efficiency in mammalian cells and controllable gene expression behavior. Thus, the mRNA nanotube provides an enabling technology toward the development of safe and stable mRNA vaccines and therapeutics.

From De Novo Design to Redesign: Harnessing Computational Protein Design for Understanding SARS-CoV-2 Molecular Mechanisms and Developing Therapeutics

The continuous emergence of novel SARS-CoV-2 variants and subvariants serves as compelling evidence that COVID-19 is an ongoing concern. The swift, well-coordinated response to the pandemic highlights how technological advancements can accelerate the detection, monitoring, and treatment of the disease. Robust surveillance systems have been established to understand the clinical characteristics of new variants, although the unpredictable nature of these variants presents significant challenges. Some variants have shown resistance to current treatments, but innovative technologies like computational protein design (CPD) offer promising solutions and versatile therapeutics against SARS-CoV-2. Advances in computing power, coupled with open-source platforms like AlphaFold and RFdiffusion (employing deep neural network and diffusion generative models), among many others, have accelerated the design of protein therapeutics with precise structures and intended functions. CPD has played a pivotal role in developing peptide inhibitors, mini proteins, protein mimics, decoy receptors, nanobodies, monoclonal antibodies, identifying drug-resistance mutations, and even redesigning native SARS-CoV-2 proteins. Pending regulatory approval, these designed therapies hold the potential for a lasting impact on human health and sustainability. As SARS-CoV-2 continues to evolve, use of such technologies enables the ongoing development of alternative strategies, thus equipping us for the "New Normal".

K-Pro: Kinetics Data on Proteins and Mutants

The study of protein folding plays a crucial role in improving our understanding of protein function and of the relationship between genetics and phenotypes. In particular, understanding the thermodynamics and kinetics of the folding process is important for uncovering the mechanisms behind human disorders caused by protein misfolding. To address this issue, it is essential to collect and curate experimental kinetic and thermodynamic data on protein folding. K-Pro is a new database designed for collecting and storing experimental kinetic data on monomeric proteins, with a two-state folding mechanism. With 1,529 records from 62 proteins corresponding to 65 structures, K-Pro contains various kinetic parameters such as the logarithm of the folding and unfolding rates, Tanford's β and the ϕ values. When available, the database also includes thermodynamic parameters associated with the kinetic data. K-Pro features a user-friendly interface that allows browsing and downloading kinetic data of interest. The graphical interface provides a visual representation of the protein and mutants, and it is cross-linked to key databases such as PDB, UniProt, and PubMed. K-Pro is open and freely accessible through https://folding.biofold.org/k-pro and supports the latest versions of popular browsers.

Stable Mammalian Serum Albumins Designed for Bacterial Expression

Albumin is the most abundant protein in the blood serum of mammals and has essential carrier and physiological roles. Albumins are also used in a wide variety of molecular and cellular experiments and in the cultivated meat industry. Despite their importance, however, albumins are challenging for heterologous expression in microbial hosts, likely due to 17 conserved intramolecular disulfide bonds. Therefore, albumins used in research and biotechnological applications either derive from animal serum, despite severe ethical and reproducibility concerns, or from recombinant expression in yeast or rice. We use the PROSS algorithm to stabilize human and bovine serum albumins, finding that all are highly expressed in E. coli. Design accuracy is verified by crystallographic analysis of a human albumin variant with 16 mutations. This albumin variant exhibits ligand binding properties similar to those of the wild type. Remarkably, a design with 73 mutations relative to human albumin exhibits over 40 °C improved stability and is stable beyond the boiling point of water. Our results suggest that proteins with many disulfide bridges have the potential to exhibit extreme stability when subjected to design. The designed albumins may be used to make economical, reproducible, and animal-free reagents for molecular and cell biology. They also open the way to high-throughput screening to study and enhance albumin carrier properties.

Transfer learning to leverage larger datasets for improved prediction of protein stability changes

Amino acid mutations that lower a protein's thermodynamic stability are implicated in numerous diseases, and engineered proteins with enhanced stability are important in research and medicine. Computational methods for predicting how mutations perturb protein stability are therefore of great interest. Despite recent advancements in protein design using deep learning, in silico prediction of stability changes has remained challenging, in part due to a lack of large, high-quality training datasets for model development. Here we introduce ThermoMPNN, a deep neural network trained to predict stability changes for protein point mutations given an initial structure. In doing so, we demonstrate the utility of a newly released mega-scale stability dataset for training a robust stability model. We also employ transfer learning to leverage a second, larger dataset by using learned features extracted from a deep neural network trained to predict a protein's amino acid sequence given its three-dimensional structure. We show that our method achieves competitive performance on established benchmark datasets using a lightweight model architecture that allows for rapid, scalable predictions. Finally, we make ThermoMPNN readily available as a tool for stability prediction and design.

A novel computationally engineered collagenase reduces the force required for tooth extraction in an ex-situ porcine jaw model

The currently employed tooth extraction methods in dentistry involve mechanical disruption of the periodontal ligament fibers, leading to inevitable trauma to the bundle bone comprising the socket walls. In our previous work, we have shown that a recombinantly expressed truncated version of clostridial collagenase G (ColG) purified from Escherichia coli efficiently reduced the force needed for tooth extraction in an ex-situ porcine jaw model, when injected into the periodontal ligament. Considering that enhanced thermostability often leads to higher enzymatic activity and to set the basis for additional rounds of optimization, we used a computational protein design approach to generate an enzyme to be more thermostable while conserving the key catalytic residues. This process generated a novel collagenase (ColG-variant) harboring sixteen mutations compared to ColG, with a nearly 4℃ increase in melting temperature. Herein, we explored the potential of ColG-variant to further decrease the physical effort required for tooth delivery using our established ex-situ porcine jaw model. An average reduction of 11% was recorded in the force applied to extract roots of mandibular split first and second premolar teeth treated with ColG-variant, relative to those treated with ColG. Our results show for the first time the potential of engineering enzyme properties for dental medicine and further contribute to minimally invasive tooth extraction.

The "Beacon" Structural Model of Protein Folding: Application for Trp-Cage in Water

Protein folding is a process in which a polypeptide must undergo folding process to obtain its three-dimensional structure. Thermodynamically, it is a process of enthalpy to overcome the loss of conformational entropy in folding. Folding is primarily related to hydrophobic interactions and intramolecular hydrogen bondings. During folding, hydrophobic interactions are regarded to be the driving forces, especially in the initial structural collapse of a protein. Additionally, folding is guided by the strong interactions within proteins, such as intramolecular hydrogen bondings related to the α-helices and β-sheets of proteins. Therefore, a protein is divided into the folding key (FK) regions related to intramolecular hydrogen bondings and the non-folding key (non-FK) regions. Various conformations are expected for FK and non-FK regions. Different from non-FK regions, it is necessary for FK regions to form the specific conformations in folding, which are regarded as the necessary folding pathways (or "beacons"). Additionally, sequential folding is expected for the FK regions, and the intermediate state is found during folding. They are reflected on the local basins in the free energy landscape (FEL) of folding. To demonstrate the structural model, molecular dynamics (MD) simulations are conducted on the folding pathway of the TRP-cage in water.

Designed active-site library reveals thousands of functional GFP variants

Mutations in a protein active site can lead to dramatic and useful changes in protein activity. The active site, however, is sensitive to mutations due to a high density of molecular interactions, substantially reducing the likelihood of obtaining functional multipoint mutants. We introduce an atomistic and machine-learning-based approach, called high-throughput Functional Libraries (htFuncLib), that designs a sequence space in which mutations form low-energy combinations that mitigate the risk of incompatible interactions. We apply htFuncLib to the GFP chromophore-binding pocket, and, using fluorescence readout, recover >16,000 unique designs encoding as many as eight active-site mutations. Many designs exhibit substantial and useful diversity in functional thermostability (up to 96 °C), fluorescence lifetime, and quantum yield. By eliminating incompatible active-site mutations, htFuncLib generates a large diversity of functional sequences. We envision that htFuncLib will be used in one-shot optimization of activity in enzymes, binders, and other proteins.

Accelerating therapeutic protein design with computational approaches toward the clinical stage

Therapeutic protein, represented by antibodies, is of increasing interest in human medicine. However, clinical translation of therapeutic protein is still largely hindered by different aspects of developability, including affinity and selectivity, stability and aggregation prevention, solubility and viscosity reduction, and deimmunization. Conventional optimization of the developability with widely used methods, like display technologies and library screening approaches, is a time and cost-intensive endeavor, and the efficiency in finding suitable solutions is still not enough to meet clinical needs. In recent years, the accelerated advancement of computational methodologies has ushered in a transformative era in the field of therapeutic protein design. Owing to their remarkable capabilities in feature extraction and modeling, the integration of cutting-edge computational strategies with conventional techniques presents a promising avenue to accelerate the progression of therapeutic protein design and optimization toward clinical implementation. Here, we compared the differences between therapeutic protein and small molecules in developability and provided an overview of the computational approaches applicable to the design or optimization of therapeutic protein in several developability issues.

Computational design and molecular dynamics simulations suggest the mode of substrate binding in ceramide synthases

Until now, membrane-protein stabilization has relied on iterations of mutations and screening. We now validate a one-step algorithm, mPROSS, for stabilizing membrane proteins directly from an AlphaFold2 model structure. Applied to the lipid-generating enzyme, ceramide synthase, 37 designed mutations lead to a more stable form of human CerS2. Together with molecular dynamics simulations, we propose a pathway by which substrates might be delivered to the ceramide synthases.

Automated optimisation of solubility and conformational stability of antibodies and proteins

Biologics, such as antibodies and enzymes, are crucial in research, biotechnology, diagnostics, and therapeutics. Often, biologics with suitable functionality are discovered, but their development is impeded by developability issues. Stability and solubility are key biophysical traits underpinning developability potential, as they determine aggregation, correlate with production yield and poly-specificity, and are essential to access parenteral and oral delivery. While advances for the optimisation of individual traits have been made, the co-optimization of multiple traits remains highly problematic and time-consuming, as mutations that improve one property often negatively impact others. In this work, we introduce a fully automated computational strategy for the simultaneous optimisation of conformational stability and solubility, which we experimentally validate on six antibodies, including two approved therapeutics. Our results on 42 designs demonstrate that the computational procedure is highly effective at improving developability potential, while not affecting antigen-binding. We make the method available as a webserver at www-cohsoftware.ch.cam.ac.uk.

AsiteDesign: a Semirational Algorithm for an Automated Enzyme Design

With advances in protein structure predictions, the number of available high-quality structures has increased dramatically. In light of these advances, structure-based enzyme engineering is expected to become increasingly important for optimizing biocatalysts for industrial processes. Here, we present AsiteDesign, a Monte Carlo-based protocol for structure-based engineering of active sites. AsiteDesign provides a framework for introducing new catalytic residues in a given binding pocket to either create a new catalytic activity or alter the existing one. AsiteDesign is implemented using pyRosetta and incorporates enhanced sampling techniques to efficiently explore the search space. The protocol was tested by designing an alternative catalytic triad in the active site of Pseudomonas fluorescens esterase (PFE). The designed variant was experimentally verified to be active, demonstrating that AsiteDesign can find alternative catalytic triads. Additionally, the AsiteDesign protocol was employed to enhance the hydrolysis of a bulky chiral substrate (1-phenyl-2-pentyl acetate) by PFE. The experimental verification of the designed variants demonstrated that F158L/F198A and F125A/F158L mutations increased the hydrolysis of 1-phenyl-2-pentyl acetate from 8.9 to 66.7 and 23.4%, respectively, and reversed the enantioselectivity of the enzyme from (R) to (S)-enantiopreference, with 32 and 55% enantiomeric excess (ee), respectively.

Domain-wise dissection of thermal stability enhancement in multidomain proteins

Stability is critical for the proper functioning of all proteins. Optimization of protein thermostability is a key step in the development of industrial enzymes and biologics. Herein, we demonstrate that multidomain proteins can be stabilized significantly using domain-based engineering followed by the recombination of the optimized domains. Domain-level analysis of designed protein variants with similar structures but different thermal profiles showed that the independent enhancement of the thermostability of a constituent domain improves the overall stability of the whole multidomain protein. The crystal structure and AlphaFold-predicted model of the designed proteins via domain-recombination provided a molecular explanation for domain-based stepwise stabilization. Our study suggests that domain-based modular engineering can minimize the sequence space for calculations in computational design and experimental errors, thereby offering useful guidance for multidomain protein engineering.

Target-template relationships in protein structure prediction and their effect on the accuracy of thermostability calculations

Improving protein thermostability has been a labor- and time-consuming process in industrial applications of protein engineering. Advances in computational approaches have facilitated the development of more efficient strategies to allow the prioritization of stabilizing mutants. Among these is FEP+, a free energy perturbation implementation that uses a thoroughly tested physics-based method to achieve unparalleled accuracy in predicting changes in protein thermostability. To gauge the applicability of FEP+ to situations where crystal structures are unavailable, here we have applied the FEP+ approach to homology models of 12 different proteins covering 316 mutations. By comparing predictions obtained with homology models to those obtained using crystal structures, we have identified that local rather than global sequence conservation between target and template sequence is a determining factor in the accuracy of predictions. By excluding mutation sites with low local sequence identity (<40%) to a template structure, we have obtained predictions with comparable performance to crystal structures (R2 of 0.67 and 0.63 and an RMSE of 1.20 and 1.16 kcal/mol for crystal structure and homology model predictions, respectively) for identifying stabilizing mutations when incorporating residue scanning into a cascade screening strategy. Additionally, we identify and discuss inherent limitations in sequence alignments and homology modeling protocols that translate into the poor FEP+ performance of a few select examples. Overall, our retrospective study provides detailed guidelines for the application of the FEP+ approach using homology models for protein thermostability predictions, which will greatly extend this approach to studies that were previously limited by structure availability.

Combinatorial assembly and design of enzymes

The design of structurally diverse enzymes is constrained by long-range interactions that are necessary for accurate folding. We introduce an atomistic and machine learning strategy for the combinatorial assembly and design of enzymes (CADENZ) to design fragments that combine with one another to generate diverse, low-energy structures with stable catalytic constellations. We applied CADENZ to endoxylanases and used activity-based protein profiling to recover thousands of structurally diverse enzymes. Functional designs exhibit high active-site preorganization and more stable and compact packing outside the active site. Implementing these lessons into CADENZ led to a 10-fold improved hit rate and more than 10,000 recovered enzymes. This design-test-learn loop can be applied, in principle, to any modular protein family, yielding huge diversity and general lessons on protein design principles.

Insights Into The Effects of Amino Acid Substitutions on The Stability of Reteplase Structure: A Molecular Dynamics Simulation Study

Background

Reteplase (recombinant plasminogen activator, r-PA) is a recombinant protein designed to imitate the endogenous tissue plasminogen activator and catalyze the plasmin production. It is known that the application of reteplase is limited by the complex production processes and protein's stability challenges. Computational redesign of proteins has gained momentum in recent years, particularly as a powerful tool for improving protein stability and consequently its production efficiency. Hence, in the current study, we implemented computational approaches to improve r-PA conformational stability, which fairly correlates with protein's resistance to proteolysis.

Objectives

The current study was developed in order to evaluate the effect of amino acid substitutions on the stability of reteplase structure using molecular dynamic simulations and computational predictions.

Materials and Methods

Several web servers designed for mutation analysis were utilized to select appropriate mutations. Additionally, the experimentally reported mutation, R103S, converting wild type r-PA into non-cleavable form, was also employed. Firstly, mutant collection, consisting of 15 structures, was constructed based on the combinations of four designated mutations. Then, 3D structures were generated using MODELLER. Finally, 17 independent 20-ns molecular dynamics (MD) simulations were conducted and different analysis were performed like root-mean-square deviation (RMSD), root-mean-square fluctuations (RMSF), secondary structure analysis, number of hydrogen bonds, principal components analysis (PCA), eigenvector projection, and density analysis.

Results

Predicted mutations successfully compensated the more flexible conformation caused by R103S substitution, so, improved conformational stability was analyzed from MD simulations. In particular, R103S/A286I/G322I indicated the best results and remarkably enhanced the protein stability.

Conclusion

The conformational stability conferred by these mutations will probably lead to more protection of r-PA in protease-rich environments in various recombinant systems and potentially enhance its production and expression level.

Exploration of Computational Aids for Effective Drug Designing and Management of Viral Diseases: A Comprehensive Review

BACKGROUND

Microbial diseases, specifically originating from viruses are the major cause of human mortality all over the world. The current COVID-19 pandemic is a case in point, where the dynamics of the viral-human interactions are still not completely understood, making its treatment a case of trial and error. Scientists are struggling to devise a strategy to contain the pandemic for over a year and this brings to light the lack of understanding of how the virus grows and multiplies in the human body.

METHODS

This paper presents the perspective of the authors on the applicability of computational tools for deep learning and understanding of host-microbe interaction, disease progression and management, drug resistance and immune modulation through in silico methodologies which can aid in effective and selective drug development. The paper has summarized advances in the last five years. The studies published and indexed in leading databases have been included in the review.

RESULTS

Computational systems biology works on an interface of biology and mathematics and intends to unravel the complex mechanisms between the biological systems and the inter and intra species dynamics using computational tools, and high-throughput technologies developed on algorithms, networks and complex connections to simulate cellular biological processes.

CONCLUSION

Computational strategies and modelling integrate and prioritize microbial-host interactions and may predict the conditions in which the fine-tuning attenuates. These microbial-host interactions and working mechanisms are important from the aspect of effective drug designing and fine- tuning the therapeutic interventions.

Hydration in Deep Eutectic Solvents Induces Non-monotonic Changes in the Conformation and Stability of Proteins

The preservation of labile biomolecules presents a major challenge in chemistry, and deep eutectic solvents (DESs) have emerged as suitable environments for this purpose. However, how the hydration of DESs impacts the behavior of proteins is often neglected. Here, we demonstrate that the amino acid environment and secondary structure of two proteins (bovine serum albumin and lysozyme) and an antibody (immunoglobulin G) in 1:2 choline chloride:glycerol and 1:2 choline chloride:urea follow a re-entrant behavior with solvent hydration. A dome-shaped transition is observed with a folded or partially folded structure at very low (<10 wt % H₂O) and high (>40 wt % H₂O) DES hydration, while protein unfolding increases between those regimes. Hydration also affects protein conformation and stability, as demonstrated for bovine serum albumin in hydrated 1:2 choline chloride:glycerol. In the neat DES, bovine serum albumin remains partially folded and unexpectedly undergoes unfolding and oligomerization at low water content. At intermediate hydration, the protein begins to refold and gradually retrieves the native monomer-dimer equilibrium. However, ca. 36 wt % H₂O is required to recover the native folding fully. The half-denaturation temperature of the protein increases with decreasing hydration, but even the dilute DESs significantly enhance the thermal stability of bovine serum albumin. Also, protein unfolding can be reversed by rehydrating the sample to the high hydration regime, also recovering protein function. This correlation provides a new perspective to understanding protein behavior in hydrated DESs, where quantifying the DES hydration becomes imperative to identifying the folding and stability of proteins.

The use of pharmacological chaperones in rare diseases caused by reduced protein stability

Rare diseases are most often caused by inherited genetic disorders that, after translation, will result in a protein with altered function. Decreased protein stability is the most frequent mechanism associated with a congenital pathogenic missense mutation and it implies the destabilization of the folded conformation in favour of unfolded or misfolded states. In the cellular context and when experimental data is available, a mutant protein with altered thermodynamic stability often also results in impaired homeostasis, with the deleterious accumulation of protein aggregates, metabolites and/or metabolic by-products. In the last decades, a significant effort has enabled the characterization of rare diseases associated to protein stability defects and triggered the development of innovative therapeutic intervention lines, say, the use of pharmacological chaperones to correct the intracellular impaired homeostasis. Here, we review the current knowledge on rare diseases caused by reduced protein stability, paying special attention to the thermodynamic aspects of the protein destabilization, also focusing on some examples where pharmacological chaperones are being tested.

Accurate protein stability predictions from homology models

Calculating changes in protein stability (ΔΔG) has been shown to be central for predicting the consequences of single amino acid substitutions in protein engineering as well as interpretation of genomic variants for disease risk. Structure-based calculations are considered most accurate, however the tools used to calculate ΔΔGs have been developed on experimentally resolved structures. Extending those calculations to homology models based on related proteins would greatly extend their applicability as large parts of e.g. the human proteome are not structurally resolved. In this study we aim to investigate the accuracy of ΔΔG values predicted on homology models compared to crystal structures. Specifically, we identified four proteins with a large number of experimentally tested ΔΔGs and templates for homology modeling across a broad range of sequence identities, and selected three methods for ΔΔG calculations to test. We find that ΔΔG-values predicted from homology models compare equally well to experimental ΔΔGs as those predicted on experimentally established crystal structures, as long as the sequence identity of the model template to the target protein is at least 40%. In particular, the Rosetta cartesian_ddg protocol is robust against the small perturbations in the structure which homology modeling introduces. In an independent assessment, we observe a similar trend when using ΔΔGs to categorize variants as low or wild-type-like abundance. Overall, our results show that stability calculations performed on homology models can substitute for those on crystal structures with acceptable accuracy as long as the model is built on a template with sequence identity of at least 40% to the target protein.

Designed High-Redox Potential Laccases Exhibit High Functional Diversity

White-rot fungi secrete an impressive repertoire of high-redox potential laccases (HRPLs) and peroxidases for efficient oxidation and utilization of lignin. Laccases are attractive enzymes for the chemical industry due to their broad substrate range and low environmental impact. Since expression of functional recombinant HRPLs is challenging, however, iterative-directed evolution protocols have been applied to improve their expression, activity, and stability. We implement a rational, stabilize-and-diversify strategy to two HRPLs that we could not functionally express. First, we use the PROSS stability-design algorithm to allow functional expression in yeast. Second, we use the stabilized enzymes as starting points for FuncLib active-site design to improve their activity and substrate diversity. Four of the FuncLib-designed HRPLs and their PROSS progenitor exhibit substantial diversity in reactivity profiles against high-redox potential substrates, including lignin monomers. Combinations of 3–4 subtle mutations that change the polarity, solvation, and sterics of the substrate-oxidation site result in orders of magnitude changes in reactivity profiles. These stable and versatile HRPLs are a step toward generating an effective lignin-degrading consortium of enzymes that can be secreted from yeast. The stabilize-and-diversify strategy can be applied to other challenging enzyme families to study and expand the utility of natural enzymes.

Dissecting the stability determinants of a challenging de novo protein fold using massively parallel design and experimentation

Designing entirely new protein structures remains challenging because we do not fully understand the biophysical determinants of folding stability. Yet, some protein folds are easier to design than others. Previous work identified the 43-residue ɑββɑ fold as especially challenging: The best designs had only a 2% success rate, compared to 39 to 87% success for other simple folds [G. J. Rocklin et al., Science 357, 168-175 (2017)]. This suggested the ɑββɑ fold would be a useful model system for gaining a deeper understanding of folding stability determinants and for testing new protein design methods. Here, we designed over 10,000 new ɑββɑ proteins and found over 3,000 of them to fold into stable structures using a high-throughput protease-based assay. NMR, hydrogen-deuterium exchange, circular dichroism, deep mutational scanning, and scrambled sequence control experiments indicated that our stable designs fold into their designed ɑββɑ structures with exceptional stability for their small size. Our large dataset enabled us to quantify the influence of universal stability determinants including nonpolar burial, helix capping, and buried unsatisfied polar atoms, as well as stability determinants unique to the ɑββɑ topology. Our work demonstrates how large-scale design and test cycles can solve challenging design problems while illuminating the biophysical determinants of folding.

De Novo Design of a Highly Stable Ovoid TIM Barrel: Unlocking Pocket Shape towards Functional Design

The ability to finely control the structure of protein folds is an important prerequisite to functional protein design. The TIM barrel fold is an important target for these efforts as it is highly enriched for diverse functions in nature. Although a TIM barrel protein has been designed de novo, the ability to finely alter the curvature of the central beta barrel and the overall architecture of the fold remains elusive, limiting its utility for functional design. Here, we report the de novo design of a TIM barrel with ovoid (twofold) symmetry, drawing inspiration from natural beta and TIM barrels with ovoid curvature. We use an autoregressive backbone sampling strategy to implement our hypothesis for elongated barrel curvature, followed by an iterative enrichment sequence design protocol to obtain sequences which yield a high proportion of successfully folding designs. Designed sequences are highly stable and fold to the designed barrel curvature as determined by a 2.1 Å resolution crystal structure. The designs show robustness to drastic mutations, retaining high melting temperatures even when multiple charged residues are buried in the hydrophobic core or when the hydrophobic core is ablated to alanine. As a scaffold with a greater capacity for hosting diverse hydrogen bonding networks and installation of binding pockets or active sites, the ovoid TIM barrel represents a major step towards the de novo design of functional TIM barrels.

Computationally Guided Redesign of a Heme-free Cytochrome with Native-like Structure and Stability

Metals can play key roles in stabilizing protein structures, but ensuring their proper incorporation is a challenge when a metalloprotein is overexpressed in a non-native cellular environment. Here, we have used computational protein design tools to redesign cytochrome b₅₆₂ (cyt b₅₆₂), which relies on the binding of its heme cofactor to achieve its proper fold, into a stable, heme-free protein. The resulting protein, ApoCyt, features only four mutations and no metal-ligand or covalent bonds, yet displays improved stability over cyt b₅₆₂. Mutagenesis studies and X-ray crystal structures reveal that the increase in stability is due to the computationally prescribed mutations, which stabilize the protein fold through a combination of hydrophobic packing interactions, hydrogen bonds, and cation-π interactions. Upon installation of the relevant mutations, ApoCyt is capable of assembling into previously reported, cytochrome-based trimeric and tetrameric assemblies, demonstrating that ApoCyt retains the structure and assembly properties of cyt b₅₆₂. The successful design of ApoCyt therefore enables further functional diversification of cytochrome-based assemblies and demonstrates that structural metal cofactors can be replaced by a small number of well-designed, non-covalent interactions.

Stabilizing proteins, simplified: A Rosetta-based webtool for predicting favorable mutations

Many proteins have low thermodynamic stability, which can lead to low expression yields and limit functionality in research, industrial and clinical settings. This article introduces two, web-based tools that use the high-resolution structure of a protein along with the Rosetta molecular modeling program to predict stabilizing mutations. The protocols were recently applied to three genetically and structurally distinct proteins and successfully predicted mutations that improved thermal stability and/or protein yield. In all three cases, combining the stabilizing mutations raised the protein unfolding temperatures by more than 20°C. The first protocol evaluates point mutations and can generate a site saturation mutagenesis heatmap. The second identifies mutation clusters around user-defined positions. Both applications only require a protein structure and are particularly valuable when a deep multiple sequence alignment is not available. These tools were created to simplify protein engineering and enable research that would otherwise be infeasible due to poor expression and stability of the native molecule.

Nonsynonymous amino acid changes in the α-chain of complement component 5 influence longitudinal susceptibility to Plasmodium falciparum infections and severe malarial anemia in kenyan children

Background: Severe malarial anemia (SMA; Hb < 5.0 g/dl) is a leading cause of childhood morbidity and mortality in holoendemic Plasmodium falciparum transmission regions such as western Kenya.

Methods: We investigated the relationship between two novel complement component 5 (C5) missense mutations [rs17216529:C>T, p(Val145Ile) and rs17610:C>T, p(Ser1310Asn)] and longitudinal outcomes of malaria in a cohort of Kenyan children (under 60 mos, n = 1,546). Molecular modeling was used to investigate the impact of the amino acid transitions on the C5 protein structure.

Results: Prediction of the wild-type and mutant C5 protein structures did not reveal major changes to the overall structure. However, based on the position of the variants, subtle differences could impact on the stability of C5b. The influence of the C5 genotypes/haplotypes on the number of malaria and SMA episodes over 36 months was determined by Poisson regression modeling. Genotypic analyses revealed that inheritance of the homozygous mutant (TT) for rs17216529:C>T enhanced the risk for both malaria (incidence rate ratio, IRR = 1.144, 95%CI: 1.059–1.236, p = 0.001) and SMA (IRR = 1.627, 95%CI: 1.201–2.204, p = 0.002). In the haplotypic model, carriers of TC had increased risk of malaria (IRR = 1.068, 95%CI: 1.017–1.122, p = 0.009), while carriers of both wild-type alleles (CC) were protected against SMA (IRR = 0.679, 95%CI: 0.542–0.850, p = 0.001).

Conclusion: Collectively, these findings show that the selected C5 missense mutations influence the longitudinal risk of malaria and SMA in immune-naïve children exposed to holoendemic P. falciparum transmission through a mechanism that remains to be defined.