Search

Allosteric regulation is a widespread strategy employed by several proteins to transduce chemical signals and perform biological functions. Metal sensor proteins are exemplary in this respect, e.g., in that they selectively bind and unbind DNA depending on the state of a distal ion coordination site. In this work, we carry out an investigation of the structural and mechanical properties of the CzrA transcription repressor through the analysis of microsecond-long molecular dynamics (MD) trajectories; the latter are processed through the mapping entropy optimisation workflow (MEOW), a recently developed information-theoretical method that highlights, in an unsupervised manner, residues of particular mechanical, functional, and biological importance. This approach allows us to unveil how differences in the properties of the molecule are controlled by the state of the zinc coordination site, with particular attention to the DNA binding region. These changes correlate with a redistribution of the conformational variability of the residues throughout the molecule, in spite of an overall consistency of its architecture in the two (ion-bound and free) coordination states. The results of this work corroborate previous studies, provide novel insight into the fine details of the mechanics of CzrA, and showcase the MEOW approach as a novel instrument for the study of allosteric regulation and other processes in proteins through the analysis of plain MD simulations.

The analysis of complex systems such as neural networks is made particularly difficult by the overwhelming number of their interacting components. In the absence of prior knowledge, identifying a small but informative subset of network nodes on which the analysis should focus is a rather challenging task. In this work, we address this problem in the context of a Hopfield model, which is observed through the lenses of low-resolution representations, or decimation mappings, consisting of subgroups of its neurons. The optimal, most informative mappings of the network are defined through a recently developed methodology, the mapping entropy optimisation workflow (MEOW), which performs an unsupervised analysis of the states sampled by the network and identifies those subgroups of spins whose configuration distribution is closest to that of the full, high-resolution model. Which neurons are retained in an optimal mapping is found to critically depend on the properties of the interaction matrix of the network and the level of detail employed to describe the system; by these means, it is thus possible to extract quantitative insight about the underlying properties of the high-resolution model through the analysis of its optimal low-resolution representations. These results show a tight and potentially fruitful relation between the level of detail at which the network is inspected and the type and amount of information that can be gathered from it, and showcase the MEOW approach as a practical, enabling tool for the study of complex systems.

Bottom-up coarse-grained (CG) models proved to be essential to complement and sometimes even replace all-atom representations of soft matter systems and biological macromolecules. The development of low-resolution models takes the moves from the reduction of the degrees of freedom employed, that is, the definition of a mapping between a system's high-resolution description and its simplified counterpart. Even in the absence of an explicit parametrisation and simulation of a CG model, the observation of the atomistic system in simpler terms can be informative: this idea is leveraged by the mapping entropy, a measure of the information loss inherent to the process of coarsening. Mapping entropy lies at the heart of the extensible coarse-graining toolbox, or EXCOGITO, developed to perform a number of operations and analyses on molecular systems pivoting around the properties of mappings. EXCOGITO can process an all-atom trajectory to compute the mapping entropy, identify the mapping that minimizes it, and establish quantitative relations between a low-resolution representation and the geometrical, structural, and energetic features of the system. Here, the software, which is available free of charge under an open-source licence, is presented and showcased to introduce potential users to its capabilities and usage. Published on the J. Chem. Inf. Model. on June 11, 2024. DOI: https://doi.org/10.1021/acs.jcim.4c00490

Coarse-grained (CG) models are simplified representations of soft matter systems that are commonly employed to overcome size and time limitations in computational studies. Many approaches have been developed to construct and parametrise such effective models for a variety of systems of natural as well as artificial origin. However, while extremely accurate in reproducing the stationary and equilibrium observables obtained with more detailed representations, CG models generally fail to preserve the original time scales of the reference system, and hence its dynamical properties. In order to improve our understanding of the impact of coarse-graining on the model system dynamics, we here formulate the Mori-Zwanzig generalised Langevin equations (GLEs) of motion of a CG model in terms of a time non-local stationary-action principle. The latter is employed in combination with a data-driven optimisation strategy to determine the parameters of the GLE. We apply this approach to a system of water molecules in standard thermodynamical conditions, showing that it can substantially improve the dynamical features of the corresponding CG model., Comment: 10 pages, 3 figures. Added references and stability analyses. Results unchanged

Simplified representations of macromolecules help in rationalising and understanding the outcome of atomistic simulations, and serve to the construction of effective, coarse-grained models. The number and distribution of coarse-grained sites bears a strict relation with the amount of information conveyed by the representation and the accuracy of the associated effective model; in this work, we investigate this relationship from the very basics: specifically, we propose a rigorous notion of scalar product among mappings, which implies a distance and a metric space of simplified representations. Making use of a Wang-Landau enhanced sampling algorithm, we exhaustively explore the space of mappings, quantifying their qualitative features in terms of their squared norm and relating them with thermodynamical properties of the underlying macromolecule. A one-to-one correspondence with an interacting lattice gas on a finite volume leads to the emergence of discontinuous phase transitions in mapping space that mark the boundaries between qualitatively different representations of the same molecule.

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.

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.

The size of chemical compound space is too large to be probed exhaustively. This leads high-throughput protocols to drastically subsample and results in sparse and non-uniform datasets. Rather than arbitrarily selecting compounds, we systematically explore chemical space according to the target property of interest. We first perform importance sampling by introducing a Markov chain Monte Carlo scheme across compounds. We then train an ML model on the sampled data to expand the region of chemical space probed. Our boosting procedure enhances the number of compounds by a factor 2 to 10, enabled by the ML model's coarse-grained representation, which both simplifies the structure-property relationship and reduces the size of chemical space. The ML model correctly recovers linear relationships between transfer free energies. These linear relationships correspond to features that are global to the dataset, marking the region of chemical space up to which predictions are reliable---a more robust alternative to the predictive variance. Bridging coarse-grained simulations with ML gives rise to an unprecedented database of drug-membrane insertion free energies for 1.3 million compounds., Comment: 9 pages, 5 figures

Unraveling the relation between the chemical structure of small drug-like compounds and their rate of passive permeation across lipid membranes is of fundamental importance for pharmaceutical applications. The elucidation of a comprehensive structure-permeability relationship expressed in terms of a few molecular descriptors is unfortunately hampered by the overwhelming number of possible compounds. In this work, we reduce a priori the size and diversity of chemical space to solve an analogous---but smoothed out---structure-property relationship problem. This is achieved by relying on a physics-based coarse-grained model that reduces the size of chemical space, enabling a comprehensive exploration of this space with greatly reduced computational cost. We perform high-throughput coarse-grained (HTCG) simulations to derive a permeability surface in terms of two simple molecular descriptors---bulk partitioning free energy and pKa. The surface is constructed by exhaustively simulating all coarse-grained compounds that are representative of small organic molecules (ranging from 30 to 160 Da) in a high-throughput scheme. We provide results for acidic, basic and zwitterionic compounds. Connecting back to the atomic resolution, the HTCG predictions for more than 500,000 compounds allow us to establish a clear connection between specific chemical groups and the resulting permeability coefficient, enabling for the first time an inverse design procedure. Our results have profound implications for drug synthesis: the predominance of commonly-employed chemical moieties narrows down the range of permeabilities., Comment: 9 pages, 4 figures (SI: 23 pages; 9 figures)

The determination of potentials of mean force for solute insertion in a membrane by means of all-atom molecular dynamics simulations is often hampered by sampling issues. A multiscale approach to conformational sampling was recently proposed by Bereau and Kremer (2016). It aims at accelerating the sampling of the atomistic conformational space by means of a systematic backmapping of coarse-grained snapshots. In this work, we first analyze the efficiency of this method by comparing its predictions for propanol insertion into a 1,2-Dimyristoyl-sn-glycero-3-phosphocholine membrane (DMPC) against reference atomistic simulations. The method is found to provide accurate results with a gain of one order of magnitude in computational time. We then investigate the role of the coarse-grained representation in affecting the reliability of the method in the case of a 1,2-Dioleoyl-sn-glycero-3-phosphocholine membrane (DOPC). We find that the accuracy of the results is tightly connected to the presence a good configurational overlap between the coarse-grained and atomistic models---a general requirement when developing multiscale simulation methods., Comment: 6 pages, 5 figures

The partitioning of small molecules in cell membranes---a key parameter for pharmaceutical applications---typically relies on experimentally-available bulk partitioning coefficients. Computer simulations provide a structural resolution of the insertion thermodynamics via the potential of mean force, but require significant sampling at the atomistic level. Here, we introduce high-throughput coarse-grained molecular dynamics simulations to screen thermodynamic properties. This application of physics based models in a large-scale study of small molecules establishes linear relationships between partitioning coefficients and key features of the potential of mean force. This allows us to predict the structure of the insertion from bulk experimental measurements for more than 400,000 compounds. The potential of mean force hereby becomes an easily accessible quantity---already recognized for its high predictability of certain properties, e.g., passive permeation. Further, we demonstrate how coarse graining helps reduce the size of chemical space, enabling a hierarchical approach to screening small molecules., Comment: 8 pages, 6 figures. Typos fixed, minor corrections

We consider a coarse-grained (CG) model with pairwise interactions, suitable to describe low-density solutions of star-branched polymers of functionality $f$. Each macromolecule is represented by a CG molecule with $(f+1)$ interaction sites, which captures the star topology. Potentials are obtained by requiring the CG model to reproduce a set of distribution functions computed in the microscopic model in the zero-density limit. Explicit results are given for $f=6,12$ and $40$. We use the CG model to compute the osmotic equation of state of the solution for concentrations $c$ such that $\Phi_p = c/c^* \lesssim 1$, where $c^*$ is the overlap concentration. We also investigate in detail the phase diagram for f=40, identifying the boundaries of the solid intermediate phase. Finally, we investigate how the polymer size changes with $c$. For $\Phi_p\lesssim 0.3$ polymers become harder as $f$ increases at fixed reduced concentration $c/c^*$. On the other hand, for $\Phi_p\gtrsim 0.3$, polymers show the opposite behavior: At fixed $\Phi_p$, the larger the value of $f$, the larger their size reduction is., Comment: 29 pages (17 pages main paper, 12 pages supplementary material), 24 figures, 8 tables. Manuscript accepted for publication in J. Chem. Phys

We discuss the reliability of integral-equation methods based on several commonly used closure relations in determining the phase diagram of coarse-grained models of soft-matter systems characterized by mutually interacting soft and hard-core particles. Specifically, we consider a set of potentials appropriate to describe a system of hard-sphere colloids and linear homopolymers in good solvent, and investigate the behavior when the soft particles are smaller than the colloids, which is the regime of validity of the coarse-grained models. Using computer-simulation results as a benchmark, we find that the hypernetted-chain approximation provides accurate estimates of thermodynamics and structure in the colloid-gas phase in which the density of colloids is small. On the other hand, all closures considered appear to be unable to describe the behavior of the mixture in the colloid-liquid phase, as they cease to converge at polymer densities significantly smaller than those at the binodal. As a consequence, integral equations appear to be unable to predict a quantitatively correct phase diagram., Comment: 16 pages, 11 figures, 3 tables

We critically discuss and review the general ideas behind single- and multi-site coarse-grained (CG) models as applied to macromolecular solutions in the dilute and semi-dilute regime. We first consider single-site models with zero-density and density-dependent pair potentials. We highlight advantages and limitations of each option in reproducing the thermodynamic behavior and the large-scale structure of the underlying reference model. As a case study we consider solutions of linear homopolymers in a solvent of variable quality. Secondly, we extend the discussion to multi-component systems presenting, as a test case, results for mixtures of colloids and polymers. Specifically, we found the CG model with zero-density potentials to be unable to predict fluid-fluid demixing in a reasonable range of densities for mixtures of colloids and polymers of equal size. For larger colloids, the polymer volume fractions at which phase separation occurs are largely overestimated. CG models with density-dependent potentials are somewhat less accurate than models with zero-density potentials in reproducing the thermodynamics of the system and, although they presents a phase separation, they significantly underestimate the polymer volume fractions along the binodal. Finally, we discuss a general multi-site strategy, which is thermodynamically consistent and fully transferable with the number of sites, and that allows us to overcome most of the limitations discussed for single-site models., Comment: 23 pages, 9 figures, 4 tables

We compare different coarse-grained models for star polymers. We find that phenomenological models inspired by the Daoud-Cotton model reproduce quite poorly the thermodynamics of these systems, even if the potential is assumed to be density dependent, as done in the analysis of experimental results. We also determine the minumum value fc of the functionality of the star polymer for which a fluid-solid transition occurs. By applying the Hansen-Verlet criterion we find 35 < fc < 40. This result is confirmed by an analysis based on the modified (reference) hypernetted chain method and is qualitatively consistent with previous work., Comment: 9 pages. In the new version, comments added and a few typos corrected. To appear in J. Chem. Phys