Evolution of cation binding in the active sites of P-loop nucleoside triphosphatases in relation to the basic catalytic mechanism
Abstract
The ubiquitous P-loop fold nucleoside triphosphatases (NTPases) are typically activated by an arginine or lysine 'finger'. Some of the apparently ancestral NTPases are, instead, activated by potassium ions. To clarify the activation mechanism, we combined comparative structure analysis with molecular dynamics (MD) simulations of Mg-ATP and Mg-GTP complexes in water and in the presence of potassium, sodium, or ammonium ions. In all analyzed structures of diverse P-loop NTPases, the conserved P-loop motif keeps the triphosphate chain of bound NTPs (or their analogs) in an extended, catalytically prone conformation, similar to that imposed on NTPs in water by potassium or ammonium ions. MD simulations of potassium-dependent GTPase MnmE showed that linking of alpha- and gamma phosphates by the activating potassium ion led to the rotation of the gamma-phosphate group yielding an almost eclipsed, catalytically productive conformation of the triphosphate chain, which could represent the basic mechanism of hydrolysis by P-loop NTPases.
Data availability
As obtained from MD simulations, we provide the structures of Mg-ATP complexes with bound K+, Na+ or NH4+ ions, as well as the structures of the G-domains of MnmE GTPases with and w/o activating potassium ions as source data files. Simulation data sets have been uploaded to Zenodo (https://zenodo.org/record/1888492#.XAasVhP7RTY).
-
Molecular Dynamics simulation data for the articleZenodo, 10.5281/zenodo.1888492.
Article and author information
Author details
Funding
Deutsche Forschungsgemeinschaft (N/A)
- Armen Y Mulkidjanian
Bundesministerium für Bildung und Forschung (Laseromix)
- Armen Y Mulkidjanian
Deutscher Akademischer Austauschdienst (Ostpartnerschaftenprogramm)
- Daria N Shalaeva
Russian Science Foundation (14-50-00029)
- Daria N Shalaeva
- Dmitry A Cherepanov
U.S. National Library of Medicine (Intramural Research Program)
- Michael Y Galperin
Lomonosov Moscow State University (RFMEFI62117X0011)
- Andrey V Golovin
Osnabrueck University, Germany (EvoCell Program)
- Armen Y Mulkidjanian
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Reviewing Editor
- Nick V Grishin, University of Texas Southwestern Medical Center, United States
Version history
- Received: April 9, 2018
- Accepted: November 26, 2018
- Accepted Manuscript published: December 11, 2018 (version 1)
- Version of Record published: December 28, 2018 (version 2)
Copyright
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Metrics
-
- 2,500
- Page views
-
- 303
- Downloads
-
- 27
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Evolutionary Biology
- Genetics and Genomics
Microbial plankton play a central role in marine biogeochemical cycles, but the timing in which abundant lineages diversified into ocean environments remains unclear. Here, we reconstructed the timeline in which major clades of bacteria and archaea colonized the ocean using a high-resolution benchmarked phylogenetic tree that allows for simultaneous and direct comparison of the ages of multiple divergent lineages. Our findings show that the diversification of the most prevalent marine clades spans throughout a period of 2.2 Ga, with most clades colonizing the ocean during the last 800 million years. The oldest clades – SAR202, SAR324, Ca. Marinimicrobia, and Marine Group II – diversified around the time of the Great Oxidation Event, during which oxygen concentration increased but remained at microaerophilic levels throughout the Mid-Proterozoic, consistent with the prevalence of some clades within these groups in oxygen minimum zones today. We found the diversification of the prevalent heterotrophic marine clades SAR11, SAR116, SAR92, SAR86, and Roseobacter as well as the Marine Group I to occur near to the Neoproterozoic Oxygenation Event (0.8–0.4 Ga). The diversification of these clades is concomitant with an overall increase of oxygen and nutrients in the ocean at this time, as well as the diversification of eukaryotic algae, consistent with the previous hypothesis that the diversification of heterotrophic bacteria is linked to the emergence of large eukaryotic phytoplankton. The youngest clades correspond to the widespread phototrophic clades Prochlorococcus, Synechococcus, and Crocosphaera, whose diversification happened after the Phanerozoic Oxidation Event (0.45–0.4 Ga), in which oxygen concentrations had already reached their modern levels in the atmosphere and the ocean. Our work clarifies the timing at which abundant lineages of bacteria and archaea colonized the ocean, thereby providing key insights into the evolutionary history of lineages that comprise the majority of prokaryotic biomass in the modern ocean.
-
- Evolutionary Biology
- Genetics and Genomics
In many species, meiotic recombination events tend to occur in narrow intervals of the genome, known as hotspots. In humans and mice, double strand break (DSB) hotspot locations are determined by the DNA-binding specificity of the zinc finger array of the PRDM9 protein, which is rapidly evolving at residues in contact with DNA. Previous models explained this rapid evolution in terms of the need to restore PRDM9 binding sites lost to gene conversion over time, under the assumption that more PRDM9 binding always leads to more DSBs. This assumption, however, does not align with current evidence. Recent experimental work indicates that PRDM9 binding on both homologs facilitates DSB repair, and that the absence of sufficient symmetric binding disrupts meiosis. We therefore consider an alternative hypothesis: that rapid PRDM9 evolution is driven by the need to restore symmetric binding because of its role in coupling DSB formation and efficient repair. To this end, we model the evolution of PRDM9 from first principles: from its binding dynamics to the population genetic processes that govern the evolution of the zinc finger array and its binding sites. We show that the loss of a small number of strong binding sites leads to the use of a greater number of weaker ones, resulting in a sharp reduction in symmetric binding and favoring new PRDM9 alleles that restore the use of a smaller set of strong binding sites. This decrease, in turn, drives rapid PRDM9 evolutionary turnover. Our results therefore suggest that the advantage of new PRDM9 alleles is in limiting the number of binding sites used effectively, rather than in increasing net PRDM9 binding. By extension, our model suggests that the evolutionary advantage of hotspots may have been to increase the efficiency of DSB repair and/or homolog pairing.