Applying linear interaction energy method for binding affinity calculations of podophyllotoxin analogues with tubulin using continuum solvent model and prediction of cytotoxic activity

Design, synthesis and in silico insights of novel 1,2,3-triazole benzenesulfonamide derivatives as potential carbonic anhydrase IX and XII inhibitors with promising anticancer activity

Novel 1,2,3-triazole benzenesulfonamide derivatives were designed as inhibitors for the tumor- related hCA IX and XII isoforms. Most of the synthesized compounds showed good inhibitory activity against hCA IX and hCA XII isoforms. Compounds 4d, 5h and 6b, exhibited remarkable activity as hCA IX inhibitors, with Ki values in the range of 0.03 to 0.06 µM, more potent than AAZ. Additionally, compounds 5b and 6d, efficiently inhibited hCA XII isoform, with Ki value of 0.02 µM, respectively, similar to AAZ. Further investigation for those potent derivatives against MCF-7, Hep-3B and WI-38 cell lines was achieved. Compounds 4d and 6d exerted dual cytotoxic activity against MCF-7 and Hep-3B cell lines, with IC50 values of 3.35 & 2.12 µM against MCF-7 cell line and 1.72 & 1.56 µM against Hep-3B cell line, with high SI values ranged from 8.92 to 17.38 on both of the cell lines. Besides, they showed a high safety profile against normal human cell line, WI-38. Moreover, compound 5h had better cytotoxic effect on MCF-7 than the reference, DOX, with IC50 value of 4.02 µM. While, compounds 5b and 6b showed higher activity against Hep-3B if compared to the reference drug, 5-FU. From ADME study, compounds 4d, 5b, 6b and 6d obeyed Lipinski's rule of five, and they might be orally active derivatives, while, compound 5h exerted less oral bioavailability than the reference standard acetazolamide. Molecular docking and MDS studies predicted the binding mode and the stability of the target compounds inside hCA IX and hCA XII active sites, especially for compounds 5b and 6b.

Efficacy of Ceftazidime and Cefepime in the Management of COVID-19 Patients: Single Center Report from Egypt

The purpose of this study was to explore the value of using cefepime and ceftazidime in treating patients with COVID-19. A total of 370 (162 males) patients, with RT-PCR-confirmed cases of COVID-19, were included in the study. Out of them, 260 patients were treated with cefepime or ceftazidime, with the addition of steroids to the treatment. Patients were divided into three groups: Group 1: patients treated with cefepime (124 patients); Group 2: patients treated with ceftazidime (136 patients); Group 3 (control group): patients treated according to the WHO guidelines and the Egyptian COVID-19 management protocol (110 patients)/ Each group was classified into three age groups: 18–30, 31–60, and >60 years. The dose of either cefepime or ceftazidime was 1000 mg twice daily for five days. Eight milligrams of dexamethasone were used as the steroidal drug. Careful follow-ups for the patients were carried out. In vitro and in silico M^pro enzyme assays were performed to investigate the antiviral potential of both antibiotics. The mean recovery time for Group 1 was 12 days, for Group 2 was 13 days, and for Group 3 (control) was 19 days. No deaths were recorded, and all patients were recovered without any complications. For Group 1, the recovery time was 10, 12, and 16 days for the age groups 18–30, 30–60, and >60 years, respectively. For Group 2, the recovery time was 11, 13, and 15 days for the age groups 18–30, 30–60, and >60 years, respectively. For Group 3 (control), the recovery time was 15, 16, and 17 days for the age groups 18–30, 30–60, and >60 years, respectively. Both ceftazidime and cefepime showed very good inhibitory activity towards SARS CoV-2′s M^pro, with IC₅₀ values of 1.81 µM and 8.53 µM, respectively. In conclusion, ceftazidime and cefepime are efficient for the management of moderate and severe cases of COVID-19 due to their potential anti-SARS CoV-2 activity and low side effects, and, hence, the currently used complex multidrug treatment protocol can be replaced by the simpler one proposed in this study.

Adequate prediction for inhibitor affinity of Aβ40 protofibril using the linear interaction energy method

The search for efficient inhibitors targeting Aβ oligomers and fibrils is an important issue in Alzheimer's disease treatment. As a consequence, an accurate and computationally cheap approach to estimate the binding affinity for many ligands interacting with Aβ peptides is very important. Here, the calculated binding free energies of 30 ligands interacting with 12Aβ11-40 peptides using the linear interaction energy (LIE) approach are found to be in good correlation with experimental data (R = 0.79). The binding affinities of these complexes are also calculated by using free energy perturbation (FEP) and molecular mechanic/Poisson-Boltzmann surface area (MM/PBSA) methods. The time-consuming FEP method provides results with similar correlation (R = 0.72), whereas MM/PBSA calculations show very low correlation with experimental data (R = 0.27). In all complexes, van der Waals interactions contribute much more than electrostatic interactions. The LIE model, which is much less time-consuming than both the FEP and MM/PBSA methods, opens the door to accurate and rapid affinity prediction of ligands with Aβ peptides and the design of new ligands.

In Silico Methods for the Discovery of Orthosteric GABAB Receptor Compounds

The GABAB receptor (GABAB-R) is a heterodimeric class C G protein-coupled receptor comprised of the GABAB1a/b and GABAB2 subunits. The endogenous orthosteric agonist γ-amino-butyric acid (GABA) binds within the extracellular Venus flytrap (VFT) domain of the GABAB1a/b subunit. The receptor is associated with numerous neurological and neuropsychiatric disorders including learning and memory deficits, depression and anxiety, addiction and epilepsy, and is an interesting target for new drug development. Ligand- and structure-based virtual screening (VS) are used to identify hits in preclinical drug discovery. In the present study, we have evaluated classical ligand-based in silico methods, fingerprinting and pharmacophore mapping and structure-based in silico methods, structure-based pharmacophores, docking and scoring, and linear interaction approximation (LIA) for their aptitude to identify orthosteric GABAB-R compounds. Our results show that the limited number of active compounds and their high structural similarity complicate the use of ligand-based methods. However, by combining ligand-based methods with different structure-based methods active compounds were identified in front of DUDE-E decoys and the number of false positives was reduced, indicating that novel orthosteric GABAB-R compounds may be identified by a combination of ligand-based and structure-based in silico methods.

Tubulin structure-based drug design for the development of novel 4β-sulfur-substituted podophyllum tubulin inhibitors with anti-tumor activity

The well-characterized anti-tubulin agent, podophyllotoxin (PTOX), with the 4′-position methoxyl group, targets the colchicines domain located between α- and β-tubulin. Two guanosine triphosphate (GTP) analogs of the tubulin-binding region were synthesized from PTOX, where a hydroxyl group was substituted with a carbon-sulfur bond. These compounds, 4-MP-PTOX and 4-TG-PTOX, reduce the dosage and greatly improve the therapeutic effect for microtubule damage in cancer cells. Here we characterize the anti-tubulin properties of these compounds. We found the stronger inhibition of tubulin polymerization (the concentration of 50% growth inhibition, GI50 < 2 μM) for compounds 4-TG-PTOX and 4-MP-PTOX, which were better than that of PTOX or colchicine. The cytotoxicity of two designed compounds on tumor cells was also significantly enhanced by comparing to those of PTOX and colchicines. The ΔH value of 4-MP-PTOX and 4-TG-PTOX binding to tubulin by isothermal titration calorimetry (ITC) was found to be −7.4 and −5.3 kcal·mol⁻¹, respectively. The wide range of enthalpy values across the series may reflect entropy/enthalpy compensation effects. Fragments 6-mercaptopurine (MP) and 6-thioguanine (TG) likely enhance the affinity of 4-MP-PTOX and 4-TG-PTOX binding to tubulin by increasing the number of binding sites. The correctness of rational drug design was strictly demonstrated by a bioactivity test.

Accessing the applicability of polarized protein-specific charge in linear interaction energy analysis

The reliability of the linear interaction energy (LIE) depends on the atomic charge model used to delineate the Coulomb interaction between the ligand and its environment. In this work, the polarized protein-specific charge (PPC) implementing a recently proposed fitting scheme has been examined in the LIE calculations of the binding affinities for avidin and β-secretase binding complexes. This charge fitting scheme, termed delta restrained electrostatic potential, bypasses the prevalent numerical difficulty of rank deficiency in electrostatic-potential-based charge fitting methods via a dual-step fitting strategy. A remarkable consistency between the predicted binding affinities and the experimental measurement has been observed. This work serves as a direct evidence of PPC's applicability in rational drug design.

Role of computational methods in pharmaceutical sciences

Over the past two decades computational methods have eased up the financial and experimental burden of early drug discovery process. The in silico methods have provided support in terms of databases, data mining of large genomes, network analysis, systems biology on the bioinformatics front and structure-activity relationship, similarity analysis, docking, and pharmacophore methods for lead design and optimization. This review highlights some of the applications of bioinformatics and chemoinformatics methods that have enriched the field of drug discovery. In addition, the review also provided insights into the use of free energy perturbation methods for efficiently computing binding energy. These in silico methods are complementary and can be easily integrated into the traditional in vitro and in vivo methods to test pharmacological hypothesis.

Rational design, synthesis and biological evaluations of amino-noscapine: a high affinity tubulin-binding noscapinoid

Noscapine and its derivatives are important microtubule-interfering agents shown to have potent anti-tumor activity. The binding free energies (ΔG (bind)) of noscapinoids computed using linear interaction energy (LIE) method with a surface generalized Born (SGB) continuum solvation model were in agreement with the experimental ΔG (bind) with average root mean square error of 0.082 kcal/mol. This LIE-SGB model guided us in designing a novel derivative of noscapine, amino-noscapine [(S)-3-((R)-9-amino-4-methoxy-6-methyl-5,6,7,8-tetrahydro [1, 3] dioxolo[4,5-g]isoquinolin-5-yl)-6,7-dimethoxy isobenzo-furan-1(3H)-one] that has higher tubulin binding activity (predicted ΔG (bind) = -6.438 kcal/mol and experimental ΔG (bind) = -6.628 kcal/mol) than noscapine, but does not significantly change the total extent of the tubulin subunit/polymer ratio. The modes of interaction of amino-noscapine with the binding pocket of tubulin involved three hydrogen bonds and are distinct compared to noscapine which involved only one hydrogen bond. Also the patterns of non-bonded interactions are albeit different between both the lignads. The 'blind docking' approach (docking of ligand with different binding sites of a protein and their evaluations) as well as the reasonable accuracy of calculating ΔG (bind) using LIE-SGB model constitutes the first evidence that this class of compounds binds to tubulin at a site overlapping with colchicine-binding site or close to it. Our results revealed that amino-noscapine has better anti-tumor activity than noscapine.