1
|
Deeks HM, Zinovjev K, Barnoud J, Mulholland AJ, van der Kamp MW, Glowacki DR. Free energy along drug-protein binding pathways interactively sampled in virtual reality. Sci Rep 2023; 13:16665. [PMID: 37794083 PMCID: PMC10551034 DOI: 10.1038/s41598-023-43523-x] [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: 06/30/2023] [Accepted: 09/25/2023] [Indexed: 10/06/2023] Open
Abstract
We describe a two-step approach for combining interactive molecular dynamics in virtual reality (iMD-VR) with free energy (FE) calculation to explore the dynamics of biological processes at the molecular level. We refer to this combined approach as iMD-VR-FE. Stage one involves using a state-of-the-art 'human-in-the-loop' iMD-VR framework to generate a diverse range of protein-ligand unbinding pathways, benefitting from the sophistication of human spatial and chemical intuition. Stage two involves using the iMD-VR-sampled pathways as initial guesses for defining a path-based reaction coordinate from which we can obtain a corresponding free energy profile using FE methods. To investigate the performance of the method, we apply iMD-VR-FE to investigate the unbinding of a benzamidine ligand from a trypsin protein. The binding free energy calculated using iMD-VR-FE is similar for each pathway, indicating internal consistency. Moreover, the resulting free energy profiles can distinguish energetic differences between pathways corresponding to various protein-ligand conformations (e.g., helping to identify pathways that are more favourable) and enable identification of metastable states along the pathways. The two-step iMD-VR-FE approach offers an intuitive way for researchers to test hypotheses for candidate pathways in biomolecular systems, quickly obtaining both qualitative and quantitative insight.
Collapse
Affiliation(s)
- Helen M Deeks
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK
| | - Kirill Zinovjev
- Departamento de Química Física, Universidad de Valencia, 46100, Burjassot, Spain
- School of Biochemistry, University of Bristol, Bristol, BS8 1TD, UK
| | - Jonathan Barnoud
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK
- CiTIUS | Centro Singular de Investigación en Tecnoloxías Intelixentes da USC, Rúa de Jenaro de la Fuente, s/n, 15705, Santiago de Compostela, A Coruña, Spain
| | - Adrian J Mulholland
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK
| | - Marc W van der Kamp
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK.
- School of Biochemistry, University of Bristol, Bristol, BS8 1TD, UK.
| | - David R Glowacki
- CiTIUS | Centro Singular de Investigación en Tecnoloxías Intelixentes da USC, Rúa de Jenaro de la Fuente, s/n, 15705, Santiago de Compostela, A Coruña, Spain.
| |
Collapse
|
2
|
Käser S, Vazquez-Salazar LI, Meuwly M, Töpfer K. Neural network potentials for chemistry: concepts, applications and prospects. DIGITAL DISCOVERY 2023; 2:28-58. [PMID: 36798879 PMCID: PMC9923808 DOI: 10.1039/d2dd00102k] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/23/2022] [Accepted: 12/20/2022] [Indexed: 12/24/2022]
Abstract
Artificial Neural Networks (NN) are already heavily involved in methods and applications for frequent tasks in the field of computational chemistry such as representation of potential energy surfaces (PES) and spectroscopic predictions. This perspective provides an overview of the foundations of neural network-based full-dimensional potential energy surfaces, their architectures, underlying concepts, their representation and applications to chemical systems. Methods for data generation and training procedures for PES construction are discussed and means for error assessment and refinement through transfer learning are presented. A selection of recent results illustrates the latest improvements regarding accuracy of PES representations and system size limitations in dynamics simulations, but also NN application enabling direct prediction of physical results without dynamics simulations. The aim is to provide an overview for the current state-of-the-art NN approaches in computational chemistry and also to point out the current challenges in enhancing reliability and applicability of NN methods on a larger scale.
Collapse
Affiliation(s)
- Silvan Käser
- Department of Chemistry, University of Basel Klingelbergstrasse 80 CH-4056 Basel Switzerland
| | | | - Markus Meuwly
- Department of Chemistry, University of Basel Klingelbergstrasse 80 CH-4056 Basel Switzerland
| | - Kai Töpfer
- Department of Chemistry, University of Basel Klingelbergstrasse 80 CH-4056 Basel Switzerland
| |
Collapse
|
3
|
Chu Q, Luo KH, Chen D. Exploring Complex Reaction Networks Using Neural Network-Based Molecular Dynamics Simulation. J Phys Chem Lett 2022; 13:4052-4057. [PMID: 35522222 DOI: 10.1021/acs.jpclett.2c00647] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Ab initio molecular dynamics (AIMD) is an established method for revealing the reactive dynamics of complex systems. However, the high computational cost of AIMD restricts the explorable length and time scales. Here, we develop a fundamentally different approach using molecular dynamics simulations powered by a neural network potential to investigate complex reaction networks. This potential is trained via a workflow combining AIMD and interactive molecular dynamics in virtual reality to accelerate the sampling of rare reactive processes. A panoramic visualization of the complex reaction networks for decomposition of a novel high explosive (ICM-102) is achieved without any predefined reaction coordinates. The study leads to the discovery of new pathways that would be difficult to uncover if established methods were employed. These results highlight the power of neural network-based molecular dynamics simulations in exploring complex reaction mechanisms under extreme conditions at the ab initio level, pushing the limit of theoretical and computational chemistry toward the realism and fidelity of experiments.
Collapse
Affiliation(s)
- Qingzhao Chu
- State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China
- Explosion Protection and Emergency Disposal Technology Engineering Research Center of the Ministry of Education, Beijing 100081, China
| | - Kai H Luo
- Department of Mechanical Engineering, University College London, Torrington Place, London WC1E 7JE, U.K
| | - Dongping Chen
- State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China
- Explosion Protection and Emergency Disposal Technology Engineering Research Center of the Ministry of Education, Beijing 100081, China
| |
Collapse
|
4
|
Ji L, Li Y, Wang J, Ning A, Zhang N, Liang S, He J, Zhang T, Qu Z, Gao J. Community Reaction Network Reduction for Constructing a Coarse-Grained Representation of Combustion Reaction Mechanisms. J Chem Inf Model 2022; 62:2352-2364. [PMID: 35442657 DOI: 10.1021/acs.jcim.2c00240] [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]
Abstract
A community-reaction network reduction (CNR) approach is presented for mechanism reduction on the basis of a network-based community detection technique, a concept related to pre-equilibrium in chemical kinetics. In this method, the detailed combustion mechanism is first transformed into a weighted network, in which communities of species that have dense inner connections under the critical ignition conditions are identified. By analyzing the community partitions in different regions, we determine the effective functional groups and driving processes. Then, a skeletal model for the overall mechanism is deduced according to the network centrality data, including transition pathway identification and reaction-path flux. The CNR method is illustrated on the hydrogen autoignition system which has been extensively investigated, and a new reduced mechanism involving seven processes is proposed. Dynamics simulations employing the present CNR model show that the computed ignition time and distribution of major species on a wide range of temperature and pressure conditions are in accord with the experiments and results from other methods.
Collapse
Affiliation(s)
- Lin Ji
- Department of Chemistry, Capital Normal University, Beijing 100048, China
| | - Yue Li
- Department of Chemistry, Capital Normal University, Beijing 100048, China
| | - Jie Wang
- Department of Chemistry, Capital Normal University, Beijing 100048, China
| | - An Ning
- Department of Chemistry, Capital Normal University, Beijing 100048, China
| | - Naixin Zhang
- Department of Chemistry, Capital Normal University, Beijing 100048, China
| | - Shengyao Liang
- Department of Chemistry, Capital Normal University, Beijing 100048, China
| | - Jiyun He
- Department of Chemistry, Capital Normal University, Beijing 100048, China
| | - Tianyu Zhang
- Department of Chemistry, Capital Normal University, Beijing 100048, China
| | - Zexing Qu
- Laboratory of Theoretical and Computational Chemistry, Jilin University, Changchun 130015, China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China.,Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
5
|
Chen X, Liu M, Gao J. CARNOT: a Fragment-Based Direct Molecular Dynamics and Virtual-Reality Simulation Package for Reactive Systems. J Chem Theory Comput 2022; 18:1297-1313. [PMID: 35129348 DOI: 10.1021/acs.jctc.1c01032] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Traditionally, the study of reaction mechanisms of complex reaction systems such as combustion has been performed on an individual basis by optimizations of transition structure and minimum energy path or by reaction dynamics trajectory calculations for one elementary reaction at a time. It is effective, but time-consuming, whereas important and unexpected processes could have been missed. In this article, we present a direct molecular dynamics (DMD) approach and a virtual-reality simulation program, CARNOT, in which plausible chemical reactions are simulated simultaneously at finite temperature and pressure conditions. A key concept of the present ab initio molecular dynamics method is to partition a large, chemically reactive system into molecular fragments that can be adjusted on the fly of a DMD simulation. The theory represents an extension of the explicit polarization method to reactive events, called ReX-Pol. We propose a highest-and-lowest adapted-spin approximation to define the local spins of individual fragments, rather than treating the entire system by a delocalized wave function. Consequently, the present ab initio DMD can be applied to reactive systems consisting of an arbitrarily varying number of closed and open-shell fragments such as free radicals, zwitterions, and separate ions found in combustion and other reactions. A graph-data structure algorithm was incorporated in CARNOT for the analysis of reaction networks, suitable for reaction mechanism reduction. Employing the PW91 density functional theory and the 6-31+G(d) basis set, the capabilities of the CARNOT program were illustrated by a combustion reaction, consisting of 28 650 atoms, and by reaction network analysis that revealed a range of mechanistic and dynamical events. The method may be useful for applications to other types of complex reactions.
Collapse
Affiliation(s)
- Xin Chen
- Peking University Shenzhen Graduate School, Shenzhen, Guangdong 581055, China.,Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong 581055, China
| | - Meiyi Liu
- Peking University Shenzhen Graduate School, Shenzhen, Guangdong 581055, China.,Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong 581055, China
| | - Jiali Gao
- Peking University Shenzhen Graduate School, Shenzhen, Guangdong 581055, China.,Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong 581055, China.,Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
6
|
Human Action Recognition in Smart Cultural Tourism Based on Fusion Techniques of Virtual Reality and SOM Neural Network. COMPUTATIONAL INTELLIGENCE AND NEUROSCIENCE 2021; 2021:3495203. [PMID: 34899891 PMCID: PMC8664517 DOI: 10.1155/2021/3495203] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/26/2021] [Revised: 11/12/2021] [Accepted: 11/18/2021] [Indexed: 11/26/2022]
Abstract
Smart cultural tourism is the development trend of the future tourism industry. Virtual reality is an important tool to realize smart tourism. The reality of virtual reality mainly comes from human-computer interaction, which is closely related to human action recognition technology. Therefore, the research takes human action recognition as the research direction, uses a self-organizing mapping network (SOM) neural network to extract the key frame of action video, combines it with multi-feature vector method to recognize human action, and compares the recognition rate and user satisfaction of different recognition methods. The results show that the recognition rate of multi-feature voting human action recognition algorithm based on SOM neural network is 93.68% on UT-Kinect action, 59.06% on MSRDailyActivity3D, and the overall action recognition time is only 3.59 s. Within six months, the total profit of human-computer interactive virtual reality tourism project with SOM neural network multi-eigenvector as the core algorithm reached 422,000 yuan, and 88% of users expressed satisfaction after use. It shows that the proposed method has a good recognition rate and can give users effective feedback in time. It is hoped that this research has a certain reference value in promoting the development of human motion recognition technology.
Collapse
|
7
|
Shannon RJ, Deeks HM, Burfoot E, Clark E, Jones AJ, Mulholland AJ, Glowacki DR. Exploring human-guided strategies for reaction network exploration: Interactive molecular dynamics in virtual reality as a tool for citizen scientists. J Chem Phys 2021; 155:154106. [PMID: 34686059 DOI: 10.1063/5.0062517] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The emerging fields of citizen science and gamification reformulate scientific problems as games or puzzles to be solved. Through engaging the wider non-scientific community, significant breakthroughs may be made by analyzing citizen-gathered data. In parallel, recent advances in virtual reality (VR) technology are increasingly being used within a scientific context and the burgeoning field of interactive molecular dynamics in VR (iMD-VR) allows users to interact with dynamical chemistry simulations in real time. Here, we demonstrate the utility of iMD-VR as a medium for gamification of chemistry research tasks. An iMD-VR "game" was designed to encourage users to explore the reactivity of a particular chemical system, and a cohort of 18 participants was recruited to playtest this game as part of a user study. The reaction game encouraged users to experiment with making chemical reactions between a propyne molecule and an OH radical, and "molecular snapshots" from each game session were then compiled and used to map out reaction pathways. The reaction network generated by users was compared to existing literature networks demonstrating that users in VR capture almost all the important reaction pathways. Further comparisons between humans and an algorithmic method for guiding molecular dynamics show that through using citizen science to explore these kinds of chemical problems, new approaches and strategies start to emerge.
Collapse
Affiliation(s)
- Robin J Shannon
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Helen M Deeks
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Eleanor Burfoot
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Edward Clark
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Alex J Jones
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | | | - David R Glowacki
- ArtSci Foundation International, 5th floor Mariner House, Bristol, BS1 4QD, United Kingdom
| |
Collapse
|
8
|
Ceriotti M, Clementi C, Anatole von Lilienfeld O. Machine learning meets chemical physics. J Chem Phys 2021; 154:160401. [PMID: 33940847 DOI: 10.1063/5.0051418] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Over recent years, the use of statistical learning techniques applied to chemical problems has gained substantial momentum. This is particularly apparent in the realm of physical chemistry, where the balance between empiricism and physics-based theory has traditionally been rather in favor of the latter. In this guest Editorial for the special topic issue on "Machine Learning Meets Chemical Physics," a brief rationale is provided, followed by an overview of the topics covered. We conclude by making some general remarks.
Collapse
Affiliation(s)
- Michele Ceriotti
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Cecilia Clementi
- Department of Physics, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | | |
Collapse
|