Molecular docking and molecular dynamics simulations studies on β-glucosidase and xylanase Trichoderma asperellum to predict degradation order of cellulosic components in oil palm leaves for nanocellulose preparation

Ammonia monooxygenase-mediated cometabolic biotransformation of volatile 4-chlorophenol in nitrifying counter-diffused biofilms: A combined molecular dynamics simulation, DFT calculation and experimental study

Ammonia monooxygenase (AMO)-mediated cometabolism of organic pollutants has been widely observed in biological nitrogen removal process. However, its molecular mechanism remains unclear, hindering its practical application. Furthermore, conventional nitrification systems encounter significant challenges such as air pollution and the loss of ammonia-oxidizing bacteria, when dealing with wastewater containing volatile organic pollutants. This study developed a nitrifying membrane-aerated biofilm reactor (MABR) to enhance the biodegradation of volatile 4-chlorophenol (4-CP). Results showed that 4-CP was primarily removed via Nitrosomonas nitrosa-mediated cometabolism in the presence of NH4+-N, supported by the increased nicotinamide adenine dinucleotide (NADH) and adenosine triphosphate (ATP) content, AMO activity and the related genes abundance. Hydroquinone, detected for the first time and produced via oxidative dechlorination, as well as 4-chlorocatechol was primary transformation products of 4-CP. Nitrosomonas nitrosa AMO structural model was constructed for the first time using homology modeling. Molecular dynamics simulation suggested that the ortho-carbon in the benzene ring of 4-CP was more prone to metabolismcompared to the ipso-carbon. Density functional theory calculation revealed that 4-CP was metabolized by AMO via H-abstraction-OH-rebound reaction, with a significantly higher rebound barrier at the ipso-carbon (16.37 kcal·mol-1) as compared to the ortho-carbon (6.7 kcal·mol-1). This study fills the knowledge gap on the molecular mechanism of AMO-mediated cometabolism of organic pollutants, providing practical and theoretical foundations for improving volatile organic pollutants removal through nitrifying MABR.

A Versatile β-Glycosidase from Petroclostridium xylanilyticum Prefers the Conversion of Ginsenoside Rb3 over Rb1, Rb2, and Rc to Rd by Its Specific Cleavage Activity toward 1,6-Glycosidic Linkages

To convert ginsenosides Rb1, Rb2, Rb3, and Rc into Rd by a single enzyme, a putative β-glycosidase (Pxbgl) from the xylan-degrading bacterium Petroclostridium xylanilyticum was identified and used. The kcat/Km value of Pxbgl for Rb3 was 18.18 ± 0.07 mM-1/s, which was significantly higher than those of Pxbgl for other ginsenosides. Pxbgl converted almost all Rb3 to Rd with a productivity of 5884 μM/h, which was 346-fold higher than that of only β-xylosidase from Thermoascus aurantiacus. The productivity of Rd from the Panax ginseng root and Panax notoginseng leaf was 146 and 995 μM/h, respectively. Mutants N293 K and I447L from site-directed mutagenesis based on bioinformatics analysis showed an increase in specific activity of 29 and 7% toward Rb3, respectively. This is the first report of a β-glycosidase that can simultaneously remove four different glycosyls at the C-20 position of natural PPD-type ginsenosides and produce Rd as the sole product from P. notoginseng leaf extracts with the highest productivity.

Trends in the development and current perspective of thermostable bacterial hemicellulases with their industrial endeavors: A review

Hemicellulases are enzymes that hydrolyze hemicelluloses, common polysaccharides in nature. Thermophilic hemicellulases, derived from microbial strains, are extensively studied as natural biofuel sources due to the complex structure of hemicelluloses. Recent research aims to elucidate the catalytic principles, mechanisms and specificity of hemicellulases through investigations into their high-temperature stability and structural features, which have applications in biotechnology and industry. This review article targets to serve as a comprehensive resource, highlighting the significant progress in the field and emphasizing the vital role of thermophilic hemicellulases in eco-friendly catalysis. The primary goal is to improve the reliability of hemicellulase enzymes obtained from thermophilic bacterial strains. Additionally, with their ability to break down lignocellulosic materials, hemicellulases hold immense potential for biofuel production. Despite their potential, the commercial viability is hindered by their high enzyme costs, necessitating the development of efficient bioprocesses involving waste pretreatment with microbial consortia to overcome this challenge.

Comparative modeling and enzymatic affinity of novel haloacid dehalogenase from Bacillus megaterium strain BHS1 isolated from alkaline Blue Lake in Turkey

This study presents the initial structural model of L-haloacid dehalogenase (DehLBHS1) from Bacillus megaterium BHS1, an alkalotolerant bacterium known for its ability to degrade halogenated environmental pollutants. The model provides insights into the structural features of DehLBHS1 and expands our understanding of the enzymatic mechanisms involved in the degradation of these hazardous pollutants. Key amino acid residues (Arg40, Phe59, Asn118, Asn176, and Trp178) in DehLBHS1 were identified to play critical roles in catalysis and molecular recognition of haloalkanoic acid, essential for efficient binding and transformation of haloalkanoic acid molecules. DehLBHS1 was modeled using I-TASSER, yielding a best TM-score of 0.986 and an RMSD of 0.53 Å. Validation of the model using PROCHECK revealed that 89.2% of the residues were located in the most favored region, providing confidence in its structural accuracy. Molecular docking simulations showed that the non-simulated DehLBHS1 preferred 2,2DCP over other substrates, forming one hydrogen bond with Arg40 and exhibiting a minimum energy of -2.5 kJ/mol. The simulated DehLBHS1 exhibited a minimum energy of -4.3 kJ/mol and formed four hydrogen bonds with Arg40, Asn176, Asp9, and Tyr11, further confirming the preference for 2,2DCP. Molecular dynamics simulations supported this preference, based on various metrics, including RMSD, RMSF, gyration, hydrogen bonding, and molecular distance. MM-PBSA calculations showed that the DehLBHS1-2,2-DCP complex had a markedly lower binding energy (-21.363 ± 1.26 kcal/mol) than the DehLBHS1-3CP complex (-14.327 ± 1.738 kcal/mol). This finding has important implications for the substrate specificity and catalytic function of DehLBHS1, particularly in the bioremediation of 2,2-DCP in contaminated alkaline environments. These results provide a detailed view of the molecular interactions between the enzyme and its substrate and may aid in the development of more efficient biocatalytic strategies for the degradation of halogenated compounds.Communicated by Ramaswamy H. Sarma.

Natural sapogenins as potential inhibitors of aquaporins for targeted cancer therapy: computational insights into binding and inhibition mechanism

Aquaporins (AQPs) are membrane proteins that facilitate the transport of water and other small molecules across biological membranes. AQPs are involved in various physiological processes and pathological conditions, including cancer, making them as potential targets for anticancer therapy. However, the development of selective and effective inhibitors of AQPs remains a challenge. In this study, we explored the possibility of using natural sapogenins, a class of plant-derived aglycones of saponins with diverse biological activities, as potential inhibitors of AQPs. We performed molecular docking, dynamics simulation and binding energy calculation to investigate the binding and inhibition mechanism of 19 sapogenins against 13 AQPs (AQP0-AQP13) that are overexpressed in various cancers. Our results showed that out of 19 sapogenins, 8 (Diosgenin, Gitogenin, Tigogenin, Ruscogenin, Yamogenin, Hecogenin, Sarsasapogenin and Smilagenin) exhibited acceptable drug-like characteristics. These sapogenin also exhibited favourable binding affinities in the range of -7.6 to -13.4 kcal/mol, and interactions within the AQP binding sites. Furthermore, MD simulations provided insights into stability and dynamics of the sapogenin-AQP complexes. Most of the fluctuations in binding pocket were observed for AQP0-Gitogenin and AQP4-Diosgenin. However, remaining protein-ligand complex showed stable root mean square deviation (RMSD) plots, strong hydrogen bonding interactions, stable solvent-accessible surface area (SASA) values and minimum distance to the receptor. These observations suggest that natural sapogenin hold promise as novel inhibitors of AQPs, offering a basis for the development of innovative therapeutic agents for cancer treatment. However, further validation of the identified compounds through experiments is essential for translating these findings into therapeutic applications.Communicated by Ramaswamy H. Sarma.

Exploring the role of microbial proteins in controlling environmental pollutants based on molecular simulation

Molecular simulation has been widely used to study microbial proteins' structural composition and dynamic properties, such as volatility, flexibility, and stability at the microscopic scale. Herein, this review describes the key elements of molecular docking and molecular dynamics (MD) simulations in molecular simulation; reviews the techniques combined with molecular simulation, such as crystallography, spectroscopy, molecular biology, and machine learning, to validate simulation results and bridge information gaps in the structure, microenvironmental changes, expression mechanisms, and intensity quantification; illustrates the application of molecular simulation, in characterizing the molecular mechanisms of interaction of microbial proteins with four different types of contaminants, namely heavy metals (HMs), pesticides, dyes and emerging contaminants (ECs). Finally, the review outlines the important role of molecular simulations in the study of microbial proteins for controlling environmental contamination and provides ideas for the application of molecular simulation in screening microbial proteins and incorporating targeted mutagenesis to obtain more effective contaminant control proteins.

Network Pharmacology, Molecular Docking, and Molecular Dynamics Simulation to Elucidate the Molecular Targets and Potential Mechanism of Phoenix dactylifera (Ajwa Dates) against Candidiasis

Candidiasis, caused by opportunistic fungal pathogens of the Candida genus, poses a significant threat to immunocompromised individuals. Natural compounds derived from medicinal plants have gained attention as potential sources of anti-fungal agents. Ajwa dates (Phoenix dactylifera L.) have been recognized for their diverse phytochemical composition and therapeutic potential. In this study, we employed a multi-faceted approach to explore the anti-candidiasis potential of Ajwa dates' phytochemicals. Utilizing network pharmacology, we constructed an interaction network to elucidate the intricate relationships between Ajwa dates phytoconstituents and the Candida-associated molecular targets of humans. Our analysis revealed key nodes in the network (STAT3, IL-2, PTPRC, STAT1, CASP1, ALB, TP53, TLR4, TNF and PPARG), suggesting the potential modulation of several crucial processes (the regulation of the response to a cytokine stimulus, regulation of the inflammatory response, positive regulation of cytokine production, cellular response to external stimulus, etc.) and fungal pathways (Th17 cell differentiation, the Toll-like receptor signaling pathway, the C-type lectin receptor signaling pathway and necroptosis). To validate these findings, molecular docking studies were conducted, revealing the binding affinities of the phytochemicals towards selected Candida protein targets of humans (ALB-rutin (-9.7 kJ/mol), STAT1-rutin (-9.2 kJ/mol), STAT3-isoquercetin (-8.7 kJ/mol), IL2-β-carotene (-8.5 kJ/mol), CASP1-β-carotene (-8.2 kJ/mol), TP53-isoquercetin (-8.8 kJ/mol), PPARG-luteolin (-8.3 kJ/mol), TNF-βcarotene (-7.7 kJ/mol), TLR4-rutin (-7.4 kJ/mol) and PTPRC-rutin (-7.0 kJ/mol)). Furthermore, molecular dynamics simulations of rutin-ALB and rutin-STAT1 complex were performed to gain insights into the stability and dynamics of the identified ligand-target complexes over time. Overall, the results not only contribute to the understanding of the molecular interactions underlying the anti-fungal potential of specific phytochemicals of Ajwa dates in humans but also provide a rational basis for the development of novel therapeutic strategies against candidiasis in humans. This study underscores the significance of network pharmacology, molecular docking and dynamics simulations in accelerating the discovery of natural products as effective anti-fungal agents. However, further experimental validation of the identified compounds is warranted to translate these findings into practical therapeutic applications.

Integrated Network Pharmacology, Molecular Docking, Molecular Simulation, and In Vitro Validation Revealed the Bioactive Components in Soy-Fermented Food Products and the Underlying Mechanistic Pathways in Lung Cancer

Globally, lung cancer remains one of the leading causes of cancer-related mortality, warranting the exploration of novel and effective therapeutic approaches. Soy-fermented food products have long been associated with potential health benefits, including anticancer properties. There is still a lack of understanding of the active components of these drugs as well as their underlying mechanistic pathways responsible for their anti-lung cancer effects. In this study, we have undertaken an integrated approach combining network pharmacology and molecular docking to elucidate the mechanism of action of soy-fermented food products against lung cancer through simulation and in vitro validation. Using network pharmacology, we constructed a comprehensive network of interactions between the identified isoflavones in soy-fermented food products and lung cancer-associated targets. Molecular docking was performed to predict the binding affinities of these compounds with key lung cancer-related proteins. Additionally, molecular simulation was utilized to investigate the stability of the compound-target complexes over time, providing insights into their dynamic interactions. Our results identified daidzein as a potential active component in soy-fermented food products with high binding affinities towards critical lung cancer targets. Molecular dynamic simulations confirmed the stability of the daidzein-MMP9 and daidzein-HSP90AA1 complexes, suggesting their potential as effective inhibitors. Additionally, in vitro validation experiments demonstrated that treatment with daidzein significantly inhibited cancer cell proliferation and suppressed cancer cell migration and the invasion of A549 lung cancer cells. Consequently, the estrogen signaling pathway was recognized as the pathway modulated by daidzein against lung cancer. Overall, the findings of the present study highlight the therapeutic potential of soy-fermented food products in lung cancer treatment and provide valuable insights for the development of targeted therapies using the identified bioactive compounds. Further investigation and clinical studies are warranted to validate these findings and translate them into clinical applications for improved lung cancer management.

Interaction of the nanobio-based reagent with sodium fluorescein and lipids via bioinformatics for forensic fingerprint visualisations

Being commonly found at crime scenes, fingerprints are crucial for human identification, attributable to their uniqueness, persistence and systematic classification of ridge patterns. In addition to latent fingerprints being invisible to the naked eye, the escalating trends of disposing forensic evidence bearing such prints in watery bodies would further complicate criminal investigations. Taking into account the toxicity of small particle reagent (SPR) commonly used in visualising latent fingerprints on wet and non-porous objects, a greener alternative using the nanobio-based reagent (NBR) has been suggested. However, NBR only applies to white and/or relatively light-coloured objects. Thus, conjugation of sodium fluorescein dye with NBR (f-NBR) may be beneficial for increasing the contrast of fingerprint on multi-colored objects. Hence, this study was aimed at investigating the possibility of such conjugation (i.e., f-NBR) as well as proposing suitable interactions between the f-NBR and lipid constituents of fingerprints (tetra-, hexa- and octadecanoic acids) via molecular docking and molecular dynamics simulations. The binding energies between CRL with its ligands were observed at -8.1, -5.0, -4.9 and -3.6 kcal/mole for sodium fluorescein, tetra-, hexa- and octadecanoic acids, respectively. Besides, the formations of hydrogen bonds observed in all complexes (ranged between 2.6 and 3.4 Å), further supported by the stabilized root mean square deviation (RMSDs) plots in MD simulations. In short, the conjugation of f-NBR was computationally feasible, and thereby merits further investigations in the laboratory.Communicated by Ramaswamy H. Sarma.

Identifying novel and potent inhibitors of EGFR protein for the drug development against the breast cancer

The epidermal growth factor receptor (EGFR) has been shown to be extremely important in numerous signaling pathways, particularly those involved in cancer progression. Many therapeutic inhibitors, consisting of both small molecules and monoclonal antibodies, have been developed to target inflammatory, triple-negative and metastatic breast cancer. With the emergence of resistance in breast cancer treatment strategies, there is a need to develop novel drug targets that not only overcome resistance, but also exhibit low toxicity and high specificity. The work presented here focuses on the identification of new inhibitors against the EGFR protein using combined computational approaches. Using a comprehensive machine learning-based virtual screening approach complemented by other computational approaches, we identified six new molecules from the ZINC database. The gold docking score of these six novel molecules is 125.95, 125.38, 123.13, 119.71, 115.64 and 113.73, respectively, while the gold score of the control group is 120.74. In addition, we also analyzed the FEC value of these compounds and found that the values of compounds 1, 2, 3 and 4 (-61.82, -63.98, -67.98 and -63.32, respectively) were higher are than those of the control group (-61.05). Furthermore, these molecules showed highly stable RMSD plots and good interaction of hydrogen bonds. The identified inhibitors provided interesting insights for understanding the electronic, hydrophobic, steric and structural requirements for EGFR inhibitory activity. Distinguishing these novel molecules could lead to the development of new drugs useful in treating breast cancer.Communicated by Ramaswamy H. Sarma.

CB1 as a novel target for Ginkgo biloba's terpene trilactone for controlling chemotherapy-induced peripheral neuropathy (CIPN)

The application of antineoplastic chemotherapeutic agents causes a common side effect known as chemotherapy-induced peripheral neuropathy (CIPN) that leads to reducing the quality of patient's life. This research involves the performance of molecular docking and molecular dynamic (MD) simulation studies to explore the impact of terpenoids of Ginkgo biloba on the targets (CB-1, TLR4, FAAH-1, COX-1, COX-2) that can significantly affect the controlling of CIPN's symptoms. According to the in-vitro and in-vivo investigations, terpenoids, particularly ginkgolides B, A, and bilobalide, can cause significant effects on neuropathic pain. The molecular docking results disclosed the tendency of our ligands to interact with mainly CB1 and FAAH-1, as well as partly with TLR4, throughout their interactions with targets. Terpene trilactone can exhibit a lower rate of binding energy than CB1's inhibitor (7dy), while being precisely located in the CB1's active site and capable of inducing stable interactions by forming hydrogen bonds. The analyses of MD simulation proved that ginkgolide B was a more suitable activator and inhibitor for CB1 and TLR4, respectively, when compared to bilobalide and ginkgolide A. Moreover, bilobalide is capable of inhibiting FAAH-1 more effectively than the two other ligands. According to the analyses of ADME, every three ligands followed the Lipinski's rule of five. Considering these facts, the exertion of three ligands is recommended for their anti-inflammatory, neuroprotective, and anti-nociception influences caused by primarily activating CB1 and inhibiting FAAH-1 and TLR4; in this regard, these compounds can stand as potential candidates for the control and treatment of CIPN's symptoms.

Biotransformation of sulfamethoxazole by microalgae: Removal efficiency, pathways, and mechanisms

Recently, the biotransformation of sulfamethoxazole (SMX) by microalgae has attracted increasing interest. In particular, cytochrome P450 (CYP450) has been suggested to be the main enzymatic contributor to this biodegradation. However, the molecular evidence of CYP450 enzymes being involved in SMX biodegradation remains relatively unclear, hindering its applicability. Herein, the biodegradation of SMX by Chlorella sorokiniana (C. sorokiniana) was investigated, and comprehensively elucidated the reaction mechanism underlying CYP450-mediated SMX metabolism. C. sorokiniana was able to efficiently remove over 80% of SMX mainly through biodegradation, in which CYP450 enzymes responded substantially to metabolize SMX in cells. Additionally, screening of transformation products (TPs) revealed that N⁴-hydroxylation-SMX (TP270) was the main TP in the SMX biodegradation pathway of microalgae. Molecular dynamics (MD) simulation suggested that the aniline of SMX was the most prone to undergo metabolism, while density functional theory (DFT) indicated that SMX was metabolized by CYP450 enzymes through H-abstraction-OH-rebound reaction. Collectively, this work reveals key details of the hydroxylamine group of SMX, elucidates the SMX biodegradation pathway involving CYP450 in microalgae in detail, and accelerates the development of using microalgae-mediated CYP450 to eliminate antibiotics.

Tyrosinase-functionalized polyhydroxyalkanoate bio-beads as a novel biocatalyst for degradation of bisphenol analogues

Bisphenol compounds are emerging contaminants of high concerns with known endocrine-disrupting effects. Biocatalysis provides a green chemistry alternative for advanced treatment in water reclamation. This study createda novel biocatalyst through genetically immobilizing the Bacillus megaterium tyrosinase enzyme (BmTyr) on the surface ofself-assembled polyhydroxyalkanoate (PHA) biopolymer beads (termed PHA-BmTyr) by using synthetic biology techniques and demonstrated one-pot in vivo production of the biocatalyst for effective degradation and detoxification of various bisphenol analogues for the first time. The degradation pathway of bisphenols was determined to be mediated by the monophenolase and diphenolase activity of BmTyr. Notably, biocatalytic bisphenol degradation by PHA-BmTyr could substantially reduce or eliminate estrogenic activity of the contaminants, and the degradation products had remarkably lower acute and chronic toxicity than their parent compounds. Furthermore, the PHA-BmTyr biocatalyst had high reusability for multiple bisphenol degradation reaction cycles and showed excellent stability that retained 100% and 86.6% of the initial activity when stored at 4 °C and room temperature, respectively for 30 days. Also, the PHA-BmTyr biocatalyst could efficiently degrade bisphenol analogues in real wastewater effluent matrix. This study provides a promising approach to develop innovative biocatalysis technologies for sustainable water reclamation.

Comparative molecular docking and molecular-dynamic simulation of wild-type- and mutant carboxylesterase with BTA-hydrolase for enhanced binding to plastic

According to the literature review, microbial degradation of polyethylene terephthalate by PETases has been detected effective and eco-friendly. However, the number of microorganisms capable of such feats is limited with some undesirable bioprospecting results. BTA-hydrolase has been already reported capable of degrading polyethylene terephthalate. Therefore, mutation by in silico site-directed mutagenesis means to introduce current isomer of PETase for polyethylene terephthalate degradative capability as a better approach to resolve this issue. This study aimed to use in silico site-directed mutagenesis to convert a carboxylesterase from Archaeoglobus fulgidus to BTA-hydrolase from Thermobifida fusca by replacing six amino acids in specific locations. This work was followed by molecular docking analysis with polyethylene terephthalate and polypropylene to compare their interactions. The best-docked enzyme-substrate complex was further subjected to molecular dynamics simulation to gauge the binding quality of the BTA-hydrolase, wild-type and mutant-carboxylesterase with only polyethylene terephthalate as a substrate. Results of molecular docking revealed lowest binding energy for the wild-type carboxylesterase-polypropylene complex (-7.5 kcal/mol). The root-mean-square deviation value was observed stable for BTA-hydrolase. Meanwhile, root-mean-square fluctuation was assessed with higher fluctuation for the mutated residue Lys178. Consequently, the Rg value for BTA-hydrolase-ligand complex (∼1.68 nm) was the lowest compared to the mutant and wild-type carboxylesterase. The collective data conveyed that mutations imparted a minimal change in the ability of the mutant carboxylesterase to bind to polyethylene terephthalate.

Computational Assessment of Botrytis cinerea Lipase for Biofuel Production

The demand for ecofriendly green catalysts for biofuel synthesis is greatly increasing with the effects of fossil fuel depletion. Fungal lipases are abundantly used as biocatalysts for the synthesis of biofuel. The use of Botrytis cinerea lipase is an excellent approach for the conversion of agroindustrial residues into biofuel. In this study, phylogenetic analyses were carried out and the physicochemical properties of B. cinerea lipase were assessed. Furthermore, the protein structure of B. cinerea lipase was predicted and refined. Putative energy-rich phytolipid compounds were explored as a substrate for the synthesis of biofuel, owing to B. cinerea lipase catalysis. Approximately 161 plant-based fatty acids were docked with B. cinerea lipase in order to evaluate their binding affinities and interactions. Among the docked fatty acids, the top ten triglycerides having the lowest number of binding affinities with B. cinerea lipase were selected, and their interactions were assessed. The top three triglycerides having the greatest number of hydrogen bonds and hydrophobic interactions were selected for simulations of 20 ns. The docking and simulations revealed that docosahexaenoic acid, dicranin, and hexadeca-7,10,13-trienoic acid had stable bonding with the B. cinerea lipase. Therefore, B. cinerea lipase has the potential to be used for the transesterification of fatty acids into biofuels, whereas docosahexaenoic acid, dicranin, and hexadeca-7,10,13-trienoic acid can be used as substrates of B. cinerea lipase for biofuel synthesis.

Unraveling the potential role of bioactive molecules produced by Trichoderma spp. as inhibitors of tomatinase enzyme having an important role in wilting disease: an in-silico approach

Tomatinase; a saponin detoxification enzyme produced by Fusarium oxysporumf.sp. lycopersici is reported as a causative agent for wilting disease in tomato crops. The disease is instigated by inhibiting the activity of α-tomatine. Trichoderma spp. widely used as biocontrol agent play an essential role in plant growth and pathogen control. In the current study, an in-silico approach using substrate docking, molecular dynamics and MM/PBSA analysis was used to evaluate the potential role of bioactive metabolites produced by Trichoderma spp. The study aims to establish the efficacy of catalytic tendency of the bioactive metabolites to combat the effect of tomatinase enzyme employing α-tomatine as the substrate. By means of the integrated molecular modeling approach; novel bioactive metabolites namely, Trichodermamide B, Trichosetin and Virone were found to be the potential inhibitors against tomatinase enzyme secreted by Fusarium oxysporum f.sp. lycopersici. Molecular dynamic (MD) simulations displayed that the screened ligands bound tomatinase during 150 ns of MD simulations. Furthermore, the (MM-PBSA) free energy calculations depicted that screened molecules possess stable and favorable energies for Trichodermamide B (-7.1 kcal/mol), Trichosetin (-7.4 kcal/mol) and Virone (-7.9 kcal/mol) thereby instigating robust binding with the enzyme's binding site. The results attained in this study, reflects that these bioactive metabolites may serve as potential substrates to control and inhibit the tomatinase enzyme; playing an integral role in combating the wilt disease.Communicated by Ramaswamy H. Sarma.

Molecular Dynamics Simulation of Nanocellulose-Stabilized Pickering Emulsions

While the economy is rapidly expanding in most emerging countries, issues coupled with a higher population has created foreseeable tension among food, water, and energy. It is crucial for more sustainable valorization of resources, for instance, nanocellulose, to address the core challenges in environmental sustainability. As the complexity of the system evolved, the timescale of project development has increased exponentially. However, research on the design and operation of integrated nanomaterials, along with energy supply, monitoring, and control infrastructure, has seriously lagged. The development cost of new materials can be significantly reduced by utilizing molecular simulation technology in the design of nanostructured materials. To realize its potential, nanocellulose, an amphiphilic biopolymer with the presence of rich -OH and -CH structural groups, was investigated via molecular dynamics simulation to reveal its full potential as Pickering emulsion stabilizer at the molecular level. This work has successfully quantified the Pickering stabilization mechanism profiles by nanocellulose, and the phenomenon could be visualized in three stages, namely the initial homogenous phase, rapid formation of micelles and coalescence, and lastly the thermodynamic equilibrium of the system. It was also observed that the high bead order was always coupled with a high volume of phase separation activities, through a coarse-grained model within 20,000 time steps. The outcome of this work would be helpful to provide an important perspective for the future design and development of nanocellulose-based emulsion products, which cater for food, cosmeceutical, and pharmaceutical industries.

Molecular modeling, docking and simulation dynamics of β-glucosidase reveals high-efficiency, thermo-stable, glucose tolerant enzyme in Paenibacillus lautus BHU3 strain

The enzyme β-glucosidase mediates the rate limiting step of conversion of cellobiose to glucose and thus plays a vital role in the process of cellulose degradation. The present study deals with analysis of the effective novel strain of Paenibacillus lautus BHU3 for identifying high-efficiency thermostable, glucose tolerant β-glucosidases. Seven counterparts with elevated T_m values ranging from 64.6 to 75.8 °C with high thermo-stability, were revealed through this analysis. The blind molecular docking of the model enzymes structures with cellobiose and pNPG gave high negative interaction energies ranging from -11.33 to -13.29 and -6.43 to -9.054 (kcal mol^-1), respectively. The enzyme WP_096774744.1 effectively formed 5 hydrogen bonds with the highest interaction energy (-13.29 kcal mol^-1) with cellobiose at its catalytic site. Molecular dynamics simulation analysis performed for the WP_096774744.1-pNPG complex predicted Glu5, Arg7, Lue68, Gly69 and Phe325 as the major contributing residues for accomplishing hydrolysis of β-1-4-linkage. Further, the molecular docking of WP_096774744.1 enzyme with glucose revealed a distinct glucose-binding site distant from the substrate-binding site, thus confirming the deficient competitive inhibition by glucose. Hence, WP_096774744.1 β-glucosidase appears to be an efficient enzyme with enhanced activity to biodegrade the cellulosic materials and highly relevant for waste management and various industrial applications.

Interaction of indole-3-acetic acid with horseradish peroxidase as a potential anticancer agent: from docking to molecular dynamics simulation

The oxidation process, catalyzed by the peroxidase enzymes, occurs in all domains of life to detoxify the hydrogen peroxide toxicity. The most well-known, applicable and vastly studied member of the peroxidases family is horseradish peroxidase (HRP), especially the isoenzyme C (HRP C). HRP (primarily HRP C) is commercially available and applicable in biotechnology and diagnosis. Recently, a novel application of HRP has been introduced in cancer therapy as the combination of HRP with indole-3-acetic acid (IAA). The anticancer activity of HRP/IAA complex is through oxidation of IAA by HRP in hypoxic tumor condition, which leads to apoptosis and cancerous cell death. However, the molecular interaction of HRP/IAA has not been elucidated. Identifying the interaction of IAA with HRP would provide a better insight into its function and applications. In this study, molecular docking and molecular dynamics (MD) simulation were applied to determine the molecular interaction of the IAA/HRP complex. The docking study represented that IAA bound at the 'exposed' heme edge of the HRP enzyme, and the IAA entrance to the enzyme was situated at the carboxymethyl side-chain of the selected structure. Our computational results showed the HRP/IAA complex structure stability. While hydrogen bond formation with ARG38 and HIS42 stabilized the substrate, hydrophobic interactions with Phe68, Gly69, Leu138, Pro139, Pro141 and Phe179 contributed to IAA/HRP complex stability. The results can help to better understand peroxidase enzyme activity and would pave the way for future development of new therapeutics with improved anticancer efficacy.Communicated by Ramaswamy H. Sarma.

Artichoke pectic oligosaccharide characterisation and virtual screening of prebiotic properties using in silico colonic fermentation

The aim of this work was to develop a comprehensive workflow to elucidate molecular features of artichoke pectic oligosaccharides (POS) contributing to high potential prebiotic activity. First, obtainment of artichoke POS by Pectinex® Ultra-Olio was optimised using an artificial neural network. Under optimal conditions (pH 6.86; 1.5 h; enzyme dose 520.5 U/g pectin) POS yield was 624 mg/g pectin. Oligosaccharide structures (M_w < 1.3 kDa) were characterised by MALDI-TOF-MS. Then, conformational analysis of glycosidic bonds was performed by replica exchange molecular dynamics simulations and interaction mechanisms between POS and several microbial glycosidases were proposed by molecular modelling. Chemical information was integrated in virtual simulations of colonic fermentation. Highest hydrolysis rate was obtained for GalA-Rha-GalA trisaccharide, while the presence of partial negative charges and high radius of gyration enhance short chain fatty acid formation in distal colon. Established structure-activity relationships could help the rational design of prebiotics and clinical trials.

Molecular docking and molecular dynamics simulation of Bacillus thuringiensis dehalogenase against haloacids, haloacetates and chlorpyrifos

The high dependency and surplus use of agrochemical products have liberated enormous quantities of toxic halogenated pollutants into the environment and threaten the well-being of humankind. Herein, this study performed molecular docking, molecular dynamic (MD) simulations, molecular mechanics-Poisson Boltzmann Surface Area (MM-PBSA) calculations on the DehH2 from Bacillus thuringiensis, to identify the order of which the enzyme degrades different substrates, haloacids, haloacetate and chlorpyrifos. The study discovered that the DehH2 favored the degradation of haloacids and haloacetates (-3.3 - 4.6 kcal/mol) and formed three hydrogen bonds with Asp125, Arg201 and Lys202. Despite the inconclusive molecular docking result, chlorpyrifos was consistently shown to be the least favored substrate of the DehH2 in MD simulations and MM-PBSA calculations. Results of MD simulations revealed the DehH2-haloacid- (RMSD 0.15 - 0.25 nm) and DehH2-haloacetates (RMSF 0.05 - 0.25 nm) were more stable, with the DehH2-L-2CP complex being the most stable while the least was the DehH2-chlorpyrifos (RMSD 0.295 nm; RMSF 0.05 - 0.59 nm). The Molecular Mechanics Poisson-Boltzmann Surface Area calculations showed the DehH2-L-2CP complex (-24.27 kcal/mol) having the lowest binding energy followed by DehH2-MCA (-22.78 kcal/mol), DehH2-D-2CP (-21.82 kcal/mol), DehH2-3CP (-21.11 kcal/mol), DehH2-2,2-DCP (-18.34 kcal/mol), DehH2-2,3-DCP (-8.34 kcal/mol), DehH2-TCA (-7.62 kcal/mol), while chlorpyrifos was unable to spontaneously bind to DehH2 (+127.16 kcal/mol). In a nutshell, the findings of this study offer valuable insights into the rational tailoring of the DehH2 for expanding its substrate specificity and catalytic activity in the near future.Communicated by Ramaswamy H. Sarma.