Study of protein complexes via homology modeling, applied to cysteine proteases and their protein inhibitors

Investigation of Multi-Subunit Mycobacterium tuberculosis DNA-Directed RNA Polymerase and Its Rifampicin Resistant Mutants

Emerging Mycobacterium tuberculosis (Mtb) resistant strains have continued to limit the efficacies of existing antitubercular therapies. More specifically, mutations in the RNA replicative machinery of Mtb, RNA polymerase (RNAP), have been widely linked to rifampicin (RIF) resistance, which has led to therapeutic failures in many clinical cases. Moreover, elusive details on the underlying mechanisms of RIF-resistance caused by Mtb-RNAP mutations have hampered the development of new and efficient drugs that are able to overcome this challenge. Therefore, in this study we attempt to resolve the molecular and structural events associated with RIF-resistance in nine clinically reported missense Mtb RNAP mutations. Our study, for the first time, investigated the multi-subunit Mtb RNAP complex and findings revealed that the mutations commonly disrupted structural-dynamical attributes that may be essential for the protein's catalytic functions, particularly at the βfork loop 2, β'zinc-binding domain, the β' trigger loop and β'jaw, which in line with previous experimental reports, are essential for RNAP processivity. Complementarily, the mutations considerably perturbed the RIF-BP, which led to alterations in the active orientation of RIF needed to obstruct RNA extension. Consequentially, essential interactions with RIF were lost due to the mutation-induced repositioning with corresponding reductions in the binding affinity of the drug observed in majority of the mutants. We believe these findings will significantly aid future efforts in the discovery of new treatment options with the potential to overcome antitubercular resistance.

Comparing sequence and structure of falcipains and human homologs at prodomain and catalytic active site for malarial peptide based inhibitor design

BACKGROUND

Falcipains are major cysteine proteases of Plasmodium falciparum involved in haemoglobin degradation and remain attractive anti-malarial drug targets. Several inhibitors against these proteases have been identified, yet none of them has been approved for malaria treatment. Other Plasmodium species also possess highly homologous proteins to falcipains. For selective therapeutic targeting, identification of sequence and structure differences with homologous human cathepsins is necessary. The substrate processing activity of these proteins is tightly controlled via a prodomain segment occluding the active site which is chopped under low pH conditions exposing the catalytic site. Current work characterizes these proteases to identify residues mediating the prodomain regulatory function for the design of peptide based anti-malarial inhibitors.

METHODS

Sequence and structure variations between prodomain regions of plasmodial proteins and human cathepsins were determined using in silico approaches. Additionally, evolutionary clustering of these proteins was evaluated using phylogenetic analysis. High quality partial zymogen protein structures were modelled using homology modelling and residue interaction analysis performed between the prodomain segment and mature domain to identify key interacting residues between these two domains. The resulting information was used to determine short peptide sequences which could mimic the inherent regulatory function of the prodomain regions. Through flexible docking, the binding affinity of proposed peptides on the proteins studied was evaluated.

RESULTS

Sequence, evolutionary and motif analyses showed important differences between plasmodial and human proteins. Residue interaction analysis identified important residues crucial for maintaining prodomain integrity across the different proteins as well as the pro-segment responsible for inhibitory mechanism. Binding affinity of suggested peptides was highly dependent on their residue composition and length.

CONCLUSIONS

Despite the conserved structural and catalytic mechanism between human cathepsins and plasmodial proteases, current work revealed significant differences between the two protein groups which may provide valuable information for selective anti-malarial inhibitor development. Part of this study aimed to design peptide inhibitors based on endogenous inhibitory portions of protease prodomains as a novel aspect. Even though peptide inhibitors may not be practical solutions to malaria at this stage, the approach followed and results offer a promising means to find new malarial inhibitors.

In Silico Analysis of Homologous Heterodimers of Cruzipain-Chagasin from Structural Models Built by Homology

The present study gives an overview of the binding energetics of the homologous heterodimers of cruzipain−chagasin based on the binding energy (ΔG_b) prediction obtained with FoldX. This analysis involves a total of 70 homologous models of the cruzipain−chagasin complex which were constructed by homology from the combinatory variation of nine papain-like cysteine peptidase structures and seven cysteine protease inhibitor structures (as chagasin-like and cystatin-like inhibitors). Only 32 systems have been evaluated experimentally, ΔG_b^experimental values previously reported. Therefore, the result of the multiple analysis in terms of the thermodynamic parameters, are shown as relative energy |ΔΔG| = |ΔG_b^fromFoldX − ΔG_b^experimental|. Nine models were identified that recorded |ΔΔG| < 1.3, five models to 2.8 > |ΔΔG| > 1.3 and the other 18 models, values of |ΔΔG| > 2.8. The energetic analysis of the contribution of ΔH and ΔS to ΔG_b to the 14-molecular model presents a ΔG_b mostly ΔH-driven at neutral pH and at an ionic strength (I) of 0.15 M. The dependence of ΔG_b(I,pH) at 298 K to the cruzipain−chagasin complex predicts a linear dependence of ΔG_b(I). The computational protocol allowed the identification and prediction of thermodynamics binding energy parameters for cruzipain−chagasin-like heterodimers.

Aminoacyl tRNA synthetases as malarial drug targets: a comparative bioinformatics study

BACKGROUND

Treatment of parasitic diseases has been challenging due to evolution of drug resistant parasites, and thus there is need to identify new class of drugs and drug targets. Protein translation is important for survival of malarial parasite, Plasmodium, and the pathway is present in all of its life cycle stages. Aminoacyl tRNA synthetases are primary enzymes in protein translation as they catalyse amino acid addition to the cognate tRNA. This study sought to understand differences between Plasmodium and human aminoacyl tRNA synthetases through bioinformatics analysis.

METHODS

Plasmodium berghei, Plasmodium falciparum, Plasmodium fragile, Plasmodium knowlesi, Plasmodium malariae, Plasmodium ovale, Plasmodium vivax, Plasmodium yoelii and human aminoacyl tRNA synthetase sequences were retrieved from UniProt database and grouped into 20 families based on amino acid specificity. These families were further divided into two classes. Both families and classes were analysed. Motif discovery was carried out using the MEME software, sequence identity calculation was done using an in-house Python script, multiple sequence alignments were performed using PROMALS3D and TCOFFEE tools, and phylogenetic tree calculations were performed using MEGA vs 7.0 tool. Possible alternative binding sites were predicted using FTMap webserver and SiteMap tool.

RESULTS

Motif discovery revealed Plasmodium-specific motifs while phylogenetic tree calculations showed that Plasmodium proteins have different evolutionary history to the human homologues. Human aaRSs sequences showed low sequence identity (below 40%) compared to Plasmodium sequences. Prediction of alternative binding sites revealed potential druggable sites in PfArgRS, PfMetRS and PfProRS at regions that are weakly conserved when compared to the human homologues. Multiple sequence analysis, motif discovery, pairwise sequence identity calculations and phylogenetic tree analysis showed significant differences between parasite and human aaRSs proteins despite functional and structural conservation. These differences may provide a basis for further exploration of Plasmodium aminoacyl tRNA synthetases as potential drug targets.

CONCLUSION

This study showed that, despite, functional and structural conservation, Plasmodium aaRSs have key differences from the human homologues. These differences in Plasmodium aaRSs can be targeted to develop anti-malarial drugs with less toxicity to the host.

PRIMO: An Interactive Homology Modeling Pipeline

The development of automated servers to predict the three-dimensional structure of proteins has seen much progress over the years. These servers make calculations simpler, but largely exclude users from the process. In this study, we present the PRotein Interactive MOdeling (PRIMO) pipeline for homology modeling of protein monomers. The pipeline eases the multi-step modeling process, and reduces the workload required by the user, while still allowing engagement from the user during every step. Default parameters are given for each step, which can either be modified or supplemented with additional external input. PRIMO has been designed for users of varying levels of experience with homology modeling. The pipeline incorporates a user-friendly interface that makes it easy to alter parameters used during modeling. During each stage of the modeling process, the site provides suggestions for novice users to improve the quality of their models. PRIMO provides functionality that allows users to also model ligands and ions in complex with their protein targets. Herein, we assess the accuracy of the fully automated capabilities of the server, including a comparative analysis of the available alignment programs, as well as of the refinement levels used during modeling. The tests presented here demonstrate the reliability of the PRIMO server when producing a large number of protein models. While PRIMO does focus on user involvement in the homology modeling process, the results indicate that in the presence of suitable templates, good quality models can be produced even without user intervention. This gives an idea of the base level accuracy of PRIMO, which users can improve upon by adjusting parameters in their modeling runs. The accuracy of PRIMO’s automated scripts is being continuously evaluated by the CAMEO (Continuous Automated Model EvaluatiOn) project. The PRIMO site is free for non-commercial use and can be accessed at https://primo.rubi.ru.ac.za/.

Targeted expression of cystatin restores fertility in cysteine protease induced male sterile tobacco plants

Fertility restoration in male sterile plants is an essential requirement for their utilization in hybrid seed production. In an earlier investigation, we have demonstrated that the targeted expression of a cysteine protease in tapetal cell layer resulted in complete male sterility in tobacco transgenic plants. In the present investigation, we have used a cystatin gene, which encodes for a cysteine protease inhibitor, from a wild peanut, Arachis diogoi and developed a plant gene based restoration system for cysteine protease induced male sterile transgenic tobacco plants. We confirmed the interaction between the cysteine protease and a cystatin of the wild peanut, A. diogoi through in silico modeling and yeast two-hybrid assay. Pollen from primary transgenic tobacco plants expressing cystatin gene under the tapetum specific promoter- TA29 restored fertility on cysteine protease induced male sterile tobacco plants developed earlier. This has confirmed the in vivo interaction of cysteine protease and cystatin in the tapetal cells, and the inactivation of cysteine protease and modulation of its negative effects on pollen fertility. Both the cysteine protease and cystatin genes are of plant origin in contrast to the analogous barnase-barstar system that deploys genes of prokaryotic origin. Because of the deployment of genes of plant origin, this system might not face biosafety problems in developing hybrids in food crops.

Plasmodium falciparum Hop: detailed analysis on complex formation with Hsp70 and Hsp90

The heat shock organizing protein (Hop) is important in modulating the activity and co-interaction of two chaperones: heat shock protein 70 and 90 (Hsp70 and Hsp90). Recent research suggested that Plasmodium falciparum Hop (PfHop), PfHsp70 and PfHsp90 form a complex in the trophozoite infective stage. However, there has been little computational research on the malarial Hop protein in complex with other malarial Hsps. Using in silico characterization of the protein, this work showed that individual domains of Hop are evolving at different rates within the protein. Differences between human Hop (HsHop) and PfHop were identified by motif analysis. Homology modeling of PfHop and HsHop in complex with their own cytosolic Hsp90 and Hsp70 C-terminal peptide partners indicated excellent conservation of the Hop concave TPR sites bound to the C-terminal motifs of partner proteins. Further, we analyzed additional binding sites between Hop and Hsp90, and showed, for the first time, that they are distinctly less conserved between human and malaria parasite. These sites are located on the convex surface of Hop TPR2, and involved in interactions with the Hsp90 middle domain. Since the convex sites are less conserved than the concave sites, it makes their potential for malarial inhibitor design extremely attractive (as opposed to the concave sites which have been the focus of previous efforts).

C1A cysteine-proteases and their inhibitors in plants

Plant cysteine-proteases (CysProt) represent a well-characterized type of proteolytic enzymes that fulfill tightly regulated physiological functions (senescence and seed germination among others) and defense roles. This article is focused on the group of papain-proteases C1A (family C1, clan CA) and their inhibitors, phytocystatins (PhyCys). In particular, the protease-inhibitor interaction and their mutual participation in specific pathways throughout the plant's life are reviewed. C1A CysProt and PhyCys have been molecularly characterized, and comparative sequence analyses have identified consensus functional motifs. A correlation can be established between the number of identified CysProt and PhyCys in angiosperms. Thus, evolutionary forces may have determined a control role of cystatins on both endogenous and pest-exogenous proteases in these species. Tagging the proteases and inhibitors with fluorescence proteins revealed common patterns of subcellular localization in the endoplasmic reticulum-Golgi network in transiently transformed onion epidermal cells. Further in vivo interactions were demonstrated by bimolecular fluorescent complementation, suggesting their participation in the same physiological processes.

Protease-associated cellular networks in malaria parasite Plasmodium falciparum

Background

Malaria continues to be one of the most severe global infectious diseases, responsible for 1-2 million deaths yearly. The rapid evolution and spread of drug resistance in parasites has led to an urgent need for the development of novel antimalarial targets. Proteases are a group of enzymes that play essential roles in parasite growth and invasion. The possibility of designing specific inhibitors for proteases makes them promising drug targets. Previously, combining a comparative genomics approach and a machine learning approach, we identified the complement of proteases (degradome) in the malaria parasite Plasmodium falciparum and its sibling species [1-3], providing a catalog of targets for functional characterization and rational inhibitor design. Network analysis represents another route to revealing the role of proteins in the biology of parasites and we use this approach here to expand our understanding of the systems involving the proteases of P. falciparum.

Results

We investigated the roles of proteases in the parasite life cycle by constructing a network using protein-protein association data from the STRING database [4], and analyzing these data, in conjunction with the data from protein-protein interaction assays using the yeast 2-hybrid (Y2H) system [5], blood stage microarray experiments [6-8], proteomics [9-12], literature text mining, and sequence homology analysis. Seventy-seven (77) out of 124 predicted proteases were associated with at least one other protein, constituting 2,431 protein-protein interactions (PPIs). These proteases appear to play diverse roles in metabolism, cell cycle regulation, invasion and infection. Their degrees of connectivity (i.e., connections to other proteins), range from one to 143. The largest protease-associated sub-network is the ubiquitin-proteasome system which is crucial for protein recycling and stress response. Proteases are also implicated in heat shock response, signal peptide processing, cell cycle progression, transcriptional regulation, and signal transduction networks.

Conclusions

Our network analysis of proteases from P. falciparum uses a so-called guilt-by-association approach to extract sets of proteins from the proteome that are candidates for further study. Novel protease targets and previously unrecognized members of the protease-associated sub-systems provide new insights into the mechanisms underlying parasitism, pathogenesis and virulence.