Spatially Resolved Simulations of Membrane Reactions and Dynamics: Multipolar Reaction DPD

Biophysical chemistry of mesoscale systems and quantitative modeling in systems biology now require a simulation methodology unifying chemical reaction kinetics with essential collective physics. This will enable the study of the collective dynamics of complex chemical and structural systems in a spatially resolved manner with a combinatorially complex variety of different system constituents. In order to allow a direct link-up with experimental data (e.g. high-throughput fluorescence images) the simulations must be constructed locally, i.e. mesoscale phenomena have to emerge from local composition and interactions that can be extracted from experimental data. Under suitable conditions, the simulation of such local interactions must lead to processes such as vesicle budding, transport of membrane-bounded compartments and protein sorting, all of which result from a sophisticated interplay between chemical and mechanical processes and require the link-up of different length scales. In this work, we show that introducing multipolar interactions between particles in dissipative particle dynamics (DPD) leads to extended membrane structures emerging in a self-organized manner and exhibiting the necessary mechanical stability for transport, correct scaling behavior, and membrane fluidity so as to provide a two-dimensional self-organizing dynamic reaction environment for kinetic studies in the context of cell biology.