Identification of stabilizing point mutations through mutagenesis of destabilized protein libraries

Molecular engineering of insulin for recombinant expression in yeast

Since the first administration of insulin to a person with diabetes in 1922, scientific contributions from academia and industry have improved insulin therapy and access. The pharmaceutical need for insulin is now more than 40 tons annually, half of which is produced by recombinant secretory expression in Saccharomyces cerevisiae. We discuss how, in this yeast species, adaptation of insulin precursors by removable structural elements is pivotal for efficient secretory expression. The technologies reviewed have been implemented at industrial scale and are seminal for the supply of human insulin and insulin analogues to people with diabetes now and in the future. Engineering of a target protein with removable structural elements may provide a general approach to yield optimisation.

A novel homozygous missense variant in TBC1D31 in a consanguineous family with congenital anomalies of the kidney and urinary tract (CAKUT)

Congenital anomalies of the kidney and urinary tract (CAKUT) is the leading cause of chronic kidney disease in the first three decades of life. Until now, more than 180 monogenic causes of isolated or syndromic CAKUT have been described. In addition, copy number variants (CNV) have also been implicated, however, all of these causative factors only explain a small fraction of patients with CAKUT, suggesting that additional yet-to-be-discovered novel genes are present. Herein, we report three siblings (two of them are monozygotic twin) of a consanguineous family with CAKUT. Whole-exome sequencing identified a homozygous variant in TBC1D31. Three dimensional protein modeling as well as molecular dynamics simulations predicted it as pathogenic. We therefore showed for the first time an association between a homozygous TBC1D31 variant with CAKUT in humans, expanding its genetic spectrum.

Combined Computational-Biochemical Approach Offers an Accelerated Path to Membrane Protein Solubilization

Membrane proteins are difficult to isolate and purify due to their dependence on the surrounding lipid membrane for structural stability. Detergents are often used to solubilize these proteins, with this approach requiring a careful balance between protein solubilization and denaturation. Determining which detergent is most appropriate for a given protein has largely been done empirically through screening, which requires large amounts of membrane protein and associated resources. Here, we describe an alternative to conventional detergent screening using a computational modeling approach to identify the most likely candidate detergents for solubilizing a protein of interest. We demonstrate our approach using ghrelin O-acyltransferase (GOAT), a member of the membrane-bound O-acyltransferase family of integral membrane enzymes that has not been solubilized or purified in active form. A computationally derived GOAT structural model provides the only structural information required for this approach. Using computational analysis of detergent ability to penetrate phospholipid bilayers and stabilize the GOAT structure, a panel of common detergents were rank-ordered for their proposed ability to solubilize GOAT. The simulations were performed at all-atom resolution for a combined simulation time of 24 μs. Independently, we biologically screened these detergents for their solubilization of fluorescently tagged GOAT constructs. We found computational prediction of protein structural stabilization was the better predictor of detergent solubilization ability, but neither approach was effective for predicting detergents that would support GOAT enzymatic function. The current rapid expansion of membrane protein computational models lacking experimental structural information and our computational detergent screening approach can greatly improve the efficiency of membrane protein detergent solubilization, supporting downstream functional and structural studies.

Biophysical Correlates of Enhanced Immunogenicity of a Stabilized Variant of the Receptor Binding Domain of SARS-CoV-2

The receptor binding domain (RBD) of SARS-CoV-2 is the primary target of neutralizing antibodies. We have previously reported the design and characterization of a mammalian cell expressed RBD derivative, mRBD1-3.2, that has higher thermal stability and greatly enhanced immunogenicity relative to the wild type mRBD. The protein is highly thermotolerant and immunogenic and is being explored for use in room temperature stable Covid-19 vaccine formulations. In the current study, we have investigated the folding pathway of both WT and stabilized RBD. It was found that chemical denaturation of RBD proceeds through a stable equilibrium intermediate. Thermal and chemical denaturation is reversible, as assayed by binding to the receptor ACE2. Unusually, in its native state, RBD binds to the hydrophobic probe ANS, and enhanced ANS binding is observed for the equilibrium intermediate state. Further characterization of the folding of mRBD1-3.2, both in solution and after reconstitution of lyophilized protein stored for a month at 37 °C, revealed a higher stability represented by higher C_m, faster refolding, slower unfolding, and enhanced resistance to proteolytic cleavage relative to WT. In contrast to WT RBD, the mutant showed decreased interaction with the hydrophobic moiety linoleic acid. Collectively, these data suggest that the enhanced immunogenicity results from reduced conformational fluctuations that likely enhance in vivo half-life as well as reduce the exposure of irrelevant non-neutralizing epitopes to the immune system.

Molecular Docking and Site-Directed Mutagenesis of GH49 Family Dextranase for the Preparation of High-Degree Polymerization Isomaltooligosaccharide

The high-degree polymerization of isomaltooligosaccharide (IMO) not only effectively promotes the growth and reproduction of Bifidobacterium in the human body but also renders it resistant to rapid degradation by gastric acid and can stimulate insulin secretion. In this study, we chose the engineered strain expressed dextranase (PsDex1711) as the research model and used the AutoDock vina molecular docking technique to dock IMO4, IMO5, and IMO6 with it to obtain mutation sites, and then studied the potential effect of key amino acids in this enzyme on its hydrolysate composition and enzymatic properties by site-directed mutagenesis method. It was found that the yield of IMO4 increased significantly to 62.32% by the mutant enzyme H373A. Saturation mutation depicted that the yield of IMO4 increased to 69.81% by the mutant enzyme H373R, and its neighboring site S374R IMO4 yield was augmented to 64.31%. Analysis of the enzymatic properties of the mutant enzyme revealed that the optimum temperature of H373R decreased from 30 °C to 20 °C, and more than 70% of the enzyme activity was maintained under alkaline conditions. The double-site saturation mutation results showed that the mutant enzyme H373R/N445Y IMO4 yield increased to 68.57%. The results suggest that the 373 sites with basic non-polar amino acids, such as arginine and histidine, affect the catalytic properties of the enzyme. The findings provide an important theoretical basis for the future marketable production of IMO4 and analysis of the structure of dextranase.

Ter-Seq: A high-throughput method to stabilize transient ternary complexes and measure associated kinetics

Regulation of biological processes by proteins often involves the formation of transient, multimeric complexes whose characterization is mechanistically important but challenging. The bacterial toxin CcdB binds and poisons DNA Gyrase. The corresponding antitoxin CcdA extracts CcdB from its complex with Gyrase through the formation of a transient ternary complex, thus rejuvenating Gyrase. We describe a high throughput methodology called Ter-Seq to stabilize probable ternary complexes and measure associated kinetics using the CcdA-CcdB-GyrA14 ternary complex as a model system. The method involves screening a yeast surface display (YSD) saturation mutagenesis library of one partner (CcdB) for mutants that show enhanced ternary complex formation. We also isolated CcdB mutants that were either resistant or sensitive to rejuvenation, and used surface plasmon resonance (SPR) with purified proteins to validate the kinetics measured using the surface display. Positions, where CcdB mutations lead to slower rejuvenation rates, are largely involved in CcdA-binding, though there were several notable exceptions suggesting allostery. Mutations at these positions reduce the affinity towards CcdA, thereby slowing down the rejuvenation process. Mutations at GyrA14-interacting positions significantly enhanced rejuvenation rates, either due to reduced affinity or complete loss of CcdB binding to GyrA14. We examined the effect of different parameters (CcdA affinity, GyrA14 affinity, surface accessibilities, evolutionary conservation) on the rate of rejuvenation. Finally, we further validated the Ter-Seq results by monitoring the kinetics of ternary complex formation for individual CcdB mutants in solution by fluorescence resonance energy transfer (FRET) studies.

Combining cysteine scanning with chemical labeling to map protein-protein interactions and infer bound structure in an intrinsically disordered region

The Mycobacterium tuberculosis genome harbours nine toxin-antitoxin (TA) systems of the mazEF family. These consist of two proteins, a toxin and an antitoxin, encoded in an operon. While the toxin has a conserved fold, the antitoxins are structurally diverse and the toxin binding region is typically intrinsically disordered before binding. We describe high throughput methodology for accurate mapping of interfacial residues and apply it to three MazEF complexes. The method involves screening one partner protein against a panel of chemically masked single cysteine mutants of its interacting partner, displayed on the surface of yeast cells. Such libraries have much lower diversity than those generated by saturation mutagenesis, simplifying library generation and data analysis. Further, because of the steric bulk of the masking reagent, labeling of virtually all exposed epitope residues should result in loss of binding, and buried residues are inaccessible to the labeling reagent. The binding residues are deciphered by probing the loss of binding to the labeled cognate partner by flow cytometry. Using this methodology, we have identified the interfacial residues for MazEF3, MazEF6 and MazEF9 TA systems of M. tuberculosis. In the case of MazEF9, where a crystal structure was available, there was excellent agreement between our predictions and the crystal structure, superior to those with AlphaFold2. We also report detailed biophysical characterization of the MazEF3 and MazEF9 TA systems and measured the relative affinities between cognate and non-cognate toxin–antitoxin partners in order to probe possible cross-talk between these systems.

Mechanistic insights into global suppressors of protein folding defects

Most amino acid substitutions in a protein either lead to partial loss-of-function or are near neutral. Several studies have shown the existence of second-site mutations that can rescue defects caused by diverse loss-of-function mutations. Such global suppressor mutations are key drivers of protein evolution. However, the mechanisms responsible for such suppression remain poorly understood. To address this, we characterized multiple suppressor mutations both in isolation and in combination with inactive mutants. We examined six global suppressors of the bacterial toxin CcdB, the known M182T global suppressor of TEM-1 β-lactamase, the N239Y global suppressor of p53-DBD and three suppressors of the SARS-CoV-2 spike Receptor Binding Domain. When coupled to inactive mutants, they promote increased in-vivo solubilities as well as regain-of-function phenotypes. In the case of CcdB, where novel suppressors were isolated, we determined the crystal structures of three such suppressors to obtain insight into the specific molecular interactions responsible for the observed effects. While most individual suppressors result in small stability enhancements relative to wildtype, which can be combined to yield significant stability increments, thermodynamic stabilisation is neither necessary nor sufficient for suppressor action. Instead, in diverse systems, we observe that individual global suppressors greatly enhance the foldability of buried site mutants, primarily through increase in refolding rate parameters measured in vitro. In the crowded intracellular environment, mutations that slow down folding likely facilitate off-pathway aggregation. We suggest that suppressor mutations that accelerate refolding can counteract this, enhancing the yield of properly folded, functional protein in vivo.

Protposer: the web server that readily proposes protein stabilizing mutations with high PPV

Protein stability is a requisite for most biotechnological and medical applications of proteins. As natural proteins tend to suffer from a low conformational stability ex vivo, great efforts have been devoted toward increasing their stability through rational design and engineering of appropriate mutations. Unfortunately, even the best currently used predictors fail to compute the stability of protein variants with sufficient accuracy and their usefulness as tools to guide the rational stabilisation of proteins is limited. We present here Protposer, a protein stabilising tool based on a different approach. Instead of quantifying changes in stability, Protposer uses structure- and sequence-based screening modules to nominate candidate mutations for subsequent evaluation by a logistic regression model, carefully trained to avoid overfitting. Thus, Protposer analyses PDB files in search for stabilization opportunities and provides a ranked list of promising mutations with their estimated success rates (eSR), their probabilities of being stabilising by at least 0.5 kcal/mol. The agreement between eSRs and actual positive predictive values (PPV) on external datasets of mutations is excellent. When Protposer is used with its Optimal kappa selection threshold, its PPV is above 0.7. Even with less stringent thresholds, Protposer largely outperforms FoldX, Rosetta and PoPMusiC. Indicating the PDB file of the protein suffices to obtain a ranked list of mutations, their eSRs and hints on the likely source of the stabilization expected. Protposer is a distinct, straightforward and highly successful tool to design protein stabilising mutations, and it is freely available for academic use at http://webapps.bifi.es/the-protposer.