DOX_BDW: Incorporating Solvation and Desolvation Effects of Cavity Water into Nonfitting Protein-Ligand Binding Affinity Prediction.

Accurate prediction of the protein-ligand binding affinity (PLBA) with an affordable cost is one of the ultimate goals in the field of structure-based drug design (SBDD), as well as a great challenge in the computational and theoretical chemistry. Herein, we have systematically addressed the complicated solvation and desolvation effects on the PLBA brought by the difference of the explicit water in the protein cavity before and after ligands bind to the protein-binding site. Based on the new solvation model, a nonfitting method at the first-principles level for the PLBA prediction was developed by taking the bridging and displaced water (BDW) molecules into account simultaneously. The newly developed method, DOX_BDW, was validated against a total of 765 noncovalent and covalent protein-ligand binding pairs, including the CASF2016 core set, Cov_2022 covalent binding testing set, and six testing sets for the hit and lead compound optimization (HLO) simulation. In all of the testing sets, the DOX_BDW method was able to produce PLBA predictions that were strongly correlated with the corresponding experimental data (R = 0.66-0.85). The overall performance of DOX_BDW is better than the current empirical scoring functions that are heavily parameterized. DOX_BDW is particularly outstanding for the covalent binding situation, implying the need for considering an electronic structure in covalent drug design. Furthermore, the method is especially recommended to be used in the HLO scenario of SBDD, where hundreds of similar derivatives need to be screened and refined. The computational cost of DOX_BDW is affordable, and its accuracy is remarkable.

Cov_FB3D: A De Novo Covalent Drug Design Protocol Integrating the BA-SAMP Strategy and Machine-Learning-Based Synthetic Tractability Evaluation.

De novo drug design actively seeks to use sets of chemical rules for the fast and efficient identification of structurally new chemotypes with the desired set of biological properties. Fragment-based de novo design tools have been successfully applied in the discovery of noncovalent inhibitors. Nevertheless, these tools are rarely applied in the field of covalent inhibitor design. Herein, we present a new protocol, called Cov_FB3D, which involves the in silico assembly of potential novel covalent inhibitors by identifying the active fragments in the covalently binding site of the target protein. In this protocol, we propose a BA-SAMP strategy, which combines the noncovalent moiety score with the X-Score as the molecular mechanism (MM) level, and the covalent candidate score with the PM7 as the QM level. The synthetic accessibility of each suggested compound could be further evaluated with machine-learning-based synthetic complexity evaluation (SCScore). An in-depth test of this protocol against the crystal structures of 15 covalent complexes consisting of BTK inhibitors, KRAS inhibitors, EGFR inhibitors, EphB1 inhibitors, MAGL inhibitors, and MAPK inhibitors revealed that most of these inhibitors could be de novo reproduced from the fragments by Cov_FB3D. The binding modes of most generated reference poses could accurately reproduce the known binding mode of most of the reference covalent adduct in the binding site (RMSD ≤ 2 Å). In particular, most of these inhibitors were ranked in the top 2%, using the BA-SAMP strategy. Notably, the novel human ALDOA inhibitor (T1) with potent inhibitory activity (0.34 ± 0.03 µM) and greater synthetic accessibility was successfully de novo designed by this protocol. The positive results confirm the abilities of Cov_FB3D protocol.

Biological evaluation and SAR analysis of novel covalent inhibitors against fructose-1,6-bisphosphatase.

Fructose-1,6-bisphosphatase (FBPase) is an attractive target for affecting the GNG pathway. In our previous study, the C128 site of FBPase has been identified as a new allosteric site, where several nitrovinyl compounds can bind to inhibit FBPase activity. Herein, a series of nitrostyrene derivatives were further synthesized, and their inhibitory activities against FBPase were investigated in vitro. Most of the prepared nitrostyrene compounds exhibit potent FBPase inhibition (IC50 < 10 µM). Specifically, when the substituents of F, Cl, OCH3, CF3, OH, COOH, or 2-nitrovinyl were installed at the R2 (meta-) position of the benzene ring, the FBPase inhibitory activities of the resulting compounds increased 4.5-55 folds compared to those compounds with the same groups at the R1 (para-) position. In addition, the preferred substituents at the R3 position were Cl or Br, thus compound HS36 exhibited the most potent inhibitory activity (IC50 = 0.15 µM). The molecular docking and site-directed mutation suggest that C128 and N125 are essential for the binding of HS36 and FBPase, which is consistent with the C128-N125-S123 allosteric inhibition mechanism. The reaction enthalpy calculations show that the order of the reactions of compounds with thiol groups at the R3 position is Cl > H > CH3. CoMSIA analysis is consistent with our proposed binding mode. The effect of compounds HS12 and HS36 on glucose production in primary mouse hepatocytes were further evaluated, showing that the inhibition was 71% and 41% at 100 µM, respectively.

Design, synthesis and algicides activities of thiourea derivatives as the novel scaffold aldolase inhibitors.

By using a new Fragment-Based Virtual Screen strategy, two series of novel FBA-II inhibitors (thiourea derivatives) were de novo discovered based on the active site of fructose-1, 6-bisphosphate aldolase from Cyanobacterial (CyFBA). In comparison, most of the N-(2-benzoylhydrazine-1-carbonothioyl) benzamide derivatives (L14∼L22) exhibit higher CyFBA-II inhibitory activities compared to N-(phenylcarbamothioyl) benzamide derivatives (L1∼L13). Especially, compound L14 not only shows higher CyFBA-II activity (Ki = 0.65 µM), but also exhibits most potent in vivo activity against Synechocystis sp. PCC 6803 (EC50 = 0.09 ppm), higher (7-fold) than that of our previous inhibitor (EC50 = 0.6 ppm). The binding modes of compound L14 and CyFBA-II were further elucidated by jointly using DOX computational protocol, MM-PBSA and site-directed mutagenesis assays. The positive results suggest that strategy adopted in this study was promising to rapidly discovery the potent inhibitors with novel scaffolds. The satisfactory algicide activities suggest that the thiourea derivatives is very likely to be a promising lead for the development of novel specific algicides to solve Cyanobacterial harmful algal blooms (CHABs).

Natural products as sources of new fungicides (IV): Synthesis and biological evaluation of isobutyrophenone analogs as potential inhibitors of class-II fructose-1,6-bisphosphate aldolase.

Several recently identified antifungal compounds share the backbone structure of acetophenones. The aim of the present study was to develop new isobutyrophenone analogs as new antifungal agents. A series of new 2,4-dihydroxy-5-methyl isobutyrophenone derivatives were prepared and characterized by 1H, 13C NMR and MS spectroscopic data. These products were evaluated for in vitro antifungal activities against seven plant fungal pathogens by the mycelial growth inhibitory rate assay. Compounds 3, 4a, 5a, 5b, 5e, 5f and 5g showed a broad-spectrum high antifungal activity. On the other hand, for the first time, these compounds were also assayed as potential inhibitors against Class II fructose-1,6-bisphosphate aldolase (Fba) from the rice blast fungus, Magnaporthe grisea. Compounds 5e and 5g were found to exhibit the inhibition constants (Ki) for 15.12 and 14.27 µM, respectively, as the strongest competitive inhibitors against Fba activity. The possible binding-modes of compounds 5e and 5g were further analyzed by molecular docking algorithms. The results strongly suggested that compound 5g could be a promising lead for the discovery of new fungicides via targeting Class II Fba.

Structure-Based Rational Design of Novel Inhibitors Against Fructose-1,6-Bisphosphate Aldolase from Candida albicans.

Class II fructose-1,6-bisphosphate aldolases (FBA-II) are attractive new targets for the discovery of drugs to combat invasive fungal infection, because they are absent in animals and higher plants. Although several FBA-II inhibitors have been reported, none of these inhibitors exhibit antifungal effect so far. In this study, several novel inhibitors of FBA-II from C. albicans (Ca-FBA-II) with potent antifungal effects were rationally designed by jointly using a specific protocols of molecular docking-based virtual screening, accurate binding-conformation evaluation strategy, synthesis and enzymatic assays. The enzymatic assays reveal that the compounds 3c, 3e-g, 3j and 3k exhibit high inhibitory activity against Ca-FBA-II (IC50 < 10 µM), and the most potential inhibitor is 3g, with IC50 value of 2.7 µM. Importantly, the compounds 3f, 3g, and 3l possess not only high inhibitions against Ca-FBA-II, but also moderate antifungal activities against C. glabrata (MIC80 = 4-64 µg/mL). The compounds 3g, 3l, and 3k in combination with fluconazole (8 µg/mL) displayed significantly synergistic antifungal activities (MIC80 < 0.0625 µg/mL) against resistant Candida strains, which are resistant to azoles drugs. The probable binding modes between 3g and the active site of Ca-FBA-II have been proposed by using the DOX (docking, ONIOM, and XO) strategy. To our knowledge, no FBA-II inhibitors with antifungal activities against wild type and resistant strains from Candida were reported previously. The positive results suggest that the strategy adopted in this study are a promising method for the discovery of novel drugs against azole-resistant fungal pathogens in the future.

Design and optimization of N-acylhydrazone pyrimidine derivatives as E. coli PDHc E1 inhibitors: Structure-activity relationship analysis, biological evaluation and molecular docking study.

By targeting the thiamin diphosphate (ThDP) binding site of Escherichia coli (E. coli) pyruvate dehydrogenase multienzyme complex E1 (PDHc E1), a series of novel 'open-chain' classes of ThDP analogs A, B, and C with N-acylhydrazone moieties was designed and synthesized to explore their activities against E. coli PHDc E1 in vitro and their inhibitory activity against microbial diseases were further evaluated in vivo. As a result, A1-23 exhibited moderate to potent inhibitory activities against E. coli PDHc E1 (IC50=0.15-23.55µM). The potent inhibitors A13, A14, A15, C2, had strong inhibitory activities with IC50 values of 0.60, 0.15, 0.39 and 0.34µM against E. coli PDHc E1 and with good enzyme-selective inhibition between microorganisms and mammals. Especially, the most powerful inhibitor A14 could 99.37% control Xanthimonas oryzae pv. Oryzae. Furthermore, the binding features of compound A14 within E. coli PDHc E1 were investigated to provide useful insights for the further construction of new inhibitor by molecular docking, site-directed mutagenesis, and enzymatic assays. The results indicated that A14 had most powerful inhibition against E. coli PDHc E1 due to the establishment of stronger interaction with Glu571, Met194, Glu522, Leu264 and Phe602 at active site of E.coli PDHc E1. It could be used as a lead compound for further optimization, and may have potential as a new microbicide.

DOX: A new computational protocol for accurate prediction of the protein-ligand binding structures.

Molecular docking techniques have now been widely used to predict the protein-ligand binding modes, especially when the structures of crystal complexes are not available. Most docking algorithms are able to effectively generate and rank a large number of probable binding poses. However, it is hard for them to accurately evaluate these poses and identify the most accurate binding structure. In this study, we first examined the performance of some docking programs, based on a testing set made of 15 crystal complexes with drug statins for the human 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR). We found that most of the top ranking HMGR-statin binding poses, predicted by the docking programs, were energetically unstable as revealed by the high theoretical-level calculations, which were usually accompanied by the large deviations from the geometric parameters of the corresponding crystal binding structures. Subsequently, we proposed a new computational protocol, DOX, based on the joint use of molecular Docking, ONIOM, and eXtended ONIOM (XO) methods to predict the accurate binding structures for the protein-ligand complexes of interest. Our testing results demonstrate that the DOX protocol can efficiently predict accurate geometries for all 15 HMGR-statin crystal complexes without exception. This study suggests a promising computational route, as an effective alternative to the experimental one, toward predicting the accurate binding structures, which is the prerequisite for all the deep understandings of the properties, functions, and mechanisms of the protein-ligand complexes.

A Rational Design, Synthesis, Biological Evaluation and Structure--Activity Relationship Study of Novel Inhibitors against Cyanobacterial Fructose-1,6-bisphosphate Aldolase.

In the present study, a series of novel maleimide derivatives were rationally designed and optimized, and their inhibitory activities against cyanobacteria class-II fructose-1,6-bisphosphate aldolase (Cy-FBA-II) and Synechocystis sp. PCC 6803 were further evaluated. The experimental results showed that the introduction of a bigger group (Br, Cl, CH3, or C6H3-o-F) on the pyrrole-2',5'-dione ring resulted in a decrease in the Cy-FBA-II inhibitory activity of the hit compounds. Generally, most of the hit compounds with high Cy-FBA-II inhibitory activities could also exhibit high in vivo activities against Synechocystis sp. PCC 6803. Especially, compound 10 not only shows a high Cy-FBA-II activity (IC50 = 1.7 µM) but also has the highest in vivo activity against Synechocystis sp. PCC 6803 (EC50 = 0.6 ppm). Thus, compound 10 was selected as a representative molecule, and its probable interactions with the surrounding important residues in the active site of Cy-FBA-II were elucidated by the joint use of molecular docking, molecular dynamics simulations, ONIOM calculations, and enzymatic assays to provide new insight into the binding mode of the inhibitors and Cy-FBA-II. The positive results indicate that the design strategy used in the present study is very likely to be a promising way to find novel lead compounds with high inhibitory activities against Cy-FBA-II in the future. The enzymatic and algal inhibition assays suggest that Cy-FBA-II is very likely to be a promising target for the design, synthesis, and development of novel specific algicides to solve cyanobacterial harmful algal blooms.

Rational design, synthesis and biological evaluation of 1,3,4-oxadiazole pyrimidine derivatives as novel pyruvate dehydrogenase complex E1 inhibitors.

On the basis of previous study on 2-methylpyrimidine-4-ylamine derivatives I, further synthetic optimization was done to find potent PDHc-E1 inhibitors with antibacterial activity. Three series of novel pyrimidine derivatives 6, 11 and 14 were designed and synthesized as potential Escherichia coli PDHc-E1 inhibitors by introducing 1,3,4-oxadiazole-thioether, 2,4-disubstituted-1,3-thiazole or 1,2,4-triazol-4-amine-thioether moiety into lead structure I, respectively. Most of 6, 11 and 14 exhibited good inhibitory activity against E. coli PHDc-E1 (IC50 0.97-19.21 µM) and obvious inhibitory activity against cyanobacteria (EC50 0.83-9.86 µM). Their inhibitory activities were much higher than that of lead structure I. 11 showed more potent inhibitory activity against both E. coli PDHc-E1 (IC50<6.62 µM) and cyanobacteria (EC50<1.63 µM) than that of 6, 14 or lead compound I. The most effective compound 11d with good enzyme-selectivity exhibited most powerful inhibitory potency against E. coli PDHc-E1 (IC50=0.97 µM) and cyanobacteria (EC50=0.83 µM). The possible interactions of the important residues of PDHc-E1 with title compounds were studied by molecular docking, site-directed mutagenesis, and enzymatic assays. The results indicated that 11d had more potent inhibitory activity than that of 14d or I due to its 1,3,4-oxadiazole moiety with more binding position and stronger interaction with Lsy392 and His106 at active site of E. coli PDHc-E1.

Synthesis and biological evaluation of novel inhibitors against 1,3,8-trihydroxynaphthalene reductase from Magnaporthe grisea.

1,3,8-Trihydroxynaphthalene reductase (3HNR) is an essential enzymes that is involved in fungal melanin biosynthesis. Based on the structural informations of active site of 3HNR, a series of ß-nitrostyrene compounds were rationally designed and synthesized. The enzymatic activities of these compounds showed that most of them exhibited high inhibitory activities (<5.0 µM) against 3HNR; compound 3-2 exhibit the highest inhibitory activity (IC50=0.29 µM). In particular, some of these compounds had moderate fungicidal activity against Magnaporthe grisea. Compound 3-4 showed high in vivo activities against M. grisea (EC50=9.5 ppm). Furthermore, compound 3-2 was selected as a representative molecule, and the probable binding mode of this compound and the surrounding residues in the active site of 3HNR was elucidated by using molecular dock. The positive results suggest that ß-nitrostyrene derivatives are most likely to be promising leads toward the discovery of novel agent of rice blast.

Clarifying and illustrating the electronic energy transfer pathways in trimeric and hexameric aggregation state of cyanobacteria allophycocyanin within the framework of Förster theory.

Within the framework of the Förster theory, the electronic excitation energy transfer pathways in the cyanobacteria allophycocyanin (APC) trimer and hexamer were studied. The associated physical quantities (i.e., excitation energy, oscillator strength, and transition dipole moments) of the phycocyanobilins (PCBs) located in APC were calculated at time-dependent density functional theory (TDDFT) level of theory. To estimate the influence of protein environment on the preceding calculated physical quantities, the long-range interactions were approximately considered with the polarizable continuum model at the TDDFT level of theory, and the short-range interaction caused by surrounding aspartate residue of PCBs were taken into account as well. The shortest energy transfer time calculated in the framework of the Förster model at TDDFT/B3LYP/6-31+G* level of theory are about 0.10 ps in the APC trimer and about 170 ps in the APC monomer, which are in qualitative agreement with the experimental finding that a very fast lifetime of 0.43-0.44 ps in APC trimers, whereas its monomers lacked any corresponding lifetime. These results suggest that the lifetime of 0.43-0.44 ps in the APC trimers determined by Sharkov et al. was most likely attributed to the energy transfer of α(1) -84 ↔ ß(3) -84 (0.23 ps), ß(1) -84 ↔ α(2) -84 (0.11 ps) or ß(2) -84 ↔ α(3) -84 (0.10 ps). So far, no experimental or theoretical energy transfer rates between two APC trimmers were reported, our calculations predict that the predominate energy transfer pathway between APC trimers is likely to occur from α(3) -84 in one trimer to α(5) -84 in an adjacent trimer with a rate of 32.51 ps.

Design, synthesis and evaluation of novel 5-phenylpyridin-2(1H)-one derivatives as potent reversible Bruton's tyrosine kinase inhibitors.

A series of novel reversible Btk inhibitors has been designed based on the structure of the recently reported preclinical drug RN486. The synthesis and SAR of these compounds are described. Among these derivatives, compound 16b was identified to be a potent and orally available reversible agent with satisfactory Btk enzymatic and cellular inhibition in vitro, as well as favorable PK properties and inhibition of arthritis in vivo.

The design, synthesis and biological evaluation of novel thiamin diphosphate analog inhibitors against the pyruvate dehydrogenase multienzyme complex E1 from Escherichia coli.

Pyruvate dehydrogenase multienzyme complex E1 (PDHc E1) is a potential target enzyme when looking for inhibitors to combat microbial disease. In this study, we designed and synthesized a series of novel thiamin diphosphate (ThDP) analogs with triazole ring and oxime ether moieties as potential inhibitors of PDHc E1. Their inhibitory activities against PDHc E1 were examined both in vitro and in vivo. Most of the tested compounds exhibited moderate inhibitory activities against PDHc E1 (IC50 = 6.1-75.5 µM). The potent inhibitors 4g, 4h and 4j, had strong inhibitory activities with IC50 values of 6.7, 6.9 and 6.1 µM against PDHc E1 in vitro and with inhibition rates of 35%, 50% and 33% at 100 µg mL(-1) against Gibberella zeae in vivo, respectively. The binding mode of 4j to PDHc E1 was analyzed by a molecular docking method. Furthermore, the possible interactions of the important residues of PDHc E1 with compound 4j were examined by site-directed mutagenesis, enzymatic assays and spectral fluorescence studies. The theoretical and experimental results are in good agreement and suggest that compound 4j could be used as a lead compound for further optimization, and may have potential as a new microbicide.

Pharmacophore-based virtual screening and experimental validation of novel inhibitors against cyanobacterial fructose-1,6-/sedoheptulose-1,7-bisphosphatase.

Cyanobacterial fructose-1,6-/sedoheptulose-1,7-bisphoshatase (cy-FBP/SBPase) is a potential enzymatic target for screening of novel inhibitors that can combat harmful algal blooms. In the present study, we targeted the substrate binding pocket of cy-FBP/SBPase. A series of novel hit compounds from the SPECs database were selected by using a pharmacophore-based virtual screening strategy. Most of the compounds tested exhibited moderate inhibitory activities (IC50 = 20.7-176.9 µM) against cy-FBP/SBPase. Compound 2 and its analogues 10 and 11 exhibited strong inhibitory activities, with IC50 values of 20.7, 13.4, and 19.0 µM against cy-FBP/SBPase in vitro and EC50 values of 12.3, 10.9, and 2.9 ppm against cyanobacteria Synechocystis PCC6803 in vivo, respectively. The compound 10 was selected in order to perform a refined docking study to investigate the rational binding mode of inhibitors with cy-FBP/SBPase. Furthermore, possible interactions of the residues with inhibitors were examined by site-directed mutagenesis, enzymatic assays, and fluorescence spectral analyses. The results provide insight into the binding mode between the inhibitors and the substrate binding pocket. The observed theoretical and experimental results are in concert, indicating that the modeling strategies and screening methods employed are appropriate to search for novel lead compounds having both structural diversity and high inhibitory activity against cy-FBP/SBPase.

Design, synthesis and molecular docking of amide and urea derivatives as Escherichia coli PDHc-E1 inhibitors.

By targeting the ThDP binding site of Escherichia coli PDHc-E1, two new 'open-chain' classes of E. coli PDHc-E1 inhibitors, amide and urea derivatives, were designed, synthesized, and evaluated. The amide derivatives of compound 6d, with 4-NO2 in the benzene ring, showed the most potent inhibition of E. coli PDHc-E1. The urea derivatives displayed more potent inhibitory activity than the corresponding amide derivatives with the same substituent. Molecular docking studies confirmed that the urea derivatives have more potency due to the two hydrogen bonds formed by two NH of urea with Glu522. The docking results also indicate it might help us to design more efficient PDHc-E1 inhibitors that could interact with Glu522.

Design, synthesis and molecular modeling of novel N-acylhydrazone derivatives as pyruvate dehydrogenase complex E1 inhibitors.

As potential inhibitors of pyruvate dehydrogenase complex E1 (PDHc-E1), a series of 19 1-((4-amino-2-methylpyrimidin-5-yl)methyl)-5-methyl-N'-(substituent)benzylidene-1H-1,2,3-triazole-4-carbohydrazide 4 has been synthesized and tested for their PDHc-E1 inhibitory activity in vitro. Some of these compounds such as 4a, 4g, 4l, 4o, 4p, and 4q were demonstrated to be effective inhibitors by the bioassay of Escherichia coli PDHc-E1. SAR analysis indicated that the PDHc-E1 inhibitory activity could be further enhanced by optimizing the substituted groups in the parent compound. Molecular modeling study with compound 4o as a model was performed to evaluate docking. The results of modeling study suggested a probable inhibition mechanism.

Role of Argonaute proteins in RNAi pathway in Plutella xylostella: A review.

Argonaute (Ago) proteins act as key elements in RNA interference (RNAi) pathway, orchestrating the intricate machinery of gene regulation within eukaryotic cells. Within the RNAi pathway, small RNA molecules, including microRNA (miRNA), small interfering RNA (siRNA), and PIWI-interacting RNA (piRNA), collaborate with Ago family member proteins such as Ago1, Ago2, and Ago3 to form the RNA-induced silencing complex (RISC). This RISC complex, in turn, either cleaves the target mRNA or inhibits the process of protein translation. The precise contributions of Ago proteins have been well-established in numerous animals and plants, although they still remain unclear in some insect species. This review aims to shed light on the specific roles played by Ago proteins within the RNAi mechanism in a destructive lepidopteran pest, the diamondback moth (Plutella xylostella). Furthermore, we explore the potential of double-stranded RNA (dsRNA)-mediated RNAi as a robust genetic tool in pest management strategies. Through an in-depth examination of Ago proteins and dsRNA-mediated RNAi, this review seeks to contribute to our understanding of innovative approaches for controlling this pest and potentially other insect species of agricultural significance.

Target Fishing Reveals a Novel Mechanism of N-Acylamino Saccharin Derivatives Targeting Glyceraldehyde-3-Phosphate Dehydrogenase toward Cyanobacterial Blooms Control.

Based on current challenges of poor targeting and limited choices in chemical control methods of cyanobacterial blooms (CBs), identifying new targets is an urgent and formidable task in the quest for target-based algaecides. This study discovered N-acylamino saccharin derivatives exhibiting potent algicidal activity. Thus, using N-acylamino saccharin as the probes, glyceraldehyde-3-phosphate dehydrogenase from cyanobacterial (CyGAPDH) was identified as a new target of algaecides through the activity-based protein profiling (ABPP) strategy for the first time. Building upon the structure of Probe2, a series of derivatives were designed and synthesized, with compound b6 demonstrating the most potent inhibitory activity against CyGAPDH and Synechocystis sp. PCC6803 (IC50 = 1.67 µM and EC50 = 1.15 µM). Furthermore, the potential covalent binding model of b6 to the cysteine residue C154 was explored through covalent possibility prediction, LC-MS experiments, substrate competitive inhibition experiments, and molecular docking. Especially, the results revealed C154 as a crucial covalent binding site, with residues T184 and R11 forming robust hydrophobic interactions and H181 establishing significant hydrogen-bonding interactions with b6, highlighting their potential as essential pharmacophores. In summary, this study not only identifies a novel target of algaecides for the control of CB but also lays the solid foundation for the development of targeted covalent algaecides.

Understanding the electronic energy transfer pathways in the trimeric and hexameric aggregation state of cyanobacteria phycocyanin within the framework of Förster theory.

In the present study, the electronic energy transfer pathways in trimeric and hexameric aggregation state of cyanobacteria C-phycocyanin (C-PC) were investigated in term of the Förster theory. The corresponding excited states and transition dipole moments of phycocyanobilins (PCBs) located into C-PC were examined by model chemistry in gas phase at time-dependent density functional theory (TDDFT), configuration interaction-singles (CIS), and Zerner's intermediate neglect of differential overlap (ZINDO) levels, respectively. Then, the long-range pigment-protein interactions were approximately taken into account by using polarizable continuum model (PCM) at TDDFT level to estimate the influence of protein environment on the preceding calculated physical quantities. The influence of the short-range interaction caused by aspartate residue nearby PCBs was examined as well. Only when the protonation of PCBs and its long- and short-range interactions were properly taken into account, the calculated energy transfer rates (1/K) in the framework of Förster model at TDDFT/B3LYP/6-31+G* level were in good agreement with the experimental results of C-PC monomer and trimer. Furthermore, the present calculated results suggested that the energy transfer pathway in C-PC monomer is predominant from ß-155 to ß-84 (1/K = 13.4 ps), however, from α-84 of one monomer to ß-84 (1/K = 0.3-0.4 ps) in a neighbor monomer in C-PC trimer. In C-PC hexamer, an additional energy flow was predicted to be from ß-155 (or α-84) in top trimer to adjacent ß-155 (or α-84) (1/K = 0.5-2.7 ps) in bottom trimer.