PCSK9 inhibitors as safer therapeutics for atherosclerotic cardiovascular disease (ASCVD): Pharmacophore design and molecular dynamics analysis.

Cardiovascular disorders are still challenging and are among the deadly diseases. As a major risk factor for atherosclerotic cardiovascular disease, dyslipidemia, and high low-density lipoprotein cholesterol in particular, can be prevented primary and secondary by lipid-lowering medications. Therefore, insights are still needed into designing new drugs with minimal side effects. Proprotein convertase subtilisin/kexin 9 (PCSK9) enzyme catalyses protein-protein interactions with low-density lipoprotein, making it a critical target for designing promising inhibitors compared to statins. Therefore, we screened for potential compounds using a redesigned PCSK9 conformational behaviour to search for a significantly extensive chemical library and investigated the inhibitory mechanisms of the final compounds using integrated computational methods, from ligand essential functional group screening to all-atoms MD simulations and MMGBSA-based binding free energy. The inhibitory mechanisms of the screened compounds compared with the standard inhibitor. K31 and K34 molecules showed stronger interactions for PCSK9, having binding energy (kcal/mol) of -33.39 and -63.51, respectively, against -27.97 of control. The final molecules showed suitable drug-likeness, non-mutagenesis, permeability, and high solubility values. The C-α atoms root mean square deviation and root mean square fluctuation of the bound-PCSK9 complexes showed stable and lower fluctuations compared to apo PCSK9. The findings present a model that unravels the mechanism by which the final molecules proposedly inhibit the PCSK9 function and could further improve the design of novel drugs against cardiovascular diseases.

Characteristic Binding Landscape of Estrogen Receptor-α36 Protein Enhances Promising Cancer Drug Design.

Breast cancer (BC) remains the most common cancer among women worldwide, and estrogen receptor-α expression is a critical diagnostic factor for BC. Estrogen receptor (ER-α36) is a dominant-negative effector of ER-α66-mediated estrogen-responsive gene pathways. ER-α36 is a novel target that mediates the non-genomic estrogen signaling pathway. However, the crystallized structure of ER-α36 remains unavailable for molecular studies. ER-positive and triple-negative BC tumors aggressively resist the FDA-approved drugs; therefore, highly potent structure-based inhibitors with preeminent benefits over toxicity will preferably replace the current BC treatment. Broussoflanol B (BFB), a B. papyrifera bark compound, exhibits potent growth inhibitory activity in ER-negative BC cells by inducing cell cycle arrest. For the first time, we unravel the comparative dynamic events of the enzymes' structures and the binding mechanisms of BFB when bound to the ER-α36 and ER-α66 ligand-binding domain using an all-atom molecular dynamics simulations approach and MM/PBSA-binding-free energy calculations. The dynamic findings have revealed that ER-α36 and ER-α66 LBD undergo timescale "coiling", opening and closing conformations favoring the high-affinity BFB-bound ER-α36 (ΔG = -52.57 kcal/mol) compared to the BFB-bound ER-α66 (ΔG = -42.41 kcal/mol). Moreover, the unbound (1.260 Å) and bound ER-α36 (1.182 Å) exhibit the highest flexibilities and atomistic motions relative to the ER-α66 systems. The RMSF (Å) of the unbound ER-α36 and ER-α66 exhibit lesser stabilities than the BFB-bound systems, resulting in higher structural flexibilities and atomistic motions than the bound variants. These findings present a model that describes the mechanisms by which the BFB compound induces downregulation-accompanied cell cycle arrest at the Gap0 and Gap1 phases.

A novel small molecule inhibitor of p38⍺ MAP kinase augments cardiomyocyte cell cycle entry in response to direct cell cycle stimulation.

BACKGROUND AND PURPOSE: Myocardial infarction (MI) is the leading cause of mortality globally due in part to the limited ability of cardiomyocytes (CMs) to regenerate. Recently, we demonstrated that overexpression of four-cell cycle factors, CDK1, CDK4, cyclin B1 and cyclin D1 (4F), induced cell division in ~20% of the post-mitotic CMs overexpressed 4F. The current study aims to identify a small molecule that augments 4F-induced CM cycle induction. EXPERIMENTAL APPROACH, KEY RESULTS: Screening of small molecules with a potential to augment 4F-induced cell-cycle induction in 60-day-old mature human induced pluripotent cardiomyocytes (hiPS-CMs) revealed N-(4,6-Dimethylpyridin-2-yl)-4-(pyridine-4-yl)piperazine-1-carbothioamide (NDPPC), which activates cell cycle progression in 4F-transduced hiPS-CMs. Autodock tool and Autodock vina computational methods showed that NDPPC has a potential interaction with the binding site at the human p38⍺ mitogen-activated protein kinase (p38⍺ MAP kinase), a critical negative regulator of the mammalian cell cycle. A p38 MAP kinase activity assay showed that NDPPC inhibits p38⍺ with 5-10 times lower IC50 compared to the other P38 isoforms in a dose-dependent manner. Overexpression of p38⍺ MAP kinase in CMs inhibited 4F cell cycle induction, and treatment with NDPPC reversed the cell cycle inhibitory effect. CONCLUSION AND IMPLICATIONS: NDPPC is a novel inhibitor for p38 MAP kinase and is a promising drug to augment CM cell cycle response to the 4F. NDPPC could become an adjunct treatment with other cell cycle activators for heart failure treatment.

Probing into the Flap-dimer Dynamics of the Mycobacterium tuberculosis Kasa Enzyme Binding Landscape Provides the Underlying Inhibitory Mechanisms of JSF-3285 and 5G.

BACKGROUND: ß-ketoacyl-ACP synthase I (KasA I) enzyme is crucial in mycolic acid synthesis via catalytic condensation reactions, hence implicated in M. tuberculosis's virulence and drug resistance. Presently, there is no known potent KasA inhibitor; thiolactomycin lacks potency. Recently reported indazole compounds JSF-3285/tr1DG167 and 5G/tr2DG167 inhibit the KasA through binding to the substrate cavity. However, the molecular mechanism is still unclear, and the unknown resistance mechanisms raise concerns about JSF-3285's novelty. METHODS: This study is the first to report the flap dimer opening and closing of the KasA pocket using combined metrics to define the symmetry impact of the flap-dimer motions and investigate the underlying inhibitory mechanism of tr1DG167 andtr2DG167 using all-atom MD simulation. RESULTS: The distance/d1 between the flap (PRO147) and dimer (LEU205) residues; TriC-α angle (θ1: PRO147-VAL83-LEU205 & θ2: PRO147-GLU199-LEU205); and the dihedral angle (Φ) were applied to investigate the flap "twisting" and dimer shift closing due to concerted motion by adjacent glycine-rich and glutamic acid-rich loops around the active site during the binding pocket's opening. The full flap-dimer of the unbound opens at 230 ns (d1 = 21.51 Å), corresponding to the largest TriC-α angle θ1 44.5° as θ2 is unreliable to describe the flap-dimer motion. The overall averages θ1 and θ2 for the bounds were ~23.13° and ~23.31°, respectively. Thus, the degree of KasA flap dimer opening is best investigated by distance and θ1. BFE (Kcal/mol) of -44.05 (tr1DG167) showed a higher affinity for the pocket than tr2DG167-KasA (-32.16). Both tr1DG167 and tr2DG167 formed hydrophobic interactions with LEU116, GLY117, ALA119, and tr1DG167 formed strong H-bonds with GLU199. The average RMSD of 2.80 Å (Apo) and RoG of 20.97 Å showed that KasA is less stable and less tightly packed without the inhibitors. CONCLUSION: These findings provide a background for a new structure-based design of novel KasA inhibitors.

Nature-Inspired O-Benzyl Oxime-Based Derivatives as New Dual-Acting Agents Targeting Aldose Reductase and Oxidative Stress.

Aldose reductase (ALR2) is the enzyme in charge of developing cellular toxicity caused by diabetic hyperglycemia, which in turn leads to the generation of reactive oxygen species triggering oxidative stress. Therefore, inhibiting ALR2 while pursuing a concomitant anti-oxidant activity through dual-acting agents is now recognized as the gold standard treatment for preventing or at least delaying the progression of diabetic complications. Herein we describe a novel series of (E)-benzaldehyde O-benzyl oximes 6a-e, 7a-e, 8a-e, and 9-11 as ALR2 inhibitors endowed with anti-oxidant properties. Inspired by the natural products, the synthesized derivatives are characterized by a different polyhydroxy substitution pattern on their benzaldehyde fragment, which proved crucial for both the enzyme inhibitory activity and the anti-oxidant capacity. Derivatives (E)-2,3,4-trihydroxybenzaldehyde O-(3-methoxybenzyl) oxime (7b) and (E)-2,3,4-trihydroxybenzaldehyde O-(4-methoxybenzyl) oxime (8b) turned out to be the most effective dual-acting products, proving to combine the best ALR2 inhibitory properties with significant anti-oxidant efficacy.

Simulation Models for Prediction of Bioavailability of Medicinal Drugs-the Interface Between Experiment and Computation.

The oral drug bioavailability (BA) problems have remained inevitable over the years, impairing drug efficacy and indirectly leading to eventual human morbidity and mortality. However, some conventional lab-based methods improve drug absorption leading to enhanced BA, and the recent experimental techniques are up-and-coming. Nevertheless, some have inherent drawbacks in improving the efficacy of poorly insoluble and low impermeable drugs. Drug BA and strategies to overcome these challenges were briefly highlighted. This review has significantly unravelled the different computational models for studying and predicting drug bioavailability. Several computational approaches provide mechanistic insights into the oral drug delivery system simulation of descriptors like solubility, permeability, transport protein-ligand interactions, and molecular structures. The in silico techniques have long been known still are just being applied to unravel drug bioavailability issues. Many publications have reported novel applications of the computational models towards achieving improved drug BA, including predicting gastrointestinal tract (GIT) drug absorption properties and passive intestinal membrane permeability, thus maximizing time and resources. Also, the classical molecular simulation models for free solvation energies of soluble-related processes such as solubilization, dissolutions, supersaturation, and precipitation have been used in virtual screening studies. A few of the tools are GastroPlusTM that supports biowaiver for drugs, mainly BCS class III and predicts drug compounds' absorption and pharmacokinetic process; SimCyp® simulator for mechanistic modelling and simulation of drug formulation processes; pharmacodynamics analysis on non-linear mixed-effects modelling; and mathematical models, predicting absorption potential/maximum absorption dose. This review provides in silico-experiment annexation in the drug bioavailability enhancement, possible insights that lead to critical opinion on the applications and reliability of the various in silico models as a growing tool for drug development and discovery, thus accelerating drug development processes.

An Allosteric Inhibitory Potential of Triterpenes from Combretum racemosum on the Structural and Functional Dynamics of Plasmodium falciparum Lactate Dehydrogenase Binding Landscape.

Multidrug resistance is a significant drawback in malaria treatment, and mutations in the active sites of the many critical antimalarial drug targets have remained challenging. Therefore, this has necessitated the global search for new drugs with new mechanisms of action. Plasmodium falciparum lactate dehydrogenase (pfLHD), a glycolytic enzyme, has emerged as a potential target for developing new drugs due to the parasite reliance on glycolysis for energy. Strong substrate-binding is required in pfLDH enzymatic catalysis; however, there is a lack of information on small molecules' inhibitory mechanism bound to the substrate-binding pocket. Therefore, this study investigated a potential allosteric inhibition of pfLDH by targeting the substrate-binding site. The structural and functional behaviour of madecassic acid (MA), the most promising among the six triterpenes bound to pfLDH, were unravelled using molecular dynamic simulations at 300 ns to gain insights into its mechanism of binding and inhibition and chloroquine as a standard drug. The docking studies identified that the substrate site has the preferred position for the compounds even in the absence of a co-factor. The bound ligands showed comparably higher binding affinity at the substrate site than at the co-factor site. Mechanistically, a characteristic loop implicated in the enzyme catalytic activity was identified at the substrate site. This loop accommodates key interacting residues (LYS174, MET175, LEU177 and LYS179) pivotal in the MA binding and inhibitory action. The MA-bound pfLHD average RMSD (1.60 Å) relative to chloroquine-bound pfLHD RMSD (2.00 Å) showed higher stability for the substrate pocket, explaining the higher binding affinity (-33.40 kcal/mol) observed in the energy calculations, indicating that MA exhibited profound inhibitory activity. The significant pfLDH loop conformational changes and the allostery substrate-binding landscape suggested inhibiting the enzyme function, which provides an avenue for designing antimalarial compounds in the future studies of pfLDH protein as a target.

Weak spots inhibition in the Mycobacterium tuberculosis antigen 85C target for antitubercular drug design through selective irreversible covalent inhibitor-SER124.

Mycobacterium tuberculosis (Mtb) encoded secreted antigen 85 enzymes (Ag85A/Ag85B/Ag85C) play that critical roles in the virulence, survival and drug-resistant TB of the pathogen. Ag85 proteins are potential antitubercular drug targets because they are essential in the catalytic synthesis of trehalose moieties and mycolic acid attachment to the Mtb cell wall. Recently, experimental protocols led to the discovery of a selective covalent Ag85 inhibitor, ß-isomer monocyclic enolphosphorus Cycliphostin (CyC8ß) compound, which targets the Ag85 serine 124 to exhibit a promising therapeutic activity. For the first time, our study unravelled the structural features among Mtb Ag85C homologs and motions and dynamics of Ag85C when the CyC8ß bound covalently and in open model conformations to the protein using bioinformatics tools and integrated Molecular dynamics simulations. Comparative Ag85C sequence analysis revealed conserved regions; 70% active site, 90% Adeniyi loop L1 and 50% loop L2, which acts as a switch between open and closed conformations. The average C-α atoms RMSD (2.05 Å) and RMSF (0.9 Å) revealed instability and high induced flexibility in the CyC8ß covalent-bound compared to the apo and open model systems, which displayed more stability and lower fluctuations. DSSP showed structural transitions of α-helices to bend and loops to 310-helices in the bound systems. SASA of CyC8ß covalent bound showed active site hydrophobic residues exposure to huge solvent. Therefore, these findings present the potential opportunity hotspots in Ag85C protein that would aid the structure-based design of novel chemical entities capable of resulting in potent antitubercular drugs.Communicated by Ramaswamy H. Sarma.

Functional Analysis of Single Nucleotide Polymorphism in ZUFSP Protein and Implication in Pathogenesis.

Researches have revealed that functional non-synonymous Single Nucleotide Polymorphism (nsSNPs) present in the Zinc-finger with UFM1-Specific Peptidase domain protein (ZUFSP) may be involved in genetic instability and carcinogenesis. For the first time, we employed in-silico approach using predictive tools to identify and validate potential nsSNPs that could be pathogenic. Our result revealed that 8 nsSNPs (rs 112738382, rs 140094037, rs 201652589, rs 201847265, rs 202076827, rs 373634906, rs 375114528, rs 772591104) are pathogenic after being subjected to rigorous filtering process. The structural impact of the nsSNPs on ZUFSP structure indicated that the nsSNPs affect the stability of the protein by lowering ZUFSP protein stability. Furthermore, conservation analysis showed that rs 201652589, rs 140094037, rs 201847265, and rs 772591104 were highly conserved. Interestingly, the protein-protein affinity between ZUFSP and Ubiquitin was altered rs 201652589, rs 140094037, rs 201847265, and rs 772591104 had a binding affinity of - 0.46, - 0.83, - 1.62, and - 1.12 kcal/mol respectively. Our study has been able to identify potential nsSNPs that could be used as genetic biomarkers for some diseases arising as a result of aberration in the ZUFSP structure, however, being a predictive study, the identified nsSNPs need to be experimentally investigated.

Leveraging on Active Site Similarities; Identification of Potential Inhibitors of Zinc-Finger and UFSP domain Protein (ZUFSP).

BACKGROUND: ZUFSP (Zinc-finger and UFSP domain protein) is a novel representative member of the recently characterized seventh class of deubiquitinating enzymes (DUBs). Due to the roles DUBs play in genetic instability, they have become a major drug target in cancer and neurodegenerative diseases. ZUFSP, being a DUB enzyme has also been implicated in genetic stability. However, no lead compound has been developed to target ZUFSP. OBJECTIVE/METHODS: Therefore, in this study, we used a combined drug repurposing, virtual screening and per-Residue Energy Decomposition (PRED) to identify ZUFSP inhibitors with therapeutic potential. 3-bromo-6-{[4-hydroxy-1-3(3-phenylbutanoyl)piperidin-4-yl]methyl}-4H,5H,6H,7H-thieno[2,3- C]pyridine-7-one (BHPTP) which is an inhibitor of USP7 was repurposed to target ZUFSP. The rationale behind this is based on the similarity of the active between USP7 and ZUFSP. RESULTS: PRED of the binding between BHPTP and ZUFSP revealed Cys223, Arg408, Met410, Asn460, and Tyr465 as the crucial residues responsible for this interaction. The pharmacophoric moieties of BHPTP responsible for this binding along with other physiochemical properties were used as a filter to retrieve potential ligands. 799 compounds were retrieved, ZINC083241427, ZINC063648749, and ZINC063648753 were selected due to the binding energy they exhibited. Cheminformatics analysis revealed that the compounds possess high membrane permeability, however, BHPTP had a low membrane permeability. Furthermore, the compounds are drug like, having obeyed Lipinski's rule of five. CONCLUSION: Taken together, findings from this study put ZINC083241427, ZINC063648749, and ZINC063648753 as potential ZUFSP inhibitor, however, more experimental validation is required to unravel the mechanism of actions of these compounds.

Delving into the Characteristic Features of "Menace" Mycobacterium tuberculosis Homologs: A Structural Dynamics and Proteomics Perspectives.

The global increase in the morbidity/mortality rate of Mycobacterial infections, predominantly renascent tuberculosis, leprosy, and Buruli ulcers have become worrisome over the years. More challenging is the incidence of resistance mediated by mutant Mycobacterium strains against front-line antitubercular drugs. Homologous to all Mycobacteria species is the GlcNAc-6-phosphate deacetylase (NagA) which catalyzes essential amino sugars synthesis required for cell wall architecture, hence, metamorphosing into an important pharmacological target for curtailing virulence and drug-resistance. This study used integrated bioinformatics methods, MD simulations, and DynaMut and PolyPhen2 to; explore unique features, monitor dynamics, and analyze the functional impact of non-synonymous single-nucleotide polymorphisms of the six NagA of most ruinous Mycobacterium species; tuberculosis (Mtb), smegmatis (MS), marinum (MM), ulcerans, africanum, and microti respectively. This approach is essential for multi-targeting and could result in the identification of potential polypharmacological antitubercular compounds. Comparative sequential analyses revealed ≤ 50% of the overall structure, including the catalytic Asp267 and reactive Cys131, remained conserved. Interestingly, MS-NagA and MM-NagA possess unique hydrophobic isoleucine (Ile) residues at their active sites in contrast to leucine (Leu) found in other variants. More so, unique to the active sites of the NagA is a 'subunit loop' that covers the active site; probably crucial in binding (entry and exit) mechanisms of targeted NagA inhibitors. Relatively, nsSNP mutations exerted a destabilizing effect on the native NagA conformation. Structural and dynamical insights provided, basically pin-pointed the "Achilles' heel" explorable for the rational drug design of target-specific 'NagA' inhibitors potent against a wide range of mycobacterial diseases.

Thompson loop: opportunities for antitubercular drug design by targeting the weak spot in demethylmenaquinone methyltransferase protein.

Drug-resistant Tuberculosis (TB) has remained the top global health challenge, with a yearly estimation of 10 million infections and 1.5 million deaths in humans. Demethylmenaquinone methyltransferase (menG) catalyzes demethylmenaquinone conversion to menaquinone (MK) that is implicated in the TB pathogenesis, hence, it has become a major drug target. DG70 is a biphenyl amide compound known to be a high binding affinity inhibitor of menG. This study investigated the structural and dynamic impacts of DG70 upon binding to menG using atom-based dynamic simulation. Our findings revealed that the modeled structure of menG possesses some Rossman-like methyltransferase characteristic features including two GXG motifs, an omega-like loop (residues 210-220) called the Thompson loop, nine α-helices, five ß-strands, etc. Furthermore, atom-based dynamic simulations revealed that the Thompson loop is critical in the therapeutic activity of DG70. The loop assumed an open conformation in the unliganded-menG structure. However, in the DG70-menG, it assumed a tightly closed conformation. This explains the high binding affinity (-32.48 kcal mol-1) observed in the energy calculations. Interestingly, these findings are further collaborated by the conformational perturbation in the menG protein. Conclusively, insights from this study, highlight the structural "Achilles heel" in menG protein which can be further leveraged by inhibitors tailored to specifically target them.

Exploring the effect of ritonavir and TMC-310911 on SARS-CoV-2 and SARS-CoV main proteases: potential from a molecular perspective.

AIM: As coronavirus (CoV) disease 2019-associated pneumonia spreads globally, there has been an urgent need to combat the spread and develop vaccines. MATERIALS & METHODS: We used an integrated computational algorithm to explore the binding mechanism of TMC-310911/ritonavir (RVT) with SARS-CoV-2 and SARS-CoV main proteases. RESULTS: RVT and TMC-310911 had favorable interactions with the proteases, and these high interactions are facilitated by some significant residues such as Asn133, Gly195 and Gln192. Our study further implicated two important rings in the structure of RVT as a possible chemical culprit in its therapeutic activity. CONCLUSION: Although there are conflicting clinical results on the therapeutic potency of RVT in the treatment of coronavirus disease 2019, our findings provided molecular insight into the binding mechanism of TMC-310911 and RVT with SARS-CoV-2 and SARS-CoV main proteases.