In silico characterization of alkaline proteases from different species of Aspergillus

Computational analysis to comprehend the structure-function properties of fibrinolytic enzymes from Bacillus spp for their efficient integration into industrial applications

Background

The fibrinolytic enzymes from Bacillus sp. are proposed as therapeutics in preventing thrombosis. Computational-based analyses of these enzymes' amino acid composition, basic physiological properties, presence of functional domain and motifs, and secondary and tertiary structure analyses can lead to developing a specific enzyme with improved catalytic activity and other properties that may increase their therapeutic potential.

Methods

The nucleotide sequences of fibrinolytic enzymes produced by the genus Bacillus and its corresponding protein sequences were retrieved from the NCBI database and aligned using the PRALINE programme. The varied physiochemical parameters and structural and functional analysis of the enzyme sequences were carried out with the ExPASy-ProtParam tool, MEME server, SOPMA, PDBsum tool, CYS-REC tool, SWISS-MODEL, SAVES servers, TMHMM program, GlobPlot, and peptide cutter software. The assessed in-silico data were compared with the published experimental results for validation.

Results

The alignment of sixty fibrinolytic serine protease enzymes (molecular mass 12-86 kDa) sequences showed 49 enzymes possess a conserved domain with a catalytic triad of Asp196, His242, and Ser569. The predicted instability and aliphatic indexes were 1.94-37.77, and 68.9-93.41, respectively, indicating high thermostability. The random coil means value suggested the predominance of this secondary structure in these proteases. A set of 50 amino acid residues representing motif 3 signifies the Peptidase S8/S53 domain that was invariably observed in 56 sequences. Additionally, 28 sequences have transmembrane helices, with two having the most disordered areas, and they pose 25 enzyme cleavage sites. A comparative analysis of the experimental work with the results of in-silico study put forward the characteristics of the enzyme sequences JF739176.1 and MF677779.1 to be considered when creating a potential mutant enzyme as these sequences are stable at high pH with thermostability and to exhibit αβ-fibrinogenase activity in both experimental and in-silico studies.

Release of CM-12 from A2-type casein by the cleavage of Ser-Leu-Xaa at the C-terminus using Aspergillus oryzae alkaline protease

BACKGROUND

Aspergillus oryzae protease can release the opioid peptide β-casomorphin-10 (CM-10, YPFPGPIPNS, 60-69) from A2-type casein. However, not only is the yield of the active peptide low, but the key enzyme involved in processing has yet to be identified.

RESULTS

A significant amount of the opioid peptide 60YPFPGPIPNSLP71 (CM-12) was produced from the A2-type casein peptide 53AQTQSLVYPFPGPIPNSLPQNIPPLTQTPV82 when the active protease in A. oryzae protease extract was fractionated with DEAE-Sepharose. The fractionated enzyme produced CM-12 from bovine A2-type casein but not from bovine A1 casein. A major protein of 34 kDa was purified and identified as an alkaline protease (Alp). Motif prediction of the Alp cleavage site using Multiple EM for Motif Elicitation analysis revealed preferable cleavage at the C-terminal end of Ser-Leu-Xaa for the release of CM-12. A2-type casein hydrolysate by Alp exhibited similar levels of opioid activity to that of synthetic CM-12 in cAMP-Glo assays with μ-opioid receptor-expressing HEK293 cells. These results suggest that CM-12 is a major opioid peptide in the casein hydrolysate.

CONCLUSION

Our findings showed that Alp fractionated from A. oryzae protease extract produced the opioid peptide CM-12 from A2-type casein as a result of preferential cleavage at the C-terminal end of Ser-Leu-Xaa and the removal of coexisting enzymes. Moreover, docking predictions suggested a stable interaction between CM-12 and the 3D structure of Alp. Casein hydrolysate with Alp-containing CM-12 has the potential for use as a bioactive peptide material with opioid activity. © 2024 Society of Chemical Industry.

Powerful cell wall biomass degradation enzymatic system from saprotrophic Aspergillus fumigatus

Cell wall biomass, Earth's most abundant natural resource, holds significant potential for sustainable biofuel production. Composed of cellulose, hemicellulose, lignin, pectin, and other polymers, the plant cell wall provides essential structural support to diverse organisms in nature. In contrast, non-plant species like insects, crustaceans, and fungi rely on chitin as their primary structural polysaccharide. The saprophytic fungus Aspergillus fumigatus has been widely recognized for its adaptability to various environmental conditions. It achieves this by secreting different cell wall biomass degradation enzymes to obtain essential nutrients. This review compiles a comprehensive collection of cell wall degradation enzymes derived from A. fumigatus, including cellulases, hemicellulases, various chitin degradation enzymes, and other polymer degradation enzymes. Notably, these enzymes exhibit biochemical characteristics such as temperature tolerance or acid adaptability, indicating their potential applications across a spectrum of industries.

In-silico characterization of GABAT protein found in gut-brain axis associated bacteria of healthy individuals and multiple sclerosis patients

Background

Multiple sclerosis (MS) is a neurodegenerative disease characterized by inflammation and demyelination of neurons. There is evidence to suggest that level of a neurotransmitter gamma-aminobutyric acid (GABA), due to the degradation by γ-aminobutyric acid transaminase (GABAT), is reduced in certain areas of the brain in MS patients. MS is always accompanied by gut bacteria dysbiosis. In healthy individuals, Faecalibacterium sp. while in MS patients A. calcoaceticus, Clostridium sp. and S. typhimurium are found abundantly. Although all these microbes produce GABAT but only in MS patients this enzyme significantly degrades GABA.

Objective

Present study is an attempt to characterize the GABAT protein sequences of these bacteria.

Methodology

Sequences of GABAT protein were retrieved from Uniprot database. Sequences were analyzed by Protparam, Gneg-mPLoc, SOSUI, PFP-FunDSeqE, Pepwheel program, PROTEUS and Alphafold and SAVES servers, MEME suite and HDOCK server.

Results

In healthy individuals gastrointestinal tract (GIT) bacteria, GABAT protein was present in inner-membrane with α helix content (61 and 62%) and β sheet content (5%), 4-helical cytokines functional domains. It has greater number of B-cell epitopes and more complex 3D configuration as compared to MS patients GIT bacterial enzymes.

Conclusion

Present study might enable us to modify the GABAT encoding gene and enzyme through site-directed mutagenesis in pathogenic bacteria thus reducing their potential of causing MS.

Exploring the association of ESR1 and ESR2 gene SNPs with polycystic ovary syndrome in human females: a comprehensive association study

BACKGROUND

Polycystic Ovary Syndrome (PCOS) affects a significant proportion of human females worldwide and is characterized by hormonal, metabolic, and reproductive dysfunctions, including infertility, irregular menstrual cycles, acanthosis nigricans, and hirsutism. Mutations in the estrogen receptor genes ESR1 and ESR2, involved in normal follicular development and ovulation, can contribute to development of the PCOS. The present study focuses on investigating the potential correlation between single nucleotide polymorphisms (SNPs) of ESR1 and ESR2 genes and the incidence of this syndrome.

METHODS

For this study, SNPs in ESR1 and ESR2 genes were retrieved from the ENSEMBL database and analyzed for their effect on mutated proteins using different bioinformatics tools including SIFT, PolyPhen, CADD, REVEL, MetaLR, I-Mutant, CELLO2GO, ProtParam, SOPMA, SWISS-MODEL and HDDOCK.

RESULTS

All the SNPs documented in the present study were deleterious. All the SNPs except rs1583384537, rs1450198518, and rs78255744 decreased protein stability. Two variants rs1463893698 and rs766843910 in the ESR2 gene altered the localization of mutated proteins i.e. in addition to the nucleus, proteins were also found in mitochondria and extracellular, respectively. SNPs rs104893956 in ESR1 and rs140630557, rs140630557, rs1596423459, rs766843910, rs1596405923, rs762454979 and rs1384121511 in ESR2 gene significantly changed the secondary structure of proteins (2D). SNPs that markedly changed 3D configuration included rs1554259481, rs188957694 and rs755667747 in ESR1 gene and rs1463893698, rs140630557, rs1596423459, rs766843910, rs1596405923, rs762454979 and rs1384121511 in ESR2 gene. Variants rs1467954450 (ESR1) and rs140630557 (ESR2) were identified to reduce the binding tendency of ESRα and β receptors with estradiol as reflected by the docking scores i.e. -164.97 and -173.23, respectively.

CONCLUSION

Due to the significant impact on the encoded proteins, these variants might be proposed as biomarkers to predict the likelihood of developing PCOS in the future and for diagnostic purposes.

In-Silico Characterization of Estrogen Reactivating β-Glucuronidase Enzyme in GIT Associated Microbiota of Normal Human and Breast Cancer Patients

Estrogen circulating in blood has been proved to be a strong biomarker for breast cancer. A β-glucuronidase enzyme (GUS) from human gastrointestinal tract (GIT) microbiota including probiotics has significant involvement in enhancing the estrogen concentration in blood through deconjugation of glucuronidated estrogens. The present project has been designed to explore GIT microbiome-encoded GUS enzymes (GUSOME) repertoire in normal human and breast cancer patients. For this purpose, a total of nineteen GUS enzymes from human GIT microbes, i.e., seven from healthy and twelve from breast cancer patients have been focused on. Protein sequences of enzymes retrieved from UniProt database were subjected to ProtParam, CELLO2GO, SOPMA (secondary structure prediction method), PDBsum (Protein Database summaries), PHYRE2 (Protein Homology/AnalogY Recognition Engine), SAVES v6.0 (Structure Validation Server), MEME version 5.4.1 (Multiple Em for Motif Elicitation), Caver Web server v 1.1, Interproscan and Predicted Antigenic Peptides tool. Analysis revealed the number of amino acids, isoelectric point, extinction coefficient, instability index and aliphatic index of GUS enzymes in the range of 586−795, 4.91−8.92, 89,980−155,075, 25.88−40.93 and 71.01−88.10, respectively. Sub-cellular localization of enzyme was restricted to cytoplasm and inner-membrane in case of breast cancer patients’ bacteria as compared to periplasmic space, outer membrane and extracellular space in normal GIT bacteria. The 2-D structure analysis showed α helix, extended strand, β turn and random coil in the range of 27.42−22.66%, 22.04−25.91%, 5.39−8.30% and 41.75−47.70%, respectively. The druggability score was found to be 0.05−0.45 and 0.06−0.80 in normal and breast cancer patients GIT, respectively. The radius, length and curvature of catalytic sites were observed to be 1.1−2.8 Å, 1.4−15.9 Å and 0.65−1.4, respectively. Ten conserved protein motifs with p < 0.05 and width 25−50 were found. Antigenic propensity-associated sequences were 20−29. Present study findings hint about the use of the bacterial GUS enzymes against breast cancer tumors after modifications via site-directed mutagenesis of catalytic sites involved in the activation of estrogens and through destabilization of these enzymes.

In vitro and in silico characterization of alkaline serine protease from Bacillus subtilis D9 recovered from Saudi Arabia

In this study, we have isolated and characterized proteolytic soil bacteria and their alkaline protease. Based on 16S rRNA sequence analysis, 12 isolates with the highest protease activity were classified as B. subtilis and B. cereus groups. B. subtilis D9 isolate showing the highest protease activity was selected for in vitro and in silico analysis for its ِِAKD9 protease. The enzyme has a molecular mass of 48 kDa, exhibiting optimal activity at 50 °C pH 9.5, and showed high stability till 65 °C and pH 8–11 for 1 h. Fe³⁺‏ stimulated, but Zn²⁺ and Hg²⁺ strongly inhibited the protease activity. Also, the maximum inhibition with PMSF indicated serine protease-type of AKD9 protease. AkD9 alkaline serine protease gene showed high sequence similarity and close phylogenetic relationship with AprX serine protease of B. subtilis isolates. Functional prediction of AKD9 resulted in the detection of subtilase domain, peptidase_S8 family, and subtilase active sites. Moreover, prediction of physicochemical properties indicated that AKD9 serine protease is hydrophilic, thermostable, and alkali-halo stable. Secondary structure prediction revealed the dominance of the coils enhances AKD9 activity and stability under saline and alkaline conditions. Based on molecular docking, AKD9 showed very promising binding affinities towards casein substrate with expected intrinsic proteolytic activities matching our obtained in vitro results. In conclusion, AKD9 alkaline serine protease seems to be a significant candidate for industrial applications because of its stability, hydrophilicity, enhanced thermostability, and alkali-halo stability.

Gene Analysis, Cloning, and Heterologous Expression of Protease from a Micromycete Aspergillus ochraceus Capable of Activating Protein C of Blood Plasma

Micromycetes are known to secrete numerous enzymes of biotechnological and medical potential. Fibrinolytic protease-activator of protein C (PAPC) of blood plasma from micromycete Aspergillus ochraceus VKM-F4104D was obtained in recombinant form utilising the bacterial expression system. This enzyme, which belongs to the proteinase-K-like proteases, is similar to the proteases encoded in the genomes of Aspergillus fumigatus ATCC MYA-4609, A. oryzae ATCC 42149 and A. flavus 28. Mature PAPC-4104 is 282 amino acids long, preceded by the 101-amino acid propeptide necessary for proper folding and maturation. The recombinant protease was identical to the native enzyme from micromycete in terms of its biological properties, including an ability to hydrolyse substrates of activated protein C (pGlu-Pro-Arg-pNA) and factor Xa (Z-D-Arg-Gly-Arg-pNA) in conjugant reactions with human blood plasma. Therefore, recombinant PAPC-4104 can potentially be used in medicine, veterinary science, diagnostics, and other applications.

In silico characterization and structural modeling of bacterial metalloprotease of family M4

BACKGROUND

The M4 family of metalloproteases is comprised of a large number of zinc-containing metalloproteases. A large number of these enzymes are important virulence factors of pathogenic bacteria and therefore potential drug targets. Whereas some enzymes have potential for biotechnological applications, the M4 family of metalloproteases is known almost exclusively from bacteria. The aim of the study was to identify the structure and properties of M4 metalloprotease proteins.

RESULTS

A total of 31 protein sequences of M4 metalloprotease retrieved from UniProt representing different species of bacteria have been characterized for various physiochemical properties. They were thermostable, hydrophillic protein of a molecular mass ranging from 38 to 66 KDa. Correlation on the basis of both enzymes and respective genes has also been studied by phylogenetic tree. B. cereus M4 metalloprotease (PDB ID: 1NPC) was selected as a representative species for secondary and tertiary structures among the M4 metalloprotease proteins. The secondary structure displaying 11 helices (H1-H11) is involved in 15 helix-helix interactions, while 4 β-sheet motifs composed of 15 β-strands in PDBsum. Possible disulfide bridges were absent in most of the cases. The tertiary structure of B. cereus M4 metalloprotease was validated by QMEAN4 and SAVES server (Ramachandran plot, verify 3D, and ERRAT) which proved the stability, reliability, and consistency of the tertiary structure of the protein. Functional analysis was done in terms of membrane protein topology, disease-causing region prediction, proteolytic cleavage sites prediction, and network generation. Transmembrane helix prediction showed absence of transmembrane helix in protein. Protein-protein interaction networks demonstrated that bacillolysin of B. cereus interacted with ten other proteins in a high confidence score. Five disorder regions were identified. Active sites analysis showed the zinc-binding residues-His-143, His-147, and Glu-167, with Glu-144 acting as the catalytic residues.

CONCLUSION

Moreover, this theoretical overview will help researchers to get a details idea about the protein structure and it may also help to design enzymes with desirable characteristics for exploiting them at industrial level or potential drug targets.

Bioinformatic mapping of a more precise Aspergillus niger degradome

Aspergillus niger has the ability to produce a large variety of proteases, which are of particular importance for protein digestion, intracellular protein turnover, cell signaling, flavour development, extracellular matrix remodeling and microbial defense. However, the A. niger degradome (the full repertoire of peptidases encoded by the A. niger genome) available is not accurate and comprehensive. Herein, we have utilized annotations of A. niger proteases in AspGD, JGI, and version 12.2 MEROPS database to compile an index of at least 232 putative proteases that are distributed into the 71 families/subfamilies and 26 clans of the 6 known catalytic classes, which represents ~ 1.64% of the 14,165 putative A. niger protein content. The composition of the A. niger degradome comprises ~ 7.3% aspartic, ~ 2.2% glutamic, ~ 6.0% threonine, ~ 17.7% cysteine, ~ 31.0% serine, and ~ 35.8% metallopeptidases. One hundred and two proteases have been reassigned into the above six classes, while the active sites and/or metal-binding residues of 110 proteases were recharacterized. The probable physiological functions and active site architectures of these peptidases were also investigated. This work provides a more precise overview of the complete degradome of A. niger, which will no doubt constitute a valuable resource and starting point for further experimental studies on the biochemical characterization and physiological roles of these proteases.

Computational Elucidation of Phylogenetic, Functional and Structural Features of Methioninase from Pseudomonas, Escherichia, Clostridium and Citrobacter Strains

BACKGROUND

L-Methioninase (EC 4.4.1.11; MGL) is a pyridoxal phosphate (PLP)-dependent enzyme that is produced by a variety of bacteria, fungi, and plants. L-methioninase, especially from Pseudomonas and Citrobacter sp., is considered as the efficient therapeutic enzyme, particularly in cancers such as glioblastomas, medulloblastoma, and neuroblastoma that are more sensitive to methionine starvation.

OBJECTIVE

The low stability is one of the main drawbacks of the enzyme; in this regard, in the current study, different features of the enzyme, including phylogenetic, functional, and structural from Pseudomonas, Escherichia, Clostridium, and Citrobacter strains were evaluated to find the best bacterial L-Methioninase.

METHODS

After the initial screening of L-Methioninase sequences from the above-mentioned bacterial strains, the three-dimensional structures of enzymes from Escherichia fergusonii, Pseudomonas fluorescens, and Clostridium homopropionicum were determined through homology modeling via GalaxyTBM server and refined by GalaxyRefine server.

RESULTS AND CONCLUSION

Afterwards, PROCHECK, verify 3D, and ERRAT servers were used for verification of the obtained models. Moreover, antigenicity, allergenicity, and physico-chemical analysis of enzymes were also carried out. In order to get insight into the interaction of the enzyme with other proteins, the STRING server was used. The secondary structure of the enzyme is mainly composed of random coils and alpha-helices. However, these outcomes should further be validated by wet-lab investigations.

In Silico Study and Optimization of Bacillus megaterium alpha-Amylases Production Obtained from Honey Sources

This study aimed to screen alpha-amylase producing microorganisms from honey as a low water activity medium, a suitable source for selecting stable and cost-beneficial bacterial enzyme production systems. Plackett-Burman method was used to select twelve effective factors including pH, inoculum size, temperature, time, corn starch, KH₂PO₄, peptone, MgSO₄, CaCl₂, NaCl, glycerin, and yeast extract concentrations on bacterial alpha-amylases production yield. The Box-Behnken method was utilized to optimize the level of selected significant factors. The stability of bacterial alpha-amylases was also determined in low pH and high-temperature conditions. In addition, in silico study was used to create the alpha-amylase structure and study the stability in high-temperature and low water available condition. Among all isolated and characterized microorganisms, Bacillus megaterium produced the highest amount of alpha-amylases. The in silico data showed the enzyme 3D structure similarity to alpha-amylase from Halothermothrix orenii and highly negative charge amino acids on its surface caused the enzyme activity and stability in low water conditions. Based on Box-Behnken results, the temperature 35 °C, pH 6 and starch 40 g/l were determined as the optimum level of significant factors to achieve the highest alpha-amylases unit (101.44 U/ml). This bacterial alpha-amylases enzyme showed stability at pH 5 and a range of temperatures from 40 to 60 °C that indicates this enzyme may possess the potential for using in industrial processes.

Computational Analysis of Therapeutic Enzyme Uricase from Different Source Organisms

Background:

Hyperuricemia and gout are the conditions, which is a response of accumulation of uric acid in the blood and urine. Uric acid is the product of purine metabolic pathway in humans. Uricase is a therapeutic enzyme that can enzymatically reduces the concentration of uric acid in serum and urine into more a soluble allantoin. Uricases are widely available in several sources like bacteria, fungi, yeast, plants and animals.

Objective:

The present study is aimed at elucidating the structure and physiochemical properties of uricase by insilico analysis.

Methods:

A total number of sixty amino acid sequences of uricase belongs to different sources were obtained from NCBI and different analysis like Multiple Sequence Alignment (MSA), homology search, phylogenetic relation, motif search, domain architecture and physiochemical properties including pI, EC, Ai, Ii, and were performed.

Results:

Multiple sequence alignment of all the selected protein sequences has exhibited distinct difference between bacterial, fungal, plant and animal sources based on the position-specific existence of conserved amino acid residues. The maximum homology of all the selected protein sequences is between 51-388. In singular category, homology is between 16-337 for bacterial uricase, 14-339 for fungal uricase, 12-317 for plants uricase, and 37-361 for animals uricase. The phylogenetic tree constructed based on the amino acid sequences disclosed clusters indicating that uricase is from different source. The physiochemical features revealed that the uricase amino acid residues are in between 300- 338 with a molecular weight as 33-39kDa and theoretical pI ranging from 4.95-8.88. The amino acid composition results showed that valine amino acid has a high average frequency of 8.79 percentage compared to different amino acids in all analyzed species.

Conclusion:

In the area of bioinformatics field, this work might be informative and a stepping-stone to other researchers to get an idea about the physicochemical features, evolutionary history and structural motifs of uricase that can be widely used in biotechnological and pharmaceutical industries. Therefore, the proposed in silico analysis can be considered for protein engineering work, as well as for gout therapy.

Analysis of the effect of temperature on protein abundance in Demodex-associated Bacillus oleronius

A potential role for bacteria in the induction of rosacea has been suggested. The aim of this work was to characterise the effect of temperature on the production of immunostimulatory proteins by Bacillus oleronius-a bacterium to which rosacea patients show sera reactivity and which was originally isolated from a Demodex mite from a rosacea patient. The affected skin of rosacea patients is at a higher temperature than unaffected skin, and it was postulated that this might alter the protein expression pattern of B. oleronius. B. oleronius growth was reduced at 37°C compared to 30°C but resulted in increased expression of the immune-reactive 62kDa protein (1.65 fold [P < 0.05]). Proteomic analysis revealed increased abundance of a wide range of proteins involved in the stress response (e.g. stress proteins [21.7-fold increase], phosphocarrier protein HPr [438.5-fold increase], 60 kDa chaperonin [12.6-fold increase]). Proteins decreased in abundance after growth at 37°C included ferredoxin (325-fold decrease) and peptidase (244-fold decrease). This work indicates that the increased skin temperature of rosacea patients may alter the growth and protein production pattern of B. oleronius and lead to the greater production of immuo-stimulatory proteins.

In vitro and in silico characterization of a novel dextranase from Pochonia chlamydosporia

The objective of this study was to purify, characterize, and phylogenetically and structurally analyze the dextranase produced by the fungus Pochonia chlamydosporia. Dextranase produced by the fungus P. chlamydosporia was purified to homogeneity in two steps, with a yield of 152%, purification factor of 6.84 and specific activity of 358.63 U/mg. Its molecular weight was estimated by SDS-PAGE at 64 kDa. The enzyme presented higher activity at 50 °C and pH 5.0, using 100 mM citrate-phosphate buffer, was inhibited by Ag1+, Hg2+, Cu2+, Mg2+, and presented KM of 23.60 µM. Mature dextranase is composed of 585 amino acids residues, with a predicted molecular weight of 64.38 kDa and pI 5.96. This dextranase showed a strong phylogenetic similarity when compared to Trichoderma harzianum dextranase. Its structure consists of two domains: the first composed by 15 β strands, and the second composed by a right-handed parallel β-helix.

Active Site Characterization of Proteases Sequences from Different Species of Aspergillus

A total of 129 proteases sequences comprising 43 serine proteases, 36 aspartic proteases, 24 cysteine protease, 21 metalloproteases, and 05 neutral proteases from different Aspergillus species were analyzed for the catalytically active site residues using MEROPS database and various bioinformatics tools. Different proteases have predominance of variable active site residues. In case of 24 cysteine proteases of Aspergilli, the predominant active site residues observed were Gln193, Cys199, His364, Asn384 while for 43 serine proteases, the active site residues namely Asp164, His193, Asn284, Ser349 and Asp325, His357, Asn454, Ser519 were frequently observed. The analysis of 21 metalloproteases of Aspergilli revealed Glu298 and Glu388, Tyr476 as predominant active site residues. In general, Aspergilli species-specific active site residues were observed for different types of protease sequences analyzed. The phylogenetic analysis of these 129 proteases sequences revealed 14 different clans representing different types of proteases with diverse active site residues.

Spatial differentiation of gene expression in Aspergillus niger colony grown for sugar beet pulp utilization

Degradation of plant biomass to fermentable sugars is of critical importance for the use of plant materials for biofuels. Filamentous fungi are ubiquitous organisms and major plant biomass degraders. Single colonies of some fungal species can colonize massive areas as large as five soccer stadia. During growth, the mycelium encounters heterogeneous carbon sources. Here we assessed whether substrate heterogeneity is a major determinant of spatial gene expression in colonies of Aspergillus niger. We analyzed whole-genome gene expression in five concentric zones of 5-day-old colonies utilizing sugar beet pulp as a complex carbon source. Growth, protein production and secretion occurred throughout the colony. Genes involved in carbon catabolism were expressed uniformly from the centre to the periphery whereas genes encoding plant biomass degrading enzymes and nitrate utilization were expressed differentially across the colony. A combined adaptive response of carbon-catabolism and enzyme production to locally available monosaccharides was observed. Finally, our results demonstrate that A. niger employs different enzymatic tools to adapt its metabolism as it colonizes complex environments.

Ecogenomics of Geminivirus from India and Neighbor Countries: An In Silico Analysis of Recombination Phenomenon

Recombination is one of the keys factor in evolutionary processes, involved in shaping the architecture of genomes and consequent phenotype. Understanding the recombination phenomenon especially among viruses will help in disease management. The present study aimed for in-silico analysis of recombination phenomenon among Begomoviruses, particularly emphasizing on viruses strains reported from India and neighboring countries. A total of 956 virus sequences have been used in this study. The Tomato yellow leaf curl China viruses, namely [Formula: see text] and [Formula: see text] were identified with the highest number of recombination event (1273). However, the Mung bean yellow mosaic India virus [Formula: see text] was found to have 1170 recombination event. The phylogenic analysis among the highly recombinant sequences was carried to get an insight of the evolution among viral sequences in this class of plant viruses. The phylogenetic analysis revealed a pattern in diversity among these virus strains and a split tree analysis showed diversity in the range of 0.049128335-10.269852. This in silico analysis may pave way for a greater understanding of recombination phenomenon in geminiviruses and it might be helpful for strategic plant viral disease management.

Ecogenomics of Geminivirus from India and neighbor countries: An in silico analysis of recombination phenomenon

Recombination is one of the keys factor in evolutionary processes, involved in shaping the architecture of genomes and consequent phenotype. Understanding the recombination phenomenon especially among viruses will help in disease management. The present study aimed for in-silico analysis of recombination phenomenon among Begomoviruses. Particularly emphasizing on viruses strains reported from India and neighboring countries. A total of 956 virus sequences have been used in this study. The Tomato yellow leaf curl China viruses, namely gi|29825986|; gi|283468151|; gi|190559151| and gi|61652782| were identified with the highest number of recombination event (1273). However, the Mung bean yellow mosaic India virus (gi|66351988|) was found to have 1170 recombination event. The phylogenic analysis among the highly recombinant sequences was carried to get an insight of the evolution among viral sequences in this class of plant viruses. The phylogenetic analysis revealed a pattern in diversity among these virus strains and a split tree analysis showed diversity in the range of 0.049128335 to 10.269852. This in silico analysis may pave way for a greater understanding of recombination phenomenon in Ggeminiviruses and it might be helpful for strategic plant viral disease management.

Ecogenomics of Geminivirus from India and neighbor countries: An in silico analysis of recombination phenomenon

Recombination is one of the keys factor in evolutionary processes, involved in shaping the architecture of genomes and consequent phenotype. Understanding the recombination phenomenon especially among viruses will help in disease management. The present study aimed for in-silico analysis of recombination phenomenon among Begomoviruses. Particularly emphasizing on viruses strains reported from India and neighboring countries. A total of 956 virus sequences have been used in this study. The Tomato yellow leaf curl China viruses, namely gi|29825986|; gi|283468151|; gi|190559151| and gi|61652782| were identified with the highest number of recombination event (1273). However, the Mung bean yellow mosaic India virus (gi|66351988|) was found to have 1170 recombination event. The phylogenic analysis among the highly recombinant sequences was carried to get an insight of the evolution among viral sequences in this class of plant viruses. The phylogenetic analysis revealed a pattern in diversity among these virus strains and a split tree analysis showed diversity in the range of 0.049128335 to 10.269852. This in silico analysis may pave way for a greater understanding of recombination phenomenon in Ggeminiviruses and it might be helpful for strategic plant viral disease management.

A genomic survey of proteases in Aspergilli

BACKGROUND

Proteases can hydrolyze peptides in aqueous environments. This property has made proteases the most important industrial enzymes by taking up about 60% of the total enzyme market. Microorganisms are the main sources for industrial protease production due to their high yield and a wide range of biochemical properties. Several Aspergilli have the ability to produce a variety of proteases, but no comprehensive comparative study has been carried out on protease productivity in this genus so far.

RESULTS

We have performed a combined analysis of comparative genomics, proteomics and enzymology tests on seven Aspergillus species grown on wheat bran and sugar beet pulp. Putative proteases were identified by homology search and Pfam domains. These genes were then clusters based on orthology and extracellular proteases were identified by protein subcellular localization prediction. Proteomics was used to identify the secreted enzymes in the cultures, while protease essays with and without inhibitors were performed to determine the overall protease activity per protease class. All this data was then integrated to compare the protease productivities in Aspergilli.

CONCLUSIONS

Genomes of Aspergillus species contain a similar proportion of protease encoding genes. According to comparative genomics, proteomics and enzymatic experiments serine proteases make up the largest group in the protease spectrum across the species. In general wheat bran gives higher induction of proteases than sugar beet pulp. Interesting differences of protease activity, extracellular enzyme spectrum composition, protein occurrence and abundance were identified for species. By combining in silico and wet-lab experiments, we present the intriguing variety of protease productivity in Aspergilli.

Cloning, Sequencing, and In Silico Analysis of β-Propeller Phytase Bacillus licheniformis Strain PB-13

β-Propeller phytases (BPPhy) are widely distributed in nature and play a major role in phytate-phosphorus cycling. In the present study, a BPPhy gene from Bacillus licheniformis strain was expressed in E. coli with a phytase activity of 1.15 U/mL and specific activity of 0.92 U/mg proteins. The expressed enzyme represented a full length ORF “PhyPB13” of 381 amino acid residues and differs by 3 residues from the closest similar existing BPPhy sequences. The PhyPB13 sequence was characterized in silico using various bioinformatic tools to better understand structural, functional, and evolutionary aspects of BPPhy class by multiple sequence alignment and homology search, phylogenetic tree construction, variation in biochemical features, and distribution of motifs and superfamilies. In all sequences, conserved sites were observed toward their N-terminus and C-terminus. Cysteine was not present in the sequence. Overall, three major clusters were observed in phylogenetic tree with variation in biophysical characteristics. A total of 10 motifs were reported with motif “1” observed in all 44 protein sequences and might be used for diversity and expression analysis of BPPhy enzymes. This study revealed important sequence features of BPPhy and pave a way for determining catalytic mechanism and selection of phytase with desirable characteristics.

Purification and biochemical characterization of a novel alkaline protease produced by Penicillium nalgiovense

Penicillium nalgiovense PNA9 produces an extracellular protease during fermentation with characteristics of growth-associated product. Enzyme purification involved ammonium sulfate precipitation, dialysis, and ultrafiltration, resulting in 12.1-fold increase of specific activity (19.5 U/mg). The protein was isolated through a series of BN-PAGE and native PAGE runs. ESI-MS analysis confirmed the molecular mass of 45.2 kDa. N-Terminal sequencing (MGFLKLLKGSLATLAVVNAGKLLTANDGDE) revealed 93 % similarity to a Penicillium chrysogenum protease, identified as major allergen. The protease exhibits simple Michaelis-Menten kinetics and K m (1.152 mg/ml), V max (0.827 mg/ml/min), and k cat (3.2 × 10(2)) (1/s) values against azocasein show that it possesses high substrate affinity and catalytic efficiency. The protease is active within 10-45 °C, pH 4.0-10.0, and 0-3 M NaCl, while maximum activity was observed at 35 °C, pH 8.0, and 0.25 M NaCl. It is active against the muscle proteins actin and myosin and inactive against myoglobin. It is highly stable in the presence of non-ionic surfactants, hydrogen peroxide, BTNB, and EDTA. Activity was inhibited by SDS, Mn(2+) and Zn(2+), and by the serine protease inhibitor PMSF, indicating the serine protease nature of the enzyme. These properties make the novel protease a suitable candidate enzyme in meat ripening and other biotechnological applications.

Serine protease identification (in vitro) and molecular structure predictions (in silico) from a phytopathogenic fungus, Alternaria solani

Serine proteases are involved in an enormous number of biological processes. The present study aims at characterizing three-dimensional (3D) molecular architecture of serine proteases from early blight pathogen, Alternaria solani that are hypothesized to be markers of phytopathogenicity. A serine protease was purified to homogeneity and MALDI-TOF-MS/MS analysis revealed that protease produced by A. solani belongs to alkaline serine proteases (AsP). AsP is made up of 403 amino acid residues with molecular weight of 42.1 kDa (Isoelectric point - 6.51) and its molecular formula was C1859 H2930 N516 O595 S4 . AsP structure model was built based on its comparative homology with serine protease using the program, MODELER. AsP had 16 β-sheets and 10 α-helices, with Ser(350) (G347-G357), Asp(158) (D158-H169), and His(193) (H193-G203) in separate turn/coil structures. Biological metal binding region situated near 6th-helix and His(193) residue is responsible for metal binding site. Also, calcium ion (Ca(2+)) is coordinated by the carboxyl groups of Lys(84), Ile(85), Lys(86), Asp(87), Phe(88), Ala(89), Ala(90) (K84-A90) for first Ca(2+) binding site and carbonyl oxygen atom of Lys(244), Gly(245), Arg(246), Thr(247), Lys(248), Lys(249), and Ala(250) (K244-A250), for second Ca(2+) binding site. Moreover, Ramachandran plot analysis of protein residues falling into most favored secondary structures were determined (83.3%). The predicted molecular 3D structural model was further verified using PROCHECK, ERRAT, and VADAR servers to confirm the geometry and stereo-chemical parameters of the molecular structural design. The functional analysis of AsP 3D molecular structure predictions familiar in the current study may provide a new perspective in the understanding and identification of antifungal protease inhibitor designing.

X-ray structure analysis and characterization of AFUEI, an elastase inhibitor from Aspergillus fumigatus

Elastase from Aspergillus sp. is an important factor for aspergillosis. AFUEI is an inhibitor of the elastase derived from Aspergillus fumigatus. AFUEI is a member of the I78 inhibitor family and has a high inhibitory activity against elastases of Aspergillus fumigatus and Aspergillus flavus, human neutrophil elastase and bovine chymotrypsin, but does not inhibit bovine trypsin. Here we report the crystal structure of AFUEI in two crystal forms. AFUEI is a wedge-shaped protein composed of an extended loop and a scaffold protein core. The structure of AFUEI shows remarkable similarity to serine protease inhibitors of the potato inhibitor I family, although they are classified into different inhibitor families. A structural comparison with the potato I family inhibitors suggests that the extended loop of AFUEI corresponds to the binding loop of the potato inhibitor I family, and AFUEI inhibits its cognate proteases through the same mechanism as the potato I family inhibitors.

In silico characterization of histidine Acid phytase sequences

Histidine acid phytases (HAPhy) are widely distributed enzymes among bacteria, fungi, plants, and some animal tissues. They have a significant role as an animal feed enzyme and in the solubilization of insoluble phosphates and minerals present in the form of phytic acid complex. A set of 50 reference protein sequences representing HAPhy were retrieved from NCBI protein database and characterized for various biochemical properties, multiple sequence alignment (MSA), homology search, phylogenetic analysis, motifs, and superfamily search. MSA using MEGA5 revealed the presence of conserved sequences at N-terminal "RHGXRXP" and C-terminal "HD." Phylogenetic tree analysis indicates the presence of three clusters representing different HAPhy, that is, PhyA, PhyB, and AppA. Analysis of 10 commonly distributed motifs in the sequences indicates the presence of signature sequence for each class. Motif 1 "SPFCDLFTHEEWIQYDYLQSLGKYYGYGAGNPLGPAQGIGF" was present in 38 protein sequences representing clusters 1 (PhyA) and 2 (PhyB). Cluster 3 (AppA) contains motif 9 "KKGCPQSGQVAIIADVDERTRKTGEAFAAGLAPDCAITVHTQADTSSPDP" as a signature sequence. All sequences belong to histidine acid phosphatase family as resulted from superfamily search. No conserved sequence representing 3- or 6-phytase could be identified using multiple sequence alignment. This in silico analysis might contribute in the classification and future genetic engineering of this most diverse class of phytase.

Virtual screening and evaluation of Ketol-Acid Reducto-Isomerase (KARI) as a putative drug target for Aspergillosis

Aspergillus is a leading causative agent for fungal morbidity and mortality in immuno-compromised patients. To identify a putative target to design or identify new antifungal drug, against Aspergillus is required. In our previous work, we have analyzed the various biochemical pathways, and we found Ketol Acid Reducto-Isomerase (KARI) an enzyme involves in the amino acid biosynthesis, could be a better target. This enzyme was found to be unique by comparing to host proteome through BLASTp analysis. A homology based model of KARI was generated by Swiss model server. The generated model had been validated by PROCHECK and WHAT IF programs. The Zinc library was generated within the limitation of the Lipinski rule of five, for docking study. Based on the dock-score six molecules have been studied for ADME/TOX analysis and subjected for pharmacophore model generation. The Zinc ID of the potential inhibitors is ZINC00720614, ZINC01068126, ZINC0923, ZINC02090678, ZINC00663057 and ZINC02284065 and found to be pharmacologically active agonist and antagonist of KARI. This study is an attempt to Insilco evaluation of the KARI as a drug target and the screened inhibitors could help in the development of the better drug against Aspergillus.