https://orcid.org/0000-0002-6336-1140
Figures
Abstract
Nitrate transport in cyanobacteria is mediated by ABC-transporter, which consists of a highly conserved ATP binding cassette (ABC) and a less conserved transmembrane domain (TMD). Under salt stress, recombinant glycinebetaine (GB) not only protected the rate of nitrate transport in transgenic Anabaena PCC 7120, rather stimulated the rate by interacting with the ABC-transporter proteins. In silico analyses revealed that nrtA protein consisted of 427 amino acids, the majority of which were hydrophobic and contained a Tat (twin-arginine translocation) signal profile of 34 amino acids (1–34). The nrtC subunit of 657 amino acids contained two hydrophobic distinct domains; the N-terminal (5–228 amino acids), which was 59% identical to nrtD (the ATP-binding subunit) and the C-terminal (268–591), 28.2% identical to nrtA, suggesting C-terminal as a solute binding domain and N-terminal as ATP binding domain. Subunit nrtD consisted of 277 amino acids and its N-terminal (21–254) was an ATP binding motif. Phylogenetic analysis revealed that nitrate-ABC-transporter proteins are highly conserved among the cyanobacterial species, though variation existed in sequences resulting in several subclades. Nostoc PCC 7120 was very close to Anabaena variabilis ATCC 29413, Anabaena sp. 4–3 and Anabaena sp. CA = ATCC 33047. On the other, Nostoc spp. NIES-3756 and PCC 7524 were often found in the same subclade suggesting more work before referring it to Anabaena PCC 7120 or Nostoc PCC 7120. The molecular interaction of nitrate with nrtA was hydrophilic, while hydrophobic with nrtC and nrtD. GB interaction with nrtACD was hydrophobic and showed higher affinity compared to nitrate.
Citation: Swapnil P, Meena M, Rai AK (2021) Molecular interaction of nitrate transporter proteins with recombinant glycinebetaine results in efficient nitrate uptake in the cyanobacterium Anabaena PCC 7120. PLoS ONE 16(11): e0257870. https://doi.org/10.1371/journal.pone.0257870
Editor: Anjani Kumar Tiwari, Babasaheb Bhimrao Ambedkar University (A Central University), INDIA
Received: August 2, 2020; Accepted: September 14, 2021; Published: November 18, 2021
Copyright: © 2021 Swapnil et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting information files.
Funding: The authors received no specific funding for this work.
Competing interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests. The authors declare that there is no conflict of interest.
Abbreviations: ABC, ATP-binding cassette; GB, Glycinebetaine; MEGA, Molecular Evolutionary Genetics Analysis; N2, Dinitrogen; NaCl, Sodium Chloride; NCBI, National Centre for Biotechnology Information; Nrt, Nitrate transporter; pBLAST, Basic Local Alignment Search Tool; PDB, Protein Data Bank; PHYLIP, PHYLogeny Inference Package; PMDB, Protein Model Data Base; ProSA, Protein Structure Analysis; TMDs, Transmembrane domains; UPGMA, Unweighted Pair Group Method with Arithmetic Mean; WT, Wild type
1. Introduction
Nitrogen is a fundamental nutrient required by plants and microorganisms in high amounts and is often a limiting factor for their growth and yield. Except for N2-fixing prokaryotes and those living in association with them, the rest of the life forms depend upon combined forms of nitrogen viz., nitrate, nitrite, ammonium, and organic nitrogen. Nitrate is the most oxidized and stable form of combined-N, and is predominantly available as an immediate source of N for most of the plants and microorganisms in terrestrial and aquatic ecosystems. In the presence of nitrate, N2-fixing cyanobacteria cease the energy effective process of N2-fixation and switch over to utilizing available nitrate. However, nitrate concentration in natural environments is relatively low (< 0.01–50 μM) [1]. Hence, the initial step in nitrate assimilation, uptake of nitrate into the cell may become limiting. Cyanobacteria possess an active nitrate transport system to meet their N requirement [2]. The high-affinity ABC-nitrate transporter in cyanobacteria belongs to the ATP-binding cassette (ABC) super-family and comprises four domains (Fig 1); two ABC proteins: (1) a cytoplasmic ATPase (nrtD) and (2) ATPase/ nitrate binding fusion protein (nrtC), and two TMDs (transmembrane domains) proteins, (3) the periplasmic nitrate binding lipoprotein (nrtA) and (4) an integral membrane permease (nrtB) [3, 4]. The TMDs possess additional short extracytoplasmic binding proteins that provide critical contact to ABCs to function [5]. After entry into the cells, nitrate is reduced to nitrite and ammonium in two sequential enzymic reductions mediated by nitrate reductase and nitrite reductase, respectively, and then incorporated into carbon skeletons through GS/GOGAT (glutamine synthetase/ glutamine 2-oxoglutarate aminotransferase) pathway yielding amino acids.
(A) The TMD proteins of the ABC-nitrate transporter, periplasmic subunit nrtA scavenges nitrate of the solution and transports across the lipid bilayer of the membrane through nrtB pore. When nitrate level rises in the cell, ABC protein, NrtC binds to nitrate and acts as negative regulator by inhibiting nitrate uptake, while nrtD binds and hydrolyzes ATP to induce conformational changes in TMDs (nrtAB). Extra cytoplasmic helices running parallel to the cytoplasmic membrane provide contact to TMDs proteins. (B) Interaction of GB with nrtACD proteins.
In cyanobacteria, nitrate assimilation is directly linked to photosynthesis and is well documented [6]. In addition to the consumption of ATP during the uptake of nitrate, 8 e- derived from photolysis of water are donated by ferredoxin to convert a molecule of nitrate into one ammonium molecule. However, most of the studies on nitrate assimilation are centered on the latter stages of nitrate assimilation, the first step being the transport of nitrate into the cell is the least studied [7]. Until recently, the mechanism of the ABC-nitrate transporter proteins in cyanobacteria has been overlooked. Salinity impairs nitrate uptake/ assimilation in cyanobacteria [6, 8]. Based on the kinetics and depressed transport of nitrate by an ionophore (monensin), sodium/nitrate symport has been proposed in Anacystis nidulans R2 [9, 10]. A mutant of Anabaena PCC 7120 defective in nitrate uptake could not take up Na+, indicated for NO3/ Na+ symport [11]. Conversely, no ABC-transporter driven by the electrochemical gradient of an ion has so far been reported. Furthermore, nitrate is noted to increase the adaptability of cyanobacterial cells to salinity by interacting with Na+ carrier and limiting Na+ influx [12–14]. To overcome the salt toxicity, halophilic and halotolerant forms synthesize/ accumulate low-molecular-weight organic compatible solutes. Glycinebetaine (N,N,N-trimethyl glycine), hereafter referred to as GB is a universal compatible solute synthesized by both prokaryotes and eukaryotes. In addition to maintaining cell osmolarity, GB is suggested to act as molecular chaperones protecting the native structure of proteins and membrane integrity [15].
In the present investigation, we used a freshwater diazotrophic cyanobacterium Anabaena PCC 7120 wild type (WT) and its transformant possessing N-methyltransferase genes (ApGSMT-DMT; glycine sarcosine methyltransferase and sarcosine dimethylglycine methyltransferase) from a halophilic cyanobacterium Aphanothece halophytica [6, 16]. The transformed cyanobacterium de novo synthesized and accumulated GB. We noted that salinity reduced the rate of nitrate transport in WT cells, while the rate increased in the ApGSMT-DMT transformed cells that synthesized recombinant GB. An attempt is made to unveil the mechanism(s) by which recombinant GB increases the efficiency of nitrate transporter system, and consequently its assimilation.
2. Materials and methods
Cyanobacterium Anabaena sp. PCC 7120 (hereafter referred to as Anabaena 7120; also known as Nostoc sp. PCC 7120) WT cells were grown and maintained in BG110 medium [17]. The transformed Anabaena 7120 containing ApGSMT-DMT genes from Aphanothece halophytica (on the expression cassette pRL488 + ApGSMT-DMT) was grown in BG110 + 0.1 M NaCl + Km10 (kanamycin 10 μg ml-1), since it was unable to grow without NaCl [6, 16]. The cultures were grown at 28 ± 1°C with an average illumination of 70 μE m2 s1 provided with daylight fluorescent lamps for 14 h a day.
2.1. Nitrate uptake
Exponential phase Anabaena 7120 WT and transgenic cells growing on N2 in the presence and absence of 0.1 M NaCl were harvested and suspended in their respective fresh media to get a cell density of 70–80 μg chl ml-1. To study the kinetics of nitrate uptake, the concentration of nitrate in the uptake mixture (nutrient solution with and without 0.1 M NaCl) was varied between 40–500 μM. The uptake experiments were performed under standard growth conditions (temperature 28 ± 1°C and light intensity 70 μE m2 s1), and samples (0.5 ml) were withdrawn at 10 min of incubation. Uptake of nitrate was determined by measuring the depletion of nitrate in the uptake mixture [18]. Since the uptake was linear up to 10 min, the initial uptake rate was determined as to nitrate taken up by the cells during the first 10 min. The assay was performed in triplicate with independent cultures.
2.2 Bioinformatics analysis
The amino acid sequences for the nitrate transporter proteins NrtA (accession no. BAB75032.1), NrtC (accession no. WP_010994786.1) and NrtD (accession no. WP_010994787) of Nostoc (Anabaena) 7120 were retrieved from National Centre for Biotechnology Information (NCBI) protein database, PDB (www.rcsb.org; http://www.ncbi.nlm.nih.gov/protein). The most suitable template for each protein was identified using BLASTP (protein BLAST, http://blast.ncbi.nlm.nih.gov/Blast.cgi) and selected for further bioinformatics approach. Protein domains, families and functional sites, patterns, and profile search were performed using PROSITE server [http://prosite.expasy.org/] [19–21]. The retrieved protein sequences were used for multiple sequence alignment using ClustalW (http://www.ebi.ac.uk/Tools/msa/clustalw2/). Sequences were aligned across the entire length or only in certain regions [22]. The top 100 sequences of each nrtA, nrtC and nrtD were taken for the phylogenetic analysis using PHYLIP (Phylogeny Inference Package) programs Seqboot, Protdist, Neighbour, and Consense to generate a phylogenetic tree [23–25]. Mega 5.0 and Tree view were used to view the phylograms [26].
2.3 Homology modeling
The nrtACD proteins of Nostoc PCC 7120 were used for template search on the NCBI PDB database for structure prediction [27, 28]. The selected templates were used for homology modeling employing Easy Modeller 4.0 [29]. It is a standalone tool with an intuitive interface that clearly defines the different steps of homology modeling for comparable protein structures based on MODELLER [30, 31]. Since no models were available for Anabaena 7120 nrtACD proteins in the PDB database, an approach was made for constructing their 3D structures using homology modeling. The validation of a generated model of nrtACD proteins was evaluated for stereochemical quality using assessment of Ramachandran Plot (RAMPAGE) server (http://mordred.boc.cam.ac.uk/~rapper/rampage.php), PROCHEK [32] and PDBSum [33–35]. The quality of the predicted model was calculated by the SAVES server (http://services.mbi.ucla.edu/SAVES/). Statistical analysis of non-bonded interactions among different atom types of the generated model was estimated using ERRAT server [36, 37], while compatibility of an atomic model (3D) with its own amino acid sequence (1D) was tested by Verify3D [37]. Potential errors in the 3D model and Z-score values were calculated and compared with a template by the ProSA-web server, which ensured overall quality as well as determined that the input structure is within the range of scores typically found for the native proteins of similar size [38]. After complete evaluation and visualization (Discovery Studio 3.0) sequence of the model is then deposited to PMDB (Protein Model Data Base) [39].
2.4 Molecular docking analysis
The chemical compounds GB (CID: 247) and nitrate (CID: 943) were retrieved using Pubchem compound database as SDF file [40] and converted to PDB using Discovery Studio Visualizer. The active sites and their residues for selected structures were identified using the MetaPocket server. The server identifies pockets on protein surfaces to predict ligand-binding sites [41, 42]. The retrieved PDB files for GB and nitrate were imported to the Patchdock server program (http://www.molegro.com/products.php) for docking with transporter proteins nrtACD models [43–45].
3. Results
Salinity (0.1 M NaCl) reduced the rate of nitrate transport into the wild-type cells of Anabaena 7120. In contrast, under identical salinity, transport rate increased into the transgenic cells (transformed with ApGSMT-DMT) that synthesized and accumulated GB (Fig 2A). Determination of kinetics of the nitrate transport (Fig 2B) showed that salinity increased the Ks value (32.3μM) in WT, while the value was lowered (Ks 17.8 μM) in ApGSMT-DMT transformed cells compared to control (Ks = 19.5 μM), i.e., transport in absence of NaCl. This indicated the increased affinity of GB synthesizing cells towards nitrate [6]. The Vmax value (26.8 μmol mg protein-1 min-1) was the highest for GB synthesizing cells, followed by control (22.6 μmol mg protein-1 min-1), and was the lowest (21.9 μmol mg protein-1 min-1) for NaCl-exposed cells. This revealed that salinization of the nutrient solution reduced the efficiency of the nitrate transport system, while GB stimulated it. Possibly, recombinant GB interacted positively with nitrate transporter proteins resulting in efficient nitrate transport. To unveil the mechanism of their interactions, we performed in silico analyses.
Exponential phase cyanobacterial cells growing on N2 in the presence and absence of 0.1 M NaCl were collected and suspended in their specific media containing varied concentrations (40–500 μ) of nitrate. Uptake was performed under standard growth conditions in light (70 μE m2 s1). Samples were withdrawn at 10 min of incubation and measured nitrate depletion in the incubation mixture. Values represent means and SDs of triplicate experiments. (A) Initial rate of nitrate transport at different substrate concentrations; (B) Hofstee plot (v versus v/[S]) to determine the kinetics.
3.1. Sequence retrieval and analysis
Prosite analysis (S1 Fig) revealed that nrtA protein was 427 amino acids long and contained a Tat (twin-arginine translocation) signal profile (PS51318; score 8.195) of 34 amino acids (1–34) necessary to initiate the transport of folded proteins [46]. The Tat signal peptide consisted of three domains, the positively charged N-terminal domain, hydrophobic domain, and C-terminal domain. The majority of the amino acids were hydrophobic in nature.
Prosite results for nrtC protein revealed that it contains 657 amino acids (S2 Fig) and two functional motifs, ABC transporter 2 (5–239 amino acids) and ABC transporter 1 (139–153 amino acids) with Walker A (GHSGCGKS, 42–49) and Walker B (LLLD, 160–163), which is responsible for ATP binding and hydrolysis. Latter, it was preceded by an alpha-helical sub-domain and a signature motif (LSGGMKQR, 139–146). Sequence alignment results revealed that the nrtC subunit contained two distinct domains; one was an N-terminal domain (5–239 amino acids) mostly composed of hydrophobic residues, and the other C-terminal domain (268–591 amino acids) possessing mostly hydrophobic amino acids. N-terminal domain (5–228 amino acids) of nrtC showed 59% identity and 77% similarity with nrtD (the ATP-binding subunit), whereas the C-terminal domain (268–591 amino acids) was 28.2% identical and 50% similar to nrtA, suggesting that the C-terminal domain of nrtC was a solute binding domain and N-terminal as ATP binding.
NrtD protein was found 277 amino acids long and contained ATP binding motif (PS50893; score 26.484) located at N-terminal (21–254 amino acids) (S3 Fig). The majority of the amino acids were hydrophobic.
3.2. Sequence alignment and phylogenetic analysis
The nrtACD protein sequences of Nostoc 7120 (target organism) were taken as input to find out the similar protein sequences in different cyanobacterial species using BLASTP. The searche against non-redundant data base and protein data bank resulted in 103, 86 and 106 homologs for nrtA, nrtC and nrtD, respectively (S1–S3 Tables), and were used for multiple sequence alignment (S4–S6 Figs). The phylogenetic analysis of nrtA, nrtC and nrtD revealed different clustering with respect to cyanobacterial nrt proteins. The nrtA cladogram (Fig 3) shows two clades in which clade 1 comprises two subclades and clade 2 comprises seven small subclades. The two subclades of clade 1 are differentiated with one having Anabaena and Nostoc species and the other containing the Microcystis and Gemnocyctis species. The subclades one, two, three and six of clad 2 showed Agrobacterium and Rhizobium, Bradyrhizobium species, Thioalkalivibrio and Psuedomonas, Thermosynechococcus and Leptolyngbya as dominating species in the phylogenetic tree, respectively. The nrtC phyllogram (Fig 4) also has two clades each having two subclades comprising of Anabaena, Nostoc, Fischerella, Planktothrix, Leptolyngbya and Microcystis species. Phylogeny of nrtD protein showed eight clades (Fig 5), where clade seven had six subclades. The phylogenetic tree was comprised of Nostoc, Anabaena, Fischerella and Calothrix species. Clade one was dominated by the presence of Nostoc and Anabaena species.
The tree was constructed by PHYLIP programs using Seqboot, Protdist, Neighbour and Consense methods. Bootstraping (1000 times) was performed to obtain support values for each branch. Values are shown at the nodes of each branch point.
The tree was constructed by PHYLIP programs using Seqboot, Protdist, Neighbour and Consense methods. Bootstraping (1000 times) was performed to obtain support values for each branch. Values are shown at the nodes of each branch point.
The tree was constructed by PHYLIP programs using Seqboot, Protdist, Neighbour and Consense methods. Bootstraping (1000 times) was performed to obtain support values for each branch. Values are shown at the nodes of each branch point.
Species like Hapalosiphon sp., Mastigocladus laminosus and Calothrix sp. had close homologues of nrtA protein, Tolipothrix sp. PCC of both the nrtC and nrtD proteins, while Fortiacontorta species appeared to have a close homologue of all three nrtACD proteins of Anabaena. Overall phylogenetic results thus indicated that the nitrate-ABC-transporter proteins (nrtACD) homology is highly conserved among the cyanobacterial species.
3.3. Modeling of nrtACD proteins of Nostoc sp. PCC 7120 and validation of 3D homology
The target proteins (nrtACD) of Nostoc sp. PCC 7120 were used for 3D structure prediction using templates (Fig 6A–6C). Details of similarity for each selected template are presented (Table 1). The results obtained by RAMPAGE for stereochemical quality of nrtA showed that 368 residues (96.1%) lies in the favored region, 15 residues (3.9%) lies in the allowed and none in the outliner region (Fig 7A). For nrtC protein, 335 residues (93.8%) lies in the favored region, 18 residues (5.0%) lies in the allowed region and 4 residues (1.1%) lies in the outliner region (Fig 7B), and for nrtD, 126 residues (94.7%) lies in the favored region, 6 residues (4.5%) lies in the allowed region and 1 residue (0.8%) lies in the outliner region (Fig 7C). PROCHEK output in the form of Ramachandran plot for nrtA showed 304 residues (91.8%) lies in the most favored region, 24 residues (7.3%) lies in the additional allowed region, 3 residues (0.9%) lies in generously allowed regions and none in the disallowed region (Fig 8A). For nrtC protein, 281 residues (90.6%) lies in the most favored region, 23 residues (7.4%) lies in the additional allowed region, 5 residues (1.6%) lies in generously allowed regions and 1 residue (0.3%) lies in the disallowed region (Fig 8B), whereas for nrtD, the 102 residues (90.3%) lies in the most favored region, 10 residues (8.8%) lies in the additional allowed region, 1 residues (0.9%) lies in generously allowed regions and none in the disallowed region (Fig 8C). ERRAT web server was used to check the overall quality factor of nrtACD proteins which was found 73.21 in nrtA (Fig 9), 68.58 in nrtC (Fig 10) and 66.94% in nrtD (Fig 11) protein models, while the Verify 3D was used to determine the compatibility of atomic model (3D) with its own amino acid sequence (1D) by assigning a structural class which is based on its position and environment (alpha, beta, loop, polar and nonpolar), and the compatibility was found more than or equal to 0.2 averaged 3D-1D score for the 95.84, 95.26 and 89.92% amino acid residue of NrtA, NrtC and NrtD, respectively. Hence, more than 80% amino acid score with the template protein demonstrated the good quality of the model proteins. The overall model quality was checked by the ProSA-web tool [35] by determining their Z-score values. Z-score value for nrtA, NrtC and NrtD was found -9.68 (Fig 12A), -8.57 (Fig 12B) and -5.2 (Fig 12C), respectively. The results taken together suggested a close similarity between the template and the predicted structures and reinforced a good quality of the modeled protein.
Models are presented as ribbon structure. Alpha helices are in red, beta sheets in cyan, turn in grey and coil in green color. (A) nrtA, (B) nrtC and (C) nrtD.
Values indicate number of the residues in favored, allowed, and outlier region. (A) nrtA, (B) nrtC and (C) nrtD.
Red regions in the graph indicate the most allowed regions and yellow the allowed region. (A) nrtA, (B) nrtC and (C) nrtD.
Values indicate over all structure quality of the traget protein.
Values indicate over all structure quality of the traget protein.
Values indicate over all structure quality of the traget protein.
Z-scores of the protein chains in PDB were determined by X-ray crystallography (light blue) and NMR spectroscopy (dark blue). The proSA web results indicated that protein structures had characteristic features of their native structures. The Z-score for all the three target proteins was highlighted as a large dot.
3.4. Molecular docking
Metapocket server predicted three possible active sites in nrtACD proteins (Table 2). Active site one in the nrtA protein (Fig 13A) and active site two in nrtC (Fig 14A) and nrtD proteins (Fig 15A) appeared as the prominent binding sites (Table 2). The polar hydrophilic and positively charged amino acids Lys251, Lys295and Arg298, neutral Ser290, and negatively charged Asp246, and hydrophobic neutral amino acid Pro48 of nrtA protein interacted in docking with nitrate (Fig 13B). The majority of the molecular interactions were hydrophilic with H-bonding between Arg298:NH1-O2, and Arg298:NH1-O2. On the other, the majority of the molecular interactions of nrtA with GB were hydrophobic. Five aliphatic hydrophobic neutral amino acid residues (Leu11, Leu66, Val13, Val185 and Val186) along with hydrophobic neutral (Pro136), polar hydrophobic neutral (Gln97), polar hydrophilic positively charged (Lys215), aliphatic neutral (Gly142) and aromatic hydrophobic neutral (Trp44) amino acid residues interacted in the absence of H-bond (Fig 13B).
(A) Green color in the sphere indicates prominent active site with which the ligand interacted. (B) 3D level interaction.
(A) Poses of docked complexes; green color in the sphere indicates prominent active site on which the ligand interacted. (B) 3D level interaction.
(A) Poses of docked complexes; green color in sphere indicates prominent active site, on which the ligand interacted. (B) 3D level interaction.
In nrtC/nitrate docking, the majority of the amino acids involved were hydrophobic and interacted in the presence of H-bond (Fig 14B). Polar hydrophilic negatively charged (sp158, aromatic hydrophobic neutral Phe81, polar hydrophobic neutral Ser82 and aliphatic hydrophobic neutral Leu85 amino acids were buried in the center with H-bond between Ser82:N-Asp158:O, and Leu85:N-Ser82:O. In the interaction of GB with nrtC, the bulk of the amino acids were hydrophobic. Aliphatic hydrophobic neutral (Leu23, 41, 4, 21, 194, Val21, Ala43, Ala44 and Ile22), aromatic hydrophobic neutral (Phe13), and polar hydrophilic charged (Asp42) and polar hydrophilic neutral (Gln45) amino acid residues bind with GB (Fig 14B). H-bond was present between Ala43:N-Asp42:OD1, Ala44:N-Gly195:O and Leu85:N-Lys83:O amino acid residues.
The results of nrtD/nitrate docking showed that the major portion of the amino acid residues that interacted with nitrate was hydrophobic and contains aliphatic hydrophobic neutral Leu17, Leu5 and Val3, polar hydrophobic neutral Cys4 and aromatic hydrophobic neutral Phe2 with hydrogen bond in between Cys4: SG-O2 (Fig 15B). Similarly, the majority of nrtD amino acids interacting with GB were hydrophobic involved aliphatic hydrophobic neutral (Leu127, Ala13 and Ile114), polar hydrophobic neutral (Asn51), aromatic hydrophobic neutral (Phe131 and Pro130) and polar hydrophilic (Gln110) with H-bond between (Tyr52:N-:Asn51:OD1, Ala113:N-:GLN110:O and Ile114:N-:Arg111:O) (Fig 15B). PatchDock server docking calculation for the ligands GB and nitrate with the receptor proteins nrtACD presented a higher score solution, larger acquired docking area and lesser atomic contact energy (ACE) for GB compared to nitrate (Table 3), revealing a higher affinity of GB towards receptor proteins compared to nitrate.
4. Discussion
Nitrate (NO3−) is the main nitrogen source for plants and their mobilization is regulated by various transporters in the plant body [47]. To transport nitrate into the cell, nrtA protein of ABC-transporter binds to nitrate and delivers it to the nrtB (Fig 1). NrtA has a high-affinity with the nitrate ABC transporter. The passage of nitrate transportation through the pore is linked to ATP binding and its consequent hydrolysis on the P-loop of the ABC domain. This brings conformational changes in the transmembrane domains and allows nitrate to cross the lipid bilayer of the membrane. Short cytoplasmic coupling helices provide critical contacts to TMD with ATP domain [4, 5, 48]. When nitrate concentration becomes high in the cell, nrtC binds to nitrate allosterically and inhibits nitrate transport [49, 50].
Salinity reduced the nitrate transport into WT cells (Fig 2), while the rate of transport was increased into the transgenic cells (synthesizing recombinant GB). This means that in addition to protecting nitrate transport from salt toxicity, recombinant GB increased the efficiency of nitrate transport system. This increased efficiency might be due to the interaction of GB with the nitrate transporter proteins. NaCl stress increases Na+ concentration in the cells, which disrupts the water structure and reduces hydrophobic interactions and hydrostatic forces within proteins [51]. Hence, NaCl-induced inhibition of nitrate transport is due to the disturbed hydrophobic–electrostatic balance necessary to maintain the protein structure [52, 53]. Compatible solute GB can interact with both the hydrophilic and hydrophobic regions of macromolecules [54] and can replace water at the surface of proteins [55] enabling the Nrt proteins to stabilize their native structure.
In silico analyses revealed that the deduced periplasmic protein nrtA in Anabaena 7120 that traps nitrate of the ambient solution was 427 amino acids long, the majority of which were hydrophobic as reported for Synechococcus sp. PCC7942 [56]. The molecular interaction of nrtA with nitrate (Fig 13B) was hydrophilic (involving Lys251, Lys295and Arg298, Ser290, Asp246, and Pro48 with H-bonding at Arg298). The protein contained a Tat signal profile of 34 amino acids (1–34) necessary to initiate the transport [46]. Tat protein transport system has been reported in most bacteria and archaea. In cyanobacteria, Tat consensus motifs have been reported in deduced protein PhoD of Aphanothece halophytica exhibiting alkaline phosphatase (assimilate phosphorous under salt stress) and phosphodiesterase activities [57]. The presence of Tat motifs also supports the capturing of nitrate by pushing nrtA off the membrane and passing to the nrtB [4]. Since one nrtA protein binds to one molecule of nitrate [58], nrtA becomes the most abundant protein in the nitrate-grown cells [59].
The N-terminal (5–239 amino acids) and C-terminal (268–591 amino acids) of nrtC protein of Anabaena 7120 contained 657 amino acids that are mostly hydrophobic. Therefore, the interaction of nrtC with nitrate was also found hydrophobic. The hydrophilic negatively charged amino acid residues (Asp158, Phe81, Ser82 and Leu85) were buried in the center with H-bond between Ser82:N-Asp158:O, and Leu85:N-Ser82:O (Fig 14B). The nrtC protein in Synechococcus sp. PCC 7942, is reported to be 659 amino acids long with the deduced amino acid sequences 1–254 for the N-terminal domain, and 279–659 for the C-terminal domain [56]. In Phormidium laminosum, the N-terminal (1–262) contained mostly hydrophobic residues, while the C-terminal (281–626) was highly hydrophilic [60]. In Anabaena 7120, the N-terminal of nrtC contained two functional motifs, ABC transporter 2 (5–239 amino acids) and ABC transporter 1 (139–153 amino acids) consisting of Walker A (GHSGCGKS, 42–49) and Walker B (LLLD, 160–163) similar to that reported in Phormidium laminosum [58], and was responsible for ATP binding and hydrolysis.
The C-terminal domain (268–591 amino acids) of nrtC in Anabaena 7120 was 28.2% identical and 50% similar to nrtA. In Synechocystis sp. PCC 6803, C-terminal domain was >30% identical and 50% similar to nrtA [4], while in Synechococcus sp. PCC 7942, the identity of C-terminal domain of nrtC (amino acid residues 279–659) to nrtA was 30% [3, 56]. This further affirmed that the C-terminal domain of nrtC protein is a solute-binding domain. The deduced protein NrtD of Anabaena 7120 was of 277 amino acids, the majority of which were hydrophobic. The interaction of nrtD with nitrate was found hydrophobic and included Leu17, Leu5 Val3, Cys4, and Phe2 with hydrogen bond between Cys4: SG-O2 (Fig 15B). The N-terminal domain (5–228 amino acids) showed 59% identity and 77% similarity with nrtD (the ATP-binding subunit). At N-terminal 21–254 amino acid sequences were ATP binding motif. In Synechococcus sp. PCC7942, the nrtD protein consisted of 274 amino acids and showed 58% identity to the N-terminal domain of nrtC [3, 56]. This protein domain of the ABC transporter is responsible to bind and hydrolyze ATP to energize the nitrate transport system.
Ammonium is preferred over nitrate when available in the solution and inhibits nitrate uptake [49, 61]. When intracellular nitrate level goes up, the C-terminal domain of nrtC, which is structurally similar to nrtA, may bind to nitrate allosterically as an effector-binding domain. Hence, nrtC role is suggested as a negative regulator involved in the ammonium-promoted inhibition of nitrate transport activity [49, 62]. The interactions of GB with deduced nrtACD proteins of Anabaena PCC 7120 were hydrophobic and showed higher affinity towards GB compared to nitrate. Further, docking of GB was more efficient (with maximum docking score) with nrtC than that of nrtA and nrtD. Among the osmoprotectants, GB is by far the most effective and the most characterized osmolyte to date [63], allowing the cells to grow luxuriantly under salt stress. The level of GB in transgenic Anabaena 7120 was low (7.5 μmol g−1 dry weight) and insufficient to contribute significantly to cellular osmotic potential, yet GB conferred considerable tolerance to salt-stress [6]. Other than the osmoprotective role, the mechanisms put forth to explain the GB enhanced stress tolerance include: (1) stabilization of proteins native structure and maintenance of membrane integrity, and (2) a regulatory function in the expression of genes i.e., induction of specific genes that led to increased stress tolerance [64, 65], including photosynthetic and respiratory electron transport [6]. It suggested that both the roles of GB, stabilization of protein native structure and regulatory function in gene expression are responsible for the efficient nitrate transport in GB synthesizing and accumulating transgenic Anabaena PCC 7120.
Salinity reduces the thermodynamic activity of water, thereby displacing the water of hydration from the protein surfaces [66] causing protein unfolding. GB induces the folding of unfolded (inactive) proteins in vitro and in vivo in a chaperone-like manner by decreasing the exposure of protein surface to solvent by drastically decreasing the water activity [66]. Experimentation regarding in situ evidence of the GB interaction with the proteins is difficult, hence scanty [67]. However, the periplasmic protein ProX of the ATP-transport system and ProU of Escherichia coli is reported to bind GB with high affinity [68, 69]. The receptor protein ProX undergoes a conformational change to form a site/cavity to accommodate GB and has an evenly negative surface potential [70]. The two conserved cis-proline align the Trp65 and Trp188 residues in a way to form an almost rectangular box in the polypeptide backbone of tryptophan residues [70, 71]. Possibly, nitrate transporter proteins nrtACD bind with GB in a similar fashion resulting in efficient nitrate transport. To our knowledge, this is the first report of GB interaction with ABC-nitrate transporter proteins and may open the way for further experimentation.
5. Conclusion
Salinity is one of the major factors limiting plant metabolic activities, growth, biomass and yield. GB is a universal osmoprotectant synthesized by prokaryotes and eukaryotes, including some crop plants to counter various stresses including salinity, and is also beneficial to humans. Nitrogen is the most limiting nutrient for plant growth and is often available as nitrate-N in the natural ecosystems, though at μM concentrations. The cyanobacterium Anabaena 7120 transported nitrate into the cells through an active process mediated by ABC-transporter proteins nrtABCD. The ApGSMT-DMT transformed cells possessed Ń-methyltransferase activities and de novo synthesized and accumulated intracellular recombinant GB. Salinity reduced the transport of nitrate into the cells. The recombinant GB not only protected the nitrate transport into the cells but increased the rate of transport under salinity. The low amount of intracellular GB (7.5 μmol g−1 dry weight) was insufficient to act as a cellular osmoticum. It is noted that the increased rate of nitrate transport in GB accumulating cells was due to the interaction of GB with ABC-nitrate transporter proteins and also by the participation of GB in the regulation of other intracellular metabolic processes. Thus, the results provided important bearing on the mechanism of GB interaction with ABC-nitrate transporter proteins.
Supporting information
S1 Fig. Prosite analysis of nrtA-nitrate transporter protein.
nrtA protein showing functional motif Tat signal profile (PS51318).
https://doi.org/10.1371/journal.pone.0257870.s001
(TIF)
S2 Fig. Prosite analysis of nrtC-nitrate transporter protein.
nrtC protein showing two functional motifs, ABC transporter 2 (5–239 amino acids) and ABC transporter 1 (139–153 amino acids) consisting Walker A (GHSGCGKS, 42–49) and Walker B (LLLD, 160–163) responsible for ATP binding and hydrolysis.
https://doi.org/10.1371/journal.pone.0257870.s002
(TIF)
S3 Fig. Prosite analysis of nrtD-nitrate transporter protein.
nrtD protein showing ATP binding motif (PS50893).
https://doi.org/10.1371/journal.pone.0257870.s003
(TIF)
S4 Fig. Multiple sequence alignment of target protein nrtA of Nostoc sp. PCC 7120.
https://doi.org/10.1371/journal.pone.0257870.s004
(PDF)
S5 Fig. Multiple sequence alignment of target protein nrtC of Nostoc sp. PCC 7120.
https://doi.org/10.1371/journal.pone.0257870.s005
(PDF)
S6 Fig. Multiple sequence alignment of target protein nrtD of Nostoc sp. PCC 7120.
https://doi.org/10.1371/journal.pone.0257870.s006
(PDF)
S1 Table. Sequence retrieval details of nrtA protein for multiple sequence alignment and phylogenetic tree construction showing accession number, proteins, organisms and name proposed.
https://doi.org/10.1371/journal.pone.0257870.s007
(DOC)
S2 Table. Sequence retrieval details of nrtC protein for multiple sequence alignment and phylogenetic tree construction showing accession number, proteins, organisms and name proposed.
https://doi.org/10.1371/journal.pone.0257870.s008
(DOC)
S3 Table. Sequence retrieval details of nrtD protein for multiple sequence alignment and phylogenetic tree construction showing accession number, proteins, organisms and name proposed.
https://doi.org/10.1371/journal.pone.0257870.s009
(DOC)
Acknowledgments
Authors acknowledge the groups whose literatures have been cited to write this manuscript.
References
- 1. Knoll AH. Cyanobacteria and earth history. The Cyanobacteria: Molecular Biology, Genomics, and Evolution. 2008; 484.
- 2. Rai AK, Kashyap AK, Gupta SL. ATP-dependent uptake of nitrate in Nostoc muscorum and inhibition by ammonium ions. Biochim Biophys Acta (BBA)-Gen Sub. 1981; 674(1):78–86. pmid:6786358
- 3. Omata T. Structure, function and regulation of the nitrate transport system of the cyanobacterium Synechococcus sp. PCC7942. Plant Cell Physiol. 1995; 36(2):207–213. pmid:7767600
- 4. Koropatkin NM, Pakrasi HB, Smith TJ. Atomic structure of a nitrate-binding protein crucial for photosynthetic productivity. Proc Natl Acad Sci USA. 2006; 103(26):9820–9825. pmid:16777960
- 5. Saurin W, Hofnung M, Dassa E. Getting in or out: Early segregation between importers and exporters in the evolution of ATP-binding cassette (ABC) transporters. J Mol Evol. 1999; 48(1):22–41. pmid:9873074
- 6. Swapnil P, Singh M, Singh S, Sharma NK, Rai AK. Recombinant glycinebetaine improves metabolic activities, ionic balance and salt tolerance in diazotrophic freshwater cyanobacteria. Algal Res. 2015; 11:194–203.
- 7. Flores E, Frías JE, Rubio LM, Herrero A. Photosynthetic nitrate assimilation in cyanobacteria. Photosynth Res. 2005; 83(2):117–133. pmid:16143847
- 8. Swapnil P, Yadav AK, Srivastav S, Sharma NK, Srikrishna S, Rai AK. Biphasic ROS accumulation and programmed cell death in a cyanobacterium exposed to salinity (NaCl and Na2SO4). Algal Res. 2017; 23:88–95.
- 9. Rodríguez R, Lara C, Guerrero MG. Nitrate transport in the cyanobacterium Anacystis nidulans R2. Kinetic and energetic aspects. Biochem J. 1992; 282(3):639–643. pmid:1554347
- 10. Rodríguez R, Guerrero MG, Lara C. Mechanism of sodium/nitrate symport in Anacystis nidulans R2. Biochim Biophys Acta (BBA)-Bioenerg. 1994; 1187(2):250–254.
- 11. Rai AK, Tiwari SP. Mutants of the cyanobacterium Anabaena sp. PCC 7120 altered in nitrate transport and reduction. Curr Microbiol. 1999; 39(5):237–243. pmid:10489430
- 12. Rai AK, Abraham G. Relationship of combined nitrogen sources to salt tolerance in freshwater cyanobacterium Anabaena doliolum. J Appl Microbiol. 1995; 78(5):501–506.
- 13. Rai AK, Tiwari SP. Response to NaCl of nitrate assimilation and nitrogenase activity of the cyanobacterium Anabaena sp. PCC 7120 and its mutants. J Appl Microbiol. 1999; 87(6):877–883. pmid:10692073
- 14. Kumari J, Rathore MS. Na+/K+-ATPase a primary membrane transporter: An overview and recent advances with special reference to algae. J Membr Biol. 2020; 253:191–204. pmid:32430620
- 15. Yancey PH. Organic osmolytes as compatible, metabolic and counteracting cytoprotectants in high osmolarity and other stresses. J Exp Biol. 2005; 208(15):2819–2830. pmid:16043587
- 16. Singh M, Sharma NK, Prasad SB, Yadav SS, Narayan G, Rai AK. The freshwater cyanobacterium Anabaena doliolum transformed with ApGSMT-DMT exhibited enhanced salt tolerance and protection to nitrogenase activity, but became halophilic. Microbiology. 2013; 159(3):641–648. pmid:23329680
- 17. Rippka R, Deruelles J, Waterbury JB, Herdman M, Stanier RY. Generic assignments, strain histories and properties of pure cultures of cyanobacteria. Microbiology. 1979; 111(1):1–61.
- 18. Cawse P. The determination of nitrate in soil solutions by ultraviolet spectrophotometry. Analyst. 1967; 92(1094):311–315.
- 19. Sigrist CJ, Cerutti L, Hulo N, Gattiker A, Falquet L, Pagni M, et al. PROSITE: a documented database using patterns and profiles as motif descriptors. Brief Bioinform. 2002; 3(3):265–274. pmid:12230035
- 20. Hulo N, Bairoch A, Bulliard V, Cerutti L, De Castro E, Langendijk-Genevaux PS, et al. The PROSITE database. Nucleic Acids Res. 2006; 34:D227–D230. pmid:16381852
- 21. Meena M, Aamir M, Vikas K, Swapnil P, Upadhyay RS. Evaluation of morpho-physiological growth parameters of tomato in response to Cd induced toxicity and characterization of metal sensitive NRAMP3 transporter protein. Environ Exp Bot. 2018; 148:144–67.
- 22. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007; 23(21):2947–2948. pmid:17846036
- 23. Felsenstein J. PHYLIP—Phylogeny inference package (Version 3.2). Cladistics. 1989; 5:164–166.
- 24. Bhaduri A, Misra R, Maji A, Bhetaria PJ, Mishra S, Arora G, et al. Mycobacterium tuberculosis cyclophilin A uses novel signal sequence for secretion and mimics eukaryotic cyclophilins for interaction with host protein repertoire. PLoS One. 2014; 9(2):e88090. pmid:24505389
- 25. Bhetaria PJ, Prajapathi M, Bhaduri A, Mandal RS, Varma A, Madan T, et al. Phylogenetic and structural analysis of polyketide synthases in Aspergilli. Evol Bioinfo. 2016; 12:101–19.
- 26. Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007; 24(8):1596–1599. pmid:17488738
- 27.
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The Protein Data Bank, 1999. In international tables for crystallography volume F: Crystallography of biological macromolecules. Springer, Netherlands. 2006; 675–684.
- 28. Sharma D, Pathak M, Sharma R, Tyagi P, Chawla R, Basu M, et al. Homology modeling and docking studies of VP24 protein of Ebola virus with an antiviral drug and its derivatives. Chem Biol Let. 2017; 4(1):27–32.
- 29. Kuntal BK, Aparoy P, Reddanna P. EasyModeller: A graphical interface to MODELLER. BMC Res Notes. 2010; 3(1):226. pmid:20712861
- 30. Šali A, Blundell TL. Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993; 234(3):779–815. pmid:8254673
- 31. Fiser A, Do RKG. Modeling of loops in protein structures. Protein Sci. 2000; 9(9):1753–1773. pmid:11045621
- 32. Laskowski RA, Rullmann JAC, MacArthur MW, Kaptein R, Thornton JM. AQUA and PROCHECK-NMR: programs for checking the quality of protein structures solved by NMR. J Biomol NMR. 1996; 8(4):477–486. pmid:9008363
- 33. Lovell SC, Word JM, Richardson JS, Richardson DC. The penultimate rotamer library. Proteins: Struct Funct Bioinf. 2000; 40(3):389–408.
- 34. Laskowski RA, Hutchinson EG, Michie AD, Wallace AC, Jones ML, Thornton JM. PDBsum: a Web-based database of summaries and analyses of all PDB structures. Trends Biochem Sci. 1997; 22(12): 488–490. pmid:9433130
- 35. Chowdhury T, Handique PJ. Homology modelling and in-silico analysis of multidrug resistant protein, PMPM of Pseudomonas aeruginosa. Biores Bull. 2014; 3:1–6.
- 36. Colovos C, Yeates TO. Verification of protein structures: patterns of nonbonded atomic interactions. Protein Sci. 1993; 2(9):1511–1519. pmid:8401235
- 37. Luthy R, Bowie JU, Eisenberg D. Assessment of protein models with three-dimensional profiles. Nature. 1992; 356(6364):83–85. pmid:1538787
- 38. Wiederstein M, Sippl MJ. ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 2007; 35:W407–W410. pmid:17517781
- 39. Castrignanò T, De Meo PDO, Cozzetto D, Talamo IG, Tramontano A. The PMDB protein model database. Nucleic Acids Res. 2006; 34:D306–D309. pmid:16381873
- 40. Kim S, Thiessen PA, Bolton EE, Chen J, Fu G, Gindulyte A, et al. PubChem substance and compound databases. Nucleic Acids Res. 2015; 44(D1): D1202–D1213. pmid:26400175
- 41. Huang B. MetaPocket: A meta approach to improve protein ligand binding site prediction. Omi A J Integr Biol. 2009; 13(4):325–330. pmid:19645590
- 42. Aamir M, Singh VK, Dubey MK, Meena M, Kashyap SP, Katari SK, et al. In silico prediction, characterization, molecular docking, and dynamic studies on fungal SDRs as novel targets for searching potential fungicides against Fusarium wilt in tomato. Front Pharmacol. 2018; 9:1038. pmid:30405403
- 43. Schneidman-Duhovny D, Inbar Y, Nussinov R, Wolfson HJ. PatchDock and SymmDock: servers for rigid and symmetric docking. Nucleic Acids Res. 2005; 33:W363–W367. pmid:15980490
- 44. Aamir M, Singh VK, Meena M, Upadhyay RS, Gupta VK, Singh S. Structural and functional insights into WRKY3 and WRKY4 transcription factors to unravel the WRKY–DNA (W-box) complex interaction in tomato (Solanum lycopersicum) (L.). A computational approach. Front Plant Sci. 2017; 8:819. pmid:28611792
- 45. Marwal A, Meena M, Gaur RK. Molecular docking studies of coronavirus proteins with medicinal plant-based phytochemicals. Def Life Sci J. 2021; 6(1):57–63.
- 46. Wickner W, Schekman R. Protein translocation across biological membranes. Science. 2005; 310(5753):1452–1456. pmid:16322447
- 47. Shitan N, Yazaki K. New insights into the transport mechanisms in plant vacuoles. Int Rev Cell Mol Biol. 2013; 305:383–433. pmid:23890387
- 48. Locher KP. Mechanistic diversity in ATP-binding cassette (ABC) transporters. Nat Struct Mol Biol. 2016; 23(6):487–494. pmid:27273632
- 49. Kobayashi M, Rodríguez R, Lara C, Omata T. Involvement of the C-terminal domain of an ATP-binding subunit in the regulation of the ABC-type nitrate/nitrite transporter of the cyanobacterium Synechococcus sp. strain PCC 7942. J Biol Chem. 1997; 272(43):27197–27201. pmid:9341163
- 50. Kobayashi M, Takatani N, Tanigawa M, Omata T. Posttranslational regulation of nitrate assimilation in the cyanobacterium Synechocystis sp. strain PCC 6803. J Bacteriol. 2005; 187(2):498–506. pmid:15629921
- 51. Ghosh T, Kalra A, Garde S. On the salt-induced stabilization of pair and many-body hydrophobic interactions. J Phys Chem B. 2005; 109(1):642–651. pmid:16851057
- 52. Pollard A, Jones RW. Enzyme activities in concentrated solutions of glycinebetaine and other solutes. Planta. 1979; 144(3):291–298. pmid:24407261
- 53.
Wyn Jones RJ, Pollard A. Proteins, enzymes and inorganic ions. In Lauchli A, Pirson A (eds), Encyclopedia of Plant Physiology. Berlin, Springer. 1983; 528–562.
- 54. Sakamoto A, Murata N. The role of glycine betaine in the protection of plants from stress: clues from transgenic plants. Plant Cell Environ. 2002; 25(2):163–171. pmid:11841661
- 55. Hasegawa PM, Bressan RA, Zhu JK, Bohnert HJ. Plant cellular and molecular responses to high salinity. Annu Rev Plant Biol. 2000; 51(1):463–499. pmid:15012199
- 56. Omata T, Andriesse X, Hirano A. Identification and characterization of a gene cluster involved in nitrate transport in the cyanobacterium Synechococcus sp. PCC7942. Mol Gen Genet MGG. 1993; 236(2):193–202. pmid:8437564
- 57. Kageyama H, Tripathi K, Rai AK, Cha-Um S, Waditee-Sirisattha R, Takabe T. An alkaline phosphatase/phosphodiesterase, PhoD, induced by salt stress and secreted out of the cells of Aphanothece halophytica, a halotolerant cyanobacterium. Appl Environ Microbiol. 2011; 77:5178–5183. pmid:21666012
- 58. Maeda SI, Omata T. Substrate-binding lipoprotein of the cyanobacterium Synechococcus sp. strain PCC 7942 involved in the transport of nitrate and nitrite. J Biol Chem. 1997; 272(5): 3036–3041. pmid:9006953
- 59. Omata T, Ohmori M, Arai N, Ogawa T. Genetically engineered mutant of the cyanobacterium Synechococcus PCC 7942 defective in nitrate transport. Proc Natl Acad Sci USA. 1989; 86(17):6612–6616. pmid:16594065
- 60. Nagore D, Llarena M, Llama MJ, Serra JL. Characterization of the N-terminal domain of NrtC, the ATP-binding subunit of ABC-type nitrate transporter of the cyanobacterium Phormidium laminosum. Biochim Biophys Acta (BBA)-Gen Sub. 2003; 1623(2):143–153. pmid:14572912
- 61. Rai AK, Singh S. Regulation of nitrate uptake in Nostoc muscorum by glutamine synthetase. FEMS Microbiol Lett. 1982; 14:303–306.
- 62. Aichi M, Yoshihara S, Yamashita M, Maeda SI, Nagai K, Omata T. Characterization of the nitrate-nitrite transporter of the major facilitator superfamily (the nrtP gene product) from the cyanobacterium Nostoc punctiforme strain ATCC 29133. Biosci Biotechnol Biochem. 2006; 70(11):2682–2689. pmid:17090931
- 63. Cayley S, Record MT. Roles of cytoplasmic osmolytes, water, and crowding in the response of Escherichia coli to osmotic stress: biophysical basis of osmoprotection by glycine betaine. Biochemistry. 2003; 42(43):12596–12609. pmid:14580206
- 64. Einset J, Nielsen E, Connolly EL, Bones A, Sparstad T, Winge P, et al. Membrane‐trafficking RabA4c involved in the effect of glycine betaine on recovery from chilling stress in Arabidopsis. Physiol Plant. 2007; 130(4):511–518.
- 65. Einset J, Winge P, Bones AM, Connolly EL. The FRO2 ferric reductase is required for glycine betaine’s effect on chilling tolerance in Arabidopsis roots. Physiol Plant. 2008; 134(2):334–341. pmid:18513375
- 66. Record MT Jr, Courtenay ES, Cayley DS, Guttman HJ. Responses of E. coli to osmotic stress: large changes in amounts of cytoplasmic solutes and water. Trends Biochem Sci. 1998; 23(4):143–148. pmid:9584618
- 67. Bourot S, Sire O, Trautwetter A, Touzé T, Wu LF, Blanco C, et al. Glycine betaine-assisted protein folding in a lysA mutant of Escherichia coli. J Biol Chem. 2000; 275(2):1050–1056. pmid:10625645
- 68. Barron A, Jung JU, Villarejo M. Purification and characterization of a glycine betaine binding protein from Escherichia coli. J Biol Chem. 1987; 262(24):11841–11846. pmid:3305496
- 69. Haardt M, Kempf B, Faatz E, Bremer E. The osmoprotectant proline betaine is a major substrate for the binding-protein-dependent transport system ProU of Escherichia coli K-12. Mol Gene Genet. 1995; 246(6):783–796. pmid:7898450
- 70. Schiefner A, Breed J, Bösser L, Kneip S, Gade J, Holtmann G, et al. Cation-π interactions as determinants for binding of the compatible solutes glycine betaine and proline betaine by the periplasmic ligand-binding protein ProX from Escherichia coli. J Biol Chem. 2004; 279(7):5588–5596. pmid:14612446
- 71. Meena M, Swapnil P. Regulation of WRKY genes in plant defense with beneficial fungus Trichoderma: Current perspectives and future prospects. Arch Phytopathology Plant Protect. 2019; 52(1–2):1–17.