1
|
Kuroda M, Kasahara Y, Hirose M, Yamaguma H, Oda M, Nagao C, Mizuguchi K. Construction of a T m-value prediction model and molecular dynamics study of AmNA-containing gapmer antisense oligonucleotide. MOLECULAR THERAPY. NUCLEIC ACIDS 2024; 35:102272. [PMID: 39176173 PMCID: PMC11339022 DOI: 10.1016/j.omtn.2024.102272] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/11/2024] [Accepted: 07/12/2024] [Indexed: 08/24/2024]
Abstract
RNase H-dependent antisense oligonucleotides (gapmer ASOs) represent a class of nucleic acid therapeutics that bind to target RNA to facilitate RNase H-mediated RNA cleavage, thereby regulating the expression of disease-associated proteins. Integrating artificial nucleic acids into gapmer ASOs enhances their therapeutic efficacy. Among these, amido-bridged nucleic acid (AmNA) stands out for its potential to confer high affinity and stability to ASOs. However, a significant challenge in the design of gapmer ASOs incorporating artificial nucleic acids, such as AmNA, is the accurate prediction of their melting temperature (T m ) values. The T m is a critical parameter for designing effective gapmer ASOs to ensure proper functioning. However, predicting accurate T m values for oligonucleotides containing artificial nucleic acids remains problematic. We developed a T m prediction model using a library of AmNA-containing ASOs to address this issue. We measured the T m values of 157 oligonucleotides through differential scanning calorimetry, enabling the construction of an accurate prediction model. Additionally, molecular dynamics simulations were used to elucidate the molecular mechanisms by which AmNA modifications elevate T m , thereby informing the design strategies of gapmer ASOs.
Collapse
Affiliation(s)
- Masataka Kuroda
- National Institutes of Biomedical Innovation, Health and Nutrition (NIBIOHN), Osaka 566-0002, Japan
- Mitsubishi Tanabe Pharma Corporation, Yokohama 227-0033, Japan
| | - Yuuya Kasahara
- National Institutes of Biomedical Innovation, Health and Nutrition (NIBIOHN), Osaka 566-0002, Japan
- Graduate School of Pharmaceutical Science, Osaka University, Osaka 565-0871, Japan
| | - Masako Hirose
- Malvern Panalytical, Spectris, Tokyo 105-0013, Japan
| | - Harumi Yamaguma
- National Institutes of Biomedical Innovation, Health and Nutrition (NIBIOHN), Osaka 566-0002, Japan
| | - Masayuki Oda
- Graduate School of Life and Environmental Sciences, Kyoto Prefectural University, Kyoto 606-8522, Japan
| | - Chioko Nagao
- Institute for Protein Research, Osaka University, Osaka 565-0871, Japan
| | - Kenji Mizuguchi
- National Institutes of Biomedical Innovation, Health and Nutrition (NIBIOHN), Osaka 566-0002, Japan
- Institute for Protein Research, Osaka University, Osaka 565-0871, Japan
| |
Collapse
|
2
|
Matsubara Y, Okabe R, Masayama R, Watanabe NM, Umakoshi H, Kasahara K, Matubayasi N. A methodology of quantifying membrane permeability based on returning probability theory and molecular dynamics simulation. J Chem Phys 2024; 161:024108. [PMID: 38984955 DOI: 10.1063/5.0214401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2024] [Accepted: 06/14/2024] [Indexed: 07/11/2024] Open
Abstract
We propose a theoretical approach to estimate the permeability coefficients of substrates (permeants) for crossing membranes from donor (D) phase to acceptor (A) phase by means of molecular dynamics (MD) simulation. A fundamental aspect of our approach involves reformulating the returning probability (RP) theory, a rigorous bimolecular reaction theory, to describe permeation phenomena. This reformulation relies on the parallelism between permeation and bimolecular reaction processes. In the present method, the permeability coefficient is represented in terms of the thermodynamic and kinetic quantities for the reactive (R) phase that exists within the inner region of a membrane. One can evaluate these quantities using multiple MD trajectories starting from phase R. We apply the RP theory to the permeation of ethanol and methylamine at different concentrations (infinitely dilute and 1 mol % conditions of permeants). Under the 1 mol% condition, the present method yields a larger permeability coefficient for ethanol (0.12 ± 0.01 cm s-1) than for methylamine (0.069 ± 0.006 cm s-1), while the values of the permeability coefficient are satisfactorily close to those obtained from the brute-force MD simulations (0.18 ± 0.03 and 0.052 ± 0.005 cm s-1 for ethanol and methylamine, respectively). Moreover, upon analyzing the thermodynamic and kinetic contributions to the permeability, we clarify that a higher concentration dependency of permeability for ethanol, as compared to methylamine, arises from the sensitive nature of ethanol's free-energy barrier within the inner region of the membrane against ethanol concentration.
Collapse
Affiliation(s)
- Yuya Matsubara
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Ryo Okabe
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Ren Masayama
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nozomi Morishita Watanabe
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Hiroshi Umakoshi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Kento Kasahara
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
3
|
Jung J, Yagi K, Tan C, Oshima H, Mori T, Yu I, Matsunaga Y, Kobayashi C, Ito S, Ugarte La Torre D, Sugita Y. GENESIS 2.1: High-Performance Molecular Dynamics Software for Enhanced Sampling and Free-Energy Calculations for Atomistic, Coarse-Grained, and Quantum Mechanics/Molecular Mechanics Models. J Phys Chem B 2024; 128:6028-6048. [PMID: 38876465 PMCID: PMC11215777 DOI: 10.1021/acs.jpcb.4c02096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2024] [Revised: 05/15/2024] [Accepted: 05/21/2024] [Indexed: 06/16/2024]
Abstract
GENeralized-Ensemble SImulation System (GENESIS) is a molecular dynamics (MD) software developed to simulate the conformational dynamics of a single biomolecule, as well as molecular interactions in large biomolecular assemblies and between multiple biomolecules in cellular environments. To achieve the latter purpose, the earlier versions of GENESIS emphasized high performance in atomistic MD simulations on massively parallel supercomputers, with or without graphics processing units (GPUs). Here, we implemented multiscale MD simulations that include atomistic, coarse-grained, and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. They demonstrate high performance and are integrated with enhanced conformational sampling algorithms and free-energy calculations without using external programs except for the QM programs. In this article, we review new functions, molecular models, and other essential features in GENESIS version 2.1 and discuss ongoing developments for future releases.
Collapse
Affiliation(s)
- Jaewoon Jung
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
| | - Kiyoshi Yagi
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
| | - Cheng Tan
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Hiraku Oshima
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan
- Graduate
School of Life Science, University of Hyogo, Harima Science Park City, Hyogo 678-1297, Japan
| | - Takaharu Mori
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
- Department
of Chemistry, Tokyo University of Science, Shinjuku-ku, Tokyo 162-8601, Japan
| | - Isseki Yu
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
- Department
of Bioinformatics, Maebashi Institute of
Technology, Maebashi, Gunma 371-0816, Japan
| | - Yasuhiro Matsunaga
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
- Graduate
School of Science and Engineering, Saitama
University, Saitama 338-8570, Japan
| | - Chigusa Kobayashi
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Shingo Ito
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
| | - Diego Ugarte La Torre
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
4
|
Lei YK, Yagi K, Sugita Y. Learning QM/MM potential using equivariant multiscale model. J Chem Phys 2024; 160:214109. [PMID: 38828815 DOI: 10.1063/5.0205123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2024] [Accepted: 05/09/2024] [Indexed: 06/05/2024] Open
Abstract
The machine learning (ML) method emerges as an efficient and precise surrogate model for high-level electronic structure theory. Its application has been limited to closed chemical systems without considering external potentials from the surrounding environment. To address this limitation and incorporate the influence of external potentials, polarization effects, and long-range interactions between a chemical system and its environment, the first two terms of the Taylor expansion of an electrostatic operator have been used as extra input to the existing ML model to represent the electrostatic environments. However, high-order electrostatic interaction is often essential to account for external potentials from the environment. The existing models based only on invariant features cannot capture significant distribution patterns of the external potentials. Here, we propose a novel ML model that includes high-order terms of the Taylor expansion of an electrostatic operator and uses an equivariant model, which can generate a high-order tensor covariant with rotations as a base model. Therefore, we can use the multipole-expansion equation to derive a useful representation by accounting for polarization and intermolecular interaction. Moreover, to deal with long-range interactions, we follow the same strategy adopted to derive long-range interactions between a target system and its environment media. Our model achieves higher prediction accuracy and transferability among various environment media with these modifications.
Collapse
Affiliation(s)
- Yao-Kun Lei
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama 351-0198, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
- RIKEN Interdisciplinary Theoretical and Mathematical Sciences Program (iTHEMS), Wako, Saitama 351-0198, Japan
| | - Kiyoshi Yagi
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama 351-0198, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama 351-0198, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
- RIKEN Interdisciplinary Theoretical and Mathematical Sciences Program (iTHEMS), Wako, Saitama 351-0198, Japan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
5
|
Ruedas R, Vuillemot R, Tubiana T, Winter JM, Pieri L, Arteni AA, Samson C, Jonic S, Mathieu M, Bressanelli S. Structure and conformational variability of the HER2-trastuzumab-pertuzumab complex. J Struct Biol 2024; 216:108095. [PMID: 38723875 DOI: 10.1016/j.jsb.2024.108095] [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: 12/02/2023] [Revised: 03/11/2024] [Accepted: 04/24/2024] [Indexed: 05/13/2024]
Abstract
Single particle analysis from cryogenic transmission electron microscopy (cryo-EM) is particularly attractive for complexes for which structure prediction remains intractable, such as antibody-antigen complexes. Here we obtain the detailed structure of a particularly difficult complex between human epidermal growth factor receptor 2 (HER2) and the antigen-binding fragments from two distinct therapeutic antibodies binding to distant parts of the flexible HER2, pertuzumab and trastuzumab (HTP). We highlight the strengths and limitations of current data processing software in dealing with various kinds of heterogeneities, particularly continuous conformational heterogeneity, and in describing the motions that can be extracted from our dataset. Our HTP structure provides a more detailed view than the one previously available for this ternary complex. This allowed us to pinpoint a previously overlooked loop in domain IV that may be involved both in binding of trastuzumab and in HER2 dimerization. This finding may contribute to explain the synergistic anticancer effect of the two antibodies. We further propose that the flexibility of the HTP complex, beyond the difficulties it causes for cryo-EM analysis, actually reflects regulation of HER2 signaling and its inhibition by therapeutic antibodies. Notably we obtain our best data with ultra-thin continuous carbon grids, showing that with current cameras their use to alleviate particle misdistribution is compatible with a protein complex of only 162 kDa. Perhaps most importantly, we provide here a dataset for such a smallish protein complex for further development of software accounting for continuous conformational heterogeneity in cryo-EM images.
Collapse
Affiliation(s)
- Rémi Ruedas
- Université Paris-Saclay, CEA, CNRS - Institute for Integrative Biology of the Cell (I2BC), 91198, Gif-sur-Yvette, France; Sanofi, Integrated Drug Discovery, 13, quai Jules Guesde 94403, Vitry-sur-Seine, France
| | - Rémi Vuillemot
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, 75005, Paris, France
| | - Thibault Tubiana
- Université Paris-Saclay, CEA, CNRS - Institute for Integrative Biology of the Cell (I2BC), 91198, Gif-sur-Yvette, France
| | - Jean-Marie Winter
- NanoImaging Core Facility, Centre de Ressources et Recherches Technologiques (C2RT), Institut Pasteur, 75015, Paris, France
| | - Laura Pieri
- Université Paris-Saclay, CEA, CNRS - Institute for Integrative Biology of the Cell (I2BC), 91198, Gif-sur-Yvette, France
| | - Ana-Andreea Arteni
- Université Paris-Saclay, CEA, CNRS - Institute for Integrative Biology of the Cell (I2BC), 91198, Gif-sur-Yvette, France
| | - Camille Samson
- Sanofi, Integrated Drug Discovery, 13, quai Jules Guesde 94403, Vitry-sur-Seine, France
| | - Slavica Jonic
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, 75005, Paris, France
| | - Magali Mathieu
- Sanofi, Integrated Drug Discovery, 13, quai Jules Guesde 94403, Vitry-sur-Seine, France
| | - Stéphane Bressanelli
- Université Paris-Saclay, CEA, CNRS - Institute for Integrative Biology of the Cell (I2BC), 91198, Gif-sur-Yvette, France.
| |
Collapse
|
6
|
Díaz Mirón G, Lien-Medrano CR, Banerjee D, Morzan UN, Sentef MA, Gebauer R, Hassanali A. Exploring the Mechanisms behind Non-aromatic Fluorescence with the Density Functional Tight Binding Method. J Chem Theory Comput 2024; 20:3864-3878. [PMID: 38634760 DOI: 10.1021/acs.jctc.4c00125] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/19/2024]
Abstract
Recent experimental findings reveal nonconventional fluorescence emission in biological systems devoid of conjugated bonds or aromatic compounds, termed non-aromatic fluorescence (NAF). This phenomenon is exclusive to aggregated or solid states and remains absent in monomeric solutions. Previous studies focused on small model systems in vacuum show that the carbonyl stretching mode along with strong interaction of short hydrogen bonds (SHBs) remains the primary vibrational mode explaining NAF in these systems. In order to simulate larger model systems taking into account the effects of the surrounding environment, in this work we propose using the density functional tight-binding (DFTB) method in combination with non-adiabatic molecular dynamics (NAMD) and the mixed quantum/molecular mechanics (QM/MM) approach. We investigate the mechanism behind NAF in the crystal structure of l-pyroglutamine-ammonium, comparing it with the related nonfluorescent amino acid l-glutamine. Our results extend our previous findings to more realistic systems, demonstrating the efficiency and robustness of the proposed DFTB method in the context of NAMD in biological systems. Furthermore, due to its inherent low computational cost, this method allows for a better sampling of the nonradiative events at the conical intersection which is crucial for a complete understanding of this phenomenon. Beyond contributing to the ongoing exploration of NAF, this work paves the way for future application of this method in more complex biological systems such as amyloid aggregates, biomaterials, and non-aromatic proteins.
Collapse
Affiliation(s)
- Gonzalo Díaz Mirón
- Condensed Matter and Statistical Physics, The Abdus Salam International Centre for Theoretical Physics, 34151 Trieste, Italy
| | - Carlos R Lien-Medrano
- Institute for Theoretical Physics and Bremen Center for Computational Materials Science, University of Bremen, 28359 Bremen, Germany
| | - Debarshi Banerjee
- Condensed Matter and Statistical Physics, The Abdus Salam International Centre for Theoretical Physics, 34151 Trieste, Italy
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), 34136 Trieste, Italy
| | - Uriel N Morzan
- Instituto de Fisica de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, C1428EGA Buenos Aires, Argentina
| | - Michael A Sentef
- Institute for Theoretical Physics and Bremen Center for Computational Materials Science, University of Bremen, 28359 Bremen, Germany
- Center for Free-Electron Laser Science (CFEL), Max Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
| | - Ralph Gebauer
- Condensed Matter and Statistical Physics, The Abdus Salam International Centre for Theoretical Physics, 34151 Trieste, Italy
| | - Ali Hassanali
- Condensed Matter and Statistical Physics, The Abdus Salam International Centre for Theoretical Physics, 34151 Trieste, Italy
| |
Collapse
|
7
|
Inoue M, Ekimoto T, Yamane T, Ikeguchi M. Computational Analysis of Activation of Dimerized Epidermal Growth Factor Receptor Kinase Using the String Method and Markov State Model. J Chem Inf Model 2024; 64:3884-3895. [PMID: 38670929 DOI: 10.1021/acs.jcim.4c00172] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/28/2024]
Abstract
Epidermal growth factor receptor (EGFR) activation is accompanied by dimerization. During the activation of the intracellular kinase domain, two EGFR kinases form an asymmetric dimer, and one side of the dimer (receiver) is activated. Using the string method and Markov state model (MSM), we performed a computational analysis of the structural changes in the activation of the EGFR dimer in this study. The string method reveals the minimum free-energy pathway (MFEP) from the inactive to active structure. The MSM was constructed from numerous trajectories of molecular dynamics simulations around the MFEP, which revealed the free-energy map of structural changes. In the activation of the receiver kinase, the unfolding of the activation loop (A-loop) is followed by the rearrangement of the C-helix, as observed in other kinases. However, unlike other kinases, the free-energy map of EGFR at the asymmetric dimer showed that the active state yielded the highest stability and revealed how interactions at the dimer interface induced receiver activation. As the H-helix of the activator approaches the C-helix of the receiver during activation, the A-loop unfolds. Subsequently, L782 of the receiver enters the pocket between the G- and H-helices of the activator, leading to a rearrangement of the hydrophobic residues around L782 of the receiver, which constitutes a structural rearrangement of the C-helix of the receiver from an outward to an inner position. The MSM analysis revealed long-time scale trajectories via kinetic Monte Carlo.
Collapse
Affiliation(s)
- Masao Inoue
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Toru Ekimoto
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Tsutomu Yamane
- HPC- and AI-driven Drug Development Platform Division, Center for Computational Science, RIKEN, 1-7-22, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Mitsunori Ikeguchi
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
- HPC- and AI-driven Drug Development Platform Division, Center for Computational Science, RIKEN, 1-7-22, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| |
Collapse
|
8
|
Okita K, Ito N, Morishita-Watanabe N, Umakoshi H, Kasahara K, Matubayasi N. Solvation dynamics on the diffusion timescale elucidated using energy-represented dynamics theory. Phys Chem Chem Phys 2024; 26:12852-12861. [PMID: 38623745 DOI: 10.1039/d4cp00235k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2024]
Abstract
Photoexcitation of a solute alters the solute-solvent interaction, resulting in the nonequilibrium relaxation of the solvation structure, often called a dynamic Stokes shift or solvation dynamics. Thanks to the local nature of the solute-solvent interaction, the characteristics of the local solvent environment dissolving the solute can be captured by the observation of this process. Recently, we derived the energy-represented Smoluchowski-Vlasov (ERSV) equation, a diffusion equation for molecular liquids, which can be used to analyze the solvation dynamics on the diffusion timescale. This equation expresses the time development for the solvent distribution on the solute-solvent pair interaction energy (energy coordinate). Since the energy coordinate can effectively treat the solvent flexibility in addition to the position and orientation, the ERSV equation can be utilized in various solvent systems. Here, we apply the ERSV equation to the solvation dynamics of 6-propionyl-2-dimethylamino naphthalene (Prodan) in water and different alcohol solvents (methanol, ethanol, and 1-propanol) for clarifying the differences of the relaxation processes among these solvents. Prodan is a solvent-sensitive fluorescent probe and is thus widely utilized for investigating heterogeneous environments. On the long timescale, the ERSV equation satisfactorily reproduces the relaxation time correlation functions obtained from the molecular dynamics (MD) simulations for these solvents. We reveal that the relaxation time coefficient on the diffusion timescale linearly correlates with the inverse of the translational diffusion coefficients for the alcohol solvents because of the Prodan-solvent energy distributions among the alcohols. In the case of water, the time coefficient deviates from the linear relationship for the alcohols due to the difference in the extent of importance of the collective motion between the water and alcohol solvents.
Collapse
Affiliation(s)
- Kazuya Okita
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
| | - Natsuumi Ito
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
| | - Nozomi Morishita-Watanabe
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
| | - Hiroshi Umakoshi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
| | - Kento Kasahara
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
| |
Collapse
|
9
|
Jung J, Tan C, Sugita Y. GENESIS CGDYN: large-scale coarse-grained MD simulation with dynamic load balancing for heterogeneous biomolecular systems. Nat Commun 2024; 15:3370. [PMID: 38643169 PMCID: PMC11032353 DOI: 10.1038/s41467-024-47654-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2023] [Accepted: 04/08/2024] [Indexed: 04/22/2024] Open
Abstract
Residue-level coarse-grained (CG) molecular dynamics (MD) simulation is widely used to investigate slow biological processes that involve multiple proteins, nucleic acids, and their complexes. Biomolecules in a large simulation system are distributed non-uniformly, limiting computational efficiency with conventional methods. Here, we develop a hierarchical domain decomposition scheme with dynamic load balancing for heterogeneous biomolecular systems to keep computational efficiency even after drastic changes in particle distribution. These schemes are applied to the dynamics of intrinsically disordered protein (IDP) droplets. During the fusion of two droplets, we find that the changes in droplet shape correlate with the mixing of IDP chains. Additionally, we simulate large systems with multiple IDP droplets, achieving simulation sizes comparable to those observed in microscopy. In our MD simulations, we directly observe Ostwald ripening, a phenomenon where small droplets dissolve and their molecules redeposit into larger droplets. These methods have been implemented in CGDYN of the GENESIS software, offering a tool for investigating mesoscopic biological processes using the residue-level CG models.
Collapse
Affiliation(s)
- Jaewoon Jung
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, 650-0047, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, 351-0198, Japan
| | - Cheng Tan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, 650-0047, Japan
| | - Yuji Sugita
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, 650-0047, Japan.
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, 351-0198, Japan.
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo, 650-0047, Japan.
| |
Collapse
|
10
|
Fujinami D, Mizui A, Miyata A, Ito S. In Vitro Characterization of an O-Specific Glycosyltransferase Involved in Flagellin Glycosylation. ACS Chem Biol 2024; 19:992-998. [PMID: 38562012 DOI: 10.1021/acschembio.4c00045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Glycosyltransferases play a fundamental role in the biosynthesis of glycoproteins and glycotherapeutics. In this study, we investigated protein glycosyltransferase FlgGT1, belonging to the GT2 family. The GT2 family includes cysteine S-glycosyltransferases involved in antimicrobial peptide biosyntheses, sharing conserved catalytic domains while exhibiting diverse C-terminal domains. Our in vitro studies revealed that FlgGT1 recognizes structural motifs rather than specific amino acid sequences when glycosylating the flagellin protein Hag. Notably, FlgGT1 is selective for serine or threonine O-glycosylation over cysteine S-glycosylation. Molecular dynamics simulations provided insights into the structural basis of FlgGT1's ability to accommodate various sugar nucleotides as donor substrates. Mutagenesis experiments on FlgGT1 demonstrated that truncating the relatively large C-terminal domain resulted in a loss of flagellin glycosylation activity. Our classification based on sequence similarity network analysis and AlphaFold2 structural predictions suggests that the acquisition of the C-terminal domain is a key evolutionary adaptation conferring distinct substrate specificities on glycosyltransferases within the GT2 family.
Collapse
Affiliation(s)
- Daisuke Fujinami
- Graduate Division of Nutritional and Environmental Sciences, University of Shizuoka, 52-1 Yada, Suruga-ku, Shizuoka 422-8526, Japan
| | - Ayumi Mizui
- Graduate Division of Nutritional and Environmental Sciences, University of Shizuoka, 52-1 Yada, Suruga-ku, Shizuoka 422-8526, Japan
| | - Azusa Miyata
- Graduate Division of Nutritional and Environmental Sciences, University of Shizuoka, 52-1 Yada, Suruga-ku, Shizuoka 422-8526, Japan
| | - Sohei Ito
- Graduate Division of Nutritional and Environmental Sciences, University of Shizuoka, 52-1 Yada, Suruga-ku, Shizuoka 422-8526, Japan
| |
Collapse
|
11
|
Watanabe K, Zhao Q, Iwatsuki R, Fukui R, Ren W, Sugita Y, Nishida N. Deciphering the Multi-state Conformational Equilibrium of HDM2 in the Regulation of p53 Binding: Perspectives from Molecular Dynamics Simulation and NMR Analysis. J Am Chem Soc 2024; 146:9790-9800. [PMID: 38549219 DOI: 10.1021/jacs.3c14383] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024]
Abstract
HDM2 negatively regulates the activity of the tumor suppressor p53. Previous NMR studies have shown that apo-HDM2 interconverts between an "open" state in which the N-terminal "lid" is disordered and a "closed" state in which the lid covers the p53-binding site in the core region. Molecular dynamics (MD) simulation studies have been performed to elucidate the conformational dynamics of HDM2, but the direct relevance of the experimental and computational analyses is unclear. In addition, how the phosphorylation of S17 in the lid contributes to the inhibition of p53 binding remains controversial. Here, we used both NMR and MD simulations to investigate the conformational dynamics of apo-HDM2. The NMR analysis revealed that apo-HDM2 exists in a fast-exchanging equilibrium within two closed states, closed 1 and closed 2, in addition to a previously demonstrated slow-exchanging "open-closed" equilibrium. MD simulations visualized two characteristic closed states, where the spatial orientation of the key residues corresponds well to the chemical shift changes of the NMR spectra. Furthermore, the phosphorylation of S17 induced an equilibrium shift toward closed 1, thereby suppressing the binding of p53 to HDM2. This study reveals a multi-state equilibrium of apo-HDM2 and provides new insights into the regulation mechanism of HDM2-p53 interactions.
Collapse
Affiliation(s)
- Kazuki Watanabe
- Graduate School of Pharmaceutical Sciences, Chiba University, 1-8-1 Inohana, Chuo-ku, Chiba 260-8675, Japan
| | - Qingci Zhao
- Graduate School of Pharmaceutical Sciences, Chiba University, 1-8-1 Inohana, Chuo-ku, Chiba 260-8675, Japan
| | - Ryosuke Iwatsuki
- Graduate School of Pharmaceutical Sciences, Chiba University, 1-8-1 Inohana, Chuo-ku, Chiba 260-8675, Japan
| | - Ryota Fukui
- Graduate School of Pharmaceutical Sciences, Chiba University, 1-8-1 Inohana, Chuo-ku, Chiba 260-8675, Japan
| | - Weitong Ren
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Hirosawa 2-1, Wako 351-0918, Saitama, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Hirosawa 2-1, Wako 351-0918, Saitama, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, 6-7-1 Minatojima-Minamimachi, Chuo-ku, Kobe 650-0047, Hyogo, Japan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, 6-7-1 Minatojima-Minamimachi, Chuo-ku, Kobe 650-0047, Hyogo, Japan
| | - Noritaka Nishida
- Graduate School of Pharmaceutical Sciences, Chiba University, 1-8-1 Inohana, Chuo-ku, Chiba 260-8675, Japan
| |
Collapse
|
12
|
Grassmann G, Miotto M, Desantis F, Di Rienzo L, Tartaglia GG, Pastore A, Ruocco G, Monti M, Milanetti E. Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments. Chem Rev 2024; 124:3932-3977. [PMID: 38535831 PMCID: PMC11009965 DOI: 10.1021/acs.chemrev.3c00550] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2023] [Revised: 02/20/2024] [Accepted: 02/21/2024] [Indexed: 04/11/2024]
Abstract
Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.
Collapse
Affiliation(s)
- Greta Grassmann
- Department
of Biochemical Sciences “Alessandro Rossi Fanelli”, Sapienza University of Rome, Rome 00185, Italy
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
| | - Mattia Miotto
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
| | - Fausta Desantis
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- The
Open University Affiliated Research Centre at Istituto Italiano di
Tecnologia, Genoa 16163, Italy
| | - Lorenzo Di Rienzo
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
| | - Gian Gaetano Tartaglia
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- Department
of Neuroscience and Brain Technologies, Istituto Italiano di Tecnologia, Genoa 16163, Italy
- Center
for Human Technologies, Genoa 16152, Italy
| | - Annalisa Pastore
- Experiment
Division, European Synchrotron Radiation
Facility, Grenoble 38043, France
| | - Giancarlo Ruocco
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- Department
of Physics, Sapienza University, Rome 00185, Italy
| | - Michele Monti
- RNA
System Biology Lab, Department of Neuroscience and Brain Technologies, Istituto Italiano di Tecnologia, Genoa 16163, Italy
| | - Edoardo Milanetti
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- Department
of Physics, Sapienza University, Rome 00185, Italy
| |
Collapse
|
13
|
Chyży P, Kulik M, Shinobu A, Re S, Sugita Y, Trylska J. Molecular dynamics in multidimensional space explains how mutations affect the association path of neomycin to a riboswitch. Proc Natl Acad Sci U S A 2024; 121:e2317197121. [PMID: 38579011 PMCID: PMC11009640 DOI: 10.1073/pnas.2317197121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2023] [Accepted: 02/15/2024] [Indexed: 04/07/2024] Open
Abstract
Riboswitches are messenger RNA (mRNA) fragments binding specific small molecules to regulate gene expression. A synthetic N1 riboswitch, inserted into yeast mRNA controls the translation of a reporter gene in response to neomycin. However, its regulatory activity is sensitive to single-point RNA mutations, even those distant from the neomycin binding site. While the association paths of neomycin to N1 and its variants remain unknown, recent fluorescence kinetic experiments indicate a two-step process driven by conformational selection. This raises the question of which step is affected by mutations. To address this, we performed all-atom two-dimensional replica-exchange molecular dynamics simulations for N1 and U14C, U14C[Formula: see text], U15A, and A17G mutants, ensuring extensive conformational sampling of both RNA and neomycin. The obtained neomycin association and binding paths, along with multidimensional free-energy profiles, revealed a two-step binding mechanism, consisting of conformational selection and induced fit. Neomycin binds to a preformed N1 conformation upon identifying a stable upper stem and U-turn motif in the riboswitch hairpin. However, the positioning of neomycin in the binding site occurs at different RNA-neomycin distances for each mutant, which may explain their different regulatory activities. The subsequent induced fit arises from the interactions of the neomycin's N3 amino group with RNA, causing the G9 backbone to rearrange. In the A17G mutant, the critical C6-A17/G17 stacking forms at a closer RNA-neomycin distance compared to N1. These findings together with estimated binding free energies coincide with experiments and elucidate why the A17G mutation decreases and U15A enhances N1 activity in response to neomycin.
Collapse
Affiliation(s)
- Piotr Chyży
- Centre of New Technologies, University of Warsaw, 02-097Warsaw, Poland
| | - Marta Kulik
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, 02-093Warsaw, Poland
| | - Ai Shinobu
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, 650-0047Kobe, Japan
| | - Suyong Re
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, 650-0047Kobe, Japan
- Artificial Intelligence Center for Health and Biomedical Research, National Institutes of Biomedical Innovation, Health, and Nutrition, 567-0085Ibaraki, Japan
| | - Yuji Sugita
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, 650-0047Kobe, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 351-0198Wako, Japan
- RIKEN Center for Computational Science, 650-0047Kobe, Japan
| | - Joanna Trylska
- Centre of New Technologies, University of Warsaw, 02-097Warsaw, Poland
| |
Collapse
|
14
|
Ito S, Sugita Y. Free-energy landscapes of transmembrane homodimers by bias-exchange adaptively biased molecular dynamics. Biophys Chem 2024; 307:107190. [PMID: 38290241 DOI: 10.1016/j.bpc.2024.107190] [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: 12/03/2023] [Revised: 01/21/2024] [Accepted: 01/21/2024] [Indexed: 02/01/2024]
Abstract
Membrane proteins play essential roles in various biological functions within the cell. One of the most common functional regulations involves the dimerization of two single-pass transmembrane (TM) helices. Glycophorin A (GpA) and amyloid precursor protein (APP) form TM homodimers in the membrane, which have been investigated both experimentally and computationally. The homodimer structures are well characterized using only four collective variables (CVs) when each TM helix is stable. The CVs are the interhelical distance, the crossing angle, and the Crick angles for two TM helices. However, conformational sampling with multi-dimensional replica-exchange umbrella sampling (REUS) requires too many replicas to sample all the CVs for exploring the conformational landscapes. Here, we show that the bias-exchange adaptively biased molecular dynamics (BE-ABMD) with the four CVs effectively explores the free-energy landscapes of the TM helix dimers of GpA, wild-type APP and its mutants in the IMM1 implicit membrane. Compared to the original ABMD, the bias-exchange algorithm in BE-ABMD can provide a more rapidly converged conformational landscape. The BE-ABMD simulations could also reveal TM packing interfaces of the membrane proteins and the dependence of the free-energy landscapes on the membrane thickness. This approach is valuable for numerous other applications, including those involving explicit solvent and a lipid bilayer in all-atom force fields or Martini coarse-grained models, and enhances our understanding of protein-protein interactions in biological membranes.
Collapse
Affiliation(s)
- Shingo Ito
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan; Computational Biophysics Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan; Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, 1-6-5 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.
| |
Collapse
|
15
|
Ikei M, Miyazaki R, Monden K, Naito Y, Takeuchi A, Takahashi YS, Tanaka Y, Murata K, Mori T, Ichikawa M, Tsukazaki T. YeeD is an essential partner for YeeE-mediated thiosulfate uptake in bacteria and regulates thiosulfate ion decomposition. PLoS Biol 2024; 22:e3002601. [PMID: 38656967 PMCID: PMC11073785 DOI: 10.1371/journal.pbio.3002601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2024] [Revised: 05/06/2024] [Accepted: 03/26/2024] [Indexed: 04/26/2024] Open
Abstract
Uptake of thiosulfate ions as an inorganic sulfur source from the environment is important for bacterial sulfur assimilation. Recently, a selective thiosulfate uptake pathway involving a membrane protein YeeE (TsuA) in Escherichia coli was characterized. YeeE-like proteins are conserved in some bacteria, archaea, and eukaryotes. However, the precise function of YeeE, along with its potential partner protein in the thiosulfate ion uptake pathway, remained unclear. Here, we assessed selective thiosulfate transport via Spirochaeta thermophila YeeE in vitro and characterized E. coli YeeD (TsuB) as an adjacent and essential protein for YeeE-mediated thiosulfate uptake in vivo. We further showed that S. thermophila YeeD possesses thiosulfate decomposition activity and that a conserved cysteine in YeeD was modified to several forms in the presence of thiosulfate. Finally, the crystal structures of S. thermophila YeeE-YeeD fusion proteins at 3.34-Å and 2.60-Å resolutions revealed their interactions. The association was evaluated by a binding assay using purified S. thermophila YeeE and YeeD. Based on these results, a model of the sophisticated uptake of thiosulfate ions by YeeE and YeeD is proposed.
Collapse
Affiliation(s)
- Mai Ikei
- Division of Biological Science, Graduate School of Science and Technology, Nara Institute of Science and Technology, Ikoma, Nara, Japan
| | - Ryoji Miyazaki
- Division of Biological Science, Graduate School of Science and Technology, Nara Institute of Science and Technology, Ikoma, Nara, Japan
| | - Keigo Monden
- Division of Biological Science, Graduate School of Science and Technology, Nara Institute of Science and Technology, Ikoma, Nara, Japan
| | - Yusuke Naito
- Division of Biological Science, Graduate School of Science and Technology, Nara Institute of Science and Technology, Ikoma, Nara, Japan
| | - Azusa Takeuchi
- Division of Biological Science, Graduate School of Science and Technology, Nara Institute of Science and Technology, Ikoma, Nara, Japan
| | - Yutaro S Takahashi
- Division of Biological Science, Graduate School of Science and Technology, Nara Institute of Science and Technology, Ikoma, Nara, Japan
| | - Yoshiki Tanaka
- Division of Biological Science, Graduate School of Science and Technology, Nara Institute of Science and Technology, Ikoma, Nara, Japan
| | - Keina Murata
- Department of Chemistry, Faculty of Science, Tokyo University of Science, Shinjuku-ku, Tokyo, Japan
| | - Takaharu Mori
- Department of Chemistry, Faculty of Science, Tokyo University of Science, Shinjuku-ku, Tokyo, Japan
| | - Muneyoshi Ichikawa
- State Key Laboratory of Genetic Engineering, Department of Biochemistry and Biophysics, School of Life Sciences, Fudan University, Shanghai, China
| | - Tomoya Tsukazaki
- Division of Biological Science, Graduate School of Science and Technology, Nara Institute of Science and Technology, Ikoma, Nara, Japan
| |
Collapse
|
16
|
Yu I, Mori T, Matsuoka D, Surblys D, Sugita Y. SPANA: Spatial decomposition analysis for cellular-scale molecular dynamics simulations. J Comput Chem 2024; 45:498-505. [PMID: 37966727 DOI: 10.1002/jcc.27260] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Revised: 10/31/2023] [Accepted: 11/02/2023] [Indexed: 11/16/2023]
Abstract
The rapid increase in computational power with the latest supercomputers has enabled atomistic molecular dynamics (MDs) simulations of biomolecules in biological membrane, cytoplasm, and other cellular environments. These environments often contain a million or more atoms to be simulated simultaneously. Therefore, their trajectory analyses involve heavy computations that can become a bottleneck in the computational studies. Spatial decomposition analysis (SPANA) is a set of analysis tools in the Generalized-Ensemble Simulation System (GENESIS) software package that can carry out MD trajectory analyses of large-scale biological simulations using multiple CPU cores in parallel. SPANA applies the spatial decomposition of a large biological system to distribute structural and dynamical analyses into individual CPU cores, which reduces the computational time and the memory size, significantly. SPANA opens new possibilities for detailed atomistic analyses of biomacromolecules as well as solvent water molecules, ions, and metabolites in MD simulation trajectories of very large biological systems containing more than millions of atoms in cellular environments.
Collapse
Affiliation(s)
- Isseki Yu
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan
- Department of Bioinformatics, Maebashi Institute of Technology, Maebashi, Gunma, Japan
| | - Takaharu Mori
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan
| | - Daisuke Matsuoka
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan
| | - Donatas Surblys
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
| |
Collapse
|
17
|
Ohno S, Manabe N, Uzawa J, Yamaguchi Y. Comparative Conformational Analysis of Acyclic Sugar Alcohols Ribitol, Xylitol and d-Arabitol by Solution NMR and Molecular Dynamics Simulations. Molecules 2024; 29:1072. [PMID: 38474585 DOI: 10.3390/molecules29051072] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2024] [Revised: 02/22/2024] [Accepted: 02/27/2024] [Indexed: 03/14/2024] Open
Abstract
Ribitol (C5H12O5) is an acyclic sugar alcohol that was recently identified in O-mannose glycan on mammalian α-dystroglycan. The conformation and dynamics of acyclic sugar alcohols such as ribitol are dependent on the stereochemistry of the hydroxyl groups; however, the dynamics are not fully understood. To gain insights into the conformation and dynamics of sugar alcohols, we carried out comparative analyses of ribitol, d-arabitol and xylitol by a crystal structure database search, solution NMR analysis and molecular dynamics (MD) simulations. The crystal structures of the sugar alcohols showed a limited number of conformations, suggesting that only certain stable conformations are prevalent among all possible conformations. The three-bond scholar coupling constants and exchange rates of hydroxyl protons were measured to obtain information on the backbone torsion angle and possible hydrogen bonding of each hydroxyl group. The 100 ns MD simulations indicate that the ribitol backbone has frequent conformational transitions with torsion angles between 180∘ and ±60∘, while d-arabitol and xylitol showed fewer conformational transitions. Taking our experimental and computational data together, it can be concluded that ribitol is more flexible than d-arabitol or xylitol, and the flexibility is at least in part defined by the configuration of the OH groups, which may form intramolecular hydrogen bonds.
Collapse
Affiliation(s)
- Shiho Ohno
- Division of Structural Biology, Institute of Molecular Biomembrane and Glycobiology, Tohoku Medical and Pharmaceutical University, 4-4-1 Komatsushima, Aoba-ku, Sendai 981-8558, Miyagi, Japan
| | - Noriyoshi Manabe
- Division of Structural Biology, Institute of Molecular Biomembrane and Glycobiology, Tohoku Medical and Pharmaceutical University, 4-4-1 Komatsushima, Aoba-ku, Sendai 981-8558, Miyagi, Japan
| | - Jun Uzawa
- Structural Glycobiology Team, RIKEN (The Institute of Physical and Chemical Research), 2-1 Hirosawa, Wako 351-0198, Saitama, Japan
| | - Yoshiki Yamaguchi
- Division of Structural Biology, Institute of Molecular Biomembrane and Glycobiology, Tohoku Medical and Pharmaceutical University, 4-4-1 Komatsushima, Aoba-ku, Sendai 981-8558, Miyagi, Japan
- Structural Glycobiology Team, RIKEN (The Institute of Physical and Chemical Research), 2-1 Hirosawa, Wako 351-0198, Saitama, Japan
| |
Collapse
|
18
|
Zhou X, Wei Z, Lu H, He J, Gao Y, Hu X, Wang C, Dong Y, Liu H. Large-Scale Molecular Dynamics Simulation Based on Heterogeneous Many-Core Architecture. J Chem Inf Model 2024; 64:851-861. [PMID: 38299978 DOI: 10.1021/acs.jcim.3c01254] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2024]
Abstract
As the application of molecular dynamics (MD) simulations continues to evolve, the demand for accelerating large-scale simulation systems and handling of enormous simulation tasks is steadily increasing. We propose a parallel acceleration method for large-scale MD simulations based on Sunway heterogeneous many-core processors. This method integrates task scheduling, simulation calculations, and data storage, effectively tackling issues related to large-scale simulations and numerous simulation tasks. The task scheduling strategy flexibly handles tasks on various scales and enables parallel execution of multiple tasks. During the simulation calculations, we ported GROMACS to the Sunway architecture and accelerated the calculation of short-range forces through a heterogeneous processor. Our method achieves approximately 10-fold acceleration and 90% scalability when executing a single simulation task. When handling numerous simulation tasks, our method achieves parallel execution of all of the tasks with 90% scalability. By employing our method, we carried out 50 ns simulations on over 3000 distinct conotoxin structures individually within just 5 h. Additionally, we evaluated more than 200 protein-ligand complexes, and the simulation efficiency significantly exceeded that of midsized to small GPU clusters.
Collapse
Affiliation(s)
- Xu Zhou
- College of Computer Science and Technology, Ocean University of China, Qingdao 266100, China
| | - Zhiqiang Wei
- College of Computer Science and Technology, Ocean University of China, Qingdao 266100, China
| | - Hao Lu
- College of Computer Science and Technology, Ocean University of China, Qingdao 266100, China
| | - Jiaqi He
- College of Computer Science and Technology, Ocean University of China, Qingdao 266100, China
| | - Yuan Gao
- College of Computer Science and Technology, Ocean University of China, Qingdao 266100, China
| | - Xiaotong Hu
- College of Computer Science and Technology, Ocean University of China, Qingdao 266100, China
| | - Cunji Wang
- College of Computer Science and Technology, Ocean University of China, Qingdao 266100, China
| | - Yujie Dong
- College of Computer Science and Technology, Ocean University of China, Qingdao 266100, China
| | - Hao Liu
- College of Computer Science and Technology, Ocean University of China, Qingdao 266100, China
| |
Collapse
|
19
|
Vuillemot R, Harastani M, Hamitouche I, Jonic S. MDSPACE and MDTOMO Software for Extracting Continuous Conformational Landscapes from Datasets of Single Particle Images and Subtomograms Based on Molecular Dynamics Simulations: Latest Developments in ContinuousFlex Software Package. Int J Mol Sci 2023; 25:20. [PMID: 38203192 PMCID: PMC10779004 DOI: 10.3390/ijms25010020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Revised: 12/16/2023] [Accepted: 12/17/2023] [Indexed: 01/12/2024] Open
Abstract
Cryo electron microscopy (cryo-EM) instrumentation allows obtaining 3D reconstruction of the structure of biomolecular complexes in vitro (purified complexes studied by single particle analysis) and in situ (complexes studied in cells by cryo electron tomography). Standard cryo-EM approaches allow high-resolution reconstruction of only a few conformational states of a molecular complex, as they rely on data classification into a given number of classes to increase the resolution of the reconstruction from the most populated classes while discarding all other classes. Such discrete classification approaches result in a partial picture of the full conformational variability of the complex, due to continuous conformational transitions with many, uncountable intermediate states. In this article, we present the software with a user-friendly graphical interface for running two recently introduced methods, namely, MDSPACE and MDTOMO, to obtain continuous conformational landscapes of biomolecules by analyzing in vitro and in situ cryo-EM data (single particle images and subtomograms) based on molecular dynamics simulations of an available atomic model of one of the conformations. The MDSPACE and MDTOMO software is part of the open-source ContinuousFlex software package (starting from version 3.4.2 of ContinuousFlex), which can be run as a plugin of the Scipion software package (version 3.1 and later), broadly used in the cryo-EM field.
Collapse
Affiliation(s)
| | | | | | - Slavica Jonic
- IMPMC-UMR 7590 CNRS, Sorbonne Université, MNHN, 75005 Paris, France
| |
Collapse
|
20
|
Ishida H, Matsumoto A, Tanaka H, Okuda A, Morishima K, Wade PA, Kurumizaka H, Sugiyama M, Kono H. Structural and Dynamic Changes of Nucleosome upon GATA3 Binding. J Mol Biol 2023; 435:168308. [PMID: 37805066 PMCID: PMC10843466 DOI: 10.1016/j.jmb.2023.168308] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Revised: 09/30/2023] [Accepted: 10/03/2023] [Indexed: 10/09/2023]
Abstract
Pioneer factors, which can directly bind to nucleosomes, have been considered to change chromatin conformations. However, the binding impact on the nucleosome is little known. Here, we show how the pioneer factor GATA3 binds to nucleosomal DNA and affects the conformation and dynamics of nucleosomes by using a combination of SAXS, molecular modeling, and molecular dynamics simulations. Our structural models, consistent with the SAXS data, indicate that only one of the two DNA binding domains, N- and C-fingers, of GATA3 binds to an end of the DNA in solution. Our MD simulations further showed that the other unbound end of the DNA increases the fluctuation and enhances the DNA dissociation from the histone core when the N-finger binds to a DNA end, a site near the entry or exit of the nucleosome. However, this was not true for the binding of the C-finger that binds to a location about 15 base pairs distant from the DNA end. In this case, DNA dissociation occurred on the bound end. Taken together, we suggest that the N-finger and C-finger bindings of GATA3 commonly enhance DNA dissociation at one of the two DNA ends (the bound end for the C-finger binding and the unbound end for the N-finger binding), leading to triggering a conformational change in the chromatin.
Collapse
Affiliation(s)
- Hisashi Ishida
- Institute for Quantum Life Science, Institutes for Quantum Science and Technology, 4-9-1 Anagawa, Inage-ku, Chiba city, Chiba 263-8555, Japan
| | - Atsushi Matsumoto
- Institute for Quantum Life Science, Institutes for Quantum Science and Technology, 4-9-1 Anagawa, Inage-ku, Chiba city, Chiba 263-8555, Japan
| | - Hiroki Tanaka
- Laboratory of Chromatin Structure and Function, Institute for Quantitative Biosciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo 113-0032, Japan; Present address: Department of Structural Virology, National Center for Global Health and Medicine, 1-21-1 Toyama, Shinjuku-ku, Tokyo 162-8655, Japan
| | - Aya Okuda
- Institute for Integrated Radiation and Nuclear Science, Kyoto University, 2-1010 Asashironishi, Kumatori, Sennan-gun, Osaka, 590-0494, Japan
| | - Ken Morishima
- Institute for Integrated Radiation and Nuclear Science, Kyoto University, 2-1010 Asashironishi, Kumatori, Sennan-gun, Osaka, 590-0494, Japan
| | - Paul A Wade
- Epigenetics and Stem Cell Biology Laboratory, National Institute of Environmental Health Sciences, Research Triangle Park, Durham, NC 27709, USA
| | - Hitoshi Kurumizaka
- Laboratory of Chromatin Structure and Function, Institute for Quantitative Biosciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo 113-0032, Japan
| | - Masaaki Sugiyama
- Institute for Integrated Radiation and Nuclear Science, Kyoto University, 2-1010 Asashironishi, Kumatori, Sennan-gun, Osaka, 590-0494, Japan
| | - Hidetoshi Kono
- Institute for Quantum Life Science, Institutes for Quantum Science and Technology, 4-9-1 Anagawa, Inage-ku, Chiba city, Chiba 263-8555, Japan; Graduate School of Science, Chiba University, 1-33 Yayoi-cho, Inage-ku, Chiba 263-8522, Japan.
| |
Collapse
|
21
|
Borges-Araújo L, Patmanidis I, Singh AP, Santos LHS, Sieradzan AK, Vanni S, Czaplewski C, Pantano S, Shinoda W, Monticelli L, Liwo A, Marrink SJ, Souza PCT. Pragmatic Coarse-Graining of Proteins: Models and Applications. J Chem Theory Comput 2023; 19:7112-7135. [PMID: 37788237 DOI: 10.1021/acs.jctc.3c00733] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.
Collapse
Affiliation(s)
- Luís Borges-Araújo
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Ilias Patmanidis
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Akhil P Singh
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
| | - Lucianna H S Santos
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
- Institut de Pharmacologie Moléculaire et Cellulaire, Université Côte d'Azur, Inserm, CNRS, 06560 Valbonne, France
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Sergio Pantano
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita, Okayama 700-8530, Japan
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Adam Liwo
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| |
Collapse
|
22
|
Kasahara K, Masayama R, Okita K, Matubayasi N. Elucidating protein-ligand binding kinetics based on returning probability theory. J Chem Phys 2023; 159:134103. [PMID: 37787130 DOI: 10.1063/5.0165692] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2023] [Accepted: 09/14/2023] [Indexed: 10/04/2023] Open
Abstract
The returning probability (RP) theory, a rigorous diffusion-influenced reaction theory, enables us to analyze the binding process systematically in terms of thermodynamics and kinetics using molecular dynamics (MD) simulations. Recently, the theory was extended to atomistically describe binding processes by adopting the host-guest interaction energy as the reaction coordinate. The binding rate constants can be estimated by computing the thermodynamic and kinetic properties of the reactive state existing in the binding processes. Here, we propose a methodology based on the RP theory in conjunction with the energy representation theory of solution, applicable to complex binding phenomena, such as protein-ligand binding. The derived scheme of calculating the equilibrium constant between the reactive and dissociate states, required in the RP theory, can be used for arbitrary types of reactive states. We apply the present method to the bindings of small fragment molecules [4-hydroxy-2-butanone (BUT) and methyl methylthiomethyl sulphoxide (DSS)] to FK506 binding protein (FKBP) in an aqueous solution. Estimated binding rate constants are consistent with those obtained from long-timescale MD simulations. Furthermore, by decomposing the rate constants to the thermodynamic and kinetic contributions, we clarify that the higher thermodynamic stability of the reactive state for DSS causes the faster binding kinetics compared with BUT.
Collapse
Affiliation(s)
- Kento Kasahara
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Ren Masayama
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Kazuya Okita
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
23
|
Miyashita O, Tama F. Advancing cryo-electron microscopy data analysis through accelerated simulation-based flexible fitting approaches. Curr Opin Struct Biol 2023; 82:102653. [PMID: 37451233 DOI: 10.1016/j.sbi.2023.102653] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2023] [Revised: 05/30/2023] [Accepted: 06/19/2023] [Indexed: 07/18/2023]
Abstract
Flexible fitting based on molecular dynamics simulation is a technique for structure modeling from cryo-EM data. It has been utilized for nearly two decades, and while cryo-EM resolution has improved significantly, it remains a powerful approach that can provide structural and dynamical insights that are not directly accessible from experimental data alone. Molecular dynamics simulations provide a means to extract atomistic details of conformational changes that are encoded in cryo-EM data and can also assist in improving the quality of structural models. Additionally, molecular dynamics simulations enable the characterization of conformational heterogeneity in cryo-EM data. We will summarize the advancements made in these techniques and highlight recent developments in this field.
Collapse
Affiliation(s)
- Osamu Miyashita
- RIKEN Center for Computational Science, 6-7-1, Minatojima-minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.
| | - Florence Tama
- RIKEN Center for Computational Science, 6-7-1, Minatojima-minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan; Department of Physics, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan; Institute of Transformative Bio-Molecules, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan.
| |
Collapse
|
24
|
Ugarte La Torre D, Takada S, Sugita Y. Extension of the iSoLF implicit-solvent coarse-grained model for multicomponent lipid bilayers. J Chem Phys 2023; 159:075101. [PMID: 37581417 DOI: 10.1063/5.0160417] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2023] [Accepted: 07/26/2023] [Indexed: 08/16/2023] Open
Abstract
iSoLF is a coarse-grained (CG) model for lipid molecules with the implicit-solvent approximation used in molecular dynamics (MD) simulations of biological membranes. Using the original iSoLF (iSoLFv1), MD simulations of lipid bilayers consisting of either POPC or DPPC and these bilayers, including membrane proteins, can be performed. Here, we improve the original model, explicitly treating the electrostatic interactions between different lipid molecules and adding CG particle types. As a result, the available lipid types increase to 30. To parameterize the potential functions of the new model, we performed all-atom MD simulations of each lipid at three different temperatures using the CHARMM36 force field and the modified TIP3P model. Then, we parameterized both the bonded and non-bonded interactions to fit the area per lipid and the membrane thickness of each lipid bilayer by using the multistate Boltzmann Inversion method. The final model reproduces the area per lipid and the membrane thickness of each lipid bilayer at the three temperatures. We also examined the applicability of the new model, iSoLFv2, to simulate the phase behaviors of mixtures of DOPC and DPPC at different concentrations. The simulation results with iSoLFv2 are consistent with those using Dry Martini and Martini 3, although iSoLFv2 requires much fewer computations. iSoLFv2 has been implemented in the GENESIS MD software and is publicly available.
Collapse
Affiliation(s)
- Diego Ugarte La Torre
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Yuji Sugita
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo, Japan
| |
Collapse
|
25
|
Jiang W. Enhanced Configurational Sampling Approaches to Alchemical Ligand Binding Free Energy Simulations: Current Status and Challenges. J Phys Chem B 2023; 127:6835-6841. [PMID: 37499215 DOI: 10.1021/acs.jpcb.3c02020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
Ligand binding free energy simulations (LB-FES) have been routine tasks in modern drug discovery campaign. A long-standing challenge for LB-FES is the difficulty in adequately sampling nontrivial environmental reorganizations in response to ligand binding. Therefore, various enhanced configurational sampling (ECS) approaches were devised to speed up fluctuations of relevant slow degrees of freedom (SDOF) and ensure simulation convergence. However, in contrast to the achievements in parametrization, software performance, and workflow automation, efficient ECS methodology suitable for high throughput screening remains in an early stage of development. Here, a review of ECS developments with LB-FES is presented, revisiting current approaches and underlining the major technical pitfalls and challenges. This Perspective focuses on alchemical LB-FES on account of their predominant role in high throughput drug screening as well as the established partnership with ECS. The critical aspects of designing ECS approaches, from both theoretical and applied perspectives, are described. This work is intended to provide a contemporary review of the scientific, technical, and practical issues associated with the accelerating convergence of alchemical LB-FES.
Collapse
Affiliation(s)
- Wei Jiang
- Computational Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Building 240, Argonne, Illinois 60439, United States
| |
Collapse
|
26
|
MacCallum JL, Hu S, Lenz S, Souza PCT, Corradi V, Tieleman DP. An implementation of the Martini coarse-grained force field in OpenMM. Biophys J 2023; 122:2864-2870. [PMID: 37050876 PMCID: PMC10398343 DOI: 10.1016/j.bpj.2023.04.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2022] [Revised: 03/05/2023] [Accepted: 04/06/2023] [Indexed: 04/14/2023] Open
Abstract
We describe a complete implementation of Martini 2 and Martini 3 in the OpenMM molecular dynamics software package. Martini is a widely used coarse-grained force field with applications in biomolecular simulation, materials, and broader areas of chemistry. It is implemented as a force field but makes extensive use of facilities unique to the GROMACS software, including virtual sites and bonded terms that are not commonly used in standard atomistic force fields. OpenMM is a flexible molecular dynamics package widely used for methods development and is competitive in speed on GPUs with other commonly used packages. OpenMM has facilities to easily implement new force field terms, external forces and fields, and other nonstandard features, which we use to implement all force field terms used in Martini 2 and Martini 3. This allows Martini simulations, starting with GROMACS topology files that are processed by custom scripts, with all the added flexibility of OpenMM. We provide a GitHub repository with test cases, compare accuracy and performance between GROMACS and OpenMM, and discuss the limitations of our implementation in terms of direct comparison with GROMACS. We describe a use case that implements the Modeling Employing Limited Data method to apply experimental constraints in a Martini simulation to efficiently determine the structure of a protein complex. We also discuss issues and a potential solution with the Martini 2 topology for cholesterol.
Collapse
Affiliation(s)
- Justin L MacCallum
- Department of Chemistry, University of Calgary, Calgary, Alberta, Canada; Centre for Molecular Simulation, University of Calgary, Calgary, Alberta, Canada.
| | - Shangnong Hu
- Department of Biological Sciences, University of Calgary, Calgary, Alberta, Canada; Centre for Molecular Simulation, University of Calgary, Calgary, Alberta, Canada
| | - Stefan Lenz
- Department of Chemistry, University of Calgary, Calgary, Alberta, Canada; Centre for Molecular Simulation, University of Calgary, Calgary, Alberta, Canada
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB - UMR 5086), CNRS & University of Lyon, Lyon, France
| | - Valentina Corradi
- Department of Biological Sciences, University of Calgary, Calgary, Alberta, Canada; Centre for Molecular Simulation, University of Calgary, Calgary, Alberta, Canada
| | - D Peter Tieleman
- Department of Biological Sciences, University of Calgary, Calgary, Alberta, Canada; Centre for Molecular Simulation, University of Calgary, Calgary, Alberta, Canada.
| |
Collapse
|
27
|
Dokainish HM, Sugita Y. Structural effects of spike protein D614G mutation in SARS-CoV-2. Biophys J 2023; 122:2910-2920. [PMID: 36397671 PMCID: PMC9671695 DOI: 10.1016/j.bpj.2022.11.025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2022] [Revised: 09/02/2022] [Accepted: 11/15/2022] [Indexed: 11/18/2022] Open
Abstract
A single mutation from aspartate to glycine at position 614 has dominated all circulating variants of the severe acute respiratory syndrome coronavirus 2. D614G mutation induces structural changes in the spike (S) protein that strengthen the virus infectivity. Here, we use molecular dynamics simulations to dissect the effects of mutation and 630-loop rigidification on S-protein structure. The introduction of the mutation orders the 630-loop structure and thereby induces global structural changes toward the cryoelectron microscopy structure of the D614G S-protein. The ordered 630-loop weakens local interactions between the 614th residue and others in contrast to disordered structures in the wild-type protein. The mutation allosterically alters global interactions between receptor-binding domains, forming an asymmetric and mobile down conformation and facilitating transitions toward up conformation. The loss of salt bridge between D614 and K854 upon the mutation generally stabilizes S-protein protomer, including the fusion peptide proximal region that mediates membrane fusion. Understanding the molecular basis of D614G mutation is crucial as it dominates in all variants of concern, including Delta and Omicron.
Collapse
Affiliation(s)
- Hisham M Dokainish
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
| | - Yuji Sugita
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan; Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan; Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo, Japan.
| |
Collapse
|
28
|
Mizutani A, Tan C, Sugita Y, Takada S. Micelle-like clusters in phase-separated Nanog condensates: A molecular simulation study. PLoS Comput Biol 2023; 19:e1011321. [PMID: 37486948 PMCID: PMC10399900 DOI: 10.1371/journal.pcbi.1011321] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Revised: 08/03/2023] [Accepted: 06/30/2023] [Indexed: 07/26/2023] Open
Abstract
The phase separation model for transcription suggests that transcription factors (TFs), coactivators, and RNA polymerases form biomolecular condensates around active gene loci and regulate transcription. However, the structural details of condensates remain elusive. In this study, for Nanog, a master TF in mammalian embryonic stem cells known to form protein condensates in vitro, we examined protein structures in the condensates using residue-level coarse-grained molecular simulations. Human Nanog formed micelle-like clusters in the condensate. In the micelle-like cluster, the C-terminal disordered domains, including the tryptophan repeat (WR) regions, interacted with each other near the cluster center primarily via hydrophobic interaction. In contrast, hydrophilic disordered N-terminal and DNA-binding domains were exposed on the surface of the clusters. Electrostatic attractions of these surface residues were responsible for bridging multiple micelle-like structures in the condensate. The micelle-like structure and condensate were dynamic and liquid-like. Mutation of tryptophan residues in the WR region which was implicated to be important for a Nanog function resulted in dissolution of the Nanog condensate. Finally, to examine the impact of Nanog cluster to DNA, we added DNA fragments to the Nanog condensate. Nanog DNA-binding domains exposed to the surface of the micelle-like cluster could recruit more than one DNA fragments, making DNA-DNA distance shorter.
Collapse
Affiliation(s)
- Azuki Mizutani
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Cheng Tan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo, Japan
| | - Yuji Sugita
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Hirosawa, Wako, Saitama, Japan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo, Japan
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| |
Collapse
|
29
|
Vuillemot R, Rouiller I, Jonić S. MDTOMO method for continuous conformational variability analysis in cryo electron subtomograms based on molecular dynamics simulations. Sci Rep 2023; 13:10596. [PMID: 37391578 PMCID: PMC10313669 DOI: 10.1038/s41598-023-37037-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2023] [Accepted: 06/14/2023] [Indexed: 07/02/2023] Open
Abstract
Cryo electron tomography (cryo-ET) allows observing macromolecular complexes in their native environment. The common routine of subtomogram averaging (STA) allows obtaining the three-dimensional (3D) structure of abundant macromolecular complexes, and can be coupled with discrete classification to reveal conformational heterogeneity of the sample. However, the number of complexes extracted from cryo-ET data is usually small, which restricts the discrete-classification results to a small number of enough populated states and, thus, results in a largely incomplete conformational landscape. Alternative approaches are currently being investigated to explore the continuity of the conformational landscapes that in situ cryo-ET studies could provide. In this article, we present MDTOMO, a method for analyzing continuous conformational variability in cryo-ET subtomograms based on Molecular Dynamics (MD) simulations. MDTOMO allows obtaining an atomic-scale model of conformational variability and the corresponding free-energy landscape, from a given set of cryo-ET subtomograms. The article presents the performance of MDTOMO on a synthetic ABC exporter dataset and an in situ SARS-CoV-2 spike dataset. MDTOMO allows analyzing dynamic properties of molecular complexes to understand their biological functions, which could also be useful for structure-based drug discovery.
Collapse
Affiliation(s)
- Rémi Vuillemot
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, CC 115, 4 Place Jussieu, 75005, Paris, France
- Department of Biochemistry and Pharmacology and Bio21 Molecular Science and Biotechnology Institute, The University of Melbourne, Melbourne, VIC, 3010, Australia
| | - Isabelle Rouiller
- Department of Biochemistry and Pharmacology and Bio21 Molecular Science and Biotechnology Institute, The University of Melbourne, Melbourne, VIC, 3010, Australia
- Australian Research Council Centre for Cryo-Electron Microscopy of Membrane Proteins, Parkville, VIC, 3052, Australia
| | - Slavica Jonić
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, CC 115, 4 Place Jussieu, 75005, Paris, France.
| |
Collapse
|
30
|
Jaffrelot Inizan T, Plé T, Adjoua O, Ren P, Gökcan H, Isayev O, Lagardère L, Piquemal JP. Scalable hybrid deep neural networks/polarizable potentials biomolecular simulations including long-range effects. Chem Sci 2023; 14:5438-5452. [PMID: 37234902 PMCID: PMC10208042 DOI: 10.1039/d2sc04815a] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2022] [Accepted: 04/03/2023] [Indexed: 07/28/2023] Open
Abstract
Deep-HP is a scalable extension of the Tinker-HP multi-GPU molecular dynamics (MD) package enabling the use of Pytorch/TensorFlow Deep Neural Network (DNN) models. Deep-HP increases DNNs' MD capabilities by orders of magnitude offering access to ns simulations for 100k-atom biosystems while offering the possibility of coupling DNNs to any classical (FFs) and many-body polarizable (PFFs) force fields. It allows therefore the introduction of the ANI-2X/AMOEBA hybrid polarizable potential designed for ligand binding studies where solvent-solvent and solvent-solute interactions are computed with the AMOEBA PFF while solute-solute ones are computed by the ANI-2X DNN. ANI-2X/AMOEBA explicitly includes AMOEBA's physical long-range interactions via an efficient Particle Mesh Ewald implementation while preserving ANI-2X's solute short-range quantum mechanical accuracy. The DNN/PFF partition can be user-defined allowing for hybrid simulations to include key ingredients of biosimulation such as polarizable solvents, polarizable counter ions, etc.… ANI-2X/AMOEBA is accelerated using a multiple-timestep strategy focusing on the model's contributions to low-frequency modes of nuclear forces. It primarily evaluates AMOEBA forces while including ANI-2X ones only via correction-steps resulting in an order of magnitude acceleration over standard Velocity Verlet integration. Simulating more than 10 μs, we compute charged/uncharged ligand solvation free energies in 4 solvents, and absolute binding free energies of host-guest complexes from SAMPL challenges. ANI-2X/AMOEBA average errors are discussed in terms of statistical uncertainty and appear in the range of chemical accuracy compared to experiment. The availability of the Deep-HP computational platform opens the path towards large-scale hybrid DNN simulations, at force-field cost, in biophysics and drug discovery.
Collapse
Affiliation(s)
- Théo Jaffrelot Inizan
- Sorbonne Université, Laboratoire de Chimie Théorique UMR 7616 CNRS Paris 75005 France
| | - Thomas Plé
- Sorbonne Université, Laboratoire de Chimie Théorique UMR 7616 CNRS Paris 75005 France
| | - Olivier Adjoua
- Sorbonne Université, Laboratoire de Chimie Théorique UMR 7616 CNRS Paris 75005 France
| | - Pengyu Ren
- Department of Biomedical Engineering, University of Texas at Austin Austin Texas USA
| | - Hatice Gökcan
- Department of Chemistry, Carnegie Mellon University Pittsburgh Pennsylvania USA
| | - Olexandr Isayev
- Department of Chemistry, Carnegie Mellon University Pittsburgh Pennsylvania USA
| | - Louis Lagardère
- Sorbonne Université, Laboratoire de Chimie Théorique UMR 7616 CNRS Paris 75005 France
- Sorbonne Université, Institut Parisien de Chimie Physique et Théorique FR 2622 CNRS Paris France
| | - Jean-Philip Piquemal
- Sorbonne Université, Laboratoire de Chimie Théorique UMR 7616 CNRS Paris 75005 France
- Department of Biomedical Engineering, University of Texas at Austin Austin Texas USA
| |
Collapse
|
31
|
Jung J, Kobayashi C, Sugita Y. Acceleration of generalized replica exchange with solute tempering simulations of large biological systems on massively parallel supercomputer. J Comput Chem 2023. [PMID: 37141320 DOI: 10.1002/jcc.27124] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2022] [Revised: 04/10/2023] [Accepted: 04/14/2023] [Indexed: 05/06/2023]
Abstract
Generalized replica exchange with solute tempering (gREST) is one of the enhanced sampling algorithms for proteins or other systems with rugged energy landscapes. Unlike the replica-exchange molecular dynamics (REMD) method, solvent temperatures are the same in all replicas, while solute temperatures are different and are exchanged frequently between replicas for exploring various solute structures. Here, we apply the gREST scheme to large biological systems containing over one million atoms using a large number of processors in a supercomputer. First, communication time on a multi-dimensional torus network is reduced by matching each replica to MPI processors optimally. This is applicable not only to gREST but also to other multi-copy algorithms. Second, energy evaluations, which are necessary for the multistate bennet acceptance ratio (MBAR) method for free energy estimations, are performed on-the-fly during the gREST simulations. Using these two advanced schemes, we observed 57.72 ns/day performance in 128-replica gREST calculations with 1.5 million atoms system using 16,384 nodes in Fugaku. These schemes implemented in the latest version of GENESIS software could open new possibilities to answer unresolved questions on large biomolecular complex systems with slow conformational dynamics.
Collapse
Affiliation(s)
- Jaewoon Jung
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Saitama, Japan
| | - Chigusa Kobayashi
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Japan
| | - Yuji Sugita
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Saitama, Japan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Japan
| |
Collapse
|
32
|
Tan C, Niitsu A, Sugita Y. Highly Charged Proteins and Their Repulsive Interactions Antagonize Biomolecular Condensation. JACS AU 2023; 3:834-848. [PMID: 37006777 PMCID: PMC10052238 DOI: 10.1021/jacsau.2c00646] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Revised: 02/09/2023] [Accepted: 02/09/2023] [Indexed: 06/19/2023]
Abstract
Biomolecular condensation is involved in various cellular processes; therefore, regulation of condensation is crucial to prevent deleterious protein aggregation and maintain a stable cellular environment. Recently, a class of highly charged proteins, known as heat-resistant obscure (Hero) proteins, was shown to protect other client proteins from pathological aggregation. However, the molecular mechanisms by which Hero proteins protect other proteins from aggregation remain unknown. In this study, we performed multiscale molecular dynamics (MD) simulations of Hero11, a Hero protein, and the C-terminal low-complexity domain (LCD) of the transactive response DNA-binding protein 43 (TDP-43), a client protein of Hero11, under various conditions to examine their interactions with each other. We found that Hero11 permeates into the condensate formed by the LCD of TDP-43 (TDP-43-LCD) and induces changes in conformation, intermolecular interactions, and dynamics of TDP-43-LCD. We also examined possible Hero11 structures in atomistic and coarse-grained MD simulations and found that Hero11 with a higher fraction of disordered region tends to assemble on the surface of the condensates. Based on the simulation results, we have proposed three possible mechanisms for Hero11's regulatory function: (i) In the dense phase, TDP-43-LCD reduces contact with each other and shows faster diffusion and decondensation due to the repulsive Hero11-Hero11 interactions. (ii) In the dilute phase, the saturation concentration of TDP-43-LCD is increased, and its conformation is relatively more extended and variant, induced by the attractive Hero11-TDP-43-LCD interactions. (iii) Hero11 on the surface of small TDP-43-LCD condensates can contribute to avoiding their fusion due to repulsive interactions. The proposed mechanisms provide new insights into the regulation of biomolecular condensation in cells under various conditions.
Collapse
Affiliation(s)
- Cheng Tan
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Ai Niitsu
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
| | - Yuji Sugita
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
33
|
Ito S, Yagi K, Sugita Y. Allosteric regulation of β-reaction stage I in tryptophan synthase upon the α-ligand binding. J Chem Phys 2023; 158:115101. [PMID: 36948822 DOI: 10.1063/5.0134117] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2023] Open
Abstract
Tryptophan synthase (TRPS) is a bifunctional enzyme consisting of α- and β-subunits that catalyzes the last two steps of L-tryptophan (L-Trp) biosynthesis. The first stage of the reaction at the β-subunit is called β-reaction stage I, which converts the β-ligand from an internal aldimine [E(Ain)] to an α-aminoacrylate [E(A-A)] intermediate. The activity is known to increase 3-10-fold upon the binding of 3-indole-D-glycerol-3'-phosphate (IGP) at the α-subunit. The effect of α-ligand binding on β-reaction stage I at the distal β-active site is not well understood despite the abundant structural information available for TRPS. Here, we investigate the β-reaction stage I by carrying out minimum-energy pathway searches based on a hybrid quantum mechanics/molecular mechanics (QM/MM) model. The free-energy differences along the pathway are also examined using QM/MM umbrella sampling simulations with QM calculations at the B3LYP-D3/aug-cc-pVDZ level of theory. Our simulations suggest that the sidechain orientation of βD305 near the β-ligand likely plays an essential role in the allosteric regulation: a hydrogen bond is formed between βD305 and the β-ligand in the absence of the α-ligand, prohibiting a smooth rotation of the hydroxyl group in the quinonoid intermediate, whereas the dihedral angle rotates smoothly after the hydrogen bond is switched from βD305-β-ligand to βD305-βR141. This switch could occur upon the IGP-binding at the α-subunit, as evidenced by the existing TRPS crystal structures.
Collapse
Affiliation(s)
- Shingo Ito
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Kiyoshi Yagi
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| |
Collapse
|
34
|
Development of hidden Markov modeling method for molecular orientations and structure estimation from high-speed atomic force microscopy time-series images. PLoS Comput Biol 2022; 18:e1010384. [PMID: 36580448 PMCID: PMC9833559 DOI: 10.1371/journal.pcbi.1010384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2022] [Revised: 01/11/2023] [Accepted: 12/20/2022] [Indexed: 12/30/2022] Open
Abstract
High-speed atomic force microscopy (HS-AFM) is a powerful technique for capturing the time-resolved behavior of biomolecules. However, structural information in HS-AFM images is limited to the surface geometry of a sample molecule. Inferring latent three-dimensional structures from the surface geometry is thus important for getting more insights into conformational dynamics of a target biomolecule. Existing methods for estimating the structures are based on the rigid-body fitting of candidate structures to each frame of HS-AFM images. Here, we extend the existing frame-by-frame rigid-body fitting analysis to multiple frames to exploit orientational correlations of a sample molecule between adjacent frames in HS-AFM data due to the interaction with the stage. In the method, we treat HS-AFM data as time-series data, and they are analyzed with the hidden Markov modeling. Using simulated HS-AFM images of the taste receptor type 1 as a test case, the proposed method shows a more robust estimation of molecular orientations than the frame-by-frame analysis. The method is applicable in integrative modeling of conformational dynamics using HS-AFM data.
Collapse
|
35
|
Okita K, Kasahara K, Matubayasi N. Diffusion theory of molecular liquids in the energy representation and application to solvation dynamics. J Chem Phys 2022; 157:244505. [PMID: 36586971 DOI: 10.1063/5.0125432] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022] Open
Abstract
The generalized Langevin equation (GLE) formalism is a useful theoretical fundament for analyzing dynamical phenomena rigorously. Despite the systematic formulation of dynamics theories with practical approximations, however, the applicability of GLE-based methods is still limited to simple polyatomic liquids due to the approximate treatment of molecular orientations involved in the static molecular liquid theory. Here, we propose an exact framework of dynamics based on the GLE formalism incorporating the energy representation theory of solution, an alternative static molecular liquid theory. A fundamental idea is the projection of the relative positions and orientations of solvents around a solute onto the solute-solvent interaction, namely the energy coordinate, enabling us to describe the dynamics on a one-dimensional coordinate. Introducing systematic approximations, such as the overdamped limit, leads to the molecular diffusion equation in the energy representation that is described in terms of the distribution function of solvents on the energy coordinate and the diffusion coefficients. The present theory is applied to the solvation dynamics triggered by the photoexcitation of benzonitrile. The long-time behavior of the solvation time correlation function is in good agreement with that obtained by the molecular dynamics simulation.
Collapse
Affiliation(s)
- Kazuya Okita
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Kento Kasahara
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
36
|
Matsunaga Y, Kamiya M, Oshima H, Jung J, Ito S, Sugita Y. Use of multistate Bennett acceptance ratio method for free-energy calculations from enhanced sampling and free-energy perturbation. Biophys Rev 2022; 14:1503-1512. [PMID: 36659993 PMCID: PMC9842838 DOI: 10.1007/s12551-022-01030-9] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2022] [Accepted: 12/01/2022] [Indexed: 12/15/2022] Open
Abstract
Multistate Bennett acceptance ratio (MBAR) works as a method to analyze molecular dynamics (MD) simulation data after the simulations have been finished. It is widely used to estimate free-energy changes between different states and averaged properties at the states of interest. MBAR allows us to treat a wide range of states from those at different temperature/pressure to those with different model parameters. Due to the broad applicability, the MBAR equations are rather difficult to apply for free-energy calculations using different types of MD simulations including enhanced conformational sampling methods and free-energy perturbation. In this review, we first summarize the basic theory of the MBAR equations and categorize the representative usages into the following four: (i) perturbation, (ii) scaling, (iii) accumulation, and (iv) full potential energy. For each, we explain how to prepare input data using MD simulation trajectories for solving the MBAR equations. MBAR is also useful to estimate reliable free-energy differences using MD trajectories based on a semi-empirical quantum mechanics/molecular mechanics (QM/MM) model and ab initio QM/MM energy calculations on the MD snapshots. We also explain how to use the MBAR software in the GENESIS package, which we call mbar_analysis, for the four representative cases. The proposed estimations of free-energy changes and thermodynamic averages are effective and useful for various biomolecular systems.
Collapse
Affiliation(s)
- Yasuhiro Matsunaga
- grid.263023.60000 0001 0703 3735Graduate School of Science and Engineering, Saitama University, Saitama, Saitama 338-8570 Japan
| | - Motoshi Kamiya
- grid.467196.b0000 0001 2285 6123Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585 Japan
| | - Hiraku Oshima
- grid.508743.dLaboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047 Japan
| | - Jaewoon Jung
- grid.474693.bComputational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047 Japan
| | - Shingo Ito
- grid.7597.c0000000094465255Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama 351-0198 Japan
| | - Yuji Sugita
- grid.508743.dLaboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047 Japan ,grid.474693.bComputational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047 Japan ,grid.7597.c0000000094465255Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama 351-0198 Japan
| |
Collapse
|
37
|
Multiple sub-state structures of SERCA2b reveal conformational overlap at transition steps during the catalytic cycle. Cell Rep 2022; 41:111760. [PMID: 36476867 DOI: 10.1016/j.celrep.2022.111760] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Revised: 09/06/2022] [Accepted: 11/09/2022] [Indexed: 12/12/2022] Open
Abstract
Sarco/endoplasmic reticulum Ca2+ ATPase (SERCA) pumps Ca2+ into the endoplasmic reticulum (ER). Herein, we present cryo-electron microscopy (EM) structures of three intermediates of SERCA2b: Ca2+-bound phosphorylated (E1P·2Ca2+) and Ca2+-unbound dephosphorylated (E2·Pi) intermediates and another between the E2P and E2·Pi states. Our cryo-EM analysis demonstrates that the E1P·2Ca2+ state exists in low abundance and preferentially transitions to an E2P-like structure by releasing Ca2+ and that the Ca2+ release gate subsequently undergoes stepwise closure during the dephosphorylation processes. Importantly, each intermediate adopts multiple sub-state structures including those like the next one in the catalytic series, indicating conformational overlap at transition steps, as further substantiated by atomistic molecular dynamic simulations of SERCA2b in a lipid bilayer. The present findings provide insight into how enzymes accelerate catalytic cycles.
Collapse
|
38
|
Harastani M, Vuillemot R, Hamitouche I, Moghadam NB, Jonic S. ContinuousFlex: Software package for analyzing continuous conformational variability of macromolecules in cryo electron microscopy and tomography data. J Struct Biol 2022; 214:107906. [PMID: 36244611 DOI: 10.1016/j.jsb.2022.107906] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2022] [Revised: 09/02/2022] [Accepted: 10/07/2022] [Indexed: 11/06/2022]
Abstract
ContinuousFlex is a user-friendly open-source software package for analyzing continuous conformational variability of macromolecules in cryo electron microscopy (cryo-EM) and cryo electron tomography (cryo-ET) data. In 2019, ContinuousFlex became available as a plugin for Scipion, an image processing software package extensively used in the cryo-EM field. Currently, ContinuousFlex contains software for running (1) recently published methods HEMNMA-3D, TomoFlow, and NMMD; (2) earlier published methods HEMNMA and StructMap; and (3) methods for simulating cryo-EM and cryo-ET data with conformational variability and methods for data preprocessing. It also includes external software for molecular dynamics simulation (GENESIS) and normal mode analysis (ElNemo), used in some of the mentioned methods. The HEMNMA software has been presented in the past, but not the software of other methods. Besides, ContinuousFlex currently also offers a deep learning extension of HEMNMA, named DeepHEMNMA. In this article, we review these methods in the context of the ContinuousFlex package, developed to facilitate their use by the community.
Collapse
Affiliation(s)
- Mohamad Harastani
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, Paris, France
| | - Rémi Vuillemot
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, Paris, France
| | - Ilyes Hamitouche
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, Paris, France
| | - Nima Barati Moghadam
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, Paris, France
| | - Slavica Jonic
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, Paris, France.
| |
Collapse
|
39
|
Yamane T, Ekimoto T, Ikeguchi M. Development of the force field for cyclosporine A. Biophys Physicobiol 2022; 19:e190045. [PMID: 36567735 PMCID: PMC9751258 DOI: 10.2142/biophysico.bppb-v19.0045] [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: 06/21/2022] [Accepted: 11/17/2022] [Indexed: 11/19/2022] Open
Abstract
Membrane permeability of cyclic peptides is an important factor in drug design. To investigate the membrane permeability of cyclic peptides using molecular dynamics (MD) simulations, the accurate force fields for unnatural amino acids present in the cyclic peptides are required. Therefore, we developed the CHARMM force fields of the unnatural amino acids present in cyclosporin A (CsA), a cyclic peptide used as an immune suppressor. Especially for N-methyl amino acids, which contribute to the membrane permeability of cyclic peptides, we developed a grid correction map (CMAP) of the energy surface using the φ and ψ dihedral angles in the main chain of CsA. To validate the developed force field, we performed MD simulations, including the generalized replica exchange with solute tempering method, of CsA in water and chloroform solvents. The conformations of CsA in water and chloroform sampled using the developed force field were consistent with those of the experimental results of the solution nuclear magnetic resonance spectroscopy.
Collapse
Affiliation(s)
- Tsutomu Yamane
- Graduate School of Medical Life Science, Yokohama City University, Yokohama, Kanagawa 230-0045, Japan,HPC- and AI-driven Drug Development Platform Division, Center for Computational Science, RIKEN, Yokohama, Kanagawa 230-0045, Japan
| | - Toru Ekimoto
- Graduate School of Medical Life Science, Yokohama City University, Yokohama, Kanagawa 230-0045, Japan
| | - Mitsunori Ikeguchi
- Graduate School of Medical Life Science, Yokohama City University, Yokohama, Kanagawa 230-0045, Japan,HPC- and AI-driven Drug Development Platform Division, Center for Computational Science, RIKEN, Yokohama, Kanagawa 230-0045, Japan
| |
Collapse
|
40
|
Hosoe Y, Miyanoiri Y, Re S, Ochi S, Asahina Y, Kawakami T, Kuroda M, Mizuguchi K, Oda M. Structural dynamics of the N‐terminal
SH2
domain of
PI3K
in its free and
CD28
‐bound states. FEBS J 2022; 290:2366-2378. [PMID: 36282120 DOI: 10.1111/febs.16666] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2022] [Revised: 08/29/2022] [Accepted: 10/24/2022] [Indexed: 11/06/2022]
Abstract
Protein conformational changes with fluctuations are fundamental aspects of protein-protein interactions (PPIs); understanding these motions is required for the rational design of PPI-regulating compounds. Src homology 2 (SH2) domains are commonly found in adapter proteins involved in signal transduction and specifically bind to consensus motifs of proteins containing phosphorylated tyrosine (pY). Here, we analysed the interaction between the N-terminal SH2 domain (nSH2) of the regulatory subunit in phosphoinositide 3-kinase (PI3K) and the cytoplasmic region of the T-cell co-receptor, CD28, using NMR and molecular dynamics (MD) simulations. First, we assigned the backbone signals of nSH2 on 1 H-15 N heteronuclear single quantum coherence spectra in the absence or presence of the CD28 phosphopeptide, SDpYMNMTPRRPG. Chemical shift perturbation experiments revealed allosteric changes at the BC loop and the C-terminal region of nSH2 upon CD28 binding. NMR relaxation experiments showed a conformational exchange associated with CD28 binding in these regions. The conformational stabilisation of the C-terminal region correlated with the regulation of PI3K catalytic function. Further, using 19 F- and 31 P-labelled CD28 phosphopeptide, we analysed the structural dynamics of CD28 and demonstrated that the aromatic ring of the pY residue fluctuated between multiple conformations upon nSH2 binding. Our MD simulations largely explained the NMR results and the structural dynamics of nSH2 and CD28 in both bound and unbound states. Notably, in addition to its major conformation, we detected a minor conformation of nSH2 in the CD28 bound state that may explain the allosteric conformational change in the BC loop.
Collapse
Affiliation(s)
- Yuhi Hosoe
- Graduate School of Life and Environmental Sciences Kyoto Prefectural University Japan
| | | | - Suyong Re
- Artificial Intelligence Center for Health and Biomedical Research National Institutes of Biomedical Innovation, Health, and Nutrition Osaka Japan
| | - Saki Ochi
- Graduate School of Life and Environmental Sciences Kyoto Prefectural University Japan
| | - Yuya Asahina
- Institute for Protein Research Osaka University Japan
| | - Toru Kawakami
- Institute for Protein Research Osaka University Japan
| | - Masataka Kuroda
- Artificial Intelligence Center for Health and Biomedical Research National Institutes of Biomedical Innovation, Health, and Nutrition Osaka Japan
- Discovery Technology Laboratories Mitsubishi Tanabe Pharma Corporation Yokohama Japan
| | - Kenji Mizuguchi
- Institute for Protein Research Osaka University Japan
- Artificial Intelligence Center for Health and Biomedical Research National Institutes of Biomedical Innovation, Health, and Nutrition Osaka Japan
| | - Masayuki Oda
- Graduate School of Life and Environmental Sciences Kyoto Prefectural University Japan
| |
Collapse
|
41
|
Hanaoka K, Iwaki S, Yagi K, Myochin T, Ikeno T, Ohno H, Sasaki E, Komatsu T, Ueno T, Uchigashima M, Mikuni T, Tainaka K, Tahara S, Takeuchi S, Tahara T, Uchiyama M, Nagano T, Urano Y. General Design Strategy to Precisely Control the Emission of Fluorophores via a Twisted Intramolecular Charge Transfer (TICT) Process. J Am Chem Soc 2022; 144:19778-19790. [PMID: 36191139 DOI: 10.1021/jacs.2c06397] [Citation(s) in RCA: 31] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Fluorogenic probes for bioimaging have become essential tools for life science and medicine, and the key to their development is a precise understanding of the mechanisms available for fluorescence off/on control, such as photoinduced electron transfer (PeT) and Förster resonance energy transfer (FRET). Here we establish a new molecular design strategy to rationally develop activatable fluorescent probes, which exhibit a fluorescence off/on change in response to target biomolecules, by controlling the twisted intramolecular charge transfer (TICT) process. This approach was developed on the basis of a thorough investigation of the fluorescence quenching mechanism of N-phenyl rhodamine dyes (commercially available as the QSY series) by means of time-dependent density functional theory (TD-DFT) calculations and photophysical evaluation of their derivatives. To illustrate and validate this TICT-based design strategy, we employed it to develop practical fluorogenic probes for HaloTag and SNAP-tag. We further show that the TICT-controlled fluorescence off/on mechanism is generalizable by synthesizing a Si-rhodamine-based fluorogenic probe for HaloTag, thus providing a palette of chemical dyes that spans the visible and near-infrared range.
Collapse
Affiliation(s)
- Kenjiro Hanaoka
- Graduate School of Pharmaceutical Sciences, Keio University, 1-5-30 Shibakoen, Minoto-ku, Tokyo105-8512, Japan
| | | | - Kiyoshi Yagi
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama351-0198, Japan
| | | | | | - Hisashi Ohno
- Graduate School of Pharmaceutical Sciences, Keio University, 1-5-30 Shibakoen, Minoto-ku, Tokyo105-8512, Japan
| | - Eita Sasaki
- Graduate School of Pharmaceutical Sciences, Keio University, 1-5-30 Shibakoen, Minoto-ku, Tokyo105-8512, Japan
| | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
42
|
Matsubara D, Kasahara K, Dokainish HM, Oshima H, Sugita Y. Modified Protein-Water Interactions in CHARMM36m for Thermodynamics and Kinetics of Proteins in Dilute and Crowded Solutions. Molecules 2022; 27:molecules27175726. [PMID: 36080494 PMCID: PMC9457699 DOI: 10.3390/molecules27175726] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Revised: 08/30/2022] [Accepted: 08/30/2022] [Indexed: 11/16/2022] Open
Abstract
Proper balance between protein-protein and protein-water interactions is vital for atomistic molecular dynamics (MD) simulations of globular proteins as well as intrinsically disordered proteins (IDPs). The overestimation of protein-protein interactions tends to make IDPs more compact than those in experiments. Likewise, multiple proteins in crowded solutions are aggregated with each other too strongly. To optimize the balance, Lennard-Jones (LJ) interactions between protein and water are often increased about 10% (with a scaling parameter, λ = 1.1) from the existing force fields. Here, we explore the optimal scaling parameter of protein-water LJ interactions for CHARMM36m in conjunction with the modified TIP3P water model, by performing enhanced sampling MD simulations of several peptides in dilute solutions and conventional MD simulations of globular proteins in dilute and crowded solutions. In our simulations, 10% increase of protein-water LJ interaction for the CHARMM36m cannot maintain stability of a small helical peptide, (AAQAA)3 in a dilute solution and only a small modification of protein-water LJ interaction up to the 3% increase (λ = 1.03) is allowed. The modified protein-water interactions are applicable to other peptides and globular proteins in dilute solutions without changing thermodynamic properties from the original CHARMM36m. However, it has a great impact on the diffusive properties of proteins in crowded solutions, avoiding the formation of too sticky protein-protein interactions.
Collapse
Affiliation(s)
- Daiki Matsubara
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe 650-0047, Hyogo, Japan
| | - Kento Kasahara
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe 650-0047, Hyogo, Japan
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka 560-8531, Osaka, Japan
| | - Hisham M. Dokainish
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako 351-0198, Saitama, Japan
| | - Hiraku Oshima
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe 650-0047, Hyogo, Japan
| | - Yuji Sugita
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe 650-0047, Hyogo, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako 351-0198, Saitama, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe 650-0047, Hyogo, Japan
- Correspondence: ; Tel.: +81-48-462-1407
| |
Collapse
|
43
|
Oide M, Sugita Y. Protein Folding Intermediates on the Dimensionality Reduced Landscape with UMAP and Native Contact Likelihood. J Chem Phys 2022; 157:075101. [DOI: 10.1063/5.0099094] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
To understand protein folding mechanisms from molecular dynamics (MD) simulations, it is important to explore not only folded/unfolded states but also representative intermediate structures on the conformational landscape. Here, we propose a novel approach to construct the landscape using the uniform manifold approximation and projection (UMAP) method, which reduces the dimensionality without losing data-point proximity. In the approach, native contact likelihood is used as feature variables rather than the conventional Cartesian coordinates or dihedral angles of protein structures. We tested the performance of UMAP for coarse-grained MD simulation trajectories of B1 domain in protein G and observed on-pathway transient structures and other metastable states on the UMAP conformational landscape. In contrast, these structures were not clearly distinguished on the dimensionality reduced landscape using principal component analysis (PCA) or time-lagged independent component analysis (tICA). This approach is also useful to obtain dynamical information through Markov State Modeling and would be applicable to large-scale conformational changes in many other biomacromolecules.
Collapse
Affiliation(s)
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN, Japan
| |
Collapse
|
44
|
Sato K, Sasaki R, Matsuda R, Nakagawa M, Ekimoto T, Yamane T, Ikeguchi M, Tabata KV, Noji H, Kinbara K. Correction to "Supramolecular Mechanosensitive Potassium Channel Formed by Fluorinated Amphiphilic Cyclophane". J Am Chem Soc 2022; 144:13983-13984. [PMID: 35861320 DOI: 10.1021/jacs.2c07339] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
45
|
Marrink SJ, Monticelli L, Melo MN, Alessandri R, Tieleman DP, Souza PCT. Two decades of Martini: Better beads, broader scope. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1620] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Siewert J. Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials University of Groningen Groningen The Netherlands
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry (MMSB ‐ UMR 5086) CNRS & University of Lyon Lyon France
| | - Manuel N. Melo
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa Oeiras Portugal
| | - Riccardo Alessandri
- Pritzker School of Molecular Engineering University of Chicago Chicago Illinois USA
| | - D. Peter Tieleman
- Centre for Molecular Simulation and Department of Biological Sciences University of Calgary Alberta Canada
| | - Paulo C. T. Souza
- Molecular Microbiology and Structural Biochemistry (MMSB ‐ UMR 5086) CNRS & University of Lyon Lyon France
| |
Collapse
|
46
|
Oshima H, Sugita Y. Modified Hamiltonian in FEP Calculations for Reducing the Computational Cost of Electrostatic Interactions. J Chem Inf Model 2022; 62:2846-2856. [PMID: 35639709 DOI: 10.1021/acs.jcim.1c01532] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The free-energy perturbation (FEP) method predicts relative and absolute free-energy changes of biomolecules in solvation and binding with other molecules. FEP is, therefore, one of the most essential tools in in silico drug design. In conventional FEP, to smoothly connect two thermodynamic states, the potential energy is modified as a linear combination of the end-state potential energies by introducing scaling factors. When the particle mesh Ewald is used for electrostatic calculations, conventional FEP requires two reciprocal-space calculations per time step, which largely decreases the computational performance. To overcome this problem, we propose a new FEP scheme by introducing a modified Hamiltonian instead of interpolation of the end-state potential energies. The scheme introduces nonuniform scaling into the electrostatic potential as used in Replica Exchange with Solute Tempering 2 (REST2) and does not require additional reciprocal-space calculations. We tested this modified Hamiltonian in FEP calculations in several biomolecular systems. In all cases, the calculated free-energy changes with the current scheme are in good agreement with those from conventional FEP. The modified Hamiltonian in FEP greatly improves the computational performance, which is particularly marked for large biomolecular systems whose reciprocal-space calculations are the major bottleneck of total computational time.
Collapse
Affiliation(s)
- Hiraku Oshima
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Integrated Innovation Building 7F, 6-7-1 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Integrated Innovation Building 7F, 6-7-1 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako-shi, Saitama 351-0198, Japan.,Computational Biophysics Research Team, RIKEN Center for Computational Science, Integrated Innovation Building 7F, 6-7-1 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
47
|
Kasahara K, Masayama R, Matsubara Y, Matubayasi N. Constructing a Memory Kernel of the Returning Probability to Efficiently Describe Molecular Binding Processes. CHEM LETT 2022. [DOI: 10.1246/cl.220236] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Kento Kasahara
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Ren Masayama
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Yuya Matsubara
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
48
|
Kohga H, Mori T, Tanaka Y, Yoshikaie K, Taniguchi K, Fujimoto K, Fritz L, Schneider T, Tsukazaki T. Crystal structure of the lipid flippase MurJ in a "squeezed" form distinct from its inward- and outward-facing forms. Structure 2022; 30:1088-1097.e3. [PMID: 35660157 DOI: 10.1016/j.str.2022.05.008] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2022] [Revised: 03/15/2022] [Accepted: 05/10/2022] [Indexed: 11/29/2022]
Abstract
The bacterial peptidoglycan enclosing the cytoplasmic membrane is a fundamental cellular architecture. The integral membrane protein MurJ plays an essential role in flipping the cell wall building block Lipid II across the cytoplasmic membrane for peptidoglycan biosynthesis. Previously reported crystal structures of MurJ have elucidated its V-shaped inward- or outward-facing forms with an internal cavity for substrate binding. MurJ transports Lipid II using its cavity through conformational transitions between these two forms. Here, we report two crystal structures of inward-facing forms from Arsenophonus endosymbiont MurJ and an unprecedented crystal structure of Escherichia coli MurJ in a "squeezed" form, which lacks a cavity to accommodate the substrate, mainly because of the increased proximity of transmembrane helices 2 and 8. Subsequent molecular dynamics simulations supported the hypothesis that the squeezed form is an intermediate conformation. This study fills a gap in our understanding of the Lipid II flipping mechanism.
Collapse
Affiliation(s)
- Hidetaka Kohga
- Nara Institute of Science and Technology, Ikoma, Nara 630-0192, Japan
| | - Takaharu Mori
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Yoshiki Tanaka
- Nara Institute of Science and Technology, Ikoma, Nara 630-0192, Japan
| | | | | | - Kei Fujimoto
- Nara Institute of Science and Technology, Ikoma, Nara 630-0192, Japan
| | - Lisa Fritz
- Institute for Pharmaceutical Microbiology, University Hospital Bonn, University of Bonn, Meckenheimer Allee 168, 53115 Bonn, Germany
| | - Tanja Schneider
- Institute for Pharmaceutical Microbiology, University Hospital Bonn, University of Bonn, Meckenheimer Allee 168, 53115 Bonn, Germany
| | - Tomoya Tsukazaki
- Nara Institute of Science and Technology, Ikoma, Nara 630-0192, Japan.
| |
Collapse
|
49
|
Derat E, Kamerlin SCL. Computational Advances in Protein Engineering and Enzyme Design. J Phys Chem B 2022; 126:2449-2451. [PMID: 35387452 DOI: 10.1021/acs.jpcb.2c01198] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Etienne Derat
- Institut Parisien de Chimie Moléculaire, UMR 8232 CNRS, Sorbonne Université, 75005 Paris, France
| | | |
Collapse
|
50
|
Tan C, Jung J, Kobayashi C, Torre DUL, Takada S, Sugita Y. Implementation of residue-level coarse-grained models in GENESIS for large-scale molecular dynamics simulations. PLoS Comput Biol 2022; 18:e1009578. [PMID: 35381009 PMCID: PMC9012402 DOI: 10.1371/journal.pcbi.1009578] [Citation(s) in RCA: 28] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Revised: 04/15/2022] [Accepted: 03/26/2022] [Indexed: 12/28/2022] Open
Abstract
Residue-level coarse-grained (CG) models have become one of the most popular tools in biomolecular simulations in the trade-off between modeling accuracy and computational efficiency. To investigate large-scale biological phenomena in molecular dynamics (MD) simulations with CG models, unified treatments of proteins and nucleic acids, as well as efficient parallel computations, are indispensable. In the GENESIS MD software, we implement several residue-level CG models, covering structure-based and context-based potentials for both well-folded biomolecules and intrinsically disordered regions. An amino acid residue in protein is represented as a single CG particle centered at the Cα atom position, while a nucleotide in RNA or DNA is modeled with three beads. Then, a single CG particle represents around ten heavy atoms in both proteins and nucleic acids. The input data in CG MD simulations are treated as GROMACS-style input files generated from a newly developed toolbox, GENESIS-CG-tool. To optimize the performance in CG MD simulations, we utilize multiple neighbor lists, each of which is attached to a different nonbonded interaction potential in the cell-linked list method. We found that random number generations for Gaussian distributions in the Langevin thermostat are one of the bottlenecks in CG MD simulations. Therefore, we parallelize the computations with message-passing-interface (MPI) to improve the performance on PC clusters or supercomputers. We simulate Herpes simplex virus (HSV) type 2 B-capsid and chromatin models containing more than 1,000 nucleosomes in GENESIS as examples of large-scale biomolecular simulations with residue-level CG models. This framework extends accessible spatial and temporal scales by multi-scale simulations to study biologically relevant phenomena, such as genome-scale chromatin folding or phase-separated membrane-less condensations.
Collapse
Affiliation(s)
- Cheng Tan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
| | - Jaewoon Jung
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan
| | - Chigusa Kobayashi
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
| | - Diego Ugarte La Torre
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Yuji Sugita
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo, Japan
- * E-mail:
| |
Collapse
|