http://orcid.org/0000-0003-3872-1869
Figures
Abstract
Proteases play an essential role in living organisms and represent one of the largest groups of industrial enzymes. The aim of this work was recombinant production and characterization of a newly identified thermostable protease 1147 from thermophilum indigenous Cohnella sp. A01. Phylogenetic tree analysis showed that protease 1147 is closely related to the cysteine proteases from DJ-1/ThiJ/PfpI superfamily, with the conserved catalytic tetrad. Structural prediction using MODELLER 9v7 indicated that protease 1147 has an overall α/β sandwich tertiary structure. The gene of protease 1147 was cloned and expressed in Escherichia coli (E. coli) BL21. The recombinant protease 1147 appeared as a homogenous band of 18 kDa in SDS-PAGE, which was verified by western blot and zymography. The recombinant protein was purified with a yield of approximately 88% in a single step using Ni-NTA affinity chromatography. Furthermore, a rapid one-step thermal shock procedure was successfully implemented to purify the protein with a yield of 73%. Using casein as the substrate, Km, and kcat, kcat/Km values of 13.72 mM, 3.143 × 10−3 (s-1), and 0.381 (M-1 S-1) were obtained, respectively. The maximum protease activity was detected at pH = 7 and 60°C with the inactivation rate constant (kin) of 2.10 × 10–3 (m-1), and half-life (t1/2) of 330.07 min. Protease 1147 exhibited excellent stability to organic solvent, metal ions, and 1% SDS. The protease activity was significantly enhanced by Tween 20 and Tween 80 and suppressed by cysteine protease specific inhibitors. Docking results and molecular dynamics (MD) simulation revealed that Tween 20 interacted with protease 1147 via hydrogen bonds and made the structure more stable. CD and fluorescence spectra indicated structural changes taking place at 100°C, very basic and acidic pH, and in the presence of Tween 20. These properties make this newly characterized protease a potential candidate for various biotechnological applications.
Citation: Tarrahimofrad H, Meimandipour A, Arjmand S, Beigi Nassiri M, Jahangirian E, Tavana H, et al. (2020) Structural and biochemical characterization of a novel thermophilic Coh01147 protease. PLoS ONE 15(6): e0234958. https://doi.org/10.1371/journal.pone.0234958
Editor: Paulo Lee Ho, Instituto Butantan, BRAZIL
Received: November 3, 2019; Accepted: June 3, 2020; Published: June 23, 2020
This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: Yes Saeed Aminzadeh, National Institute of Genetic Engineering and Biotechnology (NIGEB) grant number 971215-I-712
Competing interests: The authors have declared that no competing interests exist.
Introduction
Proteases represent a well-defined class of hydrolytic enzymes that are broadly used in numerous manufacturing processes such as detergents, leather, food, feed, and waste treatment industries. They represent one of the most important classes of commercial enzymes that account for almost a quarter of the total enzyme manufacturing worldwide [1]. Proteases act as “sharp scissors” that catalyse specific or non-specific digestive processes to produce new protein products, small peptides, or even individual amino acids. These enzymes play a vital role in a variety of biological processes, including cell growth, death, and immune responses [2].
Protease enzymes can be isolated from eukaryotic and prokaryotic sources. However, microorganisms represent the main source for the production of proteases due to their various economic and technical advantages. The culture of protease-producing microorganisms in high quantities and in a short time can be done using established fermentation methods. Furthermore, microorganisms facilitate the production of proteases with broad natural biochemical diversity and engineering of novel enzymes through genetic manipulations [3,4].
A major impediment to the commercial use of enzymes is the thermal stability constraint required for many industrial practices such as in the dairy industry and laundry detergent manufacturing [5–7]. Although several heat-stable proteases have been isolated from the thermophilus bacteria or produced through genetic engineering [8–10], there is still a strong need to identify and characterize new thermostable proteases from different microorganisms for environmental and industrial applications.
The goal of the present study was the recombinant expression of a previously isolated thermostable protease gene (Coh01144) from the thermophilic indigenous Cohnella sp. A01 [11]. The recombinant protease was purified and characterized in terms of biophysical and biochemical properties such as molecular docking and molecular dynamics simulation, thermal stability, pH profiling, kinetic parameters, and thermodynamic analysis to identify the potential industrial application of this protease.
Materials and methods
Strains, plasmids, reagents, and chemicals
E. coli strains DH5α and BL21 (DE3), and pET26b (+) vector were purchased from Invitrogen (USA). Restriction enzymes were provided by Fermentas (USA). PCR product extraction and high pure plasmid purification kits were purchased from Bioneer (Korea). High pure DNA purification kit was purchased from Roche (Germany). T4 DNA Ligase, isopropyl β-D-1-thiogalactopyranoside, and PVDF membranes were purchased from Thermo Fisher Scientific (Germany). Casein, azo-casein, 1 anilino naphthalene-8-sulfonate (ANS), and anti-polyHistidine peroxidase antibody were purchased from Sigma (USA). All other chemicals were obtained from Merck (Germany).
Structural studies
In silico structural and phylogenetic analysis.
The nucleotide sequence of Coh01147, previously obtained from Cohnella sp. A01 (Genbank Accession No: MK645908) was translated to the amino acids using Expasy translate tool (http://web.expasy.org/translate). The resulting protein sequence was used for further in silico analysis of the Cohnella sp. A01 derived protease 1147 protein. The biochemical properties of the protein, including molecular weight, isoelectric point, instability index, aliphatic coefficient, and grand average of hydropathy (GRAVY), were calculated using Expasy ProtParam tool (http://web.expasy.org/protparam) [12]. The presence of putative signal peptide and the location of its cleavage site was predicted using SignalP 5.0 server (http://www.cbs.dtu.dk/services/SignalP). pFam (http://pfam.xfam.org/search) and MEROPS (http://MEROPS.sanger.ac.uk) databases were used for sequence alignment and identifying the peptidase family. The phylogenetic relationships of the protein sequence were analysed using ClustalW in MEGA 7 programs. Multiple alignments of the resulting protein family of protease 1147 (from MEROPS) was carried out using COBALT (https://www.ncbi.nlm.nih.gov/tools/cobalt/re_cobalt.cgi), and the conserved areas were determined using DNAMAN software (www.lynnon.com). BLOSUM62 substitution matrix and a gap penalty = 12 were utilised in the search [13].
Homology modeling and validation.
BLASTP (at an E-value cutoff of < 3.5E-106 with BLOSUM62 matrix) was used for searching homologous sequences of protease 1147 in the protein data bank (PDB). The query resembled the 3D crystal structure of an intracellular protease from Pyrococcus horikoshii (PH1704, PDB ID: 1G2I_A) at a 2-Å resolution that was used as a template to construct the protease 1147 model. MODELLER 9v7 [14] was used to generate the homology model of the query. To evaluate and validate the modeled 3D structure, a Ramachandran plot was constructed using PROCHECK server (http://services.mbi.ucla.edu/PROCHECK) at the SAVE server (http://services.mbi.ucla.edu/SAVES). The ProSA web server (https://prosa.services.came.sbg.ac.at/prosa.pHp) was employed to calculate the Z-Score and evaluate the consistency between the crystal structure of the template (PH1704) and modeled protease 1147 [15]. The Chimera v1.13.1 was used to calculate the optimal energy of predicted protease 1147 structure and select the best model with the least root mean square deviation (RMSD) value by superimposing the model with the template [16, 17]. The active site of the proposed model was predicted using 3DLigandSite (http://www.sbg.bio.ic.ac.uk/3dligandsite/) and evaluated using the Pymol program (http://pymol.sourceforge.net/). The similarity of the predicted binding sites of protease 1147 and PH1704 was compared using Chimera v1.13.1.
Protein-ligand docking studies.
The modeled protease 1147 was docked with two surfactants with a positive effect on its activity (tween 20 and SDS) using Molegro Virtual Docker V.6.0 (MVD). The predicted structure was minimized, and the potential cavity (also referred to as the active site) was identified using the built-in cavity algorithm of MVD. The chemical structures of ligands were obtained from the PubChem and stored as Mol2 files using the Chimera v1.13.1. Finally, the structure of the protein was docked against the ligands using the AuotDock Vina tool. Grid size for the protein-ligand complex was set to 12 × 10 × 14 points. For each of the protease-SDS and protease-tween 20 complexes, 10 test runs were performed. Docked complex with lower binding energy was selected for MD simulation. Docking computations were performed in triplicate.
Molecular dynamics simulation.
The prepared protein and protein-ligand complexes (protease-SDS and protease-tween complexes) were used in MD simulation, performed using GROMACS v 4.6.5 with CHARMM36 all-atom force field. The CGenFF server (https://cgenff.umaryland.edu/) provides topologies and parameters of ligands compatible with the CHARMM36 all-atoms force field. Protein and protein-ligand complexes were soaked in a cubic box of water molecules, and the charges on the protein were neutralized by the addition of Na+ and Cl- ions. The energy of the system was minimised using the steepest descent algorithm to eliminate bad contact and clashes. The NVT and NPT ensembles were used during the equilibration to achieve the desired temperature (373.15 K for protease and 300 K for protease-ligand complex) and pressure (1 bar) for 100 ps and restraint forces of 1000 kJ/mol. Finally, 40 ns MD run were performed in triplicate after releasing all restraints. All bonds were constrained by the LINCS algorithm [18]. dx.doi.org/10.17504/protocols.io.bgcajsse
Cell culture and DNA extraction
The Cohnella sp. A01 cells were primarily cultured in a nutrient broth medium (0.5% peptone, 0.3% yeast extract, 0.5% sodium chloride, soluble in distilled water, pH adjusted to 6.8) overnight at 60°C on a shaker incubator (180 rpm). Of the culture medium, 2 ml (3.5 × 108 cells/ml) was added to 50 ml of the fresh nutrient broth medium and incubated in the same condition for 5 more days. The harvested cells were centrifuged at 3500 × g, and 4°C for 30 min and their genomic DNA was isolated using a high pure DNA purification kit, according to the manufacturer's instructions.
Cloning and heterologous expression of protease 1147 gene
Protease 1147 gene was amplified by PCR using specific primers (forward 5'-TACATATGAAGAAAGTCGCTTTCCTGC-3' and reverse 5'-TACTCGAGGCTCAGCTTGTTCAGCGTTTC-3') carrying NdeI and XhoI restriction enzyme recognition sites. The reverse primer was designed without stop codon to incorporate the His-tag sequence at the C-terminal of the protease. The purified PCR product was digested with the aforementioned restriction enzymes and cloned into the expression vector pET26b(+) using the T4 ligase. The recombinant plasmid pET26b(+) was transformed into E. coli DH5α through the heat shock transformation method, and the verified recombinant plasmid pET26b(+) was transformed into E. coli BL21 (DE3). A positive BL21 bacterial colony was incubated overnight at 37°C in 5 ml LB medium containing kanamycin (30 mg/ml) and inoculated to the fresh medium with the same antibiotic concentration. After the OD600 reached 0.6, the expression of the recombinant protein was initiated by the addition of 1 mM IPTG and incubating overnight at 28°C (150 rpm).
Recombinant protein purification
Ni-NTA affinity chromatography.
After the cultivation period, cells were collected by centrifugation at 3500 × g for 30 min at 4°C and resuspended in 4 ml lysis buffer (50 mM NaH2PO4, 300 mM NaCl, 10 mM Imidazole and 0.05% Tween 20 at a pH of 7). The mixture was sonicated with 4 × 40 sec pulses followed by 20 sec rest between cycles at 4°C. The crude extracts were then centrifuged at 9000 × g for 30 min at 4°C, and the resulting supernatant was subjected to Ni-NTA purification according to a previously described protocol [8]. In brief, 4 ml of sample was resuspended in binding buffer and loaded onto the Ni-NTA resin pre-equilibrated in 6 ml equilibration buffer (50 mM NaH2PO4, 300 mM NaCl, 10 mM Imidazole and 0.05% Tween 20% at pH = 7). Then, the column was washed three times with a buffer (50 mM NaH2PO4, 300 mM NaCl, 500 mM Imidazole and 0.05% Tween 20 at pH = 7) and four times with an elution buffer (50 mM NaH2PO4, 300 mM NaCl, 1200 mM Imidazole and 0.05% Tween 20 at pH = 7). The purified elution fractions were dialysed overnight in a 50 mM potassium-phosphate buffer (pH = 7) at 4°C. All the collected fractions were visualised using 12% sodium dodecyl sulfate-polyacrylamide gel electrophoresis. Protein bands were stained by Coomassie brilliant blue R250. Soluble purified protein concentration was measured using the Bradford method as described previously [19].
Single-step purification by heat shock.
Heat treatment was used as the second method for the protease 1147 purification. For this purpose, the supernatant sample obtained in the previous step was heat-treated in a hot water bath at 90°C for 15 min. Insoluble material was separated by centrifugation at 13000 × g for 10 min, and the supernatant, containing the enzyme, was analyzed on a 12% SDS–polyacrylamide. dx.doi.org/10.17504/protocols.io.bgbsjsne
Western blotting
After separation on the 12% SDS-PAGE, the proteins were transferred onto a PVDF membrane for Western blotting. The transferred membrane was blocked at room temperature for 5 h with TBST (25 mM Tris–HCl, pH 7.4, 0.14 mM NaCl, and 0.05% Tween 20) containing 5% BSA. The membrane was washed three times with TBST and incubated for 3 h in a 1:2000 diluted monoclonal anti-polyHistidin peroxidase at 37°C. After three more washes with the TBST, the target protein was visualized by developing the blot for 30 min with 3.4 mM 4-chloro-1-naphthol and 0.04% (v/v) H2O2 as a substrate.
Protease activity assay and kinetic measurements
To investigate the substrate specificity for the protease 1147, the enzyme activity was measured using a variety of commercially available protease substrates, including casein, BSA, gelatin, and azo-casein (1% w/w). Measurements of protease activity were performed according to Gulmez et al. [8] and with slight modification. In brief, to prepare the substrate solution, 250 μl of the substrate was dissolved in a 50 mM potassium phosphate buffer. The solution was mixed with 50 μl of protease 1147 enzyme and incubated at 60°C for 30 min. Of a fresh and cold 10% Trichloroacetic acid (TCA) solution, 300 μl was used to stop the reaction. After 1 h of incubation on ice, the reaction mixture was centrifuged for 10 min at 10000 × g and 4°C. The absorbance of the supernatant was measured at 280 nm.
Kinetic parameters of the purified protease 1147 were characterized in terms of Michaelis/Menten kinetic parameters (Km, Vmax, kcat, and kcat/Km) using non-linear regression by GraphPad Prism V.8 [20]. The protease activity was determined in the presence of a different concentration of casein (ranging from 10 to 70 mM) as substrate and at the optimum temperature. Assays were performed with 0.7 mg/ml of the purified enzyme and in triplicate [21].
Zymography
Casein Zymography was performed on SDS-PAGE according to the method of Garcia-Carreno et al. [22]. After electrophoresis, the gel was shaken for 1 h in 50 mM Tris-HCl buffer containing 2.5% Triton X-100 at 4°C (pH = 7) to remove SDS. To remove the Triton X-100, the gel was washed three times with 100 mM Tris–HCl buffer (pH 7). The gel was incubated for 1 h at 60°C in potassium phosphate buffer 50 mM, contained casein 1%, and subsequently stained with Coomassie Brilliant Blue R-250. Clear zones on the blue background indicated the presence of protease activity.
Effects of temperature and pH on the protease activity and stability
The effect of temperature on the protease activity was evaluated by measuring the enzyme activity at different temperatures from 10 to 100°C and with 10°C intervals. The 1% (w/w) casein was used as the substrate. To explore the temperature stability of the protein, the protease activity was measured at 60, 70, 80, and 90°C for 2 h and with 10 min intervals. All assays were performed in triplicate.
To determine the pH profile of the enzyme, the protease activity was performed by casein substrate in a pH range of 3 to 11 using 50 mM acetate (pH 3.6–5.6), 50 mM phosphate (pH 5.8–8.0), and 50 mM glycine (pH 8.6–10.6) buffers. An optimum temperature (60°C) was used for the assays. To investigate the pH stability, the residual activity was measured at pH values of 5, 9, and 11 for 2.5 h at 60°C and with 10min intervals.
Effects of metal ions, protease inhibitors, detergents, organic solvents on proteolytic activity
The effect of metal ions on protease activity was determined by measuring the enzymatic activity of the purified enzyme in the presence of 1, 2, and 5 mM concentrations of various metal ions including Fe2+, Mn2+, Mg2+, Ca2+, Co2+, Zn2+, Na+, K+, Al3+, and Li+. The effect of organic solvents was determined using the 5 and 10% (v/v) concentrations of acetone, methanol, ethanol, isopropanol, isobutanol, glycerol, DMSO, n-hexane, and chloroform as organic solvents in the reaction of enzyme activity measurements. The effect of surfactants on protease activity was evaluated using 1 and 2% concentrations of four different surfactants, i.e., Tween 20, Tween 80, Triton 100X, and SDS. To evaluate the effect of protease inhibitors, the enzyme activity was assayed in the presence of 1 and 2% concentrations of iodoacetamide (IAA), guanidinium hydrochloride (GuHCl), dithiothreitol (DTT), phenylmethylsulfonyl fluoride (PMSF), E-64 [Trans-Epoxysuccinyl-L-leucylamido (4-guanidino) butane], Leupeptin, and ethylenediaminetetraacetic acid (EDTA). The activity of the enzyme solution containing no metal ion was set as 100%, and the residual protease activity was measured.
Thermodynamic study
As described by Papamichael et al. the enzyme activation energy, and irreversible thermal inactivation can be described as the basic relations associating the changes of Gibbs free energy (ΔG), enthalpy (ΔH), and entropy (ΔS) to the values of equilibrium constants through absolute temperature and/or temperature stability [23].
The changes in the enthalpy and entropy of activation energy (ΔH‡, ΔS‡) to form enzyme/substrate complex [ES] were calculated from the rearranged Eyring absolute rate equation [23]:
(1)
where, kB is the Boltzmann constant (i.e. R/N, which is 1.38 × 10−23 J/K), N is the Avogadro’s number (6.02 × 10−23 mol-1), T is the temperature in Kelvin, ℏ is the Planck constant (6.63 × 10−34), R is the gas constant (8.314 J/K mol), ΔS‡ is the change in the entropy, and ΔH‡ is the change in the enthalpy. The ΔH‡ and ΔS‡ values were calculated from the slope and intercept, respectively, of Ln[ka/T] vs. 1/T as follows:
(2)
(3)
The Gibbs free energy of activation energy (ΔG‡) of the protease was calculated from:
(4)
or from the well-known equation:
(5)
(6)
The energy of activation (Ea‡) was calculated using the Arrhenius equation:
(7)
such that
(8)
The energy (Ea‡) involved in the activation process was calculated from the slope of a linear plot of 1/T vs. Ln [ka].
Free energy of enzyme/substrate binding [ΔG‡ E-S] was obtained from:
(9)
Free energy for transition state [ΔG‡ E-T] formation was acquired from:
(10)
Thermodynamic parameters of irreversible inactivation of protease 1147 were calculated using data obtained from the assay of temperature stability, as follows:
(11)
The values of enthalpy, entropy, and Gibbs free energy of irreversible thermal inactivation (ΔH#, ΔG#, and ΔS#) were calculated by applying Eqs 2, 3, and 4 with some modification including that in Eq 2 Ea#(in) was displaced with Ea‡, and in Eq 3, kin was used instead of kcat
The inactivation rate (kin m-1) was calculated by the following first-order expression:
(12)
which can also be expressed as:
(13)
Where, t is the incubation time, [Act]0 is the initial enzyme activity (i.e. enzyme activity at time 0), and [Act]t is the enzyme activity at the time ‘t’. The kin is the inactivation rate constant, calculated from the plots of Ln ([Act]t/[Act]0) vs. t.
The half-life (t1/2) of the enzyme is defined as the time required for the enzyme to lose one half of its initial activity, and expressed as:
(14)
The energy of activation for the inactivation process (Ea#) was calculated using the Arrhenius equation:
(15)
such that
(16)
The energy (Ea#) involved in the deactivation process was calculated from the slope of a linear plot of 1/T vs. Ln [kin].
Decimal reduction time (D value) was defined as the time required for 90% reduction in the initial enzyme activity at a specific temperature and was calculated as:
(17)
All the experiments were performed in triplicate, and the mean values were presented.
Structural analysis by fluorescence spectroscopy
Intrinsic and extrinsic fluorescence measurements were performed using Cary Eclipse spectrofluorometer (Varian Ltd., England). The measurements were made using a quartz cuvette with a 10 mm path length and at 25°C. The excitation slit width was set at 5 nm. Changes in the intrinsic fluorescence spectra were analysed to study the effects of temperature, pH, and two surfactants (SDS and Tween 20) on the tertiary structure of protease 1147. The spectra were recorded in the range of 300–400 nm by exciting the samples at 280 nm. Purified protease with the final concentration of 400 μg/ml was dialysed against 50 mM potassium phosphate buffer and used for experiments. The protease samples were incubated for 120 min at 50–100°C, at a pH range of 4–9 at 4°C, and in the presence of 1% SDS and Tween 20 (1 mM) at 4°C to study effects of temperature, pH, and the surfactants.
To study conformational changes of the recombinant protease 1147 that result in changes in its surface polarity, extrinsic fluorescent using ANS fluorescence probe was applied. The samples were excited at 375 nm, and spectra were collected from 400 to 600 nm. Effects of temperature, pH, and surfactants were assayed in the above-mentioned conditions and the presence of ANS with a final concentration of 40 μg/ml. All the experiments were performed in triplicate.
UV-far CD analysis
Circular dichroism (CD) spectrum (190–260 nm) of 0.2 mg/ml purified 1147 was recorded in the Avive-215 spectrometer (model 215, USA (equipped with a Peltier thermostat. Proteins were dialysed against 20 mM potassium phosphate buffer (pH = 7). The effects of temperature, pH, and surfactant on the protease 1147 unfolding were analyzed in the same experimental conditions used for fluorescence spectroscopy. The results were expressed as ellipticity (θ) (mdeg.cm2/dmol) using [θ] = (θ × 100 MRW)/(cl). The data obtained from three replicates and the percentage of secondary structures were calculated using CDNN 2.1 software. The results were deduced from the buffer CD signal and smoothed [24].
Results
Structural, biochemical and phylogenetic analysis of protease 1147
The DNA sequence was translated into a protein with 169 aa in length. No signal peptide was predicted by SignalP, suggesting the intracellular production of the protease 1147 (S1 Fig). Sequence alignment and phylogenetic tree analysis revealed that the protease 1147 has the highest similarity to thermostable Protease I from Alicyclobacillus macrosporangiidus and clade clustered with other PfpI families of intracellular proteases (Fig 1) that are characterised by their thermal stability and conserved catalytic tetrad (Glu77, Cys103, His104, and Gly105). The conserved position of the catalytic tetrad is shown in Fig 2A. According to these results, protease 1147 is the first reported protease from Cohnella sp. A01 that is classified in the DJ-1/ThiJ/PfpI superfamily. Amino acids involved in the catalytic tetrad of protease 1147 are shown in Fig 2B.
The unrooted neighbour-joining (NJ) tree was constructed based on the alignment of the protein sequences with the high similarities. The protease 1147 is marked with a green circle. The amino acid sequence of cysteine protease human cathepsin B was used as an out-group and marked with a red circle.
(A) Protein sequence alignment of the protease 1147 and several closely related sequences. The over residues in the boxed regions represent the conserved catalytic tetrad. (B) The amino acids involved in the catalytic tetrad of protease 1147 (Glu77 –Cys103 –His104 –Gly105) are marked in red ribbon sticks. The accession number of the sequences are represented.
The biochemical properties of protease 1147 and several related proteases were calculated by the ProtParam tool (Table A in S1 File). The predicted instability index was below 40, confirming that the protease 1147 structure is stable. The predicted aliphatic index (86.09) ascertained the high thermal stability of the protease. The negative GRAVY score denotes the hydrophilic character of the protease.
Homology modeling, docking, and molecular dynamics
Homology models were based on the crystal structure of PH1704, which had the highest sequence similarity (47%) with the protease 1147 and lowest E-value (4E-53). The ProSA-web has shown a Z-score of -8.34 that falls in the range of scores commonly found in the case of the similar native protein (Fig 3A). The ProSA results also confirmed that most of the residues have negative energy (S2 Fig).
(A) ProSA Z-score analysis for homology modeled protease 1147. The plot showed that the black spot (predicted model) with the Z-score value of -6.34 is within the range of native conformation of all proteins chains in PDB determined by X-ray crystallography or NMR spectroscopy concerning their length. (B) Ramachandran's analysis depicts that 98% of residues fall within the favoured, allowed, and acceptable regions, while the outlier region represents only 2%. (C) The 3D superimposition between the native structure of Pyrococcus horikoshii chain A (beige) and predicted protease 1147 structure (blue). The RMSD is 0.51 Å. (D) Catalytic tetrad residues are shown and labeled.
The Ramachandran plot confirmed the good quality of the final model, indicating that 74%, 20%, and 4% of residues were in the favoured, allowed, and acceptable regions, respectively (Fig 3B). The superimposed structure of native Pyrococcus horikoshii chain A and predicted protease 1147 model is shown in Fig 3C, and also, the structural comparison of protease 1147 and PH1704 indicated that the main pocket-forming amino acids (Glu77 –Cys103 –His104 –Gly105 in protease 1147 and Cys101 –His102 Gly103 in PH1704) are quite comparable. The calculated RMSD value (0.51 Å) suggests that the developed model is comparatively robust and suitable for further analysis. Catalytic tetrad residues are shown and labeled in Fig 3D.
The predicted structure of 1147, together with the amino acids forming the active site of the enzyme (Glu77 –Cys103 –His104 –Gly105), is shown in Fig 4.
(A) Predicted 3D structure of protease 1147. (B) Amino acids belonging to the catalytic tetrad (Glu77 –Cys103 –His104 –Gly105), and (C) nucleophilic elbow, enzyme cofactor, and active site flap in protease 1147 active site.
Molecular docking between the active site of protease 1147 and Tween 20 ligand predicted a hydrogen bond interaction between the oxygen atom at position 6 of Tween 20 and oxygen atom at position 6 of Arg162 with a distance of 3.10 Å. It was shown that after 40 ns of MD simulation, the predicted hydrogen bond was maintained (Fig 5A and 5B) and furthermore, a new hydrogen bond with a distance of 2.86 Å was formed between Tween 20 and Thr164 in protease (Fig 5C and 5D). No hydrogen bond was identified between the SDS ligand and the active site of protease 1147. The only predicted hydrogen bond between SDS and the protease was between oxygen atom at position 4 and Lys19 in distance of 2.62 Å (Fig 6A and 6B). The predicted hydrogen bond was lost at the end of 40 ns MD simulation, and SDS was stripped off from the protein (Fig 6C and 6D).
(A) Docked protease 1147-Tween 20 complex (before MD simulation), and (B) its 2-D display indicated an interaction between Arg162 of protease 1147 with oxygen atom at position 6 of Tween 20. (C) Docked protease 1147-Tween 20 complex (after MD simulation), and (D) its 2-D display indicated two interactions between Arg162 and Thr164 of protease 1147 with Tween20.
(A) Docked protease 1147-SDS complex indicated (before MD simulation), and (B) its 2-D display indicated an interaction between Lys19 of protease 1147 with oxygen atom at position 4. (C) Docked protease 1147-SDS complex (after MD simulation), and (D) its 2-D display indicated loss of the interaction. The results of docking runs were completely identical.
The overall RMSD of the docked protease 1147-Tween 20 complex (0.23 nm) was in the lower value compared to the protease 1147 (0.39 nm) (Fig 7A). Root Mean Square Fluctuation (RMSF) plot indicated less structural fluctuation profile in the active site of protease 1147 in complex with Tween 20 compared to the protease 1147 alone (Fig 7B). The radius of gyration (Rg) decreased in the protease 1147-Tween 20 complex, which indicates higher protease stability when it is complexed with Tween 20 (Fig 7C). The MD simulation of protease 1147 at 100°C indicated protease instability after 20 ns and loss of structure after 40 ns.
(A) RMSD graph, (B) RMSF graph, and (C) Rg Protease 1147 in the water environment (black), protease 1147-Tween20 complex (green), protease 1147-SDS complex (blue) and protease 1147 at 100°C temperature. The approximate positions of active site amino acids are shown in the RMSF plot. The simulations were run in triplicate to obtain average values. The error bars represent standard deviations and are depicted for every 10 ps.
Heterologous expression of protease 1147 in E. coli
The gene encoding protease 1147 from Cohnella sp. A01 genomic DNA was amplified using PCR and inserted into the expression plasmid pET26b(+) under the control of the T7 promoter. The accuracy of recombinant plasmid was authenticated using PCR with specific primers, double digestion, and sequencing (S3 Fig). This recombinant plasmid was used to transform E.coli BL21. The modified E.coli strain was cultivated in shake flasks. The SDS-PAGE analysis revealed intracellular accumulation of protein with a molecular weight of approximately 18 kDa and in a soluble form (Fig 8A).
(A) SDS-PAGE analysis of recombinant protease 1147. Lane 1 and 4) crude extract of recombinant BL21, lane2) crude extract of non-recombinant BL21 (negative control), lane 3) Ni-NTA purified, and lane 5) heat shock purified recombinant protease 1147. (B) Native-PAGE zymogram analysis of 6) Ni-NTA affinity column and, 7) heat shock single-step method purified protease 1147, 8) crude extract of non-recombinant BL21 (negative control). (C) Western blot analysis of (9) Ni-NTA affinity column purified protease 1147, and (10) negative control.
Purification, zymography and western blot analysis of soluble expressed protease 1147
The recombinant C-terminal His-tagged protease 1147 was purified effectively using Ni-NTA column and single-step thermal shock with overall yields of around 88% and 73%, respectively. The SDS-PAGE results of purification are shown in Fig 8A, and protein activity and efficiency of purification are summarized in Table 1. Casein substrate zymography showed the protease activity at an apparent molecular mass of 18 kDa (Fig 8B). Western blot analysis by anti-polyHistidin peroxidase antibody revealed a similar mass for the recombinant protease 1147 (Fig 8C).
Substrate specificity and kinetic measurement
Substrate specificity was inferred by comparing known protease substrates. The protease activity was higher when casein was used as substrate (100%), followed by gelatin (90%), azo-casein (25%), and BSA (10%). The kinetic parameters of protease 1147 were determined for the casein substrate by drawing the Michaelis-Menten curve (S4 Fig). The values of Km, kcat, and kcat/Km was determined 13.72±0.6 mM, 3.143 × 10−3 (s-1), and 0.381 (M-1 S-1), respectively.
High pH and thermal- stability of protease 1147
The maximum protease activity was obtained at 60°C (2.657 U/mg), and therefore, the optimal temperature was selected as 60°C (Fig 9A). The thermal stability of the recombinant enzyme is shown in Fig 9B. More than 80% of the protease activity remained after heat treatment at 70°C for 60 min. However, a sharp decrease was observed after that. The optimum pH was determined to be 7 for the protease 1147. The enzyme showed excellent pH stability in basic, neutral, and acidic pH (Fig 9C and 9D).
Effect of temperature on the recombinant protease 1147 (A) activity (B) stability. Effect of pH on the recombinant protease 1147 (C) activity (D) stability. The optimum temperature and pH were 60°C and 7, respectively. (E) Arrhenius plots for Ea‡. (F) Arrhenius plots for Ea#. (G) Thermal inactivation Ln at 60, 70, 80, 90, and 100°C.
Determination of protease 1147 activation energy and thermal inactivation parameters
The Ea‡ for protease 1147 was calculated 35.04 kJ/mol by Arrhenius equation (Fig 9E). ΔG‡, ΔH‡ and ΔS‡ at optimum temperature were 74.36 and 32.68 kJ/mole, and 125 JmolK-1, respectively (Table 2). Lower values of ΔG‡ and Ea‡ at optimum temperature revealed a higher stabilization and efficient transition state of the enzyme/substrate complex (ES‡) at 60°C.
The obtained values for ΔG‡E-T and ΔG‡E-S (-2.66 and 0.19 kJ/mol, respectively) at 60°C indicated that low amount of energy is required for protease 1147 to form the transition complex, and the negative ΔG‡E-T indicates that the reaction occurs spontaneously.
The kin values obtained for protease 1147 indicated an increasing trend with increasing the temperature from 60 to 100°C (Table 3). Using an irreversible thermal inactivation plot (Fig 9G), it is deduced that the protease 1147 has the highest t1/2 at 60°C, and increasing the temprature led to lower stability and t1/2. Furthermore, the kin, and D value were increased and decreased with increasing the temperature, respectively (Table 3). According to the results, protease 1147 can maintain its thermal stability at high temperatures.
The parameters were measured after 120 min incubation at various temperatures, and the values obtained at 60°C compared with other studied thermostable proteases.
The Ea#, for protease 1147 inactivation was calculate 50.31 kJ/mol using Arrhenius Equation (Fig 9F). ΔH# and ΔS# for inactivation of protease 1147, at optimum temperature, were 47.54 kJ/mole-1 and 154.70 JmolK-1, respectively, and showed a decreasing trend with increasing the temperatures (Table 3). The data obtained for ΔS#, ΔH#, and ΔG# corroborate the previous evidence showing the thermostable features of the enzyme. The calculated parameters of irreversible activation of protease 1147 are compared with those from other studies in Table 3.
Effects of metal ions, inhibitors, detergents, organic solvent, and surfactants on proteolytic activity of protease 1147
None of the metal ions in the selected concentrations showed a significant effect on the measured protease activity (Table B in S1 File). The effects of other molecules, including inhibitors, surfactants, and organic solvents, are summarized in Table 4. In organic solvents, β-Mercaptoethanol in lower concentrations had the highest positive effect on the activity of protease 1147. Inhibitory effects of methanol, ethanol, and isopropanol on the protease activity were observed at 3% v/v, while lower concentrations of methanol showed significant incremental effects. Also, 1% and 2% v/v of glycerol and DMSO had substantial positive effects on protease activity, respectively.
Each value represents the mean of three independent replicates with a standard deviation.
Monitoring temperature, pH, and surfactant-induced protease 1147 unfolding by fluorescence spectroscopy
The aromatic amino acids in the protein sequence are intrinsic fluorophores, which absorb light at approximately 280 nm and emit UV light. The results of fluorescence analysis in Fig 10A and 10B show that the maximum fluorescence emission was at 316 nm after excitation at 280 nm. The protease 1147 maintained its overall structure in temperatures lower than 70°C and pH of 4–8. While the addition of 1% SDS seems to cause a minor reduction in the fluorescence intensity, a substantial change was recorded after the addition of 1% Tween 20, which is due to the unfolding of protease 1147 tertiary structure (Fig 10C).
(A) Intrinsic fluorescence spectroscopy of protease 1147 in the temperature range of 40–100°C exhibited a decrease in fluorescence intensity in temperature higher than 80°C and a blue shift in 100°C. (B) Intrinsic fluorescence spectroscopy of protease 1147 in the pH range of 1–12 showed fluorescence decline in pH higher than 9 and lower than 3. A blue and red shift was observed in acidic and basic pH, respectively. (C) Intrinsic fluorescence spectroscopy of protease 1147 in the presence of 1% Tween 20 and SDS, indicated a slight reduction of fluorescence intensity for SDS and a huge increase and red shift for Tween 20. (D) Extrinsic fluorescence spectroscopy of protease 1147 in the temperature range of 40–100°C exhibited a slight increase in fluorescence strength for temperature below 80°C and a significant increase for 90 and 100°C. (E) Extrinsic fluorescence spectroscopy of protease 1147 in the pH range of 1–12 showed huge fluctuation and blue shift of fluorescence strength in strongly acidic (pH 1–3) solution. (F) Extrinsic fluorescence spectroscopy of protease 1147 in the presence of 1% Tween 20 and SDS indicated a slight increase of fluorescence intensity for SDS and a huge increase for Tween 20.
ANS fluorescence was used to analyse the exposure of hydrophobic patches in the protein. Fig 10D and 10E shows ANS spectra of protease 1147 in the various temperatures and pH. In temperatures lower than 80°C and pH between 4–12, the protease is well folded, and hydrophobic patches are not exposed. A large increase in the fluorescence strength was recorded for ANS outside of the temperature mentioned above and pH ranges, which is characteristic of protein unfolding. A striking increase in the fluorescence intensity was observed after the addition of 1% Tween 20, while 1% SDS caused slight changes in the fluorescence intensity (Fig 10F).
Determination of secondary structure changes by far-UV CD
The far-UV CD was performed to investigate protein folding and changes in the secondary structure of recombinant protease 1147 as a function of temperature, pH, and in the presence of Tween 20 and SDS. The secondary structure of purified protease 1147, at 25°C and pH of 7, was considered as control, and the result showed the typical signature of folded proteins containing less α-helix and more β-sheet structures. The far-UV CD spectra were recorded at different temperatures, pH values, and in the presence of two chemical additives, i.e., Tween 20 and SDS (Fig 11). The secondary structures were calculated using software CDNN and listed in Table 5.
Far-UV CD spectra of protease 1147, (A) at different temperatures, (B) at different pH values, and (C) in the presence of 1% Tween 20 and SDS. These spectra indicated the changes in the secondary structure content of protease 1147. The spectrum of purified protease 1147 at temperature 25°C and pH of 7 was considered as control.
The highest degree of ellipticity change and disintegration for protease 1147 structure occurred at 100°C and a wavelength of 220 nm. Furthermore, at this temperature, very intense disorganisation at 208 nm was observed compared to the control spectrum. As indicated in Table 5, the β-turn and α-helix levels increased by 0.2 nm and 0.7 nm at 100°C, respectively, while the rate of β-sheet was dropped 0.9 nm. At very acidic and basic pH values, intense disorganisation occurred, and the protease 1147 began to lose its secondary structure. The secondary structure contents were not significantly altered in the presence of 1% SDS. While Tween 20 respectively caused an increase and a decrease in the α-helix and β-sheets values.
Discussion
Because of significant industrial demands, there is a great interest in functional screening for novel proteases with improved characteristics. Due to several unique advantages such as large diversity, rapid growth, and small space required for cultivation, microbes (fungi and bacteria) are preferred sources for naturally-occurring proteases [33]. Compared to fungal proteases, bacterial proteases have higher reaction rates and better heat tolerance [34]. Several new, desired protease-encoding bacterial genes have been expressed in new hosts, with the aim of overproduction, characterization, and engineering [35, 36].
In this study, we recombinantly expressed, purified, and described the detailed structural and biochemical properties of a newly isolated protease 1147 from the thermophilic indigenous Cohnella sp. A0.1. In silico analysis of the translated sequence of protease 1147 demonstrated that it is produced intracellularly with a high aliphatic index, which is a characteristic of thermophilic bacterial proteins [37]. Similar to other related proteases that were listed in Table A in S1 File, the calculated negative GRAVY value indicated that this protein is non-polar [38]. Phylogenetic analysis reveals the high amino acid similarity between the protease 1147 and other thermostable intracellular cysteine proteases from DJ-1/ThiJ/PfpI superfamily. This superfamily has grown to a large family that has a representative in nearly all organisms. Despite closely-related 3D structures in their core fold, this superfamily has diverse functions including acting as chaperons, neuron protection, and glyoxalase III [39]. Although many members of this superfamily have been classified, few of them have been biochemically characterized [40].
The enzymatic activity of the protease 1147 was provided by the identified conserved catalytic tetrad consisting of Cys103, His104, Glu105, and Gly77, in which the thiol side chain of cysteine functions as a nucleophile during the initial stages of catalysis.
BLASTP analysis of the amino acid sequence indicated that protease 1147 exhibited the highest structural homology with the PH1704 protease from hyperthermophilic archaean Pyrococcus horikoshii, which is a member of DJ-1/ThiJ/PfpI and was identified as both aminopeptidase and endopeptidase [41]. The crystal structure of PH1704 at a 2.0-Å resolution was used as a template for the construction of the protease 1147 model. The homology model of the query was generated using MODELLER 9v7, and the quality of the modeled 3D structure was validated by the Ramachandran plot, with only 2% of residues classified as outliers. The structure was further verified using ProSA, and the calculated Z-score (-6.34) revealed a compatible value with the template and a negative balance for the potential energy of the entire structure. The RMSD value between protease 1147 and template was 0.51 Å, indicating a high structural similarity and good model structure.
The modeled building indicated that the protease 1147 has an overall α/β sandwich tertiary structure, which is present in the known structure of DJ-1/ThiJ/PfpI superfamily members [39]. Furthermore, the motif of the nucleophilic elbow, containing the conserved Cyc103 located at the very sharp turn between a β-strand and an α-helix, has been identified in the modeled structure. This motif is known to allow easy access to the cysteine by substrate [42]. The predicted 3D location and the topology of the active site in protease 1147 in the modeled structure has a strong match with known template PH1704.
Among different additives that were tested to evaluate the effect, it was shown that the protease activity was significantly enhanced in the presence of two non-ionic surfactants (Tween 20 and 80). Another ionic surfactant, SDS, at 1% concentration did not affect the enzyme activity. Interaction mechanisms of protease 1147 were investigated in the presence or absence of Tween 20 and SDS by molecular docking and MD. The results indicated that Tween 20 contacts the Arg162 in the active site of protease 1147 by hydrogen bond in aqueous solutions. The interaction remained almost stable after 40 ns of simulation and strengthened with two new contacts between the NH1 of the same Arg162 and Thr164 with the Tween 20 molecule. The only hydrogen bond between the Lys19 and SDS was diminished after the simulation time. This result explains why SDS did not affect the enzymatic activity of protease 1147.
The results of RMSD and RMSF and the reduced Rg emphasise that Tween 20 positively affects protein stability when complexed with the protease 1147. The plots of RMSD and RMSF were similar for protease 1147 alone and in complex with SDS. Tween 20 is a non-ionic surfactant commonly used in the formulation of biotherapeutic products to stabilise proteins and prevent protein aggregation and surface adsorption [43]. However, the mechanism of the stabilizing effect of Tween 20 and other non-ionic surfactants has not been fully explained. The dual nature of surfactants and their amphiphilic structure produce strong affinity for interphases, and this is the characteristic that is the basis of mechanisms by which surfactant affects the physicochemical properties of proteins [44, 45]. It is reported that, besides altering interactions of the protein with surfaces, Tween 20 acts as a chemical chaperone and aids in protein folding via direct hydrophobic interaction with proteins [46]. Other proteases with increased activity in the presence of Tween 20 are listed in Table C in S1 File.
Because protease 1147 has been obtained from thermophilic bacteria and maintained its activity up to 90°C, the MD simulation was run at 100°C to show the conformational changes. Significant differences were observed when comparing the RMSD and RMSF for the simulation at 100°C. Greater RMSD and RMSF and increased Rg observed at 100°C indicated that at this temperature, the protein started to denature, and because of a large amount of heat in the system, many bonds begin to break apart.
As shown in Fig 8, the recombinant protease 1147 was detected in the soluble form at approximately 18 kDa position and verified by zymography and Western blot analysis. The majority of bacterial proteases have a molecular weight ranging between 15 and 45 kDa [47]. The produced protein was purified using His-tag affinity column chromatography and a non-chromatographic purification method named single-step thermal shock. The latter method led to a lower purification fold and yield. However, this fast, low-cost, and convenient process makes it appropriate for use in laboratories and large-scale purification of protease 1147. The purified protease 1147, using both methods, was homogenous as it gave a single protein band on SDS-PAGE.
As we showed in this study, protease 1147 can tolerate a high range of temperatures. Thus, for the thermal shock method, we incubate the enzyme-containing supernatant at 90°C for 10 min. This method has been used successfully for purification of thermostable proteins but in lower temperatures and for a longer time (30 min at 60–70°C) [48, 49].
We subsequently characterized the purified recombinant protease 1147 biochemically, investigating the influence of pH, temperature, and different chemical additives on the enzyme activity. Casein was used as the substrate as protease 1147 displayed the highest specificity for it. The enzyme displayed excellent tolerance to high temperatures and a broad range of pH, with the highest activity at 60°C and a pH of 7. At 90°C, the protease had almost 50% activity, but the enzyme activity significantly decreased at 100°C, likely due to the denaturation of the enzyme. Furthermore, MD simulations at 100°C confirmed this result. This high degree of thermal stability may be beneficial in many industrial processes, for example, to decrease the risk of contamination at high temperatures and the cost of external cooling [50]. Concerning the influence of pH, protease 1147 shows approximately 80% and 50% of relative activity in pH values of 5 and 9, respectively, with good stability after 3 h of incubation at 60°C. The activity and stability over such a broad pH range and temperatures make protease 1147 potentially useful for different industrial applications.
Studies on the thermodynamic parameters of the enzyme have provided a basic insights into the factors that determine the enzyme stability. As we have shown here, the protease 1147 was shown to has intrinsic structural stability and maintains its activity at high temperatures.
The activation energy (Ea) is the energy that must be provided for the reaction to occur. The reduction of Ea barrier increases the fraction of reactant molecules that achieve sufficient energy to form a product. On the other hand, the maximum values of ΔH and ΔS suggested the effectiveness of the transitional state and negative ΔG measures the spontaneity of the catalytic reaction [26]. In this study, we performed a thermodynamic analysis to determine Ea‡, ΔG‡, ΔH‡, ΔS‡, ΔG‡E-T, ΔG‡E-S, and Ki for protease 1147. The low values of Ea‡, ΔG‡, besides high values of ΔH‡, ΔS‡ at the optimum temperature indicated the efficient transition state of ES‡ at this temperature. These parameters were measured for thermostable proteases in other studies such as those obtained from Bacillus stearothermophilus [25], or Aspergillus fumigatus [32]. Compared to the mentioned studies, here, we found lower ΔG‡, and Ea‡, and the lower amounts of ΔG‡E-T compared to the ΔG‡E-S in this study suggested that the reaction can occur spontaneously.
The parameters of irreversible inactivation are usually expressed in terms of kin, t1/2, and D value. Several studies have shown that t1/2 and D value are high in heat-resistant enzymes [25, 26, 51]. Compared to other heat-resistant proteases, the acquired high t1/2 and D value and notable low kin at 60°C confirmed the thermostability nature of the protease 1147. Moreover, the significant lower D value at the higher temperatures indicated its higher susceptibility to thermal inactivation in temperature above 60°C.
The activation energy for the thermal inactivation of an enzyme is a key factor in understanding its thermostability. Enzymes undergo a first-order inactivation reaction, which is responsible for irreversible denaturation. The high Ea# for inactivation means more energy is needed for denaturation of the enzyme that indicated higher thermostability. The thermodynamic parameters of irreversible inactivation, including Ea#, ΔH#, ΔG#, and ΔS# were evaluated for protease 1147 to assess the Ea#. The high value of ΔS# reflects its conformational stability and resistance towards the denaturation process [52]. The values of ΔH# and ΔS# of the protease 1147 showed a decreasing trend in higher temperatures, indicating the changes in the conformation of this enzyme toward the partially unfolded transition state [25]. The unfolding of the enzyme structure is associated with an increase in disorder or entropy led to the lower ΔS#. Furthermore, high measured amounts of ΔG# and Ea#, at the optimum temperature, revealed that protease 1147 resist against the unfolding of its transition state, and the enzyme required high inactivation energy to get denatured. However, increasing the temperature conversely affect the rigidity and thermostability of protease 1147. Overly, the findings indicate that the transition states of protease 1147 enzyme at optimum temperature was more ordered. These results, may provide a piece of evidence that the higher catalytic efficiency (kcat/Km) of the protease 1147 is due to the stability of the transition state [25].
Even though the bioinformatics analysis revealed that protease 1147 is most similar to the thermostable protease PH1704, it showed a different pattern of activity at the pH, and temperature ranges studied. It was shown that PH1704 has the highest activity at approximately 80°C and a pH of 8.5 [41]. To compare its optimum pH, temperature, and other kinetic parameters with protease 1147, a summary of characterized bacterial proteases from different sources was provided in Table D in the S1 File.
Many studies have reported that metal ions impose a stabilizing effect on enzymes and are required to maintain the enzyme structure, enhance the enzyme activity, and protect it against thermal inactivation [53, 54]. Different metal ions have been reported to affect various proteases activity [54, 55]. However, a limitation in the use of cysteine proteases is that metal ions readily suppress their activity. Therefore, these proteases require mild reductants and chelating agents and thus are not economical. As such, serine proteases that do not have this limitation may be preferred for industrial applications [56]. The effect of nine kinds of metal ions was determined at optimum pH and temperature, and the result showed that unlike many other cysteine proteases, protease 1147 is almost stable in different concentrations of metal ions with just a slight decrease in the proteolytic activity. According to the predicted structure, Zn2+ is needed as a co-factor to incorporate in the ligand-binding site of protease 1147. However, 1, 2, and 5 mM concentrations of Zn2+ did not affect the protease activity.
The protease 1147 was found to be fairly stable in the presence of studied organic solvents, and its activity enhanced by β-mercaptoethanol and lower concentrations of methanol, glycerol, and DMSO. This result is quite remarkable, considering that the technological utility of an enzyme can be greatly expanded if an enzyme can perform its activity in the presence of organic solvents rather than just an entirely aqueous reaction media [57, 58].
The protease 1147 enzyme activity was completely inhibited by IAM, E-64, and leupeptin, which are the typical cysteine protease inhibitors. However, metalloprotease inhibitor EDTA and denaturant agent GuHCl used conventionally in protease inhibitor cocktails did not have a significant inhibitory effect against protease 1147. PMSF, the serine/cysteine protease inhibitor used in a 2% concentrated solution, inhibits the protease activity up to approximately 50%. These results validated that this protease belongs to the cysteine protease family.
Fluorescence spectroscopy and CD analysis of proteins give a complete picture of the overall structure and are generally used to monitor the conformational changes of proteins with changes in the solvent composition or the environment [59]. Intrinsic fluorescence emission due to the tyrosine and tryptophan residue provides a sensitive reporter to characterize the protein. These residues are very sensitive to the polarity of the environment [60]. The intrinsic fluorescence emission of protease 1147 at a pH of 7 and 60°C displayed an emission maximum at 316 nm, suggesting that at least some of 5 tyrosine residues in the protease 1147 were exposed to the solvent. Increasing the temperature up to 80°C with pH in the range of 3 to 9 caused a small decrease in the fluorescence intensity, whereas at 100°C and in higher and lower pH values a sharp decrease in the fluorescence occurred.
The observed decrease in emission wavelength indicated that these environmental changes cause tyrosine residues to become more buried in the hydrophobic core of the enzyme, leading to the distinct structural changes. Upon the addition of 1 mM Tween 20, the emission spectra showed a 2-fold increase in intensity. This increase is likely due to the exposure of tyrosine residues. The tyrosine residues are exposed to the solvent in the unfolded state of proteins [61]. A surfactant such as Tween 20 has an anti-aggregation effect and increases the free energy of unfolding, exposing the protein side chains to the surrounding [61]. The wavelength bands of protease 1147 in the presence of 1% SDS showed an approximately 10% decrease in intensity. It was reported that a low concentration of SDS can promote protein aggregation [62], which can cause tryptophan residues to remain inaccessible. However, our results suggested negligible perturbation of protease 1147 structure in the presence of SDS 1%. Conformational changes in protease 1147 have also been studied by binding of fluorescence probe ANS. The results closely resemble those of the intrinsic fluoresce. A steep increase in fluorescence intensity was observed at 100°C, very basic and acidic pH, and in the presence of Tween 20. These results strongly suggest that protease 1147 under these conditions attains a conformation where hydrophobic patches are exposed on the surface. The protein at pH 7 is in its native compact form and thus the hydrophobic regions are out of reach. At pH 1 and 2 the highest fluorescence emission was detected, and a slight blue shift (4 nm) at pH 3 may be due to the aggregation of protein in acidic phase. The decrease in fluorescence intensity at pH 3 may be due to the presence of the molten globule (MG) intermediate state in the protease 1147. Similarly, lipase-3646 [24] and cysteine protease ZCPG [63] were shown that exhibit the behavior of the MG state at pH 3 and pH 2, respectively. The overall results were in agreement with those obtained for temperature and pH stability and the MD simulation in the presence of selected surfactants.
The far-UV CD was used to generalize the secondary structure of protease 1147 [64]. The CD spectrum shows a positive CD band near 200 nm and a negative one at 220 nm, which was consistent with the reported CD spectrum of Der p1 cysteine protease from douse dust, and papain from the latex of Carica papaya [65, 66]. The percentage of secondary structure was calculated as 17% α-helix, 31.8% β-sheet, and 20.5% β-turn. As the temperature increased, the ellipticity value considerably reduced, suggesting the loss of regular secondary structures. When the temperature was further raised to more than 80°C, the positive peak at 200 nm almost vanished. Disruption of structures in very basic and acidic pH values was also detected. At a pH of 2, a marked increase in the negative band was observed. Tween 20 causes a considerable decrease in the negative band at 200 nm, and the use of SDS shows no sign of change in the secondary structure of protease 1147. Overall, the high structural stability of protease 1147 is likely due to its high content of β-sheet, making it less susceptible to unfolding [67].
Supporting information
S1 Fig. SignalP data.
SignalP software suggested no signal peptide for protease 1147.
https://doi.org/10.1371/journal.pone.0234958.s001
(PPTX)
S2 Fig. Residues energy data.
ProSA local model quality/residue-wise energy plot shows most of the residues have negative energy.
https://doi.org/10.1371/journal.pone.0234958.s002
(PPTX)
S3 Fig. Cloning analyses.
M, DNA marker; 1, empty pET26b(+) vector; 2, recombinant pET26b(+) containing protease 1147 sequence; 3, double digestion of recombinant pET26b(+) by XhoI and NdeI restriction enzyme; 4 negative control PCR; 5, PCR from recombinant pET26b(+) vector with specific primers for protease 1147 sequence.
https://doi.org/10.1371/journal.pone.0234958.s003
(PPTX)
S4 Fig. Michaelis-Menten data.
Michaelis-Menten plot of protease 1147 activity as a function of casein concentrations.
https://doi.org/10.1371/journal.pone.0234958.s004
(PPTX)
S1 File.
Table A) Biochemical properties of protease 1147 and closely related proteases. NA, number of amino acids; MW, molecular weight; pI, theoretical isoelectric point; GRAVY, grand average of hydropathy; AI, aliphatic index; II, instability index. Table B) Effects of various metal ions on protease 1147 activity. Table C) The effect of metal ions, chemical compounds, inhibitors and surfactants on the activity of some other proteases. Table D) Type, Temperature, pH, substrate specificity and kinetic characterization of some other proteases.
https://doi.org/10.1371/journal.pone.0234958.s005
(DOCX)
Acknowledgments
The authors would like to thank the National Institute of Genetic Engineering and Biotechnology (NIGEB) for providing the research facility by grant number 971215-I-712.
References
- 1. Gupta R, Beg QK, Lorenz P. Bacterial alkaline proteases: molecular approaches and industrial applications. Applied microbiology and biotechnology. 2002;59(1):15–32. pmid:12073127.
- 2. López-Otín C, Bond JS. Proteases: multifunctional enzymes in life and disease. The Journal of biological chemistry. 2008;283(45):30433–7. pmid:18650443.
- 3. Nigam PS. Microbial enzymes with special characteristics for biotechnological applications. Biomolecules. 2013;3(3):597–611. pmid:24970183.
- 4. dos Santos Aguilar JG, Sato HH. Microbial proteases: Production and application in obtaining protein hydrolysates. Food Research International. 2018;103:253–62. https://doi.org/10.1016/j.foodres.2017.10.044 pmid:29389613
- 5. Haki GD, Rakshit SK. Developments in industrially important thermostable enzymes: a review. Bioresource Technology. 2003;89(1):17–34. https://doi.org/10.1016/S0960-8524(03)00033-6 pmid:12676497
- 6.
Margrethe Oxenboll K, Pontoppidan K, Fru F. Use of a Protease in Poultry Feed Offers Promising Environmental Benefits2011. 842 p.
- 7.
Kumar D, Thakur N, Verma R, Bhalla T. Microbial Proteases and Application as Laundry Detergent Additive2008. 661–72 p.
- 8. Gulmez C, Atakisi O, Dalginli KY, Atakisi E. A novel detergent additive: Organic solvent- and thermo-alkaline-stable recombinant subtilisin. International Journal of Biological Macromolecules. 2018;108:436–43. https://doi.org/10.1016/j.ijbiomac.2017.11.133 pmid:29174357
- 9. Ma X, Liu Y, Li Q, Liu L, Yi L, Ma L, et al. Expression, purification and identification of a thermolysin-like protease, neutral protease I, from Aspergillus oryzae with the Pichia pastoris expression system. Protein expression and purification. 2016;128:52–9. pmid:27539551.
- 10. Koma D, Yamanaka H, Moriyoshi K, Ohmoto T, Sakai K. Overexpression and characterization of thermostable serine protease in Escherichia coli encoded by the ORF TTE0824 from Thermoanaerobacter tengcongensis. Extremophiles. 2007;11(6):769–79. pmid:17657405
- 11. Aliabadi N, Aminzadeh S, Karkhane AA, Haghbeen K. Thermostable chitinase from Cohnella sp. A01: isolation and product optimization. Brazilian Journal of Microbiology. 2016;47:931–40. pmid:27528085
- 12.
Gasteiger E, Hoogland C, Gattiker A, Duvaud Se, Wilkins MR, Appel RD, et al. Protein Identification and Analysis Tools on the ExPASy Server. In: Walker JM, editor. The Proteomics Protocols Handbook. Totowa, NJ: Humana Press; 2005. p. 571–607.
- 13. Simossis VA, Heringa J. PRALINE: a multiple sequence alignment toolbox that integrates homology-extended and secondary structure information. Nucleic Acids Research. 2005;33(suppl_2):W289–W94. pmid:15980472
- 14. Šali A, Blundell TL. Comparative Protein Modelling by Satisfaction of Spatial Restraints. Journal of Molecular Biology. 1993;234(3):779–815. https://doi.org/10.1006/jmbi.1993.1626 pmid:8254673
- 15. Kandasamy S, Duraisamy S, Chinnappan S, Balakrishnan S, Thangasamy S, Muthusamy G, et al. Molecular modeling and docking of protease from Bacillus sp. for the keratin degradation. Biocatalysis and Agricultural Biotechnology. 2018;13:95–104.
- 16. Maiti R, Van Domselaar GH, Zhang H, Wishart DS. SuperPose: a simple server for sophisticated structural superposition. Nucleic Acids Research. 2004;32(suppl_2):W590–W4. pmid:15215457
- 17. Dhanavade MJ, Jalkute CB, Barage SH, Sonawane KD. Homology modeling, molecular docking and MD simulation studies to investigate role of cysteine protease from Xanthomonas campestris in degradation of Aβ peptide. Computers in Biology and Medicine. 2013;43(12):2063–70. https://doi.org/10.1016/j.compbiomed.2013.09.021 pmid:24290922
- 18. Hess B, Bekker H, Berendsen HJC, Fraaije JGEM. LINCS: A linear constraint solver for molecular simulations. Journal of Computational Chemistry. 1997;18(12):1463–72.
- 19. Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Analytical Biochemistry. 1976;72(1):248–54. https://doi.org/10.1016/0003-2697(76)90527-3.
- 20. Foukis A, Stergiou P-Y, Theodorou LG, Papagianni M, Papamichael EM. Purification, kinetic characterization and properties of a novel thermo-tolerant extracellular protease from Kluyveromyces marxianus IFO 0288 with potential biotechnological interest. Bioresource Technology. 2012;123:214–20. https://doi.org/10.1016/j.biortech.2012.06.090 pmid:22940322
- 21. Omrane Benmrad M, Moujehed E, Ben Elhoul M, Zaraî Jaouadi N, Mechri S, Rekik H, et al. A novel organic solvent- and detergent-stable serine alkaline protease from Trametes cingulata strain CTM10101. International Journal of Biological Macromolecules. 2016;91:961–72. https://doi.org/10.1016/j.ijbiomac.2016.06.025 pmid:27296442
- 22. Garciacarreno FL, Dimes LE, Haard NF. Substrate-Gel Electrophoresis for Composition and Molecular Weight of Proteinases or Proteinaceous Proteinase Inhibitors. Analytical Biochemistry. 1993;214(1):65–9. https://doi.org/10.1006/abio.1993.1457 pmid:8250256
- 23.
Papamichael EM, Stamatis H, Stergiou P-Y, Foukis A, Gkini OA. Chapter 3—Enzyme Kinetics and Modeling of Enzymatic Systems. In: Singh RS, Singhania RR, Pandey A, Larroche C, editors. Advances in Enzyme Technology: Elsevier; 2019. p. 71–104.
- 24. Golaki BP, Aminzadeh S, Karkhane AA, Yakhchali B, Farrokh P, Jazii FR, et al. The Effect of pH on Globular State of Lipase-3646; an Appropriate Model for Molten Globule Investigations. The Protein Journal. 2015;34(4):267–74. pmid:26239273
- 25. Abdel-Naby MA, Ahmed SA, Wehaidy HR, El-Mahdy SA. Catalytic, kinetic and thermodynamic properties of stabilized Bacillus stearothermophilus alkaline protease. International Journal of Biological Macromolecules. 2017;96:265–71. https://doi.org/10.1016/j.ijbiomac.2016.11.094 pmid:27899299
- 26. Sharma AK, Kikani BA, Singh SP. Biochemical, thermodynamic and structural characteristics of a biotechnologically compatible alkaline protease from a haloalkaliphilic, Nocardiopsis dassonvillei OK-18. International Journal of Biological Macromolecules. 2020;153:680–96. https://doi.org/10.1016/j.ijbiomac.2020.03.006 pmid:32145232
- 27. Ahmed SA, Saleh SAA, Abdel-Hameed SAM, Fayad AM. Catalytic, kinetic and thermodynamic properties of free and immobilized caseinase on mica glass-ceramics. Heliyon. 2019;5(5). pmid:31193050
- 28. Gohel SD, Singh SP. Characteristics and thermodynamics of a thermostable protease from a salt-tolerant alkaliphilic actinomycete. International Journal of Biological Macromolecules. 2013;56:20–7. https://doi.org/10.1016/j.ijbiomac.2013.01.028 pmid:23391710
- 29. Wehaidy HR, Abdel-Naby MA, Shousha WG, Elmallah MIY, Shawky MM. Improving the catalytic, kinetic and thermodynamic properties of Bacillus subtilis KU710517 milk clotting enzyme via conjugation with polyethylene glycol. International Journal of Biological Macromolecules. 2018;111:296–301. https://doi.org/10.1016/j.ijbiomac.2017.12.125 pmid:29309864
- 30. Gohel SD, Singh SP. Thermodynamics of a Ca2+-dependent highly thermostable alkaline protease from a haloalkliphilic actinomycete. International Journal of Biological Macromolecules. 2015;72:421–9. https://doi.org/10.1016/j.ijbiomac.2014.08.008 pmid:25150113
- 31. Guleria S, Walia A, Chauhan A, Shirkot CK. Purification and characterization of detergent stable alkaline protease from Bacillus amyloliquefaciens SP1 isolated from apple rhizosphere. Journal of Basic Microbiology. 2016;56(2):138–52. pmid:26375163
- 32. Hernández-Martínez R, Gutiérrez-Sánchez G, Bergmann CW, Loera-Corral O, Rojo-Domínguez A, Huerta-Ochoa S, et al. Purification and characterization of a thermodynamic stable serine protease from Aspergillus fumigatus. Process Biochemistry. 2011;46(10):2001–6. https://doi.org/10.1016/j.procbio.2011.07.013.
- 33. Sharma KM, Kumar R, Panwar S, Kumar A. Microbial alkaline proteases: Optimization of production parameters and their properties. J Genet Eng Biotechnol. 2017;15(1):115–26. pmid:30647648; PubMed Central PMCID: PMC6296574.
- 34. Hamza TA. Bacterial Protease Enzyme: Safe and Good Alternative for Industrial and Commercial Use. International Journal of Chemical and Biomolecular Science. 2017;3(1):1–10.
- 35. Joshi S, Satyanarayana T. Characteristics and applications of a recombinant alkaline serine protease from a novel bacterium Bacillus lehensis. Bioresour Technol. 2013;131:76–85. pmid:23340105.
- 36. Toymentseva AA, Danilova V, Tihonova AO, Sharipova MR, Balaban NP. Purification of recombinant extracellular proteases from Bacillus pumilus for ß-amyloid peptide cleavage. Russian Journal of Bioorganic Chemistry. 2016;42(1):53–8.
- 37. Ikai A. Thermostability and Aliphatic Index of Globular Proteins. The Journal of Biochemistry. 1980;88(6):1895–8. pmid:7462208
- 38. Kyte J, Doolittle RF. A simple method for displaying the hydropathic character of a protein. Journal of Molecular Biology. 1982;157(1):105–32. https://doi.org/10.1016/0022-2836(82)90515-0 pmid:7108955
- 39. Kim HJ, Lee KY, Kwon AR, Lee BJ. Structural and functional studies of SAV0551 from Staphylococcus aureus as a chaperone and glyoxalase III. Biosci Rep. 2017;37(6). pmid:29046369; PubMed Central PMCID: PMC5691139.
- 40. Wei Y, Ringe D, Wilson MA, Ondrechen MJ. Identification of functional subclasses in the DJ-1 superfamily proteins. PLoS Comput Biol. 2007;3(1):e10. pmid:17257049.
- 41. Zhan D, Bai A, Yu L, Han W, Feng Y. Characterization of the PH1704 protease from Pyrococcus horikoshii OT3 and the critical functions of Tyr120. PloS one. 2014;9(9):e103902. pmid:25192005
- 42. Udatha DB, Madsen KM, Panagiotou G, Olsson L. Multiple nucleophilic elbows leading to multiple active sites in a single module esterase from Sorangium cellulosum. J Struct Biol. 2015;190(3):314–27. pmid:25907516.
- 43. Kerwin BA. Polysorbates 20 and 80 used in the formulation of protein biotherapeutics: structure and degradation pathways. Journal of pharmaceutical sciences. 2008;97(8):2924–35. pmid:17973307.
- 44. Ruiz-Peña M, Oropesa-Nuñez R, Pons T, Louro SRW, Pérez-Gramatges A. Physico-chemical studies of molecular interactions between non-ionic surfactants and bovine serum albumin. Colloids and Surfaces B: Biointerfaces. 2010;75(1):282–9. pmid:19782541
- 45. Garidel P, Hoffmann C, Blume A. A thermodynamic analysis of the binding interaction between polysorbate 20 and 80 with human serum albumins and immunoglobulins: a contribution to understand colloidal protein stabilisation. Biophysical chemistry. 2009;143(1–2):70–8. pmid:19427732
- 46. Randolph TW, Jones LS. Surfactant-protein interactions. Pharmaceutical biotechnology. 2002;13:159–75. pmid:11987751.
- 47. Jalkute C, Waghmare S, Nadaf N, Dhanavade M, Jadhav D, Pendhari S, et al. Purification and characterization of SDS stable protease from Bacillus safensis strain CK. Biocatalysis and agricultural biotechnology. 2017;10:91–5.
- 48. Fu Z, Hamid SBA, Razak CNA, Basri M, Salleh AB, Rahman RNZA. Secretory expression in Escherichia coli and single-step purification of a heat-stable alkaline protease. Protein expression and purification. 2003;28(1):63–8. https://doi.org/10.1016/S1046-5928(02)00637-X pmid:12651108
- 49. Oka M, Yang Y-S, Nagata S, Esaki N, Tanaka H, Soda K. Overproduction of thermostable leucine dehydrogenase of Bacillus stearothermophilus and its one‐step purification from recombinant cells of Escherichia coli. Biotechnology and applied biochemistry. 1989;11(3):307–11. pmid:2503013
- 50. Amid M, Manap A, Yazid M, Zohdi NK. Purification and characterization of alkaline-thermostable protease enzyme from pitaya (Hylocereus polyrhizus) waste: a potential low cost of the enzyme. BioMed research international. 2014;2014.
- 51. Kikani BA, Kourien S, Rathod U. Stability and Thermodynamic Attributes of Starch Hydrolyzing α-Amylase of Anoxybacillus rupiensis TS-4. Starch—Stärke. 2020;72(1–2):1900105.
- 52. D'Amico S, Marx JC, Gerday C, Feller G. Activity-stability relationships in extremophilic enzymes. J Biol Chem. 2003;278(10):7891–6. Epub 2003/01/04. pmid:12511577.
- 53. Chittoor JT, Balaji L, Jayaraman G. Optimization of parameters that affect the activity of the alkaline protease from halotolerant bacterium, Bacillus acquimaris VITP4, by the application of response surface methodology and evaluation of the storage stability of the enzyme. Iranian journal of biotechnology. 2016;14(1):23. pmid:28959314
- 54. Shaikh I, Dixit P, Shaikh T. Purification and characterization of alkaline soda-bleach stable protease from Bacillus sp. APP-07 isolated from Laundromat soil. Journal of Genetic Engineering and Biotechnology. 2018;16(2):273–9. pmid:30733735
- 55. Doddapaneni KK, Tatineni R, Vellanki RN, Rachcha S, Anabrolu N, Narakuti V, et al. Purification and characterization of a solvent and detergent-stable novel protease from Bacillus cereus. Microbiological Research. 2009;164(4):383–90. pmid:17616381
- 56. Tripathi P, Tomar R, Jagannadham MV. Purification and biochemical characterisation of a novel protease streblin. Food Chemistry. 2011;125(3):1005–12.
- 57. Klibanov AM. Improving enzymes by using them in organic solvents. nature. 2001;409(6817):241. pmid:11196652
- 58. Li Q, Yi L, Marek P, Iverson BL. Commercial proteases: present and future. FEBS letters. 2013;587(8):1155–63. pmid:23318711
- 59. Reed CJ, Bushnell S, Evilia C. Circular dichroism and fluorescence spectroscopy of cysteinyl-tRNA synthetase from Halobacterium salinarum ssp. NRC-1 demonstrates that group I cations are particularly effective in providing structure and stability to this halophilic protein. PloS one. 2014;9(3):e89452. pmid:24594651
- 60. Tcherkasskaya O, Bychkova VE, Uversky VN, Gronenborn AM. Multisite Fluorescence in Proteins with Multiple Tryptophan Residues APOMYOGLOBIN NATURAL VARIANTS AND SITE-DIRECTED MUTANTS. Journal of Biological Chemistry. 2000;275(46):36285–94. pmid:10948189
- 61. Monsellier E, Bedouelle H. Quantitative measurement of protein stability from unfolding equilibria monitored with the fluorescence maximum wavelength. Protein Engineering Design and Selection. 2005;18(9):445–56.
- 62. Madsen JK, Christiansen G, Giehm L, Otzen DE. Release of Pharmaceutical Peptides in an Aggregated State: Using Fibrillar Polymorphism to Modulate Release Levels. Colloids and Interfaces. 2019;3(1):42.
- 63. Jamir K, Seshagirirao K. Fluorescence quenching, structural and unfolding studies of a purified cysteine protease, ZCPG from Zingiber montanum rhizome. International Journal of Biological Macromolecules. 2018;106:277–83. https://doi.org/10.1016/j.ijbiomac.2017.08.019 pmid:28802845
- 64. Corrêa DH, Ramos CH. The use of circular dichroism spectroscopy to study protein folding, form and function. African Journal of Biochemistry Research. 2009;3(5):164–73.
- 65. López-Rodríguez JC, Manosalva J, Cabrera-García JD, Escribese MM, Villalba M, Barber D, et al. Human glutathione-S-transferase pi potentiates the cysteine-protease activity of the Der p 1 allergen from house dust mite through a cysteine redox mechanism. Redox biology. 2019:101256. pmid:31229842
- 66. Qadeer A, Zaman M, Khan R. Inhibitory effect of post-micellar SDS concentration on thermal aggregation and activity of papain. Biochemistry (Moscow). 2014;79(8):785–96. pmid:25365488
- 67. Santucci R, Sinibaldi F, Fiorucci L. Protein folding, unfolding and misfolding: role played by intermediate States. Mini Rev Med Chem. 2008;8(1):57–62. Epub 2008/01/29. pmid:18220985.