Skip to main content

ORIGINAL RESEARCH article

Front. Mol. Biosci., 28 May 2024
Sec. Biological Modeling and Simulation

Identification of the inhibition mechanism of carbonic anhydrase II by fructooligosaccharides

Yue MuYue Mu1Qingyang Meng
Qingyang Meng2*Xinyi FanXinyi Fan2Shuyun XiShuyun Xi2Zhongli XiongZhongli Xiong3Yihua WangYihua Wang3Yanling HuangYanling Huang3Zhen Liu
Zhen Liu1*
  • 1School of Chemical Engineering, East China University of Science and Technology, Shanghai, China
  • 2Shanghai Pechoin Biotechnology Co., Ltd., Shanghai, China
  • 3Shanghai Zhengxin Biotechnology Co., Ltd., Shanghai, China

Polygonatum sibiricum (P. sibiricum), recognized as a precious nourishing Chinese traditional medicine, exhibits the pharmacological effect of anti-aging. In this work, we proposed a novel mechanism underlying this effect related to the less studied bioactive compounds fructooligosaccharides in P. sibiricum (PFOS) to identify the inhibition effect of the small glycosyl molecules on the age-related zinc metalloprotease carbonic anhydrase II (CA II). Molecular docking and molecular dynamics simulation were used to investigate the structural and energetic properties of the complex systems consisting of the CA II enzyme and two possible structures of PFOS molecules (PFOS-A and PFOS-B). The binding affinity of PFOS-A (−7.27 ± 1.02 kcal/mol) and PFOS-B (−8.09 ± 1.75 kcal/mol) shows the spontaneity of the binding process and the stability of the combination in the solvent. Based on the residue energy decomposition and nonbonded interactions analysis, the C-, D- and G-sheet fragments of the CA II were found to be crucial in binding process. Van der Waals interactions form on the hydrophobic surface of CAII mainly with 131PHE and 135VAL, while hydrogen bonds form on the hydrophilic surface mainly with 67ASN and 92GLN. The binding of PFOS results in the blocking of the zinc ions pocket and then inhibiting its catalytic activity, the stability of which has been further demonstrated by free energy landscape. These findings provide evidence of the effective inhibition of PFOS to CA II enzyme, which leads to a novel direction for exploring the mechanism of traditional Chinese medicine focused on small molecule fructooligosaccharides.

GRAPHICAL ABSTRACT

1 Introduction

Polygonatum sibiricum (P. sibiricum), known as the “GOLD” in traditional Chinese medicine, is a medicinal and food homologous plant belonging to the Liliaceae family. (Wang et al., 2019). It has been employed as a constituent of traditional medicinal formulations with multifaceted pharmacological properties for over 2000 years (An et al., 2006; Han et al., 2020; Shen et al., 2022). In recent years, considerable research efforts have been focused on investigating the potential health benefits of P. sibiricum (Wang et al., 2019; Li et al., 2021; Liu et al., 2022). Various bioactive compounds such as polysaccharides, flavonoids, and steroidal saponins have been extracted.

The anti-aging function of P. sibiricum has received considerable attention among its various pharmacological applications. Multiple mechanisms of this function have been proposed, primarily through its antioxidative properties, protection of cellular components, and improvement in metabolic functions (Cui et al., 2018). However, related research primarily revolves around the polysaccharides. Scant attention was given to other biologically active components such as oligosaccharides, which also possess potential pharmacological properties contributing to anti-aging function of P. sibiricum. A comprehensive understanding of the intricate molecular mechanism underlying this effect remains limited (Zhao et al., 2018; Zheng, 2020; Wang et al., 2022).

Notable, the significant presence of fructooligosaccharides in medicinal plants and their crucial role in the treatment of various diseases has been proposed in recent research (Chen et al., 2015). The bioactivity of fructooligosaccharides resulting from the hydrolysis of these polysaccharides in P. sibiricum has also been highlighted (Liu et al., 2004; He et al., 2021). This particular component may have the potential to protect the cellular components such as endothelial cells and mitochondria. Given their prominent role and the preliminary evidence supporting their unique contributions to the pharmacological profile of P. sibiricum, it was crucial to isolate and analyze these compounds specifically to clarify their role in anti-aging. The two possible structures of fructooligosaccharides in P. sibiricum (PFOS) was identified in previous research (Figure 1).

Figure 1
www.frontiersin.org

Figure 1. Investigation of the inhibition effect of P. sibiricum Fructooligosaccharides on carbon anhydrase II based on molecular docking and molecular dynamics simulation.

Carbonic anhydrase II (CA II) emerges as a possible target of PFOS. It shows a distribution across diverse human organs and tissues. As the first zinc metalloprotease identified in 1933, (Meldrum and Roughton, 1933), CA II represents a pivotal isoform within the spectrum of human carbonic anhydrase (CAs, CA I-CA XIV). Its significance lies in its indispensable involvement in one of the most fundamental reactions in human existence: the conversion of CO2 and H2O into bicarbonate and protons (Imtaiyaz Hassan et al., 2013; Nocentini et al., 2019; Occhipinti and Boron, 2019). Given the pervasive nature of the reactions, CA II influence a wide range of metabolic processes, including the facilitation of transmembrane CO2 and O2 transport during the respiration, the mediation of transepithelial fluid and acid-base exchange and the maintenance of the whole-body acid-base equilibrium (Alterio et al., 2012; Supuran, 2021; Angeli et al., 2023).

Actually, the inhibition of the CA II has attracted considerable attention (Nocentini et al., 2019). Based on the identification of the structure of the CA II, physiological functions, and pathological roles, it has been selected as a drug target for the treatment of diverse diseases and impairments, including glaucoma, tumorigenesis, obesity calcification (Lie-Injo et al., 1970; Supuran, 2008). Extensive efforts have been devoted to the design of the CA II inhibitors to treat these diseases and impairments by inhibiting its abnormal expression or reducing its activity (Nocentini et al., 2004; Scozzafava et al., 2006; Carta and Supuran, 2013; Supuran et al., 2015; Occhipinti and Boron, 2019; Ghorai et al., 2020).

Currently, the most active classes of the CA II inhibitors are sulfonamides and their bioisosters, sulfamides, and sulfamates (Supuran et al., 2004; Husain and Madhesia, 2012). These inhibitors play an active role by binding with the catalytic zinc ion in the CA II or blocking the reaction pocket. However, the usage of them is still limited due to the safety concerns arising from the weak selectivity towards the CA II compared to other CA isozymes. Therefore, the search for novel compounds for the CA II inhibition has been motivated by this obstacle (Supuran, 2011; Poulsen and Davis, 2014; Carreyre et al., 2017; Supuran, 2020; Supuran, 2022a; Supuran, 2022b), as these compounds exhibit a favorable profile with lower toxicity and reduced side effects compared to synthetic small molecule inhibitors.

Intriguingly, recent scientific findings have pinpointed a correlation between age-related physiological impairments and the CA II quantity within mitochondria of specific cell types, (Pollard et al., 2016; Feng et al., 2018), which signifies an increase of the CA II activity concurrent with aging. Consequently, the novel compound PFOS may have the potential to inhibit the CA II activity, providing a plausible explanation for the observed anti-aging function of P. sibiricum.

This study aims to comprehensively elucidate the intricate mechanisms underlying the anti-aging properties of P. sibiricum, especially with respect to under-explored bioactive compounds fructooligosaccharides, and to highlight the necessity to identify novel non-toxic inhibitors for the CA II. We studied the inhibition mechanism of PFOS on the CA II for three main objectives: 1) assessment of the binding affinities; 2) comprehensive investigation of the protein-ligand interactions; 3) analysis of the intricate dynamics characteristics through molecular docking and molecular dynamics (MD) simulations.

2 Materials and methods

2.1 Molecular docking

Molecular docking is a powerful method for high-throughput virtual screening of small-molecule inhibitors and predicting the ligand-protein interactions in computer-assisted drug design. It consists of three main types: rigid docking, semi-flexible docking, and flexible docking, which are categorized in ascending order of protein flexibility. Semi-flexible docking was chosen for its balance between computational efficiency and accuracy, allowing modifications in the conformations of ligands and residues within a defined area.

The crystal structure of the CA II (PDB: 12CA) was sourced from the RCSB library and does not include post-translational modifications (Berman et al., 2000). Both PFOS-A and PFOS-B were fully optimized utilizing the Gaussian09 software (Frisch MJT et al., 2009) at the density functional theory level, employing the B3LYP function with the London-dispersion correction by Grimme and the def2-SVP basis set (Schäfer et al., 1994; Grimme et al., 2011).

Autodock Vina was utilized for the docking of PFOS with CA II (Trott and Olson, 2010). A 10 × 10 × 10 nm box, set up using default parameters, encompassed almost the entire active surface of CA II, centered on the typical binding pocket identified in previous research. This pocket, delineated by hydrophobic residues (Ile91, Val121, Phe131, Leu141, Val143, and Leu198) and hydrophilic residues (Asn62, His64, Asn67, and Gln92), is crucial for rapid catalytic conversion of CO2 to bicarbonate. The energy range was set to 5.0 kcal/mol, and the exhaustiveness parameter was configured as 22 (Agarwal and Smith, 2023). The optimal pose with the best scores within multiple iterations of the docking process was selected as the initial structure for the subsequent MD simulation, according to the affinity assessment through scoring function and binding mode analysis (Tanchuk et al., 2016).

The frontier molecular orbital analysis was conducted to assess the docking reactivity of these molecules due to the significance of HOMO and LUMO (EHOMO and ELUMO) energy gaps in chemical reactivity. Various physical and chemical parameters, including chemical hardness (η), softness (s), and electrophilicity index (w), were calculated and analyzed.

2.2 Molecular dynamics simulation

Based on the initial structures obtained from the molecular docking results, MD simulations were conducted using the GROMACS 2016 software package, (Abraham et al., 2015), and the proteins and small molecules were parameterized by the all-atom CHARMM36 and CHARMM general force field (CGenFF), respectively.

The complex was put in a 10 × 10 × 10 nm periodic cubic box and then filled with the TIP3P water molecules. (Mark and Nilsson, 2001). Then, the chloride ions were added to neutralize the system charge. During the energy minimization process, the unreasonable interactions and atomic collisions were removed using the steepest decent and conjugated gradient methods. Next, the temperature of each system was incrementally increased from 0 K to 300 K during the 1 ns simulation using an NVT ensemble. The temperature of the protein-ligand and solvent-ion groups was controlled separately using a Berendsen thermostat to obtain more reliable results. After that, a Parrinello–Rahman barostat was used to maintain the pressure of the systems at 0.1 MPa in NPT equilibration for 1 ns. Finally, MD simulations were conducted for 100 ns, and the last 10 ns trajectory was taken every 10 ps for subsequent analysis.

During the simulation, the length of the hydrogen bonds was constrained by LINear Constraint Solver algorithm. The cut-off for short-range nonbonded interactions was set at 1.0 nm, while for long-range electrostatic interactions, the Particle Mesh Ewald algorithm was employed without an explicit cut-off. The determination of cut-offs was guided by empirical practices widely adopted in previous studies (Essmann et al., 1995; Kabedev et al., 2023; Sun et al., 2023). Newton’s classical equations of motion were integrated using the Verlet leapfrog algorithm. (Hockney et al., 1974). The stability of the systems was checked by the root mean square deviation (RMSD) of the protein backbone.

2.3 Binding free energy

To further assess the interaction between PFOS and CAII, the binding energy was calculated and partitioned into different parts according to the type of interaction by the MM-PBSA method implemented in the gmx_MMPBSA software package (Miller et al., 2012; Suri and Naik, 2015). It should be noted that the binding free energies calculated by the MM-PBSA method is system-dependent and sensitive to the chosen parameters, which is somewhat less accurate compared to the free energy perturbation or thermodynamic integration (Breznik et al., 2023). However, the computationally efficient MM-PBSA method has proven to be reliable in the field of computational drug design and evaluation (Genheden and Ryde, 2015; Wen et al., 2019).

This calculation was conducted using a set of snapshots extracted from the final 10 ns trajectory of the MD simulations, sampled at 10 ps intervals. The binding free energy calculations in the MM-PBSA method can be concisely represented by the following equations:

ΔG=GcomplexGreceptor+Gligand=ΔHTΔS(1)
ΔG=ΔEgas+ΔGsolTΔS(2)

where,

ΔEgas=ΔEbinded+ΔEnonbinded(3)
ΔEbinded=ΔEbind+ΔEangle+ΔEdihedral(4)
ΔEnonbinded=ΔEele+ΔEvdw(5)

and,

ΔGsol=ΔGpolar+ΔGnonpolar=ΔGPB+ΔGSA(6)

where,

ΔGSA=NPTENSION×ΔSASA+NPOFFSET=γΔSASA(7)

Within these equations, Gcomplex, Greceptor, and Gligand are the free energy of the complex system, CA II and small molecule inhibitors, respectively, collectively contributing to the binding free energy (∆G). This free energy can be further partitioned into the sum of the binding enthalpy (∆H) and the conformational entropy change upon inhibitor binding (-T∆S), as depicted in Eq. 1.

The enthalpic term, denoted as ∆H, can be formulated as the sum of ∆Egas and ∆Gsol, as elaborated in Eq. 2. ∆Egas accounts for the energy changes computed through molecular mechanics in the gas phase, while ∆Gsol pertains to solvent-related energy contributions.

In Eqs 35, the term ∆Egas is divided into two separate contributions: ∆Ebinded and ∆Enonbinded. ∆Ebinded encompasses energy contributions arising from binding, angle, and dihedral alterations, while ∆Enonbinded can be further partitioned into the ∆Eele and the ∆Evdw.

The computation of the solvation-free energy involves two components: ∆Gpolar, assessed via the Poisson–Boltzmann approach, and the ∆Gnonpolar, estimated through a SASA-only model (Eq. 6), assuming proportionality to the solvent accessible surface (∆SASA) outlined in Eq. 7.

Notably, the entropy change (∆S) effect is ignored due to computational constraints and the substantial uncertainty it introduces into the results, stemming from considerable errors during the calculation process. Several studies have shown that the MM-PBSA method can give high calculation accuracy without accounting for ∆S (Qin et al., 2010; Hou et al., 2011; Guo et al., 2021).

2.4 Simulation analysis

The trajectory obtained from the MD simulation was analyzed using auxiliary tools provided within the GROMACS 2016 package. The interaction between small molecules and CA II receptors is analyzed by Protein-Ligand Interaction Profiler (PLIP). Furthermore, the occupancies of the intermolecular hydrogen bonds within the PFOS systems were calculated using the Visual Molecular Dynamics (VMD) 1.9.3 software. With the focus on direct interaction between PFOS and CAII, the trajectory used in hydrogen binding analysis have been pre-processed to remove the waters. Hydrogen bonds were identified based on established criteria: an acceptor-hydrogen-donor angle exceeding 135° and an acceptor-hydrogen atom distance less than 3.5 Å.

In addition, a principal component analysis was conducted to delineate the dominant motion of the CA II and the correlated motions among the atoms extracted from the MD trajectories within the PFOSs system (Ma et al., 2017; John et al., 2020). Notably, an intermolecular motion was excluded by eliminating protein translation and rotation via the auxiliary tools available in the GROMACS suite. The eigenvectors, which signify the directions of the atomic motions, along with their corresponding eigenvalues that indicated their magnitudes, were calculated by evaluating the covariance matrix of the α-carbon atoms of the CA II protein. This process aimed to elucidate the impact of ligands on the conformational distribution of the protein (Mu et al., 2005; Papaleo et al., 2009). Subsequently, the first two principal components (PC1 and PC2) were utilized as reaction coordinates to construct the free energy landscape, and the free energy values were calculated with Eq. 8.

G=kBT×lnP(8)

Where kB represents Boltzmann constant, T stands for the simulation temperature, and P is the probability of the conformational distribution.

3 Results and discussion

3.1 Docking of CA II and PFOS

The initial binding structures of CA II with each PFOS compound were obtained through molecular docking. Firstly, to validate the reliability of the docking method, redocking experiments were performed on the complex of CA II with N-[(4-methylphenyl)methyl]-thiophene-2,5-disulfonamid (PDB: 1BNA, inhibitor: AL5) (Boriack-Sjodin et al., 1998). The results demonstrate that the conformation of the complex obtained by docking method matches with the native crystal structure located in the active pocket of the CA II, as shown in Supplementary Figure S1, with an RMSD of 0.78 Å. Furthermore, the interactions of the docked pose, assessed using the PLIP program, (Adasme et al., 2021), are consistent with the research conducted by Sjodin, (Boriack-Sjodin et al., 1998), further validating the accuracy of the docking method. The correlative analysis is shown in Supplementary Figure S2. Having validated the docking method, the two PFOS molecules were then docked to the CA II using the same docking method (Figure 2). The conformations with the best score for both structures are in the pocket adjacent to the zinc ions.

Figure 2
www.frontiersin.org

Figure 2. Initial structures of the PFOS-A and PFOS-B bound to CA II. The skeleton of the PFOS-A is in blue, PFOS-B is in green, and the bound residues are in purple.

The docking results reveal that both PFOS structures effectively bind to the active site on CA II as classical inhibitor AL5 used for docking verification. The chemical softness of the small molecules underscores their binding reactivity, with values of 7.58 for PFOS-A and 7.96 for PFOS-B, since a higher value of chemical softness typically correlates with increased molecular reactivity. The additional details on other chemical reactivity parameters are provided in Supplementary Table S1.

Analysis based on PLIP reveals that PFOS-A predominantly binds to residues ASN-62, HIS-64, ASN-67, ASP-72, PHE-131, GLY-132, THR-200, and PRO-201, while PFOS-B exhibits binding interactions with TRP-5, ASN-62, HIS-64, ASN-67, GLN-92, HIS-94 and THR-200. These residues are primarily located on the hydrophilic surface of the protein, attributed to the presence of multiple hydroxyl groups in these glycosyl molecules. However, the neglect of the structural flexibility inherent in the majority of protein structures in the docking method is a formidable obstacle to accurate binding analysis. Nevertheless, the docking conformations can still serve as valuable initial input structures for subsequent MD simulations and as a reference for the binding analysis. The off-target effects of PFOS have also been preliminarily assessed by the molecular docking method (Supplementary Figure S3). It should be noted that the activity of the non-target protein CA I may be slightly affected.

3.2 RMSD

With the initial structures obtained by molecular docking, 100 ns MD simulations of the systems were performed to get a more detailed and accurate understanding of their behavior. The RMSD was calculated for the resides within 10 Å of the ligand to assess the structural stability throughout the 100 ns simulations. As illustrated in Figure 3A, both PFOS complexes reached equilibrium after 50 ns, exhibiting the RMSD values of less than 2.8 Å for the majority of the equilibration stage, which are more stable than the CA II without inhibitor (averaging 3.0 Å in the last 50 ns). This observation suggests that the presence of the PFOS small molecules stabilized the core backbone of the CA II. Notably, during the equilibration period, the PFOS-B system demonstrated slightly higher stability with lower RMSD values than the PFOS-A system, while the opposite trend was observed during the early stage. This discrepancy may be attributed to the disruption or instability of the hydrogen bonds present in the initial docking structure of the PFOS-A system. Nevertheless, since all three systems reached equilibrium, further analysis can be conducted based on the obtained MD simulation trajectory.

Figure 3
www.frontiersin.org

Figure 3. (A) The RMSDs of residues reside within 10 Å of the ligand in the CA II, PFOS-A, and PFOS-B systems. (B) The RMSF plots of all the systems during 100 ns MD simulation and the β-sheet region are named as reported by Eriksson (Eriksson et al., 1988); (C) The secondary structure of the CA II. The β-sheet region is shown in blue, and the helix region is shown in red. (D) RMSF of residues located in the active site of CA II.

3.3 RMSF

The Root mean square fluctuation (RMSF) was further calculated to measure the per-residue flexibility of the CA II, PFOS-A, and PFOS-B complex systems. As illustrated in Figure 3B, the closely overlapping profiles of PFOS complexes and CA II suggest that the inhibitor has little effect on the regional flexibility of the protein. The stable regions mainly occurred at the β-sheets with the RMSF less than 0.8 Å, while the helix domain showed higher RMSF values up to 2.5 Å.

The secondary structure of CA II is shown in Figure 3C, which consists of 10 stranded twisted β-sheet and 7 helices and the zinc ions presented in the center bottom of the protein chain surrounded by the β-sheet regions. According to the secondary structure, all helix regions with high RMSF values were far away from the binding pocket and had little effect on the binding process. On the contrary, the β-sheet, especially the C-sheet, D-sheet, and G-sheet segments, showed relatively lower RMSF values in complex with PFOS as inhibitors compared to the CA II system. These segments are identified as the key sites for the combination of PFOS. In addition, systems with PFOS-B as the inhibitor are found to have higher residue stability than PFOS-A, which is in accordance with the previous RMSD analysis.

Further RMSF analysis focused on residues within the active site of CA II (Figure 3D). Upon binding with PFOS, the flexibility of these residues generally decreased, particularly for residues 131PHE and 135VAL. These two residues have also been identified as key contributors to the protein-ligand interaction in subsequent energy analysis to accommodate the hydroxyl group of PFOS.

3.4 Binding free energy

It has been shown that the PFOS effectively stabilize the CA II protein via binding with some key residues situated within the β-sheet region. However, the lack of a zinc-binding group distinguishes them from most existing CA II-targeting inhibitors. It is difficult for PFOS structures to bind directly to the catalytic zinc ion, which plays a critical role in the catalytic function of the enzyme. Hence, it is crucial to validate the exact binding affinities of these two small molecules with CA II, wherein the MM-PBSA method was employed. This calculation relied on the snapshots extracted from the final 10 ns of the equilibrium stage in the MD simulation of each respective system. The same methodology was also applied to the experimentally tested inhibitor AL5, employed in the docking validation, serving as the benchmark (Table 1).

Table 1
www.frontiersin.org

Table 1. Binding energies in PFOS-A, PFOS-B, and CA II systems by MM-PBSA method (kcal/mol).

The two PFOS structures show comparable binding affinity to the CA II protein (Table 1). The PFOS-B (−8.09 kcal/mol) shows a slightly lower ΔG compared to PFOS-A (−7.27 kcal/mol), a trend consistent with the RMSD and RMSF (Figure 3) results. From a thermodynamic perspective, these negative values corroborate the spontaneity of the binding process and the stability of the combination in the solvent. This level of binding stability aligns with findings in similar research, where a small glycosyl molecule interacts with the human sweet taste receptor, exhibiting a ΔG of −8.53 ± 1.78 kcal/mol, as validated through MD simulations and experimental data. (Mahalapbutr et al., 2020). However, the biological activity of PFOS-A and PFOS-B is relatively lower when compared to typical CA II inhibitors exemplified by the benchmark compound AL5. This difference in activity can be attributed to the absence of the zinc-binding group in the PFOS structure, which can form much more favorable electronic interactions to maintain the combination of the complex.

To further investigate the potential inhibition mechanism of PFOS on CA II, common features of the two structures are first identified by energy decomposition analysis. The total binding energy can be decomposed into different components including electrostatic energy, van der Waals energy, polar solvation-free energy, and non-polar solvation energy. Thus, relative importance of each interaction within the binding process can be assessed.

The dominant propellant for binding for both the PFOS-A and PFOS-B systems is the total nonpolar binding energy, which can be further separated into van der Waals energy and non-polar solvation-free energy (Table 1). Both these components contribute positively to the binding process. In addition, negative electrostatic energy was observed in all systems, with significantly lower values in both PFOS systems compared to AL5. This discrepancy can be attributed to the abundant hydroxyl groups in PFOS-A and PFOS-B (Figure 1), which exhibit substantial potential for establishing hydrogen bonding interactions. However, a significant unfavorable influence was detected concerning the total polar free energy, primarily attributed to the desolvation process, chiefly driven by the considerable polarity exhibited by the PFOS molecules. The formation of intramolecular hydrogen bonds by PFOS molecules may also impose limitations on their interactions with CA II molecules. The slightly weaker binding free energy of PFOS-A observed during dynamic simulations will be explained by following energy decomposition.

3.5 Per-residue energy decomposition

In order to elucidate the detailed interaction between PFOS molecules and the CA II protein, the identification of key residues critical for the binding process was pursued through energy decomposition analysis. The analysis encompassed residues within a 10 Å radius of the PFOS molecules. The ten residues exhibiting the most substantial contributions to the binding energy were ranked and presented in Figure 4. Key residues identified by contributions exceeding the threshold of −0.5 kcal/mol were summarized in Table 2 and Table 3. The detailed energy decomposition result is shown in Supplementary Tables S2, S3.

Figure 4
www.frontiersin.org

Figure 4. Total binding free energy of PFOS-A and PFOS-B in complex with CA II protein. Representative 3D structures illustrating the binding region of the small molecules on the CA II protein were taken from the last 10 ns MD snapshot with the lowest energy and color-coded according to their binding free energy values.

Table 2
www.frontiersin.org

Table 2. Key residues energy decomposition in PFOS-A system (kcal/mol).

Table 3
www.frontiersin.org

Table 3. Key residues energy decomposition in PFOS-B system (kcal/mol).

PFOS-A structure interacted with 5 key residues: TRP-5, PRO-202, PHE-20, PHE-131, and VAL-135 (Figure 4). Meanwhile, the PFOS-B structure interacted with 9 pivotal residues: PHE-131, VAL-135, LEU-198, PRO-202, LEU-204, TYR-7, THR-200, LEU-141, and TRP-5, which is nearly twice the key residues that participated in the combination of PFOS-A system. Based on the analysis of the docking conformation of the two systems, the reason for the fewer key residues in PFOS-A may be attributed to its stick line backbone structure in the combination process, which makes it less possible to access the residues located at the entrance of the docking pocket to form van der Waals interaction with them. The fusiform structure of PFOS-B makes it more suitable to block the docking pocket with more accessible residues. In addition, the van der Waals energy contribution of the key residues is significantly larger compared to the electrostatic energy, which is consistent with the previous analysis in Table 1. The TRP-5 in the PFOS-A system and the PHE-131 in the PFOS-B system exhibited the most favorable ΔGtotal values. For both PFOS structures, residues PHE-131, VAL-13 were found to play a significant role in the van der Waals interactions at the hydrophobic side of the binding pocket entrance of the CA II protein. These findings are consistent with the prior studies suggesting the importance of these residues in the binding of the hydrophobic tail in Zinc binding group inhibitors. (Alterio et al., 2012).

3.6 Binding mode analysis

The electrostatic interaction also plays an important role in inducing the complexation of PFOS molecules and CA II, in addition to van der Waals and the non-polar solvation effect. Further characterization of this kind of interaction anchoring the PFOS small molecules into the binding pocket was conducted. The occupancy of the direct intermolecular hydrogen bonds was statistically analyzed over the last ten MD simulations using the PLIP and Visual Molecular Dynamics (VMD) program. (Humphrey et al., 1996). The percentage occupancy of hydrogen bonds in the PFOS-A and PFOS-B systems are presented in Supplementary Tables S4, S5. The stable hydrogen bonding has been summarized according to the standard for stable hydrogen bonding, that the hydrogen bonding donor-acceptor interaction occupancy is greater than 10% (Table 4; Table 5). As expected, both PFOS systems showed the formation of hydrogen bonds with polar residues within the C-sheet, D-sheet and G-sheet of the CA II protein. However, it can be seen that the number and stability of these hydrogen bond formations are different for the two systems.

Table 4
www.frontiersin.org

Table 4. Main hydrogen bond occupancies in PFOS-A system.

Table 5
www.frontiersin.org

Table 5. Main hydrogen bond occupancies in PFOS-B system.

In PFOS-A system, three main hydrogen bonds formed, including H-N2 of ASN-67 with O4 and H-N2 of TRP-5 with adjacent O72 and O70, respectively (Figure 5). This observation is corroborated by the aforementioned energy decomposition analysis, which revealed the negative electrostatic energy contributions from these residues, particularly notable for 5TRP with a value of −1.94 kcal/mol. Conversely, in the PFOS-B and CA II complex, the formation of hydrogen bonds was chiefly facilitated by PRO-201, GLN-92, and GLN-136. Notably, more hydrogen bonds are present in system PFOS-B compared to PFOS-A, suggesting that PFOS-B has a stronger electrostatic interaction, in agreement with free energy calculations. However, the most stable H-bond was detected in the PFOS-A system, which persisted for 39.6% of the simulation time. While, the majority of the other hydrogen bonds endured for less than 20% of the simulation duration. In addition, the overall analysis of the H-bonds occupation data in the table indicates general consistency with the results of molecular docking, which also shows the H-bonds formation of HIS64, ASN67, THR200, and PRO201 in PFOS-A system and ASN62, ASN67, GLN92 and GLN 136 in PFOS-B system.

Figure 5
www.frontiersin.org

Figure 5. Percentage of intermolecular H-bond occupation of the CA II residues contributing to PFOS-A and PFOS-B. The 3D representative structures illustrating the H-binding between each small molecule and the CA II protein were extracted from the last snapshot of the MD simulation. The black dashed line indicates H-bond formation in the simulation process.

Considering the presence of multiple hydroxyl groups in PFOS, the formation of hydrogen bond largely depends on the structure and conformation of PFOS, which may be dynamic. Whereas, the ASN67, GLN92, GLN135 and PRO201 showed as the most potential residue to form hydrogen bond with PFOS.

3.7 Inhibition mechanism

Based on the above analysis, it can be concluded that the van der Waals interaction domain governs the binding of PFOS molecules, while the electrostatic energy contribution resulting from the formation of hydrogen bonds also has a crucial influence on the binding of PFOS molecules to the hydrophilic residue.

By integrating the comprehensive analyses of both van der Waals and electrostatic interactions, the elucidation of the binding site and interaction mechanism between PFOS small molecules and the CA II protein has been elucidated in Figure 6A. Both PFOS structures effectively bind to the entrance of the pocket that accommodates the Zn ion within CA II. While, the double-branched fusiform PFOS-B structure is better at blocking the pocket. The small molecules bind to the hydrophobic surface of the protein through van der Waals interactions and form hydrogen bonds with the hydrophilic surface through electrostatic interactions. This dual interaction strategy effectively contributes to the binding of PFOS, which blocks the carbon dioxide conversion pathway and then results in the inhibition of CA II activity.

Figure 6
www.frontiersin.org

Figure 6. (A) Plot of the CA II and PFOS/CA II complexes [right: PFOS-A system (top panel) and PFOS-B system (bottom panel)] demonstrating the interaction mechanism in the binding pocket. (B) The free energy landscape of CA II, PFOS-A, and PFOS-B systems. PC1 and PC2 represent principal components 1 and 2, respectively.

Based on the characterization of the binding pocket, the free energy landscapes (FEL) of the three systems, CA II, PFOS-A, and PFOS-B, were analyzed to further investigate the impact of the small molecules on the overall conformational dynamics of the CA II protein and also furnish a qualitative description of their inhibitory potential. The FEL was constructed by projecting the 1000 snapshots from the MD simulation onto the first two principal components, which were employed as the coordinates, as illustrated in Figure 6B.

In the FEL of the PFOS-B system, a distinct blue region indicates the high probability of the existence of a stable conformation. In contrast, the two prominent conformational regions were identified for CA II in the absence of an inhibitor molecule, indicating the stabilizing effect of PFOS-B binding. In the case of the PFOS-A system, its conformational space is characterized by a large region containing two distinct energy traps. This means that the PFOS-A system has a comparatively lower stability. In addition, mirroring the results of the RMSF analysis, the FEL of the PFOS-A system shows numerous smaller basins that are absent from the FEL graph of the PFOS-B system, implying a higher flexibility of the complex binding with PFOS-A.

4 Conclusion

In this work, the inhibition effect of PFOS on CA II protein activity was investigated through semi-flexible molecular docking, MD simulations, and MM-PBSA methods. Although the biological activity of PFOS is relatively weaker compared to typical CA II inhibitors, both structures of PFOS bind to the active pocket of CA II with high affinity (PFOS-A: −7.27 ± 1.02 kcal/mol, PFOS-B: −8.09 ± 1.75 kcal/mol). The RMSD and RMSF analysis affirmed the stabilizing effect of PFOS on the CA II structure. These results confirm the spontaneity of the binding process and the stability of the combination in the solvent, indicating the inhibitory potential of PFOS molecules.

The analysis of the binding modes revealed the importance of the three key components of CAII: the C-, D-, and G-sheets, on which a dual interaction strategy was identified by energy decomposition and hydrogen bond occupation analysis. In the binding process between PFOS and CAII, van der Waals interactions form on the hydrophobic surface of CAII mainly with 131PHE and 135VAL, while hydrogen bonds form on the hydrophilic surface mainly with 67ASN, 92GLN, 136GLN and 201PRO. This dual interaction strategy can effectively contribute to the binding of PFOS, resulting in the block of the channel where the CO2 conversion reaction takes place and thus inhibiting the activity of CA II. Finally, the FELs of the three systems were analyzed to further determine the effect of the PFOS on the overall conformational dynamics of CA II proteins, demonstrating that both PFOS structures were effective in stabilizing CA II protein into one conformation with low free energy. The inhibitory effect of fructooligosaccharides on the activity of the aging associated protein CA II may provide a mechanistic basis for the anti-aging activity of P. sibiricum.

This study provides novel mechanistic insights into the anti-aging effects of P. sibiricum, paving the way for further research and development of P. sibiricum products. Further elucidation of the pharmacological mechanisms of other small-molecule bioactive compounds in traditional Chinese medicine is also of great importance, which would facilitate the application of traditional Chinese medicine.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Author contributions

YM: Investigation, Writing–original draft. QM: Conceptualization, Writing–review and editing. XF: Visualization, Writing–review and editing. SX: Formal Analysis, Writing–review and editing. ZX: Methodology, Software, Writing–review and editing. YW: Methodology, Software, Writing–review and editing. YH: Resources, Writing–review and editing. ZL: Project administration, Writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The study was funded through the Shanghai Pujiang Program (18PJ1402500) and the Open Project of the State Key Laboratory of Chemical Engineering (SKL-ChE-18C01) at the East China University of Science and Technology.

Acknowledgments

We acknowledge the supercomputer at the East China University of Science and Technology for its generous computing resources.

Conflict of interest

Authors QM, XF, and SX were employed by the Shanghai Pechoin Biotechnology Co., Ltd. Authors ZX, YW, and YH were employed by the Shanghai Zhengxin Biotechnology Co., Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2024.1398603/full#supplementary-material

References

Abraham, M. J., Murtola, T., Schulz, R., Páll, S., Smith, J. C., Hess, B., et al. (2015). GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 1-2, 19–25. doi:10.1016/j.softx.2015.06.001

CrossRef Full Text | Google Scholar

Adasme, M. F., Linnemann, K. L., Bolz, S. N., Kaiser, F., Salentin, S., Haupt, V. J., et al. (2021). PLIP 2021: expanding the scope of the protein-ligand interaction profiler to DNA and RNA. Nucleic Acids Res. 49 (W1), W530–W534. doi:10.1093/nar/gkab294

PubMed Abstract | CrossRef Full Text | Google Scholar

Agarwal, R., and Smith, J. C. (2023). Speed vs accuracy: effect on ligand pose accuracy of varying box size and exhaustiveness in AutoDock Vina. Mol. Inf. 42 (2), e2200188. doi:10.1002/minf.202200188

CrossRef Full Text | Google Scholar

Alterio, V., Di Fiore, A., D'Ambrosio, K., Supuran, C. T., and De Simone, G. (2012). Multiple binding modes of inhibitors to carbonic anhydrases: how to design specific drugs targeting 15 different isoforms? Chem. Rev. 112 (8), 4421–4468. doi:10.1021/cr200176r

PubMed Abstract | CrossRef Full Text | Google Scholar

An, J., Liu, J.-Z., Wu, C.-F., Li, J., Dai, L., Van, D. E., et al. (2006). Anti-HIV I/II activity and molecular cloning of a novel mannose/sialic acid-binding lectin from rhizome of Polygonatum cyrtonema Hua. Acta Biochim. Biophys. Sin. (Shanghai). 38 (2), 70–78. doi:10.1111/j.1745-7270.2006.00140.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Angeli, A., Kartsev, V., Petrou, A., Lichitsky, B., Komogortsev, A., Geronikaki, A., et al. (2023). Substituted furan sulfonamides as carbonic anhydrase inhibitors: synthesis, biological and in silico studies. Bioorg Chem. 138, 106621. doi:10.1016/j.bioorg.2023.106621

PubMed Abstract | CrossRef Full Text | Google Scholar

Berman, H. M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T. N., Weissig, H., et al. (2000). The protein data bank. Nucleic Acids Res. 28 (1), 235–242. doi:10.1093/nar/28.1.235

PubMed Abstract | CrossRef Full Text | Google Scholar

Boriack-Sjodin, P. A., Zeitlin, S., Chen, H. H., Crenshaw, L., Gross, S., Dantanarayana, A., et al. (1998). Structural analysis of inhibitor binding to human carbonic anhydrase II. Protein Sci. 7 (12), 2483–2489. doi:10.1002/pro.5560071201

PubMed Abstract | CrossRef Full Text | Google Scholar

Breznik, M., Ge, Y., Bluck, J. P., Briem, H., Hahn, D. F., Christ, C. D., et al. (2023). Prioritizing small sets of molecules for synthesis through in-silico tools: a comparison of common ranking methods. ChemMedChem. 18 (1), e202200425. doi:10.1002/cmdc.202200425

PubMed Abstract | CrossRef Full Text | Google Scholar

Carreyre, H., Carre, G., Ouedraogo, M., Vandebrouck, C., Bescond, J., Supuran, C. T., et al. (2017). Bioactive natural product and superacid chemistry for lead compound identification: a case study of selective hCA III and L-Type Ca2+ current inhibitors for hypotensive agent discovery. Molecules 22 (6), 915–911. doi:10.3390/molecules22060915

PubMed Abstract | CrossRef Full Text | Google Scholar

Carta, F., and Supuran, C. T. (2013). Diuretics with carbonic anhydrase inhibitory action: a patent and literature review (2005 - 2013). Expert Opin. Ther. Pat. 23 (6), 681–691. doi:10.1517/13543776.2013.780598

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X. Y., Wang, R. F., and Liu, B. (2015). An update on oligosaccharides and their esters from traditional Chinese medicines: chemical structures and biological activities. Evid. Based Complement. Altern. Med. 2015, 512675. doi:10.1155/2015/512675

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, X., Wang, S., Cao, H., Guo, H., Li, Y., Xu, F., et al. (2018). A review: the bioactivities and pharmacological applications of polygonatum sibiricum polysaccharides. Molecules 23 (5), 1170. doi:10.3390/molecules23051170

PubMed Abstract | CrossRef Full Text | Google Scholar

Eriksson, A. E., Jones, T. A., and Liljas, A. (1988). Refined structure of human carbonic anhydrase II at 2.0 A resolution. Proteins 4 (4), 274–282. doi:10.1002/prot.340040406

PubMed Abstract | CrossRef Full Text | Google Scholar

Essmann, U., Perera, L., Berkowitz, M. L., Darden, T., Lee, H., and Pedersen, L. G. (1995). A smooth particle mesh Ewald method. J. Chem. Phys. 103 (19), 8577–8593. doi:10.1063/1.470117

CrossRef Full Text | Google Scholar

Feng, S., Wang, S., Wang, Y., Yang, Q., Wang, D., and Li, H. (2018). Identification and expression of carbonic anhydrase 2, myosin regulatory light chain 2 and selenium-binding protein 1 in zebrafish Danio rerio: implication for age-related biomarkers. Gene Expr. Patterns 29, 47–58. doi:10.1016/j.gep.2018.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Frisch Mjt, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., Scalmani, G., et al. (2009) Gaussian 09, revision D.01. Wallingford, CT: Gaussian Inc.

Google Scholar

Genheden, S., and Ryde, U. (2015). The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin. Drug Discov. 10 (5), 449–461. doi:10.1517/17460441.2015.1032936

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghorai, S., Pulya, S., Ghosh, K., Panda, P., Ghosh, B., and Gayen, S. (2020). Structure-activity relationship of human carbonic anhydrase-II inhibitors: detailed insight for future development as anti-glaucoma agents. Bioorg Chem. 95, 103557. doi:10.1016/j.bioorg.2019.103557

PubMed Abstract | CrossRef Full Text | Google Scholar

Grimme, S., Ehrlich, S., and Goerigk, L. (2011). Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 32 (7), 1456–1465. doi:10.1002/jcc.21759

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Y., Liang, J., Liu, B., and Jin, Y. (2021). Molecular mechanism of food-derived polyphenols on PD-L1 dimerization: a molecular dynamics simulation study. Int. J. Mol. Sci. 22 (20), 10924. doi:10.3390/ijms222010924

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, C., Zhu, Y., Yang, Z., Fu, S., Zhang, W., and Liu, C. (2020). Protective effect of Polygonatum sibiricum against cadmium-induced testicular injury in mice through inhibiting oxidative stress and mitochondria-mediated apoptosis. J. Ethnopharmacol. 261, 113060. doi:10.1016/j.jep.2020.113060

PubMed Abstract | CrossRef Full Text | Google Scholar

He, L., Yan, B., Yao, C., Chen, X., Li, L., Wu, Y., et al. (2021). Oligosaccharides from polygonatum cyrtonema hua: structural characterization and treatment of LPS-induced peritonitis in mice. Carbohydr. Polym. 255, 117392. doi:10.1016/j.carbpol.2020.117392

PubMed Abstract | CrossRef Full Text | Google Scholar

Hockney, R. W., Goel, S. P., and Eastwood, J. W. (1974). Quiet high-resolution computer models of a plasma. J. Comput. Phys. 14 (2), 148–158. doi:10.1016/0021-9991(74)90010-2

CrossRef Full Text | Google Scholar

Hou, T., Wang, J., Li, Y., and Wang, W. (2011). Assessing the performance of the MM/PBSA and MM/GBSA methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations. J. Chem. Inf. Model. 51 (1), 69–82. doi:10.1021/ci100275a

PubMed Abstract | CrossRef Full Text | Google Scholar

Humphrey, W., Dalke, A., and Schulten, K. (1996). VMD: visual molecular dynamics. J. Mol. Graph 14 (1), 33–38. doi:10.1016/0263-7855(96)00018-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Husain, A., and Madhesia, D. (2012). Heterocyclic compounds as carbonic anhydrase inhibitor. J. Enzyme Inhib. Med. Chem. 27 (6), 773–783. doi:10.3109/14756366.2011.617882

PubMed Abstract | CrossRef Full Text | Google Scholar

Imtaiyaz Hassan, M., Shajee, B., Waheed, A., Ahmad, F., and Sly, W. S. (2013). Structure, function and applications of carbonic anhydrase isozymes. Bioorg Med. Chem. 21 (6), 1570–1582. doi:10.1016/j.bmc.2012.04.044

PubMed Abstract | CrossRef Full Text | Google Scholar

John, A., Vetrivel, U., Sivashanmugam, M., and Natarajan, S. K. (2020). Microsecond simulation of the proteoglycan-like region of carbonic anhydrase IX and design of chemical inhibitors targeting pH homeostasis in cancer cells. ACS Omega 5 (8), 4270–4281. doi:10.1021/acsomega.9b04203

PubMed Abstract | CrossRef Full Text | Google Scholar

Kabedev, A., Bergström, C. A. S., and Larsson, P. (2023). Molecular dynamics study on micelle-small molecule interactions: developing a strategy for an extensive comparison. J. Computer-Aided Mol. Des. 38 (1), 5. doi:10.1007/s10822-023-00541-1

CrossRef Full Text | Google Scholar

Li, X. L., Ma, R. H., Zhang, F., Ni, Z. J., Thakur, K., Wang, S., et al. (2021). Evolutionary research trend of Polygonatum species: a comprehensive account of their transformation from traditional medicines to functional foods. Crit. Rev. Food Sci. Nutr. 63 (19), 3803–3820. doi:10.1080/10408398.2021.1993783

PubMed Abstract | CrossRef Full Text | Google Scholar

Lie-Injo, L. E., Lopez, C. G., and Hart, P. L. (1970). Erythrocyte carbonic anhydrase activity in health and disease. Clin. Chim. Acta 29 (3), 541–550. doi:10.1016/0009-8981(70)90026-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, F., Liu, Y., Meng, Y., Yang, M., and He, K. (2004). Structure of polysaccharide from Polygonatum cyrtonema Hua and the antiherpetic activity of its hydrolyzed fragments. Antivir. Res. 63 (3), 183–189. doi:10.1016/j.antiviral.2004.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, S., Jia, Q. J., Peng, Y. Q., Feng, T. H., Hu, S. T., Dong, J. E., et al. (2022). Advances in mechanism research on polygonatum in prevention and treatment of diabetes. Front. Pharmacol. 13, 758501. doi:10.3389/fphar.2022.758501

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, H., Li, A., and Gao, K. (2017). Network of conformational transitions revealed by molecular dynamics simulations of the carbonic anhydrase II apo-enzyme. ACS Omega 2 (11), 8414–8420. doi:10.1021/acsomega.7b01414

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahalapbutr, P., Lee, V. S., and Rungrotmongkol, T. (2020). Binding hotspot and activation mechanism of maltitol and lactitol toward the human sweet taste receptor. J. Agric. Food Chem. 68 (30), 7974–7983. doi:10.1021/acs.jafc.0c02580

PubMed Abstract | CrossRef Full Text | Google Scholar

Mark, P., and Nilsson, L. (2001). Structure and dynamics of the TIP3P, SPC, and SPC/E water models at 298 K. J. Phys. Chem. A 105 (43), 9954–9960. doi:10.1021/jp003020w

CrossRef Full Text | Google Scholar

Meldrum, N. U., and Roughton, F. J. (1933). Carbonic anhydrase. Its preparation and properties. J. Physiol. 80 (2), 113–142. doi:10.1113/jphysiol.1933.sp003077

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, B. R., McGee, T. D., Swails, J. M., Homeyer, N., Gohlke, H., and Roitberg, A. E. (2012). MMPBSA.py: an efficient program for end-state free energy calculations. J. Chem. Theory Comput. 8 (9), 3314–3321. doi:10.1021/ct300418h

PubMed Abstract | CrossRef Full Text | Google Scholar

Mu, Y., Nguyen, P. H., and Stock, G. (2005). Energy landscape of a small peptide revealed by dihedral angle principal component analysis. Proteins Struct. Funct. Bioinforma. 58 (1), 45–52. doi:10.1002/prot.20310

PubMed Abstract | CrossRef Full Text | Google Scholar

Nocentini, A., Donald, W. A., and Supuran, C. T. (2004) Carbonic anhydrases: catalytic and inhibition mechanisms, distribution and physiological roles. CRC Press LLC, 1–23.

CrossRef Full Text | Google Scholar

Nocentini, A., Donald, W. A., and Supuran, C. T. (2019). “Human carbonic anhydrases,” in Carbonic anhydrases, 151–185.

CrossRef Full Text | Google Scholar

Occhipinti, R., and Boron, W. F. (2019). Role of carbonic anhydrases and inhibitors in acid-base physiology: insights from mathematical modeling. Int. J. Mol. Sci. 20 (15), 3841. doi:10.3390/ijms20153841

PubMed Abstract | CrossRef Full Text | Google Scholar

Papaleo, E., Mereghetti, P., Fantucci, P., Grandori, R., and De Gioia, L. (2009). Free-energy landscape, principal component analysis, and structural clustering to identify representative conformations from molecular dynamics simulations: the myoglobin case. J. Mol. Graph. Model. 27 (8), 889–899. doi:10.1016/j.jmgm.2009.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Pollard, A., Shephard, F., Freed, J., Liddell, S., and Chakrabarti, L. (2016). Mitochondrial proteomic profiling reveals increased carbonic anhydrase II in aging and neurodegeneration. Aging (Albany NY) 8 (10), 2425–2436. doi:10.18632/aging.101064

PubMed Abstract | CrossRef Full Text | Google Scholar

Poulsen, S.-A., and Davis, R. A. (2014). Natural products that inhibit carbonic anhydrase. Subcell. Biochem. 75, 325–347. doi:10.1007/978-94-007-7359-2_16

PubMed Abstract | CrossRef Full Text | Google Scholar

Qin, B., Jiang, X., Lu, H., Tian, X., Barbault, F., Huang, L., et al. (2010). Diarylaniline derivatives as a distinct class of HIV-1 non-nucleoside reverse transcriptase inhibitors. J. Med. Chem. 53 (13), 4906–4916. doi:10.1021/jm1002952

PubMed Abstract | CrossRef Full Text | Google Scholar

Schäfer, A., Huber, C., and Ahlrichs, R. (1994). Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr. J. Chem. Phys. 100 (8), 5829–5835. doi:10.1063/1.467146

CrossRef Full Text | Google Scholar

Scozzafava, A., Mastrolorenzo, A., and Supuran, C. T. (2006). Carbonic anhydrase inhibitors and activators and their use in therapy. Expert Opin. Ther. Pat. 16 (12), 1627–1664. doi:10.1517/13543776.16.12.1627

CrossRef Full Text | Google Scholar

Shen, F., Xie, P., Li, C., Bian, Z., Wang, X., Peng, D., et al. (2022). Polysaccharides from polygonatum cyrtonema hua reduce depression-like behavior in mice by inhibiting oxidative stress-calpain-1-NLRP3 signaling Axis. Oxid. Med. Cell Longev. 2022, 2566917. doi:10.1155/2022/2566917

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, J., Liu, X., Zhang, S., Li, M., Zhang, Q., and Chen, J. (2023). Molecular insights and optimization strategies for the competitive binding of engineered ACE2 proteins: a multiple replica molecular dynamics study. Phys. Chem. Chem. Phys. 25 (41), 28479–28496. doi:10.1039/d3cp03392a

PubMed Abstract | CrossRef Full Text | Google Scholar

Supuran, C. T. (2008). Carbonic anhydrases: novel therapeutic applications for inhibitors and activators. Nat. Rev. Drug Discov. 7 (2), 168–181. doi:10.1038/nrd2467

PubMed Abstract | CrossRef Full Text | Google Scholar

Supuran, C. T. (2011). Carbonic anhydrase inhibition with natural products: novel chemotypes and inhibition mechanisms. Mol. Divers. 15 (2), 305–316. doi:10.1007/s11030-010-9271-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Supuran, C. T. (2020). Coumarin carbonic anhydrase inhibitors from natural sources. J. Enzyme Inhib. Med. Chem. 35 (1), 1462–1470. doi:10.1080/14756366.2020.1788009

PubMed Abstract | CrossRef Full Text | Google Scholar

Supuran, C. T. (2021). Emerging role of carbonic anhydrase inhibitors. Clin. Sci. (Lond). 135 (10), 1233–1249. doi:10.1042/CS20210040

PubMed Abstract | CrossRef Full Text | Google Scholar

Supuran, C. T. (2022a). Carbonic anhydrase inhibitors from marine natural products. Mar. Drugs 20 (11), 721. doi:10.3390/md20110721

PubMed Abstract | CrossRef Full Text | Google Scholar

Supuran, C. T. (2022b). Anti-obesity carbonic anhydrase inhibitors: challenges and opportunities. J. Enzyme Inhib. Med. Chem. 37 (1), 2478–2488. doi:10.1080/14756366.2022.2121393

PubMed Abstract | CrossRef Full Text | Google Scholar

Supuran, C. T., Capasso, C., and De Simone, G. (2015) Carbonic anhydrase II as target for drug design. Elsevier B.V., 51–90.

CrossRef Full Text | Google Scholar

Supuran, C. T., Casini, A., and Scozzafava, A. (2004) Development of sulfonamide carbonic anhydrase inhibitors. CRC Press LLC, 67–147.

CrossRef Full Text | Google Scholar

Suri, C., and Naik, P. K. (2015). Combined molecular dynamics and continuum solvent approaches (MM-PBSA/GBSA) to predict noscapinoid binding to γ-tubulin dimer. SAR QSAR Environ. Res. 26 (6), 507–519. doi:10.1080/1062936X.2015.1070200

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanchuk, V. Y., Tanin, V. O., Vovk, A. I., and Poda, G. (2016). A new, improved hybrid scoring function for molecular docking and scoring based on AutoDock and AutoDock Vina. Chem. Biol. Drug Des. 87 (4), 618–625. doi:10.1111/cbdd.12697

PubMed Abstract | CrossRef Full Text | Google Scholar

Trott, O., and Olson, A. J. (2010). AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 31 (2), 455–461. doi:10.1002/jcc.21334

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, C., Peng, D., Zhu, J., Zhao, D., Shi, Y., Zhang, S., et al. (2019a). Transcriptome analysis of Polygonatum cyrtonema Hua: identification of genes involved in polysaccharide biosynthesis. Plant Methods 15, 65. doi:10.1186/s13007-019-0441-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Li, G., Zhang, X., Wang, Y., Qiang, Y., Wang, B., et al. (2022). Structural characterization and antioxidant activity of Polygonatum sibiricum polysaccharides. Carbohydr. Polym. 291, 119524. doi:10.1016/j.carbpol.2022.119524

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W.-X., Zhang, X., Dabu, X.-L.-T., He, J., Yang, S.-C., Chen, J.-W., et al. (2019b). Analysis of chemical constituents from Polygonatum cyrtonema after “Nine-Steam-Nine-Bask” processing. Phytochem. Lett. 29, 35–40. doi:10.1016/j.phytol.2018.11.004

CrossRef Full Text | Google Scholar

Wen, C., Wang, H., Wu, X., He, L., Zhou, Q., Wang, F., et al. (2019). ROS-mediated inactivation of the PI3K/AKT pathway is involved in the antigastric cancer effects of thioredoxin reductase-1 inhibitor chaetocin. Cell Death Dis. 10 (11), 809. doi:10.1038/s41419-019-2035-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, P., Zhao, C., Li, X., Gao, Q., Huang, L., Xiao, P., et al. (2018). The genus Polygonatum: a review of ethnopharmacology, phytochemistry and pharmacology. J. Ethnopharmacol. 214, 274–291. doi:10.1016/j.jep.2017.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, S. (2020). Protective effect of Polygonatum sibiricum Polysaccharide on D-galactose-induced aging rats model. Sci. Rep. 10 (1), 2246. doi:10.1038/s41598-020-59055-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Polygonatum sibiricum, small molecule inhibitor, molecular dynamics simulation, molecular docking, inhibition mechanism

Citation: Mu Y, Meng Q, Fan X, Xi S, Xiong Z, Wang Y, Huang Y and Liu Z (2024) Identification of the inhibition mechanism of carbonic anhydrase II by fructooligosaccharides. Front. Mol. Biosci. 11:1398603. doi: 10.3389/fmolb.2024.1398603

Received: 10 March 2024; Accepted: 06 May 2024;
Published: 28 May 2024.

Edited by:

Marco Agostino Deriu, Polytechnic University of Turin, Italy

Reviewed by:

Padhmanand Sudhakar, Kumaraguru College of Technology, India
Yunhui Ge, Alkermes Inc, United States

Copyright © 2024 Mu, Meng, Fan, Xi, Xiong, Wang, Huang and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Qingyang Meng, bWVuZ3F5QHBlY2hvaW4uY29t; Zhen Liu, bGl1emhlbkBlY3VzdC5lZHUuY24=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.