Source author record

Michele Vendruscolo

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

ResearcherUnclaimed source record

Biological Physics cond-mat.soft physics.comp-ph Biomolecules cond-mat.stat-mech Molecular Networks physics.chem-ph physics.data-an physics.flu-dyn Quantitative Methods

Catalog footprint

What is connected

5works

10topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2020arXiv

A method of incorporating rate constants as kinetic constraints in molecular dynamics simulations

From the point of view of statistical mechanics, a full characterisation of a molecular system requires the experimental determination of its possible states, their populations and the respective interconversion rates. Well-established methods can incorporate in molecular dynamics simulations experimental information about states using structural restraints, and about populations using thermodynamic restraints. However, it is still unclear how to include experimental knowledge of interconversion rates. Here we introduce a method of imposing known rate constants as constraints in molecular dynamics simulations, which is based on a combination of the maximum entropy and maximum caliber principles. Starting from an existing ensemble of trajectories, obtained from either molecular dynamics or enhanced trajectory sampling, this method provides a minimally perturbed path distribution consistent with the kinetic constraints, as well as a modified free energy and committor landscape. We illustrate the application of the method to simple toy systems, as well as to all atom molecular simulations of peptide association and folding. We find that by combining experimental rate coefficient data and molecular dynamics simulations we are able to determine new transition states, reaction mechanisms and free energies. For instance, in the case of chignolin protein folding we find that imposing a slower folding rate shifts the transition state to more native like conformations, while it increases the stability of the unfolded region.

preprint2015arXiv

Metainference: A Bayesian Inference Method for Heterogeneous Systems

Modelling a complex system is almost invariably a challenging task. The incorporation of experimental observations can be used to improve the quality of a model, and thus to obtain better predictions about the behavior of the corresponding system. This approach, however, is affected by a variety of different errors, especially when a system populates simultaneously an ensemble of different states and experimental data are measured as averages over such states. To address this problem we present a Bayesian inference method, called metainference, that is able to deal with errors in experimental measurements as well as with experimental measurements averaged over multiple states. To achieve this goal, metainference models a finite sample of the distribution of models using a replica approach, in the spirit of the replica-averaging modelling based on the maximum entropy principle. To illustrate the method we present its application to a heterogeneous model system and to the determination of an ensemble of structures corresponding to the thermal fluctuations of a protein molecule. Metainference thus provides an approach to model complex systems with heterogeneous components and interconverting between different states by taking into account all possible sources of errors.

preprint2015arXiv

Particle-based simulations of steady-state mass transport at high Péclet numbers

Conventional approaches for simulating steady-state distributions of particles under diffusive and advective transport at high Péclet numbers involve solving the diffusion and advection equations in at least two dimensions. Here, we present an alternative computational strategy by combining a particle-based rather than a field-based approach with the initialisation of particles in proportion to their flux. This method allows accurate prediction of the steady state and is applicable even at high Péclet numbers where traditional particle-based Monte-Carlo methods starting from randomly initialised particle distributions fail. We demonstrate that generating a flux of particles according to a predetermined density and velocity distribution at a single fixed time and initial location allows for accurate simulation of mass transport under flow. Specifically, upon initialisation in proportion to their flux, these particles are propagated individually and detected by summing up their Monte-Carlo trajectories in predefined detection regions. We demonstrate quantitative agreement of the predicted concentration profiles with the results of experiments performed with fluorescent particles in microfluidic channels under continuous flow. This approach is computationally advantageous and readily allows non-trivial initial distributions to be considered. In particular, this method is highly suitable for simulating advective and diffusive transport in microfluidic devices.

preprint2014arXiv

Analysis of the hierarchical structure of the B. subtilis transcriptional regulatory network

The transcriptional regulation of gene expression is orchestrated by complex networks of interacting genes. Increasing evidence indicates that these transcriptional regulatory networks (TRNs) in bacteria have an inherently hierarchical architecture, although the design principles and the specific advantages offered by this type of organization have not yet been fully elucidated. In this study, we focussed on the hierarchical structure of the TRN of the gram-positive bacterium Bacillus subtilis and performed a comparative analysis with the TRN of the gram-negative bacterium Escherichia coli. Using a graph-theoretic approach, we organized the transcription factors (TFs) and sigma-factors in the TRNs of B. subtilis and E. coli into three hierarchical levels (Top, Middle and Bottom) and studied several structural and functional properties across them. In addition to many similarities, we found also specific differences, explaining the majority of them with variations in the distribution of sigma-factors across the hierarchical levels in the two organisms. We then investigated the control of target metabolic genes by transcriptional regulators to characterize the differential regulation of three distinct metabolic subsystems (catabolism, anabolism and central energy metabolism). These results suggest that the hierarchical architecture that we observed in B. subtilis represents an effective organization of its TRN to achieve flexibility in the response to diverse stimuli.

preprint2011arXiv

Nucleated polymerisation in the presence of pre-formed seed filaments

We revisit the classical problem of nucleated polymerisation and derive a range of exact results describing polymerisation in systems intermediate between the well-known limiting cases of a reaction starting from purely soluble material and for a reaction where no new growth nuclei are formed.