Source author record

Gianni De Fabritiis

Gianni De Fabritiis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2022arXiv

TorchMD-NET: Equivariant Transformers for Neural Network based Molecular Potentials

The prediction of quantum mechanical properties is historically plagued by a trade-off between accuracy and speed. Machine learning potentials have previously shown great success in this domain, reaching increasingly better accuracy while maintaining computational efficiency comparable with classical force fields. In this work we propose TorchMD-NET, a novel equivariant transformer (ET) architecture, outperforming state-of-the-art on MD17, ANI-1, and many QM9 targets in both accuracy and computational efficiency. Through an extensive attention weight analysis, we gain valuable insights into the black box predictor and show differences in the learned representation of conformers versus conformations sampled from molecular dynamics or normal modes. Furthermore, we highlight the importance of datasets including off-equilibrium conformations for the evaluation of molecular potentials.

preprint2020arXiv

AdaptiveBandit: A multi-armed bandit framework for adaptive sampling in molecular simulations

Sampling from the equilibrium distribution has always been a major problem in molecular simulations due to the very high dimensionality of conformational space. Over several decades, many approaches have been used to overcome the problem. In particular, we focus on unbiased simulation methods such as parallel and adaptive sampling. Here, we recast adaptive sampling schemes on the basis of multi-armed bandits and develop a novel adaptive sampling algorithm under this framework, \UCB. We test it on multiple simplified potentials and in a protein folding scenario. We find that this framework performs similarly or better in every type of test potentials compared to previous methods. Furthermore, it provides a novel framework to develop new sampling algorithms with better asymptotic characteristics.

preprint2020arXiv

Coarse Graining Molecular Dynamics with Graph Neural Networks

Coarse graining enables the investigation of molecular dynamics for larger systems and at longer timescales than is possible at atomic resolution. However, a coarse graining model must be formulated such that the conclusions we draw from it are consistent with the conclusions we would draw from a model at a finer level of detail. It has been proven that a force matching scheme defines a thermodynamically consistent coarse-grained model for an atomistic system in the variational limit. Wang et al. [ACS Cent. Sci. 5, 755 (2019)] demonstrated that the existence of such a variational limit enables the use of a supervised machine learning framework to generate a coarse-grained force field, which can then be used for simulation in the coarse-grained space. Their framework, however, requires the manual input of molecular features upon which to machine learn the force field. In the present contribution, we build upon the advance of Wang et al.and introduce a hybrid architecture for the machine learning of coarse-grained force fields that learns their own features via a subnetwork that leverages continuous filter convolutions on a graph neural network architecture. We demonstrate that this framework succeeds at reproducing the thermodynamics for small biomolecular systems. Since the learned molecular representations are inherently transferable, the architecture presented here sets the stage for the development of machine-learned, coarse-grained force fields that are transferable across molecular systems.

preprint2020arXiv

SkeleDock: A Web Application for Scaffold Docking in PlayMolecule

SkeleDock is a scaffold docking algorithm which uses the structure of a protein-ligand complex as a template to model the binding mode of a chemically similar system. This algorithm was evaluated in the D3R Grand Challenge 4 pose prediction challenge, where it achieved competitive performance. Furthermore, we show that, if crystallized fragments of the target ligand are available, SkeleDock can outperform rDock docking software at predicting the binding mode. This article also addresses the capacity of this algorithm to model macrocycles and deal with scaffold hopping. SkeleDock can be accessed at https://playmolecule.org/SkeleDock/.

preprint2013arXiv

Identification of slow molecular order parameters for Markov model construction

A goal in the kinetic characterization of a macromolecular system is the description of its slow relaxation processes, involving (i) identification of the structural changes involved in these processes, and (ii) estimation of the rates or timescales at which these slow processes occur. Most of the approaches to this task, including Markov models, Master-equation models, and kinetic network models, start by discretizing the high-dimensional state space and then characterize relaxation processes in terms of the eigenvectors and eigenvalues of a discrete transition matrix. The practical success of such an approach depends very much on the ability to finely discretize the slow order parameters. How can this task be achieved in a high-dimensional configuration space without relying on subjective guesses of the slow order parameters? In this paper, we use the variational principle of conformation dynamics to derive an optimal way of identifying the "slow subspace" of a large set of prior order parameters - either generic internal coordinates (distances and dihedral angles), or a user-defined set of parameters. It is shown that a method to identify this slow subspace exists in statistics: the time-lagged independent component analysis (TICA). Furthermore, optimal indicators-order parameters indicating the progress of the slow transitions and thus may serve as reaction coordinates-are readily identified. We demonstrate that the slow subspace is well suited to construct accurate kinetic models of two sets of molecular dynamics simulations, the 6-residue fluorescent peptide MR121-GSGSW and the 30-residue natively disordered peptide KID. The identified optimal indicators reveal the structural changes associated with the slow processes of the molecular system under analysis.

preprint2004arXiv

Energy controlled insertion of polar molecules in dense fluids

We present a method to search low energy configurations of polar molecules in the complex potential energy surfaces associated with dense fluids. The search is done in the configurational space of the translational and rotational degrees of freedom of the molecule, combining steepest-descent and Newton-Raphson steps which embed information on the average sizes of the potential energy wells obtained from prior inspection of the liquid structure. We perform a molecular dynamics simulation of a liquid water shell which demonstrates that the method enables fast and energy-controlled water molecule insertion in aqueous environments. The algorithm finds low energy configurations of incoming water molecules around three orders of magnitude faster than direct random insertion. This method is an important step towards dynamic simulations of open systems and it may also prove useful for energy-biased ensemble average calculations of the chemical potential.