1
|
Das M, Gangopadhyay D, Hudecová J, Kessler J, Kapitán J, Bouř P. Monitoring Conformation and Protonation States of Glutathione by Raman Optical Activity and Molecular Dynamics. Chempluschem 2023; 88:e202300219. [PMID: 37283530 DOI: 10.1002/cplu.202300219] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Revised: 06/07/2023] [Accepted: 06/07/2023] [Indexed: 06/08/2023]
Abstract
Glutathione (GSH) is a common antioxidant and its biological activity depends on the conformation and protonation state. We used molecular dynamics, Raman and Raman optical activity (ROA) spectroscopies to investigate GSH structural changes in a broad pH range. Factor analysis of the spectra provided protonation constants (2.05, 3.45, 8.62, 9.41) in good agreement with previously published values. Following the analysis, spectra of differently protonated forms were obtained by extrapolation. The complete deprotonation of the thiol group above pH 11 was clearly visible in the spectra; however, many spectral features did not change much with pH. Experimental spectra at various pH values were decomposed into the simulated ones, which allowed us to study the conformer populations and quality of molecular dynamics (MD). According to this combined ROA/MD analysis conformation of the GSH backbone is affected by the pH changes only in a limited way. The combination of ROA with the computations thus has the potential to improve the MD force field and obtain more accurate populations of the conformer species. The methodology can be used for any molecule, but for a more detailed insight better computational techniques are needed in the future.
Collapse
Affiliation(s)
- Moumita Das
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 16610, Prague 6, Czech Republic
- Department of Analytical Chemistry Faculty of Chemical Engineering, University of Chemistry and Technology, Technická 5, 16628, Prague 6, Czech Republic
| | - Debraj Gangopadhyay
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 16610, Prague 6, Czech Republic
| | - Jana Hudecová
- Department of Optics, Palacký University, 17. listopadu 12, 77146, Olomouc, Czech Republic
| | - Jiří Kessler
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 16610, Prague 6, Czech Republic
| | - Josef Kapitán
- Department of Optics, Palacký University, 17. listopadu 12, 77146, Olomouc, Czech Republic
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 16610, Prague 6, Czech Republic
- Department of Analytical Chemistry Faculty of Chemical Engineering, University of Chemistry and Technology, Technická 5, 16628, Prague 6, Czech Republic
| |
Collapse
|
2
|
Kunche L, Natarajan U. Conformations and Solvation of Synthetic Polymers in Water by Generalized Born Implicit-Solvent Molecular Dynamics Simulations: Stereoisomers of Poly(acrylic acid) and Poly(methacrylic acid). J Phys Chem B 2023; 127:1244-1253. [PMID: 36705523 DOI: 10.1021/acs.jpcb.2c06658] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
We present the GB-OBC model as an approach for implicit-solvent MD simulations of a synthetic macromolecule in water. The model is tested and found to be successful in reproducing the chain dimensions and predicting the free energy of solvation of carboxylic acid vinyl polymers. The influence of stereochemistry and the hydrophobic nature of the polymer was investigated as a function of chain length (20 < N < 600) for poly(acrylic acid) (PAA) and poly(methacrylic acid) (PMA). The dimensionless parameters of the GB-OBC model were parameterized to be applicable to PAA and PMA. Scaling relations for chain dimensions obtained using implicit-solvent MD simulations in this study are in good agreement with those from experiments, theory of solvated chains in good solvents, and all-atom MD simulations in explicit water. Results show that ⟨Rg2⟩/NL2 is greater for the atactic chain as compared to the isotactic chain, for PAA as well as PMA. ⟨Rg2⟩/NL2 values of chains attain constancy in water for N = 200, with the values being greater for PMA. The PMA chain is conformationally more perturbed than the PAA chain, for both isotactic and atactic stereochemistry. The solvation free energy ΔGhyd of PAA and PMA in water is negative for all chain lengths (N = 20-600) and becomes more favorable with an increase in molecular weight. The ΔGhyd values for isotactic and atactic chains are identical at lower values of N but differ slightly for N > 300. Irrespective of the hydrophobic nature of the polymer, the atactic chain is thermodynamically more soluble in water as compared to the isotactic chain. The isotactic chain is less hydrophilic as compared to the atactic chain due to the closer proximity of the COOH groups along the backbone. This implicit solvent method is an effective way to accurately simulate the configurational properties and solvation of synthetic polymers in water.
Collapse
Affiliation(s)
- Lakshmikumar Kunche
- Macromolecular Modeling and Simulation Lab, Department of Chemical Engineering, Indian Institute of Technology (IIT) Madras, Chennai600036, India
| | - Upendra Natarajan
- Macromolecular Modeling and Simulation Lab, Department of Chemical Engineering, Indian Institute of Technology (IIT) Madras, Chennai600036, India
| |
Collapse
|
3
|
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]
|
4
|
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
|
5
|
Bekker GJ, Kamiya N. Thermal Stability Estimation of Single Domain Antibodies Using Molecular Dynamics Simulations. Methods Mol Biol 2023; 2552:151-163. [PMID: 36346591 DOI: 10.1007/978-1-0716-2609-2_8] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
In this chapter, we describe a protocol to estimate the thermal stability of single domain antibodies (sdAbs) using molecular dynamics (MD) simulations. This method measures the Q-value, the fraction of the native contacts, along the trajectory of high-temperature MD simulations starting from the experimental X-ray structure. We show a good correlation between the Q-value and the experimental melting temperature (Tm) in seven sdAbs. Assessing the Q-value on a per-residue level enabled us to identify residues that contribute to the instability and thus demonstrate which residues could be mutated to improve the stability and have later been validated by experiments. Our protocol extends beyond the application on sdAbs, as it is also suitable for other proteins and to determine the interfacial stability between protein and ligand.
Collapse
Affiliation(s)
- Gert-Jan Bekker
- Institute for Protein Research, Osaka University, Osaka, Japan
| | - Narutoshi Kamiya
- Graduate School of Information Science, University of Hyogo, Kobe, Hyogo, Japan.
| |
Collapse
|
6
|
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
|
7
|
Yamamoto S, Ishiro S, Kessler J, Bouř P. Intense chiral signal from α-helical poly-L-alanine observed in low-frequency Raman optical activity. Phys Chem Chem Phys 2021; 23:26501-26509. [PMID: 34806737 DOI: 10.1039/d1cp04401j] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Raman optical activity (ROA) spectral features reliably indicate the structure of peptides and proteins, but the signal is often weak. However, we observed significantly enhanced low-frequency bands for α-helical poly-L-alanine (PLA) in solution. The biggest ROA signal at ∼100 cm-1 is about 10 times stronger than higher-frequency bands described previously, which facilitates the detection. The low-frequency bands of PLA were compared to those of α-helical proteins. For PLA, density functional simulations well reproduced the experimental spectra and revealed that about 12 alanine residues within two turns of the α-helix generate the strong ROA band. Averaging based on molecular dynamics (MD) provided an even more realistic spectrum compared to the static model. The low-frequency bands could be largely related to a collective motion of the α-helical backbone, partially modulated by the solvent. Helical and intermolecular vibrational coordinates have been introduced and the helical unwinding modes were assigned to the strongest ROA signal at 101-128 cm-1. Further analysis indicated that the helically arranged amide and methyl groups are important for the strong chiral signal of PLA, while the local chiral centers CαH contribute in a minor way only. The strong low-frequency ROA can thus provide precious information about the motions of the peptide backbone and facilitate future protein studies.
Collapse
Affiliation(s)
- Shigeki Yamamoto
- Department of Chemistry, Graduate School of Science, Osaka University, Osaka 560-0043, Japan.
| | - Shota Ishiro
- Department of Chemistry, Graduate School of Science, Osaka University, Osaka 560-0043, Japan.
| | - Jiří Kessler
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo náměstí 2, 16610 Prague, Czech Republic.
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo náměstí 2, 16610 Prague, Czech Republic.
| |
Collapse
|
8
|
Ikebe J, Suzuki M, Komori A, Kobayashi K, Kameda T. Enzyme modification using mutation site prediction method for enhancing the regioselectivity of substrate reaction sites. Sci Rep 2021; 11:19004. [PMID: 34602611 PMCID: PMC8488038 DOI: 10.1038/s41598-021-98433-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2021] [Accepted: 09/02/2021] [Indexed: 11/16/2022] Open
Abstract
Enzymes with low regioselectivity of substrate reaction sites may produce multiple products from a single substrate. When a target product is produced industrially using these enzymes, the production of non-target products (byproducts) causes adverse effects such as increased processing costs for purification and the amount of raw material. Thus it is required the development of modified enzymes to reduce the amount of byproducts’ production. In this paper, we report a method called mutation site prediction for enhancing the regioselectivity of substrate reaction sites (MSPER). MSPER takes conformational data for docking poses of an enzyme and a substrate as input and automatically generates a ranked list of mutation sites to destabilize docking poses for byproducts while maintaining those for target products in silico. We applied MSPER to the enzyme cytochrome P450 CYP102A1 (BM3) and the two substrates to enhance the regioselectivity for four target products with different reaction sites. The 13 of the total 14 top-ranked mutation sites predicted by MSPER for the four target products succeeded in selectively enhancing the regioselectivity up to 6.4-fold. The results indicate that MSPER can distinguish differences of substrate structures and the reaction sites, and can accurately predict mutation sites to enhance regioselectivity without selection by directed evolution screening.
Collapse
Affiliation(s)
- Jinzen Ikebe
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-ku, Tokyo, 135-0064, Japan
| | - Munenori Suzuki
- KNC Bio-Research Center, KNC Laboratories Co., Ltd., 1-1-1 Murotani, Nishi-ku, Kobe, Hyogo, 651-2241, Japan
| | - Aya Komori
- KNC Bio-Research Center, KNC Laboratories Co., Ltd., 1-1-1 Murotani, Nishi-ku, Kobe, Hyogo, 651-2241, Japan
| | - Kaito Kobayashi
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-ku, Tokyo, 135-0064, Japan
| | - Tomoshi Kameda
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-ku, Tokyo, 135-0064, Japan.
| |
Collapse
|
9
|
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
|
10
|
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
|
11
|
Higo J, Kawabata T, Kusaka A, Kasahara K, Kamiya N, Fukuda I, Mori K, Hata Y, Fukunishi Y, Nakamura H. Molecular Interaction Mechanism of a 14-3-3 Protein with a Phosphorylated Peptide Elucidated by Enhanced Conformational Sampling. J Chem Inf Model 2020; 60:4867-4880. [PMID: 32910853 DOI: 10.1021/acs.jcim.0c00551] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Enhanced conformational sampling, a genetic-algorithm-guided multidimensional virtual-system coupled molecular dynamics, can provide equilibrated conformational distributions of a receptor protein and a flexible ligand at room temperature. The distributions provide not only the most stable but also semistable complex structures and propose a ligand-receptor binding process. This method was applied to a system consisting of a receptor protein, 14-3-3ε, and a flexible peptide, phosphorylated myeloid leukemia factor 1 (pMLF1). The results present comprehensive binding pathways of pMLF1 to 14-3-3ε. We identified four thermodynamically stable clusters of MLF1 on the 14-3-3ε surface and free-energy barriers among some clusters. The most stable cluster includes two high-density spots connected by a narrow corridor. When pMLF1 passes the corridor, a salt-bridge relay (switching) related to the phosphorylated residue of pMLF1 occurs. Conformations in one high-density spot are similar to the experimentally determined complex structure. Three-dimensional distributions of residues in the intermolecular interface rationally explain the binding constant changes resulting from the alanine mutation experiment for the residues. We also performed a simulation of nonphosphorylated peptide and 14-3-3ε, which demonstrated that the complex structure was unstable, suggesting that phosphorylation of the peptide is crucially important for binding to 14-3-3ε.
Collapse
Affiliation(s)
- Junichi Higo
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Takeshi Kawabata
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka 565-0871, Japan
| | - Ayumi Kusaka
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka 565-0871, Japan
| | - Kota Kasahara
- College of Life Sciences, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
| | - Narutoshi Kamiya
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Ikuo Fukuda
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Kentaro Mori
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,Maizuru College, 234 Shiroya, Maizuru, Kyoto 625-8511 Japan
| | - Yutaka Hata
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima Minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yoshifumi Fukunishi
- Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 2-3-26, Aomi, Koto-ku, Tokyo 135-0064, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamada-oka, Suita, Osaka 565-0871, Japan
| |
Collapse
|
12
|
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
|
13
|
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] [MESH Headings] [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 EngineeringFaculty of Engineering, Kitami Institute of Technology, 165 Koen‐choKitamiHokkaido090‐8507Japan
| | - Ayumi Kusaka
- Institute for Protein ResearchOsaka University, 3‐2 YamadaokaSuitaOsaka565‐0871Japan
| | - Colin K. Kitakawa
- Graduate School of ScienceOsaka University, 1‐1 MachikaneyamachoToyonakaOsaka560‐0043Japan
| | - Jinta Onari
- Graduate School of ScienceOsaka University, 1‐1 MachikaneyamachoToyonakaOsaka560‐0043Japan
| | - Shusuke Yamanaka
- Graduate School of ScienceOsaka University, 1‐1 MachikaneyamachoToyonakaOsaka560‐0043Japan
| | - Haruki Nakamura
- Institute for Protein ResearchOsaka University, 3‐2 YamadaokaSuitaOsaka565‐0871Japan
| | - Yu Takano
- Institute for Protein ResearchOsaka University, 3‐2 YamadaokaSuitaOsaka565‐0871Japan
- Graduate School of Information SciencesHiroshima City University, 3‐4‐1 Ozuka‐Higashi Asa‐Minami‐KuHiroshima731‐3194Japan
| |
Collapse
|
14
|
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
|
15
|
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
|
16
|
Kasahara K, Sakuraba S, Fukuda I. Enhanced Sampling of Molecular Dynamics Simulations of a Polyalanine Octapeptide: Effects of the Periodic Boundary Conditions on Peptide Conformation. J Phys Chem B 2018; 122:2495-2503. [PMID: 29439570 DOI: 10.1021/acs.jpcb.7b10830] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
We investigate the problem of artifacts caused by the periodic boundary conditions (PBC) used in molecular simulation studies. Despite the long history of simulations with PBCs, the existence of measurable artifacts originating from PBCs applied to inherently nonperiodic physical systems remains controversial. Specifically, these artifacts appear as differences between simulations of the same system but with different simulation-cell sizes. Earlier studies have implied that, even in the simple case of a small model peptide in water, sampling inefficiency is a major obstacle to understanding these artifacts. In this study, we have resolved the sampling issue using the replica exchange molecular dynamics (REMD) enhanced-sampling method to explore PBC artifacts. Explicitly solvated zwitterionic polyalanine octapeptides with three different cubic-cells, having dimensions of L = 30, 40, and 50 Å, were investigated to elucidate the differences with 64 replica × 500 ns REMD simulations using the AMBER parm99SB force field. The differences among them were not large overall, and the results for the L = 30 and 40 Å simulations in the conformational free energy landscape were found to be very similar at room temperature. However, a small but statistically significant difference was seen for L = 50 Å. We observed that extended conformations were slightly overstabilized in the smaller systems. The origin of these artifacts is discussed by comparison to an electrostatic calculation method without PBCs.
Collapse
Affiliation(s)
- Kota Kasahara
- College of Life Sciences , Ritsumeikan University , 1-1-1 Noji-higashi , Kusatsu , Shiga 525-8577 , Japan
| | - Shun Sakuraba
- Graduate School of Frontier Sciences , The University of Tokyo , 5-1-5 Kashiwanoha , Kashiwa , Chiba 277-8561 , Japan
| | - Ikuo Fukuda
- Institute for Protein Research , Osaka University , 3-2 Yamada-oka , Suita , Osaka 565-0871 , Japan
| |
Collapse
|
17
|
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
|
18
|
Accessibility of the histone H3 tail in the nucleosome for binding of paired readers. Nat Commun 2017; 8:1489. [PMID: 29138400 PMCID: PMC5686127 DOI: 10.1038/s41467-017-01598-x] [Citation(s) in RCA: 54] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2017] [Accepted: 10/03/2017] [Indexed: 12/03/2022] Open
Abstract
Combinatorial polyvalent contacts of histone-binding domains or readers commonly mediate localization and activities of chromatin-associated proteins. A pair of readers, the PHD fingers of the protein CHD4, has been shown to bivalently recognize histone H3 tails. Here we describe a mechanism by which these linked but independent readers bind to the intact nucleosome core particle (NCP). Comprehensive NMR, chemical reactivity, molecular dynamics, and fluorescence analyses point to the critical roles of intra-nucleosomal histone-DNA interactions that reduce the accessibility of H3 tails in NCP, the nucleosomal DNA, and the linker between readers in modulating nucleosome- and/or histone-binding activities of the readers. We show that the second PHD finger of CHD4 initiates recruitment to the nucleosome, however both PHDs are required to alter the NCP dynamics. Our findings reveal a distinctive regulatory mechanism for the association of paired readers with the nucleosome that provides an intricate balance between cooperative and individual activities of the readers. The chromatin remodeller CHD4 contains two PHD finger reader domains that have been shown to bivalently recognize H3 histone tails. Here, the authors describe a mechanism by which the PHD fingers bind to the intact nucleosome core particle, revealing both cooperative and individual interactions.
Collapse
|
19
|
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
|
20
|
Harada R, Takano Y, Shigeta Y. Common folding processes of mini-proteins: Partial formations of secondary structures initiate the immediate protein folding. J Comput Chem 2017; 38:790-797. [DOI: 10.1002/jcc.24748] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2016] [Revised: 12/20/2016] [Accepted: 01/12/2017] [Indexed: 11/10/2022]
Affiliation(s)
- Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba; 1-1-1 Tennodai Tsukuba Ibaraki 305-8577 Japan
| | - Yu Takano
- Department of Biomedical Information Sciences; Graduate School of Information Sciences, Hiroshima City University; 3-4-1 Ozuka-Higashi, Asa-Minami-Ku Hiroshima 731-3194 Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba; 1-1-1 Tennodai Tsukuba Ibaraki 305-8577 Japan
| |
Collapse
|
21
|
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
|
22
|
Dasgupta B, Nakamura H, Higo J. Flexible binding simulation by a novel and improved version of virtual-system coupled adaptive umbrella sampling. Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2016.09.059] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
23
|
Kamiya N, Mashimo T, Takano Y, Kon T, Kurisu G, Nakamura H. Elastic properties of dynein motor domain obtained from all-atom molecular dynamics simulations. Protein Eng Des Sel 2016; 29:317-325. [PMID: 27334455 PMCID: PMC4955872 DOI: 10.1093/protein/gzw022] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2016] [Accepted: 05/09/2016] [Indexed: 01/09/2023] Open
Abstract
Dyneins are large microtubule motor proteins that convert ATP energy to mechanical power. High-resolution crystal structures of ADP-bound cytoplasmic dynein have revealed the organization of the motor domain, comprising the AAA(+) ring, the linker, the stalk/strut and the C sequence. Recently, the ADP.vanadate-bound structure, which is similar to the ATP hydrolysis transition state, revealed how the structure of dynein changes upon ATP binding. Although both the ADP- and ATP-bound state structures have been resolved, the dynamic properties at the atomic level remain unclear. In this work, we built two models named 'the ADP model' and 'the ATP model', where ADP and ATP are bound to AAA1 in the AAA(+) ring, respectively, to observe the initial procedure of the structural change from the unprimed to the primed state. We performed 200-ns molecular dynamics simulations for both models and compared their structures and dynamics. The motions of the stalk, consisting of a long coiled coil with a microtubule-binding domain, significantly differed between the two models. The elastic properties of the stalk were analyzed and compared with the experimental results.
Collapse
Affiliation(s)
- Narutoshi Kamiya
- Advanced Institute for Computational Science, RIKEN, QBiC Building B, 6-2-4, Furuedai, Suita, Osaka 565-0874, Japan
| | - Tadaaki Mashimo
- Technology Research Association for Next Generation Natural Products Chemistry (N2PC), 2-3-26, Aomi, Koto-ku, Tokyo 135-0064, Japan.,IMSBIO Co. Ltd, Owl Tower, 4-21-1, Higashi-Ikebukuro, Toshima-ku, Tokyo 170-0013, Japan
| | - Yu Takano
- Department of Biomedical Information Sciences, Graduate School of Information Sciences, Hiroshima City University, 3-4-1, Ozuka-Higashi, Asa-Minami-ku, Hiroshima 731-3194, Japan
| | - Takahide Kon
- Department of Biological Sciences, Graduate School of Science, and Faculty of Science, Osaka University, 1-1, Machikaneyama-cho, Toyonaka, Osaka 560-0043, Japan
| | - Genji Kurisu
- Institute for Protein Research, Osaka University, 3-2, Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2, Yamadaoka, Suita, Osaka 565-0871, Japan
| |
Collapse
|
24
|
H3 Histone Tail Conformation within the Nucleosome and the Impact of K14 Acetylation Studied Using Enhanced Sampling Simulation. PLoS Comput Biol 2016; 12:e1004788. [PMID: 26967163 PMCID: PMC4788430 DOI: 10.1371/journal.pcbi.1004788] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2015] [Accepted: 02/03/2016] [Indexed: 11/19/2022] Open
Abstract
Acetylation of lysine residues in histone tails is associated with gene transcription. Because histone tails are structurally flexible and intrinsically disordered, it is difficult to experimentally determine the tail conformations and the impact of acetylation. In this work, we performed simulations to sample H3 tail conformations with and without acetylation. The results show that irrespective of the presence or absence of the acetylation, the H3 tail remains in contact with the DNA and assumes an α-helix structure in some regions. Acetylation slightly weakened the interaction between the tail and DNA and enhanced α-helix formation, resulting in a more compact tail conformation. We inferred that this compaction induces unwrapping and exposure of the linker DNA, enabling DNA-binding proteins (e.g., transcription factors) to bind to their target sequences. In addition, our simulation also showed that acetylated lysine was more often exposed to the solvent, which is consistent with the fact that acetylation functions as a post-translational modification recognition site marker. Post-translational modification (PTM) of histone tails is an important component of epigenetics. Acetylation of histone tails generally functions to activate gene expression, though the molecular mechanism is not well understood. We used enhanced sampling simulation to examine the impact of acetylation on the structure of the histone H3 tail within the nucleosome. The results suggest acetylation makes the H3 tail conformation more compact and enhances dissociation of nucleosomal DNA from the histone core. Further, the acetylated lysine was more exposed to the solvent, which is consistent with its role as a PTM recognition site marker. These findings increase our understanding of the impact of PTM on nucleosome stability and dynamics and on the higher order structure of chromatin.
Collapse
|
25
|
Ikebe J, Umezawa K, Higo J. Enhanced sampling simulations to construct free-energy landscape of protein-partner substrate interaction. Biophys Rev 2016; 8:45-62. [PMID: 28510144 PMCID: PMC5425738 DOI: 10.1007/s12551-015-0189-z] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2015] [Accepted: 12/07/2015] [Indexed: 01/08/2023] Open
Abstract
Molecular dynamics (MD) simulations using all-atom and explicit solvent models provide valuable information on the detailed behavior of protein-partner substrate binding at the atomic level. As the power of computational resources increase, MD simulations are being used more widely and easily. However, it is still difficult to investigate the thermodynamic properties of protein-partner substrate binding and protein folding with conventional MD simulations. Enhanced sampling methods have been developed to sample conformations that reflect equilibrium conditions in a more efficient manner than conventional MD simulations, thereby allowing the construction of accurate free-energy landscapes. In this review, we discuss these enhanced sampling methods using a series of case-by-case examples. In particular, we review enhanced sampling methods conforming to trivial trajectory parallelization, virtual-system coupled multicanonical MD, and adaptive lambda square dynamics. These methods have been recently developed based on the existing method of multicanonical MD simulation. Their applications are reviewed with an emphasis on describing their practical implementation. In our concluding remarks we explore extensions of the enhanced sampling methods that may allow for even more efficient sampling.
Collapse
Affiliation(s)
- Jinzen Ikebe
- Molecular Modeling and Simulation Group, Japan Atomic Energy Agency, 8-1-7 Umemidai, Kizugawa, Kyoto, 619-0215, Japan
| | - Koji Umezawa
- Department of Pure and Applied Physics, Waseda University, Okubo 3-4-1, Shinjuku-Ku, Tokyo, 169-8555, Japan
| | - Junichi Higo
- Institute for Protein Research, Osaka University, Suita, Osaka, 565-0871, Japan.
| |
Collapse
|
26
|
Terakawa T, Higo J, Takada S. Multi-scale ensemble modeling of modular proteins with intrinsically disordered linker regions: application to p53. Biophys J 2015; 107:721-729. [PMID: 25099811 DOI: 10.1016/j.bpj.2014.06.026] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2014] [Revised: 05/28/2014] [Accepted: 06/18/2014] [Indexed: 10/24/2022] Open
Abstract
In eukaryotic proteins, intrinsically disordered regions (IDRs) are ubiquitous and often exist in linker regions that flank the functional domains of modular proteins, regulating their functions. For detailed structural ensemble modeling of IDRs, we propose a multiscale method for IDRs that possess significant long-range order in modular proteins and apply it to the eukaryotic transcription factor p53 as an example. First, we performed all-atom (AA) molecular dynamics (MD) simulations of the explicitly solvated p53 linker region, without experimental restraint terms, finding fractional long-range contacts within the linker. Second, we fed this AA MD ensemble into a coarse-grained (CG) model, finding an optimal set of contact potentials. The optimized CG MD simulations reproduced the contact probability map from the AA MD simulations. Finally, we performed the CG MD simulation of the tetrameric p53 fragments including the core domains, the linker, and the tetramerization domain. Using the obtained ensemble, we theoretically calculated the small angle x-ray scattering (SAXS) profile of this fragment. The obtained SAXS profile agrees well with the experiment. We also found that the long-range contacts in the p53 linker region are required to reproduce the experimental SAXS profile. The developed framework in which we calculate the long-range contact probability map from the AA MD simulation and incorporate it to the CG model can be applied to broad range of IDRs.
Collapse
Affiliation(s)
- Tsuyoshi Terakawa
- Department of Biophysics, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwakecho, Sakyo, Kyoto, 606-8502, Japan
| | - Junichi Higo
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwakecho, Sakyo, Kyoto, 606-8502, Japan.
| |
Collapse
|
27
|
Higo J, Dasgupta B, Mashimo T, Kasahara K, Fukunishi Y, Nakamura H. Virtual-system-coupled adaptive umbrella sampling to compute free-energy landscape for flexible molecular docking. J Comput Chem 2015; 36:1489-501. [DOI: 10.1002/jcc.23948] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2015] [Revised: 04/21/2015] [Accepted: 04/24/2015] [Indexed: 01/20/2023]
Affiliation(s)
- Junichi Higo
- Institute for Protein Research, Osaka University; 3-2 Yamadaoka, Suita Osaka 565-0871 Japan
| | - Bhaskar Dasgupta
- 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
- Information, Mathematical Science and Bioinformatics Co., Ltd.; 4-21-1, Higashiikebukuro Toshima-ku Tokyo 170-0013 Japan
| | - Kota Kasahara
- Institute for Protein Research, Osaka University; 3-2 Yamadaoka, Suita Osaka 565-0871 Japan
| | - Yoshifumi Fukunishi
- Molecular Profiling Research Center for Drug Discovery (molprof), National Institute of Advanced Industrial Science and Technology (AIST); 2-3-26 Aomi Koto-ku Tokyo 135-0064 Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University; 3-2 Yamadaoka, Suita Osaka 565-0871 Japan
| |
Collapse
|
28
|
Štěpánek P, Bouř P. Multi-scale modeling of electronic spectra of three aromatic amino acids: importance of conformational averaging and explicit solute-solvent interactions. Phys Chem Chem Phys 2014; 16:20639-49. [PMID: 25158079 DOI: 10.1039/c4cp02668c] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Electronic transitions in the ultraviolet and visible spectral range can reveal a wealth of information about biomolecular geometry and interactions, such as those involved in protein folding. However, the modeling that provides the necessary link between spectral shapes and the structure is often difficult even for seemingly simple systems. To understand as to how conformational equilibria and solute-solvent interaction influence spectral intensities, we collected absorption (UV-vis), electronic circular dichroism (ECD), and magnetic circular dichroism (MCD) spectra of phenylalanine (Phe), tyrosine (Tyr) and tryptophan (Trp) zwitterions in aqueous solutions, and compared them with quantum-chemical simulations. These aromatic amino acids provide a relatively strong signal in the accessible wavelength range. At the same time, they allow for a relatively accurate modeling. Energies and intensities of spectral bands were reproduced by the time-dependent density functional theory (TD DFT). The solvent was approximated by a continuum as well as clusters containing solvent molecules from the first hydration sphere. The ECD signal was found to be strongly dependent on molecular conformation, and the dependence was much weaker in UV-vis and MCD spectra. All spectral intensities, however, were significantly affected by the solvent approximation; especially for ECD and MCD the usual polarizable continuum solvent model did not yield satisfactory spectral shapes. On the other hand, averaging of the clusters obtained from molecular dynamics simulations provided an unprecedented agreement with the experiment. Proper modeling of the interactions with the environment thus makes the information about the molecular structure, as obtained from the electronic spectra, more accurate and reliable.
Collapse
Affiliation(s)
- Petr Štěpánek
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo náměstí 2, 166 10 Prague, Czech Republic.
| | | |
Collapse
|
29
|
Nishikawa Y, Oyama T, Kamiya N, Kon T, Toyoshima YY, Nakamura H, Kurisu G. Structure of the entire stalk region of the Dynein motor domain. J Mol Biol 2014; 426:3232-3245. [PMID: 25058684 DOI: 10.1016/j.jmb.2014.06.023] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2014] [Revised: 06/12/2014] [Accepted: 06/24/2014] [Indexed: 01/10/2023]
Abstract
Dyneins are large microtubule-based motor complexes that power a range of cellular processes including the transport of organelles, as well as the beating of cilia and flagella. The motor domain is located within the dynein heavy chain and comprises an N-terminal mechanical linker element, a central ring of six AAA+ modules of which four bind or hydrolyze ATP, and a long stalk extending from the AAA+ring with a microtubule-binding domain (MTBD) at its tip. A crucial mechanism underlying the motile activity of cytoskeletal motor proteins is precise coupling between the ATPase and track-binding activities. In dynein, a stalk region consisting of a long (~15nm) antiparallel coiled coil separates these two activities, which must facilitate communication between them. This communication is mediated by a small degree of helix sliding in the coiled coil. However, no high-resolution structure is available of the entire stalk region including the MTBD. Here, we have reported the structure of the entire stalk region of mouse cytoplasmic dynein in a weak microtubule-binding state, which was determined using X-ray crystallography, and have compared it with the dynein motor domain from Dictyostelium discoideum in a strong microtubule-binding state and with a mouse MTBD with its distal portion of the coiled coil fused to seryl-tRNA synthetase from Thermus thermophilus. Our results strongly support the helix-sliding model based on the complete structure of the dynein stalk with a different form of coiled-coil packing. We also propose a plausible mechanism of helix sliding together with further analysis using molecular dynamics simulations. Our results present the importance of conserved proline residues for an elastic motion of stalk coiled coil and imply the manner of change between high-affinity state and low-affinity state of MTBD.
Collapse
Affiliation(s)
- Yosuke Nishikawa
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan; Department of Macromolecular Science, Graduate School of Science, Osaka University, Toyonaka, Osaka 560-0043, Japan
| | - Takuji Oyama
- 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
| | - Takahide Kon
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Yoko Y Toyoshima
- Department of Life Sciences, Graduate School of Arts and Sciences, University of Tokyo, Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Genji Kurisu
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan; Department of Macromolecular Science, Graduate School of Science, Osaka University, Toyonaka, Osaka 560-0043, Japan.
| |
Collapse
|
30
|
Kessler J, Keiderling TA, Bouř P. Arrangement of Fibril Side Chains Studied by Molecular Dynamics and Simulated Infrared and Vibrational Circular Dichroism Spectra. J Phys Chem B 2014; 118:6937-45. [DOI: 10.1021/jp502178d] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jiří Kessler
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 166
10 Prague, Czech Republic
- Department
of Physical and Macromolecular Chemistry, Faculty of Science, Charles University, Hlavova 8, 128 40 Prague, Czech Republic
| | - Timothy A. Keiderling
- Department
of Chemistry, University of Illinois at Chicago, 845 West Taylor
Street, Chicago, Illinois 60607-7061, United States
| | - Petr Bouř
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 166
10 Prague, Czech Republic
| |
Collapse
|
31
|
Kessler J, Bouř P. Molecular dynamics with helical periodic boundary conditions. J Comput Chem 2014; 35:1552-9. [DOI: 10.1002/jcc.23653] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2014] [Revised: 05/20/2014] [Accepted: 05/26/2014] [Indexed: 11/07/2022]
Affiliation(s)
- Jiří Kessler
- Institute of Organic Chemistry and Biochemistry; Academy of Sciences; Flemingovo náměstí 2 166 10 Prague Czech Republic
- Department of Physical and Macromolecular Chemistry; Faculty of Science, Charles University; Hlavova 8 128 40 Prague Czech Republic
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry; Academy of Sciences; Flemingovo náměstí 2 166 10 Prague Czech Republic
| |
Collapse
|
32
|
Free-Energy Landscape of Intrinsically Disordered Proteins Investigated by All-Atom Multicanonical Molecular Dynamics. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2014; 805:331-51. [DOI: 10.1007/978-3-319-02970-2_14] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
33
|
Ikebe J, Sakuraba S, Kono H. Adaptive lambda square dynamics simulation: an efficient conformational sampling method for biomolecules. J Comput Chem 2013; 35:39-50. [PMID: 24166005 DOI: 10.1002/jcc.23462] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2013] [Revised: 08/23/2013] [Accepted: 09/20/2013] [Indexed: 11/10/2022]
Abstract
A novel, efficient sampling method for biomolecules is proposed. The partial multicanonical molecular dynamics (McMD) was recently developed as a method that improved generalized ensemble (GE) methods to focus sampling only on a part of a system (GEPS); however, it was not tested well. We found that partial McMD did not work well for polylysine decapeptide and gave significantly worse sampling efficiency than a conventional GE. Herein, we elucidate the fundamental reason for this and propose a novel GEPS, adaptive lambda square dynamics (ALSD), which can resolve the problem faced when using partial McMD. We demonstrate that ALSD greatly increases the sampling efficiency over a conventional GE. We believe that ALSD is an effective method and is applicable to the conformational sampling of larger and more complicated biomolecule systems.
Collapse
Affiliation(s)
- Jinzen Ikebe
- Molecular Modeling and Simulation Group, Japan Atomic Energy Agency, 8-1-7 Umemidai, Kizugawa, Kyoto, 619-0215, Japan
| | | | | |
Collapse
|
34
|
Jiang F, Han W, Wu YD. The intrinsic conformational features of amino acids from a protein coil library and their applications in force field development. Phys Chem Chem Phys 2013; 15:3413-28. [PMID: 23385383 DOI: 10.1039/c2cp43633g] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The local conformational (φ, ψ, χ) preferences of amino acid residues remain an active research area, which are important for the development of protein force fields. In this perspective article, we first summarize spectroscopic studies of alanine-based short peptides in aqueous solution. While most studies indicate a preference for the P(II) conformation in the unfolded state over α and β conformations, significant variations are also observed. A statistical analysis from various coil libraries of high-resolution protein structures is then summarized, which gives a more coherent view of the local conformational features. The φ, ψ, χ distributions of the 20 amino acids have been obtained from a protein coil library, considering both backbone and side-chain conformational preferences. The intrinsic side-chain χ(1) rotamer preference and χ(1)-dependent Ramachandran plot can be generally understood by combining the interaction of the side-chain Cγ/Oγ atom with two neighboring backbone peptide groups. Current all-atom force fields such as AMBER ff99sb-ILDN, ff03 and OPLS-AA/L do not reproduce these distributions well. A method has been developed by combining the φ, ψ plot of alanine with the influence of side-chain χ(1) rotamers to derive the local conformational features of various amino acids. It has been further applied to improve the OPLS-AA force field. The modified force field (OPLS-AA/C) reproduces experimental (3)J coupling constants for various short peptides quite well. It also better reproduces the temperature-dependence of the helix-coil transition for alanine-based peptides. The new force field can fold a series of peptides and proteins with various secondary structures to their experimental structures. MD simulations of several globular proteins using the improved force field give significantly less deviation (RMSD) to experimental structures. The results indicate that the local conformational features from coil libraries are valuable for the development of balanced protein force fields.
Collapse
Affiliation(s)
- Fan Jiang
- Laboratory of Computational Chemistry and Drug Design, Laboratory of Chemical Genomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | | | | |
Collapse
|
35
|
Sakae Y, Okamoto Y. Amino-acid-dependent main-chain torsion-energy terms for protein systems. J Chem Phys 2013; 138:064103. [PMID: 23425457 DOI: 10.1063/1.4774159] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Many commonly used force fields for protein systems such as AMBER, CHARMM, GROMACS, OPLS, and ECEPP have amino-acid-independent force-field parameters for main-chain torsion-energy terms. Here, we propose a new type of amino-acid-dependent torsion-energy terms in the force fields. As an example, we applied this approach to AMBER ff03 force field and determined new amino-acid-dependent parameters for ψ (N-C(α)-C-N) and ζ (C(β)-C(α)-C-N) angles for each amino acid by using our optimization method, which is one of the knowledge-based approach. In order to test the validity of the new force-field parameters, we then performed folding simulations of α-helical and β-hairpin peptides, using the optimized force field. The results showed that the new force-field parameters gave structures more consistent with the experimental implications than the original AMBER ff03 force field.
Collapse
Affiliation(s)
- Yoshitake Sakae
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan
| | | |
Collapse
|
36
|
Li Y, Gao Y, Zhang X, Wang X, Mou L, Duan L, He X, Mei Y, Zhang JZH. A coupled two-dimensional main chain torsional potential for protein dynamics: generation and implementation. J Mol Model 2013; 19:3647-57. [PMID: 23765039 DOI: 10.1007/s00894-013-1879-8] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2013] [Accepted: 05/01/2013] [Indexed: 11/29/2022]
Abstract
Main chain torsions of alanine dipeptide are parameterized into coupled 2-dimensional Fourier expansions based on quantum mechanical (QM) calculations at M06 2X/aug-cc-pvtz//HF/6-31G** level. Solvation effect is considered by employing polarizable continuum model. Utilization of the M06 2X functional leads to precise potential energy surface that is comparable to or even better than MP2 level, but with much less computational demand. Parameterization of the 2D expansions is against the full main chain torsion space instead of just a few low energy conformations. This procedure is similar to that for the development of AMBER03 force field, except unique weighting factor was assigned to all the grid points. To avoid inconsistency between quantum mechanical calculations and molecular modeling, the model peptide is further optimized at molecular mechanics level with main chain dihedral angles fixed before the calculation of the conformational energy on molecular mechanical level at each grid point, during which generalized Born model is employed. Difference in solvation models at quantum mechanics and molecular mechanics levels makes this parameterization procedure less straightforward. All force field parameters other than main chain torsions are taken from existing AMBER force field. With this new main chain torsion terms, we have studied the main chain dihedral distributions of ALA dipeptide and pentapeptide in aqueous solution. The results demonstrate that 2D main chain torsion is effective in delineating the energy variation associated with rotations along main chain dihedrals. This work is an implication for the necessity of more accurate description of main chain torsions in the future development of ab initio force field and it also raises a challenge to the development of quantum mechanical methods, especially the quantum mechanical solvation models.
Collapse
Affiliation(s)
- Yongxiu Li
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy and Department of Physics and Institute of Theoretical and Computational Science, East China Normal University, Shanghai, 200062, China
| | | | | | | | | | | | | | | | | |
Collapse
|
37
|
Higo J, Umezawa K, Nakamura H. A virtual-system coupled multicanonical molecular dynamics simulation: Principles and applications to free-energy landscape of protein–protein interaction with an all-atom model in explicit solvent. J Chem Phys 2013; 138:184106. [DOI: 10.1063/1.4803468] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
|
38
|
Sakae Y, Okamoto Y. Improvement of the backbone-torsion-energy term in the force field for protein systems by the double Fourier series expansion. MOLECULAR SIMULATION 2013. [DOI: 10.1080/08927022.2012.705432] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
39
|
Toward Structure Prediction for Short Peptides Using the Improved SAAP Force Field Parameters. J CHEM-NY 2013. [DOI: 10.1155/2013/407862] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
|
40
|
Chen M, Cuendet MA, Tuckerman ME. Heating and flooding: A unified approach for rapid generation of free energy surfaces. J Chem Phys 2012; 137:024102. [DOI: 10.1063/1.4733389] [Citation(s) in RCA: 58] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
|
41
|
Higo J, Ikebe J, Kamiya N, Nakamura H. Enhanced and effective conformational sampling of protein molecular systems for their free energy landscapes. Biophys Rev 2012; 4:27-44. [PMID: 22347892 PMCID: PMC3271212 DOI: 10.1007/s12551-011-0063-6] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2011] [Accepted: 11/23/2011] [Indexed: 11/29/2022] Open
Abstract
Protein folding and protein-ligand docking have long persisted as important subjects in biophysics. Using multicanonical molecular dynamics (McMD) simulations with realistic expressions, i.e., all-atom protein models and an explicit solvent, free-energy landscapes have been computed for several systems, such as the folding of peptides/proteins composed of a few amino acids up to nearly 60 amino-acid residues, protein-ligand interactions, and coupled folding and binding of intrinsically disordered proteins. Recent progress in conformational sampling and its applications to biophysical systems are reviewed in this report, including descriptions of several outstanding studies. In addition, an algorithm and detailed procedures used for multicanonical sampling are presented along with the methodology of adaptive umbrella sampling. Both methods control the simulation so that low-probability regions along a reaction coordinate are sampled frequently. The reaction coordinate is the potential energy for multicanonical sampling and is a structural identifier for adaptive umbrella sampling. One might imagine that this probability control invariably enhances conformational transitions among distinct stable states, but this study examines the enhanced conformational sampling of a simple system and shows that reasonably well-controlled sampling slows the transitions. This slowing is induced by a rapid change of entropy along the reaction coordinate. We then provide a recipe to speed up the sampling by loosening the rapid change of entropy. Finally, we report all-atom McMD simulation results of various biophysical systems in an explicit solvent.
Collapse
Affiliation(s)
- Junichi Higo
- Institute for Protein Research, Osaka University, Suita, Osaka 565-0871 Japan
| | - Jinzen Ikebe
- Institute for Protein Research, Osaka University, Suita, Osaka 565-0871 Japan
| | - Narutoshi Kamiya
- Institute for Protein Research, Osaka University, Suita, Osaka 565-0871 Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, Suita, Osaka 565-0871 Japan
| |
Collapse
|
42
|
Conformational Ensembles of an Intrinsically Disordered Protein pKID with and without a KIX Domain in Explicit Solvent Investigated by All-Atom Multicanonical Molecular Dynamics. Biomolecules 2012; 2:104-21. [PMID: 24970129 PMCID: PMC4030872 DOI: 10.3390/biom2010104] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2011] [Revised: 02/11/2012] [Accepted: 02/12/2012] [Indexed: 12/13/2022] Open
Abstract
The phosphorylated kinase-inducible activation domain (pKID) adopts a helix–loop–helix structure upon binding to its partner KIX, although it is unstructured in the unbound state. The N-terminal and C-terminal regions of pKID, which adopt helices in the complex, are called, respectively, αA and αB. We performed all-atom multicanonical molecular dynamics simulations of pKID with and without KIX in explicit solvents to generate conformational ensembles. Although the unbound pKID was disordered overall, αA and αB exhibited a nascent helix propensity; the propensity of αA was stronger than that of αB, which agrees with experimental results. In the bound state, the free-energy landscape of αB involved two low free-energy fractions: native-like and non-native fractions. This result suggests that αB folds according to the induced-fit mechanism. The αB-helix direction was well aligned as in the NMR complex structure, although the αA helix exhibited high flexibility. These results also agree quantitatively with experimental observations. We have detected that the αB helix can bind to another site of KIX, to which another protein MLL also binds with the adopting helix. Consequently, MLL can facilitate pKID binding to the pKID-binding site by blocking the MLL-binding site. This also supports experimentally obtained results.
Collapse
|
43
|
Hudecová J, Hopmann KH, Bouř P. Correction of Vibrational Broadening in Molecular Dynamics Clusters with the Normal Mode Optimization Method. J Phys Chem B 2011; 116:336-42. [DOI: 10.1021/jp208785a] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Jana Hudecová
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo nám. 2, 166 10 Prague, Czech Republic
- Institute of Physics, Charles University, Ke Karlovu 5, 12116 Prague 2, Czech Republic
| | - Kathrin H. Hopmann
- Department of Chemistry, Centre for Theoretical and Computational Chemistry, University of Tromsø, N-9037 Tromsø, Norway
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo nám. 2, 166 10 Prague, Czech Republic
| |
Collapse
|
44
|
Higo J, Nishimura Y, Nakamura H. A free-energy landscape for coupled folding and binding of an intrinsically disordered protein in explicit solvent from detailed all-atom computations. J Am Chem Soc 2011; 133:10448-58. [PMID: 21627111 DOI: 10.1021/ja110338e] [Citation(s) in RCA: 91] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
The N-terminal repressor domain of neural restrictive silencer factor (NRSF) is an intrinsically disordered protein (IDP) that binds to the paired amphipathic helix (PAH) domain of mSin3. An NMR experiment revealed that the minimal binding unit of NRSF is a 15-residue segment that adopts a helical structure upon binding to a cleft of mSin3. We computed a free-energy landscape of this system by an enhanced conformational sampling method, all-atom multicanonical molecular dynamics. The simulation started from a configuration where the NRSF segment was fully disordered and distant from mSin3 in explicit solvent. In the absence of mSin3, the disordered NRSF segment thermally fluctuated between hairpins, helices, and bent structures. In the presence of mSin3, the segment bound to mSin3 by adopting the structures involved in the isolated state, and non-native and native complexes were formed. The free-energy landscape comprised three superclusters, and free-energy barriers separated the superclusters. The native complex was located at the center of the lowest free-energy cluster. When NRSF landed in the largest supercluster, the generated non-native complex moved on the landscape to fold into the native complex, by increasing the interfacial hydrophobic contacts and the helix content. When NRSF landed in other superclusters, the non-native complex overcame the free-energy barriers between the various segment orientations in the binding cleft of mSin3. Both population-shift and induced-fit (or induced-folding) mechanisms work cooperatively in the coupled folding and binding. The diverse structural adaptability of NRSF may be related to the hub properties of the IDP.
Collapse
Affiliation(s)
- Junichi Higo
- Institute for Protein Research, Osaka University, Suita, Osaka 565-0871, Japan.
| | | | | |
Collapse
|
45
|
Ikebe J, Standley DM, Nakamura H, Higo J. Ab initio simulation of a 57-residue protein in explicit solvent reproduces the native conformation in the lowest free-energy cluster. Protein Sci 2011; 20:187-96. [PMID: 21082745 DOI: 10.1002/pro.553] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
An enhanced conformational sampling method, multicanonical molecular dynamics (McMD), was applied to the ab intio folding of the 57-residue first repeat of human glutamyl- prolyl-tRNA synthetase (EPRS-R1) in explicit solvent. The simulation started from a fully extended structure of EPRS-R1 and did not utilize prior structural knowledge. A canonical ensemble, which is a conformational ensemble thermodynamically probable at an arbitrary temperature, was constructed by reweighting the sampled structures. Conformational clusters were obtained from the canonical ensemble at 300 K, and the largest cluster (i.e., the lowest free-energy cluster), which contained 34% of the structures in the ensemble, was characterized by the highest similarity to the NMR structure relative to all alternative clusters. This lowest free-energy cluster included native-like structures composed of two anti-parallel α-helices. The canonical ensemble at 300 K also showed that a short Gly-containing segment, which adopts an α-helix in the native structure, has a tendency to be structurally disordered. Atomic-level analyses demonstrated clearly that inter-residue hydrophobic interactions drive the helix formation of the Gly-containing segment, and that increasing the hydrophobic contacts accompanies exclusion of water molecules from the vicinity of this segment. This study has shown, for the first time, that the free-energy landscape of a structurally well-ordered protein of about 60 residues is obtainable with an all atom model in explicit water without prior structural knowledge.
Collapse
Affiliation(s)
- Jinzen Ikebe
- Graduate School of Frontier Biosciences, Osaka University, Open Laboratories for Advanced Bioscience and Biotechnology, Suita, Osaka 565-0874, Japan
| | | | | | | |
Collapse
|
46
|
Ikebe J, Umezawa K, Kamiya N, Sugihara T, Yonezawa Y, Takano Y, Nakamura H, Higo J. Theory for trivial trajectory parallelization of multicanonical molecular dynamics and application to a polypeptide in water. J Comput Chem 2010; 32:1286-97. [PMID: 21425286 DOI: 10.1002/jcc.21710] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2010] [Revised: 09/09/2010] [Accepted: 10/14/2010] [Indexed: 11/10/2022]
Abstract
Trivial trajectory parallelization of multicanonical molecular dynamics (TTP-McMD) explores the conformational space of a biological system with multiple short runs of McMD starting from various initial structures. This method simply connects (i.e., trivially parallelizes) the short trajectories and generates a long trajectory. First, we theoretically prove that the simple trajectory connection satisfies a detailed balance automatically. Thus, the resultant long trajectory is regarded as a single multicanonical trajectory. Second, we applied TTP-McMD to an alanine decapeptide with an all-atom model in explicit water to compute a free-energy landscape. The theory imposes two requirements on the multiple trajectories. We have demonstrated that TTP-McMD naturally satisfies the requirements. The TTP-McMD produces the free-energy landscape considerably faster than a single-run McMD does. We quantitatively showed that the accuracy of the computed landscape increases with increasing the number of multiple runs. Generally, the free-energy landscape of a large biological system is unknown a priori. The current method is suitable for conformational sampling of such a large system to reduce the waiting time to obtain a canonical ensemble statistically reliable.
Collapse
Affiliation(s)
- Jinzen Ikebe
- Graduate School of Frontier Biosciences, Osaka University, Osaka, Japan
| | | | | | | | | | | | | | | |
Collapse
|
47
|
Sakae Y, Okamoto Y. Optimisation of OPLS–UA force-field parameters for protein systems using protein data bank. MOLECULAR SIMULATION 2010. [DOI: 10.1080/08927022.2010.507247] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
48
|
Hudecová J, Kapitán J, Baumruk V, Hammer RP, Keiderling TA, Bouř P. Side Chain and Flexibility Contributions to the Raman Optical Activity Spectra of a Model Cyclic Hexapeptide. J Phys Chem A 2010; 114:7642-51. [DOI: 10.1021/jp104744a] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jana Hudecová
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo nám. 2, Prague 6, 16610, Czech Republic, Institute of Physics, Charles University, Ke Karlovu 5, 12116, Prague 2, Czech Republic, Department of Chemistry, Louisiana State University, 232 Choppin Hall, Baton Rouge, Louisiana 70803-1804, and Department of Chemistry, University of Illinois at Chicago, 845 West Taylor Street, Chicago Illinois 60607-7061
| | - Josef Kapitán
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo nám. 2, Prague 6, 16610, Czech Republic, Institute of Physics, Charles University, Ke Karlovu 5, 12116, Prague 2, Czech Republic, Department of Chemistry, Louisiana State University, 232 Choppin Hall, Baton Rouge, Louisiana 70803-1804, and Department of Chemistry, University of Illinois at Chicago, 845 West Taylor Street, Chicago Illinois 60607-7061
| | - Vladimír Baumruk
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo nám. 2, Prague 6, 16610, Czech Republic, Institute of Physics, Charles University, Ke Karlovu 5, 12116, Prague 2, Czech Republic, Department of Chemistry, Louisiana State University, 232 Choppin Hall, Baton Rouge, Louisiana 70803-1804, and Department of Chemistry, University of Illinois at Chicago, 845 West Taylor Street, Chicago Illinois 60607-7061
| | - Robert P. Hammer
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo nám. 2, Prague 6, 16610, Czech Republic, Institute of Physics, Charles University, Ke Karlovu 5, 12116, Prague 2, Czech Republic, Department of Chemistry, Louisiana State University, 232 Choppin Hall, Baton Rouge, Louisiana 70803-1804, and Department of Chemistry, University of Illinois at Chicago, 845 West Taylor Street, Chicago Illinois 60607-7061
| | - Timothy A. Keiderling
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo nám. 2, Prague 6, 16610, Czech Republic, Institute of Physics, Charles University, Ke Karlovu 5, 12116, Prague 2, Czech Republic, Department of Chemistry, Louisiana State University, 232 Choppin Hall, Baton Rouge, Louisiana 70803-1804, and Department of Chemistry, University of Illinois at Chicago, 845 West Taylor Street, Chicago Illinois 60607-7061
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo nám. 2, Prague 6, 16610, Czech Republic, Institute of Physics, Charles University, Ke Karlovu 5, 12116, Prague 2, Czech Republic, Department of Chemistry, Louisiana State University, 232 Choppin Hall, Baton Rouge, Louisiana 70803-1804, and Department of Chemistry, University of Illinois at Chicago, 845 West Taylor Street, Chicago Illinois 60607-7061
| |
Collapse
|
49
|
Sakae Y, Okamoto Y. Controlling the secondary-structure-forming tendencies of proteins by a backbone torsion-energy term. MOLECULAR SIMULATION 2010. [DOI: 10.1080/08927020903124601] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
50
|
Verifying trivial parallelization of multicanonical molecular dynamics for conformational sampling of a polypeptide in explicit water. Chem Phys Lett 2009. [DOI: 10.1016/j.cplett.2009.03.077] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|