Next Article in Journal
Ginsenoside Rg3 Prevents Oncogenic Long Noncoding RNA ATXN8OS from Inhibiting Tumor-Suppressive microRNA-424-5p in Breast Cancer Cells
Previous Article in Journal
Immune Monitoring upon Treatment with Biologics in Sjögren’s Syndrome: The What, Where, When, and How
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Biochemical Characterization and Functional Analysis of Heat Stable High Potential Protease of Bacillus amyloliquefaciens Strain HM48 from Soils of Dachigam National Park in Kashmir Himalaya

1
Department of Environmental Science, University of Kashmir, Srinagar, Jammu & Kashmir 190006, India
2
Division of Basic Sciences and Humanities, Faculty of Agriculture, SKUAST-Kashmir, Sopore 193201, Jammu & Kashmir, India
3
Centre of Research for Development, University of Kashmir, Srinagar, Jammu & Kashmir 190006, India
*
Authors to whom correspondence should be addressed.
Biomolecules 2021, 11(1), 117; https://doi.org/10.3390/biom11010117
Submission received: 25 October 2020 / Revised: 27 November 2020 / Accepted: 4 December 2020 / Published: 18 January 2021
(This article belongs to the Section Enzymology)

Abstract

:
A novel temperature stable alkaline protease yielding bacteria was isolated from the soils of Dachigam National Park, which is known to be inhabited by a wide variety of endemic plant and animal species of Western Himalaya. This high-potential protease producing isolate was characterized and identified as Bacillus amyloliquefaciens strain HM48 by morphological, Gram’s staining and biochemical techniques followed by molecular characterization using 16S rRNA approach. The extracellular protease of B. amyloliquefaciens HM48 was purified by precipitating with ammonium sulfate (80%), followed by dialysis and Gel filtration chromatography increasing its purity by 5.8-fold. The SDS–PAGE analysis of the purified enzyme confirmed a molecular weight of about ≈25 kDa. The enzyme displayed exceptional activity in a broad temperature range (10–90 °C) at pH 8.0, retaining its maximum at 70 °C, being the highest reported for this proteolytic Bacillus sp., with KM and Vmax of 11.71 mg/mL and 357.14 µmol/mL/min, respectively. The enzyme exhibited remarkable activity and stability against various metal ions, surfactants, oxidizing agent (H2O2), organic solvents and displayed outstanding compatibility with widely used detergents. This protease showed effective wash performance by exemplifying complete blood and egg-yolk stains removal at 70 °C and efficiently disintegrated chicken feathers making it of vital importance for laundry purpose and waste management. For functional analysis, protease gene amplification of strain HM48 yielded a nucleotide sequence of about 700 bp, which, when checked against the available sequences in NCBI, displayed similarity with subtilisin-like serine protease of B. amyloliquefaciens. The structure of this protease and its highest-priority substrate β-casein was generated through protein modeling. These protein models were validated through futuristic algorithms following which protein–protein (protease from HM48 and β-casein) docking was performed. The interaction profile of these proteins in the docked state with each other was also generated, shedding light on their finer details. Such attributes make this thermally stable protease novel and suitable for high-temperature industrial and environmental applications.

1. Introduction

Dachigam National Park (henceforth, DNP) situated in the Zabarwan range of Himalayas lies to the north-east part of Srinagar, the state capital of Jammu and Kashmir, India. The park is ecologically significant as it is known to harbor a unique variety of plants and animals, including the critically endangered Kashmir red stag, locally known as “Hangul”, endangered Himalayan black bear, Himalayan brown bear, Himalayan gray langur, musk deer and snow leopard. The national park has been extensively studied for its phytodiversity and soil physicochemical characteristics [1], as well as its endangered animals [2,3,4], but is poorly understood as far as its soil bacterial studies are concerned. The molecular level study on bacterial diversity in forest ecosystems of Kashmir revealed many novel bacterial lineages [5]. However, the molecular characterization of soil bacterial populations from lower DNP, particularly for their enzyme bioprospection, has not been investigated yet. Although the great Himalayan mountain range is expected to contain a wide variety of proteolytic bacteria that could remain active over a wide range of temperatures, limited studies mainly on the cold-active proteases from psychrophilic bacteria [6,7] have been conducted. Thus, there is a need to prospect heat-stable proteases for high-temperature operations where heat acts as a limiting factor for carrying out the reactions efficiently.
Protease enzymes being among the three major industrial enzyme groups, make up for around 65% of the gross enzyme trading globally [8] and are envisaged in the development of several eco-friendly bio-remedial technologies [9,10]. Although a wide variety of soil microbes produce alkaline proteases [11,12], one of the main drawbacks of these proteases is their unstable nature in alkaline pH at higher temperatures [13] as these conditions are deleterious for most of the enzymes. Thus, there is a need to prospect new heat-stable alkaline proteases having novel properties owing to growing enzyme demand with better stability [14].
Proteases also called proteinases and peptidases [15], are a class of enzymes catalyzing the total protein hydrolysis into peptides and amino acids by cleaving their peptide bonds [16]. Gupta et al. [17] recognized protease as the most frequently used biocatalysts commercially, that occupy a vital place on account of their utilization in numerous physiologic and mercantile operations [18] like pharmaceuticals, fine chemical production, bioethanol production, and detergent industries [19]. Microbial proteases in comparison to the plant and animal sources of proteases [17] are more advantageous by virtue of their fast growth, vast diversity, ease of cultivation and genetic manipulations [20]. Approximately two-thirds of the marketable proteases are obtained from several species of bacteria, fungi, and yeasts, of which bacteria-derived proteolytic enzymes, particularly from Bacillus species, are industrially significant, accounting for about 35% of microbial enzymes [10]. However, for successful industrial application, proteases must be stable and function at elevated temperatures as well as pH [21].
Of various protease types-acidic, neutral and alkaline [17], proteases belonging to the alkaline type being active and stable at elevated pH are the most frequently used commercial enzymes [22]. Alkaline proteases, EC.3.4.21-24;99, are such peptidases that mostly operate in a range of neutral to basic pH having either a metallo-type or serine-type center and are studied extensively due to their widespread utilization in detergent, food, leather and pharmaceutical applications [23]. Bacteria secreting alkaline proteases are widespread in natural [24,25,26] as well as in human-made environments [27]. However, genus Bacillus is one of the most indispensable group and has been extensively exploited for the alkaline protease enzyme production mainly on account of their easy isolation from varied habitats (alkalescent waters, deep sea, and soil), the potentiality to flourish on proteinaceous substrates [28], a rapid rate of growth, chemo-organotrophic growth characteristics, extracellular enzyme secretions, and safe handling [29].
Thus, the present study is a pioneer attempt intended to isolate and identify a potent protease enzyme-producing bacterial strain from the soils of lower DNP followed by extraction, purification, characterization and functional analysis of its extracellular protease through protein–protein docking to assess its properties for commercial exploitation and environmental waste management. To the best of our knowledge, there is no documentation regarding the production, characterization and substrate-binding protein–protein docking study of alkaline protease by this Bacillus sp. from this part of high-altitude Himalaya.

2. Materials and Methods

Study Site: DNP located in the North-Western Himalayan range is lying between 34°04′ N–34°11′ N and 74°54′ E–75°09′ E and stretches across an altitude of some 1677–4270 m (Figure 1). Officially the national park is sprawling over a roughly rectangular area of 141 km2 [30]. DNP has two ranges: lower (26 km2) and upper Dachigam (115 km2), preferably divided on the basis of the forest types, altitudinal range, and movements of Hangul deer. The present study was conducted in the lower region of DNP (also known as Lower Dachigam), comprising about one-third of its western-most portion, and consisting of a deep gorge cut across by the river Dagwan (and its tributaries) that originates in the alpine Marsar Lake (c. 4200 m) situated in Upper Dachigam (also known as the Dagwan Valley).

2.1. Sample Collection, Isolation and Preservation of Bacteria

Soil samples were brought seasonally for two years from selected sampling sites in lower DNP with a soil corer up to a depth of 15 cm [31], and standard serial dilution spread-plate technique was used for the isolation of culturable bacteria. The pure colonies of isolates were retrieved via subculturing and preserved into slants composed of nutrient-agar, which were then kept at 4 °C for future use [32].

2.2. Preliminary Screening of Isolated Bacterial Strains for the Protease Activity

All isolated bacterial strains were initially screened (in triplicates) for proteolytic activity onto the skim-milk agar plates [33]. Strain HM48, was selected for further study, as it depicted the highest hydrolytic zone of clearance around its colonies.

2.3. Qualitative Screening for Protease Secretion on Different Media

The qualitative enzyme screening of strain HM48 was carried out at pH 7.0 and 8.0 on different basal media (in triplicates) containing 1.0% (w/v) protein source of each casein, gelatin and skim milk [34]. After the incubation period, plates were analyzed for their hydrolytic zones around the bacterial colonies.

2.4. Identification of Strain, HM48

The strain, HM48 was identified using macromorphological [35], Gram staining, biochemical [36], and molecular approaches. This strain was also screened for the production of other hydrolytic enzymes (lipase and amylase) on respective media [37]. Antibiotic susceptibility tests using concentrations as per Clinical and Laboratory Standards Institute (CLSI) [38] were also performed for strain, HM48.
For molecular identification, the retrieved sequence of 16S ribosomal RNA from SciGenom, Labs, (Kerala) obtained after purification and sequencing of polymerase chain reaction (PCR) amplicon using universal bacterial primers [39] was blasted (BLASTn) determining its phylogenetic neighbors from the nucleotide database of National Center for Biotechnology Information, NCBI [40]. The evolutionary phylogenetic tree was made through the neighbor-joining method using maximum composite likelihood as a correction factor having a bootstrap value of 1000 replicates [41], and alignment was done by Clustal-W using MEGA 7.

2.5. Enzyme Production

2.5.1. Media for Protease Production and Culture Conditions

Inoculation of the prepared basal media (50 mL) was done using 5% of 24 h bacterial broth culture having OD550 ≈ 0.2 and incubated in an orbital shaker at the rate of 120 rpm for 24 h at 37 ± 2 °C [34].

2.5.2. Extraction of the Crude Proteases

After incubation, the inoculated basal media was centrifuged at 10,000 rpm for 15 min, and after discarding the cell pellet, the supernatant containing crude proteases was preserved at 4 °C and used for supplemental study [26].

2.5.3. Protein Estimation and Protease Activity Assay

The concentration of protein (mg/mL) was estimated following Lowry et al. (1951) and determined as a measure of its absorbance at 700 nm after each step of the enzyme purification process, using bovine serum albumin (BSA) as a reference [42]. Protease activity of the enzyme (U/mL) was assessed as per standard protocol by Cupp-Enyard, [43] with modifications, as pH and temperature of the assay were optimized wherein the enzyme was incubated at different pH (6–12) with temperature ranging 10–90 °C. Thus, at each purification process as well as for characterization of enzyme, the protease activity was measured at a standardized pH, 8 and temperature, 70 °C during this study.
The activity was assessed by incubating 5 mL of casein (0.65% w/v) prepared in 0.1 M Tris-HCl buffer (pH 8.0) with 1 mL of the enzyme at 70 °C for 30 min. This reaction was stopped by adding 5 mL of 10% (w/v) trichloroacetic acid (TCA) and incubated for 10 min at 37 °C followed by 15 min centrifugation at 10,000 rpm (4 °C). To 1 mL of supernatant collected in a fresh tube, 2.5 mL sodium carbonate (0.4 M) and 0.5 mL Folin–Ciocâlteu reagent (1 N) were added, gently mixed and incubated at 37 °C for 30 min. The absorbance was measured at 660 nm. For calculations, one-unit enzyme activity was interpreted as the amount of enzyme required to hydrolyze the casein to release 1 μmol tyrosine mL−1 min−1 under the standard conditions of the assay. A standard tyrosine curve was used to determine the number of tyrosine equivalents released. The specific activity of the sample was calculated by dividing the enzyme units (U mL−1) with the overall protein content (mg mL−1). Thus, the sample activity was denoted in units (U mL−1) while the specific activity was designated as enzyme activity per mg of protein and expressed as (U mg−1).

2.6. Enzyme Purification

2.6.1. Ammonium Sulfate Precipitation of Crude Proteases

The crude protease was precipitated from the supernatant using solid ammonium sulfate (40–80%) as per standard saturation levels [44]. The precipitated protein was resuspended in a small amount of 30 mM buffer (Tris-HCl, pH 7) to obtain a concentrated enzyme solution. This suspension was dialyzed with the same buffer for 12–18 h at 4 °C to remove ammonium sulfate [45]. The dialyzed sample fraction is termed as partially purified protease enzyme (PPPE).

2.6.2. Gel Filtration Column Chromatography and Gel Electrophoresis (SDS–PAGE)

The PPPE extract was then purified by passing through a Tris-buffer (pH 7)equilibrated Sephadex G-50 column (Sigma-Aldrich). The PPPE was eluted using the same buffer, and 3 mL fractions (about 25 in number) at 1 mL min−1 flow rate were collected, and their absorbance was read at 280 nm. Those fractions that exhibited maximum protease activity were pooled, concentrated, and kept for subsequent analysis.
Gel electrophoretic analysis was executed according to Laemmli [46], and the approximate molecular weight of the protease from strain HM48 was evaluated against a standard protein ladder (Spectra, Thermo Scientific).

2.7. Enzyme Characterization

2.7.1. Effect of Temperature and pH on Activity/Stability of Protease from Strain, HM48

The effect of temperature on enzyme activity was analyzed in a varied temperature range (10–90 °C). The thermal stability, on the other hand, was checked by pre-incubating the enzyme for 30 and 60 min at different temperatures (50, 55, 60, 65 and 70 °C) followed by measuring their percent residual activities (i.e., activity of sample after incubation/activity of sample before incubation ×100).
pH effects on the enzymatic activity were assessed using 0.1 M buffer systems comprising of phosphate buffer (pH 6–7), Tris-HCl buffer (pH 8–9) and glycine–NaOH buffer (pH 10–12), whereas the pH stability analysis was conducted by prior incubation of protease enzyme for 1 h at 35 ± 2 °C in pH buffers (6–12), and residual activities were analyzed.

2.7.2. Influence of Different Substrates, Metal Ions and EDTA on Enzyme Activity of Strain, HM48

The impact of various substrates on the activity of the enzyme was investigated employing substrates (1 mg/mL), e.g., BSA, casein, gelatin and skim milk. The reaction mixtures were kept for 30 min at 70 °C, after which activity was evaluated.
The effects of several metal ions (CdCl2, CuSO4, FeCl3, HgCl2, MgCl2, and MnCl2) on the enzyme catalytic behavior was examined by its prior incubation with a specified metal ion at 1, 5 and 10 mM. Furthermore, the metal ion with which the enzyme showed maximum increase in its activity at 10 mM final concentration was examined further with 1, 5 and 10 mM of ethylene-diamine tetra acetic acid, EDTA (a chelating agent) so as to detect its impact on enzyme catalytic efficiency.

2.7.3. Impact of Various Surfactants and Oxidizing Agent on Enzyme Efficiency of Strain, HM48

The effect of 0.5 and 1.0% concentrations of surfactants (sodium dodecyl sulfate, SDS in w/v; triton X-100 and tween-80 in v/v) and an oxidizing agent (hydrogen peroxide, H2O2 in v/v) on the enzyme activity was assessed, whereby reactions were kept for incubation at 35 ± 2 °C for 30 min using casein as a substrate and the residual activities were measured.

2.7.4. Effect of Different Organic Solvents on Activity/Stability of Protease from Strain, HM48

The influence of organic solvents upon catalytic activity of the enzyme at 70 °C was determined using an organic solvent, namely benzene, toluene and xylene. However, their influence on protease stability was measured by pre-incubating each of these solvents with the enzyme at 35 ± 2 °C for 30 min. After this, the residual activity of the enzyme was determined at 70 °C following the optimized assay.

2.7.5. Effect of Varied Casein Concentrations (Enzyme Kinetics) on Enzymatic Activity of Strain, HM48

Enzyme activity was analyzed using varying casein concentrations (0.12–2.25 mg/mL) at its optimum reaction conditions (pH, 8 and temperature, 70 °C), after which the enzyme kinetic parameter values (KM and Vmax) were computed employing a Lineweaver–Burk double reciprocal plot [47].

2.8. Protease Gene Amplification of Strain, HM48

Extracted template DNA of strain, HM48 was used for the amplification process using a primer set [34] synthesized by Sigma-Aldrich. The standardized cycling parameters (for 50 μL reaction) included 5 min initial denaturation at 94 °C followed by 35 PCR cycles, each comprising of 30 s denaturation at 94 °C, 30 s annealing at 52 °C, 60 s extension at 72 °C and 10 min final extension at 72 °C with a post-hold at 4 °C. The electrophoresed amplicon was sent to SciGenom Labs for purification and Sanger sequencing. Protein blast (BLASTp) was utilized for identifying the retrieved sequence after determining the phylogenetic relations from the nucleotide databases in NCBI [40]. The neighbor-joining approach was used to develop an evolutionary phylogenetic relationship [41], and alignment was done by Clustal-W using MEGA 7.

2.8.1. Model Generation

The structure from the amino acid sequence of strain, HM48 subtilisin-like serine protease, was generated using GalaxyTBM, of GalaxyWeb, in two stages. First, the extra reliable core structure was built through a template-based approach, following which unreliable local regions (URLs) were identified and subsequently remodeled through the thermodynamic hypothesis-based ab initio method [48]. The docking of the modeled structure was accomplished with the substrate, casein, with which this enzyme displayed maximum activity to carry out its active site analysis. For this purpose, a sequence of bovine beta (β) casein (Bos taurus) retrieved from UniProt (KB-P02666) was subjected to GalaxyTBM.

2.8.2. Improvement of Model Quality

The best model of strain, HM48 subtilisin-like serine protease, was further refined using a highly reliable refinement software, namely GalaxyRefine, which uses a sequential approach involving side-chain rebuilding and repacking. Besides the defined functions, GalaxyRefine performs structural relaxation through its molecular dynamic simulations [49]. The GalaxyRefine models show improvement both in local as well as in global structure quality, which is seen, especially when the models to be refined have been generated through modeling software’s based on futuristic or advanced algorithms [49,50].

2.8.3. Model Validation through Different Approaches

As a thumb rule, the models need to be certified prior to downstream studies to restrain errors due to defective modeling. The protease model of strain, HM48 subtilisin-like serine, was tested by a variety of methods. Its stereochemical quality was tested by PROCHECK, and the compatibility of the model (3D) to its primary structure (1D) was checked by Verify3D [51,52]. Further validation was performed by protein structure analysis (ProSA) web, an easy interface for the ProSA program [53]. Validation of bovine β-casein was also done through PROCHECK and ProSA-web.

2.8.4. Protein–Protein Docking and Interaction Profile Generation

A web-based protein–protein docking server HawkDock was used for docking β-casein to strain, HM48 subtilisin-like serine protease. Molecular mechanics/generalized Born surface area (MM/GBSA) was employed to re-rank the top 10 models because of its effective re-ranking ability [54,55,56,57]. From this docked complex, the interaction profile was generated using a multifunctional web server PDBsum [58].

2.9. Enzyme Application

2.9.1. Evaluation of Enzyme as a Detergent Additive: Enzyme Compatibility with Commercial Detergents

The detergent compatibility test with some commercially available detergents like Ghadi (RSPL Ltd., India), Surf Excel and Wheel (Hindustan Unilever Ltd., India) and Tide (Procter and Gamble Ltd., India) was carried out by incubating each of the detergent solutions (7 mg mL−1) with the enzyme in 1:2 ratio at the test temperature of 70 °C. The endogenic protease in the tested detergents was deactivated at 80 °C for about an hour prior to their use, after which the residual enzyme activity was estimated.

2.9.2. Wash Performance Analysis: Blood and Egg Yolk Stain Removal

The wash performance of the enzyme from HM48 for stain removal on white cotton cloth pieces (5 × 5 cm2) was studied using detergent solutions (7 mg mL−1). The blood and egg yolk stained cloth (SC) pieces were separately dried in a hot oven for about 5 min (95–100 °C), and two sets (in triplicates) of below treatments each for human blood and egg yolk stained cloth pieces were prepared:
(i)
20 mL Surf Excel solution + SC piece
(ii)
20 mL Tide solution + SC piece
(iii)
20 mL Surf Excel solution + 400 μL of enzyme + SC piece
(iv)
20 mL Tide solution + 400 μL of enzyme + SC piece
(v)
20 mL tap water (T.W.) + 400 μL of enzyme + SC piece
(vi)
20 mL tap water (T.W.) + each of blood and egg yolk SC piece as control, respectively
The treatment plates were incubated at 70 °C for 30 min and at regular intervals visually checked for the removal of stains. After the incubation period, the stained cloth pieces were thoroughly rinsed, air-dried and examined for blood and egg yolk stains. The untreated blood and egg yolk stained cloth was taken as control [26].

2.9.3. Application of Purified Enzyme in Waste Management: Degradation of Chicken Feather

The whole chicken feathers procured from a chicken farm were washed properly with tap water to remove blood, after which they were carefully rinsed with distilled water and autoclaved. The feathers were allowed to air-dry overnight, following which the disintegration test of the whole chicken feather was executed by incubating it with the purified enzyme at a standardized test temperature of 70 °C for about 4 h [59].

3. Results

3.1. Isolation and Preliminary Screening of Bacterial Strains for the Protease Activity

Out of 92 isolated morphologically different strains of bacteria, 39 bacteria showed the proteolytic activity upon preliminary screening, of which isolate HM48 formed a remarkable hydrolytic zone of clearance around its colonies and thus was chosen for further study (Figure 2).

3.2. Qualitative Screening and Identification of Strain, HM48

3.2.1. Qualitative Screening on Different Media

The qualitative screening of isolate HM48 at pH 7.0 and 8.0 on different media revealed that it formed prominent hydrolytic zones around its colonies (Figure 3) at pH 7.0 on media containing casein (1%) as a protein source (Table 1).

3.2.2. Morphological Identification

The macromorphological colony and Gram staining characteristics of strain HM48 are given in Table 2.

3.2.3. Biochemical Identification

The biochemical identification results revealed that isolate, HM48 was able to utilize monosaccharides such as D and L-arabinose, fructose, galactose, mannose, rhamnose, sorbose, and xylose except for dextrose. Moreover, it was capable of utilizing monosaccharide derivative, esculin, but malonate and salicin were not utilized by it. This strain was also found to be inefficient in utilizing several di- and tri-saccharides except for cellobiose, citrate, inositol, inulin, mannitol, sodium gluconate, sorbitol, sucrose, and trehalose. The antibiotic susceptibility tests revealed that this strain was extremely sensitive to all the tested antibiotics (Table 3). In addition, isolate HM48 was able to produce hydrolytic enzymes, lipase and amylase as depicted by clear zones around its colonies.

3.2.4. Molecular Identification

The unidirectional 16S rRNA gene analysis, a highly conserved gene (Figure 4A) which is widely used for prokaryotic species identification, retrieved a nucleotide sequence of 885 bp that was submitted to GenBank having accession number MN006180. The sequence was checked on the NCBI database using nucleotide BLASTn depicting 100% similarity with the Bacillus amyloliquefaciens species, and its evolutionary proximity with them is clearly visible in Figure 4B.

3.3. Purification of Enzyme and Molecular Weight Determination of Protease from Strain, HM48

The supernatant was utilized as a crude enzyme source, which was then precipitated using 40–80% standard ammonium sulfate saturation levels. However, 80% saturation turned out to be the best percentage for precipitating the crude proteases. This precipitated enzyme was subjected to subsequent purification via dialysis and gel filtration chromatography. The specific activity of the culture supernatant (10.54 U mg−1) was enhanced to 22.81 U mg−1 and 38.54 U mg−1 by ammonium sulfate precipitation and dialysis, respectively, after which the column chromatographic purification raised it to 61.05 U mg−1 resulting in a purification fold of 5.8 with a yield of about 36.6% (Table 4). As per the SDS–PAGE, the estimated molecular weight of enzyme from strain, HM48, was around ≈25 kDa as indicated by the manifestation of one band in the gel (Figure 5).

3.4. Enzyme Characterization

3.4.1. Effect of Temperature on Activity/Stability of Protease from Strain, HM48

The study on enzyme optimal temperature of strain, HM48 revealed that this enzyme was active over varied temperatures (10–90 °C), peaking at 70 °C, thereby supporting the novel heat-stable attribute of the enzyme (Figure 6A). A continuous increase was observed from 10 °C up to a temperature of 70 °C after-which any further increase in temperature decreased the activity substantially. However, at temperatures 10 °C and 90 °C, the enzyme exhibited 25.1% and 39.61% relative activity, respectively, with regard to 100% relative activity at an optimum temperature of 70 °C.
Thermal stability studies conducted in the temperature range of 50 to 70 °C (Figure 6B) reflected the stable nature of enzyme at low temperature (up to 50 °C), retaining 92.15% of its activity within 30 min and 84.25% post 60 min incubation at 50 °C. However, this enzyme swiftly lost most of its activity above 50 °C, suggesting the thermal inactivation of the enzyme by heat as the incubation time increases. Even at its optimum activity temperature of 70 °C, the enzyme maintained only 17.74 and 9.29% of its activity after 30 and 60 min, respectively.

3.4.2. Effect of pH on Activity/Stability of Protease Enzyme from Strain, HM48

Ideal pH of protease was assessed using various buffers, results of which revealed the activeness of the enzyme over a wide pH range (6–12), reaching its peak at pH 8.0 with relative activity of 100%. Enzyme activity was observed to elevate proportionally from pH 6.0 to pH 8.0, after which any further deviation in the pH value reduced the activity significantly. Moreover, with the change in pH towards the acidic end (pH 6.0), a progressive decrease in the relative activity (%) was observed. The protease enzyme from strain HM48 showed an activity of 81.8%, 92.5% and 73.1% at pH 7.0, 9.0, and 10, respectively, in relation to the activity at pH 8.0. However, with a shift in pH towards alkalinity, the enzyme showed an activity of 61.8% (pH 11) and 56.1% (pH 12). On the other hand, the pH-dependent stability profile of this protease revealed its stability from pH 8–11, preserving its maximal residual activity at pH 8.0 while as considering pH 6.0 and 12, only 22.20% and 35.08% of activity (as regards to the activity at pH 8.0) was measured (Figure 6C).

3.4.3. Effect of Different Substrates on Enzyme Activity of Strain, HM48

Of all tested substrates (BSA, casein, gelatin, and skim milk), the enzyme hydrolyzed almost all proteins making it a distinct feature of this novel enzyme from strain, HM48. However, the substrate with which the enzyme depicted its highest activity (92.87%) was casein (Figure 6D).

3.4.4. Influence of Various Metal Ions and EDTA on Enzyme Activity of Strain, HM48

The effects of these ions tested at different concentrations are given in Figure 7A–C. Among the metal ions, Mn2+ at all concentrations regulated the enzyme activity positively while remaining ions such as Cd2+, Cu2+, Fe2+, and Hg2+ regulated it negatively. In contrast, the addition of 1.0, 5.0, and 10 mM final concentration of Mn2+ ion to the enzyme enhanced its protease activity by 66.48, 87.85 and 181.33%, respectively. The Hg2+ ions, on the contrary, significantly reduced it from 48.17% at 1.0 mM to 7.99% at 5.0 mM, reaching its maximum inhibition (6.99%) at 10 mM concentrations. However, Mg2+ was found to be inhibitory at 1- and 10 mM concentrations, while it tended to increase the enzyme activity by 5.23% at 5 mM concentration. During the present study, the inhibiting effect of EDTA was noticed, as the percent enzyme activity was 91.21% at 1.0 mM EDTA concentration, which was subsequently reduced to 78.27 and 18.19% at 5.0 and 10 mM concentrations of EDTA, respectively, thereby suggesting the requirement of metal ion (like Mn2+) for optimal activeness of this protease (Figure 7D).

3.4.5. Effect of Various Surfactants and Oxidizing Agent on Enzyme Efficiency of Strain, HM48

This was mainly done to investigate whether the enzyme is compatible and stable in the presence of different surfactants and oxidizing agents for assessing its applications in industries, particularly for laundry purposes. During the present investigation, apart from SDS, the alkaline protease enzyme was quite stable with almost all the assayed surfactants. At both concentrations, SDS, an anionic surfactant, inhibited the activity of this protease, retaining about 43.76 and 95.31% activity at 0.5 and 1.0% concentration, respectively. Triton-X 100 (a nonionic detergent) enhanced the activity from 4.84% at 0.5% and 7.30 at 1.0% concentrations. Moreover, Tween-80, a nonionic surfactant, was found to enhance and then inhibit the enzyme activity from lower (0.5%) to higher concentration (1%). On the other hand, the oxidizing agent, H2O2, depicted an increase of 29.07% in enzyme activity at 0.5% concentration, which was drastically reduced to 8.53% at 1.0% concentration (Figure 8A,B).

3.4.6. Effect of Various Organic Solvents on Activity/Stability of Protease from Strain, HM48

The influence of different organic solvents exhibited activity was marginally affected in the presence of xylene (96.78%), while toluene (86.4%) and benzene (85.49%) had a moderate effect on it. Furthermore, the enzyme stability studies revealed that the enzyme showed remarkable stability with each of the tested organic solvents, retaining its maximum residual activity with xylene (87.53%) followed by toluene (79.94%) and benzene (71.54%), as shown in Figure 8C.

3.4.7. Effect of Varied Casein Concentrations (Enzyme Kinetics) on Activity of Protease from Strain, HM48

The enzyme activity ascertained at standardized optimum reaction conditions with varying concentrations of substrate (casein) depicted its maximum at 2 mg/mL casein concentration (Figure 8D). A swift increase in the reaction rate was noticed, particularly from 0.5 to 2.00 mg/mL, after which the rate was found to get stabilized. Moreover, the values of enzyme kinetic parameters for this alkaline proteolytic enzyme as determined by Lineweaver–Burk plot were computed to be 11.71 mg mL−1 (KM) and 357.14 µmol mL−1 min−1 (Vmax), respectively (Figure 8E).

3.5. Identification of Protease Gene

As shown in Figure 9A, amplification of this gene yielded a product of about 700 bp, which when checked against available sequences in NCBI via BLASTp displayed a similarity of 99.03% with subtilisin-like serine alkaline protease of Bacillus amyloliquefaciens. The retrieved nucleotide sequence was deposited under accession number MT333265 in GenBank, NCBI after-which this gene sequence was translated and modeled, followed by molecular docking with a specific substrate. The tree of phylogeny was constructed employing MEGA 7 (Figure 9B).

3.5.1. Model Generation, Improvement and Validation

The protein model of Bacillus amyloliquefaciens HM48 subtilisin-like serine protease was generated through GalaxyTBM, which was improved by GalaxyRefine (Figure 10A), and rendering was done using UCSF Chimera. Ramachandran plot of this modeled protein contained 92.9% residues in the most favored region, thus reflecting its stability (Figure 10B). The quality was further validated by Verify 3D results (Figure 10C) as 100% residues exhibited averaged score (3D-1D) ≥ 0.2 as well as ProSA-web with Z-score calculation of −7.34 (Figure 11A). The bovine β-casein model also passed the Ramachandran plot with 95.1% residues in the favorable region (Figure 11B) and ProSA-web analysis with −3.67 Z-score (Figure 11C).

3.5.2. Protein–Protein Docking

Among the top 10 docked complexes, the complex with the binding free energy of −48.43 (kcal/mol) seemed to be logical based on experimental evidence (Figure 12A). Though this evidence is sufficient to prove that β-casein binds to Bacillus amyloliquefaciens HM48 subtilisin-like serine protease at the active site, but it is still debatable whether the region of β-casein interacting with the active site of this extracellular serine protease is correct or not. Interaction profile study showed that 22 residues of Bacillus amyloliquefaciens HM48 subtilisin-like serine protease interact with 20 residues of β-casein. Nine hydrogen bonds were detected between the two proteins (Figure 12B,C).

3.6. Application of Enzyme from Strain, HM48

3.6.1. Evaluation of Enzyme as a Detergent Additive: Enzyme-Detergent Compatibility and Wash Performance Analysis

The detergent affinity tests of the enzyme (Figure 13A) suggested that this alkaline protease was most compatible with Tide, retaining 63.28% of activity even when exposed to a temperature of 70 °C for 30 min thus clearly reflecting its potential use in the production of commercial detergents. These results were confirmed by complete blood, and egg yolk stain removal (Figure 13B) with Tide conducted as a part of this study.

3.6.2. Evaluation of Protease Enzyme for Application in Waste Management: Chicken Feather Degradation

The feather disintegration study suggested that the alkaline protease enzyme from strain HM48 hydrolyzed a considerable amount of protein and thus was efficient in degrading the chicken feathers to a great extent (Figure 13C).

4. Discussion

In soil, the principal sources of proteases vary from microorganisms [60], plants [61], and animal excreta [62] to the decomposition of organic matter, and dry-wet deposition [63] and the chief microbial proteases (serine-alkaline proteases, subtilisin’s and subtilisin-like proteases) are mainly produced by Bacillus and Pseudomonas strains of bacteria [64]. Microbial proteases mostly have their optimum catalytic pH in the 3.0–12 range [61]; however, those isolated from soil microorganisms retain their maximum activities at optimum temperatures between 40 and 60 °C within pH 8–9 [65]. The present study was thus aimed at identifying a high pH and thermo-tolerant protease producing bacterial strain from the soils of lower Dachigam National Park, Kashmir, through various culture-based analysis and protease characterization assays to assess its bioprospecting potential. The strain, HM48, was identified as Bacillus sp. by colony characterization as per Bergey’s manual of systematic bacteriology [35] and the same was confirmed by 16S rRNA analysis as it showed the highest similarity (100% identity) with species databases of Bacillus amyloliquefaciens in NCBI. The strain, HM48 was found to be extremely sensitive to the tested antibiotics, including cephalosporin’s (cefoxitin, ceftazidime, cephalothin), penicillin’s (amoxicillin/clavulanic acid or amoxiclav, oxacillin and penicillin-G), vancomycin and teicoplanin that cause inhibition in the cell wall formation by preventing the transpeptidation reaction during peptidoglycan synthesis. It was also sensitive towards various common aminoglycoside (gentamicin) and macrolide antibiotics (clindamycin and erythromycin) that can cause the inhibition of peptide chain elongation during protein synthesis. Furthermore, this strain was found to be sensitive to fluoroquinolone antibiotics such as ofloxacin, leading to inhibition of nucleic acid synthesis by blocking the replication and transcription of DNA [66].
In enzyme studies, purification and characterization of the enzyme are of vital importance for a better understanding of their nature and efficient application in industry as well as waste processing practices. During this study, the protease enzyme was purified via several purification processes, enhancing its purification by 5.8-fold and increasing its specific activity to 61.05 Umg−1, as also reported by Ahmad et al. [8]. A study by Thebti et al. [67] suggested that ideal specific activity is usually attained at elevated ammonium sulfate precipitations that go up to 90% sometimes, and a similar finding was made in this study whereby 80% of the precipitation was found to yield the optimal specific activity [68]. Interestingly, after each successive purification, a decrease in the total activity, as well as the total protein content was observed, which is consistent with other works on proteases [69]. Likewise, enzyme recovery yield was also observed to decrease with every purification step on account of either the autolysis of the enzyme or removal of certain low specific activity proteases [70]. Although processes of protease purification may vary, in most of the studies, precipitating (ammonium sulfate) and dialyzing enzyme along with size exclusion column chromatography [71] are mostly used either separately or in combination so as to enhance the purity of the enzyme [68].
Based on SDS–PAGE, the molecular weight (MW) of the purified enzyme of B. amyloliquefaciens HM48 was around ≈25 kDa, which was undoubtedly in compliance with the study by Sai-ut et al. [72]. Many previously reported literatures also suggest low molecular weight proteases (<30 kDa) from the genus, Bacillus including B. megaterium (MW-28 kDa, [70]), B. cereus S8 (MW-21.8 kDa, [73]), and B. subtilis APMSU6 (MW-18.3 kDa, [74]).
Studies on the temperature profile of enzymes by B. amyloliquefaciens HM48 revealed the flexibility of the enzyme remaining active across varied temperatures (10–90 °C). The enzyme activity enhanced continuously from 10 °C attaining its maximum at 70 °C, above which it decreased drastically, and these results were much higher than previously reported 50 °C (10,72] and 60 °C [26,75] for the protease enzyme activity from several species of B. amyloliquefaciens. Enzymes, when exposed to higher temperatures for a prolonged time, may undergo conformational changes leading to their denaturation [76]. It is a well-known fact that the thermostability of an enzyme is established by various factors, such as distinctive composition and sequence of amino acids, ionic interactions between them, degree of hydrogen bonds and reduction in the hydrophobic surface area [77]. Various research works have extensively assessed the means of stabilizing enzymes whereby they corroborated the significance of these salt bridges in the thermal stability of proteins [78,79], and the same has been authenticated by the studies which suggest that mesophilic proteins have less number of salt bridges in comparison to their thermophilic counterparts [80]. Enzyme thermostability is also affected by several intrinsic factors, including molecular interactions concerning hydrophobic, hydrogen, ionic and metal binding, as well as extrinsic conditions, including environmental factors such as activators, cofactors, substrate, protein concentration, and also the presence of specific salts [81]. Thermal stability studies of B. amyloliquefaciens HM48 reflected that although the enzyme maintained 92.15% and 82.45% of activity post 30 and 60 min incubation (50 °C), respectively, it constantly decreased up to its optimum of 70 °C, which could be due to the self-hydrolyzing of amino acids in the protein chains. Therefore, it is evident that the optimum activity temperature of B. amyloliquefaciens HM48 was beyond many previously recorded Bacillus derived proteases [82,83], suggesting that this alkaline protease could be exploited for high-temperature (up to 70 °C) industrial applications.
Proteases, based on active pH optimum, are grouped into acidic, neutral, and basic (or alkaline) proteases [84]. Of all these, alkaline proteases are particularly important as they remain active in neutral to alkaline pH range and are widely exploited in cosmetic, detergent and leather industries [23]. However, their exploitation in these industries is mostly limited due to their inability to withstand elevated pH and temperature, along with detergent additives and several organic solvents [85]. Thus, exploration of novel bacteria producing highly active, heat-stable alkaline protease has gained much research attention [8,86]. During the present study, the enzyme from B. amyloliquefaciens HM48 was stable with pH buffers ranging from acidic to alkaline (pH 6–12), thereby reflecting that the enzyme has a wide pH spectrum reaching its maximal activity at pH 8.0, and an analogous observation has been recorded by Guleria et al. [75]. However, the enzyme under examination had a much wider alkaline pH range than many previously reported B. amyloliquefaciens species having pH in the range of 5.0–11 [27], 4.0–10 [10], 4.0–10 [72], 6.0–10 [83] and 6.0–11 [26]. On the other hand, the pH-dependent enzyme stability profile of B. amyloliquefaciens HM48 (pH 8–11) also did not display much variation in the residual activity of the enzyme, thereby suggesting stability of the enzyme in alkaline pH range. Therefore, proteases from strain HM48, like some other reported microbial proteases, were found to have profound activity mainly in the alkaline range [87] that could pave the way for its exploitation in industries.
Casein is a commonly utilized substrate in protease studies [88] as it is subjected to proteolysis with almost all proteolytic enzymes without any need of denaturing initially mainly due to its intricate arrangement and arbitrary structure [89]. During the study, substrate, casein showed the highest enzyme activity among different tested substrates, which are in accordance with findings of Mothe and Sultanpuram [83]. Several types of research also emphasized casein, being the ideal substrate [82] as milk proteins like casein, is a great source of amino acids for the growth of high cell masses [90]. Alkaline proteases require several metal ions to retain their stable nature at a higher temperature, which generally functions as a cofactor by acting as an ion or salt bridge between two amino acids, thereby preserving the exact conformation (active forms) of enzyme molecules [91]. Among the tested metal ions, only Mn2+ ion positively regulated the activity of the protease enzyme from strain HM48 at all concentrations [10], as manganese ions are reported to stabilize the molecular structure of enzymes in some Bacillus cereus derived proteases [92]. However, Hg2+ markedly reduced the enzyme activity from lower (1 mM) to highest (10 mM) concentrations, either by binding of Hg2+ to the critical side-chains of protein changing its active site structure or by causing general unfolding [93] and a similar trend has been advocated by Ibrahim et al. [94]. Various researchers have reported the positive as well as negative effects of metal ions like Cd2+, Cu2+, Fe2+, and Hg2+, Mg2+ and Mn2+ ions on enzyme activities from several Bacillus species [86,95]. EDTA, a chelating agent, is a water softener widely utilized as a detergent additive for removing various stains. A study on the inhibitory action of EDTA on the enzymatic activity of B. amyloliquefaciens HM48 concluded that it lost most of its activity when treated with different concentrations of EDTA, reaching its lowest with 10 mM EDTA concentration, thus suggesting that the enzyme could belong to serine alkaline proteases and our findings were in compliance with EL-Eskafy et al. [96].
High-quality detergent stable proteases must be congenial with most of the detergent additives, including bleach, oxidizing agent, surfactant, etc. [17]. Although anionic surfactant, SDS has been reported to have a negative influence on the protease enzyme activity [75,94], enzymes under the current study retained > 95% of the activity with SDS (1%). Both of the nonionic surfactants (triton-X100 and tween-80), however, were found to increase the activity from 4–12%, respectively. Mainly, proteases belonging to the genus Bacillus are generally unstable in the presence of oxidizing agents like H2O2 [97], yet the enzyme from B. amyloliquefaciens HM48 depicted an increase of 8–29% in its activity against H2O2. Thus, one of the promising aspects of this protease from B. amyloliquefaciens HM48 was its affinity with tested surfactants and oxidizing agents, as only a few wild microorganisms having compatibility, particularly with SDS and oxidizing agents, H2O2 are known [68,98]. Interestingly, the enzyme under investigation exhibited remarkable activity and stability, particularly towards xylene, retaining about > 96% relative activity and > 87% of residual activity, as organic solvents are known to stabilize the protease enzymes mainly by replacing few water molecules with organic ones resulting in stabilization of its conformation [99].
The enzyme kinetic attributes (KM and Vmax) of B. amyloliquefaciens HM48 were ascertained adopting Lineweaver–Burk double reciprocal plot. These values mostly rely on the enzyme source (microorganism) and the type of substrate, thereby indicating the sensitive nature of the enzyme against its substrate, which is believed to increase as the KM value increases while the Vmax value decreases [100]. Values of KM (11.71 mg mL−1) and Vmax (357.14 µmol mL−1 min−1) were computed with good correlation (R2 = 0.9986), which was higher than many other previously reported values of KM and Vmax for several bacterial species as 3.846 mg mL−1 and 76.923 U mL−1 min−1 [101]; 3.99 mg mL−1 and 41.49 U mL−1 [102].
As per MEROPS database on peptidases (http://merops.sanger.ac.uk/), members of subtilisin-like proteases (subtilases) broadly categorized as S8 sub-family of the serine proteases super-family, possess immensely preserved arrangement of amino acids in the active site whereby neighboring as well as oppositely charged pair residues interact with each other forming hydrogen bonds or salt bridges, that affect a variety of their structural and functional traits. The model of B. amyloliquefaciens HM48 subtilisin-like serine protease showed sound stereo-chemical quality as authenticated by the Ramachandran plot with 92.9% residues in the most favored region while only 6.0% residues occurred in the additional allowed region [103]. The high quality of this model was also certified by Verify 3D results in which 100% residues exhibited averaged score (3D-1D) ≥ 0.2, the minimum criteria being only 80% [104], thus indicating the compatibility between 3D structure and primary structure (amino-acid sequence). Evaluation of B. amyloliquefaciens HM48 subtilisin-like serine protease refined model through ProSA-web (Z-score calculation) further supported its correctness as the Z-score observed was well within the range of native proteins possessing similar size [105]. The bovine β-casein model also passed the stereochemical quality test successfully through the Ramachandran plot and ProSA-web, thus validating its accuracy [53].
The compatibility test with locally available commercial detergents suggested that the enzyme was most compatible with Tide detergent, retaining 63.28% of its activity, and the same was validated by the wash performance analysis of the enzyme. Furthermore, from the results of chicken feather disintegration, it could be inferred that this enzyme can be exploited for converting waste chicken feathers into protein hydrolysate, thereby helping in proper chicken feather waste management from local poultry houses. Moreover, this protein hydrolysate may contain beneficial proteins that can be used in the production of organic fertilizers and animal fodder [59,106].

5. Conclusions

This study reports isolation, selection, characterization and functional analysis of novel alkaline protease enzyme-producing B. amyloliquefaciens HM48 from soils of lower Dachigam National Park, Kashmir Himalaya. The purification and characterization studies of the extracted protease exhibited its active and stable nature over a wide range of pH and temperature with an optimum at 8.0 and 70 °C, respectively, even in the presence of several metal ions, inhibitors, surfactants, oxidizing agent and organic solvents. This stable temperature enzyme displayed compatibility with almost all of the tested commercial detergents after pre-incubating it with them at a high-temperature (70 °C). Such characteristics advocate it as a potential alternative, especially in industrial applications requiring elevated temperatures for desired results. The improved wash performance displayed by the complete removal of blood and egg stains with the enzyme-detergent formulation reflects that this enzyme may find its application as a detergent additive. On the other hand, a disintegration study of chicken feathers suggested successful exploitation of this novel heat-stable protease for environmentally sustainable and efficient management of waste chicken feathers, the protein hydrolysate of which could be used for largescale production of organic fertilizers and animal fodder. Lastly, the protein–protein modeling of B. amyloliquefaciens HM48 helped us in understanding its catalytic behavior whereby the results revealed its subtilisin-like serine nature. The protein–protein docking analysis of this enzyme with its ideal substrate, casein, reflected its active site, which could pave the way for a better understanding of its nature for different prospective industrial and environmental exploitations.

Author Contributions

Conceptualization, H.M., A.J. and B.A.G..; methodology, H.M., B.A.G., S.F. and S.A.G..; software, S.A.G.; formal analysis, H.M., B.A.G., A.J. and S.A.G.; investigation: H.M. and B.A.G.; resources, A.J., B.A.G. and R.N.; writing—original draft, H.M.; writing—review and editing: H.M., B.A.G., A.J., and S.A.G.; supervision: A.J. and B.A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research work received no external funding and was supported through a departmental research scholarship by the Department of Environmental Science, University of Kashmir (Grant No. F(DFPMS-SCH) KU/17). However, Dr. Shabir Ahmad Ganai sincerely thank DST-SERB for financial assistance in the form of a startup grant bearing File No: YSS/2015/001267.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

Our heartfelt gratitude is due to the Office of the Principal Chief Conservator of Forests, Chief Wildlife Warden, Wildlife Warden Central Division, Srinagar and Regional Wildlife Warden, Kashmir Region, Srinagar (Department of Wildlife Protection, Govt. of Jammu and Kashmir) for providing the required sampling permission (No. WLP/Res/2017-18/493-96) and necessary arrangements for safe sampling. We would like to cordially thank the technical and administrative staff of the Department of Environmental Science and Centre of Research for Development, the University of Kashmir, for their invaluable support and for providing us with the necessary infrastructure to carry out the research work. The authors would like to acknowledge Dr. Kataria ([email protected]), School of Food Science and Nutrition, the University of Leeds, for proofreading the manuscript.

Conflicts of Interest

The authors declare no conflict of interest and harmoniously approved this work for publication.

References

  1. Yaqoob, A.; Yunus, M.; Bhat, G.; Singh, D.P. Phytodiversity and Seasonal Variations in the Soil Characteristics of Shrublands of Dachigam National Park, Jammu and Kashmir, India. Clim. Chang. Environ. Sustain. 2015, 3, 137–143. [Google Scholar] [CrossRef]
  2. Thakur, M.; Sharma, L.K.; Charoo, S.A.; Sathyakumar, S. Conflict Bear Translocation: Investigating Population Genetics and Fate of Bear Translocation in Dachigam National Park, Jammu and Kashmir, India. PLoS ONE 2015, 10, e0132005. [Google Scholar] [CrossRef] [Green Version]
  3. Mir, Z.R.; Noor, A.; Habib, B.; Gopi, G.V. Seasonal population density and winter survival strategies of endangered Kashmir gray langur (Semnopithecus ajax) in Dachigam National Park, Kashmir, India. SpringerPlus 2015, 4, 562. [Google Scholar] [CrossRef] [Green Version]
  4. Ahmad, K.; Qureshi, Q.; Agoramoorthy, G.; Nigam, P. Habitat use patterns and food habits of the Kashmir red deer or Hangul (Cervus elaphus hanglu) in Dachigam National Park, Kashmir, India. Ethol. Ecol. Evol. 2015, 28, 85–101. [Google Scholar] [CrossRef]
  5. Ahmad, N.; Johri, S.; Abdin, M.Z.; Qazi, G.N. Molecular characterization of bacterial population in the forest soil of Kashmir, India. World J. Microbiol. Biotechnol. 2008, 25, 107–113. [Google Scholar] [CrossRef]
  6. Furhan, J.; Awasthi, P.; Sharma, S. Biochemical characterization and homology modelling of cold-active alkophilic protease from Northwestern Himalayas and its application in detergent industry. Biocatal. Agric. Biotechnol. 2019, 17, 726–735. [Google Scholar] [CrossRef]
  7. Yadav, A.N.; Sachan, S.G.; Verma, P.; Kaushik, R.; Saxena, A.K. Cold active hydrolytic enzymes production by psychrotrophic Bacilli isolated from three sub-glacial lakes of NW Indian Himalayas. J. Basic Microbiol. 2016, 56, 294–307. [Google Scholar] [CrossRef]
  8. Ahmad, W.; Tayyab, M.; Aftab, M.N.; Hashmi, A.S.; Ahmad, M.D.; Firyal, S.; Wasim, M.; Awan, A.R. Optimization of Conditions for the Higher Level Production of Protease: Characterization of Protease from Geobacillus SBS-4S. Waste Biomass Valoriz. 2020, 11, 6613–6623. [Google Scholar] [CrossRef]
  9. Jisha, V.N.; Smitha, R.B.; Pradeep, S.; Sreedevi, S.; Unni, K.N.; Sajith, S.; Priji, P.; Josh, M.S.; Benjamin, S. Versatility of microbial proteases. Adv. Enzym. Res. 2013, 1, 39–51. [Google Scholar] [CrossRef]
  10. Wang, J.; Xu, A.; Wan, Y.; Li, Q. Purification and Characterization of a New Metallo-Neutral Protease for Beer Brewing from Bacillus amyloliquefaciens SYB-001. Appl. Biochem. Biotechnol. 2013, 170, 2021–2033. [Google Scholar] [CrossRef]
  11. Si, J.-B.; Jang, E.-J.; Charalampopoulos, D.; Wee, Y.-J. Purification and Characterization of Microbial Protease Produced Extracellularly from Bacillus subtilis FBL-1. Biotechnol. Bioprocess Eng. 2018, 23, 176–182. [Google Scholar] [CrossRef] [Green Version]
  12. Jadhav, H.P.; Sonawane, M.S.; Khairnar, M.H.; Sayyed, R.Z. Production of alkaline protease by rhizospheric Bacillus cereus HP_RZ17 and Paenibacillus xylanilyticus HP_RZ19. Environ. Sustain. 2020, 3, 5–13. [Google Scholar] [CrossRef]
  13. Griffin, H.L.; Greene, R.V.; Cotta, M. Isolation and characterization of an alkaline protease from the marine shipworm bacterium. Curr. Microbiol. 1992, 24, 111–117. [Google Scholar] [CrossRef]
  14. Iyer, P.V.; Ananthanarayan, L. Enzyme stability and stabilization—Aqueous and non-aqueous environment. Process. Biochem. 2008, 43, 1019–1032. [Google Scholar] [CrossRef]
  15. Qureshi, A.S.; Bhutto, M.A.; Khushk, I.; Dahot, M.U. Optimization of cultural conditions for protease production by Bacillus subtilis EFRL 01. Afr. J. Biotechnol. 2011, 10, 5173–5181. [Google Scholar] [CrossRef]
  16. Sumantha, A.; Larroche, C.; Pandey, A. Microbiology and industrial biotechnology of food-grade proteases: A perspective. Food Technol. Biotechnol. 2006, 44, 211. [Google Scholar]
  17. Gupta, R.; Beg, Q.; Lorenz, P. Bacterial alkaline proteases: Molecular approaches and industrial applications. Appl. Microbiol. Biotechnol. 2002, 59, 15–32. [Google Scholar] [CrossRef]
  18. Annamalai, N.; Rajeswari, M.V.; Balasubramanian, T. Extraction, purification and application of thermostable and halostable alkaline protease from Bacillus alveayuensis CAS 5 using marine wastes. Food Bioprod. Process. 2014, 92, 335–342. [Google Scholar] [CrossRef]
  19. Enshasy, H.E.; Enein, A.; Helmy, S.; Azaly, Y.E. Optimization of the industrial production of alkaline protease by Bacillus licheniformis in different production scales. Aust. J. Basic Appl. Sci. 2008, 2, 583–593. [Google Scholar]
  20. Wilson, P. Production and characterisation of protease enzyme produced by a novel moderate thermophilic bacterium (EP1001) isolated from an alkaline hot spring, Zimbabwe. Afr. J. Microbiol. Res. 2012, 6, 5542–5551. [Google Scholar] [CrossRef]
  21. Manachini, P.L.; Fortina, M.G. Production in sea-water of thermostable alkaline proteases by a halotolerant strain of Bacillus licheniformis. Biotechnol. Lett. 1998, 20, 565–568. [Google Scholar] [CrossRef]
  22. Denizci, A.A.; Kazan, D.; Abeln, E.; Erarslan, A. Newly isolated Bacillus clausii GMBAE 42: An alkaline protease producer capable to grow under higly alkaline conditions. J. Appl. Microbiol. 2004, 96, 320–327. [Google Scholar] [CrossRef] [Green Version]
  23. Haddar, A.; Agrebi, R.; Bougatef, A.; Hmidet, N.; Sellami-Kamoun, A.; Nasri, M. Two detergent stable alkaline serine-proteases from Bacillus mojavensis A21: Purification, characterization and potential application as a laundry detergent additive. Bioresour. Technol. 2009, 100, 3366–3373. [Google Scholar] [CrossRef]
  24. Chen, T.-L.; Chou, Y.-J.; Chen, W.-M.; Arun, B.; Young, C.-C. Tepidimonas taiwanensis sp. nov., a novel alkaline-protease-producing bacterium isolated from a hot spring. Extremophiles 2005, 10, 35–40. [Google Scholar] [CrossRef]
  25. El Hidri, D.; Guesmi, A.; Najjari, A.; Cherif, H.; Ettoumi, B.; Hamdi, C.; Boudabous, A.; Cherif, A. Cultivation-Dependant Assessment, Diversity, and Ecology of Haloalkaliphilic Bacteria in Arid Saline Systems of Southern Tunisia. BioMed Res. Int. 2013, 2013, 1–15. [Google Scholar] [CrossRef] [Green Version]
  26. Nassar, F.R.; Abdelhafez, A.A.; El-Tayeb, T.S.; Abu-Hussein, S.H. Purification, Characterization and Applications of Proteases Produced by Bacillus amyloliquefaciens 35s Isolated from Soil of the Nile Delta of Egypt. Br. Microbiol. Res. J. 2015, 6, 286–302. [Google Scholar] [CrossRef]
  27. Cho, S.-J.; Oh, S.-H.; Pridmore, R.D.; Juillerat, M.A.; Lee, C.-H. Purification and Characterization of Proteases fromBacillus amyloliquefaciensIsolated from Traditional Soybean Fermentation Starter. J. Agric. Food Chem. 2003, 51, 7664–7670. [Google Scholar] [CrossRef]
  28. Longo, M.A.; Novella, I.S.; Garcia, L.A.; Diaz, M. Comparison of Bacillus subtilis and Serratia marcescens as protease producers under different operating conditions. J. Biosci. Bioeng. 1999, 88, 35–40. [Google Scholar] [CrossRef]
  29. Cui, H.; Yang, M.; Wang, L.; Xian, C.J. Identification of a New Marine Bacterial Strain SD8 and Optimization of Its Culture Conditions for Producing Alkaline Protease. PLoS ONE 2015, 10, e0146067. [Google Scholar] [CrossRef]
  30. Jammu and Kashmir Forest Department, JKFD. Handbook of Forest Statistics; Jammu and Kashmir Forest Department (JKFD): Srinagar, Jammu and Kashmir, India, 2011.
  31. Anderson, J.M.; Ingram, J.S.I. Tropical Soil Biology and Fertility: A Handbook of Methods; CAB International: Wallingford, WA, USA, 1993. [Google Scholar]
  32. Kumar, R.; Acharya, C.; Joshi, S.R. Isolation and analyses of uranium tolerant Serratia marcescens strains and their utilization for aerobic uranium U(VI) bioadsorption. J. Microbiol. 2011, 49, 568–574. [Google Scholar] [CrossRef]
  33. Baghel, V.S.; Tripathi, R.D.; Ramteke, P.W.; Gopal, K.; Dwivedi, S.; Jain, R.K.; Rai, U.N.; Singh, S.N. Psychrotrophic proteolytic bacteria from cold environment of Gangotri glacier, Western Himalaya, India. Enzym. Microb. Technol. 2005, 36, 654–659. [Google Scholar] [CrossRef]
  34. Baweja, M.; Tiwari, R.; Singh, P.K.; Nain, L.; Shukla, P. An Alkaline Protease from Bacillus pumilus MP 27: Functional Analysis of Its Binding Model toward Its Applications As Detergent Additive. Front. Microbiol. 2016, 7, 1195. [Google Scholar] [CrossRef] [Green Version]
  35. Bergey, D.H.; Holt, J.G. Bergey’s Manual of Determinative Bacteriology; Williams & Wilkins: Baltimore, MD, USA, 1994. [Google Scholar]
  36. Pammi, N.; Chaitanya, K.; Mahmood, S.K. Serratia marcescens OU50T sp. nov., a cellulose and pha producing bacterium isolated from polluted water. Int. J. Environ. Biol. 2015, 5, 32–36. [Google Scholar]
  37. Sivanandhini, T.; Subbaiya, R.; Gopinath, M.; Angrasan, J.K.V.; Kabilan, T.; Selvam, M.M. An investigation on morphological characterization of Actinomycetes isolated from marine sediments. Res. J. Pharm. Biol. Chem. Sci. 2015, 6, 1234–1243. [Google Scholar]
  38. Saggu, S.K.; Mishra, P.C. Characterization of thermostable alkaline proteases from Bacillus infantis SKS1 isolated from garden soil. PLoS ONE 2017, 12, e0188724. [Google Scholar] [CrossRef] [Green Version]
  39. Gürtler, V.; Stanisich, V.A. New approaches to typing and identification of bacteria using the 16S-23S rDNA spacer region. Microbiology 1996, 142, 3–16. [Google Scholar] [CrossRef] [Green Version]
  40. Johnson, M.; Zaretskaya, I.; Raytselis, Y.; Merezhuk, Y.; McGinnis, S.; Madden, T.L. NCBI BLAST: A better web interface. Nucleic Acids Res. 2008, 36, W5–W9. [Google Scholar] [CrossRef]
  41. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  42. Lowry, O.I.; Rosebrough, N.J.; Farr, A.L.; Randall, R.J. Protein measurement with folin phenol reagent. J. Biol. Chem. 1951, 193, 265–275. [Google Scholar]
  43. Cupp-Enyard, C. Sigma’s Non-specific Protease Activity Assay—Casein as a Substrate. J. Vis. Exp. 2008, 17, e899. [Google Scholar] [CrossRef]
  44. Green, A.A.; Hughes, W.L. [10] Protein fractionation on the basis of solubility in aqueous solutions of salts and organic solvents. Methods Enzymol. 1955, 1, 67–90. [Google Scholar] [CrossRef]
  45. Sadia, A.A.; Shelly, D.; Mohan, P. Purification and characterization of a novel protease from Bacillus strain Sal1. Afr. J. Biotechnol. 2009, 8, 3603–3609. [Google Scholar] [CrossRef]
  46. Laemmli, U.K. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature 1970, 227, 680–685. [Google Scholar] [CrossRef]
  47. Lineweaver, H.; Burk, D. The Determination of Enzyme Dissociation Constants. J. Am. Chem. Soc. 1934, 56, 658–666. [Google Scholar] [CrossRef]
  48. Ko, J.; Park, H.; Seok, C. GalaxyTBM: Template-based modeling by building a reliable core and refining unreliable local regions. BMC Bioinform. 2012, 13, 198. [Google Scholar] [CrossRef] [Green Version]
  49. Heo, L.; Park, H.; Seok, C. GalaxyRefine: Protein structure refinement driven by side-chain repacking. Nucleic Acids Res. 2013, 41, W384–W388. [Google Scholar] [CrossRef] [Green Version]
  50. Farooq, Z.; Abdullah, E.; Banday, S.; Ganai, S.A.; Rashid, R.; Mushtaq, A.; Rashid, S.; Altaf, M. Vigilin protein Vgl1 is required for heterochromatin-mediated gene silencing in Schizosaccharomyces pombe. J. Biol. Chem. 2019, 294, 18029–18040. [Google Scholar] [CrossRef]
  51. Lüthy, R.; Bowie, J.U.; Eisenberg, D. Assessment of protein models with three-dimensional profiles. Nat. Cell Biol. 1992, 356, 83–85. [Google Scholar] [CrossRef]
  52. Laskowski, R.A.; MacArthur, M.W.; Moss, D.S.; Thornton, J.M. PROCHECK: A program to check the stereochemical quality of protein structures. J. Appl. Crystallogr. 1993, 26, 283–291. [Google Scholar] [CrossRef]
  53. Wiederstein, M.; Sippl, M.J. ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 2007, 35 (Suppl. 2), W407–W410. [Google Scholar] [CrossRef] [Green Version]
  54. Hou, T.; Qiao, X.; Zhang, W.; Xu, X. Empirical Aqueous Solvation Models Based on Accessible Surface Areas with Implicit Electrostatics. J. Phys. Chem. B 2002, 106, 11295–11304. [Google Scholar] [CrossRef]
  55. Zacharias, M. Protein-protein docking with a reduced protein model accounting for side-chain flexibility. Protein Sci. 2003, 12, 1271–1282. [Google Scholar] [CrossRef]
  56. Feng, T.; Chen, F.; Kang, Y.; Sun, H.; Liu, H.; Li, D.; Zhu, F.; Hou, T.-J. HawkRank: A new scoring function for protein–protein docking based on weighted energy terms. J. Cheminform. 2017, 9, 1–15. [Google Scholar] [CrossRef] [Green Version]
  57. Weng, G.; Wang, E.; Wang, Z.; Liu, H.; Zhu, F.; Li, D.; Hou, T.-J. HawkDock: A web server to predict and analyze the protein–protein complex based on computational docking and MM/GBSA. Nucleic Acids Res. 2019, 47, W322–W330. [Google Scholar] [CrossRef]
  58. Laskowski, R.A.; Jabłońska, J.; Pravda, L.; Vařeková, R.S.; Thornton, J.M. PDBsum: Structural summaries of PDB entries. Protein Sci. 2018, 27, 129–134. [Google Scholar] [CrossRef]
  59. Bhunia, B.; Basak, B.; Mandal, T.; Bhattacharya, P.; Dey, A. Effect of pH and temperature on stability and kinetics of novel extracellular serine alkaline protease (70 kDa). Int. J. Biol. Macromol. 2013, 54, 1–8. [Google Scholar] [CrossRef]
  60. Jedynak, L.; Kowalska, J.; Harasimowicz, J.; Golimowskim, J. Speciation analysis of arsenic in terrestrial plants from arsenic contaminated area. Sci. Total Environ. 2009, 407, 945–952. [Google Scholar] [CrossRef]
  61. Singh, S.K.; Singh, S.K.; Tripathi, V.R.; Khare, S.K.; Garg, S.K. A novel psychrotrophic, solvent tolerant Pseudomonas putida SKG-1 and solvent stability of its psychro-thermoalkalistable protease. Process. Biochem. 2011, 46, 1430–1435. [Google Scholar] [CrossRef]
  62. Sánchez-Ramos, I.; Hernandez, C.A.; Castanera, P.; Ortego, F. Proteolytic activities in body and faecal extracts of the storage mite, Acarus farris. Med. Vet. Entomol. 2004, 18, 378–386. [Google Scholar] [CrossRef]
  63. Milton, D.K.; Chawla, R.K. Cotton dust contains proteolytic and elastolytic enzymes not inhibited by alpha-1-proteinase inhibitor. Am. J. Ind. Med. 1986, 9, 247–260. [Google Scholar] [CrossRef]
  64. Rudenskaya, G.E.; Pupov, D.V. Cysteine proteinases of microorganisms and viruses. Biochemistry (Moscow) 2011, 73, 1–13. [Google Scholar] [CrossRef]
  65. Shankar, S.; Rao, M.; Laxman, R.S. Purification and characterization of an alkaline protease by a new strain of Beauveria sp. Process. Biochem. 2011, 46, 579–585. [Google Scholar] [CrossRef]
  66. Willey, J.M.; Sherwood, L.; Woolverton, C.J. Prescott’s Microbiology; McGraw-Hill: New York, NY, USA, 2011; Volume 7. [Google Scholar]
  67. Thebti, W.; Riahi, Y.; Belhadj, O. Purification and Characterization of a New Thermostable, Haloalkaline, Solvent Stable, and Detergent Compatible Serine Protease fromGeobacillus toebiiStrain LBT 77. BioMed Res. Int. 2016, 2016, 1–8. [Google Scholar] [CrossRef] [Green Version]
  68. Kamran, A.; Rehman, H.U.; Qader, S.A.U.; Baloch, A.H.; Kamal, M. Purification and characterization of thiol dependent, oxidation-stable serine alkaline protease from thermophilic Bacillus sp. J. Genet. Eng. Biotechnol. 2015, 13, 59–64. [Google Scholar] [CrossRef] [Green Version]
  69. Hakim, A.; Bhuiyan, F.R.; Iqbal, A.; Emon, T.H.; Ahmed, J.; Azad, A.K. Production and partial characterization of dehairing alkaline protease from Bacillus subtilis AKAL7 and Exiguobacterium indicum AKAL11 by using organic municipal solid wastes. Heliyon 2018, 4, e00646. [Google Scholar] [CrossRef]
  70. Ramadan, M.F.; Mahmoud, M.G.; El Shebwy, K.; Abdel-Aziz, M.S. Purification and characterization of two thermostable protease fractions from Bacillus megaterium. J. Genet. Eng. Biotechnol. 2013, 11, 103–109. [Google Scholar] [CrossRef] [Green Version]
  71. Thakur, N.; Kumar, A.; Sharma, A.; Bhalla, T.C.; Kumar, D. Purification and characterization of alkaline, thermostable and organic solvent stable protease from a mutant of Bacillus sp. Biocatal. Agric. Biotechnol. 2018, 16, 217–224. [Google Scholar] [CrossRef]
  72. Sai-Ut, S.; Benjakul, S.; Sumpavapol, P.; Kishimura, H. Purification and Characterization of Extracellular Gelatinolytic Protease fromBacillus Amyloliquefaciens H11. J. Food Biochem. 2015, 39, 119–128. [Google Scholar] [CrossRef]
  73. Lakshmi, B.K.; Kumar, D.M.; Hemalatha, K.P. Purification and characterization of alkaline protease with novel properties from Bacillus cereus strain S8. J. Genet. Eng. Biotechnol. 2018, 16, 295–304. [Google Scholar] [CrossRef]
  74. Maruthiah, T.; Esakkiraj, P.; Prabakaran, G.; Palavesam, A.; Immanuel, G. Purification and characterization of moderately halophilic alkaline serine protease from marine Bacillus subtilis AP-MSU 6. Biocatal. Agric. Biotechnol. 2013, 2, 116–119. [Google Scholar] [CrossRef]
  75. Guleria, S.; Walia, A.; Chauhan, A.; Shirkot, C.K. Purification and characterization of detergent stable alkaline protease fromBacillus amyloliquefaciensSP1 isolated from apple rhizosphere. J. Basic Microbiol. 2015, 56, 138–152. [Google Scholar] [CrossRef]
  76. Klomklao, S.; Benjakul, S.; Kishimura, H.; Chaijan, M. Extraction, purification and properties of trypsin inhibitor from Thai mung bean (Vigna radiata (L.) R. Wilczek). Food Chem. 2011, 129, 1348–1354. [Google Scholar] [CrossRef]
  77. Sinha, R.; Khare, S.K. Thermostable Proteases. In Thermophilic Microbes in Environmental and Industrial Biotechnology; Satyanarayana, T., Littlechild, J., Kawarabayasi, Y., Eds.; Springer: Dordrecht, The Netherlands, 2013; pp. 859–880. [Google Scholar]
  78. Kumar, S.; Nussinov, R. Salt bridge stability in monomeric proteins 1 1Edited by J. M. Thornton. J. Mol. Biol. 1999, 293, 1241–1255. [Google Scholar] [CrossRef]
  79. Ozawa, T.; Hakamada, Y.; Hatada, Y.; Kobayashi, T.; Shirai, T.; Ito, S. Thermostabilization by replacement of specific residues with lysine in a Bacillus alkaline cellulase: Building a structural model and implications of newly formed double intrahelical salt bridges. Protein Eng. 2001, 14, 501–504. [Google Scholar] [CrossRef] [Green Version]
  80. Szilágyi, A.; Závodszky, P. Structural differences between mesophilic, moderately thermophilic and extremely thermophilic protein subunits: Results of a comprehensive survey. Structure 2000, 8, 493–504. [Google Scholar] [CrossRef] [Green Version]
  81. Kumar, S.; Dangi, A.K.; Shukla, P.; Baishya, D.; Khare, S.K. Thermozymes: Adaptive strategies and tools for their biotechnological applications. Bioresour. Technol. 2019, 278, 372–382. [Google Scholar] [CrossRef]
  82. Farhadian, S.; Asoodeh, A.; Lagzian, M. Purification, biochemical characterization and structural modeling of a potential htrA-like serine protease from Bacillus subtilis DR8806. J. Mol. Catal. B Enzym. 2015, 115, 51–58. [Google Scholar] [CrossRef]
  83. Mothe, T.; Reddy, S.V. Production, purification and characterization of a thermotolerant alkaline serine protease from a novel species Bacillus caseinilyticus. 3 Biotech 2016, 6, 53. [Google Scholar] [CrossRef] [Green Version]
  84. Theron, L.W.; Divol, B. Microbial aspartic proteases: Current and potential applications in industry. Appl. Microbiol. Biotechnol. 2014, 98, 8853–8868. [Google Scholar] [CrossRef]
  85. Olajuyigbe, F.M.; Falade, A.M. Purification and partial characterization of serine alkaline metalloprotease from Bacillus brevis MWB-01. Bioresour. Bioprocess. 2014, 1, 8. [Google Scholar] [CrossRef]
  86. Gulmus, E.O.; Görmez, A. Characterization and biotechnological application of protease from thermophilic Thermomonas haemolytica. Arch. Microbiol. 2019, 202, 153–159. [Google Scholar] [CrossRef]
  87. Durham, D.R.; Stewart, D.B.; Stellwag, E.J. Novel alkaline- and heat-stable serine proteases from alkalophilic Bacillus sp. strain GX6638. J. Bacteriol. 1987, 169, 2762–2768. [Google Scholar] [CrossRef] [Green Version]
  88. Rick, W. Methoden der Enzymatischen Analyse, 3rd ed.; Bergmeyer, H.U., Ed.; Verlag Chemic: Weinheim, Germany, 1974; pp. 1046–1056. [Google Scholar]
  89. Reimerdes, E.H.; Klostermeyer, H. substrates. In Methods in Enzymology; Academic Press: Cambridge, MA, USA, 1976; Volume 45, pp. 26–28. [Google Scholar]
  90. Thomas, T.D.; Pritchard, G.G. Proteolytic enzymes of dairy starter cultures. FEMS Microbiol. Lett. 1987, 46, 245–268. [Google Scholar] [CrossRef]
  91. Kumar, C.; Takagi, H. Microbial alkaline proteases. Biotechnol. Adv. 1999, 17, 561–594. [Google Scholar] [CrossRef]
  92. Rao, C.S.; Sathish, T.; Ravichandra, P.; Prakasham, R. Characterization of thermo- and detergent stable serine protease from isolated Bacillus circulans and evaluation of eco-friendly applications. Process. Biochem. 2009, 44, 262–268. [Google Scholar] [CrossRef]
  93. Bernhoft, R.A. Mercury Toxicity and Treatment: A Review of the Literature. J. Environ. Public Heal. 2011, 2012, 1–10. [Google Scholar] [CrossRef]
  94. Ibrahim, A.S.S.; Elbadawi, Y.B.; El-Tayeb, M.A.; Al-Maary, K.S.; Maany, D.A.F.; Ibrahim, S.S.S.; Elagib, A.A. Alkaline serine protease from the new halotolerant alkaliphilic Salipaludibacillus agaradhaerens strain AK-R: Purification and properties. 3 Biotech 2019, 9, 1–11. [Google Scholar] [CrossRef]
  95. Hammami, A.; Hamdi, M.; Abdelhedi, O.; Jridi, M.; Nasri, M.; Bayoudh, A. Surfactant- and oxidant-stable alkaline proteases from Bacillus invictae: Characterization and potential applications in chitin extraction and as a detergent additive. Int. J. Biol. Macromol. 2017, 96, 272–281. [Google Scholar] [CrossRef]
  96. El-Eskafy, H.; Abdel-Hamid, M.; Abbas, R.; Hamza, H.; Zanaty, A.E. Production of Thermostable Protease from “Bacillus amyloliquefaciens”: An Egyptian Marine Isolate. J. Agric. Chem. Biotechnol. 2016, 7, 227–234. [Google Scholar] [CrossRef] [Green Version]
  97. Anwar, A.; Saleemuddin, M. Alkaline protease from Spilosoma obliqua: Potential applications in bio-formulations. Biotechnol. Appl. Biochem. 2000, 31, 85–89. [Google Scholar] [CrossRef]
  98. Takami, H.; Kobayashi, T.; Aono, R.; Horikoshi, K. Molecular cloning, nucleotide sequence and expression of the structural gene for a thermostable alkaline protease from Bacillus sp. no. AH-101. Appl. Microbiol. Biotechnol. 1992, 38, 101–108. [Google Scholar] [CrossRef]
  99. Ghorbel-Frikha, B.; Sellami-Kamoun, A.; Fakhfakh, N.; Haddar, A.; Manni, L.; Nasri, M. Production and purification of a calcium-dependent protease from Bacillus cereus BG1. J. Ind. Microbiol. Biotechnol. 2005, 32, 186–194. [Google Scholar] [CrossRef]
  100. Horn, S.J.; Sørlie, M.; Vaaje-Kolstad, G.; Norberg, A.L.; Synstad, B.; Vårum, K.M.; Eijsink, V.G.H. Comparative studies of chitinases A, B and C fromSerratia marcescens. Biocatal. Biotransform. 2006, 24, 39–53. [Google Scholar] [CrossRef]
  101. Emran, M.A.; Ismail, S.A.; Hashem, A.M. Production of detergent stable thermophilic alkaline protease by Bacillus licheniformis ALW1. Biocatal. Agric. Biotechnol. 2020, 26, 101631. [Google Scholar] [CrossRef]
  102. Gaonkar, S.K.; Furtado, I.J. Characterization of Extracellular Protease from the Haloarcheon Halococcus sp. Strain GUGFAWS-3 (MF425611). Curr. Microbiol. 2020, 77, 1024–1034. [Google Scholar] [CrossRef]
  103. Hooft, R.W.; Sander, C.; Vriend, G. Objectively judging the quality of a protein structure from a Ramachandran plot. Bioinformatics 1997, 13, 425–430. [Google Scholar] [CrossRef] [Green Version]
  104. Jain, C.K.; Gupta, M.; Prasad, Y.; Wadhwa, G.; Sharma, S.K. Homology modeling and protein engineering of alkane monooxygenase in Burkholderia thailandensis MSMB121: In silico insights. J. Mol. Model. 2014, 20, 1–11. [Google Scholar] [CrossRef]
  105. Prajapat, R.; Marwal, A.; Gaur, R.K. Recognition of Errors in the Refinement and Validation of Three-Dimensional Structures of AC1 Proteins of Begomovirus Strains by Using ProSA-Web. J. Viruses 2014, 2014, 1–6. [Google Scholar] [CrossRef]
  106. Fakhfakh, N.; Haddar, A.; Hmidet, N.; Frikha, F.; Nasri, M. Application of statistical experimental design for optimization of keratinases production by Bacillus pumilus A1 grown on chicken feather and some biochemical properties. Process. Biochem. 2010, 45, 617–626. [Google Scholar] [CrossRef]
Figure 1. Map showing lower Dachigam National Park.
Figure 1. Map showing lower Dachigam National Park.
Biomolecules 11 00117 g001
Figure 2. Representative image showing a preliminary screening of some isolated bacterial strains for protease on skim-milk agar.
Figure 2. Representative image showing a preliminary screening of some isolated bacterial strains for protease on skim-milk agar.
Biomolecules 11 00117 g002
Figure 3. Qualitative screening of bacterial strain, HM48 for protease production using different substrates (1% w/v). Hydrolytic clearance zone on (a) casein; (b) gelatin, and (c) skim-milk agar plates.
Figure 3. Qualitative screening of bacterial strain, HM48 for protease production using different substrates (1% w/v). Hydrolytic clearance zone on (a) casein; (b) gelatin, and (c) skim-milk agar plates.
Biomolecules 11 00117 g003
Figure 4. Molecular identification of B. amyloliquefaciens, HM48 based on its 16S rRNA gene sequence. (A) Representative electrophoretic image of amplicon; (B) neighbor-joining phylogenetic tree.
Figure 4. Molecular identification of B. amyloliquefaciens, HM48 based on its 16S rRNA gene sequence. (A) Representative electrophoretic image of amplicon; (B) neighbor-joining phylogenetic tree.
Biomolecules 11 00117 g004
Figure 5. Representative image of SDS–PAGE of (A) crude and (B) purified protease enzyme from B. amyloliquefaciens, HM48 (lane 1—crude enzyme; lane 2—protein marker; lane 3—purified enzyme).
Figure 5. Representative image of SDS–PAGE of (A) crude and (B) purified protease enzyme from B. amyloliquefaciens, HM48 (lane 1—crude enzyme; lane 2—protein marker; lane 3—purified enzyme).
Biomolecules 11 00117 g005
Figure 6. Evaluation of activity/stability of protease from B. amyloliquefaciens HM48 with varying (A-enzyme activity), (B-enzyme stability) temperature, (C) pH, and (D) substrates.
Figure 6. Evaluation of activity/stability of protease from B. amyloliquefaciens HM48 with varying (A-enzyme activity), (B-enzyme stability) temperature, (C) pH, and (D) substrates.
Biomolecules 11 00117 g006
Figure 7. Effect of metal ions, at (A) 1.0, (B) 5.0 and (C) 10 mM final concentrations and (D) EDTA (1, 5 and 10 mM) in the presence of 10 mM final concentration of metal ion, MnCl2 on the activity of protease from B. amyloliquefaciens HM48.
Figure 7. Effect of metal ions, at (A) 1.0, (B) 5.0 and (C) 10 mM final concentrations and (D) EDTA (1, 5 and 10 mM) in the presence of 10 mM final concentration of metal ion, MnCl2 on the activity of protease from B. amyloliquefaciens HM48.
Biomolecules 11 00117 g007
Figure 8. Evaluation of activity of protease from B. amyloliquefaciens HM48 in the presence of various surfactants and an oxidizing agent at (A) 0.5%, and (B) 1.0% final concentration; (C) different organic solvents; (D) varying casein concentration and (E) Lineweaver–Burk plot.
Figure 8. Evaluation of activity of protease from B. amyloliquefaciens HM48 in the presence of various surfactants and an oxidizing agent at (A) 0.5%, and (B) 1.0% final concentration; (C) different organic solvents; (D) varying casein concentration and (E) Lineweaver–Burk plot.
Biomolecules 11 00117 g008
Figure 9. Molecular identification of B. amyloliquefaciens HM48 based on its protease gene sequence. (A) Representative electrophoretic image of amplicon; (B) neighbor-joining phylogenetic tree.
Figure 9. Molecular identification of B. amyloliquefaciens HM48 based on its protease gene sequence. (A) Representative electrophoretic image of amplicon; (B) neighbor-joining phylogenetic tree.
Biomolecules 11 00117 g009
Figure 10. (A) Refined model of B. amyloliquefaciens HM48 subtilisin-like protease. The model was generated through GalaxyTBM, refined through GalaxyRefine and rendering was done with the help of UCSF Chimera; (B) Ramachandran plot of GalaxyRefine model of B. amyloliquefaciens HM48 subtilisin-like protease. As per the plot statistics, the model quality is quite good and is thus highly reliable for further studies; (C) verify 3D analysis of GalaxyRefine model of B. amyloliquefaciens HM48 subtilisin-like protease. One hundred percent of the residues displayed average score greater than or equal to 0.2. The rule of thumb is that a minimum of 80% of the residues should possess an averaged 3D-1D score of ≥ 0.2. This further indicates the compatibility between 3D structure and primary structure (amino acid residue sequence).
Figure 10. (A) Refined model of B. amyloliquefaciens HM48 subtilisin-like protease. The model was generated through GalaxyTBM, refined through GalaxyRefine and rendering was done with the help of UCSF Chimera; (B) Ramachandran plot of GalaxyRefine model of B. amyloliquefaciens HM48 subtilisin-like protease. As per the plot statistics, the model quality is quite good and is thus highly reliable for further studies; (C) verify 3D analysis of GalaxyRefine model of B. amyloliquefaciens HM48 subtilisin-like protease. One hundred percent of the residues displayed average score greater than or equal to 0.2. The rule of thumb is that a minimum of 80% of the residues should possess an averaged 3D-1D score of ≥ 0.2. This further indicates the compatibility between 3D structure and primary structure (amino acid residue sequence).
Biomolecules 11 00117 g010
Figure 11. (A) ProSA-web analysis of the refined model of B. amyloliquefaciens HM48 subtilisin-like protease. Z-score of this model (−7.34) predicted by the ProSA program was found to be well within the range typically observed for proteins (native) of similar size. Z-scores lying outside the range typical for native proteins specify inaccurate structures; (B) Ramachandran plot of the refined model of bovine β-casein authenticates the high stereochemical quality of this model. For a good quality model above 90% of the residues shall fall in most favored regions, and in our case, 95.1% of the residues were found in the defined regions, clearly indicating the quality of the model to be good; (C) ProSA-web z-score plot of β-casein (bovine) model. Z-score is indicative of the general model quality, and if this score of the model falls within the range typically displayed by native proteins bearing similar size, the model is expected to be correct, and if the Z-score of the model violates the defined range, then such model is probably erroneous. From the plot, it is quite evident that the Z-score of the β-casein model having 224 amino acid residues is −3.67 and falls within the confines of the reference range.
Figure 11. (A) ProSA-web analysis of the refined model of B. amyloliquefaciens HM48 subtilisin-like protease. Z-score of this model (−7.34) predicted by the ProSA program was found to be well within the range typically observed for proteins (native) of similar size. Z-scores lying outside the range typical for native proteins specify inaccurate structures; (B) Ramachandran plot of the refined model of bovine β-casein authenticates the high stereochemical quality of this model. For a good quality model above 90% of the residues shall fall in most favored regions, and in our case, 95.1% of the residues were found in the defined regions, clearly indicating the quality of the model to be good; (C) ProSA-web z-score plot of β-casein (bovine) model. Z-score is indicative of the general model quality, and if this score of the model falls within the range typically displayed by native proteins bearing similar size, the model is expected to be correct, and if the Z-score of the model violates the defined range, then such model is probably erroneous. From the plot, it is quite evident that the Z-score of the β-casein model having 224 amino acid residues is −3.67 and falls within the confines of the reference range.
Biomolecules 11 00117 g011
Figure 12. (A) Docked complex of β-casein and B. amyloliquefaciens HM48 subtilisin-like protease rendered with the help of UCSF Chimera software. Green colored region indicates β-casein; the remaining color designates the B. amyloliquefaciens HM48 subtilisin-like serine protease; (B) residues of β-casein (chain A) interacting with B. amyloliquefaciens HM48 subtilisin-like protease (chain B). Twenty residues of this casein interact with 22 residues of HM48 subtilisin-like protease; (C) detailed account of interactions between β-casein residues (chain A) and B. amyloliquefaciens HM48 subtilisin-like protease (chain B). Following docking with HawkDock, the interaction profile was generated with the help of PDBsum.
Figure 12. (A) Docked complex of β-casein and B. amyloliquefaciens HM48 subtilisin-like protease rendered with the help of UCSF Chimera software. Green colored region indicates β-casein; the remaining color designates the B. amyloliquefaciens HM48 subtilisin-like serine protease; (B) residues of β-casein (chain A) interacting with B. amyloliquefaciens HM48 subtilisin-like protease (chain B). Twenty residues of this casein interact with 22 residues of HM48 subtilisin-like protease; (C) detailed account of interactions between β-casein residues (chain A) and B. amyloliquefaciens HM48 subtilisin-like protease (chain B). Following docking with HawkDock, the interaction profile was generated with the help of PDBsum.
Biomolecules 11 00117 g012
Figure 13. Evaluation of protease from B. amyloliquefaciens HM48 for industrial and environmental application. (A) Detergent compatibility; (B) wash performance analysis (blood and egg-yolk stain removal, and (C) waste management (chicken feather degradation).
Figure 13. Evaluation of protease from B. amyloliquefaciens HM48 for industrial and environmental application. (A) Detergent compatibility; (B) wash performance analysis (blood and egg-yolk stain removal, and (C) waste management (chicken feather degradation).
Biomolecules 11 00117 g013
Table 1. Qualitative screening of strain HM48 for its proteolytic activity on different proteinaceous substrates (pH 7.0).
Table 1. Qualitative screening of strain HM48 for its proteolytic activity on different proteinaceous substrates (pH 7.0).
Proteinaceous SubstrateMean Diameter of Colony, d (in mm)Mean Diameter of Hydrolytic Zone, D (in mm)D/d (in mm)
Casein3289.33
Gelatin10333.3
Skim milk7263.71
Table 2. Morphological characteristics of strain, HM48.
Table 2. Morphological characteristics of strain, HM48.
SizeShapeMarginElevationTextureAppearanceColorTransparencyGram’s ReactionCell ShapeCell Arrangement
ModerateIrregularUndulateUmbonateRoughDullCreamOpaquePositiveBacilliStreptobacilli
Table 3. Carbohydrate utilization and antibiotic susceptibility tests of strain, HM48.
Table 3. Carbohydrate utilization and antibiotic susceptibility tests of strain, HM48.
CarbohydrateResultAntibioticConcentration *Result
Dextrose utilizationPositiveAmoxiclav (amoxicillin/clavulanic acid) (AMC)30 mcg (20/10 mcg)Sensitive
Esculin hydrolysisPositiveCefoxitin (Cephoxitin) (CX)30 mcgSensitive
Cellobiose utilizationPositiveCeftazidime (CAZ)30 mcgSensitive
Citrate utilizationPositiveCephalothin (CEP)30 mcgSensitive
Inositol utilizationPositiveClindamycin (CD)2 mcgSensitive
Inulin utilizationPositiveErythromycin (E)15 mcgSensitive
Mannitol utilizationPositiveGentamicin (GEN)10 mcgSensitive
Sodium gluconate utilizationPositiveOfloxacin (OF)5 mcgSensitive
Sorbitol utilizationPositiveOxacillin (OX)1 mcgSensitive
Sucrose utilizationPositivePenicillin-G (P)10 unitsSensitive
Trehalose utilizationPositiveTeicoplanin (TEI)30 mcgSensitive
Vancomycin (VA)30 mcgSensitive
* Concentration of antibiotics as per Clinical and Laboratory Standards (CLSI, https://clsi.org/).
Table 4. Purification process of enzyme produced by B. amyloliquefaciens HM48.
Table 4. Purification process of enzyme produced by B. amyloliquefaciens HM48.
Purification ProcessTotal Activity (U mL−1)Total Protein (mg mL−1)Specific Activity (U mg−1)Yield (%)Purification Fold
Culture supernatant313.1829.7010.541001
Ammonium sulfate precipitation215.639.4522.8168.92.2
Dialysis160.064.1538.5451.13.7
Gel filtration chromatography114.491.8861.0536.65.8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mushtaq, H.; Jehangir, A.; Ganai, S.A.; Farooq, S.; Ganai, B.A.; Nazir, R. Biochemical Characterization and Functional Analysis of Heat Stable High Potential Protease of Bacillus amyloliquefaciens Strain HM48 from Soils of Dachigam National Park in Kashmir Himalaya. Biomolecules 2021, 11, 117. https://doi.org/10.3390/biom11010117

AMA Style

Mushtaq H, Jehangir A, Ganai SA, Farooq S, Ganai BA, Nazir R. Biochemical Characterization and Functional Analysis of Heat Stable High Potential Protease of Bacillus amyloliquefaciens Strain HM48 from Soils of Dachigam National Park in Kashmir Himalaya. Biomolecules. 2021; 11(1):117. https://doi.org/10.3390/biom11010117

Chicago/Turabian Style

Mushtaq, Hina, Arshid Jehangir, Shabir Ahmad Ganai, Saleem Farooq, Bashir Ahmad Ganai, and Ruqeya Nazir. 2021. "Biochemical Characterization and Functional Analysis of Heat Stable High Potential Protease of Bacillus amyloliquefaciens Strain HM48 from Soils of Dachigam National Park in Kashmir Himalaya" Biomolecules 11, no. 1: 117. https://doi.org/10.3390/biom11010117

APA Style

Mushtaq, H., Jehangir, A., Ganai, S. A., Farooq, S., Ganai, B. A., & Nazir, R. (2021). Biochemical Characterization and Functional Analysis of Heat Stable High Potential Protease of Bacillus amyloliquefaciens Strain HM48 from Soils of Dachigam National Park in Kashmir Himalaya. Biomolecules, 11(1), 117. https://doi.org/10.3390/biom11010117

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop