Arthitaya Meeprasert

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 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 IC₅₀ value of laninamivir.

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.

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.

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.

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.

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.