Computer modelling reveals new conformers of the ATP binding loop of Na+/K+-ATPase involved in the transphosphorylation process of the sodium pump

Department of Biophysics, 2nd Faculty of Medicine, Charles University Prague, Prague, Czech Republic
Laboratory of Tissue Engineering, Institute of Experimental Medicine, Academy of Sciences of the Czech Republic, Prague, Czech Republic
Department of Medical Chemistry and Clinical Biochemistry, 2nd Faculty of Medicine, Charles University Prague, Prague, Czech Republic
Small Animal Clinic, Faculty of Veterinary Medicine, University of Veterinary and Pharmaceutical Science, Brno, Czech Republic
Institute of Biochemistry and Endocrinology, University of Giessen, Giessen, Germany
DOI
10.7287/peerj.preprints.2812v1
Subject Areas
Biochemistry, Bioinformatics, Biophysics, Computational Biology, Computational Science
Keywords
M4M5 loop, Open and Closed conformations, Hinge movement., Na+/K+-ATPase phosphorylation
Copyright
© 2017 Tejral et al.
Licence
This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, reproduction and adaptation in any medium and for any purpose provided that it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ Preprints) and either DOI or URL of the article must be cited.
Cite this article
Tejral G, Sopko B, Necas A, Schoner W, Amler E. 2017. Computer modelling reveals new conformers of the ATP binding loop of Na+/K+-ATPase involved in the transphosphorylation process of the sodium pump. PeerJ Preprints 5:e2812v1

Abstract

Hydrolysis of ATP by Na+/K+-ATPase, a P-Type ATPase, catalyzing active Na+ and K+ transport through cellular membranes leads transiently to a phosphorylation of its catalytical α-subunit. Surprisingly, 3-dimensional molecular structure analysis of P-type ATPases reveals that binding of ATP to the N-domain connected by a hinge to the P-domain is much too far away from the Asp369 to allow the transfer of ATP’s terminal phosphate to its aspartyl-phosphorylation site. In order to get information how the transfer of the γ‑phosphate group of ATP to the Asp369 is achieved, analogous molecular modeling of the M4-M5 loop of ATPase was performed using the crystal data of Na+/K+-ATPase of different species. Analogous molecular modeling of the cytoplasmic loop between Thr338 and Ile760 of the α2-subunit of Na+/K+-ATPase and the analysis of distances between the ATP binding site and phosphorylation site revealed the existence of 2 ATP binding sites in the open conformation, the first one close to Phe475 in the N-domain, the other one close to Asp369 in the P-domain. However, binding of Mg2+•ATP to any of these sites in the “open conformation” may not lead to phosphorylation of Asp369. Additional conformations of the cytoplasmic loop were found wobbling between “open conformation” <==> “semi-open conformation <==> “closed conformation” in the absence of 2Mg2+•ATP. The cytoplasmic loop’s conformational change to the “semi-open conformation” -- characterized by a hydrogen bond between Arg543 and Asp611 -- triggers by binding of 2Mg2+•ATP to a single ATP site and conversion to the “closed conformation” the phosphorylation of Asp369 in the P-domain, and hence the start of Na+/K+-activated ATP hydrolysis.

Author Comment

This is a submission to PeerJ for review.

Supplemental Information

Figure S1

The multialignment of the chosen target (P50993, AT1A2_HUMAN) sequence and the two templates (3B8E, 3KDP) for open conformation was prepared by MODELLER program (salign module).

DOI: 10.7287/peerj.preprints.2812v1/supp-1

Figure S2

The multialignment of the chosen target (P50993, AT1A2_HUMAN) sequence and the three templates (3WGU, 3WGV, 4HQJ) for closed conformation was prepared by MODELLER program (salign module).

DOI: 10.7287/peerj.preprints.2812v1/supp-2

Figure S3

The Berendsen coupling method ( Berendsen, Postma et al. 1984 ) was employed for the temperature and pressure coupling of a system to reflect the reference temperature of 300K and the pressure of 1bar. The leap-frog integration with 104 steps was used for stabilization, with integration step of 1fs, corresponding to 10ps simulation time to reach the equilibrium of the rectangular box. This stabilized rectangular box was used for the main thirty simulations with 5x106 steps (2 fs single step), corresponding to 10ns for each stabilization using the same simulation parameters as for the box stabilization.

DOI: 10.7287/peerj.preprints.2812v1/supp-3

Figure S4

The Berendsen coupling method ( Berendsen, Postma et al. 1984 ) was employed for the temperature and pressure coupling of a system to reflect the reference temperature of 300K and the pressure of 1bar. The leap-frog integration with 104 steps was used for stabilization, with integration step of 1fs, corresponding to 10ps simulation time to reach the equilibrium of the rectangular box. This stabilized rectangular box was used for the main thirty simulations with 5x106 steps (2 fs single step), corresponding to 10ns for each stabilization using the same simulation parameters as for the box stabilization.

DOI: 10.7287/peerj.preprints.2812v1/supp-4

Figure S5

The Berendsen coupling method ( Berendsen, Postma et al. 1984 ) was employed for the temperature and pressure coupling of a system to reflect the reference temperature of 300K and the pressure of 1bar. The leap-frog integration with 104 steps was used for stabilization, with integration step of 1fs, corresponding to 10ps simulation time to reach the equilibrium of the rectangular box. This stabilized rectangular box was used for the main thirty simulations with 5x106 steps (2 fs single step), corresponding to 10ns for each stabilization using the same simulation parameters as for the box stabilization.

DOI: 10.7287/peerj.preprints.2812v1/supp-5

Figure S6

The Berendsen coupling method ( Berendsen, Postma et al. 1984 ) was employed for the temperature and pressure coupling of a system to reflect the reference temperature of 300K and the pressure of 1bar. The leap-frog integration with 104 steps was used for stabilization, with integration step of 1fs, corresponding to 10ps simulation time to reach the equilibrium of the rectangular box. This stabilized rectangular box was used for the main thirty simulations with 5x106 steps (2 fs single step), corresponding to 10ns for each stabilization using the same simulation parameters as for the box stabilization.

DOI: 10.7287/peerj.preprints.2812v1/supp-6