4Q7Q is a protein found in Chitinophaga pinensis, a soil bacterium in the sphingobacterial family. Its structure has been previously characterized and exists in Protein Data Bank. Its function, however, has not.
After structural and sequential analysis via various databases including BLAST, Pfam, Dali, PyMOL, and ProMOL, we initially predicted that 4Q7Q is a hydrolase. More specifically, we hypothesized that 4Q7Q is a lipase, an enzyme that can hydrolyze lipids to form fatty acids and a glycerol molecule.
Characteristics of 4Q7Q
Figure 1: Plasmid Map of pMCSG73, the expression vector for 4Q7Q.
4Q7Q is composed of two chains; one chain can be seen . The colors indicate the translation direction of the peptide sequence from the N to C terminus; red represents the N-terminus while dark blue represents the C-terminus. Based on this structural model, we can see that 4Q7Q is an alpha-beta superfold; there are beta sheets (represented by the straighter strands) sandwiched between the alpha helices (represented by the coiled strands).
Using SnapGene, we analyzed the plasmid pMCSG73, which is the expression vector used to produce 4Q7Q. We found that the plasmid's open reading frame (ORF), the sequence of DNA that will be eventually translated, is 2237 base pairs in length. The size of 4Q7Q is 87.1 kDa.
Overall, as seen in Figure 1, the ORF is composed of NusA (a transcription elongation factor, two tags- 6XHis and Strep-Tag II, and a TVMV and TEV site. The plasmid map indicates the presence of a lac operon composed of three genes (lacZ, lacY, and lacA), a promoter (region that initiates transcription), terminator (region that marks the end of transcription), regulator (which codes for a repressor), and operator (repressor binds to it to prevent gene expression).
In silico Analysis
Figure 2: Alignment of 4Q7Q's putative catalytic triad (shown in green) and 3LIP's catalytic triad (shown in blue). The RMS is 2.257.
Figure 3: Alignment of 4Q7Q's putative catalytic triad (shown in green) and 1TAH's catalytic triad (shown in pink). The RMS is 2.205.
Figure 4: Alignment of 4Q7Q's putative catalytic triad (shown in green) and 1BWR's catalytic triad (shown in pink). The RMS is 2.049.
We initially analyzed 4Q7Q through the protein structure databases BLAST, Pfam, and Dali. The top hit was 4M8K, a GDSL-like lipase, a type of a lipase that has a flexible active site and therefore broad substrate specificity. Through BLAST, we found that 4M8K and 4Q7Q had a 36% sequence identity, with an E value of 0.002, indicating that it is a significant match. Since we can use the principle of homology to predict the function of an unknown protein, we first hypothesized that 4Q7Q was too a lipase.
Through analyzing the sequence of 4Q7Q in SnapGene and then analyzing the 3D structure in PyMOL, we hypothesized that a possible catalytic triad of 4Q7Q is Ser164, Asp193, and His196. We believe that this group of amino acids are involved in the active site of 4Q7Q and therefore affects how the protein works. As seen in this , all three amino acids are close in proximity to one another and are brought together in a single orientation.
We also performed further analysis in PyMOL and ProMOL which involved the homology of active sites. Top hits included 3LIP, a lipase found in Burkholderia cepacia (a human pathogen that can cause pneumonia), 1TAH, a lipase found in Burkholderia glumae (a soil bacterium), and 1BWR, a hydrolase found in Bos taurus (cattle). We aligned the putative catalytic triad of 4Q7Q with each of the catalytic triads of these known proteins.
3LIP has two chains. As seen in Figure 2, when aligning the catalytic triad of 3LIP (Asp264, Ser87, His286) to the putative catalytic triad of 4Q7Q, the RMS is 2.257.
1TAH has four chains. As seen in Figure 3, when aligning the catalytic triad of 1TAH (Asp263, Ser87, His285) to the putative catalytic triad of 4Q7Q, the RMS is 2.205.
1BWR has one chain. As seen in Figure 4, when aligning the catalytic triad of 1BWR (Asp192, Ser47, His195) to the putative catalytic triad of 4Q7Q, the RMS is 2.049.
Compiling all of the data together, we can see that 1BWR's catalytic triad is most structurally similar to the putative catalytic triad of 4Q7Q due to the lower RMS value measured. Therefore, we hypothesized that 4Q7Q is most likely a hydrolase; through experiments, we can investigate further if 4Q7Q is specifically a lipase.
Bacterial Transformation and Plasmid Purification
Figure 5: 4Q7Q bacterial transformation results.
Before characterizing the function of 4Q7Q, we first needed to synthesize the protein through first transcribing 4Q7Q's DNA to amplify it and then translating it to express it. First, 4Q7Q's DNA was transcribed using its expression vector, the plasmid pMCS573. Since transformation must occur within a cell, the plasmid was transformed into DH5α cells using New England Biolabs protocol.
After transformation, DH5α cells were lysed and spread on plates containing LB and ampicillin. Since 4Q7Q's plasmid is ampicillin resistant, only the bacteria that have been transformed with 4Q7Q's plasmid will grow on the plates. 4Q7Q's plasmid was then purified using Zyppy Plasmid Miniprep Kit.
However, although DH5α cells maximize the efficiency of transformations, they do not contain T7 polymerase, which is essential for protein expression. Therefore, the purified plasmid underwent another bacterial transformation into BL21 (DE3) cells that do contain T7 polymerase using New England Biolabs protocol.
Protein Expression and Measuring Plasmid Concentration
Figure 6: Top: how lac operon works in the absence of IPTG. Bottom: how the lac operon works in the presence of IPTG.
After the second bacterial transformation, BL21 cells were lysed and spread on LB+amp plates. A bacterial colony was then selected from a plate and suspended in liquid culture. After incubation overnight, the OD260, which essentially measures the concentration of the plasmid DNA, was measured with a biophotometer.
From our OD260, we calculated that our plasmid DNA concentration was 28.7 (μg/ mL). This concentration is low, probably due to the large size of our ORF.
The sample was then induced with IPTG. As seen in Figure 6, IPTG is a reagent that prevents the repressor from binding to the operator to allow expression to occur.
Protein Purification
Figure 7: SDS-PAGE gel. Bands are present at 87.1 kDa.
4Q7Q was then purified using the HisPUR Ni-NTA Purification kit. A nickel column as well as equilibration, wash, and elution buffers were used. We then tested for expression using SDS-PAGE, an electrophoresis method that separates proteins by mass in a polyacrylamide gel. BioRad's mini protean tetra protocol was utilized for SDS-PAGE.
Lanes in our gel contain samples of the cell extract, flow through, the third step of washing our protein through the column, as well as all three elutions. We want to see everything– our protein and all extra nucleic acids– in the cell extract lane. In the flow-through lane, we do not want to see the protein being expressed. In the wash lanes, we want to be able to see a little bit of the protein; each wash step should show more protein, as the buffer removes more and more unwanted nucleic acids. We are most interested in seeing our protein eluted in the three elution lanes. As seen in Figure 7, there are bands in the three elution lanes corresponding to our protein's weight: 87.1 kDa, confirming that we have successfully expressed 4Q7Q.
Measuring Protein Concentration
We measured the absorbance of a purified 4Q7Q sample (not used in the gel) at 280 nm in a spectrometer. We found that the A280 is 0.886.
In order to calculate the concentration of 4Q7Q, we had to use Beer's Law, shown in figure 8. Using SnapGene, we found that ε for 4Q7Q was 45840 M-1cm-1. The cuvette in our spectrometer was 1 cm in width, therefore path length was also 1 cm.
Through plugging in these values, we found that the concentration of our purified protein is 19.3μ.
pNPB Lipase Assay
We performed a lipase assay using p-nitrophenyl butyrate,
