BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Raffaello Potestio

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

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softBiological PhysicsBiomoleculesphysics.chem-phQuantitative MethodsMachine Learningphysics.comp-ph

Catalog footprint

What is connected

13works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

13 published item(s)

preprint2023arXiv

Fast, accurate, and system-specific variable-resolution modelling of proteins

In recent years, a few multiple-resolution modelling strategies have been proposed, in which functionally relevant parts of a biomolecule are described with atomistic resolution, while the remainder of the system is concurrently treated using a coarse-grained model. In most cases, the parametrisation of the latter requires lengthy reference all-atom simulations and/or the usage of off-shelf coarse-grained force fields, whose interactions have to be refined to fit the specific system under examination. Here, we overcome these limitations through a novel multi-resolution modelling scheme for proteins, dubbed coarse-grained anisotropic network model for variable resolution simulations, or CANVAS. This scheme enables the user-defined modulation of the resolution level throughout the system structure; a fast parametrisation of the potential without the necessity of reference simulations; and the straightforward usage of the model on the most commonly used molecular dynamics platforms. The method is presented and validated on two case studies, the enzyme adenylate kinase and the therapeutic antibody pembrolizumab, by comparing results obtained with the CANVAS model against fully atomistic simulations. The modelling software, implemented in python, is made freely available for the community on a collaborative github repository.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Information-theoretical measures identify accurate low-resolution representations of protein configurational space

A steadily growing computational power is employed to perform molecular dynamics simulations of biological macromolecules, which represents at the same time an immense opportunity and a formidable challenge. In fact, large amounts of data are produced, from which useful, synthetic, and intelligible information has to be extracted to make the crucial step from knowing to understanding. Here we tackled the problem of coarsening the conformational space sampled by proteins in the course of molecular dynamics simulations. We applied different schemes to cluster the frames of a dataset of protein simulations; we then employed an information-theoretical framework, based on the notion of resolution and relevance, to gauge how well the various clustering methods accomplish this simplification of the configurational space. Our approach allowed us to identify the level of resolution that optimally balances simplicity and informativeness; furthermore, we found that the most physically accurate clustering procedures are those that induce an ultrametric structure of the low-resolution space, consistently with the hypothesis that the protein conformational landscape has a self-similar organisation. The proposed strategy is general and its applicability extends beyond that of computational biophysics, making it a valuable tool to extract useful information from large datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accelerating the identification of informative reduced representations of proteins with deep learning for graphs

The limits of molecular dynamics (MD) simulations of macromolecules are steadily pushed forward by the relentless developments of computer architectures and algorithms. This explosion in the number and extent (in size and time) of MD trajectories induces the need of automated and transferable methods to rationalise the raw data and make quantitative sense out of them. Recently, an algorithmic approach was developed by some of us to identify the subset of a protein's atoms, or mapping, that enables the most informative description of it. This method relies on the computation, for a given reduced representation, of the associated mapping entropy, that is, a measure of the information loss due to the simplification. Albeit relatively straightforward, this calculation can be time consuming. Here, we describe the implementation of a deep learning approach aimed at accelerating the calculation of the mapping entropy. The method relies on deep graph networks, which provide extreme flexibility in the input format. We show that deep graph networks are accurate and remarkably efficient, with a speedup factor as large as $10^5$ with respect to the algorithmic computation of the mapping entropy. Applications of this method, which entails a great potential in the study of biomolecules when used to reconstruct its mapping entropy landscape, reach much farther than this, being the scheme easily transferable to the computation of arbitrary functions of a molecule's structure.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An information theory-based approach for optimal model reduction of biomolecules

In the theoretical modelling of a physical system a crucial step consists in the identification of those degrees of freedom that enable a synthetic, yet informative representation of it. While in some cases this selection can be carried out on the basis of intuition and experience, a straightforward discrimination of the important features from the negligible ones is difficult for many complex systems, most notably heteropolymers and large biomolecules. We here present a thermodynamics-based theoretical framework to gauge the effectiveness of a given simplified representation by measuring its information content. We employ this method to identify those reduced descriptions of proteins, in terms of a subset of their atoms, that retain the largest amount of information from the original model; we show that these highly informative representations share common features that are intrinsically related to the biological properties of the proteins under examination, thereby establishing a bridge between protein structure, energetics, and function.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ligand-protein interactions in lysozyme investigated through a dual-resolution model

A fully atomistic modelling of biological macromolecules at relevant length- and time-scales is often cumbersome or not even desirable, both in terms of computational effort required and it a posteriori analysis. This difficulty can be overcome with the use of multi-resolution models, in which different regions of the same system are concurrently described at different levels of detail. In enzymes, computationally expensive atomistic detail is crucial in the modelling of the active site in order to capture e.g. the chemically subtle process of ligand binding. In contrast, important yet more collective properties of the remainder of the protein can be reproduced with a coarser description. In the present work, we demonstrate the effectiveness of this approach through the calculation of the binding free energy of hen egg white lysozyme (HEWL) with the inhibitor di-N-acetylchitotriose. Particular attention is posed to the impact of the mapping, i.e. the selection of atomistic and coarse-grained residues, on the binding free energy. It is shown that, in spite of small variations of the binding free energy with respect to the active site resolution, the separate contributions coming from different energetic terms (such as electrostatic and van der Waals interactions) manifest a stronger dependence on the mapping, thus pointing to the existence of an optimal level of intermediate resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Open-Boundary Hamiltonian adaptive resolution. From grand canonical to non-equilibrium molecular dynamics simulations

We propose an open-boundary molecular dynamics method in which an atomistic system is in contact with an infinite particle reservoir at constant temperature, volume and chemical potential. In practice, following the Hamiltonian adaptive resolution strategy, the system is partitioned into a domain of interest and a reservoir of non-interacting, ideal gas, particles. An external potential, applied only in the interfacial region, balances the excess chemical potential of the system. To ensure that the size of the reservoir is infinite, we introduce a particle insertion/deletion algorithm to control the density in the ideal gas region. We show that it is possible to study non-equilibrium phenomena with this open-boundary molecular dynamics method. To this aim, we consider a prototypical confined liquid under the influence of an external constant density gradient. The resulting pressure-driven flow across the atomistic system exhibits a velocity profile consistent with the corresponding solution of the Navier-Stokes equation. In contrast to available computational methods in which external forces drive the system far from equilibrium, this approach conserves momentum and closely resembles experimental conditions. The presented method can be used to study various direct and indirect out-of-equilibrium conditions in complex molecular systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A multi-resolution model to capture both global fluctuations of an enzyme and molecular recognition in the ligand-binding site

In multi-resolution simulations, different system components are simultaneously modelled at different levels of resolution, these being smoothly coupled together. In the case of enzyme systems, computationally expensive atomistic detail is needed in the active site to capture the chemistry of substrate binding. Global properties of the rest of the protein also play an essential role, determining the structure and fluctuations of the binding site; however, these can be modelled on a coarser level. Similarly, in the most computationally efficient scheme only the solvent hydrating the active site requires atomistic detail. We present a methodology to couple atomistic and coarse-grained protein models, while solvating the atomistic part of the protein in atomistic water. This allows a free choice of which protein and solvent degrees of freedom to include atomistically, without loss of accuracy in the atomistic description. This multi-resolution methodology can successfully model stable ligand binding, and we further confirm its validity via an exploration of system properties relevant to enzymatic function. In addition to a computational speedup, such an approach can allow the identification of the essential degrees of freedom playing a role in a given process, potentially yielding new insights into biomolecular function.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Chirality modifies the interaction between knots

In this study we consider an idealization of a typical optical tweezers experiment involving a semiflexible double-knotted polymer, with steric hindrance and persistence length matching those of dsDNA in high salt concentration, under strong stretching. Using exhaustive Molecular Dynamics simulations we show that not only does a double-knotted dsDNA filament under tension possess a free energy minimum when the two knots are intertwined, but also that the depth of this minimum depends on the relative chirality of the two knots. We rationalize this dependence of the effective interaction on the chirality in terms of a competition between chain entropy and bending energy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Advantages and challenges in coupling an ideal gas to atomistic models in adaptive resolution simulations

In adaptive resolution simulations, molecular fluids are modeled employing different levels of resolution in different subregions of the system. When traveling from one region to the other, particles change their resolution on the fly. One of the main advantages of such approaches is the computational efficiency gained in the coarse-grained region. In this respect the best coarse-grained system to employ in the low resolution region would be the ideal gas, making intermolecular force calculations in the coarse-grained subdomain redundant. In this case, however, a smooth coupling is challenging due to the high energetic imbalance between typical liquids and a system of non-interacting particles. In the present work, we investigate this approach, using as a test case the most biologically relevant fluid, water. We demonstrate that a successful coupling of water to the ideal gas can be achieved with current adaptive resolution methods, and discuss the issues that remain to be addressed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

From classical to quantum and back: Hamiltonian coupling of classical and Path Integral models of atoms

In computer simulations, quantum delocalization of atomic nuclei can be modeled making use of the Path Integral (PI) formulation of quantum statistical mechanics. This approach, however, comes with a large computational cost. By restricting the PI modeling to a small region of space, this cost can be significantly reduced. In the present work we derive a Hamiltonian formulation for a bottom-up, theoretically solid simulation protocol that allows molecules to change their resolution from quantum-mechanical to classical and vice versa on the fly, while freely diffusing across the system. This approach renders possible simulations of quantum systems at constant chemical potential. The validity of the proposed scheme is demonstrated by means of simulations of low temperature parahydrogen. Potential future applications include simulations of biomolecules, membranes, and interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Monte Carlo adaptive resolution simulation of multicomponent molecular liquids

Complex soft matter systems can be efficiently studied with the help of adaptive resolution simulation methods, concurrently employing two levels of resolution in different regions of the simulation domain. The non-matching properties of high- and low-resolution models, however, lead to thermodynamic imbalances between the system's subdomains. Such inhomogeneities can be healed by appropriate compensation forces, whose calculation requires nontrivial iterative procedures. In this work we employ the recently developed Hamiltonian Adaptive Resolution Simulation method to perform Monte Carlo simulations of a binary mixture, and propose an efficient scheme, based on Kirkwood Thermodynamic Integration, to regulate the thermodynamic balance of multi-component systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hamiltonian adaptive resolution simulation for molecular liquids

Adaptive resolution schemes allow the simulation of a molecular fluid treating simultaneously different subregions of the system at different levels of resolution. In this work we present a new scheme formulated in terms of a global Hamiltonian. Within this approach equilibrium states corresponding to well defined statistical ensembles can be generated making use of all standard Molecular Dynamics or Monte Carlo methods. Models at different resolutions can thus be coupled, and thermodynamic equilibrium can be modulated keeping each region at desired pressure or density without disrupting the Hamiltonian framework.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Random Matrix approach to collective behavior and bulk universality in protein dynamics

Covariance matrices of amino acid displacements, commonly used to characterize the large-scale movements of proteins, are investigated through the prism of Random Matrix Theory. Bulk universality is detected in the local spacing statistics of noise-dressed eigenmodes, which is well described by a Brody distribution with parameter $β\simeq 0.8$. This finding, supported by other consistent indicators, implies a novel quantitative criterion to single out the collective degrees of freedom of the protein from the majority of high-energy, localized vibrations.

0 citations0 reviewsNo code linkedOpen access