MD-TASK: a software suite for analyzing molecular dynamics trajectories.

SUMMARY: Molecular dynamics (MD) determines the physical motions of atoms of a biological macromolecule in a cell-like environment and is an important method in structural bioinformatics. Traditionally, measurements such as root mean square deviation, root mean square fluctuation, radius of gyration, and various energy measures have been used to analyze MD simulations. Here, we present MD-TASK, a novel software suite that employs graph theory techniques, perturbation response scanning, and dynamic cross-correlation to provide unique ways for analyzing MD trajectories. AVAILABILITY AND IMPLEMENTATION: MD-TASK has been open-sourced and is available for download from https://github.com/RUBi-ZA/MD-TASK , implemented in Python and supported on Linux/Unix. CONTACT: o.tastanbishop@ru.ac.za.

Allosteric Modulation of Human Hsp90α Conformational Dynamics.

Central to Hsp90's biological function is its ability to interconvert between various conformational states. Drug targeting of Hsp90's regulatory mechanisms, including its modulation by cochaperone association, presents as an attractive therapeutic strategy for Hsp90 associated pathologies. In this study, we utilized homology modeling techniques to calculate full-length structures of human Hsp90α in closed and partially open conformations and used these structures as a basis for several molecular dynamics based analyses aimed at elucidating allosteric mechanisms and modulation sites in human Hsp90α. Atomistic simulations demonstrated that bound adenosine triphosphate (ATP) stabilizes the dimer by "tensing" each protomer, while adenosine diphosphate (ADP) and apo configurations "relax" the complex by increasing global flexibility, the former case resulting in a fully open "v-like" conformation. Dynamic residue network analysis revealed regions of the protein involved in intraprotein communication and identified several key communication hubs that correlate with known functional sites. Pairwise comparison of betweenness centrality, shortest path, and residue fluctuations revealed that a proportional relationship exists between the latter two measurables and an inverse relationship between these two and betweenness centrality. This analysis showed how protein flexibility, degree of compactness, and the distance cutoff used for network construction influence the correlations between these metrics. These findings are novel and suggest shortest path and betweenness centrality to be more relevant quantities to follow for detecting functional residues in proteins compared to residue fluctuations. Perturbation response scanning analysis identified several potential residue sites capable of modulating conformational change in favor of interstate conversion. For the ATP-bound open conformation, these sites were found to overlap with known Aha1 and client binding sites, demonstrating how naturally occurring forces associated with cofactor binding could allosterically modulate conformational dynamics.

Perturbation-Response Scanning Reveals Key Residues for Allosteric Control in Hsp70.

Hsp70 molecular chaperones play an important role in maintaining cellular homeostasis, and are implicated in a wide array of cellular processes, including protein recovery from aggregates, cross-membrane protein translocation, and protein biogenesis. Hsp70 consists of two domains, a nucleotide binding domain (NBD) and a substrate binding domain (SBD), each of which communicates via an allosteric mechanism such that the protein interconverts between two functional states, an ATP-bound open conformation and an ADP-bound closed conformation. The exact mechanism for interstate conversion is not as yet fully understood. However, the ligand-bound states of the NBD and SBD as well as interactions with cochaperones such as DnaJ and nucleotide exchange factor are thought to play crucial regulatory roles. In this study, we apply the perturbation-response scanning (PRS) method in combination with molecular dynamics simulations as a computational tool for the identification of allosteric hot residues in the large multidomain Hsp70 protein. We find evidence in support of the hypothesis that substrate binding triggers ATP hydrolysis and that the ADP-substrate complex favors interstate conversion to the closed state. Furthermore, our data are in agreement with the proposal that there is an allosterically active intermediate state between the open and closed states and vice versa, as we find evidence that ATP binding to the closed structure and peptide binding to the open structure allosterically "activate" the respective complexes. We conclude our analysis by showing how our PRS data fit the current opinion on the Hsp70 conformational cycle and present several allosteric hot residues that may provide a platform for further studies to gain additional insight into Hsp70 allostery.

Crosstalk analysis of ring resonator switches for all-optical routing.

Optical switches based on ring resonator cavities were fabricated by a silicon photonics foundry process and analyzed for optical crosstalk at various data rates and channel spacings. These devices were compared to commercial bandpass filters and at 20Gb/s, 0.5dB power penalty is observed due to spectral filtering for bit error ratio threshold of 1 × 10-9. Concurrent modulation at 20Gb/s with a channel spacing as narrow as 40GHz shows error-free transmission with 1dB power penalty as compared to wider channel spacing for the ring-based switch.

Modulation of Human Hsp90α Conformational Dynamics by Allosteric Ligand Interaction at the C-Terminal Domain.

Recent years have seen heat shock protein 90 kDa (Hsp90) attract significant interest as a viable drug target, particularly for cancer. To date, designed inhibitors that target the ATPase domain demonstrate potent anti-proliferative effects, but have failed clinical trials due to high levels of associated toxicity. To circumvent this, the focus has shifted away from the ATPase domain. One option involves modulation of the protein through allosteric activation/inhibition. Here, we propose a novel approach: we use previously obtained information via residue perturbation scanning coupled with dynamic residue network analysis to identify allosteric drug targeting sites for inhibitor docking. We probe the open conformation of human Hsp90α for druggable sites that overlap with these allosteric control elements, and identify three putative natural compound allosteric modulators: Cephalostatin 17, 20(29)-Lupene-3ß-isoferulate and 3'-Bromorubrolide F. We assess the allosteric potential of these ligands by examining their effect on the conformational dynamics of the protein. We find evidence for the selective allosteric activation and inhibition of Hsp90's conformational transition toward the closed state in response to ligand binding and shed valuable insight to further the understanding of allosteric drug design and Hsp90's complex allosteric mechanism of action.

Binding and entry of a non-enveloped T=4 insect RNA virus is triggered by alkaline pH.

Tetraviruses are small, non-enveloped, RNA viruses that exclusively infect lepidopteran insects. Their particles comprise 240 copies of a single capsid protein precursor (CP), which undergoes autoproteolytic cleavage during maturation. The molecular mechanisms of capsid assembly and maturation are well understood, but little is known about the viral infectious lifecycle due to a lack of tissue culture cell lines that are susceptible to tetravirus infection. We show here that binding and entry of the alphatetravirus, Helicoverpa armigera stunt virus (HaSV), is triggered by alkaline pH. At pH 9.0, wild-type HaSV virus particles undergo conformational changes that induce membrane-lytic activity and binding to Spodoptera frugiperda Sf9 cells. Binding is followed by entry and infection, with virus replication complexes detected by immunofluorescence microscopy within 2h post-infection and the CP after 12h. HaSV particles produced in S. frugiperda Sf9 cells are infectious. Helicoverpa armigera larval virus biofeed assays showed that pre-treatment with the V-ATPase inhibitor, Bafilomycin A1, resulted in a 50% decrease in larval mortality and stunting, while incubation of virus particles at pH 9.0 prior to infection restored infectivity. Together, these data show that HaSV, and likely other tetraviruses, requires the alkaline environment of the lepidopteran larval midgut for binding and entry into host cells.

Human, vector and parasite Hsp90 proteins: A comparative bioinformatics analysis.

The treatment of protozoan parasitic diseases is challenging, and thus identification and analysis of new drug targets is important. Parasites survive within host organisms, and some need intermediate hosts to complete their life cycle. Changing host environment puts stress on parasites, and often adaptation is accompanied by the expression of large amounts of heat shock proteins (Hsps). Among Hsps, Hsp90 proteins play an important role in stress environments. Yet, there has been little computational research on Hsp90 proteins to analyze them comparatively as potential parasitic drug targets. Here, an attempt was made to gain detailed insights into the differences between host, vector and parasitic Hsp90 proteins by large-scale bioinformatics analysis. A total of 104 Hsp90 sequences were divided into three groups based on their cellular localizations; namely cytosolic, mitochondrial and endoplasmic reticulum (ER). Further, the parasitic proteins were divided according to the type of parasite (protozoa, helminth and ectoparasite). Primary sequence analysis, phylogenetic tree calculations, motif analysis and physicochemical properties of Hsp90 proteins suggested that despite the overall structural conservation of these proteins, parasitic Hsp90 proteins have unique features which differentiate them from human ones, thus encouraging the idea that protozoan Hsp90 proteins should be further analyzed as potential drug targets.

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.

SANCDB: a South African natural compound database.

BACKGROUND: Natural products (NPs) are important to the drug discovery process. NP research efforts are expanding world-wide and South Africa is no exception to this. While freely-accessible small molecule databases, containing compounds isolated from indigenous sources, have been established in a number of other countries, there is currently no such online database in South Africa. DESCRIPTION: The current research presents a South African natural compound database, named SANCDB. This is a curated and fully-referenced database containing compound information for 600 natural products extracted directly from journal articles, book chapters and theses. There is a web interface to the database, which is simple and easy to use, while allowing for compounds to be searched by a number of different criteria. Being fully referenced, each compound page contains links to the original referenced work from which the information was obtained. Further, the website provides a submission pipeline, allowing researchers to deposit compounds from their own research into the database. CONCLUSIONS: SANCDB is currently the only web-based NP database in Africa. It aims to provide a useful resource for the in silico screening of South African NPs for drug discovery purposes. The database is supported by a submission pipeline to allow growth by entries from researchers. As such, we currently present SANCDB the starting point of a platform for a community-driven, curated database to further natural products research in South Africa. SANCDB is freely available at https://sancdb.rubi.ru.ac.za/.