Antiviral Mechanisms of N-Phenyl Benzamides on Coxsackie Virus A9.

Enteroviruses are one of the most abundant groups of viruses infecting humans, and yet there are no approved antivirals against them. To find effective antiviral compounds against enterovirus B group viruses, an in-house chemical library was screened. The most effective compounds against Coxsackieviruses B3 (CVB3) and A9 (CVA9) were CL212 and CL213, two N-phenyl benzamides. Both compounds were more effective against CVA9 and CL213 gave a better EC50 value of 1 µM with high a specificity index of 140. Both drugs were most effective when incubated directly with viruses suggesting that they mainly bound to the virions. A real-time uncoating assay showed that the compounds stabilized the virions and radioactive sucrose gradient as well as TEM confirmed that the viruses stayed intact. A docking assay, taking into account larger areas around the 2-and 3-fold axes of CVA9 and CVB3, suggested that the hydrophobic pocket gives the strongest binding to CVA9 but revealed another binding site around the 3-fold axis which could contribute to the binding of the compounds. Together, our data support a direct antiviral mechanism against the virus capsid and suggest that the compounds bind to the hydrophobic pocket and 3-fold axis area resulting in the stabilization of the virion.

Structure-based discovery and in vitro validation of inhibitors of chloride intracellular channel 4 protein.

The use of computer-aided methods have continued to propel accelerated drug discovery across various disease models, interestingly allowing the specific inhibition of pathogenic targets. Chloride Intracellular Channel Protein 4 (CLIC4) is a novel class of intracellular ion channel highly implicated in tumor and vascular biology. It regulates cell proliferation, apoptosis and angiogenesis; and is involved in multiple pathologic signaling pathways. Absence of specific inhibitors however impedes its advancement to translational research. Here, we integrate structural bioinformatics and experimental research approaches for the discovery and validation of small-molecule inhibitors of CLIC4. High-affinity allosteric binders were identified from a library of 1615 Food and Drug Administration (FDA)-approved drugs via a high-performance computing-powered blind-docking approach, resulting in the selection of amphotericin B and rapamycin. NMR assays confirmed the binding and conformational disruptive effects of both drugs while they also reversed stress-induced membrane translocation of CLIC4 and inhibited endothelial cell migration. Structural and dynamics simulation studies further revealed that the inhibitory mechanisms of these compounds were hinged on the allosteric modulation of the catalytic glutathione (GSH)-like site loop and the extended catalytic ß loop which may elicit interference with the catalytic activities of CLIC4. Structure-based insights from this study provide the basis for the selective targeting of CLIC4 to treat the associated pathologies.

Preliminary Structure-Activity Relationship Study of the MMV Pathogen Box Compound MMV675968 (2,4-Diaminoquinazoline) Unveils Novel Inhibitors of Trypanosoma brucei brucei.

New drugs are urgently needed for the treatment of human African trypanosomiasis (HAT). In line with our quest for novel inhibitors of trypanosomes, a small library of analogs of the antitrypanosomal hit (MMV675968) available at MMV as solid materials was screened for antitrypanosomal activity. In silico exploration of two potent antitrypanosomal structural analogs (7-MMV1578647 and 10-MMV1578445) as inhibitors of dihydrofolate reductase (DHFR) was achieved, together with elucidation of other antitrypanosomal modes of action. In addition, they were assessed in vitro for tentative inhibition of DHFR in a crude trypanosome extract. Their ADMET properties were also predicted using dedicated software. Overall, the two diaminoquinazoline analogs displayed approximately 40-fold and 60-fold more potency and selectivity in vitro than the parent hit, respectively (MMV1578445 (10): IC50 = 0.045 µM, SI = 1737; MMV1578467 (7): IC50 = 0.06 µM; SI = 412). Analogs 7 and 10 were also strong binders of the DHFR enzyme in silico, in all their accessible protonation states, and interacted with key DHFR ligand recognition residues Val32, Asp54, and Ile160. They also exhibited significant activity against trypanosome protein isolate. MMV1578445 (10) portrayed fast and irreversible trypanosome growth arrest between 4-72 h at IC99. Analogs 7 and 10 induced in vitro ferric iron reduction and DNA fragmentation or apoptosis induction, respectively. The two potent analogs endowed with predicted suitable physicochemical and ADMET properties are good candidates for further deciphering their potential as starting points for new drug development for HAT.

Virtual screening and in vitro validation identifies the first reported inhibitors of Salmonella enterica HPPK.

HPPK, which directly precedes DHPS in the folate biosynthetic pathway, is a promising but chronically under-exploited anti-microbial target. Here we report the identification of new S. enterica HPPK inhibitors, offering potential for new resistance circumventing S. enterica therapies as well as avenues for diversifying the current HPPK inhibitor space.

In Silico Investigation of Potential Applications of Gamma Carbonic Anhydrases as Catalysts of CO2 Biomineralization Processes: A Visit to the Thermophilic Bacteria Persephonella hydrogeniphila, Persephonella marina, Thermosulfidibacter takaii, and Thermus thermophilus.

Carbonic anhydrases (CAs) have been identified as ideal catalysts for CO2 sequestration. Here, we report the sequence and structural analyses as well as the molecular dynamics (MD) simulations of four γ-CAs from thermophilic bacteria. Three of these, Persephonella marina, Persephonella hydrogeniphila, and Thermosulfidibacter takaii originate from hydrothermal vents and one, Thermus thermophilus HB8, from hot springs. Protein sequences were retrieved and aligned with previously characterized γ-CAs, revealing differences in the catalytic pocket residues. Further analysis of the structures following homology modeling revealed a hydrophobic patch in the catalytic pocket, presumed important for CO2 binding. Monitoring of proton shuttling residue His69 (P. marina γ-CA numbering) during MD simulations of P. hydrogeniphila and P. marina's γ-CAs (γ-PhCA and γ-PmCA), showed a different behavior to that observed in the γ-CA of Escherichia coli, which periodically coordinates Zn2+. This work also involved the search for hotspot residues that contribute to interface stability. Some of these residues were further identified as key in protein communication via betweenness centrality metric of dynamic residue network analysis. T. takaii's γ-CA showed marginally lower thermostability compared to the other three γ-CA proteins with an increase in conformations visited at high temperatures being observed. Hydrogen bond analysis revealed important interactions, some unique and others common in all γ-CAs, which contribute to interface formation and thermostability. The seemingly thermostable γ-CA from T. thermophilus strangely showed increased unsynchronized residue motions at 423 K. γ-PhCA and γ-PmCA were, however, preliminarily considered suitable as prospective thermostable CO2 sequestration agents.

Mechanism of Action of Non-Synonymous Single Nucleotide Variations Associated with α-Carbonic Anhydrase II Deficiency.

Human carbonic anhydrase II (CA-II) is a Zinc (Zn 2 + ) metalloenzyme responsible for maintenance of acid-base balance within the body through the reversible hydration of CO 2 to produce protons (H + ) and bicarbonate (BCT). Due to its importance, alterations to the amino acid sequence of the protein as a result of single nucleotide variations (nsSNVs) have detrimental effects on homeostasis. Six pathogenic CA-II nsSNVs, K18E, K18Q, H107Y, P236H, P236R and N252D were identified, and variant protein models calculated using homology modeling. The effect of each nsSNV was analyzed using motif analysis, molecular dynamics (MD) simulations, principal component (PCA) and dynamic residue network (DRN) analysis. Motif analysis identified 11 functionally important motifs in CA-II. RMSD data indicated subtle SNV effects, while PCA analysis revealed that the presence of BCT results in greater conformational sampling and free energy in proteins. DRN analysis showed variant allosteric effects, and the average betweenness centrality (BC) calculations identified Glu117 as the most important residue for communication in CA-II. The presence of BCT was associated with a reduction to Glu117 usage in all variants, suggesting implications for Zn 2 + dissociation from the CA-II active site. In addition, reductions to Glu117 usage are associated with increases in the usage of the primary and secondary Zn 2 + ligands; His94, His96, His119 and Asn243 highlighting potential compensatory mechanisms to maintain Zn 2 + within the active site. Compared to traditional MD simulation investigation, DRN analysis provided greater insights into SNV mechanism of action, indicating its importance for the study of missense mutation effects in proteins and, in broader terms, precision medicine related research.

South African Abietane Diterpenoids and Their Analogs as Potential Antimalarials: Novel Insights from Hybrid Computational Approaches.

The hemoglobin degradation process in Plasmodium parasites is vital for nutrient acquisition required for their growth and proliferation. In P. falciparum, falcipains (FP-2 and FP-3) are the major hemoglobinases, and remain attractive antimalarial drug targets. Other Plasmodium species also possess highly homologous proteins to FP-2 and FP-3. Although several inhibitors have been designed against these proteins, none has been commercialized due to associated toxicity on human cathepsins (Cat-K, Cat-L and Cat-S). Despite the two enzyme groups sharing a common structural fold and catalytic mechanism, distinct active site variations have been identified, and can be exploited for drug development. Here, we utilize in silico approaches to screen 628 compounds from the South African natural sources to identify potential hits that can selectively inhibit the plasmodial proteases. Using docking studies, seven abietane diterpenoids, binding strongly to the plasmodial proteases, and three additional analogs from PubChem were identified. Important residues involved in ligand stabilization were identified for all potential hits through binding pose analysis and their energetic contribution determined by binding free energy calculations. The identified compounds present important scaffolds that could be further developed as plasmodial protease inhibitors. Previous laboratory assays showed the effect of the seven diterpenoids as antimalarials. Here, for the first time, we demonstrate that their possible mechanism of action could be by interacting with falcipains and their plasmodial homologs. Dynamic residue network (DRN) analysis on the plasmodial proteases identified functionally important residues, including a region with high betweenness centrality, which had previously been proposed as a potential allosteric site in FP-2.

Homology modeling and docking of AahII-Nanobody complexes reveal the epitope binding site on AahII scorpion toxin.

Scorpion envenoming and its treatment is a public health problem in many parts of the world due to highly toxic venom polypeptides diffusing rapidly within the body of severely envenomed victims. Recently, 38 AahII-specific Nanobody sequences (Nbs) were retrieved from which the performance of NbAahII10 nanobody candidate, to neutralize the most poisonous venom compound namely AahII acting on sodium channels, was established. Herein, structural computational approach is conducted to elucidate the Nb-AahII interactions that support the biological characteristics, using Nb multiple sequence alignment (MSA) followed by modeling and molecular docking investigations (RosettaAntibody, ZDOCK software tools). Sequence and structural analysis showed two dissimilar residues of NbAahII10 CDR1 (Tyr27 and Tyr29) and an inserted polar residue Ser30 that appear to play an important role. Indeed, CDR3 region of NbAahII10 is characterized by a specific Met104 and two negatively charged residues Asp115 and Asp117. Complex dockings reveal that NbAahII17 and NbAahII38 share one common binding site on the surface of the AahII toxin divergent from the NbAahII10 one's. At least, a couple of NbAahII10 - AahII residue interactions (Gln38 - Asn44 and Arg62, His64, respectively) are mainly involved in the toxic AahII binding site. Altogether, this study gives valuable insights in the design and development of next generation of antivenom.

The generation and characterisation of neutralising antibodies against the Theiler's murine encephalomyelitis virus (TMEV) GDVII capsid reveals the potential binding site of the host cell co-receptor, heparan sulfate.

The early stages of picornavirus capsid assembly and the host factors involved are poorly understood. Since the localisation of viral proteins in infected cells can provide information on their function, antibodies against purified Theiler's murine encephalomyelitis virus (TMEV) GDVII capsids were generated by immunisation of rabbits. The resultant anti-TMEV capsid antibodies recognised a C-terminal region of VP1 but not VP2 or VP3 by Western analysis. Examination of the sites of TMEV capsid assembly by indirect immunofluorescence and confocal microscopy showed that at 5h post infection, capsid signal was diffusely cytoplasmic with strong perinuclear staining and moved into large punctate structures from 6 to 8h post infection. A plaque reduction neutralisation assay showed that the anti-TMEV capsid antibodies but not anti-VP1 antibodies could neutralise viral infection in vitro. The VP1 C-terminal residues recognised by the anti-TMEV capsid antibodies were mapped to a loop on the capsid surface near to the putative receptor binding pocket. In silico docking experiments showed that the known TMEV co-receptor, heparan sulfate, interacts with residues of VP1 in the putative receptor binding pocket, residues of VP3 in the adjacent pit and residues of the adjoining VP1 C-terminal loop which is recognised by the anti-TMEV capsid antibodies. These findings suggest that the anti-TMEV capsid antibodies neutralise virus infection by preventing heparan sulfate from binding to the capsid. The antibodies produced in this study are an important tool for further investigating virus-host cell interactions essential to picornavirus assembly.

Subcellular localisation of Theiler's murine encephalomyelitis virus (TMEV) capsid subunit VP1 vis-á-vis host protein Hsp90.

The VP1 subunit of the picornavirus capsid is the major antigenic determinant and mediates host cell attachment and virus entry. To investigate the localisation of Theiler's murine encephalomyelitis virus (TMEV) VP1 during infection, a bioinformatics approach was used to predict a surface-exposed, linear epitope region of the protein for subsequent expression and purification. This region, comprising the N-terminal 112 amino acids of the protein, was then used for rabbit immunisation, and the resultant polyclonal antibodies were able to recognise full length VP1 in infected cell lysates by Western blot. Following optimisation, the antibodies were used to investigate the localisation of VP1 in relation to Hsp90 in infected cells by indirect immunofluorescence and confocal microscopy. At 5h post infection, VP1 was distributed diffusely in the cytoplasm with strong perinuclear staining but was absent from the nucleus of all cells analysed. Dual-label immunofluorescence using anti-TMEV VP1 and anti-Hsp90 antibodies indicated that the distribution of both proteins colocalised in the cytoplasm and perinuclear region of infected cells. This is the first report describing the localisation of TMEV VP1 in infected cells, and the antibodies produced provide a valuable tool for investigating the poorly understood mechanisms underlying the early steps of picornavirus assembly.

JMS: An Open Source Workflow Management System and Web-Based Cluster Front-End for High Performance Computing.

Complex computational pipelines are becoming a staple of modern scientific research. Often these pipelines are resource intensive and require days of computing time. In such cases, it makes sense to run them over high performance computing (HPC) clusters where they can take advantage of the aggregated resources of many powerful computers. In addition to this, researchers often want to integrate their workflows into their own web servers. In these cases, software is needed to manage the submission of jobs from the web interface to the cluster and then return the results once the job has finished executing. We have developed the Job Management System (JMS), a workflow management system and web interface for high performance computing (HPC). JMS provides users with a user-friendly web interface for creating complex workflows with multiple stages. It integrates this workflow functionality with the resource manager, a tool that is used to control and manage batch jobs on HPC clusters. As such, JMS combines workflow management functionality with cluster administration functionality. In addition, JMS provides developer tools including a code editor and the ability to version tools and scripts. JMS can be used by researchers from any field to build and run complex computational pipelines and provides functionality to include these pipelines in external interfaces. JMS is currently being used to house a number of bioinformatics pipelines at the Research Unit in Bioinformatics (RUBi) at Rhodes University. JMS is an open-source project and is freely available at https://github.com/RUBi-ZA/JMS.

Sequence and domain conservation of the coelacanth Hsp40 and Hsp90 chaperones suggests conservation of function.

Molecular chaperones and their associated co-chaperones play an important role in preserving and regulating the active conformational state of cellular proteins. The chaperone complement of the Indonesian Coelacanth, Latimeria menadoensis, was elucidated using transcriptomic sequences. Heat shock protein 90 (Hsp90) and heat shock protein 40 (Hsp40) chaperones, and associated co-chaperones were focused on, and homologous human sequences were used to search the sequence databases. Coelacanth homologs of the cytosolic, mitochondrial and endoplasmic reticulum (ER) homologs of human Hsp90 were identified, as well as all of the major co-chaperones of the cytosolic isoform. Most of the human Hsp40s were found to have coelacanth homologs, and the data suggested that all of the chaperone machinery for protein folding at the ribosome, protein translocation to cellular compartments such as the ER and protein degradation were conserved. Some interesting similarities and differences were identified when interrogating human, mouse, and zebrafish homologs. For example, DnaJB13 is predicted to be a non-functional Hsp40 in humans, mouse, and zebrafish due to a corrupted histidine-proline-aspartic acid (HPD) motif, while the coelacanth homolog has an intact HPD. These and other comparisons enabled important functional and evolutionary questions to be posed for future experimental studies.

Plasmodium falciparum Hsp70-x: a heat shock protein at the host-parasite interface.

Plasmodium falciparum 70 kDa heat shock proteins (PfHsp70s) are expressed at all stages of the pathogenic erythrocytic phase of the malaria parasite life cycle. There are six PfHsp70s, all of which have orthologues in other plasmodial species, except for PfHsp70-x which is unique to P. falciparum. This research highlights a number of original results obtained by a detailed bioinformatics analysis of the protein. Large-scale sequence analysis indicated the presence of an extended transit peptide sequence of PfHsp70-x which potentially directs it to the endoplasmic reticulum (ER). Further analysis showed that PfHsp70-x does not have an ER-retention sequence, suggesting that the protein transits through the ER and is secreted into the parasitophorous vacuole or beyond into the erythrocyte cytosol. These results are consistent with experimental findings. Next, possible interactions between PfHsp70-x and exported P. falciparum Hsp40s or host erythrocyte Hsp40 were interrogated by modelling and docking. Docking results indicated that interaction between PfHsp70-x and each of the Hsp40s, regardless of biological feasibility, seems equally likely. This suggests that J domain might not provide the specificity in the formation of unique Hsp70-Hsp40 complexes, but that the specificity might be provided by other domains of Hsp40s. By studying different structural conformations of PfHsp70-x, it was shown that Hsp40s can only bind when PfHsp70-x is in a certain conformation. Additionally, this work highlighted the possible dependence of the substrate-binding domain residues on the orientation of the α-helical lid for formation of the substrate-binding pocket.

Analysis of protein thermostability enhancing factors in industrially important thermus bacteria species.

Elucidation of evolutionary factors that enhance protein thermostability is a critical problem and was the focus of this work on Thermus species. Pairs of orthologous sequences of T. scotoductus SA-01 and T. thermophilus HB27, with the largest negative minimum folding energy (MFE) as predicted by the UNAFold algorithm, were statistically analyzed. Favored substitutions of amino acids residues and their properties were determined. Substitutions were analyzed in modeled protein structures to determine their locations and contribution to energy differences using PyMOL and FoldX programs respectively. Dominant trends in amino acid substitutions consistent with differences in thermostability between orthologous sequences were observed. T. thermophilus thermophilic proteins showed an increase in non-polar, tiny, and charged amino acids. An abundance of alanine substituted by serine and threonine, as well as arginine substituted by glutamine and lysine was observed in T. thermophilus HB27. Structural comparison showed that stabilizing mutations occurred on surfaces and loops in protein structures.

The PINIT domain of PIAS3: structure-function analysis of its interaction with STAT3.

The protein inhibitor of activated signal transducer and activator of transcription 3 (PIAS3) regulates the transcriptional activity of signal transducer and activator of transcription 3 (STAT3) which regulates transcription of genes involved in cell growth, proliferation and apoptosis. The conserved proline, isoleucine, asparagine, isoleucine, threonine (PINIT) domain of PIAS3 is thought to promote STAT3-PIAS3 interaction. The (His)(7) -PINIT domain (PIAS3(85-272) ) was heterologously expressed and purified to homogeneity by nickel affinity and size exclusion chromatography, and shown to be a folded monomer in solution. Using surface plasmon resonance spectroscopy (SPR) the PINIT domain (PIAS3(85-272) ) alone was shown to specifically bind to STAT3 in a concentration dependent manner. L97A, R99N and R99Q mutations of the PINIT domain were found to abrogate binding to STAT3, suggesting that these residues were part of a potential binding surface. An homology model for the PINIT domain was calculated to analyse the potential locations of L97 and R99 in the structure, and to evaluate the potential role of these residues in interactions with STAT3.