Standard binding free-energy calculation of glycophorin a dimer in non-isotropic media sheds light on the transmembrane alpha-helices association mechanism.

Abstract

The transmembrane (TM) protein domain recognition and association processes happening in lipid bilayers are pivotal for understanding membrane protein biogenesis. The most straightforward system to study the transmembrane protein binding mechanism is the glycophorin A (GpA) dimeric TM domain. The GpA is a glycoprotein ubiquitous to the human erythrocyte membrane, which forms a non-covalent dimer through the reversible association of its membrane-spanning domain sequences, creating a crossing angle between two alpha-helices. In our study, the membrane protein complex is primarily located in the phosphatidylcholine (POPC) lipid bilayer, surrounded by water. Applying a rigorous theoretical framework of the standard binding free-energy calculation that combines molecular dynamics and potential-of-mean-force calculations, the so-called “geometrical route”, we investigated this complex in non-isotropic media. Our analysis shows that at large separation distances, the helices are pushed together by the lipids, which triggers the interhelical interactions to occur. This observation is in agreement with the “two-stage” model of integral membrane protein folding. Additionally, we observed that the monomers, while reversibly associated, go through a meta-stable state, where the GpA alpha-helices form a different shape compared to the initial one. From the theoretical point of view, we suggest that for the membrane proteins in the non-isotropic environment, the standard binding free-energy estimation via the geometrical route should be analyzed from another perspective compared to those of protein-protein complexes in water media. We revised the separation PMF evaluation, assuming that the lipid bilayer forms a pseudo-cylindrical zone for the reversible separation of the dimer. The methodology employed herein successfully addresses the challenging transmembrane protein-protein binding problem while offering promising perspectives for the binding affinity characterization of different protein complexes in lipid bilayers.