Weakly polar interactions in proteins

The catalytic role of tyrosine 254 in flavocytochrome b2 (L-lactate dehydrogenase from baker's yeast). Comparison between the Y254F and Y254L mutant proteins

Flavocytochrome b2 catalyses the oxidation of L-lactate to pyruvate in yeast mitochondrial intermembrane space. Its flavoprotein domain is a member of a family of FMN-dependent 2-hydroxy-acid-oxidizing enzymes. Numerous solution studies suggest that the first step of the reaction consists of proton abstraction from lactate C2, leading to a carbanion that subsequently yields electrons to FMN. The crystal structure suggests that the enzyme base is His373, and that Tyr254 may be hydrogen bonded to the substrate hydroxyl. Studies carried out with the Y254F mutant [Dubois, J., Chapman, S.K., Mathews, F.S., Reid, G.A. & Lederer, F. (1990) Biochemistry 29, 6393-6400] showed that Tyr254 does not act as a base but stabilizes the transition state. As the mutation did not induce any change in substrate affinity, the question of the existence of the hydrogen bond in the Michaelis complex remained open. Similar results with glycolate oxidase, mutated at the same position, led to the suggestion that these enzymes actually operate via a hydride transfer mechanism [Macheroux, P., Kieweg, V., Massey, V., Soderlind, E., Stenberg, K. & Lindqvist, Y. (1993) Eur. J. Biochem. 213, 1047-1054]. In the present work, we have re-investigated the matter by analysing the properties of a Y254L mutant flavocytochrome b2, as well as the behaviour of the Y254F enzyme with two substrates other than lactate, and a series of inhibitors. The Y254L protein is less efficient with L-lactate than the wild-type enzyme by a factor of 500, but the substrate affinity is unchanged. In contrast, L-phenyllactate and mandelate, poor substrates (the latter acting more as an inhibitor), exhibit an increased affinity. In addition, the Y254L mutant enzyme is more efficient with phenyllactate than lactate as a substrate. In order to rationalize these observations, we have modelled phenyllactate and mandelate in the active site, using previously described modelling experiments with lactate as a starting point. The results indicate that mandelate cannot bind in an orientation allowing proton abstraction by His373, due to steric interference by the side chains of Ala198 and Leu230. It might possibly adopt a binding mode as proposed previously for lactate, which leads to a hydride transfer and with which the 198 and 230 side chains do not interfere. However, other researchers [Sinclair, R., Reid, G.A. & Chapman, S.K. (1998) Biochem. J. 333, 117-120] showed that A198G, L230A and A198G/L230A mutant enzymes exhibit a strongly improved mandelate dehydrogenase activity. These results indicate that relief of the steric crowding facilitates catalysis by enabling a better mandelate orientation at the active site, suggesting that its productive binding mode is similar to that proposed for lactate in the carbanion mechanism. The modelling studies therefore support the hypothesis of a carbanion mechanism for all substrates. In addition, we present the effect of the two mutations at position 254 on the binding of a number of competitive inhibitors (such as sulfite, D-lactate, propionate) and of inhibitors that are known to bind at the active site both when the flavin is oxidized and when it is in the semiquinone state (propionate, oxalate and L-lactate at high concentrations). Unexpectedly, the results indicate that the integrity of Tyr254 is necessary for the binding of these inhibitors at the semiquinone stage.

More hydrogen bonds for the (structural) biologist

Why does a given protein structure form and why is this structure stable? These fundamental biochemical questions remain fascinating and challenging problems because the physical bases of the forces that govern protein structure, stability and folding are still not well understood. Now, a general concept of hydrogen bonding in proteins is emerging. This concept involves not only N-H and O-H donor groups, but also C-H, and not only N and O as acceptor groups, but also pi-systems. We postulate that the incorporation of the entirety of these interactions leads to a more complete description of the problem, and that this could provide new perspectives and possibly new answers.

Non-hydrogen bond interactions involving the methionine sulfur atom

Of all the nonbonded interactions, hydrogen bond, because of its geometry involving polar atoms, is the most easily recognizable. Here we characterize two interactions involving the divalent sulfur of methionine (Met) residues that do not need any participation of proton. In one an oxygen atom of the main-chain carbonyl group or a carboxylate side chain is used. In another an aromatic atom interacting along the face of the ring is utilized. In these, the divalent sulfur behaves as an electrophile and the other electron-rich atom, a nucleophile. The stereochemistry of the interaction is such that the nucleophile tends to approach approximately along the extension of one of the covalent bonds to S. The nitrogen atom of histidine side chain is extensively used in these nonbonded contacts. There is no particular geometric pattern in the interaction of S with the edge of an aromatic ring, except when an N-H group in involved, which is found within 40 degrees from the perpendicular to the sulfide plane, thus defining the geometry of hydrogen bond interaction involving the sulfur atom. As most of the Met residues which partake in such stereospecific interactions are buried, these would be important for the stability of the protein core, and their incorporation in the binding site would be useful for molecular recognition and optimization of the site's affinity for partners (especially containing aromatic and heteroaromatic groups). Mutational studies aimed at replacing Met by other residues would benefit from the delineation of these interactions.

Folate-binding triggers the activation of folylpolyglutamate synthetase

Folic acid is an essential vitamin for normal cell growth, primarily through its central role in one-carbon metabolism. Folate analogs (antifolates) are targeted at the same reactions and are widely used as therapeutic drugs for cancer and bacterial infections. Effective retention of folates in cells and the efficacy of antifolate drugs both depend upon the addition of a polyglutamate tail to the folate or antifolate molecule by the enzyme folylpolyglutamate synthetase (FPGS). The reaction mechanism involves the ATP-dependent activation of the free carboxylate group on the folate molecule to give an acyl phosphate intermediate, followed by attack by the incoming L-glutamate substrate. FPGS shares a number of structural and mechanistic details with the bacterial cell wall ligases MurD, MurE and MurF, and these enzymes, along with FPGS, form a subfamily of the ADP-forming amide bond ligase family. High-resolution crystallographic analyses of binary and ternary complexes of Lactobacillus casei FPGS reveal that binding of the first substrate (ATP) is not sufficient to generate an active enzyme. However, binding of folate as the second substrate triggers a large conformational change that activates FPGS and allows the enzyme to adopt a form that is then able to bind the third substrate, L-glutamate, and effect the addition of a polyglutamate tail to the folate.

Direct identification of human oxytocin receptor-binding domains using a photoactivatable cyclic peptide antagonist: comparison with the human V1a vasopressin receptor

Understanding of the molecular determinants responsible for antagonist binding to the oxytocin receptor should provide important insights that facilitate rational design of potential therapeutic agents for the treatment of preterm labor. To study ligand/receptor interactions, we used a novel photosensitive radioiodinated antagonist of the human oxytocin receptor, d(CH(2))(5) [Tyr(Me)(2),Thr(4),Orn(8),Phe(3(125)I,4N(3))-NH(2)9]vasotocin. This ligand had an equivalent high affinity for human oxytocin and V(1a) vasopressin receptors expressed in Chinese hamster ovary cells. Taking advantage of this dual specificity, we conducted photoaffinity labeling experiments on both receptors. Photolabeled oxytocin and V(1a) receptors appeared as a unique protein band at 70-75 kDa and two labeled protein bands at 85-90 and 46 kDa, respectively. To identify contact sites between the antagonist and the receptors, the labeled 70-75- and the 46-kDa proteins were cleaved with CNBr and digested with Lys-C and Arg-C endoproteinases. The fragmentation patterns allowed the identification of a covalently labeled region in the oxytocin receptor transmembrane domain III consisting of the residues Leu(114)-Val(115)-Lys(116). Analysis of contact sites in the V(1a) receptor led to the identification of the homologous region consisting of the residues Val(126)-Val(127)-Lys(128). Binding domains were confirmed by mutation of several CNBr cleavage sites in the oxytocin receptor and of one Lys-C cleavage site in the V(1a) receptor. The results are in agreement with previous experimental data and three-dimensional models of agonist and antagonist binding to members of the oxytocin/vasopressin receptor family.

The interrelationships of side-chain and main-chain conformations in proteins

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.

Structural determinants of cold adaptation and stability in a large protein

The heat-labile alpha-amylase from an antarctic bacterium is the largest known protein that unfolds reversibly according to a two-state transition as shown by differential scanning calorimetry. Mutants of this enzyme were produced, carrying additional weak interactions found in thermostable alpha-amylases. It is shown that single amino acid side chain substitutions can significantly modify the melting point T(m), the calorimetric enthalpy Delta H(cal), the cooperativity and reversibility of unfolding, the thermal inactivation rate constant, and the kinetic parameters k(cat) and K(m). The correlation between thermal inactivation and unfolding reversibility displayed by the mutants also shows that stabilizing interactions increase the frequency of side reactions during refolding, leading to intramolecular mismatches or aggregations typical of large proteins. Although all mutations were located far from the active site, their overall trend is to decrease both k(cat) and K(m) by rigidifying the molecule and to protect mutants against thermal inactivation. The effects of these mutations indicate that the cold-adapted alpha-amylase has lost a large number of weak interactions during evolution to reach the required conformational plasticity for catalysis at low temperatures, thereby producing an enzyme close to the lowest stability allowing maintenance of the native conformation.

Synthetic peptides as inactivators of multimeric enzymes: inhibition of Plasmodium falciparum triosephosphate isomerase by interface peptides

Synthetic peptides corresponding to two distinct segments of the subunit interface of the homodimeric enzyme triosephosphate isomerase (residues 9-18, ANWKCNGTLE, peptide I; residues 68-79, KFGNGSYTGEVS, peptide II) from Plasmodium falciparum (PfTIM) have been investigated for their ability to act as inhibitors by interfering with the quaternary structure of the enzyme. An analog of peptide II containing cysteine at the site corresponding to position 74 and tyrosine at position 69 in the protein sequence KYGNGSCTGEVS (peptide III) was also investigated. A substantial fall in enzyme activity was observed following incubation of the enzyme with peptide II, whereas peptide I did not show any appreciable inhibition. The inhibitory effect was more pronounced on two mutants of PfTIM (Y74C and Y74G), both of which have reduced stability compared to the wild-type protein due to an interface cavity. The IC50 value determined for peptide II is in the range of 0.6-0.8 microM. This study suggests that interface peptides of oligomeric enzymes can be used to inhibit dimeric enzymes by disrupting their native multimeric states and may provide lead structures for potential inhibitor design.

trans-3-Chloroacrylic acid dehalogenase from Pseudomonas pavonaceae 170 shares structural and mechanistic similarities with 4-oxalocrotonate tautomerase

The genes (caaD1 and caaD2) encoding the trans-3-chloroacrylic acid dehalogenase (CaaD) of the 1,3-dichloropropene-utilizing bacterium Pseudomonas pavonaceae 170 were cloned and heterologously expressed in Escherichia coli and Pseudomonas sp. strain GJ1. CaaD is a protein of 50 kDa that is composed of alpha-subunits of 75 amino acid residues and beta-subunits of 70 residues. It catalyzes the hydrolytic cleavage of the beta-vinylic carbon-chlorine bond in trans-3-chloroacrylic acid with a turnover number of 6.4 s(-1). On the basis of sequence similarity, oligomeric structure, and subunit size, CaaD appears to be related to 4-oxalocrotonate tautomerase (4-OT). This tautomerase consists of six identical subunits of 62 amino acid residues and catalyzes the isomerization of 2-oxo-4-hexene-1,6-dioate, via hydroxymuconate, to yield 2-oxo-3-hexene-1,6-dioate. In view of the oligomeric architecture of 4-OT, a trimer of homodimers, CaaD is postulated to be a hexameric protein that functions as a trimer of alpha beta-dimers. The sequence conservation between CaaD and 4-OT and site-directed mutagenesis experiments suggested that Pro-1 of the beta-subunit and Arg-11 of the alpha-subunit are active-site residues in CaaD. Pro-1 could act as the proton acceptor/donor, and Arg-11 is probably involved in carboxylate binding. Based on these findings, a novel dehalogenation mechanism is proposed for the CaaD-catalyzed reaction which does not involve the formation of a covalent enzyme-substrate intermediate.

Spiral, herringbone, and triple-decker silver(I) complexes of benzopyrene derivatives assembled through eta(2)-coordination

Three novel silver(I) complexes with benzopyrene derivatives were synthesized and characterized in this paper. Treatment of AgClO(4)*H(2)O with 7-methylbenzo[a]pyrene (L(1)) afforded [Ag(2)(L(1))(toluene)(0.5)(ClO(4))(2)](n)() (1) which exhibits a 2-D sheet structure with double-stranded helical motifs. Reaction of AgCF(3)SO(3) with dibenzo[b,def ]chrysene (L(2)) gave rise to an unprecedented cocrystallization structure, ([Ag(2)(L(2))(CF(3)SO(3))(2)][Ag(2)(toluene)(2)(CF(3)SO(3))(2)])(n)() (2), formed by a 2-D neutral lamellar polymer and a 1-D neutral rodlike one. The ligand benzo[e]pyrene (L(3)) coordinated to silver(I) ions generating a closed triple-decker tetranuclear complex [Ag(4)(L(3))(4)(p-xylene)(ClO(4))(4)] (3) which can be regarded as a stacking polymer owing to existing intermolecular pi-pi stack interactions. The structural diversity of the silver(I) coordination polymers with polycyclic aromatic hydrocarbons is not only related to the stacking patterns of free polycyclic aromatic hydrocarbons in the crystalline state, but also the geometric shapes of the molecules for these free ligands. In addition, the coordination of solvents to metal ions plays a crucial role in the formation of the unprecedented coordination polymeric architectures. The ESR spectroscopic results, conductivity, and synthesis properties are also discussed.

Significance of aromatic-backbone amide interactions in protein structure

Weakly polar interactions between aromatic rings of amino acids and hydrogens of backbone amides (Ar-HN) have been shown to support local structures in proteins. Their role in secondary structures, however, has not been elucidated. To investigate the relationship between Ar-HN interaction and the stability of local and secondary structures of polypeptides and to improve the prediction of this interaction based on amino acid sequence, the structures of 560 nonhomologous proteins, from the Protein Data Bank, were searched for Ar-HN interactions between the aromatic ring of each Phe, Tyr, and Trp residue at position i and the backbone amide group of any residue, except Pro, at the positions i, i - 1, i - 2, i - 3, i + 1, i + 2, and i + 3. Ar-HN interactions were identified by calculating the chemical shift of the amide hydrogen caused by the proximal aromatic ring. Ar(i)-HN(i + 1, i + 2 and i + 3) interactions were more common (7.10%, 2.08%, and 0.54%, respectively) than were Ar(i)-HN(i - 1, i - 2, and i - 3) interactions (0.66%, <0.1%, and 0.18%, respectively). The value of the chi(1) torsion angle of the aromatic residue in position i depended on the direction of the Ar-HN interaction. The position of the aromatic ring in Ar(i)-HN(i + 1, i + 2, and i + 3) interactions was mostly trans, in Ar(i)-HN(i - 1, i - 2, and i - 3) interactions mainly gauche(-), and in Ar(i)-HN(i) interactions mostly gauche(+). The analyses of the secondary structures of the protein fragments containing Ar-HN interactions showed that Ar-HN interactions were in all types of secondary structures. Search results suggest that Ar-HN interactions have a stabilizing effect on all types of secondary structures.

Structure and function of a novel purine specific nucleoside hydrolase from Trypanosoma vivax

The purine salvage pathway of parasitic protozoa is currently considered as a target for drug development because these organisms cannot synthesize purines de novo. Insight into the structure and mechanism of the involved enzymes can aid in the development of potent inhibitors, leading to new curative drugs. Nucleoside hydrolases are key enzymes in the purine salvage pathway of Trypanosomatidae, and they are especially attractive because they have no equivalent in mammalian cells. We cloned, expressed and purified a nucleoside hydrolase from Trypanosoma vivax. The substrate activity profile establishes the enzyme to be a member of the inosine-adenosine-guanosine-preferring nucleoside hydrolases (IAG-NH). We solved the crystal structure of the enzyme at 1.6 A resolution using MAD techniques. The complex of the enzyme with the substrate analogue 3-deaza-adenosine is presented. These are the first structures of an IAG-NH reported in the literature. The T. vivax IAG-NH is a homodimer, with each subunit consisting of ten beta-strands, 12 alpha-helices and three small 3(10)-helices. Six of the eight strands of the central beta-sheet form a motif resembling the Rossmann fold. Superposition of the active sites of this IAG-NH and the inosine-uridine-preferring nucleoside hydrolase (IU-NH) of Crithidia fasciculata shows the molecular basis of the different substrate specificity distinguishing these two classes of nucleoside hydrolases. An "aromatic stacking network" in the active site of the IAG-NH, absent from the IU-NH, imposes the purine specificity. Asp10 is the proposed general base in the reaction mechanism, abstracting a proton from a nucleophilic water molecule. Asp40 (replaced by Asn39 in the IU-NH) is positioned appropriately to act as a general acid and to protonate the purine leaving group. The second general acid, needed for full enzymatic activity, is probably part of a flexible loop located in the vicinity of the active site.

Distinct cysteine sulfhydryl environments detected by analysis of Raman S-hh markers of Cys-->Ser mutant proteins

Very little is known about the character or functional relevance of hydrogen-bonded cysteine sulfhydryl (S-H) groups in proteins. The Raman S-H band is a unique and sensitive probe of the local S-H environment. Here, we report the use of Raman spectroscopy combined with site-specific mutagenesis to document the existence of five distinguishable hydrogen-bonded states of buried cysteine sulfhydryl groups in a native protein. The 666 residue subunit of the Salmonella typhimurium bacteriophage P22 tailspike contains eight cysteine residues distributed through the elongated structure. The tailspike cysteine residues display an unusual Raman S-H band complex (2500-2600 cm(-1) interval) indicative of diverse S-H hydrogen-bonding interactions in the native trimeric structure. To resolve specific Cys contributions to the complex Raman band we characterized a set of tailspike proteins with each cysteine replaced by a serine. The mutant proteins, once folded, were structurally and functionally indistinguishable from wild-type tailspikes, except for their Raman S-H signatures. Comparison of the Raman spectra of the mutant and wild-type proteins reveals the following hydrogen-bond classes for cysteine sulfhydryl groups. (i) Cys613 forms the strongest S-H...X bond of the tailspike, stronger than any heretofore observed for a protein. (ii) Cys267, Cys287 and Cys458 form robust S-H...X bonds. (iii) Moderate S-H...X bonding occurs for Cys169 and Cys635. (iv) Cys290 and Cys496 form weak hydrogen bonds. (v) It is remarkable that Cys287 contributes two Raman S-H markers, indicating the population of two distinct hydrogen-bonding states. The sum of the S-H Raman signatures of all eight mutants accurately reproduces the composite Raman band of the wild-type tailspike. The diverse cysteine states may be an outcome of the folding and assembly pathway of the tailspike, which though lacking disulfide bonds in the native state, utilizes transient disulfide bonds in the maturation pathway. This Raman study represents the first detailed assessment of local S-H hydrogen bonding in a native protein and provides information not obtainable directly by other structural probes. The method employed here should be applicable to a wide range of cysteine-containing proteins.

C-H...pi-interactions in proteins

A non-redundant set of 1154 protein structures from the Protein Data Bank was examined with respect to close interactions between C-H-donor and pi-acceptor groups. A total of 31,087 interactions were found to satisfy our selection criteria. Their geometric parameters suggest that these interactions can be classified as weak hydrogen bonds.A set of 12 interaction classes were defined based on the division of the donors into three groups and the acceptors into four groups. These classes were examined separately, and the respective interactions described in detail in each class. Most prominent were interactions between aliphatic C-H donors and aromatic pi-acceptors and interactions between aromatic C-H donors and aromatic pi-acceptors. About three-quarters of the Trp-rings, half of all Phe and Tyr-rings and a quarter of all His-rings were found to be involved as acceptors in C-H...pi-interactions. On the donor side, a preference for aromatic C-H groups was observed, but also for the aliphatic side-chains of the long, extended amino acid residues Lys, Arg and Met, and the Pro ring. The average distance between the C-donor and the center-of-mass of the pi-acceptor was observed to be significantly longer in the 174 protein structures determined at >2.5 A resolution. Also, the distribution is significantly wider. This resolution dependence suggests that the force fields commonly used for the refinement of protein structures may not be adequate. C-H...pi-interactions involving aromatic groups either as donor or as acceptor groups are found mostly in the interior of the protein. The more hydrophilic the participating groups are, the closer to the surface are the interactions located. About 40 % of all C-H...pi-interactions occur between amino acid residue side-chains that are separated by nine or less residues in sequence. Dependent on the interaction class, different preferences for secondary structure, residue type and side-chain conformation were observed. It is likely that the C-H...pi-interactions contribute significantly to the overall stability of a protein.

Understanding the role of internal lysine residues of serum albumins in conformational stability and bilirubin binding

The role of internal lysine residues of different serum albumins, viz. from human, rabbit, goat, sheep and buffalo (HSA, RbSA, GSA, SSA and BuSA), in conformational stability and bilirubin binding was investigated after blocking them using acetylation, succinylation and guanidination reactions. No significant change in the secondary structure was noticed whereas the tertiary structure of these proteins was slightly altered upon acetylation or succinylation as revealed by circular dichroism (CD), fluorescence and gel filtration results. Guanidination did not affect the native protein conformation to a measurable extent. Scatchard analysis, CD and absorption spectroscopic results showed marked reductions (5-21-fold decrease in K(a) and approximately 50% decrease in the CD Cotton effect intensity) in the affinity of albumins for bilirubin upon acetylation or succinylation whereas guanidination produced a small change. Interestingly, monosignate CD spectra of bilirubin complexed with GSA, SSA and BuSA were transformed to bisignate CD spectra upon acetylation or succinylation of internal lysine residues whereas spectra remained bisignate in the case of bilirubin bound to acetylated or succinylated derivatives of HSA and RbSA. When probed by CD spectroscopy, bilirubin bound to acetylated or succinylated derivatives of GSA and SSA rapidly switched over to native albumins and not vice versa. These results suggested that salt linkage(s) contributed by internal lysine residue(s) play an important role in the high-affinity binding of bilirubin to albumin and provide stability to the native three-dimensional conformation of the bound pigment. Chloroform severely decreased the intensity of both positive and negative CD Cotton effects of bilirubin complexed with acetylated or succinylated derivatives of all albumins which otherwise increased significantly in the case of bilirubin complexed with native and guanidinated albumin derivatives, except the bilirubin-RbSA complex which showed a small decrease in intensity. These results suggest that the presence of salt linkage(s) in bilirubin-albumin complexation is(are) crucial to bring about effective and efficient stereochemical changes in the bound pigment by co-binding of chloroform which seems to have at least one conserved binding site on these albumins that is shared with bilirubin.

Localization of the C-terminus of rat glutathione S-transferase A1-1: crystal structure of mutants W21F and W21F/F220Y

Twelve C-terminal residues of human glutathione S-transferase A1-1 form a helix in the presence of glutathione-conjugate, or substrate alone, and partly cover the active site. According to X-ray structures, the helix is disordered in the absence of glutathione, but it is not known if it is helical and delocalized, or in a random-coil conformation. Mutation to a tyrosine of residue 220 within this helix was previously shown to affect the pK(a) of Tyr-9 at the active site, in the apo form of the enzyme, and it was proposed that an on-face hydrogen bond between Tyr-220 and Tyr-9 provided a means for affecting this pK(a). In the current study, X-ray structures of the W21F and of the C-terminal mutation, W21F/F220Y, with glutathione sulfonate bound, show that the C-terminal helix is disordered (or delocalized) in the W21F crystal but is visible and ordered in a novel location, a crystal packing crevice, in one of three monomers in the W21F/F220Y crystal, and the proposed hydrogen bond is not formed. Fluorescence spectroscopy studies using an engineered F222W mutant show that the C-terminus remains delocalized in the absence of glutathione or when only the glutathione binding site is occupied, but is ordered and localized in the presence of substrate or conjugate, consistent with these and previous crystallographic studies. Proteins 2001;42:192-200.

Hydrogen bonds with pi-acceptors in proteins: frequencies and role in stabilizing local 3D structures

A comprehensive structural analysis of X--H...pi hydrogen bonding in proteins is performed based on 592 published high-resolution crystal structures (< or = 1.6 A). All potential donors and acceptors are considered, including acidic C--H groups. The sample contains 1311 putative X--H...pi hydrogen bonds with N--H, O--H or S--H donors, that is about one per 10.8 aromatic residues. By far the most efficient pi-acceptor is the side-chain of Trp, which accepts one X--H...pi hydrogen bond per 5.7 residues. The focus of the analysis is on recurrent structural patterns involving regular secondary structure elements. Numerous examples are found where peptide X--H...pi interactions are functional in stabilization of helix termini, strand ends, strand edges, beta-bulges and regular turns. Side-chain X--H...pi hydrogen bonds are formed in considerable numbers in alpha-helices and beta-sheets. Geometrical data on various types of X--H...pi hydrogen bonds are given.

The X-ray structure of the FMN-binding protein AtHal3 provides the structural basis for the activity of a regulatory subunit involved in signal transduction

BACKGROUND

The Arabidopsis thaliana HAL3 gene product encodes for an FMN-binding protein (AtHal3) that is related to plant growth and salt and osmotic tolerance. AtHal3 shows sequence homology to ScHal3, a regulatory subunit of the Saccharomyces cerevisae serine/threonine phosphatase PPz1. It has been proposed that AtHal3 and ScHal3 have similar roles in cellular physiology, as Arabidopsis transgenic plants that overexpress AtHal3 and yeast cells that overexpress ScHal3 display similar phenotypes of improved salt tolerance. The enzymatic activity of AtHal3 has not been investigated. However, the AtHal3 sequence is homologous to that of EpiD, a flavoprotein from Staphylococcus epidermidis that recognizes a peptidic substrate and subsequently catalyzes the alpha, beta-dehydrogenation of its C-terminal cysteine residue.

RESULTS

The X-ray structure of AtHal3 at 2 A resolution reveals that the biological unit is a trimer. Each protomer adopts an alpha/beta Rossmann fold consisting of a six-stranded parallel beta sheet flanked by two layers of alpha helices. The FMN-binding site of AtHal3 contains all the structural requirements of the flavoenzymes that catalyze dehydrogenation reactions. Comparison of the amino acid sequences of AtHal3, ScHal3 and EpiD reveals that a significant number of residues involved in trimer formation, the active site, and FMN binding are conserved. This observation suggests that ScHal3 and EpiD might also be trimers, having a similar structure and function to AtHal3.

CONCLUSIONS

Structural comparisons of AtHal3 with other FMN-binding proteins show that AtHal3 defines a new subgroup of this protein family that is involved in signal transduction. Analysis of the structure of AtHal3 indicates that this protein is designed to interact with another cellular component and to subsequently catalyze the alpha,beta-dehydrogenation of a peptidyl cysteine. Structural data from AtHal3, together with physiological and biochemical information from ScHal3 and EpiD, allow us to propose a model for the recognition and regulation of AtHal3/ScHal3 cellular partners.

Structural basis for the feedback regulation of Escherichia coli pantothenate kinase by coenzyme A

Pantothenate kinase (PanK) is a key regulatory enzyme in the coenzyme A (CoA) biosynthetic pathway and catalyzes the phosphorylation of pantothenic acid to form phosphopantothenate. CoA is a feedback inhibitor of PanK activity by competitive binding to the ATP site. The structures of the Escherichia coli enzyme, in complex with a nonhydrolyzable analogue of ATP, 5'-adenylimido-diphosphate (AMPPNP), or with CoA, were determined at 2.6 and 2.5 A, respectively. Both structures show that two dimers occupy an asymmetric unit; each subunit has a alpha/beta mononucleotide-binding fold with an extensive antiparallel coiled coil formed by two long helices along the dimerization interface. The two ligands, AMPPNP and CoA, associate with PanK in very different ways, but their phosphate binding sites overlap, explaining the kinetic competition between CoA and ATP. Residues Asp(127), His(177), and Arg(243) are proposed to be involved in catalysis, based on modeling of the pentacoordinate transition state. The more potent inhibition by CoA, compared with the CoA thioesters, is explained by a tight interaction of the CoA thiol group with the side chains of aromatic residues, which is predicted to discriminate against the CoA thioesters. The PanK structure provides the framework for a more detailed understanding of the mechanism of catalysis and feedback regulation of PanK.

Crystallographic and energetic analysis of binding of selected anions to the yellow variants of green fluorescent protein

The fluorescence emission of yellow fluorescent proteins (YFPs) has been shown to respond rapidly and reversibly to changes in the concentration of some small anions such as halides; this allows for the use of YFPs as genetically encodable Cl(-) sensors that may be targeted to specific organelles in living cells. Fluorescence is suppressed due to protonation of the chromophore upon anion binding, with a stronger level of interaction at low pH values. At pH 6.0, the apparent dissociation constant (K(app)) for Cl(-) is 32 mM for YFP and 22 mM for YFP-H148Q, whereas at pH 7.5, K(app) is 777 mM and 154 mM, respectively. In the cytosol, YFP-H148Q appears most promising as a halide sensor due to its high degree of sensitivity towards I(-) (K(app)=23 mM at pH 7.5). To aid in the design of variants with improved levels of specificity and affinity for Cl(-), we solved apo and I(-)-bound crystal structures of YFP-H148Q to 2.1 A resolution. The halide-binding site is found near van der Waals contact with the chromophore imidazolinone oxygen atom, in a small buried cavity adjacent to Arg96, which provides electrostatic stabilization. The halide ion is hydrogen bonded to the phenol group of T203Y, consistent with a mutational analysis that indicates that T203Y is indispensible for tight binding. A series of conformational changes occurs in the amphiphilic site upon anion binding, which appear to be propagated to the beta-bulge region around residue 148 on the protein surface. Anion binding raises the chromophore pK(a) values, since delocalization of the phenolate negative charge over the chromophore skeleton is suppressed. Extraction of microscopic binding constants for the linked equilibrium between anion and proton binding indicates that anion selectivity by YFP is related to hydration forces. Specific suggestions to improve Cl(-) binding to YFP-H148Q based on size and hydration energy are proposed.

Modeling of the TCR-MHC-peptide complex

The methodology for generating a homology model of the T1 TCR-PbCS-K(d) class I major histocompatibility complex (MHC) class I complex is presented. The resulting model provides a qualitative explanation of the effect of over 50 different mutations in the region of the complementarity determining region (CDR) loops of the T cell receptor (TCR), the peptide and the MHC's alpha(1)/alpha(2) helices. The peptide is modified by an azido benzoic acid photoreactive group, which is part of the epitope recognized by the TCR. The construction of the model makes use of closely related homologs (the A6 TCR-Tax-HLA A2 complex, the 2C TCR, the 14.3.d TCR Vbeta chain, the 1934.4 TCR Valpha chain, and the H-2 K(b)-ovalbumine peptide), ab initio sampling of CDR loops conformations and experimental data to select from the set of possibilities. The model shows a complex arrangement of the CDR3alpha, CDR1beta, CDR2beta and CDR3beta loops that leads to the highly specific recognition of the photoreactive group. The protocol can be applied systematically to a series of related sequences, permitting the analysis at the structural level of the large TCR repertoire specific for a given peptide-MHC complex.

Molecular and structural basis of the specificity of a neutralizing acetylcholine receptor-mimicking antibody, using combined mutational and molecular modeling analyses

The antagonist activity of short-chain toxins from snake venoms toward the nicotinic acetylcholine receptor (nAChR) is neutralized upon binding to a toxin-specific monoclonal antibody called Malpha2-3 (1). To establish the molecular basis of this specificity, we predicted from both mutational analyses and docking procedures the structure of the Malpha2-3-toxin complex. From knowledge of the functional paratope and epitope, and using a double-mutation cycle procedure, we gathered evidence that Asp(31) in complementarity determining region 1H is close to, and perhaps interacts with, Arg(33) in the antigen. The use of this pair of proximate residues during the selection procedure yielded three models based on docking calculations. The selected models predicted the proximity of Tyr(49) and/or Tyr(50) in the antibody to Lys(47) in the toxin. This was experimentally confirmed using another round of double-mutation cycles. The two models finally selected were submitted to energy minimization in a CHARMM22 force field, and were characterized by a root mean square deviation of 7.0 +/- 2.9 A. Both models display most features of antibody-antigen structures. Since Malpha2-3 also partially mimics some binding properties of nAChR, these structural features not only explain its fine specificity of recognition, but may also further clarify how toxins bind to nAChR.

A new allosteric site in glycogen phosphorylase b as a target for drug interactions

BACKGROUND

In muscle and liver, glycogen concentrations are regulated by the coordinated activities of glycogen phosphorylase (GP) and glycogen synthase. GP exists in two forms: the dephosphorylated low-activity form GPb and the phosphorylated high-activity form GPa. In both forms, allosteric effectors can promote equilibrium between a less active T state and a more active R state. GP is a possible target for drugs that aim to prevent unwanted glycogen breakdown and to stimulate glycogen synthesis in non-insulin-dependent diabetes. As a result of a data bank search, 5-chloro-1H-indole-2-carboxylic acid (1-(4-fluorobenzyl)-2-(4-hydroxypiperidin-1-yl)-2-oxoethy l)amide, CP320626, was identified as a potent inhibitor of human liver GP. Structural studies have been carried out in order to establish the mechanism of this unusual inhibitor.

RESULTS

The structure of the cocrystallised GPb-CP320626 complex has been determined to 2.3 A resolution. CP320626 binds at a site located at the subunit interface in the region of the central cavity of the dimeric structure. The site has not previously been observed to bind ligands and is some 15 A from the AMP allosteric site and 33 A from the catalytic site. The contacts between GPb and CP320626 comprise six hydrogen bonds and extensive van der Waals interactions that create a tight binding site in the T-state conformation of GPb. In the R-state conformation of GPa these interactions are significantly diminished.

CONCLUSIONS

CP320626 inhibits GPb by binding at a new allosteric site. Although over 30 A from the catalytic site, the inhibitor exerts its effects by stabilising the T state at the expense of the R state and thereby shifting the allosteric equilibrium between the two states. The new allosteric binding site offers a further recognition site in the search for improved GP inhibitors.

Novel aromatic residues in transmembrane domains IV and V involved in agonist binding at alpha(1a)-adrenergic receptors

We examined the role that aromatic residues located in the transmembrane helices of the alpha(1a)-adrenergic receptor play in promoting antagonist binding. Since alpha(1)-antagonists display low affinity binding at beta(2)-adrenergic receptors, two phenylalanine residues, Phe-163 and Phe-187, of the alpha(1a)-AR were mutated to the corresponding beta(2)-residue. Neither F163Q nor F187A mutations of the alpha(1a) had any effect on the affinity of the alpha(1)-antagonists. However, the affinity of the endogenous agonist epinephrine was reduced 12.5- and 8-fold by the F163Q and F187A mutations, respectively. An additive loss in affinity (150-fold) for epinephrine was observed at an alpha(1a) containing both mutations. The loss of agonist affinity scenario could be reversed by a gain of affinity with mutation of the corresponding residues in the beta(2) to the phenylalanine residues in the alpha(1a). We propose that both Phe-163 and Phe-187 are involved in independent aromatic interactions with the catechol ring of agonists. The potency but not the efficacy of epinephrine in stimulating phosphatidylinositol hydrolysis was reduced 35-fold at the F163Q/F187A alpha(1a) relative to the wild type receptor. Therefore, Phe-163 and Phe-187 represent novel binding contacts in the agonist binding pocket of the alpha(1a)-AR, but are not involved directly in receptor activation.

Structure-function analysis of alpha-helix H4 using PSE-4 as a model enzyme representative of class A beta-lactamases

We extracted maximum information for structure-function analysis of the PSE-4 class A beta-lactamase by random replacement mutagenesis of three contiguous codons in the H4 alpha-helix at amino acid positions Ala125, Thr126, Met127, Thr128 and Thr129. These positions were predicted to interact with suicide mechanism-based inhibitors when examining the PSE-4 three-dimensional model. Structure-function studies on positions 125-129 indicated that in PSE-4 these amino acids have a role distinct from those in TEM-1, in tolerating substitutions at Ala125 and being invariant at Met127. The importance of Met127 was suspected to be implicated in a structural role in maintaining the integrity of the H4 alpha-helix structure together, thus maintaining the important Ser130-Asp131-Asn132 motif positioned towards the active site. At the structural level, the H4 region was analyzed using energy minimization of the H4 regions of the PSE-4 YAM mutant and compared with wild-type PSE-4. The Tyr 125 of the mutant YAM formed an edge to face pi-pi interaction with Phe 124 which also interacts with the Trp 210 with the same interactions. Antibiotic susceptibilities showed that amino acid changes in the the H4 alpha-helix region of PSE-4 are particularly sensitive to mechanism based-inhibitors. However, kinetic analysis of PSE-4 showed that the two suicide inhibitors belonging to the penicillanic acid sulfone class, sulbactam and tazobactam, were less affected by changes in the H4 alpha-helix region than clavulanic acid, an inhibitor of the oxypenam class. The analysis of H4 alpha-helix in PSE-4 suggests its importance in interactions with the three clinically useful inhibitors and in general to all class A enzymes.

Structural characterization of protein-denaturant interactions: crystal structures of hen egg-white lysozyme in complex with DMSO and guanidinium chloride

A variety of physico-chemical methods employ chemical denaturants to unfold proteins, and study different biophysical processes involved therein. Chemical denaturants are believed to induce unfolding by stabilizing the unfolded state of proteins over the folded state, either macroscopically or through specific interactions. In order to characterize the nature of specific interactions between proteins and denaturants, we have solved crystal structures of hen egg-white lysozyme complexed with denaturants, and report here dimethyl sulfoxide and guanidinium chloride complexes. The dimethyl sulfoxide molecules and guanidinium ions were seen to bind the protein at specific sites and were involved in characteristic interactions. They share a major binding site between them, the C site in the sugar binding cleft of the enzyme. Although the overall conformations of the complexes were very similar to the native structure, spectacular conformational changes were seen to occur locally. Temperature factors were also seen to drop dramatically in the local regions close to the denaturant binding sites. An interesting observation of the present study was the generation of a sodium ion binding site in hen egg-white lysozyme in the presence of denaturants, which was hitherto unknown in any of the other lysozyme structures solved so far. Loss of some of the crucial side chain-main chain interactions may form the initial events in lysozyme unfolding.

Synthesis of thiamacrocycles and conformational studies on their precursors

1-(4-Nitrophenyl)-7-phenylthioheptane (1) and -9-phenylthiononane (2) have been synthesized and their conformations studied in solution and in the solid state. MMX calculations suggest that the global energy minimum structures are bent in the gas phase, probably owing to edge-to-face intramolecular attractive interaction between the electron rich and the electron poor terminal aryl groups. These conformations were confirmed in solution using 2D NOESY NMR. In the solid state, 1 and 2 exist in the staggered, linear conformation, stacked head-to-tail, with the plane of the nitro group being tilted above the plane of the benzene ring. It appears that the crystal lattice forces overcome the weak edge-to-face intramolecular aromatic interactions that dominate in the gas phase and in solution. The corresponding azides were treated with trifluoromethanesulfonic acid to generate the nitrenium ions, which underwent intramolecular ring-closure to give the corresponding 17- and 19-membered ring thiamacrocycles in modest yields. These results support the suggestion that MMX calculations on appropriate model compounds may be useful in predicting which precursors will lead to macrocycles and which will not.

Aromatic Ring Currents at a Protein Surface: Use of (1)H-NMR Chemical Shifts to Refine the Structure of a Naked β Sheet

Abstract The naked β sheet, a newly recognized motif of protein structure, exhibits ordered surfaces in the absence of a conventional hydrophobic core. A model is provided by an archaeal Zn ribbon homologous to eukaryotic RNA polymerase II subunit 9 (RPB9). This subunit, which regulates transcriptional start-site selection and downstream pausing, contains Zn(2+)-binding motifs similar to those of general transcription factors TFIIB and TFIIS. Interestingly, distance-geometry yields two models of the archaeal Zn ribbon differing in the orientation of a conserved tyrosine side chain on the well-ordered surface of the naked β-sheet. The models are equally consistent with conventional restraints and otherwise contain indistinguishable structural features, including a tetrahedral Cys(4) Zn(2+)-binding sites, four antiparallel β-strands, and disordered loop. Due to the change in tyrosine orientation and correlated changes in the configuration of neighboring side chains, the two models predict inequivalent patterns of aromatic ring-current shifts. The observed secondary shifts of adjoining resonances are shown to be consistent with one model but not the other. In the consistent model the surface of the β-sheet contains successive aromatic edge-to-face contacts in accord with semi-classical and ab initio potentials. We speculate that the aromatic-rich surface of the hyperthermophilic RPB9 domain contributes its thermodynamic stability and provides a nucleic-acid-binding site in the eukaryotic and archaeal transcriptional machinery. The present study demonstrates how the reduced dimensionality of a surface can lead to ambiguities in the interpretation of nuclear Overhauser enhancements. The results illustrate the utility of chemical shifts at such a surface in the cross-validation of a high-resolution solution structure.

Crystal structures of intermediates in the dehalogenation of haloalkanoates by L-2-haloacid dehalogenase

The L-2-haloacid dehalogenase from the 1,2-dichloroethane-degrading bacterium Xanthobacter autotrophicus GJ10 catalyzes the hydrolytic dehalogenation of small L-2-haloalkanoates to their corresponding D-2-hydroxyalkanoates, with inversion of the configuration at the C(2) atom. The structure of the apoenzyme at pH 8 was refined at 1.5-A resolution. By lowering the pH, the catalytic activity of the enzyme was considerably reduced, allowing the crystal structure determination of the complexes with L-2-monochloropropionate and monochloroacetate at 1.7 and 2.1 A resolution, respectively. Both complexes showed unambiguous electron density extending from the nucleophile Asp(8) to the C(2) atom of the dechlorinated substrates corresponding to a covalent enzyme-ester reaction intermediate. The halide ion that is cleaved off is found in line with the Asp(8) Odelta1-C(2) bond in a halide-stabilizing cradle made up of Arg(39), Asn(115), and Phe(175). In both complexes, the Asp(8) Odelta2 carbonyl oxygen atom interacts with Thr(12), Ser(171), and Asn(173), which possibly constitute the oxyanion hole in the hydrolysis of the ester bond. The carboxyl moiety of the substrate is held in position by interactions with Ser(114), Lys(147), and main chain NH groups. The L-2-monochloropropionate CH(3) group is located in a small pocket formed by side chain atoms of Lys(147), Asn(173), Phe(175), and Asp(176). The size and position of the pocket explain the stereospecificity and the limited substrate specificity of the enzyme. These crystallographic results demonstrate that the reaction of the enzyme proceeds via the formation of a covalent enzyme-ester intermediate at the nucleophile Asp(8).

Structure-assisted design of mechanism-based irreversible inhibitors of human rhinovirus 3C protease with potent antiviral activity against multiple rhinovirus serotypes

Human rhinoviruses, the most important etiologic agents of the common cold, are messenger-active single-stranded monocistronic RNA viruses that have evolved a highly complex cascade of proteolytic processing events to control viral gene expression and replication. Most maturation cleavages within the precursor polyprotein are mediated by rhinovirus 3C protease (or its immediate precursor, 3CD), a cysteine protease with a trypsin-like polypeptide fold. High-resolution crystal structures of the enzyme from three viral serotypes have been used for the design and elaboration of 3C protease inhibitors representing different structural and chemical classes. Inhibitors having alpha,beta-unsaturated carbonyl groups combined with peptidyl-binding elements specific for 3C protease undergo a Michael reaction mediated by nucleophilic addition of the enzyme's catalytic Cys-147, resulting in covalent-bond formation and irreversible inactivation of the viral protease. Direct inhibition of 3C proteolytic activity in virally infected cells treated with these compounds can be inferred from dose-dependent accumulations of viral precursor polyproteins as determined by SDS/PAGE analysis of radiolabeled proteins. Cocrystal-structure-assisted optimization of 3C-protease-directed Michael acceptors has yielded molecules having extremely rapid in vitro inactivation of the viral protease, potent antiviral activity against multiple rhinovirus serotypes and low cellular toxicity. Recently, one compound in this series, AG7088, has entered clinical trials.

Understanding the structure and function of catalases: clues from molecular evolution and in vitro mutagenesis

This review gives an overview about the structural organisation of different evolutionary lines of all enzymes capable of efficient dismutation of hydrogen peroxide. Major potential applications in biotechnology and clinical medicine justify further investigations. According to structural and functional similarities catalases can be divided in three subgroups. Typical catalases are homotetrameric haem proteins. The three-dimensional structure of six representatives has been resolved to atomic resolution. The central core of each subunit reveals a characteristic "catalase fold", extremely well conserved among this group. In the native tetramer structure pairs of subunits tightly interact via exchange of their N-terminal arms. This pseudo-knot structures implies a highly ordered assembly pathway. A minor subgroup ("large catalases") possesses an extra flavodoxin-like C-terminal domain. A > or = 25 A long channel leads from the enzyme surface to the deeply buried active site. It enables rapid and selective diffusion of the substrates to the active center. In several catalases NADPH is tightly bound close to the surface. This cofactor may prevent and reverse the formation of compound II, an inactive reaction intermediate. Bifunctional catalase-peroxidase are haem proteins which probably arose via gene duplication of an ancestral peroxidase gene. No detailed structural information is currently available. Even less is know about manganese catalases. Their di-manganese reaction centers may be evolutionary.

Pairwise interactions between neuronal alpha7 acetylcholine receptors and alpha-conotoxin ImI

The present work uses alpha-conotoxin ImI (CTx ImI) to probe the neurotransmitter binding site of neuronal alpha7 acetylcholine receptors. We identify key residues in alpha7 that contribute to CTx ImI affinity, and use mutant cycles analysis to identify pairs of residues that stabilize the receptor-conotoxin complex. We first mutated key residues in the seven known loops of alpha7 that converge at the subunit interface to form the ligand binding site. The mutant subunits were expressed in 293 HEK cells, and CTx ImI binding was measured by competition against the initial rate of 125I-alpha-bungarotoxin binding. The results reveal a predominant contribution by Tyr-195 in alpha7, accompanied by smaller contributions by Thr-77, Tyr-93, Asn-111, Gln-117, and Trp-149. Based upon our previous identification of bioactive residues in CTx ImI, we measured binding of receptor and toxin mutations and analyzed the results using thermodynamic mutant cycles. The results reveal a single dominant interaction between Arg-7 of CTx ImI and Tyr-195 of alpha7 that anchors the toxin to the binding site. We also find multiple weak interactions between Asp-5 of CTx ImI and Trp-149, Tyr-151, and Gly-153 of alpha7, and between Trp-10 of CTx ImI and Thr-77 and Asn-111 of alpha7. The overall results establish the orientation of CTx ImI as it bridges the subunit interface and demonstrate close approach of residues on opposing faces of the alpha7 binding site.

Relating structure to thermodynamics: the crystal structures and binding affinity of eight OppA-peptide complexes

The oligopeptide-binding protein OppA provides a useful model system for studying the physical chemistry underlying noncovalent interactions since it binds a variety of readily synthesized ligands. We have studied the binding of eight closely related tripeptides of the type Lysine-X-Lysine, where X is an abnormal amino acid, by isothermal titration calorimetry (ITC) and X-ray crystallography. The tripeptides fall into three series of ligands, which have been designed to examine the effects of small changes to the central side chain. Three ligands have a primary amine as the second side chain, two have a straight alkane chain, and three have ring systems. The results have revealed a definite preference for the binding of hydrophobic residues over the positively charged side chains, the latter binding only weakly due to unfavorable enthalpic effects. Within the series of positively charged groups, a point of lowest affinity has been identified and this is proposed to arise from unfavorable electrostatic interactions in the pocket, including the disruption of a key salt bridge. Marked entropy-enthalpy compensation is found across the series, and some of the difficulties in designing tightly binding ligands have been highlighted.