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 (ΔGbind) 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.

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.

Binding characteristics of pyrrole-scaffold hepatitis B virus capsid inhibitors and identification of novel potent compounds.

Hepatitis B virus (HBV) capsid assembly modulators (CAMs) are currently being evaluated in clinical trials as potential curative therapies for HBV. This study used in silico computational modeling to provide insights into the binding characteristics between the HBV core protein and two pyrrole-scaffold inhibitors, JNJ-6379 and GLP-26, both in the CAM-Normal (CAM-N) series. Molecular dynamics simulations showed that the pyrrole inhibitors displayed similar general binding-interaction patterns to NVR 3-778, another CAM-N, with hydrophobic interactions serving as the major driving force. However, consistent with their higher potency, the pyrrole inhibitors exhibited stronger nonpolar interactions with key residues in a solvent-accessible region as compared to NVR 3-778. This feature was facilitated by distinct hydrogen bond interactions of the pyrrole scaffold inhibitors with the residue 140 in chain B of the HBV core protein (L140B). Based on these findings, novel CAM-N compounds were designed to mimic the interaction with L140B residue while maximizing nonpolar interactions in the solvent-accessible region. Several 1H-pyrrole-2-carbonyl substituted pyrrolidine-based compounds with various hydrophobic side chains were synthesized and evaluated. Through analyses of the structure-activity and structure-druggability relations of a series of compounds, CU15 emerged as the most promising lead CAM-N compound, exhibiting sub-nanomolar potency and good pharmacokinetic profiles. Overall, the study demonstrated a practical approach to leverage computational methods for understanding key target binding features for rationale-based design, and for guiding the identification of novel compounds.

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.

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.

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.

Structural insight of DNA topoisomerases I from camptothecin-producing plants revealed by molecular dynamics simulations.

DNA topoisomerase I (Top1) catalyzes changes in DNA topology by cleaving and rejoining one strand of the double stranded (ds)DNA. Eukaryotic Top1s are the cellular target of the plant-derived anticancer indole alkaloid camptothecin (CPT), which reversibly stabilizes the Top1-dsDNA complex. However, CPT-producing plants, including Camptotheca acuminata, Ophiorrhiza pumila and Ophiorrhiza liukiuensis, are highly resistant to CPT because they possess point-mutated Top1. Here, the adaptive convergent evolution is reported between CPT production ability and mutations in their Top1, as a universal resistance mechanism found in all tested CPT-producing plants. This includes Nothapodytes nimmoniana, one of the major sources of CPT. To obtain a structural insight of the resistance mechanism, molecular dynamics simulations of CPT- resistant and -sensitive plant Top1s complexed with dsDNA and topotecan (a CPT derivative) were performed, these being compared to that for the CPT-sensitive human Top1. As a result, two mutations, Val617Gly and Asp710Gly, were identified in O. pumila Top1 and C. acuminata Top1, respectively. The substitutions at these two positions, surprisingly, are the same as those found in a CPT derivative-resistant human colon adenocarcinoma cell line. The results also demonstrated an increased linker flexibility of the CPT-resistant Top1, providing an additional explanation for the resistance mechanism found in CPT-producing plants. These mutations could reflect the long evolutionary adaptation of CPT-producing plant Top1s to confer a higher degree of resistance.

Molecular Dynamics Simulation Reveals the Selective Binding of Human Leukocyte Antigen Alleles Associated with Behçet's Disease.

Behçet's disease (BD), a multi-organ inflammatory disorder, is associated with the presence of the human leukocyte antigen (HLA) HLA-B*51 allele in many ethnic groups. The possible antigen involvement of the major histocompatibility complex class I chain related gene A transmembrane (MICA-TM) nonapeptide (AAAAAIFVI) has been reported in BD symptomatic patients. This peptide has also been detected in HLA-A*26:01 positive patients. To investigate the link of BD with these two specific HLA alleles, molecular dynamics (MD) simulations were applied on the MICA-TM nonapeptide binding to the two BD-associated HLA alleles in comparison with the two non-BD-associated HLA alleles (B*35:01 and A*11:01). The MD simulations were applied on the four HLA/MICA-TM peptide complexes in aqueous solution. As a result, stabilization for the incoming MICA-TM was found to be predominantly contributed from van der Waals interactions. The P2/P3 residue close to the N-terminal and the P9 residue at the C-terminal of the MICA-TM nonapeptide served as the anchor for the peptide accommodated at the binding groove of the BD associated HLAs. The MM/PBSA free energy calculation predicted a stronger binding of the HLA/peptide complexes for the BD-associated HLA alleles than for the non-BD-associated ones, with a ranked binding strength of B*51:01 > B*35:01 and A*26:01 > A*11:01. Thus, the HLAs associated with BD pathogenesis expose the binding efficiency with the MICA-TM nonapeptide tighter than the non-associated HLA alleles. In addition, the residues 70, 73, 99, 146, 147 and 159 of the two BD-associated HLAs provided the conserved interaction for the MICA-TM peptide binding.

In silico screening for potent inhibitors against the NS3/4A protease of hepatitis C virus.

Hepatitis C virus (HCV) infections are a serious viral health problem globally, causing liver cirrhosis and inflammation that can develop to hepatocellular carcinoma and death. Since the HCV NS3/4A protease complex cleaves the scissile peptide bond in the viral encoded polypeptide to release the non-structural proteins during the viral replication process, this protease is then an important target for drug design. The computer-aided drug design and screening targeted at NS3/4A protease of HCV were reviewed. In addition, using steered molecular dynamics simulations, potent inhibitors of the NS3/4A complex were searched for by screening the ZINC database based upon the hypothesis that a high rupture force indicates a high binding efficiency. Nine top-hit compounds (59500093, 59784724, 13527817, 26660256, 29482733, 25977181, 28005928, 13527826 and 13527826) were found that had the same or a greater maximum rupture force (and so assumed binding strength and inhibitory potency) than the four current drugs and so are potential candidates as anti- HCV chemotherapeutic agents. In addition, van der Waals interactions were found to be the main contribution in stabilizing the ligand- NS3/4A complex.

Binding pattern of the long acting neuraminidase inhibitor laninamivir towards influenza A subtypes H5N1 and pandemic H1N1.

Influenza A H5N1 and pH1N1 viruses have broadly emerged and become widespread in various countries around the world. Oseltamivir, the most commonly used antiviral drug against the seasonal and pandemic influenza viruses, is targeted at the viral neuraminidase (NA), but some isolates of this virus have become highly resistant to this drug. The novel long-acting drug, laninamivir, was recently developed to inhibit influenza A and B viruses of either the wild-type (WT) or the oseltamivir resistant mutant of NA. To understand the high efficiency of laninamivir, all-atom molecular dynamics simulations were performed on the WT and H274Y mutant of H5N1 and pH1N1 NAs with laninamivir bound. As a result, the novel drug was found to directly interact with 11 binding residues mainly through salt bridge and hydrogen bond formation (as also seen by electrostatic contribution). These are comprised of 7 of the catalytic residues (R118, D151, R152, R224, E276, R292 and R371), and 4 of the framework residues (E119, W178, E227 and E277). Laninamivir showed a similar binding pattern to all four NAs, but strong hydrogen bonding interactions were only found in the WT strain, with a slightly lowered contribution at some drug contact residues being observed in the H274Y mutation. This is in good agreement with the experimental data that the H274Y mutant has a small increase (1.3-7.5-fold, which was not statistically significant) in the IC50 value of laninamivir.

Molecular insight into the specific binding of ADP-ribose to the nsP3 macro domains of chikungunya and Venezuelan equine encephalitis viruses: molecular dynamics simulations and free energy calculations.

The outbreaks of chikungunya (CHIKV) and venezuelan equine encephalitis (VEEV) viral infections in humans have emerged or re-emerged in various countries of "Africa and southeast Asia", and "central and south America", respectively. At present, no drug or vaccine is available for the treatment and therapy of both viral infections, but the non-structural protein, nsP3, is a potential target for the design of potent inhibitors that fit at the adenosine-binding site of its macro domain. Here, so as to understand the fundamental basis of the particular interactions between the ADP-ribose bound to the nsP3 amino acid residues at the binding site, molecular dynamics simulations were applied. The results show that these two nsP3 domains share a similar binding pattern for accommodating the ADP-ribose. The ADP-ribose phosphate unit showed the highest degree of stabilization through hydrogen bond interactions with the nsP3 V33 residue and the consequent amino acid residues 110-114. The adenine base of ADP-ribose was specifically recognized by the conserved nsP3 residue D10. Additionally, the ribose and the diphosphate units were found to play more important roles in the CHIKV nsP3-ADP-ribose complex, while the ter-ribose was more important in the VEEV complex. The slightly higher binding affinity of ADP-ribose toward the nsP3 macro domain of VEEV, as predicted by the simulation results, is in good agreement with previous experimental data. These simulation results provide useful information to further assist in drug design and development for these two important viruses.