Pentagonal Model of Membrane Channel of Na,K-ATPase. Computer Simulations of Structure and Electrostatic Properties

Abstract

Computational modeling of the membrane channel of a sodium pump (Na,K-ATPase) is performed and the role of selected amino acids in binding of sodium ions is discussed. The channel is build as a pentameric 10-helix bundle. The transmembrane a-helices are determined from hydropathy calculations. The spatial arrangement of transmembrane a-helices is chosen according to the size of a pore, intersegment loops geometry, and orientation hydrophobicities of transmembrane segments. The latter property provides the numerical estimate of the distribution of the hydrophobic properties at the helical wheels. The model system involves the peptide part and 150 water molecules that soak the pore. The channel structure is submitted to geometry minimization and molecular dynamics relaxation. The relative stability of the channel states with the negatively charged acidic residues belonging to the pore interior decrease in the order Glu-334 > Asp-810 > Glu-785 > Asp-814. The estimated binding energies of 1-3 Na+ ions with the channel with the ionized Glu-334 and Glu-785 amino acids are in the range allowing the exothermic complexation.