Targeting HCV polymerase: a structural and dynamic perspective into the mechanism of selective covalent inhibition

Ion-Mediated Protein Stabilization on Nanoscopic Surfaces

The emergence of nanoparticles in biomedical applications has made their interactions with proteins inevitable. Nanoparticles conjugated with proteins and peptide-based constructs form an integral part of nanotherapeutics and have recently shown promise in treating a myriad of diseases. The proper functioning of proteins is critical to achieve their biological functions. However, interface issues result in the denaturation of proteins, and the loss of orientation and steric hindrance can adversely affect the function of the conjugate. Furthermore, surface-induced denaturation also triggers protein aggregation, resulting in amyloid-like species. Understanding the mechanistic underpinnings of protein-nanoparticle interactions and controlling their interfacial characteristics are critical and challenging due to the complex nature of the conjugates. In this milieu, we demonstrate that ionic liquids can be suitable candidates for stabilizing protein-nanoparticle interactions by virtue of their excellent protein-preserving properties. We also probe the previously unexplored mechanism of ion-mediated stabilization of the protein molecules on the nanoparticle surface. The protein-nanoparticle conjugates consist of lysozyme and choline-based ionic liquids characterized by optical and electron microscopy techniques combined with surface-sensitive plasmon-enhanced Raman spectroscopy. Furthermore, atomistic molecular dynamics simulations of the conjugates delineate interfacial interactions of the protein molecules and the modulation by the ions, particularly the conformational changes and the dynamic correlation when the protein and specific ionic liquid molecules are adsorbed on the nanoparticle surface. The combined experimental and computational studies showed the synergistic behavior of the ions of the ionic liquids, specifically the orientation and coverage of the anions aided by the cations to control the surface interactions and hence the overall protein stability. These studies pave the way for using ionic liquids, particularly their biocompatible counterparts in nanoparticle-based complexes, as stabilizing agents for biomedical applications.

Molecular Dynamics Study on Selected Bioactive Phytochemicals as Potential Inhibitors of HIV-1 Subtype C Protease

Acquired immunodeficiency syndrome (AIDS), one of the deadliest global diseases, is caused by the Human Immunodeficiency Virus (HIV). To date, there are no known conventional drugs that can cure HIV/AIDS, and this has prompted continuous scientific efforts in the search for novel and potent anti-HIV therapies. In this study, molecular dynamics simulation (MDS) and computational techniques were employed to investigate the inhibitory potential of bioactive compounds from selected South African indigenous plants against HIV-1 subtype C protease (HIVpro). Of the eight compounds (CMG, MA, UA, CA, BA, UAA, OAA and OA) evaluated, only six (CMG (-9.9 kcal/mol), MA (-9.3 kcal/mol), CA (-9.0 kcal/mol), BA (-8.3 kcal/mol), UAA (-8.5 kcal/mol), and OA (-8.6 kcal/mol)) showed favourable activities against HIVpro and binding landscapes like the reference FDA-approved drugs, Lopinavir (LPV) and Darunavir (DRV), with CMG and MA having the highest binding affinities. Using the structural analysis (root-mean-square deviation (RMSD), fluctuation (RMSF), and radius of gyration (RoG) of the bound complexes with HIVpro after 350 ns, structural evidence was observed, indicating that the six compounds are potential lead candidates for inhibiting HIVpro. This finding was further corroborated by the structural analysis of the enzyme-ligand complexe systems, where structural mechanisms of stability, flexibility, and compactness of the study metabolites were established following binding with HIVpro. Furthermore, the ligand interaction plots revealed that the metabolites interacted hydrophobically with the active site amino residues, with identification of other key residues implicated in HIVpro inhibition for drug design. Overall, this is the first computational report on the anti-HIV-1 activities of CMG and MA, with efforts on their in vitro and in vivo evaluations underway. Judging by the binding affinity, the degree of stability, and compactness of the lead metabolites (CMG, MA, CA, BA, OA, and UAA), they could be concomitantly explored with conventional HIVpro inhibitors in enhancing their therapeutic activities against the HIV-1 serotype.

A molecular dynamics perspective into estrogen receptor inhibition by selective flavonoids as alternative therapeutic options

Zearalenone is an estrogenic mycotoxin which is a common food contaminant and has been implicated in increasing the incidence of carcinogenesis and other reproductive health ailments through the estrogen receptor alpha (ERα) pathway. Competitive ERα blockers such as 4-Hydroxytamoxifen (OHT), are synthetic FDA approved drugs which, albeit being an effective anticancer agent, induces life altering side effects. For this reason, there is an increased interest in the use of naturally occurring medicinal plant products such as flavonoids. This study aimed to identity flavonoid ERα inhibitors and provide insights into the mechanism of inhibition using computational techniques. The Molecular Mechanics/Generalized Born Surface Area calculations revealed that quercetrin, hesperidin, epigallocatechin 3-gallate and kaempferol 7-O-glucoside out of 14 flavonoids had higher binding affinity for ERα than OHT. The structural analysis revealed that the binding of the compounds to the receptor lead to dynamic alterations, which induced conformational shift in the structure and orientation of the receptor resulting in stabilised, compact and low energy systems. The results of this study provide imperative information that supports the use of flavonoids in the inhibition of ERα to prevent or ameliorate the consequential adverse effects associated with zearalenone exposure.Communicated by Ramaswamy H. Sarma.

Deciphering the Mechanism of Binding Selectivity of Chlorofluoroacetamide-Based Covalent Inhibitors toward L858R/T790M Resistance Mutation

Covalent modification of the oncogenic mutant epidermal growth factor receptor (EGFR) by small molecules is an efficient strategy for achieving an enhanced and sustained pharmacological effect in the treatment of non-small-cell lung cancer. NSP-037 (18), an irreversible inhibitor of the L858R/T790M double-mutant EGFR (EGFR^DM) using α-chlorofluoroacetamide (CFA) as a novel warhead, has seven times the inhibition selectivity for EGFR^DM over the wild type (EGFR^WT), as compared to clinically approved osimertinib (7). Here, we employ multiple computational approaches to elucidate the mechanism underlining this improved selectivity, as well as the effect of CFA on the selectivity enhancement of inhibitor 18 over 7. We find that EGFR^DM undergoes significantly larger conformational changes than EGFR^WT upon binding to 18. The conformational stability of the diamine side chain and the CFA motif of 18 in the orthosteric site of EGFR^DM is identified as key for the disparate binding mechanism and inhibitory prowess of 18 with respect to EGFR^WT and EGFR^DM and 18's higher selectivity than 7. The binding free energy of the 18-bound complexes is -6.38 kcal/mol greater than that of the 7-bound complexes, explaining the difference in selectivity of these inhibitors. Further, free energy decomposition analysis indicates that the electrostatic contribution of key residues plays an important role in the 18-bound complexes. QM/MM calculations show that the most favored mechanism for the Cys797 alkylation reaction is the direct displacement mechanism through a CFA-based inhibitor, producing a reaction with the lowest energy barrier and most stable product.

Identification of promising nutraceuticals against filarial immune-modulatory proteins: insights from in silico and ex vivo studies

Lymphatic filariasis is a neglected tropical disease affecting over 863 million people in 47 countries of the world. The anti-filarial drugs, diethylcarbamazine, albendazole, and ivermectin, are effective only at the larval stages and have proven completely ineffective as adulticides. Besides this, a long-term use of these drugs is associated with several side effects including drug toxicity. Nutraceuticals have emerged as better alternatives for long term treatments due to their safety and lesser side effects. In the present work, we have used drug docking analysis and molecular dynamics simulation approaches to explore the effect of anti-inflammatory nutraceuticals against the immune-modulatory proteins of filarial worms. The filarial proteins enolase, ES-62 precursor, serpin, and cystatin, which are highly efficient in host immune modulation were targeted with more than 50 nutraceuticals. In the in silico study nutraceuticals such as naringin, β-carotene, and emodin showed higher binding efficacy and lower dissociation constant as compared to anti-filarial drugs. Molecular dynamics simulation results showed that immune-modulatory proteins formed highly stable complexes with naringin, β-carotene, and emodin over the entire MD simulation run. The nutraceutical emodin formed the most stable system in silico and hence its effect was investigated on adult filarial parasites under ex vivo conditions too. Emodin significantly affected the motility, viability, ROS production, and genomic DNA fragmentation of filarial parasites. Further in vivo and in vitro studies will help in understanding the mechanism of action of emodin at the molecular level and would help in the development of more effective anti-filarial drugs.

Here in drug docking analysis, molecular dynamics simulations and ex vivo approaches were used to demonstrate the anti-filarial effects of nutraceuticals against immune modulatory proteins of lymphatic filarial parasites.

Acquired HIV-1 Protease Conformational Flexibility Associated with Lopinavir Failure May Shape the Outcome of Darunavir Therapy after Antiretroviral Therapy Switch

Understanding the underlying molecular interaction during a therapy switch from lopinavir (LPV) to darunavir (DRV) is essential to achieve long-term virological suppression. We investigated the kinetic and structural characteristics of multidrug-resistant South African HIV-1 subtype C protease (HIV-1 PR) during therapy switch from LPV to DRV using enzyme activity and inhibition assay, fluorescence spectroscopy, and molecular dynamic simulation. The HIV-1 protease variants were from clinical isolates with a combination of drug resistance mutations; MUT-1 (M46I, I54V, V82A, and L10F), MUT-2 (M46I, I54V, L76V, V82A, L10F, and L33F), and MUT-3 (M46I, I54V, L76V, V82A, L90M, and F53L). Enzyme kinetics analysis shows an association between increased relative resistance to LPV and DRV with the progressive decrease in the mutant HIV-1 PR variants' catalytic efficiency. A direct relationship between high-level resistance to LPV and intermediate resistance to DRV with intrinsic changes in the three-dimensional structure of the mutant HIV-1 PR as a function of the multidrug-resistance mutation was observed. In silico analysis attributed these structural adjustments to the multidrug-resistance mutations affecting the LPV and DRV binding landscape. Though DRV showed superiority to LPV, as a lower concentration was needed to inhibit the HIV-1 PR variants, the inherent structural changes resulting from mutations selected during LPV therapy may dynamically shape the DRV treatment outcome after the therapy switch.

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 (CyC_8β) 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 CyC_8β 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 CyC_8β 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 3₁₀-helices in the bound systems. SASA of CyC_8β 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.

Molecular dynamic mechanism(s) of inhibition of bioactive antiviral phytochemical compounds targeting cytochrome P450 3A4 and P-glycoprotein

P-glycoprotein (ABCB1) and cytochrome P450 3A4 (CYP3A4) metabolize almost all known human immunodeficiency virus' protease inhibitor drugs (PIs). Over induction of these proteins' activities has been linked to rapid metabolism of PIs which are then pumped out of the circulatory system, eventually leading to drug-resistance in HIV-positive patients. This study aims to determine, with the use of computational tools, the inhibitory potential of four phytochemical compounds (PCs) (epigallocatechin gallate (EGCG), kaempferol-7-glucoside (K7G), luteolin (LUT) and ellagic acid (EGA)) in inhibiting the activities of these drug-metabolizing proteins. The comparative analysis of the MM/GBSA results revealed that the binding affinity (ΔG_bind) of EGCG and K7G for CYP3A4 and ABCB1 are higher than LUT and EGA and fall between the ΔG_bind of the inhibitors of CYP3A4 and ABCB1 (Ritonavir (strong inhibitor) and Lopinavir (moderate inhibitor)). The structural analysis (RMSD, RMSF, RoG and protein-ligand interaction plots) also confirmed that EGCG and K7G showed similar inhibitory activities with the inhibitors. The study has shown that EGCG and K7G have inhibitory activities against the two proteins and assumes they could decrease intracellular efflux of PIs, consequently increasing the optimal concentration of PIs in the systemic circulation.Communicated by Ramaswamy H. Sarma.

Understanding the Hsp90 N-terminal Dynamics: Structural and Molecular Insights into the Therapeutic Activities of Anticancer Inhibitors Radicicol (RD) and Radicicol Derivative (NVP-YUA922)

Heat shock protein 90 (Hsp90) is a crucial component in carcinogenesis and serves as a molecular chaperone that facilitates protein maturation whilst protecting cells against temperature-induced stress. The function of Hsp90 is highly dependent on adenosine triphosphate (ATP) binding to the N-terminal domain of the protein. Thus, inhibition through displacement of ATP by means of competitive binding with a suitable organic molecule is considered an attractive topic in cancer research. Radicicol (RD) and its derivative, resorcinylic isoxazole amine NVP-AUY922 (NVP), have shown promising pharmacodynamics against Hsp90 activity. To date, the underlying binding mechanism of RD and NVP has not yet been investigated. In this study, we provide a comprehensive understanding of the binding mechanism of RD and NVP, from an atomistic perspective. Density functional theory (DFT) calculations enabled the analyses of the compounds' electronic properties and results obtained proved to be significant in which NVP was predicted to be more favorable with solvation free energy value of -23.3 kcal/mol and highest stability energy of 75.5 kcal/mol for a major atomic delocalization. Molecular dynamic (MD) analysis revealed NVP bound to Hsp90 (NT-NVP) is more stable in comparison to RD (NT-RD). The Hsp90 protein exhibited a greater binding affinity for NT-NVP (-49.4 ± 3.9 kcal/mol) relative to NT-RD (-28.9 ± 4.5 kcal/mol). The key residues influential in this interaction are Gly 97, Asp 93 and Thr 184. These findings provide valuable insights into the Hsp90 dynamics and will serve as a guide for the design of potent novel inhibitors for cancer treatment.

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

Graphical superimposed snapshots of the Thompson novel loop (yellow) of menG protein: apo (A) and bound (B) systems. The loop switches between open and closed conformations; critical for therapeutic activity.

A Mechanistic Study of a Potent and Selective Epidermal Growth Factor Receptor Inhibitor against the L858R/T790M Resistance Mutation

Covalent targeting is a promising strategy for increasing the potency and selectivity of potential drug candidates. This therapeutic approach was recently reported for the epidermal growth factor receptor (EGFR), wherein a covalent binder, 20g [N-(3-{7-[2-methoxy-4-(4-methylpiperazin-1-yl)phenylamino]-3,4-dihydro-3-isopropyl-2,4-dioxopyrimido[4,5-d]pyrimidin-1(2H)-yl}phenyl)acrylamide], demonstrated significant selectivity and inhibitory activity toward the EGFR L858R/T790M double mutant (EGFR^DM) relative to the EGFR wild-type form (EGFR^WT). The enhanced therapeutic potency of 20g against EGFR^DM is 263 times greater than that against EGFR^WT, which necessitates a rational explanation for the underlying selective and inhibitory mechanisms. In this work, we investigate the differential binding modes of 20g with EGFR^WT and EGFR^DM using molecular dynamics simulations coupled with free energy calculations and further identify key residues involved in the selective targeting, binding, and inhibitory mechanisms mediated by 20g. We find that systematic orientational and conformational changes in the α-loop, p-loop, active loop, and αC-helix are responsible for the disparate binding mechanisms and inhibitory prowess of 20g with respect to EGFR^WT and EGFR^DM. The calculated binding free energies show good correlation with the experimental biological activity. The total binding free energy difference between EGFR^WT-20g and EGFR^DM-20g is -11.47 kcal/mol, implying that 20g binds more strongly to EGFR^DM. This enhanced binding affinity of 20g for EGFR^DM is a result of a large increase in the van der Waals and electrostatic interactions with three critical residues (Met790, Gln791, and Met793) that are chiefly responsible for the high-affinity interactions mediated by 20g with EGFR^DM relative to EGFR^WT.