1
|
Bao H, Wang W, Sun H, Chen J. The switch states of the GDP-bound HRAS affected by point mutations: a study from Gaussian accelerated molecular dynamics simulations and free energy landscapes. J Biomol Struct Dyn 2024; 42:3363-3381. [PMID: 37216340 DOI: 10.1080/07391102.2023.2213355] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2023] [Accepted: 05/04/2023] [Indexed: 05/24/2023]
Abstract
Point mutations play a vital role in the conformational transformation of HRAS. In this work, Gaussian accelerated molecular dynamics (GaMD) simulations followed by constructions of free energy landscapes (FELs) were adopted to explore the effect of mutations D33K, A59T and L120A on conformation states of the GDP-bound HRAS. The results from the post-processing analyses on GaMD trajectories suggest that mutations alter the flexibility and motion modes of the switch domains from HRAS. The analyses from FELs show that mutations induce more disordered states of the switch domains and affect interactions of GDP with HRAS, implying that mutations yield a vital effect on the binding of HRAS to effectors. The GDP-residue interaction network revealed by our current work indicates that salt bridges and hydrogen bonding interactions (HBIs) play key roles in the binding of GDP to HRAS. Furthermore, instability in the interactions of magnesium ions and GDP with the switch SI leads to the extreme disorder of the switch domains. This study is expected to provide the energetic basis and molecular mechanism for further understanding the function of HRAS.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Huayin Bao
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
2
|
Yu Z, Wang Z, Cui X, Cao Z, Zhang W, Sun K, Hu G. Conformational States of the GDP- and GTP-Bound HRAS Affected by A59E and K117R: An Exploration from Gaussian Accelerated Molecular Dynamics. Molecules 2024; 29:645. [PMID: 38338389 PMCID: PMC10856033 DOI: 10.3390/molecules29030645] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2023] [Revised: 01/01/2024] [Accepted: 01/25/2024] [Indexed: 02/12/2024] Open
Abstract
The HRAS protein is considered a critical target for drug development in cancers. It is vital for effective drug development to understand the effects of mutations on the binding of GTP and GDP to HRAS. We conducted Gaussian accelerated molecular dynamics (GaMD) simulations and free energy landscape (FEL) calculations to investigate the impacts of two mutations (A59E and K117R) on GTP and GDP binding and the conformational states of the switch domain. Our findings demonstrate that these mutations not only modify the flexibility of the switch domains, but also affect the correlated motions of these domains. Furthermore, the mutations significantly disrupt the dynamic behavior of the switch domains, leading to a conformational change in HRAS. Additionally, these mutations significantly impact the switch domain's interactions, including their hydrogen bonding with ligands and electrostatic interactions with magnesium ions. Since the switch domains are crucial for the binding of HRAS to effectors, any alterations in their interactions or conformational states will undoubtedly disrupt the activity of HRAS. This research provides valuable information for the design of drugs targeting HRAS.
Collapse
Affiliation(s)
- Zhiping Yu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China; (Z.Y.); (Z.C.)
| | - Zhen Wang
- Pingyin People’s Hospital, Jinan 250400, China; (Z.W.); (X.C.)
| | - Xiuzhen Cui
- Pingyin People’s Hospital, Jinan 250400, China; (Z.W.); (X.C.)
| | - Zanxia Cao
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China; (Z.Y.); (Z.C.)
| | - Wanyunfei Zhang
- School of Science, Xi’an Polytechnic University, Xi’an 710048, China; (W.Z.); (K.S.)
| | - Kunxiao Sun
- School of Science, Xi’an Polytechnic University, Xi’an 710048, China; (W.Z.); (K.S.)
| | - Guodong Hu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China; (Z.Y.); (Z.C.)
| |
Collapse
|
3
|
Wang Y, Yang F, Yan D, Zeng Y, Wei B, Chen J, He W. Identification Mechanism of BACE1 on Inhibitors Probed by Using Multiple Separate Molecular Dynamics Simulations and Comparative Calculations of Binding Free Energies. Molecules 2023; 28:4773. [PMID: 37375328 DOI: 10.3390/molecules28124773] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2023] [Revised: 06/12/2023] [Accepted: 06/12/2023] [Indexed: 06/29/2023] Open
Abstract
β-amyloid cleaving enzyme 1 (BACE1) is regarded as an important target of drug design toward the treatment of Alzheimer's disease (AD). In this study, three separate molecular dynamics (MD) simulations and calculations of binding free energies were carried out to comparatively determine the identification mechanism of BACE1 for three inhibitors, 60W, 954 and 60X. The analyses of MD trajectories indicated that the presence of three inhibitors influences the structural stability, flexibility and internal dynamics of BACE1. Binding free energies calculated by using solvated interaction energy (SIE) and molecular mechanics generalized Born surface area (MM-GBSA) methods reveal that the hydrophobic interactions provide decisive forces for inhibitor-BACE1 binding. The calculations of residue-based free energy decomposition suggest that the sidechains of residues L91, D93, S96, V130, Q134, W137, F169 and I179 play key roles in inhibitor-BACE1 binding, which provides a direction for future drug design toward the treatment of AD.
Collapse
Affiliation(s)
- Yiwen Wang
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Aeronautics, Shandong Jiaotong University, Jinan 250357, China
| | - Fen Yang
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
| | - Dongliang Yan
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Yalin Zeng
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
| | - Benzheng Wei
- Center for Medical Artificial Intelligence, Shandong University of Traditional Chinese Medicine, Qingdao 266112, China
| | - Jianzhong Chen
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Weikai He
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Aeronautics, Shandong Jiaotong University, Jinan 250357, China
| |
Collapse
|
4
|
Bassani D, Moro S. Past, Present, and Future Perspectives on Computer-Aided Drug Design Methodologies. Molecules 2023; 28:molecules28093906. [PMID: 37175316 PMCID: PMC10180087 DOI: 10.3390/molecules28093906] [Citation(s) in RCA: 19] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2023] [Revised: 04/28/2023] [Accepted: 05/02/2023] [Indexed: 05/15/2023] Open
Abstract
The application of computational approaches in drug discovery has been consolidated in the last decades. These families of techniques are usually grouped under the common name of "computer-aided drug design" (CADD), and they now constitute one of the pillars in the pharmaceutical discovery pipelines in many academic and industrial environments. Their implementation has been demonstrated to tremendously improve the speed of the early discovery steps, allowing for the proficient and rational choice of proper compounds for a desired therapeutic need among the extreme vastness of the drug-like chemical space. Moreover, the application of CADD approaches allows the rationalization of biochemical and interactive processes of pharmaceutical interest at the molecular level. Because of this, computational tools are now extensively used also in the field of rational 3D design and optimization of chemical entities starting from the structural information of the targets, which can be experimentally resolved or can also be obtained with other computer-based techniques. In this work, we revised the state-of-the-art computer-aided drug design methods, focusing on their application in different scenarios of pharmaceutical and biological interest, not only highlighting their great potential and their benefits, but also discussing their actual limitations and eventual weaknesses. This work can be considered a brief overview of computational methods for drug discovery.
Collapse
Affiliation(s)
- Davide Bassani
- Pharmaceutical Research & Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., 4070 Basel, Switzerland
- Molecular Modeling Section (MMS), Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Via Marzolo 5, 35131 Padova, Italy
| | - Stefano Moro
- Molecular Modeling Section (MMS), Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Via Marzolo 5, 35131 Padova, Italy
| |
Collapse
|
5
|
Pan Y, Xie N, Zhang X, Yang S, Lv S. Computational Insights into the Dynamic Structural Features and Binding Characteristics of Recombinase UvsX Compared with RecA. Molecules 2023; 28:molecules28083363. [PMID: 37110596 PMCID: PMC10144138 DOI: 10.3390/molecules28083363] [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: 03/03/2023] [Revised: 04/06/2023] [Accepted: 04/07/2023] [Indexed: 04/29/2023] Open
Abstract
RecA family recombinases are the core enzymes in the process of homologous recombination, and their normal operation ensures the stability of the genome and the healthy development of organisms. The UvsX protein from bacteriophage T4 is a member of the RecA family recombinases and plays a central role in T4 phage DNA repair and replication, which provides an important model for the biochemistry and genetics of DNA metabolism. UvsX shares a high degree of structural similarity and function with RecA, which is the most deeply studied member of the RecA family. However, the detailed molecular mechanism of UvsX has not been resolved. In this study, a comprehensive all-atom molecular dynamics simulation of the UvsX protein dimer complex was carried out in order to investigate the conformational and binding properties of UvsX in combination with ATP and DNA, and the simulation of RecA was synchronized with the property comparison learning for UvsX. This study confirmed the highly conserved molecular structure characteristics and catalytic centers of RecA and UvsX, and also discovered differences in regional conformation, volatility and the ability to bind DNA between the two proteins at different temperatures, which would be helpful for the subsequent understanding and application of related recombinases.
Collapse
Affiliation(s)
- Yue Pan
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Ningkang Xie
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Xin Zhang
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Shuo Yang
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Shaowu Lv
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China
- Bioarchaeology Laboratory, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| |
Collapse
|
6
|
Shi S, Zheng L, Ren Y, Wang Z. Impacts of Mutations in the P-Loop on Conformational Alterations of KRAS Investigated with Gaussian Accelerated Molecular Dynamics Simulations. Molecules 2023; 28:molecules28072886. [PMID: 37049650 PMCID: PMC10095679 DOI: 10.3390/molecules28072886] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2023] [Revised: 03/19/2023] [Accepted: 03/20/2023] [Indexed: 04/14/2023] Open
Abstract
G12 mutations heavily affect conformational transformation and activity of KRAS. In this study, Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the GDP-bound wild-type (WT), G12A, G12D, and G12R KRAS to probe mutation-mediated impacts on conformational alterations of KRAS. The results indicate that three G12 mutations obviously affect the structural flexibility and internal dynamics of the switch domains. The analyses of the free energy landscapes (FELs) suggest that three G12 mutations induce more conformational states of KRAS and lead to more disordered switch domains. The principal component analysis shows that three G12 mutations change concerted motions and dynamics behavior of the switch domains. The switch domains mostly overlap with the binding region of KRAS to its effectors. Thus, the high disorder states and concerted motion changes of the switch domains induced by G12 mutations affect the activity of KRAS. The analysis of interaction network of GDP with KRAS signifies that the instability in the interactions of GDP and magnesium ion with the switch domain SW1 drives the high disordered state of the switch domains. This work is expected to provide theoretical aids for understanding the function of KRAS.
Collapse
Affiliation(s)
- Shuhua Shi
- School of Science, Shandong Jianzhu University, Jinan 250101, China
| | - Linqi Zheng
- School of Science, Shandong Jianzhu University, Jinan 250101, China
| | - Yonglian Ren
- School of Science, Shandong Jianzhu University, Jinan 250101, China
| | - Ziyu Wang
- School of Science, Shandong Jianzhu University, Jinan 250101, China
- Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
| |
Collapse
|
7
|
Wang L, Wang Y, Yu Y, Liu D, Zhao J, Zhang L. Deciphering Selectivity Mechanism of BRD9 and TAF1(2) toward Inhibitors Based on Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations. Molecules 2023; 28:molecules28062583. [PMID: 36985555 PMCID: PMC10052767 DOI: 10.3390/molecules28062583] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Revised: 03/06/2023] [Accepted: 03/09/2023] [Indexed: 03/14/2023] Open
Abstract
BRD9 and TAF1(2) have been regarded as significant targets of drug design for clinically treating acute myeloid leukemia, malignancies, and inflammatory diseases. In this study, multiple short molecular dynamics simulations combined with the molecular mechanics generalized Born surface area method were employed to investigate the binding selectivity of three ligands, 67B, 67C, and 69G, to BRD9/TAF1(2) with IC50 values of 230/59 nM, 1400/46 nM, and 160/410 nM, respectively. The computed binding free energies from the MM-GBSA method displayed good correlations with that provided by the experimental data. The results indicate that the enthalpic contributions played a critical factor in the selectivity recognition of inhibitors toward BRD9 and TAF1(2), indicating that 67B and 67C could more favorably bind to TAF1(2) than BRD9, while 69G had better selectivity toward BRD9 over TAF1(2). In addition, the residue-based free energy decomposition approach was adopted to calculate the inhibitor–residue interaction spectrum, and the results determined the gatekeeper (Y106 in BRD9 and Y1589 in TAF1(2)) and lipophilic shelf (G43, F44, and F45 in BRD9 and W1526, P1527, and F1528 in TAF1(2)), which could be identified as hotspots for designing efficient selective inhibitors toward BRD9 and TAF1(2). This work is also expected to provide significant theoretical guidance and insightful molecular mechanisms for the rational designs of efficient selective inhibitors targeting BRD9 and TAF1(2).
Collapse
|
8
|
Lu J, Wang C, Ma Y, Liu K, Fu X, Xing S. Exploration of the Product Specificity of chitosanase CsnMY002 and Mutants Using Molecular Dynamics Simulations. Molecules 2023; 28:molecules28031048. [PMID: 36770713 PMCID: PMC9920700 DOI: 10.3390/molecules28031048] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2022] [Revised: 01/06/2023] [Accepted: 01/17/2023] [Indexed: 01/22/2023] Open
Abstract
Chitosanase CsnMY002 is a new type of enzyme isolated from Bacillus subtilis that is used to prepare chitosan oligosaccharide. Although mutants G21R and G21K could increase Chitosan yield and thus increase the commercial value of the final product, the mechanism by which this happens is not known. Herein, we used molecular dynamics simulations to explore the conformational changes in CsnMY002 wild type and mutants when they bind substrates. The binding of substrate changed the conformation of protein, stretching and deforming the active and catalytic region. Additionally, the mutants caused different binding modes and catalysis, resulting in different degrees of polymerization of the final Chitooligosaccharide degradation product. Finally, Arg37, Ile145 ~ Gly148 and Trp204 are important catalytic residues of CsnMY002. Our study provides a basis for the engineering of chitosanases.
Collapse
Affiliation(s)
- Jianzhang Lu
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China
- University of Edinburgh Institute (ZJU-UoE Institute), Zhejiang University, Haining 314400, China
| | - Chu Wang
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China
| | - Yingying Ma
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130012, China
| | - Kaifeng Liu
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130012, China
| | - Xueqi Fu
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China
| | - Shu Xing
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China
- Correspondence:
| |
Collapse
|
9
|
Wang L, Lu D, Wang Y, Xu X, Zhong P, Yang Z. Binding selectivity-dependent molecular mechanism of inhibitors towards CDK2 and CDK6 investigated by multiple short molecular dynamics and free energy landscapes. J Enzyme Inhib Med Chem 2023; 38:84-99. [DOI: 10.1080/14756366.2022.2135511] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Affiliation(s)
- Lifei Wang
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Dan Lu
- Department of Physics, Jiangxi Agricultural University, Nanchang, PR China
| | - Yan Wang
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Xiaoyan Xu
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Peihua Zhong
- College of Computer Information and Engineering, Jiangxi Agriculture University, Nanchang, PR China
| | - Zhiyong Yang
- Department of Physics, Jiangxi Agricultural University, Nanchang, PR China
| |
Collapse
|
10
|
Bao HY, Wang W, Sun HB, Chen JZ. Binding modes of GDP, GTP and GNP to NRAS deciphered by using Gaussian accelerated molecular dynamics simulations. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2023; 34:65-89. [PMID: 36762439 DOI: 10.1080/1062936x.2023.2165542] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/23/2022] [Accepted: 12/31/2022] [Indexed: 06/18/2023]
Abstract
Probing binding modes of GDP, GTP and GNP to NRAS are of significance for understanding the regulation mechanism on the activity of RAS proteins. Four separate Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the apo, GDP-, GTP- and GNP-bound NRAS. Dynamics analyses suggest that binding of three ligands highly affects conformational states of the switch domains from NRAS, which disturbs binding of NRAS to its effectors. The analyses of free energy landscapes (FELs) indicate that binding of GDP, GTP and GNP induces more energetic states of NRAS compared to the apo NRAS but the presence of GNP makes the switch domains more ordered than binding of GDP and GNP. The information of interaction networks of ligands with NRAS reveals that the π-π interaction of residue F28 and the salt bridge interactions of K16 and D119 with ligands stabilize binding of GDP, GTP and GNP to NRAS. Meanwhile magnesium ion plays a bridge role in interactions of ligands with NRAS, which is favourable for associations of GDP, GTP and GNP with NRAS. This work is expected to provide useful information for deeply understanding the function and activity of NRAS.
Collapse
Affiliation(s)
- H Y Bao
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - J Z Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
11
|
Chen J, Zeng Q, Wang W, Sun H, Hu G. Decoding the Identification Mechanism of an SAM-III Riboswitch on Ligands through Multiple Independent Gaussian-Accelerated Molecular Dynamics Simulations. J Chem Inf Model 2022; 62:6118-6132. [PMID: 36440874 DOI: 10.1021/acs.jcim.2c00961] [Citation(s) in RCA: 49] [Impact Index Per Article: 24.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
S-Adenosyl-l-methionine (SAM)-responsive riboswitches play a central role in the regulation of bacterial gene expression at the level of transcription attenuation or translation inhibition. In this study, multiple independent Gaussian-accelerated molecular dynamics simulations were performed to decipher the identification mechanisms of SAM-III (SMK) on ligands SAM, SAH, and EEM. The results reveal that ligand binding highly affects the structural flexibility, internal dynamics, and conformational changes of SAM-III. The dynamic analysis shows that helices P3 and P4 as well as two junctions J23 and J24 of SAM-III are highly susceptible to ligand binding. Analyses of free energy landscapes suggest that ligand binding induces different free energy profiles of SAM-III, which leads to the difference in identification sites of SAM-III on ligands. The information on ligand-nucleotide interactions not only uncovers that the π-π, cation-π, and hydrogen bonding interactions drive identification of SAM-III on the three ligands but also reveals that different electrostatic properties of SAM, SAH, and EEM alter the active sites of SAM-III. Meanwhile, the results also verify that the adenine group of SAM, SAH, and EEM is well recognized by conserved nucleotides G7, A29, U37, A38, and G48. We expect that this study can provide useful information for understanding the applications of SAM-III in chemical, synthetic RNA biology, and biomedical fields.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan250357, China
| | - Qingkai Zeng
- School of Science, Shandong Jiaotong University, Jinan250357, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan250357, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan250357, China
| | - Guodong Hu
- Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou253023, China
| |
Collapse
|