1
|
Okamoto C, Ando K. Molecular dynamics simulation analysis of structural dynamic cross correlation induced by odorant hydrogen-bonding in mouse eugenol ol- factory receptor. Biophys Physicobiol 2024; 21:e210007. [PMID: 38803338 PMCID: PMC11128758 DOI: 10.2142/biophysico.bppb-v21.0007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Accepted: 01/16/2024] [Indexed: 05/29/2024] Open
Abstract
Structural fluctuations and dynamic cross-correlations in the mouse eugenol olfactory receptor (Olfr73) were studied by molecular dynamics (MD) simulation to characterize the dynamic response of the protein upon ligand binding. The initial structure was generated by the artificial intelligence tool AlphaFold2 due to the current lack of experimental data. We focused on the hydrogen (H) bond of the odorant eugenol to Ser113, Asn207, and Tyr260 of the receptor protein, the importance of which has been suggested by previous experimental studies. The H-bond was not observed in docking simulations, but in subsequent MD simulations the H-bond to Ser113 was formed in 2-4 ns. The lifetime of the H-bond was in the range of 1-20 ns. On the trajectory with the most stable (20 ns) H-bond, the structural fluctuation of the α-carbon atoms of the receptor main chain was studied by calculating the root mean square fluctuations, the dynamic cross-correlation map, and the time-dependent dynamic cross-correlation. The analysis suggested a correlation transfer pathway Ser113 → Phe182 → (Leu259 or Tyr260) → Tyr291 induced by the ligand binding with a time scale of 4-6 ns.
Collapse
Affiliation(s)
- Chisato Okamoto
- Department of Information and Sciences, Tokyo Woman’s Christian University, Tokyo 167-8585, Japan
| | - Koji Ando
- Department of Information and Sciences, Tokyo Woman’s Christian University, Tokyo 167-8585, Japan
| |
Collapse
|
2
|
Fukuda I, Moritsugu K, Higo J, Fukunishi Y. A cutoff-based method with charge-distribution-data driven pair potentials for efficiently estimating electrostatic interactions in molecular systems. J Chem Phys 2023; 159:234116. [PMID: 38112509 DOI: 10.1063/5.0172270] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2023] [Accepted: 11/27/2023] [Indexed: 12/21/2023] Open
Abstract
We introduce a simple cutoff-based method for precise electrostatic energy calculations in the molecular dynamics (MD) simulations of point-particle systems. Our method employs a theoretically derived smooth pair potential function to define electrostatic energy, offering stability and computational efficiency in MD simulations. Instead of imposing specific physical conditions, such as dielectric environments or charge neutrality, we focus on the relationship represented by a single summation formula of charge-weighted pair potentials. This approach allows an accurate energy approximation for each particle, enabling a straightforward error analysis. The resulting particle-dependent pair potential captures the charge distribution information, making it suitable for heterogeneous systems and ensuring an enhanced accuracy through distant information inclusion. Numerical investigations of the Madelung constants of crystalline systems validate the method's accuracy.
Collapse
Affiliation(s)
- Ikuo Fukuda
- Graduate School of Science, Osaka Metropolitan University, 1-2 Gakuen-cho, Naka-ku, Sakai, Osaka 599-8231, Japan
- Japan Biological Informatics Consortium, 2-4-32 Aomi, Koto-ku, Tokyo, 135-8073, Japan
| | - Kei Moritsugu
- Graduate School of Science, Osaka Metropolitan University, 1-2 Gakuen-cho, Naka-ku, Sakai, Osaka 599-8231, Japan
| | - Junichi Higo
- Graduate School of Information Science, University of Hyogo, 7-1-28 Minatojima Minamimachi, Chuo-ku, Kobe, Hyogo650-0047, Japan
- Research Organization of Science and Technology, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, 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
| |
Collapse
|
3
|
Fukuda I, Nakamura H. Non-Ewald methods for evaluating the electrostatic interactions of charge systems: similarity and difference. Biophys Rev 2022; 14:1315-1340. [PMID: 36659982 PMCID: PMC9842848 DOI: 10.1007/s12551-022-01029-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2022] [Accepted: 11/30/2022] [Indexed: 01/13/2023] Open
Abstract
In molecular simulations, it is essential to properly calculate the electrostatic interactions of particles in the physical system of interest. Here we consider a method called the non-Ewald method, which does not rely on the standard Ewald method with periodic boundary conditions, but instead relies on the cutoff-based techniques. We focus on the physicochemical and mathematical conceptual aspects of the method in order to gain a deeper understanding of the simulation methodology. In particular, we take into account the reaction field (RF) method, the isotropic periodic sum (IPS) method, and the zero-multipole summation method (ZMM). These cutoff-based methods are based on different physical ideas and are completely distinguishable in their underlying concepts. The RF and IPS methods are "additive" methods that incorporate information outside the cutoff region, via dielectric medium and isotropic boundary condition, respectively. In contrast, the ZMM is a "subtraction" method that tries to remove the artificial effects, generated near the boundary, from the cutoff sphere. Nonetheless, we find physical and/or mathematical similarities between these methods. In particular, the modified RF method can be derived by the principle of neutralization utilized in the ZMM, and we also found a direct relationship between IPS and ZMM.
Collapse
Affiliation(s)
- Ikuo Fukuda
- Graduate School of Information Science, University of Hyogo, 7-1-28 Minatojima, Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047 Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871 Japan
| |
Collapse
|
4
|
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
|
5
|
Harada R, Yamaguchi K, Shigeta Y. Enhanced Conformational Sampling Method Based on Anomaly Detection Parallel Cascade Selection Molecular Dynamics: ad-PaCS-MD. J Chem Theory Comput 2020; 16:6716-6725. [PMID: 32926622 DOI: 10.1021/acs.jctc.0c00697] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In the present study, we propose a rare-event sampling method called anomaly detection parallel cascade selection molecular dynamics (ad-PaCS-MD). The original PaCS-MD was designed to generate conformational transition pathways from a given reactant to a product when the latter is known a priori. As an extension of the original method, ad-PaCS-MD has been designed to efficiently search transition pathways from a given reactant without referring to a given product. In ad-PaCS-MD, rarely occurring but essential states (configurations) of proteins for the transitions are identified based on the degrees of an anomaly. In more detail, ad-PaCS-MD adopts an algorithm called an anomaly detection generative adversarial network (anoGAN) as a measure for detecting rarely occurring states to be resampled. Here, the essential configurations with higher degrees of the anomaly are selected with anoGAN and intensively resampled by restarting short-time MD simulations from the selected configurations. By repeating the detections and resampling of configurations with the higher degrees of the anomaly, ad-PaCS-MD automatically and efficiently promotes the rare events and gives a wide range of the free energy landscape by combining with the Markov state model construction. As demonstrations, open-closed transitions of two globular proteins (T4 lysozyme and maltose-binding protein) were promoted with ad-PaCS-MD by referring only to the given starting configurations. In each demonstration, ad-PaCS-MD promoted the large-amplitude open-closed transitions with nanosecond-order simulation times. In conclusion, our demonstrations showed a higher conformational sampling efficiency for ad-PaCS-MD than conventional MD (CMD) because CMD required computational costs of more than microsecond-order simulation times to promote the rare events.
Collapse
Affiliation(s)
- Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan
| | - Kota Yamaguchi
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan
| |
Collapse
|
6
|
Whalley T, Dolton G, Brown PE, Wall A, Wooldridge L, van den Berg H, Fuller A, Hopkins JR, Crowther MD, Attaf M, Knight RR, Cole DK, Peakman M, Sewell AK, Szomolay B. GPU-Accelerated Discovery of Pathogen-Derived Molecular Mimics of a T-Cell Insulin Epitope. Front Immunol 2020; 11:296. [PMID: 32184781 PMCID: PMC7058665 DOI: 10.3389/fimmu.2020.00296] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Accepted: 02/05/2020] [Indexed: 01/09/2023] Open
Abstract
The strong links between (Human Leukocyte Antigen) HLA, infection and autoimmunity combine to implicate T-cells as primary triggers of autoimmune disease (AD). T-cell crossreactivity between microbially-derived peptides and self-peptides has been shown to break tolerance and trigger AD in experimental animal models. Detailed examination of the potential for T-cell crossreactivity to trigger human AD will require means of predicting which peptides might be recognised by autoimmune T-cell receptors (TCRs). Recent developments in high throughput sequencing and bioinformatics mean that it is now possible to link individual TCRs to specific pathologies for the first time. Deconvolution of TCR function requires knowledge of TCR specificity. Positional Scanning Combinatorial Peptide Libraries (PS-CPLs) can be used to predict HLA-restriction and define antigenic peptides derived from self and pathogen proteins. In silico search of the known terrestrial proteome with a prediction algorithm that ranks potential antigens in order of recognition likelihood requires complex, large-scale computations over several days that are infeasible on a personal computer. We decreased the time required for peptide searching to under 30 min using multiple blocks on graphics processing units (GPUs). This time-efficient, cost-effective hardware accelerator was used to screen bacterial and fungal human pathogens for peptide sequences predicted to activate a T-cell clone, InsB4, that was isolated from a patient with type 1 diabetes and recognised the insulin B-derived epitope HLVEALYLV in the context of disease-risk allele HLA A*0201. InsB4 was shown to kill HLA A*0201+ human insulin producing β-cells demonstrating that T-cells with this specificity might contribute to disease. The GPU-accelerated algorithm and multispecies pathogen proteomic databases were validated to discover pathogen-derived peptide sequences that acted as super-agonists for the InsB4 T-cell clone. Peptide-MHC tetramer binding and surface plasmon resonance were used to confirm that the InsB4 TCR bound to the highest-ranked peptide agonists derived from infectious bacteria and fungi. Adoption of GPU-accelerated prediction of T-cell agonists has the capacity to revolutionise our understanding of AD by identifying potential targets for autoimmune T-cells. This approach has further potential for dissecting T-cell responses to infectious disease and cancer.
Collapse
Affiliation(s)
- Thomas Whalley
- Cardiff University School of Medicine, Cardiff, United Kingdom.,Systems Immunity Research Institute, Cardiff University, Cardiff, United Kingdom
| | - Garry Dolton
- Cardiff University School of Medicine, Cardiff, United Kingdom
| | - Paul E Brown
- Zeeman Institute for Systems Biology and Infectious Disease Epidemiology Research, University of Warwick Coventry, Coventry, United Kingdom
| | - Aaron Wall
- Cardiff University School of Medicine, Cardiff, United Kingdom
| | - Linda Wooldridge
- Faculty of Health Sciences, University of Bristol, Bristol, United Kingdom
| | - Hugo van den Berg
- Mathematics Institute, University of Warwick, Coventry, United Kingdom
| | - Anna Fuller
- Cardiff University School of Medicine, Cardiff, United Kingdom
| | - Jade R Hopkins
- Cardiff University School of Medicine, Cardiff, United Kingdom
| | | | - Meriem Attaf
- Cardiff University School of Medicine, Cardiff, United Kingdom
| | - Robin R Knight
- Peter Gorer Department of Immunobiology, Guy's Hospital, London, United Kingdom
| | - David K Cole
- Cardiff University School of Medicine, Cardiff, United Kingdom
| | - Mark Peakman
- Peter Gorer Department of Immunobiology, Guy's Hospital, London, United Kingdom
| | - Andrew K Sewell
- Cardiff University School of Medicine, Cardiff, United Kingdom.,Systems Immunity Research Institute, Cardiff University, Cardiff, United Kingdom
| | - Barbara Szomolay
- Cardiff University School of Medicine, Cardiff, United Kingdom.,Systems Immunity Research Institute, Cardiff University, Cardiff, United Kingdom
| |
Collapse
|
7
|
OKAMOTO C, ANDO K. Molecular Dynamics Simulation and Binding Analysis of Olfactory Receptor Protein mOR-EG and Odorant Molecule Eugenol. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2020. [DOI: 10.2477/jccj.2021-0005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Chisato OKAMOTO
- Department of Mathematical Sciences, Tokyo Woman's Christian University, 2-6-1 Zenpuku-ji,Suginami, Tokyo 167-8585
| | - Koji ANDO
- Department of Mathematical Sciences, Tokyo Woman's Christian University, 2-6-1 Zenpuku-ji,Suginami, Tokyo 167-8585
| |
Collapse
|
8
|
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
|
9
|
Sánchez-Monroy X, Torres-Arenas J, Gil-Villegas A. Theoretical equations of state for a charged fluid. J Chem Phys 2019; 150:144507. [PMID: 30981249 DOI: 10.1063/1.5063577] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this article, we present a molecular thermodynamic study of a system of N particles contained within a volume V and interacting via a hard-core pair potential with an attractive interaction according to the Wolf model for charged systems. This variable-range potential is characterized by three parameters: the repulsive hard-core diameter σ, the energy-well depth ϵ, and the inverse range α; a fourth parameter of the model is a cut-off distance xc that depends on α according to the relation xc = 2/α. Two equations of state (EOSs) are presented and derived from thermodynamic perturbation theory and Monte Carlo (MC) simulation data. The first EOS is given by the standard Zwanzig's high-temperature expansion of the Helmholtz free energy, where the first three perturbation terms a1, a2, and a3 were obtained from MC simulations in the canonical ensemble (NVT) and parameterized as functions of α and the reduced density of particles ρ* = Nσ3/V. The second EOS was obtained from the discrete perturbation theory applied to a discrete representation of the Wolf potential. Results for pressures, internal energies, and isochoric heat capacities are compared to the MC computer simulation data of the Wolf system, including vapor-liquid coexistence curves, for different values of α. Overall, both EOSs give a very good representation of the thermodynamic properties of the Wolf fluid when 0.3 ≤ α ≤ 1.0 and 0.05 ≤ ρ* ≤ 0.8. Since the Yukawa fluid can reproduce information of screened ionic interactions, we discuss the equivalence between the Wolf and Yukawa fluids in the context of equivalent systems in liquid theory.
Collapse
Affiliation(s)
- X Sánchez-Monroy
- División de Ciencias e Ingenierías, Campus León, Universidad de Guanajuato, Loma del Bosque 103, Lomas del Campestre, 37150 León Guanajuato, Mexico
| | - J Torres-Arenas
- División de Ciencias e Ingenierías, Campus León, Universidad de Guanajuato, Loma del Bosque 103, Lomas del Campestre, 37150 León Guanajuato, Mexico
| | - A Gil-Villegas
- División de Ciencias e Ingenierías, Campus León, Universidad de Guanajuato, Loma del Bosque 103, Lomas del Campestre, 37150 León Guanajuato, Mexico
| |
Collapse
|
10
|
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
|
11
|
Bekker GJ, Ma B, Kamiya N. Thermal stability of single-domain antibodies estimated by molecular dynamics simulations. Protein Sci 2018; 28:429-438. [PMID: 30394618 PMCID: PMC6319760 DOI: 10.1002/pro.3546] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2018] [Revised: 10/03/2018] [Accepted: 10/28/2018] [Indexed: 12/20/2022]
Abstract
Single‐domain antibodies (sdAbs) function like regular antibodies, however, consist of only one domain. Because of their low molecular weight, sdAbs have advantages with respect to production and delivery to their targets and for applications such as antibody drugs and biosensors. Thus, sdAbs with high thermal stability are required. In this work, we chose seven sdAbs, which have a wide range of melting temperature (Tm) values and known structures. We applied molecular dynamics (MD) simulations to estimate their relative stability and compared them with the experimental data. High‐temperature MD simulations at 400 K and 500 K were executed with simulations at 300 K as a control. The fraction of native atomic contacts, Q, measured for the 400 K simulations showed a fairly good correlation with the Tm values. Interestingly, when the residues were classified by their hydrophobicity and size, the Q values of hydrophilic residues exhibited an even better correlation, suggesting that stabilization is correlated with favorable interactions of hydrophilic residues. Measuring the Q value on a per‐residue level enabled us to identify residues that contribute significantly to the instability and thus demonstrating how our analysis can be used in a mutant case study.
Collapse
Affiliation(s)
- Gert-Jan Bekker
- Institute for Protein Research, Osaka University, Suita, Osaka, Japan
| | - Benson Ma
- School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia, 30332
| | - Narutoshi Kamiya
- Graduate School of Simulation Studies, University of Hyogo, Kobe, Hyogo, Japan
| |
Collapse
|
12
|
Numoto N, Kamiya N, Bekker GJ, Yamagami Y, Inaba S, Ishii K, Uchiyama S, Kawai F, Ito N, Oda M. Structural Dynamics of the PET-Degrading Cutinase-like Enzyme from Saccharomonospora viridis AHK190 in Substrate-Bound States Elucidates the Ca 2+-Driven Catalytic Cycle. Biochemistry 2018; 57:5289-5300. [PMID: 30110540 DOI: 10.1021/acs.biochem.8b00624] [Citation(s) in RCA: 44] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A cutinase-type polyesterase from Saccharomonospora viridis AHK190 (Cut190) has been shown to degrade the inner block of polyethylene terephthalate. A unique feature of Cut190 is that its function and stability are regulated by Ca2+ binding. Our previous crystal structure analysis of Cut190S226P showed that one Ca2+ binds to the enzyme, which induces large conformational changes in several loop regions to stabilize an open conformation [Miyakawa, T., et al. (2015) Appl. Microbiol. Biotechnol. 99, 4297]. In this study, to analyze the substrate recognition mechanism of Cut190, we determined the crystal structure of the inactive form of a Cut190 mutant, Cut190*S176A, in complex with calcium ions and/or substrates. We found that three calcium ions bind to Cut190*S176A, which is supported by analysis using native mass spectrometry experiments and 3D Reference Interaction Site Model calculations. The complex structures with the two substrates, monoethyl succinate and monoethyl adipate (engaged and open forms), presumably correspond to the pre- and post-reaction states, as the ester bond is close to the active site and pointing outward from the active site, respectively, for the two complexes. Ca2+ binding induces the pocket to open, enabling the substrate to access the pocket more easily. Molecular dynamics simulations suggest that a post-reaction state in the engaged form presumably exists between the experimentally observed forms, indicating that the substrate would be cleaved in the engaged form and then requires the enzyme to change to the open form to release the product, a process that Ca2+ can greatly accelerate.
Collapse
Affiliation(s)
- Nobutaka Numoto
- Medical Research Institute , Tokyo Medical and Dental University , 1-5-45 Yushima , Bunkyo-ku, Tokyo 113-8510 , Japan
| | - Narutoshi Kamiya
- Institute for Protein Research , Osaka University , 3-2 Yamadaoka , Suita, Osaka 565-0871 , Japan.,Graduate School of Simulation Studies , University of Hyogo , 7-1-28 Minatojima Minami-machi , Chuo-ku, Kobe , Hyogo 650-0047 , Japan
| | - Gert-Jan Bekker
- Institute for Protein Research , Osaka University , 3-2 Yamadaoka , Suita, Osaka 565-0871 , Japan
| | - Yuri Yamagami
- Graduate School of Life and Environmental Sciences , Kyoto Prefectural University , 1-5 Hangi-cho, Shimogamo , Sakyo-ku, Kyoto , Kyoto 606-8522 , Japan
| | - Satomi Inaba
- Graduate School of Life and Environmental Sciences , Kyoto Prefectural University , 1-5 Hangi-cho, Shimogamo , Sakyo-ku, Kyoto , Kyoto 606-8522 , Japan.,Research & Utilization Division , Japan Synchrotron Radiation Research Institute , 1-1-1 Kouto , Sayo , Hyogo 679-5198 , Japan
| | - Kentaro Ishii
- Exploratory Research Center on Life and Living Systems , National Institutes of Natural Sciences , 5-1 Higashiyama, Myodaiji , Okazaki , Aichi 444-8787 , Japan
| | - Susumu Uchiyama
- Exploratory Research Center on Life and Living Systems , National Institutes of Natural Sciences , 5-1 Higashiyama, Myodaiji , Okazaki , Aichi 444-8787 , Japan.,Graduate School of Engineering , Osaka University , 2-1 Yamadaoka , Suita, Osaka 565-0871 , Japan
| | - Fusako Kawai
- Center for Fiber and Textile Science , Kyoto Institute of Technology , Matsugasaki , Sakyo-ku, Kyoto , Kyoto 606-8585 , Japan
| | - Nobutoshi Ito
- Medical Research Institute , Tokyo Medical and Dental University , 1-5-45 Yushima , Bunkyo-ku, Tokyo 113-8510 , Japan
| | - Masayuki Oda
- Graduate School of Life and Environmental Sciences , Kyoto Prefectural University , 1-5 Hangi-cho, Shimogamo , Sakyo-ku, Kyoto , Kyoto 606-8522 , Japan
| |
Collapse
|
13
|
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
|
14
|
Oda M, Inaba S, Kamiya N, Bekker GJ, Mikami B. Structural and thermodynamic characterization of endo-1,3-β-glucanase: Insights into the substrate recognition mechanism. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2018; 1866:415-425. [DOI: 10.1016/j.bbapap.2017.12.004] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/25/2017] [Revised: 11/24/2017] [Accepted: 12/11/2017] [Indexed: 11/25/2022]
|
15
|
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
|
16
|
Nobile MS, Cazzaniga P, Tangherloni A, Besozzi D. Graphics processing units in bioinformatics, computational biology and systems biology. Brief Bioinform 2017; 18:870-885. [PMID: 27402792 PMCID: PMC5862309 DOI: 10.1093/bib/bbw058] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2016] [Indexed: 01/18/2023] Open
Abstract
Several studies in Bioinformatics, Computational Biology and Systems Biology rely on the definition of physico-chemical or mathematical models of biological systems at different scales and levels of complexity, ranging from the interaction of atoms in single molecules up to genome-wide interaction networks. Traditional computational methods and software tools developed in these research fields share a common trait: they can be computationally demanding on Central Processing Units (CPUs), therefore limiting their applicability in many circumstances. To overcome this issue, general-purpose Graphics Processing Units (GPUs) are gaining an increasing attention by the scientific community, as they can considerably reduce the running time required by standard CPU-based software, and allow more intensive investigations of biological systems. In this review, we present a collection of GPU tools recently developed to perform computational analyses in life science disciplines, emphasizing the advantages and the drawbacks in the use of these parallel architectures. The complete list of GPU-powered tools here reviewed is available at http://bit.ly/gputools.
Collapse
Affiliation(s)
- Marco S Nobile
- Department of Informatics, Systems and Communication, University of Milano-Bicocca, Milano, Italy
- SYSBIO.IT Centre of Systems Biology, Milano, Italy
| | - Paolo Cazzaniga
- Department of Human and Social Sciences, University of Bergamo, Bergamo, Italy
- SYSBIO.IT Centre of Systems Biology, Milano, Italy
| | - Andrea Tangherloni
- Department of Informatics, Systems and Communication, University of Milano-Bicocca, Milano, Italy
| | - Daniela Besozzi
- Department of Informatics, Systems and Communication, University of Milano-Bicocca, Milano, Italy
- SYSBIO.IT Centre of Systems Biology, Milano, Italy
- Corresponding author. Daniela Besozzi, Department of Informatics, Systems and Communication, University of Milano-Bicocca, Milano, Italy and SYSBIO.IT Centre of Systems Biology, Milano, Italy. Tel.: +39 02 6448 7874. E-mail:
| |
Collapse
|
17
|
Kasahara K, Shiina M, Fukuda I, Ogata K, Nakamura H. Molecular mechanisms of cooperative binding of transcription factors Runx1-CBFβ-Ets1 on the TCRα gene enhancer. PLoS One 2017; 12:e0172654. [PMID: 28231333 PMCID: PMC5322934 DOI: 10.1371/journal.pone.0172654] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2016] [Accepted: 02/07/2017] [Indexed: 11/22/2022] Open
Abstract
Ets1 is an essential transcription factor (TF) for several important physiological processes, including cell proliferation and differentiation. Its recognition of the enhancer region of the TCRα gene is enhanced by the cooperative binding of the Runx1–CBFβ heterodimer, with the cancelation of phosphorylation-dependent autoinhibition. The detailed mechanism of this interesting cooperativity between Ets1 and the Runx1–CBFβ heterodimer is still largely unclear. Here, we investigated the molecular mechanisms of this cooperativity, by using molecular dynamics simulations. Consequently, we detected high flexibility of the loop region between the HI2 and H1 helices of Ets1. Upon Runx1–CBFβ heterodimer binding, this loop transiently adopts various sub-stable conformations in its interactions with the DNA. In addition, a network analysis suggested an allosteric pathway in the molecular assembly and identified some key residues that coincide with previous experimental studies. Our simulations suggest that the cooperative binding of Ets1 and the Runx1–CBFβ heterodimer alters the DNA conformation and induces sub-stable conformations of the HI2–H1 loop of Ets1. This phenomenon increases the flexibility of the regulatory module, including the HI2 helix, and destabilizes the inhibitory form of this module. Thus, we hypothesize that this effect facilitates Ets1–DNA binding and prevents the phosphorylation-dependent DNA binding autoinhibition.
Collapse
Affiliation(s)
- Kota Kasahara
- College of Life Sciences, Ritsumeikan University, Kusatsu, Shiga, Japan
- * E-mail:
| | - Masaaki Shiina
- Graduate School of Medicine, Yokohama City University, Kanazawa-ku, Yokohama, Kanagawa, Japan
| | - Ikuo Fukuda
- Institute for Protein Research, Osaka University, Suita, Osaka, Japan
| | - Kazuhiro Ogata
- Graduate School of Medicine, Yokohama City University, Kanazawa-ku, Yokohama, Kanagawa, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, Suita, Osaka, Japan
| |
Collapse
|
18
|
Harada R, Shigeta Y. Efficient Conformational Search Based on Structural Dissimilarity Sampling: Applications for Reproducing Structural Transitions of Proteins. J Chem Theory Comput 2017; 13:1411-1423. [PMID: 28170260 DOI: 10.1021/acs.jctc.6b01112] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Structural Dissimilarity Sampling (SDS) is proposed as an efficient conformational search method to promote structural transitions essential for the biological functions of proteins. In SDS, initial structures are selected based on structural dissimilarity, and conformational resampling is repeated. Conformational resampling is performed as follows: (I) arrangement of initial structures for a diverse distribution at the edge of a conformational subspace and (II) promotion of the structural transitions with multiple short-time molecular dynamics (MD) simulations restarting from the diversely distributed initial structures. Cycles of (I) and (II) are repeated to intensively promote structural transitions because conformational resampling from the initial structures would quickly expand conformational distributions toward unvisited conformational subspaces. As a demonstration, SDS was first applied to maltodextrin binding protein (MBP) in explicit water to reproduce structural transitions between the open and closed states of MBP. Structural transitions of MBP were successfully reproduced with SDS in nanosecond-order simulation times. Starting from both the open and closed forms, SDS successfully reproduced the structural transitions within 25 cycles (a total of 250 ns of simulation time). For reference, a conventional long-time (500 ns) MD simulation under NPT (300 K and 1 bar) starting from the open form failed to reproduce the structural transition. In addition to the open-closed motions of MBP, SDS was applied to folding processes of the fast-folding proteins (chignolin, Trp-cage, and villin) and successfully sampled their native states. To confirm how the selections of initial structures affected conformational sampling efficiency, numbers of base sets for characterizing structural dissimilarity of initial structures were addressed in distinct trials of SDS. The parameter searches showed that the conformational sampling efficiency was relatively insensitive with respect to the numbers of base sets, indicating the robustness of SDS for actual applications.
Collapse
Affiliation(s)
- Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba , Tennodai 1-1-1, Tsukuba, Ibaraki 305-8577, Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba , Tennodai 1-1-1, Tsukuba, Ibaraki 305-8577, Japan
| |
Collapse
|
19
|
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
|
20
|
Nishigami H, Kamiya N, Nakamura H. Revisiting antibody modeling assessment for CDR-H3 loop. Protein Eng Des Sel 2016; 29:477-484. [PMID: 27515703 PMCID: PMC5081041 DOI: 10.1093/protein/gzw028] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2016] [Revised: 06/14/2016] [Accepted: 06/14/2016] [Indexed: 02/05/2023] Open
Abstract
The antigen-binding site of antibodies, also known as complementarity-determining region (CDR), has hypervariable sequence properties. In particular, the third CDR loop of the heavy chain, CDR-H3, has such variability in its sequence, length, and conformation that ordinary modeling techniques cannot build a high-quality structure. At Stage 2 of the Second Antibody Modeling Assessment (AMA-II) held in 2013, the model structures of the CDR-H3 loops were submitted by the seven modelers and were critically assessed. After our participation in AMA-II, we rebuilt one of the long CDR-H3 loops with 13 residues (A52 antibody) by a more precise method, using enhanced conformational sampling with the explicit water model, as compared to our previous method employed at AMA-II. The current stable models obtained from the free energy landscape at 300 K include structures similar to the X-ray crystal structures. Those models were not built in our previous work at AMA-II. The current free energy landscape suggested that the CDR-H3 loop structures in the crystal are not stable in solution, but they are stabilized by the crystal packing effect.
Collapse
Affiliation(s)
- Hiroshi Nishigami
- Institute for Protein Research, Osaka University, 3-2, Yamadaoka, Suita, Osaka 565-0871, Japan
- Present address: Graduate School of Life Science, University of Hyogo, 3-2-1, Koto, Kamigori, Akoh, Hyogo 678-1297, Japan
| | - Narutoshi Kamiya
- Institute for Protein Research, Osaka University, 3-2, Yamadaoka, Suita, Osaka 565-0871, Japan
- Advanced Institute for Computational Science, RIKEN, QBiC Building B, 6-2-4, Furuedai, Suita, Osaka 565-0874, Japan
- Present address: Graduate School of Simulation Studies, University of Hyogo, 7-1-28, Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2, Yamadaoka, Suita, Osaka 565-0871, Japan
| |
Collapse
|
21
|
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]
|
22
|
Jung J, Naurse A, Kobayashi C, Sugita Y. Graphics Processing Unit Acceleration and Parallelization of GENESIS for Large-Scale Molecular Dynamics Simulations. J Chem Theory Comput 2016; 12:4947-4958. [DOI: 10.1021/acs.jctc.6b00241] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jaewoon Jung
- RIKEN Theoretical Molecular Science Laboratory, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- RIKEN Advanced Institute for Computational Science, 7-1-26 minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 640-0047, Japan
| | - Akira Naurse
- NVIDIA, 2-11-7, Akasaka,
Minato-ku, Tokyo 107-0052, Japan
| | - Chigusa Kobayashi
- RIKEN Advanced Institute for Computational Science, 7-1-26 minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 640-0047, Japan
| | - Yuji Sugita
- RIKEN Theoretical Molecular Science Laboratory, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- RIKEN Advanced Institute for Computational Science, 7-1-26 minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 640-0047, Japan
- RIKEN iTHES, 2-1 Hirosawa,
Wako, Saitama 351-0198, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN Quantitative Biology Center (QBiC), 6-7-1 minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
23
|
Shimba N, Kamiya N, Nakamura H. Model Building of Antibody–Antigen Complex Structures Using GBSA Scores. J Chem Inf Model 2016; 56:2005-2012. [DOI: 10.1021/acs.jcim.6b00066] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Noriko Shimba
- Device
Research Laboratory, Advanced Research Division, Panasonic Corporation, 3-4 Hikaridai, Seika-cho, Soraku-gun, Kyoto 619-0237, Japan
| | - Narutoshi Kamiya
- Advanced
Institute for Computational Science, RIKEN, QBiC Building B, 6-2-4 Furuedai, Suita, Osaka 565-0874, Japan
- 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
|
Kasahara K, Ma B, Goto K, Dasgupta B, Higo J, Fukuda I, Mashimo T, Akiyama Y, Nakamura H. myPresto/omegagene: a GPU-accelerated molecular dynamics simulator tailored for enhanced conformational sampling methods with a non-Ewald electrostatic scheme. Biophys Physicobiol 2016; 13:209-216. [PMID: 27924276 PMCID: PMC5060096 DOI: 10.2142/biophysico.13.0_209] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2016] [Accepted: 08/08/2016] [Indexed: 12/01/2022] Open
Abstract
Molecular dynamics (MD) is a promising computational approach to investigate dynamical behavior of molecular systems at the atomic level. Here, we present a new MD simulation engine named "myPresto/omegagene" that is tailored for enhanced conformational sampling methods with a non-Ewald electrostatic potential scheme. Our enhanced conformational sampling methods, e.g., the virtual-system-coupled multi-canonical MD (V-McMD) method, replace a multi-process parallelized run with multiple independent runs to avoid inter-node communication overhead. In addition, adopting the non-Ewald-based zero-multipole summation method (ZMM) makes it possible to eliminate the Fourier space calculations altogether. The combination of these state-of-the-art techniques realizes efficient and accurate calculations of the conformational ensemble at an equilibrium state. By taking these advantages, myPresto/omegagene is specialized for the single process execution with Graphics Processing Unit (GPU). We performed benchmark simulations for the 20-mer peptide, Trp-cage, with explicit solvent. One of the most thermodynamically stable conformations generated by the V-McMD simulation is very similar to an experimentally solved native conformation. Furthermore, the computation speed is four-times faster than that of our previous simulation engine, myPresto/psygene-G. The new simulator, myPresto/omegagene, is freely available at the following URLs: http://www.protein.osaka-u.ac.jp/rcsfp/pi/omegagene/ and http://presto.protein.osaka-u.ac.jp/myPresto4/.
Collapse
Affiliation(s)
- Kota Kasahara
- College of Life Sciences, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan
| | - Benson Ma
- College of Engineering, University of Illinois, Urbana-Champaign, United States
| | - Kota Goto
- School of Computing, Tokyo Institute of Technology, Tokyo 152-8550, Japan
| | - Bhaskar Dasgupta
- Institute for Protein Research, Osaka University, Suita, Osaka 565-0871, Japan; Technology Research Association for Next Generation Natural Products Chemistry, Tokyo 135-0064, Japan
| | - Junichi Higo
- Institute for Protein Research, Osaka University, Suita, Osaka 565-0871, Japan
| | - Ikuo Fukuda
- Institute for Protein Research, Osaka University, Suita, Osaka 565-0871, Japan
| | - Tadaaki Mashimo
- Technology Research Association for Next Generation Natural Products Chemistry, Tokyo 135-0064, Japan
| | - Yutaka Akiyama
- School of Computing, Tokyo Institute of Technology, Tokyo 152-8550, Japan; Molecular Profiling Research Center for Drug Discovery, National Institute of Advanced Industrial Science and Technology, Tokyo 135-0064, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, Suita, Osaka 565-0871, Japan
| |
Collapse
|
25
|
Bekker GJ, Nakamura H, Kinjo AR. Molmil: a molecular viewer for the PDB and beyond. J Cheminform 2016; 8:42. [PMID: 27570544 PMCID: PMC5002144 DOI: 10.1186/s13321-016-0155-1] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2016] [Accepted: 08/19/2016] [Indexed: 11/17/2022] Open
Abstract
We have developed a new platform-independent web-based molecular viewer using JavaScript and WebGL. The molecular viewer, Molmil, has been integrated into several services offered by Protein Data Bank Japan and can be easily extended with new functionality by third party developers. Furthermore, the viewer can be used to load files in various formats from the user’s local hard drive without uploading the data to a server. Molmil is available for all platforms supporting WebGL (e.g. Windows, Linux, iOS, Android) from http://gjbekker.github.io/molmil/. The source code is available at http://github.com/gjbekker/molmil under the LGPLv3 licence.
Collapse
Affiliation(s)
- Gert-Jan Bekker
- Laboratory of Protein Informatics, Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871 Japan ; Graduate School of Frontier Biosciences, Osaka University, Suita, Osaka 565-0871 Japan
| | - Haruki Nakamura
- Laboratory of Protein Informatics, Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871 Japan
| | - Akira R Kinjo
- Laboratory of Protein Informatics, Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871 Japan
| |
Collapse
|
26
|
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
|
27
|
Harada R, Takano Y, Shigeta Y. TaBoo SeArch Algorithm with a Modified Inverse Histogram for Reproducing Biologically Relevant Rare Events of Proteins. J Chem Theory Comput 2016; 12:2436-45. [DOI: 10.1021/acs.jctc.6b00082] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Ryuhei Harada
- Department
of Physics, Graduate School of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8571, Japan
- Center
for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan
- Computational
Engineering Applications Unit, RIKEN, Advanced Center for Computing and Communication, 2-1, Hirosawa, Wako-shi, Saitama 351-0198, Japan
| | - Yu Takano
- Department
of Biomedical Information Sciences, Hiroshima City University, 3-4-1
Ozuka-Higashi, Asa-Minami-Ku, Hiroshima 731-3194, Japan
| | - Yasuteru Shigeta
- Department
of Physics, Graduate School of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8571, Japan
- Center
for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan
| |
Collapse
|
28
|
Kasahara K, Mohan N, Fukuda I, Nakamura H. mDCC_tools: characterizing multi-modal atomic motions in molecular dynamics trajectories. ACTA ACUST UNITED AC 2016; 32:2531-3. [PMID: 27153575 PMCID: PMC4978922 DOI: 10.1093/bioinformatics/btw129] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2015] [Accepted: 03/03/2016] [Indexed: 11/14/2022]
Abstract
UNLABELLED We previously reported the multi-modal Dynamic Cross Correlation (mDCC) method for analyzing molecular dynamics trajectories. This method quantifies the correlation coefficients of atomic motions with complex multi-modal behaviors by using a Bayesian-based pattern recognition technique that can effectively capture transiently formed, unstable interactions. Here, we present an open source toolkit for performing the mDCC analysis, including pattern recognitions, complex network analyses and visualizations. We include a tutorial document that thoroughly explains how to apply this toolkit for an analysis, using the example trajectory of the 100 ns simulation of an engineered endothelin-1 peptide dimer. AVAILABILITY AND IMPLEMENTATION The source code is available for free at http://www.protein.osaka-u.ac.jp/rcsfp/pi/mdcctools/, implemented in C ++ and Python, and supported on Linux. CONTACT kota.kasahara@protein.osaka-u.ac.jp SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Kota Kasahara
- Institute for Protein Research, Osaka University, Suita, Osaka 565-0871, Japan
| | - Neetha Mohan
- Institute for Protein Research, Osaka University, Suita, Osaka 565-0871, Japan
| | - Ikuo Fukuda
- 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
|
29
|
Wang H, Nakamura H, Fukuda I. A critical appraisal of the zero-multipole method: Structural, thermodynamic, dielectric, and dynamical properties of a water system. J Chem Phys 2016; 144:114503. [DOI: 10.1063/1.4943956] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Affiliation(s)
- Han Wang
- CAEP Software Center for High Performance Numerical Simulation, Huayuan Road 6, 100088 Beijing, China and Zuse Institute Berlin (ZIB), Berlin, Germany
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Ikuo Fukuda
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
- RIKEN (The Institute of Physical and Chemical Research), 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| |
Collapse
|
30
|
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
|
31
|
Harada R, Takano Y, Baba T, Shigeta Y. Simple, yet powerful methodologies for conformational sampling of proteins. Phys Chem Chem Phys 2016; 17:6155-73. [PMID: 25659594 DOI: 10.1039/c4cp05262e] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Several biological functions, such as molecular recognition, enzyme catalysis, signal transduction, allosteric regulation, and protein folding, are strongly related to conformational transitions of proteins. These conformational transitions are generally induced as slow dynamics upon collective motions, including biologically relevant large-amplitude fluctuations of proteins. Although molecular dynamics (MD) simulation has become a powerful tool for extracting conformational transitions of proteins, it might still be difficult to reach time scales of the biological functions because the accessible time scales of MD simulations are far from biological time scales, even if straightforward conventional MD (CMD) simulations using massively parallel computers are employed. Thus, it is desirable to develop efficient methods to achieve canonical ensembles with low computational costs. From this perspective, we review several enhanced conformational sampling techniques of biomolecules developed by us. In our methods, multiple independent short-time MD simulations are employed instead of single straightforward long-time CMD simulations. Our basic strategy is as follows: (i) selection of initial seeds (initial structures) for the conformational sampling in restarting MD simulations. Here, the seeds should be selected as candidates with high potential to transit. (ii) Resampling from the selected seeds by initializing velocities in restarting short-time MD simulations. A cycle of these simple protocols might drastically promote the conformational transitions of biomolecules. (iii) Once reactive trajectories extracted from the cycles of short-time MD simulations are obtained, a free energy profile is evaluated by means of umbrella sampling (US) techniques with the weighted histogram analysis method (WHAM) as a post-processing technique. For the selection of the initial seeds, we proposed four different choices: (1) Parallel CaScade molecular dynamics (PaCS-MD), (2) Fluctuation Flooding Method (FFM), (3) Outlier FLOODing (OFLOOD) method, and (4) TaBoo SeArch (TBSA) method. We demonstrate applications of our methods to several biological systems, such as domain motions of proteins with large-amplitude fluctuations, conformational transitions upon ligand binding, and protein folding/refolding to native structures of proteins. Finally, we show the conformational sampling efficiencies of our methods compared with those by CMD simulations and other previously developed enhanced conformational sampling methods.
Collapse
Affiliation(s)
- Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8571, Japan.
| | | | | | | |
Collapse
|
32
|
Fukunishi Y, Mashimo T, Misoo K, Wakabayashi Y, Miyaki T, Ohta S, Nakamura M, Ikeda K. Miscellaneous Topics in Computer-Aided Drug Design: Synthetic Accessibility and GPU Computing, and Other Topics. Curr Pharm Des 2016; 22:3555-68. [PMID: 27075578 PMCID: PMC5080912 DOI: 10.2174/1381612822666160414142547] [Citation(s) in RCA: 6] [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: 03/14/2016] [Accepted: 04/12/2016] [Indexed: 01/12/2023]
Abstract
BACKGROUND Computer-aided drug design is still a state-of-the-art process in medicinal chemistry, and the main topics in this field have been extensively studied and well reviewed. These topics include compound databases, ligand-binding pocket prediction, protein-compound docking, virtual screening, target/off-target prediction, physical property prediction, molecular simulation and pharmacokinetics/pharmacodynamics (PK/PD) prediction. Message and Conclusion: However, there are also a number of secondary or miscellaneous topics that have been less well covered. For example, methods for synthesizing and predicting the synthetic accessibility (SA) of designed compounds are important in practical drug development, and hardware/software resources for performing the computations in computer-aided drug design are crucial. Cloud computing and general purpose graphics processing unit (GPGPU) computing have been used in virtual screening and molecular dynamics simulations. Not surprisingly, there is a growing demand for computer systems which combine these resources. In the present review, we summarize and discuss these various topics of drug design.
Collapse
Affiliation(s)
- 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.
| | | | | | | | | | | | | | | |
Collapse
|
33
|
Advances in Computational Techniques to Study GPCR–Ligand Recognition. Trends Pharmacol Sci 2015; 36:878-890. [DOI: 10.1016/j.tips.2015.08.006] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2015] [Revised: 07/15/2015] [Accepted: 08/07/2015] [Indexed: 12/16/2022]
|
34
|
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
|
35
|
Shimoda T, Suzuki S, Ohue M, Ishida T, Akiyama Y. Protein-protein docking on hardware accelerators: comparison of GPU and MIC architectures. BMC SYSTEMS BIOLOGY 2015; 9 Suppl 1:S6. [PMID: 25707855 PMCID: PMC4331681 DOI: 10.1186/1752-0509-9-s1-s6] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Background The hardware accelerators will provide solutions to computationally complex problems in bioinformatics fields. However, the effect of acceleration depends on the nature of the application, thus selection of an appropriate accelerator requires some consideration. Results In the present study, we compared the effects of acceleration using graphics processing unit (GPU) and many integrated core (MIC) on the speed of fast Fourier transform (FFT)-based protein-protein docking calculation. The GPU implementation performed the protein-protein docking calculations approximately five times faster than the MIC offload mode implementation. The MIC native mode implementation has the advantage in the implementation costs. However, the performance was worse with larger protein pairs because of memory limitations. Conclusion The results suggest that GPU is more suitable than MIC for accelerating FFT-based protein-protein docking applications.
Collapse
|
36
|
Kasahara K, Fukuda I, Nakamura H. A novel approach of dynamic cross correlation analysis on molecular dynamics simulations and its application to Ets1 dimer-DNA complex. PLoS One 2014; 9:e112419. [PMID: 25380315 PMCID: PMC4224484 DOI: 10.1371/journal.pone.0112419] [Citation(s) in RCA: 79] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2014] [Accepted: 10/14/2014] [Indexed: 11/18/2022] Open
Abstract
The dynamic cross correlation (DCC) analysis is a popular method for analyzing the trajectories of molecular dynamics (MD) simulations. However, it is difficult to detect correlative motions that appear transiently in only a part of the trajectory, such as atomic contacts between the side-chains of amino acids, which may rapidly flip. In order to capture these multi-modal behaviors of atoms, which often play essential roles, particularly at the interfaces of macromolecules, we have developed the "multi-modal DCC (mDCC)" analysis. The mDCC is an extension of the DCC and it takes advantage of a Bayesian-based pattern recognition technique. We performed MD simulations for molecular systems modeled from the (Ets1)2-DNA complex and analyzed their results with the mDCC method. Ets1 is an essential transcription factor for a variety of physiological processes, such as immunity and cancer development. Although many structural and biochemical studies have so far been performed, its DNA binding properties are still not well characterized. In particular, it is not straightforward to understand the molecular mechanisms how the cooperative binding of two Ets1 molecules facilitates their recognition of Stromelysin-1 gene regulatory elements. A correlation network was constructed among the essential atomic contacts, and the two major pathways by which the two Ets1 molecules communicate were identified. One is a pathway via direct protein-protein interactions and the other is that via the bound DNA intervening two recognition helices. These two pathways intersected at the particular cytosine bases (C110/C11), interacting with the H1, H2, and H3 helices. Furthermore, the mDCC analysis showed that both pathways included the transient interactions at their intermolecular interfaces of Tyr396-C11 and Ala327-Asn380 in multi-modal motions of the amino acid side chains and the nucleotide backbone. Thus, the current mDCC approach is a powerful tool to reveal these complicated behaviors and scrutinize intermolecular communications in a molecular system.
Collapse
Affiliation(s)
- Kota Kasahara
- Institute for Protein Research, Osaka University, Suita, Osaka, Japan
- * E-mail:
| | - Ikuo Fukuda
- Institute for Protein Research, Osaka University, Suita, Osaka, Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, Suita, Osaka, Japan
| |
Collapse
|
37
|
Fukuda I, Kamiya N, Nakamura H. The zero-multipole summation method for estimating electrostatic interactions in molecular dynamics: analysis of the accuracy and application to liquid systems. J Chem Phys 2014; 140:194307. [PMID: 24852538 DOI: 10.1063/1.4875693] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
In the preceding paper [I. Fukuda, J. Chem. Phys. 139, 174107 (2013)], the zero-multipole (ZM) summation method was proposed for efficiently evaluating the electrostatic Coulombic interactions of a classical point charge system. The summation takes a simple pairwise form, but prevents the electrically non-neutral multipole states that may artificially be generated by a simple cutoff truncation, which often causes large energetic noises and significant artifacts. The purpose of this paper is to judge the ability of the ZM method by investigating the accuracy, parameter dependencies, and stability in applications to liquid systems. To conduct this, first, the energy-functional error was divided into three terms and each term was analyzed by a theoretical error-bound estimation. This estimation gave us a clear basis of the discussions on the numerical investigations. It also gave a new viewpoint between the excess energy error and the damping effect by the damping parameter. Second, with the aid of these analyses, the ZM method was evaluated based on molecular dynamics (MD) simulations of two fundamental liquid systems, a molten sodium-chlorine ion system and a pure water molecule system. In the ion system, the energy accuracy, compared with the Ewald summation, was better for a larger value of multipole moment l currently induced until l ≲ 3 on average. This accuracy improvement with increasing l is due to the enhancement of the excess-energy accuracy. However, this improvement is wholly effective in the total accuracy if the theoretical moment l is smaller than or equal to a system intrinsic moment L. The simulation results thus indicate L ∼ 3 in this system, and we observed less accuracy in l = 4. We demonstrated the origins of parameter dependencies appearing in the crossing behavior and the oscillations of the energy error curves. With raising the moment l we observed, smaller values of the damping parameter provided more accurate results and smoother behaviors with respect to cutoff length were obtained. These features can be explained, on the basis of the theoretical error analyses, such that the excess energy accuracy is improved with increasing l and that the total accuracy improvement within l ⩽ L is facilitated by a small damping parameter. Although the accuracy was fundamentally similar to the ion system, the bulk water system exhibited distinguishable quantitative behaviors. A smaller damping parameter was effective in all the practical cutoff distance, and this fact can be interpreted by the reduction of the excess subset. A lower moment was advantageous in the energy accuracy, where l = 1 was slightly superior to l = 2 in this system. However, the method with l = 2 (viz., the zero-quadrupole sum) gave accurate results for the radial distribution function. We confirmed the stability in the numerical integration for MD simulations employing the ZM scheme. This result is supported by the sufficient smoothness of the energy function. Along with the smoothness, the pairwise feature and the allowance of the atom-based cutoff mode on the energy formula lead to the exact zero total-force, ensuring the total-momentum conservations for typical MD equations of motion.
Collapse
Affiliation(s)
- Ikuo Fukuda
- 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
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| |
Collapse
|
38
|
Viveros-Méndez PX, Gil-Villegas A, Aranda-Espinoza S. Monte Carlo computer simulation of sedimentation of charged hard spherocylinders. J Chem Phys 2014; 141:044905. [PMID: 25084954 DOI: 10.1063/1.4890819] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this article we present a NVT Monte Carlo computer simulation study of sedimentation of an electroneutral mixture of oppositely charged hard spherocylinders (CHSC) with aspect ratio L/σ = 5, where L and σ are the length and diameter of the cylinder and hemispherical caps, respectively, for each particle. This system is an extension of the restricted primitive model for spherical particles, where L/σ = 0, and it is assumed that the ions are immersed in an structureless solvent, i.e., a continuum with dielectric constant D. The system consisted of N = 2000 particles and the Wolf method was implemented to handle the coulombic interactions of the inhomogeneous system. Results are presented for different values of the strength ratio between the gravitational and electrostatic interactions, Γ = (mgσ)/(e(2)/Dσ), where m is the mass per particle, e is the electron's charge and g is the gravitational acceleration value. A semi-infinite simulation cell was used with dimensions Lx ≈ Ly and Lz = 5Lx, where Lx, Ly, and Lz are the box dimensions in Cartesian coordinates, and the gravitational force acts along the z-direction. Sedimentation effects were studied by looking at every layer formed by the CHSC along the gravitational field. By increasing Γ, particles tend to get more packed at each layer and to arrange in local domains with an orientational ordering along two perpendicular axis, a feature not observed in the uncharged system with the same hard-body geometry. This type of arrangement, known as tetratic phase, has been observed in two-dimensional systems of hard-rectangles and rounded hard-squares. In this way, the coupling of gravitational and electric interactions in the CHSC system induces the arrangement of particles in layers, with the formation of quasi-two dimensional tetratic phases near the surface.
Collapse
Affiliation(s)
- P X Viveros-Méndez
- Unidad Académica de Física, Universidad Autónoma de Zacatecas, Calzada Solidaridad esq. Paseo, La Bufa s/n, 98060 Zacatecas, Zacatecas, México
| | - Alejandro Gil-Villegas
- Departamento de Ingeniería Física, División de Ciencias e Ingenierías, Campus León, Universidad de Guanajuato, Loma del Bosque 103, Lomas del Campestre, 37150 León, Guanajuato, México
| | - S Aranda-Espinoza
- Unidad Académica de Física, Universidad Autónoma de Zacatecas, Calzada Solidaridad esq. Paseo, La Bufa s/n, 98060 Zacatecas, Zacatecas, México
| |
Collapse
|
39
|
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: 21] [Impact Index Per Article: 2.1] [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
|