1
|
Zhang P, Gardini AT, Xu X, Parrinello M. Intramolecular and Water Mediated Tautomerism of Solvated Glycine. J Chem Inf Model 2024. [PMID: 38620066 DOI: 10.1021/acs.jcim.4c00273] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2024]
Abstract
Understanding tautomerism and characterizing solvent effects on the dynamic processes pose significant challenges. Using enhanced-sampling molecular dynamics based on state-of-the-art deep learning potentials, we investigated the tautomeric equilibria of glycine in water. We observed that the tautomerism between neutral and zwitterionic glycine can occur through both intramolecular and intermolecular proton transfers. The latter proceeds involving a contact anionic-glycine-hydronium ion pair or separate cationic-glycine-hydroxide ion pair. These pathways with comparable barriers contribute almost equally to the reaction flux.
Collapse
Affiliation(s)
- Pengchao Zhang
- Center for Combustion Energy, Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, China
- Atomistic Simulations, Italian Institute of Technology, Genova 16152, Italy
| | - Axel Tosello Gardini
- Atomistic Simulations, Italian Institute of Technology, Genova 16152, Italy
- Department of Materials Science, Università di Milano-Bicocca, 20126 Milano, Italy
| | - Xuefei Xu
- Center for Combustion Energy, Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, Genova 16152, Italy
| |
Collapse
|
2
|
Müllender L, Rizzi A, Parrinello M, Carloni P, Mandelli D. Effective data-driven collective variables for free energy calculations from metadynamics of paths. PNAS Nexus 2024; 3:pgae159. [PMID: 38665160 PMCID: PMC11044970 DOI: 10.1093/pnasnexus/pgae159] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/05/2023] [Accepted: 04/04/2024] [Indexed: 04/28/2024]
Abstract
A variety of enhanced sampling (ES) methods predict multidimensional free energy landscapes associated with biological and other molecular processes as a function of a few selected collective variables (CVs). The accuracy of these methods is crucially dependent on the ability of the chosen CVs to capture the relevant slow degrees of freedom of the system. For complex processes, finding such CVs is the real challenge. Machine learning (ML) CVs offer, in principle, a solution to handle this problem. However, these methods rely on the availability of high-quality datasets-ideally incorporating information about physical pathways and transition states-which are difficult to access, therefore greatly limiting their domain of application. Here, we demonstrate how these datasets can be generated by means of ES simulations in trajectory space via the metadynamics of paths algorithm. The approach is expected to provide a general and efficient way to generate efficient ML-based CVs for the fast prediction of free energy landscapes in ES simulations. We demonstrate our approach with two numerical examples, a 2D model potential and the isomerization of alanine dipeptide, using deep targeted discriminant analysis as our ML-based CV of choice.
Collapse
Affiliation(s)
- Lukas Müllender
- Department of Applied Physics, Science for Life Laboratory, KTH Royal Institute of Technology, SE-171 21 Solna, Sweden
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
- Department of Physics, RWTH Aachen University, 52062 Aachen, Germany
| | - Andrea Rizzi
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
- Atomistic Simulations, Italian Institute of Technology, 16163 Genova, Italy
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, 16163 Genova, Italy
| | - Paolo Carloni
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
- Department of Physics, RWTH Aachen University, 52062 Aachen, Germany
- Universitätsklinikum, RWTH Aachen University, 52062 Aachen, Germany
| | - Davide Mandelli
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
| |
Collapse
|
3
|
Ray D, Parrinello M. Data-driven classification of ligand unbinding pathways. Proc Natl Acad Sci U S A 2024; 121:e2313542121. [PMID: 38412121 DOI: 10.1073/pnas.2313542121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2023] [Accepted: 01/26/2024] [Indexed: 02/29/2024] Open
Abstract
Studying the pathways of ligand-receptor binding is essential to understand the mechanism of target recognition by small molecules. The binding free energy and kinetics of protein-ligand complexes can be computed using molecular dynamics (MD) simulations, often in quantitative agreement with experiments. However, only a qualitative picture of the ligand binding/unbinding paths can be obtained through a conventional analysis of the MD trajectories. Besides, the higher degree of manual effort involved in analyzing pathways limits its applicability in large-scale drug discovery. Here, we address this limitation by introducing an automated approach for analyzing molecular transition paths with a particular focus on protein-ligand dissociation. Our method is based on the dynamic time-warping algorithm, originally designed for speech recognition. We accurately classified molecular trajectories using a very generic descriptor set of contacts or distances. Our approach outperforms manual classification by distinguishing between parallel dissociation channels, within the pathways identified by visual inspection. Most notably, we could compute exit-path-specific ligand-dissociation kinetics. The unbinding timescale along the fastest path agrees with the experimental residence time, providing a physical interpretation to our entirely data-driven protocol. In combination with appropriate enhanced sampling algorithms, this technique can be used for the initial exploration of ligand-dissociation pathways as well as for calculating path-specific thermodynamic and kinetic properties.
Collapse
Affiliation(s)
- Dhiman Ray
- Simulations Research Line, Italian Institute of Technology, Via Enrico Melen 83, Genova GE 16152, Italy
| | - Michele Parrinello
- Simulations Research Line, Italian Institute of Technology, Via Enrico Melen 83, Genova GE 16152, Italy
| |
Collapse
|
4
|
Yang M, Trizio E, Parrinello M. Structure and polymerization of liquid sulfur across the λ-transition. Chem Sci 2024; 15:3382-3392. [PMID: 38425540 PMCID: PMC10902632 DOI: 10.1039/d3sc06282a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Accepted: 01/18/2024] [Indexed: 03/02/2024] Open
Abstract
The anomalous λ-transition of liquid sulfur, which is supposed to be related to the transformation of eight-membered sulfur rings into long polymeric chains, has attracted considerable attention. However, a detailed description of the underlying dynamical polymerization process is still missing. Here, we study the structures and the mechanism of the polymerization processes of liquid sulfur across the λ-transition as well as its reverse process of formation of the rings. We do so by performing ab initio-quality molecular dynamics simulations thanks to a combination of machine learning potentials and state-of-the-art enhanced sampling techniques. With our approach, we obtain structural results that are in good agreement with the experiments and we report precious dynamical insights into the mechanisms involved in the process.
Collapse
Affiliation(s)
- Manyi Yang
- Atomistic Simulations, Italian Institute of Technology 16156 Genova Italy
| | - Enrico Trizio
- Atomistic Simulations, Italian Institute of Technology 16156 Genova Italy
- Department of Materials Science, Università di Milano-Bicocca 20126 Milano Italy
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology 16156 Genova Italy
| |
Collapse
|
5
|
Ruiz Munevar M, Rizzi V, Portioli C, Vidossich P, Cao E, Parrinello M, Cancedda L, De Vivo M. Cation Chloride Cotransporter NKCC1 Operates through a Rocking-Bundle Mechanism. J Am Chem Soc 2024; 146:552-566. [PMID: 38146212 PMCID: PMC10786066 DOI: 10.1021/jacs.3c10258] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2023] [Revised: 12/05/2023] [Accepted: 12/06/2023] [Indexed: 12/27/2023]
Abstract
The sodium, potassium, and chloride cotransporter 1 (NKCC1) plays a key role in tightly regulating ion shuttling across cell membranes. Lately, its aberrant expression and function have been linked to numerous neurological disorders and cancers, making it a novel and highly promising pharmacological target for therapeutic interventions. A better understanding of how NKCC1 dynamically operates would therefore have broad implications for ongoing efforts toward its exploitation as a therapeutic target through its modulation. Based on recent structural data on NKCC1, we reveal conformational motions that are key to its function. Using extensive deep-learning-guided atomistic simulations of NKCC1 models embedded into the membrane, we captured complex dynamical transitions between alternate open conformations of the inner and outer vestibules of the cotransporter and demonstrated that NKCC1 has water-permeable states. We found that these previously undefined conformational transitions occur via a rocking-bundle mechanism characterized by the cooperative angular motion of transmembrane helices (TM) 4 and 9, with the contribution of the extracellular tip of TM 10. We found these motions to be critical in modulating ion transportation and in regulating NKCC1's water transporting capabilities. Specifically, we identified interhelical dynamical contacts between TM 10 and TM 6, which we functionally validated through mutagenesis experiments of 4 new targeted NKCC1 mutants. We conclude showing that those 4 residues are highly conserved in most Na+-dependent cation chloride cotransporters (CCCs), which highlights their critical mechanistic implications, opening the way to new strategies for NKCC1's function modulation and thus to potential drug action on selected CCCs.
Collapse
Affiliation(s)
- Manuel
José Ruiz Munevar
- Laboratory
of Molecular Modelling & Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, Genoa 16163, Italy
| | - Valerio Rizzi
- Biomolecular
& Pharmaceutical Modelling Group, Université
de Genève, Rue Michel-Servet 1, Geneva CH-1211 4, Switzerland
| | - Corinne Portioli
- Laboratory
of Nanotechnology for Precision Medicine, Istituto Italiano di Tecnologia, Via Morego 30, Genoa 16163, Italy
- Laboratory
of Brain Development and Disease, Istituto
Italiano di Tecnologia, Via Morego 30, Genoa 16163, Italy
| | - Pietro Vidossich
- Laboratory
of Molecular Modelling & Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, Genoa 16163, Italy
| | - Erhu Cao
- Department
of Biochemistry, University of Utah School
of Medicine, Salt Lake City, Utah 84112-5650, United States
| | - Michele Parrinello
- Laboratory
of Atomistic Simulations, Istituto Italiano
di Tecnologia, Via Morego 30, Genoa 16163, Italy
| | - Laura Cancedda
- Laboratory
of Brain Development and Disease, Istituto
Italiano di Tecnologia, Via Morego 30, Genoa 16163, Italy
| | - Marco De Vivo
- Laboratory
of Molecular Modelling & Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, Genoa 16163, Italy
| |
Collapse
|
6
|
Bonati L, Polino D, Pizzolitto C, Biasi P, Eckert R, Reitmeier S, Schlögl R, Parrinello M. The role of dynamics in heterogeneous catalysis: Surface diffusivity and N 2 decomposition on Fe(111). Proc Natl Acad Sci U S A 2023; 120:e2313023120. [PMID: 38060558 PMCID: PMC10723053 DOI: 10.1073/pnas.2313023120] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2023] [Accepted: 10/18/2023] [Indexed: 12/17/2023] Open
Abstract
Dynamics has long been recognized to play an important role in heterogeneous catalytic processes. However, until recently, it has been impossible to study their dynamical behavior at industry-relevant temperatures. Using a combination of machine learning potentials and advanced simulation techniques, we investigate the cleavage of the N[Formula: see text] triple bond on the Fe(111) surface. We find that at low temperatures our results agree with the well-established picture. However, if we increase the temperature to reach operando conditions, the surface undergoes a global dynamical change and the step structure of the Fe(111) surface is destabilized. The catalytic sites, traditionally associated with this surface, appear and disappear continuously. Our simulations illuminate the danger of extrapolating low-temperature results to operando conditions and indicate that the catalytic activity can only be inferred from calculations that take dynamics fully into account. More than that, they show that it is the transition to this highly fluctuating interfacial environment that drives the catalytic process.
Collapse
Affiliation(s)
- Luigi Bonati
- Atomistic Simulations, Italian Institute of Technology, Genova16152, Italy
| | - Daniela Polino
- Department of Innovative Technologies, University of Applied Sciences and Arts of Southern Switzerland, Lugano6962, Switzerland
| | - Cristina Pizzolitto
- Basic Research, Research and Development Division, Casale SA, Lugano6900, Switzerland
| | - Pierdomenico Biasi
- Basic Research, Research and Development Division, Casale SA, Lugano6900, Switzerland
| | - Rene Eckert
- BU Catalysts, R&D Syngas Applications, Clariant Produkte (Deutschland) GmbH, Munich83052, Germany
| | - Stephan Reitmeier
- BU Catalysts, R&D Syngas Applications, Clariant Produkte (Deutschland) GmbH, Munich83052, Germany
| | - Robert Schlögl
- Department of Inorganic Chemistry, Fritz-Haber Institute of the Max-Planck-Society, Berlin14195, Germany
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, Genova16152, Italy
| |
Collapse
|
7
|
Rizzi A, Carloni P, Parrinello M. Free energies at QM accuracy from force fields via multimap targeted estimation. Proc Natl Acad Sci U S A 2023; 120:e2304308120. [PMID: 37931103 PMCID: PMC10655219 DOI: 10.1073/pnas.2304308120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2023] [Accepted: 09/25/2023] [Indexed: 11/08/2023] Open
Abstract
Accurate predictions of ligand binding affinities would greatly accelerate the first stages of drug discovery campaigns. However, using highly accurate interatomic potentials based on quantum mechanics (QM) in free energy methods has been so far largely unfeasible due to their prohibitive computational cost. Here, we present an efficient method to compute QM free energies from simulations using cheap reference potentials, such as force fields (FFs). This task has traditionally been out of reach due to the slow convergence of computing the correction from the FF to the QM potential. To overcome this bottleneck, we generalize targeted free energy methods to employ multiple maps-implemented with normalizing flow neural networks (NNs)-that maximize the overlap between the distributions. Critically, the method requires neither a separate expensive training phase for the NNs nor samples from the QM potential. We further propose a one-epoch learning policy to efficiently avoid overfitting, and we combine our approach with enhanced sampling strategies to overcome the pervasive problem of poor convergence due to slow degrees of freedom. On the drug-like molecules in the HiPen dataset, the method accelerates the calculation of the free energy difference of switching from an FF to a DFTB3 potential by three orders of magnitude compared to standard free energy perturbation and by a factor of eight compared to previously published nonequilibrium calculations. Our results suggest that our method, in combination with efficient QM/MM calculations, may be used in lead optimization campaigns in drug discovery and to study protein-ligand molecular recognition processes.
Collapse
Affiliation(s)
- Andrea Rizzi
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich52428, Germany
- Atomistic Simulations, Italian Institute of Technology, Genova16163, Italy
| | - Paolo Carloni
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich52428, Germany
- Department of Physics and Universitätsklinikum, RWTH Aachen University, Aachen52074, Germany
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, Genova16163, Italy
| |
Collapse
|
8
|
Abstract
Metadynamics is a popular enhanced sampling algorithm for computing the free energy landscape of rare events by using molecular dynamics simulation. Ten years ago, Tiwary and Parrinello introduced the infrequent metadynamics approach for calculating the kinetics of transitions across free energy barriers. Since then, metadynamics-based methods for obtaining rate constants have attracted significant attention in computational molecular science. Such methods have been applied to study a wide range of problems, including protein-ligand binding, protein folding, conformational transitions, chemical reactions, catalysis, and nucleation. Here, we review the principles of elucidating kinetics from metadynamics-like approaches, subsequent methodological developments in this area, and successful applications on chemical, biological, and material systems. We also highlight the challenges of reconstructing accurate kinetics from enhanced sampling simulations and the scope of future developments.
Collapse
Affiliation(s)
- Dhiman Ray
- Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, 16152 Genova, Italy
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, 16152 Genova, Italy
| |
Collapse
|
9
|
Bonati L, Trizio E, Rizzi A, Parrinello M. A unified framework for machine learning collective variables for enhanced sampling simulations: mlcolvar. J Chem Phys 2023; 159:014801. [PMID: 37409767 DOI: 10.1063/5.0156343] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2023] [Accepted: 06/13/2023] [Indexed: 07/07/2023] Open
Abstract
Identifying a reduced set of collective variables is critical for understanding atomistic simulations and accelerating them through enhanced sampling techniques. Recently, several methods have been proposed to learn these variables directly from atomistic data. Depending on the type of data available, the learning process can be framed as dimensionality reduction, classification of metastable states, or identification of slow modes. Here, we present mlcolvar, a Python library that simplifies the construction of these variables and their use in the context of enhanced sampling through a contributed interface to the PLUMED software. The library is organized modularly to facilitate the extension and cross-contamination of these methodologies. In this spirit, we developed a general multi-task learning framework in which multiple objective functions and data from different simulations can be combined to improve the collective variables. The library's versatility is demonstrated through simple examples that are prototypical of realistic scenarios.
Collapse
Affiliation(s)
- Luigi Bonati
- Atomistic Simulations, Italian Institute of Technology, 16156 Genova, Italy
| | - Enrico Trizio
- Atomistic Simulations, Italian Institute of Technology, 16156 Genova, Italy
- Department of Materials Science, Università di Milano-Bicocca, 20126 Milano, Italy
| | - Andrea Rizzi
- Atomistic Simulations, Italian Institute of Technology, 16156 Genova, Italy
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, 16156 Genova, Italy
| |
Collapse
|
10
|
Ray D, Trizio E, Parrinello M. Deep learning collective variables from transition path ensemble. J Chem Phys 2023; 158:2891484. [PMID: 37212403 DOI: 10.1063/5.0148872] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2023] [Accepted: 05/03/2023] [Indexed: 05/23/2023] Open
Abstract
The study of the rare transitions that take place between long lived metastable states is a major challenge in molecular dynamics simulations. Many of the methods suggested to address this problem rely on the identification of the slow modes of the system, which are referred to as collective variables. Recently, machine learning methods have been used to learn the collective variables as functions of a large number of physical descriptors. Among many such methods, Deep Targeted Discriminant Analysis has proven to be useful. This collective variable is built from data harvested from short unbiased simulations in the metastable basins. Here, we enrich the set of data on which the Deep Targeted Discriminant Analysis collective variable is built by adding data from the transition path ensemble. These are collected from a number of reactive trajectories obtained using the On-the-fly Probability Enhanced Sampling flooding method. The collective variables thus trained lead to more accurate sampling and faster convergence. The performance of these new collective variables is tested on a number of representative examples.
Collapse
Affiliation(s)
- Dhiman Ray
- Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, Genoa GE 16153, Italy
| | - Enrico Trizio
- Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, Genoa GE 16153, Italy
- Department of Materials Science, Università di Milano-Bicocca, Milano 20126, Italy
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, Genoa GE 16153, Italy
| |
Collapse
|
11
|
Raucci U, Sanchez DM, Martínez TJ, Parrinello M. Enhanced Sampling Aided Design of Molecular Photoswitches. J Am Chem Soc 2022; 144:19265-19271. [PMID: 36222799 DOI: 10.1021/jacs.2c04419] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Advances in the evolving field of atomistic simulations promise important insights for the design and fundamental understanding of novel molecular photoswitches. Here, we use state-of-the-art enhanced simulation techniques to unravel the complex, multistep chemistry of donor-acceptor Stenhouse adducts (DASAs). Our reaction discovery workflow consists of enhanced sampling for efficient chemical space exploration, refinement of newly observed pathways with more accurate ab initio electronic structure calculations, and structural modifications to introduce design principles within future generations of DASAs. We showcase our discovery workflow by not only recovering the full photoswitching mechanism of DASA but also predicting a plethora of new plausible thermal pathways and suggesting a way for their experimental validation. Furthermore, we illustrate the tunability of these newly discovered reactions, leading to a potential avenue for controlling DASA dynamics through multiple external stimuli. Overall, these insights could offer alternative routes to increase the efficiency and control of DASA's photoswitching mechanism, providing new elements to design more complex light-responsive materials.
Collapse
Affiliation(s)
| | - David M Sanchez
- Department of Chemistry, Stanford University, Stanford, California94305, United States.,SLAC National Accelerator Laboratory, Stanford PULSE Institute, Menlo Park, California94025, United States
| | - Todd J Martínez
- Department of Chemistry, Stanford University, Stanford, California94305, United States.,SLAC National Accelerator Laboratory, Stanford PULSE Institute, Menlo Park, California94025, United States
| | | |
Collapse
|
12
|
Abstract
We introduce a novel enhanced sampling approach named on-the-fly probability enhanced sampling (OPES) flooding for calculating the kinetics of rare events from atomistic molecular dynamics simulation. This method is derived from the OPES approach [Invernizzi and Parrinello, J. Phys. Chem. Lett. 2020, 11, 7, 2731-2736], which has been recently developed for calculating converged free energy surfaces for complex systems. In this paper, we describe the theoretical details of the OPES flooding technique and demonstrate the application on three systems of increasing complexity: barrier crossing in a two-dimensional double-well potential, conformational transition in the alanine dipeptide in the gas phase, and the folding and unfolding of the chignolin polypeptide in an aqueous environment. From extensive tests, we show that the calculation of accurate kinetics not only requires the transition state to be bias-free, but the amount of bias deposited should also not exceed the effective barrier height measured along the chosen collective variables. In this vein, the possibility of computing rates from biasing suboptimal order parameters has also been explored. Furthermore, we describe the choice of optimum parameter combinations for obtaining accurate results from limited computational effort.
Collapse
Affiliation(s)
- Dhiman Ray
- Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, 16152 Genova, Italy
| | - Narjes Ansari
- Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, 16152 Genova, Italy
| | - Valerio Rizzi
- Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, 16152 Genova, Italy.,School of Pharmaceutical Sciences and Institute of Pharmaceutical Sciences of Western Switzerland (ISPSO), University of Geneva, Rue Michel Servet 1, 1211 Genève 4, Switzerland
| | | | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, 16152 Genova, Italy
| |
Collapse
|
13
|
Novelli P, Bonati L, Pontil M, Parrinello M. Characterizing Metastable States with the Help of Machine Learning. J Chem Theory Comput 2022; 18:5195-5202. [PMID: 35920063 DOI: 10.1021/acs.jctc.2c00393] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Present-day atomistic simulations generate long trajectories of ever more complex systems. Analyzing these data, discovering metastable states, and uncovering their nature are becoming increasingly challenging. In this paper, we first use the variational approach to conformation dynamics to discover the slowest dynamical modes of the simulations. This allows the different metastable states of the system to be located and organized hierarchically. The physical descriptors that characterize metastable states are discovered by means of a machine learning method. We show in the cases of two proteins, chignolin and bovine pancreatic trypsin inhibitor, how such analysis can be effortlessly performed in a matter of seconds. Another strength of our approach is that it can be applied to the analysis of both unbiased and biased simulations.
Collapse
Affiliation(s)
- Pietro Novelli
- Computational Statistics and Machine Learning, Italian Institute of Technology, Via Enrico Melen 83, 16142 Genoa, Italy
| | - Luigi Bonati
- Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, 16142 Genoa, Italy
| | - Massimiliano Pontil
- Computational Statistics and Machine Learning, Italian Institute of Technology, Via Enrico Melen 83, 16142 Genoa, Italy.,Department of Computer Science, University College London, London WC1E 6BT, United Kingdom
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, 16142 Genoa, Italy
| |
Collapse
|
14
|
Bjelobrk Z, Rajagopalan AK, Mendels D, Karmakar T, Parrinello M, Mazzotti M. Solubility of Organic Salts in Solvent-Antisolvent Mixtures: A Combined Experimental and Molecular Dynamics Simulations Approach. J Chem Theory Comput 2022; 18:4952-4959. [PMID: 35833664 PMCID: PMC9367008 DOI: 10.1021/acs.jctc.2c00304] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We combine molecular dynamics simulations with experiments to estimate solubilities of an organic salt in complex growth environments. We predict the solubility by simulations of the growth and dissolution of ions at the crystal surface kink sites at different solution concentrations. Thereby, the solubility is identified as the solution's salt concentration, where the energy of the ion pair dissolved in solution equals the energy of the ion pair crystallized at the kink sites. The simulation methodology is demonstrated for the case of anhydrous sodium acetate crystallized from various solvent-antisolvent mixtures. To validate the predicted solubilities, we have measured the solubilities of sodium acetate in-house, using an experimental setup and measurement protocol that guarantees moisture-free conditions, which is key for a hygroscopic compound like sodium acetate. We observe excellent agreement between the experimental and the computationally evaluated solubilities for sodium acetate in different solvent-antisolvent mixtures. Given the agreement and the rich data the simulations produce, we can use them to complement experimental tasks, which in turn will reduce time and capital in the design of complicated industrial crystallization processes of organic salts.
Collapse
Affiliation(s)
- Zoran Bjelobrk
- Institute of Energy and Process Engineering, ETH Zürich, Zürich CH-8092, Switzerland
| | - Ashwin Kumar Rajagopalan
- Department of Chemical Engineering, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Dan Mendels
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Tarak Karmakar
- Department of Chemistry, Indian Institute of Technology, Delhi, Hauz Khas, New Delhi 110016, India
| | - Michele Parrinello
- Istituto Italiano di Tecnologia (IIT), Via Morego, 30, Genova 16163, Italy
| | - Marco Mazzotti
- Institute of Energy and Process Engineering, ETH Zürich, Zürich CH-8092, Switzerland
| |
Collapse
|
15
|
Abstract
![]()
In adaptive-bias
enhanced sampling methods, a bias potential is
added to the system to drive transitions between metastable states.
The bias potential is a function of a few collective variables and
is gradually modified according to the underlying free energy surface.
We show that when the collective variables are suboptimal, there is
an exploration–convergence tradeoff, and one must choose between
a quickly converging bias that will lead to fewer transitions or a
slower to converge bias that can explore the phase space more efficiently
but might require a much longer time to produce an accurate free energy
estimate. The recently proposed on-the-fly probability enhanced sampling
(OPES) method focuses on fast convergence, but there are cases where
fast exploration is preferred instead. For this reason, we introduce
a new variant of the OPES method that focuses on quickly escaping
metastable states at the expense of convergence speed. We illustrate
the benefits of this approach in prototypical systems and show that
it outperforms the popular metadynamics method.
Collapse
|
16
|
Abstract
The long time scale associated with ligand residence times renders their computation challenging. Therefore, the influence of factors like solvation and steric hindrance on residence times is not fully understood. Here, we demonstrate in a set of model host-guest systems that the recently developed Gaussian mixture based enhanced sampling allows residence times to be computed and enables an understanding of their unbinding mechanism. We observe that guest unbinding often proceeds via a series of intermediate states that can be labeled by the number of water molecules present in the binding cavity. In several cases the residence time is correlated to the water trapping times in the cavity.
Collapse
Affiliation(s)
- Jayashrita Debnath
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland.,Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy
| | | |
Collapse
|
17
|
Abstract
Over the last few decades, enhanced sampling methods have been continuously improved. Here, we exploit this progress and propose a modular workflow for blind reaction discovery and determination of reaction paths. In a three-step strategy, at first we use a collective variable derived from spectral graph theory in conjunction with the explore variant of the on-the-fly probability enhanced sampling method to drive reaction discovery runs. Once different chemical products are determined, we construct an ad-hoc neural network-based collective variable to improve sampling, and finally we refine the results using the free energy perturbation theory and a more accurate Hamiltonian. We apply this strategy to both intramolecular and intermolecular reactions. Our workflow requires minimal user input and extends the power of ab initio molecular dynamics to explore and characterize the reaction space.
Collapse
Affiliation(s)
- Umberto Raucci
- Italian Institute of Technology, Via E. Melen 83, 16152, Genova, Italy
| | - Valerio Rizzi
- Italian Institute of Technology, Via E. Melen 83, 16152, Genova, Italy
| | | |
Collapse
|
18
|
Myung CW, Hirshberg B, Parrinello M. Prediction of a Supersolid Phase in High-Pressure Deuterium. Phys Rev Lett 2022; 128:045301. [PMID: 35148160 DOI: 10.1103/physrevlett.128.045301] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/17/2021] [Revised: 08/20/2021] [Accepted: 12/15/2021] [Indexed: 06/14/2023]
Abstract
Supersolid is a mysterious and puzzling state of matter whose possible existence has stirred a vigorous debate among physicists for over 60 years. Its elusive nature stems from the coexistence of two seemingly contradicting properties, long-range order and superfluidity. We report computational evidence of a supersolid phase of deuterium under high pressure (p>800 GPa) and low temperature (T<1.0 K). In our simulations, that are based on bosonic path integral molecular dynamics, we observe a highly concerted exchange of atoms while the system preserves its crystalline order. The exchange processes are favored by the soft core interactions between deuterium atoms that form a densely packed metallic solid. At the zero temperature limit, Bose-Einstein condensation is observed as the permutation probability of N deuterium atoms approaches 1/N with a finite superfluid fraction. Our study provides concrete evidence for the existence of a supersolid phase in high-pressure deuterium and could provide insights on the future investigation of supersolid phases in real materials.
Collapse
Affiliation(s)
- Chang Woo Myung
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lenseld Road, Cambridge, CB2 1EW, United Kingdom
| | - Barak Hirshberg
- School of Chemistry, Tel Aviv University, Tel Aviv 6997801, Israel
| | | |
Collapse
|
19
|
|
20
|
Abstract
The development of enhanced sampling methods has greatly extended the scope of atomistic simulations, allowing long-time phenomena to be studied with accessible computational resources. Many such methods rely on the identification of an appropriate set of collective variables. These are meant to describe the system's modes that most slowly approach equilibrium under the action of the sampling algorithm. Once identified, the equilibration of these modes is accelerated by the enhanced sampling method of choice. An attractive way of determining the collective variables is to relate them to the eigenfunctions and eigenvalues of the transfer operator. Unfortunately, this requires knowing the long-term dynamics of the system beforehand, which is generally not available. However, we have recently shown that it is indeed possible to determine efficient collective variables starting from biased simulations. In this paper, we bring the power of machine learning and the efficiency of the recently developed on the fly probability-enhanced sampling method to bear on this approach. The result is a powerful and robust algorithm that, given an initial enhanced sampling simulation performed with trial collective variables or generalized ensembles, extracts transfer operator eigenfunctions using a neural network ansatz and then accelerates them to promote sampling of rare events. To illustrate the generality of this approach, we apply it to several systems, ranging from the conformational transition of a small molecule to the folding of a miniprotein and the study of materials crystallization.
Collapse
Affiliation(s)
- Luigi Bonati
- Department of Physics, Eidgenössische Technische Hochschule (ETH) Zürich, 8092 Zürich, Switzerland;
- Atomistic Simulations, Italian Institute of Technology, 16163 Genova, Italy
| | | | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, 16163 Genova, Italy;
| |
Collapse
|
21
|
Abstract
We present an approach that extends the theory of targeted free energy perturbation (TFEP) to calculate free energy differences and free energy surfaces at an accurate quantum mechanical level of theory from a cheaper reference potential. The convergence is accelerated by a mapping function that increases the overlap between the target and the reference distributions. Building on recent work, we show that this map can be learned with a normalizing flow neural network, without requiring simulations with the expensive target potential but only a small number of single-point calculations, and, crucially, avoiding the systematic error that was found previously. We validate the method by numerically evaluating the free energy difference in a system with a double-well potential and by describing the free energy landscape of a simple chemical reaction in the gas phase.
Collapse
Affiliation(s)
- Andrea Rizzi
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
- Atomistic Simulations, Italian Institute of Technology, Via Morego 30, Genova 16163, Italy
| | - Paolo Carloni
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
- Molecular Neuroscience and Neuroimaging (INM-11), Forschungszentrum Jülich GmbH, Jülich 52428, Germany
- Department of Physics and Universitätsklinikum, RWTH Aachen University, Aachen 52074, Germany
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, Via Morego 30, Genova 16163, Italy
| |
Collapse
|
22
|
Abstract
The determination of efficient collective variables is crucial to the success of many enhanced sampling methods. As inspired by previous discrimination approaches, we first collect a set of data from the different metastable basins. The data are then projected with the help of a neural network into a low-dimensional manifold in which data from different basins are well-discriminated. This is here guaranteed by imposing that the projected data follows a preassigned distribution. The collective variables thus obtained lead to an efficient sampling and often allow reducing the number of collective variables in a multibasin scenario. We first check the validity of the method in two-state systems. We then move to multistep chemical processes. In the latter case, at variance with previous approaches, one single collective variable suffices, leading not only to computational efficiency but also to a very clear representation of the reaction free-energy profile.
Collapse
Affiliation(s)
- Enrico Trizio
- Atomistic Simulations, Istituto Italiano di Tecnologia, 16163 Genova, Italy
- Department of Materials Science, Università di Milano-Bicocca, 20126 Milano, Italy
| | - Michele Parrinello
- Atomistic Simulations, Istituto Italiano di Tecnologia, 16163 Genova, Italy
| |
Collapse
|
23
|
Ansari N, Rizzi V, Carloni P, Parrinello M. Water-Triggered, Irreversible Conformational Change of SARS-CoV-2 Main Protease on Passing from the Solid State to Aqueous Solution. J Am Chem Soc 2021; 143:12930-12934. [PMID: 34398611 PMCID: PMC8386029 DOI: 10.1021/jacs.1c05301] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2021] [Indexed: 11/30/2022]
Abstract
The main protease from SARS-CoV-2 is a homodimer. Yet, a recent 0.1-ms-long molecular dynamics simulation performed by D. E. Shaw's research group shows that it readily undergoes a symmetry-breaking event on passing from the solid state to aqueous solution. As a result, the subunits present distinct conformations of the binding pocket. By analyzing this long simulation, we uncover a previously unrecognized role of water molecules in triggering the transition. Interestingly, each subunit presents a different collection of long-lived water molecules. Enhanced sampling simulations performed here, along with machine learning approaches, further establish that the transition to the asymmetric state is essentially irreversible.
Collapse
Affiliation(s)
- Narjes Ansari
- Italian
Institute of Technology, Via E. Melen 83, 16152 Genova, Italy
| | - Valerio Rizzi
- Italian
Institute of Technology, Via E. Melen 83, 16152 Genova, Italy
| | - Paolo Carloni
- Computational
Biomedicine, Institute for Advanced Simulation (IAS-5) and Institute
of Neuroscience and Medicine (INM-9), and JARA-Institute “Molecular
Neuroscience and Neuroimaging” (INM-11), Forschungszentrum Jülich, Jülich 52425, Germany
- Physics
Department, RWTH Aachen University, Aachen 52074, Germany
| | | |
Collapse
|
24
|
Mandelli D, Parrinello M. A modified nudged elastic band algorithm with adaptive spring lengths. J Chem Phys 2021; 155:074103. [PMID: 34418926 DOI: 10.1063/5.0059593] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a modified version of the nudged elastic band (NEB) algorithm to find minimum energy paths connecting two known configurations. We show that replacing the harmonic band-energy term with a discretized version of the Onsager-Machlup action leads to a NEB algorithm with adaptive spring lengths that automatically increase the resolution of the minimum energy path around the saddle point of the potential energy surface. The method has the same computational cost per optimization step of the standard NEB algorithm and does not introduce additional parameters. We present applications to the isomerization of alanine dipeptide, the elimination of hydrogen from ethane, and the healing of a 5-77-5 defect in graphene.
Collapse
Affiliation(s)
- D Mandelli
- Atomistic Simulations, Italian Institute of Technology, Via Morego, 30 16163 Genova, Italy
| | - M Parrinello
- Atomistic Simulations, Italian Institute of Technology, Via Morego, 30 16163 Genova, Italy
| |
Collapse
|
25
|
Yang M, Karmakar T, Parrinello M. Liquid-Liquid Critical Point in Phosphorus. Phys Rev Lett 2021; 127:080603. [PMID: 34477397 DOI: 10.1103/physrevlett.127.080603] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/04/2021] [Revised: 07/07/2021] [Accepted: 07/19/2021] [Indexed: 06/13/2023]
Abstract
The study of liquid-liquid phase transitions has attracted considerable attention. One interesting example of this phenomenon is phosphorus, for which the existence of a first-order phase transition between a low density insulating molecular phase and a conducting polymeric phase has been experimentally established. In this Letter, we model this transition by an ab initio quality molecular dynamics simulation and explore a large portion of the liquid section of the phase diagram. We draw the liquid-liquid coexistence curve and determine that it terminates into a second-order critical point. Close to the critical point, large coupled structure and electronic structure fluctuations are observed.
Collapse
Affiliation(s)
- Manyi Yang
- Italian Institute of Technology, Via Melen 83, 16152 Genova, Italy
| | - Tarak Karmakar
- Italian Institute of Technology, Via Melen 83, 16152 Genova, Italy
| | | |
Collapse
|
26
|
Grifoni E, Piccini G, Lercher JA, Glezakou VA, Rousseau R, Parrinello M. Confinement effects and acid strength in zeolites. Nat Commun 2021; 12:2630. [PMID: 33976197 PMCID: PMC8113345 DOI: 10.1038/s41467-021-22936-0] [Citation(s) in RCA: 46] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2020] [Accepted: 03/30/2021] [Indexed: 02/03/2023] Open
Abstract
Chemical reactivity and sorption in zeolites are coupled to confinement and-to a lesser extent-to the acid strength of Brønsted acid sites (BAS). In presence of water the zeolite Brønsted acid sites eventually convert into hydronium ions. The gradual transition from zeolite Brønsted acid sites to hydronium ions in zeolites of varying pore size is examined by ab initio molecular dynamics combined with enhanced sampling based on Well-Tempered Metadynamics and a recently developed set of collective variables. While at low water content (1-2 water/BAS) the acidic protons prefer to be shared between zeolites and water, higher water contents (n > 2) invariably lead to solvation of the protons within a localized water cluster adjacent to the BAS. At low water loadings the standard free energy of the formed complexes is dominated by enthalpy and is associated with the acid strength of the BAS and the space around the site. Conversely, the entropy increases linearly with the concentration of waters in the pores, favors proton solvation and is independent of the pore size/shape.
Collapse
Affiliation(s)
- Emanuele Grifoni
- grid.5801.c0000 0001 2156 2780Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, Lugano, Ticino Switzerland ,grid.29078.340000 0001 2203 2861Institute of Computational Science, Università della Svizzera italiana (USI), Via Giuseppe Buffi 13, Lugano, Ticino Switzerland ,grid.6093.cPresent Address: Scuola Normale Superiore, Piazza dei Cavalieri, Pisa, Italy
| | - GiovanniMaria Piccini
- grid.5801.c0000 0001 2156 2780Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, Lugano, Ticino Switzerland ,grid.29078.340000 0001 2203 2861Institute of Computational Science, Università della Svizzera italiana (USI), Via Giuseppe Buffi 13, Lugano, Ticino Switzerland ,grid.451303.00000 0001 2218 3491Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland, WA USA
| | - Johannes A. Lercher
- grid.451303.00000 0001 2218 3491Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland, WA USA ,grid.6936.a0000000123222966Department Chemie and Catalysis Research Center, TU München, Lichtenbergstr. 4, Garching, Germany
| | - Vassiliki-Alexandra Glezakou
- grid.451303.00000 0001 2218 3491Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland, WA USA
| | - Roger Rousseau
- grid.451303.00000 0001 2218 3491Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland, WA USA
| | - Michele Parrinello
- grid.5801.c0000 0001 2156 2780Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, Lugano, Ticino Switzerland ,grid.29078.340000 0001 2203 2861Institute of Computational Science, Università della Svizzera italiana (USI), Via Giuseppe Buffi 13, Lugano, Ticino Switzerland ,grid.25786.3e0000 0004 1764 2907Italian Institute of Technology, Via Morego 30, Genova, Italy
| |
Collapse
|
27
|
Ahlawat P, Hinderhofer A, Alharbi EA, Lu H, Ummadisingu A, Niu H, Invernizzi M, Zakeeruddin SM, Dar MI, Schreiber F, Hagfeldt A, Grätzel M, Rothlisberger U, Parrinello M. A combined molecular dynamics and experimental study of two-step process enabling low-temperature formation of phase-pure α-FAPbI 3. Sci Adv 2021; 7:eabe3326. [PMID: 33893100 PMCID: PMC8064632 DOI: 10.1126/sciadv.abe3326] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/19/2020] [Accepted: 03/05/2021] [Indexed: 05/23/2023]
Abstract
It is well established that the lack of understanding the crystallization process in a two-step sequential deposition has a direct impact on efficiency, stability, and reproducibility of perovskite solar cells. Here, we try to understand the solid-solid phase transition occurring during the two-step sequential deposition of methylammonium lead iodide and formamidinium lead iodide. Using metadynamics, x-ray diffraction, and Raman spectroscopy, we reveal the microscopic details of this process. We find that the formation of perovskite proceeds through intermediate structures and report polymorphs found for methylammonium lead iodide and formamidinium lead iodide. From simulations, we discover a possible crystallization pathway for the highly efficient metastable α phase of formamidinium lead iodide. Guided by these simulations, we perform experiments that result in the low-temperature crystallization of phase-pure α-formamidinium lead iodide.
Collapse
Affiliation(s)
- Paramvir Ahlawat
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | | | - Essa A Alharbi
- Laboratory of Photonics and Interfaces, Institute of Chemical Sciences and Engineering, EPFL, CH-1015 Lausanne, Switzerland
| | - Haizhou Lu
- Laboratory of Photonics and Interfaces, Institute of Chemical Sciences and Engineering, EPFL, CH-1015 Lausanne, Switzerland
- Laboratory of Photomolecular Science, Institute of Chemical Sciences Engineering, EPFL, CH-1015 Lausanne, Switzerland
| | - Amita Ummadisingu
- Laboratory of Photonics and Interfaces, Institute of Chemical Sciences and Engineering, EPFL, CH-1015 Lausanne, Switzerland
| | - Haiyang Niu
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8092 Zürich, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera italiana, Via G. Buffi 13, 6900 Lugano, Switzerland
| | - Michele Invernizzi
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8092 Zürich, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera italiana, Via G. Buffi 13, 6900 Lugano, Switzerland
- Italian Institute of Technology, Via Morego 30, 16163 Genova, Italy
| | - Shaik Mohammed Zakeeruddin
- Laboratory of Photonics and Interfaces, Institute of Chemical Sciences and Engineering, EPFL, CH-1015 Lausanne, Switzerland
| | - M Ibrahim Dar
- Laboratory of Photonics and Interfaces, Institute of Chemical Sciences and Engineering, EPFL, CH-1015 Lausanne, Switzerland
- Cavendish Laboratory, Department of Physics, University of Cambridge, CB3 0HE, United Kingdom
| | - Frank Schreiber
- Institut für Angewandte Physik, Universität Tübingen, 72076 Tübingen, Germany.
| | - Anders Hagfeldt
- Laboratory of Photomolecular Science, Institute of Chemical Sciences Engineering, EPFL, CH-1015 Lausanne, Switzerland.
- Department of Chemistry, Ångström Laboratory, Uppsala University, Box 523, SE-751 20 Uppsala, Sweden
| | - Michael Grätzel
- Laboratory of Photonics and Interfaces, Institute of Chemical Sciences and Engineering, EPFL, CH-1015 Lausanne, Switzerland.
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland.
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8092 Zürich, Switzerland.
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera italiana, Via G. Buffi 13, 6900 Lugano, Switzerland
- Italian Institute of Technology, Via Morego 30, 16163 Genova, Italy
| |
Collapse
|
28
|
Affiliation(s)
- Tarak Karmakar
- Institute of Computational Sciences, Faculty of Informatics, Universit della Svizzera italiana, Lugano, Switzerland
- Department of Chemistry and Applied Biosciences, ETH Zurich, Zurich, Switzerland
- Italian Institute of Technology, Genova, Italy
| | - Michele Invernizzi
- Institute of Computational Sciences, Faculty of Informatics, Universit della Svizzera italiana, Lugano, Switzerland
- Italian Institute of Technology, Genova, Italy
- Department of Physics, ETH Zurich, Zurich, Switzerland
| | - Valerio Rizzi
- Institute of Computational Sciences, Faculty of Informatics, Universit della Svizzera italiana, Lugano, Switzerland
- Department of Chemistry and Applied Biosciences, ETH Zurich, Zurich, Switzerland
- Italian Institute of Technology, Genova, Italy
| | - Michele Parrinello
- Institute of Computational Sciences, Faculty of Informatics, Universit della Svizzera italiana, Lugano, Switzerland
- Department of Chemistry and Applied Biosciences, ETH Zurich, Zurich, Switzerland
- Italian Institute of Technology, Genova, Italy
| |
Collapse
|
29
|
Yang M, Bonati L, Polino D, Parrinello M. Using metadynamics to build neural network potentials for reactive events: the case of urea decomposition in water. Catal Today 2021. [DOI: 10.1016/j.cattod.2021.03.018] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
|
30
|
Invernizzi M, Parrinello M. Correction to "Rethinking Metadynamics: From Bias Potentials to Probability Distributions". J Phys Chem Lett 2021; 12:912. [PMID: 33439651 DOI: 10.1021/acs.jpclett.0c03541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
|
31
|
Abstract
One of the main applications of atomistic computer simulations is the calculation of ligand binding free energies. The accuracy of these calculations depends on the force field quality and on the thoroughness of configuration sampling. Sampling is an obstacle in simulations due to the frequent appearance of kinetic bottlenecks in the free energy landscape. Very often this difficulty is circumvented by enhanced sampling techniques. Typically, these techniques depend on the introduction of appropriate collective variables that are meant to capture the system's degrees of freedom. In ligand binding, water has long been known to play a key role, but its complex behaviour has proven difficult to fully capture. In this paper we combine machine learning with physical intuition to build a non-local and highly efficient water-describing collective variable. We use it to study a set of host-guest systems from the SAMPL5 challenge. We obtain highly accurate binding free energies and good agreement with experiments. The role of water during the binding process is then analysed in some detail.
Collapse
Affiliation(s)
- Valerio Rizzi
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092, Zurich, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, Via G. Buffi 13, 6900, Lugano, Switzerland
| | - Luigi Bonati
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, Via G. Buffi 13, 6900, Lugano, Switzerland
- Department of Physics, ETH Zurich, 8092, Zurich, Switzerland
| | - Narjes Ansari
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092, Zurich, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, Via G. Buffi 13, 6900, Lugano, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092, Zurich, Switzerland.
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, Via G. Buffi 13, 6900, Lugano, Switzerland.
- Italian Institute of Technology, Via Morego 30, 16163, Genova, Italy.
| |
Collapse
|
32
|
Abstract
We present an ab initio molecular dynamics (MD) investigation of the tautomeric equilibrium for the aqueous solutions of glycine and acetone under realistic experimental conditions. Metadynamics is used to accelerate proton migration among tautomeric centers. Due to the formation of complex water-ion structures involved in the proton dynamics in the aqueous environment, standard enhanced sampling approaches may face severe limitations in providing a general description of the phenomenon. Recently, we have developed a set of collective variables (CVs) designed to study protons transfer reactions in complex condensed systems [Grifoni, E. Proc. Natl. Acad. Sci. U.S.A. 2019, 116, 4054 4057]. In this work, we applied this approach to study proton dissociation dynamics leading to tautomeric interconversion of biologically and chemically relevant prototypical systems, namely, glycine and acetone in water. Although relatively simple from a chemical point of view, the results show that even for these small systems, complex reaction pathways and nontrivial conversion dynamics are observed. The generality of our method allows obtaining these results without providing any prior information on the dissociation dynamics but only the atomic species that can exchange protons in the process. Our results agree with literature estimates and demonstrate the general applicability of this method in the study of tautomeric reactions.
Collapse
Affiliation(s)
- Emanuele Grifoni
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland.,Institute of Computational Science, Università della Svizzera italiana (USI), Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland
| | - GiovanniMaria Piccini
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland.,Institute of Computational Science, Università della Svizzera italiana (USI), Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland.,Institute of Computational Science, Università della Svizzera italiana (USI), Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland.,Italian Institute of Technology, Via Morego 30, 16163 Genova, Italy
| |
Collapse
|
33
|
Piaggi PM, Valsson O, Parrinello M. Erratum: Enhancing Entropy and Enthalpy Fluctuations to Drive Crystallization in Atomistic Simulations [Phys. Rev. Lett. 119, 015701 (2017)]. Phys Rev Lett 2020; 125:159902. [PMID: 33095644 DOI: 10.1103/physrevlett.125.159902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/09/2020] [Indexed: 06/11/2023]
Abstract
This corrects the article DOI: 10.1103/PhysRevLett.119.015701.
Collapse
|
34
|
Ansari N, Karmakar T, Parrinello M. Molecular Mechanism of Gas Solubility in Liquid: Constant Chemical Potential Molecular Dynamics Simulations. J Chem Theory Comput 2020; 16:5279-5286. [PMID: 32551636 DOI: 10.1021/acs.jctc.0c00450] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Accurate prediction of gas solubility in a liquid is crucial in many areas of chemistry, and a detailed understanding of the molecular mechanism of the gas solvation continues to be an active area of research. Here, we extend the idea of the constant chemical potential molecular dynamics (CμMD) approach to the calculation of the gas solubility in the liquid under constant gas chemical potential conditions. As a representative example, we utilize this method to calculate the isothermal solubility of carbon dioxide in water. Additionally, we provide microscopic insight into the mechanism of solvation that preferentially occurs in areas of the surface where the hydrogen network is broken.
Collapse
Affiliation(s)
- Narjes Ansari
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland.,Facoltà di informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, CH-6900 Lugano, Switzerland
| | - Tarak Karmakar
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland.,Facoltà di informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, CH-6900 Lugano, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland.,Facoltà di informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, CH-6900 Lugano, Switzerland.,Italian Institute of Technology, Via Morego 30, 16163 Genova, Italy
| |
Collapse
|
35
|
Capelli R, Lyu W, Bolnykh V, Meloni S, Olsen JMH, Rothlisberger U, Parrinello M, Carloni P. Accuracy of Molecular Simulation-Based Predictions of koff Values: A Metadynamics Study. J Phys Chem Lett 2020; 11:6373-6381. [PMID: 32672983 DOI: 10.1021/acs.jpclett.0c00999] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
The koff values of ligands unbinding to proteins are key parameters for drug discovery. Their predictions based on molecular simulation may under- or overestimate experiment in a system- and/or technique-dependent way. Here we use an established method-infrequent metadynamics, based on the AMBER force field-to compute the koff of the ligand iperoxo (in clinical use) targeting the muscarinic receptor M2. The ligand charges are calculated by either (i) the Amber standard procedure or (ii) B3LYP-DFT. The calculations using (i) turn out not to provide a reasonable estimation of the transition-state free energy. Those using (ii) differ from experiment by 2 orders of magnitude. On the basis of B3LYP DFT QM/MM simulations, we suggest that the observed discrepancy in (ii) arises, at least in part, from the lack of electronic polarization and/or charge transfer in biomolecular force fields. These issues might be present in other systems, such as DNA-protein complexes.
Collapse
Affiliation(s)
- Riccardo Capelli
- Computational Biomedicine Section, IAS-5/INM-9, Forschungzentrum Jülich, Wilhelm-Johnen-straße, D-52425 Jülich, Germany
- JARA-HPC, Forschungszentrum Jülich, D-54245 Jülich, Germany
| | - Wenping Lyu
- Computational Biomedicine Section, IAS-5/INM-9, Forschungzentrum Jülich, Wilhelm-Johnen-straße, D-52425 Jülich, Germany
| | - Viacheslav Bolnykh
- Computational Biomedicine Section, IAS-5/INM-9, Forschungzentrum Jülich, Wilhelm-Johnen-straße, D-52425 Jülich, Germany
- Laboratory of Computational Chemistry and Biochemistry, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Simone Meloni
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università degli Studi di Ferrara, Via Luigi Borsari 46, I-44121, Ferrara, Italy
| | - Jógvan Magnus Haugaard Olsen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zürich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana (USI), Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
- Istituto Italiano di Tecnologia, Via Morego 30, I-16163 Genova, Italy
| | - Paolo Carloni
- Computational Biomedicine Section, IAS-5/INM-9, Forschungzentrum Jülich, Wilhelm-Johnen-straße, D-52425 Jülich, Germany
- JARA-Institute INM-11: Molecular Neuroscience and Neuroimaging, Forschungzentrum Jülich, Wilhelm-Johnen-straße, D-52425 Jülich, Germany
| |
Collapse
|
36
|
Mandelli D, Hirshberg B, Parrinello M. Metadynamics of Paths. Phys Rev Lett 2020; 125:026001. [PMID: 32701329 DOI: 10.1103/physrevlett.125.026001] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/17/2020] [Revised: 05/09/2020] [Accepted: 06/12/2020] [Indexed: 05/27/2023]
Abstract
We present a method to sample reactive pathways via biased molecular dynamics simulations in trajectory space. We show that the use of enhanced sampling techniques enables unconstrained exploration of multiple reaction routes. Time correlation functions are conveniently computed via reweighted averages along a single trajectory and kinetic rates are accessed at no additional cost. These abilities are illustrated analyzing a model potential and the umbrella inversion of NH_{3} in water. The algorithm allows a parallel implementation and promises to be a powerful tool for the study of rare events.
Collapse
Affiliation(s)
- Davide Mandelli
- Atomistic Simulations, Italian Institute of Technology, via Morego 30, 16163 Genova, Italy
| | - Barak Hirshberg
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland
- Institute of Computational Sciences, Università della Svizzera Italiana, 6900 Lugano, Switzerland
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, via Morego 30, 16163 Genova, Italy
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland
- Institute of Computational Sciences, Università della Svizzera Italiana, 6900 Lugano, Switzerland
| |
Collapse
|
37
|
Abstract
We introduce an enhanced sampling method that is based on constructing a model probability density from which a bias potential is derived. The model relies on the fact that in a physical system most of the configurations visited can be grouped into isolated metastable islands. With each island we associate a distribution that is fitted to a Gaussian mixture. The different distributions are linearly combined together with coefficients that are computed self-consistently. This leads to an integrated procedure for discovering new metastable states, exploring reaction pathways, computing free energy differences, and estimating reaction rates.
Collapse
Affiliation(s)
- Jayashrita Debnath
- Department of Chemistry and Applied Biosciences, ETH Zürich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Universitá della Svizzera Italiana (USI), Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zürich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Universitá della Svizzera Italiana (USI), Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
- Instituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| |
Collapse
|
38
|
Bottaro S, Nichols PJ, Vögeli B, Parrinello M, Lindorff-Larsen K. Integrating NMR and simulations reveals motions in the UUCG tetraloop. Nucleic Acids Res 2020; 48:5839-5848. [PMID: 32427326 PMCID: PMC7293013 DOI: 10.1093/nar/gkaa399] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2020] [Revised: 04/03/2020] [Accepted: 05/17/2020] [Indexed: 12/21/2022] Open
Abstract
We provide an atomic-level description of the structure and dynamics of the UUCG RNA stem-loop by combining molecular dynamics simulations with experimental data. The integration of simulations with exact nuclear Overhauser enhancements data allowed us to characterize two distinct states of this molecule. The most stable conformation corresponds to the consensus three-dimensional structure. The second state is characterized by the absence of the peculiar non-Watson-Crick interactions in the loop region. By using machine learning techniques we identify a set of experimental measurements that are most sensitive to the presence of non-native states. We find that although our MD ensemble, as well as the consensus UUCG tetraloop structures, are in good agreement with experiments, there are remaining discrepancies. Together, our results show that (i) the MD simulation overstabilize a non-native loop conformation, (ii) eNOE data support its presence with a population of ≈10% and (iii) the structural interpretation of experimental data for dynamic RNAs is highly complex, even for a simple model system such as the UUCG tetraloop.
Collapse
Affiliation(s)
- Sandro Bottaro
- Atomistic Simulations Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
- Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Parker J Nichols
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver, Anschutz Medical Campus, Aurora, CO, USA
| | - Beat Vögeli
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver, Anschutz Medical Campus, Aurora, CO, USA
| | - Michele Parrinello
- Atomistic Simulations Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
39
|
Hirshberg B, Invernizzi M, Parrinello M. Path integral molecular dynamics for fermions: Alleviating the sign problem with the Bogoliubov inequality. J Chem Phys 2020; 152:171102. [DOI: 10.1063/5.0008720] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Barak Hirshberg
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland
- Institute of Computational Sciences, Università della Svizzera italiana, via G. Buffi 13, 6900 Lugano, Switzerland
| | - Michele Invernizzi
- Institute of Computational Sciences, Università della Svizzera italiana, via G. Buffi 13, 6900 Lugano, Switzerland
- National Centre for Computational Design and Discovery of Novel Materials MARVEL, Università della Svizzera italiana, via G. Buffi 13, 6900 Lugano, Switzerland
- Department of Physics, ETH Zurich, 8092 Zurich, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland
- Institute of Computational Sciences, Università della Svizzera italiana, via G. Buffi 13, 6900 Lugano, Switzerland
- Atomistic Simulations, Italian Institute of Technology, Via Morego 30, 16163 Genova, Italy
| |
Collapse
|
40
|
Polino D, Grifoni E, Rousseau R, Parrinello M, Glezakou VA. How Collective Phenomena Impact CO2 Reactivity and Speciation in Different Media. J Phys Chem A 2020; 124:3963-3975. [DOI: 10.1021/acs.jpca.9b11744] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Daniela Polino
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana,Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| | - Emanuele Grifoni
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana,Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| | - Roger Rousseau
- Pacific Northwest National Laboratory, 902 Battelle Blvd., P.O. Box 999, MSIN K1-83, Richland, Washington 99352, United States
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana,Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
- Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy
| | - Vassiliki-Alexandra Glezakou
- Pacific Northwest National Laboratory, 902 Battelle Blvd., P.O. Box 999, MSIN K1-83, Richland, Washington 99352, United States
| |
Collapse
|
41
|
Abstract
Designing an appropriate set of collective variables is crucial to the success of several enhanced sampling methods. Here we focus on how to obtain such variables from information limited to the metastable states. We characterize these states by a large set of descriptors and employ neural networks to compress this information in a lower-dimensional space, using Fisher's linear discriminant as an objective function to maximize the discriminative power of the network. We test this method on alanine dipeptide, using the nonlinearly separable data set composed by atomic distances. We then study an intermolecular aldol reaction characterized by a concerted mechanism. The resulting variables are able to promote sampling by drawing nonlinear paths in the physical space connecting the fluctuations between metastable basins. Lastly, we interpret the behavior of the neural network by studying its relation to the physical variables. Through the identification of its most relevant features, we are able to gain chemical insight into the process.
Collapse
Affiliation(s)
- Luigi Bonati
- Department of Physics, ETH Zurich, 8092 Zurich, Switzerland
- Institute of Computational Sciences, Università della Svizzera italiana, via Buffi 13, 6900 Lugano, Switzerland
| | - Valerio Rizzi
- Institute of Computational Sciences, Università della Svizzera italiana, via Buffi 13, 6900 Lugano, Switzerland
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland
| | - Michele Parrinello
- Institute of Computational Sciences, Università della Svizzera italiana, via Buffi 13, 6900 Lugano, Switzerland
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland
- Italian Institute of Technology, Via Morego 30, 16163 Genova, Italy
| |
Collapse
|
42
|
Abstract
Metadynamics is an enhanced sampling method of great popularity, based on the on-the-fly construction of a bias potential that is a function of a selected number of collective variables. We propose here a change in perspective that shifts the focus from the bias to the probability distribution reconstruction while retaining some of the key characteristics of metadynamics, such as flexible on-the-fly adjustments to the free energy estimate. The result is an enhanced sampling method that presents a drastic improvement in convergence speed, especially when dealing with suboptimal and/or multidimensional sets of collective variables. The method is especially robust and easy to use and in fact requires only a few simple parameters to be set, and it has a straightforward reweighting scheme to recover the statistics of the unbiased ensemble. Furthermore, it gives more control of the desired exploration of the phase space since the deposited bias is not allowed to grow indefinitely and it does not push the simulation to uninteresting high free energy regions. We demonstrate the performance of the method in a number of representative examples.
Collapse
Affiliation(s)
- Michele Invernizzi
- Department of Physics, ETH Zurich, c/o Università della Svizzera Italiana, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
- Facoltà di Informatica, Institute of Computational Science, National Center for Computational Design and Discovery of Novel Materials (MARVEL), Università della Svizzera Italiana, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o Università della Svizzera italiana, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
- Italian Institute of Technology, Via Morego 30, 16163 Genova, Italy
- Facoltà di Informatica, Institute of Computational Science, Università della Svizzera Italiana, 6900 Lugano, Switzerland
| |
Collapse
|
43
|
Ponte F, Piccini G, Sicilia E, Parrinello M. A metadynamics perspective on the reduction mechanism of the Pt(IV) asplatin prodrug. J Comput Chem 2019; 41:290-294. [PMID: 31691997 DOI: 10.1002/jcc.26100] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Revised: 10/16/2019] [Accepted: 10/17/2019] [Indexed: 11/07/2022]
Abstract
Enhanced sampling molecular dynamics has been used to model the reduction mechanism of the antitumoral Asplatin Pt(IV) complex, c,c,t-[PtCl2(NH3)2(OH)(aspirin)] in the presence of l-ascorbic acid as reducing agent. In order to overcome the timescale problem, characteristic of many chemical reactions, we enhanced the sampling of the free energy landscape using Metadynamics. To achieve such a goal, the selection of adequate collective variables is crucial for the application of the method. Recently, a new method called Multi-Class Harmonic Linear Discriminant Analysis (MC-HLDA) has been proposed as a tool for constructing collective variables (CVs) for complex chemical processes. The method reduces the dimensionality of the variable space by generating appropriate linear combinations of several relevant chemical descriptors. The aim of this work is to assess the ability and performance of this method in describing the fundamental features of complex chemical reactions such as the Asplatin reduction mechanism in a compact, simple, and physically transparent manner. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Fortuna Ponte
- Dipartimento di Chimica e Tecnologie Chimiche, Università della Calabria, Ponte P. Bucci, Cubo 14C, Arcavacata di Rende, 87030, Italy
| | - GiovanniMaria Piccini
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, Lugano, 6900, Switzerland.,Facoltà di Informatica, Istituto di Scienze Computazionali, Università della SvizzeraItaliana (USI), Via Giuseppe Buffi 13, Lugano, 6900, Switzerland
| | - Emilia Sicilia
- Dipartimento di Chimica e Tecnologie Chimiche, Università della Calabria, Ponte P. Bucci, Cubo 14C, Arcavacata di Rende, 87030, Italy
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, Lugano, 6900, Switzerland.,Facoltà di Informatica, Istituto di Scienze Computazionali, Università della SvizzeraItaliana (USI), Via Giuseppe Buffi 13, Lugano, 6900, Switzerland.,Istituto Italiano di Tecnologia, Via Morego 30, Genova, 16163, Italy
| |
Collapse
|
44
|
Karmakar T, Piaggi PM, Parrinello M. Molecular Dynamics Simulations of Crystal Nucleation from Solution at Constant Chemical Potential. J Chem Theory Comput 2019; 15:6923-6930. [DOI: 10.1021/acs.jctc.9b00795] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Tarak Karmakar
- Department of Chemistry and Applied Biosciences, ETH Zürich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computationali, Università della Svizzera Italiana, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
| | - Pablo M. Piaggi
- Department of Chemistry and Applied Biosciences, ETH Zürich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computationali, Università della Svizzera Italiana, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zürich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computationali, Università della Svizzera Italiana, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
| |
Collapse
|
45
|
Polino D, Parrinello M. Kinetics of Aqueous Media Reactions via Ab Initio Enhanced Molecular Dynamics: The Case of Urea Decomposition. J Phys Chem B 2019; 123:6851-6856. [PMID: 31286763 DOI: 10.1021/acs.jpcb.9b05271] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Aqueous solutions provide a medium for many important reactions in chemical synthesis, industrial processes, environmental chemistry, and biological functions. It is an accepted fact that aqueous solvents can be direct participants in the reaction process and not act only as simple passive dielectrics. Assisting water molecules and proton wires are thus essential for the efficiency of many reactions. Here, we study the decomposition of urea into ammonia and isocyanic acid by means of enhanced ab initio molecular dynamics simulations. We highlight the role of the solvent molecules and their interactions with the reactants providing a proper description of the reaction mechanism and how the water hydrogen-bond network affects the reaction dynamics. Reaction free energy and rates have been calculated taking into account this important effect.
Collapse
Affiliation(s)
- Daniela Polino
- Department of Chemistry and Applied Biosciences , ETH Zurich , c/o USI Campus, Via Giuseppe Buffi 13 , CH-6900 Lugano , Switzerland.,Facoltà di Informatica, Istituto di Scienze Computazionali , Università della Svizzera Italiana , Via Giuseppe Buffi 13 , CH-6900 Lugano , Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences , ETH Zurich , c/o USI Campus, Via Giuseppe Buffi 13 , CH-6900 Lugano , Switzerland.,Facoltà di Informatica, Istituto di Scienze Computazionali , Università della Svizzera Italiana , Via Giuseppe Buffi 13 , CH-6900 Lugano , Switzerland.,Istituto Italiano di Tecnologia , Via Morego 30 , 16163 Genova , Italy
| |
Collapse
|
46
|
Rizzi V, Mendels D, Sicilia E, Parrinello M. Blind Search for Complex Chemical Pathways Using Harmonic Linear Discriminant Analysis. J Chem Theory Comput 2019; 15:4507-4515. [DOI: 10.1021/acs.jctc.9b00358] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Valerio Rizzi
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera italiana (USI), Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
| | - Dan Mendels
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera italiana (USI), Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Emilia Sicilia
- Dipartimento di Chimica e Tecnologie Chimiche, Università della Calabria, 87036 Rende CS, Italy
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera italiana (USI), Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
- Istituto Italiano di Tecnologia (IIT), Via Morego, 30, 16163 Genova GE, Italy
| |
Collapse
|
47
|
Abstract
Free energy sampling methods allow studying the full dynamics of activated processes. Unfortunately, the affordable accuracy of the potential describing the energy and forces of the system is usually rather low. Here we introduce a new method that by combining metadynamics and free energy perturbation allows calculating accurate quantum chemical free energies for chemical reactions. To prove the effectiveness of this new approach we study the SN2 reaction of CH3F + Cl- → CH3Cl + F- in vacuo and solvated by water. Comparisons are made with harmonic transition-state theory to show how this method could provide accurate equilibrium and rate constants for complex systems.
Collapse
Affiliation(s)
- GiovanniMaria Piccini
- Department of Chemistry and Applied Biosciences , ETH Zurich , c/o USI Campus, Via Giuseppe Buffi 13 , CH-6900 Lugano , Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali , Università della SvizzeraItaliana (USI) , Via Giuseppe Buffi 13 , CH-6900 Lugano , Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences , ETH Zurich , c/o USI Campus, Via Giuseppe Buffi 13 , CH-6900 Lugano , Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali , Università della SvizzeraItaliana (USI) , Via Giuseppe Buffi 13 , CH-6900 Lugano , Switzerland
- Istituto Italiano di Tecnologia , Via Morego 30 , 16163 Genova , Italy
| |
Collapse
|
48
|
Abstract
From the Ising model and the Lennard-Jones fluid to water and the iron-carbon system, phase diagrams are an indispensable tool to understand phase equilibria. Despite the effort of the simulation community, the calculation of a large portion of a phase diagram using computer simulation is still today a significant challenge. Here, we propose a method to calculate phase diagrams involving liquid and solid phases by the reversible transformation of the liquid and the solid. To this end, we introduce an order parameter that breaks the rotational symmetry and we leverage our recently introduced method to sample the multithermal-multibaric ensemble. In this way, in a single molecular dynamics simulation, we are able to compute the liquid-solid coexistence line for entire regions of the temperature and pressure phase diagram. We apply our approach to the bcc-liquid phase diagram of sodium and the fcc-bcc-liquid phase diagram of aluminum.
Collapse
Affiliation(s)
- Pablo M Piaggi
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
| |
Collapse
|
49
|
Niu H, Yang YI, Parrinello M. Temperature Dependence of Homogeneous Nucleation in Ice. Phys Rev Lett 2019; 122:245501. [PMID: 31322390 DOI: 10.1103/physrevlett.122.245501] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/30/2019] [Revised: 05/06/2019] [Indexed: 06/10/2023]
Abstract
Ice nucleation is a process of great relevance in physics, chemistry, technology, and environmental sciences; much theoretical effort has been devoted to its understanding, but it still remains a topic of intense research. We shed light on this phenomenon by performing atomistic based simulations. Using metadynamics and a carefully designed set of collective variables, reversible transitions between water and ice are able to be simulated. We find that water freezes into a stacking disordered structure with the all-atom transferable intermolecular potential with 4 points/ice (TIP4P/ice) model, and the features of the critical nucleus of nucleation at the microscopic level are revealed. We have also estimated the ice nucleation rates along with other nucleation parameters at different undercoolings. Our results are in agreement with recent experimental and other theoretical works, and they confirm that nucleation is preceded by a large increase in tetrahedrally coordinated water molecules.
Collapse
Affiliation(s)
- Haiyang Niu
- Department of Chemistry and Applied Biosciences, ETH Zurich c/o USI Campus, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
- Facoltà di Informatica, Instituto di Scienze Computationali, and National Center for Computational Design and Discovery of Novel Materials MARVEL, Università della Svizzera Italiana, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
| | - Yi Isaac Yang
- Department of Chemistry and Applied Biosciences, ETH Zurich c/o USI Campus, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
- Facoltà di Informatica, Instituto di Scienze Computationali, and National Center for Computational Design and Discovery of Novel Materials MARVEL, Università della Svizzera Italiana, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich c/o USI Campus, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
- Facoltà di Informatica, Instituto di Scienze Computationali, and National Center for Computational Design and Discovery of Novel Materials MARVEL, Università della Svizzera Italiana, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
- Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy
| |
Collapse
|
50
|
Abstract
Determining the complete set of ligands' binding-unbinding pathways is important for drug discovery and for rational interpretation of mutation data. Here we have developed a metadynamics-based technique that addresses this issue and allows estimating affinities in the presence of multiple escape pathways. Our approach is shown on a lysozyme T4 variant in complex with a benzene molecule. The calculated binding free energy is in agreement with experimental data. Remarkably, not only were we able to find all the previously identified ligand binding pathways, but also we identified three pathways previously not identified as such. These results were obtained at a small computational cost, making this approach valuable for practical applications, such as screening of small compound libraries.
Collapse
Affiliation(s)
- Riccardo Capelli
- Computational Biomedicine (INM-9/IAS-5) , Forschungszentrum Jülich , Wilhelm-Johnen-Straße , D-52425 Jülich , Germany
- JARA-HPC, Forschungszentrum Jülich , D-54245 Jülich , Germany
| | - Paolo Carloni
- Computational Biomedicine (INM-9/IAS-5) , Forschungszentrum Jülich , Wilhelm-Johnen-Straße , D-52425 Jülich , Germany
- Department of Physics , RWTH Aachen University , D-52078 Aachen , Germany
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences , ETH Zürich , c/o USI Campus, Via Giuseppe Buffi 13 , CH-6900 Lugano , Ticino Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali , Università della Svizzera Italiana (USI) , Via Giuseppe Buffi 13 , CH-6900 Lugano , Ticino Switzerland
- Istituto Italiano di Tecnologia , Via Morego 30 , I-16163 Genova , Italy
| |
Collapse
|