Molecular dynamics (MD) simulations have been used to unmask details of specific interactions of anionic phospholipids with intersubunit binding sites on the surface of the bacterial potassium channel KcsA. Crystallographic data on a diacyl glycerol fragment at this site were used to model phosphatidylethanolamine (PE), or phosphatidylglycerol (PG), or phosphatidic acid (PA) at the intersubunit binding sites. Each of these models of a KcsA-lipid complex was embedded in phosphatidyl choline bilayer and explored in a 20 ns MD simulation. H-bond analysis revealed that in terms of lipid-protein interactions PA > PG ≫ PE and revealed how anionic lipids (PG and PA) bind to a site provided by two key arginine residues (R64 and R89) at the interface between adjacent subunits. A 27 ns simulation was performed in which KcsA (without any lipids initially modeled at the R64/R89 sites) was embedded in a PE/PG bilayer. There was a progressive specific increase over the course of the simulation in the number of H-bonds of PG with KcsA. Furthermore, two specific PG binding events at R64/R89 sites were observed. The phosphate oxygen atoms of bound PG formed H-bonds to the guanidinium group of R89, whereas the terminal glycerol H-bonded to R64. Overall, this study suggests that simulations can help identify and characterize sites for specific lipid interactions on a membrane protein surface.
- Anions, Bacterial Proteins, Binding Sites, Cell Membrane, Computer Simulation, Crystallography, X-Ray, Databases, Protein, Escherichia coli Proteins, Hydrogen Bonding, Ions, Lipid Bilayers, Lipids, Models, Biological, Models, Molecular, Molecular Conformation, Mutation, Oxygen, Phosphates, Phosphatidic Acids, Phosphatidylethanolamines, Phosphatidylglycerols, Phospholipids, Potassium Channels, Potassium Channels, Voltage-Gated, Pressure, Protein Folding, Streptomyces lividans, Temperature, Time Factors