In silico characterization of a novel dehalogenase (DehHX) from the halophile Pseudomonas halophila HX isolated from Tuz Gölü Lake, Turkey: insights into a hypersaline-adapted dehalogenase

Biological and molecular approaches of the degradation or decolorization potential of the hypersaline Lake Tuz Bacillus megaterium H2 isolate

The presence of synthetic dyes in water bodies and soil is one of the major issues affecting the global ecology, possibly impacting societal well-being adversely due to the colorants' recalcitrance and toxicity. Herein, the study spectrophotometrically monitored the ability of the Bacillus megaterium H2 azoreductase (AzrBmH2) to degrade four synthetic dyes, reactive blue 4, remazol brilliant red, thymol blue, and methyl red, followed by in-silico assessment using GROMACS. We found that the bacterium degraded as much as 60% of all four synthetic dyes at various tested concentrations. The genome analysis revealed five different azoreductase genes, which were then modeled into the AzrBmH21, AzrBmH22/3, and AzrBmH24/5 templates. The AzrBmH2-substrate complexes showed binding energies with all the dyes of between -10.6 to -6.9 kcal/mol and formed 4-6 hydrogen bonds with the predicted catalytic binding residues (His10, Glu 14, Ser 58, Met 99, Val 107, His 183, Asn184 and Gln 191). In contrast, the lowest binding energies were observed for the AzrBmH21-substrates (-10.6 to -7.9). Molecular dynamic simulations revealed that the AzrBmH21-substrate complexes were more stable (RMSD 0.2-0.25 nm, RMSF 0.05 - 0.3 nm) and implied strong bonding with the dyes. The Molecular Mechanics Poisson-Boltzmann Surface Area results also mirrored this outcome, showing the lowest azoreductase-dye binding energy in the order of AzrBmH21-RB4 (-78.18 ± 8.92 kcal/mol), AzrBmH21-RBR (-67.51 ± 7.74 kcal/mol), AzrBmH21-TB (-46.62 ± 5.23 kcal/mol) and AzrBmH21-MR (-40.78 ± 7.87 kcal/mol). In short, the study demonstrated the ability of the B. megaterium H2 to efficiently decolorize the above-said synthetic dyes, conveying the bacterium's promising use for large-scale dye remediation.Communicated by Ramaswamy H. Sarma.

In silico analysis of a putative dehalogenase from the genome of halophilic bacterium Halomonas smyrnensis AAD6T

Microbial-assisted removal of natural or synthetic pollutants is the prevailing green, low-cost technology to treat polluted environments. However, the challenge with enzyme-assisted bioremediation is the laborious nature of dehalogenase-producing microorganisms' bioprospecting. This bottleneck could be circumvented by in-silico analysis of certain microorganisms' whole-genome sequences to predict their protein functions and enzyme versatility for improved biotechnological applications. Herein, this study performed structural analysis on a dehalogenase (DehHsAAD6) from the genome of Halomonas smyrnensis AAD6 by molecular docking and molecular dynamic (MD) simulations. Other bioinformatics tools were also employed to identify substrate preference (haloacids and haloacetates) of the DehHsAAD6. The DehHsAAD6 preferentially degraded haloacids and haloacetates (-3.2-4.8 kcal/mol) and which formed three hydrogen bonds with Tyr12, Lys46, and Asp182. MD simulations data revealed the higher stability of DehHsAAD6-haloacid- (RMSD 0.22-0.3 nm) and DehHsAAD6-haloacetates (RMSF 0.05-0.14 nm) complexes, with the DehHsAAD6-L-2CP complex being the most stable. The detail of molecular docking calculations ranked complexes with the lowest binding free energies as: DehHsAAD6-L-2CP complex (-4.8 kcal/mol) = DehHsAAD6-MCA (-4.8 kcal/mol) < DehHsAAD6-TCA (-4.5 kcal/mol) < DehHsAAD6-2,3-DCP (-4.1 kcal/mol) < DehHsAAD6-D-2CP (-3.9 kcal/mol) < DehHsAAD6-2,2-DCP (-3.5 kcal/mol) < DehHsAAD6-3CP (-3.2 kcal/mol). In a nutshell, the study findings offer valuable perceptions into the elucidation of possible reaction mechanisms of dehalogenases for extended substrate specificity and higher catalytic activity.Communicated by Ramaswamy H. Sarma.

Molecular interactions of trichoderma β-1,4-glucosidase (ThBglT12) with mycelial cell wall components of phytopathogenic Macrophomina phaseolina

Efficacy of a β-1,4-glucosidase from Trichoderma harzianum T12 (ThBglT12) in disrupting the cell wall of the phytopathogenic fungus M. phaseolina (Macrophomina phaseolina) was studied, as the underlying molecular mechanisms of cell wall recognition remains elusive. In this study, the binding location identified by a consensus of residues predicted by COACH tool, blind docking, and multiple sequence alignment revealed that molecular recognition by ThBglT12 occurred through interactions between the α-1,3-glucan, β-1,3-glucan, β-1,3/1,4-glucan, and chitin components of M. phaseolina, with corresponding binding energies of -7.4, -7.6, -7.5 and -7.8 kcal/mol. The residue consensus verified the participation of Glu172, Tyr304, Trp345, Glu373, Glu430, and Trp431 in the active site pocket of ThBglT12 to bind the ligands, of which Trp345 was the common interacting residue. Root mean square deviation (RMSD), root mean square fluctuation (RMSF), total energy, and minimum distance calculation from molecular dynamics (MD) simulation further confirmed the stability and the closeness of the binding ligands into the ThBglT12 active site pocket. The h-bond occupancy by Glu373 and Trp431 instated the role of the nucleophile for substrate recognition and specificity, crucial for cleaving the β-1,4 linkage. Further investigation showed that the proximity of Glu373 to the anomeric carbon of β-1,3/1,4-glucan (3.5 Å) and chitin (5.5 Å) indicates the nucleophiles' readiness to form enzyme-substrate intermediates. Plus, the neighboring water molecule appeared to be correctly positioned and oriented towards the anomeric carbon to hydrolyze the β-1,3/1,4-glucan and chitin, in less than 4.0 Å. In a nutshell, the study verified that the ThBglT12 is a good alternative fungicide to inhibit the growth of M. phaseolina.Communicated by Ramaswamy H. Sarma.

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.

In silico assessment of dehalogenase from Bacillus thuringiensis H2 in relation to its salinity-stability and pollutants degradation

Increased scientific interest has led to the rise in biotechnological uses of halophilic and halotolerant microbes for hypersaline wastewater bioremediation. Hence, this study performed molecular docking, molecular dynamic (MD) simulations, and validation by Molecular Mechanic Poisson-Boltzmann Surface Area (MM-PBSA) calculations on the DehH2 from Bacillus thuringiensis H2. We aimed to identify the interactions of DehH2 with substrates haloacids, haloacetates, and chlorpyrifos under extreme salinity (35% NaCl). MD simulations revealed that DehH2 preferentially degraded haloacids and haloacetates (-6.3 to -4.7 kcal/mol) by forming three or four hydrogen bonds to the catalytic triad, Asp125, Arg201, and Lys202. Conversely, chlorpyrifos was the least preferred substrate in both MD simulations and MM-PBSA calculations. MD simulation results ranked the DehH2-L-2CP complex (RMSD □0.125-0.23 nm) as the most stable while the least was the DehH2-chlorpyrifos complex (RMSD 0.32 nm; RMSF 0.0 - 0.29). The order of stability was as follows: DehH2-L-2CP > DehH2-MCA > DehH2-D-2CP > DehH2-3CP > DehH2-2,2-DCP > DehH2-2,3-DCP > DehH2-TCA > DehH2-chlorpyrifos. The MM-PBSA calculations further affirmed the DehH2-L-2CP complex's highest stability with the lowest binding energy of -45.14 kcal/mol, followed closely by DehH2-MCA (-41.21 kcal/mol), DehH2-D-2CP (-31.59 kcal/mol), DehH2-3CP (-30.75 kcal/mol), DehH2-2,2- DCP (-29.72 kcal/mol), DehH2-2,3-DCP (-22.20 kcal/mol) and DehH2-TCA (-18.46 kcal/mol). The positive binding energy of the DehH2-chlorpyrifos complex (+180.57 kcal/mol) proved the enzyme's non-preference for the substrate. The results ultimately illustrated the unique specificity of the DehH2 to degrade the above-said pollutants under a hypersaline condition.Communicated by Ramaswamy H. Sarma.

Genomic analysis of a functional haloacid-degrading gene of Bacillus megaterium strain BHS1 isolated from Blue Lake (Mavi Gölü, Turkey)

Purpose

Bacillus megaterium strain BHS1, isolated from an alkaline water sample taken from Mavi Gölü (Blue Lake, Turkey), can grow on minimal medium containing 2,2-dichloropropionic acid. We characterized this bacterium at the genomic level.

Methods

The HiSeq platform was used to carry out genome sequencing, de novo assembly, and scaffolding with strain BHS1. Next, genome data were analyzed to demarcate DNA regions containing protein-coding genes and determine the function of certain BHS1 genes. Finally, results from a colorimetric chloride ion–release assay demonstrated that strain BHS1 produces dehalogenase.

Results

De novo assembly of the BHS1 genomic sequence revealed a genome size of ~ 5.37 Mb with an average G+C content of 38%. The predicted nuclear genome harbors 5509 protein-coding genes, 1353 tRNA genes, 67 rRNA genes, and 6 non-coding (mRNA) genes. Genomic mapping of strain BHS1 revealed its amenability to synthesize two families of dehalogenases (Cof-type haloacid dehalogenase IIB family hydrolase and haloacid dehalogenase type II), suggesting that these enzymes can participate in the catabolism of halogenated organic acids. The mapping identified seven Na+/H+ antiporter subunits that are vital for adaptation of the bacterium to an alkaline environment. Apart from a pairwise analysis to the well-established L-2-haloacid dehalogenases, whole-cell analysis strongly suggested that the haloacid dehalogenase type II might act stereospecifically on L-2-chloropropionic acid, D,L-2-chloropropionic acid, and 2,2-dichloropropionic acid. Whole-cell studies confirmed the utilization of these three substrates and the gene’s role in dehalogenation.

Conclusions

To our knowledge, this is the first report of the full genome sequence for strain BHS1, which enabled the characterization of selected genes having specific metabolic activities and their roles in the biodegradation of halogenated compounds.

Dehalogenase-producing halophiles and their potential role in bioremediation

This review aims to briefly describe the potential role of dehalogenase-producing halophilic bacteria in decontamination of organohalide pollutants. Hypersaline habitats pose challenges to life because of low water activity (water content) and is considered as the largest and ultimate sink for pollutants due to naturally and anthropogenic activities in which a substantial amount of ecological contaminants are organohalides. Several such environments appear to host and support substantial diversity of extremely halophilic and halotolerant bacteria as well as halophilic archaea. Biodegradation of several toxic inorganic and organic compounds in both aerobic and anaerobic conditions are carried out by halophilic microbes. Therefore, remediation of polluted marine/hypersaline environments are the main scorching issues in the field of biotechnology. Although many microbial species are reported as effective pollutants degrader, but little has been isolated from marine/hypersaline environments. Therefore, more novel microbial species with dehalogenase-producing ability are still desired.

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.

Biodegradation of phenolic compounds in high saline wastewater by biofilms adhering on aerated membranes

High saline phenolic wastewater is a typical toxic and refractory industrial wastewater. A single membrane-aerated biofilm reactor (MABR) was used to treat wastewater containing phenol, p-nitrophenol and hydroquinone under increasing phenolic loading and salinity conditions. More than 95 % of phenolic compounds were removed, and a removal efficiency of 8.9 g/m² d for total phenolic (TP) contents was achieved under conditions with 32 g/L of salt and 763 mg/L of influent TP contents. The microbial diversity, structure and function of a biofilm exposed to different conditions were investigated by high-throughput 16S rRNA gene sequencing and metagenomics. Salinity and specific TP loading substantially affected the bacterial community. Gammaproteobacteria, Actinobacteria and Betaproteobacteria contributed more to initial phenolic compound degradation than other classes, with Pseudomonas and Rhodococcus as the main contributing genera. The key phenolic-degrading genes of different metabolic pathways were explored, and their relative abundance was strengthened with increasing phenolic loading and salinity. The diverse cooperation and competition patterns of these microorganisms further promoted the high removal efficiency of multiple phenolic contaminants in the biofilms. These results demonstrate the feasibility of MABR for degrading multiple phenolic compounds in high saline wastewater.

Alternative Bioremediation Agents against Haloacids, Haloacetates and Chlorpyrifos Using Novel Halogen-Degrading Bacterial Isolates from the Hypersaline Lake Tuz

The indiscriminate use of chemical pesticides alongside the expansion of large-scale industries globally can critically jeopardize marine ecology and the well-being of mankind. This is because the agricultural runoffs and industrial effluents eventually enter waterways before flowing into highly saline environments i.e., oceans. Herein, the study assessed two novel bacterial isolates, Bacillus subtilis strain H1 and Bacillus thuringiensis strain H2 from the hypersaline Lake Tuz in Turkey to degrade recalcitrant haloalkanoic acids, haloacetates and chlorpyrifos, and consequently, identify their optimal pollutant concentrations, pH and temperature alongside salt-tolerance thresholds. Bacillus strains H1 and H2 optimally degraded 2,2-dichloropropionic acid (2,2-DCP) under similar incubation conditions (pH 8.0, 30 °C), except the latter preferred a higher concentration of pollutants as well as salinity at 30 mM and 35%, respectively, while strain H1 grew well on 20 mM at <30%. While both isolates could degrade all substrates used, the dehalogenase gene from strain H1 could not be amplified. Capacity of the H2 bacterial isolate to degrade 2,2-DCP was affirmed by the detection of the 795 bp putative halotolerant dehalogenase gene after a successful polymerase chain reaction (PCR) amplification. Hence, the findings envisage the potential of both isolates as bio-degraders of recalcitrant halogenated compounds and those of the same chemical family as chlorpyrifos, in saline environments.

Haloadaptation: insights from comparative modeling studies between halotolerant and non-halotolerant dehalogenases

Halophiles are extremophilic microorganisms that grow optimally at high salt concentrations by producing a myriad of equally halotolerant enzymes. Structural haloadaptation of these enzymes adept to thriving under high-salt environments, though are not fully understood. Herein, the study attempts an in silico investigation to identify and comprehend the evolutionary structural adaptation of a halotolerant dehalogenase, DehHX (GenBank accession number: KR297065) of the halotolerant Pseudomonas halophila, over its non-halotolerant counterpart, DehMX1 (GenBank accession number KY129692) produced by Pseudomonas aeruginosa. GC content of the halotolerant DehHX DNA sequence was distinctively higher (58.9%) than the non-halotolerant dehalogenases (55% average GC). Its acidic residues, Asp and Glu were 8.27% and 12.06%, respectively, compared to an average 5.5% Asp and 7% Glu, in the latter; but lower contents of basic and hydrophobic residues in the DehHX. The secondary structure of DehHX interestingly revealed a lower incidence of α-helix forming regions (29%) and a higher percentage of coils (57%), compared to 49% and 29% in the non-halotolerant homologues, respectively. Simulation models showed the DehHX is stable under a highly saline environment (25% w/v) by adopting a highly negative-charged surface with a concomitant weakly interacting hydrophobic core. The study thus, established that a halotolerant dehalogenase undergoes notable evolutionary structural changes related to GC content over its non-halotolerant counterpart, in order to adapt and thrive under highly saline environments.Communicated by Ramaswamy H. Sarma.

Toxic metals in Perna viridis mussel and surface seawater in Pasir Gudang coastal area, Malaysia, and its health implications

Contamination of toxic metals in P. viridis mussels has been prevalently reported; hence, health risk assessment for consuming this aquaculture product as well as the surrounding surface seawater at its harvesting sites appears relevant. Since Kampung Pasir Puteh, Pasir Gudang is the major harvesting site in Malaysia, and because the last heavy metal assessment was done in 2009, its current status remains unclear. Herein, flame atomic absorption spectrometry and flow injection mercury/hydride system were used to determine the concentrations of Pb, Cd, Cu and total Hg in P. viridis mussels and surface seawater (January-March 2015), respectively. Significantly higher concentrations of these metals were found in P. viridis mussels (p < 0.05) than that of surface seawater samples. The concentrations for Pb (4.27-6.55 μg/g) and Cd (1.55-2.21 μg/g) in P. viridis mussels exceeded the maximum permitted proportion prescribed by the Malaysian law. The concentrations of all metals in surface seawater also violated the Malaysia Marine Water Quality Criteria and Standards. Significant (p < 0.05) and high strength of association (r = 0.787) observed between Pb concentration in P. viridis mussel with the surface seawater indicates its possible application for inferring Pb concentrations in the mussel. Since both the calculated target hazard quotient and hazard index for Pb and Cd exceeded 1, the possible detrimental health impacts on human for consuming P. viridis mussels from this rearing site cannot be ignored. Hence, promoting continuous monitoring programmes and developing efficient toxic metal removal techniques prior to entering the market are required.