Vibrational and electronic investigations, thermodynamic parameters, HOMO and LUMO analysis on Lornoxicam by density functional theory

Synthesis, characterization, molecular modeling studies, and biological evaluation of metal piroxicam complexes (M = Ni(II), Pt(IV), Pd(II), Ag(I)) as antibacterial and anticancer agents

Four piroxicam metal complexes; NiL2 , PtL2 , PdL2 , and AgL were synthesized and characterized by different techniques with enhanced antibacterial and anticancer activity. Regarding in vitro antimicrobial activity, complex NiL2 displayed potent antibacterial effect against Escherichia coli and Pseudomonas aeruginosa that was 1.9-folds higher than piroxicam (minimum inhibitory concentration [MIC] = 31.85, 65.32 µM), respectively. In case of G+ve bacteria, complex PtL2 had potent activity on Staphylococcus aureus which was 2.1-folds higher than piroxicam (MIC = 43.12 µM), while activity of complex AgL against Enterococcus faecalis was threefolds higher than piroxicam (MIC = 74.57 µM. Complexes PtL2 and PdL2 exhibited higher inhibition of DNA gyrase than piroxicam (IC50  = 6.21 µM) in the range of 1.9-1.7-folds. The in vitro antiproliferative activity depicted that all investigated complexes showed better cytotoxic effect than piroxicam, specifically Pt and Pd complexes which had lower IC50 values than piroxicam on human liver cancer cell line HepG2 by 1.8 and 1.7-folds, respectively. While Pd and Ag complexes showed 2 and 1.6-folds better effect on human colon cancer cell line HT-29 compared with piroxicam. Molecular modeling studies including docking on Stranded DNA Duplex (1juu) and DNA gyrase enzyme (1kzn) that gave good insight about interaction of complexes with target molecules, calculation of electrostatic potential map and global reactivity descriptors were performed.

Exploring natural products library to identify potential inhibitors targeting isoniazid-resistant tuberculosis

Mycobacterium tuberculosis (MTB) causing tuberculosis (TB) infection is a leading source of illness and death in developing nations, and the emergence of drug-resistant TB remains a significant global threat and a challenge in treating the disease. Mutations in the inhA and katG genes are connected to the principal molecular mechanism of isoniazid (INH) resistance, and continuous treatment of INH for more than a decade led to the evolution of INH resistant-TB (inhR-TB). Structure-based drug discovery approaches on traditional natural compounds are the contemporary source to identify significant lead molecules. This work focuses on discovering effective small compounds from natural compound libraries and applying pharmacophore-based virtual screening to filter out the molecules. The best-identified hit complexes were used for molecular dynamics simulations (MDS) to observe their stability and compactness. A three-dimensional e-pharmacophore hypothesis and screening generated 62 hits based on phase fitness scores from the pharmacophore-based virtual screening. Molecular docking experiments in Maestro's GLIDE module indicated that ZINC000002383126 and ASN22022 may be potential inhibitors of inhA and katG (native, inhA mutants S94A, Y158A, Y158F and Y158S and D137S, Y229F, S315T, W321F, and R418L mutants of katG). In addition, MDS analysis indicated that the native and mutant docked complexes of inhA and katG had good stability and remained compact in the binding pocket of the targets. In vitro studies can further validate the compounds that can act as INH competitive inhibitors.Communicated by Ramaswamy H. Sarma.

Isolation and In Silico Inhibitory Potential against SARS-CoV-2 RNA Polymerase of the Rare Kaempferol 3-O-(6″-O-acetyl)-Glucoside from Calligonum tetrapterum

The phytochemical constituents of Calligonum tetrapterum Jaub. & Spach (Family Polygonaceae) were studied for the first time. The study resulted in the isolation of the rare flavonol glycoside, kaempferol 3-O-(6″-O-acetyl)-glucoside,(K3G-A). The potential inhibitive activity of K3G-A toward SARS-CoV-2 was investigated utilizing several in silico approaches. First, molecular fingerprints and structural similarity experiments were carried out for K3G-A against nine co-crystallized ligands of nine proteins of SARS-CoV-2 to reveal if there is a structural similarity with any of them. The conducted studies showed the high similarity of K3G-A and remdesivir, the co-crystallized ligand of SARS-CoV-2 RNA-dependent RNA polymerase (PDB ID: 7BV2), RdRp. To validate these findings, a DFT study was conducted and confirmed the proposed similarity on the electronic and orbital levels. The binding of K3G-A against RdRp was confirmed through molecular docking studies exhibiting a binding energy of −27.43 kcal/mol, which was higher than that of remdesivir. Moreover, the RdRp-K3G-A complex was subjected to several MD studies at 100 ns that authenticated the accurate mode of binding and the correct dynamic behavior. Finally, in silico ADMET and toxicity evaluation of K3G-A was conducted and denoted the safety and the drug-likeness of K3G-A. In addition to K3G-A, two other metabolites were isolated and identified to be kaempferol (K) and β-sitosterol (β-S).

Isolation and In Silico SARS-CoV-2 Main Protease Inhibition Potential of Jusan Coumarin, a New Dicoumarin from Artemisia glauca

A new dicoumarin, jusan coumarin, (1), has been isolated from Artemisia glauca aerial parts. The chemical structure of jusan coumarin was estimated, by 1D, 2D NMR as well as HR-Ms spectroscopic methods, to be 7-hydroxy-6-methoxy-3-[(2-oxo-2H-chromen-6-yl)oxy]-2H-chromen-2-one. As the first time to be introduced in nature, its potential against SARS-CoV-2 has been estimated using various in silico methods. Molecular similarity and fingerprints experiments have been utilized for 1 against nine co-crystallized ligands of COVID-19 vital proteins. The results declared a great similarity between Jusan Coumarin and X77, the ligand of COVID-19 main protease (PDB ID: 6W63), Mpro. To authenticate the obtained outputs, a DFT experiment was achieved to confirm the similarity of X77 and 1. Consequently, 1 was docked against Mpro. The results clarified that 1 bonded in a correct way inside Mpro active site, with a binding energy of -18.45 kcal/mol. Furthermore, the ADMET and toxicity profiles of 1 were evaluated and showed the safety of 1 and its likeness to be a drug. Finally, to confirm the binding and understand the thermodynamic characters between 1 and Mpro, several molecular dynamics (MD) simulations studies have been administered. Additionally, the known coumarin derivative, 7-isopentenyloxycoumarin (2), has been isolated as well as β-sitosterol (3).

Jusanin, a New Flavonoid from Artemisia commutata with an In Silico Inhibitory Potential against the SARS-CoV-2 Main Protease

A new flavonoid, Jusanin, (1) has been isolated from the aerial parts of Artemisia commutata. The chemical structure of Jusanin has been elucidated using 1D, 2D NMR, and HR-Ms spectroscopic methods to be 5,2',4'-trihydroxy-6,7,5'-trimethoxyflavone. Being new in nature, the inhibition potential of 1 has been estimated against SARS-CoV-2 using different in silico techniques. Firstly, molecular similarity and fingerprint studies have been conducted for Jusanin against co-crystallized ligands of eight different SARS-CoV-2 essential proteins. The studies indicated the similarity between 1 and X77, the co-crystallized ligand SARS-CoV-2 main protease (PDB ID: 6W63). To confirm the obtained results, a DFT study was carried out and indicated the similarity of (total energy, HOMO, LUMO, gap energy, and dipole moment) between 1 and X77. Accordingly, molecular docking studies of 1 against the target enzyme have been achieved and showed that 1 bonded correctly in the protein's active site with a binding energy of -19.54 Kcal/mol. Additionally, in silico ADMET in addition to the toxicity evaluation of Jusanin against seven models have been preceded and indicated the general safety and the likeness of Jusanin to be a drug. Finally, molecular dynamics simulation studies were applied to investigate the dynamic behavior of the Mpro-Jusanin complex and confirmed the correct binding at 100 ns. In addition to 1, three other metabolites have been isolated and identified to be сapillartemisin A (2), methyl-3-[S-hydroxyprenyl]-cumarate (3), and β-sitosterol (4).

In Silico Screening of Semi-Synthesized Compounds as Potential Inhibitors for SARS-CoV-2 Papain-like Protease: Pharmacophoric Features, Molecular Docking, ADMET, Toxicity and DFT Studies

Papain-like protease is an essential enzyme in the proteolytic processing required for the replication of SARS-CoV-2. Accordingly, such an enzyme is an important target for the development of anti-SARS-CoV-2 agents which may reduce the mortality associated with outbreaks of SARS-CoV-2. A set of 69 semi-synthesized molecules that exhibited the structural features of SARS-CoV-2 papain-like protease inhibitors (PLPI) were docked against the coronavirus papain-like protease (PLpro) enzyme (PDB ID: (4OW0). Docking studies showed that derivatives 34 and 58 were better than the co-crystallized ligand while derivatives 17, 28, 31, 40, 41, 43, 47, 54, and 65 exhibited good binding modes and binding free energies. The pharmacokinetic profiling study was conducted according to the four principles of the Lipinski rules and excluded derivative 31. Furthermore, ADMET and toxicity studies showed that derivatives 28, 34, and 47 have the potential to be drugs and have been demonstrated as safe when assessed via seven toxicity models. Finally, comparing the molecular orbital energies and the molecular electrostatic potential maps of 28, 34, and 47 against the co-crystallized ligand in a DFT study indicated that 28 is the most promising candidate to interact with the target receptor (PLpro).

In Silico Exploration of Potential Natural Inhibitors against SARS-Cov-2 nsp10

In continuation of our previous effort, different in silico selection methods were applied to 310 naturally isolated metabolites that exhibited antiviral potentialities before. The applied selection methods aimed to pick the most relevant inhibitor of SARS-CoV-2 nsp10. At first, a structural similarity study against the co-crystallized ligand, S-Adenosyl Methionine (SAM), of SARS-CoV-2 nonstructural protein (nsp10) (PDB ID: 6W4H) was carried out. The similarity analysis culled 30 candidates. Secondly, a fingerprint study against SAM preferred compounds 44, 48, 85, 102, 105, 182, 220, 221, 282, 284, 285, 301, and 302. The docking studies picked 48, 182, 220, 221, and 284. While the ADMET analysis expected the likeness of the five candidates to be drugs, the toxicity study preferred compounds 48 and 182. Finally, a density-functional theory (DFT) study suggested vidarabine (182) to be the most relevant SARS-Cov-2 nsp10 inhibitor.

Docking-Based Virtual Screening and Molecular Dynamics Simulations of Quercetin Analogs as Enoyl-Acyl Carrier Protein Reductase (InhA) Inhibitors of Mycobacterium tuberculosis

The emergence of multidrug-resistant Mycobacterium tuberculosis (MTB) has become a major problem in treating tuberculosis (TB) and shows the need to develop new and efficient drugs for better TB control. This study aimed to use in silico techniques to discover potential inhibitors to the Enoyl-[acyl-carrier-protein] reductase (InhA), which controls mycobacterial cell wall construction. Initially, 391 quercetin analogs present in the KNApSAck_3D database were selected, filters were sequentially applied by docking-based virtual screening. After recategorizing the variables (bond energy prediction and molecular interaction, including hydrogen bond and hydrophobic bond), compounds C00013874, C00006532, and C00013887 were selected as hit ligands. These compounds showed great hydrophobic contributions, and for each hit ligand, 100 ns of molecular dynamic simulations were performed, and the binding free energy was calculated. C00013874 demonstrated the greatest capacity for the InhA enzyme inhibition with ΔGbind = −148.651 kcal/mol compare to NAD (native ligand) presented a ΔGbind = −87.570 kcal/mol. These data are preliminary studies and might be a suitable candidate for further experimental analysis.

New Insights into the Structure and Reactivity of Uracil Derivatives in Different Solvents-A Computational Study

Ab initio calculations were carried out to understand the reactivity and stability of some uracil derivatives, cytosine, 1-methyl cytosine, and cytidine in solvents, water, dimethyl sulfoxide (DMSO), n-octanol, and chloroform. Geometries were fully optimized at MP2 and B3LYP using the 6-31+G(d,p) basis set by applying the Solvation Model on Density (SMD) in solvent systems. The syn conformer of cytidine (cytidine II) is the most stable conformer in the gas phase, while the anticonformer (cytidine IV) is most stable in all of the solvents. Solvation free energy and polarizability values in different solvents decrease in the order water > DMSO > n-octanol > chloroform, while dipole moment, first-order hyperpolarizability, and HOMO-LUMO energy gap values follow the order of polar protic solvent (water and n-octanol) > polar aprotic solvent (DMSO) > nonpolar solvent (chloroform). The solvation free energy, dipole moment, polarizability, and first-order hyperpolarizability values also follow the order of cytosine > 1-methyl cytosine > cytidine. To illustrate that the molecular properties correlate well with the reactivity of the molecules, ab initio calculations were carried out for the reaction of uracil derivatives with Br₂ in the gas phase, water, DMSO, n-octanol, and chloroform. All ground and transition state geometries were fully optimized at B3LYP/6-31+G(d,p), and energies were also calculated at G3MP2 for cytosine and 1-methyl cytosine. For cytosine and 1-methyl cytosine, Gibbs energies of activation decrease with the polarity of the solvent that is chloroform > n-octanol > DMSO > water, while the Gibbs energies of activation for the reaction with cytidine decrease in the order of water > DMSO > n-octanol > chloroform. These results suggest that solvent polarity is very important for the stability and reactivity of uracil derivatives. Hydrogen bonding may also play an important role mainly for cytidine. Free energies of activation decrease with the size of the molecule, i.e., cytosine > 1-methyl cytosine > cytidine.

Physical and chemical stability of a generic etoposide formulation as an alternative to etoposide phosphate

The generic Mylan® etoposide (ETP) has been investigated as an alternative for Etopophos®, in part due to a global shortage of the latter. The generic alternative is different both in its formulation and in its very limited stability (6 h at 25 °C against 4 days for Etopophos®) once reconstituted in ready-to-use chloride or glucose solutions. Its intrinsic stability has been thoroughly studied under various conditions. Two degradation products resulting from hydrolysis were characterized by LC-HR-MSⁿ and supported by density functional theory calculations of the frontier molecular orbitals energies, molecular electrostatic potential mapping, and Mulliken charge analysis. Chemical degradation increases with temperature and can be fitted to a zero order kinetic model with a half-life of 119 days and a kinetic constant of 0.0028 mM day^-1. Precipitation was only observed in solutions at 5 °C and -20 °C indicating that at these temperatures the reconstituted solutions are thermodynamically metastable. In conclusion, ETP at concentrations of 0.68 and 1 mM prepared and stored at 25 °C under good manufacturing practices remained unchanged over a period of 21 days irrespective of the nature of the solvents or the type of container.

Computational investigations of physicochemical, pharmacokinetic, toxicological properties and molecular docking of betulinic acid, a constituent of Corypha taliera (Roxb.) with Phospholipase A2 (PLA2)

BACKGROUND

Betulinic acid (BA) is a natural triterpenoid compound and exhibits a wide range of biological and medicinal properties including anti-inflammatory activity. Therefore, this theoretical investigation is performed to evaluate (a) physicochemical properties such as acid dissociation constant (pKa), distribution coefficient (logD), partition coefficient (logP), aqueous solubility (logS), solvation free energy, dipole moment, polarizability, hyperpolarizability and different reactivity descriptors, (b) pharmacokinetic properties like human intestinal absorption (HIA), cellular permeability, skin permeability (P_Skin), plasma protein binding (PPB), penetration of the blood brain barrier (BBB), (c) toxicological properties including mutagenicity, carcinogenicity, risk of inhibition of hERG gene and (d) molecular mechanism of anti-inflammatory action which will aid the development of analytical method and the synthesis of BA derivatives.

METHODS

The physicochemical properties were calculated using MarvinSketch 15.6.29 and Gaussian 09 software package. The pharmacokinetic and toxicological properties were calculated on online server PreADMET. Further, the molecular docking study was conducted on AutoDock vina in PyRx 0.8.

RESULTS

The aqueous solubility increased with increasing pH due to the ionization of BA leading to decrease in distribution coefficient. The solvation energies in water, dimethyl sulfoxide (DMSO), acetonitrile, n-octanol, chloroform and carbon tetrachloride were - 41.74 kJ/mol, - 53.80 kJ/mol, - 66.27 kJ/mol, - 69.64 kJ/mol, - 65.96 kJ/mol and - 60.13 kJ/mol, respectively. From the results of polarizability and softness, it was clear that BA is less stable and hence, kinetically more reactive in water. BA demonstrated good human intestinal absorption (HIA) and moderate cellular permeability. Further, BA also exhibited positive CNS activity due to high permeability through BBB. The toxicological study revealed that BA was a mutagenic compound but noncarcinogenic in mice model. Moreover, molecular docking study of BA with PLA2 revealed that BA interacts with GLY22 & GLY29 through hydrogen bond formation and LEU2, PHE5, HIS6, ALA17, ALA18, HIS47 and TYR51 through different types of hydrophobic interactions. The binding affinity of BA was - 41.00 kJ/mol which is comparable to the binding affinity of potent inhibitor 6-Phenyl-4(R)-(7-Phenyl-heptanoylamino)-hexanoic acid (BR4) (- 33.89 kJ/mol).

CONCLUSIONS

Our computed properties may assist the development of analytical method to assay BA or to develop BA derivatives with better pharmacokinetic and toxicological profile.