eISSN: 2353-9461
ISSN: 0860-7796
BioTechnologia
Current issue Archive About the journal Editorial board Abstracting and indexing Subscription Contact Instructions for authors Publication charge Ethical standards and procedures
Editorial System
Submit your Manuscript
SCImago Journal & Country Rank
2/2023
vol. 104
 
Share:
Share:
RESEARCH PAPERS

Identification of anti-cyanobacterial leads targeting carbonic anhydrase from phytochemical database using in silico approach

Archana Padhiary
1
,
Showkat A. Mir
1
,
Sheary S. Tete
1
,
Iswar Baitharu
2
,
Binata Nayak
1

  1. School of Life Sciences, Sambalpur University, Jyotivihar, Burla, 768019, Odisha, India
  2. Department of Environmental Sciences, Sambalpur University, Jyotivihar, Burla, 768019, Odisha, India
BioTechnologia vol. 104(2) ∙ pp. 121–136 ∙ 2023
Online publish date: 2023/06/26
Article files
Get citation
 
PlumX metrics:
 

Introduction

Cyanobacteria are recognized as the earliest photosynthetic prokaryotes that existed on Earth. They commonly assimilate CO2 from atmospheric gases, which amounts to approximately 3 × 1014 g of C (carbon) (Garcia-Pichel et al., 2003) via), through photosynthetic machinery (Arrigo, 2005). The carbon concentrating mechanism (CCM) is responsible for fixing atmospheric carbon dioxide via carbonic anhydrase (CA), a zinc (Zn) metalloenzyme that is universally present among plants, microbes, and animals (Mondal et al., 2016; DiMario et al., 2017; Aspatwar et al., 2018). CA and inorganic carbon transporters are essential for the CCM in cyanobacteria, which increases the CO2 concentration in the vicinity of Rubisco (Price et al., 2008). Cyanobacteria containing defective CAs displayed abnormal growth (Price et al., 1989). CAs are encoded by three evolutionarily unrelated gene families, α CAs, β CAs, and γCAs. The α CAs are mainly found in vertebrates, plants, algae, and eubacteria; β CAs in bacteria, algae, and fungi; and γ CAs in archaea and bacteria. Additionally, marine diatoms have δ and ζ CAs, while protozoa express η CAs (Chegwidden et al., 2000; Supuran, 2010,2008; Frost and McKenna, 2013; Imtaiyaz et al., 2013). Zn is the metal ion present in the active site of almost all CAs, coordinated by three amino acid residues and one water molecule, which ionizes into an OHion. A hydrogen bond links this OH ion to the amino acids present in the active centre. α-CAs and Γ CAs have HIS as the conserved amino acid, whereas β CAs have CYS2HIS(X) residues arranged in a tetrahedral formation, where X represents an Asp residue or any substitutable ligand (Iverson et al., 2000; Strop et al., 2001; Supuran, 2004; Sawaya et al., 2006). CA stimulates the fixation of CO2, which is carried out by a nucleophilic attack with the hydroxide ion bound to the Zn atom (Sayre, 2010).

The abundant growth of various algae and cyanobacteria in freshwater and marine water bodies frequently poses consequential threats as they are often associated with harmful toxins (Sivonen and Jones, 2009). Cyanobacterial blooms are known to contain hepatotoxins (microcystin andnodularin), cyclic heptapeptides, neurotoxins (anatoxin-a and homoanatoxin-a), and alkaloid cytotoxin (cylindrospermopsin). When the cyanobacteria mature, they lyse and release these toxins into their surroundings (Dittmann and Wiegand, 2006; Sivonen and Jones, 2009). Biofilm-producing cyanobacteria colonize various historical monuments, caves, natural/ artificial stones (lime mortar, carbonate, and sandstone), and pre-historic sites via an adhesion mechanism (Albertano, 2012; Popović et al., 2018). These biofilms are stress-tolerant and can adapt to high temperatures, desiccation, humidity, and various other factors. They also secrete biogenic pigments that cause staining and deterioration of stone surfaces (Albertano, 2012). Thus, microbial toxins significantly reduce the aesthetic value of prehistoric monuments, attributable to their acidolytic and oxido-reductive potential leading to erosion (Gaylarde, 2020). Although pesticides can inhibit cyanobacterial biomass (Nirmal Kumar et al., 2011; Debnath et al., 2012; Ni et al., 2014), continuous usage of these chemicals has harmful eco-toxicological impacts on the environment, along with various health issues and the death of non-targeted microorganisms (Singh et al., 2020). Herbicides such as diuron release toxic substances upon degradation, which are highly persistent and affect aquatic species (Field et al., 2003). Endothall is toxic to cyanobacteria, yet its application as an algaecide is not preferable due to its toxicity on zooplankton (Holdren, 2001; Osano et al., 2002; Giacomazzi and Cochet, 2004).

This study aimed to create a phytochemical database to discover lead molecules that could act as a possible inhibitor, targeting CA without adverse effects on nontargeted inhabitants. CA is the predominant enzyme that helps in increasing the biomass of cyanobacteria, as suggested by various studies (Jiang et al., 2014; Mondal et al., 2016; Tsuzuki et al., 2019). Molecular docking via in silico is a suitable approach to carry out this study (Gupta et al., 2013; Suganya et al., 2017). This method is a convenient tool to predict the binding affinity of the leads to the active site of the desired protein based on strong binding docking scores (Muhammad and Fatima, 2015; Kwatra et al., 2021). Furthermore, density functional theory (DFT) calculations provide higher occupied molecular orbital (HOMO), lower occupied molecular orbital (LUMO) and frontal orbital (FO), which give a preliminary idea about the reactivity and stability of the molecules (Parr et al., 1922). Additionally, the pk CSM web server was used to calculate the absorption, digestion, metabolism, excretion, and toxicity (ADMET) of the compounds, while SwissADME and Molsoft L.L.C online web servers were used to determine the pharmacokinetic properties of small molecules (Molsoft LLC., 2007; Pires et al., 2015).

Materials and methods

Preparation of phytochemical database

In this study, native medicinal plants belonging to Odisha (listed in Supplementary Table S1) were selected to identify lead molecules by targeting the CA (a metalloprotein enzyme) of cyanobacteria. Phytochemicals extracted from these medicinal plants (Supplementary Table S1) were searched in PubChem (https://pubchem.ncbi.nlm.nih.gov/), while those without any previous PubChem history were searched in the scientific literature (Venkatesh et al., 2014; Arunachalam et al., 2015; Panigrahi et al., 2015; Sajini et al., 2019; Jena et al., 2020). The two-dimensional structure of these phytochemicals was drawn using ChemDraw pro-8.0.1 and then converted into a three-dimensional (3D) structure using Chem3D Ultra 8.0. The structures were optimized by adding partial charges using the protonated 3D module, and their energies were minimized using the molecular operating environment (MOE) 09 software (El-Azab et al., 2010; El-Deeb, 2010). The NPACT (naturally occurring plant-based anticancer compound-activity-target) database, which consists of 1 800 molecules, was also incorporated with the phytochemicals present in Odisha. Finally, a phytochemical database containing approximately 3K molecules was created (Supplementary Table S1) by sorting them one by one in the molecular database file using MOE software.

In silico ADMET investigations

The energy-minimized molecules were converted into Simplified Molecular Input Line Entry System (SMILES) format and submitted to online web servers such as pk CSM (http://biosig.unimelb.edu.au/pkcsm/), Swiss-ADME (http://www.swissadme.ch/), and Molsoft L.L.C. (http://molsoft.com/) to predict the ADMET properties of the molecules present in the prepared database. These search engines predicted the pharmacokinetics of the molecules using the bioavailability radar, which considers six physicochemical properties, including lipophilicity, size, polarity, solubility, flexibility, and saturation, to detect toxicity (Molsoft LLC, 2007; Pires et al., 2015). The Lipinski filter, pioneer rule-of-five, was also used in this tool to predict drug likeness (Tripathi et al., 2019). Further studies were conducted on the molecules that fulfilled the above-mentioned criteria.

DFT optimizations

The lead candidates obtained from the ADMET filter were sketched using ChemDraw pro-8.0, and DFT optimizations were performed using the Becke 3-Parameter, Lee-Yang-Parr/Gaussian set (B3YLP/G*) level of theory. The HOMO and LUMO frontal orbitals were determined using the ORCA program 4.1.1. (Neese, 2012, 2017). The band gap in eV provided a preliminary idea about the nature of the molecule, such as its reactivity, hardness, softness, electronegativity, and global electrophilicity. These parameters were calculated using the equations given by Parr (1922).

Preparation of receptor protein

The 3D structure of beta-CA (Synechocystis sp. PCC 6803) was obtained from the online protein database (https://www.rcbs.org) with PDB-ID 5SWC, having a molecular weight of 162.76 kDa, an atom count of 11 602, and a resolution of 1.45 Å. The β-CA is a unique hexameric protein arranged in dimers of trimers in a three-fold symmetry axis. Refinement of this receptor protein was carried out by removing water molecules. The 3D protonation module was used to add nonpolar hydrogen atoms. Further, minimization was conducted by the MMFF94x force field, and the gradient was set at 0.05 using MOE 09 software. The binding site was predicted using ArgusLab 4.1 (Thompson, 2004).

Molecular docking simulation

To analyse the leading inhibitors of the CA enzyme, molecular docking was conducted using non-toxic leads from the 3K phytochemicals obtained from the PubChem database (Supplementary Table S1). The site finder module of MOE09 was used to identify the binding pocket, including the Zn301 atom at the base of the binding site. The docking scores of the ligands to the receptor were calculated using the following computations: rescoring 1: London-dG, placement: triangle matcher, retain: 10, refinement: forcefield, and rescoring 2: London-dG. The ligand molecules with the minimum binding scores were chosen based on their most relevant interaction with the target. The LigX module of MOE was used to identify the interactions with the receptor target.

Results

ADMET filter

The ADMET properties of all the phytochemicals were evaluated, and out of 3K molecules, 50 molecules passed all the criteria, having suitable ADMET profiles, accepting the Lipinski rule of 5, Veber rule, Muegge rule, Ghose rule, and Egan rule. Additionally, these phytochemicals did not exhibit herG-I/II cardiac muscle toxicities, AMES toxicity, blood-brain barrier, hepatotoxicity, and skin sensitization (refer to Fig. 1 and Table 1).

Fig. 1

ADMET analysis of a phytochemical database containing 3K molecules filtered using online search engines

/f/fulltexts/BTA/50640/BTA-104-2-50640-g001_min.jpg
Table 1

ADMET properties of non-toxic molecules obtained from online search engines

MoleculeDrug-likenessMWGI absorptionBBB permeantCNS permeationLipinski vioationsGhose violationsVeber violationsEgan violationsMuegge violationsBioavailabilityAMES toxicityHepatotoxicityHerG I/IIRenal OCT2 substrateSkin sensitization
alpha-tocopherol succinate, DL-0.93530.79lowno−2241120.85nonononono
Methyl alpha-D-glucopyranoside0.12194.18lowno−3.622010020.55nonononono
Prunin0.83434.4lowno−4.053101120.55nonononono
Astragalin0.67448.38lowno−3.908201130.17nonononono
Apigenin-7-glucuronide0.29446.36lowno−3.793201130.11nonononono
Chlorogenic acid0.79354.31lowno−3.856111120.11nonononono
Apigenin0.39270.24highno−2.061000000.55nonononono
2,5-dihydroxybenzoic acid0.30154.12highno−3.283030010.56nonononono
Mycophenolic acid0.69320.34highno−3.013000000.56nonononono
Plumbagic acid0.66224.21highno−2.806000000.56nonononono
Luteolin-7-O-glucoside0.60448.38lowno−3.930201130.17nonononono
Quercetin-3-O-rhamnoside0.82448.38lowno−4.156201130.17nonononono
2-isopropylmalic acid0.30176.17lowno−3.454000010.56nonononono
D-ribonic acid0.60166.13lowno−3.676020020.56nonononono
Methyl galactopyranoside0.12194.18lowno−3.622010020.55nonononono
Galactitol0.10182.17lowno−4.130120030.55nonononono
Galactaric acid0.59210.14lowno−4.224121130.11nonononono
Naproxen0.58230.26highyes−1.951000000.85nonononono
Propionic acid, 3-benzoylamino-3-(4-propoxyphenyl)-0.65327.38highyes−2.632000000.56nonononono
Auriculatin0.63420.46highno−1.763000010.55nonononono
Pseudosolasodine diacetate0.73499.74lowno−1.826120030.55nonononono
3,9-Epoxypregn-16-en-20-one, 3-methoxy-7,11,18-triacetoxy-0.62518.6lowno−3.099211140.17nonononono
Apigenin-7-glucuronide0.29446.36lowno−3.793201130.11nonononono
Mangiferin0.25422.34lowno−4.211211130.17nonononono
Orientin0.59448.38lowno−4.018211130.17nonononono
L-Ascorbyl 2,6-Dipalmitate0.50652.95lowno−2.942241130.56nonononono
Dihydromyricetin0.27320.25lowno−3.586101110.55nonononono
Turanose0.06342.3lowno−4.811211140.17nonononono
Catechin0.64290.27Highno−3.298000000.55nonononono
(−)Cycloxanthochymol0.83602.81lowno−2.168140120.56nonononono
(−)Garcinialiptone A1.17600.8lowno−2.529141120.55nonononono
(−)-kurarinone1.18438.52highno−2.910000010.55nonononono
(−)-Sonderianol0.32284.4highyes−1.789000000.55nonononono
(2S)-2-methoxykurarinone1.05452.55highno−2.424000010.55nonononono
2_4-dihydroxychalcone0.30240.26highyes−2.056000000.55nonononono
2-deoxycucurbitacin D0.42500.68highno−2.809130000.55nonononono
3_5_7-trihydroxyflavonol0.31272.26highno−2.390000000.55nonononono
3-O-Acetylpadmatin0.84360.32highno−3.663000000.55nonononono
4-bromoflavone0.34301.14highyes−1.311000000.55nonononono
5-O-caffeoylshikimic acid0.29336.3lowno−3.738001100.56nonononono
7-OH-flavanone0.74240.26highyes−1.715000000.55nonononono
Maslinic acid0.55472.71highno−1.400130110.56nonononono
Hypophyllanthin0.81430.49highyes−3.144000000.55nonononono
Morin0.46302.24highno−3.407000000.55nonononono
Pongapinone B0.96396.44highyes−2.728000000.55nonononono
Pristimerin0.67466.66lowno−1.020130110.55nonononono
Rocagloic acid0.52476.53highno−3.223000000.56nonononono
Eupafolin0.44316.27highno−3.491000000.55nonononono
Geniposide0.52388.37lowno−4.670011120.11nonononono
Genistein0.44270.24highno−2.066000000.55nonononono
Ursonic acid0.67454.7lowno−1.060130110.85nonononono

Protein structure analysis of CA

The beta-CA of Synechocystis sp. PCC 6803 is arranged in dimeric form as A–B, C–D, and E–F chains (refer to Fig. 2). The alignment tool (BLOSUM40) of the MOE was used to determine the structural identity between the chains present in the X-crystallographic structure. Chain A was set as the reference template, and the remaining chains were superposed onto the reference template. All six chains of CA showed similarities in amino acid position, including loops, helices, alpha, and beta sheets. The superposed structures with the reference template (chain A) showed a pair wise root mean square deviations (RMSD) of each chain below 1Å (refer to Fig. 3). Finally, chains A–B and single chain A were selected for docking.

Fig. 2

The complete set of cyanobacterial carbonic anhydrases consists of six chains, with each chain represented by a distinct color: chain A (red), chain B (green), chain C (blue), chain D (violet), chain E (yellow), and chain F (cyan)

/f/fulltexts/BTA/50640/BTA-104-2-50640-g002_min.jpg
Fig. 3

(A) The reference template sequence was chain A, and the remaining chains B–F were superposed onto this reference; (B) RMSD was recorded

/f/fulltexts/BTA/50640/BTA-104-2-50640-g003_min.jpg

The RMSD of superposed chains B, C, D, E, and F onto the reference chain A were found to be 0.660, 0.692, 0.406, 0.549, and 0.784 Å, respectively.

Molecular docking analysis

Docking simulations of the retrieved 50 non-toxic phytochemicals (Table 1) were carried out at the active site containing the amino acids 103 ALA, 41 ASP, 39 CYS, 101 CYS, 84 ALA, 102 GLY, and 98 HIS along with Zn301 bound to CYS 101, CYS 39, and HIS 98 with the homolog chains A–B and single chain A. Based on the interaction with chain A and chains A–B, the non-toxic phytochemicals with top binding score (S score), drug-likeness score, and Zn interaction status were considered for this study (Table 2 and Table 4).

Table 2

Phytochemicals showing their binding (S-value) and drug-likeness scores with chain A

PhytochemicalPubChem IDS-ValueDrug-likenessInteraction
Mycophenolic acid446541−11.420.69ASP 41, GLY 102, LYS 105, Zn (CYS 39, HIS 98, CYS 101)
Alpha-tocopherol succinate92257186−9.970.93ASP 41, Zn (CYS 39, HIS 98, CYS 101)
Methyl alpha-D-glucopyranoside64947−8.110.12ASP 41, SER 42, GLY 102, Zn (CYS 39, HIS 98, CYS 101)
Apigenin-7-glucuronide5319484−10.380.29ASP 41, GLU 187, ARG 43, GLY 102, Zn (CYS 39, HIS 98, CYS 101)
2,5-Dihydroxybenzoic acid3469−9.570.30ASP 41, Zn (CYS 39, HIS 98, CYS 101)
Orientin5281675−12.010.59ASP 41, GLY 102, LYS 105, LYS 109, Zn (CYS39. HIS 98, CYS 101)
Methyl galactopyranoside76935−8.090.12ASP 41, SER 42, GLY 102, Zn (CYS 39, HIS 98, CYS 101)
Galactitol11850−8.820.10ASP 41, SER 42, GLY 102, Zn (CYS 39, HIS 98, CYS 101)

With single-chain A, most of the phytochemicals exhibited an interaction with Zn301 (CYS 101, CYS 39, and HIS 98). Ligands which demonstrated interaction with Zn301, also interacted with ASP 41 and other amino acids. Alpha-tocopherol succinate and 2,5-dihydroxybenzoic acid only interacted with ASP 41. Similarly, methyl alpha-D-glucopyranoside and methyl galactopyranoside demonstrated similar interaction with ASP 41, SER 42, and GLY 102 due to their structural similarity. The interaction of apigenin-7-glucuronide was observed with ASP 41, GLU 187, ARG 43, and GLY 102. Similarly, Orientin exhibited an interaction with ASP 41, GLY 102, LYS 105, and LYS 109, while dihydromyricetin interacted with ASP 41 and LYS 105. Galactitol showed interaction with ASP 41, SER 42, and GLY 102. Table 2 and Figure 4 depict the interactions, and Table 3 provides information on the distances of the Zn-interacted molecules.

Table 3

The interactions and distances displayed by leads with zinc and other amino acids with chain A

PhytochemicalLigandReceptorResidueChainTypeScore [%]Distance
Alpha-tocopherol succinateH 3584OD 707ASP 415swc 1H-don90.11.24
Zn 3496SG 684CYS 395swc 1ionic67.62.29
Zn 3496NE 1533HIS 985swc 1ionic78.62.06
Zn 3496SG 1574CYS 1015swc 1ionic71.72.31
O 3583ZN 3496ZN 3015swc 2ionic100.02.28
Methyl alpha-D-glucopyranosideH 3517OD 707ASP 415swc 1H-don86.71.57
O 3501N 712SER 425swc 1H-acc15.63.29
O 3516N 1579GLY 1025swc 1H-acc33.02.96
Zn 3496SG 684CYS 395swc 1ionic67.62.29
Zn 3496NE 1533HIS 985swc 1ionic78.62.06
Zn 3496SG 1574CYS 1015swc 1ionic71.72.31
O 3516ZN 3496ZN 3015swc 2ionic100.02.22
OrientinH 3538OD 707ASP 415swc 1H-don81.51.67
H 3546O 1582GLY 1025swc 1H-don40.41.68
O 3537N 1579GLY 1025swc 1H-acc34.22.98
O 3510NZ 1621LYS 1055swc 1H-acc35.52.45
O 3547NZ 1688LYS 1095swc 1H-acc87.02.70
Zn 3496SG 684CYS 395swc 1ionic67.62.29
Zn 3496NE 1533HIS 985swc 1ionic78.62.06
Zn 3496SG 1574CYS 1015swc 1ionic71.42.31
O 3537ZN 3496ZN 3015swc 2ionic100.02.36
Mycophenolic acidH 3537OD 707ASP 415swc 1H-don86.01.23
O 3536N 1579GLY 1025swc 1H-acc11.52.99
O 3508NZ 1621LYS 1055swc 1H-acc55.92.43
O 3538NZ 1621LYS 1055swc 1H-acc54.62.87
Zn 3496SG 684CYS 395swc 1Ionic67.62.29
Zn 3496NE 1533HIS 985swc 1Ionic78.62.06
Zn 3496SG 1574CYS 1015swc 1Ionic71.42.31
O 3536ZN 3496ZN 3015swc 2ionic100.02.40
Dihydro-myricetinH 3510OD 707ASP 415swc 1H-don73.61.31
O 3516N 1579LYS 1055swc 1H-acc38.82.48
O 3528NZ 1621LYS 1055swc 1H-acc24.83.02
Zn 3496NZ 1621CYS 395swc 1ionic67.62.29
Zn 3496SG 684HIS 985swc 1ionic78.62.06
Zn 3496NE 1533CYS 1015swc 1ionic71.42.31
O 3509SG 1574ZN 3015swc 2ionic100.02.48
Methyl galacto-pyranosideH 3517OD 707ASP 415swc 1H-don86.81.56
O 3501N 712SER 425swc 1H-acc16.13.28
O 3516N 1579GLY 1025swc 1H-acc31.92.96
Zn 3496SG 684CYS 395swc 1ionic67.62.29
Zn 3496NE 1533HIS 985swc 1ionic78.62.06
Zn 3496SG 15794CYS 1015swc 1ionic71.42.31
O 3526ZN 3496ZN 3015swc 2ionic100.02.22
GalactitolH 3522OD 707ASP 415swc 1H-don88.21.55
O 3511OG 717SER 425swc 1H-don19.92.98
O 3511N 712SER 425swc 1H-acc72.52.70
O 3511OG 717SER 425swc 1H-acc19.92.79
O 3519N 1579GLY 1025swc 1H-acc33.32.75
O 3521N 1579GLY 1025swc 1H-acc16.53.07
Zn 3496SG 684CYS 395swc 1Ionic67.62.29
Zn 3496NE 1533HIS 985swc 1ionic78.62.06
Zn 3496SG 15794CYS 1015swc 1ionic71.42.31
O 3521ZN 3496ZN 3015swc 2ionic100.02.07
Fig. 4

The figure shows the molecular interaction between zinc ion in the binding pocket of beta-carbonic anhydrase chain A and the following compounds: (A) Alpha-tocopherol succinate, (B) Methyl alpha-D-glucopyranoside, (C) Apigenin-7-glucuronide, (D) 2,5-Dihydroxybenzoic acid, (E) Orientin, (F) Mycophenolic acid, (G) Dihydromyricetin, and (H) Galactitol

/f/fulltexts/BTA/50640/BTA-104-2-50640-g004_min.jpg

In the active site of chains, A–B, the ligands that demonstrated interaction with Zn301 were alpha-tocopherol succinate and mycophenolic acid through CYS 101, HIS 98, and CYS 39. Alpha-tocopherol succinate interacted with GLY102 of chain A and GLN 30 of chain B, while mycophenolic acid interacted with both ASP 41 and LYS 105 of chain A. Non zinc ligands interacted with amino acids in chain A (ASP A41, ARG A43, LYS A105, LYS A109, HIS A100, and GLY A102) and chain B (GLN B30, TYR B83, and ALA B87) as shown in Table 4 and Figure 5. Table 5 provides information on the distance score % of the Zn-interacted molecules.

Table 4

Phytochemicals showing their binding scores (S-value) and drug-likeness with chain A–B

PhytochemicalPubChem IDS-ValueDrug-likenessInteraction
Alpha-tocopherol succinate92257186−9.230.93GLY A102, GLN B30, Zn (CYS 101, HIS 98, CYS 39)
Hypophyllanthin182140−11.060.81ARG A43, LYS A105, LYS A109
Pseudosolasodine diacetate313026−11.670.73LYS A105, HIS A100, ARG A43 GLN B30
Methyl galactopyranoside76935−12.630.12LYS A105, ARG A43, GLY A102, HIS A100, TYR B83, GLN B30, ASP A41
2,5 dihydroxybenzoic acid3469−12.470.30TYR B83, GLN B30, HIS A100, ARG A43, HIS B28
Dihydromyricetin161557−13.580.27LYS A109, LYS A105, GLN B30
Mycophenolic acid446541−14.410.69ASP A41, LYS A105, Zn (CYS 101, HIS 98, CYS 39)
3,9-Epoxypregn-16-en-20-one, 3-methoxy-7,11,18-triacetoxy-539331−15.590.62TYR B83, ALA B87, HIS A100, GLN B30, ARG A43, TYR B83, LYS A105, LYS A109
Table 5

The interactions and distances displayed by leads with zinc and other amino acids with beta-carbonic anhydrase chain A–B

PhytochemicalLigandReceptorResidueChainTypeScore [%]Distance
Alpha-tocopherol succinateO 7070N 4014GLN 305swc 2H-acc14.92.89
O7079N 1579GLY 1025swc 1H-acc18.52.68
Zn 6992SG 684CYS 395swc 1ionic67.62.29
Zn 6992NE 1533HIS 985swc 1ionic78.62.06
Zn 6992SG 1574CYS 1015swc 1ionic71.42.31
O 7079ZN 6992ZN 3015swc 3ionic100.02.54
Mycophenolic acidH 7033OD 707ASP 415swc 1H-don84.71.20
O 7004NZ 1621LYS 1055swc 1H-acc68.92.46
O 7034NZ 1621LYS 1055swc 1H-acc62.12.54
Zn 6992SG 684CYS 395swc 1ionic67.62.29
Zn 6992NE 1533HIS 985swc 1ionic78.62.06
Zn 6992SG 1574CYS 1015swc 1ionic71.42.31
O 7031ZN 6992ZN 3015swc 3ionic100.02.62
O 7032ZN 6992ZN 3015swc 3ionic100.02.20
Fig. 5

The interaction shown by (A) Alpha-tocopherol succinate and (B) Mycophenolic acid with zinc ion in the binding pocket of beta-carbonic anhydrase chains A–B

/f/fulltexts/BTA/50640/BTA-104-2-50640-g005_min.jpg

Molecular docking simulations revealed that alphatocopherol succinate and mycophenolic acid interacted with Zn301 and neighboring amino acids present in the binding site. These findings suggest that alpha-tocopherol succinate and mycophenolic acid could be effective phytochemicals that inhibit the cellular functioning of CAs found in cyanobacteria.

Global electrophilicity parameters of lead phytochemicals

The frontier molecular orbitals provide data on the chemical reactivity and physical parameters of molecules (Sheikhshoaie et al., 2014). The HOMO acts as an electron donor, and the LUMO acts as an electron acceptor. The energy gap between HOMO and LUMO determines the stability of structures, including their chemical reactivity and kinetic stability (Kim et al., 2013). The computed data on the HOMO, LUMO, energy gap (ΔE), ionization potential, electron affinity (A), chemical potential (μ), electron negativity (χ), hardness (η), softness (σ), and electrophilicity (ω) of the molecules are presented in Table 6 and Figure 6.

Table 6

Global electrophilicity parameters of alpha-tocopherol succinate, methyl alpha-D-glucopyranoside, apigenin-7-glucuronide, mycophenolic acid, 2,5-dihydroxybenzoic acid, orientin, dihydromyricetin, methyl galactopyranoside, and galactitol

Phyto chemicalsHOMOLUMOEnergy gapIonization potentialElectron affinityChemical potentialElectron negativityHardnessSoftnessElectrophilicity
Alpha-tocopherol succinate−5.83−0.5685.2625.820.568−3.1983.1982.6310.3801.948
Methyl alpha-D-glucopyranoside−7.12−0.0257.0947.11-0.02−3.5723.5723.5470.2811.798
Apigenin-7-Glucuronide−6.091−1.9694.1226.0911.968−4.0304.0302.0610.4823.940
Mycophenolic acid−5.99−1.284.7105.9891.279−3.6343.6342.3550.4242.805
2,5-Dihydroxy-benzoic acid−6.05−1.7784.2726.0501.778−3.9143.9142.1360.4983.586
Orientin−6.041−2.0114.0296.0402.010−4.0264.2062.0140.4624.022
Dihydromyricetin−6.258−1.8374.4216.2571.837−4.0474.0472.2100.4523.705
Methyl galactopyranoside−7.12−0.0257.0947.11−0.02−3.5723.5723.5470.2811.798
Galactitol−7.3140.0227.3357.3140.022−3.6463.6463.6670.2721.812
Fig. 6

The energy shifting of HOMO and LUMO due to the occupancy of atoms present in the molecular scaffold of the following compounds: (A) Alpha-tocopherol succinate (B) Methyl alpha D-glucopyranoside (C) Apigenin-7-Glucuronide, (D) Mycophenolic acid (E) 2,5-Dihydroxybenzoic acid (F) Orientin (G) Dihydromyricetin, and (H) Galactitol

/f/fulltexts/BTA/50640/BTA-104-2-50640-g006_min.jpg

In alpha-tocopherol succinate, the energy shift was assigned to the succinate moiety, in LUMO to the tocopherol ring in HOMO. In apigenin-7-O-glucuronide, the energy shift was assigned to the apigenin moiety for both HOMO and LUMO. In mycophenolic acid, the energy shift was assigned to the phthalide intermediate in LUMO to 2-(but-2-enyl)-3-methoxy-methylphenol in HUMO. In 2,5-dihydroxy benzoic acid, the energy shift LUMO was assigned to the benzoic acid and HOMO to the hydroquinone moiety. In orientin, both LUMO and HOMO energy shift was assigned to the 2-(3,4-di-hydroxyphenyl)-5,7-dihydroxy-8-methyl-4H-chromen-4-one moiety. In dihydromyricetin, the energy shift in LUMO was assigned to the 3,5,7-trihydroxy-2,3-dihydrochrome-4-one moiety, and HOMO was assigned to benzene 1,2,3-triol. In galactitol, LUMO energy shift covered a vast surface of the galactitol molecule, whereas HOMO distribution was symmetrical. In methyl alpha-D-glucopyranoside and methyl galactopyranoside, LUMO was assigned to the methyl position, and HOMO was assigned to the overall molecule, as shown in Figure 6. Table 6 presents the band gaps of alpha-tocopherol succinate, apigenin-7-O-glucuronide, mycophenolic acid, 2-5-dihydroxybenzoic acid, orientin, dihydromyricetin, and galactitol as 5.262, 4.122,4.710, 4.272, 4.029, 4.421, and 7.335 eV, respectively, and the electrophilicity indices as 1.948, 3.940, 2.805,3.586, 4.022, 3.705, and 1.812 eV, respectively. Due to their structural similarity, methyl alpha-D-glucopyranoside and methyl galactopyranoside had an energy gap of 7.094 eV and an electrophilicity index of 1.798 eV. Overall, it was observed that all the molecules had high energy band gaps and a low global softness index. This suggests that the molecules are effective due to their stability and less reactivity.

Discussion

Cyanobacteria are primitive organisms that exist in various habitats, from hot springs to cold arctic regions. These bodies are gradually becoming eutrophied due to anthropogenic and industrial activities. Underdeveloped nations with limited technology sources rely directly on open water systems, lakes, ponds, rivers, oceans, and sanitary households for their day-to-day work. These water systems often contain cyanobacteria biofilms that secrete toxins, thereby deteriorating water quality on a large scale and leading to major health issues (O’Neil et al., 2012; Paerl et al., 2013). Therefore, measures should be taken to suppress harmful cyanobacterial growth. Naturally occurring nontoxic phytochemicals are generally secondary metabolites (terpenes, flavonoids, flavones, alkaloids, and polyphenols) that carry out several metabolic activities related to defense, stress tolerance, and interactions (Leitzmann, 2016). They are known as “nutraceuticals” due to their physiological benefits, such as being anticancer, antioxidants, anti-inflammatory, antiviral, and anti-bacterial agents (Brindha, 2016; Winter, 2017; Velmurugan et al., 2018; Ranjan et al., 2019). Natural compounds produced from plants, such as stilbenes, nostocarboline, berberine, and phenylpropanoid glucosides, have been reported as potent algaecides (Becher et al., 2005; Mizuno et al., 2008; Jančula et al., 2010; Zuo et al., 2018). Similarly, phytochemicals isolated from plants like Pistacia integerrima, brown alga Dictyopteris hoytii, flavonoids from Luffa acutangular (L.) Roxb, Murraya paniculata, and Padina pavonica,have potent CA inhibitory activity (Sarikaya et al., 2011; Sahin et al., 2012; Karioti et al., 2016; Chanda et al., 2019; Sangkaew et al., 2020; Irfan et al., 2021; Rafiq et al., 2021).

As the application of synthetic compounds to algal blooms can adversely affect aquatic organisms, we have chosen nontoxic phytochemical leads for this study that can inhibit blooms by targeting cyanobacterial beta-CAs without affecting non-targeted organisms.

A phytochemical database of 3K molecules was created, considering the medicinal plants of Odisha via a previously reported literature survey. The toxicity of phytochemicals in the database was investigated, and non-toxic phytochemicals were determined using multi software packages (Molsoft LLC., 2007; Pires et al., 2015). Fifty molecules (Table 1) from 3K molecules followed drug likenesses and ADMET properties (Pires et al., 2015; Tripathi et al., 2019). Among them, the two lead molecules, alpha-tocopherol succinate and mycophenolic acid, showed drug-likeness scores of 0.93 and 0.69, respectively, out of 1 (Molsoft LLC., 2007). The stability and reactivity of the lead molecules were determined by DFT calculations. The HOMO–LUMO energy gap was identified using B3YLP/G* level of theory, from which the ionization potential, electron affinity, chemical potential, electron negativity, hardness, softness, and electrophilicity were calculated. Hard molecules possess large energy gaps and are highly stable, whereas soft molecules have a smaller energy gap and are highly polarizable (Pilli et al., 2015). The energy gap, hardness, softness, and electrophilicity parameters of alpha-tocopherol succinate were 5.262, 2.631, 0.380, and 1.948 eV, respectively, and for mycophenolic acid, they were 4.710, 2.355, 0.424, and 2.805 eV, respectively, which supports the findings by Pilli et al. (2015).

These findings support a previous study demonstrating that phytochemicals play a pivotal role in inhibiting the biomass of cyanobacterial species, as demonstrated by both in silico and in vitro studies (Mir and Nayak, 2022).

The beta-CA of cyanobacteria Synechocystis sp. PCC 6803, with the cofactor Zn (Zn301) present in the active site, was chosen as the target protein. Zn functions as a catalyst and regulates enzymatic reactions (Escudero-Almanza et al., 2012). Initially, single chain A and later double chains A–B of beta-CA were used for docking to identify interactions with Zn2+. Alpha-tocopherol succinate and mycophenolic acid were found to be potent phytochemical leads as they exhibited interactions with Zn2+ and its coordinate amino acids CYS 101, CYS 39, and HIS 98 in both single and double chain. Previous reports have shown that phytochemicals such as pistacide-A and pistacide-B obtained from P. integerrima and lacceroic acid from brown alga D. hoytii is potent inhibitors of CA-II (Irfan et al., 2021; Rafiq et al., 2021). These inhibitors showed binding scores of −3.45, 3.64, and −9.53 Kcal/mol, respectively, against the reference molecule acetazolamide (−9.47 Kcal/mol) (Irfan et al., 2021; Rafiq et al., 2021). Similarly, triazine substituted sulfamerazine-derived compound 5 (4-ethoxy) and compound 8 (4-bromo) have also been identified as potent human CA I and II inhibitors, respectively (Bilginer et al., 2021). Compound 5 had a binding score of −7.36 Kcal/mol with human CAI (hCAI) and compound 8 (−8.26 Kcal/mol) with hCAII (Bilginer et al., 2021), respectively.

In this study, the identified most potent phytochemicals, alpha-tocopherol succinate, and mycophenolic acid, possessed binding scores of −9.23 and −14.41 Kcal/mol, respectively, which were lower than those of the reported CA II inhibitors (pistacide-A, pistacide-B, and lacceroic acid) and hCAI/hCAII inhibitors (compounds 5 and 8). This confirms the strong interaction of the lead phytochemicals with the active site of beta-CAs.

Several phytochemicals have been reported to act as anti-carbonic anhydrase agents against CA-II, CA-V, hCA-IX, and hCA-XII, which are associated with abnormalities such as tumor progression, cancer, obesity, and Alzheimer’s disease (Rahman et al., 2015; Queen et al., 2018; Fois et al., 2020; Mahmud et al., 2021). However, to our knowledge, no studies have been conducted targeting cyanobacterial beta-CA. Therefore, in this study, we have considered phytochemicals that showed robust bindings with the amino acids and co-factor present in the binding site of beta-CAs. These phytochemicals have the potential to inhibit the cellular functioning of CA, thus restricting the allocation of bicarbonates to the RuBisCo site and curtailing the formation of cyanobacterial biomass. Further studies such as molecular dynamics simulations and in-vitro experiments will be carried out to validate the present in-silico findings.

Conclusion

In rural areas and economically disadvantaged countries, people rely on open-water systems for their daily needs. Unfortunately, anthropogenic activities, such as the disposal of leached micronutrients and industrial effluents into aquatic environments, lead to cyanobacterial blooms. These harmful cyanobacteria release cyanotoxins into the water, which can cause major health issues like hepatotoxicity and immunotoxicity when ingested orally. As such, it is important to identify novel anti-cyanobacterial agents that do not have any negative side effects on aquatic organisms. In this study, alphatocopherol succinate and mycophenolic acid phytochemicals were found to exhibit binding interactions with the co-factor Zn in both single chain A and complexed chains A–B. These lead molecules interact with the Zn-binding site present in the enzyme and may inhibit the functioning of the beta-CA complex, thereby restricting biomass generation. Considering the present and future perspectives of water scarcity, this approach may be a possible way to protect the biotic environment, water bodies, and ancestral monuments. Further in vitro experiments on lead phytochemicals are necessary to validate the present in silico study.

Acknowledgments

The author SAM would like to acknowledge Ibrahim Khalifa from Benha University, Egypt, for providing us to access MOE tools. The authors SAM and BN also acknowledge partial financial support from the DST Government of India, under SHRI grant.

Conflict of interest

The authors declare no conflict of interest

References

1 

Arrigo K.R. (2005) Marine microorganisms and global nutrient cycles. Nature 437(7057): 349–355.

2 

Arunachalam K.D., Kuruva J.K., Hari S., Annamalai S.K., Baskaran K.V. (2015) HPTLC Finger print analysis and phytochemical investigation of Morinda tinctoria Roxb. Leaf extracts by HPLC and GS MS. Int. J. Pharm. Pharm. Sci. 7(2): 360–366.

3 

Aspatwar A., Haapanen S., Parkkila S. (2018) An update on the metabolic roles of carbonic anhydrases in the model alga Chlamydomonas reinhardtii. Metabolites 8(1): 22.

4 

Becher P.G., Beuchat J., Gademann K., Jüttner F. (2005) Nostocarboline: isolation and synthesis of a new cholinesterase inhibitor from Nostoc 78-12A. J. Nat. Prod. 68(12): 1793–1795.

5 

Brindha P. (2016) Role of phytochemicals as immunomodulatory agents: a review. Inter. J. Green Pharm. 10(1). https://doi.org/10.22377/ijgp.v10i1.600

6 

Chanda J., Mukherjee P.K., Biswas R., Biswas, S., Tiwari A.K., Pargaonkar A. (2019) UPLC QTOF MS analysis of a carbonic anhydrase inhibiting extract and fractions of Luffa acutangula (L.) Roxb (ridge gourd). Phytochem. Anal. 30: 148–155.

7 

Chegwidden W.R., Carter N.D., Edwards Y.H. (2000) The carbonic anhydrases: new horizons. Birkhauser Spring Sci. Business Media.

8 

Debnath M., Mandal N.C., Ray S. (2012) Effect of fungicides and insecticides on growth and enzyme activity of four cyanobacteria. Ind. J. Micro. 52(2): 275–280.

9 

DiMario R.J., Clayton H., Mukherjee A., Ludwig M., Moroney J.V. (2017) Plant carbonic anhydrases: structures, locations, evolution, and physiological roles. Mol. Plant. 10(1): 30–46.

10 

Dittmann E., Wiegand C. (2006) Cyanobacterial toxins – occurrence, biosynthesis and impact on human affairs. Mol. Nutr. Food Res. 50(1): 7–17.

11 

El-Azab A.S., Al-Omar M.A., Alaa A.M., Abdel-Aziz N.I., Magda A.A., Aleisa A.M., Abdel-Hamide S.G. (2010) Design, synthesis and biological evaluation of novel quinazoline derivatives as potential antitumor agents: molecular docking study. Eur.J. Med. Chem. 45(9): 4188–4198.

12 

El-Deeb I.M., Bayoumi S.M., El-Sherbeny M.A., Alaa A.M. (2010) Synthesis and antitumor evaluation of novel cyclic aryl sulfonyl ureas: ADME-T and pharmacophore prediction. Eur. J. Med. Chem. 45(6): 2516–2530.

13 

Escudero-Almanza D.J., Ojeda-Barrios D.L., Hernández-Rodríguez O.A., Sánchez Chávez E., Ruíz-Anchondo T., Sida-Arreola J.P. (2012) Carbonic anhydrase and zinc in plant physiology. Chile J. Agr. Res. 72(1): 140–146.

14 

Field J.A., Reed R.L., Sawyer T.E., Griffith S.M., Wigington Jr. P.J. (2003) Diuron occurrence and distribution in soil and surface and ground water associated with grass seed production. J. Env. Qual. 32(1): 171–179.

15 

Fois B., Distinto S., Meleddu R., Deplano S., Maccioni E., Floris C., Cottiglia F. et al. (2020) Coumarins from Magydaris pastinacea as inhibitors of the tumour-associated carbonic anhydrases IX and XII: isolation, biological studies and in silico evaluation. J. Enz. Inh. Med. Chem. 35(1): 539–548.

16 

Frost S.C., Mckenna R. (2013) Carbonic anhydrase: mechanism, regulation, links to disease, and industrial applications. Dordrecht Spring Sci. & Business Media.

17 

Garcia-Pichel F., Belnap J., Neuer S., Schanz F. (2003) Estimates of global cyanobacterial biomass and its distribution. Algal. Stud. 109: 213–227.

18 

Gaylarde C. (2020) Influence of environment on microbial colonization of historic stone buildings with emphasis on cyanobacteria. Heritage 3(4): 1469–1482.

19 

Giacomazzi S., Cochet N. (2004) Environmental impact of diuron transformation: a review. Chemosphere 56(11): 1021–1032.

20 

Gupta S., Kapoor P., Chaudhary K., Gautam A., Kumar R., Open-Source Drug Discovery Consortium, Raghava G.P. (2013) In silico approach for predicting toxicity of peptides and proteins. PloS one 8(9): e73957.

21 

Hewett-Emmett D., Tashian R.E. (1996) Functional diversity, conservation, and convergence in the evolution of the a-, b-, and g-carbonic anhydrase gene families. Mol. Phyl. Evol. 5: 50–77.

22 

Holdren C., Jones W., Taggart J. (2001) Managing lakes and reservoirs. North American Lake Management Society, Madison, Wisconsin.

23 

Imtaiyaz Hassan M., Shajee B., Waheed A. et al. (2013) Structure, function and applications of carbonic anhydrase isozymes. Bioorg. Med. Chem. 21(6): 1570–1582.

24 

Irfan A., Imran M., Sumrra S.H., Qaisar M.N., Khalid N., Basra M.A.R., Al-Sehemi A.G. et al. (2021) Exploration of carbonic anhydrase inhibition of bioactive metabolites from Pistacia integerrima by molecular docking and first-principles investigations. J. Saudi Chem. Soc. 25(10): 101324.

25 

Iverson T.M., Alber B.E., Kisker, C., Ferry J.G., Rees D.C. (2000) A closer look at the active site of gamma class carbonic anhydrases: high resolution crystallographic studies of the carbonic anhydrase from Methanosarcina thermophila. Biochemistry 39: 9222–9231.

26 

Jančula D., Gregorová J., Maršálek B. (2010) Algicidal and cyanocidal effects of selected isoquinoline alkaloids. Aqua. Res. 41(4): 598–601.

27 

Jena R., Rath D., Rout S.S., Kar D.M. (2020) A review on genus Millettia: traditional uses, phytochemicals and pharmacological activities. Saudi Pharm. J. 28(12): 1686.

28 

Jiang C., Tholen D., Xu J.M., Xin C., Zhang H., Zhu X., Zhao Y. (2014) Increased expression of mitochondria-localized carbonic anhydrase activity resulted in an increased biomass accumulation in Arabidopsis thaliana. J. Plant Biol. 57(6): 366–374.

29 

Karioti A., Carta F., Supuran C.T. (2016) Phenols and polyphenols as carbonic anhydrase inhibitors. Molecules 21(12): 1649.

30 

Kim B.G., Ma X., Chen C., Ie Y., Coir E.W., Hashemi H., Kim J. et al. (2013) Energy level modulation of HOMO, LUMO, and band gap in conjugated polymers for organic photovoltaic applications. Adv. Funct. Mat. 23(4): 439–445.

31 

Kwatra B., Khatun A., Bhowmik R., Rehman S. (2021) In silico-modelling of phytochemicals in septic arthritis. Pharma. Innov. 10(3): 14–21.

32 

Leitzmann C. (2016) Characteristics and health benefits of phytochemicals. Compl. Med. Res. 23(2): 69–74.

33 

Mahmud S., Rahman E., Nain, Z., Billah M., Karmakar S., Mohanto S.C., Saleh M.A. et al. (2021) Computational discovery of plant-based inhibitors against human carbonic anhydrase IX and molecular dynamics simulation. J. Biomol. Str. Dyn. 39(8): 2754–2770.

34 

Mir S.A., Nayak B. (2022) In silico analysis of binding stability of quercetin with CmpA and in vitro growth inhibition study of Cyanobacterial species using Azadirachta indica extracts. Chem. Africa: 1–11.

35 

Mizuno C.S., Schrader K.K.,Rimando A.M. (2008) Algicidal activity of stilbene analogues. J. Agr. Food Chem. 56(19): 9140–9145.

36 

Molsoft LLC (2007) 3366, North Torrey Pines Court, Suite 300, La Jolla, CA 92037, USA. https://www.molsoft.com/

37 

Mondal M., Khanra S., Tiwari O.N., Gayen K., Halder G.N. (2016) Role of carbonic anhydrase on the way to biological carbon capture through microalgae — a mini review. Env. Prog. Sust. Ener. 35(6): 1605–1615.

38 

Muhammad S.A., Fatima N. (2015) In silico analysis and molecular docking studies of potential angiotensin-converting enzyme inhibitor using quercetin glycosides. Pharma. Magazine 11(42): 123.

39 

Neese F. (2017) Software update: the ORCA program system, ver. 4.0. Wiley Interdiscip. Rev. Comp. Mol. Sci. 8: e1327.

40 

Neese F. (2012) The ORCA program system. Wiley Interdiscip. Rev. Comp. Mol. Sci. 2(1): 73–78.

41 

Ni Y., Lai J., Wan J., Chen L., (2014) Photosynthetic responses and accumulation of mesotrione in two freshwater algae. Env. Sci. Proc. Imp.16: 2288–2294.

42 

Nirmal Kumar J.I., Bora A., Amb M.K., Kumar R.N. (2011) An evaluation of pesticide stress induced proteins in three cyanobacterial species Anabaena fertilissima, Aulosira fertilissima and Westiellopsis prolifica using SDS-PAGE. Adv. Env. Biol. 5: 739–745.

43 

O’Neil J.M., Davis T.W., Burford M.A.,Gobler C.J. (2012) The rise of harmful cyanobacteria blooms: the potential roles of eutrophication and climate change. Harm. Algae 14: 313–334.

44 

Osano O., Admiraal W., Klamer H.J.C., Pastor D., Bleeker E.A.J. (2002) Comparative toxic and genotoxic effects of chloroacetanilides, formamidines and their degradation products on Vibrio fischeri and Chironomus riparius. Env. Pol. 119(2): 195–202.

45 

Paerl H.W., Otten T.G. (2013) Harmful cyanobacterial blooms: causes, consequences, and controls. Micro. Eco. 65(4): 995–1010.

46 

Panigrahi G.K., Ch R., Mudiam M.K., Vashishtha V.M., Raisuddin S., Das M. (2015) Activity-guided chemo toxic profiling of Cassia occidentalis (CO) seeds: detection of toxic compounds in body fluids of CO-exposed patients and experimental rats. Chem. Res. Toxic. 28(6): 1120–1132.

47 

Parr R.G., Szentpaly L.V., Liu S. (1922) Electrophilicity Index. J. Am. Chem. Soc. 121: 1922–1924.

48 

Pilli S.R., Banerjee T., Mohanty K. (2015) HOMO–LUMO energy interactions between endocrine disrupting chemicals and ionic liquids using the density functional theory: evaluation and comparison. J. Mol. Liq. 207: 112–124.

49 

Pires D.E., Blundell T.L., Ascher D.B. (2015) pkCSM: predicting small-molecule pharmacokinetic and toxicity properties using graph-based signatures. J. Med. Chem. 58(9): 4066–4072.

50 

Popović S., Stupar M., Unković N., Subakov Simić, G., Ljaljević Grbić M. (2018) Diversity of terrestrial cyanobacteria colonizing selected stone monuments in Serbia. Stud. Cons. 63(5): 292–302.

51 

Price G.D., Badger M.R. (1989) Isolation and characterization of high CO2-requiring-mutants of the cyanobacterium Synechococcus PCC7942. Plant Physiol. 91: 514–525.

52 

Price G.D., Badger M.R., Woodger F.J., Long B.M. (2008) Advances in understanding the cyano–bacterial CO2_concentrating mechanism (CCM): functional components, Ci transporters, diversity, genetic regulation and prospects for engineering into plants. J. Exp. Bot. 59: 1441–1461.

53 

Queen A., Khan, P., Azam A., Hassan M.I. (2018) Understanding the role and mechanism of carbonic anhydrase V in obesity and its therapeutic implications. Curr. Protein Peptide Sci. 19(9): 909–923.

54 

Rafiq K., Khan A., Ur Rehman N., Halim S.A., Khan M., Ali L., Al-Harrasi A. (2021) New carbonic anhydrase-II inhibitors from marine macro brown alga Dictyopteris hoytii supported by in silico studies. Molecules 26(23): 7074.

55 

Rahman T.U., Khattak K.F., Liaqat W., Zaman K., Musharraf S.G. (2015) Characterization of one novel flavone and four new source compounds from the bark of Millettia ovalifolia and in-vitro inhibition of carbonic anhydrase-II by the novel flavonoid. Rec. Nat. Prod. 9(4): 553.

56 

Ranjan A., Ramachandran S., Gupta N., Kaushik I., Wright S., Srivastava S., Srivastava S.K. et al. (2019) Role of phytochemicals in cancer prevention. Int. J. Mol. Sci. 20(20): 4981.

57 

Sahin H., Can Z., Yildiz O., Kolayli S., Innocenti A., Scozzafava G., Supuran C.T. (2012) Inhibition of carbonic anhydrase isozymes I and II with natural products extracted from plants, mushrooms and honey. J. Enzyme Inhib. Med. Chem. 27(3): 395–402.

58 

Sajini R., Prema S., Chitra K. (2019) Phytoconstituents, pharmacological activities of Marsilea minuta l. (Marsileaceae) – an overview. Int. J. Pharm. Sci. Res. 10(4): 1582–1587.

59 

Sangkaew A., Samritsakulchai N., Sanachai K., Rungrotmongkol T., Chavasiri W., Yompakdee C. (2020) Two flavonoid-based compounds from Murraya paniculata as novel human carbonic anhydrase isozyme II inhibitors detected by a resazurin yeast-based assay. J. Microbiol. Biotechnol. 30(4): 552–560.

60 

Sarikaya S.B.Ö., Topal F., Şentürk M., Gülçin I., Supuran C.T. (2011) In vitro inhibition of α-carbonic anhydrase isozymes by some phenolic compounds. Bioorg. Med. Chem. Lett. 21: 4259–4262.

61 

Sawaya M.R., Cannon G.C., Heinhorst S., Tanaka S., Williams E.B., Yeates T.O., Kerfeld C.A. (2006) The structure of beta-carbonic anhydrase from the carboxysomal shell reveals a distinct subclass with one active site for the price of two. J. Biol. Chem. 281: 7546–7555.

62 

Sayre R. (2010) Microalgae: the potential for carbon capture. Oxford J. Sci. Math. Biosc. 60: 722–727.

63 

Sheikhshoaie I., Ebrahimipour S.Y., Sheikhshoaie M., Rudbari H.A., Khaleghi M., Bruno G. (2014) Combined experimental and theoretical studies on the X-ray crystal structure, FT-IR, 1H NMR, 13C NMR, UV–Vis spectra, NLO behavior and antimicrobial activity of 2-hydroxyacetophenone benzoylhydrazone. Spectrochim. Acta A: Mol. Biomol. Spectr. 124: 548–555.

64 

Singh V.K., Singh S.K., Singh P.K., Verma H., Pandey K.D., Singh P.K., Kumar A. (2020) Impact of pesticides applications on the growth and function of cyanobacteria. [in:] Advances in Cyanobacterial Biology. Academic Press: 151–162.

65 

Sivonen K., Jones G. (2009) Cyanobacterial toxins. Encycl. Microbiol. 290: 307.

66 

Strop P., Smith K.S., Iverson T.M., FerryJ.G., Rees D.C. (2001) Crystal structure of the “cab"_type beta class carbonic anhydrase from the archaeon Methanobacterium thermoautotrophicum. J. Biol. Chem. 276: 10299–10305.

67 

Suganya S., Nandagopal B., Anbarasu A. (2017) Natural inhibitors of HMG CoA reductase – an in silico approach through molecular docking and simulation studies. J. Cell Biochem. 118(1): 52–57.

68 

Supuran C.T. (2008) Carbonic anhydrases – an overview. Curr. Pharm. Des. 14: 603–614.

69 

Supuran C.T. (2010) Carbonic anhydrase inhibitors. Bioorg. Med. Chem. Lett. 20(12): 3467–3474.

70 

Supuran C.T., Scozzafava A., Conway J. (2004) Carbonic anhydrase: its inhibitors and activators. Vol. 1. CRC Press.

71 

Thompson M.A. (2004) ArgusLab 4.0. 1. Planaria Software LLC, Seattle, WA, 98155. http://www.ArgusLab.com

72 

Tripathi P., Ghosh S., Talapatra S.N. (2019) Bioavailability prediction of phytochemicals present in Calotropis procera (Aiton) R. Br. by using Swiss-ADME tool. World Sci. News 131: 147–163.

73 

Tsuzuki M., Okada K., Isoda H., Hirano M., Odaka T., Saijo H., Fujiwara S. (2019) Physiological properties of photoauto trophic microalgae and cyanobacteria relevant to industrial biomass production. Marine Biotech. 21(3): 406–415.

74 

Velmurugan B.K., Rathinasamy B., Lohanathan B.P., Thiyagarajan V., Weng C.F. (2018) Neuroprotective role of phytochemicals. Molecules 23(10): 2485.

75 

Venkatesh R., Vidya R., Kalaivani K. (2014) Gas chromatography and mass spectrometry analysis of solanum villosum (mill) (solanaceae). Int. J. Pharm. Sci. Res. 5(12): 5283.

76 

Winter A.N., Brenner M.C., Punessen N. et al. (2017) Comparison of the neuroprotective and anti-inflammatory effects of the anthocyanin metabolites, protocatechuic acid and 4-hydroxybenzoic acid. Oxid. Med. Cell Longev. 2017: 6297080.

77 

Zuo Z., Chen Z., Zhang R., Gao Y., Wang Y. (2018) A potential algaecide from the pruning wastes of grape (Vitis vinifera) – stems and leaves. Environ. Eng. Manag. J. 17(12): 2895–2904.

Copyright: © 2023 Institute of Bioorganic Chemistry, Polish Academy of Sciences This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs (CC BY-NC-ND) (https://creativecommons.org/licenses/by-nc-nd/3.0/legalcode),.allowing third parties to download and share its works but not commercially purposes or to create derivative works.
 
Stosujemy się do standardu HONcode dla wiarygodnej informacji zdrowotnej This site complies with the HONcode standard for trustworthy health information: verify here