In silico insights on diverse interacting partners and phosphorylation sites of respiratory burst oxidase homolog (Rbohs) gene families from Arabidopsis and rice
BMC Plant Biology volume 18, Article number: 161 (2018)
NADPH oxidase (Nox) is a critical enzyme involved in the generation of apoplastic superoxide (O2 −), a type of reactive oxygen species (ROS) and hence regulate a wide range of biological functions in many organisms. Plant Noxes are the homologs of the catalytic subunit from mammalian NADPH oxidases and are known as respiratory burst oxidase homologs (Rbohs). Previous studies have highlighted their versatile roles in tackling different kind of stresses and in plant growth and development. In the current study, potential interacting partners and phosphorylation sites were predicted for Rboh proteins from two model species (10 Rbohs from Arabidopsis thaliana and 9 from Oryza sativa japonica). The present work is the first step towards in silico prediction of interacting partners and phosphorylation sites for Rboh proteins from two plant species.
In this work, an extensive range of potential partners (unique and common), leading to diverse functions were revealed from interaction networks and gene ontology classifications, where majority of AtRbohs and OsRbohs play role in stress-related activities, followed by cellular development. Further, 68 and 38 potential phosphorylation sites were identified in AtRbohs and OsRbohs, respectively. Their distribution, location and kinase specificities were also predicted and correlated with experimental data as well as verified with the other EF-hand containing proteins within both genomes.
Analysis of regulatory mechanisms including interaction with diverse partners and post-translational modifications like phosphorylation have provided insights regarding functional multiplicity of Rbohs. The bioinformatics-based workflow in the current study can be used to get insights for interacting partners and phosphorylation sites from Rbohs of other plant species.
Plants have developed various mechanisms to protect themselves against different stresses whether abiotic or biotic. One of them is the generation of reactive oxygen species (ROS) such as superoxide (O2−), singlet oxygen (1O2), and hydrogen peroxide (H2O2). Membrane-localized NADPH oxidases (Noxes) are the major source of ROS production in plants and transfer electrons from cytosolic NADPH/NADH to apoplastic oxygen which leads to ROS. They are the homolog of the mammalian NADPH oxidase catalytic subunit known as gp91phox . Unlike animals, plant NADPH oxidase consists of two major structural elements: Respiratory burst oxidase homologue (Rboh) and Rop (Rho-like protein; a Rac homologue of plants). The first plant NADPH oxidase was identified in Oryza sativa, known as OsRbohA  and subsequently, more Rbohs were discovered in other plant species including dicots, monocots and lower plants . Rboh proteins contain two Ca2+-binding EF-hand motifs in the N-terminal region, six transmembrane helices, FAD and NADPH binding domains in the C-terminal. The currently available crystal structure of OsRbohB N-terminal region (138–313 amino acid residues) has indicated the presence of two additional EF-hand-like motifs (EF-like 1 and EF-like 2) . Rbohs are known to perform versatile functions in the plant reproduction, growth, development, and responses to abiotic and biotic stresses [1, 4]. Recently published in silico studies on the gene structures, regulatory elements, physico-chemical characterization, topology analysis, phylogenetics and structural analysis of Rbohs have provided critical insights into their diversity and hints to design functional genomics experiments [5,6,7]. Further, few experimental studies have revealed the interaction of Rbohs with various regulatory components for their functioning which involve Ca2+, calcium-dependent protein kinases (CDPKs), Ca2+/CaM-dependent protein kinase (CCaMK), Rop, extracellular ATP (eATP), phospholipase Dα1(PLD α1) and its lipid product phosphatidic acid (PA), mitogen activated protein kinase (MAPK), Nt14–3–3 h/omega1 (a member of 14–3–3 protein family) and nitric oxide [1, 8,9,10]. These interactions may be mediated via complex signaling networks and however, the knowledge regarding the connectivity of Rbohs with these components is still a subject worth investigation.
In addition to interacting partners, functioning of Rbohs through post-translational modification like phosphorylation via various types of protein kinases such as calcium-dependent protein kinase (CDPK), Ca2+ /CaM-dependent protein kinase (CCaMK), mitogen activated protein kinase (MAPK), BIK1 (receptor-like cytoplasmic kinase) has also been observed [1, 11, 12]. However, the various aspects such as potential phosphorylation sites (serine, threonine and tyrosine), their distribution, location and kinase specificities requires extensive experimental studies.
The experimental methods are very time-consuming and expensive, hence the currently available in silico approaches provide alternative cost-effective possibilities to explore the possible interacting associates and phosphorylation sites for Rbohs. In the present study, protein-protein interaction network analysis revealed potential interacting partners for Rbohs from two model plants (10 Rbohs from Arabidopsis and 9 from rice). Further, the potential phosphorylation sites were also elucidated including their distribution, location and kinase specificities and hence correlated with the experimental information wherever available as well as verified with the other EF-hand containing proteins. To the best of our knowledge, this is the first study documenting the potential interacting partners and phosphorylation sites for Rbohs in an extensive manner.
Accession numbers for Arabidopsis and rice Rboh proteins were taken from a recent study from our lab . 10 Rboh sequences for Arabidopsis and 9 for rice were retrieved from UniProt (http://www.uniprot.org/) in FASTA format.
Analysis of protein-protein interaction network
Protein-protein interaction network studies were conducted using STRING v 9.1 (http://string-db.org/). The STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) database retrieves the physical as well as functional interactions among proteins by integrating the information from neighbourhood, gene fusion, co-occurrence, co-expression, experiments, databases, text-mining and homology . The functional interactions were analyzed by using medium confidence score, ranging from 0.4 to 1.0. Interactions with score < 0.4, 0.4 to 0.7 and > 0.7 are considered as low, high and highest confidence respectively. Three different options were used for finding number of interactors: no more than 10, 20 and 50 interactors, and the corresponding confidence scores ranged from 0.865 to 0.99, 0.8 to 0.99 and 0.659 to 0.99, respectively. More information about the functional partners was retrieved from UniProt (http://www.uniprot.org/).
Prediction of phosphorylation sites and kinase specificity
Two programs were used to predict the putative phosphorylation sites in AtRbohs and OsRbohs: Musite (http://musite.net/)  and PlantPhos (http://csb.cse.yzu.edu.tw/PlantPhos/) . We selected two models in Musite for predicting phosphorylating serine and threonine residues at 95% specificity level: General phospho-serine/threonine (A. thaliana) and General phospho-serine/threonine (Green Plants) for AtRbohs and OsRbohs, respectively. General phospho-tyrosine (Green Plants) model was used for predicting phosphorylating tyrosine residues among 19 Rbohs. Default options were employed in PlantPhos. To find any kinase specificity for the predicted sites, NetPhosK 1.0 (http://www.cbs.dtu.dk/services/NetPhosK/)  and KinasePhos 2.0 (http://kinasephos.mbc.nctu.edu.tw/)  programs were used. NetPhosK without ESS filtering method with other default options and KinasePhos with 95% prediction specificity and no specific kinase options were used.
Extraction of experimentally verified phosphorylated sites
RLIMS-P (Rule-based Literature Mining System for Protein Phosphorylation (http://research.bioinformatics.udel.edu/rlimsp/)  was used to extract any phosphorylation related information regarding Rbohs and kinases in the literature. This was also further verified with manual search.
Genome-wide analysis of EF-hand containing proteins in Arabidopsis and rice
Information regarding EF-hand containing proteins in Arabidopsis and rice were retrieved from two previous studies [19, 20]. Gene IDs and Locus IDs were converted to UniProt IDs using UniProt ID Mapping tool, TAIR and RAP-DB. Duplicates were removed manually while analysing the sequences.
Amino acids occurrence percentage and distribution of EF-hand containing proteins in Arabidopsis and rice
Amino acid occurrence percentage and distribution were computed using Residue Frequency Summarizer tool (http://omics.pnl.gov/software/amino-acid-residue-frequency-summarizer) and EMBOSS Pepstats program (http://emboss.sourceforge.net/), respectively.
In the present study, 19 Rboh proteins (10 from Arabidopsis and 9 from rice) were retrieved (Additional file 1) and their interaction partners were determined. Further, phosphorylation sites were predicted, their distribution, location and kinase specificity were analyzed, which was correlated with the available experimental information as well as verified with genome-wide analysis of the other EF-hand containing proteins in both species.
Analysis of protein-protein interaction network for Arabidopsis and rice Rbohs
In order to find the interaction among Rbohs and with other proteins within the respective plant species, a combined interaction network of 10 AtRbohs was constructed, which revealed that out of 10 AtRbohs, only four (AtRbohA, AtRbohB, AtRbohD and AtRbohF) were interacting with CDPKs (calcium dependent protein kinases) and two (AtRbohB and AtRbohF) with OST1 (open stomata 1) (Fig. 1a). AtRbohC was the only one who showed interaction with unique functional partners (not interacting with partners of any other AtRboh). No direct interactions among AtRbohs were observed. Few Rbohs (AtRbohE, AtRbohG, AtRbohH, AtRbohI and AtRbohJ) did not appear to interact with any partners. To get further hints, we also employed two other options “no more than 20 and 50 interactors” in the STRING database. As soon as we increased the number of interactors, we observed few partners for AtRbohI and AtRbohJ (Fig. 1b). AtRbohG which was not showing any interaction, appeared to interact with few partners common with AtRbohA (Fig. 1c). We also generated 10 interaction networks for each of the AtRbohs (Fig. 2). Six Rbohs; AtRbohA, AtRbohC, AtRbohD, AtRbohF, AtRbohG and AtRbohJ appeared to interact with 10 partners, AtRbohA, AtRbohD and AtRbohF had six common partners, AtRbohD and AtRbohF had eight, AtRbohE and AtRbohH had two; and AtRbohH and AtRbohJ had only one common partner. However, AtRbohC showed interaction with 10 unique partners, which was also obtained from combined network. AtRbohG and AtRbohJ have 9 unique partners. Four Rbohs; AtRbohB, AtRbohE, AtRbohH and AtRbohI were interacting with 8, 7, 3 and 6 partners. To get further insights, we also used two other options “no more than 20 and 50 interactors” in the STRING database. Four Rbohs (AtRbohB, AtRbohE, AtRbohH and AtRbohI) showed similar partners as observed in option using “no more than 10 interactors”. However, for another four Rbohs (AtRbohA, AtRbohC, AtRbohD and AtRbohF) more interactions were observed and for rest two (AtRbohJ and AtRbohG) additional one and two interactions were observed, respectively (Additional files 2 and 3). The details of their functional partners were further verified with UniProt (Table 1). Table 2 shows the various unique and common functional partners among AtRbohs in color coding, where unique partners are uncoloured. Further, they were grouped into various functional categories (Additional file 4).
Unlike AtRbohs, the combined network of OsRbohs involved all 9 Rbohs (Additional file 5a). As we go on increasing the number of interactors, we observed more partners for OsRbohs (Additional file 5b, c). However, like AtRbohs no direct interaction among OsRbohs was noticed. Individual networks revealed that 8 OsRbohs have overlapping functional partners (Additional file 6) except OsRbohA which has many unique partners. Unlike AtRbohs, increase in number of interactors for individual Rbohs resulted in more partners for all 9 OsRbohs (Additional files 7 and 8). The details of their functional partners were further verified with UniProt (Additional file 9). Various unique and common functional partners among OsRbohs are shown in Additional file 10. Further, they were grouped into various functional categories (Additional file 11). The various functional categories of potential interaction partners for AtRbohs and OsRbohs are represented in pie chart (Fig. 3).
Analysis of phosphorylation sites for Arabidopsis and rice Rbohs
Potential serine (S), threonine (T) and tyrosine (Y) phosphorylation sites were identified among 19 Rboh proteins (Table 3). Total number of phosphorylation sites predicted in AtRbohs and OsRbohs at 95% specificity level were 68 and 38, respectively. Maximum number of phosphorylation sites for AtRbohs were observed in AtRbohD whereas, in case of rice, it was OsRbohC. Further, the location of predicted sites within the Rboh proteins were identified and both the predicted and few experimentally verified sites were mapped on the multiple sequence alignment of 19 Rbohs (Additional file 12). Most of the potential phosphorylation sites were found in the N-terminal region upstream of EF-hands. For example, equivalent serine residues corresponding to S-148 and S-163 from AtRbohD were conserved in four (AtRbohA, AtRbohC, AtRbohE and AtRbohF) and all AtRbohs, respectively. Further, equivalent serine residues for S-174 from AtRbohF were found conserved among all AtRbohs while S-318 and S-322 from AtRbohC were conserved in 6 AtRbohs (AtRbohA, AtRbohD, AtRbohE, AtRbohF, AtRbohG and AtRbohI) and 3 AtRbohs (AtRbohA, AtRbohB and AtRbohD), respectively. In addition, the kinase specificities/preferences of the putative phosphorylation sites were also computed. NetPhosK and KinasePhos analysis indicated that AtRboh and OsRboh proteins possess a broad range of phosphorylation sites (Table 3). To find the abundance of S, T and Y residues in different regions of 10 AtRbohs and 9 OsRbohs, all Rbohs was divided into four regions: full N-terminal, upstream of EF-hands, EF-hands and downstream of EF-hands. The distribution patterns of S, T and Y for 9 OsRbohs and 10 AtRbohs in each region were computed (Fig. 4a, b). To verify the obtained results, the Gene IDs and Locus IDs for other EF-hand containing proteins of Arabidopsis and rice were retrieved from literature and converted into protein IDs. Overall, 230 Arabidopsis and 211 rice EF-hand containing proteins were obtained. Further, distributions of S, T and Y were analyzed for the representative proteins including 10 from Arabidopsis and 4 from rice (Fig. 5). Similarly other important residues which include lysine (K), arginine (R), proline (P) and cysteine (C) were also analyzed in OsRbohs and AtRbohs (Fig. 6a, b) as well as in the representative proteins (Fig. 7).
Plant NADPH oxidases, also known as respiratory burst oxidase homologs (Rbohs), are critical players in the production of reactive oxygen species (ROS) and play diverse roles [1, 4, 8]. In the recent past, there has been rapidly growing interest to study different aspects of Rbohs using in silico approaches [5,6,7, 21]. Analysis of regulatory mechanisms, for e.g. interaction with different partners and post-translational modification such as phosphorylation are believed to provide vital clues towards functional multiplicity of Rboh proteins. In the present study, protein-protein interaction network analysis revealed potential interacting partners for Arabidopsis and rice Rbohs. In addition, the potential phosphorylation sites were identified with their distribution, location and kinase specificities as well as correlated with experimental data wherever available. This was further verified with the other EF-hand containing proteins throughout the both genomes.
Analysis of protein-protein interaction network for Arabidopsis and rice Rbohs
Rbohs are known to mediate diverse functions and are implicated in many signal transduction pathways involving interactions with different components involved in their regulation such as Ca2+, protein kinases (CDPKs; Ser/Thr protein kinases having a Ca2+ − binding calmodulin-like domain, OST1, CBL/CIPK; calcineurin B-Like calcium sensors–interacting protein kinase) and Rop; Rac of plants . These interactions may be mediated in the form of well-regulated networks and the information regarding the connectivity of Rbohs with these components is still the subject of investigation. Therefore, it is imperative to decipher these interactions among Rbohs and with other genes/proteins. Hence, we performed protein-protein interaction network analysis using the STRING database to find any functional associations among Rbohs and specific Rboh with other proteins. The STRING stands for Search Tool for the Retrieval of Interacting Genes/Proteins that retrieves the physical and functional interactions among proteins by integrating the information from a variety of evidence types (neighbourhood, gene fusion, co-occurrence, co-expression, experiments, databases, text-mining and homology) . In the present work, the total number of partners showing interactions with AtRbohs and OsRbohs were 59 and 19, respectively. Rboh proteins with a high degree of connectivity i.e. connected to many other proteins would act as a central point involved in supervising communication in a network. On the other hand, the nodes that are not connected indicates no interaction with other proteins and hence, involved in an indirect role. A combined STRING network analysis of 10 AtRbohs predicted that various types of CDPKs were common partners for four Rbohs (AtRbohA, AtRbohB, AtRbohD and AtRbohF). However, individual network analysis of each AtRboh revealed more details of partners. AtRbohA and AtRbohD were predicted to interact with 10 CDPKs, AtRbohB with one CDPK and AtRbohF with nine CDPKs. AtRbohA, AtRbohD and AtRbohF had six common CDPKs indicating the diversity as well as overlapping in their functions. Except AtRbohA whose function is still unknown, AtRbohD and AtRbohF are known to play multiple roles . The reason might be the differences in the substrate specificities of AtCDPKs which may lead to their functional multiplicity [22, 23]. However, there are no experimentally validated interactions of CDPK with AtRbohA, AtRbohB and AtRbohF, but with that of AtRbohD have been recently reported . It has been suggested that the flg22 induced activation of CPK5 leading to the phosphorylation of AtRbohD. In the predicted networks, CPK5 was found to interact with AtRbohD at confidence score of 0.8 when we used “no more than 20 and 50 interactors” options. In addition, CPK5 was also observed to interact with AtRbohA and AtRbohF using same options. Another protein kinase OST1 was common for two Rbohs (AtRbohB and AtRbohF) with better confidence score (0.99) for AtRbohF than AtRbohB (0.408). Although its interaction with AtRbohF is known in stomatal closure , however no experimentally recorded interactions with AtRbohB are available so far. Another class of protein kinases consists of Ca2+-regulated CBL-interacting protein kinases (CIPKs) that are known to be activated upon interaction with Calcineurin B-like (CBL) Ca2+ sensor proteins. Several diverse functions have been documented for various CBL-CIPK pairs . In the present study, only AtRbohG was predicted to interact with CIPK25. A recent report has suggested the role of CIPK25 in biotic stress and energy sensing [27, 28]. However, another study has reported the interaction of AtRbohF with CBL1/9-CIPK26 complex and its involvement in the regulation of ABA responses .
In addition to protein kinases, AtRbohs also showed interaction with different Rac proteins. Rac are known to display diverse array of functions in the plants . In our study, the interactions between AtRbohA and AtRbohB with Rac6 and AtRboh C with Rop2, Rac3 and Rac5 were observed.
Among 10 AtRbohs, AtRbohC appeared to have unique functional partners, all of which are involved in root development. The role of AtRbohC in root hair formation is well known , however no experimental reports are available regarding its interaction partners. On the other hand, AtRbohH was observed to interact with only three partners (ATRDH2, AT2G35040 and GLP4). An earlier study has documented the role of ATRDH2 (also known as Rhodanese Homologue 2 or STR2; Sulfurtransferase 2) in embryo and seed development . On the other hand, GLP4 (Germin-Like Protein 4) is implicated in the plant defense, auxin-induced cell growth and exhibiting superoxide dismutase activity [33,34,35]. Although, the role of AT2G35040 has not been elucidated, it may belong to AICARFT/IMPCHase bienzyme family as evident from our UniProt analysis. Two partners (AT2G35040 and GLP4) of AtRbohH are common with AtRbohE and AtRbohI, and one (AT2G35040) with AtRbohJ. Previous reports have indicated AtRbohH and AtRbohJ as pollen-specific , which may be related to one common partner (AT2G35040) among them, however, function of AtRbohE is still unknown. In addition to two common partners (AT2G35040 and GLP4) with AtRbohH, AtRbohE has one common partner (CYP707A4) with two other Rbohs (AtRbohG and AtRbohI). CYP707A4 is an abscisic acid 8′-hydroxylase 4 which is involved in ABA catabolism during drought conditions . Further, AtRbohE appeared to interact with four unique partners (DGL, ATEXPA2, AT1G17680 and GA20OX3). Two (DGL and GA20OX3) are involved in biosynthesis of hormones, where DGL (DONGLE) possess galactolipase activity and involves in jasmonic acid biosynthesis  while GA20OX3 (Gibberellin 20 oxidase 3) in gibberellin acid biosynthesis . ATEXPA2 (α-expansin 2) plays role in cell wall loosening and development processes , while AT1G17680 is still uncharacterized. Further, AtRbohG along with AtRbohI and AtRbohJ were also observed to interact with many unique partners.
In case of rice Rbohs, except OsRbohA which has five unique and five common partners, the observed frequency of common partners was higher among other eight Rbohs. The unique interaction partners for OsRbohA involves immutans (chloroplastic alternative oxidase), three superoxide dismutases (FeSOD1, FeSOD2 and MnSOD) and protein kinase WNK4. Immutans are known in ROS-related damage prevention under high light stress , FeSOD1 and FeSOD2 in drought stress; and MnSOD in drought stress and embryogenesis [42, 43]. However, the role of WNK4 is still unknown. Another protein, two pore calcium channel protein1 (TPC1) was observed to interact with OsRbohB only. The role of TPC1 in the regulation of growth and development is well documented . Other partner (4339304) encoding Rac protein expressing under salt stress, was found to interact with two rice Rbohs (OsRbohA and OsRbohB) only . It is interesting to note that one partner with gene ID: LOC_Os04g31290.1 (basic helix-loop-helix (bHLH) DNA-binding domain containing protein) was found to interact with all nine rice Rbohs and acting as a central hub. However, no information is available for its functional annotation yet. Functional characterization of such hub will lead to addition of knowledge in the area.
Overall, the present interaction study and gene ontology classifications have provided insights into the interaction of AtRbohs and OsRbohs with a wide range of potential partners which may be critical for their diverse functions. The observed high frequency of common and total partners in AtRbohs as compared to OsRbohs may indicate more complex interactions in AtRbohs. Further, pie distribution indicated that the majority of AtRbohs play role in stress-related activities, followed by cellular development. Similar kind of trend was also obtained for OsRbohs. These observations justifies the versatility of functions played by Rbohs as evident from literature .
Analysis of phosphorylation sites for Arabidopsis and rice Rbohs
Protein phosphorylation and dephosphorylation are among the most crucial post-translational modifications, which play important role in a broad range of regulatory signaling cascades in plants. The phosphorylation of specific sites in proteins may result in conformational changes in protein structure which may lead to changes in enzyme activity, substrate specificity, biological role, intracellular localization, protein stability etc. Serine, threonine and tyrosine residues are the important amino acids which can be phosphorylated. Few studies in the past have provided hints for the regulation of Rbohs through phosphorylation via Ca2+ and different types of protein kinases such as calcium-dependent protein kinase (CDPK), Ca2+ /CaM-dependent protein kinase (CCaMK), mitogen activated protein kinase (MAPK), BIK1 (receptor-like cytoplasmic kinase) etc. [1, 9, 11, 12]. However, this area need to be fully investigated, though it is very time-consuming and expensive to identify a broad range of phosphorylation sites experimentally. Hence, in silico prediction of phosphorylation sites provides an alternative approach . In the present study, potential serine (S), threonine (T) and tyrosine (Y) phosphorylation sites were predicted. Their distribution and location were studied by mapping them on the multiple sequence alignment of 19 Rbohs. It was interesting to note that the experimentally verified sites for few AtRbohs were also obtained as potential phosphorylation sites in our prediction outputs and hence, provided hints for conservation and variability within other AtRbohs. Earlier evidences have indicated the diverse roles of AtRbohD and AtRbohF in plants involving growth and development, abiotic and biotic stresses . AtRbohD was found to be phosphorylated at S-8, S-39, S-148, S-152, S-163, S-343 and S-347 in response to pathogen elicitors [24, 47, 48]. In the present study, the conservation of equivalent serine residues corresponding to S-148 in four and S-163 in all AtRbohs, suggest N-terminal phosphorylation-mediated regulation in them. Further, S-133 and S-148 from AtRbohD corresponds to S-82 and S-97 from S. tuberosum StRbohB, which have been identified as potential phosphorylation sites for StCDPK4 and StCDPK5 . An earlier study showed that S-13 and S-174 from AtRbohF are phosphorylated by OST1 protein kinase, where OST1 is known to phosphorylate S/T from the motif [LIMVF]XRXXS/T [25, 50]. In the present work, the conservation of S-174 among all AtRbohs as well as that of arginine (R) at the − 3 position relative to S-174 may indicate that they can be phosphorylated by OST1 kinase. S-174 from AtRbohF corresponds to that of S-163 from AtRbohD. Some line of evidences have reported the functional redundancy among AtRbohD and AtRbohF, however, AtRbohD is mostly responsible for ROS in plant-pathogen interactions while AtRbohF in ABA signaling [51, 52]. This might be due to variation among few phosphorylation sites, which lead to differential regulation and function. Besides AtRbohD and AtRbohF, the function of AtRbohC has also been elucidated which is involved in root-hair development . An earlier study has identified S-318 and S-322 as possible sites of phosphorylation in AtRbohC . In the present study, the conservation of equivalent serine residues corresponding to S-318 in six and S-322 in three AtRbohs, may suggest their potential role in root development. The experimental information regarding phosphorylation for other AtRbohs and OsRbohs is still lacking. However, our mapping of the predicted potential phosphorylation sites among all 19 Rbohs indicates their abundance in the N-terminal region with higher phosphorylation of serine residues as compared to threonine and tyrosine. Further, we computed the distribution patterns of all S, T and Y residues in different regions of 19 Rbohs. It was observed that the frequency of occurrence of these residues is highest in the N-terminal upstream EF-hand region from both AtRbohs and OsRbohs. To know whether this concept holds good for other EF-hand containing proteins, genome-wide analysis of EF-hand containing proteins from Arabidopsis and rice was carried out. The distribution pattern of S, T and Y were analyzed for few representative proteins from Arabidopsis and rice. Interestingly, we obtained similar kind of pattern for the selected proteins, which verified our results and hence may point towards the critical role of the N-terminal upstream EF-hand region in contributing functional multiplicity to AtRbohs and OsRbohs. In addition to S, T and Y, other important residues including lysine (K), arginine (R), proline (P) and cysteine (C) were also studied. These residues are known to play crucial roles in transmembrane proteins such as lysine and arginine provide assistance in anchoring the transmembrane orientations, proline in stress tolerance and cysteine in membrane localization [54,55,56,57]. Like S, T and Y, the frequency of occurrence of these residues (K, R, P and C) was found highest in the N-terminal upstream EF-hand region for AtRbohs, OsRbohs and representative proteins.
The current comprehensive in silico study provides a necessary clue that may be the N-terminal amino acid residues from Rbohs and their phosphorylation are very critical for regulating various biological functions in a plant. A recent study on Rbohs sequence and structural analysis has also provided hints towards the role of the N-terminal and its variability for their functional diversity . It will be further interesting to test experimentally the interaction with the potential partners, the interaction sites, any overlapping regions, the role of correct intracellular location of the partner, the phosphorylation ability and preferences of the predicted residues among various Rbohs.
In the present work, in silico approaches were followed to comprehensively deduce the possible interacting partners and phosphorylation sites of Rboh gene family from two model plants (Arabidopsis and rice). The study elucidates an extensive range of potential partners revealed from interaction networks and gene ontology classifications, which may be responsible for their functional multiplicity. Further insights were also obtained from the prediction of potential phosphorylation sites as well as their distribution, location and kinase specificities. These results were correlated with experimental data as well as verified with the other EF-hand containing proteins. However, more and more inputs from the experimental work will further strengthen our assumptions and pave the way to modulate plant species to address the future challenges, for instance, crops with better stress adaptibility.
Kaur G, Sharma A, Guruprasad K, Pati PK. Versatile roles of plant NADPH oxidases and emerging concepts. Biotechnol Adv. 2014;32(3):551–63.
Groom QJ, Torres MA, Fordham-Skelton AP, Hammond-Kosack KE, Robinson NJ, Jones JD. rbohA, a rice homologue of the mammalian gp91phox respiratory burst oxidase gene. Plant J. 1996;10(3):515–22.
Oda T, Hashimoto H, Kuwabara N, Akashi S, Hayashi K, Kojima C, Wong HL, Kawasaki T, Shimamoto K, Sato M, et al. Structure of the N-terminal regulatory domain of a plant NADPH oxidase and its functional implications. J Biol Chem. 2010;285(2):1435–45.
Jiménez-Quesada MJ, Traverso JÁ, Alché JD. NADPH oxidase-dependent superoxide production in plant reproductive tissues. Front Plant Sci. 2016;7:359.
Kaur G, Pati PK. In silico physicochemical characterization and topology analysis of respiratory burst oxidase homolog (Rboh) proteins from Arabidopsis and rice. Bioinformation. 2018;14(2):93–100.
Kaur G, Guruprasad K, Temple BRS, Shirvanyants DG, Dokholyan NV, Pati PK. Structural complexity and functional diversity of plant NADPH oxidases. Amino Acids. 2018;50(1):79–94.
Kaur G, Pati PK. Analysis of cis-acting regulatory elements of respiratory burst oxidase homolog (Rboh) gene families in Arabidopsis and rice provides clues for their diverse functions. Comput Biol Chem. 2016;62:104–18.
Baxter A, Mittler R, Suzuki N. ROS as key players in plant stress signalling. J Exp Bot. 2014;65(5):1229–40.
Kimura S, Kaya H, Kawarazaki T, Hiraoka G, Senzaki E, Michikawa M, Kuchitsu K. Protein phosphorylation is a prerequisite for the Ca2+−dependent activation of Arabidopsis NADPH oxidases and may function as a trigger for the positive feedback regulation of Ca2+ and reactive oxygen species. Biochim Biophys Acta. 2012;1823(2):398–405.
Potocký M, Pejchar P, Gutkowska M, Jiménez-Quesada MJ, Potocká A, de Dios AJ, Kost B, Žárský V. NADPH oxidase activity in pollen tubes is affected by calcium ions, signaling phospholipids and Rac/Rop GTPases. J Plant Physiol. 2012;169(16):1654–63.
Kadota Y, Sklenar J, Derbyshire P, Stransfeld L, Asai S, Ntoukakis V, Jones JDG, Shirasu K, Menke F, Jones A. Direct regulation of the NADPH oxidase RBOHD by the PRR-associated kinase BIK1 during plant immunity. Mol Cell. 2014;54(1):43–55.
Li L, Li M, Yu L, Zhou Z, Liang X, Liu Z, Cai G, Gao L, Zhang X, Wang Y. The FLS2-associated kinase BIK1 directly phosphorylates the NADPH oxidase RbohD to control plant immunity. Cell Host Microbe. 2014;15(3):329–38.
Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, Lin J, Minguez P, Bork P, von Mering C. STRING v9. 1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013;41(D1):D808–15.
Gao J, Thelen JJ, Dunker AK, Xu D. Musite, a tool for global prediction of general and kinase-specific phosphorylation sites. Mol Cell Proteomics. 2010;9(12):2586–600.
Lee T-Y, Bretaña NA, Lu C-T. PlantPhos: using maximal dependence decomposition to identify plant phosphorylation sites with substrate site specificity. BMC Bioinformatics. 2011;12(1):261.
Blom N, Sicheritz-Pontén T, Gupta R, Gammeltoft S, Brunak S. Prediction of post-translational glycosylation and phosphorylation of proteins from the amino acid sequence. Proteomics. 2004;4(6):1633–49.
Huang H-D, Lee T-Y, Tzeng S-W, Horng J-T. KinasePhos: a web tool for identifying protein kinase-specific phosphorylation sites. Nucleic Acids Res. 2005;33(suppl 2):W226–9.
Hu Z-Z, Narayanaswamy M, Ravikumar KE, Vijay-Shanker K, Wu CH. Literature mining and database annotation of protein phosphorylation using a rule-based system. Bioinformatics. 2005;21(11):2759–65.
Boonburapong B, Buaboocha T. Genome-wide identification and analyses of the rice calmodulin and related potential calcium sensor proteins. BMC Plant Biol. 2007;7(1):4.
Day IS, Reddy VS, Shad Ali G, Reddy AS. Analysis of EF-hand-containing proteins in Arabidopsis. Genome Biol. 2002;3(10):1–0056.
Jiménez-Quesada MJ, Traverso JÁ, Zafra A, Jimenez-Lopez JC, Carmona R, Claros MG, de Dios Alché J. Identification and in silico Analysis of NADPH Oxidase Homologues Involved in Allergy from an Olive Pollen Transcriptome. In: International Conference on Bioinformatics and Biomedical Engineering. Cham: Springer; 2015 p. 450–59.
Curran A, Chang F, Chang C-L, Garg S, Miguel RM, Barron YD, Li Y, Romanowsky S, Cushman JC, Gribskov M. Calcium-dependent protein kinases from Arabidopsis show substrate specificity differences in an analysis of 103 substrates. Frontiers Plant Sci. 2011;2:36.
Boudsocq M, Sheen J. CDPKs in immune and stress signaling. Trends Plant Sci. 2013;18(1):30–40.
Dubiella U, Seybold H, Durian G, Komander E, Lassig R, Witte CP, Schulze WX, Romeis T. Calcium-dependent protein kinase/NADPH oxidase activation circuit is required for rapid defense signal propagation. P Natl Acad Sci USA. 2013;110(21):8744–9.
Sirichandra C, Gu D, Hu HC, Davanture M, Lee S, Djaoui M, Valot B, Zivy M, Leung J, Merlot S, et al. Phosphorylation of the Arabidopsis AtrbohF NADPH oxidase by OST1 protein kinase. FEBS Lett. 2009;583(18):2982–6.
Luan S. The CBL–CIPK network in plant calcium signaling. Trends Plant Sci. 2009;14(1):37–42.
Lee SC, Mustroph A, Sasidharan R, Vashisht D, Pedersen O, Oosumi T, Voesenek LACJ, Bailey-Serres J. Molecular characterization of the submergence response of the Arabidopsis thaliana ecotype Columbia. New Phytol. 2011;190(2):457–71.
Siemens J, Keller I, Sarx J, Kunz S, Schuller A, Nagel W, Schmulling T, Parniske M, Ludwig-Muller J. Transcriptome analysis of Arabidopsis clubroots indicate a key role for cytokinins in disease development. Mol Plant Microbe In. 2006;19(5):480–94.
Drerup MM, Schlucking K, Hashimoto K, Manishankar P, Steinhorst L, Kuchitsu K, Kudla J. The Calcineurin B-like calcium sensors CBL1 and CBL9 together with their interacting protein kinase CIPK26 regulate the Arabidopsis NADPH oxidase RBOHF. Mol Plant. 2013;6(2):559–69.
Craddock C, Lavagi I, Yang ZB. New insights into rho signaling from plant ROP/Rac GTPases. Trends Cell Biol. 2012;22(9):492–501.
Foreman J, Demidchik V, Bothwell JH, Mylona P, Miedema H, Torres MA, Linstead P, Costa S, Brownlee C, Jones JD, et al. Reactive oxygen species produced by NADPH oxidase regulate plant cell growth. Nature. 2003;422(6930):442–6.
Mao GH, Wang RG, Guan YF, Liu YD, Zhang SQ. Sulfurtransferases 1 and 2 play essential roles in embryo and seed development in Arabidopsis thaliana. J Biol Chem. 2011;286(9):7548–57.
Christensen AB, Thordal-Christensen H, Zimmermann G, Gjetting T, Lyngkjær MF, Dudler R, Schweizer P. The germinlike protein GLP4 exhibits superoxide dismutase activity and is an important component of quantitative resistance in wheat and barley. Mol Plant Micro Interact. 2004;17(1):109–17.
Davidson RM, Reeves PA, Manosalva PM, Leach JE. Germins: a diverse protein family important for crop improvement. Plant Sci. 2009;177(6):499–510.
Yin K, Han X, Xu Z, Xue H. Arabidopsis GLP4 is localized to the Golgi and binds auxin in vitro. Acta Biochim Biophys Sin. 2009;41(6):478–87.
Potocký M, Jones MA, Bezvoda R, Smirnoff N, Zarsky V. Reactive oxygen species produced by NADPH oxidase are involved in pollen tube growth. New Phytol. 2007;174(4):742–51.
Kushiro T, Okamoto M, Nakabayashi K, Yamagishi K, Kitamura S, Asami T, Hirai N, Koshiba T, Kamiya Y, Nambara E. The Arabidopsis cytochrome P450 CYP707A encodes ABA 8′-hydroxylases: key enzymes in ABA catabolism. EMBO J. 2004;23(7):1647–56.
Hyun Y, Choi S, Hwang H-J, Yu J, Nam S-J, Ko J, Park J-Y, Seo YS, Kim EY, Ryu SB. Cooperation and functional diversification of two closely related galactolipase genes for jasmonate biosynthesis. Dev Cell. 2008;14(2):183–92.
Han FM, Zhu BG. Evolutionary analysis of three gibberellin oxidase genesin rice, Arabidopsis, and soybean. Gene. 2011;473(1):23–35.
Sampedro J, Cosgrove DJ. The expansin superfamily. Genome Biol. 2005;6(12):242.
Rizhsky L, Hallak-Herr E, Van Breusegem F, Rachmilevitch S, Barr JE, Rodermel S, Inzé D, Mittler R. Double antisense plants lacking ascorbate peroxidase and catalase are less sensitive to oxidative stress than single antisense plants lacking ascorbate peroxidase or catalase. Plant J. 2002;32(3):329–42.
Song ZZ, Yang SY, Zuo J, Su YH. Over-expression of ApKUP3 enhances potassium nutrition and drought tolerance in transgenic rice. Biol Plant. 2014;58(4):649–58.
Zi J, Zhang J, Wang Q, Zhou B, Zhong J, Zhang C, Qiu X, Wen B, Zhang S, Fu X. Stress responsive proteins are actively regulated during rice (Oryza sativa) embryogenesis as indicated by quantitative proteomics analysis. PLoS One. 2013;8(9):e74229.
Kurusu T, Sakurai Y, Miyao A, Hirochika H, Kuchitsu K. Identification of a putative voltage-gated Ca2+−permeable channel (OsTPC1) involved in Ca2+ influx and regulation of growth and development in rice. Plant Cell Physiol. 2004;45(6):693–702.
Walia H, Wilson C, Ismail AM, Close TJ, Cui XP. Comparing genomic expression patterns across plant species reveals highly diverged transcriptional dynamics in response to salt stress. BMC Genomics. 2009;10:398.
Trost B, Kusalik A. Computational prediction of eukaryotic phosphorylation sites. Bioinformatics. 2011;27(21):2927–35.
Benschop JJ, Mohammed S, O'Flaherty M, Heck AJ, Slijper M, Menke FL. Quantitative phosphoproteomics of early elicitor signaling in Arabidopsis. Mol Cell Proteomics. 2007;6(7):1198–214.
Nühse TS, Bottrill AR, Jones AM, Peck SC. Quantitative phosphoproteomic analysis of plasma membrane proteins reveals regulatory mechanisms of plant innate immune responses. Plant J. 2007;51(5):931–40.
Kobayashi M, Ohura I, Kawakita K, Yokota N, Fujiwara M, Shimamoto K, Doke N, Yoshioka H. Calcium-dependent protein kinases regulate the production of reactive oxygen species by potato NADPH oxidase. Plant Cell. 2007;19(3):1065–80.
Vlad F, Turk BE, Peynot P, Leung J, Merlot S. A versatile strategy to define the phosphorylation preferences of plant protein kinases and screen for putative substrates. Plant J. 2008;55(1):104–17.
Kwak JM, Mori IC, Pei ZM, Leonhardt N, Torres MA, Dangl JL, Bloom RE, Bodde S, Jones JD, Schroeder JI. NADPH oxidase AtrbohD and AtrbohF genes function in ROS-dependent ABA signaling in Arabidopsis. EMBO J. 2003;22(11):2623–33.
Torres MA, Dangl JL, Jones JD. Arabidopsis gp91phox homologues AtrbohD and AtrbohF are required for accumulation of reactive oxygen intermediates in the plant defense response. P Natl Acad Sci USA. 2002;99(1):517–22.
Takeda S, Gapper C, Kaya H, Bell E, Kuchitsu K, Dolan L. Local positive feedback regulation determines cell shape in root hair cells. Science. 2008;319(5867):1241–4.
Fukuda M, Mikoshiba K. The N-terminal cysteine cluster is essential for membrane targeting of B/K protein. Biochem J. 2001;360:441–8.
Gleason NJ, Vostrikov VV, Greathouse DV, Koeppe RE. Buried lysine, but not arginine, titrates and alters transmembrane helix tilt. P Natl Acad Sci USA. 2013;110(5):1692–5.
Hristova K, Wimley WC. A look at arginine in membranes. J Membr Biol. 2011;239(1–2):49–56.
Szabados L, Savouré A. Proline: a multifunctional amino acid. Trends Plant Sci. 2010;15(2):89–97.
The authors acknowledge the financial support received under Innovation in Science Pursuit for Inspired Research (INSPIRE) Programme, Department of Science and Technology (DST), Government of India [Grant No. DST/INSPIRE Fellowship/2010]. The financial assistance from Department of Biotechnology (DBT), Government of India (Grant No. BT/PR13965/BRB/10/883/2010) is also acknowledged.
Availability of data and materials
The datasets generated and/or analyzed during this study are included in this article [and its supplementary files].
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table showing Arabidopsis and rice Rboh protein sequences retrieved from UniProt. (PDF 13 kb)
Individual network diagram of few AtRboh proteins showing potential interacting partners in evidence view using no more than 20 interactors option. Different coloured lines indicate types of evidence for association. The thickness of each line indicates the strength of the association. (PDF 549 kb)
Individual network diagram of few AtRboh proteins showing potential interacting partners in evidence view using no more than 50 interactors option. Different coloured lines indicate types of evidence for association. The thickness of each line indicates the strength of the association. (PDF 1926 kb)
Table showing functional categorization of identified interaction partners among AtRbohs. (PDF 303 kb)
Combined network diagram of 9 OsRboh proteins showing potential interacting partners in evidence view with (a) 10 interactors (b) 20 interactors and (c) 50 interactors. Different coloured lines indicate types of evidence for association. The thickness of each line indicates the strength of the association. (PDF 811 kb)
Individual network diagrams of 9 OsRboh proteins showing potential interacting partners in evidence view with 10 interactors. Different coloured lines indicate types of evidence for association. The thickness of each line indicates the strength of the association. (PDF 1978 kb)
Individual network diagram of OsRboh proteins showing potential interacting partners in evidence view using no more than 20 interactors option. Different coloured lines indicate types of evidence for association. The thickness of each line indicates the strength of the association. (PDF 855 kb)
Individual network diagram of OsRboh proteins showing potential interacting partners in evidence view using no more than 50 interactors option. Different coloured lines indicate types of evidence for association. The thickness of each line indicates the strength of the association. (PDF 491 kb)
Table showing details of identified potential functional partners of OsRboh proteins. (PDF 259 kb)
Table showing various unique and common functional partners among OsRbohs in color coding, where unique partners are uncoloured. (PDF 101 kb)
Table showing functional categorization of identified interaction partners among OsRbohs. (PDF 196 kb)
Mapping of phosphorylation sites on multiple sequence alignment of AtRboh and OsRboh protein sequences. 19 Rboh sequences were aligned with Clustal Omega (http://www.ebi.ac.uk/Tools/msa/clustalo). Gaps were indicated with dashes in the sequences. CDPK-binding motifs, two EF-hand-likes, two EF- hands, six ransmembrane spanning domains (TMD I-VI), two FAD-binding and four NADPH-binding sites were shown in black boxes. Reported (also retrieved in prediction) S/T/Y sites involved in phosphorylation were highlighted in 654 green while predicted sites were indicated in cyan. (PDF 62 kb)
About this article
Cite this article
Kaur, G., Pati, P.K. In silico insights on diverse interacting partners and phosphorylation sites of respiratory burst oxidase homolog (Rbohs) gene families from Arabidopsis and rice. BMC Plant Biol 18, 161 (2018). https://doi.org/10.1186/s12870-018-1378-2