β-glucosidases from Aspergillus unguis NII 08123: Molecular characterization, properties and applications

Lignocellulosic biomass is probably the best alternative resource for biofuel production and it is composed mainly of cellulose, hemicelluloses and lignin. Cellulose is the most abundant among the three and conversion of cellulose to glucose is catalyzed by the enzyme cellulase. Cellulases are groups of enzymes act synergistically upon cellulose to produce glucose and comprise of endoglucanase, cellobiohydrolase and β-glucosidase. β -glucosidase assumes great importance due to the fact that it is the rate limiting enzyme. Endoglucanases (EG) produces nicks in the cellulose polymer exposing reducing and non reducing ends, cellobiohydrolases (CBH) acts upon the reducing or non reducing ends to liberate cellobiose units, and β - glucosidases (BGL) cleaves the cellobiose to liberate glucose completing the hydrolysis. . β -glucosidases undergo feedback inhibition by their own product- β glucose, and cellobiose which is their substrate. Few filamentous fungi produce glucose tolerant β - glucosidases which can overcome this inhibition by tolerating the product concentration to a particular threshold. The present study had targeted a filamentous fungus producing glucose tolerant β - glucosidase which was identified by morphological as well as molecular method. The fungus showed 99% similarity to Aspergillus unguis strain which comes under the Aspergillus nidulans group where most of the glucose tolerant β -glucosidase belongs. The culture was designated the strain number NII 08123 and was deposited in the NII culture collection at CSIR-NIIST. β -glucosidase multiplicity is a common occurrence in fungal world and in A.unguis this was demonstrated using zymogram analysis. A total 5 extracellular isoforms were detected in fungus and the expression levels of these five isoforms varied based on the carbon source available in the medium. Three of these 5 isoforms were expressed in higher levels as identified by the increased fluorescence (due to larger amounts of MUG breakdown by enzyme action) and was speculated to contribute significantly to the total _- β glucosidase activity. These isoforms were named as BGL 1, BGL3 and BGL 5. Among the three, BGL5 was demonstrated to be the glucose tolerant β -glucosidase and this was a low molecular weight protein. Major fraction was a high molecular weight protein but with lesser tolerance to glucose. BGL 3 was between the two in both activity and glucose tolerance.121 Glucose tolerant .β -glucosidase was purified and characterized and kinetic analysis showed that the glucose inhibition constant (Ki) of the protein is 800mM and Km and Vmax of the enzyme was found to be 4.854 mM and 2.946 mol min-1mg protein-1respectively. The optimumtemperature was 60°C and pH 6.0. The molecular weight of the purified protein was ~10kDa in both SDS as well as Native PAGE indicating that the glucose tolerant BGL is a monomeric protein.The major β -glucosidase, BGL1 had a pH and temperature optima of 5.0 and 60 °C respectively. The apparent molecular weight of the Native protein is 240kDa. The Vmax and Km was 78.8 mol min-1mg protein-1 and 0.326mM respectively. Degenerate primers were designed for glycosyl hydrolase families 1, 3 and 5 and the BGL genes were amplified from genomic DNA of Aspergillus unguis. The sequence analyses performed on the amplicons results confirmed the presence of all the three genes. Amplicon with a size of ~500bp was sequenced and which matched to a GH1 –BGL from Aspergillus oryzae. GH3 degenerate primers producing amplicons were sequenced and the sequences matched to β - glucosidase of GH3 family from Aspergillus nidulans and Aspergillus acculateus. GH5 degenerate primers also gave amplification and sequencing results indicated the presence of GH5 family BGL gene in the Aspergillus unguis genomic DNA.From the partial gene sequencing results, specific as well as degenerate primers were designed for TAIL PCR. Sequencing results of the 1.0 Kb amplicon matched Aspergillus nidulans β -glucosidase gene which belongs to the GH1 family. The sequence mainly covered the N-Terminal region of the matching peptide. All the three BGL proteins ie. BGL1, BGL3 and BGL5 were purified by chromatography an electro elution from Native PAGE gels and were subjected to MALDI-TOF mass spectrometric analysis. The results showed that BGL1 peptide mass matched to . β -glucosidase-I of Aspergillus flavus which is a 92kDa protein with 69% protein coverage. The glucose tolerant β -glucosidase BGL5 mass matched to the catalytic C-terminal domain of β -glucosidase-F from Emericella nidulans, but the protein coverage was very low compared to the size of the Emericella nidulans protein. While comparing the size of BGL5 from Aspergillus unguis, the protein sequence coverage is more than 80%. BGL F is a glycosyl hydrolase family 3 protein.The properties of BGL5 seem to be very unique, in that it is a GH3 β -glucosidase with a very low molecular weight of ~10kDa and at the same time having catalytic activity and glucose 122 tolerance which is as yet un-described in GH β -glucosidases. The occurrence of a fully functional 10kDA protein with glucose tolerant BGL activity has tremendous implications both from the points of understanding the structure function relationships as well as for applications of BGL enzymes. BGL-3 showed similarity to BGL1 of Aspergillus aculateus which was another GH3 β -glucosidase. It may be noted that though PCR could detect GH1, GH3 and GH5 β-glucosidases in the fungus, the major isoforms BGL1 BGL3 and BGL5 were all GH3 family enzymes. This would imply that β-glucosidases belonging to other families may also co-exist in the fungus and the other minor isoforms detected in zymograms may account for them. In biomass hydrolysis, GT-BGL containing BGL enzyme was supplemented to cellulase and the performances of blends were compared with a cocktail where commercial β- glucosidase was supplemented to the biomass hydrolyzing enzyme preparation. The cocktail supplemented with A unguis BGL preparation yielded 555mg/g sugar in 12h compared to the commercial enzyme preparation which gave only 333mg/g in the same period and the maximum sugar yield of 858 mg/g was attained in 36h by the cocktail containing A. unguis BGL. While the commercial enzyme achieved almost similar sugar yield in 24h, there was rapid drop in sugar concentration after that, indicating probably the conversion of glucose back to di-or oligosaccharides by the transglycosylation activity of the BGl in that preparation. Compared this, the A.unguis enzyme containing preparation supported peak yields for longer duration (upto 48h) which is important for biomass conversion to other products since the hydrolysate has to undergo certain unit operations before it goes into the next stage ie – fermentation in any bioprocesses for production of either fuels or chemicals.. Most importantly the Aspergillus unguis BGL preparation yields approximately 1.6 fold increase in the sugar release compared to the commercial BGL within 12h of time interval and 2.25 fold increase in the sugar release compared to the control ie. Cellulase without BGL supplementation. The current study therefore leads to the identification of a potent new isolate producing glucose tolerant β - glucosidase. The organism identified as Aspergillus unguis comes under the Aspergillus nidulans group where most of the GT-BGL producers belong and the detailed studies showed that the glucose tolerant β -glucosidase was a very low molecular weight protein which probably belongs to the glycosyl hydrolase family 3. Inhibition kinetic studies helped to understand the Ki and it is the second highest among the nidulans group of Aspergilli. This has promoted us for a detailed study regarding the mechanism of glucose tolerance. The proteomic 123 analyses clearly indicate the presence of GH3 catalytic domain in the protein. Since the size of the protein is very low and still its active and showed glucose tolerance it is speculated that this could be an entirely new protein or the modification of the existing β -glucosidase with only the catalytic domain present in it. Hydrolysis experiments also qualify this BGL, a suitable candidate for the enzyme cocktail development for biomass hydrolysis

Practical applications of quantum molecular similarity measures (QMSM): programs and examples

Es descriu l'aproximació de Capes Atòmiques dins de la teoria de la Semblança Molecular Quàntica. Partint només de dades teòriques, s'ha trobat una relació entre estructura molecular i activitat biològica per a diversos conjunts de molècules. Es descriuen els aspectes teòrics de la Semblança Molecular Quàntica i alguns exemples d'aplicació

Phylogeny, function, and evolution of the cupins, a structurally conserved, functionally diverse superfamily of proteins

The cupin superfamily is a group of functionally diverse proteins that are found in all three kingdoms of life, Archaea, Eubacteria, and Eukaryota. These proteins have a characteristic signature domain comprising two histidine- containing motifs separated by an intermotif region of variable length. This domain consists of six beta strands within a conserved beta barrel structure. Most cupins, such as microbial phosphomannose isomerases (PMIs), AraC- type transcriptional regulators, and cereal oxalate oxidases (OXOs), contain only a single domain, whereas others, such as seed storage proteins and oxalate decarboxylases (OXDCs), are bi-cupins with two pairs of motifs. Although some cupins have known functions and have been characterized at the biochemical level, the majority are known only from gene cloning or sequencing projects. In this study, phylogenetic analyses were conducted on the conserved domain to investigate the evolution and structure/function relationships of cupins, with an emphasis on single- domain plant germin-like proteins (GLPs). An unrooted phylogeny of cupins from a wide spectrum of evolutionary lineages identified three main clusters, microbial PMIs, OXDCs, and plant GLPs. The sister group to the plant GLPs in the global analysis was then used to root a phylogeny of all available plant GLPs. The resulting phylogeny contained three main clades, classifying the GLPs into distinct subfamilies. It is suggested that these subfamilies correlate with functional categories, one of which contains the bifunctional barley germin that has both OXO and superoxide dismutase (SOD) activity. It is proposed that GLPs function primarily as SODs, enzymes that protect plants from the effects of oxidative stress. Closer inspection of the DNA sequence encoding the intermotif region in plant GLPs showed global conservation of thymine in the second codon position, a character associated with hydrophobic residues. Since many of these proteins are multimeric and enzymatically inactive in their monomeric state, this conservation of hydrophobicity is thought to be associated with the need to maintain the various monomer- monomer interactions. The type of structure-based predictive analysis presented in this paper is an important approach for understanding gene function and evolution in an era when genomes from a wide range of organisms are being sequenced at a rapid rate.

Microbial relatives of the seed storage proteins of higher plants: conservation of structure and diversification of function during evolution of the cupin superfamily.

This review summarizes the recent discovery of the cupin superfamily (from the Latin term "cupa," a small barrel) of functionally diverse proteins that initially were limited to several higher plant proteins such as seed storage proteins, germin (an oxalate oxidase), germin-like proteins, and auxin-binding protein. Knowledge of the three-dimensional structure of two vicilins, seed proteins with a characteristic beta-barrel core, led to the identification of a small number of conserved residues and thence to the discovery of several microbial proteins which share these key amino acids. In particular, there is a highly conserved pattern of two histidine-containing motifs with a varied intermotif spacing. This cupin signature is found as a central component of many microbial proteins including certain types of phosphomannose isomerase, polyketide synthase, epimerase, and dioxygenase. In addition, the signature has been identified within the N-terminal effector domain in a subgroup of bacterial AraC transcription factors. As well as these single-domain cupins, this survey has identified other classes of two-domain bicupins including bacterial gentisate 1, 2-dioxygenases and 1-hydroxy-2-naphthoate dioxygenases, fungal oxalate decarboxylases, and legume sucrose-binding proteins. Cupin evolution is discussed from the perspective of the structure-function relationships, using data from the genomes of several prokaryotes, especially Bacillus subtilis. Many of these functions involve aspects of sugar metabolism and cell wall synthesis and are concerned with responses to abiotic stress such as heat, desiccation, or starvation. Particular emphasis is also given to the oxalate-degrading enzymes from microbes, their biological significance, and their value in a range of medical and other applications.

A quantum chemical study on a set of non-imidazole H(3) antihistamine molecules

Molecular orbital calculations were carried out on a set of 28 non-imidazole H(3) antihistamine compounds using the Hartree-Fock method in order to investigate the possible relationships between electronic structural properties and binding affinity for H3 receptors (pK(i)). It was observed that the frontier effective-for-reaction molecular orbital (FERMO) energies were better correlated with pK(i) values than highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy values. Exploratory data analysis through hierarchical cluster (HCA) and principal component analysis (PCA) showed a separation of the compounds in two sets, one grouping the molecules with high pK(i) values, the other gathering low pK(i) value compounds. This separation was obtained with the use of the following descriptors: FERMO energies (epsilon(FERMO)), charges derived from the electrostatic potential on the nitrogen atom (N(1)), electronic density indexes for FERMO on the N(1) atom (Sigma((FERMO))c(i)(2)). and electrophilicity (omega`). These electronic descriptors were used to construct a quantitative structure-activity relationship (QSAR) model through the partial least-squares (PLS) method with three principal components. This model generated Q(2) = 0.88 and R(2) = 0.927 values obtained from a training set and external validation of 23 and 5 molecules, respectively. After the analysis of the PLS regression equation and the values for the selected electronic descriptors, it is suggested that high values of FERMO energies and of Sigma((FERMO))c(i)(2), together with low values of electrophilicity and pronounced negative charges on N(1) appear as desirable properties for the conception of new molecules which might have high binding affinity. 2010 Elsevier Inc. All rights reserved.

In Vitro Phytotoxicity and Antioxidant Activity of Selected Flavonoids

The knowledge of flavonoids involved in plant-plant interactions and their mechanisms of action are poor and, moreover, the structural characteristics required for these biological activities are scarcely known. The objective of this work was to study the possible in vitro phytotoxic effects of 27 flavonoids on the germination and early radical growth of Raphanus sativus L. and Lepidium sativum L., with the aim to evaluate the possible structure/activity relationship. Moreover, the antioxidant activity of the same compounds was also evaluated. Generally, in response to various tested flavonoids, germination was only slightly affected, whereas significant differences were observed in the activity of the various tested flavonoids against radical elongation. DPPH test confirms the antioxidant activity of luteolin, quercetin, catechol, morin, and catechin. The biological activity recorded is discussed in relation to the structure of compounds and their capability to interact with cell structures and physiology. No correlation was found between phytotoxic and antioxidant activities.

Analysis of Bothrops jararacussu venomous gland transcriptome focusing on structural and functional aspects: I - gene expression profile of highly expressed phospholipases A(2)

Snake venom glands are a rich source of bioactive molecules such as peptides, proteins and enzymes that show important pharmacological activity leading to in local and systemic effects as pain, edema, bleeding and muscle necrosis. Most studies on pharmacologically active peptides and proteins from snake venoms have been concerned with isolation and structure elucidation through methods of classical biochemistry. As an attempt to examine the transcripts expressed in the venom gland of Bothrops jararacussu and to unveil the toxicological and pharmacological potential of its products at the molecular level, we generated 549 expressed sequence tags (ESTs) from a directional cDNA library. Sequences obtained from single-pass sequencing of randomly selected cDNA clones could be identified by similarities searches on existing databases, resulting in 197 sequences with significant similarity to phospholipase A(2) (PLA(2)), of which 83.2% were Lys49-PLA(2) homologs (BOJU-1), 0.1% were basic Asp49-PLA(2)s (BOJU-II) and 0.6% were acidic Asp49-PLA(2)s (BOJU-III). Adjoining this very abundant class of proteins we found 88 transcripts codifying for putative sequences of metalloproteases, which after clustering and assembling resulted in three full-length sequences: BOJUMET-I, BOJUMET-II and BOJUMET-III; as well as 25 transcripts related to C-type lectin like protein including a full-length cDNA of a putative galactose binding C-type lectin and a cluster of eight serine-proteases transcripts including a full-length cDNA of a putative serine protease. Among the full-length sequenced clones we identified a nerve growth factor (Bj-NGF) with 92% identity with a human NGF (NGHUBM) and an acidic phospholipase A2 (BthA-I-PLA(2)) displaying 85-93% identity with other snake venom toxins. Genetic distance among PLA(2)s from Bothrops species were evaluated by phylogenetic analysis. Furthermore, analysis of full-length putative Lys49-PLA(2) through molecular modeling showed conserved structural domains, allowing the characterization of those proteins as group II PLA(2)s. The constructed cDNA library provides molecular clones harboring sequences that can be used to probe directly the genetic material from gland venom of other snake species. Expression of complete cDNAs or their modified derivatives will be useful for elucidation of the structure-function relationships of these toxins and peptides of biotechnological interest. (C) 2004 Elsevier SAS. All rights reserved.

At the interface: Crystal structures of phospholipases A(2)

The protein content of many snake venoms often includes one or more phospholipases A(2) (PLA(2)). In recent years a growing number of venoms from snakes of Agkistrodon, Bothrops and Trimeresurus species have been shown to contain a catalytically inactive PLA(2)-homologue in which the highly conserved aspartic acid at position 49 (Asp49) is substituted by lysine (Lys49). Although demonstrating little or no catalytic activity, these Lys49-PLA(2)s disrupt membranes by a Ca2+-independent mechanism of action. In addition, this family of PLA(2)s demonstrates myotoxic and cytolytic pharmacological activities, however the structural bases underlying these functional properties are poorly understood. Through the application of X-ray crystallography in combination with biophysical and bioinformatics techniques, we are studying structure/function relationships of Lys49-PLA(2)s. We here present results of a systematic X-ray crystallographic and amino acid sequence analysis study of Lys49-PLA(2)s and propose a model to explain the Ca2+ independent membrane damaging activity. (C) 1998 Elsevier B.V. Ltd. All rights reserved.

A partial least squares and principal component regression study of quinone compounds with trypanocidal activity

A quantitative structure-activity relationship (QSAR) study of 19 quinone compounds with trypanocidal activity was performed by Partial Least Squares (PLS) and Principal Component Regression (PCR) methods with the use of leave-one-out crossvalidation procedure to build the regression models. The trypanocidal activity of the compounds is related to their first cathodic potential (Ep(c1)). The regression PLS and PCR models built in this study were also used to predict the Ep(c1) of six new quinone compounds. The PLS model was built with three principal components that described 96.50% of the total variance and present Q(2) = 0.83 and R-2 = 0.90. The results obtained with the PCR model were similar to those obtained with the PLS model. The PCR model was also built with three principal components that described 96.67% of the total variance with Q(2) = 0.83 and R-2 = 0.90. The most important descriptors for our PLS and PCR models were HOMO-1 (energy of the molecular orbital below HOMO), Q4 (atomic charge at position 4), MAXDN (maximal electrotopological negative difference), and HYF (hydrophilicity index).

Novel allosteric conformation of human HB revealed by the hydration and anion effects on O-2 binding

The effect of anions on the stability of different functional conformations of Hb is examined through the determination of the dependence of O-2 affinity on water activity (a(w)). The control of a(w) is effected by varying the sucrose osmolal concentration in the bathing solution according to the osmotic stress method. Thus, the hydration change following Hb oxygenation is determined as a function of Cl- and of DPG concentration. We find that only similar to 25 additional water molecules bind to human Hb during the deoxy-to-oxy conformation transition in the absence of anions, in contrast with similar to 72 that bind in the presence of more than 50 mM Cl- or more than 15 mu M DPG. We demonstrate that the increase in the hydration change linked with oxygenation is coupled with anion binding to the deoxy-Hb. Hence, we propose that the deoxy-Hb coexists in two allosteric conformations which depend on whether anion is bound or not: the tense T-state, with low oxygen affinity and anion bound, or a new allosteric P-state, with intermediate oxygen affinity and free of bound anions. The intrinsic oxygen affinity of this unforeseen P-state and the differential binding of Cl-, DPG, and H2O between states P and T and P and R are characteristics which are consistent with those expected for a putative intermediate allosteric state of Hb. These findings represent a new opportunity to explore the structure-function relationships of hemoglobin regulation.

SEQUENCE OF A CDNA-ENCODING BOTHROPSTOXIN-I, A MYOTOXIN FROM THE VENOM OF BOTHROPS-JARARACUSSU

With the aim of further understanding the structure/function relationships in the membrane-damaging activity of the Lys(49) phospholipase A(2) (Lys(49)-PLA(2)) sub-family, we used PCR (polymerase chain reaction) on total venom gland cDNAs from Bothrops jararacussu with degenerate oligodeoxyribonucleotides encoding the N- and C-termini of myotoxin II, a Lys(49)-PLA(2) from Bothrops asper. A 350-bp cDNA coding for bothropstoxin I (BtxtxI) was amplified. Sequencing of the amplified fragment shows that BtxtxI has a Lys(49), and comparison with the known structure of myotoxin II showed that the amino acids involved in the formation of a novel dimeric structure in this protein were also conserved.

Functional and structural analysis of two fibrinogen-activating enzymes isolated from the venoms of Crotalus durissus terrificus and Crotalus durissus collilineatus

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Structure of thrombin complexed with selective non-electrophilic inhibitors having cyclohexyl moieties at P1

The crystal structures of five new non-electrophilic β-strand-templated thrombin active-site inhibitors have been determined bound to the enzyme. Four co-crystallize with hirugen and inhibitor isomorphously to produce thrombin-hirugen crystals (monoclinic, space group C2), while one co-crystallizes in the hexagonal system, space group P65. A 1,4-substituted cyclohexyl moiety is conserved at the P1 position of all the inhibitors, along with a fused hetero-bicyclic five- and six-membered ring that occupies the P2 site. Amino, amidino and aminoimidazole groups are attached to the cyclohexyl ring for recognition at the S1 specificity site, while benzylsulfonyl and diphenyl groups enhance the binding at the S3 subsite. The cyclohexyl groups at the P1 positions of three of the inhibitors appear to be in the energetically favored chair conformation, while the imidazole-substituted cyclohexyl rings are in a boat conformation. Somewhat unexpectedly, the two cyclohexyl-aminoimidazole groups bind differently in the specificity site; the unique binding of one is heretofore unreported. The other inhibitors generally mimic arginyl binding at S1. This group of inhibitors combines the nonelectrophilicity and selectivity of DAPA-like compounds and the more optimal binding features of the S1-S3 sites of thrombin for peptidic molecules, which results in highly potent (binding constants 12 nM-16 pM, one being 1.1 μM) and selective (ranging from 140 to 20 000 times more selective compared with trypsin) inhibitors of thrombin. The binding modes of these novel inhibitors are correlated with their binding constants, as is their selectivity, in order to provide further insight for the design of therapeutic antithrombotic agents that inhibit thrombin directly at the active site.

Conformational basis for the biological activity of TOAC-labeled angiotensin II and bradykinin: Electron paramagnetic resonance, circular dichroism, and fluorescence studies

N-Terminally and internally labeled analogues of the hormones angiotensin (AII, DRVYIHPF) and bradykinin (BK, RPPGFSPFR) were synthesized containing the paramagnetic amino acid 2,2,6,6-tetramethylpiperidine-1-oxyl-4-amino-4- carboxylic acid (TOAC). TOAC replaced Asp 1 (TOAC 1-AII) and Val 3 (TOAC 3-AII) in AII and was inserted prior to Arg 1 (TOAC 0-BK) and replacing Pro 3 (TOAC 3-BK) in BK. The peptide conformational properties were examined as a function of trifluoroethanol (TFE) content and pH. Electron paramagnetic resonance spectra were sensitive to both variables and showed that internally labeled analogues yielded rotational correlation times (TC) considerably larger than N-terminally labeled ones, evincing the greater freedom of motion of the N-terminus. In TFE, τ C increased due to viscosity effects. Calculation of τ Cpeptide/τ CTOAC ratios indicated that the peptides acquired more folded conformations. Circular dichroism spectra showed that, except for TOAC 1-AII in TFE, the N-terminally labeled analogues displayed a conformational behavior similar to that of the parent peptides. In contrast, under all conditions, the TOAC 3 derivatives acquired more restricted conformations. Fluorescence spectra of All and its derivatives were especially sensitive to the ionization of Tyr 4. Fluorescence quenching by the nitroxide moiety was much more pronounced for TOAC 3-AII The conformational behavior of the TOAC derivatives bears excellent correlation with their biological activity, since, while the N-terminally labeled peptides were partially active, their internally labeled counterparts were inactive [Nakaie, C. R., et al., Peptides 2002, 23, 65-70]. The data demonstrate that insertion of TOAC in the middle of the peptide chain induces conformational restrictions that lead to loss of backbone flexibility, not allowing the peptides to acquire their receptor-bound conformation. © 2004 Wiley Periodicals, Inc.

Conformation and lytic activity of eumenine mastoparan: A new antimicrobial peptide from wasp venom

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)