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
|
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
|
3
|
Acharya S, Bagchi B. Diffusion in a two-dimensional energy landscape in the presence of dynamical correlations and validity of random walk model. Phys Rev E 2023; 107:024127. [PMID: 36932553 DOI: 10.1103/physreve.107.024127] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2022] [Accepted: 01/26/2023] [Indexed: 06/18/2023]
Abstract
Diffusion in a multidimensional energy surface with minima and barriers is a problem of importance in statistical mechanics and it also has wide applications, such as protein folding. To understand it in such a system, we carry out theory and simulations of a tagged particle moving on a two-dimensional periodic potential energy surface, both in the presence and absence of noise. Langevin dynamics simulations at multiple temperatures are carried out to obtain the diffusion coefficient of a solute particle. Friction is varied from zero to large values. Diffusive motion emerges in the limit of a long time, even in the absence of noise. Noise destroys the correlations and increases diffusion at small friction. Diffusion thus exhibits a nonmonotonic friction dependence at the intermediate value of the damping, ultimately converging to our theoretically predicted value. The latter is obtained using the well-established relationship between diffusion and random walk. An excellent agreement is obtained between theory and simulations in the high-friction limit but not so in the intermediate regime. We explain the deviation in the low- to intermediate-friction regime using the modified random walk theory. The rate of escape from one cell to another is obtained from the multidimensional rate theory of Langer. We find that enhanced dimensionality plays an important role. To quantify the effects of noise on the potential-imposed coherence on the trajectories, we calculate the Lyapunov exponent. At small friction values, the Lyapunov exponent mimics the friction dependence of the rate.
Collapse
Affiliation(s)
- Subhajit Acharya
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012, India
| | - Biman Bagchi
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012, India
| |
Collapse
|
4
|
Understanding How Cells Probe the World: A Preliminary Step towards Modeling Cell Behavior? Int J Mol Sci 2023; 24:ijms24032266. [PMID: 36768586 PMCID: PMC9916635 DOI: 10.3390/ijms24032266] [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: 12/19/2022] [Revised: 01/16/2023] [Accepted: 01/20/2023] [Indexed: 01/26/2023] Open
Abstract
Cell biologists have long aimed at quantitatively modeling cell function. Recently, the outstanding progress of high-throughput measurement methods and data processing tools has made this a realistic goal. The aim of this paper is twofold: First, to suggest that, while much progress has been done in modeling cell states and transitions, current accounts of environmental cues driving these transitions remain insufficient. There is a need to provide an integrated view of the biochemical, topographical and mechanical information processed by cells to take decisions. It might be rewarding in the near future to try to connect cell environmental cues to physiologically relevant outcomes rather than modeling relationships between these cues and internal signaling networks. The second aim of this paper is to review exogenous signals that are sensed by living cells and significantly influence fate decisions. Indeed, in addition to the composition of the surrounding medium, cells are highly sensitive to the properties of neighboring surfaces, including the spatial organization of anchored molecules and substrate mechanical and topographical properties. These properties should thus be included in models of cell behavior. It is also suggested that attempts at cell modeling could strongly benefit from two research lines: (i) trying to decipher the way cells encode the information they retrieve from environment analysis, and (ii) developing more standardized means of assessing the quality of proposed models, as was done in other research domains such as protein structure prediction.
Collapse
|
5
|
Cordier BA, Sawaya NPD, Guerreschi GG, McWeeney SK. Biology and medicine in the landscape of quantum advantages. J R Soc Interface 2022; 19:20220541. [PMID: 36448288 PMCID: PMC9709576 DOI: 10.1098/rsif.2022.0541] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Accepted: 11/04/2022] [Indexed: 12/03/2022] Open
Abstract
Quantum computing holds substantial potential for applications in biology and medicine, spanning from the simulation of biomolecules to machine learning methods for subtyping cancers on the basis of clinical features. This potential is encapsulated by the concept of a quantum advantage, which is contingent on a reduction in the consumption of a computational resource, such as time, space or data. Here, we distill the concept of a quantum advantage into a simple framework to aid researchers in biology and medicine pursuing the development of quantum applications. We then apply this framework to a wide variety of computational problems relevant to these domains in an effort to (i) assess the potential of practical advantages in specific application areas and (ii) identify gaps that may be addressed with novel quantum approaches. In doing so, we provide an extensive survey of the intersection of biology and medicine with the current landscape of quantum algorithms and their potential advantages. While we endeavour to identify specific computational problems that may admit practical advantages throughout this work, the rapid pace of change in the fields of quantum computing, classical algorithms and biological research implies that this intersection will remain highly dynamic for the foreseeable future.
Collapse
Affiliation(s)
- Benjamin A. Cordier
- Department of Medical Informatics and Clinical Epidemiology, Oregon Health and Science University, Portland, OR 97202, USA
| | | | | | - Shannon K. McWeeney
- Department of Medical Informatics and Clinical Epidemiology, Oregon Health and Science University, Portland, OR 97202, USA
- Knight Cancer Institute, Oregon Health and Science University, Portland, OR 97202, USA
- Oregon Clinical and Translational Research Institute, Oregon Health and Science University, Portland, OR 97202, USA
| |
Collapse
|
6
|
Mascart C, Scarella G, Reynaud-Bouret P, Muzy A. Scalability of Large Neural Network Simulations via Activity Tracking With Time Asynchrony and Procedural Connectivity. Neural Comput 2022; 34:1915-1943. [PMID: 35896155 DOI: 10.1162/neco_a_01524] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2021] [Accepted: 04/27/2022] [Indexed: 11/04/2022]
Abstract
We present a new algorithm to efficiently simulate random models of large neural networks satisfying the property of time asynchrony. The model parameters (average firing rate, number of neurons, synaptic connection probability, and postsynaptic duration) are of the order of magnitude of a small mammalian brain or of human brain areas. Through the use of activity tracking and procedural connectivity (dynamical regeneration of synapses), computational and memory complexities of this algorithm are proved to be theoretically linear with the number of neurons. These results are experimentally validated by sequential simulations of millions of neurons and billions of synapses running in a few minutes using a single thread of an equivalent desktop computer.
Collapse
Affiliation(s)
| | - Gilles Scarella
- Université Côte d'Azur, CNRS, I3S, France.,Université Côte d'Azur, CNRS, LJAD, 06103 Nice, France
| | | | | |
Collapse
|
7
|
Wang Y, Zhang C, Tang K, Wang X. En route for molecular dynamics simulation of a living cell. FUNDAMENTAL RESEARCH 2022. [DOI: 10.1016/j.fmre.2022.06.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022] Open
|
8
|
Gupta C, Sarkar D, Tieleman DP, Singharoy A. The ugly, bad, and good stories of large-scale biomolecular simulations. Curr Opin Struct Biol 2022; 73:102338. [PMID: 35245737 DOI: 10.1016/j.sbi.2022.102338] [Citation(s) in RCA: 28] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Revised: 11/29/2021] [Accepted: 01/24/2022] [Indexed: 12/20/2022]
Abstract
Molecular modeling of large biomolecular assemblies exemplifies a disruptive area holding both promises and contentions. Propelled by peta and exascale computing, several simulation methodologies have now matured into user-friendly tools that are successfully employed for modeling viruses, membranous nano-constructs, and key pieces of the genetic machinery. We present three unifying biophysical themes that emanate from some of the most recent multi-million atom simulation endeavors. Despite connecting molecular changes with phenotypic outcomes, the quality measures of these simulations remain questionable. We discuss the existing and upcoming strategies for constructing representative ensembles of large systems, how new computing technologies will boost this area, and make a point that integrative modeling guided by experimental data is the future of biomolecular computations.
Collapse
Affiliation(s)
- Chitrak Gupta
- School of Molecular Sciences, Center for Applied Structural Discovery, Arizona State University at Tempe, Tempe, AZ, 85282, USA; Biodesign Institute, Tempe, AZ, 85281, USA. https://twitter.com/ChitrakGupta2
| | - Daipayan Sarkar
- School of Molecular Sciences, Center for Applied Structural Discovery, Arizona State University at Tempe, Tempe, AZ, 85282, USA; MSU-DOE Plant Research Laboratory, Michigan State University, East Lansing, MI, 48824-1319, USA. https://twitter.com/17Dsarkar
| | - D Peter Tieleman
- Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada.
| | - Abhishek Singharoy
- School of Molecular Sciences, Center for Applied Structural Discovery, Arizona State University at Tempe, Tempe, AZ, 85282, USA; Biodesign Institute, Tempe, AZ, 85281, USA.
| |
Collapse
|
9
|
Bongrand P. Is There a Need for a More Precise Description of Biomolecule Interactions to Understand Cell Function? Curr Issues Mol Biol 2022; 44:505-525. [PMID: 35723321 PMCID: PMC8929073 DOI: 10.3390/cimb44020035] [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: 11/25/2021] [Revised: 01/15/2022] [Accepted: 01/17/2022] [Indexed: 11/16/2022] Open
Abstract
An important goal of biological research is to explain and hopefully predict cell behavior from the molecular properties of cellular components. Accordingly, much work was done to build extensive “omic” datasets and develop theoretical methods, including computer simulation and network analysis to process as quantitatively as possible the parameters contained in these resources. Furthermore, substantial effort was made to standardize data presentation and make experimental results accessible to data scientists. However, the power and complexity of current experimental and theoretical tools make it more and more difficult to assess the capacity of gathered parameters to support optimal progress in our understanding of cell function. The purpose of this review is to focus on biomolecule interactions, the interactome, as a specific and important example, and examine the limitations of the explanatory and predictive power of parameters that are considered as suitable descriptors of molecular interactions. Recent experimental studies on important cell functions, such as adhesion and processing of environmental cues for decision-making, support the suggestion that it should be rewarding to complement standard binding properties such as affinity and kinetic constants, or even force dependence, with less frequently used parameters such as conformational flexibility or size of binding molecules.
Collapse
Affiliation(s)
- Pierre Bongrand
- Lab Adhesion and Inflammation (LAI), Inserm UMR 1067, Cnrs UMR 7333, Aix-Marseille Université UM 61, Marseille 13009, France
| |
Collapse
|
10
|
Puech PH, Bongrand P. Mechanotransduction as a major driver of cell behaviour: mechanisms, and relevance to cell organization and future research. Open Biol 2021; 11:210256. [PMID: 34753321 PMCID: PMC8586914 DOI: 10.1098/rsob.210256] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Accepted: 10/18/2021] [Indexed: 01/04/2023] Open
Abstract
How do cells process environmental cues to make decisions? This simple question is still generating much experimental and theoretical work, at the border of physics, chemistry and biology, with strong implications in medicine. The purpose of mechanobiology is to understand how biochemical and physical cues are turned into signals through mechanotransduction. Here, we review recent evidence showing that (i) mechanotransduction plays a major role in triggering signalling cascades following cell-neighbourhood interaction; (ii) the cell capacity to continually generate forces, and biomolecule properties to undergo conformational changes in response to piconewton forces, provide a molecular basis for understanding mechanotransduction; and (iii) mechanotransduction shapes the guidance cues retrieved by living cells and the information flow they generate. This includes the temporal and spatial properties of intracellular signalling cascades. In conclusion, it is suggested that the described concepts may provide guidelines to define experimentally accessible parameters to describe cell structure and dynamics, as a prerequisite to take advantage of recent progress in high-throughput data gathering, computer simulation and artificial intelligence, in order to build a workable, hopefully predictive, account of cell signalling networks.
Collapse
Affiliation(s)
- Pierre-Henri Puech
- Lab Adhesion and Inflammation (LAI), Inserm UMR 1067, CNRS UMR 7333, Aix-Marseille Université UM61, Marseille, France
| | - Pierre Bongrand
- Lab Adhesion and Inflammation (LAI), Inserm UMR 1067, CNRS UMR 7333, Aix-Marseille Université UM61, Marseille, France
| |
Collapse
|
11
|
|
12
|
Figon F, Casas J. The integrative biology of pigment organelles, a quantum chemical approach. Integr Comp Biol 2021; 61:1490-1501. [PMID: 33940609 DOI: 10.1093/icb/icab045] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Coloration is a complex phenotypic trait involving both physical and chemical processes at a multiscale level, from molecules to tissues. Pigments, whose main property is to absorb specific wavelengths of visible light, are usually deposited in specialized organelles or complex matrices comprising proteins, metals, ions and redox compounds, among others. By modulating electronic properties and stability, interactions between pigments and these molecular actors can lead to color tuning. Furthermore, pigments are not only important for visual effects but also provide other critical functions, such as detoxification and antiradical activity. Hence, integrative studies of pigment organelles are required to understand how pigments interact with their cellular environment. In this review, we show how quantum chemistry, a computational method that models the molecular and optical properties of pigments, has provided key insights into the mechanisms by which pigment properties, from color to reactivity, are modulated by their organellar environment. These results allow to rationalize and to predict the way pigments behave in supramolecular complexes, up to the complete modelling of pigment organelles. We also discuss the main limitations of quantum chemistry, emphasizing the need for carrying experimental work with identical vigor. We finally suggest that taking into account the ecology of pigments (i.e. how they interact with these various other cellular components and at higher organizational levels) will lead to a greater understanding of how and why animals are vividly and variably colored, two fundamental questions in organismal biology.
Collapse
Affiliation(s)
- Florent Figon
- Institut de Recherche sur la Biologie de l'Insecte, UMR 7261, CNRS-Université de Tours, 37200 Tours, France
| | - Jérôme Casas
- Institut de Recherche sur la Biologie de l'Insecte, UMR 7261, CNRS-Université de Tours, 37200 Tours, France
| |
Collapse
|
13
|
|
14
|
A three dimensional computer model of urothelium and bladder cancer initiation, progress and collective invasion. INFORMATICS IN MEDICINE UNLOCKED 2021. [DOI: 10.1016/j.imu.2021.100750] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022] Open
|