Search

Leukemia is one of the ten types of cancer that causes the biggest death in the world. Compared to other types of cancer, leukemia has a low life expectancy, so an early diagnosis of the cancer is necessary. A new strategy has been developed to identify various leukemia biomarkers by making blood cancer biosensors, especially by developing nanomaterial applications so that they can improve the performance of the biosensor. Although many biosensors have been developed, the detection of leukemia by using nanomaterials with electrochemical and optical methods is still less carried out compare to other types of cancer biosensors. Even the acoustic and calorimetric testing methods for the detection of leukemia by utilizing nanomaterials have not yet been carried out. Most of the reviewed works reported the use of gold nanoparticles and electrochemical characterization methods for leukemia detection with the object of study being conventional cancer cells. In order to be used clinically by the community, future research must be carried out with a lot of patient blood objects, develop non-invasive leukemia detection, and be able to detect all types of blood cancer specifically with one biosensor. This can lead to a fast and accurate diagnosis thus allowing for early treatment and easy periodic condition monitoring for various types of leukemia based on its biomarker and future design controlable via internet of things (IoT) so that why would be monitoring real times.

Abstract New coumarin derivatives were designed using a 2-(2-oxo-2H-chromen-4-yl)acetic acid scaffold conjugated with amino acid esters or tyramine. The anti-tyrosinase and anti-lipid peroxidation activities of the synthesized compounds were investigated. Coumarin derivatives 7,9, 11–13, 15–18 showed strong anti-lipid peroxidation activity. Compound 13 exhibited uncompetitive tyrosinase inhibitory activity with an IC50 value of 68.86 µM. Compound 14 (% activity = 123.41) showed stronger tyrosinase activating activity than 8-methoxypsolaren (8-MOP, % activity = 109.46). In silico studies revealed different poses between the inhibitors and activators near the tyrosinase catalytic site. Compounds 13 (25–50 μM) and 14 (25–100 μM) did not show cytotoxicity against B16F10 cells. In contrast to the tyrosinase inhibition assay, compound 13 (50 μM) suppressed melanogenesis in B16F10 cells with two times higher potency than KA (100 μM). Compound 14 at 100 μM showed melanogenesis enhancement in B16F10 cells in a dose-dependent manner, however, inferior to the 8-MOP. Based on the findings, compound 13 and 14 offer potential for development as skin-lightening agents and vitiligo therapy agents, respectively.

Abstract The prevalence of HIV-1 infection continues to pose a significant global public health issue, highlighting the need for antiretroviral drugs that target viral proteins to reduce viral replication. One such target is HIV-1 protease (PR), responsible for cleaving viral polyproteins, leading to the maturation of viral proteins. While darunavir (DRV) is a potent HIV-1 PR inhibitor, drug resistance can arise due to mutations in HIV-1 PR. To address this issue, we developed a novel approach using the fragment molecular orbital (FMO) method and structure-based drug design to create DRV analogs. Using combinatorial programming, we generated novel analogs freely accessible via an on-the-cloud mode implemented in Google Colab, Combined Analog generator Tool (CAT). The designed analogs underwent cascade screening through molecular docking with HIV-1 PR wild-type and major mutations at the active site. Molecular dynamics (MD) simulations confirmed the assess ligand binding and susceptibility of screened designed analogs. Our findings indicate that the three designed analogs guided by FMO, 19–0–14–3, 19–8–10–0, and 19–8–14–3, are superior to DRV and have the potential to serve as efficient PR inhibitors. These findings demonstrate the effectiveness of our approach and its potential to be used in further studies for developing new antiretroviral drugs.

Abstract The Musashi (MSI) family of RNA-binding proteins, comprising the two homologs Musashi-1 (MSI1) and Musashi-2 (MSI2), typically regulates translation and is involved in cell proliferation and tumorigenesis. MSI proteins contain two ribonucleoprotein-like RNA-binding domains, RBD1 and RBD2, that bind single-stranded RNA motifs with a central UAG trinucleotide with high affinity and specificity. The finding that MSI also promotes the replication of Zika virus, a neurotropic Flavivirus, has triggered further investigations of the biochemical principles behind MSI–RNA interactions. However, a detailed molecular understanding of the specificity of MSI RBD1/2 interaction with RNA is still missing. Here, we performed computational studies of MSI1–RNA association complexes, investigating different RNA pentamer motifs using molecular dynamics simulations with binding free energy calculations based on the solvated interaction energy method. Simulations with Alphafold2 suggest that predicted MSI protein structures are highly similar to experimentally determined structures. The binding free energies show that two out of four RNA pentamers exhibit a considerably higher binding affinity to MSI1 RBD1 and RBD2, respectively. The obtained structural information on MSI1 RBD1 and RBD2 will be useful for a detailed functional and mechanistic understanding of this type of RNA–protein interactions.

Abstract Anti-interferon gamma autoantibodies (anti-IFN-γ autoAbs) neutralize the IFN-γ-mediated functions, contributing to immunodeficiency. A particular autoAb in patient serum had been previously demonstrated to recognize the same determinant on IFN-γ as the neutralizing anti-IFN-γ monoclonal antibody clone B27 (B27 mAb). This study explored the epitope recognized by B27 mAb. The specific peptide sequence recognized by B27 mAb, TDFLRMMLQEER, was retrieved from a phage display random peptide library. Sequence alignment and homology modeling demonstrated that the queried phage peptide sequence and structure were similar to amino acids at position 27–40 (TLFLGILKNWKEES) of the human IFN-γ. This determinant resides in the contact surface of IFN-γ and interferon gamma receptor 1. To elucidate the crucial amino acids, mutations were introduced by substituting T27 and T27F29L30 with alanine or deleting the amino acid residues T27–L33. The binding of B27 mAb to IFN-γ T27A using western blotting was lesser than that to wild-type. The interaction with triple mutant and T27–L33 deletion mutant using western blotting and sandwich ELISA was abolished. The finding demonstrated that T27, F29, and L30 are critical residues in the B27 antigenic determinant. Identification of the functional domain of IFN-γ decrypted the relevance of neutralizing autoAb in adult-onset immunodeficiency.

Abstract Cholangiocarcinoma (CCA) is a lethal cancer with rapid progression and poor survival. Novel and more effective therapies than those currently available are, therefore, urgently needed. Our research group previously reported the combination of gemcitabine and cytotoxic T lymphocytes to be more effective than single-agent treatment for the elimination of CCA cells. However, gemcitabine treatment of CCA cells upregulates the expression of an immune checkpoint protein (programmed death-ligand 1 [PD-L1]) that consequently inhibits the cytotoxicity of T lymphocytes. To overcome this challenge and take advantage of PD-L1 upregulation upon gemcitabine treatment, we generated recombinant PD-L1xCD3 bispecific T cell engagers (BiTEs) to simultaneously block PD-1/PD-L1 signaling and recruit T lymphocytes to eliminate CCA cells. Two recombinant PD-L1xCD3 BiTEs (mBiTE and sBiTE contain anti-PD-L1 scFv region from atezolizumab and from a published sequence, respectively) were able to specifically bind to both CD3 on T lymphocytes, and to PD-L1 overexpressed after gemcitabine treatment on CCA (KKU213A, KKU055, and KKU100) cells. mBiTE and sBiTE significantly enhanced T lymphocyte cytotoxicity against CCA cells, especially after gemcitabine treatment, and their magnitudes of cytotoxicity were positively associated with the levels of PD-L1 expression. Our findings suggest combination gemcitabine and PD-L1xCD3 BiTE as a potential alternative therapy for CCA.

SARS-CoV-2 virus continues to evolve and mutate causing most of the mutated variants resist to many of the therapeutic monoclonal antibodies (mAbs). Despite several mAbs retained neutralizing capability for Omicron BA.1 and BA.2, reduction in neutralization potency was reported. Hence, effort of searching for mAb that is broader in neutralization breadth without losing the neutralizing ability is continued. MW06 was reported with capability in neutralizing most of the variants of concern (VOC) and it binds to the conserved region (left flank) near epitope mAb sotrovimab (S309). In this study, binding affinity of mAb MW06 and its cocktail formulation with MW05 for receptor binding domain (RBD) SARS-CoV-2 virus was investigated under molecular dynamics simulations (MDs). Binding free energies computed by Molecular Mechanics Generalised Born Surface Area (MM-GBSA) algorithm predicted the binding affinity of MW06 for RBD BA.1 (−53 kcal/mol) as strong as RBD wildtype (−58 kcal/mol) while deterioration was observed for RBD BA.2 (−43 kcal/mol). Alike S309 and MW06, simulated cocktail mAb (MW05 and MW06)-RBD interactions suggested the neutralizing capability of the cocktail formulation for RBD BA.1 and BA.2 reduced. Meanwhile, residue pairs that favour the communication between the mAb and RBD have been identified by decomposing the free energy per pairwise residue basis. Apart from understanding the effects of mutation occurred in the RBD region on human angiotensin-converting enzyme 2 (hACE2) binding, impact of heavily mutated RBD on mAb-RBD interactions was investigated in this study as well. In addition to energetic profile obtained from MDs, plotting the dynamics cross-correlation map of the mAb-RBD complex under elastic network model (ENM) was aimed to understand the cross-correlations between residue fluctuations. It allows simple and rapid analysis on the motions or dynamics of the protein residues of mAbs and RBD in complex. Protein residues having correlated motions are normally part of the structural domains of the protein and their respective motions and protein function are related. Motion of mutated RBD residues and mAb residues was less correlated while their respective interactions energy computed to be higher. The combined techniques of MDs and ENM offered simplicity in understanding dynamics and energy contribution that explain binding affinity of mAb-RBD complexes.

Abstract Owing to their ability to maintain a thermodynamically stable fold at extremely high temperatures, thermophilic proteins (TTPs) play a critical role in basic research and a variety of applications in the food industry. As a result, the development of computation models for rapidly and accurately identifying novel TTPs from a large number of uncharacterized protein sequences is desirable. In spite of existing computational models that have already been developed for characterizing thermophilic proteins, their performance and interpretability remain unsatisfactory. We present a novel sequence-based thermophilic protein predictor, termed SCMTPP, for improving model predictability and interpretability. First, an up-to-date and high-quality dataset consisting of 1853 TPPs and 3233 non-TPPs was compiled from published literature. Second, the SCMTPP predictor was created by combining the scoring card method (SCM) with estimated propensity scores of g-gap dipeptides. Benchmarking experiments revealed that SCMTPP had a cross-validation accuracy of 0.883, which was comparable to that of a support vector machine-based predictor (0.906–0.910) and 2–17% higher than that of commonly used machine learning models. Furthermore, SCMTPP outperformed the state-of-the-art approach (ThermoPred) on the independent test dataset, with accuracy and MCC of 0.865 and 0.731, respectively. Finally, the SCMTPP-derived propensity scores were used to elucidate the critical physicochemical properties for protein thermostability enhancement. In terms of interpretability and generalizability, comparative results showed that SCMTPP was effective for identifying and characterizing TPPs. We had implemented the proposed predictor as a user-friendly online web server at http://pmlabstack.pythonanywhere.com/SCMTPP in order to allow easy access to the model. SCMTPP is expected to be a powerful tool for facilitating community-wide efforts to identify TPPs on a large scale and guiding experimental characterization of TPPs.

Hand, foot and mouth disease (HFMD) is a childhood disease caused by enterovirus A71 (EV-A71) and coxsackievirus A16 (CV-A16). Capsid loops are important epitopes for EV-A71 and CV-A16. Seven chimeric EV-A71 (ChiE71) involving VP1 BC (45.5% similarity), DE, EF, GH and HI loops, VP2 EF loop and VP3 GH loop (91.3% similarity) were substituted with corresponding CV-A16 loops. Only ChiE71-1-BC, ChiE71-1-EF, ChiE71-1-GH and ChiE71-3-GH were viable. EV-A71 and CV-A16 antiserum neutralized ChiE71-1-BC and ChiE71-1-EF. Mice immunized with inactivated ChiE71 elicited high IgG, IFN-γ, IL-2, IL-4 and IL-10. Neonatal mice receiving passive transfer of WT EV-A71, ChiE71-1-EF and ChiE71-1-BC immune sera had 100%, 80.0% and no survival, respectively, against lethal challenges with EV-A71, suggesting that the substituted CV-A16 loops disrupted EV-A71 immunogenicity. Passive transfer of CV-A16, ChiE71-1-EF and ChiE71-1-BC immune sera provided 40.0%, 20.0% and 42.9% survival, respectively, against CV-A16. One-day-old neonatal mice immunized with WT EV-A71, ChiE71-1-BC, ChiE71-1-EF and CV-A16 achieved 62.5%, 60.0%, 57.1%, and no survival, respectively, after the EV-A71 challenge. Active immunization using CV-A16 provided full protection while WT EV-A71, ChiE71-1-BC and ChiE71-1-EF immunization showed partial cross-protection in CV-A16 lethal challenge with survival rates of 50.0%, 20.0% and 40%, respectively. Disruption of a capsid loop could affect virus immunogenicity, and future vaccine design should include conservation of the enterovirus capsid loops.

Leukemia, a condition characterized by the abnormal proliferation of blood cells, poses significant challenges in cancer treatment. Thymoquinone (TQ), a bioactive compound derived from black seed, has demonstrated anticancer properties, including telomerase inhibition and the induction of apoptosis. However, TQ’s poor solubility and limited bioavailability hinder its clinical application. This study explored the use of Sulfobutylether-β-cyclodextrin (SBE-β-CD), a cyclodextrin derivative, to enhance the solubility and stability of TQ for leukemia treatment. SBE-β-CD offers low hemolytic activity and has been successfully employed in controlled drug release systems. The study investigated the formation of inclusion complexes between TQ and SBE-β-CD and evaluated their effects on leukemia cell growth and telomerase activity. The results indicated that the TQ/SBE-β-CD complex exhibited improved solubility and enhanced cytotoxic effects against K-562 leukemia cells compared to TQ alone, suggesting the potential of SBE-β-CD as a drug delivery system for TQ. The annexin V-FITC assay demonstrated increased apoptosis, while the qPCR quantification assay revealed reduced telomerase activity in leukemia cells treated with TQ/SBE-β-CD, supporting its anti-leukemic potential. The molecular docking analysis indicated a strong binding affinity between TQ and telomerase. However, further research is needed to optimize the apoptotic effects and minimize necrosis induction. In conclusion, TQ/SBE-β-CD shows promise as a novel strategy for leukemia treatment by inhibiting telomerase and enhancing the cytotoxic effects of TQ, offering a potential solution to overcome the limitations of TQ’s poor solubility and bioavailability.

Abstract As anticancer peptides (ACPs) have attracted great interest for cancer treatment, several approaches based on machine learning have been proposed for ACP identification. Although existing methods have afforded high prediction accuracies, however such models are using a large number of descriptors together with complex ensemble approaches that consequently leads to low interpretability and thus poses a challenge for biologists and biochemists. Therefore, it is desirable to develop a simple, interpretable and efficient predictor for accurate ACP identification as well as providing the means for the rational design of new anticancer peptides with promising potential for clinical application. Herein, we propose a novel flexible scoring card method (FSCM) making use of propensity scores of local and global sequential information for the development of a sequence-based ACP predictor (named iACP-FSCM) for improving the prediction accuracy and model interpretability. To the best of our knowledge, iACP-FSCM represents the first sequence-based ACP predictor for rationalizing an in-depth understanding into the molecular basis for the enhancement of anticancer activities of peptides via the use of FSCM-derived propensity scores. The independent testing results showed that the iACP-FSCM provided accuracies of 0.825 and 0.910 as evaluated on the main and alternative datasets, respectively. Results from comparative benchmarking demonstrated that iACP-FSCM could outperform seven other existing ACP predictors with marked improvements of 7% and 17% for accuracy and MCC, respectively, on the main dataset. Furthermore, the iACP-FSCM (0.910) achieved very comparable results to that of the state-of-the-art ensemble model AntiCP2.0 (0.920) as evaluated on the alternative dataset. Comparative results demonstrated that iACP-FSCM was the most suitable choice for ACP identification and characterization considering its simplicity, interpretability and generalizability. It is highly anticipated that the iACP-FSCM may be a robust tool for the rapid screening and identification of promising ACPs for clinical use.

A series of C4-substituted tertiary nitrogen-bearing 2′-hydroxychalcones were designed and synthesised based on a previous mixed type acetylcholinesterase inhibitor. Majority of the 2′-hydroxychalcone analogues displayed a better inhibition against acetylcholinesterase (AChE) than butyrylcholinesterase (BuChE). Among them, compound 4c was identified as the most potent AChE inhibitor (IC50: 3.3 µM) and showed the highest selectivity for AChE over BuChE (ratio >30:1). Molecular docking studies suggested that compound 4c interacts with both the peripheral anionic site (PAS) and catalytic anionic site (CAS) regions of AChE. ADMET analysis confirmed the therapeutic potential of compound 4c based on its blood–brain barrier penetrating. Overall, the results suggest that this 2′-hydroxychalcone deserves further investigation into the therapeutic lead for Alzheimer’s disease (AD).

Phosphodiesterase 1B (PDE1B) and PDE10A are dual-specificity PDEs that hydrolyse both cyclic adenosine monophosphate and cyclic guanosine monophosphate, and are highly expressed in the striatum. Several reports have suggested that PDE10A inhibitors may present a promising approach for the treatment of positive symptoms of schizophrenia, whereas PDE1B inhibitors may present a novel mechanism to modulate cognitive deficits. Previously, we have reported a novel dual inhibitor of PDE1B and PDE10A, compound 2 [(3-fluorophenyl)(2-methyl-2,3-dihydro-4H-benzo[b][1,4]oxazin-4-yl)methanone] which has shown inhibitory activity for human recombinant PDE1B and PDE10A in vitro. In the present study, the safety profile of compound 2 has been evaluated in rats in the acute oral toxicity study, as well as; the antipsychotic-like effects in the rat model of schizophrenia. Compound 2 was tolerated up to 1 g/kg when administered at a single oral dose. Additionally, compound 2 has strongly suppressed ketamine-induced hyperlocomotion, which presented a model for the positive symptoms of schizophrenia. It has also shown an ability to attenuate social isolation induced by chronic administration of ketamine and enhanced recognition memory of rats ​in the novel object recognition test. Altogether, our results suggest that compound 2 represents a promising therapy for the treatment of the three symptomatic domains of schizophrenia.

There is an increasing demand in developing new, effective, and affordable anti-cancer against colon and rectal. In this study, our aim is to identify the potential anthraquinone compounds from the root bark of Morinda citrifolia to be tested in vitro against colorectal cancer cell lines. Eight potential anthraquinone compounds were successfully isolated, purified and tested for both in-silico and in-vitro analyses. Based on the in-silico prediction, two anthraquinones, morindone and rubiadin, exhibit a comparable binding affinity towards multitargets of β-catenin, MDM2-p53 and KRAS. Subsequently, we constructed a 2D interaction analysis based on the above results and it suggests that the predicted anthraquinones from Morinda citrifolia offer an attractive starting point for potential antiproliferative agents against colorectal cancer. In vitro analyses further indicated that morindone and damnacanthal have significant cytotoxicity effect and selectivity activity against colorectal cancer cell lines.

SARS-CoV-2 causes the current global pandemic coronavirus disease 2019. Widely-available effective drugs could be a critical factor in halting the pandemic. The main protease (3CLpro) plays a vital role in viral replication; therefore, it is of great interest to find inhibitors for this enzyme. We applied the combination of virtual screening based on molecular docking derived from the crystal structure of the peptidomimetic inhibitors (N3, 13b, and 11a), and experimental verification revealed FDA-approved drugs that could inhibit the 3CLpro of SARS-CoV-2. Three drugs were selected using the binding energy criteria and subsequently performed the 3CLpro inhibition by enzyme-based assay. In addition, six common drugs were also chosen to study the 3CLpro inhibition. Among these compounds, lapatinib showed high efficiency of 3CLpro inhibition (IC50 value of 35 ± 1 μM and Ki of 23 ± 1 μM). The binding behavior of lapatinib against 3CLpro was elucidated by molecular dynamics simulations. This drug could well bind with 3CLpro residues in the five subsites S1', S1, S2, S3, and S4. Moreover, lapatinib's key chemical pharmacophore features toward SAR-CoV-2 3CLpro shared important HBD and HBA with potent peptidomimetic inhibitors. The rational design of lapatinib was subsequently carried out using the obtained results. Our discovery provides an effective repurposed drug and its newly designed analogs to inhibit SARS-CoV-2 3CLpro.

Computational and experimental approaches were adopted to utilize a chromophore diglycolic functionalized fluorescein derivative as a Ca2+ receptor. Fluorescein diglycolic acid (Fl-DGA, 1) was synthesized and used in multivariate determination of Ca2+ and K+. Full-structure computation shows that the complexes of 1 and Ca2+ have comparable energies regardless of additional interaction with lactone moiety. The initial formation of diglycolic-Ca2+ complex followed by macrocyclization is thermodynamically disfavored. A U-shaped pre-organized 1 allows Ca2+ to interact simultaneously with diglycolic and lactone motifs. Both motifs actively participate in Ca2+ recognition and the eleven methylene units in the undecyl arm provides excellent flexibility for reorganization and optimum interaction. Principal component analysis (PCA) of computational molecular properties reveals a simple method in evaluating motifs for cation recognition. Fragment models support full-structure results that negative charge causes significant structural changes, but do not reproduce the full extent of C-O bond breaking observed in the latter. Experimental optical responses show that 1 is selective towards Ca2+ and discriminates against K+ and Mg2+. PCA of emission intensities affords distinct clusters of 0.01, 0.1 and 1 mM Ca2+ and K+, and suggests applicability of this technique for simultaneous determination of cationic plant macronutrients in precision agriculture and a wide variety of other applications.

Inhibition of phosphodiesterase 4 (PDE4) is a promising therapeutic approach for the treatment of inflammatory pulmonary disorders, i.e. asthma and chronic obstructive pulmonary disease. However, the treatment with non-selective PDE4 inhibitors is associated with side effects such as nausea and vomiting. Among the subtypes of PDE4 inhibited by these inhibitors, PDE4B is expressed in immune, inflammatory and airway smooth muscle cells, whereas, PDE4D is expressed in the area postrema and nucleus of the solitary tract. Thus, PDE4D inhibition is responsible for the emetic response. In this regard, a selective PDE4B inhibitor is expected to be a potential drug candidate for the treatment of inflammatory pulmonary disorders. Therefore, a shared feature pharmacophore model was developed and used as a query for the virtual screening of Maybridge and SPECS databases. A number of filters were applied to ensure only compounds with drug-like properties were selected. Accordingly, nine compounds have been identified as final hits, where HTS04529 showed the highest affinity and selectivity for PDE4B over PDE4D in molecular docking. The docked complexes of HTS04529 with PDE4B and PDE4D were subjected to molecular dynamics simulations for 100ns to assess their binding stability. The results showed that HTS04529 was bound tightly to PDE4B and formed a more stable complex with it than with PDE4D.

Given the fact that lithium aluminum hydride (LiAlH4) can exist in distinct crystalline structures under different conditions, in this study, we aim to theoretically investigate the structural properties and the pressure-induced phase transformations of its 13 closely related crystal structures by means of the density functional theory (DFT). The present study reveals that the phase transformation of LiAlH4 from the most stable form (α-phase) to the second most stable form (β-phase) occurs at approximately 3.3 GPa, corresponding to a volume collapse of ∼14% and a reduction of 22% in the crystal volume. Due to the relatively higher hydrogen weight content, β-LiAlH4 becomes a potentially attractive candidate for solid-state hydrogen storage at moderate pressures. The two most stable forms, i.e., the structures with the (i) P21/c (α-LiAlH4) and (ii) I41/a (β-LiAlH4) space groups, have been selected so that their structural and electronic properties can be discussed in greater detail. Our study also shows that the numerical results are greatly influenced by the choice of the DFT methods used, such as the exchange-correlation functionals and optimization schemes.

Protein-protein interaction plays an essential role in almost all cellular processes and biological functions. Coupling molecular dynamics (MD) simulations and nanoparticle tracking analysis (NTA) assay offered a simple, rapid, and direct approach in monitoring the protein-protein binding process and predicting the binding affinity. Our case study of designed ankyrin repeats proteins (DARPins)—AnkGAG1D4 and the single point mutated AnkGAG1D4-Y56A for HIV-1 capsid protein (CA) were investigated. As reported, AnkGAG1D4 bound with CA for inhibitory activity; however, it lost its inhibitory strength when tyrosine at residue 56 AnkGAG1D4, the most key residue was replaced by alanine (AnkGAG1D4-Y56A). Through NTA, the binding of DARPins and CA was measured by monitoring the increment of the hydrodynamic radius of the AnkGAG1D4-gold conjugated nanoparticles (AnkGAG1D4-GNP) and AnkGAG1D4-Y56A-GNP upon interaction with CA in buffer solution. The size of the AnkGAG1D4-GNP increased when it interacted with CA but not AnkGAG1D4-Y56A-GNP. In addition, a much higher binding free energy (∆GB) of AnkGAG1D4-Y56A (−31 kcal/mol) obtained from MD further suggested affinity for CA completely reduced compared to AnkGAG1D4 (−60 kcal/mol). The possible mechanism of the protein-protein binding was explored in detail by decomposing the binding free energy for crucial residues identification and hydrogen bond analysis.

Extracellular signal-regulated kinases 1 and 2 (ERK1/2) play key roles in promoting cell survival and proliferation through the phosphorylation of various substrates. Remarkable antitumour activity is found in many inhibitors that act upstream of the ERK pathway. However, drug-resistant tumour cells invariably emerge after their use due to the reactivation of ERK1/2 signalling. ERK1/2 inhibitors have shown clinical efficacy as a therapeutic strategy for the treatment of tumours with mitogen-activated protein kinase (MAPK) upstream target mutations. These inhibitors may be used as a possible strategy to overcome acquired resistance to MAPK inhibitors. Here, we report a class of repeat proteins—designed ankyrin repeat protein (DARPin) macromolecules targeting ERK2 as inhibitors. The structural basis of ERK2–DARPin interactions based on molecular dynamics (MD) simulations was studied. The information was then used to predict stabilizing mutations employing a web-based algorithm, MAESTRO. To evaluate whether these design strategies were successfully deployed, we performed all-atom, explicit-solvent molecular dynamics (MD) simulations. Two mutations, Ala → Asp and Ser → Leu, were found to perform better than the original sequence (DARPin E40) based on the associated energy and key residues involved in protein-protein interaction. MD simulations and analysis of the data obtained on these mutations supported our predictions.

Abstract Background Computational analysis of protein-protein interaction provided the crucial information to increase the binding affinity without a change in basic conformation. Several docking programs were used to predict the near-native poses of the protein-protein complex in 10 top-rankings. The universal criteria for discriminating the near-native pose are not available since there are several classes of recognition protein. Currently, the explicit criteria for identifying the near-native pose of ankyrin-protein complexes (APKs) have not been reported yet. Results In this study, we established an ensemble computational model for discriminating the near-native docking pose of APKs named “AnkPlex”. A dataset of APKs was generated from seven X-ray APKs, which consisted of 3 internal domains, using the reliable docking tool ZDOCK. The dataset was composed of 669 and 44,334 near-native and non-near-native poses, respectively, and it was used to generate eleven informative features. Subsequently, a re-scoring rank was generated by AnkPlex using a combination of a decision tree algorithm and logistic regression. AnkPlex achieved superior efficiency with ≥1 near-native complexes in the 10 top-rankings for nine X-ray complexes compared to ZDOCK, which only obtained six X-ray complexes. In addition, feature analysis demonstrated that the van der Waals feature was the dominant near-native pose out of the potential ankyrin-protein docking poses. Conclusion The AnkPlex model achieved a success at predicting near-native docking poses and led to the discovery of informative characteristics that could further improve our understanding of the ankyrin-protein complex. Our computational study could be useful for predicting the near-native poses of binding proteins and desired targets, especially for ankyrin-protein complexes. The AnkPlex web server is freely accessible at http://ankplex.ams.cmu.ac.th .

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.

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.

The data is obtained from exploring the modulatory activities of bioflavonoids on P-glycoprotein function by ligand-based approaches. Multivariate Linear-QSAR models for predicting the induced/inhibitory activities of the flavonoids were created. Molecular descriptors were initially used as independent variables and a dependent variable was expressed as pFAR. The variables were then used in MLR analysis by stepwise regression calculation to build the linear QSAR data. The entire dataset consisted of 23 bioflavonoids was used as a training set. Regarding the obtained MLR QSAR model, R of 0.963, R2=0.927, Radj2=0.900, SEE=0.197, F=33.849 and q2=0.927 were achieved. The true predictabilities of QSAR model were justified by evaluation with the external dataset (Table 4). The pFARs of representative flavonoids were predicted by MLR QSAR modelling. The data showed that internal and external validations may generate the same conclusion.

β-Glucosidases and β-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 β-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 −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 β-glycosidases.

The title compound, C12H11N3O2, is a second monoclinic polymorph (P21, with Z′ = 4) of the previously reported monoclinic (P21/c, with Z′ = 2) form [Akhmad Aznan et al. (2010). Acta Cryst. E66, o2400]. Four independent molecules comprise the asymmetric unit, which have the common features of a syn disposition of the pyridine N atom and the toluene ring, and an intramolecular amine–nitro N—H...O hydrogen bond. The differences between molecules relate to the dihedral angles between the rings which range from 2.92 (19) to 26.24 (19)°. The geometry-optimized structure [B3LYP level of theory and 6–311 g+(d,p) basis set] has the same features except that the entire molecule is planar. In the crystal, the three-dimensional architecture is consolidated by a combination of C—H...O, C—H...π, nitro-N—O...π and π–π interactions [inter-centroid distances = 3.649 (2)–3.916 (2) Å].

Combined computational and experimental strategies for the systematic design of chemical sensor arrays using carbonitrile neutral receptors are presented. Binding energies of acetonitrile, n-pentylcarbonitrile and malononitrile with Ca(II), Mg(II), Be(II) and H+ have been investigated with the B3LYP, G3, CBS-QB3, G4 and MQZVP methods, showing a general trend H+ > Be(II) > Mg(II) > Ca(II). Hydrogen bonding, donor-acceptor and cation-lone pair electron simple models were employed in evaluating the performance of computational methods. Mg(II) is bound to acetonitrile in water by 12.5 kcal/mol, and in the gas phase the receptor is more strongly bound by 33.3 kcal/mol to Mg(II) compared to Ca(II). Interaction of bound cations with carbonitrile reduces the energies of the MOs involved in the proposed σ-p conjugated network. The planar malononitrile-Be(II) complex possibly involves a π-network with a cationic methylene carbon. Fabricated potentiometric chemical sensors show distinct signal patterns that can be exploited in sensor array applications.

Background The 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. Methods 3D 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. Results We 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. Discussion The 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.

Burkholderia pseudomallei, the causative agent of melioidosis, is intrinsically resistant to many conventional antibiotics. Therefore, alternative antimicrobial agents such as antimicrobial peptides (AMPs) are extensively studied to combat this issue. Our study aims to identify and understand the mode of action of the potential AMP(s) that are effective against B. pseudomallei in both planktonic and biofilm state as well as to predict the possible binding targets on using in vitro and in silico approaches. In the in vitro study, 11 AMPs were tested against 100 B. pseudomallei isolates for planktonic cell susceptibility, where LL-37, and PG1, demonstrated 100.0% susceptibility and TP1 demonstrated 83% susceptibility. Since the B. pseudomallei activity was reported on LL-37 and PG1, TP1 was selected for further investigation. TP1 inhibited B. pseudomallei cells at 61.69 μM, and membrane blebbing was observed using scanning electron microscopy. Moreover, TP1 inhibited B. pseudomallei cell growth, reaching bactericidal endpoint within 2 h post exposure as compared to ceftazidime (CAZ) (8 h). Furthermore, TP1 was shown to suppress the growth of B. pseudomallei cells in biofilm state at concentrations above 221 μM. However, TP1 was cytotoxic to the mammalian cell lines tested. In the in silico study, molecular docking revealed that TP1 demonstrated a strong interaction to the common peptide or inhibitor binding targets for lipopolysaccharide of Escherichia coli, as well as autolysin, pneumolysin, and pneumococcal surface protein A (PspA) of Streptococcus pneumoniae. Homology modelled B. pseudomallei PspA protein (YDP) also showed a favourable binding with a strong electrostatic contribution and nine hydrogen bonds. In conclusion, TP1 demonstrated a good potential as an anti-B. pseudomallei agent.

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.

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ïve HIV-1 infected individuals in northern Thailand during 2001⁻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.

A new macromolecule possessing two dansyl moieties and based on 2-[4-(2-aminoethylthio)butylthio]ethanamine was prepared as a fluorescent sensor and its mercury sensing properties toward various transition metal, alkali, and alkali earth ions were investigated. The designed compound exhibited pronounced Hg2+-selective ON-OFF type fluorescence switching upon binding. The new compoundprovided highly selective sensing to Hg2+ in acetonitrile-water solvent mixtures with a detection limit of 2.49 x 10-7 M or 50 ppb. The molecular modeling results indicated that ions-recognition of the sensor originated from a self assembly process of the reagentand Hg2+ to form a helical wrapping structure with the favorable electrostatic interactions of Hg2+coordinated with sulfur, oxygen, nitrogen atoms and aromatic moieties.

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.