A rigid network of long-range contacts increases thermostability in a mutant endoglucanase

Genetic heterogeneity and respiratory chain enzyme analysis in pediatric Indian patients with mitochondrial disorder: Report of novel variants in POLG1 gene and their functional implication using molecular dynamic simulation

Mitochondrial disorders are a heterogeneous group of disorders caused by mutations in the mitochondrial DNA or in nuclear genes encoding the mitochondrial proteins and subunits. Polymerase Gamma (POLG) is a nuclear gene and mutation in the POLG gene are one of the major causes of inherited mitochondrial disorders. In this study, 15 pediatric patients, with a wide spectrum of clinical phenotypes were screened using blood samples (n = 15) and muscle samples (n = 4). Respiratory chain enzyme analysis in the muscle samples revealed multi-complex deficiencies with Complex I deficiency present in (1/4) patients, Complex II (2/4), Complex III (3/4) and Complex IV (2/4) patients. Multiple large deletions were observed in 4/15 patients using LR-PCR. Whole exome sequencing (WES) revealed a compound heterozygous mutation consisting of a POLG1 novel variant (NP_002684.1:p.Trp261X) and a missense variant (NP_002684.1:p. Leu304Arg) in one patient and another patient harboring a novel homozygous POLG1 variant (NP_002684.1:p. Phe750Val). These variants (NP_002684.1:p. Leu304Arg) and (NP_002684.1:p. Phe750Val) and their interactions with DNA were modelled using molecular docking and molecular dynamics (MD) simulation studies. The protein conformation was analyzed as root mean square deviation (RMSD), root mean square fluctuation (RMSF) which showed local fluctuations in the mutants compared to the wildtype. However, Solvent Accessible Surface Area (SASA) significantly increased for NP_002684.1:p.Leu304Arg and decreased in NP_002684.1:p.Phe750Val mutants. Further, Contact Order analysis indicated that the Aromatic-sulfur interactions were destabilizing in the mutants. Overall, these in-silico analysis has revealed a destabilizing mutations suggesting pathogenic variants in POLG1 gene.

Computational analysis of protein stability and allosteric interaction networks in distinct conformational forms of the SARS-CoV-2 spike D614G mutant: reconciling functional mechanisms through allosteric model of spike regulation

In this study, we used an integrative computational approach to examine molecular mechanisms underlying functional effects of the D614G mutation by exploring atomistic modeling of the SARS-CoV-2 spike proteins as allosteric regulatory machines. We combined coarse-grained simulations, protein stability and dynamic fluctuation communication analysis with network-based community analysis to examine structures of the native and mutant SARS-CoV-2 spike proteins in different functional states. Through distance fluctuations communication analysis, we probed stability and allosteric communication propensities of protein residues in the native and mutant SARS-CoV-2 spike proteins, providing evidence that the D614G mutation can enhance long-range signaling of the allosteric spike engine. By combining functional dynamics analysis and ensemble-based alanine scanning of the SARS-CoV-2 spike proteins we found that the D614G mutation can improve stability of the spike protein in both closed and open forms, but shifting thermodynamic preferences towards the open mutant form. Our results revealed that the D614G mutation can promote the increased number of stable communities and allosteric hub centers in the open form by reorganizing and enhancing the stability of the S1-S2 inter-domain interactions and restricting mobility of the S1 regions. This study provides atomistic-based view of allosteric communications in the SARS-CoV-2 spike proteins, suggesting that the D614G mutation can exert its primary effect through allosterically induced changes on stability and communications in the residue interaction networks.Communicated by Ramaswamy H. Sarma.

A low CO2-responsive mutant of Setaria viridis reveals that reduced carbonic anhydrase limits C4 photosynthesis

In C4 species, β-carbonic anhydrase (CA), localized to the cytosol of the mesophyll cells, accelerates the interconversion of CO2 to HCO3-, the substrate used by phosphoenolpyruvate carboxylase (PEPC) in the first step of C4 photosynthesis. Here we describe the identification and characterization of low CO2-responsive mutant 1 (lcr1) isolated from an N-nitroso-N-methylurea- (NMU) treated Setaria viridis mutant population. Forward genetic investigation revealed that the mutated gene Sevir.5G247800 of lcr1 possessed a single nucleotide transition from cytosine to thymine in a β-CA gene causing an amino acid change from leucine to phenylalanine. This resulted in severe reduction in growth and photosynthesis in the mutant. Both the CO2 compensation point and carbon isotope discrimination values of the mutant were significantly increased. Growth of the mutants was stunted when grown under ambient pCO2 but recovered at elevated pCO2. Further bioinformatics analyses revealed that the mutation has led to functional changes in one of the conserved residues of the protein, situated near the catalytic site. CA transcript accumulation in the mutant was 80% lower, CA protein accumulation 30% lower, and CA activity ~98% lower compared with the wild type. Changes in the abundance of other primary C4 pathway enzymes were observed; accumulation of PEPC protein was significantly increased and accumulation of malate dehydrogenase and malic enzyme decreased. The reduction of CA protein activity and abundance in lcr1 restricts the supply of bicarbonate to PEPC, limiting C4 photosynthesis and growth. This study establishes Sevir.5G247800 as the major CA allele in Setaria for C4 photosynthesis and provides important insights into the function of CA in C4 photosynthesis that would be required to generate a rice plant with a functional C4 biochemical pathway.

Genome-Wide Mutagenesis of Hepatitis C Virus Reveals Ability of Genome To Overcome Detrimental Mutations

To gain insight into the impact of mutations on the viability of the hepatitis C virus (HCV) genome, we created a set of full-genome mutant libraries, differing from the parent sequence as well as each other, by using a random mutagenesis approach; the proportion of mutations increased across these libraries with declining template amount or dATP concentration. The replication efficiencies of full-genome mutant libraries ranged between 71 and 329 focus-forming units (FFU) per 10⁵ Huh7.5 cells. Mutant libraries with low proportions of mutations demonstrated low replication capabilities, whereas those with high proportions of mutations had their replication capabilities restored. Hepatoma cells transfected with selected mutant libraries, with low (4 mutations per 10,000 bp copied), moderate (33 mutations), and high (66 mutations) proportions of mutations, and their progeny were subjected to serial passage. Predominant virus variants (mutants) from these mutant libraries (Mutant^l, Mutant^m, and Mutant^h, respectively) were evaluated for changes in growth kinetics and particle-to-FFU unit ratio, virus protein expression, and modulation of host cell protein synthesis. Mutant^m and Mutant^l variants produced >3.0-log-higher extracellular progeny per ml than the parent, and Mutant^h produced progeny at a rate 1.0-log lower. More than 80% of the mutations were in a nonstructural part of the mutant genomes, the majority were nonsynonymous, and a moderate to large proportion were in the conserved regions. Our results suggest that the HCV genome has the ability to overcome lethal/deleterious mutations because of the high reproduction rate but highly selects for random, beneficial mutations.IMPORTANCE Hepatitis C virus (HCV) in vivo displays high genetic heterogeneity, which is partly due to the high reproduction and random substitutions during error-prone genome replication. It is difficult to introduce random substitutions in vitro because of limitations in inducing mutagenesis from the 5' end to the 3' end of the genome. Our study has overcome this limitation. We synthesized full-length genomes with few to several random mutations in the background of an HCV clone that can recapitulate all steps of the life cycle. Our study provides evidence of the capability of the HCV genome to overcome deleterious mutations and remain viable. Mutants that emerged from the libraries had diverse phenotype profiles compared to the parent, and putative adaptive mutations mapped to segments of the conserved nonstructural genome. We demonstrate the potential utility of our system for the study of sequence variation that ensures the survival and adaptation of HCV.

Recyclable magnetic carboxymethyl chitosan/calcium alginate - cellulase bioconjugates for corn stalk hydrolysis

The use of cellulase hydrolysis of straw to produce fermentable sugars has many application prospects. However, cellulase is very expensive, which hampers its industrial applications. To improve cellulase's catalytic activity and reduce the enzyme cost, magnetite carboxymethyl chitosan/calcium alginate - cellulase bioconjugate (MCCCB) was synthesized via an improved hydrothermal method, molecular self-assembly technology, physical absorption, embedding and covalent bonding. Its loading capacity was 3.95mg/mL, and the catalytic activity increased to 267.18%. We decreased the release rate, improved the reusability, and enhanced the stability of MCCCB. Corn stalk hydrolysis also greatly improved, and the overall yield of fermentable sugars increased by 698.26%. All of these results indicate that MCCCB could significantly improve the efficiency of cellulase, greatly reduce the cost of enzyme, and effectively promote the production of fermentable sugars.

Application of Rigidity Theory to the Thermostabilization of Lipase A from Bacillus subtilis

Protein thermostability is a crucial factor for biotechnological enzyme applications. Protein engineering studies aimed at improving thermostability have successfully applied both directed evolution and rational design. However, for rational approaches, the major challenge remains the prediction of mutation sites and optimal amino acid substitutions. Recently, we showed that such mutation sites can be identified as structural weak spots by rigidity theory-based thermal unfolding simulations of proteins. Here, we describe and validate a unique, ensemble-based, yet highly efficient strategy to predict optimal amino acid substitutions at structural weak spots for improving a protein’s thermostability. For this, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state has been found as the cause for thermostability. When applied prospectively to lipase A from Bacillus subtilis, we achieved both a high success rate (25% over all experimentally tested mutations, which raises to 60% if small-to-large residue mutations and mutations in the active site are excluded) in predicting significantly thermostabilized lipase variants and a remarkably large increase in those variants’ thermostability (up to 6.6°C) based on single amino acid mutations. When considering negative controls in addition and evaluating the performance of our approach as a binary classifier, the accuracy is 63% and increases to 83% if small-to-large residue mutations and mutations in the active site are excluded. The gain in precision (predictive value for increased thermostability) over random classification is 1.6-fold (2.4-fold). Furthermore, an increase in thermostability predicted by our approach significantly points to increased experimental thermostability (p < 0.05). These results suggest that our strategy is a valuable complement to existing methods for rational protein design aimed at improving thermostability.

Protein thermostability is a crucial factor for biotechnological enzyme applications. However, performance studies of computational approaches for predicting effects of mutations on protein (thermo)stability have suggested that there is still room for improvement. We describe and validate a novel and unique strategy to predict optimal amino acid substitutions at structural weak spots. At variance with other rational approaches, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state is the underlying cause for thermostability. When applied prospectively on lipase LipA from Bacillus subtilis, a high success rate in predicting thermostabilized lipase variants and a remarkably large increase in their thermostability is achieved. This demonstrates the value of the novel strategy, which extends the existing portfolio of methods for rational protein design.

Polarity Alteration of a Calcium Site Induces a Hydrophobic Interaction Network and Enhances Cel9A Endoglucanase Thermostability

Structural calcium sites control protein thermostability and activity by stabilizing native folds and changing local conformations. Alicyclobacillus acidocaldarius survives in thermal-acidic conditions and produces an endoglucanase Cel9A (AaCel9A) which contains a calcium-binding site (Ser465 to Val470) near the catalytic cleft. By superimposing the Ca(2+)-free and Ca(2+)-bounded conformations of the calcium site, we found that Ca(2+) induces hydrophobic interactions between the calcium site and its nearby region by driving a conformational change. The hydrophobic interactions at the high-B-factor region could be enhanced further by replacing the surrounding polar residues with hydrophobic residues to affect enzyme thermostability and activity. Therefore, the calcium-binding residue Asp468 (whose side chain directly ligates Ca(2+)), Asp469, and Asp471 of AaCel9A were separately replaced by alanine and valine. Mutants D468A and D468V showed increased activity compared with those of the wild type with 0 mM or 10 mM Ca(2+) added, whereas the Asp469 or Asp471 substitution resulted in decreased activity. The D468A crystal structure revealed that mutation D468A triggered a conformational change similar to that induced by Ca(2+) in the wild type and developed a hydrophobic interaction network between the calcium site and the neighboring hydrophobic region (Ala113 to Ala117). Mutations D468V and D468A increased 4.5°C and 5.9°C, respectively, in melting temperature, and enzyme half-life at 75°C increased approximately 13 times. Structural comparisons between AaCel9A and other endoglucanases of the GH9 family suggested that the stability of the regions corresponding to the AaCel9A calcium site plays an important role in GH9 endoglucanase catalysis at high temperature.

Structural Rigidity and Protein Thermostability in Variants of Lipase A from Bacillus subtilis

Understanding the origin of thermostability is of fundamental importance in protein biochemistry. Opposing views on increased or decreased structural rigidity of the folded state have been put forward in this context. They have been related to differences in the temporal resolution of experiments and computations that probe atomic mobility. Here, we find a significant (p = 0.004) and fair (R² = 0.46) correlation between the structural rigidity of a well-characterized set of 16 mutants of lipase A from Bacillus subtilis (BsLipA) and their thermodynamic thermostability. We apply the rigidity theory-based Constraint Network Analysis (CNA) approach, analyzing directly and in a time-independent manner the statics of the BsLipA mutants. We carefully validate the CNA results on macroscopic and microscopic experimental observables and probe for their sensitivity with respect to input structures. Furthermore, we introduce a robust, local stability measure for predicting thermodynamic thermostability. Our results complement work that showed for pairs of homologous proteins that raising the structural stability is the most common way to obtain a higher thermostability. Furthermore, they demonstrate that related series of mutants with only a small number of mutations can be successfully analyzed by CNA, which suggests that CNA can be applied prospectively in rational protein design aimed at higher thermodynamic thermostability.

A novel mechanism of protein thermostability: a unique N-terminal domain confers heat resistance to Fe/Mn-SODs

Superoxide dismutases (SODs), especially thermostable SODs, are widely applied in medical treatments, cosmetics, food, agriculture, and other industries given their excellent antioxidant properties. A novel thermostable cambialistic SOD from Geobacillus thermodenitrificans NG80-2 exhibits maximum activity at 70°C and high thermostability over a broad range of temperatures (20–80°C). Unlike other reported SODs, this enzyme contains an extra repeat-containing N-terminal domain (NTD) of 244 residues adjacent to the conserved functional SODA domain. Deletion of the NTD dramatically decreased its optimum active temperature (OAT) to 30°C and also impaired its thermostability. Conversely, appending the NTD to a mesophilic counterpart from Bacillus subtilis led to a moderately thermophilic enzyme (OAT changed from 30 to 55°C) with improved heat resistance. Temperature-dependant circular dichroism analysis revealed the enhanced conformational stability of SODs fused with this NTD. Furthermore, the NTD also contributes to the stress resistance of host proteins without altering their metal ion specificity or oligomerisation form except for a slight effect on their pH profile. We therefore demonstrate that the NTD confers outstanding thermostability to the host protein. To our knowledge, this is the first discovery of a peptide capable of remarkably improving protein thermostability and provides a novel strategy for bioengineering thermostable SODs.

Fast closure of long loops at the initiation of the folding transition of globular proteins studied by time-resolved FRET-based methods

The protein folding problem would be considered “solved” when it will be possible to “read genes”, i.e., to predict the native fold of proteins, their dynamics, and the mechanism of fast folding based solely on sequence data. The long-term goal should be the creation of an algorithm that would simulate the stepwise mechanism of folding, which constrains the conformational space and in which random search for stable interactions is possible. Here, we focus attention on the initial phases of the folding transition starting with the compact disordered collapsed ensemble, in search of the initial sub-domain structural biases that direct the otherwise stochastic dynamics of the backbone. Our studies are designed to test the “loop hypothesis”, which suggests that fast closure of long loop structures by non-local interactions between clusters of mainly non-polar residues is an essential conformational step at the initiation of the folding transition of globular proteins. We developed and applied experimental methods based on time-resolved resonance excitation energy transfer (trFRET) measurements combined with fast mixing methods and studied the initial phases of the folding of

Endoglucanases: insights into thermostability for biofuel applications

Obtaining bioethanol from cellulosic biomass involves numerous steps, among which the enzymatic conversion of the polymer to individual sugar units has been a main focus of the biotechnology industry. Among the cellulases that break down the polymeric cellulose are endoglucanases that act synergistically for subsequent hydrolytic reactions. The endoglucanases that have garnered relatively more attention are those that can withstand high temperatures, i.e., are thermostable. Although our understanding of thermostability in endoglucanases is incomplete, some molecular features that are responsible for increased thermostability have been recently identified. This review focuses on the investigations of endoglucanases and their implications for biofuel applications.