Membrane Thickness Strain from Protein Inclusions: A Multiscale Simulation and X-Ray Scattering Study of Proteoliposomes

Abstract

Integral membrane proteins remodel the surrounding lipid bilayer, but quantifying the resulting deformations and linking them to protein density in the membrane has remained challenging. Here, we introduce an integrative methodology that combines all-atom molecular dynamics (MD) simulations with multiscale small-angle X-ray scattering (SAXS) analysis to connect membrane strain to the protein/lipid ratio in proteoliposomes. Using outer membrane phospholipase A (OmpLA) reconstituted into lipid bilayers with both increased and decreased hydrophobic thickness, we systematically probe the effects of positive and negative hydrophobic mismatch. MD simulations demonstrate that OmpLA causes anisotropic, oscillatory thickness deformations extending up to eight times the radius of the first lipid shell surrounding the protein, yet the net change in average membrane thickness remains below 1%. Through our multiscale SAXS analysis, we quantitatively extract structural parameters, ranging from proteoliposome size to internal membrane architecture, using constrained Bayesian inference, with priors derived from MD findings. Specifically, we determine the protein/lipid molar ratio and average membrane strain, revealing excellent agreement between experiment and simulation. In thinner bilayers, substantial protein loss limits the analysis, highlighting the role of bilayer stability in sample preparation. Moreover, the predominance of OmpLA monomers in the thicker membranes is consistent with weak, membrane-mediated repulsive interactions between protein inclusions. Collectively, this integrative approach establishes a framework for quantifying protein-lipid interactions across molecular and mesoscale dimensions.