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

Unveiling the synergy: a combined experimental and theoretical study of β-cyclodextrin with melatonin

Melatonin (MT) is a vital hormone controlling biorhythms, and optimizing its release in the human body is crucial. To address MT's unfavorable pharmacokinetics, we explored the inclusion complexes of MT with β-cyclodextrin (β-CD). Nano spray drying was applied to efficiently synthesize these complexes in three molar ratios (MT : β-CD = 1 : 1, 2 : 1, and 1 : 2), reducing reagent use and expediting inclusion. The complex powders were characterized through thermal analyses (TGA and DSC), Fourier transform infrared spectroscopy (FTIR), and in vitro MT release measurements via high-performance liquid chromatography (HPLC). In parallel, computational studies were conducted, examining the stability of MT : β-CD complexes by means of unbiased semi-empirical conformational searches refined by DFT, which produced a distribution of MT : β-CD binding enthalpies. Computational findings highlighted that these complexes are stabilized by specific hydrogen bonds and non-specific dispersive forces, with stronger binding in the 1 : 1 complex, which was corroborated by in vitro release data. Furthermore, the alignment between simulated and experimental FTIR spectra demonstrated the quality of both the structural model and computational methodology, which was crucial to enhance our comprehension of optimizing MT's release for therapeutic applications.

Mechanochemical Approach to Obtaining a Multicomponent Fisetin Delivery System Improving Its Solubility and Biological Activity

In this study, binary amorphous solid dispersions (ASDs, fisetin-Eudragit®) and ternary amorphous solid inclusions (ASIs, fisetin-Eudragit®-HP-β-cyclodextrin) of fisetin (FIS) were prepared by the mechanochemical method without solvent. The amorphous nature of FIS in ASDs and ASIs was confirmed using XRPD (X-ray powder diffraction). DSC (Differential scanning calorimetry) confirmed full miscibility of multicomponent delivery systems. FT-IR (Fourier-transform infrared analysis) confirmed interactions that stabilize FIS's amorphous state and identified the functional groups involved. The study culminated in evaluating the impact of amorphization on water solubility and conducting in vitro antioxidant assays: 2,2-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid)-ABTS, 2,2-diphenyl-1-picrylhydrazyl-DPPH, Cupric Reducing Antioxidant Capacity-CUPRAC, and Ferric Reducing Antioxidant Power-FRAP and in vitro neuroprotective assays: inhibition of acetylcholinesterase-AChE and butyrylcholinesterase-BChE. In addition, molecular docking allowed for the determination of possible bonds and interactions between FIS and the mentioned above enzymes. The best preparation turned out to be ASI_30_EPO (ASD fisetin-Eudragit® containing 30% FIS in combination with HP-β-cyclodextrin), which showed an improvement in apparent solubility (126.5 ± 0.1 µg∙mL-1) and antioxidant properties (ABTS: IC50 = 10.25 µg∙mL-1, DPPH: IC50 = 27.69 µg∙mL-1, CUPRAC: IC0.5 = 9.52 µg∙mL-1, FRAP: IC0.5 = 8.56 µg∙mL-1) and neuroprotective properties (inhibition AChE: 39.91%, and BChE: 42.62%).

Interactions between β-cyclodextrin as a carrier for anti-cancer drug delivery: a molecular dynamics simulation study

A series of molecular dynamics simulations were performed on 5-fluorouracil (5-Fu), Alendronate (Ald), and Temozolomide (TMZ) anticancer drugs in the presence and absence of β-cyclodextrin (βCD) as a carrier. Thermodynamic investigations showed that the van der Waals interaction energy was dominant in loading all drugs inside the βCD cavity. The sum of the interaction energies illustrated that the highest affinity was related to Ald (-136.5 kJ/mol), which in turn was due to the presence of bulky and charged atoms of phosphorus and oxygen, although TMZ (-115.92 kJ/mol) showed a very high affinity as well. At the same time, the hydrogen bond analysis also represented that Ald had the most hydrogen bond (1.97) with the highest half-life (3.13 ps) with βCD. Investigation of the root mean fluctuation (RMSF) indicated that all the drugs had a relatively rigid structure and maintain this rigidity during loading in the βCD cavity, and in the meantime, Ald was slightly more flexible than 5-Fu and TMZ. The area of ​the primary hydroxyl rim decreased in all drug-containing systems, which in turn was caused by the attractive interaction of drugs with oxygens in the primary hydroxyl rim. Especially for those drugs that were able to penetrate to the end of the primary hydroxyl rim of the βCD, that means TMZ and 5-Fu. Meanwhile, due to the lack of Ald penetration to the end of the primary hydroxyl rim, the area change in the Ald-containing system was less than in the two others.Communicated by Ramaswamy H. Sarma.

Current Status of Quantum Chemical Studies of Cyclodextrin Host-Guest Complexes

This article aims to review the application of various quantum chemical methods (semi-empirical, density functional theory (DFT), second order Møller-Plesset perturbation theory (MP2)) in the studies of cyclodextrin host-guest complexes. The details of applied approaches such as functionals, basis sets, dispersion corrections or solvent treatment methods are analyzed, pointing to the best possible options for such theoretical studies. Apart from reviewing the ways that the computations are usually performed, the reasons for such studies are presented and discussed. The successful applications of theoretical calculations are not limited to the determination of stable conformations but also include the prediction of thermodynamic properties as well as UV-Vis, IR, and NMR spectra. It has been shown that quantum chemical calculations, when applied to the studies of CD complexes, can provide results unobtainable by any other methods, both experimental and computational.

Complexation of Lignin Dimers with β-Cyclodextrin and Binding Stability Analysis by ESI-MS, Isothermal Titration Calorimetry, and Molecular Dynamics Simulations

Lignin derived from lignocellulosic biomass is the largest source of renewable bioaromatics present on earth and requires environmentally sustainable separation strategies to selectively obtain high-value degradation products. Applications of supramolecular interactions have the potential to isolate lignin compounds from biomass degradation fractions by the formation of variable inclusion complexes with cyclodextrins (CDs). CDs are commonly used as selective adsorbents for many applications and can capture guest molecules in their internal hydrophobic cavity. The strength of supramolecular interactions between CDs and lignin model compounds that represent potential lignocellulosic biomass degradation products can be characterized by assessing the thermodynamics of binding stability. Consequently, the inclusion interactions of β-CD and lignin model compounds G-(β-O-4')-G, G-(β-O-4')-truncG (guaiacylglycerol-β-guaiacyl ether), and G-(β-β')-G (pinoresinol) were investigated empirically by electrospray ionization mass spectrometry and isothermal titration calorimetry, complemented by molecular dynamics (MD) simulations. Empirical results indicate that there are substantial differences in binding stability dependent on the linkage type. The lignin model β-β' dimer showed more potential bound states including 1:1, 2:1, and 1:2 (guest:host) complexation and, based on binding stability determinations, was consistently the most energetically favorable guest. Empirical results are supported by MD simulations that reveal that the capture of G-(β-β')-G by β-CD is promising with a 66% probability of being bound for G-(β-O-4')-truncG compared to 88% for G-(β-β')-G (unbiased distance trajectory and explicit counting of bound states). These outcomes indicate CDs as a promising material to assist in separations of lignin oligomers from heterogeneous mixtures for the development of environmentally sustainable isolations of lignin compounds from biomass fractions.

NS3/4A helicase inhibitory alkaloids from Aptenia cordifolia as HCV target

Chemical investigation of Aptenia cordifolia roots extract, using chromatographic and spectroscopic techniques, resulted in isolation and identification of eight known compounds. The basic ethyl acetate fraction (alkaloidal fraction) afforded O-methylsceletenone, epinine, 4-methoxy phenethylamine, and N-methyl tyramine while, the acidic ethyl acetate fraction (non-alkaloidal fraction) afforded only cis-N-coumaroyl tyramine. Moreover, the petroleum ether fraction afforded capric acid, tricosanol, and a mixture of β-sitosterol & stigma sterol. Upon screening of anti HCV activity of these three fractions, only the basic ethyl acetate fraction had high activity against HCV with an IC50 value equal to 2.4 μg mL-1 which provoked us to carry out structure based in silico virtual screening on the drug targets of HCV of isolated alkaloidal compounds as well as the previously dereplicated alkaloids through metabolomics from the antiviral active fraction. The tortuosamine compound exhibited the strongest binding to the active site of NS3/4A helicase with a binding affinity (-7.1 kcal mol-1) which is very close to the native ligand (-7.7 kcal mol-1).

Illustration of a computational pipeline for evaluating cyclodextrin host-guest complex formation through conformational capture of bullvalene

Cyclodextrins have a diverse range of applications, including as supramolecular hosts, as enzyme active-site analogs, in improving drug solubility and delivery, and in molecular selection. We have investigated their ability to form stable complexes with bullvalenes, unusual organic cage molecules that spontaneously interconvert between numerous degenerate isomers. The shape-shifting nature of substituted bullvalenes raises the potential for dynamic adaptive binding to biological targets. We tested whether β- and γ-cyclodextrins can capture particular bullvalene isomers and whether the preferred binding mode(s) differ between isomers. We first applied our computational host-guest interaction potential energy profiling to determine the best binding mode(s) of unsubstituted bullvalene and each isomer of methylenehydroxybullvalene to β- and γ-cyclodextrin. Subsequent molecular dynamics simulations of the predicted host-guest complexes showed that while unsubstituted bullvalene has a single, albeit ill-defined, binding mode with either cyclodextrin, each isomer of methylenehydroxybullvalene has two possible modes of binding to β-cyclodextrin but only a single, nebulous mode of binding to γ-cyclodextrin. Experimental determination of the binding free energy of each methylenehydroxybullvalene-cyclodextrin complex showed that methylenehydroxybullvalene is more likely to bind to β-cyclodextrin than to γ-cyclodextrin, despite its smaller cavity. Together, our results suggest that β-cyclodextrin, but not γ-cyclodextrin, shows promise for conformational capture of mono-substituted bullvalenes. More broadly, our computational pipeline should prove useful for rapid characterization of cyclodextrin host-guest complexes, avoiding the need for costly synthesis of guest molecules that are unlikely to bind stably, as well as providing detailed atomic-level insight into the nature of complexation.

Effect of Water Microsolvation on the Excited-State Proton Transfer of 3-Hydroxyflavone Enclosed in γ-Cyclodextrin

The effect of microsolvation on excited-state proton transfer (ESPT) reaction of 3-hydroxyflavone (3HF) and its inclusion complex with γ-cyclodextrin (γ-CD) was studied using computational approaches. From molecular dynamics simulations, two possible inclusion complexes formed by the chromone ring (C-ring, Form I) and the phenyl ring (P-ring, Form II) of 3HF insertion to γ-CD were observed. Form II is likely more stable because of lower fluctuation of 3HF inside the hydrophobic cavity and lower water accessibility to the encapsulated 3HF. Next, the conformation analysis of these models in the ground (S₀) and the first excited (S₁) states was carried out by density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations, respectively, to reveal the photophysical properties of 3HF influenced by the γ-CD. The results show that the intermolecular hydrogen bonding (interHB) between 3HF and γ-CD, and intramolecular hydrogen bonding (intraHB) within 3HF are strengthened in the S₁ state confirmed by the shorter interHB and intraHB distances and the red-shift of O–H vibrational modes involving in the ESPT process. The simulated absorption and emission spectra are in good agreement with the experimental data. Significantly, in the S₁ state, the keto form of 3HF is stabilized by γ-CD, explaining the increased quantum yield of keto emission of 3HF when complexing with γ-CD in the experiment. In the other word, ESPT of 3HF is more favorable in the γ-CD hydrophobic cavity than in aqueous solution.

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.

Multiple Virtual Screening Strategies for the Discovery of Novel Compounds Active Against Dengue Virus: A Hit Identification Study

Dengue infection is caused by a mosquito-borne virus, particularly in children, which may even cause death. No effective prevention or therapeutic agents to cure this disease are available up to now. The dengue viral envelope (E) protein was discovered to be a promising target for inhibition in several steps of viral infection. Structure-based virtual screening has become an important technique to identify first hits in a drug screening process, as it is possible to reduce the number of compounds to be assayed, allowing to save resources. In the present study, pharmacophore models were generated using the common hits approach (CHA), starting from trajectories obtained from molecular dynamics (MD) simulations of the E protein complexed with the active inhibitor, flavanone (FN5Y). Subsequently, compounds presented in various drug databases were screened using the LigandScout 4.2 program. The obtained hits were analyzed in more detail by molecular docking, followed by extensive MD simulations of the complexes. The highest-ranked compound from this procedure was then synthesized and tested on its inhibitory efficiency by experimental assays.

Detailed chemical characterization and molecular modeling of serotonin inclusion complex with unmodified β-cyclodextrin

In this study, we analyzed the capability of unmodified β-cyclodextrin (β-CD) to form the stable complex with serotonin hydrochloride (SER), as an important neurotransmitter in the brain. The stable β-CD: SER formulation was prepared and characterized using spectroscopic, thermal, molecular docking, and molecular dynamics techniques, revealing the phenomenon of H-bond formations and the domination of hydrophobic forces between the host molecule and its guest via the amine group of SER and the narrow side of β-CD. The complexation mechanism was mainly enthalpy-driven, representing the improvement in SER photo-stability. Overall, the results highlighted the possibility to use this formulation with improved stability in clinical practice for treatment and prevention of various depressive conditions, such as anxiety disorders.

Effect of the Incorporation of Functionalized Cyclodextrins in the Liposomal Bilayer

Liposomes loaded with drug–cyclodextrin complexes are widely used as drug delivery systems, especially for species with low aqueous solubility and stability. Investigation of the intimate interactions of macrocycles with liposomes are essential for formulation of efficient and stable drug-in-cyclodextrin-in-liposome carriers. In this work, we reported the preparation of unilamellar vesicles of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) embedded with native β-cyclodextrin and two synthetic derivatives: heptakis(2,3,6-tri-O-methyl)-β-cyclodextrin (TMCD) and heptakis(2,3-di-O-acetyl)-β-cyclodextrin (DACD). We then studied the effect of these macrocycles on the liposomal size, membrane viscosity, and liposomal stability at different temperatures and concentrations. We observed that TMCD and DACD affected vesicle size and the change of size was related to CD concentration. Irrespective of its nature, the macrocycle established interactions with the phospholipidic head groups, preventing cyclodextrins to diffuse into the lipid bilayer, as confirmed by molecular dynamics simulations. Such supramolecular structuring improves liposome stability making these colloid systems promising carriers for biologically active compounds.

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.

Interactions of HLA-DR and Topoisomerase I Epitope Modulated Genetic Risk for Systemic Sclerosis

The association of systemic sclerosis with anti-Topoisomerase 1 antibody (ATASSc) with specific alleles of human leukocyte antigen (HLA)-DR has been observed among various ethnics. The anti-Topoisomerase 1 antibody is a common autoantibody in SSc with diffuse cutaneous scleroderma, which is one of the clinical subtypes of SSc. On the other hand, an immunodominant peptide of topoisomerase 1 (Top1) self-protein (residues 349-368) was reported to have strong association with ATASSc. In this study, molecular dynamics simulation was performed on the complexes of Top1 peptide with various HLA-DR subtypes divided into ATASSc-associated alleles (HLA-DRB1*08:02, HLA-DRB1*11:01 and HLA-DRB1*11:04), suspected allele (HLA-DRB5*01:02), and non-associated allele (HLA-DRB1*01:01). The unique interaction for each system was compared to the others in terms of dynamical behaviors, binding free energies and solvation effects. Our results showed that three HLA-DR/Top1 complexes of ATASSc association mostly exhibited high protein stability and increased binding efficiency without solvent interruption, in contrast to non-association. The suspected case (HLA-DRB5*01:02) binds Top1 as strongly as the ATASSc association case, which implied a highly possible risk for ATASSc development. This finding might support ATASSc development mechanism leading to a guideline for the treatment and avoidance of pathogens like Top1 self-peptide risk for ATASSc.

Cardol triene inhibits dengue infectivity by targeting kl loops and preventing envelope fusion

Dengue virus causes a global burden that specific chemotherapy has not been established. A previous report suggested that anacardic acid inhibited hepatitis C virus infection. Here, we explored structure activity relationship of anacardic acid, cardanol, and cardol homologues with anti-DENV cellular infectivities. Cardol triene showed the highest therapeutic index at 29.07 with the CC₅₀ and EC₅₀ of 207.30 ± 5.24 and 7.13 ± 0.72 µM, respectively. Moreover, we observed that the more unsaturated the hydrocarbon tail, the higher the CC₅₀s in all head groups. High CC₅₀s were also found in HepG-2, THP-1, and HEK-293 cell lines where cardol triene CC₅₀s were 140.27 ± 8.44, 129.77 ± 12.08, and 92.80 ± 3.93 µM, respectively. Cardol triene expressed pan-dengue inhibition with the EC₅₀s of 5.35 to 8.89 µM and kl loops of dengue envelope proteins were major targets. The strong binding energy at T48, E49, A50, P53, K128, V130, L135, M196, L198, Q200, W206, L207, I270, and L277 prevented cellular pH-dependent fusion. Zika virus kl loops were aligned in the closed position preventing cardol triene to bind and inhibit fusion and infectivity. This study showed for the first time that cardol triene had a potential for further development as anti-dengue inhibitors.

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.

Deciphering the Odorant Binding, Activation, and Discrimination Mechanism of Dhelobp21 from Dastarus Helophoroides

Odorant-binding proteins (OBPs) play a pivotal role in transporting odorants through the sensillar lymph of insect chemosensory sensilla and increasing the sensitivity of the olfactory system. To address the ligand binding, activation, and release mechanisms of OBPs, we performed a set of conventional molecular dynamics simulations for binding of the odorant-binding protein DhelOBP21 from Dastarcus helophoroides with 18 ligands (1-NPN and 17 volatiles), as well as four constant-pH molecular dynamics simulations. We found that the open pocket DhelOBP21 at pH 5.0 could bind volatiles and form a closed pocket complex via transformation of its N-terminus into regular Helix at pH 7.0 and vice versa. Moreover, the discrimination of volatiles (selectivity and promiscuity) was determined by the characteristics of both the volatiles and the ‘essential’ and ‘selective’ amino acid residues in OBP binding pockets, rather than the binding affinity of the volatiles. This study put forward a new hypothesis that during the binding of volatiles there are two transitions for the DhelOBP21 amino-terminus: pH- and odorant binding-dependent random-coil-to-helix. Another important finding is providing a framework for the exploration of the complete coil-to-helix transition process and theoretically analyzing its underlying causes at molecular level.

A novel flavanone derivative inhibits dengue virus fusion and infectivity

Dengue infection is a global burden affecting millions of world population. Previous studies indicated that flavanones were potential dengue virus inhibitors. We discovered that a novel flavanone derivative, 5-hydroxy-7-methoxy-6-methylflavanone (FN5Y), inhibited DENV2 pH-dependent fusion in cell-based system with strong binding efficiency to DENV envelope protein at K (P83, L107, K128, L198), K' (T48, E49, A50, L198, Q200, L277), X' (Y138, V354, I357), and Y' (V97, R99, N103, K246) by molecular dynamic simulation. FN5Y inhibited DENV2 infectivity with EC₅₀s (and selectivity index) of 15.99 ± 5.38 (>6.25), and 12.31 ± 1.64 (2.23) μM in LLC/MK2 and Vero cell lines, respectively, and inhibited DENV4 at 11.70 ± 6.04 (>8.55) μM. CC50s in LLC/MK2, HEK-293, and HepG2 cell lines at 72 h were higher than 100 μM. Time-of-addition study revealed that the maximal efficacy was achieved at early after infection corresponded with pH-dependent fusion. Inactivating the viral particle, interfering with cellular receptors, inhibiting viral protease, or the virus replication complex were not major targets of this compound. FN5Y could become a potent anti-flaviviral drug and can be structurally modified for higher potency using simulation to DENV envelope as a molecular target.

Theoretical and Experimental Studies on Inclusion Complexes of Pinostrobin and β-Cyclodextrins

Pinostrobin (PNS) belongs to the flavanone subclass of flavonoids which shows several biological activities such as anti-inflammatory, anti-cancerogenic, anti-viral and anti-oxidative effects. Similar to other flavonoids, PNS has a quite low water solubility. The purpose of this work is to improve the solubility and the biological activities of PNS by forming inclusion complexes with β-cyclodextrin (βCD) and its derivatives, heptakis-(2,6-di-O-methyl)-β-cyclodextrin (2,6-DMβCD) and (2-hydroxypropyl)-β-cyclodextrin (HPβCD). The A_L-type diagram of the phase solubility studies of PNS exhibited the formed inclusion complexes with the 1:1 molar ratio. Inclusion complexes were prepared by the freeze-drying method and were characterized by differential scanning calorimetry (DSC). Two-dimensional nuclear magnetic resonance (2D-NMR) and steered molecular dynamics (SMD) simulation revealed two different binding modes of PNS, i.e., its phenyl- (P-PNS) and chromone- (C-PNS) rings preferably inserted into the cavity of βCD derivatives whilst only one orientation of PNS, where the C-PNS ring is inside the cavity, was detected in the case of the parental βCD. All PNS/βCDs complexes had a higher dissolution rate than free PNS. Both PNS and its complexes significantly exerted a lowering effect on the IL-6 secretion in LPS-stimulated macrophages and showed a moderate cytotoxic effect against MCF-7 and HeLa cancer cell lines in vitro.

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.

Host-guest interactions between sildenafil and cyclodextrins: Spectrofluorometric study and molecular dynamic modeling

Sildenafil (SF) was included in cyclodextrins (CD) to enhance its solubility. Spectrofluorometry was used to confirm the complexation constant (stability constant). The orientations of SF inside the β-CDs and γ-CDs were fully illustrated. Molecular dynamics simulations were performed on two inclusion complexes (β-CD/SF and γ-CD/SF) in the aqueous system. The polar methylpiperazine group was found to locate inside the β-CD cavity, both in the wide and narrow side and was positioned 2°A from the center. In contrast, the methylpiperazine group did not fit well within the γ-CD cavity. Moreover, these results also confirmed hydrogen bonding that the highest number of bonding formed between the polar methylpiperazine sulfonyl structure and the hydroxyl group of β-CD. The simulated binding free energy of the methylpiperazine-β-CD inclusion complex (-6.01kcal/mol), featured a large contribution from electrostatic and van der Waals forces, which was the most stable complex. The association constant of β-CD/SF (12.3) was higher than γ-CD/SF (3.3) and confirmed with in silico measurements of binding free energy. In summary, SF forms a stable complex with β-CD.

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.

Physical properties and biological activities of hesperetin and naringenin in complex with methylated β-cyclodextrin

The aim of this work is to improve physical properties and biological activities of the two flavanones hesperetin and naringenin by complexation with β-cyclodextrin (β-CD) and its methylated derivatives (2,6-di-O-methyl-β-cyclodextrin, DM-β-CD and randomly methylated-β-CD, RAMEB). The free energies of inclusion complexes between hesperetin with cyclodextrins (β-CD and DM-β-CD) were theoretically investigated by molecular dynamics simulation. The free energy values obtained suggested a more stable inclusion complex with DM-β-CD. The vdW force is the main guest–host interaction when hesperetin binds with CDs. The phase solubility diagram showed the formation of a soluble complex of A_L type, with higher increase in solubility and stability when hesperetin and naringenin were complexed with RAMEB. Solid complexes were prepared by freeze-drying, and the data from differential scanning calorimetry (DSC) confirmed the formation of inclusion complexes. The data obtained by the dissolution method showed that complexation with RAMEB resulted in a better release of both flavanones to aqueous solution. The flavanones-β-CD/DM-β-CD complexes demonstrated a similar or a slight increase in anti-inflammatory activity and cytotoxicity towards three different cancer cell lines. The overall results suggested that solubilities and bioactivities of both flavanones were increased by complexation with methylated β-CDs.

Co-solvation effect on the binding mode of the α-mangostin/β-cyclodextrin inclusion complex

Cyclodextrins (CDs) have been extensively utilized as host molecules to enhance the solubility, stability and bioavailability of hydrophobic drug molecules through the formation of inclusion complexes. It was previously reported that the use of co-solvents in such studies may result in ternary (host:guest:co-solvent) complex formation. The objective of this work was to investigate the effect of ethanol as a co-solvent on the inclusion complex formation between α-mangostin (α-MGS) and β-CD, using both experimental and theoretical studies. Experimental phase-solubility studies were carried out in order to assess complex formation, with the mechanism of association being probed using a mathematical model. It was found that α-MGS was poorly soluble at low ethanol concentrations (0–10% v/v), but higher concentrations (10–40% v/v) resulted in better α-MGS solubility at all β-CD concentrations studied (0–10 mM). From the equilibrium constant calculation, the inclusion complex is still a binary complex (1:1), even in the presence of ethanol. The results from our theoretical study confirm that the binding mode is binary complex and the presence of ethanol as co-solvent enhances the solubility of α-MGS with some effects on the binding affinity with β-CD, depending on the concentration employed.