Development of Metal Complexes for Treatment of Coronaviruses

Coronavirus disease (SARS-CoV-2) is a global epidemic. This pandemic, which has been linked to high rates of death, has forced some countries throughout the world to implement complete lockdowns in order to contain the spread of infection. Because of the advent of new coronavirus variants, it is critical to find effective treatments and vaccines to prevent the virus’s rapid spread over the world. In this regard, metal complexes have attained immense interest as antibody modifiers and antiviral therapies, and they have a lot of promise towards SARS-CoV-2 and their suggested mechanisms of action are discussed, i.e., a new series of metal complexes’ medicinal vital role in treatment of specific proteins or SARS-CoV-2 are described. The structures of the obtained metal complexes were fully elucidated by different analytical and spectroscopic techniques also. Molecular docking and pharmacophore studies presented that most of complexes studied influenced good binding affinity to the main protease SARS-CoV-2, which also was attained as from the RCSB pdb (Protein Data Bank) data PDB ID: 6 W41, to expect the action of metal complexes in contradiction of COVID-19. Experimental research is required to determine the pharmacokinetics of most of the complexes analyzed for the treatment of SARS-CoV-2-related disease. Finally, the toxicity of a metal-containing inorganic complex will thus be discussed by its capability to transfer metals which may bind with targeted site.     

In Vitro and In Silico Evaluation of New 1,3,4-Oxadiazole Derivatives of Pyrrolo[3,4-d]pyridazinone as Promising Cyclooxygenase Inhibitors

Since long-term use of classic NSAIDs can cause severe side effects related mainly to the gastroduodenal tract, discovery of novel cyclooxygenase inhibitors with a safe gastric profile still remains a crucial challenge. Based on the most recent literature data and previous own studies, we decided to modify the structure of already reported 1,3,4-oxadiazole based derivatives of pyrrolo[3,4-d]pyridazinone in order to obtain effective COX inhibitors. Herein we present the synthesis, biological evaluation and molecular docking studies of 12 novel compounds with disubstituted arylpiperazine pharmacophore linked in a different way with 1,3,4-oxadiazole ring. None of the obtained molecules show cytotoxicity on NHDF and THP-1 cell lines and, therefore, all were qualified for further investigation. In vitro cyclooxygenase inhibition assay revealed almost equal activity of new derivatives towards both COX-1 and COX-2 isoenzymes. Moreover, all compounds inhibit COX-2 isoform better than Meloxicam which was used as reference. Anti-inflammatory activity was confirmed in biological assays according to which title molecules are able to reduce induced inflammation within cells. Molecular docking studies were performed to describe the binding mode of new structures to cyclooxygenase. Investigated derivatives take place in the active site of COX, very similar to Meloxicam. For some compounds, promising druglikeness was calculated using in silico predictions.     

蛋清中ACE抑制肽的筛选及其作用机制

表征来自卵清蛋白中的血管紧张素转换酶（angiotensin converting enzyme,ACE）抑制肽。基于在线程序,采用胃肠道蛋白酶对卵清蛋白进行模拟酶切。对获得的三肽毒性、溶解性和ADME性质进行预测,以ACE为蛋白靶标,筛选出能够与ACE紧密结合的三肽,对这些三肽的活性进行验证。结果表明：获得了一种新的ACE抑制三肽CIK,其IC50值为（161±0.06）μmol/L。CIK-ACE复合物通过5 个常规氢键,2 个碳氢键和4 个盐桥达到稳定状态。CIK最佳的对接构型可以与ACE的氨基酸残基Gln281、His353、Ala354、Glu376、Val380、His383、Glu384、Lys511、His513、Tyr523和Phe527形成相互作用。最佳的药效团模型由3 个氢键受体,一个疏水中心和一个正电荷中心组成,三肽CIK可与其中4 个特征相匹配。     

Identification of Selective BRD9 Inhibitor via Integrated Computational Approach

Bromodomain-containing protein 9 (BRD9), a member of the bromodomain and extra terminal domain (BET) protein family, works as an epigenetic reader. BRD9 has been considered an essential drug target for cancer, inflammatory diseases, and metabolic disorders. Due to its high similarity among other isoforms, no effective treatment of BRD9-associated disorders is available. For the first time, we performed a detailed comparative analysis among BRD9, BRD7, and BRD4. The results indicate that residues His42, Gly43, Ala46, Ala54, Val105, and Leu109 can confer the BRD9 isoform selectivity. The predicted crucial residues were further studied. The pharmacophore model’s features were precisely mapped with some key residues including, Gly43, Phe44, Phe45, Asn100, and Tyr106, all of which play a crucial role in BRD9 inhibition. Docking-based virtual screening was utilized with the consideration of the conserved water network in the binding cavity to identify the potential inhibitors of BRD9. In this workflow, 714 compounds were shortlisted. To attain selectivity, 109 compounds were re-docked to BRD7 for negative selection. Finally, four compounds were selected for molecular dynamics studies. Our studies pave the way for the identification of new compounds and their role in causing noticeable, functional differences in isoforms and between orthologues.     

In Silico Identification of Multi-Target Ligands as Promising Hit Compounds for Neurodegenerative Diseases Drug Development

The conventional treatment of neurodegenerative diseases (NDDs) is based on the “one molecule—one target” paradigm. To combat the multifactorial nature of NDDs, the focus is now shifted toward the development of small-molecule-based compounds that can modulate more than one protein target, known as “multi-target-directed ligands” (MTDLs), while having low affinity for proteins that are irrelevant for the therapy. The in silico approaches have demonstrated a potential to be a suitable tool for the identification of MTDLs as promising drug candidates with reduction in cost and time for research and development. In this study more than 650,000 compounds were screened by a series of in silico approaches to identify drug-like compounds with predicted activity simultaneously towards three important proteins in the NDDs symptomatic treatment: acetylcholinesterase (AChE), histone deacetylase 2 (HDAC2), and monoamine oxidase B (MAO-B). The compounds with affinities below 5.0 µM for all studied targets were additionally filtered to remove known non-specifically binding or unstable compounds. The selected four hits underwent subsequent refinement through in silico blood-brain barrier penetration estimation, safety evaluation, and molecular dynamics simulations resulting in two hit compounds that constitute a rational basis for further development of multi-target active compounds against NDDs.     

Computer-Aided Design of Orally Bioavailable Pyrrolidine Carboxamide Inhibitors of Enoyl-Acyl Carrier Protein Reductase of Mycobacterium tuberculosis with Favorable Pharmacokinetic Profiles

We have carried out a computational structure-based design of new potent pyrrolidine carboxamide (PCAMs) inhibitors of enoyl-acyl carrier protein reductase (InhA) of Mycobacterium tuberculosis (MTb). Three-dimensional (3D) models of InhA-PCAMx complexes were prepared by in situ modification of the crystal structure of InhA-PCAM1 (Protein Data Bank (PDB) entry code: 4U0J), the reference compound of a training set of 20 PCAMs with known experimental inhibitory potencies (IC₅₀^exp). First, we built a gas phase quantitative structure-activity relationships (QSAR) model, linearly correlating the computed enthalpy of the InhA-PCAM complex formation and the IC₅₀^exp. Further, taking into account the solvent effect and loss of inhibitor entropy upon enzyme binding led to a QSAR model with a superior linear correlation between computed Gibbs free energies (ΔΔG_com) of InhA-PCAM complex formation and IC₅₀^exp (pIC₅₀^exp = −0.1552·ΔΔG_com + 5.0448, R² = 0.94), which was further validated with a 3D-QSAR pharmacophore model generation (PH4). Structural information from the models guided us in designing of a virtual combinatorial library (VL) of more than 17 million PCAMs. The VL was adsorption, distribution, metabolism and excretion (ADME) focused and reduced down to 1.6 million drug like orally bioavailable analogues and PH4 in silico screened to identify new potent PCAMs with predicted IC₅₀^pre reaching up to 5 nM. Combining molecular modeling and PH4 in silico screening of the VL resulted in the proposed novel potent antituberculotic agent candidates with favorable pharmacokinetic profiles.     

In Silico Drug Repositioning to Target the SARS-CoV-2 Main Protease as Covalent Inhibitors Employing a Combined Structure-Based Virtual Screening Strategy of Pharmacophore Models and Covalent Docking

The epidemic caused by the SARS-CoV-2 coronavirus, which has spread rapidly throughout the world, requires urgent and effective treatments considering that the appearance of viral variants limits the efficacy of vaccines. The main protease of SARS-CoV-2 (M^pro) is a highly conserved cysteine proteinase, fundamental for the replication of the coronavirus and with a specific cleavage mechanism that positions it as an attractive therapeutic target for the proposal of irreversible inhibitors. A structure-based strategy combining 3D pharmacophoric modeling, virtual screening, and covalent docking was employed to identify the interactions required for molecular recognition, as well as the spatial orientation of the electrophilic warhead, of various drugs, to achieve a covalent interaction with Cys145 of M^pro. The virtual screening on the structure-based pharmacophoric map of the SARS-CoV-2 M^pro in complex with an inhibitor N3 (reference compound) provided high efficiency by identifying 53 drugs (FDA and DrugBank databases) with probabilities of covalent binding, including N3 (Michael acceptor) and others with a variety of electrophilic warheads. Adding the energy contributions of affinity for non-covalent and covalent docking, 16 promising drugs were obtained. Our findings suggest that the FDA-approved drugs Vaborbactam, Cimetidine, Ixazomib, Scopolamine, and Bicalutamide, as well as the other investigational peptide-like drugs (DB04234, DB03456, DB07224, DB7252, and CMX-2043) are potential covalent inhibitors of SARS-CoV-2 M^pro.     

Identification of New KRAS G12D Inhibitors through Computer-Aided Drug Discovery Methods

Owing to several mutations, the oncogene Kirsten rat sarcoma 2 viral oncogene homolog (KRAS) is activated in the majority of cancers, and targeting it has been pharmacologically challenging. In this study, using an in silico approach comprised of pharmacophore modeling, molecular docking, and molecular dynamics simulations, potential KRAS G12D inhibitors were investigated. A ligand-based common feature pharmacophore model was generated to identify the framework necessary for effective KRAS inhibition. The chemical features in the selected pharmacophore model comprised two hydrogen bond donors, one hydrogen bond acceptor, two aromatic rings and one hydrophobic feature. This model was used for screening in excess of 214,000 compounds from InterBioScreen (IBS) and ZINC databases. Eighteen compounds from the IBS and ten from the ZINC database mapped onto the pharmacophore model and were subjected to molecular docking. Molecular docking results highlighted a higher affinity of four hit compounds towards KRAS G12D in comparison to the reference inhibitor, BI-2852. Sequential molecular dynamics (MD) simulation studies revealed all four hit compounds them possess higher KRAS G12D binding free energy and demonstrate stable polar interaction with key residues. Further, Principal Component Analysis (PCA) analysis of the hit compounds in complex with KRAS G12D also indicated stability. Overall, the research undertaken provides strong support for further in vitro testing of these newly identified KRAS G12D inhibitors, particularly Hit1 and Hit2.