|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 281, Issue 3, 1532-1546, January 20, 2006
Formation of a Distinctive Complex between the Inducible Bacterial Lysine Decarboxylase and a Novel AAA+ ATPase* 1 2![]() ![]() 4
From the
Received for publication, October 13, 2005
AAA+ ATPases are ubiquitous proteins that employ the energy obtained from ATP hydrolysis to remodel proteins, DNA, or RNA. The MoxR family of AAA+ proteins is widespread throughout bacteria and archaea but is largely uncharacterized. Limited work with specific members has suggested a potential role as molecular chaperones involved in the assembly of protein complexes. As part of an effort aimed at determining the function of novel AAA+ chaperones in Escherichia coli, we report the characterization of a representative member of the MoxR family, YieN, which we have renamed RavA (regulatory ATPase variant A). We show that the ravA gene exists on an operon with another gene encoding a protein, YieM, of unknown function containing a Von Willebrand Factor Type A domain. RavA expression is under the control of the S transcription factor, and its levels increase toward late log/early stationary phase, consistent with its possible role as a general stress-response protein. RavA functions as an ATPase and forms hexameric oligomers. Importantly, we demonstrate that RavA interacts strongly with inducible lysine decarboxylase (LdcI or CadA) forming a large cage-like structure consisting of two LdcI decamers linked by a maximum of five RavA oligomers. Surprisingly, the activity of LdcI does not appear to be affected by binding to RavA in a number of in vitro and in vivo assays, however, complex formation results in the stimulation of RavA ATPase activity. Data obtained suggest that the RavA-LdcI interaction may be important for the regulation of RavA activity against its targets.
The AAA+ (ATPases associated with various cellular activities) superfamily of proteins comprises a large and functionally diverse class of P-loop NTPases. Members of this superfamily are involved in a wide range of cellular processes, including protein refolding/degradation, transcriptional regulation, ribosome and organelle biogenesis, DNA repair and replication, and molecular transport events (1, 2). Proteins belonging to this superfamily contain one or more copies of the AAA+ module, which is a 200- to 250-amino acid region containing several distinct conserved motifs (1, 3, 4), including the Walker A and Walker B consensus sequences. AAA+ modules are responsible for ATP binding and hydrolysis, the energy of which is then harnessed for use in molecular rearrangement events (5). AAA+ proteins oligomerize, typically forming ring-shaped hexamers (6). Recent phylogenetic analyses using sequence and structural information have identified numerous families of AAA+ proteins (1, 2, 7, 8). One major but poorly characterized family is the MoxR family, which is widespread throughout bacteria and archaea, with members being represented in all major lineages (1, 2, 8). The exact functional role of MoxR proteins is unclear, although limited work with members of this family suggests a chaperone-like role in the assembly/activation of specific protein complexes. The first characterized member of this family is MoxR in Paracoccus denitrificans, which was shown to be important in the biogenesis of methanol dehydrogenase. Cells in which the moxR gene was disrupted produced wild-type levels of methanol dehydrogenase, but, surprisingly, the enzyme was not functional (9). Another MoxR-related protein, NorQ/NirQ, was shown to be important in the biogenesis of nitric-oxide reductase. As with MoxR and methanol dehydrogenase, disruptions of the norQ/nirQ gene resulted in the production of non-functional nitric-oxide reductase enzyme in Pseudomonas stutzeri, P. denitrificans, and Rhodobacter sphaeroides (1012). Also, the CbbQ gene product, a MoxR-related protein, was shown to alter Rubisco4 conformation and stimulate its activity upon overexpression of recombinant forms of the proteins together in an Escherichia coli system (13, 14). Yet another MoxR family member, GvpN, has been shown to be important in gas vesicle biogenesis in certain organisms (1518). The exact role and mechanism of function of MoxR proteins in all of these processes are not clear.
Genes encoding MoxR proteins frequently co-occur in close proximity to genes encoding proteins containing Von Willebrand Factor Type A (VWA) domains, suggesting that these gene products may be functionally linked. For example, the MoxR family members, norQ and cbbQ genes, are encoded on the same operon as, and immediately upstream of, the norD and cbbO genes, respectively, which encode VWA proteins. Both NorD and CbbO were shown to be important in the activation of functional nitric-oxide reductase and Rubisco, respectively (1113). The VWA domain is a metal-binding domain often mediating protein-protein interactions (19). The metal is typically magnesium. Metal-binding occurs through a non-contiguous metal ion-dependent adhesion site, which is important for binding protein ligands (20). VWA proteins have been well studied in eukaryotes, where they are involved in a range of processes, including cell adhesion, transport, the complement system, proteolysis, transcription, DNA repair, and ribosome biogenesis. Bacterial and archaeal VWA proteins are not as well characterized, although limited work has identified members involved in surface adhesion, serum opacity, fibrinogen binding, and metal insertion into protoporphyrin IX (2123). Thus, to carry out their varied remodeling functions, many MoxR AAA+ proteins probably work in conjunction with VWA domain-containing proteins. As part of a research project aimed at identifying novel AAA+ chaperones in E. coli, we became interested in the uncharacterized protein YieN (SwissProt number P31473 [GenBank] ), which defines a MoxR subfamily (2). We have renamed this protein RavA (regulatory ATPase variant A). Here, we show that RavA is a cytoplasmic, hexameric AAA+ protein that functions with a VWA domain-containing protein, YieM (SwissProt number P03818), which we termed ViaA (VWA interacting with AAA+ ATPase). Detailed information is presented on RavA gene organization, expression patterns, protein subcellular localization, enzymatic activity, and oligomerization properties. Furthermore, we show that RavA strongly interacts with the inducible lysine decarboxylase (LdcI or CadA, SwissProt number P0A9H3), an enzyme important in bacterial acid stress response, to form a distinctive 5-fold symmetric cage-like complex as determined by electron microscopy.
BioinformaticsA total of 156 proteins belonging to the MoxR AAA+ family, as defined by the Clusters of Orthologous Groups (MoxR COG0714), were selected for analysis (24). These sequences included those currently identified in the COG data base as belonging to COG0714, as well as additional sequences extracted using the BLASTP program from 427 microbial genomes in the NCBI sequence data base (25). All additional sequences were analyzed using the COGNITOR program to ensure that they belonged to COG0714 (26). A multiple sequence alignment was then constructed using the ClustalW program (27). The phylogenetic tree was generated using the PROTDIST and FITCH programs included in the PHYLIP package (28). A PMB model of substitution (29) and global rearrangement were used. Default settings were used for all other parameters. 36 sequences identified as belonging to the RavA subfamily of MoxR AAA+ were selected and aligned using ClustalW. This alignment was then used to identify residues conserved in >90% of sequences. Conservation was based upon the following amino acid groupings: ILVM, KR, DE, ST, YFW, and NQ. Ungrouped amino acids included Ala, Cys, Pro, His, and Gly.
Bacterial StrainsWild-type MG1655 strain was purchased from ATCC (ATCC Number 47076). MG1655
The RavA-TAP-encoding strain was constructed using the Lambda Red system in a manner similar to that described for the construction of the knock-out strains. Specific details on strain construction, as well as on the pull-down and identification of bands by mass spectrometry, are provided elsewhere (33). MG1655 Gene CloningravA (yieN_ECOLI), ldcI (cadA/dcly_ECOLI), and ldcC (dclz_ECOLI) genes were cloned from E. coli K12 MG1655, whereas viaA was cloned from E. coli O157:H7. It should be noted that the amino acid sequences of ViaA from E. coli K12 MG1655 and O157:H7 strains are identical. Primers used were as follows: RAVA_forward (5'-atcacgatcatatggctcaccctcatttattag-3'), RAVA_reverse (5'-tacgtaggatccttagcattgttgtgcctggcg-3'), VIAA_forward (5'-aactcgatcatatgctaacgctggatacg-3'), VIAA_reverse (5'-ctatggatccttatcgccgccagcgtctgagc-3'), LDCI_forward (5'-atgctaccatatgaacgttattgcaatattg-3'), LDCI_ reverse (5'-ctggatccttatttcttgctttcttctttcaatacc-3'), LDCC_forward (5'-tctaccatgggcatgaacatcattgccattatg-3'), and LDCC_reverse (5'-cacctcgagttatcccgccatttttaggac-3'). ravA was cloned into p11 plasmid (34), a modified pET15b, which adds His6 followed by a tobacco etch virus-cut site to the N terminus of the expressed protein (HV tag). viaA and ldcI were cloned into pET3a (Novagen), which expresses untagged proteins. ldcC was cloned into pET16b such that it was expressed untagged, with the exception of additional N-terminal Gly and Met residues. All constructs were verified by DNA sequencing.
Protein Expression and PurificationProteins were expressed in BL21(DE3) pLysS (Stratagene). RavA, ViaA, and LdcI protein expression was induced by addition of 0.4 mM isopropyl 1-thio- RavA was purified using Qiagen nickel-nitrilotriacetic acid-agarose beads according to the manufacturer's protocols. The HV tag was removed by digestion with tobacco etch virus protease at 1:10 molar ratio of tobacco etch virus to protein for 5 h at 4°C. The protein was further purified using a Mono S HR 5/5 column. ViaA was precipitated from cell lysate using 30% ammonium sulfate. Pellets were resuspended and exchanged into a no-salt buffer using PD10 desalting columns. The protein was further purified on a Heparin HP5 column. Cells overexpressing LdcI were lysed in buffer LI1 (25 mM Tris-HCl, pH 7.5, 300 mM NaCl, 0.1 mM pyridoxal 5'-phosphate, 5% glycerol, and 1 mM DTT). Cell debris was removed by centrifugation. The supernatant was heated at 70 °C for 5 min resulting in the precipitation of the bulk of E. coli proteins, however, most of LdcI remained soluble. Precipitate was removed by centrifugation at 30,000 x g for 15 min. LdcI in the supernatant was further purified by ion exchange chromatography using a Mono Q HR 5/5 column and then by size exclusion chromatography using a Superdex 200 HR 10/30 column. LdcC was precipitated from cell lysate using 30% ammonium sulfate, and the pellets were resuspended and exchanged into low salt buffer by dialysis. The protein was further purified using a Mono Q HR 5/5 column. The final yields of purified proteins per 1 liter of culture were as follows: 6 mg for RavA, 2 mg for ViaA, 10 mg for LdcI, and 5 mg for LdcC. Proteins were judged to be >95% pure by SDS-PAGE analysis. All chromatography columns used were from Amersham Biosciences. Protein concentration was determined using the Bio-Rad Protein Assay. Measurement of RavA ATPase ActivityRavA ATPase activity was measured by determining the free phosphate released after ATP hydrolysis using a colorimetric assay similar to that reported previously (35). Briefly, RavA was placed in buffer RA1 (5 mM MgCl2, 0.02% Triton X-100, and 1 mM DTT) containing various concentrations of ATP substrate. RA1 solutions of different pH were generated using 10 mM concentrations of citrate, acetate, MES, HEPES, or TAPS buffers. Where applicable, BSA, ViaA, LdcI, and LdcC were included in the buffer. Assays were performed in 100-µl volumes at 37 °C in Eppendorf tubes or in a quartz plate. At different time points, 10-µl aliquots of the assay mixture were removed and added to 200 µl of RA2 (1 mM malachite green, 8.5 mM ammonium molybdate, and 1 M HCl). Color development was allowed to proceed for 1 min and then was stopped by the addition of 25 µl of 37% citric acid. Absorbance was measured at 660 nm and converted to moles of phosphate produced using a KH2PO4 standard curve.
Preparation of Cell Extracts for LdcI Activity AssaysWT and Measurement of LdcI Lysine Decarboxylase ActivityThe ability of LdcI to generate cadaverine from lysine was measured as previously described with slight modification (36). Purified LdcI was dissolved in buffer LA1 (2.5 mM MgCl2, 1 mM ATP, 0.1 mM PLP, and 4 mM lysine). LA1 solutions of different pH levels were generated using acetate, HEPES, or TAPS buffers at a concentration of 50100 mM. Total cell protein (25 µg/assay) was dissolved in buffer LA2 (25 mM MES, pH 6.0, 0.1 mM PLP, and 4 mM lysine). Assays were conducted in 100-µl volumes at 37 °C using either a quartz plate or Eppendorf tubes. At specific time points, 15-µl aliquots were removed and added to an equal volume of 1 M Na2CO3 stop solution. Subsequently, 15 µl of 10 mM 2,4,6-trinitrobenzenesulfonic acid (TNBS) was added, and samples were heated at 40 °C for 5 min to allow the formation of cadaverine-TNBS and lysine-TNBS adducts. Cadaverine-TNBS (N,N'-bistrinitrophenylcadaverine) was extracted using 500 µl of toluene, and the absorbance at 340 nm was measured. An extinction coefficient of 2.5 x 104 M1 cm1 was used (36) to convert absorbance into moles of cadaverine produced.
Growth in LdcI Indicator MediaE. coli K12 MG1655 WT,
Extreme Acid Shock AssaysExtreme acid shock assays were adapted from previously described protocols (37). Briefly, E. coli K12 MG1655 WT and
Northern Blot AnalysisE. coli MG1655 cultures were grown in LB media at 37 °C to an A600 of 0.60.8. Total RNA was isolated using the single step acid guanidinium thiocyanate-phenol-chloroform extraction method as described (38). Briefly, 1.5-ml samples of culture were centrifuged and resuspended in 500 µl of buffer NB1 (4 M guanidinium thiocyanate, 25 mM sodium citrate, 0.5% sodium lauryl sarcosinate, and 0.1 M 20 µg of total RNA was separated on a 1.2% agarose gel containing 6.7% formaldehyde. RNA was transferred to Hybond N+ Nylon membrane (Amersham Biosciences) and was covalently cross-linked to the membrane by UV irradiation using a Stratalinker UV Crosslinker (Stratagene). Probe generation and Northern blotting were carried out using the AlkPhos Direct Labeling and Detection System (Amersham Biosciences). For RavA, the DNA probe used was complementary to bp 6231497 of the ravA gene, whereas for ViaA, the DNA probe used was complementary to bp 7191452 of the viaA gene. RT-PCRRT-PCR analysis was performed using the Qiagen OneStep RT-PCR kit as described in the manufacturer's instructions. 2 µgof E. coli K12 MG1655 total RNA was used as a template. Amplification was carried out using primers RAVA_forward2 (5'-tatcacggcgccatggctcaccctcatttattag-3') and VIAA_reverse and primers RAVA_internal_forward (5'-gcgaattccatatgcaacaaattgatgtattgatgaccg-3') and VIAA_internal_reverse (5'-cgttcgatcactttttcacg-3'). Western Blot AnalysisTypically, cell pellets were lysed by sonication, and lysates were centrifuged to remove cellular debris. Supernatant was transferred to a fresh tube, and protein concentration was determined using the Bio-Rad Protein Assay. 540 µg of proteins was separated on a 10% SDS-PAGE gel and transferred to Biotrace polyvinylidene difluoride membrane using a Bio-Rad semidry transfer apparatus as per the manufacturer's instructions. Membranes were then incubated with rabbit polyclonal anti-RavA or anti-LdcI antibodies generated at the Division of Comparative Medicine, University of Toronto. Membranes were then washed and incubated with protein LA-peroxidase. ECL detection reagent (Amersham Biosciences) was used to visualize RNA bands on a Bioflex chemiluminescence detection film (Clonex). Subcellular Fractionation of E. coli200 ml of cells grown at 37 °C until an A600 of 2.5 was divided into 2 x 100 ml and centrifuged at 4000 x g for 10 min at 4 °C. Cell pellets were washed once with buffer SF1 (10 mM Tris-HCl, pH 7.1, and 30 mM NaCl) and then centrifuged again. One pellet was used for the isolation of perisplasmic and cytoplasmic proteins, and the other pellet was used for the isolation of membrane proteins. Periplasmic and cytoplasmic proteins were isolated using a protocol based on a previously described osmotic shock method (39), whereas total membrane proteins were isolated as previously described (40).
Size Exclusion ChromatographyA calibrated Superose 6 HR 10/30 column connected to an AKTA fast-protein liquid chromatography system (Amersham Biosciences) was used for size exclusion chromatography. Samples were typically run in buffer SE1 (25 mM Tris-HCl, pH 7.5, 300 mM NaCl, 2.5 mM MgCl2, and 1 mM DTT) in the presence or absence of 1 mM ATP. 0.1 mM PLP was added to the buffer when LdcI or LdcC were included in the sample. The column was calibrated using molecular mass standards (Sigma): thyroglobulin (669 kDa), apoferritin (443 kDa), Analytical Ultracentrifugation Analysis of RavASedimentation equilibrium experiments on RavA were carried out at the Ultracentrifugation Service Facility in the Department of Biochemistry at the University of Toronto. 20 µM RavA protein in Buffer AU1 (25 mM Tris, pH 7.5, 300 mM NaCl, 1 mM DTT, and 10% glycerol) containing 0.2 mM ATP was spun at 5000 and 8000 rpm, at 4 °C, in an An-60 Ti rotor in a Beckman Optima XL-A analytical ultracentrifuge. Absorbance was recorded at 230 and 280 nm. Data analysis was performed using the Origin MicroCal XL-A/CL-I Data Analysis Software Package Version 4.0.
Mapping of RavA Domains7.5 µg of purified RavA protein was dissolved in 140 µl of buffer RD1 (10 mM Tris, pH 7.5, 100 mM KCl, 10 mM MgCl2, 2.5 mM DTT, and 5 mM ATP) and digested with 0.075 µgof trypsin (Roche Applied Science, 1418475) at 4 °C. At specific time points aliquots of digested protein were removed and separated using SDS-PAGE. Sample preparation of bands for mass spectrometry was adapted from a previously described method (41). Briefly, the samples were processed in 500 µl of PCR grade tubes. The gel band was cubed and washed sequentially in high-performance liquid chromatography grade water (VWR, EM-WX0004-1) and 100 mM ammonium bicarbonate (Sigma, A6141). The gel cubes were dehydrated between the following steps with 100% acetonitrile (VWR, EM-AX0151-1): gel cubes were reduced with 1 mM DTT (Sigma D9163) at 50 °C for 30 min and alkylated with 55 mM iodoacetamide (Sigma, I1149) at room temperature for 20 min in the dark. Gel cubes were dehydrated with acetonitrile washes and SpeedVac-treated to remove acetonitrile. Protein was digested by swelling the cubes with 50 µl of 0.1 ng/µl trypsin in 50 mM ammonium carbonate and 0.1% calcium chloride overnight. The tryptic peptides were extracted from the gel slices by shaking with two washes of 150 µl of 100 mM ammonium carbonate buffer at 37 °C for 45 min. Liquid from the overnight digestion and two extractions were pooled on ice and acidified with 2 µl of 1% glacial acetic acid (VWR EM-AX0073-6). Peptides were desalted and concentrated using C18 resin (Sigma, H8261) by a batch method. Final volume was 10 µl of extracted tryptic peptides in 65% acetonitrile and 1% acetic acid. A 1-µl aliquot of tryptic peptide sample was mixed with 1 µl of a saturated High mass samples were analyzed in a Bruker Reflex III (Bruker Daltonics) in high mass linear mode with a high mass detector. Carbonic anhydrase (Sigma, C3934) and lysozyme (Sigma, L7651) were used as external calibrants. Tryptic peptides were analyzed in a Bruker Reflex III (Bruker Daltonics) in reflector mode. Samples were internally calibrated using trypsin peptide peaks. Mass mapping was done using PAWS (Genomic Solutions) against the known sequence of the protein. The tryptic peptide fingerprint of a gel band was used to determine an approximate region for the partial proteolysis fragment from the full sequence. The observed high mass was then used to define the exact cut sites for the partial proteolysis fragments.
Electron MicroscopyFor preparation of negatively stained protein samples, 3 µl of protein solutions was applied to glow-discharged continuous carbon copper grids and left for 1 min before washing with two drops of 1% (w/v) uranyl acetate. The grids were observed with a Philips CM 12 transmission electron microscope with an LaB6 filament at 80 kV. Images were recorded on Kodak SO-163 films at 35,000x magnification in the case of RavA, LdcI, and ATP
Image AnalysisFor LdcI, 3631 subframes of 88 x 88 pixels containing single particles were extracted interactively from several micrographs at similar defocus of
For ADP·RavA·LdcI and ATP
RavA Represents a Distinct Subfamily Within the MoxR AAA+ FamilyInitially, phylogenetic analysis of the MoxR AAA+ family was carried out to determine the major subfamilies within this family. A phylogenetic tree was generated using 156 MoxR AAA+ sequences, which belong to MoxR Clusters of Orthologous Groups 0714 (COG0714) (National Center for Biotechnology Information, National Institutes of Health) (24). The tree reveals the presence of six distinct subfamilies: MoxR-Proper, CGN, RavA, YehL, APE0892, and PA2707 (Fig. 1A). The first five of these groups are consistent with those reported earlier (2). None of the proteins of the RavA, YehL, APE0892, and PA2707 subfamilies have yet been characterized; however, limited work with members of the MoxR-Proper and CGN subfamilies has suggested that the MoxR AAA+ proteins function as molecular chaperones. Closer examination of the 36 genes in the RavA subfamily reveals that 33 of these occur immediately adjacent to genes encoding proteins containing VWA domains (Fig. 1B). The co-occurrence of MoxR AAA+-containing genes with those encoding proteins containing VWA domains is a common feature within the MoxR family as a whole (data not shown) but appears with particularly high frequency among the RavA subfamily. This suggests a functional linkage between the RavA AAA+ proteins and their VWA-containing counterparts. In the case of E. coli RavA, the VWA-containing protein is YieM, which we have renamed ViaA. Examination of a ClustalW (27) alignment of all 36 RavA sequences revealed the presence of a number of highly conserved signature sequences in the E. coli K12 MG1655 RavA sequence (marked in red in Fig. 2). Upon comparison to common AAA+ signature sequences (1), it is observed that many of these conserved residues distinguish RavA proteins from other AAA+ proteins. At the N terminus, the Box II motif contains two highly conserved residues, an Arg and hydrophobic residue (Leu in the RavA sequence). The Walker A motif contains a number of highly conserved residues, but most of these are common to all AAA+ proteins. One feature of note, however, is that members of the RavA subfamily favor an alanine residue in place of the more commonly observed glycine of the classical "GK(S/T)" motif. A series of highly conserved residues also exist in the Box IV and Box IV' motifs. Based upon comparison to published alignments, at least part of the sequence in this region may lie on a small "insert" in the second helix of the AAA+ core module (2). The Walker B motif also contains a number of highly conserved residues. Most of these are common to all AAA+ proteins, however, the conserved tryptophan/phenylalanine residue is particularly unique to RavA and may have important functional significance to RavA proteins. A unique, highly conserved LP motif also exists immediately prior to the Walker B motif.
Within Box VI, which follows the Walker B, members of the MoxR family possess a highly conserved LLXX
A conserved LXAXYDRXXXR motif exists in Box VII-Box VII' region. This region contains the conserved arginine finger of AAA+ proteins, believed to be responsible for mediating intersubunit communication in response to ATP hydrolysis (5). The first Arg residue in the conserved motif may correspond to the arginine finger. Some additional conserved hydrophobic residues also exist in Box VII' and in the sequence following Box VII''. The Sensor 2 region contains a short stretch of conserved residues, SDRR. One of the arginines is expected to correspond to the conserved Sensor 2 residue present throughout most of the AAA+ superfamily. This residue is believed to be responsible for sensing ATP and mediating movements between the AAA+ core domain and small
RavA Is a Soluble Cytoplasmic ATPase Consisting of Two Domains To experimentally determine the domains present in RavA structure, we expressed and purified the protein from E. coli. The purified protein was then treated with limited amounts of sequencing grade trypsin for varying amounts of time (Fig. 3A). The boundaries of the fragments generated were then determined using mass spectrometry (Fig. 3B). Based on this analysis, the protein can be divided into two main domains: a highly conserved N-terminal domain (Met1Arg347), which corresponds to the AAA+ module, and a poorly conserved C-terminal domain (Arg348Cys498) that is predicted to be predominantly
Subsequently, the subcellular localization of RavA was determined. Wild-type and ravA cells grown at 37 °C to early stationary phase were lysed and separated into cytoplasmic, periplasmic, and crude membrane fractions. Western blot analysis of these fractions using RavA antibodies (Fig. 3C) clearly showed that RavA is predominantly cytoplasmic. This is consistent with the lack of any predicted signal sequences or transmembrane regions in RavA. A faint band observed in the membrane fractions appears to be nonspecific, as it appears in both the WT and ravA strains. As shown in Fig. 3D, purified RavA at 0.4 µM has the highest ATPase activity at neutral to alkaline pH, with the highest activity around pH 7.5. Below pH 7.5, the activity appears to drop rapidly and is minimal below pH 6.0. Complete Michaelis-Menten kinetic analysis was performed at pH 7.5 (Fig. 3D, inset). The kcat and Km values for RavA ATPase are in a similar range to those obtained for other AAA+ proteins such as HslU (49), NtrC (50), and Lon (51). Furthermore, RavA also has GTPase activity. GTP hydrolysis appears to be much slower than ATP hydrolysis at lower enzyme/substrate concentrations but reaches comparable levels with increasing amounts of enzyme/substrate (data not shown). This effect may be the result of a lower ability of GTP versus ATP to promote RavA oligomerization (data not shown).
ravA and viaA Genes Form an Operon under the Direct Control of S The ravA and viaA gene sequences encode predicted proteins of similar length consisting of 498 and 483 amino acid residues, respectively. The genes reside on the E. coli chromosome close to the origin of replication (OriC) (52), 2200 bp clockwise from OriC, and the genes are transcribed in the counterclockwise direction. Interestingly, the 5'-end of the ViaA coding region overlaps the 3'-end of the RavA coding sequence by seven bases (Fig. 4A). Such an overlap might hint to the presence of an ancient gene whose sequence was a fusion of the ravA and viaA gene sequences.
To determine whether ravA and viaA are on the same operon, Northern blot analysis was carried out on E. coli MG1655 wild-type (WT), MG1655
These data suggest that ravA and viaA exist as a single operon, consistent with the lack of a predicted viaA promoter sequence. The full-length 3100 bases transcript appears to be inherently unstable, however, and undergoes at least two distinct internal cleavage events shown schematically in Fig. 4E. The first of these cleavages produces the 1900-base fragment observed in the ravA blot (Fig. 4C, band 2 in the first lane of panel I) and the 1200-base fragment observed in the viaA blot (Fig. 4C, band 2 in the first lane of panel II). The second cleavage produces the 1100-base fragment observed in the ravA blot (Fig. 4C, band 3 in the first lane of panel I) and the 2000-base fragment observed in the viaA blot (Fig. 4C, band 1 in the first lane of panel II). Whether or not these cleavage events have a regulatory role is currently under investigation.
Northern blot analysis of MG1655
Analysis of the sequence immediately upstream of the ravA gene reveals the presence of a highly conserved S promoter consensus sequence (Fig. 4A) (53). This consensus sequence has been shown to be present in genes that are up-regulated in stationary phase. Consistent with this, Western blot analysis of MG1655 cells shows a significant induction of RavA levels as cells approach late log/early stationary phase (Fig. 4F, left panel). Such an induction does not occur in an MG1655 rpoS strain, although constant levels of RavA protein appear to be present at all times (Fig. 4F, right panel). The data suggest that RavA induction is partially dependent on the stationary phase sigma factor S, likely in a direct manner. The presence of low, constitutive levels of RavA protein in the rpoS strain suggests that the ravA promoter may display conditional selectivity allowing its expression to be regulated by other sigma factor(s) (e.g. 70). Such conditional promoter selectivity has been reported for other E. coli genes, for instance dps and yaiA, whose expression may be regulated by either S or 70, depending on cellular conditions (53).
RavA Forms a Hexameric OligomerMany AAA+ proteins typically function as oligomeric rings (6); therefore, the nature of the oligomeric state of RavA was investigated using size exclusion chromatography and analytical ultracentrifugation. Size exclusion chromatography analysis using 3 µM RavA revealed that the protein forms oligomeric species and that this oligomerization is enhanced in the presence of nucleotide (Fig. 5A). The largest species formed runs near the 443-kDa marker. To more accurately confirm the size of this species, analytical ultracentrifugation sedimentation equilibrium studies were performed in the presence of ATP. The best fit to the data set corresponds to a monodisperse species of Negative stain electron microscopy of RavA in the presence of ADP, ATP, or non-hydrolyzable ATP analogues invariably showed a mixture of different oligomeric states of the protein, ranging from monomers to bigger oligomers. The top-view appearance of the largest oligomers clearly favored their hexameric composition (Fig. 5C). However, due to the inhomogeneity of RavA oligomeric states on the grid and the resulting difficulty in unambiguous identification of corresponding side views, no attempts have been made to obtain a three-dimensional reconstruction of the RavA hexamer (however, see below). RavA Binds to the Inducible Lysine Decarboxylase LdcI/CadAIn an effort to determine the function of RavA, experiments were carried out in an attempt to identify its interacting partners. E. coli DY330 strain (30, 31), in which a TAP tag (54) was fused C-terminally to the endogenous ravA chromosomal gene, was used for pull-down experiments (33). When RavA-TAP was isolated from cells grown to stationary phase in TB-rich media at 37 °C, the inducible cytoplasmic lysine decarboxylase (LdcI) was found to be strongly bound to RavA (Fig. 6A). LdcI (also termed CadA, SwissProt name DCLY_ECOLI, SwissProt number P0A9H3) is a pyridoxal phosphate (PLP, a vitamin B6 derivative)-dependent decameric enzyme (55, 56), which plays a major role in the acid stress response in bacteria (5759). LdcI decarboxylates a lysine to produce the polyamine cadaverine, which is then pumped out of the cell by an inner membrane lysine/cadaverine antiporter protein, CadB (59, 60). This process removes a proton from the cell cytoplasm and reduces the acidification of the intracellular environment. The genes cadB and ldcI/cadA are on the same operon (60, 61), and the operon is activated by an inner membrane anchored transcription activator, CadC, that binds a Pcad element in the promoter region of cadBA (61). The regulation of cadBA transcription is rather complex. The operon is induced by low external pH, by the presence of exogenous lysine, and under anaerobic conditions (55, 62). The binding of RavA to LdcI leads us to speculate that RavA might also be involved in the bacterial stress response or that the interaction between the two proteins is regulatory in nature.
Analysis of the binding of RavA to LdcI was carried out using size exclusion chromatography (Fig. 6, B and C). All experiments shown in Fig. 6B were performed in the presence of 1 mM ATP. 18 µM RavA protein migrates as a single species with an apparent molecular mass of 440 kDa representing the hexameric species (Fig. 6B, panel 1, top panel). 10 µM LdcI (Fig. 6B, panel 2) also elutes as a single species, with an apparent molecular mass of 600 kDa probably representing the decameric species that is known to be formed by LdcI (also see below) (55, 56). A mixture of 10 µM LdcI with excess RavA (18 µM) results in a large shift in the migration of both species indicating a strong interaction (Fig. 6B, panel 3). A range of higher order oligomers are formed, with the predominant complex (boxed) migrating at or near the void volume of the column, suggesting a size of 50 MDa or greater. The predominant complex formed changes if 18 µM RavA is incubated with 2-fold excess of LdcI (40 µM). The formed complex elutes within the column volume, and its molecular mass is 35 MDa as estimated by sedimentation equilibrium experiments (Fig. 6B, left box in panel 4 and data not shown). We term this complex the fraction three complex (FTC), which, as discussed below, might be the physiologically relevant RavA·LdcI complex. E. coli contains another protein that is a very close paralogue of LdcI. This protein is the constitutive lysine decarboxylase LdcC (SwissProt name DCLZ_ECOLI, SwissProt number P52095 [GenBank] ). LdcI and LdcC are 69% identical and 84% similar as determined using pairwise BLAST P analysis. LdcC is expressed in E. coli at very low levels, and its exact role in the cell remains unclear (63, 64). Like LdcI, LdcC is also a PLP-dependent enzyme that decarboxylates lysine and possibly assembles into a decamer (63, 64) (Fig. 6B, panel 5). However, the enzymatic activity profile and stability of LdcC differ from that of LdcI (63, 64). Due to the high sequence similarity between LdcC and LdcI, we expected that LdcC might also bind RavA, however, we were surprised to find that there is no significant interaction between RavA and LdcC (Fig. 6B, panel 6). The specific interaction of RavA with LdcI and not LdcC possibly indicates that RavA is a stress response protein.
The size of the LdcI oligomer (decamer) is unaffected by the presence or absence of ATP (Fig. 6C, panel 1, upper panel) or by LdcI concentration (compare elution profile of 10 µM LdcI in Fig. 6B, panel 2, to the elution profile of 1.5 µM LdcI in Fig. 6C, panel 1). As shown in Fig. 6C, the size of the RavA·LdcI complex is dependent upon RavA oligomerization. 3 µM RavA in the absence of ATP predominantly elutes at Similar interactions were observed between RavA and LdcI when ADP was added to the buffer rather than ATP. Furthermore, the inclusion of lysine or cadaverine or the removal of the LdcI cofactor PLP from the buffer did not affect the interaction between RavA and LdcI (data not shown). These observations clearly indicate that the interaction between RavA and LdcI is highly robust and functionally significant.
RavA·LdcI Forms a Very Large Distinctive Cage-like Complex as Visualized by Negative Stain Electron MicroscopyA 20 Å resolution structure of the D5-symmetric LdcI oligomer (Fig. 7, AC) was obtained as described under "Materials and Methods." The existence of a 5-fold symmetry axis in the particle was corroborated by multivariate statistical analysis. Fig. 7B shows the pronounced 5-fold symmetrical shape of the top-view class average and the 2-fold symmetry of the side-view class average. Some dimeric dissociation products could always be detected in the background of the micrographs (Fig. 7A), in agreement with the early publication of Sabo et al. (55). The random orientation of single particles on the carbon film made this specimen particularly suitable for a three-dimensional reconstruction by angular reconstitution approach (Fig. 7C). The decamer of LdcI revealed itself as two pentameric rings, tightly stacked together back-to-back and turned
Grids prepared by mixing ADP·RavA or ATP
Hence, the RavA·LdcI complex is a unique complex in which five RavA oligomers (at most) bridge two LdcI decamers resulting in a distinctive cage-like structure. Although the volume of the densities attributed to RavA would be consistent with its inclusion in the complex as a hexameric species, the current quality of the reconstruction does not allow us to claim this with absolute certainty. Remarkably, no differences could be observed between the RavA·LdcI complex in ADP and the one in ATP S (data not shown), at least at the current resolution estimated to be 30 Å. We believe that the complex seen in Fig. 7 corresponds to the FTC in Fig. 6.
Binding of RavA to LdcI Affects RavA Activity but Not LdcI Activity The formation of a tight complex between RavA and LdcI either indicates that RavA affects LdcI function or vice versa. Because RavA is an AAA+ protein, we expected that it would remodel LdcI in a nucleotide-dependent manner and, hence, affect the activity of the enzyme. LdcI activity was measured at a range of different pH values in the presence and absence of RavA protein (Fig. 8A). LdcI activity is greater at neutral pH than at alkaline pH, with highest activities observed at pH 6.0 and 7.0. At pH 8.0 and 9.0, LdcI activity is virtually undetectable. Surprisingly, the activity of the LdcI enzyme was unaffected by the presence of RavA under any of the conditions tested even at higher concentrations of RavA and LdcI or under different incubation conditions (data not shown). The inclusion of ViaA protein in the assay also failed to produce any alteration in the activity of the LdcI enzyme (data not shown).
Furthermore, LdcI activity in whole cell extracts was also not affected by deletion of the ravA gene (Fig. 8B), in which the enzyme was induced by growth in unbuffered LB media containing excess glucose. The utilization of glucose as a carbon source results in the decrease in the pH of the media as cells ferment the glucose and produce lactic and acetic acids (65). The addition of purified RavA and nucleotides to the extract also did not affect LdcI activity (data not shown). The lack of effect of RavA on LdcI activity in vivo was further confirmed by comparing the ability of WT and
The ability of WT and The effect of LdcI on RavA ATPase was measured. RavA ATPase activity was assessed using 1 mM ATP substrate in the presence of ViaA, LdcI, LdcC, and BSA (Fig. 8F). None of the added proteins appear to possess significant ATPase activity of their own (columns 14). The activity of 0.25 µM RavA alone is shown (column 5). Addition of 0.5 µM BSA or 0.5 µM LdcC did not have any significant effect on RavA ATPase activity (columns 6 and 7) consistent with the lack of binding of BSA or LdcC to RavA. Addition of 0.5 µM LdcI results in 1.4-fold stimulation of RavA ATPase (column 8). Increasing the concentration of LdcI did not result in further stimulation of RavA ATPase (data not shown). This observed increase in activity can be explained in several ways. One possibility is that RavA may exist in a more active form upon binding LdcI; alternatively, LdcI may help stimulate RavA hexamer formation in the presence of nucleotide by acting as a scaffold, with a concomitant increase in the overall ATPase activity of the population. Addition of ViaA results in a 1.9-fold stimulation of RavA ATPase activity (column 9). RavA activity in the presence of both 0.5 µM LdcI and 0.25 µM ViaA was stimulated by 2.4-fold (column 10). Because the fold stimulation obtained is not an average of 1.4 and 1.9, this might indicate the possible formation of a ternary complex of RavA·ViaA·LdcI. Alternatively, in the presence of LdcI and ViaA, RavA exists in a different conformation or oligomeric state than in the presence of LdcI or ViaA alone. It should be noted that the stimulations observed in the above experiments are reproducible with different preparations of proteins. Based on the above data, our current speculation is that LdcI functions to regulate RavA and that the RavA·LdcI complex is formed under certain conditions to modulate RavA function. Although we have identified a putative regulator of RavA, efforts are currently underway to determine the exact role of RavA and ViaA in the cell and the substrates targeted by this AAA+ and VWA system.
RavA·ViaA Might Play a Role in Metal InsertionOur analysis of RavA reveals that it is a member of a distinct subfamily of the MoxR AAA+ proteins whose function remains poorly characterized. We were able to divide the MoxR proteins into six subfamilies (Fig. 1A). Intriguingly, despite the distant relationship between members of the MoxR-proper and CGN subfamilies, there appears to be a remarkable similarity in the nature of their function. The parallels between the importance of the CGN family members NorQ/NirQ and CbbQ for the production/activation of nitric-oxide reductase and Rubisco, respectively, and the importance of MoxR for the production of active methanol dehydrogenase are striking. Even the GvpN protein, although not important for the biogenesis of an enzyme, appears to be crucial for the production of the large proteinaceous structure that comprises gas vesicles in certain cyanobacteria and halophilic archaea. Thus, it does not seem unreasonable to presume that RavA proteins, and indeed members of the other uncharacterized MoxR subfamilies, may play similar roles. The co-occurrence of genes encoding proteins with VWA domains with MoxR AAA+ genes is also notable. We have demonstrated that RavA and its corresponding VWA protein, ViaA, comprise an operon (Fig. 4). This is consistent with results reported for the NorQ/NirQ and CbbQ proteins, whose genes are also part of an operon containing immediately adjacent VWA protein-encoding genes. These VWA domain-containing proteins were also shown to be important in the biogenesis of active nitric-oxide reductase and Rubisco. Although the exact role of ViaA is unclear, one possibility is that ViaA and RavA function in a role similar to that of the metal-chelatases, a family of AAA+ proteins to which the MoxR family has been shown to be closely related (2, 66). In these metal-chelatases, the AAA+ components work in conjunction with VWA proteins to mediate the insertion of Co+2 or Mg+2 into porphyrin rings as part of the synthesis of cobalamin or (bacterio)chlorophyll, respectively (47). In the case of Mg-chelatase, one of the subunits, BchD, actually contains an AAA+ and VWA domain fused together on a single polypeptide, the only example of such an arrangement known in bacteria and archaea (19). It is intriguing to speculate that RavA and ViaA may, therefore, play a role in metal insertion, a process that may be important in the activation of their substrate(s). Interestingly, both the nitric-oxide reductase for which NorQ/NirQ are important and the Rubisco for which CbbQ is important depend upon metal ions for their function. Efforts are currently underway to explore this possibility. Possible Implications of LdcI Binding to RavAThe interaction of RavA with LdcI is interesting, however, its exact role is puzzling. RavA does not appear to exert any effect upon the activity of LdcI under a range of different in vitro conditions (Fig. 8A). In addition, cells in which the RavA gene has been disrupted display normal LdcI induction and activity and do not appear to be compromised in their response to a variety of acid stresses (Fig. 8, BD). RavA levels are also not induced by acid, and the protein is present under growth conditions in which the LdcI enzyme is not present (Fig. 8E). The interaction of LdcI with RavA does result in an increase in RavA ATPase activity, however, suggesting that the complex is indeed functionally important. A further increase is observed in the presence of ViaA, suggesting ternary complex formation. The interaction also appears to be highly specific to LdcI, because RavA does not interact with the closely related LdcC protein. One possibility is that the complex formed between RavA and LdcI represents a means of regulating RavA activity in response to acid stress conditions. Thus RavA can be envisioned to carry out an independent function of some form, possibly helping to mediate the folding/assembly of particular enzyme complexes during stationary phase. Upon exposure to acid stress conditions, levels of LdcI increase dramatically, leading to the formation of the RavA·LdcI complex. This complex may represent a way of sequestering RavA protein from its natural substrates in response to acid stress. Alternatively, the complex may provide a means of maintaining RavA function under acid stress conditions, something more in line with the observed stimulation of RavA activity. By acting as a "scaffold," LdcI may stabilize RavA oligomers under acid stress conditions. It is also possible that the complex formed may lead to a "gain of function," directing RavA toward entirely new substrate proteins via localization and/or conformational changes. Implications of the RavA·LdcI Complex StructureEM analysis has provided us with an exceptional view of both the LdcI enzyme alone and in complex with RavA. The structure of the LdcI particle, deduced from the presented electron microscopy data (Fig. 7, AC), confirms the suggestion of Sabo et al. (55) that the LdcI oligomer is composed of five stable dimers. The characterization of LdcI in that work had indeed indicated that the protein was dimeric at pH 8.0 and at low ionic strength, whereas lowering of the pH to 7.0 and a concomitant increase in ionic strength favored an association of dimers into decamers. The dimer-decamer equilibrium was shown to be dependent on protein concentration. These observations explain some presence of dimers in our experimental conditions. It should be noted that, based on a visual analysis of individual side views, Sabo et al. (55), proposed the orientation of the dimer to be roughly perpendicular to the plane of the ring. On the contrary, our three-dimensional reconstruction unambiguously demonstrates a pronounced rotation of the monomer inside the dimeric building block of the LdcI double-toroid. The cage-like structure of the imaged RavA·LdcI complex makes it tempting to speculate that it may provide a space in which another protein or proteins might be confined. Proteins held in such a position might then be subject to remodeling by the RavA oligomers. Such a complex could thus represent a novel chaperone structure assembled under acid stress conditions when levels of LdcI enzyme increase dramatically.
In summary, considerable information has been gathered on the RavA protein and its interaction with ViaA and LdcI. Its regulation by
* This work was supported in part by a grant from the Natural Sciences and Engineering Research Council of Canada (to W. A. H.). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
1 Recipient of a graduate fellowship from the Natural Sciences and Engineering Research Council of Canada.
3 Supported by the Ontario Genomics Institute and Genome Canada. 2 To whom questions regarding EM methodology should be addressed: IVMS, c/o EMBL Grenoble Outstation 38042 Grenoble Cedex 9, France. E-mail: gutsche{at}emblgrenoble.fr. 4 A Canadian Institutes of Health Research New Investigator. To whom correspondence should be addressed: 1 King's College Circle, Medical Sciences Bldg., Dept. of Biochemistry, University of Toronto, Toronto, Ontario M5S 1A8, Canada. Tel.: 416-946-7141; Fax: 416-978-8548; E-mail: walid.houry{at}utoronto.ca.
4 The abbreviations used are: Rubisco, ribulose-bisphosphate carboxylase/oxygenase; VWA, Von Willebrand Factor Type A; LdcI, inducible lysine decarboxylase; DTT, dithiothreitol; MES, 4-morpholineethanesulfonic acid; TAPS, 3-{[2-hydroxy-1,1-bis(hydroxymethyl)ethyl]amino}-1-propanesulfonic acid; BSA, bovine serum albumin; WT, wild-type; PLP, pyridoxal phosphate; TNBS, 2,4,6-trinitrobenzenesulfonic acid; RT, reverse transcription; EM, electron microscopy; ATP
5 D. M. Belmap, J. F. Conway, and J. B. Heymann (2004) www.niams.nih.gov/rcn/labbranch/lsbr/software/pft2_em3dr2/pft2_em3dr2.htm.
We thank Usheer Kanjee, Asad Merchant, Bharat Sharma, and Philip Wong for contributing to different aspects of this project. J. S. and W. A. H. thank Dr. L. Aravind at the National Center for Biotechnology Information, National Institutes of Health, Bethesda, MD for his help with the bioinformatic analysis.
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||