Toggle switch residues control allosteric transitions in bacterial adhesins by participating in a concerted repacking of the protein core

Critical molecular events that control conformational transitions in most allosteric proteins are ill-defined. The mannose-specific FimH protein of Escherichia coli is a prototypic bacterial adhesin that switches from an 'inactive' low-affinity state (LAS) to an 'active' high-affinity state (HAS) conformation allosterically upon mannose binding and mediates shear-dependent catch bond adhesion. Here we identify a novel type of antibody that acts as a kinetic trap and prevents the transition between conformations in both directions. Disruption of the allosteric transitions significantly slows FimH's ability to associate with mannose and blocks bacterial adhesion under dynamic conditions. FimH residues critical for antibody binding form a compact epitope that is located away from the mannose-binding pocket and is structurally conserved in both states. A larger antibody-FimH contact area is identified by NMR and contains residues Leu-34 and Val-35 that move between core-buried and surface-exposed orientations in opposing directions during the transition. Replacement of Leu-34 with a charged glutamic acid stabilizes FimH in the LAS conformation and replacement of Val-35 with glutamic acid traps FimH in the HAS conformation. The antibody is unable to trap the conformations if Leu-34 and Val-35 are replaced with a less bulky alanine. We propose that these residues act as molecular toggle switches and that the bound antibody imposes a steric block to their reorientation in either direction, thereby restricting concerted repacking of side chains that must occur to enable the conformational transition. Residues homologous to the FimH toggle switches are highly conserved across a diverse family of fimbrial adhesins. Replacement of predicted switch residues reveals that another E. coli adhesin, galactose-specific FmlH, is allosteric and can shift from an inactive to an active state. Our study shows that allosteric transitions in bacterial adhesins depend on toggle switch residues and that an antibody that blocks the switch effectively disables adhesive protein function.

Structure of 3-mercaptopropionic acid dioxygenase with a substrate analog reveals bidentate substrate binding at the iron center

Thiol dioxygenases are a subset of nonheme iron oxygenases that catalyze the formation of sulfinic acids from sulfhydryl-containing substrates and dioxygen. Among this class, cysteine dioxygenases (CDOs) and 3-mercaptopropionic acid dioxygenases (3MDOs) are the best characterized, and the mode of substrate binding for CDOs is well understood. However, the manner in which 3-mercaptopropionic acid (3MPA) coordinates to the nonheme iron site in 3MDO remains a matter of debate. A model for bidentate 3MPA coordination at the 3MDO Fe-site has been proposed on the basis of computational docking, whereas steady-state kinetics and EPR spectroscopic measurements suggest a thiolate-only coordination of the substrate. To address this gap in knowledge, we determined the structure of Azobacter vinelandii 3MDO (Av3MDO) in complex with the substrate analog and competitive inhibitor, 3-hydroxypropionic acid (3HPA). The structure together with DFT computational modeling demonstrates that 3HPA and 3MPA associate with iron as chelate complexes with the substrate-carboxylate group forming an additional interaction with Arg168 and the thiol bound at the same position as in CDO. A chloride ligand was bound to iron in the coordination site assigned as the O2-binding site. Supporting HYSCORE spectroscopic experiments were performed on the (3MPA/NO)-bound Av3MDO iron nitrosyl (S = 3/2) site. In combination with spectroscopic simulations and optimized DFT models, this work provides an experimentally verified model of the Av3MDO enzyme-substrate complex, effectively resolving a debate in the literature regarding the preferred substrate-binding denticity. These results elegantly explain the observed 3MDO substrate specificity, but leave unanswered questions regarding the mechanism of substrate-gated reactivity with dioxygen.

RMSD analysis of structures of the bacterial protein FimH identifies five conformations of its lectin domain

FimH is a bacterial adhesin protein located at the tip of Escherichia coli fimbria that functions to adhere bacteria to host cells. Thus, FimH is a critical factor in bacterial infections such as urinary tract infections and is of interest in drug development. It is also involved in vaccine development and as a model for understanding shear-enhanced catch bond cell adhesion. To date, over 60 structures have been deposited in the Protein Data Bank showing interactions between FimH and mannose ligands, potential inhibitors, and other fimbrial proteins. In addition to providing insights about ligand recognition and fimbrial assembly, these structures provide insights into conformational changes in the two domains of FimH that are critical for its function. To gain further insights into these structural changes, we have superposed FimH's mannose binding lectin domain in all these structures and categorized the structures into five groups of lectin domain conformers using RMSD as a metric. Many structures also include the pilin domain, which anchors FimH to the fimbriae and regulates the conformation and function of the lectin domain. For these structures, we have also compared the relative orientations of the two domains. These structural analyses enhance our understanding of the conformational changes associated with FimH ligand binding and domain-domain interactions, including its catch bond behavior through allosteric action of force in bacterial adhesion.

The Evolution of SlyA/RovA Transcription Factors from Repressors to Countersilencers in Enterobacteriaceae

Gene duplication and subsequent evolutionary divergence have allowed conserved proteins to develop unique roles. The MarR family of transcription factors (TFs) has undergone extensive duplication and diversification in bacteria, where they act as environmentally responsive repressors of genes encoding efflux pumps that confer resistance to xenobiotics, including many antimicrobial agents. We have performed structural, functional, and genetic analyses of representative members of the SlyA/RovA lineage of MarR TFs, which retain some ancestral functions, including repression of their own expression and that of divergently transcribed multidrug efflux pumps, as well as allosteric inhibition by aromatic carboxylate compounds. However, SlyA and RovA have acquired the ability to countersilence horizontally acquired genes, which has greatly facilitated the evolution of Enterobacteriaceae by horizontal gene transfer. SlyA/RovA TFs in different species have independently evolved novel regulatory circuits to provide the enhanced levels of expression required for their new role. Moreover, in contrast to MarR, SlyA is not responsive to copper. These observations demonstrate the ability of TFs to acquire new functions as a result of evolutionary divergence of both cis-regulatory sequences and in trans interactions with modulatory ligands.IMPORTANCE Bacteria primarily evolve via horizontal gene transfer, acquiring new traits such as virulence and antibiotic resistance in single transfer events. However, newly acquired genes must be integrated into existing regulatory networks to allow appropriate expression in new hosts. This is accommodated in part by the opposing mechanisms of xenogeneic silencing and countersilencing. An understanding of these mechanisms is necessary to understand the relationship between gene regulation and bacterial evolution. Here we examine the functional evolution of an important lineage of countersilencers belonging to the ancient MarR family of classical transcriptional repressors. We show that although members of the SlyA lineage retain some ancestral features associated with the MarR family, their cis-regulatory sequences have evolved significantly to support their new function. Understanding the mechanistic requirements for countersilencing is critical to understanding the pathoadaptation of emerging pathogens and also has practical applications in synthetic biology.

Identifying G protein-coupled receptor dimers from crystal packings

Dimers of G protein-coupled receptors (GPCRs) are believed to be important for signaling with their associated G proteins. Low-resolution electron microscopy has shown rhodopsin dimers in native retinal membranes, and CXCR4 dimers have been found in several different crystal structures. Evidence for dimers of other GPCRs is more indirect. An alternative to computational modeling studies is to search for parallel dimers in the packing environments of the reported crystal structures of GPCRs. Two major structural types of GPCR dimers exist (as predicted by others), but there is considerable structural variation within each cluster. The different structural variants described here might reflect different functional properties and should provide a range of model structures for computational and experimental examination.

A Streptavidin Binding Site Mutation Yields an Unexpected Result: An Ionized Asp128 Residue Is Not Essential for Strong Biotin Binding

We report a detailed study of a point mutation of the crucial binding site residue, D128, in the biotin-streptavidin complex. The conservative substitution, D128N, preserves the detailed structure observed for the wild-type complex but has an only minimal impact on biotin binding, even though previous experimental and computational studies suggested that a charged D128 residue was crucial for high-affinity binding. These results show clearly that the fundamental basis for streptavidin's extremely strong biotin binding affinity is more complex than assumed and illustrate some of the challenges that may arise when analyzing extremely strong ligand-protein binding interactions.

Comparison of epsilon- and delta-class glutathione S-transferases: the crystal structures of the glutathione S-transferases DmGSTE6 and DmGSTE7 from Drosophila melanogaster

Cytosolic glutathione transferases (GSTs) comprise a large family of enzymes with canonical structures that diverge functionally and structurally among mammals, invertebrates and plants. Whereas mammalian GSTs have been characterized extensively with regard to their structure and function, invertebrate GSTs remain relatively unstudied. The invertebrate GSTs do, however, represent potentially important drug targets for infectious diseases and agricultural applications. In addition, it is essential to fully understand the structure and function of invertebrate GSTs, which play important roles in basic biological processes. Invertebrates harbor delta- and epsilon-class GSTs, which are not found in other organisms. Drosophila melanogaster GSTs (DmGSTs) are likely to contribute to detoxication or antioxidative stress during development, but they have not been fully characterized. Here, the structures of two epsilon-class GSTs from Drosophila, DmGSTE6 and DmGSTE7, are reported at 2.1 and 1.5 Å resolution, respectively, and are compared with other GSTs to identify structural features that might correlate with their biological functions. The structures of DmGSTE6 and DmGSTE7 are remarkably similar; the structures do not reveal obvious sources of the minor functional differences that have been observed. The main structural difference between the epsilon- and delta-class GSTs is the longer helix (A8) at the C-termini of the epsilon-class enzymes.

Crystal and molecular structure of the analgesic tetrapeptide, L-Phe-L-Leu-L-Pro-L-Ser

The tetrapeptide, L-Phe-L-Leu-L-Pro-L-Ser (FLPS), alleviates pain in a rat model of post-surgery pain. The crystal structure of the tetrapeptide is solved at high resolution (0.54 Å). The asymmetric unit contains two FLPS molecules, one Zn ion, and four molecules of water with a formula of [Zn(C23H33N4O6)2(H2O)4]. Each Zn ion is octahedrally coordinated with Phe and Ser residues from four peptide molecules [2N+4O]. The linking of Phe and Ser residues of one FLPS molecule to three other FLPS molecules by Zn ion forms a complex consisting of chains of metal ions and FLPS molecules oriented along the b axis. Analysis of molecular packing reveals the coexistence of two FLPS conformers in the same crystal. The crystallographic parameters for [Zn(C23H33N4O6)2(H2O)4] are as follows: space group P21 21 21 , a = 9.8698(2) Å, α = 90°, b = 20.1844(4) Å, β = 90°, c = 25.9302(6) Å, γ = 90°. Volume = 5165.71(19) Å(3), Z = 4, density (calc) = 1.364 Mg/cm(3), and agreement factor R1 = 4.13%.

The structure of rice weevil pectin methylesterase

Rice weevils (Sitophilus oryzae) use a pectin methylesterase (EC 3.1.1.11), along with other enzymes, to digest cell walls in cereal grains. The enzyme is a right-handed β-helix protein, but is circularly permuted relative to plant and bacterial pectin methylesterases, as shown by the crystal structure determination reported here. This is the first structure of an animal pectin methylesterase. Diffraction data were collected to 1.8 Å resolution some time ago for this crystal form, but structure solution required the use of molecular-replacement techniques that have been developed and similar structures that have been deposited in the last 15 years. Comparison of the structure of the rice weevil pectin methylesterase with that from Dickeya dandantii (formerly Erwinia chrysanthemi) indicates that the reaction mechanisms are the same for the insect, plant and bacterial pectin methylesterases. The similarity of the structure of the rice weevil enzyme to the Escherichia coli lipoprotein YbhC suggests that the evolutionary origin of the rice weevil enzyme was a bacterial lipoprotein, the gene for which was transferred to a primitive ancestor of modern weevils and other Curculionidae. Structural comparison of the rice weevil pectin methylesterase with plant and bacterial enzymes demonstrates that the rice weevil protein is circularly permuted relative to the plant and bacterial molecules.

E2~Ub conjugates regulate the kinase activity of Shigella effector OspG during pathogenesis

Pathogenic bacteria introduce effector proteins directly into the cytosol of eukaryotic cells to promote invasion and colonization. OspG, a Shigella spp. effector kinase, plays a role in this process by helping to suppress the host inflammatory response. OspG has been reported to bind host E2 ubiquitin-conjugating enzymes activated with ubiquitin (E2~Ub), a key enzyme complex in ubiquitin transfer pathways. A co-crystal structure of the OspG/UbcH5c~Ub complex reveals that complex formation has important ramifications for the activity of both OspG and the UbcH5c~Ub conjugate. OspG is a minimal kinase domain containing only essential elements required for catalysis. UbcH5c~Ub binding stabilizes an active conformation of the kinase, greatly enhancing OspG kinase activity. In contrast, interaction with OspG stabilizes an extended, less reactive form of UbcH5c~Ub. Recognizing conserved E2 features, OspG can interact with at least ten distinct human E2s~Ub. Mouse oral infection studies indicate that E2~Ub conjugates act as novel regulators of OspG effector kinase function in eukaryotic host cells.

Structural consequences of cutting a binding loop: two circularly permuted variants of streptavidin

Circular permutation of streptavidin was carried out in order to investigate the role of a main-chain amide in stabilizing the high-affinity complex of the protein and biotin. Mutant proteins CP49/48 and CP50/49 were constructed to place new N-termini at residues 49 and 50 in a flexible loop involved in stabilizing the biotin complex. Crystal structures of the two mutants show that half of each loop closes over the binding site, as observed in wild-type streptavidin, while the other half adopts the open conformation found in the unliganded state. The structures are consistent with kinetic and thermodynamic data and indicate that the loop plays a role in enthalpic stabilization of the bound state via the Asn49 amide-biotin hydrogen bond. In wild-type streptavidin, the entropic penalties of immobilizing a flexible portion of the protein to enhance binding are kept to a manageable level by using a contiguous loop of medium length (six residues) which is already constrained by its anchorage to strands of the β-barrel protein. A molecular-dynamics simulation for CP50/49 shows that cleavage of the binding loop results in increased structural fluctuations for Ser45 and that these fluctuations destabilize the streptavidin-biotin complex.

Second-contact shell mutation diminishes streptavidin-biotin binding affinity through transmitted effects on equilibrium dynamics

We report a point mutation in the second contact shell of the high-affinity streptavidin-biotin complex that appears to reduce binding affinity through transmitted effects on equilibrium dynamics. The Y54F streptavidin mutation causes a 75-fold loss of binding affinity with 73-fold faster dissociation, a large loss of binding enthalpy (ΔΔH = 3.4 kcal/mol at 37 °C), and a small gain in binding entropy (TΔΔS = 0.7 kcal/mol). The removed Y54 hydroxyl is replaced by a water molecule in the bound structure, but there are no observable changes in structure in the first contact shell and no additional changes surrounding the mutation. Molecular dynamics simulations reveal a large increase in the atomic fluctuation amplitudes for W79, a key biotin contact residue, compared to the fluctuation amplitudes in the wild-type. The increased W79 atomic fluctuation amplitudes are caused by loss of water-mediated hydrogen bonds between the Y54 hydroxyl group and peptide backbone atoms in and near W79. We propose that the increased atomic fluctuation amplitudes diminish the integrity of the W79-biotin interaction and represents a loosening of the "tryptophan collar" that is critical to the slow dissociation and high affinity of streptavidin-biotin binding. These results illustrate how changes in protein dynamics distal to the ligand binding pocket can have a profound impact on ligand binding, even when equilibrium structure is unperturbed.

Streptavidin and its biotin complex at atomic resolution

Atomic resolution crystallographic studies of streptavidin and its biotin complex have been carried out at 1.03 and 0.95 Å, respectively. The wild-type protein crystallized with a tetramer in the asymmetric unit, while the crystals of the biotin complex contained two subunits in the asymmetric unit. Comparison of the six subunits shows the various ways in which the protein accommodates ligand binding and different crystal-packing environments. Conformational variation is found in each of the polypeptide loops connecting the eight strands in the β-sandwich subunit, but the largest differences are found in the flexible binding loop (residues 45-52). In three of the unliganded subunits the loop is in an `open' conformation, while in the two subunits binding biotin, as well as in one of the unliganded subunits, this loop `closes' over the biotin-binding site. The `closed' loop contributes to the protein's high affinity for biotin. Analysis of the anisotropic displacement parameters included in the crystallographic models is consistent with the variation found in the loop structures and the view that the dynamic nature of the protein structure contributes to the ability of the protein to bind biotin so tightly.

The bacterial fimbrial tip acts as a mechanical force sensor

The subunits that constitute the bacterial adhesive complex located at the tip of the fimbria form a hook-chain that acts as a rapid force-sensitive anchor at high flow.

There is increasing evidence that the catch bond mechanism, where binding becomes stronger under tensile force, is a common property among non-covalent interactions between biological molecules that are exposed to mechanical force in vivo. Here, by using the multi-protein tip complex of the mannose-binding type 1 fimbriae of Escherichia coli, we show how the entire quaternary structure of the adhesive organella is adapted to facilitate binding under mechanically dynamic conditions induced by flow. The fimbrial tip mediates shear-dependent adhesion of bacteria to uroepithelial cells and demonstrates force-enhanced interaction with mannose in single molecule force spectroscopy experiments. The mannose-binding, lectin domain of the apex-positioned adhesive protein FimH is docked to the anchoring pilin domain in a distinct hooked manner. The hooked conformation is highly stable in molecular dynamics simulations under no force conditions but permits an easy separation of the domains upon application of an external tensile force, allowing the lectin domain to switch from a low- to a high-affinity state. The conformation between the FimH pilin domain and the following FimG subunit of the tip is open and stable even when tensile force is applied, providing an extended lever arm for the hook unhinging under shear. Finally, the conformation between FimG and FimF subunits is highly flexible even in the absence of tensile force, conferring to the FimH adhesin an exploratory function and high binding rates. The fimbrial tip of type 1 Escherichia coli is optimized to have a dual functionality: flexible exploration and force sensing. Comparison to other structures suggests that this property is common in unrelated bacterial and eukaryotic adhesive complexes that must function in dynamic conditions.

Noncovalent biological interactions are commonly subjected to mechanical force, particularly when they are involved in adhesion or cytoskeletal movements. While one might expect mechanical force to break these interactions, some of them form so-called catch bonds that lock on harder under force, like a nanoscale finger-trap. In this study, we show that the catch-bond forming adhesive protein FimH, which is located at the tip of E. coli fimbriae, allows bacteria to bind to urinary epithelial cells in a shear-dependent manner; that is, they bind at high but not at low flow. We show that isolated fimbrial tips, consisting of elongated protein complexes with FimH at the apex, reproduce this behavior in vitro. Our molecular dynamics simulations of the fimbrial tip structure show that FimH is shaped like a hook that is normally rigid but opens under force, causing structural changes that lead to firm anchoring of the bacteria on the surface. In contrast, the more distal adaptor proteins of the fimbrial tip create a flexible connection of FimH to the rigid fimbria, enhancing the ability of the adhesin to move into position and form bonds with mannose on the surface. We suggest that the entire tip complex forms a hook-chain, ideal for rapid and stable anchoring in flow.

A distal point mutation in the streptavidin-biotin complex preserves structure but diminishes binding affinity: experimental evidence of electronic polarization effects?

We have identified a distal point mutation in streptavidin that causes a 1000-fold reduction in biotin binding affinity without disrupting the equilibrium complex structure. The F130L mutation creates a small cavity occupied by a water molecule; however, all neighboring side chain positions are preserved, and protein-biotin hydrogen bonds are unperturbed. Molecular dynamics simulations reveal a reduced mobility of biotin binding residues but no observable destabilization of protein-ligand interactions. Our combined structural and computational studies suggest that the additional water molecule may affect binding affinity through an electronic polarization effect that impacts the highly cooperative hydrogen bonding network in the biotin binding pocket.

Donor strand exchange and conformational changes during E. coli fimbrial formation

Fimbriae and pili are macromolecular structures on the surface of Gram negative bacteria that are important for cellular adhesion. A 2.7Å resolution crystal structure of a complex of Escherichia coli fimbrial proteins containing FimH, FimG, FimF, and FimC provides the most complete model to date for the arrangement of subunits assembled in the native structure. The first three proteins form the tip of the fimbriae while FimC is the chaperone protein involved in the usher/chaperone assembly process. The subunits interact through donor strand complementation where a β-strand from a subunit completes the β-sandwich structure of the neighboring subunit or domain closer to the tip of the fimbria. The function of FimC is to provide a surrogate donor strand before delivery of each subunit to the FimD usher and the growing fimbria. Comparison of the subunits in this structure and their chaperone-bound complexes show that the two FimH domains change their relative orientation and position in forming the tip structure. Also, the non-chaperone subunits undergo a conformational change in their first β-strand when the chaperone is replaced by the native donor strand. Some residues move as much as 14Å in the process. This structural shift has not been noted in structural studies of other bacterial adhesion sub-structures assembled via donor strand complementation. The domains undergo a significant structural change in the donor strand binding groove during fimbrial assembly, and this likely plays a role in determining the specificity of subunit-subunit interactions among the fimbrial proteins.

Atomic resolution studies of carbonic anhydrase II

Carbonic anhydrase has been well studied structurally and functionally owing to its importance in respiration. A large number of X-ray crystallographic structures of carbonic anhydrase and its inhibitor complexes have been determined, some at atomic resolution. Structure determination of a sulfonamide-containing inhibitor complex has been carried out and the structure was refined at 0.9 A resolution with anisotropic atomic displacement parameters to an R value of 0.141. The structure is similar to those of other carbonic anhydrase complexes, with the inhibitor providing a fourth nonprotein ligand to the active-site zinc. Comparison of this structure with 13 other atomic resolution (higher than 1.25 A) isomorphous carbonic anhydrase structures provides a view of the structural similarity and variability in a series of crystal structures. At the center of the protein the structures superpose very well. The metal complexes superpose (with only two exceptions) with standard deviations of 0.01 A in some zinc-protein and zinc-ligand bond lengths. In contrast, regions of structural variability are found on the protein surface, possibly owing to flexibility and disorder in the individual structures, differences in the chemical and crystalline environments or the different approaches used by different investigators to model weak or complicated electron-density maps. These findings suggest that care must be taken in interpreting structural details on protein surfaces on the basis of individual X-ray structures, even if atomic resolution data are available.

Substrate specificity combined with stereopromiscuity in glutathione transferase A4-4-dependent metabolism of 4-hydroxynonenal

Conjugation to glutathione (GSH) by glutathione transferase A4-4 (GSTA4-4) is a major route of elimination for the lipid peroxidation product 4-hydroxynonenal (HNE), a toxic compound that contributes to numerous diseases. Both enantiomers of HNE are presumed to be toxic, and GSTA4-4 has negligible stereoselectivity toward them, despite its high catalytic chemospecificity for alkenals. In contrast to the highly flexible, and substrate promiscuous, GSTA1-1 isoform that has poor catalytic efficiency with HNE, GSTA4-4 has been postulated to be a rigid template that is preorganized for HNE metabolism. However, the combination of high substrate chemoselectivity and low substrate stereoselectivity is intriguing. The mechanism by which GSTA4-4 achieves this combination is important, because it must metabolize both enantiomers of HNE to efficiently detoxify the biologically formed mixture. The crystal structures of GSTA4-4 and an engineered variant of GSTA1-1 with high catalytic efficiency toward HNE, cocrystallized with a GSH-HNE conjugate analogue, demonstrate that GSTA4-4 undergoes no enantiospecific induced fit; instead, the active site residue Arg15 is ideally located to interact with the 4-hydroxyl group of either HNE enantiomer. The results reveal an evolutionary strategy for achieving biologically useful stereopromiscuity toward a toxic racemate, concomitant with high catalytic efficiency and substrate specificity toward an endogenously formed toxin.

Structural analysis of a glutathione transferase A1-1 mutant tailored for high catalytic efficiency with toxic alkenals

The specificity of human glutathione transferase (GST) A1-1 is drastically altered to favor alkenal substrates in the GIMFhelix mutant designed to mimic first-sphere interactions utilized by GSTA4-4. This redesign serves as a model for improving our understanding of the structural determinants that contribute to the distinct specificities of alpha class GSTs. Herein we report the first crystal structures of GIMFhelix, both in complex with GSH and in apo form at 1.98 and 2.38 A resolution. In contrast to the preorganized hydrophobic binding pocket that accommodates alkenals in GSTA4-4, GSTA1-1 includes a dynamic alpha9 helix that undergoes a ligand-dependent localization to complete the active site. Comparisons of the GIMFhelix structures with previously reported structures show a striking similarity with the GSTA4-4 active site obtained within an essentially GSTA1-1 scaffold and reveal the alpha9 helix assumes a similar localized structure regardless of active site occupancy in a manner resembling that of GSTA4-4. However, we cannot fully account for all the structural elements important in GSTA4-4 within the mutant's active site. The contribution of Phe10 to the Tyr212-Phe10-Phe220 network prevents complete C-terminal closure and demonstrates that the presence of Phe10 within the context of a GSTA4-4-like active site may ultimately hinder Phe220, a key C-terminal residue, from effectively contributing to the active site. In total, these results illustrate the remaining structural differences presumably reflected in the previously reported catalytic efficiencies of GIMFhelix and GSTA4-4 and emphasize the F10P mutation as being necessary to completely accomplish the transformation to a highly specific GST from the more promiscuous GSTA1-1 enzyme.

Dynamics of the streptavidin-biotin complex in solution and in its crystal lattice: distinct behavior revealed by molecular simulations

We present a 250 ns simulation of the wild-type, biotin-liganded streptavidin tetramer in the solution phase and compare the trajectory to two previously published simulations of the protein in its crystal lattice. By performing both types of simulations, we are able to interpret the protein's behavior in solution in the context of its X-ray structure. We find that the rate of conformational sampling is increased in solution over the lattice environment, although the relevant conformational space in solution is also much larger, as indicated by overall fluctuations in the positions of backbone atoms. We also compare the distributions of chi1 angles sampled by side chains exposed to solvent in the lattice and in the solution phase, obtaining overall good agreement between the distributions obtained in our most rigorous lattice simulation and the crystallographic chi1 angles. We observe changes in the chi1 distributions in the solution phase, and note an apparent progression of the distributions as the environment changes from a tightly packed lattice filled with crystallization media to a bath of pure water. Finally, we examine the interaction of biotin and streptavidin in each simulation, uncovering a possible alternate conformation of the biotin carboxylate tail. We also note that a hydrogen bond observed to break transiently in previous solution-phase simulations is predominantly broken in this much longer solution-phase trajectory; in the lattice simulations, the lattice environment appears to help maintain the hydrogen bond, but more sampling will be needed to confirm whether the simulation model truly gives good agreement with the X-ray data in the lattice simulations. We expect that pairing solution-phase biomolecular simulations with crystal lattice simulations will help to validate simulation models and improve the interpretation of experimentally determined structures.

Simulations of a protein crystal: explicit treatment of crystallization conditions links theory and experiment in the streptavidin-biotin complex

A 250 ns molecular dynamics simulation of the biotin-liganded streptavidin crystal lattice, including cryoprotectant molecules and crystallization salts, is compared to a 250 ns simulation of the lattice solvated with pure water. The simulation using detailed crystallization conditions preserves the initial X-ray structure better than the simulation using pure water, even though the protein molecules display comparable mobility in either simulation. Atomic fluctuations computed from the simulation with crystallization conditions closely reproduce fluctuations derived from experimental temperature factors (correlation coefficient of 0.88, omitting two N-terminal residues with very high experimental B-factors). In contrast, fluctuations calculated from the simulation with pure water were less accurate, particularly for two of the streptavidin loops exposed to solvent in the crystal lattice. Finally, we obtain good agreement between the water and cryoprotectant densities obtained from the simulated crystallization conditions and the electron density due to solvent molecules in the X-ray structure. Our results suggest that detailed lattice simulations with realistic crystallization conditions can be used to assess potential function parameters, validate simulation protocols, and obtain valuable insights that solution-phase simulations do not easily provide. We anticipate that this will prove to be a powerful strategy for molecular dynamics simulations of biomolecules.

Alternative models for two crystal structures of bovine rhodopsin

The space-group symmetry of two crystal forms of rhodopsin (PDB codes 1gzm and 2j4y; space group P3(1)) can be re-interpreted as hexagonal (space group P6(4)). Two molecules of the G protein-coupled receptor are present in the asymmetric unit in the trigonal models. However, the noncrystallographic twofold axes parallel to the c axis can be treated as crystallographic symmetry operations in the hexagonal space group. This halves the asymmetric unit and makes all of the protein molecules equivalent in these structures. Corrections for merohedral twinning were also applied in the refinement in the higher symmetry space group for one of the structures (2j4y).

Crystal structure of the BARD1 ankyrin repeat domain and its functional consequences

BARD1 is the constitutive nuclear partner to the breast and ovarian cancer-specific tumor suppressor BRCA1. Together, they form a heterodimeric complex responsible for maintaining genomic stability through nuclear functions involving DNA damage signaling and repair, transcriptional regulation, and cell cycle control. We report the 2.0A structure of the BARD1 ankyrin repeat domain. The structure includes four ankyrin repeats with a non-canonical C-terminal capping ankyrin repeat and a well ordered extended loop preceding the first repeat. Conserved surface features show an acidic patch and an acidic pocket along the surface typically used by ankyrin repeat domains for binding cognate proteins. We also demonstrate that two reported mutations, N470S and V507M, in the ankyrin repeat domain do not result in observable structural defects. These results provide a structural basis for exploring the biological function of the ankyrin repeat domain and for modeling BARD1 isoforms.

Reprint of "Crystal packing analysis of Rhodopsin crystals" [J. Struct. Biol. 158 (2007) 455-462]

Oligomerization has been proposed as one of several mechanisms to regulate the activity of G protein-coupled receptors (GPCRs), but little is known about the structure of GPCR oligomers. Crystallographic analyses of two new crystal forms of rhodopsin reveal an interaction surface which may be involved in the formation of functional dimers or oligomers. New crystallization conditions lead to the formation of two crystal forms with similar rhodopsin-rhodopsin interactions, but changes in the crystal lattice are induced by the addition of different surfactant additives. However, the intermolecular interactions between rhodopsin molecules in these crystal structures may reflect the contacts necessary for the maintenance of dimers or oligomers in rod outer segment membranes. Similar contacts may assist in the formation of dimers or oligomers in other GPCRs as well. These new dimers are compared with other models proposed by crystallography or EM and AFM studies. The inter-monomer surface contacts are different for each model, but several of these models coincide in implicating helix I, II, and H-8 as contributors to the main contact surface stabilizing the dimers.

Binding of Dr adhesins of Escherichia coli to carcinoembryonic antigen triggers receptor dissociation

Carcinoembryonic antigen (CEA)-related cell adhesion molecules (CEACAMs) are host receptors for the Dr family of adhesins of Escherichia coli. To define the mechanism for binding of Dr adhesins to CEACAM receptors, we carried out structural studies on the N-terminal domain of CEA and its complex with the Dr adhesin. The crystal structure of CEA reveals a dimer similar to other dimers formed by receptors with IgV-like domains. The structure of the CEA/Dr adhesin complex is proposed based on NMR spectroscopy and mutagenesis data in combination with biochemical characterization. The Dr adhesin/CEA interface overlaps appreciably with the region responsible for CEA dimerization. Binding kinetics, mutational analysis and spectroscopic examination of CEA dimers suggest that Dr adhesins can dissociate CEA dimers prior to the binding of monomeric forms. Our conclusions include a plausible mechanism for how E. coli, and perhaps other bacterial and viral pathogens, exploit CEACAMs. The present structure of the complex provides a powerful tool for the design of novel inhibitory strategies to treat E. coli infections.

An alternate description of two crystal structures of phospholipase A2 from Bungarus caeruleus

Reinterpretations of the space-group symmetry are reported for two crystal structures of phospholipase A(2) isoforms (PDB codes 1u4j and 1g2x). The two structures reported in space groups R3 and C2 are isomorphous with a third isoform with space group R32 (PDB code 1fe5). The original structure reports were interpreted in terms of different oligomeric forms of the isoforms, but these conclusions are not supported by the isomorphous structures.

Crystal packing analysis of Rhodopsin crystals

Oligomerization has been proposed as one of several mechanisms to regulate the activity of G protein-coupled receptors (GPCRs), but little is known about the structure of GPCR oligomers. Crystallographic analyses of two new crystal forms of rhodopsin reveal an interaction surface which may be involved in the formation of functional dimers or oligomers. New crystallization conditions lead to the formation of two crystal forms with similar rhodopsin-rhodopsin interactions, but changes in the crystal lattice are induced by the addition of different surfactant additives. However, the intermolecular interactions between rhodopsin molecules in these crystal structures may reflect the contacts necessary for the maintenance of dimers or oligomers in rod outer segment membranes. Similar contacts may assist in the formation of dimers or oligomers in other GPCRs as well. These new dimers are compared with other models proposed by crystallography or EM and AFM studies. The inter-monomer surface contacts are different for each model, but several of these models coincide in implicating helix I, II, and H-8 as contributors to the main contact surface stabilizing the dimers.

Improvements in G protein-coupled receptor purification yield light stable rhodopsin crystals

G protein-coupled receptors (GPCRs) represent the largest family of transmembrane signaling proteins and are the target of approximately half of all therapeutic agents. Agonist ligands bind their cognate GPCRs stabilizing the active conformation that is competent to bind G proteins, thus initiating a cascade of intracellular signaling events leading to modification of the cell activity. Despite their biomedical importance, the only known GPCR crystal structures are those of inactive rhodopsin forms. In order to understand how GPCRs are able to transduce extracellular signals across the plasma membrane, it is critical to determine the structure of these receptors in their ligand-bound, active state. Here, we report a novel combination of purification procedures that allowed the crystallization of rhodopsin in two new crystal forms and can be applicable to the purification and crystallization of other membrane proteins. Importantly, these new crystals are stable upon photoactivation and the preliminary X-ray diffraction analysis of both photoactivated and ground state rhodopsin crystals are also reported.

Crystal structure of a photoactivated deprotonated intermediate of rhodopsin

The changes that lead to activation of G protein-coupled receptors have not been elucidated at the structural level. In this work we report the crystal structures of both ground state and a photoactivated deprotonated intermediate of bovine rhodopsin at a resolution of 4.15 A. In the photoactivated state, the Schiff base linking the chromophore and Lys-296 becomes deprotonated, reminiscent of the G protein-activating state, metarhodopsin II. The structures reveal that the changes that accompany photoactivation are smaller than previously predicted for the metarhodopsin II state and include changes on the cytoplasmic surface of rhodopsin that possibly enable the coupling to its cognate G protein, transducin. Furthermore, rhodopsin forms a potentially physiologically relevant dimer interface that involves helices I, II, and 8, and when taken with the prior work that implicates helices IV and V as the physiological dimer interface may account for one of the interfaces of the oligomeric structure of rhodopsin seen in the membrane by atomic force microscopy. The activation and oligomerization models likely extend to the majority of other G protein-coupled receptors.

Crystal structure and mutational analysis of the DaaE adhesin of Escherichia coli

DaaE is a member of the Dr adhesin family of Escherichia coli, members of which are associated with diarrhea and urinary tract infections. A receptor for Dr adhesins is the cell surface protein, decay-accelerating factor (DAF). We have carried out a functional analysis of Dr adhesins, as well as mutagenesis and crystallographic studies of DaaE, to obtain detailed molecular information about interactions of Dr adhesins with their receptors. The crystal structure of DaaE has been solved at 1.48 A resolution. Trimers of the protein are found in the crystal, as has been the case for other Dr adhesins. Naturally occurring variants and directed mutations in DaaE have been generated and analyzed for their ability to bind DAF. Mapping of the mutation sites onto the DaaE molecular structure shows that several of them contribute to a contiguous surface that is likely the primary DAF-binding site. The DAF-binding properties of purified fimbriae and adhesin proteins from mutants and variants correlated with the ability of bacteria expressing these proteins to bind to human epithelial cells in culture. DaaE, DraE, AfaE-III, and AfaE-V interact with complement control protein (CCP) domains 2-4 of DAF, and analysis of the ionic strength dependence of their binding indicates that the intermolecular interactions are highly electrostatic in nature. The adhesins AfaE-I and NfaE-2 bind to CCP-3 and CCP-4 of DAF, and electrostatic interactions contribute significantly less to these interactions. These observations are consistent with structural predictions for these Dr variants and also suggest a role for the positively charged region linking CCP-2 and CCP-3 of DAF in electrostatic Dr adhesin-DAF interactions.

The high-resolution structure of (+)-epi-biotin bound to streptavidin

(+)-Epi-biotin differs from (+)-biotin in the configuration of the chiral center at atom C2. This could lead to a difference in the mode of binding of (+)-epi-biotin to streptavidin, a natural protein receptor for (+)-biotin. Diffraction data were collected to a maximum of 0.85 Angstrom resolution for structural analysis of the complex of streptavidin with a sample of (+)-epi-biotin and refinement was carried out at both 1.0 and 0.85 Angstrom resolution. The structure determination shows a superposition of two ligands in the binding site, (+)-biotin and (+)-epi-biotin. The molecules overlap in the model for the complex except for the position of S1 in the tetrahydrothiophene ring. Differences in the conformation of the ring permits binding of each molecule to streptavidin with little observable difference in the protein structures at this high resolution.

Functional and structural characterization of rhodopsin oligomers

A major question in G protein-coupled receptor signaling concerns the quaternary structure required for signal transduction. Do these transmembrane receptors function as monomers, dimers, or larger oligomers? We have investigated the oligomeric state of the model G protein-coupled receptor rhodopsin (Rho), which absorbs light and initiates a phototransduction-signaling cascade that forms the basis of vision. In this study, different forms of Rho were isolated using gel filtration techniques in mild detergents, including n-dodecyl-beta-D-maltoside, n-tetradecyl-beta-D-maltoside, and n-hexadecyl-beta-D-maltoside. The quaternary structure of isolated Rho was determined by transmission electron microscopy, demonstrating that in micelles containing n-dodecyl-beta-D-maltoside, Rho exists as a mixture of monomers and dimers whereas in n-tetradecyl-beta-D-maltoside and n-hexadecyl-beta-D-maltoside Rho forms higher ordered structures. Especially in n-hexadecyl-beta-D-maltoside, most of the particles are present in tightly packed rows of dimers. The oligomerization of Rho seems to be important for interaction with its cognate G protein, transducin. Although the activated Rho (Meta II) monomer or dimers are capable of activating the G protein, transducin, the activation process is much faster when Rho exists as organized dimers. Our studies provide direct comparisons between signaling properties of Meta II in different quaternary complexes.

Crystallographic Analysis of a Full-length Streptavidin with Its C-terminal Polypeptide Bound in the Biotin Binding Site

The structure of a full-length streptavidin has been determined at 1.7 A resolution and shows that the 20 residue extension at the C terminus forms a well-ordered polypeptide loop on the surface of the tetramer. Residues 150-153 of the extension are bound to the ligand-binding site, possibly competing with exogenous ligands. The binding mode of these residues is compared with that of biotin and peptidic ligands. The observed structure helps to rationalize the observations that full-length mature streptavidin binds biotinylated macromolecules with reduced affinity.

Cooperative hydrogen bond interactions in the streptavidin-biotin system

The thermodynamic and structural cooperativity between the Ser45- and D128-biotin hydrogen bonds was measured by calorimetric and X-ray crystallographic studies of the S45A/D128A double mutant of streptavidin. The double mutant exhibits a binding affinity approximately 2x10(7) times lower than that of wild-type streptavidin at 25 degrees C. The corresponding reduction in binding free energy (DeltaDeltaG) of 10.1 kcal/mol was nearly completely due to binding enthalpy losses at this temperature. The loss of binding affinity is 11-fold greater than that predicted by a linear combination of the single-mutant energetic perturbations (8.7 kcal/mol), indicating that these two mutations interact cooperatively. Crystallographic characterization of the double mutant and comparison with the two single mutant structures suggest that structural rearrangements at the S45 position, when the D128 carboxylate is removed, mask the true energetic contribution of the D128-biotin interaction. Taken together, the thermodynamic and structural analyses support the conclusion that the wild-type hydrogen bond between D128-OD and biotin-N2 is thermodynamically stronger than that between S45-OG and biotin-N1.

Anatomy of atrans–cispeptide transition during least-squares refinement of rubrerythrin

A detailed view is presented of the effects of one round of 20 cycles of restrained least-squares refinement of rubreythrin in which a trans peptide between Gly78 and Ile79 converts to a cis conformation automatically. While the omega angle for the peptide changes by nearly 180 degrees , the maximum shift in any atomic position is 1.32 Angstroms. The peptide converts by passing through a non-ideal structure containing a nearly linear C-N-C(alpha) bond angle. The overall motion is not possible for real or virtual molecular models with ideal bond lengths, angles and torsion angles. Strengthening the stereochemical bond length and bond angle restraints halts the structural change.

Rhodopsin: A Structural Primer for G-Protein Coupled Receptors

An overview of the rhodopsin crystal structure provides a structural basis for understanding the structures and functions of other G-protein coupled receptors (GPCRs). All of the structural details observed to date for rhodopsin will not necessarily carry over to other GPCRs, but major features such as the arrangement of the seven transmembrane helices, the retinal/ligand binding site, the D(E)RY and NPXXY sequence and structural motifs, and the bent helices are likely characteristics of the GPCRs most closely related to rhodopsin. A general view of these structural features is presented here.

Functional characterization of rhodopsin monomers and dimers in detergents

Rhodopsin (Rho) is a G protein-coupled receptor that initiates phototransduction in rod photoreceptors. High expression levels of Rho in the disc membranes of rod outer segments and the propensity of Rho to form higher oligomeric structures are evident from atomic force microscopy, transmission electron microscopy, and chemical cross-linking experiments. To explore the structural and functional properties of Rho in n-dodecyl-beta-maltoside, frequently used to purify heterologously expressed Rho and its mutants, we used gel filtration techniques, blue native gel electrophoresis, and functional assays. Here, we show that in micelles containing n-dodecyl-beta-maltoside at concentrations greater than 3 mM, Rho is present as a single monomer per detergent micelle. In contrast, in 12 mM 3-[(3-cholamidopropyl)dimethyl-ammonio]-1-propanesulfonate (CHAPS), micelles contain mostly dimeric Rho. The cognate G protein transducin (Gt) appears to have a preference for binding to the Rho dimer, and the complexes fall apart in the presence of guanosine 5'-3-O-(thio)triphosphate. Cross-linked Rho dimers release the chromophore at a slower rate than monomers and are much more resistant to heat denaturation. Both Rho(*) monomers and dimers are capable of activating Gt, and both of them are phosphorylated by Rho kinase. Rho expressed in HEK293 cells is also readily cross-linked by a bifunctional reagent. These studies provide an explanation of how detergent influences the oligomer-dimermonomer equilibrium of Rho and describe the functional characterization of Rho monomers and dimers in detergent.

A naturally occurring mutation of the opsin gene (T4R) in dogs affects glycosylation and stability of the G protein-coupled receptor

Rho (rhodopsin; opsin plus 11-cis-retinal) is a prototypical G protein-coupled receptor responsible for the capture of a photon in retinal photoreceptor cells. A large number of mutations in the opsin gene associated with autosomal dominant retinitis pigmentosa have been identified. The naturally occurring T4R opsin mutation in the English mastiff dog leads to a progressive retinal degeneration that closely resembles human retinitis pigmentosa caused by the T4K mutation in the opsin gene. Using genetic approaches and biochemical assays, we explored the properties of the T4R mutant protein. Employing immunoaffinity-purified Rho from affected RHO(T4R/T4R) dog retina, we found that the mutation abolished glycosylation at Asn(2), whereas glycosylation at Asn(15) was unaffected, and the mutant opsin localized normally to the rod outer segments. Moreover, we found that T4R Rho(*) lost its chromophore faster as measured by the decay of meta-rhodopsin II and that it was less resistant to heat denaturation. Detergent-solubilized T4R opsin regenerated poorly and interacted abnormally with the G protein transducin (G(t)). Structurally, the mutation affected mainly the "plug" at the intradiscal (extracellular) side of Rho, which is possibly responsible for protecting the chromophore from the access of bulk water. The T4R mutation may represent a novel molecular mechanism of degeneration where the unliganded form of the mutant opsin exerts a detrimental effect by losing its structural integrity.

Evolutionary analysis of rhodopsin and cone pigments: connecting the three-dimensional structure with spectral tuning and signal transfer

Extensive sequence data and structural sampling of expressed proteins from different species lead to the idea that entire molecules or specific domain folds belong to large superfamilies of proteins. A subset of G protein-coupled receptors, one of the largest families involved in cellular signaling, rod and cone opsins are involved in phototransduction in photoreceptor cells. Here, the evolutionary analysis of opsin sequences and structures predicts key residues involved in the transmission of the signal from the binding site of the chromophore to the cytoplasmic surface and residues that are involved in the spectral tuning of opsins to short wavelengths of light.

Structural studies of hydrogen bonds in the high-affinity streptavidin-biotin complex: mutations of amino acids interacting with the ureido oxygen of biotin

An elaborate hydrogen-bonding network contributes to the tight binding of biotin to streptavidin. The specific energetic contributions of hydrogen bonds to the biotin ureido oxygen have previously been investigated by mapping the equilibrium and activation thermodynamic signatures of N23A, N23E, S27A, Y43A and Y43F site-directed mutants [Klumb et al. (1998), Biochemistry, 37, 7657-7663]. The crystal structures of these variants in the unbound and biotin-bound states provide structural insight into the energetic alterations and are described here. High (1.5-2.2 A) to atomic resolution (1.14 A) structures were obtained and structural models were refined to R values ranging from 0.12 to 0.20. The overall folding of streptavidin as described previously has not changed in any of the mutant structures. Major deviations such as side-chain shifts of residues in the binding site are observed only for the N23A and Y43A mutations. In none of the mutants is a systematic shift of biotin observed when one of the hydrogen-bonding partners to the ureido oxygen of biotin is removed. Recent thermodynamic studies report increases of DeltaDeltaG(o) of 5.0-14.6 kJ mol(-1) for these mutants with respect to the wild-type protein. The decreasing stabilities of the complexes of the mutants are discussed in terms of their structures.

Structural characterization and comparison of RGD cell-adhesion recognition sites engineered into streptavidin

The RGD (arginine-glycine-aspartic acid) sequence is found in several important extracellular matrix proteins and serves as an adhesion ligand for members of the integrin family of cell-surface receptors. This sequence and flanking residues from fibronectin or osteopontin have been engineered into an accessible surface loop of streptavidin to create two new streptavidin variants (FN-SA or OPN-SA, respectively) that bind cells through the alpha(v)beta(3) and/or alpha(5)beta(1) integrin receptors. Their crystal structures confirm the design and construction of the mutants and provide evidence about the conformational dynamics of the RGD loops. The loops in the isomorphous crystal structures are involved in crystal-packing interactions and this stabilizes their structures. Even so, the loop in OPN-SA is slightly disordered and two of the residues are not seen in difference electron-density maps. Comparison with other experimentally determined structures of RGD loops in cell-adhesion molecules shows that these loops occupy a large subset of conformational space. This is consistent with the view that RGD loops, at least those involved in cell adhesion, sample a number of structures dynamically, a few of which display high affinity for appropriate receptors.

G protein-coupled receptor rhodopsin: a prospectus

Rhodopsin is a retinal photoreceptor protein of bipartite structure consisting of the transmembrane protein opsin and a light-sensitive chromophore 11-cis-retinal, linked to opsin via a protonated Schiff base. Studies on rhodopsin have unveiled many structural and functional features that are common to a large and pharmacologically important group of proteins from the G protein-coupled receptor (GPCR) superfamily, of which rhodopsin is the best-studied member. In this work, we focus on structural features of rhodopsin as revealed by many biochemical and structural investigations. In particular, the high-resolution structure of bovine rhodopsin provides a template for understanding how GPCRs work. We describe the sensitivity and complexity of rhodopsin that lead to its important role in vision.

Crystal structures and increased stabilization of the protein G variants with switched folding pathways NuG1 and NuG2

We recently described two protein G variants (NuG1 and NuG2) with redesigned first hairpins that were almost twice as stable, folded 100-fold faster, and had a switched folding mechanism relative to the wild-type protein. To test the structural accuracy of our design algorithm and to provide insights to the dramatic changes in the kinetics and thermodynamics of folding, we have now determined the crystal structures of NuG1 and NuG2 to 1.8 A and 1.85 A, respectively. We find that they adopt hairpin structures that are closer to the computational models than to wild-type protein G; the RMSD of the NuG1 hairpin to the design model and the wild-type structure are 1.7 A and 5.1 A, respectively. The crystallographic B factor in the redesigned first hairpin of NuG1 is systematically higher than the second hairpin, suggesting that the redesigned region is somewhat less rigid. A second round of structure-based design yielded new variants of NuG1 and NuG2, which are further stabilized by 0.5 kcal/mole and 0.9 kcal/mole.

Crystal structure of rhodopsin: a template for cone visual pigments and other G protein-coupled receptors

The crystal structure of rhodopsin has provided the first three-dimensional molecular model for a G-protein-coupled receptor (GPCR). Alignment of the molecular model from the crystallographic structure with the helical axes seen in cryo-electron microscopic (cryo-EM) studies provides an opportunity to investigate the properties of the molecule as a function of orientation and location within the membrane. In addition, the structure provides a starting point for modeling and rational experimental approaches of the cone pigments, the GPCRs in cone cells responsible for color vision. Homology models of the cone pigments provide a means of understanding the roles of amino acid sequence differences that shift the absorption maximum of the retinal chromophore in the environments of different opsins.

1.3-A resolution structure of human glutathione S-transferase with S-hexyl glutathione bound reveals possible extended ligandin binding site

Cytosolic glutathione S-transferases (GSTs) play a critical role in xenobiotic binding and metabolism, as well as in modulation of oxidative stress. Here, the high-resolution X-ray crystal structures of homodimeric human GSTA1-1 in the apo form and in complex with S-hexyl glutathione (two data sets) are reported at 1.8, 1.5, and 1.3A respectively. At this level of resolution, distinct conformations of the alkyl chain of S-hexyl glutathione are observed, reflecting the nonspecific nature of the hydrophobic substrate binding site (H-site). Also, an extensive network of ordered water, including 75 discrete solvent molecules, traverses the open subunit-subunit interface and connects the glutathione binding sites in each subunit. In the highest-resolution structure, three glycerol moieties lie within this network and directly connect the amino termini of the glutathione molecules. A search for ligand binding sites with the docking program Molecular Operating Environment identified the ordered water network binding site, lined mainly with hydrophobic residues, suggesting an extended ligand binding surface for nonsubstrate ligands, the so-called ligandin site. Finally, detailed comparison of the structures reported here with previously published X-ray structures reveal a possible reaction coordinate for ligand-dependent conformational changes in the active site and the C-terminus.