Enzyme Commission Research Articles

Accurately predicting enzyme functions through geometric graph learning on ESMFold-predicted structures

Enzymes are crucial in numerous biological processes, with the Enzyme Commission (EC) number being a commonly used method for defining enzyme function. However, current EC number prediction technologies have not fully recognized the importance of enzyme active sites and structural characteristics. Here, we propose GraphEC, a geometric graph learning-based EC number predictor using the ESMFold-predicted structures and a pre-trained protein language model. Specifically, we first construct a model to predict the enzyme active sites, which is utilized to predict the EC number. The prediction is further improved through a label diffusion algorithm by incorporating homology information. In parallel, the optimum pH of enzymes is predicted to reflect the enzyme-catalyzed reactions. Experiments demonstrate the superior performance of our model in predicting active sites, EC numbers, and optimum pH compared to other state-of-the-art methods. Additional analysis reveals that GraphEC is capable of extracting functional information from protein structures, emphasizing the effectiveness of geometric graph learning. This technology can be used to identify unannotated enzyme functions, as well as to predict their active sites and optimum pH, with the potential to advance research in synthetic biology, genomics, and other fields.

GloEC: a hierarchical-aware global model for predicting enzyme function.

The annotation of enzyme function is a fundamental challenge in industrial biotechnology and pathologies. Numerous computational methods have been proposed to predict enzyme function by annotating enzyme labels with Enzyme Commission number. However, the existing methods face difficulties in modelling the hierarchical structure of enzyme label in a global view. Moreover, they haven't gone entirely to leverage the mutual interactions between different levels of enzyme label. In this paper, we formulate the hierarchy of enzyme label as a directed enzyme graph and propose a hierarchy-GCN (Graph Convolutional Network) encoder to globally model enzyme label dependency on the enzyme graph. Based on the enzyme hierarchy encoder, we develop an end-to-end hierarchical-aware global model named GloEC to predict enzyme function. GloEC learns hierarchical-aware enzyme label embeddings via the hierarchy-GCN encoder and conducts deductive fusion of label-aware enzyme features to predict enzyme labels. Meanwhile, our hierarchy-GCN encoder is designed to bidirectionally compute to investigate the enzyme label correlation information in both bottom-up and top-down manners, which has not been explored in enzyme function prediction. Comparative experiments on three benchmark datasets show that GloEC achieves better predictive performance as compared to the existing methods. The case studies also demonstrate that GloEC is capable of effectively predicting the function of isoenzyme. GloEC is available at: https://github.com/hyr0771/GloEC.

FURNA: A database for functional annotations of RNA structures.

Despite the increasing number of 3D RNA structures in the Protein Data Bank, the majority of experimental RNA structures lack thorough functional annotations. As the significance of the functional roles played by noncoding RNAs becomes increasingly apparent, comprehensive annotation of RNA function is becoming a pressing concern. In response to this need, we have developed FURNA (Functions of RNAs), the first database for experimental RNA structures that aims to provide a comprehensive repository of high-quality functional annotations. These include Gene Ontology terms, Enzyme Commission numbers, ligand-binding sites, RNA families, protein-binding motifs, and cross-references to related databases. FURNA is available at https://seq2fun.dcmb.med.umich.edu/furna/ to enable quick discovery of RNA functions from their structures and sequences.

A complementary metaproteomic approach to interrogate microbiome cultivated from clinical colon biopsies.

The human gut microbiome plays a vital role in preserving individual health and is intricately involved in essential functions. Imbalances or dysbiosis within the microbiome can significantly impact human health and are associated with many diseases. Several metaproteomics platforms are currently available to study microbial proteins within complex microbial communities. In this study, we attempted to develop an integrated pipeline to provide deeper insights into both the taxonomic and functional aspects of the cultivated human gut microbiomes derived from clinical colon biopsies. We combined a rapid peptide search by MSFragger against the Unified Human Gastrointestinal Protein database and the taxonomic and functional analyses with Unipept Desktop and MetaLab-MAG. Across seven samples, we identified and matched nearly 36,000 unique peptides to approximately 300 species and 11 phyla. Unipept Desktop provided gene ontology, InterPro entries, and enzyme commission number annotations, facilitating the identification of relevant metabolic pathways. MetaLab-MAG contributed functional annotations through Clusters of Orthologous Genes and Non-supervised Orthologous Groups categories. These results unveiled functional similarities and differences among the samples. This integrated pipeline holds the potential to provide deeper insights into the taxonomy and functions of the human gut microbiome for interrogating the intricate connections between microbiome balance and diseases.

CASTpFold: Computed Atlas of Surface Topography of the universe of protein Folds.

Geometric and topological properties of protein structures, including surface pockets, interior cavitiesand cross channels, are of fundamental importance for proteins to carry out their functions. Computed Atlas of Surface Topography of proteins(CASTp) is a widely used web server for locating, delineating, and measuring these geometric and topological properties of protein structures. Recent developments in AI-based protein structure prediction such as AlphaFold2 (AF2) have significantly expanded our knowledge on protein structures. Here we present CASTpFold, a continuation of CASTp that provides accurate and comprehensive identifications and quantifications of protein topography. It now provides (i) results on an expanded database of proteins, including the Protein Data Bank(PDB) and non-singleton representative structures of AlphaFold2 structures, covering 183 million AF2 structures; (ii) functional pockets prediction with corresponding Gene Ontology (GO) terms or Enzyme Commission (EC) numbers for AF2-predicted structuresand (iii) pocket similarity search function for surface and protein-protein interface pockets. The CASTpFold web server is freely accessible at https://cfold.bme.uic.edu/castpfold/.

IfDEEPre: large protein language-based deep learning enables interpretable and fast predictions of enzyme commission numbers

Abstract Accurate understanding of the biological functions of enzymes is vital for various tasks in both pathologies and industrial biotechnology. However, the existing methods are usually not fast enough and lack explanations on the prediction results, which severely limits their real-world applications. Following our previous work, DEEPre, we propose a new interpretable and fast version (ifDEEPre) by designing novel self-guided attention and incorporating biological knowledge learned via large protein language models to accurately predict the commission numbers of enzymes and confirm their functions. Novel self-guided attention is designed to optimize the unique contributions of representations, automatically detecting key protein motifs to provide meaningful interpretations. Representations learned from raw protein sequences are strictly screened to improve the running speed of the framework, 50 times faster than DEEPre while requiring 12.89 times smaller storage space. Large language modules are incorporated to learn physical properties from hundreds of millions of proteins, extending biological knowledge of the whole network. Extensive experiments indicate that ifDEEPre outperforms all the current methods, achieving more than 14.22% larger F1-score on the NEW dataset. Furthermore, the trained ifDEEPre models accurately capture multi-level protein biological patterns and infer evolutionary trends of enzymes by taking only raw sequences without label information. Meanwhile, ifDEEPre predicts the evolutionary relationships between different yeast sub-species, which are highly consistent with the ground truth. Case studies indicate that ifDEEPre can detect key amino acid motifs, which have important implications for designing novel enzymes. A web server running ifDEEPre is available at https://proj.cse.cuhk.edu.hk/aihlab/ifdeepre/ to provide convenient services to the public. Meanwhile, ifDEEPre is freely available on GitHub at https://github.com/ml4bio/ifDEEPre/.

REME: an integrated platform for reaction enzyme mining and evaluation.

A key challenge in pathway design is finding proper enzymes that can be engineered to catalyze a non-natural reaction. Although existing tools can identify potential enzymes based on similar reactions, these tools encounter several issues. Firstly, the calculated similar reactions may not even have the same reaction type. Secondly, the associated enzymes are often numerous and identifying the most promising candidate enzymes is difficult due to the lack of data for evaluation. Thirdly, existing web tools do not provide interactive functions that enable users to fine-tune results based on their expertise. Here, we present REME(https://reme.biodesign.ac.cn/), the first integrated web platform for reaction enzyme mining and evaluation. Combining atom-to-atom mapping, atom type change identification, and reaction similarity calculation enables quick ranking and visualization of reactions similar to an objective non-natural reaction. Additional functionality enables users to filter similar reactions by their specified functional groups and candidate enzymes can be further filtered (e.g.by organisms) or expanded by Enzyme Commission number (EC) or sequence homology. Afterward, enzyme attributes (such as kcat, Km, optimal temperature and pH) can be assessed with deep learning-based methods, facilitating the swift identification of potential enzymes that can catalyze the non-natural reaction.

Bulk and single cells transcriptomes with experimental validation identify USP18 as a novel glioma prognosis and proliferation indicator.

The mechanism by which ubiquitin-specific protease 18 (USP18) (enzyme commission: 3.4.19.12) inhibition in cancer promotes cell pyroptosis via the induction of interferon (IFN)-stimulated genes has been recently demonstrated. It is also known that USP18 influences the epithelial-mesenchymal transition of glioma cells. In the present study, the upregulation of USP18 in glioma was revealed through bulk transcriptome analysis, which was associated with poor prognosis in patients with glioma. Furthermore, USP18 levels affected the response to immunotherapy in patients with glioma. Single-cell transcriptome and enrichment analyses demonstrated that USP18 was associated with type 1 IFN responses in glioma T cells. To demonstrate the effect of USP18 expression levels on glioma cells, USP18 expression was knocked down in U251 and U87MG ATCC cell lines. A subsequent Cell Counting Kit-8 assay revealed that glioma cell viability was significantly decreased 4 days after USP18 knockdown. In addition, the knockdown of USP18 expression significantly inhibited the clonogenicity of U251 and U87MG ATCC cells. In conclusion, the present study demonstrated that knockdown of USP18 expression inhibited the proliferation of glioma cells, which may be mediated by the effect of USP18 on the IFN-I response.

Macroalgal microbiomes unveil a valuable genetic resource for halogen metabolism

BackgroundMacroalgae, especially reds (Rhodophyta Division) and browns (Phaeophyta Division), are known for producing various halogenated compounds. Yet, the reasons underlying their production and the fate of these metabolites remain largely unknown. Some theories suggest their potential antimicrobial activity and involvement in interactions between macroalgae and prokaryotes. However, detailed investigations are currently missing on how the genetic information of prokaryotic communities associated with macroalgae may influence the fate of organohalogenated molecules.ResultsTo address this challenge, we created a specialized dataset containing 161 enzymes, each with a complete enzyme commission number, known to be involved in halogen metabolism. This dataset served as a reference to annotate the corresponding genes encoded in both the metagenomic contigs and 98 metagenome-assembled genomes (MAGs) obtained from the microbiome of 2 red (Sphaerococcus coronopifolius and Asparagopsis taxiformis) and 1 brown (Halopteris scoparia) macroalgae. We detected many dehalogenation-related genes, particularly those with hydrolytic functions, suggesting their potential involvement in the degradation of a wide spectrum of halocarbons and haloaromatic molecules, including anthropogenic compounds. We uncovered an array of degradative gene functions within MAGs, spanning various bacterial orders such as Rhodobacterales, Rhizobiales, Caulobacterales, Geminicoccales, Sphingomonadales, Granulosicoccales, Microtrichales, and Pseudomonadales. Less abundant than degradative functions, we also uncovered genes associated with the biosynthesis of halogenated antimicrobial compounds and metabolites.ConclusionThe functional data provided here contribute to understanding the still largely unexplored role of unknown prokaryotes. These findings support the hypothesis that macroalgae function as holobionts, where the metabolism of halogenated compounds might play a role in symbiogenesis and act as a possible defense mechanism against environmental chemical stressors. Furthermore, bacterial groups, previously never connected with organohalogen metabolism, e.g., Caulobacterales, Geminicoccales, Granulosicoccales, and Microtrichales, functionally characterized through MAGs reconstruction, revealed a biotechnologically relevant gene content, useful in synthetic biology, and bioprospecting applications.7KFSwtia72NY8FA2HCcej3Video

PredictEFC: a fast and efficient multi-label classifier for predicting enzyme family classes

BackgroundEnzymes play an irreplaceable and important role in maintaining the lives of living organisms. The Enzyme Commission (EC) number of an enzyme indicates its essential functions. Correct identification of the first digit (family class) of the EC number for a given enzyme is a hot topic in the past twenty years. Several previous methods adopted functional domain composition to represent enzymes. However, it would lead to dimension disaster, thereby reducing the efficiency of the methods. On the other hand, most previous methods can only deal with enzymes belonging to one family class. In fact, several enzymes belong to two or more family classes.ResultsIn this study, a fast and efficient multi-label classifier, named PredictEFC, was designed. To construct this classifier, a novel feature extraction scheme was designed for processing functional domain information of enzymes, which counting the distribution of each functional domain entry across seven family classes in the training dataset. Based on this scheme, each training or test enzyme was encoded into a 7-dimenion vector by fusing its functional domain information and above statistical results. Random k-labelsets (RAKEL) was adopted to build the classifier, where random forest was selected as the base classification algorithm. The two tenfold cross-validation results on the training dataset shown that the accuracy of PredictEFC can reach 0.8493 and 0.8370. The independent test on two datasets indicated the accuracy values of 0.9118 and 0.8777.ConclusionThe performance of PredictEFC was slightly lower than the classifier directly using functional domain composition. However, its efficiency was sharply improved. The running time was less than one-tenth of the time of the classifier directly using functional domain composition. In additional, the utility of PredictEFC was superior to the classifiers using traditional dimensionality reduction methods and some previous methods, and this classifier can be transplanted for predicting enzyme family classes of other species. Finally, a web-server available at http://124.221.158.221/ was set up for easy usage.

Insights into the inner workings of transformer models for protein function prediction.

We explored how explainable artificial intelligence (XAI) can help to shed light into the inner workings of neural networks for protein function prediction, by extending the widely used XAI method of integrated gradients such that latent representations inside of transformer models, which were finetuned to Gene Ontology term and Enzyme Commission number prediction, can be inspected too. The approach enabled us to identify amino acids in the sequences that the transformers pay particular attention to, and to show that these relevant sequence parts reflect expectations from biology and chemistry, both in the embedding layer and inside of the model, where we identified transformer heads with a statistically significant correspondence of attribution maps with ground truth sequence annotations (e.g. transmembrane regions, active sites) across many proteins. Source code can be accessed at https://github.com/markuswenzel/xai-proteins.

Parallel convolutional contrastive learning method for enzyme function prediction.

The function labeling of enzymes has a wide range of application value in the medical field, industrial biology and other fields. Scientists define enzyme categories by enzyme commission (EC) numbers. At present, although there are some tools for enzyme function prediction, their effects have not reached the application level. To improve the precision of enzyme function prediction, we propose a parallel convolutional contrastive learning (PCCL) method to predict enzyme functions. First, we use the advanced protein language model ESM-2 to preprocess the protein sequences. Second, PCCL combines convolutional neural networks (CNNs) and contrastive learning to improve the prediction precision of multifunctional enzymes. Contrastive learning can make the model better deal with the problem of class imbalance. Finally, the deep learning framework is mainly composed of three parallel CNNs for fully extracting sample features. we compare PCCL with state-of-art enzyme function prediction methods based on three evaluation metrics. The performance of our model improves on both two test sets. Especially on the smaller test set, PCCL improves the AUC by 2.57%. The source code can be downloaded from https://github.com/biomg/PCCL.

Evidential deep learning for trustworthy prediction of enzyme commission number.

The rapid growth of uncharacterized enzymes and their functional diversity urge accurate and trustworthy computational functional annotation tools. However, current state-of-the-art models lack trustworthiness on the prediction of the multilabel classification problem with thousands of classes. Here, we demonstrate that a novel evidential deep learning model (named ECPICK) makes trustworthy predictions of enzyme commission (EC) numbers with data-driven domain-relevant evidence, which results in significantly enhanced predictive power and the capability to discover potential new motif sites. ECPICK learns complex sequential patterns of amino acids and their hierarchical structures from 20million enzyme data. ECPICK identifies significant amino acids that contribute to the prediction without multiple sequence alignment. Our intensive assessment showed not only outstanding enhancement of predictive performance on the largest databases of Uniprot, Protein Data Bank (PDB) and Kyoto Encyclopedia of Genes and Genomes (KEGG), but also a capability to discover new motif sites in microorganisms. ECPICK is a reliable EC number prediction tool to identify protein functions of an increasing number of uncharacterized enzymes.

Increased bacterial taxonomic and functional diversity is associated with impaired rotavirus vaccine immunogenicity in infants from India and Malawi

The immunogenicity and effectiveness of oral rotavirus vaccines (ORVs) against severe rotavirus-associated gastroenteritis are impaired in low- and middle-income countries (LMICs) where the burden of disease is highest. Determining risk factors for impaired ORV response may help identify strategies to enhance vaccine effectiveness. In this study, we use metagenomic sequencing to provide a high-resolution taxonomic analysis of stool samples collected at 6 weeks of age (coinciding with the first ORV dose) during a prospective study of ORV immunogenicity in India and Malawi. We then analyse the functional capacity of the developing microbiome in these cohorts. Microbiome composition differed significantly between countries, although functional capacity was more similar than taxonomic composition. Our results confirm previously reported findings that the developing microbiome is more diverse in taxonomic composition in ORV non-seroconverters compared with seroconverters, and we additionally demonstrate a similar pattern in functional capacity. Although taxonomic or functional feature abundances are poor predictors of ORV response, we show that skews in the direction of associations within these microbiome data can be used to identify consistent markers of ORV response across LMIC infant cohorts. We also highlight the systemic under-representation of reference genes from LMICs that limit functional annotation in our study (7% and 13% annotation at pathway and enzyme commission level, respectively). Overall, higher microbiome diversity in early life may act as marker for impaired ORV response in India and Malawi, whilst a holistic perspective of functional capacity may be hidden in the “dark matter” of the microbiome.

Functional annotation of enzyme-encoding genes using deep learning with transformer layers

Functional annotation of open reading frames in microbial genomes remains substantially incomplete. Enzymes constitute the most prevalent functional gene class in microbial genomes and can be described by their specific catalytic functions using the Enzyme Commission (EC) number. Consequently, the ability to predict EC numbers could substantially reduce the number of un-annotated genes. Here we present a deep learning model, DeepECtransformer, which utilizes transformer layers as a neural network architecture to predict EC numbers. Using the extensively studied Escherichia coli K-12 MG1655 genome, DeepECtransformer predicted EC numbers for 464 un-annotated genes. We experimentally validated the enzymatic activities predicted for three proteins (YgfF, YciO, and YjdM). Further examination of the neural network’s reasoning process revealed that the trained neural network relies on functional motifs of enzymes to predict EC numbers. Thus, DeepECtransformer is a method that facilitates the functional annotation of uncharacterized genes.

RefSeq and the prokaryotic genome annotation pipeline in the age of metagenomes.

The Reference Sequence (RefSeq) project at the National Center for Biotechnology Information (NCBI) contains over 315000 bacterial and archaeal genomes and 236 million proteins with up-to-date and consistent annotation. In the past 3 years, we have expanded the diversity of the RefSeq collection by including the best quality metagenome-assembled genomes (MAGs) submitted to INSDC (DDBJ, ENA and GenBank), while maintaining its quality by adding validation checks. Assemblies are now more stringently evaluated for contamination and for completeness of annotation prior to acceptance into RefSeq. MAGs now account for over 17000 assemblies in RefSeq, split over 165 orders and 362 families. Changes in the Prokaryotic Genome Annotation Pipeline (PGAP), which is used to annotate nearly all RefSeq assemblies include better detection of protein-coding genes. Nearly 83% of RefSeq proteins are now named by a curated Protein Family Model, a 4.7% increase in the past three years ago. In addition to literature citations, Enzyme Commission numbers, and gene symbols, Gene Ontology terms are now assigned to 48% of RefSeq proteins, allowing for easier multi-genome comparison. RefSeq is found at https://www.ncbi.nlm.nih.gov/refseq/. PGAP is available as a stand-alone tool able to produce GenBank-ready files at https://github.com/ncbi/pgap.

"Revealing the genetic arsenal of Bacillus firmus TNAU1: Unleashing nematicidal and plant growth promotion traits"

Bacillus firmus is a promising nematicidal bacterium that also confers plant growth promotion, nutrient acquisition and antimicrobial activities. In the present study, attempts were made for whole genome sequencing of B. firmus TNAU1 isolated from the rhizosphere region of cucumber. The results confirmed that the isolate was B. firmus which contained a single circular genome of 5,370,919 bp (5.3 Mb) in length with an average GC content of 41.16 %. Annotation was performed using PATRIC tool to explore the subsystem analysis, specialty genes, Enzyme Commission (EC) numbers, Gene Ontology (GO) mappings, KEGG pathways of B. firmus TNAU1. A total of 5804 coding sequence (CDS) genes that encode proteins were predicted based on the KEGG and clusters of orthologous groups (COGs) database. Comparative genomic analysis was performed to discover the unique features and their relationship with other available B. firmus strains. Blastp analysis revealed that TNAU 1 genome assemblies had 26 homologs corresponding to nematode-virulent proteases and various genes responsible for plant growth promotion which includes indole acetic acid (IAA) production, nitrate transport, phosphorus and potassium solubilization. Four secondary metabolite biosynthetic clusters encoding antibiotics and siderophore production were confirmed through comparative genomics. B. firmus TNAU1 was evaluated for its antagonistic activity against root-knot nematode Meloidogyne incognita of tomato. The results exhibited that B. firmus inhibited the egg hatching and mortality of juveniles by 96% over control under in vitro conditions and greenhouse experiments confirmed that the TNAU1 strain not only had the potential to suppress the soil and root population of M. incognita (70% reduction over control), but also triggered the plant growth parameters of tomato.

ProtEC: A Transformer Based Deep Learning System for Accurate Annotation of Enzyme Commission Numbers.

The advancements in next-generation sequencing technologies have given rise to large-scale, open-source protein databases consisting of hundreds of millions of sequences. However, to make these sequences useful in biomedical applications, they need to be painstakingly annotated by curating them from literature. To counteract this problem, many automated annotation algorithms have been developed over the years including deep learning models, especially in recent times. In this work, we propose a transformer-based deep-learning model that can predict the Enzyme Commission numbers of an enzyme from full-scale sequences with state-of-the-art accuracy compared to other recent machine learning annotation algorithms. The system does especially well on clustered split dataset which consists of training and testing samples derived from different distributions that are structurally dissimilar from each other. This proves that the model is able to understand deep patterns within the sequences and can accurately identify the motifs responsible for the different enzyme commission numbers. Moreover, the algorithm is able to retain similar accuracy even when the training size is significantly reduced, and also, the model accuracy is independent of the sequence length making it suitable for a wide range of applications consisting of varying sequence structures.

Predicting enzymatic function of protein sequences with attention.

There is a growing number of available protein sequences, but only a limited amount has been manually annotated. For example, only 0.25% of all entries of UniProtKB are reviewed by human annotators. Further developing automatic tools to infer protein function from sequence alone can alleviate part of this gap. In this article, we investigate the potential of Transformer deep neural networks on a specific case of functional sequence annotation: the prediction of enzymatic classes. We show that our EnzBert transformer models, trained to predict Enzyme Commission (EC) numbers by specialization of a protein language model, outperforms state-of-the-art tools for monofunctional enzyme class prediction based on sequences only. Accuracy is improved from 84% to 95% on the prediction of EC numbers at level two on the EC40 benchmark. To evaluate the prediction quality at level four, the most detailed level of EC numbers, we built two new time-based benchmarks for comparison with state-of-the-art methods ECPred and DeepEC: the macro-F1 score is respectively improved from 41% to 54% and from 20% to 26%. Finally, we also show that using a simple combination of attention maps is on par with, or better than, other classical interpretability methods on the EC prediction task. More specifically, important residues identified by attention maps tend to correspond to known catalytic sites. Quantitatively, we report a max F-Gain score of 96.05%, while classical interpretability methods reach 91.44% at best. Source code and datasets are respectively available at https://gitlab.inria.fr/nbuton/tfpc and https://doi.org/10.5281/zenodo.7253910.

Measuring β-galactosidase activity in opaque dairy solutions under optimum conditions for galactooligosaccharide synthesis by isothermal titration calorimetry

The dairy industry uses enzymes to make cheese, alter product flavor, and eliminate lactose. The activities of these enzymes have been measured in clear buffered solutions, but because of the limitations of spectrophotometric methods, enzyme activities have not been measured in opaque or colored dairy products where they are used. Isothermal titration calorimetry (ITC) can be used to determine reaction kinetics in opaque and colored solutions by measuring the heat rate (thermal power) from enzyme-catalyzed reactions as a function of time. This study used ITC to measure β-galactosidase activity in opaque solutions of milk, sweet whey, sweet whey permeate, acid whey, and acid whey permeate with 2 β-galactosidase (Enzyme Commission number 3.2.1.23) isozymes derived from Aspergillus oryzae and Kluyveromyces lactis. The components of the dairy fluids alter the enzyme kinetics and reaction thermodynamics, and the reactions catalyzed by the 2 homologues differ as shown by differing thermodynamic profiles. The study demonstrates that ITC can be used to measure enzyme activity in opaque and colored dairy fluids and identify reactions by their thermodynamic properties.