Sigma-Class Glutathione Transferases (GSTσ): A New Target with Potential for Helminth Control.

Glutathione transferases (GSTs EC 2.5.1.18) are critical components of phase II metabolism, instrumental in xenobiotics' metabolism. Their primary function involves conjugating glutathione to both endogenous and exogenous toxic compounds, which increases their solubility and enables their ejection from cells. They also play a role in the transport of non-substrate compounds and immunomodulation, aiding in parasite establishment within its host. The cytosolic GST subfamily is the most abundant and diverse in helminths, and sigma-class GST (GSTσ) belongs to it. This review focuses on three key functions of GSTσ: serving as a detoxifying agent that provides drug resistance, functioning as an immune system modulator through its involvement in prostaglandins synthesis, and acting as a vaccine antigen.

Biochemical characterization and gene structure analysis of the 24-kDa glutathione transferase sigma from Taenia solium.

Taenia solium can cause human taeniasis and/or cysticercosis. The latter can in some instances cause human neurocysticercosis which is considered a priority in disease-control strategies and the prevention of mental health problems. Glutathione transferases are crucial for the establishment and long-term survival of T. solium; therefore, we structurally analyzed the 24-kDa glutathione transferase gene (Ts24gst) of T. solium and biochemically characterized its product. The gene promoter showed potential binding sites for transcription factors and xenobiotic regulatory elements. The gene consists of a transcription start site, four exons split by three introns, and a polyadenylation site. The gene architecture is conserved in cestodes. Recombinant Ts24GST (rTs24GST) was active and dimeric. Anti-rTs24GST serum showed slight cross-reactivity with human sigma-class GST. A 3D model of Ts24GST enabled identification of putative residues involved in interactions of the G-site with GSH and of the H-site with CDNB and prostaglandin D2. Furthermore, rTs24GST showed optimal activity at 45 °C and pH 9, as well as high structural stability in a wide range of temperatures and pHs. These results contribute to the better understanding of this parasite and the efforts directed to fight taeniasis/cysticercosis.

Circular dichroism spectroscopic assessment of structural changes upon protein thermal unfolding at contrasting pH: Comparison with molecular dynamics simulations.

In most instances, the usual fastness of protein unfolding events hinders determining changes in secondary structures associated with this process because these determinations rely on the recording of high-resolution circular dichroism (CD) spectra. In this work, far-UV CD spectra, recorded at ten-minute intervals, were used to evaluate the time course followed by four classes of secondary structures in the slow temperature-induced unfolding of yeast triosephosphate isomerase (yTIM) under distinct pH conditions. CONTIN-LL and SELCON3 algorithms were used for the deconvolution of spectra. Both algorithms furnished helix and unordered structure contents that changed according to first-order kinetics, agreeing with the behavior shown by CD data at specific wavelengths. Analyses of unfolded yTIM spectra, using a dataset that includes spectra of unfolded proteins and either one of the two algorithms, clearly showed a more unordered protein structure at high pH; this finding was corroborated with analysis of the difference spectra. Molecular dynamics (MD) simulations performed with AMBER and OPLS force fields resulted in more extensive loss of helices and gain in coils at high pH, in agreement with spectroscopic results. However, structural differences between low- and high-pH unfolded yTIM were relatively small. Comparison of results from CD and MD thus point to the need of fine-tuning of MD procedures.

Non-conventional interactions of N3 inhibitor with the main protease of SARS-CoV and SARS-CoV-2.

The extensive spread of COVID-19 in every continent shows that SARS-CoV-2 virus has a higher transmission rate than SARS-CoV virus which emerged in 2002. This results in a global pandemic that is difficult to control. In this investigation, we analyze the interaction of N3 inhibitor and the main protease of SARS-CoV and SARS-CoV-2 by quantum chemistry calculations. Non-covalent interactions involved in these systems were studied using a model of 469 atoms. Density Functional Theory and Quantum Theory of Atoms in Molecules calculations lead us to the conclusion that non-conventional hydrogen bonds are important to describe attractive interactions in these complexes. The energy of these non-conventional hydrogen bonds represents more than a half of the estimated interaction energy for non-covalent contacts. This means that hydrogen bonds are crucial to correctly describe the bonds between inhibitors and the main proteases. These results could be useful for the design of new drugs, since non-covalent interactions are related to possible mechanisms of action of molecules used against these viruses.

Main interactions of dopamine and risperidone with the dopamine D2 receptor.

Psychosis is one of the psychiatric disorders that is controlled by dopaminergic drugs such as antipsychotics that have affinity for the dopamine D2 receptor (DRD2). In this investigation we perform quantum chemical calculations of two molecules [dopamine and risperidone] within a large cavity of DRD2 that represents the binding site of the receptor. Dopamine is an endogenous neurotransmitter and risperidone is a second-generation antipsychotic. Non-covalent interactions of dopamine and risperidone with DRD2 are analyzed using the Quantum Theory of Atoms in Molecules (QTAIM) and the Non-Covalent Interaction index (NCI). The QTAIM results show that these molecules strongly interact with the receptor. There are 22 non-covalent interactions for dopamine and 54 for risperidone. The electron density evaluated at each critical binding point is small in both systems but it is higher for dopamine than for risperidone, indicating that the interactions of DRD2 with the first are stronger than with the second molecule. However, the binding energy is higher for risperidone (-72.6 kcal mol-1) than for dopamine (-22.8 kcal mol-1). Thus, the strength of the binding energy is due to the number of contacts rather than the strength of the interactions themselves. This could be related to the ability of risperidone to block DRD2 and may explain the efficacy of this drug for controlling the symptoms of schizophrenia, but likewise its secondary effects.

Discovery of a new non-substrate inhibitor of the 26.5 kDa glutathione transferase from Taenia solium by virtual screening.

The inappropriate use of anthelmintics, such as praziquantel and albendazole, has generated resistance and the need to develop new drugs. Glutathione transferases, GSTs, are bisubstrate dimeric enzymes that constitute the main detoxification mechanism against electrophiles, drugs and oxidative damage in Taenia solium. Therefore, GSTs are important targets for the development of new anthelmintics. In this work, we reported a successful virtual screen aimed at the identification of novel inhibitors of a 26.5 kDa GST from T. solium (TsGST26). We found that a compound, i7, able to inhibit selectively TsGST26 concerning human GSTs, showing a non-competitive inhibition mechanism towards substrate glutathione with a Ki (GSH) of 55.7 µM and mixed inhibition towards the electrophilic substrate 1-chloro-2,4-dinitrobenzene with a Ki (CDNB) of 8.64 µM. These results are in agreement with those of docking simulations, which showed i7 binds a site adjacent to the electrophilic site and furthest from the glutathione site.

Correction to: Stabilized Human Cystatin C, Variant L47C/G69C, is a Better Reporter than the Wild-type Inhibitor for Characterizing the Thermodynamics of Binding to Cysteine Proteases.

The original publication of this article contained a number of grammatical errors. Unfortunately, an incorrect version of the file that did not include some final language editing was inadvertently published online. The original article has been corrected.

Stabilized Human Cystatin C Variant L47C/G69C Is a Better Reporter Than the Wild-Type Inhibitor for Characterizing the Thermodynamics of Binding to Cysteine Proteases.

Human cystatin C (HCC) binds and inhibits all types of cysteine proteases from the papain family, including cathepsins (a group of enzymes that participate in a variety of physiological processes), which are some of its natural targets. The affinities of diverse proteases for HCC, expressed as equilibrium binding constants (Kb), range from 106 to 1014 M-1. Isothermal titration calorimetry (ITC) is one of the most useful techniques to characterize the thermodynamics of molecular associations, making it possible to dissect the binding free energy into its enthalpic and entropic components. This information, together with the structural changes that occur during the different associations, could enable better understanding of the molecular basis of affinity. Notwithstanding the high sensitivity of modern calorimeters, ITC requires protein concentrations in at least the 10-100 µM range to obtain reliable data, and it is known that HCC forms oligomers in this concentration range. We present herein a comparative study of the structural, thermal stability, and oligomerization properties of HCC and its stabilized variant (sHCC) L47C/G69C (which possesses an additional disulfide bridge) as well as their binding thermodynamics to the protease chymopapain, analyzed by ITC. The results show that, because sHCC remains monomeric, it is a better reporter than wild-type HCC to characterize the thermodynamics of binding to cysteine proteases.

Coordination numbers in hydrated Cu(II) ions.

The potential energy surface of [Cu(H2O)n]2+ clusters with n = 12, 16, and 18 was explored by using a modified version of the simulated annealing method. Such exploration was carried out by using the PM7 semiempirical method to obtain around 100,000 isomers, which provide candidates to be optimized with PBE0-D3, M06-2X, and BHLYP exchange-correlation functionals coupled with the 6-311++G** basis set. These methods based on the Kohn-Sham approach delivered isomers with coordination numbers of 4, 5, and 6. The analysis used to obtain coordination numbers was based on geometrical parameters and the quantum theory of atoms in molecules (QTAIM) approach. Our methodology found only one isomer with fourfold coordination and its probabilities to appear in these clusters are quite small for high temperatures. The procedure used in this article predicts important populations of fivefold and sixfold coordination clusters, in fact, the fivefold coordination dominates for PBE0-D3 and BHLYP methods, although the sixfold coordination starts to be important when the number of water molecules is increased. The nature of axial and equatorial contacts is discussed in the context of the QTAIM and the noncovalent interaction index (NCI), which gives a clear classification of such orientations. Also, these methods suggest a partial covalent interaction between the Cu2+ and water molecules in both positions; equatorial and axial.

Exploring interfacial water trapping in protein-ligand complexes with multithermal titration calorimetry.

In this work, we examine the hypothesis about how trapped water molecules at the interface between triosephosphate isomerase (TIM) and either of two phosphorylated inhibitors, 2-phosphoglycolate (2PG) or phosphoglycolohydroxamate (PGH), can explain the anomalous highly negative binding heat capacities (ΔCp,b) of both complexes, TIM-2PG and TIM-PGH. We performed fluorimetric titrations of the enzyme with PGH inhibitor under osmotic stress conditions, using various concentrations of either osmolyte: sucrose, ethylene glycol or glycine betaine. We also analyze the binding processes under various stressor concentrations using a novel calorimetric methodology that allows ΔCp,b determinations in single experiments: Multithermal Titration Calorimetry. The binding constant of the TIM-PGH complex decreased gradually with the concentration of all osmolytes, but at diverse extents depending on the osmolyte nature. According to the osmotic stress theory, this decrease indicates that the number of water molecules associated with the enzyme increases with inhibitor binding, i.e. some solvent molecules became trapped. Additionally, the binding heat capacities became less negative at higher osmolyte concentrations, their final values depending on the osmolyte. These effects were also observed in the TIM-2PG complex using sucrose as stressor. Our results strongly suggest that some water molecules became immobilized when the TIM-inhibitor complexes were formed. A computational analysis of the hydration state of the binding site of TIM in both its free state and its complexed form with 2PG or PGH, based on molecular dynamics (MD) simulations in explicit solvent, showed that the binding site effectively immobilized additional water molecules after binding these inhibitors.

Complex kinetics and residual structure in the thermal unfolding of yeast triosephosphate isomerase.

BACKGROUND: Saccharomyces cerevisiae triosephosphate isomerase (yTIM) is a dimeric protein that shows noncoincident unfolding and refolding transitions (hysteresis) in temperature scans, a phenomenon indicative of the slow forward and backward reactions of the native-unfolded process. Thermal unfolding scans suggest that no stable intermediates appear in the unfolding of yTIM. However, reported evidence points to the presence of residual structure in the denatured monomer at high temperature. RESULTS: Thermally denatured yTIM showed a clear trend towards the formation of aggregation-prone, ß-strand-like residual structure when pH decreased from 8.0 to 6.0, even though thermal unfolding profiles retained a simple monophasic appearance regardless of pH. However, kinetic studies performed over a relatively wide temperature range revealed a complex unfolding mechanism comprising up to three observable phases, with largely different time constants, each accompanied by changes in secondary structure. Besides, a simple sequential mechanism is unlikely to explain the observed variation of amplitudes and rate constants with temperature. This kinetic complexity is, however, not linked to the appearance of residual structure. Furthermore, the rate constant for the main unfolding phase shows small, rather unvarying values in the pH region where denatured yTIM gradually acquires a ß-strand-like conformation. It appears, therefore, that the residual structure has no influence on the kinetic stability of the native protein. However, the presence of residual structure is clearly associated with increased irreversibility. CONCLUSIONS: The slow temperature-induced unfolding of yeast TIM shows three kinetic phases. Rather than a simple sequential pathway, a complex mechanism involving off-pathway intermediates or even parallel pathways may be operating. ß-strand-type residual structure, which appears below pH 8.0, is likely to be associated with increased irreversible aggregation of the unfolded protein. However, this denatured form apparently accelerates the refolding process.

Effects of pH on the association between the inhibitor cystatin and the proteinase chymopapain.

Cysteine proteinases are involved in many aspects of physiological regulation. In humans, some cathepsins have shown another function in addition to their role as lysosomal proteases in intracellular protein degradation; they have been implicated in the pathogenesis of several heart and blood vessel diseases and in cancer development. In this work, we present a fluorometric and computational study of the binding of one representative plant cysteine proteinase, chymopapain, to one of the most studied inhibitors of these proteinases: chicken cystatin. The binding equilibrium constant, Kb, was determined in the pH range between 3.5 and 10.0, revealing a maximum in the affinity at pH 9.0. We constructed an atomic model for the chymopapain-cystatin dimer by docking the individual 3D protein structures; subsequently, the model was refined using a 100 ns NPT molecular dynamics simulation in explicit water. Upon scrutiny of this model, we identified 14 ionizing residues at the interface of the complex using a cutoff distance of 5.0 Å. Using the pKa values predicted with PROPKA and a modified proton-linkage model, we performed a regression analysis on our data to obtain the composite pKavalues for three isoacidic residues. We also calculated the electrostatic component of the binding energy (ΔGb,elec) at different pH values using an implicit solvent model and APBS software. The pH profile of this calculated energy compares well with the experimentally obtained binding energy, ΔGb. We propose that the residues that form an interchain ionic pair, Lys139A from chymopapain and Glu19B from cystatin, as well as Tyr61A and Tyr67A from chymopapain are the main residues responsible for the observed pH dependence in the chymopapain- cystatin affinity.

Surface Tension of Organic Liquids Using the OPLS/AA Force Field.

Molecular dynamics simulations are performed to obtain the surface tension of 61 organic liquids using the OPLS/AA (all-atom optimized potential for liquid simulations). The force field parameters are the same as those recently used (Caleman et al. J. Chem. Theory Comput.2012, 8, 61) to determine several thermodynamic properties of 146 organic liquids. The correct evaluation of surface tension using slab simulations of liquids requires one to properly take into account the long-range interactions (Trukhymchuk and Alejandre J. Chem. Phys.1999, 111, 8510). In addition, the liquid density from slab simulations has to be the same as that obtained in liquid simulations at constant temperature and pressure. The new results of surface tensions from this work improve those reported by Caleman et al. The OPLS/AA force field gives good surface tensions compared with experimental data for most of the systems studied in this work, although it was developed to simulate liquids.

Additive effect of single amino acid replacements on the kinetic stability of ß-glucosidase B.

Previously, we applied in vitro evolution to generate the thermoresistant triple mutant H62R/N223Y/M319I of ß-glucosidase B (BglB) from Paenibacillus polymyxa. In order to dissect the energetic contributions to protein stabilization achieved by these mutations, we measured the kinetic constants of the heat denaturation of wild type BglB, the triple mutant and the three single mutants (H62R, N223Y, M319I) by circular dichroism at various temperatures. Our results show that all four mutants delayed the denaturation process. Based on the Transition State theory, the increase of the activation barrier for the thermal denaturation of the triple mutant (ΔΔG ( N→TS )) is equivalent to that produced by the sum of the contributions from the three single mutants, whose C ( ß ) s are located at least 18 Å apart. This analysis provides a formal demonstration of the generally accepted idea that protein thermal stability can be increased through sequential addition of individual mutations. Each of the mutations described here contribute in part to the overall effect, which in this case affects the unfolding barrier.

Thermal denaturation of ß-glucosidase B from Paenibacillus polymyxa proceeds through a Lumry-Eyring mechanism.

ß-glucosidase B (BglB), 1,4-ß-D: -glucanohydrolase, is an enzyme with various technological applications for which some thermostable mutants have been obtained. Because BglB denatures irreversibly with heating, the stabilities of these mutants are assessed kinetically. It, therefore, becomes relevant to determine whether the measured rate constants reflect one or several elementary kinetic steps. We have analyzed the kinetics of heat denaturation of BglB from Paenibacillus polymyxa under various conditions by following the loss of secondary structure and enzymatic activity. The denaturation is accompanied by aggregation and an initial reversible step at low temperatures. At T ≥ T ( m ), the process follows a two-state irreversible mechanism for which the kinetics does not depend on the enzyme concentration. This behavior can be explained by a Lumry-Eyring model in which the difference between the rates of the irreversible and the renaturation steps increases with temperature. Accordingly, at high scan rates (≥1 °C min(-1)) or temperatures (T ≥ T ( m )), the measurable activation energy involves only the elementary step of denaturation.

Binding thermodynamics of phosphorylated inhibitors to triosephosphate isomerase and the contribution of electrostatic interactions.

Electrostatic interactions have a central role in some biological processes, such as recognition of charged ligands by proteins. We characterized the binding energetics of yeast triosephosphate isomerase (TIM) with phosphorylated inhibitors 2-phosphoglycollate (2PG) and phosphoglycolohydroxamate (PGH). We determined the thermodynamic parameters of the binding process (K(b), ΔG(b), ΔH(b), ΔS(b) and ΔC(p)) with different concentrations of NaCl, using fluorimetric and calorimetric titrations in the conventional mode of ITC and a novel method, multithermal titration calorimetry (MTC), which enabled us to measure ΔC(p) in a single experiment. We ruled out specific interactions of Na(+) and Cl(-) with the native enzyme and did not detect significant linked protonation effects upon the binding of inhibitors. Increasing ionic strength (I) caused K(b), ΔG(b) and ΔH(b) to become less favorable, while ΔS(b) became less unfavorable. From the variation of K(b) with I, we determined the electrostatic contribution of TIM-2PG and TIM-PGH to ΔG(b) at I=0.06 M and 25 °C to be 36% and 26%, respectively. The greater affinity of PGH for TIM is due to a more favorable ΔH(b) compared to 2PG (by 19-24 kJ mol(-1) at 25 °C). This difference is compatible with PGH establishing up to five more hydrogen bonds with TIM. Both binding ΔC(p)s were negative, and less negative with increasing ionic strength. ΔC(p)s at I=0.06 M were much more negative than predicted by surface area models. Water molecules trapped in the interface when ligands bind to protein could explain the highly negative ΔCps. Thermodynamic binding functions for TIM-2PG changed more with ionic strength than those for TIM-PGH. This greater dependence is consistent with linked, but compensated, protonation equilibriums yielding the dianionic species of 2PG that binds to TIM, process that is not required for PGH.

Effect of a specific inhibitor on the unfolding and refolding kinetics of dimeric triosephosphate isomerase: establishing the dimeric and similarly structured nature of the main transition states on the forward and backward reactions.

2-Phosphoglycolate (PGA), a strong competitive inhibitor of the dimeric enzyme triosephosphate isomerase (TIM), brings about a large decrease in the unfolding rate constant of the protein. The data set of rate constants versus ligand concentration may be equally well explained by regarding either a monomeric or a dimeric transition state (TS). However, if the thermodynamics for binding of PGA to native TIM is taken into account, it becomes clear that a dimeric TS is the right assumption. Furthermore, by studying the effect of the ligand on the second-order refolding reaction, we found results indicating similar PGA-binding affinities to be present in the transition states for the rate-limiting steps of the forward and backward reactions. Most likely, therefore, both TS resemble each other in respect to the active site architecture. It should be mentioned, however, that our data do not rule out the possible occurrence of an unstable, (partially) folded monomeric intermediate, which would rapidly interconvert with the unfolded monomer.

Effect of a new ionic pair on the unfolding activation barrier of beta-glucosidase B.

Thermal unfolding kinetics of beta-glucosidase B from Paenibacillus polymyxa and its thermoresistant mutant H62R were determined from far-UV circular dichroism (CD) measurements at different temperatures. The unfolding of both enzymes followed simple two-state kinetics. The new ionic pair formed between Arg62 and Glu429 in the H62R variant did not change substantially the enzyme structure as judged by far-UV CD and fluorescence spectra, but produced an increase in the unfolding activation barrier of 0.95 +/- 0.10 kcal mol(-1), in good agreement with the energetic contribution reported for surface salt bridges in proteins. Eyring's analysis of the unfolding kinetic constants showed that the activation enthalpies for thermal denaturation of both enzymes were essentially the same. Thus, the greater kinetic stability rendered by the salt bridge seems to be due to a reduction in the activation entropy.

Multithermal titration calorimetry: a rapid method to determine binding heat capacities.

Herein a new method that allows binding DeltaCp to be determined with a single experiment is presented. Multithermal titration calorimetry (MTC) is a simple extension of isothermal titration calorimetry (ITC) that explicitly takes into account the thermal dependences of DeltaH and the binding constant. Experimentally, this is accomplished by performing a single stepwise titration with ITC equipment, allowing temperature re-adjustments of the system at intermediate states of the titration process. Thus, from the resulting multitherm, DeltaCp can also be determined. The experimental feasibility of MTC was tested by using the well-characterized lysozyme-chitotriose complex as a model system.

Conserved cysteine 126 in triosephosphate isomerase is required not for enzymatic activity but for proper folding and stability.

In triosephosphate isomerase, Cys126 is a conserved residue located close to the catalytic glutamate, Glu165. Although it has been mentioned that Cys126 and other nearby residues are required to maintain the active site geometry optimal for catalysis, no evidence supporting this idea has been reported to date. In this work, we studied the catalytic and stability properties of mutants C126A and C126S of Saccharomyces cerevisiae TIM (wtTIM). None of these amino acid replacements induced significant changes in the folding of wtTIM, as indicated by spectroscopic studies. C126S and C126A have K(M) and k(cat) values that are concomitantly reduced by only 4-fold and 1.5-fold, respectively, compared to those of wtTIM; in either case, however, the catalytic efficiency (k(cat)/K(M)) of the enzyme is barely affected. The affinity of mutated TIMs for the competitive inhibitor 2-phosphoglycolate augmented also slightly. In contrast, greater susceptibility to thermal denaturation resulted from mutation of Cys126, especially when it was changed to Ser. By using values of the rate constants for unfolding and refolding, we estimated that, at 25 degrees C, C126A and C126S are less stable than wtTIM by about 5.0 and 9.0 kcal mol(-)(1), respectively. Moreover, either of these mutations slows down the folding rate by a factor of 10 and decreases the recovery of the active enzyme after thermal unfolding. Thus, Cys126 is required for proper stability and efficient folding of TIM rather than for enzymatic catalysis.