Molecular modelling and de novo fragment-based design of potential inhibitors of beta-tubulin gene of Necator americanus from natural products.

The emergence of drug resistance against the known hookworm drugs namely albendazole and mebendazole and their reduced efficacies necessitate the need for new drugs. Chemically diverse natural products present plausible templates to augment hookworm drug discovery. The present work utilized pharmacoinformatics techniques to predict African natural compounds ZINC95486082, ZINC95486052 and euphohelionon as potential inhibitory molecules of the hookworm Necator americanus ß tubulin gene. A library of 3390 compounds was screened against a homology-modelled structure of ß tubulin. The docking results obtained from AutoDock Vina was validated with an acceptable area under the curve (AUC) of 0.714 computed from the receiver operating characteristic (ROC) curve. The three selected compounds had favourable binding affinities and were predicted to form no interactions with the resistance-associated mutations Phe167, Glu198 and Phe200. The compounds were predicted as anthelmintics using a Bayesian-based technique and were pharmacologically profiled to be druglike. Further molecular dynamics simulations and MM-PBSA calculations showed the compounds as promising anthelmintic drug leads. Novel critical residues comprising Leu246, Asn247 and Asn256 were also predicted for binding. Euphohelionon was selected as a template for the de novo fragment-based design of five compounds labelled A1, A2, A3, A4 and A5; with four of them having SAscore values below 6, denoting easy synthesis. All the five de novo molecules docked firmly in the binding pocket of the ß tubulin with no binding interactions with the three known resistance mutation residues. Binding energies of -8.2, -7.6, -7.3, -7.2 and -6.8 kcal/mol were obtained for A1, A2, A3, A4 and A5, respectively. The identified compounds can serve as treasure troves from which future potent anthelmintics can be designed. The current study strives to assuage the hookworm disease burden, especially making available molecules with the potential to circumvent the chemoresistance.

Computational Study on Potential Novel Anti-Ebola Virus Protein VP35 Natural Compounds.

Ebola virus (EBOV) is one of the most lethal pathogens that can infect humans. The Ebola viral protein VP35 (EBOV VP35) inhibits host IFN-α/ß production by interfering with host immune responses to viral invasion and is thus considered as a plausible drug target. The aim of this study was to identify potential novel lead compounds against EBOV VP35 using computational techniques in drug discovery. The 3D structure of the EBOV VP35 with PDB ID: 3FKE was used for molecular docking studies. An integrated library of 7675 African natural product was pre-filtered using ADMET risk, with a threshold of 7 and, as a result, 1470 ligands were obtained for the downstream molecular docking using AutoDock Vina, after an energy minimization of the protein via GROMACS. Five known inhibitors, namely, amodiaquine, chloroquine, gossypetin, taxifolin and EGCG were used as standard control compounds for this study. The area under the curve (AUC) value, evaluating the docking protocol obtained from the receiver operating characteristic (ROC) curve, generated was 0.72, which was considered to be acceptable. The four identified potential lead compounds of NANPDB4048, NANPDB2412, ZINC000095486250 and NANPDB2476 had binding affinities of -8.2, -8.2, -8.1 and -8.0 kcal/mol, respectively, and were predicted to possess desirable antiviral activity including the inhibition of RNA synthesis and membrane permeability, with the probable activity (Pa) being greater than the probable inactivity (Pi) values. The predicted anti-EBOV inhibition efficiency values (IC50), found using a random forest classifier, ranged from 3.35 to 11.99 µM, while the Ki values ranged from 0.97 to 1.37 µM. The compounds NANPDB4048 and NANPDB2412 had the lowest binding energy of -8.2 kcal/mol, implying a higher binding affinity to EBOV VP35 which was greater than those of the known inhibitors. The compounds were predicted to possess a low toxicity risk and to possess reasonably good pharmacological profiles. Molecular dynamics (MD) simulations of the protein-ligand complexes, lasting 50 ns, and molecular mechanisms Poisson-Boltzmann surface area (MM-PBSA) calculations corroborated the binding affinities of the identified compounds and identified novel critical interacting residues. The antiviral potential of the molecules could be confirmed experimentally, while the scaffolds could be optimized for the design of future novel anti-EBOV chemotherapeutics.

Phosphorylated VP30 of Marburg Virus Is a Repressor of Transcription

The filoviruses Marburg virus (MARV) and Ebola virus (EBOV) cause hemorrhagic fever in humans and nonhuman primates, with high case fatality rates. MARV VP30 is known to be phosphorylated and to interact with nucleoprotein (NP), but its role in regulation of viral transcription is disputed. Here, we analyzed phosphorylation of VP30 by mass spectrometry, which resulted in identification of multiple phosphorylated amino acids. Modeling the full-length three-dimensional structure of VP30 and mapping the identified phosphorylation sites showed that all sites lie in disordered regions, mostly in the N-terminal domain of the protein. Minigenome analysis of the identified phosphorylation sites demonstrated that phosphorylation of a cluster of amino acids at positions 46 through 53 inhibits transcription. To test the effect of VP30 phosphorylation on its interaction with other MARV proteins, coimmunoprecipitation analyses were performed. They demonstrated the involvement of VP30 phosphorylation in interaction with two other proteins of the MARV ribonucleoprotein complex, NP and VP35. To identify the role of protein phosphatase 1 (PP1) in the identified effects, a small molecule, 1E7-03, targeting a noncatalytic site of the enzyme that previously was shown to increase EBOV VP30 phosphorylation was used. Treatment of cells with 1E7-03 increased phosphorylation of VP30 at a cluster of phosphorylated amino acids from Ser-46 to Thr-53, reduced transcription of MARV minigenome, enhanced binding to NP and VP35, and dramatically reduced replication of infectious MARV particles. Thus, MARV VP30 phosphorylation can be targeted for development of future antivirals such as PP1-targeting compounds. IMPORTANCE The largest outbreak of MARV occurred in Angola in 2004 to 2005 and had a 90% case fatality rate. There are no approved treatments available for MARV. Development of antivirals as therapeutics requires a fundamental understanding of the viral life cycle. Because of the close similarity of MARV to another member of Filoviridae family, EBOV, it was assumed that the two viruses have similar mechanisms of regulation of transcription and replication. Here, characterization of the role of VP30 and its phosphorylation sites in transcription of the MARV genome demonstrated differences from those of EBOV. The identified phosphorylation sites appeared to inhibit transcription and appeared to be involved in interaction with both NP and VP35 ribonucleoproteins. A small molecule targeting PP1 inhibited transcription of the MARV genome, effectively suppressing replication of the viral particles. These data demonstrate the possibility developing antivirals based on compounds targeting PP1.

Determination of structural models of the complex between the cytoplasmic domain of erythrocyte band 3 and ankyrin-R repeats 13-24.

The adaptor protein ankyrin-R interacts via its membrane binding domain with the cytoplasmic domain of the anion exchange protein (AE1) and via its spectrin binding domain with the spectrin-based membrane skeleton in human erythrocytes. This set of interactions provides a bridge between the lipid bilayer and the membrane skeleton, thereby stabilizing the membrane. Crystal structures for the dimeric cytoplasmic domain of AE1 (cdb3) and for a 12-ankyrin repeat segment (repeats 13-24) from the membrane binding domain of ankyrin-R (AnkD34) have been reported. However, structural data on how these proteins assemble to form a stable complex have not been reported. In the current studies, site-directed spin labeling, in combination with electron paramagnetic resonance (EPR) and double electron-electron resonance, has been utilized to map the binding interfaces of the two proteins in the complex and to obtain inter-protein distance constraints. These data have been utilized to construct a family of structural models that are consistent with the full range of experimental data. These models indicate that an extensive area on the peripheral domain of cdb3 binds to ankyrin repeats 18-20 on the top loop surface of AnkD34 primarily through hydrophobic interactions. This is a previously uncharacterized surface for binding of cdb3 to AnkD34. Because a second dimer of cdb3 is known to bind to ankyrin repeats 7-12 of the membrane binding domain of ankyrin-R, the current models have significant implications regarding the structural nature of a tetrameric form of AE1 that is hypothesized to be involved in binding to full-length ankyrin-R in the erythrocyte membrane.

A new model defines the minimal set of polymorphism in HLA-DQ and -DR that determines susceptibility and resistance to autoimmune diabetes.

BACKGROUND: The mechanism underlying autoimmune diabetes has been difficult to define. There is a strong genetic contribution and numerous studies associate the major histocompatibility complex, especially the class II region, with predisposition or resistance. However, how these molecules are implicated remains obscure. PRESENTATION OF THE HYPOTHESIS: We have supplemented structural analysis with computational biophysical and sequence analyses and propose an heuristic for distinguishing between human leukocyte antigen molecules that predispose to insulin dependent diabetes mellitus and those that are protective. Polar residues at both beta37 and beta9 suffice to distinguish accurately between class II alleles that predispose to type 1 diabetes and those that do not. The electrostatic potential within the peptide binding pocket exerts a strong influence on diabetogenic epitopes with basic residues. Diabetes susceptibility alleles are predicted to bind autoantigens strongly with tight affinity, prolonged association and altered cytokine expression profile. Protective alleles bind moderately, and neutral alleles poorly or not at all. Non-Asp beta57 is a modifier that supplements disease risk but only in the presence of the polymorphic, polar pair at beta9 and beta37. The nature of beta37 determines resistance on one hand, and susceptibility or dominant protection on the other. CONCLUSION: The proposed ideas are illustrated with structural, functional and population studies from the literature. The hypothesis, in turn, rationalizes their results. A plausible mechanism of immune mediated diabetes based on binding affinity and peptide kinetics is discussed. The number of the polymorphic markers present correlates with onset of disease and severity. The molecular elucidation of disease susceptibility and resistance paves the way for risk prediction, treatment and prevention of disease based on analogue peptides. REVIEWERS: This article was reviewed by Eugene V. Koonin, Michael Lenardo, Hossam Ashour, and Bhagirath Singh. For the full reviews, please go to the Reviewers' comments section.

Flanking p10 contribution and sequence bias in matrix based epitope prediction: revisiting the assumption of independent binding pockets.

BACKGROUND: Eluted natural peptides from major histocompatibility molecules show patterns of conserved residues. Crystallographic structures show that the bound peptide in class II major histocompatibility complex adopts a near uniform polyproline II-like conformation. This way allele-specific favoured residues are able to anchor into pockets in the binding groove leaving other peptide side chains exposed for recognition by T cells. The anchor residues form a motif. This sequence pattern can be used to screen large sequences for potential epitopes. Quantitative matrices extend the motif idea to include the contribution of non-anchor peptide residues. This report examines two new matrices that extend the binding register to incorporate the polymorphic p10 pocket of human leukocyte antigen DR1. Their performance is quantified against experimental binding measurements and against the canonical nine-residue register matrix. RESULTS: One new matrix shows significant improvement over the base matrix; the other does not. The new matrices differ in the sequence of the peptide library. CONCLUSION: One of the extended quantitative matrices showed significant improvement in prediction over the original nine residue matrix and over the other extended matrix. Proline in the sequence of the peptide library of the better performing matrix presumably stabilizes the peptide conformation through neighbour interactions. Such interactions may influence epitope prediction in this test of quantitative matrices. This calls into question the assumption of the independent contribution of individual binding pockets.

Crystallographic structure of the human leukocyte antigen DRA, DRB3*0101: models of a directional alloimmune response and autoimmunity.

We describe structural studies of the human leukocyte antigen DR52a, HLA-DRA/DRB3*0101, in complex with an N-terminal human platelet integrin alphaII(B)betaIII glycoprotein peptide which contains a Leu/Pro dimorphism. The 33:Leu dimorphism is the epitope for the T cell directed response in neonatal alloimmune thrombocytopenia and post-transfusion purpura in individuals with the alphaII(B)betaIII 33:Pro allele, and defines the unidirectional alloimmune response. This condition is always associated with DR52a. The crystallographic structure has been refined to 2.25 A. There are two alphabeta heterodimers to the asymmetric unit in space group P4(1)2(1)2. The molecule is characterized by two prominent hydrophobic pockets at either end of the peptide binding cleft and a deep, narrower and highly charged P4 opening underneath the beta 1 chain. Further, the peptide in the second molecule displays a sharp upward turn after pocket P9. The structure reveals the role of pockets and the distinctive basic P4 pocket, shared by DR52a and DR3, in selecting their respective binding peptide repertoire. We observe an interesting switch in a residue from the canonically assigned pocket 6 seen in prior class II structures to pocket 4. This occludes the P6 pocket helping to explain the distinctive "1-4-9" peptide binding motif. A beta57 Asp-->Val substitution abrogates the salt-bridge to alpha76 Arg and along with a hydrophobic beta37 is important in shaping the P9 pocket. DRB3*0101 and DRB1*0301 belong to an ancestral haplotype and are associated with many autoimmune diseases linked to antigen presentation, but whereas DR3 is susceptible to type 1 diabetes DR52a is not. This dichotomy is explored for clues to the disease.