Structural and dynamical thermostability of psychrophilic enzyme at various temperatures: Molecular dynamics simulations of tryptophan synthase

A Molecular Dynamics Simulation Study of the Effects of βGln114 Mutation on the Dynamic Behavior of the Catalytic Site of the Tryptophan Synthase

L-tryptophan (l-Trp), a vital amino acid for the survival of various organisms, is synthesized by the enzyme tryptophan synthase (TS) in organisms such as eubacteria, archaebacteria, protista, fungi, and plantae. TS, a pyridoxal 5'-phosphate (PLP)-dependent enzyme, comprises α and β subunits that typically form an α2β2 tetramer. The enzyme's activity is regulated by the conformational switching of its α and β subunits between the open (T state) and closed (R state) conformations. Many microorganisms rely on TS for growth and replication, making the enzyme and the l-Trp biosynthetic pathway potential drug targets. For instance, Mycobacterium tuberculosis, Chlamydiae bacteria, Streptococcus pneumoniae, Francisella tularensis, Salmonella bacteria, and Cryptosporidium parasitic protozoa depend on l-Trp synthesis. Antibiotic-resistant salmonella strains have emerged, underscoring the need for novel drugs targeting the l-Trp biosynthetic pathway, especially for salmonella-related infections. A single amino acid mutation can significantly impact enzyme function, affecting stability, conformational dynamics, and active or allosteric sites. These changes influence interactions, catalytic activity, and protein-ligand/protein-protein interactions. This study focuses on the impact of mutating the βGln114 residue on the catalytic and allosteric sites of TS. Extensive molecular dynamics simulations were conducted on E(PLP), E(AEX1), E(A-A), and E(C3) forms of TS using the WT, βQ114A, and βQ114N versions. The results show that both the βQ114A and βQ114N mutations increase protein backbone root mean square deviation fluctuations, destabilizing all TS forms. Conformational and hydrogen bond analyses suggest the significance of βGln114 drifting away from cofactor/intermediates and forming hydrogen bonds with water molecules necessary for l-Trp biosynthesis. The βQ114A mutation creates a gap between βAla114 and cofactor/intermediates, hindering hydrogen bond formation due to short side chains and disrupting β-sites. Conversely, the βQ114N mutation positions βAsn114 closer to cofactor/intermediates, forming hydrogen bonds with O3 of cofactors/intermediates and nearby water molecules, potentially disrupting the l-Trp biosynthetic mechanism.

Helping proteins come in from the cold: 5 burning questions about cold-active enzymes

Enzymes from psychrophilic (cold-loving) organisms have attracted considerable interest over the past decades for their potential in various low-temperature industrial processes. However, we still lack large-scale commercialization of their activities. Here, we review their properties, limitations and potential. Our review is structured around answers to 5 central questions: 1. How do cold-active enzymes achieve high catalytic rates at low temperatures? 2. How is protein flexibility connected to cold-activity? 3. What are the sequence-based and structural determinants for cold-activity? 4. How does the thermodynamic stability of psychrophilic enzymes reflect their cold-active capabilities? 5. How do we effectively identify novel cold-active enzymes, and can we apply them in an industrial context? We conclude that emerging screening technologies combined with big-data handling and analysis make it reasonable to expect a bright future for our understanding and exploitation of cold-active enzymes.

Investigation of thermal stability characteristic in family A DNA polymerase - A theoretical study

DNA polymerases create complementary DNA strands in living cells and are crucial to genome transmission and maintenance. These enzymes possess similar human right-handed folds which contain thumb, fingers, and palm subdomains and contribute to polymerization activities. These enzymes are classified into seven evolutionary families, A, B, C, D, X, Y, and RT, based on amino acid sequence analysis and biochemical characteristics. Family A DNA polymerases exist in an extended range of organisms including mesophilic, thermophilic, and hyper-thermophilic bacteria, participate in DNA replication and repair, and have a broad application in molecular biology and biotechnology. In this study, we attempted to detect factors that play a role in the thermostability properties of this family member despite their remarkable similarities in structure and function. For this purpose, similarities and differences in amino acid sequences, structure, and dynamics of these enzymes have been inspected. Our results demonstrated that thermophilic and hyper-thermophilic enzymes have more charged, aromatic, and polar residues than mesophilic ones and consequently show further electrostatic and cation-pi interactions. In addition, in thermophilic enzymes, aliphatic residues tend to position in buried states more than mesophilic enzymes. These residues within their aliphatic parts increase hydrophobic core packing and therefore enhance the thermostability of these enzymes. Furthermore, a decrease in thermophilic cavities volumes assists in the protein compactness enhancement. Moreover, molecular dynamic simulation results revealed that increasing temperature impacts mesophilic enzymes further than thermophilic ones that reflect on polar and aliphatic residues surface area and hydrogen bonds changes.

Mining of key genes for cold adaptation from Pseudomonas fragi D12 and analysis of its cold-adaptation mechanism

The psychrotroph Pseudomonas fragi D12, which grew strongly under low temperatures, was screened from tundra soil collected from the permanent alpine zone on Changbai Mountain. To mine the genes critical for cold tolerance and to investigate the cold-adaptation mechanism, whole-genome sequencing, comparative genomic analysis, and transcriptome analysis were performed with P. fragi. A total of 124 potential cold adaptation genes were identified, including nineteen unique cold-adaptive genes were detected in the genome of P. fragi D12. Three unique genes associated with pili protein were significantly upregulated at different degrees of low temperature, which may be the key to the strong low-temperature adaptability of P. fragi D12. Meanwhile, we were pleasantly surprised to find that Pseudomonas fragi D12 exhibited different cold-adaptation mechanisms under different temperature changes. When the temperature declined from 30°C to 15°C, the response included maintenance of the fluidity of cell membranes, increased production of extracellular polymers, elevation in the content of compatibility solutes, and reduction in the content of reactive oxygen species, thereby providing a stable metabolic environment. When the temperature decreased from 15°C to 4°C, the response mainly included increases in the expression of molecular chaperones and transcription factors, enabling the bacteria to restore normal transcription and translation. The response mechanism of P. fragi D12 to low-temperature exposure is discussed. The results provide new ideas for the cold-adaptation mechanism of cold-tolerant microorganisms.

Allosteric regulation of β-reaction stage I in tryptophan synthase upon the α-ligand binding

Tryptophan synthase (TRPS) is a bifunctional enzyme consisting of α- and β-subunits that catalyzes the last two steps of L-tryptophan (L-Trp) biosynthesis. The first stage of the reaction at the β-subunit is called β-reaction stage I, which converts the β-ligand from an internal aldimine [E(Ain)] to an α-aminoacrylate [E(A-A)] intermediate. The activity is known to increase 3-10-fold upon the binding of 3-indole-D-glycerol-3'-phosphate (IGP) at the α-subunit. The effect of α-ligand binding on β-reaction stage I at the distal β-active site is not well understood despite the abundant structural information available for TRPS. Here, we investigate the β-reaction stage I by carrying out minimum-energy pathway searches based on a hybrid quantum mechanics/molecular mechanics (QM/MM) model. The free-energy differences along the pathway are also examined using QM/MM umbrella sampling simulations with QM calculations at the B3LYP-D3/aug-cc-pVDZ level of theory. Our simulations suggest that the sidechain orientation of βD305 near the β-ligand likely plays an essential role in the allosteric regulation: a hydrogen bond is formed between βD305 and the β-ligand in the absence of the α-ligand, prohibiting a smooth rotation of the hydroxyl group in the quinonoid intermediate, whereas the dihedral angle rotates smoothly after the hydrogen bond is switched from βD305-β-ligand to βD305-βR141. This switch could occur upon the IGP-binding at the α-subunit, as evidenced by the existing TRPS crystal structures.

Computational Analysis on the Allostery of Tryptophan Synthase: Relationship between α/β-Ligand Binding and Distal Domain Closure

Tryptophan synthase (TRPS) is a bifunctional enzyme consisting of α and β-subunits and catalyzes the last two steps of l-tryptophan (L-Trp) biosynthesis, namely, cleavage of 3-indole-d-glycerol-3'-phosphate (IGP) into indole and glyceraldehyde-3-phosphate (G3P) in the α-subunit, and a pyridoxal phosphate (PLP)-dependent reaction of indole and l-serine (L-Ser) to produce L-Trp in the β-subunit. Importantly, the IGP binding at the α-subunit affects the β-subunit conformation and its ligand-binding affinity, which, in turn, enhances the enzymatic reaction at the α-subunit. The intersubunit communications in TRPS have been investigated extensively for decades because of the fundamental and pharmaceutical importance, while it is still difficult to answer how TRPS allostery is regulated at the atomic detail. Here, we investigate the allosteric regulation of TRPS by all-atom classical molecular dynamics (MD) simulations and analyze the potential of mean-force (PMF) along conformational changes of the α- and β-subunits. The present simulation has revealed a widely opened conformation of the β-subunit, which provides a pathway for L-Ser to enter into the β-active site. The IGP binding closes the α-subunit and induces a wide opening of the β-subunit, thereby enhancing the binding affinity of L-Ser to the β-subunit. Structural analyses have identified critical hydrogen bonds (HBs) at the interface of the two subunits (αG181-βS178, αP57-βR175, etc.) and HBs between the β-subunit (βT110 - βH115) and a complex of PLP and L-Ser (an α-aminoacrylate intermediate). The former HBs regulate the allosteric, β-subunit opening, whereas the latter HBs are essential for closing the β-subunit in a later step. The proposed mechanism for how the interdomain communication in TRPS is realized with ligand bindings is consistent with the previous experimental data, giving a general idea to interpret the allosteric regulations in multidomain proteins.

Insights of conformational dynamics on catalytic activity in the computational stability design of Bacillus subtilis LipA

In protein engineering, the contributions of individual mutations to designed combinatorial mutants are unpredictable. Screening designed mutations that affect enzyme catalytic activity enables evolutions towards efficient activities. Here, Bacillus subtilis LipA (BSLA) was selected as a model protein for thermostabilization designs, and the circular dichroism measurements showed six combinatorial designs with improved stability (from 5.81 °C to 13.61 °C). Based on molecular dynamic simulations, the conformational dynamics of the mutants revealed that mutations alter the populations of conformational states and the increased ensembles of inactive conformations might lead to a reduction in activity. We further demonstrated that the mutations responsible for the reduced enzyme catalytic activity involved a short dynamic correlation path to disturbing the equilibrium conformation of active sites. By removing N82V, which had a close dynamic correlation to the active sites in mutant D3, the redesigned mutant RD3 had an increased activity of 57.6%. By combining computational simulation with experimental verification, this work established that essential sites to counteract the activity-stability trade-off in multipoint combinatorial mutants could be computationally predicted and thus provide a possible strategy by which to indirectly or directly guide protein design.

Development of Fluorescence Polarization Immunoassay With scFv to Detect Fumonisin Bs in Maize and Simultaneous Study of Their Molecular Recognition Mechanism

In this study, a fluorescence polarization immunoassay (FPIA) was developed based on the single-chain variable fragments (scFvs) for fumonisin B_s (FB_s). The scFvs were prepared from FB_s-specific monoclonal antibody secreting hybridomas (4F5 and 4B9). The established FPIA could determine the sum of fumonisin B₁ (FB₁) and fumonisin B₂ (FB₂) within a short time. The IC₅₀ of FPIA for the detection of FB₁ and FB₂ were 29.36 ng/ml and 1,477.82 ng/ml with 4F5 scFv, and 125.16 ng/ml and 30.44 ng/ml with 4B9 scFv, so the 4B9 scFv was selected for detection of FB₁ and FB₂ in maize samples with a limit of detection of 441.54 μg/kg and 344.933 μg/kg. The recoveries ranged from 84.7 to 104.1% with a coefficient of variation less than 14.1% in spiked samples, and the result of the FPIA method was in good consistency with that of HPLC-MS/MS. To supply a better understanding of the immunoassay results, the interactions mechanism of scFvs-FB_s was further revealed by the homology modelling, molecular docking, and molecular dynamic simulation. It was indicated that six complementarity-determining regions (CDRs) were involved in 4B9 scFv recognition, forming a narrow binding cavity, and FB₁/FB₂ could be inserted into this binding cavity stably through strong hydrogen bonds and other interactions. While in 4F5 scFv, only the FB₁ stably inserted in the binding pocket formed by four CDRs through strong hydrogen bonds, and FB₂ did not fit the binding cavity due to the lack of hydroxyl at C10, which is the key recognition site of 4F5 scFv. Also, the binding energy of FB₂-4B9 scFv complex is higher than the FB₂-4F5 scFv complex. This study established a FPIA method with scFv for the detection of FB₁ and FB₁ in maize, and systematically predicted recognition mechanism of FB_s and scFvs, which provided a reference for the better understanding of the immunoassay mechanism.

Computational design of a cutinase for plastic biodegradation by mining molecular dynamics simulations trajectories

Polyethylene terephthalate (PET) has caused serious environmental concerns but could be degraded at high temperature. Previous studies show that cutinase from Thermobifida fusca KW3 (TfCut2) is capable of degrading and upcycling PET but is limited by its thermal stability. Nowadays, Popular protein stability modification methods rely mostly on the crystal structures, but ignore the fact that the actual conformation of protein is complex and constantly changing. To solve these problems, we developed a computational approach to design variants with enhanced protein thermal stability by mining Molecular Dynamics simulation trajectories using Machine Learning methods (MDL). The optimal classification accuracy and the optimal Pearson correlation coefficient of MDL model were 0.780 and 0.716, respectively. And we successfully designed variants with high ΔT_m values using MDL method. The optimal variant S121P/D174S/D204P had the highest ΔT_m value of 9.3 °C, and the PET degradation ratio increased by 46.42-fold at 70℃, compared with that of wild type TfCut2. These results deepen our understanding on the complex conformations of proteins and may enhance the plastic recycling and sustainability at glass transition temperature.

Thermostable Cellulases / Xylanases From Thermophilic and Hyperthermophilic Microorganisms: Current Perspective

The bioconversion of lignocellulose into monosaccharides is critical for ensuring the continual manufacturing of biofuels and value-added bioproducts. Enzymatic degradation, which has a high yield, low energy consumption, and enhanced selectivity, could be the most efficient and environmentally friendly technique for converting complex lignocellulose polymers to fermentable monosaccharides, and it is expected to make cellulases and xylanases the most demanded industrial enzymes. The widespread nature of thermophilic microorganisms allows them to proliferate on a variety of substrates and release substantial quantities of cellulases and xylanases, which makes them a great source of thermostable enzymes. The most significant breakthrough of lignocellulolytic enzymes lies in lignocellulose-deconstruction by enzymatic depolymerization of holocellulose into simple monosaccharides. However, commercially valuable thermostable cellulases and xylanases are challenging to produce in high enough quantities. Thus, the present review aims at giving an overview of the most recent thermostable cellulases and xylanases isolated from thermophilic and hyperthermophilic microbes. The emphasis is on recent advancements in manufacturing these enzymes in other mesophilic host and enhancement of catalytic activity as well as thermostability of thermophilic cellulases and xylanases, using genetic engineering as a promising and efficient technology for its economic production. Additionally, the biotechnological applications of thermostable cellulases and xylanases of thermophiles were also discussed.

Nobiletin, a hexamethoxyflavonoid from citrus pomace, attenuates G1 cell cycle arrest and apoptosis in hypoxia-induced human trophoblast cells of JEG-3 and BeWo via regulating the p53 signaling pathway

Background

Hypoxia is associated with abnormal cell apoptosis in trophoblast cells, which causes fetal growth restriction and related placental pathologies. Few effective methods for the prevention and treatment of placenta-related diseases exist. Natural products and functional foods have always been a rich source of potential anti-apoptotic drugs. Nobiletin (NOB), a hexamethoxyflavonoid derived from the citrus pomace, shows an anti-apoptotic activity, which is a non-toxic constituent of dietary phytochemicals approved by the Food and Drug Administration. However, their effects on hypoxia-induced human trophoblast cells have not been fully studied.

Objective

The aim of this study was to investigate the protective effects of NOB on hypoxia-induced apoptosis of human trophoblast JEG-3 and BeWo cells, and their underlying mechanisms.

Design

First, the protective effect of NOB on hypoxia-induced apoptosis of JEG-3 and BeWo cells was studied. Cell viability and membrane integrity were determined by CCK-8 assay and lactate dehydrogenase activity, respectively. Real Time Quantitative PCR (RT-qPCR) and Western blot analysis were used to detect the mRNA and protein levels of HIF1α. Propidium iodide (PI)-labeled flow cytometry was used to detect cell cycle distribution. Cell apoptosis was detected by flow cytometry with Annexin V-FITC and PI double staining, and the expression of apoptosis marker protein cl-PARP was detected by Western blot analysis. Then, the molecular mechanism of NOB against apoptosis was investigated. Computer molecular docking and dynamics were used to simulate the interaction between NOB and p53 protein, and this interaction was verified in vitro by Ultraviolet and visible spectrum (UV-visible spectroscopy), fluorescence spectroscopy and circular dichroism. Furthermore, the changes in the expression of p53 signaling pathway genes and proteins were detected by RT-qPCR and Western blot analysis, respectively.

Results

Hypoxia treatment resulted in a decreased cell viability and cell membrane integrity in JEG-3 and BeWo cell lines, and an increased expression of HIF1α, cell cycle arrest in the G1 phase, and massive cell apoptosis, which were alleviated after NOB treatment. Molecular docking and dynamics simulations found that NOB spontaneously bonded to human p53 protein, leading to the change of protein conformation. The intermolecular interaction between NOB and human p53 protein was further confirmed by UV-visible spectroscopy, fluorescence spectroscopy and circular dichroism. After the treatment of 100 μM NOB, a down-regulation of mRNA and protein levels of p53 and p21 and an up-regulation of BCL2/BAX mRNA and protein ratio were observed in JEG-3 cells; however, there was also a down-regulation of mRNA and protein levels observed for p53 and p21 in BeWo cells after the treatment of NOB. The BCL2/BAX ratio of BeWo cells did not change after the treatment of 100 μM NOB.

Conclusion

NOB attenuated hypoxia-induced apoptosis in JEG-3 and BeWo cell lines and might be a potential functional ingredient to prevent pregnancy-related diseases caused by hypoxia-induced apoptosis. These findings would also suggest the exploration and utilization of citrus resources, and the development of citrus industry.

Application of bioinformatics and molecular dynamics simulation approaches for identification of fibroblast growth factor 10 analogues with potentially improved thermostability

Fibroblast growth factor 10 functions as a paracrine mesenchymal molecule to initiate signalling pathways regarding to cellular development and health. However, the low thermal stability restricts it's functionality in the human body and the shelf-life of FGF10-based formulations. The current study aimed to employ rational design and bioinformatics approaches to identify some point mutations which may improve the thermal stability of FGF10. Bioinformatics analyses resulted in N105D, C106F, K144R, K153M and I156R as the potential stability conferring mutations. The identified mutants were subjected to MD simulation indicating that all mutations are both structurally and energetically favoured. Finally, the effects of the identified mutations on receptor binding of FGF10 were predicted and the results showed that K144R and K153M mutations may increase the binding affinity relative to the wild type. The findings of the current study propose potentially improved FGF10 analogues for further experimental investigations.

Improving the thermostability and acid resistance of Rhizopus oryzae α-amylase by using multiple sequence alignment based site-directed mutagenesis

Higher thermostability or acid resistance for fungal α-amylase will help to improve the sugar-making process and cut down the production costs. Here, the thermostability or acid resistance of Rhizopus oryzae α-amylase (ROAmy) was significantly enhanced by site-directed evolution based on multiple sequence alignment (MSA) method. For instance, compared with the wild-type ROAmy, the optimum temperature of mutants G136D and A144Y was increased from 50 to 55 °C, whereas for mutants V174R and I276P, the optimum temperature was increased from 50 to 60 °C. The optimum pH of mutants G136D and A144Y shifted from 5.5 to 5.0, whereas for mutants V174R and T253E, the optimum pH changed from 5.5 to 4.5. The results showed that mutant V174R had a 2.52-fold increase in half-life at 55 °C, a 2.55-fold increase in half-life at pH 4.5, and a 1.61-fold increase in catalytic efficiency (k_cat /K_m ) on soluble starch. The three-dimensional model simulation revealed that changes of hydrophilicity, hydrogen bond, salt bridge, or rigidity observed in mutants might mainly account for the improvement of thermostability and acid resistance. The mutants with improved catalytic properties attained in this work may render an accessible and operable approach for directed evolution of fungal α-amylase aimed at interesting functions.