Anti-Aging Potential of the Two Major Flavonoids Occurring in Asian Water Lily Using In Vitro and In Silico Molecular Modeling Assessments

Our previous study investigated the major flavonoids and antioxidant potential of Asian water lily (Nymphaea lotus L., family Nymphaeaceae) stamens and perianth extracts. Quercetin-3-O-rhamnoside (Que-3-Rha) and kaempferol-3-O-galactoside (Kae-3-Gal) were reported as the two most prominent flavonoids found in these extracts. Many flavonoids have been reported on the skin anti-aging effect that are useful for cosmeceutical/phytopharmaceutical application. However, Que-3-Rha and Kae-3-Gal occurring in this medicinal plant have not yet been evaluated for their ability to inhibit skin-aging enzymes. Therefore, this study aimed (1) to assess the enzyme inhibitory activity of Que-3-Rha and Kae-3-Gal, and (2) to conduct molecular modeling of these compounds against critical enzymes involved in skin aging such as collagenase, elastase, and tyrosinase. In vitro enzymatic assays demonstrated that both of the two most prominent flavonoids exhibited moderate to good inhibitory activity toward these enzymes. These experimental findings were supported by molecular docking analysis, which indicated that Que-3-Rha and Kae-3-Gal showed superior binding affinity to the target enzymes compared to the positive controls. Additionally, computational predictions suggested favorable skin permeability and no severe toxicity for both compounds. The results from molecular dynamic (MD) simulation revealed that all the complexes remained stable during the 200 ns MD simulation. Structural analyses and binding free energy calculations also supported the inhibitory potential of these two flavonoids against skin-aging enzymes. In conclusion, this study provides valuable insights into the anti-aging potential of the two major flavonoids occurring in this medicinal plant, paving the way for further development of cosmeceutical/phytopharmaceutical products targeting skin aging.

Efficiency of membrane fusion inhibitors on different hemagglutinin subtypes: insight from a molecular dynamics simulation perspective

The challenge in vaccine development, along with drug resistance issues, has encouraged the search for new anti-influenza drugs targeting different viral proteins. Hemagglutinin (HA) glycoprotein, crucial in the viral replication cycle, has emerged as a promising therapeutic target. CBS1117 and JNJ4796 were reported to exhibit similar potencies against infectious group 1 influenza, which included H1 and H5 HAs; however, their potencies were significantly reduced against group 2 HA. This study aims to explore the molecular binding mechanisms and group specificity of these fusion inhibitors against both group 1 (H5) and group 2 (H3) HA influenza viruses using molecular dynamics simulations. CBS1117 and JNJ4796 exhibit stronger interactions with key residues within the H5 HA binding pocket compared to H3-ligand complexes. Hydrogen bonding and hydrophobic interactions involving residues, such as H381, Q401, T3251 (H5-CBS1117), T3181 (H5-JNJ4796), W212, I452, V482, and V522 predominantly contribute to stabilizing H5-ligand systems. In contrast, these interactions are notably weakened in H3-inhibitor complexes. Predicted protein-ligand binding free energies align with experimental data, indicating CBS1117 and JNJ4796's preference for heterosubtypic group 1 HA binding. Understanding the detailed atomistic mechanisms behind the varying potencies of these inhibitors against the two HA groups can significantly contribute to the development and optimization of effective HA fusion inhibitors. To accomplish this, the knowledge of the transition of HA from its pre- to post-fusion states, the molecular size of ligands, and their potential binding regions, could be carefully considered.Communicated by Ramaswamy H. Sarma.

Naturally occurring quercetin and myricetin as potent inhibitors for human ectonucleotide pyrophosphatase/phosphodiesterase 1

Ecto-nucleotide pyrophosphatases/phosphodiesterases 1 (ENPP1) is a key enzyme in purinergic signaling pathways responsible for cell-to-cell communications and regulation of several fundamental pathophysiological processes. In this study, Kyoto Green, a rapid chemical sensor of pyrophosphate, was employed to screen for effective ENPP1 inhibitors among five representative flavonoids (quercetin, myricetin, morin, kaempferol, and quercetin-3-glucoside), five nucleosides (adenosine, guanosine, inosine, uridine, and cytidine), and five deoxynucleosides (2'- and 3'-deoxyadenosine, 2'-deoxyguanosine, 2'-deoxyinosine, and 2'-deoxyuridine). Conventional colorimetric, fluorescence, and bioluminescence assays revealed that ENPP1 was effectively inhibited by quercetin (Ki ~ 4 nM) and myricetin (Ki ~ 32 nM) when ATP was used as a substrate at pH 7.4. In silico analysis indicated that the presence of a chromone scaffold, particularly one containing a hydroxyl group at the 3' position on the B ring, may promote binding to the active site pocket of ENPP1 and enhance inhibition. This study demonstrated that the naturally derived quercetin and myricetin could effectively inhibit ENPP1 enzymatic activity and may offer health benefits in arthritis management.

Molecular binding of different classes of organophosphates to methyl parathion hydrolase from Ochrobactrum species

Methyl parathion hydrolase (MPH) is an enzyme of the metallo-β-lactamase superfamily, which hydrolyses a wide range of organophosphates (OPs). Recently, MPH has attracted attention as a promising enzymatic bioremediator. The crystal structure of MPH enzyme shows a dimeric form, with each subunit containing a binuclear metal ion center. MPH also demonstrates metal ion-dependent selectivity patterns. The origins of these patterns remain unclear but are linked to open questions about the more general role of metal ions in functional evolution and divergence within enzyme superfamilies. We aimed to investigate and compare the binding of different OP pesticides to MPH with cobalt(II) metal ions. In this study, MPH was modeled from Ochrobactrum sp. with different OP pesticides bound, including methyl paraoxon and dichlorvos and profenofos. The docked structures for each substrate optimized by DFT calculation were selected and subjected to atomistic molecular dynamics simulations for 500 ns. It was found that alpha metal ions did not coordinate with all the pesticides. Rather, the pesticides coordinated with less buried beta metal ions. It was also observed that the coordination of beta metal ions was perturbed to accommodate the pesticides. The binding free energy calculations and structure-based pharmacophore model revealed that all the three substrates could bind well at the active site. However, profenofos exhibit a stronger binding affinity to MPH in comparison to the other two substrates. Therefore, our findings provide molecular insight on the binding of different OP pesticides which could help us design the enzyme for OP pesticides degradation.

Exploring Major Flavonoid Phytochemicals from Nelumbo nucifera Gaertn. as Potential Skin Anti-Aging Agents: In Silico and In Vitro Evaluations

Nelumbo nucifera Gaertn., an aquatic medicinal plant (Nelumbonaceae family), has a history of use in traditional medicine across various regions. Our previous study demonstrated the skin anti-aging potential of its stamen ethanolic extract by effectively inhibiting collagenase and tyrosinase enzymes. While the major constituents of this extract are well documented, there is a lack of research on the individual compounds' abilities to inhibit skin aging enzymes. Therefore, this study aimed to evaluate the anti-aging potential of the primary flavonoids found in N. nucifera using both in silico and in vitro approaches. Our initial step involved molecular docking to identify compounds with the potential to inhibit collagenase, elastase, and tyrosinase. Among the seven flavonoids studied, kaempferol-3-O-robinobioside (Kae-3-Rob) emerged as the most promising candidate, exhibiting the highest docking scores for three skin aging-related enzymes. Subsequent enzyme-based inhibition assays confirmed that Kae-3-Rob displayed robust inhibitory activity against collagenase (58.24 ± 8.27%), elastase (26.29 ± 7.16%), and tyrosinase (69.84 ± 6.07%). Furthermore, we conducted extensive 200-ns molecular dynamics (MD) simulations, revealing the stability of the complexes formed between Kae-3-Rob and each enzyme along the MD simulation time. MM/PBSA-based binding free energy calculations indicated the considerably stronger binding affinity of Kae-3-Rob for collagenase and tyrosinase compared to elastase, which was related to the greater percentage of hydrogen bond occupations. These computational findings were consistent with the relatively high inhibitory activity of Kae-3-Rob against collagenase and tyrosinase observed in our in vitro experiment. In conclusion, the results obtained from this comprehensive study suggest that Kae-3-Rob, a key flavonoid from N. nucifera, holds significant potential as a source of bioactive compounds for anti-aging cosmeceutical and other phytopharmaceutical application.

Structural Dynamics of the Precatalytic State of Human Cytochrome c upon T28C, G34C, and A50C Mutations: A Molecular Dynamics Simulation Perspective

The native structure of cytochrome c (cytc) contains hexacoordinate heme iron with His18 and Met80 residues ligated at the axial sites. Mutations of cytc at Ω-loops have been investigated in modulating the peroxidase activity and, hence, related to the initiation of the apoptotic pathway. Our previous experimental data reported on the peroxidase activity of the cysteine-directed mutants at different parts of the Ω-loop of human cytc (hCytc), that is, T28C, G34C, and A50C. In this work, we performed 1 μs molecular dynamics (MD) simulations to elucidate the detailed structural and dynamic changes upon these mutations, particularly at the proximal Ω-loop. The structures of hCytc were modeled in the hexacoordinated form, which was referred to as the "precatalytic state". The results showed that the structural features of the G34C mutant were more distinctive than those of other mutants. G34C mutation caused local destabilization and flexibility at the proximal Ω-loop (residues 12-28) and an extended distance between this Ω-loop region and heme iron. Besides, analysis of the orientation of the Arg38 side chain of the G34C mutant revealed the Arg38 conformer facing away from the heme iron. The obtained MD results also suggested structural diversity of the precatalytic states for the three hCytc mutants, specifically the effect of G34C mutation on the flexibility of the proximal Ω-loops. Therefore, our MD simulations combined with previous experimental data provide detailed insights into the structural basis of hCytc that could contribute to its pro-apoptotic function.

In Vitro and In Silico Study on the Molecular Encapsulation of α-Tocopherol in a Large-Ring Cyclodextrin

α-tocopherol is the physiologically most active form of vitamin E, with numerous biological activities, such as significant antioxidant activity, anticancer capabilities, and anti-aging properties. However, its low water solubility has limited its potential use in the food, cosmetic, and pharmaceutical industries. One possible strategy for addressing this issue is the use of a supramolecular complex with large-ring cyclodextrins (LR-CDs). In this study, the phase solubility of the CD26/α-tocopherol complex was investigated to assess the possible ratios between host and guest in the solution phase. Next, the host-guest association of the CD26/α-tocopherol complex at different ratios of 1:2, 1:4, 1:6, 2:1, 4:1, and 6:1 was studied by all-atom molecular dynamics (MD) simulations. At 1:2 ratio, two α-tocopherol units interact spontaneously with CD26, forming an inclusion complex, as supported by the experimental data. In the 2:1 ratio, a single α-tocopherol unit was encapsulated by two CD26 molecules. In comparison, increasing the number of α-tocopherol or CD26 molecules above two led to self-aggregation and consequently limited the solubility of α-tocopherol. The computational and experimental results indicate that a 1:2 ratio could be the most suitable stoichiometry to use in the CD26/α-tocopherol complex to improve α-tocopherol solubility and stability in inclusion complex formation.

Molecular insights on the conformational dynamics of a P76C mutant of human cytochrome c and the enhancement on its peroxidase activity

In apoptotic pathway, the interaction of Cytochrome c (Cytc) with cardiolipin in vivo is a key process to induce peroxidase activity of Cytc and trigger the release of Cytc in the inner mitochondria into cytosol. The peroxidase active form of Cytc occurs due to local conformational changes that support the opening of the heme crevice and the loss of an axial ligand between Met80 and heme Fe. Structural adjustments at the Ω-loop segments of Cytc are required for such process. To study the role of the distal Ω-loop segments comprising residues 71-85 in human Cytc (hCytc), we investigated a cysteine mutation at Pro76, one of the highly conserved residues in this loop. The effect of P76C mutant was explored by the combination of experimental characterizations and molecular dynamics (MD) simulations. The peroxidase activity of the P76C mutant was found to be significantly increased by ∼13 folds relative to the wild type. Experimental data on global denaturation, alkaline transition, heme bleaching, and spin-labeling Electron Spin Resonance were in good agreement with the enhancement of peroxidase activity. The MD results of hCytc in the hexacoordinate form suggest the important changes in P76C mutant occurred due to the unfolding at the central Ω-loop (residues 40-57), and the weakening of H-bond between Tyr67 and Met80. Whereas the experimental data implied that the P76C mutant tend to be in equilibrium between the pentacoordinate and hexacoordinate forms, the MD and experimental information are complementary and were used to support the mechanisms of peroxidase active form of hCytc.

Bibenzyl analogue DS-1 inhibits MDM2-mediated p53 degradation and sensitizes apoptosis in lung cancer cells

BACKGROUND

Lung cancer is a leading fatal malignancy due to the high incidence of treatment failure. Dysfunction of the tumor suppressor p53 contributes to cancer initiation, progression, and therapeutic resistance. Targeting MDM2, a negative regulator of p53, has recently attracted interest in cancer drug research as it may restore tumor suppressive function.

PURPOSE

The present study aimed to investigate the effect of 3,4-dihydroxy-5,4'-dimethoxybibenzyl (DS-1) on targeting MDM2 and restoring p53 function in lung cancer cells.

METHODS

The efficacy of DS-1 alone or in combination with cisplatin in lung cancer cells was determined by MTT, nuclear staining, and annexin V/PI assay. The expression of apoptosis-related proteins was determined by western blot analysis. To evaluate the role of DS-1 on the stabilization and degradation of p53, cycloheximide chasing assay and immunoprecipitation were conducted, and the active form of p53 was investigated by immunofluorescent staining assay. To confirm and demonstrate the site interaction between DS-1 and the MDM2 protein, in silico computational analysis was performed.

RESULTS

DS-1 exhibited a cytotoxic effect and sensitized lung cancer cells to cisplatin-induced apoptosis. DS-1 caused a significant increase in the cellular level of p53 protein, while the active form of p53 (phosphorylation at Ser15) was unaltered. DS-1 treatment in combination with cisplatin could enhance activated p-p53 (Ser15) and p53 downstream signaling (Bax, Bcl-2, and Akt), leading to a higher level of apoptosis. Immunoprecipitation analysis revealed that DS-1 decreased the p53-ubiquitin complex, a prerequisite step in p53 proteasomal degradation. Molecular docking simulation further evidenced that DS-1 interacts with MDM2 within the p53-binding domain by carbon-hydrogen bond interaction at Lys27, π-alkyl interactions at Ile37 and Leu30, and van der Waals interactions at Ile75, Val51, Val69, Phe67, Met38, Tyr43, Gly34, and Phe31. Treatment by DS-1 and cisplatin in patient-derivated primary lung cancer cells showed consistent effects by increasing cisplatin sensitivity.

CONCLUSIONS

Our findings provide evidence that DS-1 is an MDM2 inhibitor and its underlying mechanism involves MDM2 binding and p53 induction, which may benefit the development of this compound for lung cancer treatment.

Exploring of paritaprevir and glecaprevir resistance due to A156T mutation of HCV NS3/4A protease: molecular dynamics simulation study

Hepatitis C virus (HCV) NS3/4A serine protease is a promising drug target for the discovery of anti-HCV drugs. However, its amino acid mutations, particularly A156T, commonly lead to rapid emergence of drug resistance. Paritaprevir and glecaprevir, the newly FDA-approved HCV drugs, exhibit distinct resistance profiles against the A156T mutation of HCV NS3/4A serine protease. To illustrate their different molecular resistance mechanisms, molecular dynamics simulations and binding free energy calculations were carried out on the two compounds complexed with both wild-type (WT) and A156T variants of HCV NS3/4A protease. QM/MM-GBSA-based binding free energy calculations revealed that the binding affinities of paritaprevir and glecaprevir towards A156T NS3/4A were significantly reduced by ∼4 kcal/mol with respect to their WT complexes, which were in line with the experimental resistance folds. Moreover, the relatively weak intermolecular interactions with amino acids such as H57, R155, and T156 of NS3 protein, the steric effect and the destabilized protein binding surface, which is caused by the loss of salt bridge between R123 and D168, are the main contributions for the higher fold-loss in potency of glecaprevir due to A156T mutation. An insight into the difference of molecular mechanism of drug resistance against the A156T substitution among the two classes of serine protease inhibitors could be useful for further optimization of new generation HCV NS3/4A inhibitors with enhanced inhibitory potency.Communicated by Ramaswamy H. Sarma.

Structural Basis of Specific Glucoimidazole and Mannoimidazole Binding by Os3BGlu7

β-Glucosidases and β-mannosidases hydrolyze substrates that differ only in the epimer of the nonreducing terminal sugar moiety, but most such enzymes show a strong preference for one activity or the other. Rice Os3BGlu7 and Os7BGlu26 β-glycosidases show a less strong preference, but Os3BGlu7 and Os7BGlu26 prefer glucosides and mannosides, respectively. Previous studies of crystal structures with glucoimidazole (GIm) and mannoimidazole (MIm) complexes and metadynamic simulations suggested that Os7BGlu26 hydrolyzes mannosides via the B_2,5 transition state (TS) conformation preferred for mannosides and glucosides via their preferred ⁴H₃/⁴E TS conformation. However, MIm is weakly bound by both enzymes. In the present study, we found that MIm was not bound in the active site of crystallized Os3BGlu7, but GIm was tightly bound in the -1 subsite in a ⁴H₃/⁴E conformation via hydrogen bonds with the surrounding residues. One-microsecond molecular dynamics simulations showed that GIm was stably bound in the Os3BGlu7 active site and the glycone-binding site with little distortion. In contrast, MIm initialized in the B_2,5 conformation rapidly relaxed to a E₃/⁴H₃ conformation and moved out into a position in the entrance of the active site, where it bound more stably despite making fewer interactions. The lack of MIm binding in the glycone site in protein crystals and simulations implies that the energy required to distort MIm to the B_2,5 conformation for optimal active site residue interactions is sufficient to offset the energy of those interactions in Os3BGlu7. This balance between distortion and binding energy may also provide a rationale for glucosidase versus mannosidase specificity in plant β-glycosidases.

Delivery of Alpha-Mangostin Using Cyclodextrins through a Biological Membrane: Molecular Dynamics Simulation

α-Mangostin (MGS) exhibits various pharmacological activities, including antioxidant, anticancer, antibacterial, and anti-inflammatory properties. However, its low water solubility is the major obstacle for its use in pharmaceutical applications. To increase the water solubility of MGS, complex formation with beta-cyclodextrins (βCDs), particularly with the native βCD and/or its derivative 2,6-dimethyl-β-CD (DMβCD) is a promising technique. Although there have been several reports on the adsorption of βCDs on the lipid bilayer, the release of the MGS/βCDs inclusion complex through the biological membrane remains unclear. In this present study, the release the MGS from the two different βCDs (βCD and DMβCD) across the lipid bilayer was investigated. Firstly, the adsorption of the free MGS, free βCDs, and inclusion complex formation was studied by conventional molecular dynamics simulation. The MGS in complex with those two βCDs was able to spontaneously release free MGS into the inner membrane. However, both MGS and DMβCD molecules potentially permeated into the deeper region of the interior membrane, whereas βCD only adsorbed at the outer membrane surface. The interaction between secondary rim of βCD and the 1-palmitoeyl-2-oleoyl-glycero-3-phosphocholine (POPC) phosphate groups showed the highest number of hydrogen bonds (up to 14) corresponding to the favorable location of βCD on the POPC membrane. Additionally, the findings suggested that electrostatic energy was the main driving force for βCD adsorption on the POPC membrane, while van der Waals interactions played a predominant role in DMβCD adsorption. The release profile of MGS from the βCDs pocket across the lipid bilayer exhibited two energy minima along the reaction coordinate associated with the permeation of the MGS molecule into the deeper region of the POPC membrane.

Why Are Lopinavir and Ritonavir Effective against the Newly Emerged Coronavirus 2019? Atomistic Insights into the Inhibitory Mechanisms

Since the emergence of a novel coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first reported from Wuhan, China, neither a specific vaccine nor an antiviral drug against SARS-CoV-2 has become available. However, a combination of two HIV-1 protease inhibitors, lopinavir and ritonavir, has been found to be effective against SARS-CoV, and both drugs could bind well to the SARS-CoV 3C-like protease (SARS-CoV 3CLpro). In this work, molecular complexation between each inhibitor and SARS-CoV-2 3CLpro was studied using all-atom molecular dynamics simulations, free energy calculations, and pair interaction energy analyses based on MM/PB(GB)SA and FMO-MP2/PCM/6-31G* methods. Both anti-HIV drugs interacted well with the residues at the active site of SARS-CoV-2 3CLpro. Ritonavir showed a somewhat higher number atomic contacts, a somewhat higher binding efficiency, and a somewhat higher number of key binding residues compared to lopinavir, which correspond with the slightly lower water accessibility at the 3CLpro active site. In addition, only ritonavir could interact with the oxyanion hole residues N142 and G143 via the formation of two hydrogen bonds. The interactions in terms of electrostatics, dispersion, and charge transfer played an important role in the drug binding. The obtained results demonstrated how repurposed anti-HIV drugs could be used to combat COVID-19.

Understanding of the drug resistance mechanism of hepatitis C virus NS3/4A to paritaprevir due to D168N/Y mutations: A molecular dynamics simulation perspective

Hepatitis C virus (HCV) NS3/4A protease is an attractive target for the development of antiviral therapy. However, the evolution of antiviral drug resistance is a major problem for treatment of HCV infected patients. Understanding of drug-resistance mechanisms at molecular level is therefore very important for the guidance of further design of antiviral drugs with high efficiency and specificity. Paritaprevir is a potent inhibitor against HCV NS3/4A protease genotype 1a. Unfortunately, this compound is highly susceptible to the substitution at D168 in the protease. In this work, molecular dynamics simulations of paritaprevir complexed with wild-type (WT) and two mutated strains (D168 N and D168Y) were carried out. Due to such mutations, the inhibitor-protein hydrogen bonding between them was weakened and the salt-bridge network among residues R123, R155 and D168 responsible for inhibitor binding was disrupted. Moreover, the per-residue free energy decomposition suggested that the main contributions from key residues such as Q80, V132, K136, G137 and R155 were lost in the D168 N/Y mutations. These lead to a lower binding affinity of paritaprevir for D168 N/Y variants of the HCV NS3/4A protease, consistent with the experimental data. This detailed information could be useful for further design of high potency anti-HCV NS3/4A inhibitors.

The reaction mechanism of Zika virus NS2B/NS3 serine protease inhibition by dipeptidyl aldehyde: a QM/MM study

Zika virus (ZIKV) infection has become a global public health problem, associated with microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there are no commercially available anti-ZIKV drugs. The viral protease NS2B/NS3, which is involved in viral replication and maturation, is a potential drug target. Peptidomimetic aldehyde inhibitors bind covalently to the catalytic S135 of the NS3 protease. Here, we apply hybrid quantum mechanics/molecular mechanics (QM/MM) free-energy simulations at the PDDG-PM3/ff14SB level to investigate the inhibition mechanism of the ZIKV protease by a dipeptidyl aldehyde inhibitor (acyl-KR-aldehyde). The results show that proton transfer from the catalytic S135 to H51 occurs in concert with nucleophilic addition on the aldehyde warhead by S135. The anionic covalent complex between the dipeptidyl aldehyde and the ZIKV protease is analogous to the tetrahedral intermediate for substrate hydrolysis. Spontaneous protonation by H51 forms the hemiacetal. In addition, we use correlated ab initio QM/MM potential energy path calculations at levels up to LCCSD(T)/(aug)-cc-pVTZ to obtain accurate potential energy profiles of the reaction, which also support a concerted mechanism. These results provide detailed insight into the mechanism of ZIKV protease inhibition by a peptidyl aldehyde inhibitor, which will guide in the design of inhibitors.

Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations Support a Concerted Reaction Mechanism for the Zika Virus NS2B/NS3 Serine Protease with Its Substrate

Zika virus (ZIKV) is mainly transmitted to humans by Aedes species mosquitoes and is associated with serious pathological disorders including microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there is no vaccine or anti-ZIKV drug available for preventing or controlling ZIKV infection. An attractive drug target for ZIKV treatment is a two-compartment (NS2B/NS3) serine protease that processes viral polyprotein during infection. Here, conventional molecular dynamics simulations of the ZIKV protease in complex with peptide substrate (TGKRS) sequence at the C-terminus of NS2B show that the substrate is in the active conformation for the cleavage reaction by ZIKV protease. Hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling simulations (PM6/ff14SB) of acylation results reveal that proton transfer from S135 to H51 and nucleophilic attack on the substrate by S135 are concerted. The rate-limiting step involves the formation of a tetrahedral intermediate. In addition, the single-point energy QM/MM calculations, precisely at the level of coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), were performed to correct the potential energy profiles for the first step of the acylation process. The average computed activation barrier at this level of theory is 16.3 kcal mol^-1. Therefore, the computational approaches presented here are helpful for further designing of NS2B/NS3 inhibitors based on transition-state analogues.

Low susceptibility of asunaprevir towards R155K and D168A point mutations in HCV NS3/4A protease: A molecular dynamics simulation

Hepatitis C has become an important health problem that requires expensive treatment and leads to liver tumorigenesis. Hepatitis C virus (HCV), which is the main cause of hepatitis C, has a high mutation rate due to the lack of proofreading activity of the RNA polymerase enzyme. The NS3/4A serine protease is an important target for anti-HCV drug discovery and development because of its crucial role in the cleavage of the polypeptides involved in viral replication. In the present study, all-atom molecular dynamics simulation was performed to elucidate the effect of the single point mutations R155K and D168A in the HCV genotype 1 NS3/4A protease on the structural dynamics, molecular interactions and susceptibility of asunaprevir (ASV), a second-generation NS3/4A protease inhibitor. Principal component analysis indicated that these two mutations converted the direction of motion of residues 123, 155 and 168 in the binding pocket to significantly point outwards from ASV, resulting in a loss of the hydrogen bond network of residues R123···R155···D168. The free energy calculations based on different semiempirical QM/MM-GBSA methods revealed that the binding affinity of ASV with the two mutant forms of the NS3/4A protease was significantly decreased in the order of wild-type < R155K < D168A. This work provided useful structural information regarding the atomistic understanding of acquired drug resistance against ASV caused by the R155K and D168A mutations.

Computational screening of fatty acid synthase inhibitors against thioesterase domain

Thioesterase (TE) domain of fatty acid synthase (FAS) is an attractive therapeutic target for design and development of anticancer drugs. In this present work, we search for the potential FAS inhibitors of TE domain from the ZINC database based on similarity search using three natural compounds as templates, including flavonoids, terpenoids, and phenylpropanoids. Molecular docking was used to predict the interaction energy of each screened ligand compared to the reference compound, which is methyl γ-linolenylfluorophosphonate (MGLFP). Based on this computational technique, rosmarinic acid and its eight analogs were observed as a new series of potential FAS inhibitors, which showed a stronger binding affinity than MGLFP. Afterward, nine docked complexes were studied by molecular dynamics simulations for investigating protein-ligand interactions and binding free energies using MM-PB(GB)SA, MM-3DRISM-KH, and QM/MM-GBSA methods. The binding free energy calculation indicated that the ZINC85948835 (R34) displayed the strongest binding efficiency against the TE domain of FAS. There are eight residues (S2308, I2250, E2251, Y2347, Y2351, F2370, L2427, and E2431) mainly contributed for the R34 binding. Moreover, R34 could directly form hydrogen bonds with S2308, which is one of the catalytic triad of TE domain. Therefore, our finding suggested that R34 could be a potential candidate as a novel FAS-TE inhibitor for further drug design.

Molecular insights into inclusion complexes of mansonone E and H enantiomers with various β-cyclodextrins

The structural dynamics and stability of inclusion complexes of mansonone E (ME) and H (MH) including their stereoisomers with various βCDs (methylated- and hydroxypropylated-βCDs) were investigated by classical molecular dynamics (MD) simulations and binding free energy calculations. The simulation results revealed that mansonones are able to form inclusion complexes with βCDs. The guest molecules are not completely inserted into the host cavity, their preferably positions are nearby the secondary rim with the oxane ring dipping into the hydrophobic inner cavity. The encapsulation process leads to a higher rigidity of the βCDs enhancing the intramolecular hydrogen bond formation ability and decreasing the chance of glucopyranose rotation. According to the MM-PBSA binding free energy calculation, all considered inclusion complexes are stable and the binding energies are mainly caused by van der Waals interactions. Moreover, the free energy calculations showed significant differences in the complexation energies for the stereoisomers, which could enable the separation of the isomers by analytical techniques for further pharmaceutical applications.

Metadynamics supports molecular dynamics simulation-based binding affinities of eucalyptol and beta-cyclodextrin inclusion complexes

The development of various molecular dynamics methods enables the detailed investigation of association processes, like host–guest complexes, including their dynamics and, additionally, the release of the guest compound.

Susceptibility of inhibitors against 3C protease of coxsackievirus A16 and enterovirus A71 causing hand, foot and mouth disease: A molecular dynamics study

Hand foot and mouth disease (HFMD) epidemic has occurred in many countries. Coxsackievirus A16 (CV-A16) and Enterovirus A71 (EV-A71) are the main causes of HFMD. Up to now, there are no anti-HFMD drugs available. Rupintrivir, a broad-spectrum inhibitor, is a drug candidate for HFMD treatment, while other HFMD inhibitors designed from several studies have a relatively low efficiency. Therefore, in this work we aim to study the binding mechanisms of rupintrivir and a peptidic α,β-unsaturated ethyl ester (SG85) against both CV-A16 and EV-A71 3C proteases (3Cpro) using all-atoms molecular dynamics simulation. The obtained results indicate that SG85 shows a stronger binding affinity than rupintrivir against CV-A16. Both inhibitors exhibit a comparable affinity against EV-A71 3Cpro. The molecular information of the binding of the two inhibitors to the proteases will be elucidated. Thus, it is implied that these two compounds may be used as leads for further anti-HFMD drug design and development.

Screening of hepatitis C NS5B polymerase inhibitors containing benzothiadiazine core: a steered molecular dynamics

Hepatic C virus (HCV) is a global health problem, resulting in liver cirrhosis and inflammation that can develop to hepatocellular carcinoma and fatality. The NS5B polymerase of HCV plays an important role in viral RNA replication process, making it an attractive therapeutic target for design and development of anti-HCV drugs. To search new potent compounds against the HCV NS5B polymerase, the molecular docking and the steered molecular dynamics (SMD) simulation techniques were performed. The potential potent inhibitors of the NS5B polymerase were screened out from the ZINC database using structural similarity search and molecular docking technique. Five top-hit compounds (the ZINC compounds 49888724, 49054741, 49777239, 49793673, and 49780355) were then studied by the SMD simulations based on the hypothesis that a high rupture force relates to a high binding efficiency. The results demonstrated that the ZINC compound 49888724 had a greater maximum rupture force, reflecting a good binding strength and inhibitory potency than known inhibitors and the rest four ZINC compounds. Therefore, our finding indicated that the ZINC compound 49888724 is a potential candidate to be a novel NS5B inhibitor for further design. Besides, the van der Waals interaction could be considered as the main contribution for stabilizing the NS5B-ligand complex.

Binding mode and free energy prediction of fisetin/β-cyclodextrin inclusion complexes

In the present study, our aim is to investigate the preferential binding mode and encapsulation of the flavonoid fisetin in the nano-pore of β-cyclodextrin (β-CD) at the molecular level using various theoretical approaches: molecular docking, molecular dynamics (MD) simulations and binding free energy calculations. The molecular docking suggested four possible fisetin orientations in the cavity through its chromone or phenyl ring with two different geometries of fisetin due to the rotatable bond between the two rings. From the multiple MD results, the phenyl ring of fisetin favours its inclusion into the β-CD cavity, whilst less binding or even unbinding preference was observed in the complexes where the larger chromone ring is located in the cavity. All MM- and QM-PBSA/GBSA free energy predictions supported the more stable fisetin/β-CD complex of the bound phenyl ring. Van der Waals interaction is the key force in forming the complexes. In addition, the quantum mechanics calculations with M06-2X/6-31G(d,p) clearly showed that both solvation effect and BSSE correction cannot be neglected for the energy determination of the chosen system.