Do Homologous Thermophilic–Mesophilic Proteins Exhibit Similar Structures and Dynamics at Optimal Growth Temperatures? A Molecular Dynamics Simulation Study

Molecular dynamics simulations at high temperatures of the Aeropyrum pernix L7Ae thermostable protein: Insight into the unfolding pathway

Most life forms on earth live at temperatures below 50 °C. Within these organisms are proteins that form the three-dimensional structures essential to their biological activity and function. However, some thermophilic life forms can resist higher temperatures and have corresponding adaptations to preserve protein function at these high temperatures. Among the structural factors responsible for this resistance of thermophilic proteins to high temperatures is the presence of additional hydrogen bonds in the thermophilic proteins, which means that the structure of the protein is more resistant to unfolding. Similarly, thermostable proteins are rich in structure-stabilizing salt bridges and/or disulfide bridges. In this context, we perform multiple replica molecular dynamics simulations at different temperatures on the Aeropyrum pernix (L7Ae) protein (from the crenarchaeal species A. pernix), known for its high melting temperature, and this in the aim to elucidate the structural factors responsible for its high thermostability. The results reveal that between the most sensitive regions of the protein to the increase of temperature are the loops L1, and L5, which surround the hydrophobic core region of the protein, besides the loop L9, and the C-terminal α5 region. This latter is the longer alpha helix of the protein secondary structure motifs and it is the first to be denaturated at 450 K, while the rest of the protein secondary structure motifs at this temperature were intact. The mechanism of unfolding that follows this protein at 550 K is similar to other thermophile proteins found in literature, with the opening of the loops that surround the hydrophobic core of the protein. So, the latter is completely exposed to the solvent, and partially denatured. The total denaturation process of the protein takes an average time of 40 ns to be achieved. Our investigation also shows that all the calculated salt bridges, with distances less than or equal to 6 A°, are on the periphery part of the protein, exposed to the solvent. However, the hydrophobic core of the protein is not involved in the formation of salt bridges, but rather with formation of some important hydrogen bondings that still persist even at 450 K. So, optimizing hydrogen bonding, near or within the core region, at high temperatures is a strategy that follows this thermostable protein to protect its hydrophobic core from denaturation, and ensure the thermal stability of the protein.

Molecular Dynamics Simulations of HPr Proteins from a Thermophilic and a Mesophilic Organism: A Comparative Thermal Study

The histidine-containing phosphocarrier (HPr) is a monomeric protein conserved in Gram-positive bacteria, which may be of mesophilic or thermophilic nature. In particular, the HPr protein from the thermophilic organism B. stearothermophilus is a good model system for thermostability studies, since experimental data, such as crystal structure and thermal stability curves, are available. However, its unfolding mechanism at higher temperatures is yet unclear at a molecular level. Therefore, in this work, we researched the thermal stability of this protein using molecular dynamics simulations, subjecting it to five different temperatures during a time span of 1 μs. The analyses of the structural parameters and molecular interactions were compared with those of the mesophilic homologue HPr protein from B. subtilis. Each simulation was run in triplicate using identical conditions for both proteins. The results showed that the two proteins lose stability as the temperature increases, but the mesophilic structure is more affected. We found that the salt bridge network formed by the triad of Glu3-Lys62-Glu36 residues and the salt bridge made up of Asp79-Lys83 ion pair are key factors to keep stable the thermophilic protein, maintaining the hydrophobic core protected and the structure packed. In addition, these molecular interactions neutralize the negative surface charge, acting as "natural molecular staples".

In-silico design of peptides for inhibition of HLA-A*03-KLIETYFSK complex as a new drug design for treatment of multiples sclerosis disease

Multiple sclerosis is recognized as a chronic inflammatory disease. Human leukocyte antigen (HLA) plays an important role in initiating adaptive immune responses. HLA class I is present in almost all nucleated cells and presents the cleaved endogenous peptide antigens to cytotoxic T cells. HLA-A*03 is one of the HLA class I alleles, which is reported as substantially related HLA to MS disease. In 2011, the structure of the HLA-A*03 in complex was identified with an immunodominant proteolipid protein (PLP) epitope (KLIETYFSK). This complex has been reported as an important autoantigen-presenting complex in MS pathogenesis. In this study, new peptides were designed to bind to this complex that may prevent specific pathogenic cytotoxic T cell binding to this autoantigen-presenting complex and CNS demyelination. Herein, 14 new helical peptides containing 19 amino acids were designed and their structures were predicted using the PEP-FOLD server. The binding of each designed peptide to the mentioned complex was then performed. A mutation approach was used by the BeAtMuSiC server to improve the binding affinity of the designed peptide. In each position, amino acid substitutions leading to an increase in the binding affinity of the peptide to the mentioned complex were determined. Finally, the resulting complexes were simulated for 40 ns using AMBER18 software. The results revealed that out of 14 designed peptides, "WRYWWKDWAKQFRQFYRWF" peptide exhibited the highest affinity for binding to the mentioned complex. This peptide can be considered as a potential drug to control multiple sclerosis disease in patients carrying the HLA-A*03 allele.

Temperature-resistant and solvent-tolerant lipases as industrial biocatalysts: Biotechnological approaches and applications

The development of new biocatalytic systems to replace the chemical catalysts, with suitable characteristics in terms of efficiency, stability under high temperature reactions and in the presence of organic solvents, reusability, and eco-friendliness is considered a very important step to move towards the green processes. From this basis, the use of lipase as a catalyst is highly desired for many industrial applications because it offers the reactions in which could be used, stability in harsh conditions, reusability and a greener process. Therefore, the introduction of temperature-resistant and solvent-tolerant lipases have become essential and ideal for industrial applications. Temperature-resistant and solvent-tolerant lipases have been involved in many large-scale applications including biodiesel, detergent, food, pharmaceutical, organic synthesis, biosensing, pulp and paper, textile, animal feed, cosmetics, and leather industry. So, the present review provides a comprehensive overview of the industrial use of lipase. Moreover, special interest in biotechnological and biochemical techniques for enhancing temperature-resistance and solvent-tolerance of lipases to be suitable for the industrial uses.

A novel hyperthermophilic methylglyoxal synthase: molecular dynamic analysis on the regional fluctuations

Two putative methylglyoxal synthases, which catalyze the conversion of dihydroxyacetone phosphate to methylglyoxal, from Oceanithermus profundus DSM 14,977 and Clostridium difficile 630 have been characterized for activity and thermal stability. The enzyme from O. profundus was found to be hyperthermophilic, with the optimum activity at 80 °C and the residual activity up to 59% after incubation of 15 min at 95 °C, whereas the enzyme from C. difficile was mesophilic with the optimum activity at 40 °C and the residual activity less than 50% after the incubation at 55 °C or higher temperatures for 15 min. The structural analysis of the enzymes with molecular dynamics simulation indicated that the hyperthermophilic methylglyoxal synthase has a rigid protein structure with a lower overall root-mean-square-deviation value compared with the mesophilic or thermophilic counterparts. In addition, the simulation results identified distinct regions with high fluctuations throughout those of the mesophilic or thermophilic counterparts via root-mean-square-fluctuation analysis. Specific molecular interactions focusing on the hydrogen bonds and salt bridges in the distinct regions were analyzed in terms of interatomic distances and positions of the individual residues with respect to the secondary structures of the enzyme. Key interactions including specific salt bridges and hydrogen bonds between a rigid beta-sheet core and surrounding alpha helices were found to contribute to the stabilisation of the hyperthermophilic enzyme by reducing the regional fluctuations in the protein structure. The structural information and analysis approach in this study can be further exploited for the engineering and industrial application of the enzyme.

On the Thermal Stability of O6-Methylguanine-DNA Methyltransferase from Archaeon Pyrococcus kodakaraensis by Molecular Dynamics Simulations

We have employed molecular dynamics simulations to analyze the thermal stability of the O⁶-methylguanine-DNA methyltransferase (MGMT) protein, both hyperthermophilic archaeon Pyrococcus kodakaraensis (Pk-MGMT) and its mesophilic homologue pair, obtained from enterobacterium Escherichia coli (AdaC). This theoretical study was done at three different temperatures: 302, 371, and 450 K. The molecular dynamics has been performed in explicit aqueous solvent during a period of time of 95 ns, including periodic boundary conditions and constant pressure. The same procedure has been used for both proteins, and each simulation has been carried out by triplicate. Hence, we performed 18 simulations. In this way, we have done different analyses to explore the factors that may affect the thermal stability of Pk-MGMT. The structural behavior was analyzed using indicators such as root-mean-square deviation, radius of gyration, solvent-accessible surface area, hydrogen bonds, native contacts, secondary structure, and salt bridge formation. The results showed that when the temperature increases, the global atomic fluctuations increase too, which suggests that both proteins lose thermal stability, but as expected, this fact is highlighted in AdaC. Moreover, the contacts of the native state in AdaC are considerably lower than those found in Pk-MGMT at 450 K. Also, the structural studies showed that conserved and nonconserved salt bridges kept close contacts with the Pk-MGMT protein at high temperatures. These interaction types act as molecular staples and are mainly responsible to provide thermostability to the hyperthermophilic protein.

In silico study of the association of the HLA-A*31:01 allele (human leucocyte antigen allele 31:01) with neuroantigenic epitopes of PLP (proteolipid protein), MBP (myelin basic protein) and MOG proteins (myelin oligodendrocyte glycoprotein) for studying the multiple sclerosis disease pathogenesis

The main pathologic hallmark of multiple sclerosis is a demyelinating plaque that contains a prominent immunologic response dominated by T cells of the immune system. PLP (proteolipid protein), MPB (myelin basic protein), and Myelin oligodendrocyte glycoprotein (MOG) proteins are important autoantigens for the demyelinating of CNS in multiple sclerosis. There is good evidence indicating that T CD8⁺ cells and MHC class I molecules play an important role in this disease. The HLA-A*31:01 allele of MHC class I is a member of HLA-A3 superfamily and there is no clear report concerning the relationship of this allele with MS. Feeling this gap, we studied the possible association of the HLA-A*31:01 with MS by prediction of neuroantigenic epitopes of human MBP, PLP, and MOG proteins of myelin sheath using in silico methods. PLP did not show any neuroantigenic epitope, but the two epitopes of MBP and seven epitopes of MOG for HLA-A*31:01 were determined via bioinformatics servers. In silico study of the nine epitope showed that MOG195-204 (LIICYNWLHR) peptide of the membrane-associated/cytoplasmic part of human MOG has suitable binding affinity to the HLA-A*31:01 allele as a potential neuroantigenic epitope. Further investigations of this peptide revealed that the binding of C-terminal residue of this peptide has a more significant effect on binding to this allele than the N-terminal part of the peptide. Altogether, this combination of "LIICYNWLHR/A*31:01 allele "may play an important role in MS pathogenesis and this complex is suggested for further studies such as T cell receptor.Communicated by Ramaswamy H. Sarma.

Improving the solubility, activity, and stability of reteplase using in silico design of new variants

Reteplase (recombinant plasminogen activator, r-PA) is a thrombolytic agent recombined from tissue-type plasminogen activator (t-PA), which has several prominent features such as strong thrombolytic ability and E. coli expressibility. Despite these outstanding features, it demonstrates reduced fibrin binding affinity, reduced stimulation of protease activity, and lower solubility, hence higher aggregation propensity, compared to t-PA. The present study was devoted to design r-PA variants with comparable structural stability, enhanced biological activity, and high solubility. For this purpose, computational molecular modeling techniques were utilized. The supercharging technique was applied for r-PA to designing new species of the protein. Based on the results from in silico evaluation of selected mutations in comparison to the wild-type r-PA, the designed supercharged mutant (S7 variant) exhibited augmented stability, decreased solvation energy, as well as enhanced binding affinity to fibrin. The data also implied increased plasminogen cleavage activity of the new variant. These findings have implications to therapies which involve removal of intravascular blood clots, including the treatment of acute myocardial infarction.

Temperature dependent dynamics in highly homologous adenylate kinases

To correlate the structural features of enzymes to temperature adaptation, we studied psychrophile, mesophile, and thermophile adenylate kinases as model enzymes using bioinformatics and computational tools. Phylogenetic analysis revealed that mesophile and thermophile variants are clustered in one stem of phylogenetic tree and are close to contemporary time, while psychrophile enzyme is more close to their common ancestor. This finding is in good agreement with the process of environmental changes from ice age toward current warm conditions on the earth. We also performed Molecular Dynamics simulation at corresponding temperatures of all enzyme variants including 308, 318, and 328 K. It was found that mesophile enzyme has no distinct deviation of Root Mean Square Deviation (RMSD) and Radius of Gyration (Rg) values from equilibrium states at operating temperature of thermophile enzyme as well as its own optimum temperature. However, psychrophile enzyme undergoes more fluctuations with higher amplitude of change; particularly at 328 K. It was also found that initial increasing of RMSD and Rg for Psychrophile enzyme at all temperatures is occurred gradually; while, the increment of this structural parameters for thermophile enzyme at 328 K is occurred in a highly cooperative and switching manner demonstrating snap structural change of thermophile enzyme in its own temperature. By analysis of Root Mean Square Fluctuation values at different temperatures, we identified two flexible fragments in adenylate kinases so that different dynamic behavior of these regions in mesophile enzyme against operating temperatures of psychrophile and thermophile variants is critical in compensation of flexibility challenges at respective temperatures.

Computational comparison of two new fusion proteins for multiple sclerosis

Multiple sclerosis (MS), as one of the human autoimmune diseases, demyelinates the neurons of the central nervous system (CNS). Activation of the T cells which target the CNS antigens is the first autoimmune event in MS. Myelin oligodendrocyte glycoprotein (MOG) and myelin basic protein (MBP) are two proteins of the myelin sheath and have been shown to be among the high antigens contributing to the pathogenesis of MS. Production of the drugs with high specificity for the immune system diseases is a concern for various researchers. Therefore, tolerogenic vaccines are considered as a new strategy for the treatment of MS by presenting specific antigens. This study aimed to design and compare two fusion proteins by a combination of two neuroantigens linked to interleukin-16 (IL-16) (MOG-Linker-MBP-IL16 and MBP-Linker-MOG-IL16) as vaccines for MS. In this study, at first two models MOG (aa 11-30) linked to MBP (aa 13-32) was made by Modeler 9.10 and simulated for 20 ns via Gromacs 5.1.1 package. Then simulated antigen domains connected to the N-terminal domain of IL-16 and obtained structures simulated for 50 ns. The results revealed that both constructs had stable structures and the linker could keep two antigenic fragments separate enough, preventing undesired interactions. While MOG-Linker-MBP-IL16 showed better solubility, more accessible surface areas, more flexibility of its IL-16 domain, and better functionality of its IL-16 domain as well as more specific cleavage of its related epitopes after endocytosis lead to a better presentation of its antigenic property.

Discovery of an intermolecular disulfide bond required for the thermostability of a heterodimeric protein from the thermophile Hydrogenobacter thermophilus

Factors that increase protein thermostability are of considerable interest in both scientific and industrial fields. Disulfide bonds are one of such factors that increase thermostability, but are rarely found in intracellular proteins because of the reducing environment of the cytosol. Here, we report the first example of an intermolecular disulfide bond between heteromeric subunits of a novel-type phosphoserine phosphatase from a thermophilic bacterium Hydrogenobacter thermophilus, which contributes to the protein thermostability at the physiological temperature. Comparison of remaining soluble proteins between wild-type and cysteine-deleted mutant using SDS-PAGE revealed that the disulfide bond increases the thermostability of the whole protein by tightly connecting a subunit with low solubility to the partner with higher solubility. Furthermore, it was strongly suggested that the disulfide bond is formed and contributes to the stability in vivo. This finding will open new avenues for the design of proteins with increased thermostability.

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.

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.

Interface matters: the stiffness route to stability of a thermophilic tetrameric malate dehydrogenase

In this work we investigate by computational means the behavior of two orthologous bacterial proteins, a mesophilic and a thermophilic tetrameric malate dehydrogenase (MalDH), at different temperatures. Namely, we quantify how protein mechanical rigidity at different length- and time-scales correlates to protein thermophilicity as commonly believed. In particular by using a clustering analysis strategy to explore the conformational space of the folded proteins, we show that at ambient conditions and at the molecular length-scale the thermophilic variant is indeed more rigid that the mesophilic one. This rigidification is the result of more efficient inter-domain interactions, the strength of which is further quantified via ad hoc free energy calculations. When considered isolated, the thermophilic domain is indeed more flexible than the respective mesophilic one. Upon oligomerization, the induced stiffening of the thermophilic protein propagates from the interface to the active site where the loop, controlling the access to the catalytic pocket, anchors down via an extended network of ion-pairs. On the contrary in the mesophilic tetramer the loop is highly mobile. Simulations at high temperature, could not re-activate the mobility of the loop in the thermophile. This finding opens questions on the similarities of the binding processes for these two homologues at their optimal working temperature and suggests for the thermophilic variant a possible cooperative role of cofactor/substrate.

How conformational flexibility stabilizes the hyperthermophilic elongation factor G-domain

Proteins from thermophilic organisms are stable and functional well above ambient temperature. Understanding the molecular mechanism underlying such a resistance is of crucial interest for many technological applications. For some time, thermal stability has been assumed to correlate with high mechanical rigidity of the protein matrix. In this work we address this common belief by carefully studying a pair of homologous G-domain proteins, with their melting temperatures differing by 40 K. To probe the thermal-stability content of the two proteins we use extensive simulations covering the microsecond time range and employ several different indicators to assess the salient features of the conformational landscape and the role of internal fluctuations at ambient condition. At the atomistic level, while the magnitude of fluctuations is comparable, the distribution of flexible and rigid stretches of amino-acids is more regular in the thermophilic protein causing a cage-like correlation of amplitudes along the sequence. This caging effect is suggested to favor stability at high T by confining the mechanical excitations. Moreover, it is found that the thermophilic protein, when folded, visits a higher number of conformational substates than the mesophilic homologue. The entropy associated with the occupation of the different substates and the thermal resilience of the protein intrinsic compressibility provide a qualitative insight on the thermal stability of the thermophilic protein as compared to its mesophilic homologue. Our findings potentially open the route to new strategies in the design of thermostable proteins.