Insights Into The Effects of Amino Acid Substitutions on The Stability of Reteplase Structure: A Molecular Dynamics Simulation Study

Overview

Journal Iran J Biotechnol

Specialty Biotechnology

Date 2023 Feb 22

PMID 36811105

Authors

Kaveh Haji-Allahverdipoor

Mokhtar Jalali Javaran

Sajad Rashidi Monfared

Mohamad Bagher Khadem-Erfan

Bahram Nikkhoo

Zhila Bahrami Rad

Habib Eslami

Sherko Nasseri

Affiliations

Soon will be listed here.

Abstract

Background: Reteplase (recombinant plasminogen activator, r-PA) is a recombinant protein designed to imitate the endogenous tissue plasminogen activator and catalyze the plasmin production. It is known that the application of reteplase is limited by the complex production processes and protein's stability challenges. Computational redesign of proteins has gained momentum in recent years, particularly as a powerful tool for improving protein stability and consequently its production efficiency. Hence, in the current study, we implemented computational approaches to improve r-PA conformational stability, which fairly correlates with protein's resistance to proteolysis.

Objectives: The current study was developed in order to evaluate the effect of amino acid substitutions on the stability of reteplase structure using molecular dynamic simulations and computational predictions.

Materials And Methods: Several web servers designed for mutation analysis were utilized to select appropriate mutations. Additionally, the experimentally reported mutation, R103S, converting wild type r-PA into non-cleavable form, was also employed. Firstly, mutant collection, consisting of 15 structures, was constructed based on the combinations of four designated mutations. Then, 3D structures were generated using MODELLER. Finally, 17 independent 20-ns molecular dynamics (MD) simulations were conducted and different analysis were performed like root-mean-square deviation (RMSD), root-mean-square fluctuations (RMSF), secondary structure analysis, number of hydrogen bonds, principal components analysis (PCA), eigenvector projection, and density analysis.

Results: Predicted mutations successfully compensated the more flexible conformation caused by R103S substitution, so, improved conformational stability was analyzed from MD simulations. In particular, R103S/A286I/G322I indicated the best results and remarkably enhanced the protein stability.

Conclusion: The conformational stability conferred by these mutations will probably lead to more protection of r-PA in protease-rich environments in various recombinant systems and potentially enhance its production and expression level.

Citing Articles

Direct Interaction of Minocycline to p47phox Contributes to its Attenuation of TNF-α-Mediated Neuronal PC12 Cell Death: Experimental and Simulation Validation.

Eslami H, Rokhzadi K, Basiri M, Esmaeili-Mahani S, Mahmoodi Z, Haji-Allahverdipoor K Cell Biochem Biophys. 2024; 82(2):1261-1277.

PMID: 38739323 DOI: 10.1007/s12013-024-01279-9.

References

Hidalgo D, Abdoli-Nasab M, Jalali-Javaran M, Bru-Martinez R, Cusido R, Corchete P . Biotechnological production of recombinant tissue plasminogen activator protein (reteplase) from transplastomic tobacco cell cultures. Plant Physiol Biochem. 2017; 118:130-137. DOI: 10.1016/j.plaphy.2017.06.013. View

Cai Y, Bao J, Lao X, Zheng H, Chen J, Yu R . Computational design, functional analysis and antigenic epitope estimation of a novel hybrid of 12 peptides of hirudin and reteplase. J Mol Model. 2015; 21(9):229. DOI: 10.1007/s00894-015-2774-2. View

Goldenzweig A, Fleishman S . Principles of Protein Stability and Their Application in Computational Design. Annu Rev Biochem. 2018; 87:105-129. DOI: 10.1146/annurev-biochem-062917-012102. View

Wijesinha-Bettoni R, Alexeev Y, Johnson P, Marsh J, Sancho A, Abdullah S . The structural characteristics of nonspecific lipid transfer proteins explain their resistance to gastroduodenal proteolysis. Biochemistry. 2010; 49(10):2130-9. DOI: 10.1021/bi901939z. View

Monzon A, Zea D, Marino-Buslje C, Parisi G . Homology modeling in a dynamical world. Protein Sci. 2017; 26(11):2195-2206. PMC: 5654859. DOI: 10.1002/pro.3274. View

Dyson H, Wright P, Scheraga H . The role of hydrophobic interactions in initiation and propagation of protein folding. Proc Natl Acad Sci U S A. 2006; 103(35):13057-61. PMC: 1559752. DOI: 10.1073/pnas.0605504103. View

Gong L, Liu M, Zeng T, Shi X, Yuan C, Andreasen P . Crystal Structure of the Michaelis Complex between Tissue-type Plasminogen Activator and Plasminogen Activators Inhibitor-1. J Biol Chem. 2015; 290(43):25795-804. PMC: 4646234. DOI: 10.1074/jbc.M115.677567. View

Abdoli Nasab M, Javaran M, Cusido R, Palazon J . Purification of recombinant tissue plasminogen activator (rtPA) protein from transplastomic tobacco plants. Plant Physiol Biochem. 2016; 108:139-144. DOI: 10.1016/j.plaphy.2016.06.029. View

Magliery T . Protein stability: computation, sequence statistics, and new experimental methods. Curr Opin Struct Biol. 2015; 33:161-8. PMC: 4858314. DOI: 10.1016/j.sbi.2015.09.002. View

10.

Pires D, Ascher D, Blundell T . mCSM: predicting the effects of mutations in proteins using graph-based signatures. Bioinformatics. 2013; 30(3):335-42. PMC: 3904523. DOI: 10.1093/bioinformatics/btt691. View

11.

Lindorff-Larsen K, Piana S, Palmo K, Maragakis P, Klepeis J, Dror R . Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins. 2010; 78(8):1950-8. PMC: 2970904. DOI: 10.1002/prot.22711. View

12.

van der Spoel D, Lindahl E, Hess B, Groenhof G, Mark A, Berendsen H . GROMACS: fast, flexible, and free. J Comput Chem. 2005; 26(16):1701-18. DOI: 10.1002/jcc.20291. View

13.

Schiermeyer A, Schinkel H, Apel S, Fischer R, Schillberg S . Production of Desmodus rotundus salivary plasminogen activator alpha1 (DSPAalpha1) in tobacco is hampered by proteolysis. Biotechnol Bioeng. 2005; 89(7):848-58. DOI: 10.1002/bit.20410. View

14.

Kulshreshtha S, Chaudhary V, Goswami G, Mathur N . Computational approaches for predicting mutant protein stability. J Comput Aided Mol Des. 2016; 30(5):401-12. DOI: 10.1007/s10822-016-9914-3. View

15.

Pace C, Scholtz J, Grimsley G . Forces stabilizing proteins. FEBS Lett. 2014; 588(14):2177-84. PMC: 4116631. DOI: 10.1016/j.febslet.2014.05.006. View

16.

Webb B, Sali A . Comparative Protein Structure Modeling Using MODELLER. Curr Protoc Bioinformatics. 2016; 54:5.6.1-5.6.37. PMC: 5031415. DOI: 10.1002/cpbi.3. View

17.

Kalhor H, Sadeghi S, Marashiyan M, Enssi M, Kalhor R, Ganji M . In silico mutagenesis in recombinant human keratinocyte growth factor: Improvement of stability and activity in addition to decrement immunogenicity. J Mol Graph Model. 2020; 97:107551. DOI: 10.1016/j.jmgm.2020.107551. View

18.

Goulay R, Naveau M, Gaberel T, Vivien D, Parcq J . Optimized tPA: A non-neurotoxic fibrinolytic agent for the drainage of intracerebral hemorrhages. J Cereb Blood Flow Metab. 2017; 38(7):1180-1189. PMC: 6434446. DOI: 10.1177/0271678X17719180. View

19.

Sadeghi S, Poorebrahim M, Rahimi H, Karimipoor M, Azadmanesh K, Khorramizadeh M . In silico studying of the whole protein structure and dynamics of Dickkopf family members showed that N-terminal domain of Dickkopf 2 in contrary to other Dickkopfs facilitates its interaction with low density lipoprotein receptor related protein 5/6. J Biomol Struct Dyn. 2018; 37(10):2564-2580. DOI: 10.1080/07391102.2018.1491891. View

20.

Worth C, Preissner R, Blundell T . SDM--a server for predicting effects of mutations on protein stability and malfunction. Nucleic Acids Res. 2011; 39(Web Server issue):W215-22. PMC: 3125769. DOI: 10.1093/nar/gkr363. View