A selectivity study of benzenesulfonamide derivatives on human carbonic anhydrase II/IX by 3D-QSAR, Molecular Docking and Molecular Dynamics Simulation

Caffeic acid and ferulic acid can improve toxicological damage caused by iron overload mediated by carbonic anhydrase inhibition

The iron ion is an essential element for most forms of life, however, it can damage biological systems when found in free form. Chelation therapy is very important, but it is precarious. Caffeic and ferulic acid are antioxidant compounds with many properties described in research such as anti-inflammatory, antiobesogenic, antithrombotic, vasodilator, and anti-tumor. The aim of the study was to evaluate presenting an in silico approach on the toxicity and bioavailability of caffeic and ferulic acid, subsequently, evaluating them in an iron overload model in vivo and providing a pharmacophoric model through molecular docking. The predictive in silico test did not show relevant toxicity of the compounds, therefore, the in vivo test was performed. The rats received dextran iron and the test groups received caffeic and ferulic acid orally for six weeks. Biochemical, hematological parameters, and tissue oxidative stress marker were analyzed. The experimental model showed increased serum iron levels and changes in several serum parameters such as glucose (215.8 ± 20.3 mg/dL), ALT (512.2 ± 128.7 U/L), creatine kinase (186.8 ± 30.1 U/L), and creatine kinase isoform MB (373.3 ± 69.7 U/L). Caffeic acid and, to a lessed degree, ferullic acid, attenuated the effects of iron overload on the rat serum biochemical parameters. Docking showed a pharmacophoric model where carbonic anhydrase interacted with the test molecules and caffeic acid showed less energy expenditure in this interaction. The results illustrate a new therapeutic action of phenolic compounds on iron overload. The possible interference of carbonic anhydrase in iron metabolism needs to be elucidated.

Synthesis, Characterization, Cytotoxicity, Cellular Imaging, Molecular Docking, and ADMET Studies of Piperazine-Linked 1,8-Naphthalimide-Arylsulfonyl Derivatives

To reduce the mortality and morbidity associated with cancer, new cancer theranostics are in high demand and are an emerging area of research. To achieve this goal, we report the synthesis and characterization of piperazine-linked 1,8-naphthalimide-arylsulfonyl derivatives (SA1-SA7). These compounds were synthesized in good yields following a two-step protocol and characterized using multiple analytical techniques. In vitro cytotoxicity and fluorescent cellular imaging of the compounds were assessed against non-cancerous fibroblast (3T3) and breast cancer (4T1) cell lines. Although the former study indicated the safe nature of the compounds (viability = 82-95% at 1 μg/mL), imaging studies revealed that the designed probes had good membrane permeability and could disperse in the whole cell cytoplasm. In silico studies, including molecular docking, molecular dynamics (MD) simulation, and ADME/Tox results, indicated that the compounds had the ability to target CAIX-expressing cancers. These findings suggest that piperazine-linked 1,8-naphthalimide-arylsulfonyl derivatives are potential candidates for cancer theranostics and a valuable backbone for future research.

3D-QSAR studies on the structure-bitterness analysis of citrus flavonoids

Despite their important bioactivities, the unpleasant bitter taste of citrus derived flavonoids limits their applications in the food industry, and the structure-bitterness relationship of flavonoids is still far from clear. In this study, 26 flavonoids were characterized by their bitterness threshold and their common skeleton using sensory evaluation and molecular superposition, respectively. The quantitative conformational relationship of the structure-bitterness of flavonoids was explored using 3D-QSAR based on comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA). The results showed that increases of a hydrogen bond donor at A-5 or B-3', a bulky group at A-8, or an electron-withdrawing group at B-4' would enhance the bitterness of flavonoids. The bitterness of some flavonoids was predicted and evaluated, and the results were similar to the bitter intensity of the counterparts from the 3D-QSAR and contour plots, confirming the validation of 3D-QSAR. This study explains the theory of the structure-bitterness relationship of flavonoids, by showing potential information for understanding the bitterness in citrus flavonoids and developing a debittering process.

Identification of the structural features of quinazoline derivatives as EGFR inhibitors using 3D-QSAR modeling, molecular docking, molecular dynamics simulations and free energy calculations

Epidermal growth factor receptor (EGFR) is a promising target for the treatment of different types of malignant tumors. Therefore, a combined molecular modeling study was performed on a series of quinazoline derivatives as EGFR inhibitors. The optimum ligand-based CoMFA and CoMSIA models showed reliable and satisfactory predictability (with R²_cv=0.681, R²_ncv=0.844, R²_pred=0.8702 and R²_cv=0.643, R²_ncv=0.874, R²_pred=0.6423). The derived contour maps provide structural features to improve inhibitory activity. Furthermore, the contour maps, molecular docking, and molecular dynamics (MD) simulations have good consistency, illustrating that the derived models are reliable. In addition, MD simulations and binding free energy calculations were also carried out to understand the conformational fluctuations at the binding pocket of the receptor. The results indicate that hydrogen bond, hydrophobic and electrostatic interactions play significant roles on activity and selectivity. Furthermore, amino acids Val31, Lys50, Thr95, Leu149 and Asp160 are considered as essential residues to participate in the ligand-receptor interactions. Overall, this work would offer reliable theoretical basis for future structural modification, design and synthesis of novel EGFR inhibitors with good potency.Communicated by Ramaswamy H. Sarma.

Impact of SNPs, off-targets, and passive permeability on efficacy of BCL6 degrading drugs assigned by virtual screening and 3D-QSAR approach

B-cell lymphoma 6 (BCL6) regulates various genes and is reported to be overexpressed in lymphomas and other malignancies. Thus, BCL6 inhibition or its tagging for degradation would be an amenable therapeutic approach. A library of 2500 approved drugs was employed to find BCL6 inhibitory molecules via virtual screening. Moreover, the 3D core structure of 170 BCL6 inhibitors was used to build a 3D QSAR model and predict the biological activity. The SNP database was analyzed to study the impact on the destabilization of BCL6/drug interactions. Structural similarity search and molecular docking analyses were used to assess the interaction between possible off-targets and BCL6 inhibitors. The tendency of drugs for passive membrane permeability was also analyzed. Lifitegrast (DB11611) had favorable binding properties and biological activity compared to the BI-3802. Missense SNPs were located at the essential interaction sites of the BCL6. Structural similarity search resulted in five BTB-domain containing off-target proteins. BI-3802 and Lifitegrast had similar chemical behavior and binding properties against off-target candidates. More interestingly, the binding affinity of BI-3802 (against off-targets) was higher than Lifitegrast. Energetically, Lifitegrast was less favorable for passive membrane permeability. The interaction between BCL6 and BI-3802 is more prone to SNP-derived variations. On the other hand, higher nonspecific binding of BI-3802 to off-target proteins could bring about higher undesirable properties. It should also be noted that energetically less desirable passive membrane translocation of Lifitegrast would demand drug delivery vehicles. However, further empirical evaluation of Lifitegrast would unveil its true potential.

In the footsteps of an inhibitor unbinding from the active site of human carbonic anhydrase II

The crystal structure of human carbonic anhydrase (HCA) II bound to an inhibitor molecule, 6-hydroxy-2-thioxocoumarin (FC5), shows FC5 to be located in a hydrophobic pocket at the active site. The present work employs classical molecular dynamics (MD) simulation to follow the FC5 molecule for 1 μs as it unbinds from its binding location, adopts the path of substrate/product diffusion (path 1) to leave the active site at around 75 ns. It is then found to undergo repeated binding and unbinding at different locations on the surface of the enzyme in water. Several transient excursions through different regions of the enzyme are also observed prior to its exit from the active site. These transient paths are combined with functionally relevant cavities/channels to enlist five additional pathways (path 2-6). Pathways 1-6 are subsequently explored using steered MD and umbrella sampling simulations. A free energy barrier of 0.969 kcal mol^-1 is encountered along path 1, while barriers in the range of 0.57-2.84 kcal mol^-1 are obtained along paths 2, 4 and 5. We also analyze in detail the interaction between FC5 and the enzyme along each path as the former leaves the active site of HCA II. Our results indicate path 1 to be the major exit pathway for FC5, although competing contributions may also come from the paths 2, 4 and 5.Communicated by Ramaswamy H. Sarma.

In silico design novel (5-imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine derivatives as inhibitors for glycogen synthase kinase 3 based on 3D-QSAR, molecular docking and molecular dynamics simulation

Glycogen Synthase Kinase 3 (GSK-3) is a member of cellular kinase with various functions, such as glucose regulation, cellular differentiation, neuronal function and cell apoptosis. It has been proved as an important therapeutic target in type 2 diabetes mellitus and Alzheimer's disease. To better understand their structure-activity relationships and mechanism of action, an integrated computational study, including three dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, and molecular dynamics (MD), was performed on 79 (5-Imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine GSK-3 inhibitors. In this paper, we constructed 3D-QSAR using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) method. The results showed that the CoMFA model (q² = 0.743,  r² = 0.980) and the CoMSIA model (q² = 0.813,  r² = 0.976) had stable and reliable predictive ability. The electrostatic and H-bond donor fields play important roles in the models. The contour maps of the model visually showed the relationship between the activity of compounds and their three-dimensional structure. Molecular docking was used to identify the key amino acid residues at the active site of GSK-3 and explore its binding mode with ligands. Based on 3D-QSAR models, contour maps and the binding feature between GSK-3 and inhibitor, we designed 10 novel compounds with good potential activity and ADME/T profile. Molecular dynamics simulation results validated that Ile62, Val70 and Lys85 located in the active site play a key role for GSK-3 complexed with inhibitors. These results might provide important information for designing GSK-3 inhibitors with high activity.

Computationally exploring novel xanthine oxidase inhibitors using docking-based 3D-QSAR, molecular dynamics, and virtual screening

Computationally exploring novel potential xanthine oxidase inhibitors using a systematic modeling study.