RESUMEN

Glioblastoma (GB) is an aggressive brain malignancy characterized by its invasive nature. Current treatment has limited effectiveness, resulting in poor patients' prognoses. ß-Amino carbonyl (ß-AC) compounds have gained attention due to their potential anticancerous properties. In vitro assays were performed to evaluate the effects of an in-house synthesized ß-AC compound, named SHG-8, upon GB cells. Small RNA sequencing (sRNA-seq) and biocomputational analyses investigated the effects of SHG-8 upon the miRNome and its bioavailability within the human body. SHG-8 exhibited significant cytotoxicity and inhibition of cell migration and proliferation in U87MG and U251MG GB cells. GB cells treated with the compound released significant amounts of reactive oxygen species (ROS). Annexin V and acridine orange/ethidium bromide staining also demonstrated that the compound led to apoptosis. sRNA-seq revealed a shift in microRNA (miRNA) expression profiles upon SHG-8 treatment and significant upregulation of miR-3648 and downregulation of miR-7973. Real-time polymerase chain reaction (RT-qPCR) demonstrated a significant downregulation of CORO1C, an oncogene and a player in the Wnt/ß-catenin pathway. In silico analysis indicated SHG-8's potential to cross the blood-brain barrier. We concluded that SHG-8's inhibitory effects on GB cells may involve the deregulation of various miRNAs and the inhibition of CORO1C.

RESUMEN

Since the discovery of the anti-influenza drugs oseltamivir and zanamivir using computer-aided drug design methods, there have been significant applications of molecular modelling methodologies applied to influenza A virus drug discovery, such as molecular dynamics (MD) simulation, molecular docking, and virtual screening (VS). In this review, we provide a brief general introduction to molecular modelling in the context of drug discovery and then focus on the advances and impact of integrating these methods with specific reference to potential influenza A antiviral drug targets.

Asunto(s)

RESUMEN

The nuclear export protein (NEP) of the influenza A virus exports viral ribonucleoproteins to the host cell cytoplasm following nuclear transcription. In this work conservation analysis of 3000 protein sequences and molecular modelling of full-length NEP identified ligand binding sites overlapping with high sequence conservation. Two binding hot spots were identified close to the first nuclear export signal and several hot spots overlapped with highly conserved amino acids such as Arg42, Asp43, Lys39, Ile80, Gln101 and Val109. Virtual screening with ~43,000 compounds against a binding site showed affinities of up to -8.95 kcal/mol, while ~1700 approved drugs showed affinities of up to -8.31 kcal/mol. A drug-like compounds predicted was ZINC01564229 that could be used as probe to investigate NEP function or as a new drug lead. The approved drugs Nandrolone phenylpropionate and Estropipate were predicted to bind with high affinity and may be investigated for repurposing as anti-influenza drugs. IMPORTANCE: The influenza A virus causes respiratory illness in humans and farm animals annually across the world. Antigenic shifts and drifts in the surface proteins lead to genome diversity and unpredictable pandemics and epidemics. The high evolution rate of the RNA genome can also limit the effectiveness of antivirals and is the cause of emerging resistance. From a human health perspective, it is important that compounds identified as potential influenza replication inhibitors remain effective long-term. This work presents results which are based on computational predictions that reveal interactions between available compounds and regions of the influenza A nuclear export protein which display high conservation. Due to a low probability of highly conserved regions undergoing genomic changes, these compounds may serve as ideal leads for new antivirals.

Asunto(s)

RESUMEN

Genes coding for nucleotide-binding leucine-rich repeat (LRR) receptors (NLRs) control resistance against intracellular (cell-penetrating) pathogens. However, evidence for a role of genes coding for proteins with LRR domains in resistance against extracellular (apoplastic) fungal pathogens is limited. Here, the distribution of genes coding for proteins with eLRR domains but lacking kinase domains was determined for the Brassica napus genome. Predictions of signal peptide and transmembrane regions divided these genes into 184 coding for receptor-like proteins (RLPs) and 121 coding for secreted proteins (SPs). Together with previously annotated NLRs, a total of 720 LRR genes were found. Leptosphaeria maculans-induced expression during a compatible interaction with cultivar Topas differed between RLP, SP and NLR gene families; NLR genes were induced relatively late, during the necrotrophic phase of pathogen colonization. Seven RLP, one SP and two NLR genes were found in Rlm1 and Rlm3/Rlm4/Rlm7/Rlm9 loci for resistance against L. maculans on chromosome A07 of B. napus. One NLR gene at the Rlm9 locus was positively selected, as was the RLP gene on chromosome A10 with LepR3 and Rlm2 alleles conferring resistance against L. maculans races with corresponding effectors AvrLm1 and AvrLm2, respectively. Known loci for resistance against L. maculans (extracellular hemi-biotrophic fungus), Sclerotinia sclerotiorum (necrotrophic fungus) and Plasmodiophora brassicae (intracellular, obligate biotrophic protist) were examined for presence of RLPs, SPs and NLRs in these regions. Whereas loci for resistance against P. brassicae were enriched for NLRs, no such signature was observed for the other pathogens. These findings demonstrate involvement of (i) NLR genes in resistance against the intracellular pathogen P. brassicae and a putative NLR gene in Rlm9-mediated resistance against the extracellular pathogen L. maculans.

Asunto(s)

RESUMEN

The influenza A basic polymerase protein 2 (PB2) functions as part of a heterotrimer to replicate the viral RNA genome. To investigate novel PB2 antiviral target sites, this work identified evolutionary conserved regions across the PB2 protein sequence amongst all sub-types and hosts, as well as ligand binding hot spots which overlap with highly conserved areas. Fifteen binding sites were predicted in different PB2 domains; some of which reside in areas of unknown function. Virtual screening of ~50,000 drug-like compounds showed binding affinities of up to -10.3kcal/mol. The highest affinity molecules were found to interact with conserved residues including Gln138, Gly222, Ile529, Asn540 and Thr530. A library containing 1738 FDA approved drugs was screened additionally and revealed Paliperidone as a top hit with a binding affinity of -10kcal/mol. Predicted ligands are ideal leads for new antivirals as they were targeted to evolutionary conserved binding sites.

Asunto(s)

RESUMEN

To determine if mutations introduced into phospholemman (PLM) could increase the level of PLM-Na,K-ATPase (NKA) binding, we performed scanning mutagenesis of the transmembrane domain of PLM and measured Förster resonance energy transfer (FRET) between each mutant and NKA. We observed an increased level of binding to NKA for several PLM mutants compared to that of the wild type (WT), including L27A, L30A, and I32A. In isolated cardiomyocytes, overexpression of WT PLM increased the amplitude of the Ca2+ transient compared to the GFP control. The Ca2+ transient amplitude was further increased by L30A PLM overexpression. The L30A mutation also delayed Ca2+ extrusion and increased the duration of cardiomyocyte contraction. This mimics aspects of the effect of cardiac glycosides, which are known to increase contractility through inhibition of NKA. No significant differences between WT and L30A PLM-expressing myocytes were observed after treatment with isoproterenol, suggesting that the superinhibitory effects of L30A are reversible with ß-adrenergic stimulation. We also observed a decrease in the extent of PLM tetramerization with L30A compared to WT using FRET, suggesting that L30 is an important residue for mediating PLM-PLM binding. Molecular dynamics simulations revealed that the potential energy of the L30A tetramer is greater than that of the WT, and that the transmembrane α helix is distorted by the mutation. The results identify PLM residue L30 as an important determinant of PLM tetramerization and of functional inhibition of NKA by PLM.

Asunto(s)

RESUMEN

Virtual screening is used in biomedical research to predict the binding affinity of a large set of small organic molecules to protein receptor targets. This report shows the development and evaluation of a novel yet straightforward attempt to improve this ranking in receptor-based molecular docking using a receptor-decoy strategy. This strategy includes defining a decoy binding site on the receptor and adjusting the ranking of the true binding-site virtual screen based on the decoy-site screen. The results show that by docking against a receptor-decoy site with Autodock Vina, improved Receiver Operator Characteristic Enrichment (ROCE) was achieved for 5 out of fifteen receptor targets investigated, when up to 15 % of a decoy site rank list was considered. No improved enrichment was seen for 7 targets, while for 3 targets the ROCE was reduced. The extent to which this strategy can effectively improve ligand prediction is dependent on the target receptor investigated.

Asunto(s)

RESUMEN

Sequence variations in the binding sites of influenza A proteins are known to limit the effectiveness of current antiviral drugs. Clinically, this leads to increased rates of virus transmission and pathogenicity. Potential influenza A inhibitors are continually being discovered as a result of high-throughput cell based screening studies, whereas the application of computational tools to aid drug discovery has further increased the number of predicted inhibitors reported. This review brings together the aspects that relate to the identification of influenza A drug target sites and the findings from recent antiviral drug discovery strategies.

Asunto(s)

RESUMEN

The molecular dynamics (MD) simulation of membrane proteins requires the setup of an accurate representation of lipid bilayers. This chapter describes the setup of a lipid bilayer system from scratch using generally available tools, starting with a definition of the lipid molecule POPE, generation of a lipid bilayer, energy minimization, MD simulation, and data analysis. The data analysis includes the calculation of area and volume per lipid, deuterium order parameters, self-diffusion constant, and the electron density profile.

Asunto(s)

RESUMEN

Cellulose-paper-based colorimetric bioassays may be used at the point of sampling without sophisticated equipment. This study reports the development of a colorimetric bioassay based on cellulose that can detect pathogen DNA. The detection was based on covalently attached single-stranded DNA probes and visual analysis. A cellulose surface functionalized with tosyl groups was prepared by the N,N-dimethylacetamide-lithium chloride method. Tosylation of cellulose was confirmed by scanning electron microscopy, Fourier transform infrared spectroscopy and elemental analysis. Sulfhydryl-modified oligonucleotide probes complementary to a segment of the DNA sequence IS6110 of Mycobacterium tuberculosis were covalently immobilized on the tosylated cellulose. On hybridization of biotin-labelled DNA oligonucleotides with these probes, a colorimetric signal was obtained with streptavidin-conjugated horseradish peroxidase catalysing the oxidation of tetramethylbenzamidine by H2O2. The colour intensity was significantly reduced when the bioassay was subjected to DNA oligonucleotide of randomized base composition. Initial experiments have shown a sensitivity of 0.1 µM. A high probe immobilization efficiency (more than 90 %) was observed with a detection limit of 0.1 µM, corresponding to an absolute amount of 10 pmol. The detection of M. tuberculosis DNA was demonstrated using this technique coupled with PCR for biotinylation of the DNA. This work shows the potential use of tosylated cellulose as the basis for point-of-sampling bioassays.

Asunto(s)

RESUMEN

G-quadruplexes are higher-order nucleic acid structures formed of square-planar arrangements of four guanine bases held together by Hoogsteen-type hydrogen bonds. Stacks of guanine tetrads are stabilised by intercalating potassium ions. FXYD1 encodes for phospholemman, a regulatory subunit of the cardiac Na(+)/K(+)-ATPase. Computational sequence analysis of FXYD1 pre-mRNA predicted the formation of stable intramolecular G-quadruplexes in human and orthologue sequences. Multiple sequence alignment indicated that G-rich sequences are conserved in evolution suggesting a potential role of G-quadruplexes in FXYD1 gene expression. The existence of a non-functional alternative splicing product indicated that the G-quadruplex formation may control alternative splicing. Quadruplex formation of human and bovine oligonucleotides was confirmed in vitro by native polyacrylamide gel electrophoresis and intrinsic fluorescence emission spectroscopy. Taking together the evolutionary conservation of G-quadruplex forming sequences with the confirmation of G-quadruplex formation in vitro by two FXYD1 homologues the results point to a potential role of these structures in regulating the expression of FXYD1 and thus regulate indirectly the activity of the cardiac Na(+)/K(+)-ATPase.

Asunto(s)

RESUMEN

The nucleoprotein (NP) of the influenza A virus encapsidates the viral RNA and participates in the infectious life cycle of the virus. The aims of this study were to find the degree of conservation of NP among all virus subtypes and hosts and to identify conserved binding sites, which may be utilised as potential drug target sites. The analysis of conservation based on 4430 amino acid sequences identified high conservation in known functional regions as well as novel highly conserved sites. Highly variable clusters identified on the surface of NP may be associated with adaptation to different hosts and avoidance of the host immune defence. Ligand binding potential overlapping with high conservation was found in the tail-loop binding site and near the putative RNA binding region. The results provide the basis for developing antivirals that may be universally effective and have a reduced potential to induce resistance through mutations.

Asunto(s)

RESUMEN

The etiology of Alzheimer's disease is thought to be linked to interactions between amyloid-ß (Aß) and neural cell membranes, causing membrane disruption and increased ion conductance. The effects of Aß on lipid behavior have been characterized experimentally, but structural and causal details are lacking. We used atomistic molecular dynamics simulations totaling over 6 µs in simulation time to investigate the behavior of Aß(42) in zwitterionic and anionic lipid bilayers. We simulated transmembrane ß-sheets (monomer and tetramer) resulting from a global optimization study and a helical structure obtained from an NMR study. In all simulations Aß(42) remained embedded in the bilayer. It was found that the surface charge and the lipid tail type are determinants for transmembrane stability of Aß(42) with zwitterionic surfaces and unsaturated lipids promoting stability. From the considered structures, the ß-sheet tetramer is most stable as a result of interpeptide interactions. We performed an in-depth analysis of the translocation of water in the Aß(42)-bilayer systems. We observed that this process is generally fast (within a few nanoseconds) yet generally slower than in the peptide-free bilayers. It is mainly governed by the lipid type, simulation temperature and Aß(42) conformation. The rate limiting step is the permeation through the hydrophobic core, where interactions between Aß(42) and permeating H(2)O molecules slow the translocation process. The ß-sheet tetramer allows more water molecules to pass through the bilayer compared to monomeric Aß, allowing us to conclude that the experimentally observed permeabilization of membranes must be due to membrane-bound Aß oligomers, and not monomers.

Asunto(s)

RESUMEN

The inhibitor κB kinase-ß (IKK-ß) phosphorylates the NF-κB inhibitor protein IκB leading to the translocation of the transcription factor NF-κB to the nucleus. The transcription factor NF-κB and consequently IKK-ß are central to signal transduction pathways of mammalian cells. The purpose of this research was to develop a 3D structural model of the IKK-ß kinase domain with its ATP cofactor and investigate its dynamics and ligand binding potential. Through a combination of comparative modelling and simulated heating/annealing molecular dynamics (SAMD) simulation in explicit water the model accuracy could be substantially improved compared to comparative modelling on its own as shown by model validation measures. The structure revealed the details of ATP/Mg(2+) binding indicating hydrophobic interactions with the adenine base and a significant contribution of Mg(2+) as a bridge between ATP phosphate groups and negatively charged side chains. The molecular dynamics trajectories of the ATP-bound and free enzyme showed two conformations in each case, which contributed to the majority of the trajectory. The ATP-free enzyme revealed a novel binding site distant from the ATP binding site that was not encountered in the ATP bound enzyme. Based on the overall structural flexibility, it is suggested that a truncated version of the kinase domain from Ala14 to Leu265 should be subjected to crystallisation trials. The 3D structure of this enzyme will enable rational design of new ligands and analysis of protein-protein interactions. Furthermore, our results may provide a new impetus for wet-lab based structural investigation focussing on a truncated kinase domain.

Asunto(s)

RESUMEN

In order to exploit the advantages of receptor-based virtual screening, namely time/cost saving and specificity, it is important to rely on algorithms that predict a high number of active ligands at the top ranks of a small molecule database. Towards that goal consensus methods combining the results of several docking algorithms were developed and compared against the individual algorithms. Furthermore, a recently proposed rescoring method based on drug efficiency indices was evaluated. Among AutoDock Vina 1.0, AutoDock 4.2 and GemDock, AutoDock Vina was the best performing single method in predicting high affinity ligands from a database of known ligands and decoys. The rescoring of predicted binding energies with the water/octanol partition coefficient did not lead to an improvement averaged over ten receptor targets. Various consensus algorithms were investigated and a simple combination of AutoDock and AutoDock Vina results gave the most consistent performance that showed early enrichment of known ligands for all receptor targets investigated. In case a number of ligands is known for a specific target, every method proposed in this study should be evaluated.

Asunto(s)

RESUMEN

The non-structural protein 1 (NS1) of the influenza A virus and the NS2 protein, which is also known as nuclear export protein, play important roles in the infectious life cycle of the virus. The objective of this study was to find the degree of conservation in the NS proteins and to identify conserved sites of functional or structural importance that may be utilized as potential drug target sites. The analysis was based on 2620 amino acid sequences for the NS1 protein and 1195 sequences for the NS2 protein. The degree of conservation and potential binding sites were mapped onto the protein structures obtained from a combination of experimentally available structure fragments with predicted threading models. In addition to high conservation in protein regions of known function, novel highly conserved sites have been identified, namely Glu159, Thr171, Val192, Arg200, Glu208 and Gln218 on the NS1 protein and Ser24, Leu28, Arg66, Arg84, Ser93, Ile97 and Leu103 on the NS2 protein. Using the Q-SiteFinder binding site prediction algorithm, several highly conserved binding sites were found, including two spatially close sites on the NS1 protein, which could be targeted with a bivalent ligand that would interfere with double-stranded RNA binding. Altogether, this work reveals novel universally conserved residues that are candidates for protein-protein interactions and provide the basis for designing universal anti-influenza drugs.

Asunto(s)

RESUMEN

United-atom force fields for molecular dynamics (MD) simulations provide a higher computational efficiency, especially in lipid membrane simulations, with little sacrifice in accuracy, when compared to all-atom force fields. Excellent united-atom lipid models are available, but in combination with depreciated protein force fields. In this work, a united-atom model of the lipid 1,2-dipalmitoyl-sn-glycero-3-phosphocholine has been built with standard parameters of the force field GROMOS96 53a6 that reproduces the experimental area per lipid of a lipid bilayer within 3% accuracy to a value of 0.623 ± 0.011 nm(2) without the assumption of a constant surface area or the inclusion of surface pressure. In addition, the lateral self-diffusion constant and deuterium order parameters of the acyl chains are in agreement with experimental data. Furthermore, models for 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) result in areas per lipid of 0.625 nm(2) (DMPC), 0.693 nm(2) (POPC), and 0.700 nm(2) (POPG) from 40 ns MD simulations. Experimental lateral self-diffusion coefficients are reproduced satisfactorily by the simulation. The lipid models can form the basis for molecular dynamics simulations of membrane proteins with current and future versions of united-atom protein force fields.

RESUMEN

Elucidating the structure of transmembrane proteins domains with high-resolution methods is a difficult and sometimes impossible task. Here, we explain the method of combining a limited amount of experimental data with automated high-throughput molecular dynamics (MD) simulations of alpha-helical transmembrane bundles in an explicit lipid bilayer/water environment. The procedure uses a systematic conformational search of the helix rotation with experimentally constrained MDs simulations. The experimentally determined helix tilt and rotational angle of a labeled residue with site-specific infrared dichroism allows us to select a unique high-resolution model from a number of possible energy minima encountered in the systematic conformational search.

Asunto(s)

RESUMEN

The 17O-'diluted' glycine-14 sites in a phospholemman (PLM) transmembrane domain protein are characterized by solid-state 17O NMR spectroscopy. The PLM transmembrane domain is an alpha-helical tetramer unit of four 28-residue peptides and is rigidly embedded in a bilayer where each alpha-helix has an average tilt of 7.3 degrees against the membrane normal. The PLM sample investigated here consists of a high lipid/peptide molar ratio (25:1) with one glycine residue in each helix enriched to <40% (17)O; thus, this is a very dilute 17O-sample and is the most dilute 17O-membrane protein to date to be characterized by solid-state 17O NMR spectroscopy. Based on the spectral analysis of 17O magic angle spinning (MAS) at 14.1 and 18.8T, the PLM transmembrane domain protein consists of multiple crystallographic gly14 sites, suggesting that the tetramer protein is an asymmetric unit with either C2- or C1-rotational symmetry along the bilayer normal.

Asunto(s)

RESUMEN

The structure of ghrelin, a 28-residue octanoylated peptide hormone, is only known up to the level of primary structure identifying an active core of residues 1-5 or 1-4 including octanoyl-Ser3 as necessary to elicit receptor response. This chapter reviews the results and limitations of experimental and computer modeling studies, which have appeared in the literature. The (1)H NMR spectroscopy experimental studies revealed an unstructured and/or fast interconverting peptide at acidic pH, while molecular dynamics (MD) simulation studies at neutral pH pointed to a stable conformation over a time period of 25 ns in water and in the presence of a lipid bilayer. The significance of these findings is discussed with regards to the pH difference, the timescales accessible to simulation and NMR spectroscopy, and the limitations of computational modeling. MD simulations of ghrelin in the presence of a lipid membrane revealed that the octanoyl side chain did not insert into the lipid bilayer, but instead the peptide bound to the lipid headgroups with residues Arg15, Lys16, Glu17, and Ser18, which are located in a hairpin-like bend in the structure. The implications of these findings with regards to a recently obtained homology model of the ghrelin receptor are discussed.

Asunto(s)