Structural and dynamical thermostability of psychrophilic enzyme at various temperatures: Molecular dynamics simulations of tryptophan synthase.

Mesophilic enzymes are among the most frequently used biocatalysts, however, psychrophilic enzymes are crucially important for their use in heat-sensitive reactions. How enzymes can work efficiently at various range of temperatures is an interesting subject for researchers, and yet it is very least explored. The structural and dynamical behavior of psychrophilic enzymes and their thermostability at various temperatures can help to understand the mechanism and function at molecular level, and for this purpose the ligand-free α-subunit of Shewanella frigidimarina's tryptophan synthase (Sf-TRPS) in isolated monomeric and in hetero-αß-dimeric states was subjected to molecular dynamics (MD) simulations study. The simulation sampled a complete open conformation of Loop L6 in α-subunit with and without ß-partner, which was further investigated under three temperatures mimicking psychrophilic, mesophilic and thermophilic environment. The results indicated an imperative role of ß-subunit in the dynamics of L6 loop as well as in the thermostability of α-subunit by increasing interaction strength at the αß-interface. An interesting relation was observed between the numbers of H-bonds and residue-pairs forming salt bridges at every temperature, and the combine effect seemed to regulate the balance between protein rigidity and flexibility. The outcome of the study will help to understand the driving forces that lead to the stability of the protein at different temperature, and thereby, assist in enzyme engineering that will be beneficial from industrial point of view.

Discovery of Klotho peptide antagonists against Wnt3 and Wnt3a target proteins using combination of protein engineering, protein-protein docking, peptide docking and molecular dynamics simulations.

The Klotho is known as lifespan enhancing protein involved in antagonizing the effect of Wnt proteins. Wnt proteins are stem cell regulators, and uninterrupted exposure of Wnt proteins to the cell can cause stem and progenitor cell senescence, which may lead to aging. Keeping in mind the importance of Klotho in Wnt signaling, in silico approaches have been applied to study the important interactions between Klotho and Wnt3 and Wnt3a (wingless-type mouse mammary tumor virus (MMTV) integration site family members 3 and 3a). The main aim of the study is to identify important residues of the Klotho that help in designing peptides which can act as Wnt antagonists. For this aim, a protein engineering study is performed for Klotho, Wnt3 and Wnt3a. During the theoretical analysis of homology models, unexpected role of number of disulfide bonds and secondary structure elements has been witnessed in case of Wnt3 and Wnt3a proteins. Different in silico experiments were carried out to observe the effect of correct number of disulfide bonds on 3D protein models. For this aim, total of 10 molecular dynamics (MD) simulations were carried out for each system. Based on the protein-protein docking simulations of selected protein models of Klotho with Wnt3 and Wnt3a, different peptides derived from Klotho have been designed. Wnt3 and Wnt3a proteins have three important domains: Index finger, N-terminal domain and a patch of ∼10 residues on the solvent exposed surface of palm domain. Protein-peptide docking of designed peptides of Klotho against three important domains of palmitoylated Wnt3 and Wnt3a yields encouraging results and leads better understanding of the Wnt protein inhibition by proposed Klotho peptides. Further in vitro studies can be carried out to verify effects of novel designed peptides as Wnt antagonists.

Synthesis of novel derivatives of oxindole, their urease inhibition and molecular docking studies.

We synthesized a series of novel 5-24 derivatives of oxindole. The synthesis started from 5-chlorooxindole, which was condensed with methyl 4-carboxybezoate and result in the formation of benzolyester derivatives of oxindole which was then treated with hydrazine hydrate. The oxindole benzoylhydrazide was treated with aryl acetophenones and aldehydes to get target compounds 5-24. The synthesized compounds were evaluated for urease inhibition; the compound 5 (IC50 = 13.00 ± 0.35 µM) and 11 (IC50 = 19.20 ± 0.50 µM) showed potent activity as compared to the standard drug thiourea (IC50 = 21.00 ± 0.01 µM). Other compounds showed moderate to weak activity. All synthetic compounds were characterized by different spectroscopic techniques including (1)H NMR, (13)C NMR, IR and EI MS. The molecular interactions of the active compounds within the binding site of urease enzyme were studied through molecular docking simulations.

Natural products can be potential inhibitors of metalloproteinase II from Bacteroides fragilis to intervene colorectal cancer.

Bacteroides fragilis, a gram negative and obligate anaerobe bacterium, is a member of normal gut microbiota and facilitates many essential roles being performed in human body in normal circumstances specifically in Gastrointestinal or GI tract. Sometimes, due to genetics, epigenetics, and environmental factors, Bacteroides fragilis and their protein(s) start interacting with intestinal epithelium thus damaging the lining leading to colorectal cancers (CRC). To identify these protein(s), we incorporated a novel subtractive proteomics approach in the study. Metalloproteinase II (MPII), a Bacteroides fragilis toxin (bft), was investigated for its virulence and unique pathways to demonstrate its specificity and uniqueness in pathogenicity followed by molecular docking against a set of small drug-like natural molecules to discover potential inhibitors against the toxin. All these identified inhibitor-like molecules were analyzed for their ADMET calculations and detailed physiochemical properties to predict their druggability, GI absorption, blood brain barrier and skin permeation, and others. Resultantly, a total of ten compounds with the least binding energies were obtained and were subjected to protein-compound interaction analysis. Interaction analysis revealed the most common ligand-interacting residues in MPII are His 345, Glu 346, His 339, Gly 310, Tyr 341, Pro 340, Asp 187, Phe 309, Lys 307, Ile 185, Thr 308, and Pro 184. Therefore, top three compounds complexed with MPII having best binding energies were selected in order to analyze their trajectories. RMSD, RMSF, Rg and MMPBSA analysis revealed that all compounds showed good binding and keeping the complex stable and compact throughout the simulation time in addition to all properties and qualities of being a potential inhibitor against MPII.

Discovering Potential Bacteriocins Against Pseudomonas fragi: a Subtractive Proteomics and Molecular Dynamic Simulation Study for Food Preservation.

Food preservation is a schematic and scientific procedure employed for the maintenance and improvement of food's quality, shelf life, and nutritional value. Although, on one hand, ancient conventional methods such as freezing, pasteurization, canning, and chemical methods have the potential to lengthen the shelf life of edible substances, but on the other hand, they can deteriorate its nutritional value as well. Present research focuses on the identification of promising bacteriocins against Pseudomonas fragi via subtractive proteomics pipeline as an alternative approach for food preservation. Bacteriocins are small peptides produced by certain microbes to naturally defend themselves by destroying other closely related bacteria residing in their neighborhood. P. fragi lies among the most notable microbes responsible for the elicitation of food spoilage. Due to increasing emergence and prevalence of multidrug resistance bacteria, there is a need to unravel novel drug targets, crucially involved in food decay process. Based on subtractive scrutinization, UDP-N-acetylglucosamine O-acyltransferase (LpxA) was chosen as promising therapeutic protein target that could play a significant role in progression of food spoilage. Subtilosin A, thuricin-CD, and mutacin B-NY266 were found as the most robust inhibitors of LpxA according to the molecular docking assay results. Molecular dynamic simulations and binding energy calculations via MM/PBSA method of LpxA and three top hit docked complexes, i.e., LpxA-subtilosin A, LpxA-thuricin-CD, and LpxA-mutacin B-NY266, revealed stability throughout simulations and ensured that shortlisted bacteriocins had strong affinity for LpxA.

Computational medicinal chemistry applications to target Asian-prevalent strain of hepatitis C virus.

Hepatitis C Virus (HCV), affecting millions of people worldwide, is the leading cause of liver disorder, cirrhosis, and hepatocellular carcinoma. HCV is genetically diverse having eight genotypes and several subtypes predominant in different regions of the globe. The HCV NS3/4A protease is a primary therapeutic target for HCV with various FDA-approved antivirals and several clinical developments. However, available protease inhibitors (PIs) have lower potency against HCV genotype 3 (GT3), prevalent in South Asia. In this study, the incumbent computational tools were utilized to understand and explore interactions of the HCV GT3 receptor with the potential inhibitors after the virtual screening of one million compounds retrieved from the ZINC database. The molecular dynamics, pharmacological studies, and experimental studies uncovered the potential PIs as ZINC000224449889, ZINC000224374291, and ZINC000224374456 and the derivative of ZINC000224374456 from the ZINC library. The study revealed that these top-hit compounds exhibited good binding and better pharmacokinetics properties that might be considered the most promising compound against HCV GT3 protease. Viability test, on primary healthy Human Gingival Fibroblasts (HGFs) and cancerous AGS cell line, was also carried out to assess their safety profile after administration. In addition, Surface Plasmon Resonance (SPR) was also performed for the determination of affinity and kinetics of synthesized compounds with target proteins.

Antiprotozoal Activity of Thymoquinone (2-Isopropyl-5-methyl-1,4-benzoquinone) for the Treatment of Leishmania major-Induced Leishmaniasis: In Silico and In Vitro Studies.

Leishmaniasis, a neglected tropical parasitic disease (NTPD), is caused by various Leishmania species. It transmits through the bites of the sandfly. The parasite is evolving resistance to commonly prescribed antileishmanial drugs; thus, there is an urgent need to discover novel antileishmanial drugs to combat drug-resistant leishmaniasis. Thymoquinone (2-isopropyl-5-methyl-1,4-benzoquinone; TQ), a primary pharmacologically active ingredient of Nigella sativa (black seed) essential oil, has been reported to possess significant antiparasitic activity. Therefore, the present study was designed to investigate the in vitro and in silico antileishmanial activity of TQ against various infectious stages of Leishmania major (L. major), i.e., promastigotes and amastigotes, and its cytotoxicity against mice macrophages. In silico molecular dockings of TQ were also performed with multiple selected target proteins of L. major, and the most preferred antileishmanial drug target protein was subjected to in silico molecular dynamics (MD) simulation. The in vitro antileishmanial activity of TQ revealed that the half-maximal effective concentration (EC50), half-maximal cytotoxic concentration (CC50), and selectivity index (SI) values for promastigotes are 2.62 ± 0.12 µM, 29.54 ± 0.07 µM, and 11.27, while for the amastigotes, they are 17.52 ± 0.15 µM, 29.54 ± 0.07 µM, and 1.69, respectively. The molecular docking studies revealed that squalene monooxygenase is the most preferred antileishmanial drug target protein for TQ, whereas triosephosphate isomerase is the least preferred. The MD simulation revealed that TQ remained stable in the binding pocket throughout the simulation. Additionally, the binding energy calculations using Molecular Mechanics Generalized-Born Surface Area (MMGBSA) indicated that TQ is a moderate binder. Thus, the current study shows that TQ is a promising antileishmanial drug candidate that could be used to treat existing drug-resistant leishmaniasis.

X-ray and NMR crystallography in an enzyme active site: the indoline quinonoid intermediate in tryptophan synthase.

Chemical-level details such as protonation and hybridization state are critical for understanding enzyme mechanism and function. Even at high resolution, these details are difficult to determine by X-ray crystallography alone. The chemical shift in NMR spectroscopy, however, is an extremely sensitive probe of the chemical environment, making solid-state NMR spectroscopy and X-ray crystallography a powerful combination for defining chemically detailed three-dimensional structures. Here we adopted this combined approach to determine the chemically rich crystal structure of the indoline quinonoid intermediate in the pyridoxal-5'-phosphate-dependent enzyme tryptophan synthase under conditions of active catalysis. Models of the active site were developed using a synergistic approach in which the structure of this reactive substrate analogue was optimized using ab initio computational chemistry in the presence of side-chain residues fixed at their crystallographically determined coordinates. Various models of charge and protonation state for the substrate and nearby catalytic residues could be uniquely distinguished by their calculated effects on the chemical shifts measured at specifically (13)C- and (15)N-labeled positions on the substrate. Our model suggests the importance of an equilibrium between tautomeric forms of the substrate, with the protonation state of the major isomer directing the next catalytic step.

The role of oligomerization and cooperative regulation in protein function: the case of tryptophan synthase.

The oligomerization/co-localization of protein complexes and their cooperative regulation in protein function is a key feature in many biological systems. The synergistic regulation in different subunits often enhances the functional properties of the multi-enzyme complex. The present study used molecular dynamics and Brownian dynamics simulations to study the effects of allostery, oligomerization and intermediate channeling on enhancing the protein function of tryptophan synthase (TRPS). TRPS uses a set of α/ß-dimeric units to catalyze the last two steps of L-tryptophan biosynthesis, and the rate is remarkably slower in the isolated monomers. Our work shows that without their binding partner, the isolated monomers are stable and more rigid. The substrates can form fairly stable interactions with the protein in both forms when the protein reaches the final ligand-bound conformations. Our simulations also revealed that the α/ß-dimeric unit stabilizes the substrate-protein conformation in the ligand binding process, which lowers the conformation transition barrier and helps the protein conformations shift from an open/inactive form to a closed/active form. Brownian dynamics simulations with a coarse-grained model illustrate how protein conformations affect substrate channeling. The results highlight the complex roles of protein oligomerization and the fine balance between rigidity and dynamics in protein function.

In Vitro and In Silico Approaches for the Antileishmanial Activity Evaluations of Actinomycins Isolated from Novel Streptomyces smyrnaeus Strain UKAQ_23.

Leishmaniasis, a Neglected Tropical Parasitic Disease (NTPD), is induced by several Leishmania species and is disseminated through sandfly (Lutzomyia longipalpis) bites. The parasite has developed resistance to currently prescribed antileishmanial drugs, and it has become pertinent to the search for new antileishmanial agents. The current study aimed to investigate the in vitro and in silico antileishmanial activity of two newly sourced actinomycins, X2 and D, produced by the novel Streptomyces smyrnaeus strain UKAQ_23. The antileishmanial activity conducted on promastigotes and amastigotes of Leishmania major showed actinomycin X2 having half-maximal effective concentrations (EC50), at 2.10 ± 0.10 µg/mL and 0.10 ± 0.0 µg/mL, and selectivity index (SI) values of 0.048 and 1, respectively, while the actinomycin D exhibited EC50 at 1.90 ± 0.10 µg/mL and 0.10 ± 0.0 µg/mL, and SI values of 0.052 and 1. The molecular docking studies demonstrated squalene synthase as the most favorable antileishmanial target protein for both the actinomycins X2 and D, while the xanthine phosphoribosyltransferase was the least favorable target protein. The molecular dynamics simulations confirmed that both the actinomycins remained stable in the binding pocket during the simulations. Furthermore, the MMPBSA (Molecular Mechanics Poisson-Boltzmann Surface Area) binding energy calculations established that the actinomycin X2 is a better binder than the actinomycin D. In conclusion, both actinomycins X2 and D from Streptomyces smyrnaeus strain UKAQ_23 are promising antileishmanial drug candidates and have strong potential to be used for treating the currently drug-resistant leishmaniasis.

Synergistic regulation and ligand-induced conformational changes of tryptophan synthase.

Conformational changes of enzyme complexes are often related to regulating and creating an optimal environment for efficient chemistry. We investigated the synergistic regulation of the tryptophan synthase (TRPS) complex, studied for decades as a model of allosteric regulation and substrate channeling within protein complexes. TRPS is a bifunctional tetrameric alphabetabetaalpha enzyme complex that exhibits cooperative motions of the alpha- and beta-subunits by tightly controlled allosteric interactions. We have delineated the atomically detailed dynamics and conformational changes of TRPS in the absence and presence of substrates using molecular dynamics simulations. The computed energy and entropy associated with the protein motions also offer mechanistic insights into the conformational fluctuations and the ligand binding mechanism. The flexible alpha-L6 loop samples both open and partially closed conformations in the ligand-free state but shifts to fully closed conformations when its substrates are present. The fully closed conformations are induced by favorable protein-ligand interactions but are partly compensated by configurational entropy loss. Considerable local rearrangements exist during ligand binding processes when the system is searching for energy minima. The motion of the region that closes the beta-subunit during catalysis, the COMM domain, couples with the motion of the alpha-subunit, although the fluctuations are smaller than in the flexible loop regions. Because of multiple conformations of ligand-free TRPS in the open and partially closed states, the alpha-L6 loop fluctuations have preferential directionality, which may facilitate the fully closed conformations induced by alpha- and beta-substrates binding to both subunits. Such cooperative and directional motion may be a general feature that contributes to catalysis in many enzyme complexes.

Reactivation potency of two novel oximes (K456 and K733) against paraoxon-inhibited acetyl and butyrylcholinesterase: In silico and in vitro models.

Organophosphates (OPs) irreversibly inhibit acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) enzymes. The reactivation of these inhibited enzymes is paramount for their normal function. Present study evaluates reactivation potency of two newly developed oximes, K456 and K733, against paraoxon (POX)-inhibited human-RBC-AChE and human-plasma-BChE in comparison to reported reactivator, pralidoxime (2-PAM). In vitro studies showed higher intrinsic toxicities of both oximes than 2-PAM for AChE. No substantial reactivation of hBChE was noted by tested concentration. Contrary to 2-PAM, the in silico study predicted lower binding free energies for both oximes. However, the detailed interaction study revealed inability of oximes to interact with catalytic anionic site of AChE and hBChE in contrast to 2-PAM. Both in vitro and in silico studies conclude that K456 and K733 are unlikely to be used as reactivators of paraoxon-inhibited AChE or BChE.

Exploring structure and dynamics of the diaquotriamminezinc(II) complex by QM/MM MD simulation.

The structural and dynamical properties of the cis-(O-Zn-O angle approximately 90 degrees) and trans-(O-Zn-O angle approximately 180 degrees) isomers of the model diaquotriamminezinc(II) complex in aqueous solution have been evaluated using the hybrid quantum mechanical/molecular mechanical molecular dynamics simulation approach at ab initio Hartree-Fock level. In both complexes, the first hydration shell contains five ligands (two water and three ammonia molecules) arranged in a trigonal bipyramidal geometry. In the metastable cis-isomer two different bond lengths of 2.34 and 2.13 A are observed for the Zn-Oax and Zn-Oeq bonds, respectively. The trans-isomer shows the maximum of the Zn-O distance at 2.26 A. The Zn-N bond distances in both cases are approximately 2.12 A. A geometrical transformation of the cis-isomer into the trans-isomer was observed after 11.5 ps of simulation, and the trans-isomer then remained stable throughout the whole simulation time of 30 ps. A comparative study for both isomers has been performed in terms of radial distribution functions, coordination number distributions, angular distribution functions, tilt and theta angle distributions, ligands' mean residence time, ion-ligand stretching frequencies, and the vibrational and librational motions of water ligands. The results are compared with the data for the previously studied zinc-monoamine and -diamine complexes.

Meta-Analysis of miRNAs and Their Involvement as Biomarkers in Oral Cancers.

Oral Squamous Cell Carcinoma (OSCC) is one of the most common cancers worldwide. Recent studies have highlighted the role of miRNA in disease pathology, indicating its potential use as an early diagnostic marker. Dysregulated expression of miRNAs is known to affect cell growth, and these may function as tumor suppressors or oncogenes in various cancers. The main objective of this study was to characterize the extracellular miRNAs involved in oral cancer (OC) that can potentially be used as biomarkers of OC. A total of 318 miRNAs involved in oral carcinoma were shortlisted. Differentially expressed genes (DEGs) of oral carcinoma from reported experiments were identified. Common genes between lists of DEGs of OC of each miRNA were identified. These common genes are the targets of specific miRNA, which may be used as biomarkers of OC. A list of significant biomarkers for cancer was generated like CDH2 and CDK7, and functional enrichment analysis identified the role of miRNAs in major pathways like cell adhesion molecules pathway affected by cancer. We observed that at least 25 genes like ABCF3, ALDH2, CD163L1, and so forth are regulated by a maximum number of miRNAs; thereby, they can be used as biomarkers of OC.

High prevalence of G3 rotavirus in hospitalized children in Rawalpindi, Pakistan during 2014.

Rotavirus A species (RVA) is the leading cause of severe diarrhea among children in both developed and developing countries. Among different RVA G types, humans are most commonly infected with G1, G2, G3, G4 and G9. During 2003-2004, G3 rotavirus termed as "new variant G3" emerged in Japan that later disseminated to multiple countries across the world. Although G3 rotaviruses are now commonly detected globally, they have been rarely reported from Pakistan. We investigated the genetic diversity of G3 strains responsible RVA gastroenteritis in children hospitalized in Rawalpindi, Pakistan during 2014. G3P[8] (18.3%; n = 24) was detected as the most common genotype causing majority of infections in children less than 06 months. Phylogenetic analysis of Pakistani G3 strains showed high amino acid similarity to "new variant G3" and G3 strains reported from China, Russia, USA, Japan, Belgium and Hungary during 2007-2012. Pakistani G3 strains belonged to lineage 3 within sub-lineage 3d, containing an extra N-linked glycosylation site compared to the G3 strain of RotaTeqTM. To our knowledge, this is the first report on the molecular epidemiology of G3 rotavirus strains from Pakistan and calls for immediate response measures to introduce RV vaccine in the routine immunization program of the country on priority.

Stability of different zinc(II)-diamine complexes in aqueous solution with respect to structure and dynamics: a QM/MM MD study.

In the context of our detailed study of the chemical behavior of aquo- and ammine-Zn(II) complexes, ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations were performed at the Hartree-Fock (HF) level for the zinc(II)-diamine complexes in aqueous solution. The initial structures of cis and trans isomers of the tetraaquodiamminezinc(II) complex were found to transform into the triaquodiamminezinc(II) complex by releasing one water ligand after approximately 6 and approximately 22 ps of simulation time, respectively. The structural and dynamical properties of these three zinc complexes, i.e., cis-[Zn(NH3)2(H2O)4]2+, trans-[Zn(NH3)2(H2O)4]2+, and [Zn(NH3)2(H2O)3]2+, were analyzed in terms of radial distribution functions (RDF), coordination number distributions (CND), angular distribution functions (ADF), tilt and theta angle distributions, ligands' mean residence times (MRTs), and ion-ligand stretching frequencies. One considerably elongated Zn-O bond of 2.43 A was observed in the case of the cis isomer for one of the water ligands located in the trans position to an ammonia ligand. In the trans isomer the average Zn-O bond length was observed to be 2.23 A, while in the triaquodiamminezinc(II) complex two distinct Zn-O bonds, namely 2.12 A for the ligands in the trigonal plane and 2.26 A for axial water molecules, were observed. As both of the octahedral isomers are transformed into the pentacoordinated structure within the picosecond range, they might be regarded as "metastable species or intermediates", while the triaquodiamminezinc(II) complex is the most stable species of the zinc(II)-diamine complex in aqueous solution.

Temperature effects on the structural and dynamical properties of the Zn(II)-water complex in aqueous solution: a QM/MM molecular dynamics study.

An ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulation at double-zeta restricted Hartree-Fock (RHF) level was performed at an elevated temperature of 363 K (90 degrees C) to study the temperature effects on the structural and dynamical properties of a Zn(II)-water complex in aqueous solution. The first hydration shell, consisting of 6 water molecules at a mean Zn-O distance of 2.16 A, was found to remain stable also at 90 degrees C with respect to exchange processes. The flexible second shell contains, in average, approximately 27 water ligands. To fully characterize the hydration structure, several other parameters such as radial and angular distribution functions (RDF and ADF) and tilt- and theta-angle distributions were evaluated and compared to data obtained at 298 K (25 degrees C). Temperature effects on the dynamics of the Zn(II)-water complex were studied in terms of water reorientations, mean ligand residence times (MRTs), and number of ligand exchange processes. To get further insight into the solute dynamics, additional data, in particular, librational and vibrational motions of water ligands and Zn-O stretching frequencies, were calculated. The second shell is considerably influenced by the elevated temperature, as the ligands' mean residence time is shortened to 4 ps from the value of 10.5 observed at room temperature. The values of the QM/MM MD simulation were also compared to the results of a classical molecular dynamics (CMD) simulation with two- plus three-body potential performed at 90 degrees C, revealing that an accurate description of the second shell and the dynamics of the Zn(II) hydrate needs the inclusion of quantum mechanics in the description.

Virtual screening of eighteen million compounds against dengue virus: Combined molecular docking and molecular dynamics simulations study.

Dengue virus is a major issue of tropical and sub-tropical regions. Dengue virus has been the cause behind the major alarming epidemics in the history with mass causalities from the decades. Unavailability of on-shelf drugs for the prevention of further proliferation of virus inside the human body results in immense number of deaths each year. This issue necessitates the design of novel anti-dengue drug. The protease enzyme pathway is the critical target for drug design due to its significance in the replication, survival and other cellular activities of dengue virus. Therefore, approximately eighteen million compounds from the ZINC database have been virtually screened against nonstructural protein 3 (NS3). The incremental construction algorithm of Glide docking program has been used with its features high throughput virtual screening (HTVS), standard precision (SP), extra precision (XP) and in combination of Prime module, induced fit docking (IFD) approach has also been applied. Five top-ranked compounds were then selected from the IFD results with better predicted binding energies with the catalytic triad residues (His51, Asp75, and Ser135) that may act as potential inhibitors for the underlying target protease enzyme. The top-ranked compounds ZINC95518765, ZINC44921800, ZINC71917414, ZINC39500661, ZINC36681949 have shown the predicted binding energies of -7.55, -7.36, -8.04, -8.41, -9.18kcal/mol, respectively, forming binding interactions with three catalytically important amino acids. Top-docking poses of compounds are then used in molecular dynamics (MD) simulations. In computational studies, our proposed compounds confirm promising results against all the four serotypes of dengue virus, strengthening the opportunity of these compounds to work as potential on-shelf drugs against dengue virus. Further experimentation on the proposed compounds can result in development of strong inhibitors.

Computational Docking Study of p7 Ion Channel from HCV Genotype 3 and Genotype 4 and Its Interaction with Natural Compounds.

BACKGROUND: The current standard care therapy for hepatitis C virus (HCV) infection consists of two regimes, namely interferon-based and interferon-free treatments. The treatment through the combination of ribavirin and pegylated interferon is expensive, only mildly effective, and is associated with severe side effects. In 2011, two direct-acting antiviral (DAA) drugs, boceprevir and telaprevir, were licensed that have shown enhanced sustained virologic response (SVR) in phase III clinical trial, however, these interferon-free treatments are more sensitive to HCV genotype 1 infection. The variable nature of HCV, and the limited number of inhibitors developed thus aim in expanding the repertoire of available drug targets, resulting in targeting the virus assembly therapeutically. AIM: We conducted this study to predict the 3D structure of the p7 protein from the HCV genotypes 3 and 4. Approximately 63 amino acid residues encoded in HCV render this channel sensitive to inhibitors, making p7 a promising target for novel therapies. HCV p7 protein forms a small membrane known as viroporin, and is essential for effective self-assembly of large channels that conduct cation assembly and discharge infectious virion particles. METHOD: In this study, we screened drugs and flavonoids known to disrupt translation and production of HCV proteins, targeted against the active site of p7 residues of HCV genotype 3 (GT3) (isolatek3a) and HCV genotype 4a (GT4) (isolateED43). Furthermore, we conducted a quantitative structure-activity relationship and docking interaction study. RESULTS: The drug NB-DNJ formed the highest number of hydrogen bond interactions with both modeled p7 proteins with high interaction energy, followed by BIT225. A flavonoid screen demonstrated that Epigallocatechin gallate (EGCG), nobiletin, and quercetin, have more binding modes in GT3 than in GT4. Thus, the predicted p7 protein molecule of HCV from GT3 and GT4 provides a general avenue to target structure-based antiviral compounds. CONCLUSIONS: We hypothesize that the inhibitors of viral p7 identified in this screen may be a new class of potent agents, but further confirmation in vitro and in vivo is essential. This structure-guided drug design for both GT3 and GT4 can lead to the identification of drug-like natural compounds, confirming p7 as a new target in the rapidly increasing era of HCV.

Fungal transformation of cedryl acetate and α-glucosidase inhibition assay, quantum mechanical calculations and molecular docking studies of its metabolites.

The fungal transformation of cedryl acetate (1) was investigated for the first time by using Cunninghamella elegans. The metabolites obtained include, 10ß-hydroxycedryl acetate (3), 2α, 10ß-dihydroxycedryl acetate (4), 2α-hydroxy-10-oxocedryl acetate (5), 3α,10ß-dihydroxycedryl acetate (6), 3α,10α-dihydroxycedryl acetate (7), 10ß,14α-dihydroxy cedryl acetate (8), 3ß,10ß-cedr-8(15)-ene-3,10-diol (9), and 3α,8ß,10ß -dihydroxycedrol (10). Compounds 1, 2, and 4 showed α-glucosidase inhibitory activity, whereby 1 was more potent than the standard inhibitor, acarbose, against yeast α-glucosidase. Detailed docking studies were performed on all experimentally active compounds to study the molecular interaction and binding mode in the active site of the modeled yeast α-glucosidase and human intestinal maltase glucoamylase. All active ligands were found to have greater binding affinity with the yeast α-glucosidase as compared to that of human homolog, the intestinal maltase, by an average value of approximately -1.4 kcal/mol, however, no significant difference was observed in the case of pancreatic amylase.