Shortening a loop can increase protein native state entropy

Insertion-Deletion Events Are Depleted in Protein Regions with Predicted Secondary Structure

A fundamental goal in evolutionary biology and population genetics is to understand how selection shapes the fate of new mutations. Here, we test the null hypothesis that insertion-deletion (indel) events in protein-coding regions occur randomly with respect to secondary structures. We identified indels across 11,444 sequence alignments in mouse, rat, human, chimp, and dog genomes and then quantified their overlap with four different types of secondary structure-alpha helices, beta strands, protein bends, and protein turns-predicted by deep-learning methods of AlphaFold2. Indels overlapped secondary structures 54% as much as expected and were especially underrepresented over beta strands, which tend to form internal, stable regions of proteins. In contrast, indels were enriched by 155% over regions without any predicted secondary structures. These skews were stronger in the rodent lineages compared to the primate lineages, consistent with population genetic theory predicting that natural selection will be more efficient in species with larger effective population sizes. Nonsynonymous substitutions were also less common in regions of protein secondary structure, although not as strongly reduced as in indels. In a complementary analysis of thousands of human genomes, we showed that indels overlapping secondary structure segregated at significantly lower frequency than indels outside of secondary structure. Taken together, our study shows that indels are selected against if they overlap secondary structure, presumably because they disrupt the tertiary structure and function of a protein.

New robust subtilisins from halotolerant and halophilic Bacillaceae

The aim of the present study was the characterisation of three true subtilisins and one phylogenetically intermediate subtilisin from halotolerant and halophilic microorganisms. Considering the currently growing enzyme market for efficient and novel biocatalysts, data mining is a promising source for novel, as yet uncharacterised enzymes, especially from halophilic or halotolerant Bacillaceae, which offer great potential to meet industrial needs. Both halophilic bacteria Pontibacillus marinus DSM 16465^T and Alkalibacillus haloalkaliphilus DSM 5271^T and both halotolerant bacteria Metabacillus indicus DSM 16189 and Litchfieldia alkalitelluris DSM 16976^T served as a source for the four new subtilisins SPPM, SPAH, SPMI and SPLA. The protease genes were cloned and expressed in Bacillus subtilis DB104. Purification to apparent homogeneity was achieved by ethanol precipitation, desalting and ion-exchange chromatography. Enzyme activity could be observed between pH 5.0-12.0 with an optimum for SPPM, SPMI and SPLA around pH 9.0 and for SPAH at pH 10.0. The optimal temperature for SPMI and SPLA was 70 °C and for SPPM and SPAH 55 °C and 50 °C, respectively. All proteases showed high stability towards 5% (w/v) SDS and were active even at NaCl concentrations of 5 M. The four proteases demonstrate potential for future biotechnological applications. KEY POINTS: • Halophilic and halotolerant Bacillaceae are a valuable source of new subtilisins. • Four new subtilisins were biochemically characterised in detail. • The four proteases show potential for future biotechnological applications.

Structural Dynamics Predominantly Determine the Adaptability of Proteins to Amino Acid Deletions

The insertion or deletion (indel) of amino acids has a variety of effects on protein function, ranging from disease-forming changes to gaining new functions. Despite their importance, indels have not been systematically characterized towards protein engineering or modification goals. In the present work, we focus on deletions composed of multiple contiguous amino acids (mAA-dels) and their effects on the protein (mutant) folding ability. Our analysis reveals that the mutant retains the native fold when the mAA-del obeys well-defined structural dynamics properties: localization in intrinsically flexible regions, showing low resistance to mechanical stress, and separation from allosteric signaling paths. Motivated by the possibility of distinguishing the features that underlie the adaptability of proteins to mAA-dels, and by the rapid evaluation of these features using elastic network models, we developed a positive-unlabeled learning-based classifier that can be adopted for protein design purposes. Trained on a consolidated set of features, including those reflecting the intrinsic dynamics of the regions where the mAA-dels occur, the new classifier yields a high recall of 84.3% for identifying mAA-dels that are stably tolerated by the protein. The comparative examination of the relative contribution of different features to the prediction reveals the dominant role of structural dynamics in enabling the adaptation of the mutant to mAA-del without disrupting the native fold.

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.

Phylogenetic survey of the subtilase family and a data-mining-based search for new subtilisins from Bacillaceae

The subtilase family (S8), a member of the clan SB of serine proteases are ubiquitous in all kingdoms of life and fulfil different physiological functions. Subtilases are divided in several groups and especially subtilisins are of interest as they are used in various industrial sectors. Therefore, we searched for new subtilisin sequences of the family Bacillaceae using a data mining approach. The obtained 1,400 sequences were phylogenetically classified in the context of the subtilase family. This required an updated comprehensive overview of the different groups within this family. To fill this gap, we conducted a phylogenetic survey of the S8 family with characterised holotypes derived from the MEROPS database. The analysis revealed the presence of eight previously uncharacterised groups and 13 subgroups within the S8 family. The sequences that emerged from the data mining with the set filter parameters were mainly assigned to the subtilisin subgroups of true subtilisins, high-alkaline subtilisins, and phylogenetically intermediate subtilisins and represent an excellent source for new subtilisin candidates.

Estimating Change in Foldability Due to Multipoint Deletions in Protein Structures

Insertions/deletions of amino acids in the protein backbone potentially result in altered structural/functional specifications. They can either contribute positively to the evolutionary process or can result in disease conditions. Despite being the second most prevalent form of protein modification, there are no databases or computational frameworks that delineate harmful multipoint deletions (MPD) from beneficial ones. We introduce a positive unlabeled learning-based prediction framework (PROFOUND) that utilizes fold-level attributes, environment-specific properties, and deletion site-specific properties to predict the change in foldability arising from such MPDs, both in the non-loop and loop regions of protein structures. In the absence of any protein structure dataset to study MPDs, we introduce a dataset with 153 MPD instances that lead to native-like folded structures and 7650 unlabeled MPD instances whose effect on the foldability of the corresponding proteins is unknown. PROFOUND on 10-fold cross-validation on our newly introduced dataset reports a recall of 82.2% (86.6%) and a fall out rate (FR) of 14.2% (20.6%), corresponding to MPDs in the protein loop (non-loop) region. The low FR suggests that the foldability in proteins subject to MPDs is not random and necessitates unique specifications of the deleted region. In addition, we find that additional evolutionary attributes contribute to higher recall and lower FR. The first of a kind foldability prediction system owing to MPD instances and the newly introduced dataset will potentially aid in novel protein engineering endeavors.

Engineering a Trypsin-Resistant Thermophilic α-Galactosidase to Enhance Pepsin Resistance, Acidic Tolerance, Catalytic Performance, and Potential in the Food and Feed Industry

α-Galactosidase has potential applications, and attempts to improve proteolytic resistance of enzymes have important values. We use a novel strategy for genetic manipulation of a pepsin-sensitive region specific for a pepsin-sensitive but trypsin-resistant high-temperature-active Gal27B from Neosartorya fischeri to screen mutants with enhanced pepsin resistance. All enzymes were produced in Pichia pastoris to identify the roles of loop 4 (Gal27B-A23) and its key residue at position 156 (Gly156Arg/Pro/His) in pepsin resistance. Gal27B-A23 and Gly156Arg/Pro/His elevated pepsin resistance, thermostability, stability at low pH, activity toward raffinose (5.3-6.9-fold) and stachyose (about 1.3-fold), and catalytic efficiencies (up to 4.9-fold). Replacing the pepsin cleavage site Glu155 with Gly improved pepsin resistance but had no effect on pepsin resistance when Arg/Pro/His was at position 156. Thus, pepsin resistance could appear to occur through steric hindrance between the residue at the altered site and neighboring pepsin active site. In the presence of pepsin or trypsin, all mutations increased the ability of Gal27B to hydrolyze galactosaccharides in soybean flour (up to 9.6- and 4.3-fold, respectively) and promoted apparent metabolizable energy and nutrient digestibility in soybean meal for broilers (1.3-1.8-fold). The high activity and tolerance to heat, low pH, and protease benefit food and feed industry in a cost-effective way.

Biochemical Characterization and Mutational Analysis of a Lactone Hydrolase from Phialophora americana

The mycotoxin zearalenone (ZEN) is a secondary metabolite produced mainly by Fusarium species. ZEN poses health hazards both for humans and animals, as a major contaminant in the food and feed industries. Currently, there is no effective technique for degrading ZEN during industrial processes. In this study, we isolated and biochemically characterized a novel lactone hydrolase, ZHD607, isolated from Phialophora americana, cloned, and exogenously expressed in Pichia pastoris. ZHD607 was characterized as a mesophilic lactone hydrolase having a neutral pH and showing optimal activity at 35 °C and pH 8.0. Two mutants, ZHDM1 and I160Y, generated from ZHD607 based on structure and sequence alignment analyses, exhibited 2.9- and 3.4-fold higher activity towards ZEN than did ZHD607. Molecular dynamics simulation revealed diverse mechanisms driving this improved catalytic activity. These findings enrich our knowledge about ZHD enzyme family and represent an important step toward industrialization of ZEN-detoxifying lactone hydrolases.

Effect of disulphide loop length on mechanochemical structural stability of macromolecules

Polymer chains folded with a single disulphide loop are shown to present distinct rates of mechanochemical fragmentation.

Finding the generalized molecular principles of protein thermal stability

Are there any generalized molecular principles of thermal adaptation? Here, integrating the concepts of structural bioinformatics, sequence analysis, and classical knot theory, we develop a robust computational framework that seeks for mechanisms of thermal adaptation by comparing orthologous mesophilic-thermophilic and mesophilic-hyperthermophilic proteins of remarkable structural and topological similarities, and still leads us to context-independent results. A comprehensive analysis of 4741 high-resolution, non-redundant X-ray crystallographic structures collected from 11 hyperthermophilic, 32 thermophilic and 53 mesophilic prokaryotes unravels at least five "nearly universal" signatures of thermal adaptation, irrespective of the enormous sequence, structure, and functional diversity of the proteins compared. A careful investigation further extracts a set of amino acid changes that can potentially enhance protein thermal stability, and remarkably, these mutations are overrepresented in protein crystallization experiments, in disorder-to-order transitions and in engineered thermostable variants of existing mesophilic proteins. These results could be helpful to find a precise, global picture of thermal adaptation.

Analyzing Change in Protein Stability Associated with Single Point Deletions in a Newly Defined Protein Structure Database

Protein backbone alternation due to insertion/deletion or mutation operation often results in a change of fundamental biophysical properties of proteins. The proposed work intends to encode the protein stability changes associated with single point deletions (SPDs) of amino acids in proteins. The encoding will help in the primary screening of detrimental backbone modifications before opting for expensive in vitro experimentations. In the absence of any benchmark database documenting SPDs, we curate a data set containing SPDs that lead to both folded conformations and unfolded state. We differentiate these SPD instances with the help of simple structural and physicochemical features and eventually classify the foldability resulting out of SPDs using a Random Forest classifier and an Elliptic Envelope based outlier detector. Adhering to leave one out cross validation, the accuracy of the Random Forest classifier and the Elliptic Envelope is of 99.4% and 98.1%, respectively. The newly defined database and the delineation of SPD instances based on its resulting foldability provide a head start toward finding a solution to the given problem.

Insight into the Thermophilic Mechanism of a Glycoside Hydrolase Family 5 β-Mannanase

To study the molecular basis for thermophilic β-mannanase of glycoside hydrolase family 5, two β-mannanases, TlMan5A and PMan5A, from Talaromyces leycettanus JCM12802 and Penicillium sp. WN1 were used as models. The four residues, His112 and Phe113, located near the antiparallel β-sheet at the barrel bottom and Leu375 and Ala408 from loop 7 and loop 8 of PMan5A, were inferred to be key thermostability contributors through module substitution, truncation, and site-directed mutagenesis. The effects of these four residues on the thermal properties followed the order H112Y > A408P > L375H > F113Y and were strongly synergetic. These results were interpreted structurally using molecular dynamics (MD) simulations, which showed that improved hydrophobic interactions in the inner wall of the β-barrel and the rigidity of loop 8 were caused by the outside domain of the barrel bottom and proline, respectively. The TIM barrel bottom and four specific residues responsible for the thermostability of GH5 β-mannanases were elucidated.

An NMR Confirmation for Increased Folded State Entropy Following Loop Truncation

Previous studies conducted on flexible loop regions in proteins revealed that the energetic consequences of changing loop length predominantly arise from the entropic cost of ordering a loop during folding. However, in an earlier study of human acylphosphatase (hmAcP) using experimental and computational approaches, we showed that thermodynamic stabilization upon loop truncation can be attributed mainly to the increased entropy of the folded state. Here, using 15N NMR spectroscopy, we studied the effect of loop truncation on hmAcP backbone dynamics on the picosecond-nanosecond timescale with the aim of confirming the effect of folded state entropy on protein stability. NMR-relaxation-derived N-H squared generalized order parameters reveal that loop truncation results in a significant increase in protein conformational flexibility. Comparison of these results with previously acquired all-atom molecular dynamics simulation, analyzed here in terms of squared generalized NMR order parameters, demonstrates general agreement between the two methods. The NMR study not only provides direct evidence for the enhanced conformational entropy of the folded state of hmAcP upon loop truncation but also gives a quantitative measure of the observed effects.

The Cost of Long Catalytic Loops in Folding and Stability of the ALS-Associated Protein SOD1

A conspicuous feature of the amyotrophic lateral sclerosis (ALS)-associated protein SOD1 is that its maturation into a functional enzyme relies on local folding of two disordered loops into a catalytic subdomain. To drive the disorder-to-order transition, the protein employs a single Zn²⁺ ion. The question is then if the entropic penalty of maintaining such disordered loops in the immature apoSOD1 monomer is large enough to explain its unusually low stability, slow folding, and pathological aggregation in ALS. To find out, we determined the effects of systematically altering the SOD1-loop lengths by protein redesign. The results show that the loops destabilize the apoSOD1 monomer by ∼3 kcal/mol, rendering the protein marginally stable and accounting for its aggregation behavior. Yet the effect on the global folding kinetics remains much smaller with a transition-state destabilization of <1 kcal/mol. Notably, this 1/3 transition-state to folded-state stability ratio provides a clear-cut example of the enigmatic disagreement between the Leffler α value from loop-length alterations (typically 1/3) and the "standard" reaction coordinates based on solvent perturbations (typically >2/3). Reconciling the issue, we demonstrate that the disagreement disappears when accounting for the progressive loop shortening that occurs along the folding pathway. The approach assumes a consistent Flory loop entropy scaling factor of c = 1.48 for both equilibrium and kinetic data and has the added benefit of verifying the tertiary interactions of the folding nucleus as determined by phi-value analysis. Thus, SOD1 not only represents a case where evolution of key catalytic function has come with the drawback of a destabilized apo state but also stands out as a well-suited model system for exploring the physicochemical details of protein self-organization.

Improving the thermostability of Serratia marcescens B4A chitinase via G191V site-directed mutagenesis

Chitinases with high thermostability are important for many industrial and biotechnological applications. This study was conducted to enhance the stability of Serratia marcescens B4A chitinase by site directed mutagenesis of G191 V. Further characterization showed that the thermal stability of the mutant showed marked increase of about 5 and 15 fold at 50 and 60 °C respectively, while the optimum temperature and pH was retained. Kinetic analysis showed decreased K_m and V_max of the mutant in comparison with the wild type chitinase of about 1.3 and 3 fold, respectively. Based on structural prediction, it was speculated that this replacement shortened an important loop concomitant with the extension of adjacent β sheets. Accordingly, a higher thermostability of G191 V up to 90 °C supporting the decreased flexibility of unfolded state was also indicated. Finally, a practical proof of kinetic and thermal stabilization of chitinase was provided through decreased flexibility and entropic stabilization of its surface loops.

Identification of the critical residues responsible for differential reactivation of the triosephosphate isomerases of two trypanosomes

The reactivation of triosephosphate isomerase (TIM) from unfolded monomers induced by guanidine hydrochloride involves different amino acids of its sequence in different stages of protein refolding. We describe a systematic mutagenesis method to find critical residues for certain physico-chemical properties of a protein. The two similar TIMs of Trypanosoma brucei and Trypanosoma cruzi have different reactivation velocities and efficiencies. We used a small number of chimeric enzymes, additive mutants and planned site-directed mutants to produce an enzyme from T. brucei with 13 mutations in its sequence, which reactivates fast and efficiently like wild-type (WT) TIM from T. cruzi, and another enzyme from T. cruzi, with 13 slightly altered mutations, which reactivated slowly and inefficiently like the WT TIM of T. brucei. Our method is a shorter alternative to random mutagenesis, saturation mutagenesis or directed evolution to find multiple amino acids critical for certain properties of proteins.

Backbone Circularization Coupled with Optimization of Connecting Segment in Effectively Improving the Stability of Granulocyte-Colony Stimulating Factor

Backbone circularization of protein is a powerful method to improve its structural stability. In this paper, we presumed that a tight connection leads to much higher stability. Therefore, we designed circularized variants of a granulocyte-colony stimulating factor (G-CSF) with a structurally optimized terminal connection. To estimate the appropriate length of the connection, we surveyed the Protein Data Bank to find local structures as a model for the connecting segment. We set the library of local structures composed of "helix-loop-helix," subsequently selected entries similar to the G-CSF terminus, and finally sorted the hit structures according to the loop length. Two, five, or nine loop residues were frequently observed; thus, three circularized variants (C163, C166, and C170) were constructed, prepared, and evaluated. All circularized variants demonstrated a higher thermal stability than linear G-CSF (L175). In particular, C166 that retained five connecting residues demonstrated apparent T_m values of 69.4 °C, which is 8.7 °C higher than that of the circularized variant with no truncation (C177), indicating that the optimization of the connecting segment is effective for enhancing the overall structural stability. C166 also showed higher proteolytic stability against both endoprotease and exopeptidase than L175. We anticipate that the present study will contribute to the improvement in the general design of circularized protein and development of G-CSF biobetters.

Physical and molecular bases of protein thermal stability and cold adaptation

The molecular bases of thermal and cold stability and adaptation, which allow proteins to remain folded and functional in the temperature ranges in which their host organisms live and grow, are still only partially elucidated. Indeed, both experimental and computational studies fail to yield a fully precise and global physical picture, essentially because all effects are context-dependent and thus quite intricate to unravel. We present a snapshot of the current state of knowledge of this highly complex and challenging issue, whose resolution would enable large-scale rational protein design.