Comparative modeling of thioredoxin glutathione reductase from Schistosoma mansoni: a multifunctional target for antischistosomal therapy

Analysis of Pleurotin binding to human thioredoxin reductase using docking and molecular dynamics simulation

Thioredoxin reductase (TrxR) has been considered a potential target for cancer chemotherapy. It acts by controlling the redox homeostasis of human cells and, therefore, interfering in its function may trigger apoptosis, which is a crucial tumor suppression mechanism. Despite the great effort in the search for TrxR inhibitors, none was approved for human therapy. In the present study a virtual screening for natural organic compounds is discussed for a set of 72 compounds with known IC-50 for TrxR inhibition. The results suggest the Pleurotin, a naphthoquinone obtained from Hohenbuehelia grisea fungus, as a potential TrxR inhibitor, which acts by binding to the active site of the enzyme, between the N- and C-terminal domains. The presence of the ligand blocks the approximation of the C-terminal arm to the N-terminal, which is an essential step of the enzyme function. Besides, the two equivalent binding sites of TrxR were explored, by docking two ligands simultaneously. The results indicate that both sites have an allosteric correlation and, the presence of the ligand in one site may interfere, or even prevent, the binding of the second ligand at the other site. All these findings are quantitatively discussed based on the analysis of long molecular dynamics trajectories, which provides a full description of the ligand-receptor binding modes, average binding energies and conformational changes.Communicated by Ramaswamy H. Sarma.

Insight into the Mechanistic Basis of the Hysteretic-Like Kinetic Behavior of Thioredoxin-Glutathione Reductase (TGR)

A kinetic study of thioredoxin-glutathione reductase (TGR) from Taenia crassiceps metacestode (cysticerci) was carried out. The results obtained from both initial velocity and product inhibition experiments suggest the enzyme follows a two-site ping-pong bi bi kinetic mechanism, in which both substrates and products are bound in rapid equilibrium fashion. The substrate GSSG exerts inhibition at moderate or high concentrations, which is concomitant with the observation of hysteretic-like progress curves. The effect of NADPH on the apparent hysteretic behavior of TGR was also studied. At low concentrations of NADPH in the presence of moderate concentrations of GSSG, atypical time progress curves were observed, consisting of an initial burst-like stage, followed by a lag whose amplitude and duration depended on the concentration of both NADPH and GSSG. Based on all the kinetic and structural evidence available on TGR, a mechanism-based model was developed. The model assumes a noncompetitive mode of inhibition by GSSG in which the disulfide behaves as an affinity label-like reagent through its binding and reduction at an alternative site, leading the enzyme into an inactive state. The critical points of the model are the persistence of residual GSSG reductase activity in the inhibited GSSG-enzyme complexes and the regeneration of the active form of the enzyme by GSH. Hence, the hysteretic-like progress curves of GSSG reduction by TGR are the result of a continuous competition between GSH and GSSG for driving the enzyme into active or inactive states, respectively. By using an arbitrary but consistent set of rate constants, the experimental full progress curves were successfully reproduced in silico.

Structural insights into natural compounds as inhibitors of Fasciola gigantica thioredoxin glutathione reductase

Fascioliasis is caused by the helminth parasites of genus Fasciola. Thioredoxin glutathione reductase (TGR) is an important enzyme in parasitic helminths and plays an indispensable role in its redox biology. In the present study, we conducted a structure-based virtual screening of natural compounds against the Fasciola gigantica TGR (FgTGR). The compounds were docked against FgTGR in four sequential docking modes. The screened ligands were further assessed for Lipinski and ADMET prediction so as to evaluate drug proficiency and likeness property. After refinement, three potential inhibitors were identified that were subjected to 50 ns molecular dynamics simulation and free energy binding analyses to evaluate the dynamics of protein-ligand interaction and the stability of the complexes. Key residues involved in the interaction of the selected ligands were also determined. The results suggested that three top hits had a negative binding energy greater than GSSG (-91.479 KJ · mol^-1 ), having -152.657, -141.219, and -92.931 kJ · mol^-1 for compounds with IDs ZINC85878789, ZINC85879991, and ZINC36369921, respectively. Further analysis showed that the compound ZINC85878789 and ZINC85879991 displayed substantial pharmacological and structural properties to be a drug candidate. Thus, the present study might prove useful for the future design of new derivatives with higher potency and specificity.

Proteomic Analysis of the Schistosoma mansoni Miracidium

Despite extensive control efforts, schistosomiasis continues to be a major public health problem in developing nations in the tropics and sub-tropics. The miracidium, along with the cercaria, both of which are water-borne and free-living, are the only two stages in the life-cycle of Schistosoma mansoni which are involved in host invasion. Miracidia penetrate intermediate host snails and develop into sporocysts, which lead to cercariae that can infect humans. Infection of the snail host by the miracidium represents an ideal point at which to interrupt the parasite’s life-cycle. This research focuses on an analysis of the miracidium proteome, including those proteins that are secreted. We have identified a repertoire of proteins in the S. mansoni miracidium at 2 hours post-hatch, including proteases, venom allergen-like proteins, receptors and HSP70, which might play roles in snail-parasite interplay. Proteins involved in energy production and conservation were prevalent, as were proteins predicted to be associated with defence. This study also provides a strong foundation for further understanding the roles that neurohormones play in host-seeking by schistosomes, with the potential for development of novel anthelmintics that interfere with its various life-cycle stages.

Molecular docking to explore the possible binding mode of potential inhibitors of thioredoxin glutathione reductase

Praziquantel (PZQ) is the treatment of choice for schistosomiasis, one of the most important but neglected tropical diseases. Recently, however, Schistosoma have exhibited reduced susceptibility to PZQ, and an urgent need to develop new drugs to treat schistosomiasis has emerged. Thioredoxin glutathione reductase (TGR) plays a crucial role in the redox balance of the parasite, combining glutaredoxin (Grx), glutathione reductase and thioredoxin reductase (TR) activities. Several compounds, including oxadiazole 2-oxides, phosphinic acid amides, isoxazolones and phosphoramidites, have been identified as agents that inhibit TGR from Schistosoma mansoni (smTGR) and exhibit anti-schistosomal activity. 4-Phenyl-1,2,5-oxadiazole-3-carbonitrile-2-oxide has also been shown to be active against TGR from Schistosoma japonicum (sjTGR). The binding sites of these inhibitors, however, remain unclear. To explore the binding interactions of these compounds, we selected six compounds to dock into the NADPH binding site, the active site of the TR domain and the Grx active site of both smTGR and sjTGR using AutoDock 4.2.5.1. The results suggested that the most favoured binding site for all compounds in either sjTGR or smTGR was the oxidised glutathione-binding pocket of the TR domain. Although all of the compounds could fit into the sjTGR site, the inhibition efficiency of these compounds towards sjTGR was marginally lower than it was towards smTGR, suggesting that it would be necessary to design specific inhibitors of TGR for different Schistosoma species. The docking results showed that all compounds docking in smTGR and sjTGR adopted similar binding modes in the TR domain. Two peptide fragments from another subunit, Phe505′–Leu508′ and Pro572′–Thr577′, played a critical role in the interactions with the inhibitors. In conclusion, the present study has revealed binding mechanisms for potential inhibitors of Schistosoma TGRs and could lead to structure-based ligand design and the development of new anti-schistosomiasis drugs.

Charge group partitioning in biomolecular simulation

Molecular simulation techniques are increasingly being used to study biomolecular systems at an atomic level. Such simulations rely on empirical force fields to represent the intermolecular interactions. There are many different force fields available--each based on a different set of assumptions and thus requiring different parametrization procedures. Recently, efforts have been made to fully automate the assignment of force-field parameters, including atomic partial charges, for novel molecules. In this work, we focus on a problem arising in the automated parametrization of molecules for use in combination with the GROMOS family of force fields: namely, the assignment of atoms to charge groups such that for every charge group the sum of the partial charges is ideally equal to its formal charge. In addition, charge groups are required to have size at most k. We show NP-hardness and give an exact algorithm that solves practical problem instances to provable optimality in a fraction of a second.

Molecular dynamics simulation of the interaction between protein tyrosine phosphatase 1B and aryl diketoacid derivatives

The protein tyrosine phosphatase 1B (PTP-1B) is acknowledged as an outstanding therapeutic target for the treatment of diabetes, obesity and cancer. In this work, six aryl diketoacid compounds have been studied on the basis of molecular dynamics simulations. Hydrogen bonds, binding energies and conformation changes of the WPD loop have been analyzed. The results indicated that their activation model falls into two parts: the target region of the monomeric aryl diketoacid compounds is the active site, whereas the target region of the dimeric aryl diketoacid compounds is the WPD loop or the R loop. The van der Waals interactions exhibit stronger effects than the short-range electrostatic interactions. The van der Waals interaction energy and the IC50 values exhibit an approximately exponential relationship. Furthermore, the van der Waals interactions cooperate with the hydrogen bond interactions. This study provides a more thorough understanding of the PTP-1B inhibitor binding processes.

Thioredoxin glutathione reductase: its role in redox biology and potential as a target for drugs against neglected diseases

BACKGROUND

There are two, largely autonomous antioxidant pathways in many organisms, one based on thioredoxin and one based on glutathione, with each pathway having a unique flavoprotein oxidoreductase to maintain them in a reduced state. A recently discovered protein, thioredoxin glutathione reductase (TGR) potentially connects these two pathways. In a large group of parasitic worms, responsible for hundreds of millions of infections in humans and animals, untold morbidity and significant mortality, TGR is the sole enzyme present to maintain redox balance.

SCOPE OF REVIEW

In this review, the current understanding of the biochemical properties of TGR enzymes is compared to the related enzymes thioredoxin reductase and glutathione reductase. The role of the rare amino acid selenocysteine is discussed. An overview of the potential to target TGR for drug development against a range of parasitic worms and preliminary results to identify TGR inhibitors for schistosomiasis treatment is presented.

MAJOR CONCLUSIONS

TGR has properties that are both unique and common to other flavoprotein oxidoreductases. TGR plays a fundamentally different and essential role in the redox biology of parasitic flatworms. Therefore, TGR is a promising target for drug development for schistosomiasis and other trematode and cestode infections.

GENERAL SIGNIFICANCE

TGR may have differing functions in host organisms, but through analyses to understand its ability to reduce both glutathione and thioredoxin we can better understand the reaction mechanisms of an important class of enzymes. The unique properties of TGR in parasitic flatworms provide promising routes to develop new treatments for diseases.

HomPPI: a class of sequence homology based protein-protein interface prediction methods

BACKGROUND

Although homology-based methods are among the most widely used methods for predicting the structure and function of proteins, the question as to whether interface sequence conservation can be effectively exploited in predicting protein-protein interfaces has been a subject of debate.

RESULTS

We studied more than 300,000 pair-wise alignments of protein sequences from structurally characterized protein complexes, including both obligate and transient complexes. We identified sequence similarity criteria required for accurate homology-based inference of interface residues in a query protein sequence.Based on these analyses, we developed HomPPI, a class of sequence homology-based methods for predicting protein-protein interface residues. We present two variants of HomPPI: (i) NPS-HomPPI (Non partner-specific HomPPI), which can be used to predict interface residues of a query protein in the absence of knowledge of the interaction partner; and (ii) PS-HomPPI (Partner-specific HomPPI), which can be used to predict the interface residues of a query protein with a specific target protein.Our experiments on a benchmark dataset of obligate homodimeric complexes show that NPS-HomPPI can reliably predict protein-protein interface residues in a given protein, with an average correlation coefficient (CC) of 0.76, sensitivity of 0.83, and specificity of 0.78, when sequence homologs of the query protein can be reliably identified. NPS-HomPPI also reliably predicts the interface residues of intrinsically disordered proteins. Our experiments suggest that NPS-HomPPI is competitive with several state-of-the-art interface prediction servers including those that exploit the structure of the query proteins. The partner-specific classifier, PS-HomPPI can, on a large dataset of transient complexes, predict the interface residues of a query protein with a specific target, with a CC of 0.65, sensitivity of 0.69, and specificity of 0.70, when homologs of both the query and the target can be reliably identified. The HomPPI web server is available at http://homppi.cs.iastate.edu/.

CONCLUSIONS

Sequence homology-based methods offer a class of computationally efficient and reliable approaches for predicting the protein-protein interface residues that participate in either obligate or transient interactions. For query proteins involved in transient interactions, the reliability of interface residue prediction can be improved by exploiting knowledge of putative interaction partners.

Investigations of the catalytic mechanism of thioredoxin glutathione reductase from Schistosoma mansoni

Thioredoxin glutathione reductase from Schistosoma mansoni (SmTGR) catalyzes the reduction of both thioredoxin and glutathione disulfides (GSSG), thus playing a crucial role in maintaining redox homeostasis in the parasite. In line with this role, previous studies have demonstrated that SmTGR is a promising drug target for schistosomiasis. To aid in the development of efficacious drugs that target SmTGR, it is essential to understand the catalytic mechanism of SmTGR. SmTGR is a dimeric flavoprotein in the glutathione reductase family and has a head-to-tail arrangement of its monomers; each subunit has the components of both a thioredoxin reductase (TrxR) domain and a glutaredoxin (Grx) domain. However, the active site of the TrxR domain is composed of residues from both subunits: FAD and a redox-active Cys-154/Cys-159 pair from one subunit and a redox-active Cys-596'/Sec-597' pair from the other; the active site of the Grx domain contains a redox-active Cys-28/Cys-31 pair. Via its Cys-28/Cys-31 dithiol and/or its Cys-596'/Sec-597' thiol-selenolate, SmTGR can catalyze the reduction of a variety of substrates by NADPH. It is presumed that SmTGR catalyzes deglutathionylation reactions via the Cys-28/Cys-31 dithiol. Our anaerobic titration data suggest that reducing equivalents from NADPH can indeed reach the Cys-28/Cys-31 disulfide in the Grx domain to facilitate reductions effected by this cysteine pair. To clarify the specific chemical roles of each redox-active residue with respect to its various reactivities, we generated variants of SmTGR. Cys-28 variants had no Grx deglutathionylation activity, whereas Cys-31 variants retained partial Grx deglutathionylation activity, indicating that the Cys-28 thiolate is the nucleophile initiating deglutathionylation. Lags in the steady-state kinetics, found when wild-type SmTGR was incubated at high concentrations of GSSG, were not present in Grx variants, indicating that this cysteine pair is in some way responsible for the lags. A Sec-597 variant was still able to reduce a variety of substrates, albeit slowly, showing that selenocysteine is important but is not the sole determinant for the broad substrate tolerance of the enzyme. Our data show that Cys-520 and Cys-574 are not likely to be involved in the catalytic mechanism.

Practical considerations for building GROMOS-compatible small-molecule topologies

Molecular dynamics simulations are being applied to increasingly complex systems, including those involving small endogenous compounds and drug molecules. In order to obtain meaningful and accurate data from these simulations, high-quality topologies for small molecules must be generated in a manner that is consistent with the derivation of the force field applied to the system. Often, force fields are designed for use with macromolecules such as proteins, making their transferability to other species challenging. Investigators are increasingly attracted to automated topology generation programs, although the quality of the resulting topologies remains unknown. Here we assess the applicability of the popular PRODRG server that generates small-molecule topologies for use with the GROMOS family of force fields. We find that PRODRG does not reproduce topologies for even the most well-characterized species in the force field due to inconsistent charges and charge groups. We assessed the effects of PRODRG-derived charges on several systems: pure liquids, amino acids at a hydrophobic-hydrophilic interface, and an enzyme-cofactor complex. We found that partial atomic charges generated by PRODRG are largely incompatible with GROMOS force fields, and the behavior of these systems deviates substantially from that of simulations using GROMOS parameters. We conclude by proposing several points as "best practices" for parametrization of small molecules under the GROMOS force fields.

Schistosomiasis in the People's Republic of China: the era of the Three Gorges Dam

The potential impact of the Three Gorges Dam (TGD) on schistosomiasis transmission in China has invoked considerable global concern. The TGD will result in changes in the water level and silt deposition downstream, favoring the reproduction of Oncomelania snails. Combined with blockages of the Yangtze River's tributaries, these changes will increase the schistosomiasis transmission season within the marshlands along the middle and lower reaches of the Yangtze River. The changing schistosome transmission dynamics necessitate a comprehensive strategy to control schistosomiasis. This review discusses aspects of the epidemiology and transmission of Schistosoma japonicum in China and considers the pathology, clinical outcomes, diagnosis, treatment, immunobiology, and genetics of schistosomiasis japonica together with an overview of current progress in vaccine development, all of which will have an impact on future control efforts. The use of synchronous praziquantel (PZQ) chemotherapy for humans and domestic animals is only temporarily effective, as schistosome reinfection occurs rapidly. Drug delivery requires a substantial infrastructure to regularly cover all parts of an area of endemicity. This makes chemotherapy expensive and, as compliance is often low, a less than satisfactory control option. There is increasing disquiet about the possibility that PZQ-resistant schistosomes will develop. Consequently, as mathematical modeling predicts, vaccine strategies represent an essential component in the future control of schistosomiasis in China. With the inclusion of focal mollusciciding, improvements in sanitation, and health education into the control scenario, China's target of reducing the level of schistosome infection to less than 1% by 2015 may be achievable.

Inhibition of Schistosoma mansoni thioredoxin-glutathione reductase by auranofin: structural and kinetic aspects

Schistosomiasis is a parasitic disease affecting over 200 million people currently treated with one drug, praziquantel. A possible drug target is the seleno-protein thioredoxin-glutathione reductase (TGR), a key enzyme in the pathway of the parasite for detoxification of reactive oxygen species. The enzyme is a unique fusion of a glutaredoxin domain with a thioredoxin reductase domain, which contains a selenocysteine (Sec) as the penultimate amino acid. Auranofin (AF), a gold-containing compound already in clinical use as an anti-arthritic drug, has been shown to inhibit TGR and to substantially reduce worm burden in mice. Using x-ray crystallography we solved (at 2.5 A resolution) the structure of wild type TGR incubated with AF. The electron density maps show that the actual inhibitor is gold, released from AF. Gold is bound at three different sites not directly involving the C-terminal Sec residue; however, because the C terminus in the electron density maps is disordered, we cannot exclude the possibility that gold may also bind to Sec. To investigate the possible role of Sec in the inactivation kinetics, we tested the effect of AF on a model enzyme of the same superfamily, i.e. the naturally Sec-lacking glutathione reductase, and on truncated TGR. We demonstrate that the role of selenium in the onset of inhibition by AF is catalytic and can be mimicked by an external source of selenium (benzeneselenol). Therefore, we propose that Sec mediates the transfer of gold from its ligands in AF to the redox-active Cys couples of TGR.