- 1Council of Scientific and Industrial Research (CSIR), Institute of Microbial Technology, Chandigarh, India
- 2Council of Scientific and Industrial Research (CSIR), Indian Institute of Chemical Biology, Kolkata, India
HapR is designated as a high cell density quorum sensing master regulatory protein of Vibrio cholerae. It is a member of the TetR family protein and functions both as an activator and a repressor by directly communicating with cognate promoters, thus controlling the expression of a plethora of genes in a density-dependent manner. Molecular insights reveal the domain architecture and further unveil the significance of a cross talk between the DNA binding domain and the dimerization domain for the functionality of the wild-type protein. The DNA binding domain is made up of three α-helices, where a helix-turn-helix motif spans between the helices α2 and α3. The essentiality of the glycine-rich linker linking helices α1 and α2 came into prominence while unraveling the molecular basis of a natural non-functional variant of HapR. Subsequently, the importance of linker length was demonstrated. The present study, involving a series of biochemical analyses coupled with molecular dynamics simulation, has illustrated the indispensability of a critical arginine within the linker at position 37 contributing to HapR–DNA binding activity.
Introduction
Communication is an essential characteristic of life. Even microbes, once perceived as organisms that maintain a reclusive lifestyle, also participate in intra- and interspecies social networking to shape their behavior in diverse environmental settings. Such communication of the microbial world known as “quorum-sensing” has become a major focus of research. Currently, a large body of cumulative knowledge strongly underpins how the quorum sensing communication process allows microbes to coordinate a plethora of cellular events including pathogenesis and survival in various ecological niches (Velicer, 2003).
Vibrio cholerae causes cholera, a life-threatening diarrheal disease with worldwide distribution (Camacho et al., 2018). In recent past, an elaborate study on the quorum sensing of V. cholerae unveiled its impact on the individual and collective behavior of this bacterium. Acting in concert, a myriad of factors including small metabolites, non-coding RNAs, and regulatory proteins control the optimal performance of the quorum sensing sensory network in V. cholerae (Lenz et al., 2004; Papenfort et al., 2017). Among the various regulatory proteins, HapR is crowned as the “high cell density master regulator.” Functionally, HapR belongs to the TetR family regulatory proteins, and it modulates a large number of disparate physiological events, thus governing virulence to environmental survival of the bacterium (Ng and Bassler, 2009). Crystal structure analysis reveals a two-domain architecture where the N-terminal DNA binding domain comprised three α-helices and the remaining six α-helices are found in the C-terminal-located dimerization domain (De Silva et al., 2007). Mechanistically, HapR communicates with cognate promoters by recognizing unique DNA binding motifs (Tsou et al., 2009), thereby fulfilling its commitment as a high cell density master regulator by regulating the constellation of diverse cellular activities. The motif 1 sequence comprised a dyad symmetry and a consensus “AATAR” (where R represents A or G). On the contrary, the motif 2 sequence is devoid of any dyad symmetry, but includes highly conserved consensus “TGT” (Tsou et al., 2009). Besides HapR, the crystal structure of SmcR, a HapR homolog of Vibrio vulnificus, has also revealed a similar architecture, having an N-terminal DNA binding domain and a C-terminal dimerization domain. Like HapR, SmcR also interacts with many promoters (Kim et al., 2010).
Previously, we have shown the contribution of a conserved glycine 39 located in the DNA binding domain in HapR function (Dongre et al., 2011). It is to be noted that all residues in the glycine linker are highly conserved in HapR homologs. In a continuing effort, we have further established that the length of this glycine-rich linker is also important for HapR function as a DNA binding protein (Singh et al., 2014). Recognition that the linker composition is also linked to protein function (Van Leeuwen et al., 1997; Seedorff and Schleif, 2011) further prompted us to embark on the present work.
Our present work focuses on the contribution of each residue of the glycine-rich linker in the DNA binding activity of HapR. Combined in vitro and in vivo results clearly identified a critical arginine at position 37 (R37) imparting in the HapR-mediated cross talk with various cognate promoters. Further molecular dynamics studies suggested the role of arginine in positioning the DNA binding domain and facilitating the interaction with cognate DNA.
Materials and Methods
Bacterial Strains and Media
The bacterial strains and plasmids used in this study are listed in Supplementary Table 1. V. cholerae strains were derived from a non-O1, non-O139 strain S7, and serogroup O37. In addition, V. cholerae strains V2S and GK178 were also included in this study. The strains were maintained at −70°C in Luria–Bertani (LB) medium containing 20% glycerol. Escherichia coli BL21 (DE3, Novagen) was used for the overexpression of proteins. All strains were propagated at 37°C in liquid with agitation or on solid (1.5% agar) in Luria broth, unless mentioned otherwise. For protease assay, V. cholerae recombinant strains were grown with aeration at 37°C in tryptic soya broth without dextrose (TSB-D). The growth medium was supplemented with streptomycin (100 μg ml–1), chloramphenicol (17 μg ml–1), and ampicillin (100 μg ml–1), whenever appropriate. All antibiotics were purchased from Sigma-Aldrich. Media ingredients were purchased from Himedia and Difco.
Site-Specific Mutagenesis
All the mutants were constructed using QuickChange PCR kit (Stratagene) as per the manufacturer’s guidelines. The primers are listed in Supplementary Table 2. Positive clones were sequenced in their entirety to confirm the clones and desired mutations at the corresponding positions. The desired constructs were further transformed into protease-negative strains of V. cholerae S7 and V2S.
Protease Assay
Protease activity was measured by employing an azocasein assay described earlier (Raychaudhuri et al., 2006). Briefly, wild-type and recombinant derivatives of V. cholerae strains V2S and S7 (Supplementary Table 1) were grown in TSB-D, containing appropriate antibiotics, with agitation to stationary phase at 37°C. Of the stationary phase culture supernatant, 100 μl was incubated with 100 μl azocasein (5 mg ml–1 in 100 mM Tris, pH 8.0) for 1 h at 37°C. The reaction was stopped by the addition of 400 μl of 10% trichloroacetic acid (TCA). After centrifugation, the supernatant was transferred to 700 μl of 525 mM NaOH and the optical density was determined at 442 nm. One azocasein unit was defined as the amount of enzyme producing an increase of 0.01 optical density (OD) units per hour.
Colony Morphology and Motility Assay
Colony corrugation was examined as per published protocols (Giglio et al., 2013). Briefly, overnight-grown cultures were diluted in fresh medium and grown to an OD600 = 0.3–0.4. Five microliters of the secondary culture was spotted on LB agar containing appropriate antibiotics. The plates were incubated at 30°C and photographs were taken after 48 h. For motility assays, all recombinant strains were grown overnight on LB agar, followed by the inoculation of a single colony by toothpick on plates composed of 1% tryptone, 0.5% NaCl, and 0.3% agar (Moisi et al., 2009). The plates were incubated at 37°C and photographs were taken after 12 h.
Hemagglutination Assay
Hemagglutinin assay was done as described elsewhere (Benitez et al., 2001). The cultures were grown overnight and the cells were pelleted to collect the supernatant. The chicken red blood cells collected from a slaughterhouse were washed twice with normal saline (0.85% NaCl) and once with Krebbs-Ringer-Tris (KRT) buffer (128 mM NaCl, 5.1 mM KCl, 1.34 mM MgSO4, 2.7 mM CaCl2, 10 mM Tris, pH 7.5). The red blood cells (RBCs) were resuspended in the KRT buffer at 1.2% (vol/vol). A serial dilution of 50 μl of the supernatant in the KRT buffer was mixed with 50 μl of freshly prepared chicken erythrocytes in the wells of a 96-well polystyrene V-bottom microtiter plate. The mixture was then incubated at room temperature for 30 min and the hemagglutination monitored visually.
Protein Purification
Wild-type and mutant proteins were purified by Ni-NTA chromatography. The genes encoding these proteins were cloned into the NdeI-BamHI site of the pET15b vector (Novagen) to generate N-terminal 6X His-HapR fusion protein. All positive clones were sequenced and subsequently transformed into E. coli BL21 (DE3). After induction with 0.4 mM isopropyl ß-D-1-thiogalactopyranoside (IPTG) at 37°C, HapR proteins were purified through Qiagen Ni2+-nitrilotriacetic acid columns. All proteins were dialyzed overnight in a solution of buffer A containing 10 mM Tris, pH 7.9, 300 mM KCl, and 0.1 mM EDTA.
Electrophoretic Gel Mobility Shift Assay With Promoter Regions of hapA, aphA, and vc0900
Gel mobility shift assay was done essentially as described earlier (Dongre et al., 2011; Singh et al., 2013). Briefly, the promoter regions of vca0865(hapA), vc2647(aphA), and vc0900(cdgG), respectively, were amplified with primer pairs as listed in Supplementary Table 2 from the genomic DNA of V. cholerae strains V2 and S7. The fragments were gel purified and end labeled with γ-dATP32 using T4 polynucleotide kinase (New England Biolabs® inc). The binding reaction was carried out with 4 ng of labeled fragments in 10 mM Tris–HCl, pH 7.9, 1 mM EDTA, 1 mM dithiothreitol (DTT), 60 mM KCl, 10% glycerol, 5 μg bovine serum albumin (BSA), and 0.5 μg poly(dI-dC) in a 20 μl reaction volume for 20 min at 26°C. The reaction mixture was applied to a 5.5% native polyacrylamide gel and subjected to electrophoresis in 1 × TAE, pH 8.5, at 4°C. Gel electrophoresis was performed in a GIBCO-BRL gel apparatus (model V16-2). The gel was dried and autoradiographed to examine the shift of the band.
Western Blot Analysis
To examine HapA production, Western blot analysis was carried out with 12 h grown spent cultures of recombinant strains harboring alanine variants along with the vector and wild-type controls. The cell-free supernatants were centrifuged and concentrated to 200 μl (10-fold) through 10 kDa Amicon ultracentrifugal filters (Millipore). The protein concentration was determined by the bicinchoninic acid (BCA) method and subjected to 12% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). Protein samples were electrophoretically transferred onto an Immobilon-P polyvinylidene fluoride (PVDF) membrane (Millipore). The Hap protein was detected by probing rabbit anti-Hap serum primary antibody at a dilution of 1:2,500 and horseradish peroxidase (HRP)-conjugated goat anti-rabbit immunoglobulin G (IgG) secondary antibody at a dilution of 1:5,000. Molecular masses were deduced by referring to the protein molecular weight markers.
Western blot analysis was done to check the cellular stability of the wild-type and mutant proteins. Briefly, V. cholerae strain S7 harboring FLAG-tagged recombinant derivatives of the wild-type and mutant constructs were grown in tryptic soy broth (TSB) without glucose overnight with agitation. The protein samples were prepared from an equal number of bacterial cells and separated on 12% gel by SDS-PAGE. For immunoblotting, the proteins were transferred to Immobilon-P PVDF membrane (Millipore) at 75 mA for 1 h. The membrane was subsequently blocked in phosphate-buffered saline (PBS) with 5% skim milk at 37°C for 2 h and then shifted to 4°C overnight with shaking. Blot was then washed in PBS Tween (PBST) five times for 10 min each, incubated in monoclonal HRP-conjugated anti-flag (Sigma Aldrich) at a dilution of 1:5,000 in PBST with 2% skim milk for 1 h, and again washed in PBST five times for 10 min each. The proteins were visualized using the Luminata Forte Western HRP substrate (Millipore). Molecular masses were calculated with reference to the SDS-PAGE molecular mass standards (broad range) from NAX-GEN ALPHA PS ladder from Genetix Biotech.
In vivo Oligomerization
The V. cholerae S7 strain bearing FLAG-tagged recombinant derivatives of the wild-type and mutant constructs were grown overnight. Equal numbers of cells were pelleted in duplicate for each culture and resuspended in protein sampling buffer ± DTT (100 mM) (Yang et al., 2013). The samples were then subjected to SDS-PAGE. Protein HapR was detected by Western blot using HRP-conjugated monoclonal anti-FLAG tag antibody (Sigma).
β-Galactosidase Reporter Assay
The LacZ reporter assay was performed as reported earlier (Miller, 1972). For the β-galactosidase assay, V. cholerae strain GK178 (aphA-lacZ fusion) (Kovacikova and Skorupski, 2002) was used by propagating the cultures overnight and diluting in AKI medium (Iwanaga et al., 1986) at 1:100 dilution. The assay was carried out after harvesting the cells at 0.3–0.4 OD600nm. The cell pellets were washed with Z-buffer (60 mM Na2HPO4⋅7H2O, 40 mM NaH2PO4×H2O, 10 mM KCl, 1 mM MgSO4⋅7H2O, and 50 mM β-mercaptoethanol) and resuspended in Sigma cell lytic buffer at a dilution of 1:10. The cells were sonicated at 20% amplitude for two cycles (10 s on and 10 s off). Of the clear lysate, 150 μl was mixed with 50 μl of ortho-nitrophenyl-β-galactoside [4 mg/ml in Z-buffer (pH 7) without β-mercaptoethanol] in 96-well plates and incubated at 30°C for 45 min. The β-galactosidase activity was observed by measuring the absorbance at 420 and 550 nm to calculate the Miller units (Miller, 1972).
Molecular Dynamics Simulation
Modeling of the B-DNA was performed using 3D-DART (van Dijk and Bonvin, 2009) web server. Energy minimization of the DNA was performed using the steepest descent algorithm with AMBER ff99bsc0 force field (Wang et al., 2004). The DNA was docked into HapR using Genetic Algorithm as the search parameter in the AutoDock-4.2 package. The output of the docking run was generated with the Lamarckian GA algorithm in AutoDock (Morris et al., 2009). The docked model was parameterized using AMBER ff99bsc force field and also applied to the system. Particle mesh Ewald algorithm was used for long-range electrostatics calculations and a 1.4 nm cutoff was used for short-range non-bonded interaction calculations. Periodic boundary condition was also applied to the system. The structure was well equilibrated using a 3 nm cubic box with TIP3P water. NaCl (100 mM) was added to the system with neutralizing counterions. Energy minimization was performed using the steepest descent minimization algorithm and a two-phase equilibration was followed thereafter. In the first phase, the structure was equilibrated with 50 ns of constant volume (NVT) ensemble. In second phase, another 50 ns of constant pressure (NPT) equilibration was performed. After the equilibration cycles, production molecular dynamics (MD) was carried out for 80 ns with no position restraint applied to the protein–DNA complex. During production MD, a Nosé–Hoover thermostat (Nosé, 1984; Hoover, 1985) was used for maintaining temperature and Parrinello–Rahmanbarostat (Parrinello and Rahman, 1981) was used for isotropic regulation of pressure. All the simulations were performed using GROMACS-0.4.6.2 package (Hess et al., 2008). Trajectory analysis and visualizations were performed using Chimera (Pettersen et al., 2004).
Results
Arginine 37 of Glycine Linker Controls the DNA Binding Activity of HapR
In order to identify the critical residue(s) in HapR–DNA interaction, we replaced each residue of the glycine linker by alanine from positions 33 to 39 (33RGIGRGG39; Figure 1). As HapR acts as an activator and repressor, we have chosen two well-characterized HapR targets, namely, HapA (hemagglutinin/protease) and AphA, to evaluate the function of alanine variants. HapR induces protease (HapA) production at a high cell density (Zhu et al., 2002). To evaluate functionality, the alanine variants of the linker region were tested for their capacity to reinstate protease production in V. cholerae strain S7 harboring a non-functional HapR (Supplementary Table 1). Protease production was measured with cell-free supernatants of the recombinant strains containing various alanine derivatives by employing a chromogenic substrate assay published earlier (Benitez et al., 2001; Singh et al., 2014). Out of all the variants, only three strains, those harboring HapR-G34A, HapR-I35A, and HapR-R37A, showed a significant loss in protease activity as compared to the vector and functional HapR controls (Figure 2A). Of these, the strain expressing HapR-R37A had the most dramatic defect in protease activity. The protease result was further corroborated by Western blot of HapA protease (Figure 2B). Our data also indicated the presence of both 45 and 35 kDa matured processed forms, as evidenced earlier (Pal and Wai, 2013). It could be surmised that loss in the protease activity of HapR-R37A could be a result of the instability of the protein under in vivo condition. To rule out the possibility, FLAG epitope was inserted in the C-terminal variant of the wild-type and alanine variants of HapR and stability was examined by Western blot (Supplementary Figure 1). To ensure that FLAG insertion did not alter the activity of the wild-type and alanine variants of HapR, protease production was examined with FLAG-tagged constructs. We observed no significant alteration in protease production as compared to the non-FLAG-tagged constructs (Figure 2A and Supplementary Figure 2). Our Western blot data clearly indicated that HapR-R37A is stable and its intactness is comparable to other alanine variants and wild-type protein. In other words, the loss in activity of HapR-R37A is not linked to its in vivo stability.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/544591/fmicb-11-01949-HTML/image_m/fmicb-11-01949-g001.jpg)
Figure 1. Pictorial representation of the glycine-rich linker (R33GIGRGG39) bridging the α-helix 1 and the HTH motif (α2 and α3) in the N-terminal binding domain and a schematic diagram showing the alanine scanning library.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/544591/fmicb-11-01949-HTML/image_m/fmicb-11-01949-g002.jpg)
Figure 2. HA/P protease activity and production in recombinant derivatives of Vibrio cholerae strain S7. (A) Azocasein assay. The strains were grown for 12 h in TSB-D at 37°C (200 rpm) and the proteolytic activity was determined by assaying the digestion of azocasein. Values are the mean ± SD (error bars) plotted from experiments performed in triplicate (*p < 0.05, ***p < 0.001; unpaired two-tailed Student’s t-test). S7-HapR was taken as the positive control. (B) Western blot analysis of Hap production in the cell-free supernatants. The supernatants of the 12 h grown cultures were concentrated to 200 μl (10-fold) through 10 kDa Amicon ultracentrifugal filters (Millipore). Hap protein was detected by probing rabbit anti-Hap serum primary antibody and horseradish peroxidase (HRP)-conjugated goat anti-rabbit immunoglobulin G (IgG) secondary antibody. Molecular masses were deduced by referring to the protein molecular weight markers. The 45- and 35-kDa sizes of the protein denote the differently processed HAP protein.
To examine the generality of our observation and rule out any possible strain-specific phenomenon, all the alanine variants along with the vector and wild-type functional HapR were transformed in another protease-negative strain of V. cholerae V2S (Dongre et al., 2011; Supplementary Table 1). We observed a similar trend of protease production in this strain background as well, where HapR-R37A exhibited a non-functional trait (Supplementary Figure 3). Taken together, the results obtained from protease with cell-free culture supernatants of V. cholerae strains S7 and V2S confirmed the compromise in the functionality of the HapR-R37A variant.
HapR negatively regulates vps gene expression and suppresses biofilm formation in V. cholerae (Beyhan et al., 2007; Lim et al., 2007). In some V. cholerae strains, mutation in HapR leads to the development of a rugose colony morphology (Yildiz et al., 2004). Interestingly, out of the two HapR-negative strains used in this study, V. cholerae V2S produces rugose colony on LB agar (Figure 3A). Subsequently, all recombinant strains of V. cholerae V2S carrying alanine variants along with the vector and wild-type functional HapR were examined for rugose colony development. We observed smooth colonies in all recombinant strains except the strain having HapR-R37A, which exhibited rugose morphotype (Figure 3A). As reported, rugose strains are less motile than their smooth counterparts (Ali et al., 2002). To examine this, motility assay was performed following a published protocol (Moisi et al., 2009). We observed that the smooth variants are more motile than the rugose variant (Figure 3B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/544591/fmicb-11-01949-HTML/image_m/fmicb-11-01949-g003.jpg)
Figure 3. Colony morphology, motility phenotypes, and β-galactosidase reporter assay to examine the functionality of HapR and its alanine variants. (A) Colony corrugation was examined on Luria–Bertani (LB) agar containing appropriate antibiotics after incubation at 30°C for 48 h. (B) Motility phenotypes were analyzed on swarm plates after incubation at 37°C for 12 h. (C) The Vibrio cholerae lacZ reporter strain GK178 (aphA-LacZ fusion) and the recombinants were grown in AKI medium at 37°C for 12 h. All strains were diluted in fresh AKI media and allowed to grow up to 0.3–0.4 OD600 nm. The β-gal production was measured as a readout of o-nitrophenol, a cleaved product of substrate ONPG, and further calculated the Miller unit. The amount of Miller units denoted the levels of lacZ expression and correlates the functionality of wild-type HapR and its alanine variants. HapR was taken as the positive control. Values shown are the mean and standard error of at least two biological experiments and three technical replicates (*p < 0.05; unpaired two-tailed Student’s t-test).
HapR represses the expression of the gene encoding AphA by binding to its promoter region and occludes the engagement of VpsR, the latter being an activator of AphA (Kovacikova and Skorupski, 2002; Lin et al., 2007). To examine the interaction of HapR-R37A with the promoter region of AphA, a lacZ reporter assay was conducted in a recombinant V. cholerae strain C6706, designated as GK178, carrying chromosomal fusion of aphA–lacZ (Supplementary Table 1; Kovacikova and Skorupski, 2002). In the reporter strain, binding of functional HapR to the aphA promoter region will block lacZ expression; therefore, the presence or absence of lacZ expression is directly proportional to the function of HapR. Our reporter assay with V. cholerae strain GK178 harboring all linker alanine variants and wild-type functional HapR clearly indicated that HapR-R37A failed to inhibit lacZ expression, further indicating its inability to bind the aphA promoter region (Figure 3C).
Since arginine at position 37 (R37) turned out to be the most important residue among others in the glycine linker, we subsequently kept our experimental focus on R37 to understand its role in the HapR regulatory activity. HapR directly interacts with the promoter region of aphA and hapA, thus modulating the expression of these proteins (Zhu et al., 2002). To evaluate the binding ability of HapR-R37A with these promoter regions, gel shift assays were performed as described earlier (Dongre et al., 2011). Gel shift analysis clearly demonstrated complete failure of the binding of HapR-R37A to the promoter region of aphA and hapA over a threefold range of concentrations (Figure 4). The data are in congruence with in vivo protease, rugosity, motility, and lacZ reporter assays (Figures 2A, 3), which further bolstered the importance of arginine 37 in the DNA binding regulatory activity of HapR.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/544591/fmicb-11-01949-HTML/image_m/fmicb-11-01949-g004.jpg)
Figure 4. Gel shift experiments showing the inability of mutant HapR to bind to its cognate promoters. Electrophoretic mobility shift assay (EMSA) of the purified wild-type protein HapR and its mutant HapRR37A was performed with 32P-labeled promoter regions of aphA (VC2647), hapA (VCA0865), and cdgG (VC0900). The solid wedges represent the increasing concentrations of proteins. The results are representative of two or more experiments.
As promoter regions of aphA and hapA harbor HapR binding motif 2, we chose another cognate promoter, vc0900, having HapR binding motif 1 (Tsou et al., 2009) and examined the DNA binding ability of HapR-R37A. Our data clearly demonstrated the inability of HapR-R37A to shift the promoter regions of vc0900 (Figure 4).
In essence, combined in vivo and in vitro results aid in the identification of an essential arginine at position 37 (R37) for HapR function. Replacement of arginine with alanine (R37A) renders a completely non-functional HapR that fails to bind with the cognate promoters examined in this study.
Necessity of Positive Charge at Position 37 to Maintain HapR Regulatory Function
It is evident from the previous section that arginine 37 is indispensable for HapR function. The question arises whether residue or charge is important at position 37. It should be noted that replacement of arginine with lysine is also associated with functional modulation, further indicating the importance of residue over charge (Aukerman et al., 1991; Tang et al., 2009). To address this, arginine 37 was substituted with lysine, glutamate, and aspartate. Incidentally, a natural variant of arginine 37 has been reported where the residue is substituted with histidine (Wang et al., 2011). All variants (R37A, R37K, R37E, R37D, and R37H) were subjected to a battery of in vivo and in vitro experiments (e.g., protease, hemagglutination, rugosity, motility, and gel shift). As documented, the HapR-R37K variant is functionally equivalent, while the functionality of the HapR-R37A, HapR-R37D, and HapR-R37E variants is significantly compromised in all experimental conditions in comparison to their wild-type counterpart (HapR-R37) (Figures 5, 6), whereas HapR-R37H appears to retain partial activity. To underscore that the loss in functionality of the HapR-R37D, HapR-R37E, and HapR-R37H variants is not linked with protein stability, FLAG Western blot was performed. No significant changes were observed in protease production after insertion of the FLAG epitope at the C-terminus of these variants (Supplementary Figure 4). Western blot data also documented the in vivo stability of all proteins (Supplementary Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/544591/fmicb-11-01949-HTML/image_m/fmicb-11-01949-g005.jpg)
Figure 5. In vivo assays of the recombinant derivatives of Vibrio cholerae strains S7 and V2S bearing wild-type HapR and its R37 variants (R37A, R37K, R37D, R37H, and R37E). (A) Protease activity. All V. cholerae S7 recombinant strains were grown for 12 h in TSB-D at 37°C (200 rpm) and the proteolytic activity was determined in cell-free culture supernatants by assaying the digestion of azocasein. Values are the mean ± standard deviation (error bars) plotted from experiments performed in triplicate (**p < 0.01; unpaired two-tailed Student’s t-test; ns, non-significant). S7-HapR was taken as the positive control. (B) Hemagglutination activity. The assay was performed by making a serial twofold dilution of the culture supernatants of the overnight (12 h)-grown V. cholerae S7 recombinant variants. The samples were incubated with 1.2% washed chicken erythrocytes and the ability of HA/P to hemagglutinate the red blood cells (RBCs) was checked. The results are representative of three independent experiments. (C) Colony morphology. Colony morphotype of V. cholerae V2S recombinant strains was examined on Luria–Bertani (LB) agar containing appropriate antibiotics after incubation at 30°C for 48 h. (D) Motility phenotypes. Phenotypic analysis of the motility of V. cholerae V2S recombinant variants was done on swarm plates containing appropriate antibiotics with 12 h post-incubation at 37°C.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/544591/fmicb-11-01949-HTML/image_m/fmicb-11-01949-g006.jpg)
Figure 6. Gel shift experiments of HapR and its R37 mutants. Electrophoretic mobility shift assay (EMSA) of the purified wild-type protein HapR and its mutants HapRR37A, HapRR37K, HapRR37D, HapRR37H, and HapRR37E was performed with 32P-labeled promoter regions of (A) aphA (VC2647), (B) hapA (VCA0865), and (C) cdgG (VC0900). The solid wedges represent the increasing concentrations of proteins. The results are representative of two or more experiments.
DNA Binding Impairment Does Not Affect Oligomerization
As documented, HapR exists as a dimer (De Silva et al., 2007; Dongre et al., 2011). It could be surmised that the non-functionality of the HapR-R37A, HapR-R37H, HapR37E, and HapR-R37D variants could be a result of an altered oligomeric pattern under in vivo condition. To further investigate, the oligomeric status under in vivo condition was examined by Western blot following protocols published earlier (Yang et al., 2013). Western blot data clearly exhibited the dimeric status of HapR-R37A, HapR-R37H, HapR-R37E, and HapR-R37D under the experimental condition (Figure 7).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/544591/fmicb-11-01949-HTML/image_m/fmicb-11-01949-g007.jpg)
Figure 7. In vivo oligomerization. S7 strain of Vibrio cholerae possessing the C-terminal FLAG-tagged hapR and its R37 mutant derivatives were propagated in TSB-D until the late log phase. Equal numbers of cells were harvested and the samples were prepared from cell lysates in duplicate ± DTT (100 mM in the loading buffer). Samples were subjected to gel electrophoresis and the oligomeric state of the proteins was detected by immunoblotting with anti-FLAG antibody. The sizes ∼27 and 54 kDa depict the monomeric and dimeric states of the FLAG-tagged HapR protein, thereby suggesting that the mutations do not affect the oligomerization of the protein.
MD Simulation Unravels the Role of Arg37 (R37) in HapR-Mediated DNA Binding
De Silva et al. (2007) demonstrated that Phe55 is very crucial for DNA binding by HapR. However, in this study, using various experimental approaches, we have been able to point out that, besides F55, another residue (R37) plays a very crucial role in DNA binding. To generate a more vivid picture of the scenario, we further performed structural analysis of binding using extensive MD simulations. Earlier, the promoter region of hapA was considered in MD simulation experiments (Singh et al., 2014). In this study too, we used the same DNA sequence and modeled a B-DNA using 3D-DART web server (van Dijk and Bonvin, 2009). Energy minimization of the DNA was performed using AMBER ff99bsc0 force field (Wang et al., 2004). Further, to generate a HapR–DNA complex, a docking procedure was carried out using AutoDock-4.2 package (Morris et al., 2009). Docking results were screened to obtain a complex structure with a total energy of −9.7 kcal/mol. The docked structure was then subjected to 50 ns of NVT and NPT equilibration using AMBER force field (Wang et al., 2004) with position restrains applied to the complex. The well-equilibrated structure was further introduced to 80 ns of non-restrained MD simulation. Analyzing the final structure, we found out that F55 from both the chains of HapR was taking part in stabilizing the bound DNA by forming electrostatic interactions (De Silva et al., 2007; Singh et al., 2014). We also found out that the R37 residues from both the chains of HapR are forming a hydrogen bond with a “handshake” topology to hold the two domains of the protein in a desired distance for binding with the DNA (Figure 8). Analysis of the trajectories of the simulation also showed that the handshake is not due to any charged ions, added to neutralize the counterion during simulations, lying nearby R37 (Supplementary Figure 6). R37 is also found to form electrostatic interaction with the main chain oxygen atom of I35, although this interaction is not true for both the R37 residues. However, to substantiate our claim, we further needed to mutate the R37 and performed extensive MD simulation. For this purpose, we took the final structure produced by an earlier simulation and substituted both the R37 with alanine using PyMOL package (DeLano, 2009). We used alanine substitution to completely eliminate any possibility of ionic interactions between two of these residues and, in turn, would provide us a way to check whether R37’s “handshake” topology was affecting the DNA binding or we just hit on an artifactual phenomenon.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/544591/fmicb-11-01949-HTML/image_m/fmicb-11-01949-g008.jpg)
Figure 8. HapR with bound DNA. (A) HapR–DNA complex after 80 ns of non-restrained molecular dynamics (MD) simulation. Arginine 37 (R37) from both the chains are represented in a ball-and-stick model. (B) Close-up view shows that the orientation of R37 residues forms an ion pair possessing a “handshake” topology.
The mutant (HapR-R37A) was then further equilibrated for 80 ns of NVT and NPT simulations until the backbone root mean square division (RMSD) reached a certain amount of stability (Supplementary Figure 7). The equilibrated structure was then further subjected to 80 ns of non-restrained MD simulation. Analysis of the trajectory showed that, after a few cycles of non-restrained simulations, the two domains of HapR-R37A failed to hold their position properly for providing interaction with the DNA. This vigorous domain movement eventually leads to the disruption of the electrostatic interactions between HapR-R37A and the DNA (Figure 9, Supplementary Figure 8, and Supplementary Video 1). Further, to check whether other positively charged amino acids have a similar effect, we chose to mutate the R37 with histidine (R37H). Simulations were carried out using a similar set of parameterizations to that used before. Non-restrained simulation results show that R37H is also unable to hold the DNA as the length of the histidine side chain is not long enough to establish the handshake (Supplementary Figure 9). In another similar simulation setup, we evaluated how the complex would behave if R37 is substituted with glutamate (E), which has the negatively charged side chains. The R37E simulations also failed to hold the DNA when subjected to 80 ns of non-restrained simulation. Analysis of the simulation snapshot shows that the E is 7 Å apart from each other and, due to this, is not able to form any kind of stabilizing interactions (Supplementary Figure 10). These MD simulation studies clarify the fact that the “handshake” mechanism of R37 is highly significant to fine-tune the domain movement in a corroborated manner, which allows the DNA to bind with HapR.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/544591/fmicb-11-01949-HTML/image_m/fmicb-11-01949-g009.jpg)
Figure 9. Superimposition of HapR-DNA bound states (HapR–DNA is represented in light brown color; HapRR37A–DNA is represented in transparent cyan color). (A) Structural alignment shows that a large helix movement occurred due to the alanine substitution and which eventually led to the disruption of the protein–DNA complex. (B) 90° rotated view of the aligned structures.
Based on our investigations performed so far, it could be hypothesized that not the residue in the position but the residue carrying specific charge matters establishes a positive HapR–DNA interaction. Hence, we further substituted the Arg37 residue with lysine, a residue closely resembling the same electrostatic properties. The mutant (HapR-R37K) was then equilibrated using 80 ns of NVT and NPT simulations to obtain a stable backbone RMSD (Supplementary Figure 11). After successful equilibration, another 80 ns of non-restrained MD simulation was carried out. Analysis of the trajectories of the non-restrained run of HapR-R37K showed that the rapid domain movements were very few compared to the HapRR37A mutant. Snapshots collected from the MD run showed that the Lys37 residue formed a hydrogen bond with each other and kept the vigorous domain movement at bay to facilitate the stability of the HapR (Figure 10), thus sustaining the HapR–DNA interaction. This also explains the fact that after 80 ns of non-restrained simulations, the DNA was found to be in a bound state with HapR-R37K. This phenomenon is also in good agreement with the in vitro and in vivo experimental study described in the preceding sections. Our current investigation delineates that the R37 residue is quite indispensable in HapR–DNA interactions, however, it can be replaced with specific residues with nearly similar electrostatic and structural signatures. In essence, our experimental and computational results have led to the identification of an arginine residue at position 37 that plays a critical role in the DNA binding activity of HapR. However, a similar investigation with the crystal structures of these complexes should provide more information if performed in the future.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/544591/fmicb-11-01949-HTML/image_m/fmicb-11-01949-g010.jpg)
Figure 10. HapRR37K with bound DNA. (A) HapRR37K–DNA complex after 80 ns of non-restrained molecular dynamics (MD) simulation. Lysine 37 (K37) from both the chains are represented as sticks. (B) Close-up snapshot of the K37 residues shows that these are forming ion pairs to hold the domains, controlling the rapid uncorrelated domain movement.
Discussion
HapR contributes immensely in host and non-host life cycles of Vibrio cholerae by modulating the function of the genes responsible for disparate cellular events in a cell density-dependent manner. To achieve this feat, HapR communicates effectively with a series of cognate promoters by recognizing the binding motifs in the promoter region and acting both as an activator and a repressor. The recognition is done by a putative DNA binding domain (DBD) situated in the N-terminal region spanning over three α-helices and comprising initial 57 amino acids, where helix 3 is the recognition helix. Previously, phenylalanine at position 55 (F55) in the recognition helix was the only amino acid linked to the DNA binding ability of HapR (De Silva et al., 2007). Later, the contribution of a glycine-rich linker connecting helices 1 and 2 in the DNA binding activity was brought to light (Dongre et al., 2011; Singh et al., 2014). The current study further delves into the composition of the linker region and eventually identifies an important arginine residue at position 37 contributing toward the DNA binding ability of HapR. The next question arises on the importance of charge vs. residue at this position. Interestingly, accumulated evidences indicating that a conserved arginine replaced with lysine impairs the DNA binding ability (Aukerman et al., 1991; Tang et al., 2009) left us to examine the effect of such substitution on R37 of HapR. Subsequently, R37 was replaced with negatively charged moieties such as aspartate and glutamate. We observed no loss in the DNA binding activity of the lysine-substituted variant, but the aspartate and glutamate variants turned functionally inert under experimental conditions. The data further reinforced the importance of positive charge rather than residue at position 37 for HapR–DNA cross talk.
It should be noted that natural variants of the glycine linker have also been reported where arginine 37 and glycine 39 were found to be substituted with histidine (R37H) and aspartate (G39D), respectively (Dongre et al., 2011; Wang et al., 2011). Previously, the molecular characterization of HapR-G39D has been done (Dongre et al., 2011; Cruite et al., 2018). In this study, we confirmed the non-functionality of R37H, described in the preceding sections.
The linker contributes immensely in the function of a protein. Though the majority of linkers are rigid, glycine-rich linkers are flexible in nature and maintain the functional integrity of discrete domains (Reddy Chichili et al., 2013). A large body of evidence bolsters the various functional roles of glycine-rich linkers in diverse proteins. While the glycine linker in PAX6, a human transcriptional factor, promotes contacts with DNA (Xu et al., 1999), the glycine linker in the transmembrane glycoprotein of retroviruses participates in membrane fusion (Wilson et al., 2005). Likewise, the glycine-rich linker of human T-cell leukemia virus type 1 transmembrane is also involved in membrane fusion (Wilson et al., 2005). In the case of SRSF1 (sequence-specific RNA binding factor), the glycine-rich linker connects two RNA recognition motifs and promotes interaction with splicing enhancers (Cho et al., 2011). Notably, mutations in the Gly-rich domain of TDP-43 have led to the development of amyotrophic lateral sclerosis (Pesiridis et al., 2009).
Quorum sensing defective strains of V. cholerae harboring non-functional variants of HapR are abundant in nature, and loss of HapR function also attributes certain physiological advantages to such strains (Joelsson et al., 2007; Dongre et al., 2011; Wang et al., 2011; Singh et al., 2013). Intriguingly, functionally challenged natural variants also aid in gaining novel mechanistic insights on HapR (Dongre et al., 2011; Singh et al., 2013). Incidentally, the present work stems out from our previous findings on the molecular basis of a non-functional natural quorum sensing variant where the significance of the glycine-rich linker, for the first time in the DNA binding activity of HapR, was revealed (Dongre et al., 2011). The current work further adds on the existing knowledge on HapR function by identifying a key arginine residue in linker sequence. It is noteworthy to mention that arginine 36 (equivalent to arginine 37 of HapR) of SmcR also contributes in DNA binding (Kim et al., 2010), further emphasizing the importance of arginine 37 in addition to phenylalanine 55 (De Silva et al., 2007) in the DNA binding activity of HapR. As both arginine 37 and phenylalanine 55 are well conserved among HapR homologs in several species of Vibrio (Supplementary Figure 12), it further necessitates evaluating the contribution of both residues in the functionality of HapR homologs in different Vibrio species. This warrants further investigation.
Data Availability Statement
All datasets presented in this study are included in the article/Supplementary Material.
Author Contributions
SR conceived the idea and wrote the manuscript. SR and ME designed the experiments and analyzed all biochemical data. ME performed the experiments. AM and SD performed and analyzed the MD simulation data. HS and RS contributed to the reagents, materials, and analysis tools. All authors gave editorial inputs.
Funding
This work was partly supported by grants from the Department of Biotechnology (BT/PR3583/MED/30/659/2011-GAP/106) and the Science and Engineering Research Board (CRG/2018/000297/SERB-GAP/0185), respectively. ME and HS acknowledge University Grant Commission (UGC), SERB (DST), and Council of Scientific and Industrial Research (CSIR), respectively, for fellowships.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We gratefully acknowledge Dr. Giesla Storz, National Institutes of Health, Bethesda, MD, United States, for providing the pKK177-3RI plasmid; Drs. Gabriela Kovacikova and Karen Skorupski, Department of Microbiology and Immunology, Dartmouth Medical School, Hanover, United States, for the V. cholerae GK178 strain; and Dr. Ranjan Nandy, NICED, Kolkata, India, for the V. cholerae S7 strain. We extend our appreciation to Dr. Johnny Peterson for the anti-Hap antibody. We greatly acknowledge Drs. Mitesh Dongre and Naorem Santa Singh for their initial participation in this project. SR and ME extend gratitude to Nandita Sharma and Manpreet Kaur, lab members of SR, for critical evaluation of the entire manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.01949/full#supplementary-material
References
Ali, A., Rashid, M. H., and Karaolis, D. K. R. (2002). High-frequency rugose exopolysaccharide production by Vibrio cholerae. Appl. Environ. Microbiol. 68, 5773–5778. doi: 10.1128/AEM.68.11.5773-5778.2002
Aukerman, M. J., Schmidt, R. J., Burr, B., and Burr, F. A. (1991). An arginine to lysine substitution in the bZIP domain of an opaque-Z mutant in maize abolishes specific DNA binding. Genes Dev. 5, 310–320. doi: 10.1101/gad.5.2.310
Benitez, J. A., Silva, A. J., and Finkelstein, R. A. (2001). Environmental signals controlling production of hemagglutinin/protease in Vibrio cholerae. Infect. Immun. 69, 6549–6553. doi: 10.1128/IAI.69.10.6549-6553.2001
Beyhan, S., Bilecen, K., Salama, S. R., Casper-Lindley, C., and Yildiz, F. H. (2007). Regulation of rugosity and biofilm formation in Vibrio cholerae: comparison of VpsT and VpsR regulons and epistasis analysis of vpsT, vpsR, and hapR. J. Bacteriol. 189, 388–402. doi: 10.1128/JB.00981-06
Camacho, A., Bouhenia, M., Alyusfi, R., Alkohlani, A., Naji, M. A. M., de Radiguès, X., et al. (2018). Cholera epidemic in Yemen, 2016–18: an analysis of surveillance data. Lancet Glob. Heal. 6, e680–e690. doi: 10.1016/S2214-109X(18)30230-4
Cho, S., Hoang, A., Chakrabarti, S., Huynh, N., Huang, D.-B., and Ghosh, G. (2011). The SRSF1 linker induces semi-conservative ESE binding by cooperating with the RRMs. Nucleic Acids Res. 39, 9413–9421. doi: 10.1093/nar/gkr663
Cruite, J., Succo, P., Raychaudhuri, S., and Kull, F. J. (2018). Crystal structure of an inactive variant of the quorum-sensing master regulator HapR from the protease-deficient non-O1, non-O139 Vibrio cholerae strain V2. Acta Crystallogr. Sect. F Struct. Biol. Commun. 74, 331–336. doi: 10.1107/S2053230X18006519
De Silva, R. S., Kovacikova, G., Lin, W., Taylor, R. K., Skorupski, K., and Kull, F. J. (2007). Crystal structure of the Vibrio cholerae quorum-sensing regulatory protein HapR. J. Bacteriol. 189, 5683–5691. doi: 10.1128/JB.01807-06
Dongre, M., Singh, N. S., Dureja, C., Peddada, N., Solanki, A. K., and Ashish, et al. (2011). Evidence on how a conserved glycine in the hinge region of HapR regulates its DNA binding ability: lessons from a natural variant. J. Biol. Chem. 286, 15043–15049. doi: 10.1074/jbc.M110.209346
Giglio, K. M., Fong, J. C., Yildiz, F. H., and Sondermann, H. (2013). Structural basis for biofilm formation via the Vibrio cholerae matrix protein RbmA. J. Bacteriol. 195, 3277–3286. doi: 10.1128/JB.00374-13
Hess, B., Kutzner, C., Van Der Spoel, D., and Lindahl, E. (2008). GRGMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 4, 435–447. doi: 10.1021/ct700301q
Hoover, W. G. (1985). Canonical dynamics: equilibrium phase-space distributions. Phys. Rev. A 31, 1695–1697. doi: 10.1103/PhysRevA.31.1695
Iwanaga, M., Yamamoto, K., Higa, N., Ichinose, Y., Nakasone, N., and Tanabe, M. (1986). Culture Conditions for Stimulating Cholera Toxin Production by Vibrio cholerae 01 El Tor. Microbiol. Immunol. 30, 1075–1083. doi: 10.1111/j.1348-0421.1986.tb03037.x
Joelsson, A., Kan, B., and Zhu, J. (2007). Quorum sensing enhances the stress response in Vibrio cholerae. Appl. Environ. Microbiol. 73, 3742–3746. doi: 10.1128/AEM.02804-06
Kim, Y., Kim, B. S., Park, Y. J., Choi, W.-C., Hwang, J., Kang, B. S., et al. (2010). Crystal structure of SmcR, a quorum-sensing master regulator of Vibrio vulnificus, provides insight into its regulation of transcription. J. Biol. Chem. 285, 14020–14030. doi: 10.1074/jbc.M109.100248
Kovacikova, G., and Skorupski, K. (2002). Regulation of virulence gene expression in Vibrio cholerae by quorum sensing: HapR functions at the aphA promoter. Mol. Microbiol. 46, 1135–1147. doi: 10.1046/j.1365-2958.2002.03229.x
Lenz, D. H., Mok, K. C., Lilley, B. N., Kulkarni, R. V., Wingreen, N. S., and Bassler, B. L. (2004). The small RNA chaperone Hfq and multiple small RNAs control quorum sensing in Vibrio harveyi and Vibrio cholerae. Cell 118, 69–82. doi: 10.1016/j.cell.2004.06.009
Lim, B., Beyhan, S., and Yildiz, F. H. (2007). Regulation of Vibrio polysaccharide synthesis and virulence factor production by CdgC, a GGDEF-EAL domain protein, in Vibrio cholerae. J. Bacteriol. 189, 717–729. doi: 10.1128/JB.00834-06
Lin, W., Kovacikova, G., and Skorupski, K. (2007). The quorum sensing regulator HapR downregulates the expression of the virulence gene transcription factor AphA in Vibrio cholerae by antagonizing Lrp- and VpsR-mediated activation. Mol. Microbiol. 64, 953–967. doi: 10.1111/j.1365-2958.2007.05693.x
Miller, J. H. (ed.) (1972). “Assay of b-galactosidase,” in Experiments in Molecular Genetics (Cold Spring, NY: Cold Spring Harbor Laboratory), 352–355.
Moisi, M., Jenul, C., Butler, S. M., New, A., Tutz, S., Reidl, J., et al. (2009). A novel regulatory protein involved in motility of Vibrio cholerae. J. Bacteriol. 191, 7027–7038. doi: 10.1128/JB.00948-09
Morris, G. M., Ruth, H., Lindstrom, W., Sanner, M. F., Belew, R. K., Goodsell, D. S., et al. (2009). Software news and updates AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility. J. Comput. Chem. 30, 2785–2791. doi: 10.1002/jcc.21256
Ng, W.-L., and Bassler, B. L. (2009). Bacterial quorum-sensing network architectures. Annu. Rev. Genet. 43, 197–222. doi: 10.1146/annurev-genet-102108-134304
Nosé, S. (1984). A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 81, 511–519. doi: 10.1063/1.447334
Pal, A., and Wai, S. N. (2013). “Proteases in Vibrio cholerae and their role in pathogenesis,” in Proteases in Health and Disease, Advances in Biochemistry in Health and Disease, Vol. 7, eds S. Chakraborti and N. Dhalla (New York: Springer), 243–256. doi: 10.1007/978-1-4614-9233-7_14
Papenfort, K., Silpe, J. E., Schramma, K. R., Cong, J. P., Seyedsayamdost, M. R., and Bassler, B. L. (2017). A Vibrio cholerae autoinducer-receptor pair that controls biofilm formation. Nat. Chem. Biol. 13, 551–557. doi: 10.1038/nchembio.2336
Parrinello, M., and Rahman, A. (1981). Polymorphic transitions in single crystals: a new molecular dynamics method. J. Appl. Phys. 52, 7182–7190. doi: 10.1063/1.328693
Pesiridis, G. S., Lee, V. M. Y., and Trojanowski, J. Q. (2009). Mutations in TDP-43 link glycine-rich domain functions to amyotrophic lateral sclerosis. Hum. Mol. Genet. 18, R156–R162. doi: 10.1093/hmg/ddp303
Pettersen, E. F., Goddard, T. D., Huang, C. C., Couch, G. S., Greenblatt, D. M., Meng, E. C., et al. (2004). UCSF Chimera - A visualization system for exploratory research and analysis. J. Comput. Chem. 25, 1605–1612. doi: 10.1002/jcc.20084
Raychaudhuri, S., Jain, V., and Dongre, M. (2006). Identification of a constitutively active variant of LuxO that affects production of HA/protease and biofilm development in a non-O1, non-O139 Vibrio cholerae O110. Gene 369, 126–133. doi: 10.1016/j.gene.2005.10.031
Reddy Chichili, V. P., Kumar, V., and Sivaraman, J. (2013). Linkers in the structural biology of protein-protein interactions. Protein Sci. 22, 153–167. doi: 10.1002/pro.2206
Seedorff, J., and Schleif, R. (2011). Active role of the interdomain linker of AraC. J. Bacteriol. 193, 5737–5746. doi: 10.1128/JB.05339-11
Singh, N. S., Kachhap, S., Singh, R., Mishra, R. C., Singh, B., and Raychaudhuri, S. (2014). The length of glycine-rich linker in DNA-binding domain is critical for optimal functioning of quorum-sensing master regulatory protein HapR. Mol. Genet. Genomics 289, 1171–1182. doi: 10.1007/s00438-014-0878-5
Singh, R., Rathore, Y. S., Singh, N. S., Peddada, N., Ashish, and Raychaudhuri, S. (2013). Substitution of glutamate residue by lysine in the dimerization domain affects dna binding ability of HapR by inducing structural deformity in the DNA binding domain. PLoS One 8:e76033. doi: 10.1371/journal.pone.0076033
Tang, X., Narayanan, S., Peruzzi, G., Apara, A., Natarajan, K., Margulies, D. H., et al. (2009). A single residue, arginine 65, is critical for the functional interaction of leukocyte-associated inhibitory receptor-1 with collagens. J. Immunol. 182, 5446–5452. doi: 10.4049/jimmunol.0804052
Tsou, A. M., Cai, T., Liu, Z., Zhu, J., and Kulkarni, R. V. (2009). Regulatory targets of quorum sensing in Vibrio cholerae: evidence for two distinct HapR-binding motifs. Nucleic Acids Res. 37, 2747–2756. doi: 10.1093/nar/gkp121
van Dijk, M., and Bonvin, A. M. J. J. (2009). 3D-DART: a DNA structure modelling server. Nucleic Acids Res. 37, W235–W239. doi: 10.1093/nar/gkp287
Van Leeuwen, H. C., Strating, M. J., Rensen, M., De Laat, W., and Van Der Vliet, P. C. (1997). Linker length and composition influence the flexibility of Oct-1 DNA binding. EMBO J. 16, 2043–2053. doi: 10.1093/emboj/16.8.2043
Velicer, G. J. (2003). Social strife in the microbial world. Trends Microbiol. 11, 330–337. doi: 10.1016/S0966-842X(03)00152-5
Wang, J., Wolf, R. M., Caldwell, J. W., Kollman, P. A., and Case, D. A. (2004). Development and testing of a general Amber force field. J. Comput. Chem. 25, 1157–1174. doi: 10.1002/jcc.20035
Wang, Y., Wang, H., Cui, Z., Chen, H., Zhong, Z., Kan, B., et al. (2011). The prevalence of functional quorum-sensing systems in recently emerged Vibrio cholerae toxigenic strains. Environ. Microbiol. Rep. 3, 218–222. doi: 10.1111/j.1758-2229.2010.00212.x
Wilson, K. A., Bar, S., Maerz, A. L., Alizon, M., and Poumbourios, P. (2005). The conserved glycine-rich segment linking the N-terminal fusion peptide to the coiled coil of human T-cell leukemia virus type 1 transmembrane glycoprotein gp21 is a determinant of membrane fusion function. J. Virol. 79, 4533–4539. doi: 10.1128/jvi.79.7.4533-4539.2005
Xu, E. H., Rould, M. A., Xu, W., Epstein, J. A., Maas, R. L., and Pabo, C. O. (1999). Crystal structure of the human Pax6 paired domain-DNA complex reveals specific roles for the linker region and carboxy-terminal subdomain in DNA binding. Genes Dev. 13, 1263–1275. doi: 10.1101/gad.13.10.1263
Yang, M., Liu, Z., Hughes, C., Stern, A. M., Wang, H., Zhong, Z., et al. (2013). Bile salt-induced intermolecular disulfide bond formation activates Vibrio cholerae virulence. Proc. Natl. Acad. Sci. U.S.A. 110, 2348–2353. doi: 10.1073/pnas.1218039110
Yildiz, F. H., Liu, X. S., Heydorn, A., and Schoolnik, G. K. (2004). Molecular analysis of rugosity in a Vibrio cholerae O1 El Tor phase variant. Mol. Microbiol. 53, 497–515. doi: 10.1111/j.1365-2958.2004.04154.x
Keywords: HapR, Vibrio cholerae, quorum sensing, linker, arginine, DNA binding, molecular dynamics
Citation: Ekka M, Mondal A, Singh R, Sen H, Datta S and Raychaudhuri S (2020) Arginine 37 of Glycine Linker Dictates Regulatory Function of HapR. Front. Microbiol. 11:1949. doi: 10.3389/fmicb.2020.01949
Received: 21 March 2020; Accepted: 24 July 2020;
Published: 21 August 2020.
Edited by:
Dongsheng Zhou, Beijing Institute of Microbiology and Epidemiology, ChinaReviewed by:
Eric Krukonis, University of Detroit Mercy, United StatesSamit Watve, Tufts Medical Center, United States
Copyright © 2020 Ekka, Mondal, Singh, Sen, Datta and Raychaudhuri. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Saumya Raychaudhuri, saumya@rocketmail.com; saumya@imtech.res.in