Showing 1–10 of 10 results for author: Rizzi, V

DeepLNE++ leveraging knowledge distillation for accelerated multi-state path-like collective variables

Authors: Thorben Fröhlking, Valerio Rizzi, Simone Aureli, Francesco Luigi Gervasio

Abstract: Path-like collective variables can be very effective for accurately modeling complex biomolecular processes in molecular dynamics simulations. Recently, we introduced DeepLNE, a machine learning-based path-like CV that provides a progression variable s along the path as a non-linear combination of several descriptors, effectively approximating the reaction coordinate. However, DeepLNE is computati… ▽ More Path-like collective variables can be very effective for accurately modeling complex biomolecular processes in molecular dynamics simulations. Recently, we introduced DeepLNE, a machine learning-based path-like CV that provides a progression variable s along the path as a non-linear combination of several descriptors, effectively approximating the reaction coordinate. However, DeepLNE is computationally expensive for realistic systems needing many descriptors and limited in its ability to handle multi-state reactions. Here we present DeepLNE++, which uses a knowledge distillation approach to significantly accelerate the evaluation of DeepLNE, making it feasible to compute free energy landscapes for large and complex biomolecular systems. In addition, DeepLNE++ encodes system-specific knowledge within a supervised multitasking framework, enhancing its versatility and effectiveness. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Deep learning path-like collective variable for enhanced sampling molecular dynamics

Authors: Thorben Fröhlking, Luigi Bonati, Valerio Rizzi, Francesco Luigi Gervasio

Abstract: Several enhanced sampling techniques rely on the definition of collective variables to effectively explore free energy landscapes. Existing variables that describe the progression along a reactive pathway offer an elegant solution but face a number of limitations. In this paper, we address these challenges by introducing a new path-like collective variable called the `Deep-locally-non-linear-embed… ▽ More Several enhanced sampling techniques rely on the definition of collective variables to effectively explore free energy landscapes. Existing variables that describe the progression along a reactive pathway offer an elegant solution but face a number of limitations. In this paper, we address these challenges by introducing a new path-like collective variable called the `Deep-locally-non-linear-embedding', which is inspired by principles of the locally linear embedding technique and is trained on a reactive trajectory. The variable mimics the ideal reaction coordinate by automatically generating a non-linear combination of features through a differentiable generalized autoencoder that combines a neural network with a continuous k-nearest-neighbor selection. Among the key advantages of this method is its capability to automatically choose the metric for searching neighbors and to learn the path from state A to state B without the need to handpick landmarks a priori. We demonstrate the effectiveness of DeepLNE by showing that the progression along the path variable closely approximates the ideal reaction coordinate in toy models such as the Müller-Brown-potential and alanine dipeptide. We then use it in molecular dynamics simulations of an RNA tetraloop, where we highlight its capability to accelerate transitions and converge the free energy of folding. △ Less

Submitted 2 February, 2024; originally announced February 2024.

Rare Event Kinetics from Adaptive Bias Enhanced Sampling

Authors: Dhiman Ray, Narjes Ansari, Valerio Rizzi, Michele Invernizzi, Michele Parrinello

Abstract: We introduce a novel enhanced sampling approach named OPES flooding for calculating the kinetics of rare events from atomistic molecular dynamics simulation. This method is derived from the On-the-fly-Probability-Enhanced-Sampling (OPES) approach [Invernizzi and Parrinello, JPC Lett. 2020], which has been recently developed for calculating converged free energy surfaces for complex systems. In thi… ▽ More We introduce a novel enhanced sampling approach named OPES flooding for calculating the kinetics of rare events from atomistic molecular dynamics simulation. This method is derived from the On-the-fly-Probability-Enhanced-Sampling (OPES) approach [Invernizzi and Parrinello, JPC Lett. 2020], 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 gas phase, and the folding and unfolding of the chignolin polypeptide in 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. △ Less

Submitted 16 September, 2022; v1 submitted 8 August, 2022; originally announced August 2022.

Water regulates the residence time of Benzamidine in Trypsin

Authors: Narjes Ansari, Valerio Rizzi, Michele Parrinello

Abstract: The process of ligand-protein unbinding is crucial in biophysics. Water is an essential part of any biological system and yet, many aspects of its role remain elusive. Here, we simulate with state-of-the-art enhanced sampling techniques the binding of Benzamidine to Trypsin which is a much studied and paradigmatic ligand-protein system. We use machine learning methods to determine efficient collec… ▽ More The process of ligand-protein unbinding is crucial in biophysics. Water is an essential part of any biological system and yet, many aspects of its role remain elusive. Here, we simulate with state-of-the-art enhanced sampling techniques the binding of Benzamidine to Trypsin which is a much studied and paradigmatic ligand-protein system. We use machine learning methods to determine efficient collective coordinates for the complex non-local network of water. These coordinates are used to perform On-the-fly Probability Enhanced Sampling simulations, which we adapt to calculate also the ligand residence time. Our results, both static and dynamic, are in good agreement with experiments. We find that the presence of a water molecule located at the bottom of the binding pocket allows via a network of hydrogen bonds the ligand to be released into the solution. On a finer scale, even when unbinding is allowed, another water molecule further modulates the exit time. △ Less

Submitted 28 September, 2022; v1 submitted 12 April, 2022; originally announced April 2022.

Journal ref: Nat. Commun. 13, 5438 (2022)

The role of water in host-guest interaction

Authors: Valerio Rizzi, Luigi Bonati, Narjes Ansari, Michele Parrinello

Abstract: One of the main applications of atomistic computer simulations is the calculation of ligand binding energies. The accuracy of these calculations depends on the force field quality and on the thoroughness of configuration sampling. Sampling is an obstacle in modern simulations due to the frequent appearance of kinetic bottlenecks in the free energy landscape. Very often this difficulty is circumven… ▽ More One of the main applications of atomistic computer simulations is the calculation of ligand binding energies. The accuracy of these calculations depends on the force field quality and on the thoroughness of configuration sampling. Sampling is an obstacle in modern 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 of host-guest systems from the SAMPL5 challenge. We obtain highly accurate binding energies and good agreement with experiments. The role of water during the binding process is then analysed in some detail. △ Less

Submitted 8 January, 2021; v1 submitted 23 June, 2020; originally announced June 2020.

Journal ref: Nature Communications 12, 93 (2021)

Data-Driven Collective Variables for Enhanced Sampling

Authors: Luigi Bonati, Valerio Rizzi, Michele Parrinello

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… ▽ More 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 non-linearly separable dataset 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 non-linear 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. △ Less

Submitted 7 April, 2020; v1 submitted 16 February, 2020; originally announced February 2020.

Journal ref: J. Phys. Chem. Lett. 2020, 11, 2998-3004

Path Integral Molecular Dynamics for Bosons

Authors: Barak Hirshberg, Valerio Rizzi, Michele Parrinello

Abstract: Trapped Bosons exhibit fundamental physical phenomena and are potentially useful for quantum technologies. We present a method for simulating Bosons using path integral molecular dynamics. A main challenge for simulations is including all permutations due to exchange symmetry. We show that evaluation of the potential can be done recursively, avoiding explicit enumeration of permutations, and scale… ▽ More Trapped Bosons exhibit fundamental physical phenomena and are potentially useful for quantum technologies. We present a method for simulating Bosons using path integral molecular dynamics. A main challenge for simulations is including all permutations due to exchange symmetry. We show that evaluation of the potential can be done recursively, avoiding explicit enumeration of permutations, and scales cubically with system size. The method is applied to Bosons in a 2D trap and agrees with essentially exact results. An analysis of the role of exchange with decreasing temperature is also presented. △ Less

Submitted 22 May, 2019; originally announced May 2019.

Journal ref: Proc. Natl. Acad. Sci. U.S.A. 116 (2019) 21445-21449

Blind search for complex chemical pathways using Harmonic Linear Discriminant Analysis

Authors: Valerio Rizzi, Dan Mendels, Emilia Sicilia, Michele Parrinello

Abstract: Disentangling the mechanistic details of a chemical reaction pathway is a hard problem that often requires a considerable amount of chemical intuition and a component of luck. Experiments struggle in observing short-life metastable intermediates, while computer simulations often rely upon a good initial guess. In this work, we propose a method that, from the simulations of a reactant and a product… ▽ More Disentangling the mechanistic details of a chemical reaction pathway is a hard problem that often requires a considerable amount of chemical intuition and a component of luck. Experiments struggle in observing short-life metastable intermediates, while computer simulations often rely upon a good initial guess. In this work, we propose a method that, from the simulations of a reactant and a product state, searches for reaction mechanisms connecting the two by exploring the configuration space through metadynamics, a well known enhanced molecular dynamics method. The key quantity underlying this search is based on the use of an approach called Harmonic Linear Discriminant Analysis which allows a systematic construction of collective variables. Given the reactant and product states, we choose a set of descriptors capable of discriminating between the two states. In order not to prejudge the results, generic descriptors are introduced. The fluctuations of the descriptors in the two states are used to construct collective variables. We use metadynamics in an exploratory mode to discover the intermediates and the transition states that lead from reactant to product. The search is at first conducted at a low theory level. The calculation is then refined and the energy of the intermediates and transition states discovered during metadynamics is computed again using a higher level of theory. The method's aim is to offer a simple reaction mechanism search procedure that helps in saving time and is able to find unexpected mechanisms that defy well established chemical paradigms. We apply it to two reactions, showing that a high level of complexity can be hidden even in seemingly trivial and small systems. The method can be applied to larger systems, such as reactions in solution or catalysis. △ Less

Submitted 31 July, 2019; v1 submitted 12 April, 2019; originally announced April 2019.

The onset of dehydrogenation in solid Ammonia Borane, an ab-initio metadynamics study

Authors: Valerio Rizzi, Daniela Polino, Emilia Sicilia, Nino Russo, Michele Parrinello

Abstract: The discovery of effective hydrogen storage materials is fundamental for the progress of a clean energy economy. Ammonia borane ($\mathrm{H_3BNH_3}$) has attracted great interest as a promising candidate but the reaction path that leads from its solid phase to hydrogen release is not yet fully understood. To address the need for insights in the atomistic details of such a complex solid state proce… ▽ More The discovery of effective hydrogen storage materials is fundamental for the progress of a clean energy economy. Ammonia borane ($\mathrm{H_3BNH_3}$) has attracted great interest as a promising candidate but the reaction path that leads from its solid phase to hydrogen release is not yet fully understood. To address the need for insights in the atomistic details of such a complex solid state process, in this work we use \textit{ab-initio} molecular dynamics and metadynamics to study the early stages of AB dehydrogenation. We show that the formation of ammonia diborane ($\mathrm{H_3NBH_2(}$$μ$$\mathrm{-H)BH_3}$) leads to the release of $\mathrm{NH_4^+}$, which in turn triggers an autocatalytic $\mathrm{H_2}$ production cycle. Our calculations provide a model for how complex solid state reactions can be theoretically investigated and rely upon the presence of multiple ammonia borane molecules, as substantiated by standard quantum-mechanical simulations on a cluster. △ Less

Submitted 29 January, 2019; v1 submitted 17 September, 2018; originally announced September 2018.

Inelastic electron injection in a water chain

Authors: Valerio Rizzi, Tchavdar N. Todorov, Jorge J. Kohanoff

Abstract: Irradiation of biological matter triggers a cascade of secondary particles that interact with their surroundings, resulting in damage. Low-energy electrons are one of the main secondary species and electron-phonon interaction plays a fundamental role in their dynamics. We have developed a method to capture the electron-phonon inelastic energy exchange in real time and have used it to inject electr… ▽ More Irradiation of biological matter triggers a cascade of secondary particles that interact with their surroundings, resulting in damage. Low-energy electrons are one of the main secondary species and electron-phonon interaction plays a fundamental role in their dynamics. We have developed a method to capture the electron-phonon inelastic energy exchange in real time and have used it to inject electrons into a simple system that models a biological environment, a water chain. We simulated both an incoming electron pulse and a steady stream of electrons and found that electrons with energies just outside bands of excited molecular states can enter the chain through phonon emission or absorption. Furthermore, this phonon-assisted dynamical behaviour shows great sensitivity to the vibrational temperature, highlighting a crucial controlling factor for the injection and propagation of electrons in water. △ Less

Submitted 23 February, 2017; v1 submitted 23 December, 2016; originally announced December 2016.

Journal ref: Scientific Reports 7, Article number: 45410 (2017)