Abstract
P4-ATPase translocates lipids from the exoplasmic to the cytosolic plasma membrane leaflet to maintain lipid asymmetry distribution in eukaryotic cells. P4-ATPase is associated with severe neurodegenerative and metabolic diseases such as neurological and motor disorders. Thus, it is important to understand its transport mechanism. However, even with progress in X-ray diffraction and cryo-electron microscopy techniques, it is difficult to obtain the dynamic information of the phospholipid transport process in detail. There are still some problems required to be resolved: (1) when does the lipid transport happen? (2) How do the key residues on the transmembrane helices contribute to the free energy of important states? In this work, we explore the phospholipid transport mechanism using a coarse-grained model and binding free energy calculations. We obtained the free energy landscape by coupling the protein conformational changes and the phospholipid transport event, taking ATP8A1-CDC50 (the typical subtype of P4-ATPase) as the research object. According to the results, we found that the phospholipid would bind to the ATP8A1-CDC50 at the early stage when ATP8A1-CDC50 changes from E2P to E2Pi-PL state. We also found that the electrostatic effects play crucial roles in the phospholipid transport process. The information obtained from this work could help us in designing novel drugs for P-type flippase disorders.