1
|
Pyzer-Knapp EO, Curioni A. Advancing biomolecular simulation through exascale HPC, AI and quantum computing. Curr Opin Struct Biol 2024; 87:102826. [PMID: 38733863 DOI: 10.1016/j.sbi.2024.102826] [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: 12/15/2023] [Revised: 04/11/2024] [Accepted: 04/11/2024] [Indexed: 05/13/2024]
Abstract
Biomolecular simulation can act as both a digital microscope and a crystal ball; offering the potential for a deeper understanding of experimental observations whilst also presenting a forward-looking avenue for the in silico design and evaluation of hitherto unsynthesized compounds. Indeed, as the intricacy of our scientific inquiries has grown, so too has the computational prowess we seek to deploy in our pursuit of answers. As we enter the Exascale era, this mini-review surveys the computational landscape from both the point of view of the development of new and ever more powerful systems, and the simulations that are run on them. Moreover, as we stand on the cusp of a transformative phase in computational biology, this article offers a contemplative glance into the future, speculating on the profound implications of artificial intelligence and quantum computing for large-scale biomolecular simulations.
Collapse
|
2
|
Li J, Zhou Y, Chen SJ. Embracing exascale computing in nucleic acid simulations. Curr Opin Struct Biol 2024; 87:102847. [PMID: 38815519 PMCID: PMC11283969 DOI: 10.1016/j.sbi.2024.102847] [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: 12/17/2023] [Revised: 04/18/2024] [Accepted: 05/09/2024] [Indexed: 06/01/2024]
Abstract
This mini-review reports the recent advances in biomolecular simulations, particularly for nucleic acids, and provides the potential effects of the emerging exascale computing on nucleic acid simulations, emphasizing the need for advanced computational strategies to fully exploit this technological frontier. Specifically, we introduce recent breakthroughs in computer architectures for large-scale biomolecular simulations and review the simulation protocols for nucleic acids regarding force fields, enhanced sampling methods, coarse-grained models, and interactions with ligands. We also explore the integration of machine learning methods into simulations, which promises to significantly enhance the predictive modeling of biomolecules and the analysis of complex data generated by the exascale simulations. Finally, we discuss the challenges and perspectives for biomolecular simulations as we enter the dawning exascale computing era.
Collapse
Affiliation(s)
- Jun Li
- Department of Physics, Department of Biochemistry and Institute for Data Science and Informatics, University of Missouri, 223 Physics Bldg., Columbia, 65211, MO, USA
| | - Yuanzhe Zhou
- Department of Physics, Department of Biochemistry and Institute for Data Science and Informatics, University of Missouri, 223 Physics Bldg., Columbia, 65211, MO, USA
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry and Institute for Data Science and Informatics, University of Missouri, 223 Physics Bldg., Columbia, 65211, MO, USA.
| |
Collapse
|
3
|
Brown CM, Marrink SJ. Modeling membranes in situ. Curr Opin Struct Biol 2024; 87:102837. [PMID: 38744147 DOI: 10.1016/j.sbi.2024.102837] [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: 12/22/2023] [Revised: 03/26/2024] [Accepted: 04/23/2024] [Indexed: 05/16/2024]
Abstract
Molecular dynamics simulations of cellular membranes have come a long way-from simple model lipid bilayers to multicomponent systems capturing the crowded and complex nature of real cell membranes. In this opinionated minireview, we discuss the current challenge to simulate the dynamics of membranes in their native environment, in situ, with the prospect of reaching the level of whole cells and cell organelles using an integrative modeling framework.
Collapse
Affiliation(s)
- Chelsea M Brown
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, the Netherlands. https://twitter.com/chelseabrowncg
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, the Netherlands. s.j.marrinkrug.nl
| |
Collapse
|
4
|
Biriukov D, Vácha R. Pathways to a Shiny Future: Building the Foundation for Computational Physical Chemistry and Biophysics in 2050. ACS PHYSICAL CHEMISTRY AU 2024; 4:302-313. [PMID: 39069976 PMCID: PMC11274290 DOI: 10.1021/acsphyschemau.4c00003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/07/2024] [Revised: 03/15/2024] [Accepted: 03/18/2024] [Indexed: 07/30/2024]
Abstract
In the last quarter-century, the field of molecular dynamics (MD) has undergone a remarkable transformation, propelled by substantial enhancements in software, hardware, and underlying methodologies. In this Perspective, we contemplate the future trajectory of MD simulations and their possible look at the year 2050. We spotlight the pivotal role of artificial intelligence (AI) in shaping the future of MD and the broader field of computational physical chemistry. We outline critical strategies and initiatives that are essential for the seamless integration of such technologies. Our discussion delves into topics like multiscale modeling, adept management of ever-increasing data deluge, the establishment of centralized simulation databases, and the autonomous refinement, cross-validation, and self-expansion of these repositories. The successful implementation of these advancements requires scientific transparency, a cautiously optimistic approach to interpreting AI-driven simulations and their analysis, and a mindset that prioritizes knowledge-motivated research alongside AI-enhanced big data exploration. While history reminds us that the trajectory of technological progress can be unpredictable, this Perspective offers guidance on preparedness and proactive measures, aiming to steer future advancements in the most beneficial and successful direction.
Collapse
Affiliation(s)
- Denys Biriukov
- CEITEC
− Central European Institute of Technology, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
- National
Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
| | - Robert Vácha
- CEITEC
− Central European Institute of Technology, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
- National
Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
- Department
of Condensed Matter Physics, Faculty of Science, Masaryk University, Kotlářská 267/2, 611 37 Brno, Czech
Republic
| |
Collapse
|
5
|
Cornet J, Coulonges N, Pezeshkian W, Penissat-Mahaut M, Desgrez-Dautet H, Marrink SJ, Destainville N, Chavent M, Manghi M. There and back again: bridging meso- and nano-scales to understand lipid vesicle patterning. SOFT MATTER 2024; 20:4998-5013. [PMID: 38884641 DOI: 10.1039/d4sm00089g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2024]
Abstract
We describe a complete methodology to bridge the scales between nanoscale molecular dynamics and (micrometer) mesoscale Monte Carlo simulations in lipid membranes and vesicles undergoing phase separation, in which curving molecular species are furthermore embedded. To go from the molecular to the mesoscale, we notably appeal to physical renormalization arguments enabling us to rigorously infer the mesoscale interaction parameters from its molecular counterpart. We also explain how to deal with the physical timescales at stake at the mesoscale. Simulating the as-obtained mesoscale system enables us to equilibrate the long wavelengths of the vesicles of interest, up to the vesicle size. Conversely, we then backmap from the meso- to the nano-scale, which enables us to equilibrate in turn the short wavelengths down to the molecular length-scales. By applying our approach to the specific situation of patterning a vesicle membrane, we show that macroscopic membranes can thus be equilibrated at all length-scales in achievable computational time offering an original strategy to address the fundamental challenge of timescale in simulations of large bio-membrane systems.
Collapse
Affiliation(s)
- Julie Cornet
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, France.
| | - Nelly Coulonges
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, France.
- Institut de Pharmacologie et Biologie Structurale, Université de Toulouse, CNRS, Université Toulouse III - Paul Sabatier, 31400, Toulouse, France.
| | - Weria Pezeshkian
- Niels Bohr International Academy, Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100 Copenhagen, Denmark
| | - Maël Penissat-Mahaut
- Institut de Pharmacologie et Biologie Structurale, Université de Toulouse, CNRS, Université Toulouse III - Paul Sabatier, 31400, Toulouse, France.
| | - Hermes Desgrez-Dautet
- Laboratoire de Microbiologie et Génétique Moléculaires (LMGM), Centre de Biologie Intégrative (CBI), Université de Toulouse, CNRS, UPS, France
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | | | - Matthieu Chavent
- Institut de Pharmacologie et Biologie Structurale, Université de Toulouse, CNRS, Université Toulouse III - Paul Sabatier, 31400, Toulouse, France.
- Laboratoire de Microbiologie et Génétique Moléculaires (LMGM), Centre de Biologie Intégrative (CBI), Université de Toulouse, CNRS, UPS, France
| | - Manoel Manghi
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, France.
| |
Collapse
|
6
|
Klippenstein V, Wolf N, van der Vegt NFA. A Gauss-Newton method for iterative optimization of memory kernels for generalized Langevin thermostats in coarse-grained molecular dynamics simulations. J Chem Phys 2024; 160:204115. [PMID: 38804493 DOI: 10.1063/5.0203832] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2024] [Accepted: 05/08/2024] [Indexed: 05/29/2024] Open
Abstract
In molecular dynamics simulations, dynamically consistent coarse-grained (CG) models commonly use stochastic thermostats to model friction and fluctuations that are lost in a CG description. While Markovian, i.e., time-local, formulations of such thermostats allow for an accurate representation of diffusivities/long-time dynamics, a correct description of the dynamics on all time scales generally requires non-Markovian, i.e., non-time-local, thermostats. These thermostats typically take the form of a Generalized Langevin Equation (GLE) determined by a memory kernel. In this work, we use a Markovian embedded formulation of a position-independent GLE thermostat acting independently on each CG degree of freedom. Extracting the memory kernel of this CG model from atomistic reference data requires several approximations. Therefore, this task is best understood as an inverse problem. While our recently proposed approximate Newton scheme allows for the iterative optimization of memory kernels (IOMK), Markovian embedding remained potentially error-prone and computationally expensive. In this work, we present an IOMK-Gauss-Newton scheme (IOMK-GN) based on IOMK that allows for the direct parameterization of a Markovian embedded model.
Collapse
Affiliation(s)
- Viktor Klippenstein
- Department of Chemistry, Technical University of Darmstadt, 64287 Darmstadt, Germany
| | - Niklas Wolf
- Department of Chemistry, Technical University of Darmstadt, 64287 Darmstadt, Germany
| | - Nico F A van der Vegt
- Department of Chemistry, Technical University of Darmstadt, 64287 Darmstadt, Germany
| |
Collapse
|
7
|
Antila HS, Dixit S, Kav B, Madsen JJ, Miettinen MS, Ollila OHS. Evaluating Polarizable Biomembrane Simulations against Experiments. J Chem Theory Comput 2024; 20:4325-4337. [PMID: 38718349 PMCID: PMC11137822 DOI: 10.1021/acs.jctc.3c01333] [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: 12/06/2023] [Revised: 04/10/2024] [Accepted: 04/10/2024] [Indexed: 05/29/2024]
Abstract
Owing to the increase of available computational capabilities and the potential for providing a more accurate description, polarizable molecular dynamics force fields are gaining popularity in modeling biomolecular systems. It is, however, crucial to evaluate how much precision is truly gained with increasing cost and complexity of the simulation. Here, we leverage the NMRlipids open collaboration and Databank to assess the performance of available polarizable lipid models─the CHARMM-Drude and the AMOEBA-based parameters─against high-fidelity experimental data and compare them to the top-performing nonpolarizable models. While some improvement in the description of ion binding to membranes is observed in the most recent CHARMM-Drude parameters, and the conformational dynamics of AMOEBA-based parameters are excellent, the best nonpolarizable models tend to outperform their polarizable counterparts for each property we explored. The identified shortcomings range from inaccuracies in describing the conformational space of lipids to excessively slow conformational dynamics. Our results provide valuable insights for the further refinement of polarizable lipid force fields and for selecting the best simulation parameters for specific applications.
Collapse
Affiliation(s)
- Hanne S. Antila
- Department
of Theory and Bio-Systems, Max Planck Institute
of Colloids and Interfaces, Potsdam 14476, Germany
- Department
of Biomedicine, University of Bergen, Bergen 5020, Norway
- Computational
Biology Unit, Department of Informatics, University of Bergen, Bergen 5008, Norway
| | - Sneha Dixit
- Department
of Theory and Bio-Systems, Max Planck Institute
of Colloids and Interfaces, Potsdam 14476, Germany
| | - Batuhan Kav
- Institute
of Biological Information Processing: Structural Biochemistry (IBI-7), Forschungszentrum Jülich, Jïulich 52428, Germany
| | - Jesper J. Madsen
- Department
of Molecular Medicine, Morsani College of Medicine, University of South Florida, Tampa, Florida 33612, United States
- Center
for Global Health and Infectious Diseases Research, Global and Planetary
Health, College of Public Health, University
of South Florida, Tampa, Florida 33612, United States of America
| | - Markus S. Miettinen
- Department
of Theory and Bio-Systems, Max Planck Institute
of Colloids and Interfaces, Potsdam 14476, Germany
- Computational
Biology Unit, Department of Informatics, University of Bergen, Bergen 5008, Norway
- Department
of Chemistry, University of Bergen, Bergen 5007, Norway
| | - O. H. Samuli Ollila
- VTT Technical
Research Centre of Finland, Espoo 02044, Finland
- Institute
of Biotechnology, University of Helsinki, Helsinki 00014, Finland
| |
Collapse
|
8
|
Tiemann JKS, Szczuka M, Bouarroudj L, Oussaren M, Garcia S, Howard RJ, Delemotte L, Lindahl E, Baaden M, Lindorff-Larsen K, Chavent M, Poulain P. MDverse: Shedding Light on the Dark Matter of Molecular Dynamics Simulations. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.05.02.538537. [PMID: 37205542 PMCID: PMC10187166 DOI: 10.1101/2023.05.02.538537] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
The rise of open science and the absence of a global dedicated data repository for molecular dynamics (MD) simulations has led to the accumulation of MD files in generalist data repositories, constituting the dark matter of MD - data that is technically accessible, but neither indexed, curated, or easily searchable. Leveraging an original search strategy, we found and indexed about 250,000 files and 2,000 datasets from Zenodo, Figshare and Open Science Framework. With a focus on files produced by the Gromacs MD software, we illustrate the potential offered by the mining of publicly available MD data. We identified systems with specific molecular composition and were able to characterize essential parameters of MD simulation such as temperature and simulation length, and could identify model resolution, such as all-atom and coarse-grain. Based on this analysis, we inferred metadata to propose a search engine prototype to explore the MD data. To continue in this direction, we call on the community to pursue the effort of sharing MD data, and to report and standardize metadata to reuse this valuable matter.
Collapse
|
9
|
Wang Y, Hernandez R. Construction of Multiscale Dissipative Particle Dynamics (DPD) Models from Other Coarse-Grained Models. ACS OMEGA 2024; 9:17667-17680. [PMID: 38645334 PMCID: PMC11025104 DOI: 10.1021/acsomega.4c01868] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/26/2024] [Revised: 03/06/2024] [Accepted: 03/15/2024] [Indexed: 04/23/2024]
Abstract
We present a general scheme for converting coarse-grained models into Dissipative Particle Dynamics (DPD) models. We build the corresponding DPD models by analogy with the de novo DPD coarse-graining scheme suggested by Groot and Warren (J. Chem. Phys., 1997). Electrostatic interactions between charged DPD particles are represented though the addition of a long-range Slater Coulomb potential as suggested by González-Melchor et al. (J. Chem. Phys., 2006). The construction is illustrated by converting MARTINI models for various proteins into a DPD representation, but it not restricted to the usual potential form in the MARTINI model-viz., Lennard-Jones potentials. We further extended the DPD scheme away from the typical use of homogeneous particle sizes, therefore faithfully representing the variations in the particle sizes seen in the underlying MARTINI model. The accuracy of the resulting construction of our generalized DPD models with respect to several structural observables has been benchmarked favorably against all-atom and MARTINI models for a selected set of peptides and proteins, and variations in the scales of the coarse-graining of the water solvent.
Collapse
Affiliation(s)
- Yinhan Wang
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Rigoberto Hernandez
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
10
|
Abettan A, Nguyen MH, Ladant D, Monticelli L, Chenal A. CyaA translocation across eukaryotic cell membranes. Front Mol Biosci 2024; 11:1359408. [PMID: 38584704 PMCID: PMC10995232 DOI: 10.3389/fmolb.2024.1359408] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Accepted: 01/23/2024] [Indexed: 04/09/2024] Open
Affiliation(s)
- Amiel Abettan
- Institut Pasteur, Université de Paris Cité, CNRS UMR3528, Biochemistry of Macromolecular Interactions Unit, Paris, France
- Molecular Microbiology and Structural Biochemistry Laboratory, CNRS UMR 5086, University of Lyon, IBCP, Lyon, France
| | - Minh-Ha Nguyen
- Institut Pasteur, Université de Paris Cité, CNRS UMR3528, Biochemistry of Macromolecular Interactions Unit, Paris, France
- Université de Paris Cité, Paris, France
- Institut Pasteur, Université de Paris Cité, CNRS UMR3528, Biological NMR and HDX-MS Technological Platform, Paris, France
| | - Daniel Ladant
- Institut Pasteur, Université de Paris Cité, CNRS UMR3528, Biochemistry of Macromolecular Interactions Unit, Paris, France
- Université de Paris Cité, Paris, France
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry Laboratory, CNRS UMR 5086, University of Lyon, IBCP, Lyon, France
- Institut National de la Santé et de la Recherche Médicale (INSERM), Lyon, France
| | - Alexandre Chenal
- Institut Pasteur, Université de Paris Cité, CNRS UMR3528, Biochemistry of Macromolecular Interactions Unit, Paris, France
- Université de Paris Cité, Paris, France
| |
Collapse
|
11
|
Singh PK, Stan RC. ThermoPCD: a database of molecular dynamics trajectories of antibody-antigen complexes at physiologic and fever-range temperatures. Database (Oxford) 2024; 2024:baae015. [PMID: 38502609 PMCID: PMC10950042 DOI: 10.1093/database/baae015] [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: 10/25/2023] [Revised: 02/12/2024] [Accepted: 02/22/2024] [Indexed: 03/21/2024]
Abstract
Progression of various cancers and autoimmune diseases is associated with changes in systemic or local tissue temperatures, which may impact current therapies. The role of fever and acute inflammation-range temperatures on the stability and activity of antibodies relevant for cancers and autoimmunity is unknown. To produce molecular dynamics (MD) trajectories of immune complexes at relevant temperatures, we used the Research Collaboratory for Structural Bioinformatics (RCSB) database to identify 50 antibody:antigen complexes of interest, in addition to single antibodies and antigens, and deployed Groningen Machine for Chemical Simulations (GROMACS) to prepare and run the structures at different temperatures for 100-500 ns, in single or multiple random seeds. MD trajectories are freely available. Processed data include Protein Data Bank outputs for all files obtained every 50 ns, and free binding energy calculations for some of the immune complexes. Protocols for using the data are also available. Individual datasets contain unique DOIs. We created a web interface, ThermoPCD, as a platform to explore the data. The outputs of ThermoPCD allow the users to relate thermally-dependent changes in epitopes:paratopes interfaces to their free binding energies, or against own experimentally derived binding affinities. ThermoPCD is a free to use database of immune complexes' trajectories at different temperatures that does not require registration and allows for all the data to be available for download. Database URL: https://sites.google.com/view/thermopcd/home.
Collapse
Affiliation(s)
- Puneet K Singh
- Department of Basic Medical Science, Chonnam National University, Hwasun 58128, Republic of Korea
| | - Razvan C Stan
- Department of Basic Medical Science, Chonnam National University, Hwasun 58128, Republic of Korea
| |
Collapse
|
12
|
Beck TL, Carloni P, Asthagiri DN. All-Atom Biomolecular Simulation in the Exascale Era. J Chem Theory Comput 2024; 20:1777-1782. [PMID: 38382017 DOI: 10.1021/acs.jctc.3c01276] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2024]
Abstract
Exascale supercomputers have opened the door to dynamic simulations, facilitated by AI/ML techniques, that model biomolecular motions over unprecedented length and time scales. This new capability holds the potential to revolutionize our understanding of fundamental biological processes. Here we report on some of the major advances that were discussed at a recent CECAM workshop in Pisa, Italy, on the topic with a primary focus on atomic-level simulations. First, we highlight examples of current large-scale biomolecular simulations and the future possibilities enabled by crossing the exascale threshold. Next, we discuss challenges to be overcome in optimizing the usage of these powerful resources. Finally, we close by listing several grand challenge problems that could be investigated with this new computer architecture.
Collapse
Affiliation(s)
- Thomas L Beck
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830, United States
| | - Paolo Carloni
- INM-9/IAS-5 Computational Biomedicine, Forschungszentrum Jülich, Wilhelm-Johnen-Straße, D-54245 Jülich, Germany
- Department of Physics, RWTH Aachen University, D-52078 Aachen, Germany
| | - Dilipkumar N Asthagiri
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830, United States
| |
Collapse
|
13
|
Kiirikki AM, Antila HS, Bort LS, Buslaev P, Favela-Rosales F, Ferreira TM, Fuchs PFJ, Garcia-Fandino R, Gushchin I, Kav B, Kučerka N, Kula P, Kurki M, Kuzmin A, Lalitha A, Lolicato F, Madsen JJ, Miettinen MS, Mingham C, Monticelli L, Nencini R, Nesterenko AM, Piggot TJ, Piñeiro Á, Reuter N, Samantray S, Suárez-Lestón F, Talandashti R, Ollila OHS. Overlay databank unlocks data-driven analyses of biomolecules for all. Nat Commun 2024; 15:1136. [PMID: 38326316 PMCID: PMC10850068 DOI: 10.1038/s41467-024-45189-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2023] [Accepted: 01/17/2024] [Indexed: 02/09/2024] Open
Abstract
Tools based on artificial intelligence (AI) are currently revolutionising many fields, yet their applications are often limited by the lack of suitable training data in programmatically accessible format. Here we propose an effective solution to make data scattered in various locations and formats accessible for data-driven and machine learning applications using the overlay databank format. To demonstrate the practical relevance of such approach, we present the NMRlipids Databank-a community-driven, open-for-all database featuring programmatic access to quality-evaluated atom-resolution molecular dynamics simulations of cellular membranes. Cellular membrane lipid composition is implicated in diseases and controls major biological functions, but membranes are difficult to study experimentally due to their intrinsic disorder and complex phase behaviour. While MD simulations have been useful in understanding membrane systems, they require significant computational resources and often suffer from inaccuracies in model parameters. Here, we demonstrate how programmable interface for flexible implementation of data-driven and machine learning applications, and rapid access to simulation data through a graphical user interface, unlock possibilities beyond current MD simulation and experimental studies to understand cellular membranes. The proposed overlay databank concept can be further applied to other biomolecules, as well as in other fields where similar barriers hinder the AI revolution.
Collapse
Affiliation(s)
- Anne M Kiirikki
- University of Helsinki, Institute of Biotechnology, Helsinki, Finland
| | - Hanne S Antila
- Department of Theory and Bio-Systems, Max Planck Institute of Colloids and Interfaces, 14424, Potsdam, Germany
- Department of Biomedicine, University of Bergen, 5020, Bergen, Norway
| | - Lara S Bort
- Department of Theory and Bio-Systems, Max Planck Institute of Colloids and Interfaces, 14424, Potsdam, Germany
- University of Potsdam, Institute of Physics and Astronomy, 14476, Potsdam-Golm, Germany
| | - Pavel Buslaev
- Nanoscience Center and Department of Chemistry, University of Jyväskylä, 40014, Jyväskylä, Finland
| | - Fernando Favela-Rosales
- Departamento de Ciencias Básicas, Tecnológico Nacional de México - ITS Zacatecas Occidente, Sombrerete, 99102, Zacatecas, Mexico
| | - Tiago Mendes Ferreira
- NMR group - Institute for Physics, Martin Luther University Halle-Wittenberg, 06120, Halle (Saale), Germany
| | - Patrick F J Fuchs
- Sorbonne Université, Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules (LBM), F-75005, Paris, France
- Université Paris Cité, F-75006, Paris, France
| | - Rebeca Garcia-Fandino
- Center for Research in Biological Chemistry and Molecular Materials (CiQUS), Universidade de Santiago de Compostela, E-15782, Santiago de Compostela, Spain
| | | | - Batuhan Kav
- Institute of Biological Information Processing: Structural Biochemistry (IBI-7), Forschungszentrum Jülich, 52428, Jülich, Germany
- ariadne.ai GmbH (Germany), Häusserstraße 3, 69115, Heidelberg, Germany
| | - Norbert Kučerka
- Department of Physical Chemistry of Drugs, Faculty of Pharmacy, Comenius University Bratislava, 832 32, Bratislava, Slovakia
| | - Patrik Kula
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 542/2, CZ-16610, Prague, Czech Republic
| | - Milla Kurki
- School of Pharmacy, University of Eastern Finland, 70211, Kuopio, Finland
| | | | - Anusha Lalitha
- Institut Charles Gerhardt Montpellier (UMR CNRS 5253), Université Montpellier, Place Eugène Bataillon, 34095, Montpellier, Cedex 05, France
| | - Fabio Lolicato
- Heidelberg University Biochemistry Center, 69120, Heidelberg, Germany
- Department of Physics, University of Helsinki, FI-00014, Helsinki, Finland
| | - Jesper J Madsen
- Department of Molecular Medicine, Morsani College of Medicine, University of South Florida, 33612, Tampa, FL, USA
- Center for Global Health and Infectious Diseases Research, Global and Planetary Health, College of Public Health, University of South Florida, 33612, Tampa, FL, USA
| | - Markus S Miettinen
- Department of Theory and Bio-Systems, Max Planck Institute of Colloids and Interfaces, 14424, Potsdam, Germany
- Department of Chemistry, University of Bergen, 5007, Bergen, Norway
- Department of Informatics, Computational Biology Unit, University of Bergen, 5008, Bergen, Norway
| | - Cedric Mingham
- Hochschule Mannheim, University of Applied Sciences, 68163, Mannheim, Germany
| | - Luca Monticelli
- University of Lyon, CNRS, Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), F-69007, Lyon, France
- Institut National de la Santé et de la Recherche Médicale (INSERM), Lyon, France
| | - Ricky Nencini
- University of Helsinki, Institute of Biotechnology, Helsinki, Finland
- Division of Pharmaceutical Biosciences, Faculty of Pharmacy, University of Helsinki, 00014, Helsinki, Finland
| | - Alexey M Nesterenko
- Department of Chemistry, University of Bergen, 5007, Bergen, Norway
- Department of Informatics, Computational Biology Unit, University of Bergen, 5008, Bergen, Norway
| | - Thomas J Piggot
- Chemistry, University of Southampton, Highfield, SO17 1BJ, Southampton, UK
| | - Ángel Piñeiro
- Department of Applied Physics, Faculty of Physics, University of Santiago de Compostela, E-15782, Santiago de Compostela, Spain
| | - Nathalie Reuter
- Department of Chemistry, University of Bergen, 5007, Bergen, Norway
- Department of Informatics, Computational Biology Unit, University of Bergen, 5008, Bergen, Norway
| | - Suman Samantray
- Institute of Biological Information Processing: Structural Biochemistry (IBI-7), Forschungszentrum Jülich, 52428, Jülich, Germany
- Institute of Biotechnology, RWTH Aachen University, Worringerweg 3, 52074, Aachen, Germany
| | - Fabián Suárez-Lestón
- Center for Research in Biological Chemistry and Molecular Materials (CiQUS), Universidade de Santiago de Compostela, E-15782, Santiago de Compostela, Spain
- Department of Applied Physics, Faculty of Physics, University of Santiago de Compostela, E-15782, Santiago de Compostela, Spain
- MD.USE Innovations S.L., Edificio Emprendia, 15782, Santiago de Compostela, Spain
| | - Reza Talandashti
- Department of Chemistry, University of Bergen, 5007, Bergen, Norway
- Department of Informatics, Computational Biology Unit, University of Bergen, 5008, Bergen, Norway
| | - O H Samuli Ollila
- University of Helsinki, Institute of Biotechnology, Helsinki, Finland.
- VTT Technical Research Centre of Finland, Espoo, Finland.
| |
Collapse
|
14
|
Chen LH, Hu JN. Development of nano-delivery systems for loaded bioactive compounds: using molecular dynamics simulations. Crit Rev Food Sci Nutr 2024:1-22. [PMID: 38206576 DOI: 10.1080/10408398.2023.2301427] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
Over the past decade, a remarkable surge in the development of functional nano-delivery systems loaded with bioactive compounds for healthcare has been witnessed. Notably, the demanding requirements of high solubility, prolonged circulation, high tissue penetration capability, and strong targeting ability of nanocarriers have posed interdisciplinary research challenges to the community. While extensive experimental studies have been conducted to understand the construction of nano-delivery systems and their metabolic behavior in vivo, less is known about these molecular mechanisms and kinetic pathways during their metabolic process in vivo, and lacking effective means for high-throughput screening. Molecular dynamics (MD) simulation techniques provide a reliable tool for investigating the design of nano-delivery carriers encapsulating these functional ingredients, elucidating the synthesis, translocation, and delivery of nanocarriers. This review introduces the basic MD principles, discusses how to apply MD simulation to design nanocarriers, evaluates the ability of nanocarriers to adhere to or cross gastrointestinal mucosa, and regulates plasma proteins in vivo. Moreover, we presented the critical role of MD simulation in developing delivery systems for precise nutrition and prospects for the future. This review aims to provide insights into the implications of MD simulation techniques for designing and optimizing nano-delivery systems in the healthcare food industry.
Collapse
Affiliation(s)
- Li-Hang Chen
- SKL of Marine Food Processing & Safety Control, National Engineering Research Center of Seafood, Collaborative Innovation Center of Seafood Deep Processing, School of Food Science and Technology, Dalian Polytechnic University, Dalian, China
| | - Jiang-Ning Hu
- SKL of Marine Food Processing & Safety Control, National Engineering Research Center of Seafood, Collaborative Innovation Center of Seafood Deep Processing, School of Food Science and Technology, Dalian Polytechnic University, Dalian, China
| |
Collapse
|
15
|
Sarkar D, Lee H, Vant JW, Turilli M, Vermaas JV, Jha S, Singharoy A. Adaptive Ensemble Refinement of Protein Structures in High Resolution Electron Microscopy Density Maps with Radical Augmented Molecular Dynamics Flexible Fitting. J Chem Inf Model 2023; 63:5834-5846. [PMID: 37661856 DOI: 10.1021/acs.jcim.3c00350] [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: 09/05/2023]
Abstract
Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge-driven or artificial intelligence-driven and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary-structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models that contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly while refining them against the density maps. We introduce such an adaptive decision-making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two prototypical protein systems, adenylate kinase and carbon monoxide dehydrogenase. For these test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2-3 Å map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that are inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from the chimpanzee adenovirus. The overhead for decision making remains low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install R-MDFF on different computing environments, from local Linux-based workstations to high-performance computing environments.
Collapse
Affiliation(s)
- Daipayan Sarkar
- MSU-DOE Plant Research Laboratory, East Lansing, Michigan 48824, United States
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85281, United States
| | - Hyungro Lee
- Pacific Northwest National Laboratory, Richland, Washington 99354, United States
- Electrical & Computer Engineering, Rutgers University, New Brunswick, New Jersey 08854, United States
| | - John W Vant
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85281, United States
| | - Matteo Turilli
- Electrical & Computer Engineering, Rutgers University, New Brunswick, New Jersey 08854, United States
- Computational Science Initiative, Brookhaven National Laboratory, Upton, New York 11973, United States
| | - Josh V Vermaas
- MSU-DOE Plant Research Laboratory, East Lansing, Michigan 48824, United States
| | - Shantenu Jha
- Electrical & Computer Engineering, Rutgers University, New Brunswick, New Jersey 08854, United States
- Computational Science Initiative, Brookhaven National Laboratory, Upton, New York 11973, United States
| | - Abhishek Singharoy
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85281, United States
| |
Collapse
|
16
|
Melo MCR, Bernardi RC. Fostering discoveries in the era of exascale computing: How the next generation of supercomputers empowers computational and experimental biophysics alike. Biophys J 2023; 122:2833-2840. [PMID: 36738105 PMCID: PMC10398237 DOI: 10.1016/j.bpj.2023.01.042] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2022] [Revised: 01/24/2023] [Accepted: 01/30/2023] [Indexed: 02/05/2023] Open
Abstract
Over a century ago, physicists started broadly relying on theoretical models to guide new experiments. Soon thereafter, chemists began doing the same. Now, biological research enters a new era when experiment and theory walk hand in hand. Novel software and specialized hardware became essential to understand experimental data and propose new models. In fact, current petascale computing resources already allow researchers to reach unprecedented levels of simulation throughput to connect in silico and in vitro experiments. The reduction in cost and improved access allowed a large number of research groups to adopt supercomputing resources and techniques. Here, we outline how large-scale computing has evolved to expand decades-old research, spark new research efforts, and continuously connect simulation and observation. For instance, multiple publicly and privately funded groups have dedicated extensive resources to develop artificial intelligence tools for computational biophysics, from accelerating quantum chemistry calculations to proposing protein structure models. Moreover, advances in computer hardware have accelerated data processing from single-molecule experimental observations and simulations of chemical reactions occurring throughout entire cells. The combination of software and hardware has opened the way for exascale computing and the production of the first public exascale supercomputer, Frontier, inaugurated by the Oak Ridge National Laboratory in 2022. Ultimately, the popularization and development of computational techniques and the training of researchers to use them will only accelerate the diversification of tools and learning resources for future generations.
Collapse
Affiliation(s)
- Marcelo C R Melo
- Auburn University, Department of Physics, Auburn University, Auburn, Alabama
| | - Rafael C Bernardi
- Auburn University, Department of Physics, Auburn University, Auburn, Alabama.
| |
Collapse
|
17
|
Chaisson EH, Heberle FA, Doktorova M. Building Asymmetric Lipid Bilayers for Molecular Dynamics Simulations: What Methods Exist and How to Choose One? MEMBRANES 2023; 13:629. [PMID: 37504995 PMCID: PMC10384462 DOI: 10.3390/membranes13070629] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/07/2023] [Revised: 06/20/2023] [Accepted: 06/25/2023] [Indexed: 07/29/2023]
Abstract
The compositional asymmetry of biological membranes has attracted significant attention over the last decade. Harboring more differences from symmetric membranes than previously appreciated, asymmetric bilayers have proven quite challenging to study with familiar concepts and techniques, leaving many unanswered questions about the reach of the asymmetry effects. One particular area of active research is the computational investigation of composition- and number-asymmetric lipid bilayers with molecular dynamics (MD) simulations. Offering a high level of detail into the organization and properties of the simulated systems, MD has emerged as an indispensable tool in the study of membrane asymmetry. However, the realization that results depend heavily on the protocol used for constructing the asymmetric bilayer models has sparked an ongoing debate about how to choose the most appropriate approach. Here we discuss the underlying source of the discrepant results and review the existing methods for creating asymmetric bilayers for MD simulations. Considering the available data, we argue that each method is well suited for specific applications and hence there is no single best approach. Instead, the choice of a construction protocol-and consequently, its perceived accuracy-must be based primarily on the scientific question that the simulations are designed to address.
Collapse
Affiliation(s)
- Emily H. Chaisson
- Department of Chemistry, University of Tennessee Knoxville, Knoxville, TN 37916, USA
| | - Frederick A. Heberle
- Department of Chemistry, University of Tennessee Knoxville, Knoxville, TN 37916, USA
| | - Milka Doktorova
- Department of Molecular Physiology and Biological Physics, University of Virginia, Charlottesville, VA 22903, USA
| |
Collapse
|
18
|
Hsieh YC, Delarue M, Orland H, Koehl P. Analyzing the Geometry and Dynamics of Viral Structures: A Review of Computational Approaches Based on Alpha Shape Theory, Normal Mode Analysis, and Poisson-Boltzmann Theories. Viruses 2023; 15:1366. [PMID: 37376665 DOI: 10.3390/v15061366] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2023] [Revised: 06/05/2023] [Accepted: 06/09/2023] [Indexed: 06/29/2023] Open
Abstract
The current SARS-CoV-2 pandemic highlights our fragility when we are exposed to emergent viruses either directly or through zoonotic diseases. Fortunately, our knowledge of the biology of those viruses is improving. In particular, we have more and more structural information on virions, i.e., the infective form of a virus that includes its genomic material and surrounding protective capsid, and on their gene products. It is important to have methods that enable the analyses of structural information on such large macromolecular systems. We review some of those methods in this paper. We focus on understanding the geometry of virions and viral structural proteins, their dynamics, and their energetics, with the ambition that this understanding can help design antiviral agents. We discuss those methods in light of the specificities of those structures, mainly that they are huge. We focus on three of our own methods based on the alpha shape theory for computing geometry, normal mode analyses to study dynamics, and modified Poisson-Boltzmann theories to study the organization of ions and co-solvent and solvent molecules around biomacromolecules. The corresponding software has computing times that are compatible with the use of regular desktop computers. We show examples of their applications on some outer shells and structural proteins of the West Nile Virus.
Collapse
Affiliation(s)
- Yin-Chen Hsieh
- Institute for Arctic and Marine Biology, Department of Biosciences, Fisheries, and Economics, UiT The Arctic University of Norway, 9037 Tromso, Norway
| | - Marc Delarue
- Institut Pasteur, Université Paris-Cité and CNRS, UMR 3528, Unité Architecture et Dynamique des Macromolécules Biologiques, 75015 Paris, France
| | - Henri Orland
- Institut de Physique Théorique, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Patrice Koehl
- Department of Computer Science, University of California, Davis, CA 95616, USA
| |
Collapse
|
19
|
Sinha S, Pindi C, Ahsan M, Arantes PR, Palermo G. Machines on Genes through the Computational Microscope. J Chem Theory Comput 2023; 19:1945-1964. [PMID: 36947696 PMCID: PMC10104023 DOI: 10.1021/acs.jctc.2c01313] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/24/2023]
Abstract
Macromolecular machines acting on genes are at the core of life's fundamental processes, including DNA replication and repair, gene transcription and regulation, chromatin packaging, RNA splicing, and genome editing. Here, we report the increasing role of computational biophysics in characterizing the mechanisms of "machines on genes", focusing on innovative applications of computational methods and their integration with structural and biophysical experiments. We showcase how state-of-the-art computational methods, including classical and ab initio molecular dynamics to enhanced sampling techniques, and coarse-grained approaches are used for understanding and exploring gene machines for real-world applications. As this review unfolds, advanced computational methods describe the biophysical function that is unseen through experimental techniques, accomplishing the power of the "computational microscope", an expression coined by Klaus Schulten to highlight the extraordinary capability of computer simulations. Pushing the frontiers of computational biophysics toward a pragmatic representation of large multimegadalton biomolecular complexes is instrumental in bridging the gap between experimentally obtained macroscopic observables and the molecular principles playing at the microscopic level. This understanding will help harness molecular machines for medical, pharmaceutical, and biotechnological purposes.
Collapse
|
20
|
Pezeshkian W, Grünewald F, Narykov O, Lu S, Arkhipova V, Solodovnikov A, Wassenaar TA, Marrink SJ, Korkin D. Molecular architecture and dynamics of SARS-CoV-2 envelope by integrative modeling. Structure 2023; 31:492-503.e7. [PMID: 36870335 DOI: 10.1016/j.str.2023.02.006] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2022] [Revised: 11/15/2022] [Accepted: 02/07/2023] [Indexed: 03/06/2023]
Abstract
Despite tremendous efforts, the exact structure of SARS-CoV-2 and related betacoronaviruses remains elusive. SARS-CoV-2 envelope is a key structural component of the virion that encapsulates viral RNA. It is composed of three structural proteins, spike, membrane (M), and envelope, which interact with each other and with the lipids acquired from the host membranes. Here, we developed and applied an integrative multi-scale computational approach to model the envelope structure of SARS-CoV-2 with near atomistic detail, focusing on studying the dynamic nature and molecular interactions of its most abundant, but largely understudied, M protein. The molecular dynamics simulations allowed us to test the envelope stability under different configurations and revealed that the M dimers agglomerated into large, filament-like, macromolecular assemblies with distinct molecular patterns. These results are in good agreement with current experimental data, demonstrating a generic and versatile approach to model the structure of a virus de novo.
Collapse
Affiliation(s)
- Weria Pezeshkian
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, 9747AG Groningen, the Netherlands; Niels Bohr International Academy, Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100 Copenhagen, Denmark
| | - Fabian Grünewald
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, 9747AG Groningen, the Netherlands
| | - Oleksandr Narykov
- Department of Computer Science, Worcester Polytechnic Institute, Worcester, MA 01609, USA
| | - Senbao Lu
- Bioinformatics and Computational Biology Program, Worcester Polytechnic Institute, Worcester, MA 01609, USA
| | | | | | - Tsjerk A Wassenaar
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, 9747AG Groningen, the Netherlands; Institute for Life Science and Technology, Hanze University of Applied Sciences, 9747AS Groningen, the Netherlands
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, 9747AG Groningen, the Netherlands.
| | - Dmitry Korkin
- Department of Computer Science, Worcester Polytechnic Institute, Worcester, MA 01609, USA; Bioinformatics and Computational Biology Program, Worcester Polytechnic Institute, Worcester, MA 01609, USA.
| |
Collapse
|
21
|
Stevens JA, Grünewald F, van Tilburg PAM, König M, Gilbert BR, Brier TA, Thornburg ZR, Luthey-Schulten Z, Marrink SJ. Molecular dynamics simulation of an entire cell. Front Chem 2023; 11:1106495. [PMID: 36742032 PMCID: PMC9889929 DOI: 10.3389/fchem.2023.1106495] [Citation(s) in RCA: 35] [Impact Index Per Article: 35.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2022] [Accepted: 01/09/2023] [Indexed: 01/19/2023] Open
Abstract
The ultimate microscope, directed at a cell, would reveal the dynamics of all the cell's components with atomic resolution. In contrast to their real-world counterparts, computational microscopes are currently on the brink of meeting this challenge. In this perspective, we show how an integrative approach can be employed to model an entire cell, the minimal cell, JCVI-syn3A, at full complexity. This step opens the way to interrogate the cell's spatio-temporal evolution with molecular dynamics simulations, an approach that can be extended to other cell types in the near future.
Collapse
Affiliation(s)
- Jan A. Stevens
- Molecular Dynamics Group, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, Netherlands
| | - Fabian Grünewald
- Molecular Dynamics Group, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, Netherlands
| | - P. A. Marco van Tilburg
- Molecular Dynamics Group, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, Netherlands
| | - Melanie König
- Molecular Dynamics Group, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, Netherlands
| | - Benjamin R. Gilbert
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Champaign, IL, United States
| | - Troy A. Brier
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Champaign, IL, United States
| | - Zane R. Thornburg
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Champaign, IL, United States
| | - Zaida Luthey-Schulten
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Champaign, IL, United States
| | - Siewert J. Marrink
- Molecular Dynamics Group, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, Netherlands,*Correspondence: Siewert J. Marrink,
| |
Collapse
|
22
|
Koehl P, Akopyan A, Edelsbrunner H. Computing the Volume, Surface Area, Mean, and Gaussian Curvatures of Molecules and Their Derivatives. J Chem Inf Model 2023; 63:973-985. [PMID: 36638318 PMCID: PMC9930125 DOI: 10.1021/acs.jcim.2c01346] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
Geometry is crucial in our efforts to comprehend the structures and dynamics of biomolecules. For example, volume, surface area, and integrated mean and Gaussian curvature of the union of balls representing a molecule are used to quantify its interactions with the water surrounding it in the morphometric implicit solvent models. The Alpha Shape theory provides an accurate and reliable method for computing these geometric measures. In this paper, we derive homogeneous formulas for the expressions of these measures and their derivatives with respect to the atomic coordinates, and we provide algorithms that implement them into a new software package, AlphaMol. The only variables in these formulas are the interatomic distances, making them insensitive to translations and rotations. AlphaMol includes a sequential algorithm and a parallel algorithm. In the parallel version, we partition the atoms of the molecule of interest into 3D rectangular blocks, using a kd-tree algorithm. We then apply the sequential algorithm of AlphaMol to each block, augmented by a buffer zone to account for atoms whose ball representations may partially cover the block. The current parallel version of AlphaMol leads to a 20-fold speed-up compared to an independent serial implementation when using 32 processors. For instance, it takes 31 s to compute the geometric measures and derivatives of each atom in a viral capsid with more than 26 million atoms on 32 Intel processors running at 2.7 GHz. The presence of the buffer zones, however, leads to redundant computations, which ultimately limit the impact of using multiple processors. AlphaMol is available as an OpenSource software.
Collapse
Affiliation(s)
- Patrice Koehl
- Department
of Computer Science, University of California, Davis, California95616, United States,
| | | | | |
Collapse
|
23
|
Sarkar D, Kulke M, Vermaas JV. LongBondEliminator: A Molecular Simulation Tool to Remove Ring Penetrations in Biomolecular Simulation Systems. Biomolecules 2023; 13:biom13010107. [PMID: 36671493 PMCID: PMC9856086 DOI: 10.3390/biom13010107] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2022] [Revised: 12/30/2022] [Accepted: 01/01/2023] [Indexed: 01/07/2023] Open
Abstract
We develop a workflow, implemented as a plugin to the molecular visualization program VMD, that can fix ring penetrations with minimal user input. LongBondEliminator, detects ring piercing artifacts by the long, strained bonds that are the local minimum energy conformation during minimization for some assembled simulation system. The LongBondEliminator tool then automatically treats regions near these long bonds using multiple biases applied through NAMD. By combining biases implemented through the collective variables module, density-based forces, and alchemical techniques in NAMD, LongBondEliminator will iteratively alleviate long bonds found within molecular simulation systems. Through three concrete examples with increasing complexity, a lignin polymer, an viral capsid assembly, and a large, highly glycosylated protein aggrecan, we demonstrate the utility for this method in eliminating ring penetrations from classical MD simulation systems. The tool is available via gitlab as a VMD plugin, and has been developed to be generically useful across a variety of biomolecular simulations.
Collapse
|
24
|
Corey RA, Baaden M, Chavent M. A brief history of visualizing membrane systems in molecular dynamics simulations. FRONTIERS IN BIOINFORMATICS 2023; 3:1149744. [PMID: 37213533 PMCID: PMC10196259 DOI: 10.3389/fbinf.2023.1149744] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2023] [Accepted: 03/13/2023] [Indexed: 05/23/2023] Open
Abstract
Understanding lipid dynamics and function, from the level of single, isolated molecules to large assemblies, is more than ever an intensive area of research. The interactions of lipids with other molecules, particularly membrane proteins, are now extensively studied. With advances in the development of force fields for molecular dynamics simulations (MD) and increases in computational resources, the creation of realistic and complex membrane systems is now common. In this perspective, we will review four decades of the history of molecular dynamics simulations applied to membranes and lipids through the prism of molecular graphics.
Collapse
Affiliation(s)
- R. A. Corey
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| | - M. Baaden
- Centre Nationale de la Recherche Scientifique, Laboratoire de Biochimie Théorique, Université Paris Cité, Paris, France
| | - M. Chavent
- Institut de Pharmacologie et Biologie Structurale, CNRS, Université de Toulouse, Toulouse, France
- *Correspondence: M. Chavent,
| |
Collapse
|
25
|
Hadjidemetriou K, Kaur S, Cassidy CK, Zhang P. Mechanisms of E. coli chemotaxis signaling pathways visualized using cryoET and computational approaches. Biochem Soc Trans 2022; 50:1595-1605. [PMID: 36421737 PMCID: PMC9788364 DOI: 10.1042/bst20220191] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2022] [Revised: 11/08/2022] [Accepted: 11/11/2022] [Indexed: 11/25/2022]
Abstract
Chemotaxis signaling pathways enable bacteria to sense and respond to their chemical environment and, in some species, are critical for lifestyle processes such as biofilm formation and pathogenesis. The signal transduction underlying chemotaxis behavior is mediated by large, highly ordered protein complexes known as chemosensory arrays. For nearly two decades, cryo-electron tomography (cryoET) has been used to image chemosensory arrays, providing an increasingly detailed understanding of their structure and function. In this mini-review, we provide an overview of the use of cryoET to study chemosensory arrays, including imaging strategies, key results, and outstanding questions. We further discuss the application of molecular modeling and simulation techniques to complement structure determination efforts and provide insight into signaling mechanisms. We close the review with a brief outlook, highlighting promising future directions for the field.
Collapse
Affiliation(s)
| | - Satinder Kaur
- Diamond Light Source, Harwell Science and Innovation Campus, Didcot OX11 0DE, U.K
| | - C. Keith Cassidy
- Diamond Light Source, Harwell Science and Innovation Campus, Didcot OX11 0DE, U.K
| | - Peijun Zhang
- Diamond Light Source, Harwell Science and Innovation Campus, Didcot OX11 0DE, U.K
- Division of Structural Biology, Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, U.K
- Chinese Academy of Medical Sciences Oxford Institute, University of Oxford, Oxford OX3 7BN, U.K
| |
Collapse
|
26
|
Pantano S. Back and forth modeling through biological scales. Biochem Biophys Res Commun 2022; 633:39-41. [DOI: 10.1016/j.bbrc.2022.09.037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2022] [Accepted: 09/09/2022] [Indexed: 11/06/2022]
|
27
|
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
|
28
|
Kampfrath M, Staritzbichler R, Hernández GP, Rose AS, Tiemann JKS, Scheuermann G, Wiegreffe D, Hildebrand PW. MDsrv: visual sharing and analysis of molecular dynamics simulations. Nucleic Acids Res 2022; 50:W483-W489. [PMID: 35639717 PMCID: PMC9252803 DOI: 10.1093/nar/gkac398] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2022] [Revised: 05/02/2022] [Accepted: 05/05/2022] [Indexed: 11/24/2022] Open
Abstract
Molecular dynamics simulation is a proven technique for computing and visualizing the time-resolved motion of macromolecules at atomic resolution. The MDsrv is a tool that streams MD trajectories and displays them interactively in web browsers without requiring advanced skills, facilitating interactive exploration and collaborative visual analysis. We have now enhanced the MDsrv to further simplify the upload and sharing of MD trajectories and improve their online viewing and analysis. With the new instance, the MDsrv simplifies the creation of sessions, which allows the exchange of MD trajectories with preset representations and perspectives. An important innovation is that the MDsrv can now access and visualize trajectories from remote datasets, which greatly expands its applicability and use, as the data no longer needs to be accessible on a local server. In addition, initial analyses such as sequence or structure alignments, distance measurements, or RMSD calculations have been implemented, which optionally support visual analysis. Finally, based on Mol*, MDsrv now provides faster and more efficient visualization of even large trajectories compared to its predecessor tool NGL.
Collapse
Affiliation(s)
- Michelle Kampfrath
- Image and Signal Processing Group, Department of Computer Science, Leipzig University, Augustusplatz 10, 04109 Leipzig, Germany
| | - René Staritzbichler
- Institute for Medical Physics and Biophysics, Medical Faculty, Leipzig University, Härtelstraße 16-18, 04107 Leipzig, Germany
| | - Guillermo Pérez Hernández
- Charité Universitätsmedizin Berlin, corporate member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Institute of Medical Physics and Biophysics, Berlin, Germany
| | | | - Johanna K S Tiemann
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N., Denmark
| | - Gerik Scheuermann
- Image and Signal Processing Group, Department of Computer Science, Leipzig University, Augustusplatz 10, 04109 Leipzig, Germany
| | - Daniel Wiegreffe
- Image and Signal Processing Group, Department of Computer Science, Leipzig University, Augustusplatz 10, 04109 Leipzig, Germany
| | - Peter W Hildebrand
- Institute for Medical Physics and Biophysics, Medical Faculty, Leipzig University, Härtelstraße 16-18, 04107 Leipzig, Germany.,Charité Universitätsmedizin Berlin, corporate member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Institute of Medical Physics and Biophysics, Berlin, Germany.,Berlin Institute of Health, 10178 Berlin, Germany
| |
Collapse
|