Effect of halogen substitutions on dUMP to stability of thymidylate synthase/dUMP/mTHF ternary complex using molecular dynamics simulation

Voriconazole Eye Drops: Enhanced Solubility and Stability through Ternary Voriconazole/Sulfobutyl Ether β-Cyclodextrin/Polyvinyl Alcohol Complexes

Voriconazole (VCZ) is a broad-spectrum antifungal agent used to treat ocular fungal keratitis. However, VCZ has low aqueous solubility and chemical instability in aqueous solutions. This study aimed to develop VCZ eye drop formulations using cyclodextrin (CD) and water-soluble polymers, forming CD complex aggregates to improve the aqueous solubility and chemical stability of VCZ. The VCZ solubility was greatly enhanced using sulfobutyl ether β-cyclodextrin (SBEβCD). The addition of polyvinyl alcohol (PVA) showed a synergistic effect on VCZ/SBEβCD solubilization and a stabilization effect on the VCZ/SBEβCD complex. The formation of binary VCZ/SBEβCD and ternary VCZ/SBEβCD/PVA complexes was confirmed by spectroscopic techniques and in silico studies. The 0.5% w/v VCZ eye drop formulations were developed consisting of 6% w/v SBEβCD and different types and concentrations of PVA. The VCZ/SBEβCD systems containing high-molecular-weight PVA prepared under freeze-thaw conditions (PVA-H hydrogel) provided high mucoadhesion, sustained release, good ex vivo permeability through the porcine cornea and no sign of irritation. Additionally, PVA-H hydrogel was effective against the filamentous fungi tested. The stability study revealed that our VCZ eye drops provide a shelf-life of more than 2.5 years at room temperature, while a shelf-life of only 3.5 months was observed for the extemporaneous Vfend^® eye drops.

Early Life Stage Folic Acid Deficiency Delays the Neurobehavioral Development and Cognitive Function of Rat Offspring by Hindering De Novo Telomere Synthesis

Early life stage folate status may influence neurodevelopment in offspring. The developmental origin of health and disease highlights the importance of the period of the first 1000 days (from conception to 2 years) of life. This study aimed to evaluate the effect of early life stage folic acid deficiency on de novo telomere synthesis, neurobehavioral development, and the cognitive function of offspring rats. The rats were divided into three diet treatment groups: folate-deficient, folate-normal, and folate-supplemented. They were fed the corresponding diet from 5 weeks of age to the end of the lactation period. After weaning, the offspring rats were still fed with the corresponding diet for up to 100 days. Neurobehavioral tests, folic acid and homocysteine (Hcy) levels, relative telomere length in brain tissue, and uracil incorporation in telomere in offspring were measured at different time points. The results showed that folic acid deficiency decreased the level of folic acid, increased the level of Hcy of brain tissue in offspring, increased the wrong incorporation of uracil into telomeres, and hindered de novo telomere synthesis. However, folic acid supplementation increased the level of folic acid, reduced the level of Hcy of brain tissue in offspring, reduced the wrong incorporation of uracil into telomeres, and protected de novo telomere synthesis of offspring, which was beneficial to the development of early sensory-motor function, spatial learning, and memory in adolescence and adulthood. In conclusion, early life stage folic acid deficiency had long-term inhibiting effects on neurodevelopment and cognitive function in offspring.

Solubility enhancement of poorly water soluble domperidone by complexation with the large ring cyclodextrin

The water solubility of domperidone (DMP) could be improved by complexation with large ring cyclodextrins (LR-CDs). LR-CDs contain a relatively hydrophobic cavity that is capable of entrapping the molecules to form inclusion complexes. The complex formation capability of mixture LR-CDs having a degree of polymerization (DP) of 22-48, with DMP was investigated. The phase solubility profile of mixture LR-CD/DMP was classified as A_N-type, resulting in increased DMP solubility in water by 3-fold. Various physicochemical techniques confirmed the mixture LR-CD/DMP complex formation. Single LR-CD with DP of 26, 27, 28, 29, 30, 33 and 34 (CD26 ~ CD34) were isolated from LR-CD mixtures using ODS column for HPLC separation. The CD33/DMP complex has demonstrated the most significant improvement compared to other single LR-CD complexes with a 2.7-fold increase in DMP solubility. The molecular dynamic result revealed that DMP formed stable complexes with CD33 by positioned fully encapsulated inside the cavity and covered by 13-14 subunits of CD33.

Metastatic colorectal cancer: Advances in the folate-fluoropyrimidine chemotherapy backbone

Notwithstanding recent treatment advances in metastatic colorectal cancer (mCRC), chemotherapy with a combination of a fluoropyrimidine and a folate agent, often 5-fluorouracil (5-FU) and leucovorin, remains the backbone of treatment regimens for the majority of patients with mCRC. This is despite a recent focus on molecular-targeted treatments and patient stratification according to mutational status or expression levels of specific genes. Intracellular folate concentration was discovered to be pivotal in the cytotoxic efficacy of 5-FU, paving the way to the current standard combination therapy approach. Subsequent discovery that systemic chemotherapy agents, such as irinotecan and oxaliplatin, can further increase the efficacy of 5-FU-based treatments led to the development of several combination chemotherapy regimens, including FOLFOX, FOLFIRI and FOLFOXIRI. Subsequent efforts to optimise 5-FU-based treatments have focused on 5-FU analogues, initially capecitabine and the combination drug tegafur/gimeracil/oteracil (S-1) and then TAS-102, which has recently been evaluated in phase 3 clinical trials for refractory colorectal cancer. Further approaches taken to improve the efficacy of 5-FU chemotherapy regimens have focused on optimising the route and dosing schedules and regulating folate metabolism. Pharmacokinetic variability caused by the requirement for metabolic conversion of leucovorin has been central to recent research, and the development of agents such as arfolitixorin which bypass the need for metabolic conversion remains promising for future therapeutic candidates. In this review, we summarise the evidence leading to the current treatment regimens employing 5-FU and leucovorin, focusing on recent approaches taken to optimise and refine treatments to improve clinical outcomes in patients with mCRC.

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.

Molecular dynamics simulations of sulfone derivatives in complex with DNA topoisomerase IIα ATPase domain

Human topoisomerase II alpha (TopoIIα) is a crucial enzyme involved in maintaining genomic integrity during the process of DNA replication and mitotic division. It is a vital therapeutic target for designing novel anticancer agents in targeted cancer therapy. Sulfones, members of organosulfur compounds, have been reported to possess various biological activities such as antimicrobial, anti-inflammatory, anti-HIV, anticancer, and antimalarial properties. In the present study, a series of sulfones was selected to evaluate their inhibitory activity against TopoIIα using computational approaches. Molecular docking results revealed that several sulfone analogs bind efficiently to the ATPase domain of TopoIIα. Among them, sulfones 18a, 60a, *4 b, *8 b, *3c, and 8c exhibit higher binding affinity than the known TopoII inhibitor, salvicine. Molecular dynamics simulations and free energy calculations based on MM/PB(GB)SA method demonstrated that sulfone *8 b strongly interacts with amino acid residues in the ATP-binding pocket (E87, N91, D94, I125, I141, F142, S149, G161, and A167), driven mainly by an electrostatic attraction and a strong H-bond formation at G161 residue. Altogether, the obtained results predicted that sulfones could have a high potential to be a lead molecule for targeting TopoIIα.Communicated by Ramaswamy H. Sarma.

dUMP/F-dUMP Binding to Thymidylate Synthase: Human Versus Mycobacterium tuberculosis

Thymidylate synthase is an enzyme that catalyzes deoxythymidine monophosphate (dTMP) synthesis from substrate deoxyuridine monophosphate (dUMP). Thymidylate synthase of Mycobacterium tuberculosis (MtbThyX) is structurally distinct from its human analogue human thymidylate synthase (hThyA), thus drawing attention as an attractive drug target for combating tuberculosis. Fluorodeoxyuridylate (F-dUMP) is a successful inhibitor of both MtbThyX and hThyA, thus limited by poor selectivity. Understanding the dynamics and energetics associated with substrate/inhibitor binding to thymidylate synthase in atomic details remains a fundamental unsolved problem, which is necessary for a new selective inhibitor design. Structural studies of MtbThyX and hThyA bound substrate/inhibitor complexes not only revealed the extensive specific interaction network between protein and ligands but also opened up the possibility of directly computing the energetics of the substrate versus inhibitor recognition. Using experimentally determined structures as a template, we report extensive computer simulations (∼4.5 μs) that allow us to quantitatively estimate ligand selectivity (dUMP vs F-dUMP) by MtbThyX and hThyA. We show that MtbThyX prefers deprotonated dUMP (enolate form) as the substrate, whereas hThyA binds to the keto form of dUMP. Computed energetics clearly show that MtbThyX is less selective between dUMP and F-dUMP, favoring the latter, relative to hThyA. The simulations reveal the role of tyrosine at position 135 (Y135) of hThyA in amplifying the selectivity. The protonation state of the pyrimidine base of the ligand (i.e., keto or enolate) seems to have no role in MtbThyX ligand selectivity. A molecular gate (consists of Y108, K165, H203, and a water molecule) restricts water accessibility and offers a desolvated dry ligand-binding pocket for MtbThyX. The ligand-binding pocket of hThyA is relatively wet and exposed to bulk water.

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.

Atomistic mechanisms underlying the activation of the G protein-coupled sweet receptor heterodimer by sugar alcohol recognition

The human T1R2-T1R3 sweet taste receptor (STR) plays an important role in recognizing various low-molecular-weight sweet-tasting sugars and proteins, resulting in the release of intracellular heterotrimeric G protein that in turn leads to the sweet taste perception. Xylitol and sorbitol, which are naturally occurring sugar alcohols (polyols) found in many fruits and vegetables, exhibit the potential caries-reducing effect and are widely used for diabetic patients as low-calorie sweeteners. In the present study, computational tools were applied to investigate the structural details of binary complexes formed between these two polyols and the T1R2-T1R3 heterodimeric STR. Principal component analysis revealed that the Venus flytrap domain (VFD) of T1R2 monomer was adapted by the induced-fit mechanism to accommodate the focused polyols, in which α-helical residues 233-268 moved significantly closer to stabilize ligands. This finding likely suggested that these structural transformations might be the important mechanisms underlying polyols-STR recognitions. The calculated free energies also supported the VFD of T1R2 monomer as the preferential binding site for such polyols, rather than T1R3 region, in accord with the lower number of accessible water molecules in the T1R2 pocket. The E302 amino acid residue in T1R2 was found to be the important recognition residue for polyols binding through a strongly formed hydrogen bond. Additionally, the binding affinity of xylitol toward the T1R2 monomer was significantly higher than that of sorbitol, making it a sweeter tasting molecule.

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.

A flap motif in human serine hydroxymethyltransferase is important for structural stabilization, ligand binding, and control of product release

Human cytosolic serine hydroxymethyltransferase (hcSHMT) is a promising target for anticancer chemotherapy and contains a flexible "flap motif" whose function is yet unknown. Here, using size-exclusion chromatography, analytical ultracentrifugation, small-angle X-ray scattering (SAXS), molecular dynamics (MD) simulations, and ligand-binding and enzyme-kinetic analyses, we studied the functional roles of the flap motif by comparing WT hcSHMT with a flap-deleted variant (hcSHMT/Δflap). We found that deletion of the flap results in a mixture of apo-dimers and holo-tetramers, whereas the WT was mostly in the tetrameric form. MD simulations indicated that the flap stabilizes structural compactness and thereby enhances oligomerization. The hcSHMT/Δflap variant exhibited different catalytic properties in (6S)-tetrahydrofolate (THF)-dependent reactions compared with the WT but had similar activity in THF-independent aldol cleavage of β-hydroxyamino acid. hcSHMT/Δflap was less sensitive to THF inhibition than the WT (K_i of 0.65 and 0.27 mm THF at pH 7.5, respectively), and the THF dissociation constant of the WT was also 3-fold lower than that of hcSHMT/Δflap, indicating that the flap is important for THF binding. hcSHMT/Δflap did not display the burst kinetics observed in the WT. These results indicate that, upon removal of the flap, product release is no longer the rate-limiting step, implying that the flap is important for controlling product release. The findings reported here improve our understanding of the functional roles of the flap motif in hcSHMT and provide fundamental insight into how a flexible loop can be involved in controlling the enzymatic reactions of hcSHMT and other enzymes.

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.

Computational screening of chalcones acting against topoisomerase IIα and their cytotoxicity towards cancer cell lines

Targeted cancer therapy has become one of the high potential cancer treatments. Human topoisomerase II (hTopoII), which catalyzes the cleavage and rejoining of double-stranded DNA, is an important molecular target for the development of novel cancer therapeutics. In order to diversify the pharmacological activity of chalcones and to extend the scaffold of topoisomerase inhibitors, a series of chalcones was screened against hTopoIIα by computational techniques, and subsequently tested for their in vitro cytotoxicity. From the experimental IC₅₀ values, chalcone 3d showed a high cytotoxicity with IC₅₀ values of 10.8, 3.2 and 21.1 µM against the HT-1376, HeLa and MCF-7 cancer-derived cell lines, respectively, and also exhibited an inhibitory activity against hTopoIIα-ATPase that was better than the known inhibitor, salvicine. The observed ligand-protein interactions from a molecular dynamics study affirmed that 3d strongly interacts with the ATP-binding pocket residues. Altogether, the newly synthesised chalcone 3d has a high potential to serve as a lead compound for topoisomerase 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.

Therapeutic use of fluorinated nucleosides – progress in patents

Fluorinated nucleosides constitute a large class of chemotherapeutics approved for clinical use. The pharmacokinetic and pharmacodynamic properties of a drug can be affected, as a consequence of modulation of electronic, lipophilic and steric parameters, by the introduction of fluorine into the structure of drug-like molecule. Herein, we focus on fluorinated-nucleoside analogs, their therapeutic use and applications based on the patent literature from 2014 to 2018. We briefly discuss the clinical properties of anticancer and antiviral fluorine-containing nucleos(t)ides US FDA-approved or in development, and highlight their resistance mechanisms and limitations in the clinic. We emphasize patent inventions related to improved synthetic methods toward selected nucleos(t)ide analogs including the phosphoramidate sofosbuvir and 18F-labeled nucleosides FLT and FMAU, used as a 18F-PET tracers.

Nitric oxide promotes cancer cell dedifferentiation by disrupting an Oct4:caveolin-1 complex: A new regulatory mechanism for cancer stem cell formation

Cancer stem cells (CSCs) are unique populations of cells that can self-renew and generate different cancer cell lineages. Although CSCs are believed to be a promising target for novel therapies, the specific mechanisms by which these putative therapeutics could intervene are less clear. Nitric oxide (NO) is a biological mediator frequently up-regulated in tumors and has been linked to cancer aggressiveness. Here, we search for targets of NO that could explain its activity. We find that it directly affects the stability and function of octamer-binding transcription factor 4 (Oct4), known to drive the stemness of lung cancer cells. We demonstrated that NO promotes the CSC-regulatory activity of Oct4 through a mechanism that involves complex formation between Oct4 and the scaffolding protein caveolin-1 (Cav-1). In the absence of NO, Oct4 forms a molecular complex with Cav-1, which promotes the ubiquitin-mediated proteasomal degradation of Oct4. NO promotes Akt-dependent phosphorylation of Cav-1 at tyrosine 14, disrupting the Cav-1:Oct4 complex. Site-directed mutagenesis and computational modeling studies revealed that the hydroxyl moiety at tyrosine 14 of Cav-1 is crucial for its interaction with Oct4. Both removal of the hydroxyl via mutation to phenylalanine and phosphorylation lead to an increase in binding free energy (ΔG_bind) between Oct4 and Cav-1, destabilizing the complex. Together, these results unveiled a novel mechanism of CSC regulation through NO-mediated stabilization of Oct4, a key stem cell transcription factor, and point to new opportunities to design CSC-related therapeutics.

Polarizable Empirical Force Field for Halogen-Containing Compounds Based on the Classical Drude Oscillator

The quality of the force field is crucial to ensure the accuracy of simulations used in molecular modeling, including computer-aided drug design (CADD). To perform more accurate modeling and simulations of halogenated molecules, in this study the polarizable force field based on the classical Drude oscillator model was extended to both aliphatic and aromatic systems using halogenated ethane and benzene model compounds for the halogens F, Cl, Br, and I. The force field parameters were optimized targeting quantum mechanical dipole moments, water interactions, and molecular polarizabilities as well as experimental observables, including enthalpies of vaporization, molecular volumes, hydration free energies, and dielectric constants. The developed halogenated polarizable force field is capable of reproducing QM relative energies and geometries of both halogen bonds and halogen-hydrogen bond donor interactions at an unprecedented level due to the inclusion of a virtual particle and anisotropic atomic polarizability on the halogen and, notably, the inclusion of Lennard-Jones parameters on the halogen Drude particle. The model was validated on the basis of its ability to accurately reproduce pure solvent properties for halogenated naphthalenes and alkanes, including species analogous to those used as refrigerants. Accordingly, it is anticipated that the model will be applicable for the study of halogenated derivatives in CADD as well as in other chemical and biophysical studies.

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.

Loss of Uracil DNA Glycosylase Selectively Resensitizes p53-Mutant and -Deficient Cells to 5-FdU

Thymidylate synthase (TS) inhibitors including fluoropyrimidines [e.g., 5-Fluorouracil (5-FU) and 5-Fluorodeoxyuridine (5-FdU, floxuridine)] and antifolates (e.g., pemetrexed) are widely used against solid tumors. Previously, we reported that shRNA-mediated knockdown (KD) of uracil DNA glycosylase (UDG) sensitized cancer cells to 5-FdU. Because p53 has also been shown as a critical determinant of the sensitivity to TS inhibitors, we further interrogated 5-FdU cytotoxicity after UDG depletion with regard to p53 status. By analyzing a panel of human cancer cells with known p53 status, it was determined that p53-mutated or -deficient cells are highly resistant to 5-FdU. UDG depletion resensitizes 5-FdU in p53-mutant and -deficient cells, whereas p53 wild-type (WT) cells are not affected under similar conditions. Utilizing paired HCT116 p53 WT and p53 knockout (KO) cells, it was shown that loss of p53 improves cell survival after 5-FdU, and UDG depletion only significantly sensitizes p53 KO cells. This sensitization can also be recapitulated by UDG depletion in cells with p53 KD by shRNAs. In addition, sensitization is also observed with pemetrexed in p53 KO cells, but not with 5-FU, most likely due to RNA incorporation. Importantly, in p53 WT cells, the apoptosis pathway induced by 5-FdU is activated independent of UDG status. However, in p53 KO cells, apoptosis is compromised in UDG-expressing cells, but dramatically elevated in UDG-depleted cells. Collectively, these results provide evidence that loss of UDG catalyzes significant cell death signals only in cancer cells mutant or deficient in p53.Implications: This study reveals that UDG depletion restores sensitivity to TS inhibitors and has chemotherapeutic potential in the context of mutant or deficient p53. Mol Cancer Res; 16(2); 212-21. ©2017 AACR.

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.

Molecular dynamics investigation of stereoselective inhibition mechanism of HIF-2α/ARNT heterodimer

Hypoxia-inducible factors (HIFs) are heterodimeric transcription factors related with the onset and progression of solid tumors. Studies demonstrated a class of tetrazole containing chiral inhibitors could stereoselectively disrupt the HIF-2 dimerization and reduce the target gene expression. However, the dynamical features and structural motifs of the HIF-2 heterodimer caused by the binding of enantiomers have not been rationalized at the atomistic level. In this work, molecular dynamics (MD) simulations combined with adaptive steered MD (ASMD) simulations were used to investigate stereoselective interrupting mechanism of HIF-2. Our results decipher that the binding of ligand A (S, R)-24 begets the significant conformation changes of β-sheets and interrupts the HIF-2α/ARNT heterodimerization, which may be attributed to the disruption of the hydrogen bond and salt bridge interactions formed by the 4 foremost residues (Asp240, Arg247, Glu362, and Arg366) and the destruction of hydrophobic interactions on the binding interface. By contrast, the binding of ligand B (R, S)-24 does not disrupt protein dimerization and causes the motion of Fα helix in HIF-2α PAS-B domain to further change the major tunnel for ligand ingress and engress. The present work provides important molecular-level insight into the effect of the binding enantiomers on HIF-2 heterodimerization and bridges the gap between theory and the experimental results, which may conduce to develop highly potent antagonists for intervening the HIF-2-driven tumors.

Fluorinated nucleosides as an important class of anticancer and antiviral agents

Fluorine-containing nucleoside analogs (NAs) represent a significant class of the US FDA-approved chemotherapeutics widely used in the clinic. The incorporation of fluorine into drug-like agents modulates lipophilic, electronic and steric parameters, thus influencing pharmacodynamic and pharmacokinetic properties of drugs. Fluorine can block oxidative metabolism of drugs and the formation of undesired metabolites by changing H-bonding interactions. In this review, we focus our attention on chemical fluorination reagents and methods used in the NAs field, including positron emission tomography radiochemistry. We briefly discuss both the cellular biology and clinical properties of FDA-approved and fluorine-containing nucleoside/nucleotide analogs in development as well as common resistance mechanisms associated with their use. Finally, we emphasize pronucleotide strategies used to improve therapeutic outcome of NAs in the clinic.

Effect of D168V mutation in NS3/4A HCV protease on susceptibilities of faldaprevir and danoprevir

Hepatitis C virus (HCV) is a serious cause of liver inflammation, cirrhosis and the development of hepatocellular carcinoma. Its NS3/4A serine protease functions to cleave a specific peptide bond, which is an important step in HCV replication. Thus the NS3/4A protease has become one of the main drug-targets in the design and development of anti-HCV agents. Unfortunately, high mutation rates in HCV have been reported due to the lack of RNA proofreading activity resulting in drug resistance. Herein, all-atom molecular dynamics simulations were employed to understand and illustrate the effects of the NS3/4A D168V mutation on faldaprevir (FDV) and danoprevir (DNV) binding efficiency. The D168V mutation was shown to interrupt the hydrogen bonding network of Q80R155D168R123 embedded in the extended S2 and partial S4 subsites of the NS3 protein and as a result the R123 side chain was displaced and moved out from the binding pocket. By means of MM/PBSA and MM/GBSA binding free energy calculations, the FDV and DNV binding affinities were shown to be significantly reduced by ∼10-15 kcal mol^-1 and ∼4-9 kcal mol^-1 relative to the wild-type complexes, respectively, which somewhat agrees with the experimental resistance folds.

Simultaneous, But Not Consecutive, Combination With Folinate Salts Potentiates 5-Fluorouracil Antitumor Activity In Vitro and In Vivo

The combination of folinate salts to 5-fluoruracil (5-FU)-based schedules is an established clinical routine in the landscape of colorectal cancer treatment. The aim of this study was to investigate the pharmacological differences between the sequential administration of folinate salts (1 h before, as in clinical routine) followed by 5-FU and the simultaneous administration of both drugs. Proliferation and apoptotic assays were performed on human colon cancer cells exposed to 5-FU, calcium (CaLV), or disodium (NaLV) levofolinate or their simultaneous and sequential combination for 24 and 72 h. TYMS and SLC19A1 gene expression was performed with real-time PCR. In vivo experiments were performed in xenografted nude mice, which were treated with 5-FU escalating doses and CaLV or NaLV alone or in simultaneous and sequential combination. The simultaneous combination of folinate salts and 5-FU was synergistic (NaLV) or additive (CaLV) in a 24-h treatment in both cell lines. In contrast, the sequential combination of both folinate salts and 5-FU was antagonistic at 24 and 72 h. The simultaneous combination of 5-FU and NaLV or CaLV inhibited TYMS gene expression at 24 h, whereas the sequential combination reduced SLC19A1 gene expression. In vivo experiments confirmed the enhanced antitumor activity of the 5-FU + NaLV simultaneous combination with a good toxicity profile, whereas the sequential combination with CaLV failed to potentiate 5-FU activity. In conclusion, only the simultaneous, but not the consecutive, in vitro and in vivo combination of 5-FU and both folinate salt formulations potentiated the antiproliferative effects of the drugs.

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.

Two novel antimicrobial defensins from rice identified by gene coexpression network analyses

Defensins form an antimicrobial peptides (AMP) family, and have been widely studied in various plants because of their considerable inhibitory functions. However, their roles in rice (Oryza sativa L.) have not been characterized, even though rice is one of the most important staple crops that is susceptible to damaging infections. Additionally, a previous study identified 598 rice genes encoding cysteine-rich peptides, suggesting there are several uncharacterized AMPs in rice. We performed in silico gene expression and coexpression network analyses of all genes encoding defensin and defensin-like peptides, and determined that OsDEF7 and OsDEF8 are coexpressed with pathogen-responsive genes. Recombinant OsDEF7 and OsDEF8 could form homodimers. They inhibited the growth of the bacteria Xanthomonas oryzae pv. oryzae, X. oryzae pv. oryzicola, and Erwinia carotovora subsp. atroseptica with minimum inhibitory concentration (MIC) ranging from 0.6 to 63μg/mL. However, these OsDEFs are weakly active against the phytopathogenic fungi Helminthosporium oryzae and Fusarium oxysporum f.sp. cubense. This study describes a useful method for identifying potential plant AMPs with biological activities.

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.

Inclusion complexation of pinostrobin with various cyclodextrin derivatives

Pinostrobin (PNS) is one of the important flavonoids and can be abundantly found in the rhizomes of fingerroot (Boesenbergia rotrunda) and galangal (Alpinia galangal and Alpinia officinarum), the herbal basis of Southeast Asian cooking. Similar to other flavonoids, PNS exhibits anti-oxidative, anti-inflammatory and anti-cancer properties. However, this compound has an extremely low water solubility that limits its use in pharmaceutical applications. Beta-cyclodextrin (βCD) and its derivatives, 2,6-dimethyl-βCD (2,6-DMβCD) and the three hydroxypropyl-βCDs (2-HPβCD, 6-HPβCD and 2,6-DHPβCD), have unique properties that enhance the stability and solubility of such low-soluble guest molecules. In the present study, molecular dynamics simulations were applied to investigate the dynamics and stability of PNS inclusion complexes with βCD and its derivatives (2,6-DMβCD, 2,6-DHPβCD, 2-HPβCD and 6-HPβCD). PNS was able to form complexes with βCD and all four of its derivatives by either the chromone (C-PNS) or phenyl (P-PNS) ring dipping toward the cavity. According to the molecular mechanics-generalized Born surface area binding free energy values, the stability of the different PNS/βCD complexes was ranked as 2,6-DHPβCD>2,6-DMβCD>2-HPβCD>6-HPβCD>βCD. These theoretical results were in good agreement with the stability constants that had been determined by the solubility method.

Expression of Folate Pathway Genes in Stage III Colorectal Cancer Correlates with Recurrence Status Following Adjuvant Bolus 5-FU-Based Chemotherapy

Colorectal cancer is commonly treated with 5-fluorouracil and 5-formyltetrahydrofolate (leucovorin). Metabolic action of leucovorin requires several enzymatic steps that are dependent on expression of corresponding coding genes. To identify folate pathway genes with possible impact on leucovorin metabolism, a retrospective study was performed on 193 patients with stage III colorectal cancer. Relative expression of 22 genes putatively involved in leucovorin transport, polyglutamation and metabolism was determined in tumor and mucosa samples using quantitative real-time polymerase chain reaction. After surgery, patients received adjuvant 5-fluorouracil-based bolus chemotherapy with leucovorin during six months, and were followed for 3 to 5 years. Cox regression analysis showed that high tumoral expression of the genes SLC46A1/PCFT (proton-coupled folate transporter) and SLC19A1/RFC-1 (reduced folate carrier 1) correlated significantly (p < 0.001 and p < 0.01, respectively) with a decreased risk of recurrent disease, measured as disease-free survival (DFS). These two genes are involved in the transport of folates into the cells and each functions optimally at a different pH. We conclude that SLC46A1/PCFT and SLC19A1/RFC-1 are associated with DFS of patients with colorectal cancer and hypothesize that poor response to 5-fluorouracil plus leucovorin therapy in some patients may be linked to low expression of these genes. Such patients might need a more intensified therapeutic approach than those with high gene expression. Future prospective studies will determine if the expression of any of these genes can be used to predict response to leucovorin.

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.

A pharmacokinetic and pharmacodynamic investigation of Modufolin® compared to Isovorin® after single dose intravenous administration to patients with colon cancer: a randomized study

PURPOSE

Leucovorin is commonly used as folate supplement in 5-fluorouracil-based chemotherapy, but needs to be converted to active 5,10-methylenetetrahydrofolate (methyleneTHF) intracellularly. This provides for interindividual differences. MethyleneTHF has recently been developed into the stable, distributable drug, Modufolin®. The aim was to compare the concentration of folate metabolites in tumor, mucosa, and plasma of patients with colon cancer after administration of Modufolin® or Isovorin® (levo-leucovorin).

METHODS

Thirty-two patients scheduled for colon resection were randomized to receive Modufolin® or Isovorin® at dosage of 60 or 200 mg/m². The study drug was given as one i.v. bolus injection after anesthesia. Plasma was collected for pharmacokinetic (PK) analysis before, during, and after surgery. Tissue biopsies were collected at surgery. Folate metabolites were analyzed by LC-MS/MS.

RESULTS

MethyleneTHF and THF concentrations were significantly higher in mucosa (p < 0.01, both dosages) and tumors (p < 0.01, 200 mg/m²) after Modufolin® as compared to Isovorin® administration. The results correlated with PK observations. The Modufolin® to Isovorin® C(max) ratio for methyleneTHF was 113 at 200 mg/m² and 52 at 60 mg/m²; the AUC(last) ratios were 17 and 9, respectively. The THF plasma concentrations were also higher after Modufolin® administration (C(max) ratio 23, AUC(last) ratio 13 at 200 mg/m²; C(max) ratio 15, AUC(last) ratio 11 at 60 mg/m²).

CONCLUSION

Modufolin® administration resulted in significantly higher methyleneTHF levels than Isovorin® and may potentially increase the efficacy of 5-fluorouracil-based chemotherapy. The results encourage further evaluation of Modufolin® as a substitute to Isovorin® including the potential clinical benefits.

Key binding and susceptibility of NS3/4A serine protease inhibitors against hepatitis C virus

Hepatitis C virus (HCV) causes an infectious disease that manifests itself as liver inflammation, cirrhosis, and can lead to the development of liver cancer. Its NS3/4A serine protease is a potent target for drug design and development since it is responsible for cleavage of the scissile peptide bonds in the polyprotein important for the HCV life cycle. Herein, the ligand-target interactions and the binding free energy of the four current NS3/4A inhibitors (boceprevir, telaprevir, danoprevir, and BI201335) were investigated by all-atom molecular dynamics simulations with three different initial atomic velocities. The per-residue free energy decomposition suggests that the key residues involved in inhibitor binding were residues 41-43, 57, 81, 136-139, 155-159, and 168 in the NS3 domain. The van der Waals interactions yielded the main driving force for inhibitor binding at the protease active site for the cleavage reaction. In addition, the highest number of hydrogen bonds was formed at the reactive P1 site of the four studied inhibitors. Although the hydrogen bond patterns of these inhibitors were different, their P3 site was most likely to be recognized by the A157 backbone. Both molecular mechanic (MM)/Poisson-Boltzmann surface area and MM/generalized Born surface area approaches predicted the relative binding affinities of the four inhibitors in a somewhat similar trend to their experimentally derived biological activities.

Enhanced stability of a naringenin/2,6-dimethyl β-cyclodextrin inclusion complex: molecular dynamics and free energy calculations based on MM- and QM-PBSA/GBSA

The structure, dynamic behavior and binding affinity of the inclusion complexes between naringenin and the two cyclodextrins (CDs), β-CD and its 2,6-dimethyl derivative (DM-β-CD), were theoretically studied by multiple molecular dynamics simulations and free energy calculations. Naringenin most likely prefers to bind with CDs through the phenyl ring. Although a lower hydrogen bond formation of naringenin with the 3-hydroxyl group of DM-β-CD (relative to β-CD) was observed, the higher cavity could encapsulate almost the whole naringenin molecule. In contrast for the naringenin/β-CD complex, the phenyl ring feasibly passed through the primary rim resulting in the chromone ring binding inside instead. MM-PBSA/GBSA and QM-PBSA/GBSA binding free energies strongly suggested a greater stability of the naringenin/DM-β-CD inclusion complex. Van der Waals force played an important role as the key guest-host interaction for the complexation between naringenin and each cyclodextrin.