1
|
Sheikhzadeh A, Safaei M, Fadaei Naeini V, Baghani M, Foroutan M, Baniassadi M. Multiscale modeling of unfolding and bond dissociation of rubredoxin metalloprotein. J Mol Graph Model 2024; 129:108749. [PMID: 38442439 DOI: 10.1016/j.jmgm.2024.108749] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2023] [Revised: 02/20/2024] [Accepted: 02/20/2024] [Indexed: 03/07/2024]
Abstract
Mechanical properties of proteins that have a crucial effect on their operation. This study used a molecular dynamics simulation package to investigate rubredoxin unfolding on the atomic scale. Different simulation techniques were applied, and due to the dissociation of covalent/hydrogen bonds, this protein demonstrates several intermediate states in force-extension behavior. A conceptual model based on the cohesive finite element method was developed to consider the intermediate damages that occur during unfolding. This model is based on force-displacement curves derived from molecular dynamics results. The proposed conceptual model is designed to accurately identify bond rupture points and determine the associated forces. This is achieved by conducting a thorough comparison between molecular dynamics and cohesive finite element results. The utilization of a viscoelastic cohesive zone model allows for the consideration of loading rate effects. This rate-dependent model can be further developed and integrated into the multiscale modeling of large assemblies of metalloproteins, providing a comprehensive understanding of mechanical behavior while maintaining a reduced computational cost.
Collapse
Affiliation(s)
- Aliakbar Sheikhzadeh
- School of Mechanical Engineering, College of Engineering, University of Tehran, Tehran, Iran; Department of Chemical and Materials Engineering, University of Alberta, Edmonton, AB, T6G 2V4, Canada
| | - Mohammad Safaei
- School of Mechanical Engineering, College of Engineering, University of Tehran, Tehran, Iran
| | - Vahid Fadaei Naeini
- Division of Machine Elements, Luleå University of Technology, Luleå, SE-97187, Sweden
| | - Mostafa Baghani
- School of Mechanical Engineering, College of Engineering, University of Tehran, Tehran, Iran
| | - Masumeh Foroutan
- Department of Physical Chemistry, School of Chemistry, College of Science, University of Tehran, Tehran, Iran.
| | - Majid Baniassadi
- School of Mechanical Engineering, College of Engineering, University of Tehran, Tehran, Iran; University of Strasbourg, CNRS, ICUBE Laboratory, Strasbourg, France.
| |
Collapse
|
2
|
Ghaedizadeh S, Zeinali M, Dabirmanesh B, Rasekh B, Khajeh K, Banaei-Moghaddam AM. Rational design engineering of a more thermostable Sulfurihydrogenibium yellowstonense carbonic anhydrase for potential application in carbon dioxide capture technologies. BIOCHIMICA ET BIOPHYSICA ACTA. PROTEINS AND PROTEOMICS 2024; 1872:140962. [PMID: 37716447 DOI: 10.1016/j.bbapap.2023.140962] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/19/2023] [Revised: 08/18/2023] [Accepted: 09/06/2023] [Indexed: 09/18/2023]
Abstract
Implementing hyperthermostable carbonic anhydrases into CO2 capture and storage technologies in order to increase the rate of CO2 absorption from the industrial flue gases is of great importance from technical and economical points of view. The present study employed a combination of in silico tools to further improve thermostability of a known thermostable carbonic anhydrase from Sulfurihydrogenibium yellowstonense. Experimental results showed that our rationally engineered K100G mutant not only retained the overall structure and catalytic efficiency but also showed a 3 °C increase in the melting temperature and a two-fold improvement in the enzyme half-life at 85 °C. Based on the molecular dynamics simulation results, rearrangement of salt bridges and hydrogen interactions network causes a reduction in local flexibility of the K100G variant. In conclusion, our study demonstrated that thermostability can be improved through imposing local structural rigidity by engineering a single-point mutation on the surface of the enzyme.
Collapse
Affiliation(s)
- Shima Ghaedizadeh
- Laboratory of Genomics and Epigenomics (LGE), Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - Majid Zeinali
- Microbiology and Biotechnology Research Group, Research Institute of Petroleum Industry (RIPI), Tehran, Iran.
| | - Bahareh Dabirmanesh
- Department of Biochemistry, Faculty of Biological Sciences, Tarbiat Modares University, Tehran, Iran
| | - Behnam Rasekh
- Microbiology and Biotechnology Research Group, Research Institute of Petroleum Industry (RIPI), Tehran, Iran
| | - Khosrow Khajeh
- Department of Biochemistry, Faculty of Biological Sciences, Tarbiat Modares University, Tehran, Iran
| | - Ali Mohammad Banaei-Moghaddam
- Laboratory of Genomics and Epigenomics (LGE), Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran.
| |
Collapse
|
3
|
Zhang X, Li L, Zheng Q. Dissecting the Effect of Temperature on Hyperthermophilic Pf2001 Esterase Dimerization by Molecular Dynamics. J Chem Inf Model 2023; 63:4762-4771. [PMID: 37452749 DOI: 10.1021/acs.jcim.3c00415] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/18/2023]
Abstract
Pf2001 esterase (Pf2001) from Pyrococcus furiosus has hyperthermophilic properties and exerts a biocatalytic function in a dimeric state. Crystal structures revealed that the structural rearrangement of the cap domain is responsible for the Pf2001 dimer formation. However, the details of the cap domain remodeling and the effects of temperature on the dimerization process remain elusive at the molecular level, taking into account that experimental methods are difficult to capture the dynamic process of dimerization to some extent. Herein, four dimer models based on the monomeric crystal structure (PDB ID: 5G59) were constructed to investigate the conformational transition details and temperature effects in the dimerization by conventional molecular dynamics and accelerated molecular dynamics simulations. Our simulation results indicate that the monomer undergoes a conformational change into a "preparatory state" at high temperatures, which is more favorable for its transformation into a stable dimer. The subsequent free energy landscape analysis further identifies four intermediate states (from separated state to dimeric state) and discloses that a more accessible α-helix driven by stronger hydrophobic interactions induces a rearrangement of the cap domain, displaying a "tic-tac-toe" activation feature that is important for stabilizing the dimer interface and facilitating the formation of hydrophobic pockets. In addition, the electrostatic potential surface analysis illustrates that the weaker electrostatic repulsion (Lys and Arg) in the dimer interface at high temperatures is also a key factor for dimer stabilization. Altogether, our results can provide molecular-level insight into the dimer formation process of hyperthermophilic esterase and would be useful to understand the enzymatic specificity of α/β-hydrolase.
Collapse
Affiliation(s)
- Xue Zhang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China
| | - Lei Li
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China
| | - Qingchuan Zheng
- School of Pharmaceutical Sciences, Jilin University, Changchun 130021, China
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China
| |
Collapse
|
4
|
Okazaki K, Kobashigawa Y, Morita H, Yamauchi S, Fukuda N, Liu C, Toyota Y, Sato T, Morioka H. Molecular Dynamics-Based Design and Biophysical Evaluation of Thermostable Single-Chain Fv Antibody Mutants Derived from Pharmaceutical Antibodies. ACS OMEGA 2023; 8:22945-22954. [PMID: 37396255 PMCID: PMC10308585 DOI: 10.1021/acsomega.3c01948] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/23/2023] [Accepted: 05/26/2023] [Indexed: 07/04/2023]
Abstract
Antibody drugs are denatured under physical stress, e.g., friction, heat, and freezing, which triggers formation of aggregates and resultant allergic reactions. Design of a stable antibody is thus critical for the development of antibody drugs. Here, we obtained a thermostable single-chain Fv (scFv) antibody clone by rigidifying the flexible region. We first conducted a short molecular dynamics (MD) simulation (3 runs of 50 ns) to search for weak spots in the scFv antibody, i.e., flexible regions located outside the CDR (complementarity determining region) and the interface between the heavy-chain and light-chain variable regions. We then designed a thermostable mutant and evaluated it by means of a short MD simulation (3 runs of 50 ns) based on reductions in the root-mean-square fluctuation (RMSF) values and formation of new hydrophilic interactions around the weak spot. Finally, we designed the VL-R66G mutant by applying our strategy to scFv derived from trastuzumab. Trastuzumab scFv variants were prepared by using an Escherichia coli expression system, and the melting temperature-measured as a thermostability index-was 5 °C higher than that of the wild-type trastuzumab scFv, while the antigen-binding affinity was unchanged. Our strategy required few computational resources, and would be applicable to antibody drug discovery.
Collapse
|
5
|
Gómez-Flores AK, López-Pérez E, Alas-Guardado SJ. Molecular Dynamics Simulations of HPr Proteins from a Thermophilic and a Mesophilic Organism: A Comparative Thermal Study. Int J Mol Sci 2023; 24:ijms24119557. [PMID: 37298508 DOI: 10.3390/ijms24119557] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2023] [Revised: 05/23/2023] [Accepted: 05/24/2023] [Indexed: 06/12/2023] Open
Abstract
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".
Collapse
Affiliation(s)
- Ana K Gómez-Flores
- Departamento de Ciencias Naturales, Universidad Autónoma Metropolitana Unidad Cuajimalpa, Mexico City 05300, Mexico
| | - Edgar López-Pérez
- Posgrado en Ciencias Naturales e Ingeniería, Universidad Autónoma Metropolitana Unidad Cuajimalpa, Mexico City 05300, Mexico
| | - Salomón J Alas-Guardado
- Departamento de Ciencias Naturales, Universidad Autónoma Metropolitana Unidad Cuajimalpa, Mexico City 05300, Mexico
| |
Collapse
|
6
|
Rao D, Huo R, Yan Z, Guo Z, Liu W, Lu M, Luo H, Tao X, Yang W, Su L, Chen S, Wang L, Wu J. Multiple approaches of loop region modification for thermostability improvement of 4,6-α-glucanotransferase from Limosilactobacillus fermentum NCC 3057. Int J Biol Macromol 2023; 233:123536. [PMID: 36740130 DOI: 10.1016/j.ijbiomac.2023.123536] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2022] [Revised: 01/17/2023] [Accepted: 01/31/2023] [Indexed: 02/05/2023]
Abstract
4,6-α-glucanotransferase (4,6-α-GT), as a member of the glycoside hydrolase 70 (GH70) family, converts starch/maltooligosaccharides into α,1-6 bond-containing α-glucan and possesses potential applications in food, medical and related industries but does not satisfy the high-temperature requirement due to its poor thermostability. In this study, a 4,6-α-GT (ΔGtfB) from Limosilactobacillus fermentum NCC 3057 was used as a model enzyme to improve its thermostability. The loops of ΔGtfB as the target region were optimized using directed evolution, sequence alignment, and computer-aided design. A total of 11 positive mutants were obtained and iteratively combined to obtain a combined mutant CM9, with high resistance to temperature (50 °C). The activity of mutant CM9 was 2.08-fold the activity of the wild type, accompanied by a 5 °C higher optimal temperature, a 5.76 °C higher melting point (Tm, 59.46 °C), and an 11.95-fold longer half-life time (t1/2). The results showed that most of the polar residues in the loop region of ΔGtfB are mutated into rigid proline residues. Molecular dynamics simulation demonstrated that the root mean square fluctuation of CM9 significantly decreased by "Breathing" movement reduction of the loop region. This study provides a new strategy for improving the thermostability of 4,6-α-GT through rational loop region modification.
Collapse
Affiliation(s)
- Deming Rao
- State Key Laboratory of Food Science and Technology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; School of Biotechnology and Key Laboratory of Industrial Biotechnology Ministry of Education, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; International Joint Laboratory on Food Safety, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China
| | - Runtian Huo
- State Key Laboratory of Food Science and Technology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; School of Biotechnology and Key Laboratory of Industrial Biotechnology Ministry of Education, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; International Joint Laboratory on Food Safety, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China
| | - Zhengfei Yan
- State Key Laboratory of Food Science and Technology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; School of Biotechnology and Key Laboratory of Industrial Biotechnology Ministry of Education, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; International Joint Laboratory on Food Safety, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China
| | - Zhiyong Guo
- State Key Laboratory of Food Science and Technology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; School of Biotechnology and Key Laboratory of Industrial Biotechnology Ministry of Education, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; International Joint Laboratory on Food Safety, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China
| | - Weiqiong Liu
- State Key Laboratory of Food Science and Technology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; School of Biotechnology and Key Laboratory of Industrial Biotechnology Ministry of Education, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; International Joint Laboratory on Food Safety, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China
| | - Mengwei Lu
- State Key Laboratory of Food Science and Technology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; School of Biotechnology and Key Laboratory of Industrial Biotechnology Ministry of Education, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; International Joint Laboratory on Food Safety, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China
| | - Hui Luo
- State Key Laboratory of Food Science and Technology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; School of Biotechnology and Key Laboratory of Industrial Biotechnology Ministry of Education, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; International Joint Laboratory on Food Safety, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China
| | - Xiumei Tao
- State Key Laboratory of Food Science and Technology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; School of Biotechnology and Key Laboratory of Industrial Biotechnology Ministry of Education, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; International Joint Laboratory on Food Safety, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China
| | - Weikang Yang
- State Key Laboratory of Food Science and Technology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; School of Biotechnology and Key Laboratory of Industrial Biotechnology Ministry of Education, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; International Joint Laboratory on Food Safety, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China
| | - Lingqia Su
- State Key Laboratory of Food Science and Technology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; School of Biotechnology and Key Laboratory of Industrial Biotechnology Ministry of Education, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; International Joint Laboratory on Food Safety, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China
| | - Sheng Chen
- State Key Laboratory of Food Science and Technology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; School of Biotechnology and Key Laboratory of Industrial Biotechnology Ministry of Education, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; International Joint Laboratory on Food Safety, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China
| | - Lei Wang
- State Key Laboratory of Food Science and Technology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; School of Biotechnology and Key Laboratory of Industrial Biotechnology Ministry of Education, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; International Joint Laboratory on Food Safety, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China.
| | - Jing Wu
- State Key Laboratory of Food Science and Technology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; School of Biotechnology and Key Laboratory of Industrial Biotechnology Ministry of Education, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China; International Joint Laboratory on Food Safety, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China.
| |
Collapse
|
7
|
Rahban M, Zolghadri S, Salehi N, Ahmad F, Haertlé T, Rezaei-Ghaleh N, Sawyer L, Saboury AA. Thermal stability enhancement: Fundamental concepts of protein engineering strategies to manipulate the flexible structure. Int J Biol Macromol 2022; 214:642-654. [DOI: 10.1016/j.ijbiomac.2022.06.154] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2022] [Revised: 06/22/2022] [Accepted: 06/23/2022] [Indexed: 01/28/2023]
|
8
|
Qu Y, Davey K, Sun Y, Middelberg A, Bi J. Engineered Design of the E-Helix Structure on Ferritin Nanoparticles. ACS APPLIED BIO MATERIALS 2022; 5:3167-3179. [PMID: 35770389 DOI: 10.1021/acsabm.2c00154] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Insertion of an immunogenic epitope at the C-terminus of ferritin has shown the potential to produce a stable and efficacious vaccine. There is however limited understanding of how C-terminus insertion affects ferritin protein stability. The E-helix at the C-terminus has attracted interest because there are contradictory reports as to whether it has a role in protein stabilization. Here, we report, for the first time, combining molecular dynamics simulation (MDS) with experiment to engineer the design of the E-helix at the C-terminus of engineered human ferritin heavy chain (F1) inserted with Epstein-Barr nuclear antigen 1 (EBNA1, E1) and flexible linker (L3) residues (to afford F1L3E1). Hot spots on the E-helix of the C-terminus were predicted by MDS at aa 167 (Glu) and aa 171 (Asp). Five (5) variants of F1L3E1 were constructed by considering hot spots and alteration of electrostatic or hydrophobic interfaces, namely, (1) C1, hot spots substituted with noncharged residue Gln; (2) C2, hot spots substituted with positively charged residue Arg; (3) C3, hydrophobic residues substituted with the most hydrophobic residues Val and Ile; (4) C4, hydrophobic residues substituted with the most hydrophilic residues Gln and Asn; and (5) C5, a heptad repeat structure in the E-helix disrupted by substituting "a" and "d" heptad residues with noncharged polar residue Gln. It was found that the E-helix is essential to maintain integrated protein stability and that changing the hydrophobic interface (C3 and C4) had more significant effects on protein folding and stability than changing the electrostatic interface (C1 and C2). It was confirmed by both MDS and experiment that variants C1, C2, and C5 were able to fold to form stable conformational structures with protein surface hydrophobicity similar to that of F1L3E1. However, they are less thermally stable than F1L3E1. Significant changes in hydrophobicity drove significant protein aggregation for variants C3 and C4. It is concluded that the molecular design of the C-terminus in engineered ferritin, especially the E-helix, is important to ensure the epitope-based chimeric vaccine is safe (aggregate free) and efficacious.
Collapse
Affiliation(s)
- Yiran Qu
- School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, SA 5005, Australia
| | - Kenneth Davey
- School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, SA 5005, Australia
| | - Yan Sun
- Department of Biochemical Engineering and Key Laboratory of Systems Bioengineering of the Ministry of Education, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
| | - Anton Middelberg
- Division of Research and Innovation, The University of Adelaide, Adelaide, SA 5005, Australia
| | - Jingxiu Bi
- School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, SA 5005, Australia
| |
Collapse
|
9
|
Palaniappan C, Narayanan RC, Sekar K. Mutation-Dependent Refolding of Prion Protein Unveils Amyloidogenic-Related Structural Ramifications: Insights from Molecular Dynamics Simulations. ACS Chem Neurosci 2021; 12:2810-2819. [PMID: 34296847 DOI: 10.1021/acschemneuro.1c00142] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
The main focus of prion structural biology studies is to understand the molecular basis of prion diseases caused by misfolding, and aggregation of the cellular prion protein PrPC remains elusive. Several genetic mutations are linked with human prion diseases and driven by the conformational conversion of PrPC to the toxic PrPSc. The main goal of this study is to gain a better insight into the molecular effect of disease-associated V210I mutation on this process by molecular dynamics simulations. This inherited mutation elicited copious structural changes in the β1-α1-β2 subdomain, including an unfolding of a helix α1 and the elongation of the β-sheet. These unusual structural changes likely appeared to detach the β1-α1-β2 subdomain from the α2-α3 core, an early misfolding event necessary for the conformational conversion of PrPC to PrPSc. Ultimately, the unfolded α1 and its prior β1-α1 loop further engaged with unrestrained conformational dynamics and were widely considered as amyloidogenic-inducing traits. Furthermore, the resulting folding intermediate possesses a highly unstable β1-α1-β2 subdomain, thereby enhancing the aggregation of misfolded PrPC through intermolecular interactions between frequently refolding regions. Briefly, these remarkable changes as seen in the mutant β1-α1-β2 subdomain are consistent with previous experimental results and thus provide a molecular basis of PrPC misfolding associated with the conformational conversion of PrPC to PrPSc.
Collapse
Affiliation(s)
| | - Rahul C. Narayanan
- Department of Computational and Data Sciences, Indian Institute of Science, Bangalore 560 012, India
| | - Kanagaraj Sekar
- Department of Computational and Data Sciences, Indian Institute of Science, Bangalore 560 012, India
| |
Collapse
|
10
|
Arnittali M, Rissanou AN, Amprazi M, Kokkinidis M, Harmandaris V. Structure and Thermal Stability of wtRop and RM6 Proteins through All-Atom Molecular Dynamics Simulations and Experiments. Int J Mol Sci 2021; 22:ijms22115931. [PMID: 34073028 PMCID: PMC8199364 DOI: 10.3390/ijms22115931] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2021] [Revised: 05/23/2021] [Accepted: 05/25/2021] [Indexed: 01/07/2023] Open
Abstract
In the current work we study, via molecular simulations and experiments, the folding and stability of proteins from the tertiary motif of 4-α-helical bundles, a recurrent motif consisting of four amphipathic α-helices packed in a parallel or antiparallel fashion. The focus is on the role of the loop region in the structure and the properties of the wild-type Rop (wtRop) and RM6 proteins, exploring the key factors which can affect them, through all-atom molecular dynamics (MD) simulations and supporting by experimental findings. A detailed investigation of structural and conformational properties of wtRop and its RM6 loopless mutation is presented, which display different physical characteristics even in their native states. Then, the thermal stability of both proteins is explored showing RM6 as more thermostable than wtRop through all studied measures. Deviations from native structures are detected mostly in tails and loop regions and most flexible residues are indicated. Decrease of hydrogen bonds with the increase of temperature is observed, as well as reduction of hydrophobic contacts in both proteins. Experimental data from circular dichroism spectroscopy (CD), are also presented, highlighting the effect of temperature on the structural integrity of wtRop and RM6. The central goal of this study is to explore on the atomic level how a protein mutation can cause major changes in its physical properties, like its structural stability.
Collapse
Affiliation(s)
- Maria Arnittali
- Institute of Applied and Computational Mathematics (IACM), Foundation for Research and Technology Hellas (FORTH), IACM/FORTH, GR-71110 Heraklion, Crete, Greece; (M.A.); (V.H.)
- Department of Mathematics and Applied Mathematics, University of Crete, GR-71409 Heraklion, Crete, Greece
| | - Anastassia N. Rissanou
- Institute of Applied and Computational Mathematics (IACM), Foundation for Research and Technology Hellas (FORTH), IACM/FORTH, GR-71110 Heraklion, Crete, Greece; (M.A.); (V.H.)
- Department of Mathematics and Applied Mathematics, University of Crete, GR-71409 Heraklion, Crete, Greece
- Correspondence: ; Tel.: +30-2810-393746
| | - Maria Amprazi
- Department of Biology, University of Crete, GR-71409 Heraklion, Crete, Greece; (M.A.); (M.K.)
- Institute of Molecular Biology and Biotechnology, Foundation of Research and Technology, GR-70013 Heraklion, Crete, Greece
| | - Michael Kokkinidis
- Department of Biology, University of Crete, GR-71409 Heraklion, Crete, Greece; (M.A.); (M.K.)
- Institute of Molecular Biology and Biotechnology, Foundation of Research and Technology, GR-70013 Heraklion, Crete, Greece
| | - Vagelis Harmandaris
- Institute of Applied and Computational Mathematics (IACM), Foundation for Research and Technology Hellas (FORTH), IACM/FORTH, GR-71110 Heraklion, Crete, Greece; (M.A.); (V.H.)
- Department of Mathematics and Applied Mathematics, University of Crete, GR-71409 Heraklion, Crete, Greece
- Computation-Based Science and Technology Research Center, The Cyprus Institute, 2121 Nicosia, Cyprus
| |
Collapse
|
11
|
Qu Y, Wang L, Yin S, Zhang B, Jiao Y, Sun Y, Middelberg A, Bi J. Stability of Engineered Ferritin Nanovaccines Investigated by Combined Molecular Simulation and Experiments. J Phys Chem B 2021; 125:3830-3842. [PMID: 33825471 DOI: 10.1021/acs.jpcb.1c00276] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Human ferritin is regarded as an attractive and promising vaccine platform because of its uniform structure, good plasticity, and desirable thermal and chemical stabilities. Besides, it is biocompatible and presumed safe when used as a vaccine carrier. However, there is a lack of knowledge of how different antigen insertion sites on the ferritin nanocage impact the resulting protein stability and performance. To address this question, we selected Epstein-Barr nuclear antigen 1 as a model epitope and fused it at the DNA level with different insertion sites, namely, the N- and C-termini of ferritin, to engineer proteins E1F1 and F1E1, respectively. Protein properties including hydrophobicity and thermal, pH, and chemical stability were investigated both by molecular dynamics (MD) simulation and by experiments. Both methods demonstrate that the insertion site plays an important role in protein properties. The C-terminus insertion (F1E1) leads to a less hydrophobic surface and more tolerance to the external influence of high temperature, pH, and high concentration of chemical denaturants compared to N-terminus insertion (E1F1). Simulated protein hydrophobicity and thermal stability by MD were in high accordance with experimental results. Thus, MD simulation can be used as a valuable tool to engineer nanovaccine candidates, cutting down costs by reducing the experimental effort and accelerating vaccine design.
Collapse
Affiliation(s)
- Yiran Qu
- School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, South Australia 5005, Australia
| | - Lijie Wang
- Department of Biochemical Engineering and Key Laboratory of Systems Bioengineering of the Ministry of Education, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
| | - Shuang Yin
- School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, South Australia 5005, Australia
| | - Bingyang Zhang
- School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, South Australia 5005, Australia
| | - Yan Jiao
- School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, South Australia 5005, Australia
| | - Yan Sun
- Department of Biochemical Engineering and Key Laboratory of Systems Bioengineering of the Ministry of Education, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
| | - Anton Middelberg
- School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, South Australia 5005, Australia
| | - Jingxiu Bi
- School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, South Australia 5005, Australia
| |
Collapse
|
12
|
Xu C, Suo H, Xue Y, Qin J, Chen H, Hu Y. Experimental and theoretical evidence of enhanced catalytic performance of lipase B from Candida antarctica acquired by the chemical modification with amino acid ionic liquids. MOLECULAR CATALYSIS 2021. [DOI: 10.1016/j.mcat.2020.111355] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
13
|
Miotto M, Olimpieri PP, Di Rienzo L, Ambrosetti F, Corsi P, Lepore R, Tartaglia GG, Milanetti E. Insights on protein thermal stability: a graph representation of molecular interactions. Bioinformatics 2020; 35:2569-2577. [PMID: 30535291 PMCID: PMC6662296 DOI: 10.1093/bioinformatics/bty1011] [Citation(s) in RCA: 58] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2018] [Revised: 10/29/2018] [Accepted: 12/07/2018] [Indexed: 11/14/2022] Open
Abstract
Motivation Understanding the molecular mechanisms of thermal stability is a challenge in protein biology. Indeed, knowing the temperature at which proteins are stable has important theoretical implications, which are intimately linked with properties of the native fold, and a wide range of potential applications from drug design to the optimization of enzyme activity. Results Here, we present a novel graph-theoretical framework to assess thermal stability based on the structure without any a priori information. In this approach we describe proteins as energy-weighted graphs and compare them using ensembles of interaction networks. Investigating the position of specific interactions within the 3D native structure, we developed a parameter-free network descriptor that permits to distinguish thermostable and mesostable proteins with an accuracy of 76% and area under the receiver operating characteristic curve of 78%. Availability and implementation Code is available upon request to edoardo.milanetti@uniroma1.it Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Mattia Miotto
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, Rome, Italy.,Center for Life Nano Science@Sapienza, Instituto Italiano di Tecnologia, Viale Regina Elena, 291 Roma (RM), Italy.,Soft and Living Matter Laboratory, Institute of Nanotechnology, Consiglio Nazionale delle Ricerche, Rome, Italy
| | | | - Lorenzo Di Rienzo
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, Rome, Italy
| | - Francesco Ambrosetti
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, Rome, Italy.,Bijvoet Center for Biomolecular Research, Faculty of Science - Chemistry, Utrecht University, Padualaan 8, Utrecht, the Netherlands
| | - Pietro Corsi
- Department of Science, Università degli Studi "Roma Tre", via della Vasca Navale 84, Rome, Italy
| | - Rosalba Lepore
- Biozentrum, University of Basel, Klingelbergstrasse 50-70, CH-4056 Basel, Switzerland.,SIB Swiss Institute of Bioinformatics, Biozentrum, University of Basel, Klingelbergstrasse 50-70, CH-4056 Basel, Switzerland
| | - Gian Gaetano Tartaglia
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader St. 88, Barcelona, Spain.,Institucio' Catalana de Recerca i Estudis Avancats (ICREA), 23 Passeig Lluìs Companys, Barcelona, Spain.,Department of Biology and Biotechnology, Sapienza University of Rome, Piazzale Aldo Moro 5, Rome, Italy
| | - Edoardo Milanetti
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, Rome, Italy.,Center for Life Nano Science@Sapienza, Instituto Italiano di Tecnologia, Viale Regina Elena, 291 Roma (RM), Italy
| |
Collapse
|
14
|
López-Chávez E, Pérez-Hernández G, Aparicio F, Alas SJ. On the Thermal Stability of O 6-Methylguanine-DNA Methyltransferase from Archaeon Pyrococcus kodakaraensis by Molecular Dynamics Simulations. J Chem Inf Model 2020; 60:2138-2154. [PMID: 32250621 DOI: 10.1021/acs.jcim.0c00012] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We have employed molecular dynamics simulations to analyze the thermal stability of the O6-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.
Collapse
Affiliation(s)
- Erick López-Chávez
- Posgrado en Ciencias Naturales e Ingeniería, Unidad Cuajimalpa, Universidad Autónoma Metropolitana, Ciudad de México 05300, Mexico
| | - Gerardo Pérez-Hernández
- Departamento de Ciencias Naturales, Unidad Cuajimalpa, Universidad Autónoma Metropolitana, Ciudad de México 05300, Mexico
| | - Felipe Aparicio
- Departamento de Ciencias Naturales, Unidad Cuajimalpa, Universidad Autónoma Metropolitana, Ciudad de México 05300, Mexico
| | - Salomón J Alas
- Departamento de Ciencias Naturales, Unidad Cuajimalpa, Universidad Autónoma Metropolitana, Ciudad de México 05300, Mexico.,Departamento de Química, Unidad Iztapalapa, Universidad Autónoma Metropolitana, Ciudad de México 09340, Mexico
| |
Collapse
|
15
|
Elucidating the preference of dimeric over monomeric form for thermal stability of Thermus thermophilus isopropylmalate dehydrogenase: A molecular dynamics perspective. J Mol Graph Model 2020; 96:107530. [PMID: 31986322 DOI: 10.1016/j.jmgm.2020.107530] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2019] [Revised: 01/06/2020] [Accepted: 01/06/2020] [Indexed: 11/23/2022]
Abstract
An oligomer usually refers to a macromolecular complex formed by non-covalent interactions of monomers. Several thermophilic proteins are oligomers. The significance of oligomerization of individual proteins for stability at higher temperature is of prime importance for understanding evolution and increasing industrial productivity. The functional form of Thermus thermophilius isopropylmalate dehydrogenase (IPMDH), a widely studied protein to understand the factors affecting the thermal stability of a protein is a dimer, a simplest oligomer. To decipher the relationship between the effects of oligomerization on thermal stability of a protein, we have applied all-atom molecular mechanics approach by analyzing how temperature effects dynamics of a subunit in the presence and absence of another subunit in dimeric (SS) and monomeric forms (SA), respectively, before its denaturation begins. Comparing the difference in overall dynamic structural aspects at two different temperatures, 300 K and 337 K. Analysis of root mean square deviation (RMSD), root mean square fluctuations (RMSF) and Cα-Cα distance with an increase in temperature from 300 K to 337 K for a total of 0.2 μs reveals higher thermal stability of the dimer as compared to monomer. In contrast to dimeric form, the monomer is relatively stable at 300 K but cannot withstand the structural stability at 337 K leading to loosening of intramolecular interactions with maximum fluctuation at B23-B24 within a subunit. Energetic and structural properties indicate that B24-B24' is the major contributor to maintaining subunit-subunit interaction at 337 K. Correlation between the favorable interaction energy (IE) with the minimal perturbance in Cα atoms of domain 2 in a subunit in the presence of another subunit enhances the rigidity of the domain with subunit-subunit interaction. Overall, the study indicates that the dimeric over monomeric form enhances the protein's thermal stability and not all major subunit interacting regions contribute equally in maintaining the former.
Collapse
|
16
|
Chandrasekaran P, Santosh Kumar C, Rangachari K, Sekar K. Disassociation of β1-α1-β2 from the α2-α3 domain of prion protein (PrP) is a prerequisite for the conformational conversion of PrPC into PrPSc: Driven by the free energy landscape. Int J Biol Macromol 2019; 136:368-376. [DOI: 10.1016/j.ijbiomac.2019.06.099] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2019] [Revised: 05/29/2019] [Accepted: 06/13/2019] [Indexed: 12/16/2022]
|
17
|
Rahban M, Salehi N, Saboury AA, Hosseinkhani S, Karimi-Jafari MH, Firouzi R, Rezaei-Ghaleh N, Moosavi-Movahedi AA. Histidine substitution in the most flexible fragments of firefly luciferase modifies its thermal stability. Arch Biochem Biophys 2017; 629:8-18. [PMID: 28711358 DOI: 10.1016/j.abb.2017.07.003] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2017] [Revised: 07/08/2017] [Accepted: 07/11/2017] [Indexed: 11/26/2022]
Abstract
Molecular dynamics (MD) at two temperatures of 300 and 340 K identified two histidine residues, His461 and His489, in the most flexible regions of firefly luciferase, a light emitting enzyme. We therefore designed four protein mutants H461D, H489K, H489D and H489M to investigate their enzyme kinetic and thermodynamic stability changes. Substitution of His461 by aspartate (H461D) decreased ATP binding affinity, reduced the melting temperature of protein by around 25 °C and shifted its optimum temperature of activity to 10 °C. In line with the common feature of psychrophilic enzymes, the MD data showed that the overall flexibility of H461D was relatively high at low temperature, probably due to a decrease in the number of salt bridges around the mutation site. On the other hand, substitution of His489 by aspartate (H489D) introduced a new salt bridge between the C-terminal and N-terminal domains and increased protein rigidity but only slightly improved its thermal stability. Similar changes were observed for H489K and, to a lesser degree, H489M mutations. Based on our results we conclude that the MD simulation-based rational substitution of histidines by salt-bridge forming residues can modulate conformational dynamics in luciferase and shift its optimal temperature activity.
Collapse
Affiliation(s)
- Mahdie Rahban
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - Najmeh Salehi
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - Ali Akbar Saboury
- Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran.
| | - Saman Hosseinkhani
- Department of Biochemistry, Faculty of Biological Sciences, Tarbiat Modares University, Tehran, Iran.
| | | | - Rohoullah Firouzi
- Department of Physical Chemistry, Chemistry and Chemical Engineering Research Center of Iran, Tehran, Iran
| | | | | |
Collapse
|
18
|
Chakravorty D, Khan MF, Patra S. Multifactorial level of extremostability of proteins: can they be exploited for protein engineering? Extremophiles 2017; 21:419-444. [PMID: 28283770 DOI: 10.1007/s00792-016-0908-9] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2016] [Accepted: 12/19/2016] [Indexed: 12/20/2022]
Abstract
Research on extremostable proteins has seen immense growth in the past decade owing to their industrial importance. Basic research of attributes related to extreme-stability requires further exploration. Modern mechanistic approaches to engineer such proteins in vitro will have more impact in industrial biotechnology economy. Developing a priori knowledge about the mechanism behind extreme-stability will nurture better understanding of pathways leading to protein molecular evolution and folding. This review is a vivid compilation about all classes of extremostable proteins and the attributes that lead to myriad of adaptations divulged after an extensive study of 6495 articles belonging to extremostable proteins. Along with detailing on the rationale behind extreme-stability of proteins, emphasis has been put on modern approaches that have been utilized to render proteins extremostable by protein engineering. It was understood that each protein shows different approaches to extreme-stability governed by minute differences in their biophysical properties and the milieu in which they exist. Any general rule has not yet been drawn regarding adaptive mechanisms in extreme environments. This review was further instrumental to understand the drawback of the available 14 stabilizing mutation prediction algorithms. Thus, this review lays the foundation to further explore the biophysical pleiotropy of extreme-stable proteins to deduce a global prediction model for predicting the effect of mutations on protein stability.
Collapse
Affiliation(s)
- Debamitra Chakravorty
- Department of Biosciences and Bioengineering, Indian Institute of Technology Guwahati, Guwahati, 781039, Assam, India
| | - Mohd Faheem Khan
- Department of Biosciences and Bioengineering, Indian Institute of Technology Guwahati, Guwahati, 781039, Assam, India
| | - Sanjukta Patra
- Department of Biosciences and Bioengineering, Indian Institute of Technology Guwahati, Guwahati, 781039, Assam, India.
| |
Collapse
|
19
|
BIKKINA SWETHA, BHATI AGASTYAP, PADHI SILADITYA, PRIYAKUMAR UDEVA. Temperature Dependence of the Stability of Ion Pair Interactions, and its Implications on the Thermostability of Proteins from Thermophiles. J CHEM SCI 2017. [DOI: 10.1007/s12039-017-1231-4] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
20
|
Contribution of main chain and side chain atoms and their locations to the stability of thermophilic proteins. J Mol Graph Model 2016; 64:85-93. [DOI: 10.1016/j.jmgm.2016.01.001] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2015] [Accepted: 01/03/2016] [Indexed: 11/21/2022]
|
21
|
Sawle L, Ghosh K. Convergence of Molecular Dynamics Simulation of Protein Native States: Feasibility vs Self-Consistency Dilemma. J Chem Theory Comput 2016; 12:861-9. [DOI: 10.1021/acs.jctc.5b00999] [Citation(s) in RCA: 43] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Lucas Sawle
- Department of Physics and
Astronomy, University of Denver, Denver, Colorado 80209, United States
| | - Kingshuk Ghosh
- Department of Physics and
Astronomy, University of Denver, Denver, Colorado 80209, United States
| |
Collapse
|
22
|
Sharma R, Sastry GN. Deciphering the Dynamics of Non-Covalent Interactions Affecting Thermal Stability of a Protein: Molecular Dynamics Study on Point Mutant of Thermus thermophilus Isopropylmalate Dehydrogenase. PLoS One 2015; 10:e0144294. [PMID: 26657745 PMCID: PMC4689552 DOI: 10.1371/journal.pone.0144294] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2015] [Accepted: 11/16/2015] [Indexed: 11/18/2022] Open
Abstract
Thermus thermophilius isopropylmalate dehydrogenase catalyzes oxidative decarboxylation and dehydrogenation of isopropylmalate. Substitution of leucine to alanine at position 172 enhances the thermal stability among the known point mutants. Exploring the dynamic properties of non-covalent interactions such as saltbridges, hydrogen bonds and hydrophobic interactions to explain thermal stability of a protein is interesting in its own right. In this study dynamic changes in the non-covalent interactions are studied to decipher the deterministic features of thermal stability of a protein considering a case study of a point mutant in Thermus thermophilus isopropylmalate dehydrogenase. A total of four molecular dynamic simulations of 0.2 μs were carried out on wild type and mutant's functional dimers at 300 K and 337 K. Higher thermal stability of the mutant as compared to wild type is revealed by root mean square deviation, root mean square fluctuations and Cα-Cα distance with an increase in temperature from 300 K to 337 K. Most of the regions of wild type fluctuate higher than the corresponding regions of mutant with an increase in temperature. Cα-Cα distance analysis suggests that long distance networks are significantly affected in wild type as compared to the mutant. Short lived contacts are higher in wild type, while long lived contacts are lost at 337 K. The mutant forms less hydrogen bonds with water as compared to wild type at 337 K. In contrast to wild type, the mutant shows significant increase in unique saltbridges, hydrogen bonds and hydrophobic contacts at 337 K. The current study indicates that there is a strong inter-dependence of thermal stability on the way in which non-covalent interactions reorganize, and it is rewarding to explore this connection in single mutant studies.
Collapse
Affiliation(s)
- Reetu Sharma
- Centre for Molecular Modeling, CSIR-Indian Institute of Chemical Technology, Tarnaka, Hyderabad, 500007, India
| | - G. Narahari Sastry
- Centre for Molecular Modeling, CSIR-Indian Institute of Chemical Technology, Tarnaka, Hyderabad, 500007, India
- * E-mail:
| |
Collapse
|
23
|
Bezsudnova EY, Petrova TE, Popinako AV, Antonov MY, Stekhanova TN, Popov VO. Intramolecular hydrogen bonding in the polyextremophilic short-chain dehydrogenase from the archaeon Thermococcus sibiricus and its close structural homologs. Biochimie 2015; 118:82-9. [PMID: 26300061 DOI: 10.1016/j.biochi.2015.08.010] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2015] [Accepted: 08/19/2015] [Indexed: 01/27/2023]
Abstract
The short-chain alcohol dehydrogenase from the archaeon Thermococcus sibiricus (TsAdh319) exhibits adaptation to different kinds of stress: high temperature, high salinity, and the presence of organic solvents and denaturants. Previously a comparison of TsAdh319 with close structural homologs revealed an abnormally large number of charged residues on the surface of TsAdh319 tetramer. We further focused on the analysis of hydrogen bonding of TsAdh319 and its structural homologs from thermophilic and mesophilic organisms as a structural factor of adaptation to extreme environment. The calculation and analysis of the dynamics of hydrogen bonds of different kind were performed. In particular, the intramolecular hydrogen bonds of different kind according to their location and the type of a.a. residues involved in the bond were analyzed. TsAdh319 showed the greatest contribution of charged residues to the formation of surface hydrogen bonds, inner hydrogen bonding, and the bonds between different subunits compared to its structural homologs. Molecular dynamics simulations revealed that, of three enzyme molecules analyzed, TsAdh319 shows the least change in the number of hydrogen bonds of different kinds upon a temperature shift from 27 to 85 °C. The greatest changes were observed for a homologous enzyme from a mesophilic host. Only guanidine hydrochloride being a charged agent was able to deactivate TsAdh319. We suggest that the percentage of charged residues plays a key role in the resistance of TsAdh319 to environmental stress. The analysis shows that salt bridges in TsAdh319 serve as a universal instrument of stabilization under different extreme conditions.
Collapse
Affiliation(s)
- Ekaterina Yu Bezsudnova
- A.N. Bach Institute of Biochemistry, Research Center of Biotechnology of the Russian Academy of Sciences, Leninsky pr. 33, bld. 2, 119071, Moscow, Russian Federation.
| | - Tatiana E Petrova
- Institute of Mathematical Problems of Biology, RAS, Institutskaya str. 4, Pushchino, 142290, Russian Federation
| | - Anna V Popinako
- A.N. Bach Institute of Biochemistry, Research Center of Biotechnology of the Russian Academy of Sciences, Leninsky pr. 33, bld. 2, 119071, Moscow, Russian Federation
| | - Mikhail Yu Antonov
- M.K. Ammosov North-Eastern Federal University, Belinskiy str., 58, Suite 312, Yakutsk, 677980, Republic of Sakha (Yakutia), Russian Federation
| | - Tatiana N Stekhanova
- A.N. Bach Institute of Biochemistry, Research Center of Biotechnology of the Russian Academy of Sciences, Leninsky pr. 33, bld. 2, 119071, Moscow, Russian Federation
| | - Vladimir O Popov
- A.N. Bach Institute of Biochemistry, Research Center of Biotechnology of the Russian Academy of Sciences, Leninsky pr. 33, bld. 2, 119071, Moscow, Russian Federation; Dep. "Protein Factory", NBICS Center, National Research Centre "Kurchatov Institute", Akad. Kurchatova sqr., 1, Moscow, 123182, Russian Federation
| |
Collapse
|
24
|
Singh B, Bulusu G, Mitra A. Understanding the thermostability and activity of Bacillus subtilis lipase mutants: insights from molecular dynamics simulations. J Phys Chem B 2015; 119:392-409. [PMID: 25495458 DOI: 10.1021/jp5079554] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
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.
Collapse
Affiliation(s)
- Bipin Singh
- Center for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology Hyderabad (IIIT-H) , Gachibowli, Hyderabad, 500032, India
| | | | | |
Collapse
|
25
|
Fuciños P, Atanes E, López-López O, Solaroli M, Cerdán ME, González-Siso MI, Pastrana L, Rúa ML. Cloning, expression, purification and characterization of an oligomeric His-tagged thermophilic esterase from Thermus thermophilus HB27. Process Biochem 2014. [DOI: 10.1016/j.procbio.2014.03.006] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
|
26
|
Kalimeri M, Rahaman O, Melchionna S, Sterpone F. How conformational flexibility stabilizes the hyperthermophilic elongation factor G-domain. J Phys Chem B 2013; 117:13775-85. [PMID: 24087838 DOI: 10.1021/jp407078z] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
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.
Collapse
Affiliation(s)
- Maria Kalimeri
- Laboratoire de Biochimie Théorique, IBPC, CNRS, UPR9080, Université Paris Diderot , Sorbonne Paris Cité, France
| | | | | | | |
Collapse
|