Combined molecular docking, homology modelling and density functional theory studies to modify dioxygenase to efficiently degrade aromatic hydrocarbons

Nutrient amendments enrich microbial hydrocarbon degradation metagenomic potential in freshwater coastal wetland microcosm experiments

Biostimulating native microbes with fertilizers has proven to be a highly effective strategy to speed up biodegradation rates in microbial communities. This study investigates the genetic potential of microbes to degrade light synthetic crude oil in a freshwater coastal wetland. Experimental sediment microcosms were exposed to a variety of conditions (biological control, a light synthetic crude oil amendment, and light synthetic crude oil with nutrient amendment) and incubated for 30 days before volatile organic compounds (BTEX) were quantified and DNA was sequenced for metagenomic analysis. The resulting DNA sequences were binned into metagenome-assembled genomes (MAGs). Analyses of MAGs uncovered a 13-fold significant increase in the abundance of rate-limiting hydrocarbon degrading monooxygenases and dioxygenases, identified only in MAGs from the light synthetic crude oil with nutrient amendments. Further, complete degradation pathways for BTEX compounds were found only in MAGs resulting from the light synthetic crude with nutrient amendment. Moreover, volatile organic compounds (BTEX, cyclohexane, and naphthalene) analyses of microcosm sediments in the presence of nutrients documented that benzene was degraded below detection limits, toluene (98%) and ethylbenzene (67%) were predominantly reduced within 30 days. Results indicate that the genetic potential to degrade BTEX compounds in this freshwater wetland can be linked to the functional potential for bioremediation. BTEX compounds are typically more recalcitrant and tougher to degrade than alkanes. This study demonstrated that stimulating a microbial community with nutrients to enhance its ability to biodegrade hydrocarbons, even in a relatively nutrient-rich habitat like a freshwater wetland, is an effective remediation tactic.

IMPORTANCE

The impact of oil spills in a freshwater aquatic environment can pose dire social, economic, and ecological effects on the region. An oil spill in the Laurentian Great Lakes region has the potential to affect the drinking water of more than 30 million people. The light synthetic crude oil used in this experimental microcosm study is transported through an underground pipeline crossing the waterway between two Laurentian Great Lakes. This study collected metagenomic data (experiments in triplicate) and assessed the quantity of BTEX compounds, which connected microbial degradation function to gene potential. The resulting data documented the bioremediation capabilities of native microbes in a freshwater coastal wetland. This study also provided evidence for this region that bioremediation can be a viable remediation strategy instead of invasive physical methods.

Cultivation and characterization of functional-yet-uncultivable phenanthrene degraders by stable-isotope-probing and metagenomic-binning directed cultivation (SIP-MDC)

High-throughput identification and cultivation of functional-yet-uncultivable microorganisms is a fundamental goal in environmental microbiology. It remains as a critical challenge due to the lack of routine and effective approaches. Here, we firstly proposed an approach of stable-isotope-probing and metagenomic-binning directed cultivation (SIP-MDC) to isolate and characterize the active phenanthrene degraders from petroleum-contaminated soils. From SIP and metagenome, we assembled 13 high-quality metagenomic bins from 13C-DNA, and successfully obtained the genome of an active PHE degrader Achromobacter (genome-MB) from 13C-DNA metagenomes, which was confirmed by gyrB gene comparison and average nucleotide/amino identity (ANI/AAI), as well as the quantification of PAH dioxygenase and antibiotic resistance genes. Thereinto, we modified the traditional cultivation medium with antibiotics and specific growth factors (e.g., vitamins and metals), and separated an active phenanthrene degrader Achromobacter sp. LJB-25 via directed isolation. Strain LJB-25 could degrade phenanthrene and its identity was confirmed by ANI/AAI values between its genome and genome-MB (>99 %). Our results hinted at the feasibility of SIP-MDC to identify, isolate and cultivate functional-yet-uncultivable microorganisms (active phenanthrene degraders) from their natural habitats. Our findings developed a state-of-the-art SIP-MDC approach, expanded our knowledge on phenanthrene biodegradation mechanisms, and proposed a strategy to mine functional-yet-uncultivable microorganisms.

Fungal enzyme degradation of lignin-PLA composites: Insights from experiments and molecular docking simulations

Fungal enzymes are effective in degrading various polymeric materials. In this study, we assessed the initial degradation of composites consisting of lignin-poly(lactic acid) (PLA) with both unmodified lignin (LIG) and oxypropylated lignin (oLIG) incorporated at 10 % and 40 % weight within the PLA matrix in a fungal environment. Trametes versicolor fungi were used, and the samples were treated only for eight weeks. Although there was no significant difference in weight loss, the degradation process impacted the chemical and thermal properties of the composites, as shown by Fourier transform infrared spectroscopy (FTIR) and Differential scanning calorimetry (DSC) analyses. After the degradation process, the carbonyl index values decreased for all composites and the hydroxyl index values increased for LIG/PLA and a reverse trend was observed for oLIG/PLA composites. The first heating scan from DSC results showed that the melting peak and the cold crystallization peak disappeared after the degradation process. Microscopic analysis revealed that LIG/PLA exhibited higher roughness than oLIG/PLA. Molecular docking simulations were carried out using guaiacylglycerol-β-guaiacyl ether (GGE) and lactic acid (LA) as model compounds for lignin and PLA, respectively, with laccase (Lac) enzyme for Trametes versicolor. The docking results showed that GGE had the strongest binding interaction and affinity with Lac than lactic acid and oxypropylated GGE. The oxypropylated GGE formed a shorter hydrogen bonding with the Lac enzyme than GGE and LA. The trend associated with the degradation of composites from experimental and molecular docking findings was consistent. This combined approach provided insights into the degradation process using fungi and had the potential to be applied to different polymeric composites.

Catalytic resilience of multicomponent aromatic ring-hydroxylating dioxygenases in Pseudomonas for degradation of polycyclic aromatic hydrocarbons

Hydrophobic organic compounds, either natural or introduced through anthropogenic activities, pose a serious threat to all spheres of life, including humankind. These hydrophobic compounds are recalcitrant and difficult to degrade by the microbial system; however, microbes have also evolved their metabolic and degradative potential. Pseudomonas species have been reported to have a multipotential role in the biodegradation of aromatic hydrocarbons through aromatic ring-hydroxylating dioxygenases (ARHDs). The structural complexity of different hydrophobic substrates and their chemically inert nature demands the explicit role of evolutionary conserved multicomponent enzyme ARHDs. These enzymes catalyze ring activation and subsequent oxidation by adding two molecular oxygen atoms onto the vicinal carbon of the aromatic nucleus. This critical metabolic step in the aerobic mode of degradation of polycyclic aromatic hydrocarbons (PAHs) catalyzed by ARHDs can also be explored through protein molecular docking studies. Protein data analysis enables an understanding of molecular processes and monitoring complex biodegradation reactions. This review summarizes the molecular characterization of five ARHDs from Pseudomonas species already reported for PAH degradation. Homology modeling for the amino acid sequences encoding the catalytic α-subunit of ARHDs and their docking analyses with PAHs suggested that the enzyme active sites show flexibility around the catalytic pocket for binding of low molecular weight (LMW) and high molecular weight (HMW) PAH substrates (naphthalene, phenanthrene, pyrene, benzo[α]pyrene). The alpha subunit harbours variable catalytic pockets and broader channels, allowing relaxed enzyme specificity toward PAHs. ARHD's ability to accommodate different LMW and HMW PAHs demonstrates its 'plasticity', meeting the catabolic demand of the PAH degraders.

Bioremediation potential of laccase for catalysis of glyphosate, isoproturon, lignin, and parathion: Molecular docking, dynamics, and simulation

The present study aimed to investigate the catalytic degradation produced by laccase in the detoxification of glyphosate, isoproturon, lignin polymer, and parathion. We explored laccase-glyphosate, laccase-lignin polymer, laccase-isoproturon, and laccase-parathion using molecular docking (MD) and molecular dynamics simulation (MDS) approaches. The results suggest that laccase interacts well with glyphosate, lignin polymer, isoproturon, and parathion during biodegradation. We calculated the root mean square deviations (RMSD) of laccase-glyphosate, laccase-lignin polymer, laccase-isoproturon, and laccase-parathion as 0.24 ± 0.02, 0.59 ± 0.32, 0.43 ± 0.07, and 0.43 ± 0.06 nm, respectively. In an aqueous solution, the stability of laccase with glyphosate, lignin polymer, isoproturon, and parathion is mediated through the formation of hydrophobic interactions, hydrogen bonds, and van der Waals interactions. The presence of xenobiotic toxic compounds in the active site changed the conformation of laccase. MDS of the laccase-substrate complexes confirmed their stability during catalytic degradation. Laccase assay results confirmed that the degradation of syringol, dihydroconiferyl alcohol, guaiacol, parathion, isoproturon, and glyphosate were 100%, 99.31%, 95.69%, 60.96%, 54.51%, and 48.34% within 2 h, respectively. Taken together, we describe a novel method to understand the molecular-level biodegradation of xenobiotic compounds through laccase and its potential application in contaminant removal.

In silico analytical toolset for predictive degradation and toxicity of hazardous pollutants in water sources

Different phenolic compounds, including multimeric lignin derivatives in the β-O-4 form, are among the most prevalent compounds in wastewater, often generated from paper industries. Relatively small concentrations of lignin are hazardous to aquatic organisms and can trigger severe environmental hazards. Herein, we present a predictive toolset to insight the induced toxic hazards prediction, and their Lignin peroxidase (LiP)-assisted degradation mechanism of selected multimeric lignin model compounds. T.E.ST and Toxtree toolset were deployed for toxic hazards estimation in different endpoints. To minimize the concerning hazards, we screened multimeric compounds for binding affinity with LiP. The binding affinity was found to be significantly lower than the reference compound. An Extra precision (XP) Glide score of -6.796 kcal/mol was found for dimer (guaiacyl 4-O-5 guaiacyl) complex as lowest compared to reference compound (-4.007 kcal/mol). The active site residues ASP-153, HIP-226, VAL-227, ARG-244, GLU-215, 239, PHE-261 were identified as site-specific key binding AA residues actively involved with corresponding ligands, forming Hydrophobic, H-Bond, π-Stacking, π-π type interactions. The DESMOND-assisted molecular dynamics simulation's (MDS) trajectories of protein-ligand revealed the considerable binding behavior and attained stability and system equilibrium state. Such theoretical and predictive conclusions indicted the feasibility of LiP assisted sustainable mitigation of lignin-based compounds, and such could be used to protect the environment from the potential hazards posed by recognized similar pollutants.

Unveiling metabolic characteristics of an uncultured Gammaproteobacterium responsible for in situ PAH biodegradation in petroleum polluted soil

Exploring the metabolic characteristics of indigenous PAH degraders is critical to understanding the PAH bioremediation mechanism in the natural environment. While stable-isotopic probing (SIP) is a viable method to identify functional microorganisms in complex environments, the metabolic characteristics of uncultured degraders are still elusive. Here, we investigated the naphthalene (NAP) biodegradation of petroleum polluted soils by combining SIP, amplicon sequencing and metagenome binning. Based on the SIP and amplicon sequencing results, an uncultured Gammaproteobacterium sp. was identified as the key NAP degrader. Additionally, the assembled genome of this uncultured degrader was successfully obtained from the ¹³ C-DNA metagenomes by matching its 16S rRNA gene with the SIP identified OTU sequence. Meanwhile, a number of NAP degrading genes encoding naphthalene/PAH dioxygenases were identified in this genome, further confirming the direct involvement of this indigenous degrader in the NAP degradation. The degrader contained genes related to the metabolisms of several carbon sources, energy substances and vitamins, illuminating potential reasons for why microorganisms cannot be cultivated and finally realize their cultivation. Our findings provide novel information on the mechanisms of in situ PAH biodegradation and add to our current knowledge on the cultivation of non-culturable microorganisms by combining both SIP and metagenome binning.

Fungal enzymes for the degradation of polyethylene: Molecular docking simulation and biodegradation pathway proposal

Polyethylene (PE) is one of the most highly consumed petroleum-based polymers and its accumulation as waste causes environmental pollution. In this sense, the use of microorganisms and their enzymes represents the most ecofriendly and effective decontamination approach. In this work, molecular docking simulation for catalytic enzyme degradation of PE was carried out using individual enzymes: laccase (Lac), manganese peroxidase (MnP), lignin peroxidase (LiP) and unspecific peroxygenase (UnP). PE-binding energy, PE-binding affinity and dimensions of PE-binding sites in the enzyme cavity were calculated in each case. Four hypothetical PE biodegradation pathways were proposed using individual enzymes, and one pathway was proposed using a synergic enzyme combination. These results show that in nature, enzymes act in a synergic manner, using their specific features to undertake an extraordinarily effective sequential catalytic process for organopollutants degradation. In this process, Lac (oxidase) is crucial to provide hydrogen peroxide to the medium to ensure pollutant breakdown. UnP is a versatile enzyme that offers a promising practical application for the degradation of PE and other pollutants due to its cavity features. This is the first in silico report of PE enzymatic degradation, showing the mode of interaction of PE with enzymes as well as the degradation mechanism.

Enhanced plant-microbe remediation of PCBs in soil using enzyme modification technique combined with molecular docking and molecular dynamics

The study on the enhanced mechanisms of the enzymes involved in plant absorption, plant degradation, and microbial mineralization in the remediation of soils contaminated with polychlorinated biphenyls (PCBs) is of great significance for the application of plant-microbe combined remediation technique in PCB-contaminated soils. The present study first used a combination of molecular docking and molecular dynamics methods to calculate the effects of the plant absorption enzyme, plant degradation enzyme, and microbial mineralization enzyme on the PCBs in the soil environment. A multifunctional plant degradation enzyme was constructed with three functional roles of absorption, degradation, and mineralization using amino acid sequence recombination and site-directed mutagenesis to modify the template of plant degradation enzyme. Finally, using the Taguchi experimental design-assisted molecular dynamics simulation method, the suitable external environmental conditions of plant-microbe combined remediation of the PCB-contaminated soil were determined. In total, six multifunctional plant degradation enzymes were designed, which exhibited a significantly improved efficiency of PCB degradation. In comparison to the complex of plant absorption enzyme, plant degradation enzyme, and microorganism mineralization enzyme (6QIM-3GZX-1B85), the six multifunctional plant degradation enzymes exhibited significantly higher efficiency (2.10-2.38 times) in degrading the PCBs, with a maximum of 2.69 times under suitable external environmental conditions.

In silico approach: biological prediction of nordentatin derivatives as anticancer agent inhibitors in the cAMP pathway

A combination of computational techniques has been carried out to predict the binding of nordentatin derivatives based on pyranocoumarin semi-synthesis with the target protein from the expression of the PDE4B gene. The inhibition of the cAMP pathway is the main target of anti-cancer drugs, which is responsible for uncontrolled cell division in cancer. Modeling was done using a combination of semi-empirical methods and the density functional theory (PM3-DFT/6-31G*/B3LYP) to obtain the optimal structure of a small ligand that could be modeled. Studies on the interaction of the ligands and amino acid residues on protein targets were carried out using a combination of molecular docking and molecular dynamic simulation. Molecular docking based on functional grid scores showed a very good native ligand pose with an RMSD of 0.93 Å in determining the initial coordinates of the ligand-receptor interactions. Furthermore, the amino acid residues responsible for interaction through H-bonds were Tyr103, His104, His177, Met217, and Gln313. The binding free energy (kcal mol-1) results of the candidates were PS-1 (-36.84 ± 0.31), PS-2 (-35.34 ± 0.28), PS-3 (-26.65 ± 0.30), PS-5 (-42.66 ± 0.26), PS-7 (-35.33 ± 0.23), and PS-9 (-32.57 ± 0.20), which are smaller than that of the native ligand Z72 (-24.20 ± 0.19), and thus these have good potential as drugs that can inhibit the cAMP pathway. These results provide theoretical information for the efficient inhibition of the cAMP pathway in the future.

Molecular Modification of Fluoroquinolone-Biodegrading Enzymes Based on Molecular Docking and Homology Modelling

To improve the biodegradation efficiency of fluoroquinolone antibiotics during sewage treatment, fluoroquinolone aerobic, anaerobic and facultative degrading enzymes for fluoroquinolone degradation were modified by molecular docking and homology modelling. First, amino acid residues of the binding sites of degrading enzymes for the target fluoroquinolones ciprofloxacin (CIP), norfloxacin (NOR) and ofloxacin (OFL) were analysed by the molecular docking method. The hydrophobic amino acid residues within 5 Å of the target fluoroquinolone molecules were selected as the modification sites. The hydrophobic amino acid residues at the modified sites were replaced by the hydrophilic amino acid residues, and 150 amino acid sequence modification schemes of the degrading enzymes were designed. Subsequently, a reconstruction scheme of the degrading enzyme amino acid sequence reconstruction scheme was submitted to the SWISS-MODEL server and a selected homology modelling method was used to build a new structure of the degrading enzyme. At the same time, the binding affinities between the novel degrading enzymes and the target fluoroquinolones (represented by the docking scoring function) were evaluated by the molecular docking method. It was found that the novel enzymes can simultaneously improve the binding affinities for the three target fluoroquinolones, and the degradation ability of the six modification schemes was increased by more than 50% at the same time. Among the novel enzymes, the affinity effect of the novel anaerobic enzyme (6-1) with CIP, NOR and OFL was significantly increased, with increases of 129.24%, 165.06% and 169.59%, respectively, followed by the facultative enzyme and aerobic enzyme. In addition, the designed degrading enzymes had certain selectivity for the degradation of the target quinolone. Among the novel enzymes, the binding affinities of the novel anaerobic enzyme (6-3) and CIP, the novel aerobic enzyme (3-6) and NOR, and the novel facultative enzyme (13-6) and OFL were increased by 149.71%, 178.57% and 297.12% respectively. Calculations using the Gaussian09 software revealed that the degradation reaction barrier of the novel degrading enzyme (7-1) and CIP NOR and OFL decreased by 37.65 kcal·mol-1, 6.28 kcal·mol-1 and 6.28 kcal·mol-1, respectively, which would result in efficient degradation of the target fluoroquinolone molecules. By analysing the binding affinity of the degrading enzymes before and after the modification with methanol, it was further speculated that the degradation effect of the modified aerobic degrading enzymes on organic matter was lower than that before the modification, and the increase or decrease in the degradation effect was less than 10%. The mechanism analysis found that the interaction between the modified amino acid residues of the degrading enzymes and the fluoroquinolone molecules increased. The average distance between the amino acid residues and the fluoroquinolone molecules represented a comprehensive affinity effect, and its value was positively correlated with the degradation effect of the novel degrading enzymes.