1
|
Tahir Khan M, Dumont E, Chaudhry AR, Wei DQ. Free energy landscape and thermodynamics properties of novel mutations in PncA of pyrazinamide resistance isolates of Mycobacterium tuberculosis. J Biomol Struct Dyn 2023:1-12. [PMID: 37837425 DOI: 10.1080/07391102.2023.2268216] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2023] [Accepted: 09/29/2023] [Indexed: 10/16/2023]
Abstract
Pyrazinamide (PZA) is one of the first-line antituberculosis therapy, active against non-replicating Mycobacterium tuberculosis (Mtb). The conversion of PZA into pyrazinoic acid (POA), the active form, required the activity of pncA gene product pyrazinamidase (PZase) activity. Mutations occurred in pncA are the primary cause behind the PZA resistance. However, the resistance mechanism is important to explore using high throughput computational approaches. Here we aimed to explore the mechanism of PZA resistance behind novel P62T, L120R, and V130M mutations in PZase using 200 ns molecular dynamics (MD) simulations. MD simulations were performed to observe the structural changes for these three mutants (MTs) compared to the wild types (WT). Root means square fluctuation, the radius of gyration, free energy landscape, root means square deviation, dynamic cross-correlation motion, and pocket volume were found in variation between WT and MTs, revealing the effects of P62T, L120R, and V130M. The free energy conformational landscape of MTs differs significantly from the WT system, lowering the binding of PZA. The geometric shape complementarity of the drug (PZA) and target protein (PZase) further confirmed that P62T, L120R, and V130M affect the protein structure. These effects on PZase may cause vulnerability to convert PZA into POA.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Muhammad Tahir Khan
- Zhongjing Research and Industrialization Institute of Chinese Medicine, Zhongguancun Scientific Park, Nanyang, PR China
- Institute of Molecular Biology and Biotechnology (IMBB), The University of Lahore, Lahore, Pakistan
| | - Elise Dumont
- Université Côte d'Azur, CNRS, Institut de Chimie de Nice, UMR7272, Nice, France
- Institut Universitaire de France, Paris, France
| | | | | |
Collapse
|
2
|
Xie Q, Kasahara K, Higo J, Takahashi T. Molecular Mechanisms of Functional Modulation of Transcriptional Coactivator PC4 via Phosphorylation on Its Intrinsically Disordered Region. ACS OMEGA 2023; 8:14572-14582. [PMID: 37125110 PMCID: PMC10134458 DOI: 10.1021/acsomega.3c00364] [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: 01/17/2023] [Accepted: 04/03/2023] [Indexed: 05/03/2023]
Abstract
To investigate the effects of phosphorylation on the function of the human positive cofactor 4 (PC4), an enhanced molecular dynamics (MD) simulation was performed. The simulation system consists of the N-terminal intrinsic disordered region (IDR) of PC4 and a complex that comprises the C-terminal acidic activation domain of a herpes simplex virion protein 16 (VP16ad) and a homodimer of the C-terminal structured core domain of PC4 (PC4ctd). An earlier report of an experimental study reported that the PC4-VP16ad interaction is modulated by incremental phosphorylation of the IDR. That report also proposed a dynamic model where phosphorylated serine residues of a segment (SEAC) in the IDR contact positively charged residues (lysin and arginine) of another segment (K-rich) in the IDR. This contact formation induced by the phosphorylation results in variation of PC4-VP16ad interaction. However, this contact formation has not yet been measured directly because it is transiently formed and because the SEAC and K-rich segments are unstructured with high flexibility. We performed two simulations to mimic the incremental phosphorylation: the IDR was not phosphorylated in one simulation and only partially phosphorylated in the other. Our simulation results indicate that the phosphorylation weakens the IDR-VP16ad contact considerably with the induction of a compact structure in the IDR. This structure was stabilized by electrostatic interactions between the phosphorylated serine residues of a segment and lysine or arginine residues of another segment in the IDR, but the conformational fluctuation of this compact structure was considerably large. Consequently, the present study supports the experimentally proposed dynamic model. Results of this study can be important for computational elucidation of the functional modulation of PC4.
Collapse
Affiliation(s)
- Qilin Xie
- Graduate
School of Life Sciences, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
| | - Kota Kasahara
- College
of Life Sciences, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
| | - Junichi Higo
- Graduate
School of Information Science, University
of Hyogo, 7-1-28 minatojima
Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Research
Organization of Science and Technology, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
| | - Takuya Takahashi
- College
of Life Sciences, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
| |
Collapse
|
3
|
Fukunishi Y, Higo J, Kasahara K. Computer simulation of molecular recognition in biomolecular system: from in silico screening to generalized ensembles. Biophys Rev 2022; 14:1423-1447. [PMID: 36465086 PMCID: PMC9703445 DOI: 10.1007/s12551-022-01015-8] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Accepted: 11/06/2022] [Indexed: 11/29/2022] Open
Abstract
Prediction of ligand-receptor complex structure is important in both the basic science and the industry such as drug discovery. We report various computation molecular docking methods: fundamental in silico (virtual) screening, ensemble docking, enhanced sampling (generalized ensemble) methods, and other methods to improve the accuracy of the complex structure. We explain not only the merits of these methods but also their limits of application and discuss some interaction terms which are not considered in the in silico methods. In silico screening and ensemble docking are useful when one focuses on obtaining the native complex structure (the most thermodynamically stable complex). Generalized ensemble method provides a free-energy landscape, which shows the distribution of the most stable complex structure and semi-stable ones in a conformational space. Also, barriers separating those stable structures are identified. A researcher should select one of the methods according to the research aim and depending on complexity of the molecular system to be studied.
Collapse
Affiliation(s)
- Yoshifumi Fukunishi
- Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 2-3-26, Aomi, Koto-Ku, Tokyo, 135-0064 Japan
| | - Junichi Higo
- Graduate School of Information Science, University of Hyogo, 7-1-28 Minatojima Minamimachi, Chuo-Ku, Kobe, Hyogo 650-0047 Japan ,Research Organization of Science and Technology, Ritsumeikan University, 1-1-1 Noji-Higashi, Kusatsu, Shiga 525-8577 Japan
| | - Kota Kasahara
- College of Life Sciences, Ritsumeikan University, 1-1-1 Noji-Higashi, Kusatsu, Shiga 525-8577 Japan
| |
Collapse
|
4
|
Fly casting with ligand sliding and orientational selection supporting complex formation of a GPCR and a middle sized flexible molecule. Sci Rep 2022; 12:13792. [PMID: 35963875 PMCID: PMC9376114 DOI: 10.1038/s41598-022-17920-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2022] [Accepted: 08/02/2022] [Indexed: 11/08/2022] Open
Abstract
A GA-guided multidimensional virtual-system coupled molecular dynamics (GA-mD-VcMD) simulation was conducted to elucidate binding mechanisms of a middle-sized flexible molecule, bosentan, to a GPCR protein, human endothelin receptor type B (hETB). GA-mD-VcMD is a generalized ensemble method that produces a free-energy landscape of the ligand-receptor binding by searching large-scale motions accompanied with stable maintenance of the fragile cell-membrane structure. All molecular components (bosentan, hETB, membrane, and solvent) were represented with an all-atom model. Then sampling was conducted from conformations where bosentan was distant from the binding site in the hETB binding pocket. The deepest basin in the resultant free-energy landscape was assigned to native-like complex conformation. The following binding mechanism was inferred. First, bosentan fluctuating randomly in solution is captured using a tip region of the flexible N-terminal tail of hETB via nonspecific attractive interactions (fly casting). Bosentan then slides occasionally from the tip to the root of the N-terminal tail (ligand–sliding). During this sliding, bosentan passes the gate of the binding pocket from outside to inside of the pocket with an accompanying rapid reduction of the molecular orientational variety of bosentan (orientational selection). Last, in the pocket, ligand–receptor attractive native contacts are formed. Eventually, the native-like complex is completed. The bosentan-captured conformations by the tip-region and root-region of the N-terminal tail correspond to two basins in the free-energy landscape. The ligand-sliding corresponds to overcoming of a free-energy barrier between the basins.
Collapse
|
5
|
Higo J, Takashima H, Fukunishi Y, Yoshimori A. Generalized-ensemble method study: A helix-mimetic compound inhibits protein-protein interaction by long-range and short-range intermolecular interactions. J Comput Chem 2021; 42:956-969. [PMID: 33755222 DOI: 10.1002/jcc.26516] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2020] [Revised: 02/25/2021] [Accepted: 03/01/2021] [Indexed: 12/30/2022]
Abstract
A heterocyclic compound mS-11 is a helix-mimetic designed to inhibit binding of an intrinsic disordered protein neural restrictive silence factor/repressor element 1 silencing factor (NRSF/REST) to a receptor protein mSin3B. We apply a generalized ensemble method, multi-dimensional virtual-system coupled molecular dynamics developed by ourselves recently, to a system consisting of mS-11 and mSin3B, and obtain a thermally equilibrated distribution, which is comprised of the bound and unbound states extensively. The lowest free-energy position of mS-11 coincides with the NRSF/REST position in the experimentally-determined NRSF/REST-mSin3B complex. Importantly, the molecular orientation of mS-11 is ordering in a wide region around mSin3B. The resultant binding scenario is: When mS-11 is distant from the binding site of mSin3B, mS-11 descends the free-energy slope toward the binding site maintaining the molecular orientation to be advantageous for binding. Then, finally a long and flexible hydrophobic sidechain of mS-11 fits into the binding site, which is the lowest-free-energy complex structure inhibiting NRSF/REST binding to mSin3B.
Collapse
Affiliation(s)
- Junichi Higo
- Graduate School of Simulation Studies, University of Hyogo, Kobe, Japan
| | - Hajime Takashima
- Department of Research and Development, PRISM BioLab Co., Ltd., Fujisawa, Japan
| | - Yoshifumi Fukunishi
- Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Tokyo, Japan
| | - Atsushi Yoshimori
- Chemoinformatics & AI Research Group, Institute for Theoretical Medicine, Inc., Fujisawa, Japan
| |
Collapse
|
6
|
Difference of binding modes among three ligands to a receptor mSin3B corresponding to their inhibitory activities. Sci Rep 2021; 11:6178. [PMID: 33731831 PMCID: PMC7971087 DOI: 10.1038/s41598-021-85612-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2021] [Accepted: 02/19/2021] [Indexed: 11/08/2022] Open
Abstract
A preceding experiment suggested that a compound, which inhibits binding of the REST/NRSF segment to the cleft of a receptor protein mSin3B, can be a potential drug candidate to ameliorate many neuropathies. We have recently developed an enhanced conformational sampling method, genetic-algorithm-guided multi-dimensional virtual-system-coupled canonical molecular dynamics, and in the present study, applied it to three systems consisting of mSin3B and one of three compounds, sertraline, YN3, and acitretin. Other preceding experiments showed that only sertraline inhibits the binding of REST/NRSF to mSin3B. The current simulation study produced the spatial distribution of the compounds around mSin3B, and showed that sertraline and YN3 bound to the cleft of mSin3B with a high propensity, although acitretin did not. Further analyses of the simulation data indicated that only the sertraline-mSin3B complex produced a hydrophobic core similar to that observed in the molecular interface of the REST/NRSF-mSin3B complex: An aromatic ring of sertraline sunk deeply in the mSin3B's cleft forming a hydrophobic core contacting to hydrophobic amino-acid residues located at the bottom of the cleft. The present study proposes a step to design a compound that inhibits competitively the binding of a ligand to its receptor.
Collapse
|
7
|
Higo J, Kusaka A, Kasahara K, Kamiya N, Hayato I, Qilin X, Takahashi T, Fukuda I, Mori K, Hata Y, Fukunishi Y. GA-guided mD-VcMD: A genetic-algorithm-guided method for multi-dimensional virtual-system coupled molecular dynamics. Biophys Physicobiol 2021; 17:161-176. [PMID: 33585149 DOI: 10.2142/biophysico.bsj-2020008] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2020] [Accepted: 12/22/2020] [Indexed: 12/01/2022] Open
Abstract
We introduced a conformational sampling method in an earlier report: The multi-dimensional virtual-system coupled molecular dynamics (mD-VcMD) enhances conformational sampling of a biomolecular system by computer simulations. Herein, new sampling method, a subzone-based mD-VcMD, is presented as an extension of mD-VcMD. Then, the subzone-based method is extended further using a genetic algorithm (GA) named the GA-guided mD-VcMD. In these methods, iterative simulation runs are performed to increase the sampled region gradually. The new methods have the following benefits: (1) They are free from a production run: i.e., all snapshots from all iterations are useful for analyses. (2) They are free from fine tuning of a weight function (probability distribution function or potential of mean force). (3) A canonical ensemble (i.e., a thermally equilibrated ensemble) is generated from a simple procedure. A thermodynamic weight is assigned to each snapshot. (4) Selective sampling can be performed for particularly addressing a poorly sampled region without breaking the proportion of the canonical ensemble if the poorly sampled conformational region emerges in sampling. By applying the methods to a simple system that involves an energy barrier between potential-energy minima, we demonstrated that the new methods have considerably higher sampling efficiency than the original mD-VcMD does.
Collapse
Affiliation(s)
- Junichi Higo
- Graduate School of Simulation Studies, University of Hyogo, Kobe, Hyogo 650-0047, Japan
| | - Ayumi Kusaka
- Institute for Protein Research, Osaka University, Suita, Osaka 565-0871, Japan
| | - Kota Kasahara
- College of Life Sciences, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan
| | - Narutoshi Kamiya
- Graduate School of Simulation Studies, University of Hyogo, Kobe, Hyogo 650-0047, Japan
| | - Itaya Hayato
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan
| | - Xie Qilin
- College of Pharmaceutical Sciences, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan
| | - Takuya Takahashi
- College of Life Sciences, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan
| | - Ikuo Fukuda
- Graduate School of Simulation Studies, University of Hyogo, Kobe, Hyogo 650-0047, Japan
| | - Kentaro Mori
- Graduate School of Simulation Studies, University of Hyogo, Kobe, Hyogo 650-0047, Japan.,National Institute of Technology, Maizuru College, Maizuru, Kyoto 625-8511, Japan
| | - Yutaka Hata
- Graduate School of Simulation Studies, University of Hyogo, Kobe, Hyogo 650-0047, Japan
| | - Yoshifumi Fukunishi
- Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Koto-ku, Tokyo, 135-0064, Japan
| |
Collapse
|
8
|
Abdelsattar AS, Mansour Y, Aboul-Ela F. The Perturbed Free-Energy Landscape: Linking Ligand Binding to Biomolecular Folding. Chembiochem 2021; 22:1499-1516. [PMID: 33351206 DOI: 10.1002/cbic.202000695] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Revised: 12/19/2020] [Indexed: 12/24/2022]
Abstract
The effects of ligand binding on biomolecular conformation are crucial in drug design, enzyme mechanisms, the regulation of gene expression, and other biological processes. Descriptive models such as "lock and key", "induced fit", and "conformation selection" are common ways to interpret such interactions. Another historical model, linked equilibria, proposes that the free-energy landscape (FEL) is perturbed by the addition of ligand binding energy for the bound population of biomolecules. This principle leads to a unified, quantitative theory of ligand-induced conformation change, building upon the FEL concept. We call the map of binding free energy over biomolecular conformational space the "binding affinity landscape" (BAL). The perturbed FEL predicts/explains ligand-induced conformational changes conforming to all common descriptive models. We review recent experimental and computational studies that exemplify the perturbed FEL, with emphasis on RNA. This way of understanding ligand-induced conformation dynamics motivates new experimental and theoretical approaches to ligand design, structural biology and systems biology.
Collapse
Affiliation(s)
- Abdallah S Abdelsattar
- Center for X-Ray Determination of the Structure of Matter, Zewail City of Science and Technology, Ahmed Zewail Road, October Gardens, 12578, Giza, Egypt
| | - Youssef Mansour
- Center for X-Ray Determination of the Structure of Matter, Zewail City of Science and Technology, Ahmed Zewail Road, October Gardens, 12578, Giza, Egypt
| | - Fareed Aboul-Ela
- Center for X-Ray Determination of the Structure of Matter, Zewail City of Science and Technology, Ahmed Zewail Road, October Gardens, 12578, Giza, Egypt
| |
Collapse
|
9
|
Higo J, Kawabata T, Kusaka A, Kasahara K, Kamiya N, Fukuda I, Mori K, Hata Y, Fukunishi Y, Nakamura H. Molecular Interaction Mechanism of a 14-3-3 Protein with a Phosphorylated Peptide Elucidated by Enhanced Conformational Sampling. J Chem Inf Model 2020; 60:4867-4880. [PMID: 32910853 DOI: 10.1021/acs.jcim.0c00551] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Enhanced conformational sampling, a genetic-algorithm-guided multidimensional virtual-system coupled molecular dynamics, can provide equilibrated conformational distributions of a receptor protein and a flexible ligand at room temperature. The distributions provide not only the most stable but also semistable complex structures and propose a ligand-receptor binding process. This method was applied to a system consisting of a receptor protein, 14-3-3ε, and a flexible peptide, phosphorylated myeloid leukemia factor 1 (pMLF1). The results present comprehensive binding pathways of pMLF1 to 14-3-3ε. We identified four thermodynamically stable clusters of MLF1 on the 14-3-3ε surface and free-energy barriers among some clusters. The most stable cluster includes two high-density spots connected by a narrow corridor. When pMLF1 passes the corridor, a salt-bridge relay (switching) related to the phosphorylated residue of pMLF1 occurs. Conformations in one high-density spot are similar to the experimentally determined complex structure. Three-dimensional distributions of residues in the intermolecular interface rationally explain the binding constant changes resulting from the alanine mutation experiment for the residues. We also performed a simulation of nonphosphorylated peptide and 14-3-3ε, which demonstrated that the complex structure was unstable, suggesting that phosphorylation of the peptide is crucially important for binding to 14-3-3ε.
Collapse
Affiliation(s)
- Junichi Higo
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Takeshi Kawabata
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka 565-0871, Japan
| | - Ayumi Kusaka
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka 565-0871, Japan
| | - Kota Kasahara
- College of Life Sciences, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
| | - Narutoshi Kamiya
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Ikuo Fukuda
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Kentaro Mori
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,Maizuru College, 234 Shiroya, Maizuru, Kyoto 625-8511 Japan
| | - Yutaka Hata
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yoshifumi Fukunishi
- Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 2-3-26, Aomi, Koto-ku, Tokyo 135-0064, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka 565-0871, Japan
| |
Collapse
|
10
|
Kasahara K, Terazawa H, Itaya H, Goto S, Nakamura H, Takahashi T, Higo J. myPresto/omegagene 2020: a molecular dynamics simulation engine for virtual-system coupled sampling. Biophys Physicobiol 2020; 17:140-146. [PMID: 33240741 PMCID: PMC7671739 DOI: 10.2142/biophysico.bsj-2020013] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Accepted: 10/10/2020] [Indexed: 12/03/2022] Open
Abstract
The molecular dynamics (MD) method is a promising approach for investigating the molecular mechanisms of microscopic phenomena. In particular, generalized ensemble MD methods can efficiently explore the conformational space with a rugged free-energy surface. However, the implementation and acquisition of technical knowledge for each generalized ensemble MD method are not straightforward for end-users. Here, we present a new version of the myPresto/omegagene software, which is an MD simulation engine tailored for a series of generalized ensemble methods, which are virtual-system coupled multicanonical MD (V-McMD), virtual-system coupled adaptive umbrella sampling (V-AUS), and virtual-system coupled canonical MD (VcMD). This program has been applied in several studies analyzing free-energy landscapes of a variety of molecular systems with all-atom simulations. The updated version provides new functionality for coarse-grained simulations powered by the hydrophobicity scale method. The software package includes a step-by-step tutorial document for enhanced conformational sampling of the poly-glutamine (poly-Q) oligomer expressed as a one-bead per residue model. The myPresto/omegagene software is freely available at the following URL: https://github.com/kotakasahara/omegagene under the Apache2 license.
Collapse
Affiliation(s)
- Kota Kasahara
- College of Life Sciences, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan
| | - Hiroki Terazawa
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan
| | - Hayato Itaya
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan
| | - Satoshi Goto
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, Suita, Osaka 565-0871, Japan
| | - Takuya Takahashi
- College of Life Sciences, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan
| | - Junichi Higo
- Graduate School of Simulation Studies, University of Hyogo, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|