Computational Approach for Structural Feature Determination of Grapevine NHX Antiporters

Ayadi, Mariem; Ben Ayed, Rayda; Mzid, Rim; Aifa, Sami; Hanana, Mohsen

doi:https://doi.org/10.1155/2019/1031839

BioMed Research International

On this page

Abstract Introduction Methods Results Discussion Abbreviations Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 1031839 | https://doi.org/10.1155/2019/1031839

Computational Approach for Structural Feature Determination of Grapevine NHX Antiporters

Mariem Ayadi,^1,2Rayda Ben Ayed,²Rim Mzid,¹Sami Aifa,²and Mohsen Hanana¹

Guest Editor: Mostafa Abdelrahman

Received29 Aug 2018

Revised01 Nov 2018

Accepted08 Nov 2018

Published09 Jan 2019

Abstract

Plant NHX antiporters are responsible for monovalent cation/H⁺ exchange across cellular membranes and play therefore a critical role for cellular pH regulation, Na⁺ and K⁺ homeostasis, and salt tolerance. Six members of grapevine NHX family (VvNHX1-6) have been structurally characterized. Phylogenetic analysis revealed their organization in two groups: VvNHX1-5 belonging to group I (vacuolar) and VvNHX6 belonging to group II (endosomal). Conserved domain analysis of these VvNHXs indicates the presence of different kinds of domains. Out of these, two domains function as monovalent cation-proton antiporters and one as the aspartate-alanine exchange; the remaining are not yet with defined function. Overall, VvNHXs proteins are typically made of 11-13 putative transmembrane regions at their N-terminus which contain the consensus amiloride-binding domain in the TM domain and a cation-binding site in between the and TM domain, followed by a hydrophilic C-terminus that is the target of several and diverse regulatory posttranslational modifications. Using a combination of primary structure analysis, secondary structure alignments, and the tertiary structural models, the VvNHXs revealed mainly 18 α helices although without β sheets. Homology modeling of the 3D structure showed that VvNHX antiporters are similar to the bacterial sodium proton antiporters MjNhaP1 (Methanocaldococcus jannaschii) and PaNhaP (Pyrococcus abyssi).

1. Introduction

The NHX antiporters family mediates cation/proton exchange across different cellular membranes, using electrochemical gradients generated by proton translocating enzymes, namely, H⁺-ATPase in the plasma membrane and vacuolar ATPase and pyrophosphatase within the intracellular compartments [1, 2]. Members of this NHX family are found in all living organisms and cells, from prokaryotic to eukaryotic. In plants, they are thought to play function in several biological roles through their cation/proton exchange activity, such as cellular pH regulation, potassium homeostasis, growth, developmental processes, and osmotic stress tolerance [3–5]. In grapevine Vitis vinifera L., one member (VvNHX1) of the VvNHX family has been already characterized and reported by Hanana et al. [5]; the remaining members still to be described and studied. Although the number of Vitis vinifera protein sequences is growing considerably since the inception of genome sequencing, however, functional annotation of these proteins is often less accurate or missing. Therefore, functional studies and biological roles identification are clearly needed. Protein structure determination, whether performed through X-ray crystallography and NMR spectroscopy or computational and in silico methods, would improve our knowledge about their functions and biological roles. However, protein crystallization, mainly transmembranous, is a staggering task due to the problem of solubility and stability and is therefore time and cost consuming. Computational methods and bioinformatics tools would be an interesting alternative for protein structure prediction. Indeed, several computer tools and programs for structure prediction have been developed and are now available [6]. In this context, our work focuses on structural features determination of grapevine NHX antiporters via the use of computational and in silico methods. The objective of our study is to predict the three-dimensional structure of Vitis vinifera NHX antiporter protein family through homology modeling and examine its physicochemical properties using in silico approaches. Biocomputational analyses of the targeted proteins were performed using an array of online bioinformatics tools and databases. The homology model was developed using different software packages and the best model was selected upon evaluation.

2. Material and Methods

2.1. Database Search and Identification of Grapevine VvNHX Genes

To search for putative NHX genes in the genome of Vitis vinifera, a keyword search “sodium/hydrogen antiporter” was performed in National Center for Biotechnology Information (NCBI) database. The selected genes were used for a basic local alignment search tool (BLAST) on the grapevine genome in the Phytozome database V12 (http://www.phytozome.net/). To validate the reliability of these potential candidates, the Pfam database (http://pfam.sanger.ac.uk/search) was exploited to confirm each VvNHX candidate as a true member of NHX family.

2.2. Phylogenetic Relationships and Identification of Conserved Protein Motifs

Alignments were performed using the Clustal W module from MEGA6 software [7] which help to draw the unrooted phylogenetic tree [7] using the neighbor-joining method with statistical bootstrapping procedure involving 1000 replicates. MEME software (Multiple Expectation Maximization for Motif Elicitation) (MEME version 4.10.2) (http://meme.sdsc.edu/meme/intro.html) led us to the identification of conserved motifs and the description of their position and frequency [8]. The used parameters in the analysis were repetitions number: 3, maximum number of motifs: 10; minimum motif width: 10; and maximum motif width: 50. Structural motif annotation was used directly at the Conserved Domain Database (CDD) (http://www.ncbi.nlm.nih.gov/cdd/).

2.3. Structural Features Characterization

The analyzing tool of protein primary structures is provided by the Protpram ExPASy server (Expert Protein Analysis System; https://web.expasy.org/protparam/) [9] including MW (Molecular Weight), IEP (Isoelectric point), and composition of amino acids. The subcellular localization of VvNHX proteins was predicted with WOLF PSORT [10] (https://wolfpsort.hgc.jp/) and CELLO Prediction (http://cello.life.nctu.edu.tw/) [11]. Sequence similarity was calculated with the global alignment tool using BLOSUM50. Transmembrane (TM) helix prediction was performed by the consensus extracted from ARAMEMNON database which has been planned to collect various computational predictions for plant membrane proteins (Eiconda_v1, TmHMM, TMMOD, MemSAT_v3, Philius, Phobius, DAS-TMfilter, PredTmr-v1, Scampi, SosuiG_V1.1, THUBUP_v1, Hmmtop-V2, TmPred, and TopPred-v2) [12]. In particular, the alignment of these VvNHX proteins was drawn with PRALINE multiple sequence alignment using DSSP [13] and PSIPRED [14] to obtain secondary structure prediction. Also, the secondary structure of the protein was predicted by NPS Secondary Structure Prediction Method [15]. Structural disorder was examined using PONDR-FIT accessible from the platform of DisProt [16] using the VSL2B predictor option [17]. All curves obtained are overlapped by a tool of plot on the R graphic [18]. The 3D models of the targeted VvNHX proteins were constructed using protein structure homology model building programs PHYRE2 [19] (http://www.sbg.bio.ic.ac.uk/phyre2/html/page.cgi?id=index). The FASTA sequences of the query proteins were entered and the intensive mode was selected to attain 3D models. Molecular graphics and analyses were performed with the PyMOL Molecular Graphics System, Version 2.0 Schrödinger, LLC.

Posttranslational modification analyses were predicted using the CBS server (Center for Biological Sequence Analysis) (http://www.cbs.dtu.dk/services). Supposed phosphorylation sites were analyzed using KinasePhos (http://kinasephos.mbc.nctu.edu.tw) [20].

3. Results

3.1. Members of VvNHXs Family

6 members of VvNHX genes family were retrieved from the NCBI database in the whole grapevine genome after a procedure of several plant NHX sequences blast request and further screening for the presence of “sodium/hydrogen exchanger family” domain. Therefore, 6 nonredundant genes were found and confirmed as NHX genes in the grape genome, and then named VvNHX1 to VvNHX6 (Table 1). The prediction of their subcellular localization showed that most of VvNHXs belong to intracellular membranes.

3.2. Phylogenetic Groups and Motif Analysis of VvNHX Proteins

The phylogenetic tree (Figure 1) drawn using both Grapevine and Arabidopsis NHX protein sequences reveals the existence of two major groups of VvNHX antiporters (I and II) that gather either with AtNHX1–4 (group I) or AtNHX5/6 (group II). VvNHX6, presumed to be of endosomal location, is the only member of group II.

Figure 1

Conserved protein motifs and phylogeny of VvNHXs. Phylogeny and schematic representation of the 10 conserved motifs in NHX proteins of Grapevine and Arabidopsis. The tree was constructed with VvNHX and AtNHX proteins using the neighbor-joining algorithm of MEGA 6 [7]. The bootstrap consensus tree was inferred from 1000 replicates. Scale bar corresponds to 0.1 amino acid substitution per residue. Motif analysis was performed using MEME 4.0 software as described in Methods. Different motifs, numbered 1–10, are displayed in different colored boxes. A detailed motif introduction for NHXs proteins is shown in Table 2.

Furthermore, in order to identify conserved motifs and consensus domains constituting the NHX proteins, the online MEME Suite (v4.8.2) program was used (Figure 1). The sequence details of each motif are shown in Table 2. Analogous motifs were shared between the two different NHX groups, suggesting common conserved functions inside the NHX family. The first motif corresponding to aspartate-alanine exchange (Asp-Al_Ex) pattern is found in all groups between TM8 and TM9 regions. The conserved motifs 2 and 4 were characterized as Monovalent Cation: Proton Antiporter-1 (CPA1), which were broadly distributed in all NHX protein sequences. The CPA1 family members allow the electroneutral exchange of monovalent cations for H⁺ in a pH-dependent way [23]. Motifs 6 and 9, present in all VvNHXs, are conserved sequences in sodium/proton exchangers. The motifs 3, 5, 7, and 8 were found only within VvNHXs belonging to group I, while the motif 10 was distributed exclusively in VvNHX6 and in twice, though the function of these motifs are still unknown.

3.3. Structural Characterization of VvNHX Proteins

3.3.1. Primary Structure

The size of encoded proteins of the identified VvNHX genes ranged from 524 to 541 amino acids (aa), with an average of 535 aa, and the corresponding predicted molecular masses ranged from 58.4 to 60.1 kDa, with an average of 59.4 kDa (Table 1). The computed IEP of these proteins was ~7.13 on average, indicating that they are likely to precipitate in either acidic or basic buffers and can be maintained within a neutral buffer.

The VvNHX1-6 proteins shared 30 to 80% of similarity (Table 3). VvNHX1-5 (group I) in particular shared more than 80% of identity among VvNHX1 and 2 and more than 70% of identity among VvNHX1 and 4. However, VvNHX6 (group II) shares the lowest similarity (30%) with VvNHX3 and 5, although displaying similar lengths. This wide range of variability among VvNHX members would be synonymous of functional diversification.

Main represented aa (Table 4) of these NHX members are leucine (13%), alanine (9%), glycine (9%), valine (9%), and serine (8%). The least common aa residues were cysteine and tryptophan which accounted for ~1% of the protein’s primary structure. The low amounts of cysteine residues indicated that the probabilities of disulfide bond formation were low. Leucine, alanine and valine are hydrophobic, aliphatic and nonpolar amino acids and thus expected to be inside the protein or within lipid membranes. Glycine, the smallest aa, known for its flexibility, would allow the protein twist and rotation. Serine, a polar aa, is qualified as hydrogen bond donor. The contribution of most readily oxidable amino acids was only 10% in all VvNHXs, which could be gainful to get better stability of targeted proteins under stress. These assumed that the targeted vacuolar-NHX protein of Vitis vinifera is hydrophobic in nature. Negatively and positively charged aa number of VvNHXs ranged from 30 to 45 (Table 4). However, the residual and global charge (calculated at pH 7) of the whole NHX protein varies from negative (VvNHX2, 4, and 5), neutral (VvNHX1), positive (VvNHX3) to highly positive for VvNHX6 (Table 4), depending on the pH environment of the different extra-membrane VvNHX segments.

3.3.2. Secondary Structure

Data extracted from ARAMEMNON database [12] indicated that the amino acid sequences of VvNHXs comprise a predicted secretion pathway. Topological prediction analyses of transmembrane (TM) domains of VvNHXs showed a range from 11 to 13 TM domains, 23 residues for each (Figure 2(a); supplementary file 1). Our MEME analysis, set with the mentioned parameters, led to the identification and localization of 10 motifs within the VvNHX proteins (Figure 1). Motifs 2 and 4 are the functional CPA1 conserved domain of NHE proteins and are widely distributed in all VvNHXs proteins. All the members of VvNHXs contained a conserved sequence FFI/LY/FLLPPI in the motif 2, at a conserved position of the TM domain (Supplementary file 1), corresponding to an amiloride-binding site in mammals that is known to be an inhibitor of eukaryotic NHXs [24]. Harris and Fliegel [25] have demonstrated that sodium ions and amiloride molecule interact at unique regions of NHE-like Na⁺/H⁺ transporters. Cation-binding site was predicted to be located in between the and TM domain. The NHX TM domains are well conserved in different sequences; however, variations are seen at N- and C-termini of amino acids sequences (Figure 2(a), supplementary file 1). In contrast, VvNHX1-5 proteins (Group I) have a calmodulin (CaM) binding domain at the end of their C-termini, like the one identified in AtNHX1 [26]. This same domain was identified in VvNHX1 by Hanana et al. [5] from 489 to 524 amino acids. Helical wheel representation (Figure 2(b)) of these peptide regions displayed an amphipathic pattern. However, VvNHX6 does not share a high similarity of this region in its C-terminus, probably because of its endosomal location and its implication in other functions.

(a)

(b)

(c)

The two dimensional secondary structure of NHX protein sequences was predicted using GORIV of NPS [15]. Secondary structure prediction showed that all VvNHXs contain around 30% α-helices, 25% extended strands, and 45% random coils and loops (Table 5). In light of the above, the presence of helices in the protein makes it more flexible for folding, which might increase protein interactions [27]. A secondary factor that influences proteins’ stability is the presence of extensive hydrogen bonds.

To further analyze their structural flexibility, we carried out sequence analysis with PONDR-FIT, a computational tool for prediction of structured and unstructured regions [16]. As shown in Figure 2(c), several regions of VvNHXs are disordered. Unexpectedly, this disorder was roughly preserved all along the VvNHXs. These unstructured regions are registered along more than two-thirds of each sequence. Surprisingly, unstructured regions were observed especially in N- and C-termini, showing that the protein conservation within the NHX family is not only restricted to their primary sequence but also extended to their structural flexibility.

3.3.3. Posttranslational Modifications

Posttranslational modifications (PTMs) can adjust localization and protein activity including interactions between proteins and intrinsically disordered regions, hence making plants able to face environmental changes and constraints [28]. These modifications comprise acetylation, phosphorylation, glycosylation, sumoylation, methylation and many other types. The family of VvNHXs harbors multiple sites of activity regulation and PTMs illustrated in Table 6. Glycosylation PTM site is commonly known for biosynthetic processing of transporter proteins in yeast like the NHX1 of Saccharomyces cerevisiae which is described as a glycoprotein [29]. In average, two putative N-glycosylation sites were identified in conserved positions of the 6 sequences of VvNHXs; yet supplementary sites of N-glycosylation have been predicted elsewhere for VvNHX1, VvNHX2 and VvNHX3.

Phosphorylation, one of the most common PTMs of proteins [30] has been shown to regulate the activity of animal NHE antiporters [31]. In each member of VvNHXs, putative sites of phosphorylation were found and which ranged from 10 to 18, most of them were of serine type (Table 6; Figure 2(a)); however, only few sites of tyrosine type were detected.

Moreover, various regulatory sites in VvNHX proteins were detected including sumoylation which is a PTM associated with some cellular processes, such as protein stability, transcriptional regulation, response to stress, nuclear-cytosolic transport, and apoptosis.

3.3.4. Tertiary Structure

Protein Homology/analogY Recognition Engine (PHYRE 2) was used for the prediction of the 3D structure of VvNHX proteins. This tool develops template-based modeling in which a protein is aligned to another of known structure on the basis of patterns of evolutionary variation [19]. Based on homology search, all the VvNHX proteins show homology with the sodium proton antiporter MjNhaP1 from Methanocaldococcus jannaschii (c4czbB) with expectation (E) value of ≤ 1 x (Figure 3, Supplementary file 2). An additional homology with the sodium proton antiporter PaNhaP from Pyrococcus abyssi (c4cz8A) was also reported by PHYRE 2. The alignment of VvNHXs with the references sequences of MjNhaP1 (c4czbB) and PaNhaP (c4cz8A) was drawn with PRALINE multiple sequence alignment using DSSP and PSIPRED to perform secondary structure prediction (Figure 3). The selected eight sequences shared similar α-helical and extended strands/β-sheet content (Figure 3; Table 5). The analysis revealed that the α-helices were dominant among the secondary structures (18 helices) except VvNHX2, 4 and 5 which showed only a β-sheet structure.

Among the two structural analogs generated by PHYRE 2, c4czbB was chosen to be the appropriate template for homology modeling of the VvNHX antiporter proteins. Wöhlert and collaborators [32] have validated its structure in two different conformations at pH 8 and pH 4 (4CZ9_A and 4CZ8_A), also with bound thallium ion (4CZA_A). Although, no X-ray crystallographic structures for intracellular plant NHX antiporters are available, the intracellular NHXs mediate the H⁺-coupled transport of both Na⁺ and K⁺ [4, 33] while NhaA does not mediate K⁺ transport; meaning that the structure of NHXs must accommodate the relatively larger K⁺ ion [34].

Using the PyMol software, a comparison was made between the template and all VvNHXs (Figure 4). However, the beta strands existing in the Escherichia coli NhaA appear to be absent in the VvNHX structure. All the models of VvNHXs consisted mainly of alpha helices with high confidence levels followed by the random coil. As shown in Figure 4(a), the predicted 12 TM domains are localized mainly within alpha helices of VvNHXs except the model of VvNHX3 that harbors only 11 TM domains. Based on the longest sequence in alpha helix, the VvNHX2 protein sequence was selected to align the other sequences of NHX. Figure 4(b) shows that the TM domains are conserved along all the VvNHX; but C-terminal structures are different from each other remarkably and especially for VvNHX4, 5 and 6.

(a)

(b)

Figure 4

Ribbon representation of the predicted model of Vitis vinifera NHX exchangers. (a) α-helices and coils are colored from N- to C- terminal by red, blue, orange, yellow, purple, and cyan for, respectively, VvNHX1, VvNHX2, VvNHX3, VvNHX4, VvNHX5, and VvNHX6. The TM domains are colored in green and numbered from 1 to 12. (b) The aligned models are in the order: VvNHX2/1 (blue/ red), VvNHX2/3 (blue/ orange), VvNHX2/4 (blue/ yellow), VvNHX2/5 (blue/ purple), VvNHX2/6 (blue/ cyan), and all VvNHX studied. All TM domains are colored green. Molecular graphics were performed with the PyMOL Molecular Graphics System.

4. Discussion

Our work led to the identification of 6 NHX members in grapevine Vitis vinifera L. Based on their subcellular localization, as for Arabidopsis, VvNHXs formed two distinct phylogenetic groups: group I vacuolar (VvNHX1-5) and group II endosomal (VvNHX6).

Like other known intracellular antiporters, the deduced amino acid sequences of VvNHXs family reveal 11-13 putative TM regions, with a consensus of 12 TM.

Though their structures (secondary, tertiary, and quaternary) have yet to be elucidated, the first bacterial crystal structure of Escherichia coli Na⁺/H⁺ antiporter NhaA has been established [35, 36]. Recently, Wöhlert et al. [32] have crystallized a sodium/proton antiporter from a single-celled organism Pyrococcus abyssi called PaNhaP. They determined the substrate ion in the dimeric PaNhaP at 3.2 Å and have resolved its structure in two different conformations at pH 8 and pH 4. However, to date, no X-ray crystallographic structures for plant NHX antiporters are presented. Nonetheless, two separate topology models have been suggested for AtNHX1. The first model contained nine TM domains with a C-terminal hydrophilic in front of the vacuolar lumen [37] whereas the second topology model suggests that several TM domains of AtNHX1 preserve similar topogenic properties like human NHE1 [38].

Despite the quick progress in molecular identification and biochemical characterization of NHX antiporters, the systematic and the structural analyses of plant NHX antiporter genes have not yet been vigorously investigated. Afore mentioned, the Na⁺/H⁺ antiporter has a vital role in the exchange of Na⁺ for H⁺ across membranes, which is essential for the plant's salt tolerance due to the fact that it keeps cellular ion homeostasis.

Our framework provides new hypotheses and insights into some potential structural features of NHXs from grapevine and offers a basis for further characterizations of the 3D structure of VvNHXs.

Our study pointed out that cation-binding sites and an amiloride-binding site were conserved in all VvNHX proteins. Cation-binding sites shared a high similarity within AtNHX1 homologous. The AtNHX1 has been applied to hold the predicted cation-binding domain as mapped by analogy to other antiporters [39]. Particularly, VvNHX6 displayed the lowest predicted IEP (5.4) which could be linked to its endosomal localization; indeed, according to some authors [40–42], subcellular localization of proteins would be correlated to their IEP values, according to the intracellular pH conditions. In plant cells, there is a decreasing gradient of pH from nuclear to vacuolar compartments, meaning that endosomal compartments (ER, TGN and PVC) are closest to pH 7, matching therefore with the low value of VvNHX6 IEP. Indeed, such neutral pH would be more suitable for VvNHX6 stability and activity; whereas, the acidic vacuolar pH (~ pH 5) is more appropriate to VvNHX1-5 according to their higher IEP. The VvNHX6, like MaNHX6 which belongs to the group II of the plant NHX family [43], does not have the region of CaM-Binding domain and encodes a putative K⁺/H⁺ exchanger. Nevertheless, nor the structural nor the regulatory mechanisms have been unraveled in the plant group II NHX isoforms.

The family of VvNHX contains several sites of activity regulation and PTMs including phosphorylation and glycosylation and shares the same interacting partners. Phosphoproteomic studies of Arabidopsis, barley and rice tonoplast preparations, further suggested the regulation of NHX antiporters by phosphorylation [44]. Knowing that the S471 conserved residue at C-terminus of the rice vacuolar-NHX isoforms was phosphorylated in OsNHX3, this could be confirmed for the grape VvNHX1 and VvNHX2 which have the same type of a conserved residue at serine 462.

Herein, the 3D structure of VvNHXs has been modeled in order to provide insights on their protein architecture and transport function. These antiporters share similarity with crystallized bacterial sodium proton antiporter, despite overall low sequence identity. This hypothesis is based on the observation that cotransporter structures are significantly more conserved in species than linear sequences [45]. However, Sato and Sakaguchi [38] have demonstrated that AtNHX1 has the same topology as human NHE1.

Using a combination of primary structure analysis, secondary structure alignments, and the tertiary structural models, VvNHX revealed conserved 12 TM located within 18 alpha helices. These TM compile a hollow cylinder and embedded in the membrane to provide the antiporter for Na⁺ and H⁺ transport. Recent researches of bacterial NhaA homologs have revealed a lack of β sheet, raising interesting questions about the role of the β sheet and the importance of the dimerization in these antiporters [46]. This helical phenomenon has been recommended and seen before in other proteins [47] and might be attributed to a particular function. Results support the idea that the core of all VvNHX resembles a sodium proton antiporter. These cation/H⁺ antiporters are essential for diverse functions, such as tolerance to salt stress, K⁺ uptake into vacuoles, pH regulation, and protein targeting [1, 48–51]. Furthermore, it was observed that 3D structures of VvNHX5 (Group I) and VvNHX6 (Group II) were the most divergent especially at their C-terminal structures. In fact, Yamaguchi and coworkers [37] solved the topology of AtNHX1 in which the C-terminus appears to be involved in the determination of the ion selectivity of the transporter. Thus 3D model structures can offer identify core residues of transporters to conclude function.

Abbreviations

CPA1:	Cation: Proton Antiporter-1
MEME:	Multiple Expectation Maximization for Motif Elicitation
TM:	Transmembrane
PTM:	Posttranslational Modification.

Data Availability

The accession numbers of VvNHX sequences data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Mariem Ayadi and Mohsen Hanana designed the study. Mariem Ayadi, Rayda Ben Ayed, Sami Aifa, and Rim Mzid performed data analyses. Mariem Ayadi and Mohsen Hanana analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgments

This research work was financially supported by the Ministry of Higher Education and Scientific Research of Tunisia.

Supplementary Materials

Supplementary File 1: alignment global of VvNHX proteins (.png). Supplementary File 2: details of tertiary structure of VvNHXs proteins (.doc). (Supplementary Materials)

References

M. P. Rodríguez-Rosales, F. J. Gálvez, R. Huertas et al., “Plant NHX cation/proton antiporters,” Plant Signaling and Behavior, vol. 4, no. 4, pp. 265–276, 2009.
View at: Publisher Site | Google Scholar
W. Dong, D. Li, N. Qiu, and Y. Song, “The functions of plant cation/proton antiporters,” Biologia Plantarum, vol. 62, no. 3, pp. 421–427, 2018.
View at: Publisher Site | Google Scholar
K. Venema, F. J. Quintero, J. M. Pardo, and J. P. Donaire, “The Arabidopsis Na+/H+ exchanger AtNHX1 catalyzes low affinity Na+ and K+ transport in reconstituted liposomes,” The Journal of Biological Chemistry, vol. 277, no. 4, pp. 2413–2418, 2002.
View at: Publisher Site | Google Scholar
M. P. Apse, J. B. Sottosanto, and E. Blumwald, “Vacuolar cation/H+ exchange, ion homeostasis, and leaf development are altered in a T-DNA insertional mutant of AtNHX1, the Arabidopsis vacuolar Na+/H+ antiporter,” The Plant Journal, vol. 36, no. 2, pp. 229–239, 2003.
View at: Publisher Site | Google Scholar
M. Hanana, O. Cagnac, T. Yamaguchi, S. Hamdi, A. Ghorbel, and E. Blumwald, “A grape berry (Vitis vinifera L.) cation/proton antiporter is associated with berry ripening,” Plant & Cell Physiology (PCP), vol. 48, no. 6, pp. 804–811, 2007.
View at: Publisher Site | Google Scholar
C. L. Mills, P. J. Beuning, and M. J. Ondrechen, “Biochemical functional predictions for protein structures of unknown or uncertain function,” Computational and Structural Biotechnology Journal, vol. 13, pp. 182–191, 2015.
View at: Publisher Site | Google Scholar
K. Tamura, G. Stecher, D. Peterson, A. Filipski, and S. Kumar, “MEGA6: Molecular Evolutionary Genetics Analysis version 6.0,” Molecular Biology and Evolution, vol. 30, no. 12, pp. 2725–2729, 2013.
View at: Publisher Site | Google Scholar
T. L. Bailey, M. Boden, F. A. Buske et al., “MEME SUITE: tools for motif discovery and searching,” Nucleic Acids Research, vol. 37, no. 2, pp. W202–W208, 2009.
View at: Publisher Site | Google Scholar
E. Gasteiger, A. Gattiker, C. Hoogland, I. Ivanyi, R. D. Appel, and A. Bairoch, “ExPASy: the proteomics server for in-depth protein knowledge and analysis,” Nucleic Acids Research, vol. 31, no. 13, pp. 3784–3788, 2003.
View at: Publisher Site | Google Scholar
P. Horton, K. Park, T. Obayashi et al., “WoLF PSORT: protein localization predictor,” Nucleic Acids Research, vol. 35, no. 2, pp. W585–W587, 2007.
View at: Publisher Site | Google Scholar
C.-S. Yu, Y.-C. Chen, C.-H. Lu, and J.-K. Hwang, “Prediction of protein subcellular localization,” Proteins: Structure, Function, and Bioinformatics, vol. 64, no. 3, pp. 643–651, 2006.
View at: Publisher Site | Google Scholar
R. Schwacke, A. Schneider, E. van der Graaff et al., “ARAMEMNON, a novel database for Arabidopsis integral membrane proteins,” Plant Physiology, vol. 131, no. 1, pp. 16–26, 2003.
View at: Publisher Site | Google Scholar
W. Kabsch and C. Sander, “Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features,” Biopolymers—Peptide Science Section, vol. 22, no. 12, pp. 2577–2637, 1983.
View at: Publisher Site | Google Scholar
D. T. Jones, “Protein secondary structure prediction based on position-specific scoring matrices,” Journal of Molecular Biology, vol. 292, no. 2, pp. 195–202, 1999.
View at: Publisher Site | Google Scholar
J. Garnier, J.-F. Gibrat, and B. Robson, “GOR method for predicting protein secondary structure from amino acid sequence,” Methods in Enzymology, vol. 266, pp. 540–553, 1996.
View at: Publisher Site | Google Scholar
B. Xue, R. L. Dunbrack, R. W. Williams, A. K. Dunker, and V. N. Uversky, “PONDR-FIT: a meta-predictor of intrinsically disordered amino acids,” Biochimica et Biophysica Acta (BBA) - Proteins and Proteomics, vol. 1804, no. 4, pp. 996–1010, 2010.
View at: Publisher Site | Google Scholar
Z. Obradovic, K. Peng, S. Vucetic, P. Radivojac, and A. K. Dunker, “Exploiting heterogeneous sequence properties improves prediction of protein disorder,” Proteins: Structure, Function, and Genetics, vol. 61, no. 7, pp. 176–182, 2005.
View at: Publisher Site | Google Scholar
“R: The R Project for Statistical Computing,” 2016, https://www.r-project.org/, [Accessed 10 Oct 2016].
View at: Google Scholar
L. A. Kelley, S. Mezulis, C. M. Yates, M. N. Wass, and M. J. E. Sternberg, “The Phyre2 web portal for protein modeling, prediction and analysis,” Nature Protocols, vol. 10, no. 6, pp. 845–858, 2015.
View at: Publisher Site | Google Scholar
H. Huang, T. Lee, S. Tzeng, and J. Horng, “KinasePhos: a web tool for identifying protein kinase-specific phosphorylation sites,” Nucleic Acids Research, vol. 33, no. 2, pp. W226–W229, 2005.
View at: Publisher Site | Google Scholar
U. Omasits, C. H. Ahrens, S. Müller, and B. Wollscheid, “Protter: Interactive protein feature visualization and integration with experimental proteomic data,” Bioinformatics, vol. 30, no. 6, pp. 884–886, 2014.
View at: Publisher Site | Google Scholar
R. Gautier, D. Douguet, B. Antonny, and G. Drin, “HELIQUEST: A web server to screen sequences with specific α-helical properties,” Bioinformatics, vol. 24, no. 18, pp. 2101-2102, 2008.
View at: Publisher Site | Google Scholar
Y. Ma, J. Wang, Y. Zhong, F. Geng, G. R. Cramer, and Z. Cheng, “Subfunctionalization of cation/proton antiporter 1 genes in grapevine in response to salt stress in different organs,” Horticulture Research, vol. 2, no. 1, 2015.
View at: Publisher Site | Google Scholar
L. Counillon, A. Franchi, and J. Pouyssegur, “A point mutation of the Na+/H+ exchanger gene (NHE1) and amplification of the mutated allele confer amiloride resistance upon chronic acidosis,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 90, no. 10, pp. 4508–4512, 1993.
View at: Publisher Site | Google Scholar
C. Harris and L. Fliegel, “Amiloride and the Na(+)/H(+) exchanger protein: mechanism and significance of inhibition of the Na(+)/H(+) exchanger (review).,” International Journal of Molecular Medicine, 1999.
View at: Publisher Site | Google Scholar
T. Yamaguchi, G. S. Aharon, J. B. Sottosanto, and E. Blumwald, “Vacuolar Na+/H+ antiporter cation selectivity is regulated by calmodulin from within the vacuole in a Ca2+- and pH-dependent manner,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 102, no. 44, pp. 16107–16112, 2005.
View at: Publisher Site | Google Scholar
A. R. Curran and D. M. Engelman, “Sequence motifs, polar interactions and conformational changes in helical membrane proteins,” Current Opinion in Structural Biology, vol. 13, no. 4, pp. 412–417, 2003.
View at: Publisher Site | Google Scholar
A. Hashiguchi and S. Komatsu, “Impact of Post-Translational Modifications of Crop Proteins under Abiotic Stress,” Proteomes, vol. 4, no. 4, p. 42, 2016.
View at: Publisher Site | Google Scholar
K. M. Wells and R. Rao, “The yeast Na+/H+ exchanger Nhx1 is an N-linked glycoprotein. Topological implications,” The Journal of Biological Chemistry, vol. 276, no. 5, pp. 3401–3407, 2001.
View at: Publisher Site | Google Scholar
F. J. M. Maathuis, “Conservation of protein phosphorylation sites within gene families and across species,” Plant Signaling and Behavior, vol. 3, no. 11, pp. 1011–1013, 2008.
View at: Publisher Site | Google Scholar
M. Vila-Petroff, C. Mundiña-Weilenmann, N. Lezcano et al., “Ca2+/calmodulin-dependent protein kinase II contributes to intracellular pH recovery from acidosis via Na+/H+ exchanger activation,” Journal of Molecular and Cellular Cardiology, vol. 49, no. 1, pp. 106–112, 2010.
View at: Publisher Site | Google Scholar
D. Wöhlert, W. Kühlbrandt, and Ö. Yildiz, “Structure and substrate ion binding in the sodium/proton antiporter PaNhaP,” eLife, vol. 3, 2014.
View at: Publisher Site | Google Scholar
H.-X. Zhang and E. Blumwald, “Transgenic salt-tolerant tomato plants accumulate salt in foliage but not in fruit,” Nature Biotechnology, vol. 19, no. 8, pp. 765–768, 2001.
View at: Publisher Site | Google Scholar
E. Bassil, A. Coku, and E. Blumwald, “Cellular ion homeostasis: Emerging roles of intracellular NHX Na +/H+ antiporters in plant growth and development,” Journal of Experimental Botany, vol. 63, no. 16, pp. 5727–5740, 2012.
View at: Publisher Site | Google Scholar
K. A. Williams, U. Geldmacher-Kaufer, E. Padan, S. Schuldiner, and W. Kühlbrandt, “Projection structure of NhaA, a secondary transporter from Escherichia coli, at 4.0 Å resolution,” EMBO Journal, vol. 18, no. 13, pp. 3558–3563, 1999.
View at: Publisher Site | Google Scholar
E. Padan, T. Tzubery, K. Herz, L. Kozachkov, A. Rimon, and L. Galili, “NhaA of Escherichia coli, as a model of a pH-regulated Na +/H+antiporter,” Biochimica et Biophysica Acta (BBA) - Bioenergetics, vol. 1658, no. 1-2, pp. 2–13, 2004.
View at: Publisher Site | Google Scholar
T. Yamaguchi, M. P. Apse, H. Shi, and E. Blumwald, “Topological analysis of a plant vacuolar Na+/H+ antiporter reveals a luminal C terminus that regulates antiporter cation selectivity,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 100, no. 21, pp. 12510–12515, 2003.
View at: Publisher Site | Google Scholar
Y. Sato and M. Sakaguchi, “Topogenic properties of transmembrane segments of Arabidopsis thaliana NHX1 reveal a common topology model of the Na+/H+ exchanger family,” The Journal of Biochemistry, vol. 138, no. 4, pp. 425–431, 2005.
View at: Publisher Site | Google Scholar
M. Hernandez, D. Ghersi, and R. Sanchez, “SITEHOUND-web: a server for ligand binding site identification in protein structures,” Nucleic Acids Research, vol. 37, no. 2, pp. W413–W416, 2009.
View at: Publisher Site | Google Scholar
W. A. Bickmore and H. G. E. Sutherland, “Addressing protein localization within the nucleus,” EMBO Journal, vol. 21, no. 6, pp. 1248–1254, 2002.
View at: Publisher Site | Google Scholar
R. Schwartz, C. S. Ting, and J. King, “Whole proteome pl values correlate with subcellular localizations of proteins for organisms within the three domains of life,” Genome Research, vol. 11, no. 5, pp. 703–709, 2001.
View at: Publisher Site | Google Scholar
S. Wu, P. Wan, J. Li, D. Li, Y. Zhu, and F. He, “Multi-modality of pI distribution in whole proteome,” Proteomics, vol. 6, no. 2, pp. 449–455, 2006.
View at: Publisher Site | Google Scholar
B. Cao, D. Long, M. Zhang, C. Liu, Z. Xiang, and A. Zhao, “Molecular characterization and expression analysis of the mulberry Na+/H+ exchanger gene family,” Plant Physiology and Biochemistry, vol. 99, pp. 49–58, 2016.
View at: Publisher Site | Google Scholar
S.-A. Whiteman, L. Serazetdinova, A. M. E. Jones et al., “Identification of novel proteins and phosphorylation sites in a tonoplast enriched membrane fraction of Arabidopsis thaliana,” Proteomics, vol. 8, no. 17, pp. 3536–3547, 2008.
View at: Publisher Site | Google Scholar
L. R. Forrest, R. Krämer, and C. Ziegler, “The structural basis of secondary active transport mechanisms,” Biochimica et Biophysica Acta (BBA) - Bioenergetics, vol. 1807, no. 2, pp. 167–188, 2011.
View at: Publisher Site | Google Scholar
A. Ganoth, R. Alhadeff, and I. T. Arkin, “Computational study of the Na +/H + antiporter from Vibrio parahaemolyticus,” Journal of Molecular Modeling, vol. 17, no. 8, pp. 1877–1890, 2011.
View at: Publisher Site | Google Scholar
H. Luecke, “Atomic resolution structures of bacteriorhodopsin photocycle intermediates: The role of discrete water molecules in the function of this light-driven ion pump,” Biochimica et Biophysica Acta (BBA) - Bioenergetics, vol. 1460, no. 1, pp. 133–156, 2000.
View at: Publisher Site | Google Scholar
E. Bassil, H. Tajima, Y.-C. Liang et al., “The arabidopsis Na+/H+ antiporters NHX1 and NHX2 control vacuolar ph and K+ homeostasis to regulate growth, flower development, and reproduction,” The Plant Cell, vol. 23, no. 9, pp. 3482–3497, 2011.
View at: Publisher Site | Google Scholar
E. Bassil, M.-A. Ohto, T. Esumi et al., “The Arabidopsis intracellular Na+/H+ antiporters NHX5 and NHX6 are endosome associated and necessary for plant growth and development,” The Plant Cell, vol. 23, no. 1, pp. 224–239, 2011.
View at: Publisher Site | Google Scholar
J. M. Pardo, B. Cubero, E. O. Leidi, and F. J. Quintero, “Alkali cation exchangers: Roles in cellular homeostasis and stress tolerance,” Journal of Experimental Botany, vol. 57, no. 5, pp. 1181–1199, 2006.
View at: Publisher Site | Google Scholar
M. Reguera, E. Bassil, and E. Blumwald, “Intracellular NHX-Type cation/H+ antiporters in plants,” Molecular Plant, vol. 7, no. 2, pp. 261–263, 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Mariem Ayadi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1251

Downloads

1256

Citations