1
|
Bekker GJ, Oshima K, Araki M, Okuno Y, Kamiya N. Binding Mechanism between Platelet Glycoprotein and Cyclic Peptide Elucidated by McMD-Based Dynamic Docking. J Chem Inf Model 2024; 64:4158-4167. [PMID: 38751042 DOI: 10.1021/acs.jcim.4c00100] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/28/2024]
Abstract
The cyclic peptide OS1 (amino acid sequence: CTERMALHNLC), which has a disulfide bond between both termini cysteine residues, inhibits complex formation between the platelet glycoprotein Ibα (GPIbα) and the von Willebrand factor (vWF) by forming a complex with GPIbα. To study the binding mechanism between GPIbα and OS1 and, therefore, the inhibition mechanism of the protein-protein GPIbα-vWF complex, we have applied our multicanonical molecular dynamics (McMD)-based dynamic docking protocol starting from the unbound state of the peptide. Our simulations have reproduced the experimental complex structure, although the top-ranking structure was an intermediary one, where the peptide was bound in the same location as in the experimental structure; however, the β-switch of GPIbα attained a different conformation. Our analysis showed that subsequent refolding of the β-switch results in a more stable binding configuration, although the transition to the native configuration appears to take some time, during which OS1 could dissociate. Our results show that conformational changes in the β-switch are crucial for successful binding of OS1. Furthermore, we identified several allosteric binding sites of GPIbα that might also interfere with vWF binding, and optimization of the peptide to target these allosteric sites might lead to a more effective inhibitor, as these are not dependent on the β-switch conformation.
Collapse
Affiliation(s)
- Gert-Jan Bekker
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Kanji Oshima
- Bio-Pharma Research Laboratories, Kaneka Corporation, 1-8 Miyamae-cho, Takasago-cho, Takasago, Hyogo 676-8688, Japan
| | - Mitsugu Araki
- Graduate School of Medicine, Kyoto University, 53 Shogoin-Kawaharacho, Sakyo-ku, Kyoto 606-8507, Japan
| | - Yasushi Okuno
- Graduate School of Medicine, Kyoto University, 53 Shogoin-Kawaharacho, Sakyo-ku, Kyoto 606-8507, Japan
| | - Narutoshi Kamiya
- Graduate School of Information Science, University of Hyogo, 7-1-28 minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
2
|
Bekker GJ, Fukunishi Y, Higo J, Kamiya N. Binding Mechanism of Riboswitch to Natural Ligand Elucidated by McMD-Based Dynamic Docking Simulations. ACS OMEGA 2024; 9:3412-3422. [PMID: 38284074 PMCID: PMC10809319 DOI: 10.1021/acsomega.3c06826] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Revised: 12/16/2023] [Accepted: 12/28/2023] [Indexed: 01/30/2024]
Abstract
Flavin mononucleotide riboswitches are common among many pathogenic bacteria and are therefore considered to be an attractive target for antibiotics development. The riboswitch binds riboflavin (RBF, also known as vitamin B2), and although an experimental structure of their complex has been solved with the ligand bound deep inside the RNA molecule in a seemingly unreachable state, the binding mechanism between these molecules is not yet known. We have therefore used our Multicanonical Molecular Dynamics (McMD)-based dynamic docking protocol to analyze their binding mechanism by simulating the binding process between the riboswitch aptamer domain and the RBF, starting from the apo state of the riboswitch. Here, the refinement stage was crucial to identify the native binding configuration, as several other binding configurations were also found by McMD-based docking simulations. RBF initially binds the interface between P4 and P6 including U61 and G62, which forms a gateway where the ligand lingers until this gateway opens sufficiently to allow the ligand to pass through and slip into the hidden binding site including A48, A49, and A85.
Collapse
Affiliation(s)
- Gert-Jan Bekker
- Institute
for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, 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
| | - Junichi Higo
- Graduate
School of Information Science, University
of Hyogo, 7-1-28 minatojima
Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Narutoshi Kamiya
- Graduate
School of Information Science, University
of Hyogo, 7-1-28 minatojima
Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
3
|
Bekker GJ, Numoto N, Kawasaki M, Hayashi T, Yabuno S, Kozono Y, Shimizu T, Kozono H, Ito N, Oda M, Kamiya N. Elucidation of binding mechanism, affinity, and complex structure between mWT1 tumor-associated antigen peptide and HLA-A*24:02. Protein Sci 2023; 32:e4775. [PMID: 37661929 PMCID: PMC10510467 DOI: 10.1002/pro.4775] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Revised: 08/02/2023] [Accepted: 08/29/2023] [Indexed: 09/05/2023]
Abstract
We have applied our advanced computational and experimental methodologies to investigate the complex structure and binding mechanism of a modified Wilms' Tumor 1 (mWT1) protein epitope to the understudied Asian-dominant allele HLA-A*24:02 (HLA-A24) in aqueous solution. We have applied our developed multicanonical molecular dynamics (McMD)-based dynamic docking method to analyze the binding pathway and mechanism, which we verified by comparing the highest probability structures from simulation with our experimentally solved x-ray crystal structure. Subsequent path sampling MD simulations elucidated the atomic details of the binding process and indicated that first an encounter complex is formed between the N-terminal's positive charge of the 9-residue mWT1 fragment peptide and a cluster of negative residues on the surface of HLA-A24, with the major histocompatibility complex (MHC) molecule preferring a predominantly closed conformation. The peptide first binds to this closed MHC conformation, forming an encounter complex, after which the binding site opens due to increased entropy of the binding site, allowing the peptide to bind to form the native complex structure. Further sequence and structure analyses also suggest that although the peptide loading complex would help with stabilizing the MHC molecule, the binding depends in a large part on the intrinsic affinity between the MHC molecule and the antigen peptide. Finally, our computational tools and analyses can be of great benefit to study the binding mechanism of different MHC types to their antigens, where it could also be useful in the development of higher affinity variant peptides and for personalized medicine.
Collapse
Affiliation(s)
- Gert-Jan Bekker
- Institute for Protein Research, Osaka University, Suita, Osaka, Japan
| | - Nobutaka Numoto
- Medical Research Institute, Tokyo Medical and Dental University (TMDU), Tokyo, Japan
| | - Maki Kawasaki
- Graduate School of Life and Environmental Sciences, Kyoto Prefectural University, Kyoto, Kyoto, Japan
| | - Takahiro Hayashi
- Graduate School of Life and Environmental Sciences, Kyoto Prefectural University, Kyoto, Kyoto, Japan
| | - Saaya Yabuno
- Graduate School of Life and Environmental Sciences, Kyoto Prefectural University, Kyoto, Kyoto, Japan
| | - Yuko Kozono
- Research Institute for Biomedical Sciences, Tokyo University of Science, Noda, Chiba, Japan
| | - Takeyuki Shimizu
- Department of Immunology, Kochi Medical School, Kochi University, Nankoku-shi, Kochi, Japan
| | - Haruo Kozono
- Research Institute for Biomedical Sciences, Tokyo University of Science, Noda, Chiba, Japan
| | - Nobutoshi Ito
- Medical Research Institute, Tokyo Medical and Dental University (TMDU), Tokyo, Japan
| | - Masayuki Oda
- Graduate School of Life and Environmental Sciences, Kyoto Prefectural University, Kyoto, Kyoto, Japan
| | - Narutoshi Kamiya
- Graduate School of Information Science, University of Hyogo, Kobe, Hyogo, Japan
| |
Collapse
|
4
|
Zhou C, Li J, Wang S, Zhao J, Ai L, Chen Q, Chen Q, Zhao D. Development of Molecular Dynamics and Research Progress in the Study of Slag. MATERIALS (BASEL, SWITZERLAND) 2023; 16:5373. [PMID: 37570076 PMCID: PMC10419983 DOI: 10.3390/ma16155373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/06/2023] [Revised: 07/07/2023] [Accepted: 07/12/2023] [Indexed: 08/13/2023]
Abstract
Molecular dynamics is a method of studying microstructure and properties by calculating and simulating the movement and interaction of molecules. The molecular dynamics simulation method has become an important method for studying the structural and dynamic characteristics of slag systems and can make up for the shortcomings of existing detection methods and experiments. Firstly, this paper analyzes the development process and application fields of molecular dynamics, summarizes the general simulation steps and software algorithms of molecular dynamics simulation methods, and discusses the advantages and disadvantages of the algorithms and the common functions of the software. Secondly, the research status and application progress of molecular dynamics simulation methods in the study of phosphate, silicate, aluminate and aluminosilicate are introduced. On this basis, a method of combining molecular dynamics simulation with laboratory experiments is proposed, which will help obtain more accurate simulation results. This review provides theoretical guidance and a technical framework for the effective analysis of the microstructure of different slag systems via molecular dynamics, so as to finally meet the needs of iron and steel enterprises in producing high-quality steel grades.
Collapse
Affiliation(s)
- Chaogang Zhou
- College of Metallurgy and Energy, North China University of Science and Technology, Tangshan 063210, China; (C.Z.); (J.L.); (L.A.); (Q.C.); (Q.C.); (D.Z.)
- Tangshan Special Metallurgy and Material Preparation Laboratory, Tangshan 063210, China
| | - Jinyue Li
- College of Metallurgy and Energy, North China University of Science and Technology, Tangshan 063210, China; (C.Z.); (J.L.); (L.A.); (Q.C.); (Q.C.); (D.Z.)
- Tangshan Special Metallurgy and Material Preparation Laboratory, Tangshan 063210, China
| | - Shuhuan Wang
- College of Metallurgy and Energy, North China University of Science and Technology, Tangshan 063210, China; (C.Z.); (J.L.); (L.A.); (Q.C.); (Q.C.); (D.Z.)
- Tangshan Special Metallurgy and Material Preparation Laboratory, Tangshan 063210, China
| | - Jingjing Zhao
- College of Pharmaceutical Sciences, North China University of Science and Technology, Tangshan 063210, China;
| | - Liqun Ai
- College of Metallurgy and Energy, North China University of Science and Technology, Tangshan 063210, China; (C.Z.); (J.L.); (L.A.); (Q.C.); (Q.C.); (D.Z.)
- Tangshan Special Metallurgy and Material Preparation Laboratory, Tangshan 063210, China
| | - Qinggong Chen
- College of Metallurgy and Energy, North China University of Science and Technology, Tangshan 063210, China; (C.Z.); (J.L.); (L.A.); (Q.C.); (Q.C.); (D.Z.)
- Tangshan Special Metallurgy and Material Preparation Laboratory, Tangshan 063210, China
| | - Qiya Chen
- College of Metallurgy and Energy, North China University of Science and Technology, Tangshan 063210, China; (C.Z.); (J.L.); (L.A.); (Q.C.); (Q.C.); (D.Z.)
- Tangshan Special Metallurgy and Material Preparation Laboratory, Tangshan 063210, China
| | - Dingguo Zhao
- College of Metallurgy and Energy, North China University of Science and Technology, Tangshan 063210, China; (C.Z.); (J.L.); (L.A.); (Q.C.); (Q.C.); (D.Z.)
- Tangshan Special Metallurgy and Material Preparation Laboratory, Tangshan 063210, China
| |
Collapse
|
5
|
Bekker GJ, Araki M, Oshima K, Okuno Y, Kamiya N. Mutual induced-fit mechanism drives binding between intrinsically disordered Bim and cryptic binding site of Bcl-xL. Commun Biol 2023; 6:349. [PMID: 36997643 PMCID: PMC10063584 DOI: 10.1038/s42003-023-04720-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2022] [Accepted: 03/16/2023] [Indexed: 04/03/2023] Open
Abstract
The intrinsically disordered region (IDR) of Bim binds to the flexible cryptic site of Bcl-xL, a pro-survival protein involved in cancer progression that plays an important role in initiating apoptosis. However, their binding mechanism has not yet been elucidated. We have applied our dynamic docking protocol, which correctly reproduced both the IDR properties of Bim and the native bound configuration, as well as suggesting other stable/meta-stable binding configurations and revealed the binding pathway. Although the cryptic site of Bcl-xL is predominantly in a closed conformation, initial binding of Bim in an encounter configuration leads to mutual induced-fit binding, where both molecules adapt to each other; Bcl-xL transitions to an open state as Bim folds from a disordered to an α-helical conformation while the two molecules bind each other. Finally, our data provides new avenues to develop novel drugs by targeting newly discovered stable conformations of Bcl-xL.
Collapse
Affiliation(s)
- Gert-Jan Bekker
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan.
| | - Mitsugu Araki
- Graduate School of Medicine, Kyoto University, 53 Shogoin-Kawaharacho, Sakyo-ku, Kyoto, 606-8507, Japan
| | - Kanji Oshima
- Bio-Pharma Research Laboratories, KANEKA CORPORATION, 1-8 Miyamae-cho, Takasago-cho, Takasago, Hyogo, 676-8688, Japan
| | - Yasushi Okuno
- Graduate School of Medicine, Kyoto University, 53 Shogoin-Kawaharacho, Sakyo-ku, Kyoto, 606-8507, Japan
| | - Narutoshi Kamiya
- Graduate School of Information Science, University of Hyogo, 7-1-28 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo, 650-0047, Japan.
| |
Collapse
|
6
|
Bekker GJ, Kamiya N. Advancing the field of computational drug design using multicanonical molecular dynamics-based dynamic docking. Biophys Rev 2022; 14:1349-1358. [PMID: 36659995 PMCID: PMC9842809 DOI: 10.1007/s12551-022-01010-z] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2022] [Accepted: 10/14/2022] [Indexed: 11/20/2022] Open
Abstract
Multicanonical molecular dynamics (McMD)-based dynamic docking is a powerful tool to not only predict the native binding configuration between two flexible molecules, but it can also be used to accurately simulate the binding/unbinding pathway. Furthermore, it can also predict alternative binding sites, including allosteric ones, by employing an exhaustive sampling approach. Since McMD-based dynamic docking accurately samples binding/unbinding events, it can thus be used to determine the molecular mechanism of binding between two molecules. We developed the McMD-based dynamic docking methodology based on the powerful, but woefully underutilized McMD algorithm, combined with a toolset to perform the docking and to analyze the results. Here, we showcase three of our recent works, where we have applied McMD-based dynamic docking to advance the field of computational drug design. In the first case, we applied our method to perform an exhaustive search between Hsp90 and one of its inhibitors to successfully predict the native binding configuration in its binding site, as we refined our analysis methods. For our second case, we performed an exhaustive search of two medium-sized ligands and Bcl-xL, which has a cryptic binding site that differs greatly between the apo and holo structures. Finally, we performed a dynamic docking simulation between a membrane-embedded GPCR molecule and a high affinity ligand that binds deep within its receptor's pocket. These advanced simulations showcase the power that the McMD-based dynamic docking method has, and provide a glimpse of the potential our methodology has to unravel and solve the medical and biophysical issues in the modern world. Supplementary Information The online version contains supplementary material available at 10.1007/s12551-022-01010-z.
Collapse
Affiliation(s)
- Gert-Jan Bekker
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871 Japan
| | - Narutoshi Kamiya
- Graduate School of Information Science, University of Hyogo, 7-1-28 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047 Japan
| |
Collapse
|
7
|
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
|