Search

The present study investigates deep eutectic solvents (DESs) as potential media for enzymatic hydrolysis. A series of ternary ammonium and phosphonium-based DESs were prepared at different molar ratios by mixing with aqueous glycerol (85%). The physicochemical properties including surface tension, conductivity, density, and viscosity were measured at a temperature range of 298.15 K – 363.15 K. The eutectic points were highly influenced by the variation of temperature. The eutectic point of the choline chloride: glycerol: water (ratio of 1: 2.55: 2.28) and methyltriphenylphosphonium bromide:glycerol:water (ratio of 1: 4.25: 3.75) is 213.4 K and 255.8 K, respectively. The stability of the lipase enzyme isolated from porcine pancreas (PPL) and Rhizopus niveus (RNL) toward hydrolysis in ternary DESs medium was investigated. The PPL showed higher activity compared to the RNL in DESs. Molecular docking simulation of the selected DES with the substrate (p-nitrophenyl palmitate) toward PPL was also reported. It is worth noting that ternary DES systems would be viable lipase activators in hydrolysis reactions.

Alzheimer’s disease (AD) is considered one of the greatest global public burdens. Pathophysiology of AD is proposed to be associated with reduced levels of the neurotransmitter acetylcholine (ACh). Cholinesterase enzymes, namely acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) cleave ACh via hydrolysis. Cholinesterase inhibitors (ChEIs) are the main group of drugs currently used for the treatment of AD. Novel chromone-2-carboxamido-alkylamines (7–18) were designed, synthesized, and evaluated for cholinesterase inhibitory activity. The compounds exhibited potent AChE inhibitory activities at micromolar range (IC50 0.09–9.16 µM) and demonstrated weak BChE inhibitory activities (IC50 12.09–44.56 µM). Compound 14 (IC50 0.09 ± 0.02 µM) was the most potent AChEI in this series; it showed higher activity than the clinical used drug tacrine. Enzyme kinetic study suggested that 14 was an uncompetitive inhibitor. Molecular docking study revealed that 14 was a dual-binding site inhibitor. Compound 14 did not induce any concentration-related cytotoxic effect against SH-SY5Y cells. It also showed neuroprotective effect in the cell line. Chromone-2-carboxamido-alkylamines can be promising lead compounds for development of anti-Alzheimer’s agents.

Dengue virus (DENV) infection causes mild to severe illness in humans that can lead to fatality in severe cases. Currently, no specific drug is available for the treatment of DENV infection. Thus, the development of an anti-DENV drug is urgently required. Cordycepin (3′-deoxyadenosine), which is a major bioactive compound in Cordyceps (ascomycete) fungus that has been used for centuries in Chinese traditional medicine, was reported to exhibit antiviral activity. However, the anti-DENV activity of cordycepin is unknown. We hypothesized that cordycepin exerts anti-DENV activity and that, as an adenosine derivative, it inhibits DENV replication. To test this hypothesis, we investigated the anti-DENV activity of cordycepin in DENV-infected Vero cells. Cordycepin treatment significantly decreased DENV protein at a half-maximal effective concentration (EC50) of 26.94 μM. Moreover, DENV RNA was dramatically decreased in cordycepin-treated Vero cells, indicating its effectiveness in inhibiting viral RNA replication. Via in silico molecular docking, the binding of cordycepin to DENV non-structural protein 5 (NS5), which is an important enzyme for RNA synthesis, at both the methyltransferase (MTase) and RNA-dependent RNA polymerase (RdRp) domains, was predicted. The results of this study demonstrate that cordycepin is able to inhibit DENV replication, which portends its potential as an anti-dengue therapy.

A total of 20 of isolates of lactic acid bacteria (LAB) were selected and screened for antagonistic activity against clinical strains of 30 clinical isolates of extremely drug-resistant (XDR) Acinetobacter baumannii using the well diffusion assay method. Results showed that 50% of the highly LAB strains possessed inhibitory activity against (up to 66%) of the XDR A. baumannii strains tested. The supernatant of the twenty LAB strains was subjected to gas chromatography mass spectrometry (GCMS) revealed that the common compound found in the active isolates against XDR A. baumannii was 3-Isobutyl-2,3,6,7,8,8a-hexahydropyrrolo[1,2-a]pyrazine-1,4-dione, a known potential diketopiperazine group. The molecular docking study against potential antibacterial targets with selected ligands was performed to predict the binding mode of interactions, which is responsible for antibacterial activity. The docking analysis of the potent compounds supported the potential antibacterial activity exhibiting high inhibition constant and binding affinity in silico.

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

A series of 2-aryl-3-(arylideneamino)-1,2-dihydroquinazoline-4(3H)-ones were evaluated as inhibitors of acetylcholinesterase (AChE), butyrylcholinesterase (BuChE) and self-induced β-amyloid (Aβ) aggregation. All the compounds were found to inhibit both forms of cholinesterase (IC50 in the range 4–32 μM) with some selectivity for BuChE. Most of the compounds also showed self-induced Aβ aggregation inhibitory activities, which were comparable or higher than those obtained for reference compounds, curcumin and myricetin. Docking and molecular dynamics (MD) simulation experiments suggested that the compounds are able to disrupt the dimer form of Aβ.

Breast cancer is the second most common cancer among women worldwide. The Wnt??-catenin pathway appears to be deregulated in most cancer cells including breast cancer. The role of zerumbone, the active sesquiterpene from Zingiber zerumbet Roscoe, on the Wnt??-catenin pathway is relatively unknown, especially detailed molecular studies have yet to be published. Using the Chemistry at HARvard Macromolecular Mechanics (CHARMm) force field-based docking protocol, CDOCKER, the molecular interactions between zerumbone and key proteins of the Wnt??-catenin pathway were evaluated in this study. The results suggest that zerumbone has a strong affinity for free ?-catenin in the cytoplasm, as well as the ?-catenin?transcription factor 4 complex in the nucleus. The overall hydrophobic nature of zerumbone allowed its interaction with other hydrophobic residues, such as Trp383, while its active ?,?-unsaturated carbonyl facilitated its interaction with positively charged residues, such as Lys345, Arg386 and Asn415 in the ?-catenin binding pocket. However, the Wnt protein and its frizzled receptor showed no attraction to zerumbone.

Enterovirus A71 (EV-A71) causes hand, foot and mouth disease epidemics with neurological complications and fatalities. However, the neuropathogenesis of EV-A71 remains poorly understood. In mice, adaptation and virulence determinants have been mapped to mutations at VP2-149, VP1-145 and VP1-244. We investigate how these amino acids alter heparin-binding phenotype and shapes EV-A71 virulence in one-day old mice. We constructed six viruses with varying residues at VP1-98, VP1-145 (which are both heparin-binding determinants) and VP2-149 (based on the wild type 149K/98E/145Q, termed KEQ) to generate KKQ, KKE, KEE, IEE and IEQ variants. We demonstrated that the weak heparin-binder IEE was highly lethal in mice. The initially strong heparin-binding IEQ variant acquired an additional mutation VP1-K244E, which confers weak heparin-binding phenotype resulting in elevated viremia and increased virus antigens in mice brain, with subsequent high virulence. IEE and IEQ-244E variants inoculated into mice disseminated efficiently and displayed high viremia. Increasing polymerase fidelity and impairing recombination of IEQ attenuated the virulence, suggesting the importance of population diversity in EV-A71 pathogenesis in vivo. Combining in silico docking and deep sequencing approaches, we inferred that virus population diversity is shaped by electrostatic interactions at the five-fold axis of the virus surface. Electrostatic surface charges facilitate virus adaptation by generating poor heparin-binding variants for better in vivo dissemination in mice, likely due to reduced adsorption to heparin-rich peripheral tissues, which ultimately results in increased neurovirulence. The dynamic switching between heparin-binding and weak heparin-binding phenotype in vivo explained the neurovirulence of EV-A71., Author summary Enterovirus A71 (EV-A71) is the primary cause of hand, foot and mouth disease, and it can also infect the central nervous system and cause fatal outbreaks in young children. EV-A71 pathogenesis remains elusive. In this study, we demonstrated that EV-A71 variants with strong affinity to heparan sulfate (heparin) have a growth advantage in cell culture, but are disadvantaged in vivo. When inoculated into one-day old mice, strong heparin-binding virus variants are more likely to be adsorbed to peripheral tissues, resulting in impaired ability to disseminate, and are cleared from the bloodstream rapidly. The lower viremia level resulted in no neuroinvasion. In contrast, weak heparin-binding variants show greater levels of viremia, dissemination and subsequent neurovirulence in mice. We also provide evidence that the EV-A71 heparin-binding pattern is mediated by electrostatic surface charges on the virus capsid surface. In mice, EV-A71 undergoes adaptive mutation to acquire greater negative surface charges, thus generating new virulent variants with weak heparin-binding ability which allows greater viral spread. Our study underlines the importance of electrostatic surface charges in shaping EV-A71 virulence.

In this study, we have demonstrated an automated workflow by using KNIME Analytical Platform for modelling and predicting potential HIV-1 protease (HIVP) inhibitors. The workflow has been simplified in three easy steps i.e., 1) retrievethe database of inhibitors for the target disease from ChEMBL website and well-known drug from DrugBank database, 2) generate the descriptors and, 3) select the optimal number of features after machine learning models training. Our results have indicated that the random forest with auto prediction validation method is the most reliable with the best R2 value of 0.9394. Apparently, this workflow can be transformed easily for any other diseases and the quantitative structure-activity relationship (QSAR) model that has been developed can accurately predict in silico how chemical modifications might influence biological behaviour. Overall, the automated workflow which has been presented in this study may significantly reduce the time, cost and efforts needed to design or develop potential HIVP inhibitors.

Enterovirus A71 (EV-A71) causes hand, foot and mouth disease epidemics with neurological complications and fatalities. However, the neuropathogenesis of EV-A71 remains poorly understood. In mice, adaptation and virulence determinants have been mapped to mutations at VP1-145, VP1-244 and VP2-149. We hypothesized that heparin-binding phenotype shapes EV-A71 virulence in mice. We constructed six viruses with varying residues at VP1-98, VP1-145 (which are both heparin-binding determinants) and VP2-149 (based on the wild type 98E/145Q/149K, termed EQK) to generate KQK, KEK, EEK, EEI and EQI variants. We demonstrated that the weak heparin-binder EEI was highly lethal in mice. The initially strong heparin-binding EQI variant acquired an additional mutation VP1-K244E, which confers weak heparin-binding phenotype resulting in elevated viremia and increased brain inflammation and virus antigens in mice, with subsequent high virulence. EEI and EQI-K244E variants inoculated into mice disseminated efficiently and displayed high viremia. Increasing polymerase fidelity and impairing recombination of EQI attenuated virulence, suggesting the importance of population diversity in EV-A71 pathogenesisin vivo. Combiningin silicodocking and deep sequencing approaches, we inferred that virus population diversity is shaped by electrostatic interactions at the five-fold axis of the virus surface. Electrostatic surface charges facilitate virus adaptation by generating poor heparin-binding variants for betterin vivodissemination in mice, likely due to reduced adsorption to heparin-rich peripheral tissues, which ultimately results in increased neurovirulence. The dynamic switching from heparin-binding to weak heparin-binding phenotypein vivoexplained the neurovirulence of EV-A71.Author summaryEnterovirus A71 (EV-A71) is the primary cause of hand, foot and mouth disease, and it can also infect the central nervous system and cause fatal outbreaks in young children. EV-A71 pathogenesis remains elusive. In this study, we demonstrated that EV-A71 variants with strong affinity to heparan sulfate (heparin) have a growth advantage in tissue culture, but are disadvantageousin vivo. When inoculated into mice, strong heparin-binding virus variants are more likely to be adsorbed to peripheral tissues, resulting in impaired ability to disseminate and being cleared from the bloodstream rapidly. The lower viremia level resulted in no neuroinvasion. In contrast, weak heparin-binding variants show greater levels of viremia, dissemination and subsequent neurovirulence in mice. We also provide evidence that the ability of EV-A71 to bind heparin is mediated by electrostatic surface charges due to amino acids on the virus capsid surface. In mice, EV-A71 undergoes adaptive mutation to acquire greater negative surface charges, thus generating new virulent variants with weak heparin-binding which allows greater viral spread. Our study underlines the importance of electrostatic surface charges in shaping EV-A71 virulence.

Naturally occurring proteins are emerging as novel therapeutics in the protein-based biopharmaceutical industry for the treatment of diabetes and obesity. However, proteins are not suitable for oral delivery due to short half-life, reduced physical and chemical stability and low permeability across the membrane. Chemical modification has been identified as a formulation strategy to enhance the stability and bioavailability of protein drugs. The present study aims to study the effect of charge-specific modification of basic amino acids (Lys, Arg) and guanidination on the interaction of insulin with its receptor using molecular modelling. Our investigation revealed that the guanidination of insulin (Lys-NHC = NHNH2) enhanced and exerted stronger binding of the protein to its receptor through electrostatic interaction than native insulin (Lys-NH3+). Point mutations of Lys and Arg (R22, K29; R22K, K29; R22, K29R; R22K, K29R) were attempted and the effects on the interaction and stability between insulin/modified insulins and insulin receptor were also analyzed in this study. The findings from the study are expected to provide a better understanding of the possible mechanism of action of the modified protein at a molecular level before advancing to real experiments.

P-glycoprotein (P-gp) is a 170-kDa membrane protein. It provides a barrier function and help to excrete toxins from the body as a transporter. Some bioflavonoids have been shown to block P-gp activity. To evaluate the important amino acid residues within nucleotide binding domain 1 (NBD1) of P-gp that play a key role in molecular interactions with flavonoids using structure-based pharmacophore model. In the molecular docking with NBD1 models, a putative binding site of flavonoids was proposed and compared with the site for ATP. The binding modes for ligands were achieved using LigandScout to generate the P-gp–flavonoid pharmacophore models. The binding pocket for flavonoids was investigated and found these inhibitors compete with the ATP for binding site in NBD1 including the NBD1 amino acid residues identified by the in silico techniques to be involved in the hydrogen bonding and van der Waals (hydrophobic) interactions with flavonoids. These flavonoids occupy with the same binding site of ATP in NBD1 proffering that they may act as an ATP competitive inhibitor.

Modulation of P-glycoprotein (P-gp)’s function may lead to significant changes in the prescription drugs’ pharmacokinetic profiles and increase potential risks in occurring of drug-drug including herb-drug interactions. This computational structured-based study aimed to screen bioflavonoids which play a role in herb-drug interactions as a P-gp inhibitor utilising molecular docking analysis. 25 flavonoids were used as ligands for docking study. The mouse P-gp (code: 4Q9H) was acquired from the Protein Data Bank (PDB). Docking analysis was operated utilising AutoDock 4.2.6. A lowest estimated free energy of binding value in a cluster of each flavonoid against nucleotide-binding domain 1 (NBD1) of P-gp was analysed. The result points out the high correlation between the estimated free energies of binding and percentage of inhibitory efficiency values of flavonoids applied in the molecular modelling that is powerful and capable to predict potential P-gp interactions among flavonoids and its substrates.

Certain proteins have demonstrated proficient human immunodeficiency virus (HIV-1) life cycle disturbance. Recently, the ankyrin repeat protein targeting the HIV-1 capsid, AnkGAG1D4, showed a negative effect on the viral assembly of the HIV-1NL4-3 laboratory strain. To extend its potential for future clinical application, the activity of AnkGAG1D4 in the inhibition of other HIV-1 circulating strains was evaluated. Chimeric NL4-3 viruses carrying patient-derived Gag/PR-coding regions were generated from 131 antiretroviral drug-na&iuml, ve HIV-1 infected individuals in northern Thailand during 2001&ndash, 2012. SupT1, a stable T-cell line expressing AnkGAG1D4 and ankyrin non-binding control (AnkA32D3), were challenged with these chimeric viruses. The p24CA sequences were analysed and classified using the K-means clustering method. Among all the classes of virus classified using the p24CA sequences, SupT1/AnkGAG1D4 demonstrated significantly lower levels of p24CA than SupT1/AnkA32D3, which was found to correlate with the syncytia formation. This result suggests that AnkGAG1D4 can significantly interfere with the chimeric viruses derived from patients with different sequences of the p24CA domain. It supports the possibility of ankyrin-based therapy as a broad alternative therapeutic molecule for HIV-1 gene therapy in the future.

Zika virus (ZIKV) is an arthropod-borne flavivirus that causes newborn microcephaly and Guillian-Barre syndrome in adults. No therapeutics are available to treat ZIKV infection or other flaviviruses. In this study, we explored the inhibitory effect of glycosaminoglycans and analogues against ZIKV infection. Highly sulfated heparin, dextran sulfate and suramin significantly inhibited ZIKV infection in Vero cells. De-sulfated heparin analogues lose inhibitory effect, implying that sulfonate groups are critical for viral inhibition. Suramin, an FDA-approved anti-parasitic drug, inhibits ZIKV infection with 3-5 log10 PFU viral reduction with IC50 value of ∼2.5-5 μg/ml (1.93 μM-3.85 μM). A time-of-drug-addition study revealed that suramin remains potent even when administrated at 1-24 hpi. Suramin inhibits ZIKV infection by preventing viral adsorption, entry and replication. Molecular dynamics simulation revealed stronger interaction of suramin with ZIKV NS3 helicase than with the envelope protein. Suramin warrants further investigation as a potential antiviral candidate for ZIKV infection. Heparan sulfate (HS) is a cellular attachment receptor for multiple flaviviruses. However, no direct ZIKV-heparin interaction was observed in heparin-binding analysis, and downregulate or removal of cellular HS with sodium chlorate or heparinase I/III did not inhibit ZIKV infection. This indicates that cell surface HS is not utilized by ZIKV as an attachment receptor.

The ankyrin repeat (AR) can be used as a versatile scaffold for protein–protein interactions. It consists of a 33-residues sequence motif found in proteins with diverse functions, such as transcription initiation, cell cycle regulation, cytoskeletal integrity, ion transport, and cell–cell signaling. Using AR with high affinity for the Escherichia coli maltose binding protein (MBP) as our model system, we explored a structure-based computational protocol to probe and characterize binding affinity hot-spots at protein–protein interfaces. In this study, the long time scale dynamics simulations with GPU accelerated molecular dynamics (MD) simulations in AMBER12 have been performed to locate the hot-spots of protein-protein interaction by the analysis of the Molecular Mechanics–Poisson-Boltzmann Surface Area/Generalized Born Solvent Area (MM-PBSA/GBSA) of the MD trajectories. The two designed AR systems with different binding affinities from ELISA were simulated. Our calculations gave the absolute binding affi...

Thiazolidine-4-carboxylic acid derivative compounds with their inhibitory activity against influenza A neuraminidase (NA) were used as a data set for devel- oping the quantitative structure-activity relationship (QSAR) model. The 2D-QSAR model was developed using multiple linear regression analysis with r 2 and r 2 (CV) value of 0.98 and 0.70, respectively. The generated QSAR model has shown that the electrostatic and steric properties have the predominant influence on biological activity. 3D- QSAR was modeled with r 2 value of 0.98 and RMSD of 0.12. The pharmacophore alignments were generated and suggested that the hydrophobic, hydrogen bond donor, and hydrogen bond acceptor features on R 1 and R 2 site sub- stitutions of the core of thiazolidine are important proper- ties to enhance the activities of molecule against influenza A NA.

Manscript received May 13, 2013; revised July 11, 2013. Abstract—In vitro Moringa oleifera Lam. seedlings were successfully produced through seed culture on Murashige & Skoog (MS) agar medium containing 3% (w/v) sucrose and 0.2% (w/v) Gelrite in the absence of growth regulators under 1,500 lux of light density, 16 hour photoperiod light at temperature of 25 ± 2°C. Shoot-derived callus and rootderived callus of M. oleifera were established via culture of shoot and root on the MS medium supplemented with 0.5 mg/l of 2,4-dichlorophenoxy acetic acid (2,4-D) in the dark at 25 ± 2°C. Stem, leaf and root of native M. oleifera and M. oleifera callus were assayed for peroxidase activity using guaiacol as a substrate of the enzyme. In native plant, crude extract from root provided the highest peroxidase specific activity, followed by those from stem and leaf with the specific activity of 19.73 ± 0.18, 16.56 ± 1.43 and 13.38 ± 1.04 unit/mg protein, respectively. Crude extract of root-derived callus and shoot-derived callus of M. oleifera possessed specific activity of 167.25 ± 16.12 and 103.99 ± 10.64 unit/mg protein. These values are significantly higher than their counter parts from native M. oleifera suggesting the potential use of the callus cultures as new and improved sources of peroxidase.

In this work, molecular docking, pharmacophore modeling and molecular dynamics (MD) simulation were rendered for the mouse P-glycoprotein (P-gp) (code: 4Q9H) and bioflavonoids; amorphigenin, chrysin, epigallocatechin, formononetin and rotenone including a positive control; verapamil to identify protein-ligand interaction features including binding affinities, interaction characteristics, hot-spot amino acid residues and complex stabilities. These flavonoids occupied the same binding site with high binding affinities and shared the same key residues for their binding interactions and the binding region of the flavonoids was revealed that overlapped the ATP binding region with hydrophobic and hydrophilic interactions suggesting a competitive inhibition mechanism of the compounds. Root mean square deviations (RMSDs) analysis of MD trajectories of the protein-ligand complexes and NBD2 residues, and ligands pointed out these residues were stable throughout the duration of MD simulations. Thus, the applied preliminary structure-based molecular modeling approach of interactions between NBD2 and flavonoids may be gainful to realize the intimate inhibition mechanism of P-gp at NBD2 level and on the basis of the obtained data, it can be concluded that these bioflavonoids have the potential to cause herb-drug interactions or be used as lead molecules for the inhibition of P-gp (as anti-multidrug resistance agents) via the NBD2 blocking mechanism in future.

BackgroundThe re-emerging,Aedes spp.transmitted Chikungunya virus (CHIKV) has recently caused large outbreaks in a wide geographical distribution of the world including countries in Europe and America. Though fatalities associated with this self-remitting disease were rarely reported, quality of patients’ lives have been severely diminished by polyarthralgia recurrence. Neither effective antiviral treatment nor vaccines are available for CHIKV. Our previous in vitro screening showed that hesperetin, a bioflavonoid exhibits inhibitory effect on the virus intracellular replication. Here, we present a study using the computational approach to identify possible target proteins for future mechanistic studies of hesperetin.Methods3D structures of CHIKV nsP2 (3TRK) and nsP3 (3GPG) were retrieved from Protein Data Bank (PDB), whereas nsP1, nsP4 and cellular factor SPK2 were modeled using Iterative Threading Assembly Refinement (I-TASSER) server based on respective amino acids sequence. We performed molecular docking on hesperetin against all four CHIKV non-structural proteins and SPK2. Proteins preparation and subsequent molecular docking were performed using Discovery Studio 2.5 and AutoDock Vina 1.5.6. The Lipinski’s values of the ligand were computed and compared with the available data from PubChem. Two non-structural proteins with crystal structures 3GPG and 3TRK in complexed with hesperetin, demonstrated favorable free energy of binding from the docking study, were further explored using molecular dynamics (MD) simulations.ResultsWe observed that hesperetin interacts with different types of proteins involving hydrogen bonds, pi-pi effects, pi-cation bonding and pi-sigma interactions with varying binding energies. Among all five tested proteins, our compound has the highest binding affinity with 3GPG at −8.5 kcal/mol. The ligand used in this study also matches the Lipinski’s rule of five in addition to exhibiting closely similar properties with that of in PubChem. The docking simulation was performed to obtain a first guess of the binding structure of hesperetin complex and subsequently analysed by MD simulations to assess the reliability of the docking results. Root mean square deviation (RMSD) of the simulated systems from MD simulations indicated that the hesperetin complex remains stable within the simulation timescale.DiscussionThe ligand’s tendencies of binding to the important proteins for CHIKV replication were consistent with our previous in vitro screening which showed its efficacy in blocking the virus intracellular replication. NsP3 serves as the highest potential target protein for the compound’s inhibitory effect, while it is interesting to highlight the possibility of interrupting CHIKV replication via interaction with host cellular factor. By complying the Lipinski’s rule of five, hesperetin exhibits drug-like properties which projects its potential as a therapeutic option for CHIKV infection.

Enterovirus A71 (EV-A71) is a neurotropic enterovirus that uses heparan sulfate as an attachment receptor. The molecular determinants of EV-A71-heparan sulfate interaction are unknown. With In silico heparin docking and mutagenesis of all possible lysine residues in VP1, we identified that K162, K242 and K244 are responsible for heparin interaction and inhibition. EV-A71 mutants with K242A and K244A rapidly acquired compensatory mutations, T100K or E98A, and Q145R-T237N respectively, which restored the heparin-binding phenotype. Both VP1-98 and VP1-145 modulates heparin binding. Heparin-binding phenotype was completely abolished with VP1-E98-E145, but was restored by an E98K or E145Q substitution. During cell culture adaptation, EV-A71 rapidly acquired K98 or Q/G145 to restore the heparin-binding phenotype. Together with next-generation sequencing analysis, our results implied that EV-A71 has high genetic plasticity by modulating positively-charged residues at the five-fold axis during in vitro heparin adaptation. Our finding has impact on EV-A71 vaccine production, evolutionary studies and pathogenesis.

Aedes aegypti, Aedes albopictus and Culex quinquefasciatus are vectors of dengue fever and West Nile virus diseases. This study was conducted to determine the toxicity, mechanism of action and the binding interaction of three active phenylpropanoids from Piper sarmentosum (Piperaceae) toward late 3rd or early 4th larvae of above vectors. A bioassay guided-fractionation on the hexane extract from the roots of Piper sarmentosum led to the isolation and identification of three active phenylpropanoids; asaricin 1, isoasarone 2 and trans-asarone 3. The current study involved evaluation of the toxicity and acetylcholinesterase (AChE) inhibition of these compounds against Aedes aegypti, Aedes albopictus and Culex quinquefasciatus larvae. Asaricin 1 and isoasarone 2 were highly potent against Aedes aegypti, Aedes albopictus and Culex quinquefasciatus larvae causing up to 100% mortality at ≤ 15 μg/mL concentration. The ovicidal activity of asaricin 1, isoasarone 2 and trans-asarone 3 were evaluated through egg hatching. Asaricin 1 and isoasarone 2 showed potent ovicidal activity. Ovicidal activity for both compounds was up to 95% at 25μg/mL. Asaricin 1 and isoasarone 2 showed strong inhibition on acetylcholinesterase with relative IC50 values of 0.73 to 1.87 μg/mL respectively. These findings coupled with the high AChE inhibition may suggest that asaricin 1 and isoasarone 2 are neuron toxic compounds toward Aedes aegypti, Aedes albopictus and Culex quinquefasciatus. Further computational docking with Autodock Vina elaborates the possible interaction of asaricin 1 and isoasarone 2 with three possible binding sites of AChE which includes catalytic triads (CAS: S238, E367, H480), the peripheral sites (PAS: E72, W271) and anionic binding site (W83). The binding affinity of asaricin 1 and isoasarone 2 were relatively strong with asaricin 1 showed a higher binding affinity in the anionic pocket.

Many proteases are produced as zymogens bearing the N-terminal proregions acting both as intramolecular chaperones and as protease inhibitors. The latter role of the proregions as potent and specific inhibitors of their associated protease has been demonstrated in various peptidases and therefore has been targeted for alternative pest control. Here, we isolated amino acid sequence of Plutella xylostella midgut trypsin from the larvae of diamondback moth and tested in silico for its inhibitory activity toward propeptide models using computational modeling and docking. The propeptide models (AAAPGHR, AAAPGRR, AAAPGKR, AAPGHRI, APGHRIV, PGHRIVG, AAAAPGH, and AAAAAPG) were designed based on histidine-mutated and frame-shifted modifications of the 7-amino-acid proregion (AAAPGHR) of the Plutella xylostella trypsin. Among the eight peptides, AAAPGRR was found to give the best docking scores, showing a strong binding to the cognate enzyme. In addition, the obtained structure of trypsin-AAAPGRR complex was found to share a similar binding mode with a crystal structure of plant protease inhibitor complex. Our results may guide the experiment for the design of future peptide inhibitor with specificity and selectivity for the target enzyme.

Low-pressure plasma has been used to improve the hydrophobicity of Thai silk. In this study, Glycine–Alanine (GA) and Alanine–Glycine (AG) were chosen to represent model compounds of Bombyx mori silk. Single crystals of the simplified model compounds were characterized by polarizing microscopy and X-ray diffraction. The space groups of P212121 and P21 were found for AG and GA, respectively. The initial structures for calculation were obtained from the experimental crystal structures. Density functional theory at the BHandHLYP levels was used to investigate possible mechanisms of fluorine radicals reacting with AG and GA in the SF6 plasma treatment. The results indicate that hydrogen atoms of silk model compounds were most likely to be abstracted from the alanine residue.

The molecular structures of the binding between human immunodeficiency virus-1 protease (HIV-1PR) and various inhibitors including existing extensive natural products extracts have been investigated for anti-HIV drug development. In this study, the binding of HIV-1PR and a terpenoid from Litchi chinensis extracts (3-oxotrirucalla-7,24-dien-21-oic acid) was investigated in order to clarify the inhibition effectiveness of this compound. Molecular dynamics (MD) simulations of HIV-1PR complex with 3-oxotrirucalla-7,24-dien-21-oic acid were performed including water molecules. The MD simulation results indicated the formation of hydrogen bonds between the oxygen atoms of the inhibitor and the catalytic aspartates, which are commonly found in inhibitors–protease complexes. On the other hand, no hydrogen bonding of this particular inhibitor to the flap region was found. In addition, the radial distribution function of water oxygens around the catalytic carboxylate nitrogens of Asp29 and Asp30 suggests that at le...

Since the emergence of the severe acute respiratory syndrome (SARS) to date, neither an effective antiviral drug nor a vaccine against SARS is available. However, it was found that a mixture of two HIV-1 proteinase inhibitors, lopinavir and ritonavir, exhibited some signs of effectiveness against the SARS virus. To understand the fine details of the molecular interactions between these proteinase inhibitors and the SARS virus via complexation, molecular dynamics simulations were carried out for the SARS-CoV 3CL(pro) free enzyme (free SARS) and its complexes with lopinavir (SARS-LPV) and ritonavir (SARS-RTV). The results show that flap closing was clearly observed when the inhibitors bind to the active site of SARS-CoV 3CL(pro). The binding affinities of LPV and RTV to SARS-CoV 3CL(pro) do not show any significant difference. In addition, six hydrogen bonds were detected in the SARS-LPV system, while seven hydrogen bonds were found in SARS-RTV complex.

Bombyx mori silk treated by cold SF 6 plasma was found to show higher hydrophobic property than the original silk. In order to clarify the chemical changes occurring in this treatment, the changes of functional groups of silk surface were investigated by attenuated total reflection (ATR) infrared spectroscopy, surface charge determination, and quantum mechanical calculation. Infrared spectra of original and plasma-treated silks do not show any change in the frequency regions of amide I, II, and III bands which locate at around 1627, 1513, and 1228 cm -1 , respectively. Slight changes were detected for the peak intensities of the bands locating in the frequency region of 1000-1050 cm -1 after plasma treatment. This suggests the formation of CF groups in the Bombyx mori silk chain. The zeta potential experiment suggested that the electrostatic charges of the silk surface were not affected by the plasma treatment. In order to investigate the surface state of the plasma treated silk, the density functional calculation was performed for the model compounds with similar chemical structure as that of Bombyx mori silk. In this calculation, a fluorine radical was located at the various positions of the model compound, and the energetically most plausible structures were extracted to show the chemical reaction of CH + F - → CF + H - .

In order to understand the mechanisms of ligand binding and interaction between two commercial drugs (ligands), zanamivir and oseltamivir and H5N1 Influenza Virus Neuraminidase subtype N1, a three-dimensional model of N1-ligand (GenBank accession no. AAS654617) was initially generated by homology modeling using the 13 high-resolution X-ray structures of neuraminidase N2 and N9 as the template. With the aid of the molecular mechanics and molecular dynamics methods, the final implicit solvent refined model was obtained. It was, then, assessed by PROCHECK, PROSA and VERIFY3D. With this model, a flexible docking study was performed. The results show strong hydrogen bond interactions between the glycerol side chains of zanamivir and Arg29 of the N1. Common hydrogen bonds between the carboxyl groups and Arg279 were found for both drugs. It was also found that the Glu30, Asp62, Arg63, Arg204, Trp310, Tyr313, Glu336, Ile338, Trp348, Ala349 were observed to facilitate the enzyme-ligand non-bonding interactions as ...

The phylogenetic relationship of Coptotermes gestroi (Wasmann) and Coptotermes vastator Light (Isoptera: Rhinotermitidae) was determined using DNA sequence comparisons of mitochondrial genes. Partial sequences of the ribosomal RNA small subunit 12S, ribosomal RNA large subunit 16S, and mitochondrial COII were obtained from nine populations of C. gestroi from South East Asia (Malaysia, Singapore, Thailand, and Indonesia) and four populations of C. vastator from the Philippines and Hawaii. In addition, four populations of Coptotermes formosanus Shiraki and Globitermes sulphureus (Haviland) were used as the outgroups. Consensus sequences were obtained and aligned. C. vastator and C. gestroi are synonymous, based on high sequence homology across the 12S, 16S, and COII genes. The interspecific pairwise sequence divergence, based on Kimura 2-parameter model between C. gestroi and C. vastator, varied only up to 0.80%. Morphometric measurements of 16 characteristics revealed numerous overlaps between the examined individuals of both species. Based on the molecular phylogenetics and morphometric data, it is proposed that C. vastator is a junior synonym of C. gestroi.

We have previously designed a series of antimicrobial peptides (AMPs) and in the current study, the in vivo therapeutic efficacy and toxicity were investigated. Among all the peptides, DM3 conferred protection to a substantial proportion of the lethally infected mice caused by a strain of penicillin-resistant Streptococcus pneumoniae. Synergism was reported and therapeutic efficacy was significantly enhanced when DM3 was formulated in combination with penicillin (PEN). No toxicity was observed in mice receiving these treatments. The in silico molecular docking study results showed that, DM3 has a strong affinity towards three protein targets; autolysin and pneumococcal surface protein A (pspA). Thus AMPs could serve as supporting therapeutics in combination with conventional antibiotics to enhance treatment outcome.

Candida spp. are the most common causes of fungal infections worldwide. Among the Candida species, Candida albicans remains the predominant species that causes invasive candidiasis in most countries. In this study, we used two peptides, KABT-AMP and uperin 3.6 as templates to develop novel antifungal peptides. Their anticandidal activity was assessed using a combination of MIC, time-killing assay and biofilm reduction assay. Hybrid peptides, KU2 and KU3 containing a mixed backbone of KABT-AMP and Uperin 3.6 demonstrated the most potent anticandidal activity with MIC values ranging from 8–16 mg/L. The number of Trp residues and the amphipathic structure of peptides probably enhanced the anticandidal activity of peptides. Increasing the cationicity of the uperin 3.6 analogues resulted in reduced MIC from the range of 64–128 mg/L to 16–64 mg/L and this was also correlated with the antibiofilm activity and killing kinetics of the peptides. Peptides showed synergistic effects when used in combination with conventional antifungals. Peptides demonstrated low haemolytic activity but significant toxicity on two normal human epithelial cell lines. This study provides us with a better understanding on the structure-activity relationship and the balance between cationicity and hydrophobicity of the peptides although the therapeutic application of the peptides is limited.

Breast cancer is the second most common cancer among women worldwide. Several signaling pathways have been implicated as causative and progression agents. The tumor necrosis factor (TNF) α protein plays a dual role in promoting and inhibiting cancer depending largely on the pathway initiated by the binding of the protein to its receptor. Zerumbone, an active constituent of Zingiber zerumbet, Smith, is known to act on the tumor necrosis factor pathway upregulating tumour necrosis factor related apoptosis inducing ligand (TRAIL) death receptors and inducing apoptosis in cancer cells. Zerumbone is a sesquiterpene that is able to penetrate into the hydrophobic pockets of proteins to exert its inhibiting activity with several proteins. We found a good binding with the tumor necrosis factor, kinase κB (IKKβ) and the Nuclear factor κB (NF-κB) component proteins along the TNF pathway. Our results suggest that zerumbone can exert its apoptotic activities by inhibiting the cytoplasmic proteins. It inhibits the IKKβ kinase that activates the NF-κB and also binds to the NF-κB complex in the TNF pathway. Blocking both proteins can lead to inhibition of cell proliferating proteins to be downregulated and possibly ultimate induction of apoptosis.

From the last two decades, computational alanine scanning mutagenesis (cASM) have been successively applied to a variety of protein complexes to study the energetics and structural characteristics of hot spot residues at protein interface. The cASM combines a continuum approach to model solvent interactions with a MM-based approach to atomistically model protein-protein interactions. In the present study, this methodology was used to study the hotspot residues involved in the binding complex of Integrin linked kinase containing Ankyrin repeats with PINCH1 which have profound effect in cell migration, spreading and signalling. Molecular dynamic simulations using a continuum solvent approach (MM-PB/GBSA) were used in one of our previous study and later alanine scanning with post processing protocol was done as a downstream analysis to predict accurately the differences in binding free energies in solution between wild-type and alanine mutated complexes (ΔΔGbind).

Computational approaches have been used by the molecular biologists all around the world as a vital tool in developing and improving the functional and binding properties of proteins, particularly antibody-antigen (Ab-Ag) complexes. Based on a previous work, it has been identified that several residues located in the complementary determining region (CDR) of antibody fragment of FAB23.5 may serve as a potential point mutation modification in order to improve the binding interaction and the inhibitory activity between the single chain variable fragment (scFv) antibody and the epitope variant of HIV-1 p24 capsid protein. In this study, several binding residues from the variable chain domain of scFv anti-p24 were modified. The designed antibody along with the bound antigen was then subjected for molecular dynamics (MD) simulation by using PMEMD.CUDA under Amber 12 with NVIDIA ® Kepler ® GPUs system. GPU system has been fully utilized in this work as a mean to speed up the calculation time. From the simulated model, the newly possessed binding interaction of the modified AbAg complex were acquired by using the Molecular Mechanics/Poisson-Boltzmann surface area (MM-PBSA) and Molecular Mechanics/Generalized Born surface area (MM-GBSA), respectively.

Antimicrobial peptides (AMPs) represent promising alternatives to conventional antibiotics in order to defeat multidrug-resistant bacteria such as Streptococcus pneumoniae. In this study, thirteen antimicrobial peptides were designed based on two natural peptides indolicidin and ranalexin. Our results revealed that four hybrid peptides RN7-IN10, RN7-IN9, RN7-IN8, and RN7-IN6 possess potent antibacterial activity against 30 pneumococcal clinical isolates (MIC 7.81-15.62µg/ml). These four hybrid peptides also showed broad spectrum antibacterial activity (7.81µg/ml) against S. aureus, methicillin resistant S. aureus (MRSA), and E. coli. Furthermore, the time killing assay results showed that the hybrid peptides were able to eliminate S. pneumoniae within less than one hour which is faster than the standard drugs erythromycin and ceftriaxone. The cytotoxic effects of peptides were tested against human erythrocytes, WRL-68 normal liver cell line, and NL-20 normal lung cell line. The results revealed that none of the thirteen peptides have cytotoxic or hemolytic effects at their MIC values. The in silico molecular docking study was carried out to investigate the binding properties of peptides with three pneumococcal virulent targets by Autodock Vina. RN7IN6 showed a strong affinity to target proteins; autolysin, pneumolysin, and pneumococcal surface protein A (PspA) based on rigid docking studies. Our results suggest that the hybrid peptides could be suitable candidates for antibacterial drug development.

Computational approach was employed to evaluate the biological activity of novel cyclooxygenase-2 COX-2 inhibitor, O-desmethylquinine, in comparison to quinine as common inhibitor which can also be used an agent of antipyretic, antimalaria, analgesic and antiinflamation. The molecular models of the compound were constructed and optimized with the density function theory with at the B3LYP/6-31G (d,p) level using Gaussian 09 program. Molecular docking studies of the compounds were done to obtain the COX-2 complex structures and their binding energies were analyzed using the AutoDock Vina. The results of docking of the two ligands were comparable and cannot be differentiated from the energy scoring function with AutoDock Vina.

The ankyrin repeat (AR) protein can be used as a versatile scaffold for protein–protein interactions. It has been found that the heterotrimeric complex between integrin-linked kinase (ILK), PINCH, and parvin is an essential signaling platform, serving as a convergence point for integrin and growth-factor signaling and regulating cell adhesion, spreading, and migration. Using ILK-AR with high affinity for the PINCH1 as our model system, we explored a structure-based computational protocol to probe and characterize binding affinity hot spots at protein–protein interfaces. In this study, the long time scale dynamics simulations with GPU accelerated molecular dynamics (MD) simulations in AMBER12 have been performed to locate the hot spots of protein-protein interaction by the analysis of the Molecular Mechanics–Poisson-Boltzmann Surface Area/Generalized Born Solvent Area (MM-PBSA/GBSA) of the MD trajectories. Our calculations suggest good binding affinity of the complex and also the residues critical in the binding.

International audience; Ankyrins are cellular repeat proteins, which can be genetically modified to randomize amino-acid residues located at defined positions in each repeat unit, and thus create a potential binding surface adaptable to macromolecular ligands. From a phage-display library of artificial ankyrins, we have isolated Ank(GAG)1D4, a trimodular ankyrin which binds to the HIV-1 capsid protein N-terminal domain (NTD(CA)) and has an antiviral effect at the late steps of the virus life cycle. In this study, the determinants of the Ank(GAG)1D4-NTD(CA) interaction were analyzed using peptide scanning in competition ELISA, capsid mutagenesis, ankyrin crystallography and molecular modeling. We determined the Ank(GAG)1D4 structure at 2.2 A resolution, and used the crystal structure in molecular docking with a homology model of HIV-1 capsid. Our results indicated that NTD(CA) alpha-helices H1 and H7 could mediate the formation of the capsid-Ank(GAG)1D4 binary complex, but the interaction involving H7 was predicted to be more stable than with H1. Arginine-18 (R18) in H1, and R132 and R143 in H7 were found to be the key players of the Ank(GAG)1D4-NTD(CA) interaction. This was confirmed by R-to-A mutagenesis of NTD(CA), and by sequence analysis of trimodular ankyrins negative for capsid binding. In Ank(GAG)1D4, major interactors common to H1 and H7 were found to be S45, Y56, R89, K122 and K123. Collectively, our ankyrin-capsid binding analysis implied a significant degree of flexibility within the NTD(CA) domain of the HIV-1 capsid protein, and provided some clues for the design of new antivirals targeting the capsid protein and viral assembly.

Peroxidases (POD) and polyphenol oxidase (PPO) are enzymes that are well known to be involved in the enzymatic browning reaction of fruits and vegetables with different catalytic mechanisms. Both enzymes have some common substrates, but each also has its specific substrates. In our computational study, the amino acid sequence of grape peroxidase (ABX) was used for the construction of models employing homology modeling method based on the X-ray structure of cytosolic ascorbate peroxidase from pea (PDB ID:1APX), whereas the model of grape polyphenol oxidase was obtained directly from the available X-ray structure (PDB ID:2P3X). Molecular docking of common substrates of these two enzymes was subsequently studied. It was found that epicatechin and catechin exhibited high affinity with both enzymes, even though POD and PPO have different binding pockets regarding the size and the key amino acids involved in binding. Predicted binding modes of substrates with both enzymes were also compared. The calculated docking interaction energy of trihydroxybenzoic acid related compounds shows high affinity, suggesting specificity and potential use as common inhibitor to grape ascorbate peroxidase and polyphenol oxidase.

The reaction mechanism of creatinine-creatininase binding to form creatine as a final product has been investigated by using a combined ab initio quantum mechanical/molecular mechanical approach and classical molecular dynamics (MD) simulations. In MD simulations, an X-ray crystal structure of the creatininase/creatinine was modified for creatininase/creatinine complexes and the MD simulations were run for free creatininase and creatinine in water. MD results reveal that two X-ray water molecules can be retained in the active site as catalytic water. The binding free energy from Molecular Mechanics Poisson-Boltzmann Surface Area calculation predicted the strong binding of creatinine with Zn2+, Asp45 and Glu183. Two step mechanisms via Mn2+/Zn2+ (as in X-ray structure) and Zn2+/Zn2+ were proposed for water adding step and ring opening step with two catalytic waters. The pathway using synchronous transit methods with local density approximations with PWC functional for the fragment in the active region were obtained. Preferable pathway Zn2+/Zn2+ was observed due to lower activation energy in water adding step. The calculated energy in the second step for both systems were comparable with the barrier of 26.03 and 24.44 kcal/mol for Mn2+/Zn2+ and Zn2+/Zn2+, respectively.

As it is known that the understanding of the basic properties of the enzyme/inhibitor complex leads directly to enhancing the capability in drug designing and drug discovery. Molecular dynamics simulations have been performed to examine detailed information on the structure and dynamical properties of the HIV-1 PR complexed with saquinavir in the three protonated states, monoprotonates at Asp25 (Mono-25) and Asp25'(Mono-25') and diprotonate (Di-Pro) at both Asp25 and Asp25'. The obtained results support clinical data which reveal that Ile84 and Gly48 are two of the most frequent residues where mutation toward a protease inhibitor takes place. In contrast to the Ile84 mutation due to high displacement of Ile84 in the presence of saquinavir, source of the Gly48 mutation was observed to be due to the limited space in the HIV-1 PR pocket. The Gly48 was, on one side, found to form strong hydrogen bonds with saquinavir, while on the other side this residue was repelled by the hydrophobic Phe53 residue. In terms of inhibitor/enzyme binding, interactions between saquinavir and a catalytic triad of the HIV-1 PR were calculated using the ab initio method. The results show an order of the binding energy of Mono25Di-proMono-25', suggesting that the active site in the HIV-1 PR complexed with saquinavir is monoprotonated states on Asp25. In contrast to the binding energy, 3, 6 and 12 hydrogen bonds between saquinavir and HIV-1 PR were found for the Mono-25, Mono-25' and Di-pro states, respectively. Discrepancy between the two trends suggests us to conclude that interaction between inhibitor and catalytic residues should be used as a criteria to enhance capability in drug designing and drug screening instead of using the total inhibitor/enzyme interaction which is normally reported in the literature. In addition, the distribution and binding of water molecules, in terms of hydrogen bonding, to the donor atoms of saquinavir were investigated and discussed, referring to that which was reported experimentally.

Attempts have been made to predict the binding structures of the human immunodeficiency virus-1 protease (HIV-1Pr) with various inhibitors within the shortest simulation time consuming. The purpose here is to improve the structural prediction by using statistical approach. We use a combination of molecular docking and non-parametric binomial distribution test considering the combination of binding energy, hydrogen bonding, and hydrophobic-hydrophilic interaction in term of binding residues to select the most probable binding structure. In this study, the binding of HIV-1Pr and two inhibitors: Saquinavir and Litchi chinensis extracts (3-oxotrirucalla-7, 24-dien-21-oic acid) were investigated. Each inhibitor was positioned in the active site of HIV-1Pr in many different ways using Lamarckian genetic algorithm and then score each orientation by applying a reasonable evaluation function by AutoDock3.0 program. The results from search methods were screened out using non-parametric binomial distribution test and compared with the binding structure from explicit molecular dynamic simulation. Both complexes from statistical selected docking simulation were found to be comparable with those from X-ray diffraction analysis and explicit molecular dynamic simulation structures.