Bilayer-thickness-mediated interactions between integral membrane proteins

TL;DR

This paper develops an analytic and numerical framework to accurately predict bilayer-mediated elastic interactions between integral membrane proteins of arbitrary shapes and orientations, considering membrane elasticity and protein-induced deformations.

Contribution

It introduces a combined analytic and numerical methodology for solving membrane elastic equations, enabling precise modeling of protein interactions in complex membrane environments.

Findings

Exact solutions for cylindrical proteins with varying thickness

Perturbative solutions for non-cylindrical shapes

Validated numerical schemes for complex protein geometries

Abstract

Hydrophobic thickness mismatch between integral membrane proteins and the surrounding lipid bilayer can produce lipid bilayer thickness deformations. Experiment and theory have shown that protein-induced lipid bilayer thickness deformations can yield energetically favorable bilayer-mediated interactions between integral membrane proteins, and large-scale organization of integral membrane proteins into protein clusters in cell membranes. Within the continuum elasticity theory of membranes, the energy cost of protein-induced bilayer thickness deformations can be captured by considering compression and expansion of the bilayer hydrophobic core, membrane tension, and bilayer bending, resulting in biharmonic equilibrium equations describing the shape of lipid bilayers for a given set of bilayer-protein boundary conditions. Here we develop a combined analytic and numerical methodology for the solution of the equilibrium elastic equations associated with protein-induced lipid bilayer deformations. Our methodology allows accurate prediction of thickness-mediated protein interactions for arbitrary protein symmetries at arbitrary protein separations and relative orientations. We provide exact analytic solutions for cylindrical integral membrane proteins with constant and varying hydrophobic thickness, and develop perturbative analytic solutions for non-cylindrical protein shapes. We complement these analytic solutions, and assess their accuracy, by developing both finite element and finite difference numerical solution schemes. Taken together, the work presented here puts into place an analytic and numerical framework which allows calculation of bilayer-mediated elastic interactions between integral membrane proteins for the complicated protein shapes suggested by structural biology and at the small protein separations most relevant for the crowded membrane environments provided by living cells.