Insights into thermal stability of thermophilic nitrile hydratases by molecular dynamics simulation

Unveiling potential PET degrading eukaryotes through in silico bioprospecting of PETases

This study addresses the environmental problem of PET plastic through in silico bioprospecting for the identification and experimental validation of novel PET degrading eukaryotes through the in silico bioprospectingI of PETases, employing a methodology that combines Hidden Markov Models (HMMs), clustering techniques, molecular docking, and dynamic simulations. A total of 424 putative PETase sequences were identified from 219 eukaryotic organisms, highlighting six sequences with low affinity energies. The Aspergillus luchuensis sequence showed the lowest Gibbs free energy and exhibited stability at different temperatures in molecular dynamics assays. Experimental validation, through a plate clearance assay and HPLC, confirmed PETase activity in three wild-type fungal strains, with A. luchuensis showing the highest efficiency. The results obtained demonstrate the effectiveness of combining computational and experimental approaches as proof of concept to discover and validate eukaryotes with PET-degrading capabilities opening new perspectives for the sustainable management of this type of waste and contributing to its environmental mitigation.

Design of potential anti-cancer agents as COX-2 inhibitors, using 3D-QSAR modeling, molecular docking, oral bioavailability proprieties, and molecular dynamics simulation

Modeling the structural properties of novel morpholine-bearing 1, 5-diaryl-diazole derivatives as potent COX-2 inhibitor, two proposed models based on CoMFA and CoMSIA were evaluated by external and internal validation methods. Partial least squares analysis produced statistically significant models with Q 2 values of 0.668 and 0.652 for CoMFA and CoMSIA, respectively, and also a significant non-validated correlation coefficient R² with values of 0.882 and 0.878 for CoMFA and CoMSIA, respectively. Both models met the requirements of Golbraikh and Tropsha, which means that both models are consistent with all validation techniques. Analysis of the CoMFA and CoMSIA contribution maps and molecular docking revealed that the R1 substituent has a very significant effect on their biological activity. The most active molecules were evaluated for their thermodynamic stability by performing MD simulations for 100 ns; it was revealed that the designed macromolecular ligand complex with 3LN1 protein exhibits a high degree of structural and conformational stability. Based on these results, we predicted newly designed compounds, which have acceptable oral bioavailability properties and would have high synthetic accessibility.

Enhancing the thermostability of Streptomyces cyaneofuscatus strain Ms1 tyrosinase by multi-factors rational design and molecular dynamics simulations

This study presents a multi-factor rational design strategy combined with molecular dynamics simulation to improve the thermostability of Streptomyces cyaneofuscatus strain Ms1 tyrosinase. Candidate mutation sites were identified using Discovery Studio and FoldX software, and the double mutant G124W/G137W was obtained. The mutant was heterogeneously expressed in Escherichia coli strain Rosetta2 (DE3), and its thermostability was verified. Results indicate that the rational design method, combined with molecular dynamics simulation and protein energy calculation, improved the enzyme's thermostability more accurately and effectively. The double mutant G124W/G137W had an optimum temperature of 60°C, about 5.0°C higher than that of the wild-type TYRwt, and its activity was 171.06% higher than the wild-type TYRwt. Its thermostability was enhanced, 42.78% higher than the wild-type at 50°C. These findings suggest that the rational design strategy applied in this study can facilitate the application of industrial enzymes in the pharmaceutical industry.

Binding of Anionic Polyacrylamide with Amidase and Laccase under 298, 303, and 308 K: Docking and Molecular Dynamics Simulation Studies Combined with Experiments

Amidase and laccase play a key role in the degradation process of anionic polyacrylamide (HPAM). However, the largest challenge of HPAM enzymatic degradation is whether the enzyme can bind with a substrate for a period of time. Here, the most suitable complexes, namely, Rh Amidase-HPAM-2 and Bacillus subtilis (B. subtilis) laccase-HPAM-3, were obtained by docking, and they were carried out for molecular dynamics simulation (MDS) under 298, 303, and 308 K. MDS result analysis showed that Rh Amidase-HPAM-2 was the most stable at 298 K mainly due to a salt bridge and a hydrogen bond, and B. subtilis laccase-HPAM-3 was the most stable at 298 K mainly due to two electrostatic and hydrogen bonds. The LYS96 in Rh Amidase-HPAM-2 and LYS135 in B. subtilis laccase-HPAM-3 had been the most important in their binding process. The binding of Rh Amidase-HPAM-2 and B. subtilis laccase-HPAM-3 was optimal at 303 and 298 K, respectively. HPAM was degraded by mixed bacteria, and the optimal conditions were determined to be 308 K, initial pH = 7, and an inoculated dosage of 2 mL. Under these conditions, the degradation ratio reached 39.24%. The effect of parameters on the HPAM degradation ratio followed a decreasing order of temperature > initial pH > inoculated dosage. The HPAM codegradation mechanism was supposed by mixed bacteria according to test data. The mixed bacteria secreted both amidase and laccase, and they interacted jointly with HPAM. These results lay a theoretical foundation to design and modify the enzyme through mutation experiments in the future.

Engineering enhanced thermostability into the Geobacillus pallidus nitrile hydratase

Nitrile hydratases (NHases) are important biocatalysts for the enzymatic conversion of nitriles to industrially-important amides such as acrylamide and nicotinamide. Although thermostability in this enzyme class is generally low, there is not sufficient understanding of its basis for rational enzyme design. The gene expressing the Co-type NHase from the moderate thermophile, Geobacillus pallidus RAPc8 (NRRL B-59396), was subjected to random mutagenesis. Four mutants were selected that were 3 to 15-fold more thermostable than the wild-type NHase, resulting in a 3.4–7.6 ​kJ/mol increase in the activation energy of thermal inactivation at 63 ​°C. High resolution X-ray crystal structures (1.15–1.80 ​Å) were obtained of the wild-type and four mutant enzymes. Mutant 9E, with a resolution of 1.15 ​Å, is the highest resolution crystal structure obtained for a nitrile hydratase to date. Structural comparisons between the wild-type and mutant enzymes illustrated the importance of salt bridges and hydrogen bonds in enhancing NHase thermostability. These additional interactions variously improved thermostability by increased intra- and inter-subunit interactions, preventing cooperative unfolding of α-helices and stabilising loop regions. Some hydrogen bonds were mediated via a water molecule, specifically highlighting the significance of structured water molecules in protein thermostability. Although knowledge of the mutant structures makes it possible to rationalize their behaviour, it would have been challenging to predict in advance that these mutants would be stabilising.

•

Random mutagenesis yields a 15-fold increase in nitrile hydratase thermostability.

•

Salt bridges and hydrogen bonds improves nitrile hydratase thermostability.

•

Water-mediated hydrogen bonds improves protein thermostability.

Discovery of Novel Cyclic Salt Bridge in Thermophilic Bacterial Protease and Study of its Sequence and Structure

The plausible explanation behind the stability of thermophilic protein is still yet to be defined more clearly. Here, an in silico study has been undertaken by investigating the sequence and structure of protease from thermophilic (tPro) bacteria and mesophilic (mPro) bacteria. Results showed that charged and uncharged polar residues have higher abundance in tPro. In extreme environment, the tPro is stabilized by high number of isolated and network salt bridges. A novel cyclic salt bridge is also found in a structure of tPro. High number of metal ion-binding site also helps in protein stabilization of thermophilic protease. Aromatic-aromatic interactions also play a crucial role in tPro stabilization. Formation of long network aromatic-aromatic interactions also first time reported here. Finally, the present study provides a major insight with a newly identified cyclic salt bridge in the stability of the enzyme, which may be helpful for protein engineering. It is also used in industrial applications for human welfare.

Computational Design of Nitrile Hydratase from Pseudonocardia thermophila JCM3095 for Improved Thermostability

High thermostability and catalytic activity are key properties for nitrile hydratase (NHase, EC 4.2.1.84) as a well-industrialized catalyst. In this study, rational design was applied to tailor the thermostability of NHase from Pseudonocardia thermophila JCM3095 (PtNHase) by combining FireProt server prediction and molecular dynamics (MD) simulation. Site-directed mutagenesis of non-catalytic residues provided by the rational design was subsequentially performed. The positive multiple-point mutant, namely, M10 (αI5P/αT18Y/αQ31L/αD92H/βA20P/βP38L/βF118W/βS130Y/βC189N/βC218V), was obtained and further analyzed. The Melting temperature (T_m) of the M10 mutant showed an increase by 3.2 °C and a substantial increase in residual activity of the enzyme at elevated temperatures was also observed. Moreover, the M10 mutant also showed a 2.1-fold increase in catalytic activity compared with the wild-type PtNHase. Molecular docking and MD simulations demonstrated better substrate affinity and improved thermostability for the mutant.

Enhanced Thermostability and Enzymatic Activity of Cel6A Variants from Thermobifida fusca by Empirical Domain Engineering (Short Title: Domain Engineering of Cel6A)

Cellulases are a set of lignocellulolytic enzymes, capable of producing eco-friendly low-cost renewable bioethanol. However, low stability and hydrolytic activity limit their wide-scale applicability at the industrial scale. In this work, we report the domain engineering of endoglucanase (Cel6A) of Thermobifida fusca to improve their catalytic activity and thermal stability. Later, enzymatic activity and thermostability of the most efficient variant named as Cel6A.CBC was analyzed by molecular dynamics simulations. This variant demonstrated profound activity against soluble and insoluble cellulosic substrates like filter paper, alkali-treated bagasse, regenerated amorphous cellulose (RAC), and bacterial microcrystalline cellulose. The variant Cel6A.CBC showed the highest catalysis of carboxymethyl cellulose (CMC) and other related insoluble substrates at a pH of 6.0 and a temperature of 60 °C. Furthermore, a sound rationale was observed between experimental findings and molecular modeling of Cel6A.CBC which revealed thermostability of Cel6A.CBC at 26.85, 60.85, and 74.85 °C as well as structural flexibility at 126.85 °C. Therefore, a thermostable derivative of Cel6A engineered in the present work has enhanced biological performance and can be a useful construct for the mass production of bioethanol from plant biomass.

Recent Advances and Promises in Nitrile Hydratase: From Mechanism to Industrial Applications

Nitrile hydratase (NHase, EC 4.2.1.84) is one type of metalloenzyme participating in the biotransformation of nitriles into amides. Given its catalytic specificity in amide production and eco-friendliness, NHase has overwhelmed its chemical counterpart during the past few decades. However, unclear catalytic mechanism, low thermostablity, and narrow substrate specificity limit the further application of NHase. During the past few years, numerous studies on the theoretical and industrial aspects of NHase have advanced the development of this green catalyst. This review critically focuses on NHase research from recent years, including the natural distribution, gene types, posttranslational modifications, expression, proposed catalytic mechanism, biochemical properties, and potential applications of NHase. The developments of NHase described here are not only useful for further application of NHase, but also beneficial for the development of the fields of biocatalysis and biotransformation.

Improving the Efficiency of Ligand-Binding Protein Design with Molecular Dynamics Simulations

Custom-designed ligand-binding proteins represent a promising class of macromolecules with exciting applications toward the design of new enzymes or the engineering of antibodies and small-molecule recruited proteins for therapeutic interventions. However, several challenges remain in designing a protein sequence such that the binding site organization results in high affinity interaction with a bound ligand. Here, we study the dynamics of explicitly solvated designed proteins through all-atom molecular dynamics (MD) simulations to gain insight into the causes that lead to the low affinity or instability of most of these designs, despite the prediction of their success by the computational design methodology. Simulations ranging from 500 to 1000 ns per replicate were conducted on 37 designed protein variants encompassing two distinct folds and a range of ligand affinities, resulting in more than 180 μs of combined sampling. The simulations provide retrospective insights into the properties affecting ligand affinity that can prove useful in guiding further steps of design optimization. Features indicate that entropic components are particularly important for affinity, which are not easily incorporated in the empirical models often used in design protocols. Additionally, we demonstrate that the application of machine learning approaches built upon the output from the simulations can help discriminate between successful and failed binders, such that MD could act as a screening step in protein design, resulting in a more efficient process.

Temperature-dependent conformational dynamics of cytochrome c: Implications in apoptosis

Heat, electric shock, and burn injuries induce apoptosis by releasing cytochrome c (cyt-c) from mitochondria and by subsequently activating the death protease, caspases-3. During apoptosis, cyt-c undergoes changes in the secondary structure that have been suggested to increase its peroxidase activity. Information about these structural changes will provide better understanding of the apoptotic mechanism. Hence, temperature-dependent conformational dynamics of cyt-c has been investigated through molecular dynamics (MD) simulations to explain the structural changes and to correlate them with its apoptotic behavior. We observe that, at lower temperatures (223, 248, and 300K), the secondary structure of cyt-c, remains stable, while at higher temperatures (323, 373, 423, and 473K), the secondary structural regions change significantly. Further, our MD results indicate that these structural changes are mainly localized on α-helices, turns, β-sheets, and important loops that were involved in the stabilization of the heme conformation. This conformational transition between specific regions of secondary structure of cyt-c directly affects the electron tunneling properties of the proteins as observed experimentally. We quantify and compare these changes and explain that the temperature plays a vital role in assuring the structural stability of cyt-c and thus its functions. Our findings from this MD study reproduce experimental results at high temperatures and provide evidence for the alteration of the heme through the disruption of the H-bonding interactions between specific regions of cyt-c, thereby enhancing its peroxidase activity which plays a crucial role in the apoptotic process.

A spectroscopic and molecular dynamics simulation approach towards the stabilizing effect of ammonium-based ionic liquids on bovine serum albumin

Experimental and theoretical evidence in support of the stabilizing effect of ammonium-based ionic liquids on thermal unfolding/refolding of bovine serum albumin is provided in this article.

Construction of a subunit-fusion nitrile hydratase and discovery of an innovative metal ion transfer pattern

Metallochaperones are metal-binding proteins designed to deliver the appropriate metal to a target protein. The metal is usually transferred between different proteins. In this study, we discovered that metal was transferred between the same subunit of a mutant nitrile hydratase (NHase). Various “activator proteins” mediate the trafficking of metal ions into NHases. We constructed fusion NHases by fusing the β- and α-subunits and/or the “activator proteins” of the NHase from Pseudomonas putida. The fusion NHases exhibited higher thermostability and tolerance to high concentrations of the product amide. The mechanism of the cobalt incorporation changed from a self-subunit swapping pattern to an apoprotein-specific molecular chaperone pattern in vivo and a metallochaperone pattern in vitro. Notably, the cobalt transfer occurred between the same α-subunit in the metallochaperone pattern. These results not only demonstrated the superiority of fusion-type NHases, but also revealed an innovative metal ion transfer pattern in metalloprotein biosynthesis.

The effect of ligands on the thermal stability of sulfotransferases: a molecular dynamics simulation study

Human cytosolic sulfotransferases (hSULTs) are important phase II metabolic enzymes. They catalyze transfer of the sulfuryl-group (-SO3) from 3'-phosphoadenosine 5'-phosphosulfate (PAPS) to the hydroxyl or primary amine moieties of a large number of endogenous and xenobiotic substrates. Broad selectivity and specificity of binding and activity within the sulfortransferases family could be detected by thermal denaturation assays, which have been made more and more suitable for high throughput screening based on recent technical advances. Here molecular dynamics simulations were used to explore the effect of the cofactor (PAPS) and substrate (LCA) on the thermal stability of the enzyme. It was found that the apo-enzyme unfolded fastest upon heating. The holo-enzyme with bound substrate LCA unfolded slowest. This thermo-denaturation order is consistent with that observed in experiments. Further it was found that the cofactor and substrate will pronouncedly increase the thermal stability of the active pocket regions that interact directly with the ligands. In addition, cofactor and substrate show noticeable synergy effect on the thermal stability of the enzyme.

Understanding the thermostability and activity of Bacillus subtilis lipase mutants: insights from molecular dynamics simulations

Improving the thermostability of industrial enzymes is an important protein engineering challenge. Point mutations, induced to increase thermostability, affect the structure and dynamics of the target protein in several ways and thus can also affect its activity. There appears to be no general rules for improving the thermostabilty of enzymes without adversely affecting their enzymatic activity. We report MD simulations, of wild type Bacillus subtilis lipase (WT) and its six progressively thermostable mutants (2M, 3M, 4M, 6M, 9M, and 12M), performed at different temperatures, to address this issue. Less thermostable mutants (LTMs), 2M to 6M, show WT-like dynamics at all simulation temperatures. However, the two more thermostable mutants (MTMs) show the required flexibility at appropriate temperature ranges and maintain conformational stability at high temperature. They show a deep and rugged free-energy landscape, confining them within a near-native conformational space by conserving noncovalent interactions, and thus protecting them from possible aggregation. In contrast, the LTMs having marginally higher thermostabilities than WT show greater probabilities of accessing non-native conformations, which, due to aggregation, have reduced possibilities of reverting to their respective native states under refolding conditions. Our analysis indicates the possibility of nonadditive effects of point mutations on the conformational stability of LTMs.

The role of proline substitutions within flexible regions on thermostability of luciferase

Improving the stability of firefly luciferase has been a critical issue for its wider industrial applications. Studies about hyperthermophile proteins show that flexibility could be an effective indicator to find out weak spots to engineering thermostability of proteins. However, the relationship among flexibility, activity and stability in most of proteins is unclear. Proline is the most rigid residue and can be introduced to rigidify flexible regions to enhance thermostability of proteins. We firstly apply three different methods, molecular dynamics (MD) simulation, B-FITTER and framework rigidity optimized dynamics algorithm (FRODA) to determine the flexible regions of Photinus pyralis luciferase: Fragment 197-207; Fragment 471-481 and Fragment 487-495. Then, introduction of proline is used to rigidify these flexible regions. Two mutants D476P and H489P within most flexible regions are finally designed. In the results, H489P mutant shows improved thermostability while maintaining its catalytic efficiency compared to that of wild type luciferase. Flexibility analysis confirms that the overall rigidity and local rigidity of H489P mutant are greatly strengthened. D476P mutant shows decreased thermosatbility and the reason for this is elucidated at the molecular level. S307P mutation is randomly chosen outside the flexible regions as a control. Thermostability analysis shows that S307P mutation has decreased kinetic stability and enhanced thermodynamic stability.

Biochemical characterization of a carboxylesterase from the archaeon Pyrobaculum sp. 1860 and a rational explanation of its substrate specificity and thermostability

In this work, genome mining was used to identify esterase/lipase genes in the archaeon Pyrobaculum sp. 1860. A gene was cloned and functionally expressed in Escherichia coli as His-tagged protein. The recombinant enzyme (rP186_1588) was verified by western blotting and peptide mass fingerprinting. Biochemical characterization revealed that rP186_1588 exhibited optimum activity at pH 9.0 and 80 °C towards p-nitrophenyl acetate (K(m): 0.35 mM, k(cat): 11.65 s⁻¹). Interestingly, the purified rP186_1588 exhibited high thermostability retaining 70% relative activity after incubation at 90 °C for 6 h. Circular dichroism results indicated that rP186_1588 showed slight structure alteration from 60 to 90 °C. Structural modeling showed P186_1588 possessed a typical α/β hydrolase's fold with the catalytic triad consisting of Ser97, Asp147 and His172, and was further confirmed by site-directed mutagenesis. Comparative molecular simulations at different temperatures (300, 353, 373 and 473 K) revealed that its thermostability was associated with its conformational rigidity. The binding free energy analysis by MM-PBSA method revealed that the van der Waals interaction played a major role in p-NP ester binding for P186_1588. Our data provide insights into the molecular structures of this archaeal esterase, and may help to its further protein engineering for industrial applications.

Surface modification of polyacrylonitrile fibre by nitrile hydratase from Corynebacterium nitrilophilus

Previously, nitrile hydratase (NHase) from Corynebacterium nitrilophilus was obtained and showed potential in polyacrylonitrile (PAN) fibre modification. In the present study, the modification conditions of C. nitrilophilus NHase on PAN were investigated. In the optimal conditions, the wettability and dyeability (anionic and reactive dyes) of PAN treated by C. nitrilophilus NHase reached a similar level of those treated by alkali. In addition, the chemical composition and microscopically observable were changed in the PAN surface after NHase treatment. Meanwhile, it revealed that cutinase combined with NHase facilitates the PAN hydrolysis slightly because of the ester existed in PAN as co-monomer was hydrolyzed. All these results demonstrated that C. nitrilophilus NHase can modify PAN efficiently without textile structure damage, and this study provides a foundation for the further application of C. nitrilophilus NHase in PAN modification industry.

Enhanced thermal stability and hydrolytic ability of Bacillus subtilis aminopeptidase by removing the thermal sensitive domain in the non-catalytic region

Besides the catalytic ability, many enzymes contain conserved domains to perform some other physiological functions. However, sometimes these conserved domains were unnecessary or even detrimental to the catalytic process for industrial application of the enzymes. In this study, based on homology modeling and molecular dynamics simulations, we found that Bacillus subtilis aminopeptidase contained a thermal sensitive domain (protease-associated domain) in the non-catalytic region, and predicted that deletion of this flexible domain can enhance the structure stability. This prediction was then verified by the deletion of protease-associated domain from the wild-type enzyme. The thermal stability analysis showed that deletion of this domain improved the T50 (the temperature required to reduce initial activity by 50% in 30 min) of the enzyme from 71 °C to 77 °C. The melting temperature (Tm) of the enzyme also increased, which was measured by thermal denaturation experiments using circular dichroism spectroscopy. Further studies indicated that this deletion did not affect the activity and specificity of the enzyme toward aminoacyl-p-nitroanilines, but improved its hydrolytic ability toward a 12-aa-long peptide (LKRLKRFLKRLK) and soybean protein. These findings suggested the possibility of a simple technique for enzyme modification and the artificial enzyme obtained here was more suitable for the protein hydrolysis in food industry than the wild-type enzyme.

Engineering proteins for thermostability through rigidifying flexible sites

Engineering proteins for thermostability is an exciting and challenging field since it is critical for broadening the industrial use of recombinant proteins. Thermostability of proteins arises from the simultaneous effect of several forces such as hydrophobic interactions, disulfide bonds, salt bridges and hydrogen bonds. All of these interactions lead to decreased flexibility of polypeptide chain. Structural studies of mesophilic and thermophilic proteins showed that the latter need more rigid structures to compensate for increased thermal fluctuations. Hence flexibility can be an indicator to pinpoint weak spots for enhancing thermostability of enzymes. A strategy has been proven effective in enhancing proteins' thermostability with two steps: predict flexible sites of proteins firstly and then rigidify these sites. We refer to this approach as rigidify flexible sites (RFS) and give an overview of such a method through summarizing the methods to predict flexibility of a protein, the methods to rigidify residues with high flexibility and successful cases regarding enhancing thermostability of proteins using RFS.

Stability mechanisms of a thermophilic laccase probed by molecular dynamics

Laccases are highly stable, industrially important enzymes capable of oxidizing a large range of substrates. Causes for their stability are, as for other proteins, poorly understood. In this work, multiple-seed molecular dynamics (MD) was applied to a Trametes versicolor laccase in response to variable ionic strengths, temperatures, and glycosylation status. Near-physiological conditions provided excellent agreement with the crystal structure (average RMSD ∼0.92 Å) and residual agreement with experimental B-factors. The persistence of backbone hydrogen bonds was identified as a key descriptor of structural response to environment, whereas solvent-accessibility, radius of gyration, and fluctuations were only locally relevant. Backbone hydrogen bonds decreased systematically with temperature in all simulations (∼9 per 50 K), probing structural changes associated with enthalpy-entropy compensation. Approaching T_opt (∼350 K) from 300 K, this change correlated with a beginning “unzipping” of critical β-sheets. 0 M ionic strength triggered partial denucleation of the C-terminal (known experimentally to be sensitive) at 400 K, suggesting a general salt stabilization effect. In contrast, F⁻ (but not Cl⁻) specifically impaired secondary structure by formation of strong hydrogen bonds with backbone NH, providing a mechanism for experimentally observed small anion destabilization, potentially remedied by site-directed mutagenesis at critical intrusion sites. N-glycosylation was found to support structural integrity by increasing persistent backbone hydrogen bonds by ∼4 across simulations, mainly via prevention of F⁻ intrusion. Hydrogen-bond loss in distinct loop regions and ends of critical β-sheets suggest potential strategies for laboratory optimization of these industrially important enzymes.

Engineering of Rhodococcus cell catalysts for tolerance improvement by sigma factor mutation and active plasmid partition

Tolerance to various stresses is a key phenotype for cell catalysts, which are used widely in bioproduction of diverse valuable chemicals. Using the Rhodococcus ruber TH strain, which exhibits high nitrile hydratase activity, as the target cell catalyst for acrylamide production, we established a method to improve cell tolerance by stably introducing global transcription perturbation. The σ70 gene (sigA) of R. ruber was cloned and randomly mutated. An R. ruber TH3/pNV-sigAM library containing additional sigA mutants was constructed and used for survival selection. The TH3/M4N1-59 mutant was selected by acrylonitrile/acrylamide double stress and exhibited a 160 % extension of the half-life of nitrile hydratase upon exposure to 40 % acrylamide. A redesigned parDEM gene was introduced to Rhodococcus to accomplish stable inheritance of plasmids. A two-batch acrylonitrile hydration reaction was performed using the engineered cells as a catalyst. Compared to TH3, the acrylamide productivity of TH3/M4N1-59DEM catalysis increased by 27.8 and 37.5 % in the first and second bioreaction batches, respectively. These data suggest a novel method for increasing the bioconversion productivity of target chemicals through sigA mutation of the cell catalyst.

Comparison of the structural basis for thermal stability between archaeal and bacterial proteins

In this study, the structural basis for thermal stability in archaeal and bacterial proteins was investigated. There were many common factors that confer resistance to high temperature in both archaeal and bacterial proteins. These factors include increases in the Lys content, the bends and blanks of secondary structure, the Glu content of salt bridge; decreases in the number of main-side chain hydrogen bond and exposed surface area, and changes in the bends and blanks of amino acids. Certainly, the utilization of charged amino acids to form salt bridges is a primary factor. In both heat-resistant archaeal and bacterial proteins, most Glu and Asp participate in the formation of salt bridges. Other factors may influence either archaeal or bacterial protein thermostability, which includes the more frequent occurrence of shorter 3(10)-helices and increased hydrophobicity in heat-resistant archaeal proteins. However, there were increases in average helix length, the Glu content in salt bridges, temperature factors and decreases in the number of main-side chain hydrogen bonds, uncharged-uncharged hydrogen bonds, hydrophobicity, and buried and exposed polar surface area in heat-resistant bacterial proteins. Evidently, there are few similarities and many disparities between the heat-resistant mechanisms of archaeal and bacterial proteins.

Increased number of Arginine-based salt bridges contributes to the thermotolerance of thermotolerant acetic acid bacteria, Acetobacter tropicalis SKU1100

Thermotolerant acetic acid bacteria (AAB), Acetobacter tropicalis SKU1100, can grow above 40°C. To investigate the basis of its thermotolerance, we compared the genome of A. tropicalis SKU1100 with that of mesophilic AAB strain Acetobacter pasteurianus IFO3283-01. The comparative genomic study showed that amino acid substitutions from large to small residue and Lys to Arg occur in many orthologous genes. Furthermore, comparative modeling study was carried out with the orthologous proteins between SKU1100 and IFO3283-01 strains, indicating that the number of Arg-based salt bridges increased in protein models. Since it has been reported that Arg-based salt bridges are important factor for thermo-stability of protein structure, our results strongly suggest that the increased number of Arg-based salt bridges may contributes to the thermotolerance of A. tropicalis SKU1100 (the thermo-stability of proteins in A. tropicalis SKU1100).

Enhanced thermostability of methyl parathion hydrolase from Ochrobactrum sp. M231 by rational engineering of a glycine to proline mutation

Protein thermostability can be increased by some glycine to proline mutations in a target protein. However, not all glycine to proline mutations can improve protein thermostability, and this method is suitable only at carefully selected mutation sites that can accommodate structural stabilization. In this study, homology modeling and molecular dynamics simulations were used to select appropriate glycine to proline mutations to improve protein thermostability, and the effect of the selected mutations was proved by the experiments. The structure of methyl parathion hydrolase (MPH) from Ochrobactrum sp. M231 (Ochr-MPH) was constructed by homology modeling, and molecular dynamics simulations were performed on the modeled structure. A profile of the root mean square fluctuations of Ochr-MPH was calculated at the nanosecond timescale, and an eight-amino acid loop region (residues 186-193) was identified as having high conformational fluctuation. The two glycines nearest to this region were selected as mutation targets that might affect protein flexibility in the vicinity. The structures and conformational fluctuations of two single mutants (G194P and G198P) and one double mutant (G194P/G198P) were modeled and analyzed using molecular dynamics simulations. The results predicted that the mutant G194P had the decreased conformational fluctuation in the loop region and might increase the thermostability of Ochr-MPH. The thermostability and kinetic behavior of the wild-type and three mutant enzymes were measured. The results were consistent with the computational predictions, and the mutant G194P was found to have higher thermostability than the wild-type enzyme.

Using molecular dynamics to probe the structural basis for enhanced stability in thermal stable cytochromes P450

High-temperature molecular dynamics (MD) has been used to assess if MD can be employed as a useful tool for probing the structural basis for enhanced stability in thermal stable cytochromes P450. CYP119, the most thermal stable P450 known, unfolds more slowly during 500 K MD simulations than P450s that melt at lower temperatures, P450cam and P450cin. A comparison of the 500 K MD trajectories shows that the Cys ligand loop, a critically important structural feature just under the heme, in both P450cin and P450cam completely unfolds while this region is quite stable in CYP119. In CYP119, this region is stabilized by tight nonpolar interactions involving Tyr26 and Leu308. The corresponding residues in P450cam are Gly and Thr, respectively. The in silico generated Y26A/L308A CYP119 double mutant is substantially less stable than wild-type CYP119, and the Cys ligand loop unfolds in a manner similar to that of P450cam. The MD thus has identified a potential "hot spot" important for stability. As an experimental test of the MD results, the Y26A/L308A double mutant was prepared, and thermal melting curves show that the double mutant exhibits a melting temperature (T(m)) 16 degrees C lower than that of wild-type CYP119. Control mutations that were predicted by MD not to destabilize the protein were also generated, and the experimental melting temperature was not significantly different from that of the wild-type enzyme. Therefore, high-temperature MD is a useful tool in predicting the structural underpinnings of thermal stability in P450s.

Temperature dependence of the flexibility of thermophilic and mesophilic flavoenzymes of the nitroreductase fold

A widely held hypothesis regarding the thermostability of thermophilic proteins states asserts that, at any given temperature, thermophilic proteins are more rigid than their mesophilic counterparts. Many experimental and computational studies have addressed this question with conflicting results. Here, we compare two homologous enzymes, one mesophilic (Escherichia coli FMN-dependent nitroreductase; NTR) and one thermophilic (Thermus thermophilus NADH oxidase; NOX), by multiple molecular dynamics simulations at temperatures from 5 to 100 degrees C. We find that the global rigidity/flexibility of the two proteins, assessed by a variety of metrics, is similar on the time scale of our simulations. However, the thermophilic enzyme retains its native conformation to a much greater degree at high temperature than does the mesophilic enzyme, both globally and within the active site. The simulations identify the helix F-helix G 'arm' as the region with the greatest difference in loss of native contacts between the two proteins with increasing temperature. In particular, a network of electrostatic interactions holds helix F to the body of the protein in the thermophilic protein, and this network is absent in the mesophilic counterpart.

Structural study of carboxylesterase from hyperthermophilic bacteria Geobacillus stearothermophilus by molecular dynamics simulation

Carboxylesterases are ubiquitous enzymes with important physiological, industrial and medical applications such as synthesis and hydrolysis of stereo specific compounds, including the metabolic processing of drugs, and antimicrobial agents. Here, we have performed molecular dynamics simulations of carboxylesterase from hyperthermophilic bacterium Geobacillus stearothermophilus (GsEst) for 10ns each at five different temperatures namely at 300K, 343K, 373K, 473K and 500K. Profiles of root mean square fluctuation (RMSF) identify thermostable and thermosensitive regions of GsEst. Unfolding of GsEst initiates at the thermosensitive alpha-helices and proceeds to the thermostable beta-sheets. Five ion-pairs have been identified as critical ion-pairs for thermostability and are maintained stably throughout the higher temperature simulations. A detailed investigation of the active site residues of this enzyme suggests that the geometry of this site is well preserved up to 373K. Furthermore, the hydrogen bonds between Asp188 and His218 of the active site are stably maintained at higher temperatures imparting stability of this site. Radial distribution functions (RDFs) show similar pattern of solvent ordering and water penetration around active site residues up to 373K. Principal component analysis suggests that the motion of the entire protein as well as the active site is similar at 300K, 343K and 373K. Our study may help to identify the factors responsible for thermostability of GsEst that may endeavor to design enzymes with enhanced thermostability.