Coarse-grained molecular dynamics simulations of lipid-protein interactions in SLC4 proteins

The SLC4 family of secondary bicarbonate transporters is responsible for the transport of HCO ₃ -, CO ₃ ^2- , Cl ^- , Na ⁺ , K ⁺ , NH ₃ and H ⁺ necessary for regulation of pH and ion homeostasis. They are widely expressed in numerous tissues throughout the body and function in different cell types with different membrane properties. Potential lipid roles in SLC4 function have been reported in experimental studies, focusing mostly on two members of the family: AE1 (Cl ^- /HCO ₃ ^- exchanger) and NBCe1 (Na ⁺ -CO ₃ ^2- cotransporter). Previous computational studies of the outward facing (OF) state of AE1 with model lipid membranes revealed enhanced protein-lipid interactions between cholesterol (CHOL) and phosphatidylinositol bisphosphate (PIP2). However, the protein-lipid interactions in other members of the family and other conformation states are still poorly understood and this precludes the detailed studies of a potential regulatory role for lipids in the SLC4 family. In this work, we performed multiple 50 µs coarse-grained molecular dynamics simulations on three members of the SLC4 family with different transport modes: AE1, NBCe1 and NDCBE (a Na ⁺ -CO ₃ ^2- /Cl ^- exchanger), in model HEK293 membranes consisting of CHOL, PIP2, phosphatidylcholine (POPC), phosphatidylethanolamine (POPE), phosphatidylserine (POPS), and sphingomyelin (POSM). The recently resolved inward-facing (IF) state of AE1 was also included in the simulations. Lipid-protein contact analysis of the simulated trajectories was performed with the ProLint server, which provides a multitude of visualization tools for illustration of areas of enhanced lipid-protein contact and identification of putative lipid binding sites within the protein matrix. We observed enrichment of CHOL and PIP2 around all proteins with subtle differences in their distribution depending on the protein type and conformation state. Putative binding sites were identified for CHOL, PIP2, POPC, and POSM in the three studied proteins and their potential roles in the SLC4 transport function, conformational transition and protein dimerization were discussed.

The SLC4 protein family is involved in critical physiological processes like pH and blood pressure regulation and maintenance of ion homeostasis. Its members can be found in various tissues. A number of studies suggest possible lipid regulation of the SLC4 function. However, the protein-lipid interactions in the SLC4 family are still poorly understood. Here we make use of long coarse-grained molecular dynamics simulations to assess the protein-lipid interactions in three SLC4 proteins with different transport modes, AE1, NBCe1, and NDCBE. We identify putative lipid binding sites for several lipid types of potential mechanistic importance, discuss them in the framework of the known experimental data and provide a necessary basis for further studies on lipid regulation of SLC4 function.