Kinetics of self-assembly of inclusions due to lipid membrane thickness interactions

Self-assembly of proteins on lipid membranes underlies many important processes in cell biology, such as, exo- and endo-cytosis, assembly of viruses, etc. An attractive force that can cause self-assembly is mediated by membrane thickness interactions between proteins. The free energy profile associated with this attractive force is a result of the overlap of thickness deformation fields around the proteins which can be calculated from the solution of a boundary value problem. Yet, the time scales over which two inclusions coalesce has not been explored, even though the evolution of particle concentrations on membranes has been modeled using phase-field approaches. In this paper we compute this time scale as a function of the initial distance between two inclusions by viewing their coalescence as a first passage time problem. The mean first passage time is computed using Langevin dynamics and a partial differential equation, and both methods are found to be in excellent agreement. Inclusions of three different shapes are studied and it is found that for two inclusions separated by about hundred nanometers the time to coalescence is hundreds of milliseconds irrespective of shape. An efficient computation of the interaction energy of inclusions is central to our work. We compute it using a finite difference technique and show that our results are in excellent agreement with those from a previously proposed semi-analytical method based on Fourier-Bessel series. The computational strategies described in this paper could potentially lead to efficient methods to explore the kinetics of self-assembly of proteins on lipid membranes.