Conformational changes and reaction of clostridial glycosylating toxins

The crystal structures of the catalytic fragments of ‘lethal toxin’ from Clostridium sordellii and of ‘α-toxin’ from Clostridium novyi have been established. Almost half of the residues follow the chain fold of the glycosyl-transferase type A family of enzymes; the other half forms large α-helical protrusions that are likely to confer specificity for the respective targeted subgroup of Rho proteins in the cell. In the crystal, the active center of α-toxin contained no substrates and was disassembled, whereas that of lethal toxin, which was ligated with the donor substrate UDP-glucose and cofactor Mn^2 +, was catalytically competent. Surprisingly, the structure of lethal toxin with Ca^2 + (instead of Mn^2 +) at the cofactor position showed a bound donor substrate with a disassembled active center, indicating that the strictly octahedral coordination sphere of Mn^2 + is indispensable to the integrity of the enzyme. The homologous structures of α-toxin without substrate, distorted lethal toxin with Ca^2 + plus donor, active lethal toxin with Mn^2 + plus donor and the homologous Clostridium difficile toxin B with a hydrolyzed donor have been lined up to show the geometry of several reaction steps. Interestingly, the structural refinement of one of the three crystallographically independent molecules of Ca^2 +-ligated lethal toxin resulted in the glucosyl half-chair conformation expected for glycosyl-transferases that retain the anomeric configuration at the C1″ atom. A superposition of six acceptor substrates bound to homologous enzymes yielded the position of the nucleophilic acceptor atom with a deviation of < 1 Å. The resulting donor-acceptor geometry suggests that the reaction runs as a circular electron transfer in a six-membered ring, which involves the deprotonation of the nucleophile by the β-phosphoryl group of the donor substrate UDP-glucose.     

A maurotoxin with constrained standard disulfide bridging: innovative strategy of chemical synthesis,pharmacology, and docking on K+ channels

Maurotoxin (MTX) is a 34-residue toxin that has been isolated initially from the venom of the scorpion Scorpio maurus palmatus. It presents a large number of pharmacological targets, including small conductance Ca2+-activated and voltage-gated K+ channels. Contrary to other toxins of the alpha-KTx6 family (Pi1, Pi4, Pi7, and HsTx1), MTX exhibits a unique disulfide bridge organization of the type C1-C5, C2-C6, C3-C4, and C7-C8 (instead of the conventional C1-C5, C2-C6, C3-C7, and C4-C8, herein referred to as Pi1-like) that does not prevent its folding along the classic alpha/beta scaffold of scorpion toxins. Here, we developed an innovative strategy of chemical peptide synthesis to produce an MTX variant (MTXPi1) with a conventional pattern of disulfide bridging without any alteration of the toxin chemical structure. This strategy was used solely to address the impact of half-cystine pairings on MTX structural properties and pharmacology. The data indicate that MTXPi1 displays some marked changes in affinities toward the target K+ channels. Computed docking analyses using molecular models of both MTXPi1 and the various voltage-gated K+ channel subtypes (Shaker B, Kv1.2, and Kv1.3) were found to correlate with MTXPi1 pharmacology. A functional map detailing the interaction between MTXPi1 and Shaker B channel was generated in line with docking experiments.     

Strong effects of an individual water molecule on the rate of light-driven charge separation in the Rhodobacter sphaeroides reaction center

The role of a water molecule (water A) located between the primary electron donor (P) and first electron acceptor bacteriochlorophyll (B(A)) in the purple bacterial reaction center was investigated by mutation of glycine M203 to leucine (GM203L). The x-ray crystal structure of the GM203L reaction center shows that the new leucine residue packs in such a way that water A is sterically excluded from the complex, but the structure of the protein-cofactor system around the mutation site is largely undisturbed. The results of absorbance and resonance Raman spectroscopy were consistent with either the removal of a hydrogen bond interaction between water A and the keto carbonyl group of B(A) or a change in the local electrostatic environment of this carbonyl group. Similarities in the spectroscopic properties and x-ray crystal structures of reaction centers with leucine and aspartic acid mutations at the M203 position suggested that the effects of a glycine to aspartic acid substitution at the M203 position can also be explained by steric exclusion of water A. In the GM203L mutant, loss of water A was accompanied by an approximately 8-fold slowing of the rate of decay of the primary donor excited state, indicating that the presence of water A is important for optimization of the rate of primary electron transfer. Possible functions of this water molecule are discussed, including a switching role in which the redox potential of the B(A) acceptor is rapidly modulated in response to oxidation of the primary electron donor.     

Polypeptide Hormones and Chromatin-Associated Proteins Act as Acceptors for Cholera Toxin-Catalyzed ADP-Ribosylation

Abstract: Cholera toxin catalyzed the ADP-ribosylation of the pituitary protein hormones thyrotropin (TSH), lutropin (LH), follitropin (FSH), human chorionic gonadotropin (hCG). and corticotropin (ACTH)_1–24, and ADP-ribosylation of the basic proteins histone subfraction H₁ and protamine. Casein and phosvitin, acidic nuclear proteins, did not act as acceptors for toxin-catalyzed ADP-ribosylation. The isolated TSH A and B subunits were tested for their ADP-ribose acceptor activity. The TSH A subunit showed fourfold greater ADP-ribose acceptor activity than the TSH B subunit. The ADP-ribose acceptor protein protamine was analyzed by sodium dodecyl sulfate polyacrylamide gel electrophoresis following incubation with cholera toxin under ADP-ribosylating conditions. [³H]ADP-ribose incorporated into protein from [³H]NAD migrated with the acceptor protein protamine. In the absence of added acceptor protein, the [³H]ADP-ribose incorporated into protein migrated with the A₁ fragment of cholera toxin. Cholera toxin A and B subunits were isolated and tested for their ability to catalyze the transfer of ADP-ribose to protamine. The cholera toxin A subunit showed 50-fold greater ADP-ribosyltransferase activity than the B subunit. Our data indicate that a variety of adenohypophyseal hormones and regulatory proteins act as acceptors for toxin-catalyzed ADP-ribosylation. These studies may help in understanding the role of endogenous ADP-ribosyltransferases and the physiological effects of this modification of protein.     

Role of glutamine 148 of human 15-hydroxyprostaglandin dehydrogenase in catalytic oxidation of prostaglandin E2

NAD+-dependent 15-hydroxyprostaglandin dehydrogenase (15-PGDH), a member of the short-chain dehydrogenase/reductase (SDR) family, catalyzes the first step in the catabolic pathways of prostaglandins and lipoxins. This enzyme oxidizes the C-15 hydroxyl group of prostaglandins and lipoxins to produce 15-keto metabolites which exhibit greatly reduced biological activities. A three-dimensional (3D) structure of 15-PGDH based on the crystal structures of the levodione reductase and tropinone reductase-II was generated and used for docking study with NAD+ coenzyme and PGE2 substrate. Three well-conserved residues among SDR family which correspond to Ser-138, Tyr-151, and Lys-155 of 15-PGDH have been shown to participate in the catalytic reaction. Based on the molecular interactions observed from 3D structure of 15-PGDH, we further propose that Gln-148 in 15-PGDH is important in properly positioning the 15-hydroxyl group of PGE2 by hydrogen bonding with the side-chain oxygen atom of Gln-148. This residue is found to be less conserved and replaceable by glutamyl, histidinyl, and asparaginyl residues in SDR family. Accordingly, site-directed mutagenesis of Gln-148 of 15-PGDH to alanine, glutamic acid, histidine, and asparagine (Q148A, Q148E, Q148H, and Q148N) was carried out. The activity of mutant Q148A was not detectable, whereas those of mutants Q148E, Q148H, and Q148N were comparable to or higher than the wild type. This indicates that the side-chain oxygen or nitrogen atom at position 148 of 15-PGDH plays an important role in anchoring C-15 hydroxyl group of PGE2 through hydrogen bonding for catalytic reaction.     

Modeling the three-dimensional structures of an unbound single-chain variable fragment (scFv) and its hypothetical complex with a Corynespora cassiicola toxin, cassiicolin

Rubber trees infected with a host-specific cassiicolin toxin often experience considerable leaf fall, which in turn results in loss of crop productivity. It was recently revealed that cassiicolin-specific single-chain variable fragments (scFv) can successfully reduce the toxic effects of cassiicolin. However, the detailed mechanism of antibody action remains poorly understood. The primary sequence of the newly sequenced cassiicolin-specific scFv was highly homologous to several members of single-chain antibodies in the 14B7 family. In this study, with the aid of homology modeling, the three-dimensional structure of cassiicolin-specific scFv was elucidated, and was found to exhibit a characteristic immunoglobulin fold that mainly consists of β sheets. Additionally, molecular docking between the modeled scFv antibody and the available three-dimensional crystal structure of cassiicolin toxin was also performed. The predicted structural complex and the change in accessible surface area between the toxin and the scFv antibody upon complexation reveal the potential role of certain complementarity determining region (CDR) amino acid residues in the formation of the complex. These computational results suggest that mutagenesis experiments that are aimed at validating the model and improving the binding affinity of cassiicolin-specific scFv antibodies for the toxin should be performed.     

Crystal structure of Paracoccus denitrificans electron transfer flavoprotein: structural and electrostatic analysis of a conserved flavin binding domain

The crystal structure of electron transfer flavoprotein (ETF) from Paracoccus denitrificans was determined and refined to an R-factor of 19.3% at 2.6 A resolution. The overall fold is identical to that of the human enzyme, with the exception of a single loop region. Like the human structure, the structure of the P. denitrificans ETF is comprised of three distinct domains, two contributed by the alpha-subunit and the third from the beta-subunit. Close analysis of the structure reveals that the loop containing betaI63 is in part responsible for conferring the high specificity of AMP binding by the ETF protein. Using the sequence and structures of the human and P. denitrificans enzymes as models, a detailed sequence alignment has been constructed for several members of the ETF family, including sequences derived for the putative FixA and FixB proteins. From this alignment, it is evident that in all members of the ETF family the residues located in the immediate vicinity of the FAD cofactor are identical, with the exception of the substitution of serine and leucine residues in the W3A1 ETF protein for the human residues alphaT266 and betaY16, respectively. Mapping of ionic differences between the human and P. denitrificans ETF onto the structure identifies a surface that is electrostatically very similar between the two proteins, thus supporting a previous docking model between human ETF and pig medium-chain acyl-CoA dehydrogenase (MCAD). Analysis of the ionic strength dependence of the electron transfer reaction between either human or P. denitrificans ETF and MCAD demonstrates that the human ETF functions optimally at low ( approximately 10 mequiv) ionic strength, while P. denitrificans ETF is a better electron acceptor at higher (>75 mequiv) ionic strength. This suggests that the electrostatic surface potential of the two proteins is very different and is consistent with the difference in isoelectric points between the proteins. Analysis of the electrostatic potentials of the human and P. denitrificans ETFs reveals that the P. denitrificans ETF is more negatively charged. This excess negative charge may contribute to the difference in redox potentials between the two ETF flavoproteins and suggests an explanation for the opposing ionic strength dependencies for the reaction of MCAD with the two ETFs. Furthermore, by analysis of a model of the previously described human-P. denitrificans chimeric ETF protein, it is possible to identify one region of ETF that participates in docking with ETF-ubiquinone oxidoreductase, the physiological electron acceptor for ETF.     

Bacillus thuringiensis insecticidal Cry1Aa toxin binds to a highly conserved region of aminopeptidase N in the host insect leading to its evolutionary success.

Bacillus thuringiensis insecticidal protein, Cry1Aa toxin, binds to a specific receptor in insect midguts and has insecticidal activity. Therefore, the structure of the receptor molecule is probably a key factor in determining the binding affinity of the toxin and insect susceptibility. The cDNA fragment (PX frg1) encoding the Cry1Aa toxin-binding region of an aminopeptidase N (APN) or an APN family protein from diamondback moth, Plutella xylostella midgut was cloned and sequenced. A comparison between the deduced amino acid sequence of PX frg1 and other insect APN sequences shows that Cry1Aa toxin binds to a highly conserved region of APN family protein. In this paper, we propose a model to explain the mechanism that causes B. thuringiensis evolutionary success and differing insect susceptibility to Cry1Aa toxin.     

In silico modelling and molecular dynamics simulation studies of thiazolidine based PTP1B inhibitors

Protein tyrosine phosphatase 1B (PTP1B) has been identified as a negative regulator of insulin and leptin signalling pathway; hence, it can be considered as a new therapeutic target of intervention for the treatment of type 2 diabetes. Inhibition of this molecular target takes care of both diabetes and obesity, i.e. diabestiy. In order to get more information on identification and optimization of lead, pharmacophore modelling, atom-based 3D QSAR, docking and molecular dynamics studies were carried out on a set of ligands containing thiazolidine scaffold. A six-point pharmacophore model consisting of three hydrogen bond acceptor (A), one negative ionic (N) and two aromatic rings (R) with discrete geometries as pharmacophoric features were developed for a predictive 3D QSAR model. The probable binding conformation of the ligands within the active site was studied through molecular docking. The molecular interactions and the structural features responsible for PTP1B inhibition and selectivity were further supplemented by molecular dynamics simulation study for a time scale of 30 ns. The present investigation has identified some of the indispensible structural features of thiazolidine analogues which can further be explored to optimize PTP1B inhibitors.     

Solution structure of BmKK2, a new potassium channel blocker from the venom of chinese scorpion Buthus martensi Karsch

A natural K+ channel blocker, BmKK2 (a member of scorpion toxin subfamily alpha-KTx 14), which is composed of 31 amino acid residues and purified from the venom of the Chinese scorpion Buthus martensi Karsch, was characterized using whole-cell patch-clamp recording in rat hippocampal neurons. The three dimensional structure of BmKK2 was determined with two-dimensional NMR spectroscopy and molecular modelling techniques. In solution this toxin adopted a common alpha/beta-motif, but showed distinct local conformation in the loop between alpha-helix and beta-sheet in comparison with typical short-chain scorpion toxins (e.g., CTX and NTX). Also, the alpha helix is shorter and the beta-sheet element is smaller (each strand consisted only two residues). The unusual structural feature of BmKK2 was attributed to the shorter loop between the alpha-helix and beta-sheet and the presence of two consecutive Pro residues at position 21 and 22 in the loop. Moreover, two models of BmKK2/hKv1.3 channel and BmKK2/rSK2 channel complexes were simulated with docking calculations. The results demonstrated the existence of a alpha-mode binding between the toxin and the channels. The model of BmKK2/rSK2 channel complex exhibited favorable contacts both in electrostatic and hydrophobic, including a network of five hydrogen bonds and bigger interface containing seven pairs of inter-residue interactions. In contrast, the model of BmKK2/hKv1.3 channel complex, containing only three pairs of inter-residue interactions, exhibited poor contacts and smaller interface. The results well explained its lower activity towards Kv channel, and predicted that it may prefer a type of SK channel with a narrower entryway as its specific receptor.     

Ligand-induced conformational changes: improved predictions of ligand binding conformations and affinities

A computational docking strategy using multiple conformations of the target protein is discussed and evaluated. A series of low molecular weight, competitive, nonpeptide protein tyrosine phosphatase inhibitors are considered for which the x-ray crystallographic structures in complex with protein tyrosine phosphatase 1B (PTP1B) are known. To obtain a quantitative measure of the impact of conformational changes induced by the inhibitors, these were docked to the active site region of various structures of PTP1B using the docking program FlexX. Firstly, the inhibitors were docked to a PTP1B crystal structure cocrystallized with a hexapeptide. The estimated binding energies for various docking modes as well as the RMS differences between the docked compounds and the crystallographic structure were calculated. In this scenario the estimated binding energies were not predictive inasmuch as docking modes with low estimated binding energies corresponded to relatively large RMS differences when aligned with the corresponding crystal structure. Secondly, the inhibitors were docked to their parent protein structures in which they were cocrystallized. In this case, there was a good correlation between low predicted binding energy and a correct docking mode. Thirdly, to improve the predictability of the docking procedure in the general case, where only a single target protein structure is known, we evaluate an approach which takes possible protein side-chain conformational changes into account. Here, side chains exposed to the active site were considered in their allowed rotamer conformations and protein models containing all possible combinations of side-chain rotamers were generated. To evaluate which of these modeled active sites is the most likely binding site conformation for a certain inhibitor, the inhibitors were docked against all active site models. The receptor rotamer model corresponding to the lowest estimated binding energy is taken as the top candidate. Using this protocol, correct inhibitor binding modes could successfully be discriminated from proposed incorrect binding modes. Moreover, the ranking of the estimated ligand binding energies was in good agreement with experimentally observed binding affinities.     

X-ray studies reveal lanthanide binding sites at the A/B5 interface of E. coli heat labile enterotoxin.

The crystal structure determination of heat labile enterotoxin (LT) bound to two different lanthanide ions, erbium and samarium, revealed two distinct ion binding sites in the interface of the A subunit and the B pentamer of the toxin. One of the interface sites is conserved in the very similar cholera toxin sequence. These sites may be potential calcium binding sites. Erbium and samarium binding causes a change in the structure of LT: a rotation of the A1 subunit of up to two degrees relative to the B pentamer.     

Modeling the structure of mAb 14B7 bound to the anthrax protective antigen

The anthrax protective antigen (PA) is a key component of the tripartite anthrax toxin. Monoclonal antibody (mAb) 14B7 and its engineered, affinity-matured variants have been shown to be effective in blocking PA binding to cellular receptors and mitigating anthrax toxicity. Here, we perform computational structural modeling of the mAb 14B7-PA interaction. Our objectives are to determine the structure of the 14B7-PA complex, to deduce a structural explanation for the affinity maturation from the docking models, and to study the effect of inaccuracies in the antibody homology model on docking. We used the RosettaDock program to dock PA with the mAb 14B7 crystal structure or homology model. Our simulations generate two distinct binding orientations consistent with experimental residue mutations that diminish 14B7-PA binding. Furthermore, the models suggest new site-directed mutations to positively identify one of these two solutions as the correct 14B7-PA docking orientation. The models indicate that PA regions 648-660 and 712-720 may be important for 14B7 binding in addition to the known PA epitope, and the binding interfaces are similar to that seen in the PA complex with cellular receptor CMG2. Antibody residues involved in affinity maturation do not contact the antigen in the docking models, suggesting that affinity maturation in the 14B7 family does not result from direct enhancements of antibody-antigen contacts. Docking the homology model produces low-resolution representations of the crystal structure docking orientations, but homology model docking is frustrated by antibody H3 loop conformation errors. This work demonstrates the usefulness and limitations of computational structure prediction for the development of antibody therapeutics, and reemphasizes the need for flexible backbone docking algorithms to achieve high-resolution docking using homology models.     

Rho-glucosylating Clostridium difficile toxins A and B: new insights into structure and function

Clostridium difficile causes pseudomembranous colitis and is responsible for many cases of nosocomial antibiotic-associated diarrhea. Major virulence factors of C. difficile are the glucosylating exotoxins A and B. Both toxins enter target cells in a pH- dependent manner from endosomes by forming pores. They translocate the N-terminal catalytic domains into the cytosol of host cells and inactivate Rho guanosine triphosphatases by glucosylation. The crystal structure of the catalytic domain of toxin B was solved in a complex with uridine diphosphate, glucose, and manganese ion, exhibiting a folding of type A family glycosyltransferases. Crystallization of fragments of the C-terminus of toxin A, which is characterized by polypeptide repeats, revealed a solenoid-like structure often found in bacterial cell surface proteins. These studies, which provide new insights into structure, uptake, and function of the family of clostridial glucosylating toxins, are reviewed.     

The role of the diphthamide-containing loop within eukaryotic elongation factor 2 in ADP-ribosylation by Pseudomonas aeruginosa exotoxin A

eEF2 (eukaryotic elongation factor 2) contains a post-translationally modified histidine residue, known as diphthamide, which is the specific ADP-ribosylation target of diphtheria toxin, cholix toxin and Pseudomonas aeruginosa exotoxin A. Site-directed mutagenesis was conducted on residues within the diphthamide-containing loop (Leu693-Gly703) of eEF2 by replacement with alanine. The purified yeast eEF2 mutant proteins were then investigated to determine the role of this loop region in ADP-ribose acceptor activity of elongation factor 2 as catalysed by exotoxin A. A number of single alanine substitutions in the diphthamide-containing loop caused a significant reduction in the eEF2 ADP-ribose acceptor activities, including two strictly conserved residues, His694 and Asp696. Analysis by MS revealed that all of these mutant proteins lacked the 2'-modification on the His699 residue and that eEF2 is acetylated at Lys509. Furthermore, it was revealed that the imidazole ring of Diph699 (diphthamide at position 699) still functions as an ADP-ribose acceptor (albeit poorly), even without the diphthamide modification on the His699. Therefore, this diphthamide-containing loop plays an important role in the ADP-ribosylation of eEF2 catalysed by toxin and also for modification of His699 by the endogenous diphthamide modification machinery.     

Analyzing of expression of novel polypeptide complexes consisting of Shiga toxin B subunit and Adherence Fimbriae of Escherichia coli based on in silico modeling

Enterohemorrhagic (EHEC) and enteroaggregative (EAEC) are two pathotypes of diarrheagenic Escherichia coli. EAEC strains express adhesins called aggregate adherence fimbriae (AAFs) which the bacteria use to adhere to intestinal mucosa. EHEC virulence factor is Shiga toxin which belongs to the AB5 toxin family. B subunit, the nontoxic part of Shiga toxin (StxB), forms a homo pentamer and is responsible for binding to target cells. StxB has recently been proven to have adjuvant activity. In the current study we fused StxB encoding gene to 3' end of genes encoding two variants of AAFs, i.e., AAF/I and AAF/II. The in silico studies on tertiary structure and biochemical characteristics of Shiga toxin A subunit (StxA) revealed more resemblance to AAF/II than AAF/I. The constructs were prepared in a way that StxB could imitate its natural structure (pentamer formation) and its position (C-terminus) in the native toxin complex. The expression of these constructs showed the formation of AAF/II-B as a protein complex but with lower molecular mass than its expected size. In contrast, the AAF/I-B complex was not formed. Overall, the results of in silico studies and expression experiments together revealed that despite AAF/II-B expression, StxB failed to form pentamer. Therefore the observed protein complex has lower molecular mass. Since StxB is bound to AAF/II through disulfide bond, this bond prevents pentamer formation of StxB. However, due to the lack of disulfide bond between AAF/I and StxB, no protein complex is formed, thus StxB maintains its pentamer structure.     

A holistic molecular docking approach for predicting protein-protein complex structure

A holistic protein-protein molecular docking approach, HoDock, was established, composed of such steps as binding site prediction, initial complex structure sampling, refined complex structure sampling, structure clustering, scoring and final structure selection. This article explains the detailed steps and applications for CAPRI Target 39. The CAPRI result showed that three predicted binding site residues, A191HIS, B512ARG and B531ARG, were correct, and there were five submitted structures with a high fraction of correct receptor-ligand interface residues, indicating that this docking approach may improve prediction accuracy for protein-protein complex structures.     

How to choose templates for modeling of protein complexes: Insights from benchmarking template-based docking

Comparative docking is based on experimentally determined structures of protein-protein complexes (templates), following the paradigm that proteins with similar sequences and/or structures form similar complexes. Modeling utilizing structure similarity of target monomers to template complexes significantly expands structural coverage of the interactome. Template-based docking by structure alignment can be performed for the entire structures or by aligning targets to the bound interfaces of the experimentally determined complexes. Systematic benchmarking of docking protocols based on full and interface structure alignment showed that both protocols perform similarly, with top 1 docking success rate 26%. However, in terms of the models' quality, the interface-based docking performed marginally better. The interface-based docking is preferable when one would suspect a significant conformational change in the full protein structure upon binding, for example, a rearrangement of the domains in multidomain proteins. Importantly, if the same structure is selected as the top template by both full and interface alignment, the docking success rate increases 2-fold for both top 1 and top 10 predictions. Matching structural annotations of the target and template proteins for template detection, as a computationally less expensive alternative to structural alignment, did not improve the docking performance. Sophisticated remote sequence homology detection added templates to the pool of those identified by structure-based alignment, suggesting that for practical docking, the combination of the structure alignment protocols and the remote sequence homology detection may be useful in order to avoid potential flaws in generation of the structural templates library.     

Highlighting a π–π interaction: a protein modeling and molecular dynamics simulation study on Anopheles gambiae glutathione S-transferase 1-2

Cytosolic insect theta class glutathione S-transferases (GSTs) have not been studied completely and their physiological roles are unknown. A detailed understanding of Anopheles gambiae GST (Aggst1-2) requires an accurate structure, which has not yet been determined. A high quality model structure of Aggst1-2 was constructed using homology modeling and the ligand–protein complex was obtained by the docking method. Molecular dynamics (MD) simulations were carried out to study conformational changes and to calculate binding free energy. The results of MD simulation indicate that Aggst1-2 undergoes small conformational changes after ligands dock to the protein, which facilitate the catalytic reaction. An essential hydrogen bond was found between the sulfur atom of glutathione (GSH) and the hydrogen atom of hydroxyl group in Ser9, which was in good agreement with experimental data. A π–π interaction between Phe204 and CDNB ligand was also found. This interaction seems to be important in stabilization of the ligand. Further study of binding free energy decomposition revealed a van der Waals interaction between two ligands that may play a key role in nucleophilic addition reaction. This work will be a good starting point for further determination of the biological role of cytosolic insect theta class GSTs and will aid the design of structure-based inhibitors.     

Structure-activity relationship study on the 6-membered heteroaromatic ring system of diphenylpyrazine-type prostacyclin receptor agonists

A series of prostacyclin receptor agonists was prepared by modifying the central heteroaromatic ring of lead compound 2, and a docking study was performed to investigate their structure-activity relationships by using a homology-modeled structure of the prostacyclin receptor. Compound 2 and its derivatives could be docked to the prostacyclin receptor in two ways depending on the position of the nitrogen atom within the heteroaromatic ring. Furthermore, hydrogen bonding between the nitrogen atom in the heteroaromatic ring and the hydroxyl group of Ser20 or Tyr75 of the receptor appears to be important for the potent expression of biological activity.