Evolution of cation binding in the active sites of P-loop nucleoside triphosphatases in relation to the basic catalytic mechanism

  1. Daria N Shalaeva
  2. Dmitry A Cherepanov
  3. Michael Y Galperin
  4. Andrey V Golovin
  5. Armen Y Mulkidjanian  Is a corresponding author
  1. University of Osnabrück, Germany
  2. Lomonosov Moscow State University, Russia
  3. Russian Academy of Sciences, Russia
  4. National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, United States
8 figures, 3 tables and 2 additional files

Figures

Mg-NTP complexes and their binding in the active sites of P-loop NTPases.

Phosphate chains of NTP molecules and their analogs are colored by atoms: oxygen atoms in red, phosphorus in orange. The K+ ion is shown as a purple sphere, Na+ ion is shown as a blue sphere, Mg2+

https://doi.org/10.7554/eLife.37373.002
Figure 2 with 3 supplements
Binding of monovalent cations to the Mg-ATP in water.

The color scheme is as in Figure 1. (A) Superposition of the ATP phosphate chain conformations observed in the MD simulations in the presence of K+ ions (shown in purple); Na+ ions (shown in blue) …

https://doi.org/10.7554/eLife.37373.003
Figure 2—figure supplement 1
Radial distribution of cations in the proximity of each oxygen atom.

Radial distributions are shown for all atoms of the ATP phosphate chain. (A) Atom names are in accordance with the CHARMM naming scheme (Vanommeslaeghe et al., 2010) and the recent IUPAC …

https://doi.org/10.7554/eLife.37373.004
Figure 2—figure supplement 2
Properties of cation binding to the ATP as derived from MD simulations.

(A) Probability distribution functions for cations around the phosphate chain. We have plotted the number of atoms inside the area centered on phosphorus atoms of the ATP phosphate chain as a …

https://doi.org/10.7554/eLife.37373.005
Figure 2—figure supplement 3
Binding of monovalent cations to the Mg-GTP in water.

Distances to the AG and BG binding sites (RAG and RBG) were calculated as averages of the distances to the two corresponding oxygen atoms (see Figure 2 in the main text). The distances to the oxygen …

https://doi.org/10.7554/eLife.37373.006
Figure 3 with 1 supplement
Cation binding induces eclipsed conformation of the phosphate chain.

(A) Conformations of Mg-ATP complexes with one and two K+ ions bound as inferred from MD simulations; left structure, no K+ ion bound in the AG site; right structure, a K+ ion is bound in the AG …

https://doi.org/10.7554/eLife.37373.007
Figure 3—figure supplement 1
Coupling between cation binding in the AG site and rotation of γ-phosphate relative to α- and β-phosphates.

Data from MD simulations with restraints on the positions of K+ ions (see the text and Supplementary file 1C). The top graph shows free energy calculated from normalized probabilities of ATP …

https://doi.org/10.7554/eLife.37373.008
Figure 4 with 3 supplements
Dynamics of the phosphate chain of the Mg-ATP complex with and without monovalent cations.

Each left panel shows the PA-PG distance (upper trace) and the PB-O3B-PG angle (bottom trace) in the course of MD simulations. Thin gray lines show actual values measured from each frame of the MD …

https://doi.org/10.7554/eLife.37373.009
Figure 4—figure supplement 1
Coordination of the Mg22+ion by the oxygen atoms of the ATP phosphate chain during MD simulations.

Black vertical lines indicate borders between independent simulations, thick colored lines show moving averages of distances measured during MD simulations. Oxygen atoms are labeled as in Figure 1D. …

https://doi.org/10.7554/eLife.37373.010
Figure 4—figure supplement 2
Estimation of correlation times for the PA-PG distances.

A,B, Changes of the distance value upon MD simulations of βγ-coordinated Mg-ATP complexes with no additional monovalent cations (A) and with K+ ions (B) provided as examples. (C) Autocorrelation …

https://doi.org/10.7554/eLife.37373.011
Figure 4—figure supplement 3
Estimation of correlation times for the PB-O-PG angles.

A, B, Changes of the angle value upon MD simulations of βγ-coordinated Mg-ATP complexes with no additional monovalent cations (A) and with K+ ions (B) provided as examples. (C) Autocorrelation …

https://doi.org/10.7554/eLife.37373.012
Figure 5 with 2 supplements
Heat maps of the Mg-ATP phosphate chain conformations distribution characterized by the PA-PG distances (X-axis) and PB-O3B-PG angles (Y-axis).

Heat maps for systems with monovalent cations include only conformations of Mg-ATP complexes with at least one cation present within 4 Å radius. The color intensity is proportional to the …

https://doi.org/10.7554/eLife.37373.014
Figure 5—figure supplement 1
Heat maps of the Mg-GTP phosphate chain conformations distribution characterized by the PA-PG distances (X-axis) and PB-O3B-PG angles (Y-axis).

Heat maps for systems with monovalent cations include only conformations of Mg-GTP complexes with at least one cation present within 4 Å area, and with Mg2+ in βγ coordination. The color intensity …

https://doi.org/10.7554/eLife.37373.015
Figure 5—figure supplement 2
Phosphate chain shape of ATP and GTP analogs in the X-ray structures of P-loop NTPases.

PDB entries for structures of P-loop NTPases were extracted from InterPro database entry IPR027417 ‘P-loop containing nucleoside triphosphate hydrolase’ and filtered to contain only those X-ray …

https://doi.org/10.7554/eLife.37373.016
Figure 6 with 4 supplements
Location of positive charges around the phosphate chain of Mg-NTP complexes in solution and in protein structures.

The color scheme is as in Figure 1; dark blue spheres indicate positions of positively charged side-chain nitrogen atoms of Lys and Arg residues, P-loop regions are shown as cartoons in grey. (A) …

https://doi.org/10.7554/eLife.37373.017
Figure 6—figure supplement 1
Active sites of P-loop NTPases with established K+-dependent activity (see Supplementary file 1A for the full list and references).

Each of the proteins shown has both Asn residues that were shown to be associated with binding of monovalent cations in related proteins (Ash et al., 2012). Switch I, including the K-loop, and its …

https://doi.org/10.7554/eLife.37373.018
Figure 6—figure supplement 2
Activation of the MnmE GTPase upon dimerization.

(A) Inactive dimer of the full-length MnmE in the GTP-bound form (the structure (PDB: 3GEI) was resolved with non-hydrolyzable GTP analogs). The P-loop domain is shown in grey, the K-loop is not …

https://doi.org/10.7554/eLife.37373.019
Figure 6—figure supplement 3
Activation of the GTPase Era upon RNA binding.

(A) Inactive Era in the GDP-bound form [PDB: 3IEU] (Tu et al., 2009) in two projections. (B) Active Era in complex with nucleotides 1506–1542 of 16S rRNA and a non-hydrolyzable analog of GTP [PDB: …

https://doi.org/10.7554/eLife.37373.020
Figure 6—figure supplement 4
Positively charged moieties in the active site of RecA-like recombinases.

(A) Cation-dependent RadA recombinase from Methanococcus voltae [PDB: 2F1H] (Qian et al., 2006). (B) Cation-independent RecA recombinase from E. coli [PDB: 3CMX] (Chen et al., 2008). The protein …

https://doi.org/10.7554/eLife.37373.021
Figure 7 with 2 supplements
Molecular dynamics of MnmE GTPase.

(A) Superposition of the GTP-binding sites of the inactive, monomeric G-domain of MnmE (the 2GJ8W system, blue) and the active K+-bound dimer of G-domains (the 2GJ8K system, red); representative …

https://doi.org/10.7554/eLife.37373.024
Figure 7—figure supplement 1
Hydrogen bonds lengths during MD simulations of the G-domain of MnmE.

Distances between phosphate chain oxygen atoms and surrounding amino acid residues were measured in the course of 100-ns MD simulations. (A) inactive monomer without a full-fledged K-loop (the 3GEI …

https://doi.org/10.7554/eLife.37373.025
Figure 7—figure supplement 2
The distance between O2A and O3G in GTP bound to MnmE as inferred from MD simulations.

Plotted are data for the active dimer of G-domains with K+ ions bound (the 2GJ8K system, red and magenta for individual monomers) and inactive, monomeric G-domain of MnmE with the K+ ion replaced by …

https://doi.org/10.7554/eLife.37373.026
Effects of Na+ binding on the shape of phosphate chain in solution and in Na+-adapted P-loop NTPases.

The color scheme is as in Figure 1, except that Al and F atoms in the GDP-AlF4- complexes are colored grey and cyan, respectively. (A) Superposition of the K+-bound (solid structure) and Na+-bound …

https://doi.org/10.7554/eLife.37373.027

Tables

Table 1
Effects of monovalent cations on the shape of the triphosphate chain of the Mg-ATP complex in water, as inferred from the MD simulation data.
https://doi.org/10.7554/eLife.37373.013
Added cationConformation of the triphosphate chain of Mg-ATP*
βγ-coordinationβγ-coordination,
‘curled’ phosphate chain
αβγ-coordination
PA-PG distance, ÅPB-O3B-PG anglePA-PG distance, ÅPB-O3B-PG anglePA-PG distance, ÅPB-O3B-PG angle
None5.46 ± 0.34122.3 ± 3.5N/A4.76 ± 0.18124.9 ± 3.3
K+4.91 ± 0.24122.0 ± 3.3N/A4.32 ± 0.24128.0 ± 3.5
Na+4.69 ± 0.22122.9 ± 3.24.60 ± 0.22124.0 ± 3.34.26 ± 0.37127.7 ± 3.6
NH4+4.85 ± 0.22122.3 ± 3.34.56 ± 0.21124.6 ± 3.34.22 ± 0.16127.8 ± 3.
  1. *The conformations of the Mg-ATP complex were determined as described in the text. Mean values and standard deviations of PA-PG distance (in Å) and the PB-O3B-PG angle (in degrees) were measured over the respective parts of the simulations. Simulation periods corresponding to βγ and αβγ conformations were identified by tracking distances between Mg2+ and non-bridging oxygen atoms of the phosphate chain (Figure 4—figure supplement 1); simulation periods corresponding to the ‘curled’ conformation were identified from PA-PG distance tracks and visual inspection of the phosphate chain shape (Figure 4). Data for the αβγ coordination of the Mg-ATP complex and conformations with curled phosphate chain were calculated from simulations 1–4 in Supplementary file 1C; characterization of the βγ-coordination was based on simulations 5–8 in Supplementary file 1C, see Supplementary file 1E for further details.

Table 2
Activation mechanisms within the classes of P-loop NTPases that contain both cation-dependent and cation-independent enzymes.
https://doi.org/10.7554/eLife.37373.022
SuperfamilyFamilyActivating chargeActivation
mechanism
Kinase-GTPase division, TRAFAC class
Classic translation factor GTPasesEF-G/EF-2K+Functional interaction with ribosomal RNA/other protein(s)/other domain(s) of the same protein
(Hwang and Inouye, 2001;Moreau et al., 2008; Tomar et al., 2011; Achila et al., 2012;Fasano et al., 1982; Ebel et al., 1992; Dubnoff and Maitra, 1972;Kuhle and Ficner, 2014; Manikas et al., 2016; Daigle and Brown, 2004; Foucher et al., 2012;Rafay et al., 2012; Pérez-Arellano et al., 2013; Villarroya et al., 2008)
EF-Tu/EF-1AK+
EIF2GK+
ERF3K+
IF-2K+
LepAK+
OBG-HflX-like GTPasesHflXK+
OBGK+
NOGK+
YchF/OLA1K+
YlqF/YawG GTPasesNOG2K+
RsgAK+
TrmE-Era-EngA-EngB-Septin-like GTPasesEngA (Der)K+
EngBK+
EraK+
FeoBK+Dimerization (e.g. mRNA-associated in the case of MnmE) (Chappie et al., 2010; Koenig et al., 2008; Gasper et al., 2009)
MnmEK+
SeptinArg finger
Toc34-likeArg finger
Dynamin-like GTPaseshGBPArg finger
DynaminK+/Na+
Extended RasRas familyArg fingerInteraction with a specialized activating protein or domain(Bos et al., 2007;
Cherfils and Zeghouf, 2013)
Gα subunitsArg finger
Myosin/kinesinMyosinArg finger
KinesinArg finger
ASCE division, RecA/F1-like class
DNA-repair and recombination ATPasesRecALys fingerDNA/RNA-dependent oligomerization(Chen et al., 2008)
RadAK+
Rho helicasesRhoArg fingerInteraction with the neighboring subunit within a conformationally coupled hexamer (Komoriya et al., 2012; Walker, 1998;Senior et al., 2002Skordalakes and Berger, 2006
T3SS ATPasesYscNArg finger
FlilArg finger
F-/V-type ATPasesV-type AArg finger
F-type β
V-type B
F-type α
Table 3
Monovalent cation binding in crystal structures of P-loop NTPases.
https://doi.org/10.7554/eLife.37373.023
ProteinPDB entryBound NTP analogOccupation of the AG sitePhosphate chain shape
CationDistance to the closest O atom of PA, Å*Distance to the closest O atom of PG, Å*,†PA-PG distance, Å*PB-O3B-PG angle, degrees
TRAFAC class NTPases
GTPase MnmE(TrmE)2gj8GDP AlF4-K+2.82.65.4136.3
2gjaGDP AlF4-NH4+2.92.55.4136.9
2gj9GDP AlF4-Rb+2.92.85.5131.6
GTPase FeoB3ss8GDP AlF4-K+2.82.65.4144.9
Dynamin-like proteins2x2eGDP AlF4-Na+4.02.55.3131.2
2x2fGDP AlF4-Na+4.12.65.3133.6
3w6pGDP AlF4-Na+42.45.5135.3
3t34GDP AlF4-Na+3.82.45.6149.3
GTPase Era3r9wGNPH2O33.45.1129.2
Eukaryotic translation initiation factor eIF5B4ncnGTPNa+2.42.45.0126.6
4tmvGSPNa+2.42.8 (S)§4.9126.3
4tmwGTPNa+2.42.44.9125.9
4tmzGSPK+2.73.3 (S)§4.9122.1
RecA/F1-like class NTPases
DNA recombinase RadA3ew9ANPK+6.23.35.1124.5
2f1hANPK+6.63.55.3125.3
2fpmANPK+5.92.65.1124.2
1xu4ANPK+6.12.75.2125.0
  1. *The values were measured directly in the respective protein structures displayed in PyMOL.

     If the γ-phosphate was replaced by an AlF4- complex, the distance was measured to the closest F atom

  2.  While GTPase Era has been shown to be K+-dependent (Rafay et al., 2012; Meier et al., 2000), the crystallization solution contained no K+, only Na+, so that the likely cation-binding site is occupied by a water molecule, which forms hydrogen bonds with K+ ligands.

    § Non-hydrolyzable GTP analog GDP-monothiophosphate (GSP) contains a sulfur atom in the place of the O1G atom of γ-phosphate; this atom in involved in coordination of monovalent cations in respective structures.

Additional files

Supplementary file 1

(A) Monovalent cation requirements of P-loop GTPases and ATPases. (B) Properties of monovalent cations and their interactions with the Mg2+-ATP complex. (C) Molecular dynamics simulations performed in this work. (D) Values of dihedral angles of the phosphate chains of Mg-ATP in the presence of K+ ions. (E) Lifetimes of the βγ-conformation of Mg-ATP complex during MD simulations. (F) Characteristics of the triphosphate chain for different interactions between the Mg2+ ion and ATP. (G) Comparison of the PA-PG distance measurements of the βγ-coordinated Mg-ATP complexes. (H) Comparison of the PA-PG distance measurements of the αβγ-coordinated Mg-ATP complexes. (I) Comparison of the PA-PG distance measurements for the αβγ-coordinated and ‘curled’ βγ-coordinated Mg-ATP complexes in different systems. (J) Comparison of the PB-O3B-PG angle measurements for the βγ-coordinated Mg-ATP complexes. (K) Comparison of the PB-O3B-PG angle measurements for the αβγ-coordinated Mg-ATP complexes. (L) Comparison of the PA-PG distance measurements for the αβγ-coordinated and ‘curled’ βγ-coordinated Mg-ATP complexes.

https://doi.org/10.7554/eLife.37373.028
Transparent reporting form
https://doi.org/10.7554/eLife.37373.029

Download links