t i l e A Thermally Stable Acidic Chitinase from Paenibacillus sp. Y412MC10: Molecular Characterization and its Structural Modeling

Chitinase is an enzyme that breaks down chitin polymer’s β-1, 4-glycosidic bonds. Chitinases having high thermostability makes them suitable to use in different industrial sectors. A 1464-bp chitinase gene, from Paenibacillus sp . Y412MC10 was cloned and produced in E. coli BL21 cells using the pRSETB vector. Gel filtration and ion exchange chromatography were used to purify the recombinant 52 kDa protein. Purified chitinase showed optimum activity at 60°C and pH 5.5 with colloidal chitin breakdown. Enzyme showed stable residual activity within pH range of 4.5–6.5 and >70% thermal stability upto 60°C for 2.5 h. The activity of chitinase increased in the presence of Mn +2 , SDS, and methanol. Chitinase has K m and V max values of 2.287 mg/ml and 6.784 µM/min, respectively, towards colloidal chitin. The protein- ligand docking analysis and molecular dynamic simulation indicated that the interactions of amino acids Asn 204, Glu 277, Leu 278, Asp 312, Glu 314, Gln 372, Tyr 374, Trp 463, Trp 467 from GH18 catalytic domain with (NAG) 2 molecule are involved in the enzyme catalytic mechanism. This is the first time a thermostable chitinase from Paenibacillus sp. Y412MC10 has been cloned, expressed heterologously, and purified.


INTRODUCTION
C hitin composed of β-(1, 4) linked repeating units of N-acetyl-D-glucosamine (GlcNAc) and is widespread in nature. Chitin is a protein that is made up of β-(1, 4) connected repeating units of N-acetyl-D-glucosamine (GlcNAc). It is the second most prevalent polymer in nature, following cellulose, and can be found in a variety of creatures, including marine shells, exoskeletons, and gut linings of worms and insects, arthropods, as well as the fungal cell wall (Stoykov et al., 2015). Chitin can be found in two different crystalline forms: α-chitin, the most common kind of chitin and a rigid form and have antiparallel chain arrangement, and β-chitin, having open structure, with chains arrange in parallel fashion (Komi and Hamblin, 2016).
Tons of chitin polysaccharide are produced every year in nature, and shellfish debris (10,000 tons) are produced each year that is not decomposed properly (Rameshthangam et al., 2018;Wang et al., 2018). So, thorough decomposition of chitin waste is required and create soluble and low molecular weight useful GlcNAc. These soluble products have been utilizing in agriculture, medicine, nutrition, food etc. (Liaqat and Eltem, 2018). Traditionally, different chemical methods have been applied for the degradation which not only gives low yield but also cause environmental pollutions. Enzymatic chitin breakdown with chitinolytic enzymes is now recognized to be of major economic and environmental significance (Gao et al., 2018).
The primary degraders of chitin i.e. microorganisms, secrete several different chitin degrading enzymes (Zhang et al., 2018b). Endochitinases and exochitinases are two types of chitinase enzymes. Endochitinases break the chitin chains at random locations internally whereas exo-chitinases are further split into chitobiosidases and O n l i n e

F i r s t A r t i c l e
1-4-glucosaminidases that cleave from non-reducing and reducing end, respectively (Hamid et al., 2013;Zhou et al., 2019). Depending upon the sequence of amino acid and the primary structure, chitinolytic enzymes consist of three families; GH18, GH 19 and GH 20 (Lombard et al., 2014). Families 18 and 19 are called chitinases, because of their involvement in the breakdown of chitin polsaccharide (Vaaje-Kolstad et al., 2013). The most diverse members of GH family 18 chitinases are bacteria, fungi and insect. Bacterial chitinases have environmental adaptability, high thermostability and fast and easy expression and because of this, can be easily engineered within a labs (Le and Yang, 2019). These bacterial chitinases can be isolated from different sources i.e. shellfish waste, soil, and hot springs (Hamid et al., 2013). Chitinases producing bacteria include the genera Serratia, Vibrio, Bacillus, Chromobacterium, Pseudomonas, Paenibacillus, etc. Most of these bacteria produce chitinases that work in moderate pH and temperature (Bouacem et al., 2018). But chitinases having high thermostability and tolerability makes them suitable to work in different industrial applications (Karthik et al., 2015).
In this study, we are describing the cloning, heterologous expression and characterization of chitinase from Paenibacillus sp. Y412MC10, a novel Paenibacillus lautus strain, which is derived from Yellowstone National Park (in obsidian hot Spring) for the first time. The properties of this enzyme were evaluated to provide a valuable data for efficient utilization in the breaking of chitin.

Bioinformatics analysis
Database for Automated Carbohydrate-active enzyme ANnotation (DbCAN) was used to predict the domains (Yin et al., 2012). ClustalW was used for multiple sequence alignment (Higgins and Sharp, 1988). Putative chitinase pI value and molecular weight were determined from ExPASy (Compute pI/Mw) (Gasteiger et al., 2005). Homology for putative chitinase with the other bacterial species was determined by using pBLAST. SWISS-MODEL software was used for the prediction of 3-D structure of chitinase enzyme (Waterhouse et al., 2018).

Cloning of chitinase gene
The putative, multi-domain chitinase (Plchi1), which has 1464 bp (GenBank: ACX65998.1) and codes for a 487-amino-acid protein, was isolated from Paenibacillus sp. Y412MC10 and cloned into the pRSETB expression vector (Invitrogen, Carlsbad, CA, USA). Chitinase gene from Paenibacillus sp. Y412MC10 was synthesized from Bio Basic (1509 bp), including the BsmI and HindIII restriction sites and, the flanking region at 3' and 5' end that was also complementary to the region already present in the pRSETB expression vector to achieve restrictionfree strategy of cloning. The gene of interest was fused to the pRSET B vector, which had been linearized, using the GeneArt™ Seamless Cloning and Assembly Kit (Thermo scientific). The final construct (pRSETB-chitinase) was introduced into competent E. coli BL21 (DE3) cells (Invitrogen). DNA sequencing was used to confirm the construct.

Expression and purification of chitinase
Overnight cultures grown from -80°C stocks of E. coli BL21 (DE3) cells harboring pRSETB-chitinase construct were used to inoculate 500 ml of LB medium having ampicillin (50 µg/ml). Without induction, the culture was incubated with shaking for 20 h at 37°C. The periplasmic protein fraction was isolated by the use of osmotic shock method with some adjustments after centrifugation (6,000 rpm) at 4°C for 10 min (Tuveng et al., 2017). The cells were resuspended in 50 ml of cold spheroplast buffer (17.1 percent (wt/vol) sucrose, 0.5 M EDTA, 1 M Tris HCl pH 8, 50 mM PMSF) and cooled on ice for 10 min to prepare periplasmic extracts. The pellet was made by centrifuging the supernatant for 12 min at 4°C at 6,000 rpm before warming it to room temperature (RT). The pellet was suspended in 30 ml of cool dH 2 0, then 625 µl of 20 mM MgCl 2 was added, and cells were incubated on ice upto1 minute before pelleted down by centrifugation at 8,000 rpm for 10 min at 4°C. The supernatant, which included mature enzyme and was sterilized by filtration, was then tested on a 12% sodium dodecyl sulphate polyacrylamide gel electrophoresis (SDS-PAGE) gel (Laemmli, 1970).
Chitinase protein was isolated in two steps: first, gel filtration with a HiLoad 16/60 Superdex (GE Healthcare) equilibrated with 50 mM sodium acetate buffer solution at pH 5.5, and then anion exchange chromatography with a 5 mL of Q Sepharose FF column having pH 8.0 (50 mM Tris buffer) (GE Healthcare, Uppsala, Sweden). For protein elution, a linear salt (NaCl) gradient of 30 fractional volumes was used. SDS/PAGE at 12% validated the protein's purity. Using Amicon ultra centrifugal filters of MW cut-off (MWCO) 10,000 Da, the pure protein fractions were concentrated. The protein quantity was calculated using the Bradford test, which employed bovine serum albumin (BSA) as a standard (Bradford, 1976).

Colloidal chitin preparation
Murthy and Bleakley (2017) method was used for colloidal chitin preparation with some modifications. 5g of O n l i n e F i r s t A r t i c l e powdered chitin (crab shell) was progressively introduced into conc. HCl (60 ml) and maintained at room temperature (RT) for 1 h with continuous stirring. The mixture was then added to 200 ml of ice-cold 50 percent ethanol and stirred vigorously overnight at 4°C. Precipitates were collected by centrifugation at 10,000 rpm for 20 min at 4°C, and then rinsed with sterile water up to the point, when the pH of the colloidal chitin reached to pH 7.0. After that, it was air dried on Whatman filter paper and stored at 4°C in the dark.
Enzyme activity assay Miller (1959) revealed how the amount of N-acetylglucosamine (NAG) liberated from the substrate affects chitinase enzyme activity. In a test tube, a reaction mixture of adequately diluted enzyme (0.5 ml) and 1 percent (w/v) suspension of colloidal chitin (0.5 ml) in sodium acetate buffer (50 mM, pH 5.5) was incubated for 30 min in a water bath at 60°C along with appropriate blank. Addition of 3 mL DNS reagent in the mixture caused the stopping of the reaction, which was then boiled in a water bath for 10 min before being cooled at RT. A clear supernatant was obtained by centrifuging for 10 min at 1000 rpm. The reducing sugars liberated in the solution were quantified using a UV/visible light spectrophotometer (UV-2450; Shimadzu) set to 540nm. Chitinase of 1 unit (U) is defined as the amount of enzyme which produced reducing sugars corresponding to 1 micromole of N-acetyl glucosamine equivalents from colloidal chitin per minute under the assay conditions.

Biochemical characterization of chitinase
Using 1% chitin in colloidal form as a substrate, the effect of pH on pure chitinase activity was investigated. Sodium acetate solution (50 mM) with a pH range of 3.5-5.5, sodium phosphate buffer (50 mM) with pH range 6.0-7.5, and 50 mM glycine NaOH buffer with a pH range of 8.0-10.5 were used. The pH stability of chitinase enzyme was estimated by incubating the purified enzyme (without substrate) for 1 h at room temperature in buffer solutions ranging from 3.5 to 10.5 pH. After that, the residual enzyme activity was found using the usual procedure.
Heating a known quantity of purified chitinase enzyme at multiple temperatures (40 to 80°C) while using 1% substrate (colloidal chitin) produced in 50 mM sodium acetate buffer pH 5.5 was used to find its optimal temperature. Prior to adding the substrate, the thermal stability of isolated enzyme was tested by incubating in water bath for 1 h at various temperatures (40-80°C). Enzyme activity (residual) was further tested under standard conditions of temperature.

Substrate specificity and kinetic studies
Using 1% (w/v) of several substrates such as colloidal chitin, powdered α-chitin and β-chitin, Avicel, and carboxymethyl cellulose (CMC-Na), the specificity of substrates for recombinant chitinase was investigated. Then relative activity of enzyme was calculated by following standard protocol with the maximum activity assumed to be 100%.
Enzyme kinetics was studied by performing the experiments in triplicates with the varying concentration of colloidal chitin from 1 mg/ml to 12 mg/ml in sodium acetate buffer (pH 5.5). The enzyme's specific activity was then determined, and Michaels constants were calculated using a Lineweaver-Burk double reciprocal plot.

3D structure prediction and molecular docking study of protein
Molecular docking, a popular computer simulation approach, for determining the shape of a complex (receptorligand), in which mostly receptor is nucleic acid molecule or a protein, and the ligand is either another protein or a small molecule (Meng et al., 2011). In our study, appropriate binding of the enzyme with N-acetylglucosamine (a subunit of chitin) was confirmed by molecular docking. For this purpose, the 3D structures of N-acetylglucosamine (NAG) 2 was downloaded from the PDB database [PDB ID: 3N15] (Berman et al., 2002). Alpha Fold2 ColabFold was used to estimate the 3D structure of our enzyme (Mirdita et al., 2021). Five models were generated which were further validated through SAVES server (https://saves.mbi.ucla. edu/). The computational servers such as Ramachandran Plot (Gopalakrishnan et al., 2007), Verify 3D (Eisenberg et al., 1997) and ERRAT (Colovos and Yeates, 1993) were used for the evaluation of generated models and one model was selected for downstream analysis. CB-dock (Liu et al., 2020b) web server was used for protein-ligand interactional study which uses AutoDock vina to perform molecular docking (Trott and Olson, 2010). (NAG) 2 obtained from PDB database was taken as ligands and O n l i n e

F i r s t A r t i c l e
our protein as receptor. The best pose having the lowest binding energy was further chosen for analysis.

Molecular dynamic simulation
The molecular dynamics simulation approach was used to investigate the substrate's binding affinities with the protein at the atomic level. For this purpose, GROMACS version 5.0.5 was used to follow the dynamic pattern of docked complexes applying the CHARMM36 force field on proteins and CGenFF on drug-like compounds (Abraham et al., 2015;Yu et al., 2012). The system was equilibrated for 100 ps, at 300K followed by a 50 ns production phase. Eventually, MD trajectories were analyzed for the calculation of RMSF; root mean square fluctuations, RMSD; root mean square deviation, HB; number of hydrogen bonds and Coul-SR and LJ-SR; compound interaction energies.

Bioinformatics analysis
A full ORF encoding a 487-amino-acid protein was found in the chitinase gene. Chitinase is made up of three domains: The N-terminal ChBD3: chitin-binding domain, the Fn3D: fibronectin-type-III domain, and the C-terminal CatD: GH 18 catalytic domain, according to domain architectural study (Fig. 1A). The structure (3-D) of chitinase was predicted using software SWISS-MODEL and shown in Figure 1A.
Although, the amino acid sequences of various organisms differ, some essential amino acids in the GH18 domain are substantially conserved. These two sequences SXGG and DXXDXDXE, are considered as conserved motifs that plays an important role in binding and catalytic degradation of substrate, respectively ( Fig.  2A). Glutamic acid (E) is an acid/base, which allows oxygen atoms from sugar molecules to be protonated. The region between chitin binding domain and catalytic domain was found to be FN3 domain, which is found in many bacterial chitinases, amylases and cellulases enzymes. FN3 is an extracellular plasma protein that is involved in cell adhesion. It was considered to play a main role in the degradation of recalcitrant polysaccharides. FN3 domain shows conserved region with sequence TXYXFXXXAKDAXG, and is shown in Figure 2B. The putative chitin binding domain (ChtBD3) was compared with other characterized bacterial chitin binding domain and analysis showed that ChtBD3 is required for the specific recognition of crystalline and insoluble chitin. In the ChtBD3 of Paenibacillus and other chitinases, in chitin polymer, Tryptophan (W) and Tyrosine (Y) residues along with polar hydrogen groups have conservation in the genome and are expected to attach directly against the N-acetylglucosamine residues (pyranose ring form) (Fig.  2C). Homologues of the protein found in Protein Data Bank were determined and shown in Table I.

Cloning of chitinase gene
The safe and efficient method for production of biological products i.e. enzymes in bulk quantity is the recombination DNA technology. Considering the significance of enzyme, chitinase gene from Paenibacillus sp. Y412MC10, with open reading frame of 1464 bp, was synthesized from Bio Basic and cloned into pRSETB expression vector. BsmI and HindIII restriction enzymes were used for the confirmation of insert by digesting the purified recombinant plasmid. Fragments 2835 bp of pRSETB and chitinase gene of 1464 bp were observed, that is shown in Figure 3.

Expression and purification of recombinant chitinase
To gain a better understanding of an enzyme's mode of action and its specificity for different substrate, enzyme purification is required. The multi-domain chitinase gene was isolated from Paenibacillus sp. Y412MC10 strain and was expressed in competent cells of E. coli BL21 (DE3). Chromatography i.e. gel filtration and anion exchange were being used for the purification of protein. Chitinase is expected to have a molecular weight of 52 kDa, visualized on SDS-PAGE (Fig. 4), that is nearly close to the expected weight, calculated from EXPASY tool. The molecular weight of chitinases from diverse Paenibacillus species, to be between 37 and 154 kDa (Fu et al., 2014). As compared to the molecular weight of purified chitinase in this study, lower molecular weight chitinases have been purified from Paenibacillus pasadenensis NCIM (35 kDa) (Loni et al., 2014) and Streptomyces sp. CS501 (43 kDa) (Rahman et al., 2014). In contrast to the current work, greater molecular weight chitinase has been found in Paenibacillus barengoltzii (74 kDa) by Fu et al. (2014) and Paenibacillus ehimensis MA2012 (100 kDa) by Seo et al. (2016).

Biochemical characterization
One of the crucial variables that determines the speed at which an enzyme-catalyzed process takes place is temperature. To resist the enzymes up to thermal denaturation, the factors like ionic stability, hydrogen bonding and hydrophobic interactions within a polypeptide chain is of great importance (Younas et al., 2016). Effect of different range of temperatures (40 to 80°C) on chitinase activity were determined (Fig. 5A). With the increasing temperature from 40°C, the chitinase activity increased gradually and reached at its optimum activity at 60°C. The chitinase activity decreased clearly above 60°C. The findings of Ueda and Kurosawa (2015) for Paenibacillus thermoaerophilus chitinase are consistent with this investigation. However optimum temperature of chitinase isolated from Paenibacillus pasadenensis (Loni et al., 2014) is relatively low i.e. 37 o C while chitinase from other bacteria has a much higher optimal temperature, such as 85 o C for chitinase from Thermococcus kodakaraensis (Tanaka et al., 2001). Among the various properties of enzymes, stability of temperature is considered as a substantial tool for chitin industry and other industrial areas. Purified chitinase's thermostability was tested against the range of temperatures (40 to 80°C) (Fig. 5B). The enzyme was stable upto 60°C for 1.5 h of incubation. Afterwards, enzyme activity decreased gradually with increasing temperature but retained more than 70% activity when incubated at 40-60°C for 2.5 h. Enzyme activity suddenly decreased to 50% after 1 h incubation at 80°C. pH would be considered as an important factor in studying the tertiary and quaternary structure of enzymes. The change in pH causes the ionization of atoms and molecules of amino acids, thus change the structure and shape of protein. This change will damage the enzyme function and as a result, solubility, activity and stability of the enzyme was effected (Chaplin and Bucke, 1990). Chitinase obtained from present study had an optimum pH of 5.5, which has similar findings as Yang et al. (2016) who purified a chitinase from Paenibacillus barengoltzii. While chitinase obtained from Paenibacillus thermoaerophilus TC22-2b possessed pH optimum at 3.0 (Ueda and Kurosawa, 2015). On contrary, purified chitinase at pH 10, by Paenibacillus pasadenensis exhibited maximum activity (Loni et al., 2014). Because pure chitinase was most active in the acidic range in our study, it could be useful in the control of fungal diseases (Moore et al., 2004). When incubated in buffers consisting of various pH, chitinase maintains significant stability having ≥ 90% residual activity between pH 4.5-6.5 (Fig. 6B) and shows 100% residual activity on its optimal pH (pH 5.5). After pH 6.5, relative activity starts to decrease gradually. The pH range of some chitinases was even wider. Fu et al. (2014) purified a chitinase from Paenibacillus barengoltzii exhibiting stability from pH 4.0-9.0, while chitinase purified by Kim et al. (2017) maintained its activity in wider pH range i.e. 3-11. Metal ions operate as a cofactor in the catalytic process, which helps to stabilize the enzyme-substrate complex (Andreini et al., 2008). At various concentrations (1 mM to 5 mM), metal ions effect on the activity of purified chitinase is illustrated in Figure 7A. Ca +2 and Mn +2 increased enzyme activity (almost three times). The enhancement of chitinase activity by these metal ions in Paenibacillus chitinolyticus strain UMBR 0002 and Thermobifida fusca Tfu_0580 were reported by Liu et al. (2020a) and Yan and Fong (2018), respectively. The chitinase activity was also increased by Co +2 , Fe +2 ions but this is in the contrast with study obtained from Rahman et al. (2014) in which the enzyme is reported to be inhibited by Fe +2 and Co +2 . Various metal ions inhibited the enzyme and caused the destruction of the tertiary structure of protein which resulted in their inactivation. In the current study, In the presence of Cu +2 , chitinase activity was suppressed which is similar with results of chitinase from P. pasadenensis, in which Cu +2 metal ion hindered O n l i n e

F i r s t A r t i c l e
Molecular Characterization and Structural Modeling of Chitinase from Paenibacillus sp. 7 their enzyme activity as well (Loni et al., 2014). Deng et al. (2019) proposed that Cu +2 ion induces the autoxidation of cysteine, resulting in intramolecular disulfide bonds or sulfonic acid production, thus inhibiting chitinase function. Cd +2 has showed no effect on enzyme activity (Fig. 7A). The influence of chemical reagents at concentrations of 1 mM and 5 mM on pure chitinase activity is demonstrated in Figure 7B. It was found that certain reducing agents like SDS and DTT increased the activity of chitinase to 131.2% and 146.5%, respectively at 5mM concentration which is in opposite with the results obtained from chitinase of Paenibacillus sp. TKU052 (Doan et al., 2021). In the presence of a chelating reagent, EDTA, enzyme activity increases only little (115.6%) at 1mM concentration but this result is not consistent with chitinase from Streptomyces violascens (Gangwar et al., 2016) and Streptomyces sp. CS501 (Rahman et al., 2014). 37.8%, 18.3, 2.9% and 6.6% reduction was observed when purified enzyme was incubated in β-mercaptoethanol, urea sodium bisulphite and PMSF, respectively at 1 mM concentration and more decrease in enzyme activity i.e. 61 %, 72.8%, 18.5% and 22.2% was found at 5 mM concentration of β-mercaptoethanol, urea sodium bisulphite and PMSF, respectively.
In enzyme assays, organic solvents are utilized to make hydrophobic substrates soluble. Solvents may disturb the hydrophobic region of the enzymes so interfering enzyme structure and function. Furthermore, according to Jaouadi et al. (2013), polar solvents cause structural changes in enzymes, which can disrupt the substrate-enzyme active site interaction. Keeping in view such studies, various organic solvents to determine the impact on enzyme activity were used. The effect of several organic solvents on enzyme activity is demonstrated in Figure 7C. The activity was substantially increased in case of methanol i.e. 190%, 158.3%, respectively when used at 1% and 5% final concentration. However, chitinase activity decreased with increased percentage of solvent. Most of the organic solvents utilized in this study, such as twenty, ethanol, DMSO, isopropanol and triton X-100 had an inhibitory effect on enzyme activity, lowering it to 2.9%, 37.8%, 68.8%, 76.5%, 14.4% respectively, when used at 5% final concentration. Similar results were observed from the study of chitinase from Streptomyces sp. CS501 and Paenibacillus sp. TKU052 by Rahman et al. (2014) and Doan et al. (2021).

Substrate specificity and kinetic studies
Chitinase hydrolyzes substances with glycosidic linkages in a specific way. Specificity of various substrates for purified enzyme was determined. The residual activity of purified chitinase against various substrates are shown in Figure 8A. Colloidal chitin had the highest specific activity while considerable activity was also achieved when β-chitin substrate was used. Enzyme showed very little activity toward α-chitin as compared to β-chitin. is anti-parallel beta-sheets, strong forces of attraction between the molecules are present, due to which it is not easily hydrolyzed as compared to β-chitin (parallel betasheets) (Beier and Bertilsson, 2013;Martínez et al., 2014) and thus, showed very little activity. In the present study, enzyme shows activity towards Avicel (microcrystalline cellulose) due to the the binding of ChBD onto cellulose which can be attributed to interactions between aromatic residues in the enzyme and the pyranose ring of the polysaccharide (Kikkawa et al., 2011). But no catalytic activity was detected toward carboxymethyl cellulose because Avicel is crystalline cellulose which is insoluble in water whereas CMC-Na is soluble form of cellulose and showed no activity. In the instance of Paenibacillus chitinolyticus chitinase, a similar effect was seen (Liu et al., 2020a) and Paenibacillus sp. TKU052 (Doan et al., 2021). The catalytic efficiency of enzyme and its affinity for substrate is determined by its kinetic studies. Enzyme activity in different concentration of substrate were measured and The Michaelis constant was calculated using the Line weaver-Burk plot as shown in Fig. 8B. The plot was linear (R2= 0.990) over substrate concentration and the V max and K m values were found to be 6.784 µM/ min and 2.287 mg/ml, respectively. Fig. 9. The Ramachandran plot of predicted chitinase indicating total residues in the the most favored areas (red region) and additional additional permitted areas regions (yellow region) are 82.9 and 12.3%, respectively.

Structure prediction and docking study
To study the behavior of enzyme-substrate complex, molecular docking analysis was done in this study. AlphaFold2 ColabFold predicted five different structures of our protein. Ramachandran Plot, Verify-3D and ERRAT were used to further assess these models. One of the best predicted model among five showed different residues lies within different regions in Ramachandran plot analysis i.e. 82.9% (most favorable regions), 12.3% (additional allowed regions), 3.1% (generously allowed regions) and 1.7% (disallowed regions) (Fig. 9). Likewise, Verify-3D result showed 88.71% of the residues have averaged 3D-1D score >= 0.2, and ERRAT results predicted 86.25 the overall quality factor scores for the same model, which signify the consistency of the model predicted. The ligand-protein docking helped to analyze the molecular interactions between substrate and enzyme. CB-dock showed different binding modes of (NAG) 2 with our O n l i n e

F i r s t A r t i c l e
Molecular Characterization and Structural Modeling of Chitinase from Paenibacillus sp. 9 chitinase. The best pose of (NAG) 2 was selected which showed minimum binding energy of -8.2 (kcal/mol) and RMSD value as 0.00 (Table II). The best docked 2D and 3D structure represented in Figure 10A and B, respectively which clearly showed that the residues Asn 204, Glu 277, Leu 278, Asp 312, Glu 314, Gln 372, Tyr 374, Trp 463 and Trp 467 of chitinase are actively involved in binding with (NAG) 2 through strong hydrogen bonds.

MD simulation
During enzyme-substrate complex formation, to investigate the structural variations and to confirm whether the given complex is stable or not, RMSD plot of chitinase backbone was generated (Fig. 11A). The RMSD values in the first 5 ns showed fluctuations and reached to 0.6 Å but with the increase in the time, the RMSD values were found to be at constant rate of around 0.4 Å, which shows the overall stability of enzyme-substrate complex at this constant temperature of 310K. In docked complex of chitinase, to estimate the residual flexibility, RMSF fluctuations were computed (Fig. 11B). The RMSF is a measure of the difference between a particle's position and a reference position. Mostly the residues showed fluctuation below 0.5 Å which is also important for catalysis, apart from N-terminus, which showed a violent fluctuation between the residues window between 0-30 (up to 1.8 Å) and showed high fluctuation between residues 45-90. The overall RMSF values predicted that chitinase complexed with (NAG) 2 is in stable mode. During MD simulation, the presence of substrate causes the overall energy of the system to decrease which demonstrating that a more stabilized protein ligand in aqueous conditions was formed. The binding energies were calculated in terms of Lennard Jones Short-Range (LJ-SR) and Coul-SR interaction energies. The interaction energy of Coul-SR was -113.27 kJ/mol and LJ-SR energy value was -129.685 kJ/mol this much of lower energies indicate the stability of enzyme-substrate complex (Fig. 11C). Likewise, the consistency of hydrogen bonds (HB) up to 50 ns indicates the stability of enzyme-substrate complex (Fig. 11D). Gene cloning, expression and biochemical characterization of a novel chitinase of Paenibacillus sp. Y412MC10 were done. At pH 5.5 and 60°C, the recombinant chitinase enzyme was most active against colloidal chitin. Chitinase has a broad specificity of substrates and degrade colloidal chitin form and other crystalline substrates of chitin to produce chitooligosaccharides. Molecular dynamic simulations of chitinase in complex with (NAG) 2 for 50 ns run help in determining the active binding sites and their key residues roles in the regulation of an enzyme. There is a need to investigate its potential in biomass conversion and to study its biotechnological applications for future purposes. This is the first study that we are aware of; on the molecular characterization of thermally stable acidic recombinant chitinase from Paenibacillus sp. Y412MC10.

Statement of conflict of interest
The authors have declared no conflict of interest.