Analysis of P-Loop and its Flanking Region Subsequence of Diverse NTPases Reveals Evolutionary Selected Residues

The P-loop NTPases are involved in diverse cellular functions. Members of the P-loop NTPase superfamily are characterized by presence of a highly conserved sequence pattern GxxxxGKS/T, known as Walker A motif. This motif adopts an archetypal P-loop conformation which allows accommodation of the triphosphate moiety of a bound nucleotide. Despite the presence of Walker A as a common sequence motif, P-loop NTPases exhibit extreme sequence divergence which hampers their phylogenetic or evolutionary classification. Here, we show that P-loop and its flanking region subsequence (termed as “extended-WalkerA motif”) contain distinct signatures that can be utilized to classify NTPase domain of functionally diverse proteins. We find a clearly classified group of diverse NTPases of Conserved Domain Database such as G-proteins, Ylqf, RecA like, DExDc, AAA, CPT, NK, ABC transporter and NifH proteins.     

Intradimeric Walker A ATPases: Conserved Features of A Functionally Diverse Family

Nucleoside-triphosphate hydrolases (NTPases) are a diverse, but essential group of enzymes found in all living organisms. NTPases that have a G-X-X-X-X-G-K-[S/T] consensus sequence (where X is any amino acid), known as the Walker A or P-loop motif, constitute a superfamily of P-loop NTPases. A subset of ATPases within this superfamily contains a modified Walker A motif, X-K-G-G-X-G-K-[S/T], wherein the first invariant lysine residue is essential to stimulate nucleotide hydrolysis. Although the proteins in this subset have vastly differing functions, ranging from electron transport during nitrogen fixation to targeting of integral membrane proteins to their correct membranes, they have evolved from a shared ancestor and have thus retained common structural features that affect their functions. These commonalities have only been disparately characterized in the context of their individual proteins systems, but have not been generally annotated as features that unite the members of this family. In this review, we report an analysis based on the sequences, structures, and functions of several members in this family that highlight their remarkable similarities. A principal feature of these proteins is their dependence on homodimerization. Since their functionalities are heavily influenced by changes that happen in conserved elements at the dimer interface, we refer to the members of this subclass as intradimeric Walker A ATPases.     

A novel family of P-loop NTPases with an unusual phyletic distribution and transmembrane segments inserted within the NTPase domain

Background

Recent sequence-structure studies on P-loop-fold NTPases have substantially advanced the existing understanding of their evolution and functional diversity. These studies provide a framework for characterization of novel lineages within this fold and prediction of their functional properties.

Results

Using sequence profile searches and homology-based structure prediction, we have identified a previously uncharacterized family of P-loop NTPases, which includes the neuronal membrane protein and receptor tyrosine kinase substrate Kidins220/ARMS, which is conserved in animals, the F-plasmid PifA protein involved in phage T7 exclusion, and several uncharacterized bacterial proteins. We refer to these (predicted) NTPases as the KAP family, after Kidins220/ARMS and PifA. The KAP family NTPases are sporadically distributed across a wide phylogenetic range in bacteria but among the eukaryotes are represented only in animals. Many of the prokaryotic KAP NTPases are encoded in plasmids and tend to undergo disruption to form pseudogenes. A unique feature of all eukaryotic and certain bacterial KAP NTPases is the presence of two or four transmembrane helices inserted into the P-loop NTPase domain. These transmembrane helices anchor KAP NTPases in the membrane such that the P-loop domain is located on the intracellular side. We show that the KAP family belongs to the same major division of the P-loop NTPase fold with the AAA+, ABC, RecA-like, VirD4-like, PilT-like, and AP/NACHT-like NTPase classes. In addition to the KAP family, we identified another small family of predicted bacterial NTPases, with two transmembrane helices inserted into the P-loop domain. This family is not specifically related to the KAP NTPases, suggesting independent acquisition of the transmembrane helices.

Conclusions

We predict that KAP family NTPases function principally in the NTP-dependent dynamics of protein complexes, especially those associated with the intracellular surface of cell membranes. Animal KAP NTPases, including Kidins220/ARMS, are likely to function as NTP-dependent regulators of the assembly of membrane-associated signaling complexes involved in neurite growth and development. One possible function of the prokaryotic KAP NTPases might be in the exclusion of selfish replicons, such as viruses, from the host cells. Phylogenetic analysis and phyletic patterns suggest that the common ancestor of the animals acquired a KAP NTPase via lateral transfer from bacteria. However, an earlier transfer into eukaryotes followed by multiple losses in several eukaryotic lineages cannot be ruled out.     

Classification and evolution of P-loop GTPases and related ATPases

Sequences and available structures were compared for all the widely distributed representatives of the P-loop GTPases and GTPase-related proteins with the aim of constructing an evolutionary classification for this superclass of proteins and reconstructing the principal events in their evolution. The GTPase superclass can be divided into two large classes, each of which has a unique set of sequence and structural signatures (synapomorphies). The first class, designated TRAFAC (after translation factors) includes enzymes involved in translation (initiation, elongation, and release factors), signal transduction (in particular, the extended Ras-like family), cell motility, and intracellular transport. The second class, designated SIMIBI (after signal recognition particle, MinD, and BioD), consists of signal recognition particle (SRP) GTPases, the assemblage of MinD-like ATPases, which are involved in protein localization, chromosome partitioning, and membrane transport, and a group of metabolic enzymes with kinase or related phosphate transferase activity. These two classes together contain over 20 distinct families that are further subdivided into 57 subfamilies (ancient lineages) on the basis of conserved sequence motifs, shared structural features, and domain architectures. Ten subfamilies show a universal phyletic distribution compatible with presence in the last universal common ancestor of the extant life forms (LUCA). These include four translation factors, two OBG-like GTPases, the YawG/YlqF-like GTPases (these two subfamilies also consist of predicted translation factors), the two signal-recognition-associated GTPases, and the MRP subfamily of MinD-like ATPases. The distribution of nucleotide specificity among the proteins of the GTPase superclass indicates that the common ancestor of the entire superclass was a GTPase and that a secondary switch to ATPase activity has occurred on several independent occasions during evolution. The functions of most GTPases that are traceable to LUCA are associated with translation. However, in contrast to other superclasses of P-loop NTPases (RecA-F1/F0, AAA+, helicases, ABC), GTPases do not participate in NTP-dependent nucleic acid unwinding and reorganizing activities. Hence, we hypothesize that the ancestral GTPase was an enzyme with a generic regulatory role in translation, with subsequent diversification resulting in acquisition of diverse functions in transport, protein trafficking, and signaling. In addition to the classification of previously known families of GTPases and related ATPases, we introduce several previously undetected families and describe new functional predictions.     

X-ray structure of HPr kinase: a bacterial protein kinase with a P-loop nucleotide-binding domain

HPr kinase/phosphatase (HprK/P) is a key regulatory enzyme controlling carbon metabolism in Gram- positive bacteria. It catalyses the ATP-dependent phosphorylation of Ser46 in HPr, a protein of the phosphotransferase system, and also its dephosphorylation. HprK/P is unrelated to eukaryotic protein kinases, but contains the Walker motif A characteristic of nucleotide-binding proteins. We report here the X-ray structure of an active fragment of Lactobacillus casei HprK/P at 2.8 A resolution, solved by the multiwavelength anomalous dispersion method on a seleniated protein (PDB code 1jb1). The protein is a hexamer, with each subunit containing an ATP-binding domain similar to nucleoside/nucleotide kinases, and a putative HPr-binding domain unrelated to the substrate-binding domains of other kinases. The Walker motif A forms a typical P-loop which binds inorganic phosphate in the crystal. We modelled ATP binding by comparison with adenylate kinase, and designed a tentative model of the complex with HPr based on a docking simulation. The results confirm that HprK/P represents a new family of protein kinases, first identified in bacteria, but which may also have members in eukaryotes.     

Nucleotide Binding States of Subunit A of the A-ATP Synthase and the Implication of P-Loop Switch in Evolution

The crystal structures of the nucleotide-empty (A_E), 5′-adenylyl-β,γ-imidodiphosphate (A_PNP)-bound, and ADP (A_DP)-bound forms of the catalytic A subunit of the energy producer A₁A_O ATP synthase from Pyrococcus horikoshii OT3 have been solved at 2.47 Å and 2.4 Å resolutions. The structures provide novel features of nucleotide binding and depict the residues involved in the catalysis of the A subunit. In the A_E form, the phosphate analog SO₄²⁻ binds, via a water molecule, to the phosphate binding loop (P-loop) residue Ser238, which is also involved in the phosphate binding of ADP and 5′-adenylyl-β,γ-imidodiphosphate. Together with amino acids Gly234 and Phe236, the serine residue stabilizes the arched P-loop conformation of subunit A, as shown by the 2.4-Å structure of the mutant protein S238A in which the P-loop flips into a relaxed state, comparable to the one in catalytic β subunits of F₁F_O ATP synthases. Superposition of the existing P-loop structures of ATPases emphasizes the unique P-loop in subunit A, which is also discussed in the light of an evolutionary P-loop switch in related A₁A_O ATP synthases, F₁F_O ATP synthases, and vacuolar ATPases and implicates diverse catalytic mechanisms inside these biological motors.     

Effect of ionic strength and oxidation on the P-loop conformation of the protein tyrosine phosphatase-like phytase, PhyAsr

The protein tyrosine phosphatase (PTP)-like phytase, PhyAsr, from Selenomonas ruminantium is a novel member of the PTP superfamily, and the only described member that hydrolyzes myo-inositol-1,2,3,4,5,6-hexakisphosphate. In addition to the unique substrate specificity of PhyAsr, the phosphate-binding loop (P-loop) has been reported to undergo a conformational change from an open (inactive) to a closed (active) conformation upon ligand binding at low ionic strength. At high ionic strengths, the P-loop was observed in the closed, active conformation in both the presence and absence of ligand. To test whether the P-loop movement can be induced by changes in ionic strength, we examined the effect that ionic strength has on the catalytic efficiency of PhyAsr, and determined the structure of the enzyme at several ionic strengths. The catalytic efficiency of PhyAsr is highly sensitive to ionic strength, with a seven-fold increase in k(cat)/K(m) and a ninefold decrease in K(m) when the ionic strength is increased from 100 to 500 mm. Surprisingly, the P-loop is observed in the catalytically competent conformation at all ionic strengths, despite the absence of a ligand. Here we provide structural evidence that the ionic strength dependence of PhyAsr and the conformational change in the P-loop are not linked. Furthermore, we demonstrate that the previously reported P-loop conformational change is a result of irreversible oxidation of the active site thiolate. Finally, we rationalize the observed P-loop conformational changes observed in all oxidized PTP structures.     

Changing nucleotide specificity of the DEAD-box helicase Hera abrogates communication between the Q-motif and the P-loop

DEAD-box proteins disrupt or remodel RNA and protein/RNA complexes at the expense of ATP. The catalytic core is composed of two flexibly connected RecA-like domains. The N-terminal domain contains most of the motifs involved in nucleotide binding and serves as a minimalistic model for helicase/nucleotide interactions. A single conserved glutamine in the so-called Q-motif has been suggested as a conformational sensor for the nucleotide state. To reprogram the Thermus thermophilus RNA helicase Hera for use of oxo-ATP instead of ATP and to investigate the sensor function of the Q-motif, we analyzed helicase activity of Hera Q28E. Crystal structures of the Hera N-terminal domain Q28E mutant (TthDEAD_Q28E) in apo- and ligand-bound forms show that Q28E does change specificity from adenine to 8-oxoadenine. However, significant structural changes accompany the Q28E mutation, which prevent the P-loop from adopting its catalytically active conformation and explain the lack of helicase activity of Hera_Q28E with either ATP or 8-oxo-ATP as energy sources. 8-Oxo-adenosine, 8-oxo-AMP, and 8-oxo-ADP weakly bind to TthDEAD_Q28E but in non-canonical modes. These results indicate that the Q-motif not only senses the nucleotide state of the helicase but could also stabilize a catalytically competent conformation of the P-loop and other helicase signature motifs.     

Fitness consequences of centrality in mutualistic individual-based networks

The relationships among the members of a population can be visualized using individual networks, where each individual is a node connected to each other by means of links describing the interactions. The centrality of a given node captures its importance within the network. We hypothesize that in mutualistic networks, the centrality of a node should benefit its fitness. We test this idea studying eight individual-based networks originated from the interaction between Erysimum mediohispanicum and its flower visitors. In these networks, each plant was considered a node and was connected to conspecifics sharing flower visitors. Centrality indicates how well connected is a given E. mediohispanicum individual with the rest of the co-occurring conspecifics because of sharing flower visitors. The centrality was estimated by three network metrics: betweenness, closeness and degree. The complex relationship between centrality, phenotype and fitness was explored by structural equation modelling. We found that the centrality of a plant was related to its fitness, with plants occupying central positions having higher fitness than those occupying peripheral positions. The structural equation models (SEMs) indicated that the centrality effect on fitness was not merely an effect of the abundance of visits and the species richness of visitors. Centrality has an effect even when simultaneously accounting for these predictors. The SEMs also indicated that the centrality effect on fitness was because of the specific phenotype of each plant, with attractive plants occupying central positions in networks, in relation to the distribution of conspecific phenotypes. This finding suggests that centrality, owing to its dependence on social interactions, may be an appropriate surrogate for the interacting phenotype of individuals.     

Switch II mutants reveal coupling between the nucleotide- and actin-binding regions in myosin V

Conserved active-site elements in myosins and other P-loop NTPases play critical roles in nucleotide binding and hydrolysis; however, the mechanisms of allosteric communication among these mechanoenzymes remain unresolved. In this work we introduced the E442A mutation, which abrogates a salt-bridge between switch I and switch II, and the G440A mutation, which abolishes a main-chain hydrogen bond associated with the interaction of switch II with the γ phosphate of ATP, into myosin V. We used fluorescence resonance energy transfer between mant-labeled nucleotides or IAEDANS-labeled actin and FlAsH-labeled myosin V to examine the conformation of the nucleotide- and actin-binding regions, respectively. We demonstrate that in the absence of actin, both the G440A and E442A mutants bind ATP with similar affinity and result in only minor alterations in the conformation of the nucleotide-binding pocket (NBP). In the presence of ADP and actin, both switch II mutants disrupt the formation of a closed NBP actomyosin.ADP state. The G440A mutant also prevents ATP-induced opening of the actin-binding cleft. Our results indicate that the switch II region is critical for stabilizing the closed NBP conformation in the presence of actin, and is essential for communication between the active site and actin-binding region.     

Biochemical and X-ray crystallographic studies on shikimate kinase: the important structural role of the P-loop lysine

Shikimate kinase, despite low sequence identity, has been shown to be structurally a member of the nucleoside monophosphate (NMP) kinase family, which includes adenylate kinase. In this paper we have explored the roles of residues in the P-loop of shikimate kinase, which forms the binding site for nucleotides and is one of the most conserved structural features in proteins. In common with many members of the P-loop family, shikimate kinase contains a cysteine residue 2 amino acids upstream of the essential lysine residue; the side chains of these residues are shown to form an ion pair. The C13S mutant of shikimate kinase was found to be enzymatically active, whereas the K15M mutant was inactive. However, the latter mutant had both increased thermostability and affinity for ATP when compared to the wild-type enzyme. The structure of the K15M mutant protein has been determined at 1.8 A, and shows that the organization of the P-loop and flanking regions is heavily disturbed. This indicates that, besides its role in catalysis, the P-loop lysine also has an important structural role. The structure of the K15M mutant also reveals that the formation of an additional arginine/aspartate ion pair is the most likely reason for its increased thermostability. From studies of ligand binding it appears that, like adenylate kinase, shikimate kinase binds substrates randomly and in a synergistic fashion, indicating that the two enzymes have similar catalytic mechanisms.     

Crystal structure of the nucleotide-binding domain of the ABC-transporter haemolysin B: identification of a variable region within ABC helical domains

The ABC-transporter haemolysin B is a central component of the secretion machinery that translocates the toxin, haemolysin A, in a Sec-independent fashion across both membranes of E. coli. Here, we report the X-ray crystal structure of the nucleotide-binding domain (NBD) of HlyB. The molecule shares the common overall architecture of ABC-transporter NBDs. However, the last three residues of the Walker A motif adopt a 3(10) helical conformation, stabilized by a bound anion. In consequence, this results in an unusual interaction between the Walker A lysine residue and the Walker B glutamate residue. As these residues are normally required to be available for ATP binding, for catalysis and for dimer formation of ABC domains, we suggest that this conformation may represent a latent monomeric form of the NBD. Surprisingly, comparison of available NBD structures revealed a structurally diverse region (SDR) of about 30 residues within the helical arm II domain, unique to each of the eight NBDs analyzed. As this region interacts with the transmembrane part of ABC-transporters, the SDR helps to explain the selectivity and/or targeting of different NBDs to their cognate transmembrane domains.     

P-loop Flexibility in Na+ Channel Pores Revealed by Single- and Double-cysteine Replacements

Replacement of individual P-loop residues with cysteines in rat skeletal muscle Na⁺ channels (SkM1) caused an increased sensitivity to current blockade by Cd²⁺ thus allowing detection of residues lining the pore. Simultaneous replacement of two residues in distinct P-loops created channels with enhanced and reduced sensitivity to Cd²⁺ block relative to the individual single mutants, suggesting coordinated Cd²⁺ binding and cross-linking by the inserted sulfhydryl pairs. Double-mutant channels with reduced sensitivity to Cd²⁺ block showed enhanced sensitivity after the application of sulfhydryl reducing agents. These results allow identification of residue pairs capable of approaching one another to within less than 3.5 Å. We often observed that multiple consecutive adjacent residues in one P-loop could coordinately bind Cd²⁺ with a single residue in another P-loop. These results suggest that, on the time-scale of Cd²⁺ binding to mutant Na⁺ channels, P-loops show a high degree of flexibility.     

A synthetic hexapeptide designed to resemble a proteinaceous P-loop nest is shown to bind inorganic phosphate

The hexapeptide Ser-Gly-Ala-Gly-Lys-Thr has been synthesized and characterized. It was designed as a minimal soluble peptide that would be likely to have the phosphate-binding properties observed in the P-loops of proteins that bind the β-phosphate of GTP or ATP. The β-phosphate in such proteins is bound by a combination of the side chain ε-amino group of the lysine residue plus the concavity formed by successive main chain peptide NH groups called a nest, which is favored by the glycines. The hexapeptide is shown to bind HPO(4) (2-) strongly at neutral pH. The affinities of the various ionized species of phosphate and hexapeptide are analyzed, showing that they increase with pH. It is likely the main chain NH groups of the hexapeptide bind phosphate in much the same way as the corresponding P-loop atoms bind the phosphate ligand in proteins. Most proteinaceous P-loops are situated at the N-termini of α-helices, and this observation has frequently been considered a key aspect of these binding sites. Such a hexapeptide in isolation seems unlikely to form an α-helix, an expectation in accord with the CD spectra examined; this suggests that being at the N-terminus of an α-helix is not essential for phosphate binding. An unexpected finding about the hexapeptide-HPO(4) (2-) complex is that the side chain ε-amino group of the lysine occurs in its deprotonated form, which appears to bind HPO(4) (2-) via an N···H-O-P hydrogen bond.     

Functional characterization of EngA(MS), a P-loop GTPase of Mycobacterium smegmatis

Bacterial P-loop GTPases belong to a family of proteins that selectively hydrolyze a small molecule guanosine tri-phosphate (GTP) to guanosine di-phosphate (GDP) and inorganic phosphate, and regulate several essential cellular activities such as cell division, chromosomal segregation and ribosomal assembly. A comparative genome sequence analysis of different mycobacterial species indicates the presence of multiple P-loop GTPases that exhibit highly conserved motifs. However, an exact function of most of these GTPases in mycobacteria remains elusive. In the present study we characterized the function of a P-loop GTPase in mycobacteria by employing an EngA homologue from Mycobacterium smegmatis, encoded by an open reading frame, designated as MSMEG_3738. Amino acid sequence alignment and phylogenetic analysis suggest that MSMEG_3738 (termed as EngA(MS)) is highly conserved in mycobacteria. Homology modeling of EngA(MS) reveals a cloverleaf structure comprising of α/β fold typical to EngA family of GTPases. Recombinant EngA(MS) purified from E. coli exhibits a GTP hydrolysis activity which is inhibited by the presence of GDP. Interestingly, the EngA(MS) protein is co-eluted with 16S and 23S ribosomal RNA during purification and exhibits association with 30S, 50S and 70S ribosomal subunits. Further studies demonstrate that GTP is essential for interaction of EngA(MS) with 50S subunit of ribosome and specifically C-terminal domains of EngA(MS) are required to facilitate this interaction. Moreover, EngA(MS) devoid of N-terminal region interacts well with 50S even in the absence of GTP, indicating a regulatory role of the N-terminal domain in EngA(MS)-50S interaction.     

The ATP Sites of AAA+ Clamp Loaders Work Together as a Switch to Assemble Clamps on DNA

Clamp loaders belong to a family of proteins known as ATPases associated with various cellular activities (AAA+). These proteins utilize the energy from ATP binding and hydrolysis to perform cellular functions. The clamp loader is required to load the clamp onto DNA for use by DNA polymerases to increase processivity. ATP binding and hydrolysis are coordinated by several key residues, including a conserved Lys located within the Walker A motif (or P-loop). This residue is required for each subunit to bind ATP. The specific function of each ATP molecule bound to the Saccharomyces cerevisiae clamp loader is unknown. A series of point mutants, each lacking a single Walker A Lys residue, was generated to study the effects of abolishing ATP binding in individual clamp loader subunits. A variety of biochemical assays were used to analyze the function of ATP binding during discrete steps of the clamp loading reaction. All mutants reduced clamp binding/opening to different degrees. Decreased clamp binding activity was generally correlated with decreases in the population of open clamps, suggesting that differences in the binding affinities of Walker A mutants stem from differences in stabilization of proliferating cell nuclear antigen in an open conformation. Walker A mutations had a smaller effect on DNA binding than clamp binding/opening. Our data do not support a model in which each ATP site functions independently to regulate a different step in the clamp loading cycle to coordinate these steps. Instead, the ATP sites work in unison to promote conformational changes in the clamp loader that drive clamp loading.     

Block of sodium channels by divalent mercury: role of specific cysteinyl residues in the P-loop region

Divalent mercury (Hg(2+)) blocked human skeletal Na(+) channels (hSkM1) in a stable dose-dependent manner (K(d) = 0.96 microM) in the absence of reducing agent. Dithiothreitol (DTT) significantly prevented Hg(2+) block of hSkM1, and Hg(2+) block was also readily reversed by DTT. Both thimerosal and 2,2'-dithiodipyridine had little effect on hSkM1; however, pretreatment with thimerosal attenuated Hg(2+) block of hSkM1. Y401C+E758C rat skeletal muscle Na(+) channels (mu1) that form a disulfide bond spontaneously between two cysteines at the 401 and 758 positions showed a significantly lower sensitivity to Hg(2+) (K(d) = 18 microM). However, Y401C+E758C mu1 after reduction with DTT had a significantly higher sensitivity to Hg(2+) (K(d) = 0.36 microM) than wild-type hSkM1. Mutants C753Amu1 (K(d) = 8.47 microM) or C1521A mu1 (K(d) = 8.63 microM) exhibited significantly lower sensitivity to Hg(2+) than did wild-type hSkM1, suggesting that these two conserved cysteinyl residues of the P-loop region may play an important role in the Hg(2+) block of the hSkM1 isoform. The heart Na(+) channel (hH1) was significantly more sensitive to low-dose Hg(2+) (K(d) = 0.43 microM) than was hSkM1. The C373Y hH1 mutant exhibited higher resistance (K(d) = 1.12 microM) to Hg(2+) than did wild-type hH1. In summary, Hg(2+) probably inhibits the muscle Na(+) channels at more than one cysteinyl residue in the Na(+) channel P-loop region. Hg(2+) exhibits a lower K(d) value (<1. 23 microM) for inhibition by forming a sulfur-Hg-sulfur bridge, as compared to reaction at a single cysteinyl residue with a higher K(d) value (>8.47 microM) by forming sulfur-Hg(+) covalently. The heart Na(+) channel isoform with more than two cysteinyl residues in the P-loop region exhibits an extremely high sensitivity (K(d) < 0. 43 microM) to Hg(+), accounting for heart-specific high sensitivity to the divalent mercury.