1
|
Xu Q, Yang M, Ji J, Weng J, Wang W, Xu X. Impact of Nonnative Interactions on the Binding Kinetics of Intrinsically Disordered p53 with MDM2: Insights from All-Atom Simulation and Markov State Model Analysis. J Chem Inf Model 2024; 64:5219-5231. [PMID: 38916177 DOI: 10.1021/acs.jcim.3c01833] [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: 06/26/2024]
Abstract
Intrinsically disordered proteins (IDPs) lack a well-defined tertiary structure but are essential players in various biological processes. Their ability to undergo a disorder-to-order transition upon binding to their partners, known as the folding-upon-binding process, is crucial for their function. One classical example is the intrinsically disordered transactivation domain (TAD) of the tumor suppressor protein p53, which quickly forms a structured α-helix after binding to its partner MDM2, with clinical significance for cancer treatment. However, the contribution of nonnative interactions between the IDP and its partner to the rapid binding kinetics, as well as their interplay with native interactions, is not well understood at the atomic level. Here, we used molecular dynamics simulation and Markov state model (MSM) analysis to study the folding-upon-binding mechanism between p53-TAD and MDM2. Our results suggest that the system progresses from the nascent encounter complex to the well-structured encounter complex and finally reaches the native complex, following an induced-fit mechanism. We found that nonnative hydrophobic and hydrogen bond interactions, combined with native interactions, effectively stabilize the nascent and well-structured encounter complexes. Among the nonnative interactions, Leu25p53-Leu54MDM2 and Leu25p53-Phe55MDM2 are particularly noteworthy, as their interaction strength is close to the optimum. Evidently, strengthening or weakening these interactions could both adversely affect the binding kinetics. Overall, our findings suggest that nonnative interactions are evolutionarily optimized to accelerate the binding kinetics of IDPs in conjunction with native interactions.
Collapse
Affiliation(s)
- Qianjun Xu
- Department of Chemistry, Institute of Biomedical Sciences and Multiscale Research Institute of Complex Systems, Fudan University, Shanghai 200438, China
| | - Maohua Yang
- Department of Chemistry, Institute of Biomedical Sciences and Multiscale Research Institute of Complex Systems, Fudan University, Shanghai 200438, China
| | - Jie Ji
- Department of Chemistry, Institute of Biomedical Sciences and Multiscale Research Institute of Complex Systems, Fudan University, Shanghai 200438, China
| | - Jingwei Weng
- Department of Chemistry, Institute of Biomedical Sciences and Multiscale Research Institute of Complex Systems, Fudan University, Shanghai 200438, China
| | - Wenning Wang
- Department of Chemistry, Institute of Biomedical Sciences and Multiscale Research Institute of Complex Systems, Fudan University, Shanghai 200438, China
| | - Xin Xu
- Department of Chemistry, Institute of Biomedical Sciences and Multiscale Research Institute of Complex Systems, Fudan University, Shanghai 200438, China
| |
Collapse
|
2
|
Jin L, Liu Y, Wu Y, Huang Y, Zhang D. REST Is Not Resting: REST/NRSF in Health and Disease. Biomolecules 2023; 13:1477. [PMID: 37892159 PMCID: PMC10605157 DOI: 10.3390/biom13101477] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2023] [Revised: 09/26/2023] [Accepted: 09/29/2023] [Indexed: 10/29/2023] Open
Abstract
Chromatin modifications play a crucial role in the regulation of gene expression. The repressor element-1 (RE1) silencing transcription factor (REST), also known as neuron-restrictive silencer factor (NRSF) and X2 box repressor (XBR), was found to regulate gene transcription by binding to chromatin and recruiting chromatin-modifying enzymes. Earlier studies revealed that REST plays an important role in the development and disease of the nervous system, mainly by repressing the transcription of neuron-specific genes. Subsequently, REST was found to be critical in other tissues, such as the heart, pancreas, skin, eye, and vascular. Dysregulation of REST was also found in nervous and non-nervous system cancers. In parallel, multiple strategies to target REST have been developed. In this paper, we provide a comprehensive summary of the research progress made over the past 28 years since the discovery of REST, encompassing both physiological and pathological aspects. These insights into the effects and mechanisms of REST contribute to an in-depth understanding of the transcriptional regulatory mechanisms of genes and their roles in the development and progression of disease, with a view to discovering potential therapeutic targets and intervention strategies for various related diseases.
Collapse
Affiliation(s)
- Lili Jin
- School of Life Sciences, Liaoning University, Shenyang 110036, China
| | - Ying Liu
- Department of Stem Cells and Regenerative Medicine, Key Laboratory of Cell Biology, National Health Commission of China, and Key Laboratory of Medical Cell Biology, Ministry of Education of China, China Medical University, Shenyang 110122, China
| | - Yifan Wu
- Department of Stem Cells and Regenerative Medicine, Key Laboratory of Cell Biology, National Health Commission of China, and Key Laboratory of Medical Cell Biology, Ministry of Education of China, China Medical University, Shenyang 110122, China
| | - Yi Huang
- Department of Stem Cells and Regenerative Medicine, Key Laboratory of Cell Biology, National Health Commission of China, and Key Laboratory of Medical Cell Biology, Ministry of Education of China, China Medical University, Shenyang 110122, China
| | - Dianbao Zhang
- Department of Stem Cells and Regenerative Medicine, Key Laboratory of Cell Biology, National Health Commission of China, and Key Laboratory of Medical Cell Biology, Ministry of Education of China, China Medical University, Shenyang 110122, China
| |
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
|
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
|
5
|
Kondo HX, Nakamura H, Takano Y. Negative fragmentation approach for investigating the depolarization effect of neighboring residues on hydrogen bonds in π-helix. Chem Phys Lett 2023. [DOI: 10.1016/j.cplett.2023.140361] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
|
6
|
Quantum chemical studies on hydrogen bonds in helical secondary structures. Biophys Rev 2023; 14:1369-1378. [PMID: 36659988 PMCID: PMC9842822 DOI: 10.1007/s12551-022-01034-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Accepted: 11/30/2022] [Indexed: 01/07/2023] Open
Abstract
We present a brief review of our recent computational studies of hydrogen bonds (H-bonds) in helical secondary structures of proteins, α-helix and 310-helix, using a Negative Fragmentation Approach with density functional theory. We found that the depolarized electronic structures of the carbonyl oxygen of the ith residue and the amide hydrogen of the (i + 4)th residue cause weaker H-bond in an α-helix than in an isolated H-bond. Our calculations showed that the H-bond energies in the 310-helix were also weaker than those of the isolated H-bonds. In the 310-helices, the adjacent N-H group at the (i + 1)th residue was closer to the C=O group of the H-bond pair than the adjacent C=O group in the 310-helices, whereas the adjacent C=O group at the (i + 1)th residue was close to the H-bond acceptor in α-helices. Therefore, the destabilization of the H-bond is attributed to the depolarization caused by the adjacent residue of the helical backbone connecting the H-bond donor and acceptor. The differences in the change in electron density revealed that such depolarizations were caused by the local electronic interactions in their neighborhood inside the helical structure and redistributed the electron density. We also present the improvements in the force field of classical molecular simulation, based on our findings. Supplementary Information The online version contains supplementary material available at 10.1007/s12551-022-01034-5.
Collapse
|
7
|
Nakamura H. Some reflections on a career in science and a note of thanks to the contributors of this Special Issue. Biophys Rev 2022; 14:1223-1226. [PMID: 36659991 PMCID: PMC9842830 DOI: 10.1007/s12551-022-01035-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/08/2022] [Indexed: 12/24/2022] Open
Affiliation(s)
- Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871 Japan
| |
Collapse
|
8
|
Okuda M, Tsunaka Y, Nishimura Y. Dynamic structures of intrinsically disordered proteins related to the general transcription factor TFIIH, nucleosomes, and histone chaperones. Biophys Rev 2022; 14:1449-1472. [PMID: 36659983 PMCID: PMC9842849 DOI: 10.1007/s12551-022-01014-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Accepted: 11/06/2022] [Indexed: 11/19/2022] Open
Abstract
Advances in structural analysis by cryogenic electron microscopy (cryo-EM) and X-ray crystallography have revealed the tertiary structures of various chromatin-related proteins, including transcription factors, RNA polymerases, nucleosomes, and histone chaperones; however, the dynamic structures of intrinsically disordered regions (IDRs) in these proteins remain elusive. Recent studies using nuclear magnetic resonance (NMR), together with molecular dynamics (MD) simulations, are beginning to reveal dynamic structures of the general transcription factor TFIIH complexed with target proteins including the general transcription factor TFIIE, the tumor suppressor p53, the cell cycle protein DP1, the DNA repair factors XPC and UVSSA, and three RNA polymerases, in addition to the dynamics of histone tails in nucleosomes and histone chaperones. In complexes of TFIIH, the PH domain of the p62 subunit binds to an acidic string formed by the IDR in TFIIE, p53, XPC, UVSSA, DP1, and the RPB6 subunit of three RNA polymerases by a common interaction mode, namely extended string-like binding of the IDR on the positively charged surface of the PH domain. In the nucleosome, the dynamic conformations of the N-tails of histones H2A and H2B are correlated, while the dynamic conformations of the N-tails of H3 and H4 form a histone tail network dependent on their modifications and linker DNA. The acidic IDRs of the histone chaperones of FACT and NAP1 play important roles in regulating the accessibility to histone proteins in the nucleosome.
Collapse
Affiliation(s)
- Masahiko Okuda
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehiro-Cho, Tsurumi-Ku, Yokohama, 230-0045 Japan
| | - Yasuo Tsunaka
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehiro-Cho, Tsurumi-Ku, Yokohama, 230-0045 Japan
| | - Yoshifumi Nishimura
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehiro-Cho, Tsurumi-Ku, Yokohama, 230-0045 Japan
- Graduate School of Integrated Sciences for Life, Hiroshima University, 1-4-4 Kagamiyama, Higashi-Hiroshima, 739-8528 Japan
| |
Collapse
|
9
|
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
|
10
|
Depolarizing Effects in Hydrogen Bond Energy in 3 10-Helices Revealed by Quantum Chemical Analysis. Int J Mol Sci 2022; 23:ijms23169032. [PMID: 36012292 PMCID: PMC9409261 DOI: 10.3390/ijms23169032] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2022] [Revised: 08/10/2022] [Accepted: 08/11/2022] [Indexed: 12/18/2022] Open
Abstract
Hydrogen-bond (H-bond) energies in 310-helices of short alanine peptides were systematically examined by precise DFT calculations with the negative fragmentation approach (NFA), a modified method based on the molecular tailoring approach. The contribution of each H-bond was evaluated in detail from the 310-helical conformation of total energies (whole helical model, WH3-10 model), and the results were compared with the property of H-bond in α-helix from our previous study. The H-bond energies of the WH3-10 model exhibited tendencies different from those exhibited by the α-helix in that they depended on the helical position of the relevant H-bond pair. H-bond pairs adjacent to the terminal H-bond pairs were observed to be strongly destabilized. The analysis of electronic structures indicated that structural characteristics cause the destabilization of the H-bond in 310-helices. We also found that the longer the helix length, the more stable the H-bond in the terminal pairs of the WH3-10 model, suggesting the action of H-bond cooperativity.
Collapse
|
11
|
Sato M. Biological Significance of Intrinsically Disordered Protein Structure. CHEM-BIO INFORMATICS JOURNAL 2022. [DOI: 10.1273/cbij.22.26] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Affiliation(s)
- Mamoru Sato
- Graduate School of Medical Life Science, Yokohama City University
| |
Collapse
|
12
|
Constructing conformational library for amyloid-β42 dimers as the smallest toxic oligomers using two CHARMM force fields. J Mol Graph Model 2022; 115:108207. [DOI: 10.1016/j.jmgm.2022.108207] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2021] [Revised: 04/21/2022] [Accepted: 04/23/2022] [Indexed: 11/19/2022]
|
13
|
Ueda T, Imai S, Shimada I. Function-related dynamics of GPCRs. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2022; 336:107164. [PMID: 35168190 DOI: 10.1016/j.jmr.2022.107164] [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: 11/16/2021] [Revised: 01/29/2022] [Accepted: 02/03/2022] [Indexed: 06/14/2023]
Abstract
G protein-coupled receptors (GPCRs) include various neurotransmitters and hormones, and over 30% of modern drugs target GPCRs. The number of GPCR crystal structures has rapidly increased, and many structures of GPCRs in complexes with their binding partners are being solved by cryo-electron microscopy. However, crystallographic or cryo-electron microscopy data alone cannot fully explain the important features of GPCR signaling determined experimentally. Recent studies have suggested that GPCRs are structurally dynamic, and exchange between multiple conformations. In this respect, NMR methods provide information about the dynamics of proteins over a wide range of frequencies, in aqueous solutions at nearphysiological temperatures. Although NMR studies of GPCRs are challenging due to their innate instability and relatively large molecular weights, recent methodological advances have enabled us to observe the NMR signals of various GPCRs. These NMR studies revealed that GPCRs exist in function-related equilibria between locally different conformations that are simultaneously populated. Here we will describe solution NMR studies that have clarified the function-related conformational dynamics of two GPCRs, β2 adrenergic receptor and adenosine A2A receptor.
Collapse
Affiliation(s)
- Takumi Ueda
- Graduate School of Pharmaceutical Sciences, The University of Tokyo, Tokyo, Japan
| | - Shunsuke Imai
- Center for Biosystems Dynamics Research (BDR), RIKEN, Kanagawa, Japan
| | - Ichio Shimada
- Graduate School of Pharmaceutical Sciences, The University of Tokyo, Tokyo, Japan; Center for Biosystems Dynamics Research (BDR), RIKEN, Kanagawa, Japan.
| |
Collapse
|
14
|
Kamenik AS, Linker SM, Riniker S. Enhanced sampling without borders: on global biasing functions and how to reweight them. Phys Chem Chem Phys 2022; 24:1225-1236. [PMID: 34935813 PMCID: PMC8768491 DOI: 10.1039/d1cp04809k] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Accepted: 12/14/2021] [Indexed: 12/17/2022]
Abstract
Molecular dynamics (MD) simulations are a powerful tool to follow the time evolution of biomolecular motions in atomistic resolution. However, the high computational demand of these simulations limits the timescales of motions that can be observed. To resolve this issue, so called enhanced sampling techniques are developed, which extend conventional MD algorithms to speed up the simulation process. Here, we focus on techniques that apply global biasing functions. We provide a broad overview of established enhanced sampling methods and promising new advances. As the ultimate goal is to retrieve unbiased information from biased ensembles, we also discuss benefits and limitations of common reweighting schemes. In addition to concisely summarizing critical assumptions and implications, we highlight the general application opportunities as well as uncertainties of global enhanced sampling.
Collapse
Affiliation(s)
- Anna S Kamenik
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
| | - Stephanie M Linker
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
| | - Sereina Riniker
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
| |
Collapse
|
15
|
Bekker GJ, Kamiya N. N-Terminal-Driven Binding Mechanism of an Antigen Peptide to Human Leukocyte Antigen-A*2402 Elucidated by Multicanonical Molecular Dynamic-Based Dynamic Docking and Path Sampling Simulations. J Phys Chem B 2021; 125:13376-13384. [PMID: 34856806 DOI: 10.1021/acs.jpcb.1c07230] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have applied our advanced multicanonical molecular dynamics (McMD)-based dynamic docking methodology to investigate the binding mechanism of an HIV-1 Nef protein epitope to the Asian-dominant allele human leukocyte antigen (HLA)-A*2402. Even though pMHC complex formation [between a Major histocompatibility complex (MHC) class I molecule, which is encoded by an HLA allele, and an antigen peptide] is one of the fundamental processes of the adaptive human immune response, its binding mechanism has not yet been well studied, partially due to the high allelic variation of HLAs in the population. We have used our developed McMD-based dynamic docking method and have successfully reproduced the native complex structure, which is located near the free energy global minimum. Subsequent path sampling MD simulations elucidated the atomic details of the binding process and indicated that the peptide binding is initially driven by the highly positively charged N-terminus of the peptide that is attracted to the various negatively charged residues on the MHC molecule's surface. Upon nearing the pocket, the second tyrosine residue of the peptide anchors the peptide by strongly binding to the B-site of the MHC molecule via hydrophobic driven interactions, resulting in a very strong bound complex structure. Our methodology can be effectively used to predict the bound complex structures between MHC molecules and their antigens to study their binding mechanism in close detail, which would help with the development of new vaccines against cancers, as well as viral infections such as HIV and COVID-19.
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
|
16
|
Matsuo N, Goda N, Tenno T, Hiroaki H. Cryoprotective activities of FK20, a human genome-derived intrinsically disordered peptide against cryosensitive enzymes without a stereospecific molecular interaction. PEERJ PHYSICAL CHEMISTRY 2021. [DOI: 10.7717/peerj-pchem.20] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Background
Intrinsically disordered proteins (IDPs) have been shown to exhibit cryoprotective activity toward other cellular enzymes without any obvious conserved sequence motifs. This study investigated relationships between the physical properties of several human genome-derived IDPs and their cryoprotective activities.
Methods
Cryoprotective activity of three human-genome derived IDPs and their truncated peptides toward lactate dehydrogenase (LDH) and glutathione S-transferase (GST) was examined. After the shortest cryoprotective peptide was defined (named FK20), cryoprotective activity of all-D-enantiomeric isoform of FK20 (FK20-D) as well as a racemic mixture of FK20 and FK20-D was examined. In order to examine the lack of increase of thermal stability of the target enzyme, the CD spectra of GST and LDH in the presence of a racemic mixture of FK20 and FK20-D at varying temperatures were measured and used to estimate Tm.
Results
Cryoprotective activity of IDPs longer than 20 amino acids was nearly independent of the amino acid length. The shortest IDP-derived 20 amino acid length peptide with sufficient cryoprotective activity was developed from a series of TNFRSF11B fragments (named FK20). FK20, FK20-D, and an equimolar mixture of FK20 and FK20-D also showed similar cryoprotective activity toward LDH and GST. Tm of GST in the presence and absence of an equimolar mixture of FK20 and FK20-D are similar, suggesting that IDPs’ cryoprotection mechanism seems partly from a molecular shielding effect rather than a direct interaction with the target enzymes.
Collapse
Affiliation(s)
- Naoki Matsuo
- Graduate School of Pharmaceutical Sciences, Laboratory of Structural and Molecular Pharmacology, Nagoya University, Nagoya, AICHI, JAPAN
| | - Natsuko Goda
- Graduate School of Pharmaceutical Sciences, Laboratory of Structural and Molecular Pharmacology, Nagoya University, Nagoya, AICHI, JAPAN
| | - Takeshi Tenno
- Graduate School of Pharmaceutical Sciences, Laboratory of Structural and Molecular Pharmacology, Nagoya University, Nagoya, AICHI, JAPAN
- BeCellBar, LLC., Nagoya, Aichi, Japan
| | - Hidekazu Hiroaki
- Graduate School of Pharmaceutical Sciences, Laboratory of Structural and Molecular Pharmacology, Nagoya University, Nagoya, AICHI, JAPAN
- BeCellBar, LLC., Nagoya, Aichi, Japan
| |
Collapse
|
17
|
Ono S, Naylor MR, Townsend CE, Okumura C, Okada O, Lee HW, Lokey RS. Cyclosporin A: Conformational Complexity and Chameleonicity. J Chem Inf Model 2021; 61:5601-5613. [PMID: 34672629 DOI: 10.1021/acs.jcim.1c00771] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The chameleonic behavior of cyclosporin A (CsA) was investigated through conformational ensembles employing multicanonical molecular dynamics simulations that could sample the cis and trans isomers of N-methylated amino acids; these assessments were conducted in explicit water, dimethyl sulfoxide, acetonitrile, methanol, chloroform, cyclohexane (CHX), and n-hexane (HEX) using AMBER ff03, AMBER10:EHT, AMBER12:EHT, and AMBER14:EHT force fields. The conformational details were discussed employing the free-energy landscapes (FELs) at T = 300 K; it was observed that the experimentally determined structures of CsA were only a part of the conformational space. Comparing the ROESY measurements in CHX-d12 and HEX-d14, the major conformations in those apolar solvents were essentially the same as that in CDCl3 except for the observation of some sidechain rotamers. The effects of the metal ions on the conformations, including the cis/trans isomerization, were also investigated. Based on the analysis of FELs, it was concluded that the AMBER ff03 force field best described the experimentally derived conformations, indicating that CsA intrinsically formed membrane-permeable conformations and that the metal ions might be the key to the cis/trans isomerization of N-methylated amino acids before binding a partner protein.
Collapse
Affiliation(s)
- Satoshi Ono
- Modality Laboratories, Innovative Research Division, Mitsubishi Tanabe Pharma Corporation, 1000 Kamoshida-cho, Aoba-ku, Yokohama, Kanagawa 227-0033, Japan
| | - Matthew R Naylor
- Department of Chemistry and Biochemistry, University of California Santa Cruz, 1156 High Street, Santa Cruz, California 95064, United States
| | - Chad E Townsend
- Department of Chemistry and Biochemistry, University of California Santa Cruz, 1156 High Street, Santa Cruz, California 95064, United States
| | - Chieko Okumura
- Modality Laboratories, Innovative Research Division, Mitsubishi Tanabe Pharma Corporation, 1000 Kamoshida-cho, Aoba-ku, Yokohama, Kanagawa 227-0033, Japan
| | - Okimasa Okada
- Modality Laboratories, Innovative Research Division, Mitsubishi Tanabe Pharma Corporation, 1000 Kamoshida-cho, Aoba-ku, Yokohama, Kanagawa 227-0033, Japan
| | - Hsiau-Wei Lee
- Department of Chemistry and Biochemistry, University of California Santa Cruz, 1156 High Street, Santa Cruz, California 95064, United States
| | - R Scott Lokey
- Department of Chemistry and Biochemistry, University of California Santa Cruz, 1156 High Street, Santa Cruz, California 95064, United States
| |
Collapse
|
18
|
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
|
19
|
Kawasaki M, Oda M. DNA-binding function of c-Myb R2R3 around thermal denaturation temperature. Biophys Physicobiol 2021; 18:78-84. [PMID: 33977005 PMCID: PMC8056152 DOI: 10.2142/biophysico.bppb-v18.009] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2021] [Accepted: 03/22/2021] [Indexed: 01/23/2023] Open
Abstract
The minimum DNA-binding domain of the transcriptional factor c-Myb R2R3 remarkably fluctuates in the solution. In the present study, we evaluated the protein fluctuation of R2R3 C130I mutant, R2R3*, on its DNA-binding and folding thermodynamics. DNA-binding analysis using isothermal titration calorimetry revealed that the heat capacity change determined from the correlation between temperature and binding enthalpy change is highly negative above 35°C, indicating that the fluctuation increases with increasing temperature and elevates the conformational change on DNA binding. The results were in accordance with those of differential scanning calorimetry, which revealed that the heat capacity corresponding to thermal denaturation gradually increased above 35°C, followed by the broad transition peak. In contrast, the transition peak of R2R3* in the DNA-bound state was sharper and larger than that in the DNA-unbound state. The fluctuating form could transform into lesser fluctuating form upon DNA binding, resulting in a larger enthalpy change for denaturation of R2R3* in the DNA-bound state. It should also be noted that R2R3* could specifically bind to DNA around thermal denaturation temperature. This would be due to proteins with numerous fluctuations. Moreover, we discuss specific and non-specific DNA binding accompanied by the conformational change between well-ordered and disordered forms of R2R3* observed around the denaturation temperature.
Collapse
Affiliation(s)
- Maki Kawasaki
- Graduate School of Life and Environmental Sciences, Kyoto Prefectural University, Kyoto 606-8522, Japan
| | - Masayuki Oda
- Graduate School of Life and Environmental Sciences, Kyoto Prefectural University, Kyoto 606-8522, Japan
| |
Collapse
|
20
|
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
|
21
|
Damjanovic J, Miao J, Huang H, Lin YS. Elucidating Solution Structures of Cyclic Peptides Using Molecular Dynamics Simulations. Chem Rev 2021; 121:2292-2324. [PMID: 33426882 DOI: 10.1021/acs.chemrev.0c01087] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Protein-protein interactions are vital to biological processes, but the shape and size of their interfaces make them hard to target using small molecules. Cyclic peptides have shown promise as protein-protein interaction modulators, as they can bind protein surfaces with high affinity and specificity. Dozens of cyclic peptides are already FDA approved, and many more are in various stages of development as immunosuppressants, antibiotics, antivirals, or anticancer drugs. However, most cyclic peptide drugs so far have been natural products or derivatives thereof, with de novo design having proven challenging. A key obstacle is structural characterization: cyclic peptides frequently adopt multiple conformations in solution, which are difficult to resolve using techniques like NMR spectroscopy. The lack of solution structural information prevents a thorough understanding of cyclic peptides' sequence-structure-function relationship. Here we review recent development and application of molecular dynamics simulations with enhanced sampling to studying the solution structures of cyclic peptides. We describe novel computational methods capable of sampling cyclic peptides' conformational space and provide examples of computational studies that relate peptides' sequence and structure to biological activity. We demonstrate that molecular dynamics simulations have grown from an explanatory technique to a full-fledged tool for systematic studies at the forefront of cyclic peptide therapeutic design.
Collapse
Affiliation(s)
- Jovan Damjanovic
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, United States
| | - Jiayuan Miao
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, United States
| | - He Huang
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, United States
| | - Yu-Shan Lin
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, United States
| |
Collapse
|
22
|
Umezawa K, Kii I. Druggable Transient Pockets in Protein Kinases. Molecules 2021; 26:molecules26030651. [PMID: 33513739 PMCID: PMC7865889 DOI: 10.3390/molecules26030651] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2020] [Revised: 01/23/2021] [Accepted: 01/26/2021] [Indexed: 12/29/2022] Open
Abstract
Drug discovery using small molecule inhibitors is reaching a stalemate due to low selectivity, adverse off-target effects and inevitable failures in clinical trials. Conventional chemical screening methods may miss potent small molecules because of their use of simple but outdated kits composed of recombinant enzyme proteins. Non-canonical inhibitors targeting a hidden pocket in a protein have received considerable research attention. Kii and colleagues identified an inhibitor targeting a transient pocket in the kinase DYRK1A during its folding process and termed it FINDY. FINDY exhibits a unique inhibitory profile; that is, FINDY does not inhibit the fully folded form of DYRK1A, indicating that the FINDY-binding pocket is hidden in the folded form. This intriguing pocket opens during the folding process and then closes upon completion of folding. In this review, we discuss previously established kinase inhibitors and their inhibitory mechanisms in comparison with FINDY. We also compare the inhibitory mechanisms with the growing concept of “cryptic inhibitor-binding sites.” These sites are buried on the inhibitor-unbound surface but become apparent when the inhibitor is bound. In addition, an alternative method based on cell-free protein synthesis of protein kinases may allow the discovery of small molecules that occupy these mysterious binding sites. Transitional folding intermediates would become alternative targets in drug discovery, enabling the efficient development of potent kinase inhibitors.
Collapse
Affiliation(s)
- Koji Umezawa
- Department of Biomolecular Innovation, Institute for Biomedical Sciences, Shinshu University, 8304 Minami-Minowa, Kami-ina, Nagano 399-4598, Japan;
| | - Isao Kii
- Laboratory for Drug Target Research, Faculty & Graduate School of Agriculture, Shinshu University, 8304 Minami-Minowa, Kami-ina, Nagano 399-4598, Japan
- Correspondence: ; Tel.: +81-265-77-1521
| |
Collapse
|
23
|
Staby L, Bugge K, Falbe-Hansen RG, Salladini E, Skriver K, Kragelund BB. Connecting the αα-hubs: same fold, disordered ligands, new functions. Cell Commun Signal 2021; 19:2. [PMID: 33407551 PMCID: PMC7788954 DOI: 10.1186/s12964-020-00686-8] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2020] [Accepted: 11/18/2020] [Indexed: 11/25/2022] Open
Abstract
Background Signal fidelity depends on protein–protein interaction–‘hubs’ integrating cues from large interactomes. Recently, and based on a common secondary structure motif, the αα-hubs were defined, which are small α-helical domains of large, modular proteins binding intrinsically disordered transcriptional regulators.
Methods Comparative structural biology. Results We assign the harmonin-homology-domain (HHD, also named the harmonin N-terminal domain, NTD) present in large proteins such as harmonin, whirlin, cerebral cavernous malformation 2, and regulator of telomere elongation 1 to the αα-hubs. The new member of the αα-hubs expands functionality to include scaffolding of supra-modular complexes mediating sensory perception, neurovascular integrity and telomere regulation, and reveal novel features of the αα-hubs. As a common trait, the αα-hubs bind intrinsically disordered ligands of similar properties integrating similar cellular cues, but without cross-talk. Conclusion The inclusion of the HHD in the αα-hubs has uncovered new features, exemplifying the utility of identifying groups of hub domains, whereby discoveries in one member may cross-fertilize discoveries in others. These features make the αα-hubs unique models for decomposing signal specificity and fidelity. Using these as models, together with other suitable hub domain, we may advance the functional understanding of hub proteins and their role in cellular communication and signaling, as well as the role of intrinsically disordered proteins in signaling networks. Video Abstract
Collapse
Affiliation(s)
- Lasse Staby
- REPIN, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Katrine Bugge
- REPIN, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | | | - Edoardo Salladini
- REPIN, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Karen Skriver
- REPIN, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| | - Birthe B Kragelund
- REPIN, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| |
Collapse
|
24
|
Bugge K, Staby L, Salladini E, Falbe-Hansen RG, Kragelund BB, Skriver K. αα-Hub domains and intrinsically disordered proteins: A decisive combo. J Biol Chem 2021; 296:100226. [PMID: 33361159 PMCID: PMC7948954 DOI: 10.1074/jbc.rev120.012928] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2020] [Revised: 12/22/2020] [Accepted: 12/22/2020] [Indexed: 01/02/2023] Open
Abstract
Hub proteins are central nodes in protein-protein interaction networks with critical importance to all living organisms. Recently, a new group of folded hub domains, the αα-hubs, was defined based on a shared αα-hairpin supersecondary structural foundation. The members PAH, RST, TAFH, NCBD, and HHD are found in large proteins such as Sin3, RCD1, TAF4, CBP, and harmonin, which organize disordered transcriptional regulators and membrane scaffolds in interactomes of importance to human diseases and plant quality. In this review, studies of structures, functions, and complexes across the αα-hubs are described and compared to provide a unified description of the group. This analysis expands the associated molecular concepts of "one domain-one binding site", motif-based ligand binding, and coupled folding and binding of intrinsically disordered ligands to additional concepts of importance to signal fidelity. These include context, motif reversibility, multivalency, complex heterogeneity, synergistic αα-hub:ligand folding, accessory binding sites, and supramodules. We propose that these multifaceted protein-protein interaction properties are made possible by the characteristics of the αα-hub fold, including supersite properties, dynamics, variable topologies, accessory helices, and malleability and abetted by adaptability of the disordered ligands. Critically, these features provide additional filters for specificity. With the presentations of new concepts, this review opens for new research questions addressing properties across the group, which are driven from concepts discovered in studies of the individual members. Combined, the members of the αα-hubs are ideal models for deconvoluting signal fidelity maintained by folded hubs and their interactions with intrinsically disordered ligands.
Collapse
Affiliation(s)
- Katrine Bugge
- REPIN and The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark; Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Lasse Staby
- REPIN and The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark; Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Edoardo Salladini
- REPIN and The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Rasmus G Falbe-Hansen
- REPIN and The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Birthe B Kragelund
- REPIN and The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark; Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| | - Karen Skriver
- REPIN and The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| |
Collapse
|
25
|
Gerlach GJ, Carrock R, Stix R, Stollar EJ, Ball KA. A disordered encounter complex is central to the yeast Abp1p SH3 domain binding pathway. PLoS Comput Biol 2020; 16:e1007815. [PMID: 32925900 PMCID: PMC7514057 DOI: 10.1371/journal.pcbi.1007815] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2020] [Revised: 09/24/2020] [Accepted: 08/15/2020] [Indexed: 12/20/2022] Open
Abstract
Protein-protein interactions are involved in a wide range of cellular processes. These interactions often involve intrinsically disordered proteins (IDPs) and protein binding domains. However, the details of IDP binding pathways are hard to characterize using experimental approaches, which can rarely capture intermediate states present at low populations. SH3 domains are common protein interaction domains that typically bind proline-rich disordered segments and are involved in cell signaling, regulation, and assembly. We hypothesized, given the flexibility of SH3 binding peptides, that their binding pathways include multiple steps important for function. Molecular dynamics simulations were used to characterize the steps of binding between the yeast Abp1p SH3 domain (AbpSH3) and a proline-rich IDP, ArkA. Before binding, the N-terminal segment 1 of ArkA is pre-structured and adopts a polyproline II helix, while segment 2 of ArkA (C-terminal) adopts a 310 helix, but is far less structured than segment 1. As segment 2 interacts with AbpSH3, it becomes more structured, but retains flexibility even in the fully engaged state. Binding simulations reveal that ArkA enters a flexible encounter complex before forming the fully engaged bound complex. In the encounter complex, transient nonspecific hydrophobic and long-range electrostatic contacts form between ArkA and the binding surface of SH3. The encounter complex ensemble includes conformations with segment 1 in both the forward and reverse orientation, suggesting that segment 2 may play a role in stabilizing the correct binding orientation. While the encounter complex forms quickly, the slow step of binding is the transition from the disordered encounter ensemble to the fully engaged state. In this transition, ArkA makes specific contacts with AbpSH3 and buries more hydrophobic surface. Simulating the binding between ApbSH3 and ArkA provides insight into the role of encounter complex intermediates and nonnative hydrophobic interactions for other SH3 domains and IDPs in general.
Collapse
Affiliation(s)
- Gabriella J. Gerlach
- Department of Chemistry, Skidmore College, Saratoga Springs, New York, United States
| | - Rachel Carrock
- Department of Chemistry, Skidmore College, Saratoga Springs, New York, United States
| | - Robyn Stix
- Department of Chemistry, Skidmore College, Saratoga Springs, New York, United States
| | - Elliott J. Stollar
- School of Life Sciences, University of Liverpool, Liverpool, United Kingdom
| | - K. Aurelia Ball
- Department of Chemistry, Skidmore College, Saratoga Springs, New York, United States
| |
Collapse
|
26
|
Sala D, Cosentino U, Ranaudo A, Greco C, Moro G. Dynamical Behavior and Conformational Selection Mechanism of the Intrinsically Disordered Sic1 Kinase-Inhibitor Domain. Life (Basel) 2020; 10:life10070110. [PMID: 32664566 PMCID: PMC7399826 DOI: 10.3390/life10070110] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2020] [Revised: 07/02/2020] [Accepted: 07/08/2020] [Indexed: 01/04/2023] Open
Abstract
Intrinsically Disordered Peptides and Proteins (IDPs) in solution can span a broad range of conformations that often are hard to characterize by both experimental and computational methods. However, obtaining a significant representation of the conformational space is important to understand mechanisms underlying protein functions such as partner recognition. In this work, we investigated the behavior of the Sic1 Kinase-Inhibitor Domain (KID) in solution by Molecular Dynamics (MD) simulations. Our results point out that application of common descriptors of molecular shape such as Solvent Accessible Surface (SAS) area can lead to misleading outcomes. Instead, more appropriate molecular descriptors can be used to define 3D structures. In particular, we exploited Weighted Holistic Invariant Molecular (WHIM) descriptors to get a coarse-grained but accurate definition of the variegated Sic1 KID conformational ensemble. We found that Sic1 is able to form a variable amount of folded structures even in absence of partners. Among them, there were some conformations very close to the structure that Sic1 is supposed to assume in the binding with its physiological complexes. Therefore, our results support the hypothesis that this protein relies on the conformational selection mechanism to recognize the correct molecular partners.
Collapse
Affiliation(s)
- Davide Sala
- Dipartimento di Biotecnologie e Bioscienze, Università di Milano-Bicocca, P.zza della Scienza 2, 20126 Milano, Italy;
| | - Ugo Cosentino
- Dipartimento di Scienze dell’Ambiente e della Terra, Università di Milano-Bicocca, P.zza della Scienza 1, 20126 Milano, Italy; (U.C.); (A.R.)
| | - Anna Ranaudo
- Dipartimento di Scienze dell’Ambiente e della Terra, Università di Milano-Bicocca, P.zza della Scienza 1, 20126 Milano, Italy; (U.C.); (A.R.)
| | - Claudio Greco
- Dipartimento di Scienze dell’Ambiente e della Terra, Università di Milano-Bicocca, P.zza della Scienza 1, 20126 Milano, Italy; (U.C.); (A.R.)
- Correspondence: (C.G.); (G.M.)
| | - Giorgio Moro
- Dipartimento di Biotecnologie e Bioscienze, Università di Milano-Bicocca, P.zza della Scienza 2, 20126 Milano, Italy;
- Correspondence: (C.G.); (G.M.)
| |
Collapse
|
27
|
Ezerski JC, Zhang P, Jennings NC, Waxham MN, Cheung MS. Molecular Dynamics Ensemble Refinement of Intrinsically Disordered Peptides According to Deconvoluted Spectra from Circular Dichroism. Biophys J 2020; 118:1665-1678. [PMID: 32145192 PMCID: PMC7136346 DOI: 10.1016/j.bpj.2020.02.015] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 02/07/2020] [Accepted: 02/18/2020] [Indexed: 02/08/2023] Open
Abstract
We have developed a computational method of atomistically refining the structural ensemble of intrinsically disordered peptides (IDPs) facilitated by experimental measurements using circular dichroism spectroscopy (CD). A major challenge surrounding this approach stems from the deconvolution of experimental CD spectra into secondary structure features of the IDP ensemble. Currently available algorithms for CD deconvolution were designed to analyze the spectra of proteins with stable secondary structures. Herein, our work aims to minimize any bias from the peptide deconvolution analysis by implementing a non-negative linear least-squares fitting algorithm in conjunction with a CD reference data set that contains soluble and denatured proteins (SDP48). The non-negative linear least-squares method yields the best results for deconvolution of proteins with higher disordered content than currently available methods, according to a validation analysis of a set of protein spectra with Protein Data Bank entries. We subsequently used this analysis to deconvolute our experimental CD data to refine our computational model of the peptide secondary structure ensemble produced by all-atom molecular dynamics simulations with implicit solvent. We applied this approach to determine the ensemble structures of a set of short IDPs, that mimic the calmodulin binding domain of calcium/calmodulin-dependent protein kinase II and its 1-amino-acid and 3-amino-acid mutants. Our study offers a, to our knowledge, novel way to solve the ensemble secondary structures of IDPs in solution, which is important to advance the understanding of their roles in regulating signaling pathways through the formation of complexes with multiple partners.
Collapse
Affiliation(s)
- Jacob C Ezerski
- Department of Physics, University of Houston, Houston, Texas
| | - Pengzhi Zhang
- Department of Physics, University of Houston, Houston, Texas
| | | | - M Neal Waxham
- Department of Neurobiology and Anatomy, University of Texas, Health Science Center at Houston, Houston, Texas
| | - Margaret S Cheung
- Department of Physics, University of Houston, Houston, Texas; Center for Theoretical Biological Physics, Rice University, Houston, Texas.
| |
Collapse
|
28
|
Zou R, Zhou Y, Wang Y, Kuang G, Ågren H, Wu J, Tu Y. Free Energy Profile and Kinetics of Coupled Folding and Binding of the Intrinsically Disordered Protein p53 with MDM2. J Chem Inf Model 2020; 60:1551-1558. [DOI: 10.1021/acs.jcim.9b00920] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Rongfeng Zou
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, 10691 Stockholm, Sweden
| | - Yang Zhou
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, 10691 Stockholm, Sweden
| | - Yong Wang
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
| | - Guanglin Kuang
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, 10691 Stockholm, Sweden
| | - Hans Ågren
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, 10691 Stockholm, Sweden
- College of Chemistry and Chemical Engineering, Henan University, 475004 Kaifeng, Henan, P. R. China
| | - Junchen Wu
- Key Laboratory for Advanced Materials & Institute of Fine Chemicals, School of Chemistry and Molecular Engineering, East China University of Science and Technology, 200237 Shanghai, China
| | - Yaoquan Tu
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, 10691 Stockholm, Sweden
| |
Collapse
|
29
|
Mutual population-shift driven antibody-peptide binding elucidated by molecular dynamics simulations. Sci Rep 2020; 10:1406. [PMID: 31996730 PMCID: PMC6989527 DOI: 10.1038/s41598-020-58320-z] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Accepted: 01/14/2020] [Indexed: 11/08/2022] Open
Abstract
Antibody based bio-molecular drugs are an exciting, new avenue of drug development as an alternative to the more traditional small chemical compounds. However, the binding mechanism and the effect on the conformational ensembles of a therapeutic antibody to its peptide or protein antigen have not yet been well studied. We have utilized dynamic docking and path sampling simulations based on all-atom molecular dynamics to study the binding mechanism between the antibody solanezumab and the peptide amyloid-β (Aβ). Our docking simulations reproduced the experimental structure and gave us representative binding pathways, from which we accurately estimated the binding free energy. Not only do our results show why solanezumab has an explicit preference to bind to the monomeric form of Aβ, but that upon binding, both molecules are stabilized towards a specific conformation, suggesting that their complex formation follows a novel, mutual population-shift model, where upon binding, both molecules impact the dynamics of their reciprocal one.
Collapse
|
30
|
Higo J, Kasahara K, Wada M, Dasgupta B, Kamiya N, Hayami T, Fukuda I, Fukunishi Y, Nakamura H. Free-energy landscape of molecular interactions between endothelin 1 and human endothelin type B receptor: fly-casting mechanism. Protein Eng Des Sel 2019; 32:297-308. [PMID: 31608410 PMCID: PMC7052515 DOI: 10.1093/protein/gzz029] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2019] [Revised: 06/28/2019] [Accepted: 07/08/2019] [Indexed: 01/05/2023] Open
Abstract
The free-energy landscape of interaction between a medium-sized peptide, endothelin 1 (ET1), and its receptor, human endothelin type B receptor (hETB), was computed using multidimensional virtual-system coupled molecular dynamics, which controls the system's motions by introducing multiple reaction coordinates. The hETB embedded in lipid bilayer was immersed in explicit solvent. All molecules were expressed as all-atom models. The resultant free-energy landscape had five ranges with decreasing ET1-hETB distance: completely dissociative, outside-gate, gate, binding pocket, and genuine-bound ranges. In the completely dissociative range, no ET1-hETB interaction appeared. In the outside-gate range, an ET1-hETB attractive interaction was the fly-casting mechanism. In the gate range, the ET1 orientational variety decreased rapidly. In the binding pocket range, ET1 was in a narrow pathway with a steep free-energy slope. In the genuine-bound range, ET1 was in a stable free-energy basin. A G-protein-coupled receptor (GPCR) might capture its ligand from a distant place.
Collapse
Affiliation(s)
- Junichi Higo
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Osaka, Suita 565-0871, Japan
| | - Kota Kasahara
- College of Life Sciences, Ritsumeikan University, 1-1-1 Noji-higashi, Shiga, Kusatsu 525-8577, Japan
| | - Mitsuhito Wada
- Technology Research Association for Next Generation Natural Products Chemistry, 2-3-26, Aomi, Tokyo, Koto-ku 135-0064, Japan
| | - Bhaskar Dasgupta
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Osaka, Suita 565-0871, Japan
| | - Narutoshi Kamiya
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Osaka, Suita 565-0871, Japan
| | - Tomonori Hayami
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Osaka, Suita 565-0871, Japan
| | - Ikuo Fukuda
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yoshifumi Fukunishi
- Molecular Profiling Research Center for Drug Discovery (molprof), National Institute of Advanced Industrial Science and Technology (AIST), 2-3-26, Aomi, Tokyo, Koto-ku 135-0064, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Osaka, Suita 565-0871, Japan
| |
Collapse
|
31
|
Hayami T, Higo J, Nakamura H, Kasahara K. Multidimensional virtual-system coupled canonical molecular dynamics to compute free-energy landscapes of peptide multimer assembly. J Comput Chem 2019; 40:2453-2463. [PMID: 31282023 DOI: 10.1002/jcc.26020] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2019] [Revised: 05/26/2019] [Accepted: 06/17/2019] [Indexed: 11/05/2022]
Abstract
An enhanced-sampling method termed multidimensional virtual-system coupled canonical molecular dynamics (mD-VcMD) method is developed. In many cases, generalized-ensemble methods realizing enhanced sampling, for example, adaptive umbrella sampling, apply an effective potential, which is derived from temporarily assumed canonical distribution as a function of one or more arbitrarily defined reaction coordinates. However, it is not straightforward to estimate the appropriate canonical distribution, especially for cases applying multiple reaction coordinates. The current method, mD-VcMD, does not rely on the form of the canonical distribution. Therefore, it is practically useful to explore a high-dimensional reaction-coordinate space. In this article, formulation of mD-VcMD and its evaluation with the simple molecular models consisting of three or four alanine peptides are presented. We confirmed that mD-VcMD efficiently searched 2D and 3D reaction-coordinate spaces defined as interpeptide distances. Direct comparisons with results of long-term canonical MD simulations revealed that mD-VcMD produces correct canonical ensembles. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Tomonori Hayami
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka, 565-0871, Japan
| | - Junichi Higo
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka, 565-0871, Japan
| | - Haruki Nakamura
- 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
| |
Collapse
|
32
|
Ono S, Naylor MR, Townsend CE, Okumura C, Okada O, Lokey RS. Conformation and Permeability: Cyclic Hexapeptide Diastereomers. J Chem Inf Model 2019; 59:2952-2963. [PMID: 31042375 PMCID: PMC7751304 DOI: 10.1021/acs.jcim.9b00217] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Conformational ensembles of eight cyclic hexapeptide diastereomers in explicit cyclohexane, chloroform, and water were analyzed by multicanonical molecular dynamics (McMD) simulations. Free-energy landscapes (FELs) for each compound and solvent were obtained from the molecular shapes and principal component analysis at T = 300 K; detailed analysis of the conformational ensembles and flexibility of the FELs revealed that permeable compounds have different structural profiles even for a single stereoisomeric change. The average solvent-accessible surface area (SASA) in cyclohexane showed excellent correlation with the cell permeability, whereas this correlation was weaker in chloroform. The average SASA in water correlated with the aqueous solubility. The average polar surface area did not correlate with cell permeability in these solvents. A possible strategy for designing permeable cyclic peptides from FELs obtained from McMD simulations is proposed.
Collapse
Affiliation(s)
- Satoshi Ono
- Modality Laboratories, Innovative Research Division,
Mitsubishi Tanabe Pharma Corporation, 1000, Kamoshida-cho, Aoba-ku, Yokohama,
Kanagawa 227-0033, Japan
| | - Matthew R. Naylor
- Department of Chemistry and Biochemistry, University
of California Santa Cruz, 1156 High Street, Santa Cruz, California 95064, United
States
| | - Chad E. Townsend
- Department of Chemistry and Biochemistry, University
of California Santa Cruz, 1156 High Street, Santa Cruz, California 95064, United
States
| | - Chieko Okumura
- Modality Laboratories, Innovative Research Division,
Mitsubishi Tanabe Pharma Corporation, 1000, Kamoshida-cho, Aoba-ku, Yokohama,
Kanagawa 227-0033, Japan
| | - Okimasa Okada
- Modality Laboratories, Innovative Research Division,
Mitsubishi Tanabe Pharma Corporation, 1000, Kamoshida-cho, Aoba-ku, Yokohama,
Kanagawa 227-0033, Japan
| | - R. Scott Lokey
- Department of Chemistry and Biochemistry, University
of California Santa Cruz, 1156 High Street, Santa Cruz, California 95064, United
States
| |
Collapse
|
33
|
Kasahara K, Terazawa H, Takahashi T, Higo J. Studies on Molecular Dynamics of Intrinsically Disordered Proteins and Their Fuzzy Complexes: A Mini-Review. Comput Struct Biotechnol J 2019; 17:712-720. [PMID: 31303975 PMCID: PMC6603302 DOI: 10.1016/j.csbj.2019.06.009] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2019] [Revised: 05/29/2019] [Accepted: 06/11/2019] [Indexed: 11/19/2022] Open
Abstract
The molecular dynamics (MD) method is a promising approach toward elucidating the molecular mechanisms of intrinsically disordered regions (IDRs) of proteins and their fuzzy complexes. This mini-review introduces recent studies that apply MD simulations to investigate the molecular recognition of IDRs. Firstly, methodological issues by which MD simulations treat IDRs, such as developing force fields, treating periodic boundary conditions, and enhanced sampling approaches, are discussed. Then, several examples of the applications of MD to investigate molecular interactions of IDRs in terms of the two kinds of complex formations; coupled-folding and binding and fuzzy complex. MD simulations provide insight into the molecular mechanisms of these binding processes by sampling conformational ensembles of flexible IDRs. In particular, we focused on all-atom explicit-solvent MD simulations except for studies of higher-order assembly of IDRs. Recent advances in MD methods, and computational power make it possible to dissect the molecular details of realistic molecular systems involving the dynamic behavior of IDRs.
Collapse
Affiliation(s)
- Kota Kasahara
- College of Life Sciences, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
- Corresponding author.
| | - Hiroki Terazawa
- Graduate School of Life Sciences, 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
| | - Junichi Higo
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima-minamimachi, Chuo-ku, Kobe 650-0047, Japan
| |
Collapse
|
34
|
Kondo HX, Kusaka A, Kitakawa CK, Onari J, Yamanaka S, Nakamura H, Takano Y. Hydrogen bond donors and acceptors are generally depolarized in α-helices as revealed by a molecular tailoring approach. J Comput Chem 2019; 40:2043-2052. [PMID: 31099907 PMCID: PMC6767508 DOI: 10.1002/jcc.25859] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2019] [Revised: 04/09/2019] [Accepted: 04/29/2019] [Indexed: 11/24/2022]
Abstract
Hydrogen‐bond (H‐bond) interaction energies in α‐helices of short alanine peptides were systematically examined by precise density functional theory calculations, followed by a molecular tailoring approach. The contribution of each H‐bond interaction in α‐helices was estimated in detail from the entire conformation energies, and the results were compared with those in the minimal H‐bond models, in which only H‐bond donors and acceptors exist with the capping methyl groups. The former interaction energies were always significantly weaker than the latter energies, when the same geometries of the H‐bond donors and acceptors were applied. The chemical origin of this phenomenon was investigated by analyzing the differences among the electronic structures of the local peptide backbones of the α‐helices and those of the minimal H‐bond models. Consequently, we found that the reduced H‐bond energy originated from the depolarizations of both the H‐bond donor and acceptor groups, due to the repulsive interactions with the neighboring polar peptide groups in the α‐helix backbone. The classical force fields provide similar H‐bond energies to those in the minimal H‐bond models, which ignore the current depolarization effect, and thus they overestimate the actual H‐bond energies in α‐helices. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Hiroko X Kondo
- School of Regional Innovation and Social Design Engineering, Faculty of Engineering, Kitami Institute of Technology, 165 Koen-cho, Kitami, Hokkaido, 090-8507, Japan
| | - Ayumi Kusaka
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan
| | - Colin K Kitakawa
- Graduate School of Science, Osaka University, 1-1 Machikaneyamacho, Toyonaka, Osaka, 560-0043, Japan
| | - Jinta Onari
- Graduate School of Science, Osaka University, 1-1 Machikaneyamacho, Toyonaka, Osaka, 560-0043, Japan
| | - Shusuke Yamanaka
- Graduate School of Science, Osaka University, 1-1 Machikaneyamacho, Toyonaka, Osaka, 560-0043, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan
| | - Yu Takano
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan.,Graduate School of Information Sciences, Hiroshima City University, 3-4-1 Ozuka-Higashi Asa-Minami-Ku, Hiroshima, 731-3194, Japan
| |
Collapse
|
35
|
Efficient construction of a diverse conformational library for amyloid-β as an intrinsically disordered protein. J Mol Graph Model 2019; 88:183-193. [DOI: 10.1016/j.jmgm.2019.01.009] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2018] [Revised: 01/11/2019] [Accepted: 01/16/2019] [Indexed: 12/16/2022]
|
36
|
Iida S, Kawabata T, Kasahara K, Nakamura H, Higo J. Multimodal Structural Distribution of the p53 C-Terminal Domain upon Binding to S100B via a Generalized Ensemble Method: From Disorder to Extradisorder. J Chem Theory Comput 2019; 15:2597-2607. [DOI: 10.1021/acs.jctc.8b01042] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Shinji Iida
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Takeshi Kawabata
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Kota Kasahara
- College of Life Sciences, Ritsumeikan University, Noji-higashi 1-1-1, Kusatsu, Shiga 525-8577, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Junichi Higo
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
37
|
Higo J, Kasahara K, Nakamura H. Multi-dimensional virtual system introduced to enhance canonical sampling. J Chem Phys 2018; 147:134102. [PMID: 28987097 DOI: 10.1063/1.4986129] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
When an important process of a molecular system occurs via a combination of two or more rare events, which occur almost independently to one another, computational sampling for the important process is difficult. Here, to sample such a process effectively, we developed a new method, named the "multi-dimensional Virtual-system coupled Monte Carlo (multi-dimensional-VcMC)" method, where the system interacts with a virtual system expressed by two or more virtual coordinates. Each virtual coordinate controls sampling along a reaction coordinate. By setting multiple reaction coordinates to be related to the corresponding rare events, sampling of the important process can be enhanced. An advantage of multi-dimensional-VcMC is its simplicity: Namely, the conformation moves widely in the multi-dimensional reaction coordinate space without knowledge of canonical distribution functions of the system. To examine the effectiveness of the algorithm, we introduced a toy model where two molecules (receptor and its ligand) bind and unbind to each other. The receptor has a deep binding pocket, to which the ligand enters for binding. Furthermore, a gate is set at the entrance of the pocket, and the gate is usually closed. Thus, the molecular binding takes place via the two events: ligand approach to the pocket and gate opening. In two-dimensional (2D)-VcMC, the two molecules exhibited repeated binding and unbinding, and an equilibrated distribution was obtained as expected. A conventional canonical simulation, which was 200 times longer than 2D-VcMC, failed in sampling the binding/unbinding effectively. The current method is applicable to various biological systems.
Collapse
Affiliation(s)
- Junichi Higo
- 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
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka 565-0871, Japan
| |
Collapse
|
38
|
Kasahara K, Shiina M, Higo J, Ogata K, Nakamura H. Phosphorylation of an intrinsically disordered region of Ets1 shifts a multi-modal interaction ensemble to an auto-inhibitory state. Nucleic Acids Res 2018; 46:2243-2251. [PMID: 29309620 PMCID: PMC5861456 DOI: 10.1093/nar/gkx1297] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2017] [Revised: 12/15/2017] [Accepted: 12/19/2017] [Indexed: 12/28/2022] Open
Abstract
Multi-modal interactions are frequently observed in intrinsically disordered regions (IDRs) of proteins upon binding to their partners. In many cases, post-translational modifications in IDRs are accompanied by coupled folding and binding. From both molecular simulations and biochemical experiments with mutational studies, we show that the IDR including a Ser rich region (SRR) of the transcription factor Ets1, just before the DNA-binding core domain, undergoes multi-modal interactions when the SRR is not phosphorylated. In the phosphorylated state, the SRR forms a few specific complex structures with the Ets1 core, covering the recognition helix in the core and drastically reducing the DNA binding affinities as the auto-inhibitory state. The binding kinetics of mutated Ets1 indicates that aromatic residues in the SRR can be substituted with other hydrophobic residues for the interactions with the Ets1 core.
Collapse
Affiliation(s)
- Kota Kasahara
- College of Life Sciences, Ritsumeikan University, Noji-higashi 1-1-1, Kusatsu, Shiga 525-8577, Japan
| | - Masaaki Shiina
- Graduate School of Medicine, Yokohama City University, Fuku-ura 3–9, Kanazawa-ku, Yokohama, Kanagawa 236-0004, Japan
| | - Junichi Higo
- Institute for Protein Research, Osaka University, Yamada-oka 3-2, Suita, Osaka 565-0871, Japan
| | - Kazuhiro Ogata
- Graduate School of Medicine, Yokohama City University, Fuku-ura 3–9, Kanazawa-ku, Yokohama, Kanagawa 236-0004, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, Yamada-oka 3-2, Suita, Osaka 565-0871, Japan
| |
Collapse
|
39
|
Hayami T, Kasahara K, Nakamura H, Higo J. Molecular dynamics coupled with a virtual system for effective conformational sampling. J Comput Chem 2018; 39:1291-1299. [PMID: 29464736 DOI: 10.1002/jcc.25196] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2017] [Revised: 02/01/2018] [Accepted: 02/02/2018] [Indexed: 11/12/2022]
Abstract
An enhanced conformational sampling method is proposed: virtual-system coupled canonical molecular dynamics (VcMD). Although VcMD enhances sampling along a reaction coordinate, this method is free from estimation of a canonical distribution function along the reaction coordinate. This method introduces a virtual system that does not necessarily obey a physical law. To enhance sampling the virtual system couples with a molecular system to be studied. Resultant snapshots produce a canonical ensemble. This method was applied to a system consisting of two short peptides in an explicit solvent. Conventional molecular dynamics simulation, which is ten times longer than VcMD, was performed along with adaptive umbrella sampling. Free-energy landscapes computed from the three simulations mutually converged well. The VcMD provided quicker association/dissociation motions of peptides than the conventional molecular dynamics did. The VcMD method is applicable to various complicated systems because of its methodological simplicity. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Tomonori Hayami
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan
| | - Kota Kasahara
- College of Life Sciences, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga, 525-8577, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan
| | - Junichi Higo
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan
| |
Collapse
|
40
|
Mascarenhas NM, Terse VL, Gosavi S. Intrinsic Disorder in a Well-Folded Globular Protein. J Phys Chem B 2018; 122:1876-1884. [PMID: 29304275 DOI: 10.1021/acs.jpcb.7b12546] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The folded structure of the heterodimeric sweet protein monellin mimics single-chain proteins with topology β1-α1-β2-β3-β4-β5 (chain A: β3-β4-β5; chain B: β1-α1-β2). Furthermore, like naturally occurring single-chain proteins of a similar size, monellin folds cooperatively with no detectable intermediates. However, the two monellin chains, A and B, are marginally structured in isolation and fold only upon binding to each other. Thus, monellin presents a unique opportunity to understand the design of intrinsically disordered proteins that fold upon binding. Here, we study the folding of a single-chain variant of monellin (scMn) using simulations of an all heavy-atom structure-based model. These simulations can explain mechanistic details derived from scMn experiments performed using several different structural probes. scMn folds cooperatively in our structure-based simulations, as is also seen in experiments. We find that structure formation near the transition-state ensemble of scMn is not uniformly distributed but is localized to a hairpin-like structure which contains one strand from each chain (β2, β3). Thus, the sequence and the underlying energetics of heterodimeric monellin promote the early formation of the interchain interface (β2-β3). By studying computational scMn mutants whose "interchain" interactions are deleted, we infer that this energy distribution allows the two protein chains to remain largely disordered when this interface is not folded. From these results, we suggest that cutting the protein backbone of a globular protein between residues which lie within its folding nucleus may be one way to construct two disordered fragments which fold upon binding.
Collapse
Affiliation(s)
| | - Vishram L Terse
- Simons Centre for the Study of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research , Bangalore 560065, India
| | - Shachi Gosavi
- Simons Centre for the Study of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research , Bangalore 560065, India
| |
Collapse
|
41
|
Bernetti M, Masetti M, Pietrucci F, Blackledge M, Jensen MR, Recanatini M, Mollica L, Cavalli A. Structural and Kinetic Characterization of the Intrinsically Disordered Protein SeV NTAIL through Enhanced Sampling Simulations. J Phys Chem B 2017; 121:9572-9582. [DOI: 10.1021/acs.jpcb.7b08925] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Mattia Bernetti
- Department
of Pharmacy and Biotechnology, Alma Mater Studiorum − Università di Bologna, Via Belmeloro 6, 40126, Bologna, Italy
- CompuNet, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genova, Italy
| | - Matteo Masetti
- Department
of Pharmacy and Biotechnology, Alma Mater Studiorum − Università di Bologna, Via Belmeloro 6, 40126, Bologna, Italy
| | - Fabio Pietrucci
- Institut
de Minéralogie,
de Physique des Matériaux et de Cosmochimie, Sorbonne Universités−Université
Pierre et Marie Curie Paris 6, CNRS UMR 7590, IRD UMR 206, Museum national d’Histoire naturelle, F-75005 Paris, France
| | - Martin Blackledge
- Protein
Dynamics and Flexibility by NMR Group, Institut de Biologie Structurale, 38044 Grenoble, France
| | - Malene Ringkjobing Jensen
- Protein
Dynamics and Flexibility by NMR Group, Institut de Biologie Structurale, 38044 Grenoble, France
| | - Maurizio Recanatini
- Department
of Pharmacy and Biotechnology, Alma Mater Studiorum − Università di Bologna, Via Belmeloro 6, 40126, Bologna, Italy
| | - Luca Mollica
- CompuNet, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genova, Italy
| | - Andrea Cavalli
- Department
of Pharmacy and Biotechnology, Alma Mater Studiorum − Università di Bologna, Via Belmeloro 6, 40126, Bologna, Italy
- CompuNet, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genova, Italy
| |
Collapse
|
42
|
Zhang P, Tripathi S, Trinh H, Cheung MS. Opposing Intermolecular Tuning of Ca 2+ Affinity for Calmodulin by Neurogranin and CaMKII Peptides. Biophys J 2017; 112:1105-1119. [PMID: 28355539 PMCID: PMC5374985 DOI: 10.1016/j.bpj.2017.01.020] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2016] [Revised: 12/27/2016] [Accepted: 01/23/2017] [Indexed: 12/03/2022] Open
Abstract
We investigated the impact of bound calmodulin (CaM)-target compound structure on the affinity of calcium (Ca2+) by integrating coarse-grained models and all-atomistic simulations with nonequilibrium physics. We focused on binding between CaM and two specific targets, Ca2+/CaM-dependent protein kinase II (CaMKII) and neurogranin (Ng), as they both regulate CaM-dependent Ca2+ signaling pathways in neurons. It was shown experimentally that Ca2+/CaM (holoCaM) binds to the CaMKII peptide with overwhelmingly higher affinity than Ca2+-free CaM (apoCaM); the binding of CaMKII peptide to CaM in return increases the Ca2+ affinity for CaM. However, this reciprocal relation was not observed in the Ng peptide (Ng13–49), which binds to apoCaM or holoCaM with binding affinities of the same order of magnitude. Unlike the holoCaM-CaMKII peptide, whose structure can be determined by crystallography, the structural description of the apoCaM-Ng13–49 is unknown due to low binding affinity, therefore we computationally generated an ensemble of apoCaM-Ng13–49 structures by matching the changes in the chemical shifts of CaM upon Ng13–49 binding from nuclear magnetic resonance experiments. Next, we computed the changes in Ca2+ affinity for CaM with and without binding targets in atomistic models using Jarzynski’s equality. We discovered the molecular underpinnings of lowered affinity of Ca2+ for CaM in the presence of Ng13–49 by showing that the N-terminal acidic region of Ng peptide pries open the β-sheet structure between the Ca2+ binding loops particularly at C-domain of CaM, enabling Ca2+ release. In contrast, CaMKII peptide increases Ca2+ affinity for the C-domain of CaM by stabilizing the two Ca2+ binding loops. We speculate that the distinctive structural difference in the bound complexes of apoCaM-Ng13–49 and holoCaM-CaMKII delineates the importance of CaM’s progressive mechanism of target binding on its Ca2+ binding affinities.
Collapse
Affiliation(s)
- Pengzhi Zhang
- Department of Physics, University of Houston, Houston, Texas
| | | | - Hoa Trinh
- Department of Physics, University of Houston, Houston, Texas
| | - Margaret S Cheung
- Department of Physics, University of Houston, Houston, Texas; Center for Theoretical Biological Physics, Rice University, Houston, Texas.
| |
Collapse
|
43
|
Okuda M, Higo J, Komatsu T, Konuma T, Sugase K, Nishimura Y. Dynamics of the Extended String-Like Interaction of TFIIE with the p62 Subunit of TFIIH. Biophys J 2017; 111:950-62. [PMID: 27602723 DOI: 10.1016/j.bpj.2016.07.042] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2016] [Revised: 07/27/2016] [Accepted: 07/28/2016] [Indexed: 01/09/2023] Open
Abstract
General transcription factor II E (TFIIE) contains an acid-rich region (residues 378-393) in its α-subunit, comprising 13 acidic and two hydrophobic (Phe387 and Val390) residues. Upon binding to the p62 subunit of TFIIH, the acidic region adopts an extended string-like structure on the basic groove of the pleckstrin homology domain (PHD) of p62, and inserts Phe387 and Val390 into two shallow pockets in the groove. Here, we have examined the dynamics of this interaction by NMR and molecular dynamics (MD) simulations. Although alanine substitution of Phe387 and/or Val390 greatly reduced binding to PHD, the binding mode of the mutants was similar to that of the wild-type, as judged by the chemical-shift changes of the PHD. NMR relaxation dispersion profiles of the interaction exhibited large amplitudes for residues in the C-terminal half-string in the acidic region (Phe387, Glu388, Val390, Ala391, and Asp392), indicating a two-site binding mode: one corresponding to the final complex structure, and one to an off-pathway minor complex. To probe the off-pathway complex structure, an atomically detailed free-energy landscape of the binding mode was computed by all-atom multicanonical MD. The most thermodynamically stable cluster corresponded to the final complex structure. One of the next stable clusters was the off-pathway structure cluster, showing the reversed orientation of the C-terminal half-string on the PHD groove, as compared with the final structure. MD calculations elucidated that the C-terminal half-acidic-string forms encounter complexes mainly around the positive groove region with nearly two different orientations of the string, parallel and antiparallel to the final structure. Interestingly, the most encountered complexes exhibit a parallel-like orientation, suggesting that the string has a tendency to bind around the groove in the proper orientation with the aid of Phe387 and/or Val390 to proceed smoothly to the final complex structure.
Collapse
Affiliation(s)
- Masahiko Okuda
- Graduate School of Medical Life Science, Yokohama City University, Tsurumi-ku, Yokohama, Kanagawa, Japan
| | - Junichi Higo
- Institute for Protein Research, Osaka University, Suita, Osaka, Japan
| | - Tadashi Komatsu
- Graduate School of Medical Life Science, Yokohama City University, Tsurumi-ku, Yokohama, Kanagawa, Japan
| | - Tsuyoshi Konuma
- Department of Structural and Chemical Biology, Icahn School of Medicine at Mount Sinai, New York, New York
| | - Kenji Sugase
- Department of Molecular Engineering, Kyoto University, Kyoto, Japan
| | - Yoshifumi Nishimura
- Graduate School of Medical Life Science, Yokohama City University, Tsurumi-ku, Yokohama, Kanagawa, Japan.
| |
Collapse
|
44
|
Warren C, Shechter D. Fly Fishing for Histones: Catch and Release by Histone Chaperone Intrinsically Disordered Regions and Acidic Stretches. J Mol Biol 2017; 429:2401-2426. [PMID: 28610839 DOI: 10.1016/j.jmb.2017.06.005] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2017] [Revised: 06/05/2017] [Accepted: 06/06/2017] [Indexed: 01/21/2023]
Abstract
Chromatin is the complex of eukaryotic DNA and proteins required for the efficient compaction of the nearly 2-meter-long human genome into a roughly 10-micron-diameter cell nucleus. The fundamental repeating unit of chromatin is the nucleosome: 147bp of DNA wrapped about an octamer of histone proteins. Nucleosomes are stable enough to organize the genome yet must be dynamically displaced and reassembled to allow access to the underlying DNA for transcription, replication, and DNA damage repair. Histone chaperones are a non-catalytic group of proteins that are central to the processes of nucleosome assembly and disassembly and thus the fluidity of the ever-changing chromatin landscape. Histone chaperones are responsible for binding the highly basic histone proteins, shielding them from non-specific interactions, facilitating their deposition onto DNA, and aiding in their eviction from DNA. Although most histone chaperones perform these common functions, recent structural studies of many different histone chaperones reveal that there are few commonalities in their folds. Importantly, sequence-based predictions show that histone chaperones are highly enriched in intrinsically disordered regions (IDRs) and acidic stretches. In this review, we focus on the molecular mechanisms underpinning histone binding, selectivity, and regulation of these highly dynamic protein regions. We highlight new evidence suggesting that IDRs are often critical for histone chaperone function and play key roles in chromatin assembly and disassembly pathways.
Collapse
Affiliation(s)
- Christopher Warren
- Department of Biochemistry, Albert Einstein College of Medicine, 1300 Morris Park Avenue, Bronx, NY 10461, USA
| | - David Shechter
- Department of Biochemistry, Albert Einstein College of Medicine, 1300 Morris Park Avenue, Bronx, NY 10461, USA.
| |
Collapse
|
45
|
Drobnak I, Gradišar H, Ljubetič A, Merljak E, Jerala R. Modulation of Coiled-Coil Dimer Stability through Surface Residues while Preserving Pairing Specificity. J Am Chem Soc 2017; 139:8229-8236. [PMID: 28553984 DOI: 10.1021/jacs.7b01690] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The coiled-coil dimer is a widespread protein structural motif and, due to its designability, represents an attractive building block for assembling modular nanostructures. The specificity of coiled-coil dimer pairing is mainly based on hydrophobic and electrostatic interactions between residues at positions a, d, e, and g of the heptad repeat. Binding affinity, on the other hand, can also be affected by surface residues that face away from the dimerization interface. Here we show how design of the local helical propensity of interacting peptides can be used to tune the stabilities of coiled-coil dimers over a wide range. By designing intramolecular charge pairs, regions of high local helical propensity can be engineered to form trigger sequences, and dimer stability is adjusted without changing the peptide length or any of the directly interacting residues. This general principle is demonstrated by a change in thermal stability by more than 30 °C as a result of only two mutations outside the binding interface. The same approach was successfully used to modulate the stabilities in an orthogonal set of coiled-coils without affecting their binding preferences. The stability effects of local helical propensity and peptide charge are well described by a simple linear model, which should help improve current coiled-coil stability prediction algorithms. Our findings enable tuning the stabilities of coiled-coil-based building modules match a diverse range of applications in synthetic biology and nanomaterials.
Collapse
Affiliation(s)
- Igor Drobnak
- Department of Synthetic Biology and Immunology, National Institute of Chemistry , Hajdrihova 19, SI-1000 Ljubljana, Slovenia
| | - Helena Gradišar
- Department of Synthetic Biology and Immunology, National Institute of Chemistry , Hajdrihova 19, SI-1000 Ljubljana, Slovenia.,EN-FIST Centre of Excellence , Trg OF 13, SI-1000 Ljubljana, Slovenia
| | - Ajasja Ljubetič
- Department of Synthetic Biology and Immunology, National Institute of Chemistry , Hajdrihova 19, SI-1000 Ljubljana, Slovenia
| | - Estera Merljak
- Department of Synthetic Biology and Immunology, National Institute of Chemistry , Hajdrihova 19, SI-1000 Ljubljana, Slovenia
| | - Roman Jerala
- Department of Synthetic Biology and Immunology, National Institute of Chemistry , Hajdrihova 19, SI-1000 Ljubljana, Slovenia.,EN-FIST Centre of Excellence , Trg OF 13, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
46
|
Bekker GJ, Kamiya N, Araki M, Fukuda I, Okuno Y, Nakamura H. Accurate Prediction of Complex Structure and Affinity for a Flexible Protein Receptor and Its Inhibitor. J Chem Theory Comput 2017; 13:2389-2399. [DOI: 10.1021/acs.jctc.6b01127] [Citation(s) in RCA: 36] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Gert-Jan Bekker
- Institute
for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Narutoshi Kamiya
- Institute
for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
- Advanced Institute for Computational Science, RIKEN, 7-1-26 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Graduate
School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima
Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Mitsugu Araki
- Advanced Institute for Computational Science, RIKEN, 7-1-26 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Graduate
School of Medicine, Kyoto University, 54 Shogoin-Kawaharacho, Sakyo-ku, Kyoto 606-8507, Japan
| | - Ikuo Fukuda
- Institute
for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Yasushi Okuno
- Advanced Institute for Computational Science, RIKEN, 7-1-26 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Graduate
School of Medicine, Kyoto University, 54 Shogoin-Kawaharacho, Sakyo-ku, Kyoto 606-8507, Japan
| | - Haruki Nakamura
- Institute
for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| |
Collapse
|
47
|
Enhanced conformational sampling to visualize a free-energy landscape of protein complex formation. Biochem J 2017; 473:1651-62. [PMID: 27288028 PMCID: PMC4901360 DOI: 10.1042/bcj20160053] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2016] [Accepted: 02/29/2016] [Indexed: 01/23/2023]
Abstract
We introduce various, recently developed, generalized ensemble methods, which are useful to sample various molecular configurations emerging in the process of protein-protein or protein-ligand binding. The methods introduced here are those that have been or will be applied to biomolecular binding, where the biomolecules are treated as flexible molecules expressed by an all-atom model in an explicit solvent. Sampling produces an ensemble of conformations (snapshots) that are thermodynamically probable at room temperature. Then, projection of those conformations to an abstract low-dimensional space generates a free-energy landscape. As an example, we show a landscape of homo-dimer formation of an endothelin-1-like molecule computed using a generalized ensemble method. The lowest free-energy cluster at room temperature coincided precisely with the experimentally determined complex structure. Two minor clusters were also found in the landscape, which were largely different from the native complex form. Although those clusters were isolated at room temperature, with rising temperature a pathway emerged linking the lowest and second-lowest free-energy clusters, and a further temperature increment connected all the clusters. This exemplifies that the generalized ensemble method is a powerful tool for computing the free-energy landscape, by which one can discuss the thermodynamic stability of clusters and the temperature dependence of the cluster networks.
Collapse
|
48
|
Olson MA. On the Helix Propensity in Generalized Born Solvent Descriptions of Modeling the Dark Proteome. Front Mol Biosci 2017; 4:3. [PMID: 28197405 PMCID: PMC5281587 DOI: 10.3389/fmolb.2017.00003] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2016] [Accepted: 01/12/2017] [Indexed: 01/06/2023] Open
Abstract
Intrinsically disordered proteins that populate the so-called “Dark Proteome” offer challenging benchmarks of atomistic simulation methods to accurately model conformational transitions on a multidimensional energy landscape. This work explores the application of parallel tempering with implicit solvent models as a computational framework to capture the conformational ensemble of an intrinsically disordered peptide derived from the Ebola virus protein VP35. A recent X-ray crystallographic study reported a protein-peptide interface where the VP35 peptide underwent a folding transition from a disordered form to a helix-β-turn-helix topological fold upon molecular association with the Ebola protein NP. An assessment is provided of the accuracy of two generalized Born solvent models (GBMV2 and GBSW2) using the CHARMM force field and applied with temperature-based replica exchange dynamics to calculate the disorder propensity of the peptide and its probability density of states in a continuum solvent. A further comparison is presented of applying an explicit/implicit solvent hybrid replica exchange simulation of the peptide to determine the effect of modeling water interactions at the all-atom resolution.
Collapse
Affiliation(s)
- Mark A Olson
- Department of Cell Biology and Biochemistry, Molecular and Translational Sciences Division, United States Army Medical Research Institute for Infectious Diseases Fredrick, MD, USA
| |
Collapse
|
49
|
Higo J, Kasahara K, Dasgupta B, Nakamura H. Enhancement of canonical sampling by virtual-state transitions. J Chem Phys 2017; 146:044104. [DOI: 10.1063/1.4974087] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Junichi Higo
- 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
| | - Bhaskar Dasgupta
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka 565-0871,
Japan
- Technology Research Association for Next Generation Natural Products Chemistry, 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
|
50
|
Iida S, Mashimo T, Kurosawa T, Hojo H, Muta H, Goto Y, Fukunishi Y, Nakamura H, Higo J. Variation of free-energy landscape of the p53 C-terminal domain induced by acetylation: Enhanced conformational sampling. J Comput Chem 2016; 37:2687-2700. [PMID: 27735058 PMCID: PMC5242334 DOI: 10.1002/jcc.24494] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2016] [Revised: 07/26/2016] [Accepted: 08/31/2016] [Indexed: 12/29/2022]
Abstract
The C-terminal domain (CTD) of tumor suppressor protein p53 is an intrinsically disordered region that binds to various partner proteins, where lysine of CTD is acetylated/nonacetylated and histidine neutralized/non-neutralized. Because of the flexibility of the unbound CTD, a free-energy landscape (FEL) is a useful quantity for determining its statistical properties. We conducted enhanced conformational sampling of CTD in the unbound state via virtual system coupled multicanonical molecular dynamics, in which the lysine was acetylated or nonacetylated and histidine was charged or neutralized. The fragments were expressed by an all-atom model and were immersed in an explicit solvent. The acetylation and charge-neutralization varied FEL greatly, which might be convenient to exert a hub property. The acetylation slightly enhanced alpha-helix structures that are more compact than sheet/loop conformations. The charge-neutralization produced hairpins. Additionally, circular dichroism experiments confirmed the computational results. We propose possible binding mechanisms of CTD to partners by investigating FEL. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Shinji Iida
- Institute for Protein Research Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan
| | - Tadaaki Mashimo
- Technology Research Association for Next Generation Natural Products Chemistry, 2-3-26 Aomi, Koto-Ku, Tokyo, 135-0064, Japan
- IMSBIO Co, Ltd, Owl Tower 6F, 4-21-1, Higashi-ikebukuro, Toshima-ku, Tokyo, 170-0013, Japan
| | - Takashi Kurosawa
- Technology Research Association for Next Generation Natural Products Chemistry, 2-3-26 Aomi, Koto-Ku, Tokyo, 135-0064, Japan
- Hitachi Solutions East Japan, 21-1 Ekimaehoncho, Kawasaki-ku, Kanagawa, 210-0007, Japan
| | - Hironobu Hojo
- Institute for Protein Research Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan
| | - Hiroya Muta
- Institute for Protein Research Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan
| | - Yuji Goto
- Institute for Protein Research Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan
| | - Yoshifumi Fukunishi
- Technology Research Association for Next Generation Natural Products Chemistry, 2-3-26 Aomi, Koto-Ku, Tokyo, 135-0064, Japan
- Molecular Profiling Research Center for Drug Discovery (molprof), National Institute of Advanced Industrial Science and Technology (AIST), 2-3-36, Aomi, Koto-ku, Tokyo, 135-0064, Japan
| | - Haruki Nakamura
- Institute for Protein Research Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan
| | - Junichi Higo
- Institute for Protein Research Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan.
| |
Collapse
|