A Structural Portrait of the PDZ Domain Family

PDZ (PSD-95/Discs-large/ZO1) domains are interaction modules that typically bind to specific C-terminal sequences of partner proteins and assemble signaling complexes in multicellular organisms. We have analyzed the existing database of PDZ domain structures in the context of a specificity tree based on binding specificities defined by peptide-phage binding selections. We have identified 16 structures of PDZ domains in complex with high-affinity ligands and have elucidated four additional structures to assemble a structural database that covers most of the branches of the PDZ specificity tree. A detailed comparison of the structures reveals features that are responsible for the diverse specificities across the PDZ domain family. Specificity differences can be explained by differences in PDZ residues that are in contact with the peptide ligands, but these contacts involve both side-chain and main-chain interactions. Most PDZ domains bind peptides in a canonical conformation in which the ligand main chain adopts an extended β-strand conformation by interacting in an antiparallel fashion with a PDZ β-strand. However, a subset of PDZ domains bind peptides with a bent main-chain conformation and the specificities of these non-canonical domains could not be explained based on canonical structures. Our analysis provides a structural portrait of the PDZ domain family, which serves as a guide in understanding the structural basis for the diverse specificities across the family.     

Structural Basis of Conserved Cysteine in the Fibroblast Growth Factor Family: Evidence for a Vestigial Half-Cystine

The 22 members of the mouse/human fibroblast growth factor (FGF) family of proteins contain a conserved cysteine residue at position 83 (numbering scheme of the 140-residue form of FGF-1). Sequence and structure information suggests that this position is a free cysteine in 16 members and participates as a half-cystine in at least 3 (and perhaps as many as 6) other members. While a structural role as a half-cystine provides a stability basis for possible selective pressure, it is less clear why this residue is conserved as a free cysteine (although free buried thiols can limit protein functional half-life). To probe the structural role of the free cysteine at position 83 in FGF-1, we constructed Ala, Ser, Thr, Val, and Ile mutations and determined their effects on structure and stability. These results show that position 83 in FGF-1 is thermodynamically optimized to accept a free cysteine. A second cysteine mutation was introduced into wild-type FGF-1 at adjacent position Ala66, which is known to participate as a half-cystine with position 83 in FGF-8, FGF-19, and FGF-23. Results show that, unlike position 83, a free cysteine at position 66 destabilizes FGF-1; however, upon oxidation, a near-optimal disulfide bond is formed between Cys66 and Cys83, resulting in ∼ 14 kJ/mol of increased thermostability. Thus, while the conserved free cysteine at position 83 in the majority of the FGF proteins may have a principal role in limiting functional half-life, evidence suggests that it is a vestigial half-cystine.     

Systematic Mutational Analysis of Peptide Inhibition of the p53-MDM2/MDMX Interactions

Inhibition of the interaction between the tumor suppressor protein p53 and its negative regulators MDM2 and MDMX is of great interest in cancer biology and drug design. We previously reported a potent duodecimal peptide inhibitor, termed PMI (TSFAEYWNLLSP), of the p53-MDM2 and -MDMX interactions. PMI competes with p53 for MDM2 and MDMX binding at an affinity roughly 2 orders of magnitude higher than that of ^17-28p53 (ETFSDLWKLLPE) of the same length; both peptides adopt nearly identical α-helical conformations in the complexes, where the three highlighted hydrophobic residues Phe, Trp, and Leu dominate PMI or ^17-28p53 binding to MDM2 and MDMX. To elucidate the molecular determinants for PMI activity and specificity, we performed a systematic Ala scanning mutational analysis of PMI and ^17-28p53. The binding affinities for MDM2 and MDMX of a total of 35 peptides including 10 truncation analogs were quantified, affording a complete dissection of energetic contributions of individual residues of PMI and ^17-28p53 to MDM2 and MDMX association. Importantly, the N8A mutation turned PMI into the most potent dual-specific antagonist of MDM2 and MDMX reported to date, registering respective K_d values of 490 pM and 2.4 nM. The co-crystal structure of N8A-PMI-^25-109MDM2 was determined at 1.95 Å, affirming that high-affinity peptide binding to MDM2/MDMX necessitates, in addition to optimized intermolecular interactions, enhanced helix stability or propensity contributed by non-contact residues. The powerful empirical binding data and crystal structures present a unique opportunity for computational studies of peptide inhibition of the p53-MDM2/MDMX interactions.     

Substitution of Glu122 by Glutamine Revealed the Function of the Second Water Molecule as a Proton Donor in the Binuclear Metal Enzyme Creatininase

Creatininase is a binuclear zinc enzyme and catalyzes the reversible conversion of creatinine to creatine. It exhibits an open-closed conformational change upon substrate binding, and the differences in the conformations of Tyr121, Trp154, and the loop region containing Trp174 were evident in the enzyme-creatine complex when compared to those in the ligand-free enzyme. We have determined the crystal structure of the enzyme complexed with a 1-methylguanidine. All subunits in the complex existed as the closed form, and the binding mode of creatinine was estimated. Site-directed mutagenesis revealed that the hydrophobic residues that show conformational change upon substrate binding are important for the enzyme activity.We propose a catalytic mechanism of creatininase in which two water molecules have significant roles. The first molecule is a hydroxide ion (Wat1) that is bound as a bridge between the two metal ions and attacks the carbonyl carbon of the substrate. The second molecule is a water molecule (Wat2) that is bound to the carboxyl group of Glu122 and functions as a proton donor in catalysis. The activity of the E122Q mutant was very low and it was only partially restored by the addition of ZnCl₂ or MnCl₂. In the E122Q mutant, k_cat is drastically decreased, indicating that Glu122 is important for catalysis. X-ray crystallographic study and the atomic absorption spectrometry analysis of the E122Q mutant-substrate complex revealed that the drastic decrease of the activity of the E122Q was caused by not only the loss of one Zn ion at the Metal1 site but also a critical function of Glu122, which most likely exists for a proton transfer step through Wat2.     

Structural and Thermodynamic Characterization of Pre- and Postpolymerization States in the F4 Fimbrial Subunit FaeG

Enterotoxigenic Escherichia coli expressing F4 fimbriae are the major cause of porcine colibacillosis and are responsible for significant death and morbidity in neonatal and postweaned piglets. Via the chaperone-usher pathway, F4 fimbriae are assembled into thin, flexible polymers mainly composed of the single-domain adhesin FaeG. The F4 fimbrial system has been labeled eccentric because the F4 pilins show some features distinct from the features of pilins of other chaperone-usher-assembled structures. In particular, FaeG is much larger than other pilins (27  versus ∼ 17 kDa), grafting an additional carbohydrate binding domain on the common immunoglobulin-like core. Structural data of FaeG during different stages of the F4 fimbrial biogenesis process, combined with differential scanning calorimetry measurements, confirm the general principles of the donor strand complementation/exchange mechanisms taking place during pilus biogenesis via the chaperone-usher pathway.     

Mapping of protein-protein interaction sites by the 'absence of interference' approach

Protein-protein interactions are critical to most biological processes, and locating protein-protein interfaces on protein structures is an important task in molecular biology. We developed a new experimental strategy called the ‘absence of interference’ approach to determine surface residues involved in protein-protein interaction of established yeast two-hybrid pairs of interacting proteins. One of the proteins is subjected to high-level randomization by error-prone PCR. The resulting library is selected by yeast two-hybrid system for interacting clones that are isolated and sequenced. The interaction region can be identified by an absence or depletion of mutations. For data analysis and presentation, we developed a Web interface that analyzes the mutational spectrum and displays the mutational frequency on the surface of the structure (or a structural model) of the randomized protein†. Additionally, this interface might be of use for the display of mutational distributions determined by other types of random mutagenesis experiments. We applied the approach to map the interface of the catalytic domain of the DNA methyltransferase Dnmt3a with its regulatory factor Dnmt3L. Dnmt3a was randomized with high mutational load. A total of 76 interacting clones were isolated and sequenced, and 648 mutations were identified. The mutational pattern allowed to identify a unique interaction region on the surface of Dnmt3a, which comprises about 500-600 Å². The results were confirmed by site-directed mutagenesis and structural analysis. The absence-of-interference approach will allow high-throughput mapping of protein interaction sites suitable for functional studies and protein docking.     

Using Soft X-Rays for a Detailed Picture of Divalent Metal Binding in the Nucleosome

Divalent metals associate with DNA in a site-selective manner, which can influence nucleosome positioning, mobility, compaction, and recognition by nuclear factors. We previously characterized divalent metal binding in the nucleosome core using hard (short-wavelength) X-rays allowing high-resolution crystallographic determination of the strongest affinity sites, which revealed that Mn²⁺ associates with the DNA major groove in a sequence- and conformation-dependent manner. In this study, we obtained diffraction data with soft X-rays at the Mn²⁺ absorption edge for a core particle crystal in the presence of 10 mM MnSO₄, mimicking prevailing Mg²⁺ concentration in the nucleus. This provides an exceptional view of counterion binding in the nucleosome through identification of 45 divalent metal binding sites.In addition to that at the well-characterized major interparticle interface, only one other histone-divalent metal binding site is found, which corresponds to a symmetry-related counterpart on the ‘free’ H2B α1 helix C-terminus. This emphasizes the importance of the α-helix dipole in ion binding and suggests that the H2B motif may serve as a nucleation site in nucleosome compaction. The 43 sites associated with the DNA are characterized by (1) high-affinity direct coordination at the most electrostatically favorable major groove locations, (2) metal hydrate binding to the major groove, (3) direct coordination to phosphate groups at sites of high charge density, (4) metal hydrate binding in the minor groove, or (5) metal hydrate-divalent anion pairing. Metal hydrates are found within the minor groove only at locations displaying a narrow range of high-intermediate width and to which histone N-terminal tails are not associated or proximal. This indicates that divalent metals and histone tails can both collaborate and compete in minor groove association, which sheds light on nucleosome solubility and chromatin compaction behavior.     

The capsid of the small RNA phage PRR1 is stabilized by metal ions

Many nonenveloped virus particles are stabilized by calcium ions bound in the interfaces between the protein subunits. These ions may have a role in the disassembly process. The small RNA phages of the Leviviridae family have T = 3 quasi-symmetry and are unique among simple viruses in that they have a coat protein with a translational repressor activity and a fold that has not been observed in other viruses. The crystal structure of phage PRR1 has been determined to 3.5 Å resolution. The structure shows a tentative binding site for a calcium ion close to the quasi-3-fold axis. The RNA-binding surface used for repressor activity is mostly conserved. The structure does not show any significant differences between quasi-equivalent subunits, which suggests that the assembly is not controlled by conformational switches as in many other simple viruses.     

On the Structure of the Proton-Binding Site in the Fo Rotor of Chloroplast ATP Synthases

The recently reported crystal structures of the membrane-embedded proton-dependent c-ring rotors of a cyanobacterial F₁F_o ATP synthase and a chloroplast F₁F_o ATP synthase have provided new insights into the mechanism of this essential enzyme. While the overall features of these c-rings are similar, a discrepancy in the structure and hydrogen-bonding interaction network of the H⁺ sites suggests two distinct binding modes, potentially reflecting a mechanistic differentiation. Importantly, the conformation of the key glutamate side chain to which the proton binds is also altered. To investigate the nature of these differences, we use molecular dynamics simulations of both c-rings embedded in a phospholipid membrane. We observe that the structure of the c₁₅ ring from Spirulina platensis is unequivocally stable within the simulation time. By contrast, the proposed structure of the H⁺ site in the chloroplast c₁₄ ring changes rapidly and consistently into that reported for the c₁₅ ring, indicating that the latter represents a common binding mode. To assess this hypothesis, we have remodeled the c₁₄ ring by molecular replacement using the published structure factors. The resulting structure provides clear evidence in support of a common binding site conformation and is also considerably improved statistically. These findings, taken together with a sequence analysis of c-subunits in the ATP synthase family, indicate that the so-called proton-locked conformation observed in the c₁₅ ring may be a common characteristic not only of light-driven systems such as chloroplasts and cyanobacteria but also of a selection of other bacterial species.     

The role of RNA sequence and structure in RNA--protein interactions

We investigate the sequence and structural properties of RNA-protein interaction sites in 211 RNA-protein chain pairs, the largest set of RNA-protein complexes analyzed to date. Statistical analysis confirms and extends earlier analyses made on smaller data sets. There are 24.6% of hydrogen bonds between RNA and protein that are nucleobase specific, indicating the importance of both nucleobase-specific and -nonspecific interactions. While there is no significant difference between RNA base frequencies in protein-binding and non-binding regions, distinct preferences for RNA bases, RNA structural states, protein residues, and protein secondary structure emerge when nucleobase-specific and -nonspecific interactions are considered separately. Guanine nucleobase and unpaired RNA structural states are significantly preferred in nucleobase-specific interactions; however, nonspecific interactions disfavor guanine, while still favoring unpaired RNA structural states. The opposite preferences of nucleobase-specific and -nonspecific interactions for guanine may explain discrepancies between earlier studies with regard to base preferences in RNA-protein interaction regions. Preferences for amino acid residues differ significantly between nucleobase-specific and -nonspecific interactions, with nonspecific interactions showing the expected bias towards positively charged residues. Irregular protein structures are strongly favored in interactions with the protein backbone, whereas there is little preference for specific protein secondary structure in either nucleobase-specific interaction or -nonspecific interaction. Overall, this study shows strong preferences for both RNA bases and RNA structural states in protein-RNA interactions, indicating their mutual importance in protein recognition.     

Tertiary Motifs Revealed in Analyses of Higher-Order RNA Junctions

RNA junctions are secondary-structure elements formed when three or more helices come together. They are present in diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze currently solved 3D RNA junctions in terms of base-pair interactions and 3D configurations. First, we study base-pair interaction diagrams for solved RNA junctions with 5 to 10 helices and discuss common features. Second, we compare these higher-order junctions to those containing 3 or 4 helices and identify global motif patterns such as coaxial stacking and parallel and perpendicular helical configurations. These analyses show that higher-order junctions organize their helical components in parallel and helical configurations similar to lower-order junctions. Their sub-junctions also resemble local helical configurations found in three- and four-way junctions and are stabilized by similar long-range interaction preferences such as A-minor interactions. Furthermore, loop regions within junctions are high in adenine but low in cytosine, and in agreement with previous studies, we suggest that coaxial stacking between helices likely forms when the common single-stranded loop is small in size; however, other factors such as stacking interactions involving noncanonical base pairs and proteins can greatly determine or disrupt coaxial stacking. Finally, we introduce the ribo-base interactions: when combined with the along-groove packing motif, these ribo-base interactions form novel motifs involved in perpendicular helix-helix interactions. Overall, these analyses suggest recurrent tertiary motifs that stabilize junction architecture, pack helices, and help form helical configurations that occur as sub-elements of larger junction networks. The frequent occurrence of similar helical motifs suggest nature's finite and perhaps limited repertoire of RNA helical conformation preferences. More generally, studies of RNA junctions and tertiary building blocks can ultimately help in the difficult task of RNA 3D structure prediction.     

Energetics-based protein profiling on a proteomic scale: identification of proteins resistant to proteolysis

Native states of proteins are flexible, populating more than just the unique native conformation. The energetics and dynamics resulting from this conformational ensemble are inherently linked to protein function and regulation. Proteolytic susceptibility is one feature determined by this conformational energy landscape. As an attempt to investigate energetics of proteins on a proteomic scale, we challenged the Escherichia coli proteome with extensive proteolysis and determined which proteins, if any, have optimized their energy landscape for resistance to proteolysis. To our surprise, multiple soluble proteins survived the challenge. Maltose binding protein, a survivor from thermolysin digestion, was characterized by in vitro biophysical studies to identify the physical origin of proteolytic resistance. This experimental characterization shows that kinetic stability is responsible for the unusual resistance in maltose binding protein. The biochemical functions of the identified survivors suggest that many of these proteins may have evolved extreme proteolytic resistance because of their critical roles under stressed conditions. Our results suggest that under functional selection proteins can evolve extreme proteolysis resistance by modulating their conformational energy landscapes without the need to invent new folds, and that proteins can be profiled on a proteomic scale according to their energetic properties by using proteolysis as a structural probe.     

Determination of the transition state ensemble for the folding of ubiquitin from a combination of Phi and Psi analyses

Protein engineering techniques have emerged as powerful tools for characterizing transition states (TSs) for protein folding. Recently, the Ψ analysis, in which double-histidine mutations create the possibility of reversible crosslinking in the native state, has been proposed as an additional approach to the well-established Φ analysis. We present here a combination of these two procedures for defining the structure of the TS of ubiquitin, a small α/β protein that has been used extensively as a model system for both experimental and computational studies of the protein-folding process. We performed a series of molecular dynamics simulations in which Φ and Ψ values were used as ensemble-averaged structural restraints to determine an ensemble of structures representing the TS of ubiquitin. Although the available Ψ values for ubiquitin did not, by themselves, generate well-defined TS ensembles, the inclusion of the restricted set of zero or unity values, but not fractional ones, provided useful complementary information to the Φ analysis. Our results show that the TS of ubiquitin is formed by a relatively narrow ensemble of structures exhibiting an overall native-like topology in which the N-terminal and C-terminal regions are in close proximity.     

Mechanism for Attenuation of DNA Binding by MarR Family Transcriptional Regulators by Small Molecule Ligands

Members of the multiple antibiotic resistance regulator (MarR) family control gene expression in a variety of metabolic processes in bacteria and archaea. Hypothetical uricase regulator (HucR), which belongs to the ligand-responsive branch of the MarR family, regulates uricase expression in Deinococcus radiodurans by binding a shared promoter region between uricase and HucR genes. We show here that HucR responds only to urate and, to a lesser extent, to xanthine by attenuated DNA binding, compared to other intermediates of purine degradation. Using molecular-dynamics-guided mutational analysis, we identified the ligand-binding site in HucR. Electrophoretic mobility shift assays and intrinsic Trp fluorescence have identified W20 from the N-terminal helix and R80 from helix 3, which serves as a scaffold for the DNA recognition helix, as being essential for ligand binding. Using structural data combined with in silico and in vitro analyses, we propose a mechanism for the attenuation of DNA binding in which a conformational change initiated by charge repulsion due to a bound ligand propagates to DNA recognition helices. This mechanism may apply generally to MarR homologs that bind anionic phenolic ligands.     

Structure of Bacteriophage T4 Endonuclease II Mutant E118A, a Tetrameric GIY-YIG Enzyme

Coliphage T4 endonuclease II (EndoII), encoded by gene denA, is a small (16 kDa, 136 aa) enzyme belonging to the GIY-YIG family of endonucleases, which lacks a C-terminal domain corresponding to that providing most of the binding energy in the structurally characterized GIY-YIG endonucleases, I-TevI and UvrC. In vivo, it is involved in degradation of host DNA, permitting scavenging of host-derived nucleotides for phage DNA synthesis. EndoII primarily catalyzes single-stranded nicking of DNA; 5- to 10-fold less frequently double-stranded breaks are produced. The Glu118Ala mutant of EndoII was crystallized in space group P2₁ with four monomers in the asymmetric unit. The fold of the EndoII monomer is similar to that of the catalytic domains of UvrC and I-TevI. In contrast to these enzymes, EndoII forms a striking X-shaped tetrameric structure composed as a dimer of dimers, with a protruding hairpin domain not present in UvrC or I-TevI providing most of the dimerization and tetramerization interfaces. A bound phosphate ion in one of the four active sites of EndoII likely mimics the scissile phosphate in a true substrate complex. In silico docking experiments showed that a protruding loop containing a nuclease-associated modular domain 3 element is likely to be involved in substrate binding, as well as residues forming a separate nucleic acid binding surface adjacent to the active site. The positioning of these sites within the EndoII primary dimer suggests that the substrate would bind to a primary EndoII dimer diagonally over the active sites, requiring significant distortion of the enzyme or the substrate DNA, or both, for simultaneous nicking of both DNA strands. The scarcity of potential nucleic acid binding residues between the active sites indicates that EndoII may bind its substrate inefficiently across the two sites in the dimer, offering a plausible explanation for the catalytic preponderance of single-strand nicks. Mutations analyzed in earlier functional studies are discussed in their structural context.     

Crystal structure of the C-terminal domain of a flagellar hook-capping protein from Xanthomonas campestris

The crystal structure of the C-terminal domain of a hook-capping protein FlgD from the plant pathogen Xanthomonas campestris (Xc) has been determined to a resolution of ca 2.5 Å using X-ray crystallography. The monomer of whole FlgD comprises 221 amino acids with a molecular mass of 22.7 kDa, but the flexible N-terminus is cleaved for up to 75 residues during crystallization. The final structure of the C-terminal domain reveals a novel hybrid comprising a tudor-like domain interdigitated with a fibronectin type III domain. The C-terminal domain of XcFlgD forms three types of dimers in the crystal. In agreement with this, analytical ultracentrifugation and gel filtration experiments reveal that they form a stable dimer in solution. From these results, we propose that the Xc flagellar hook cap protein FlgD comprises two individual domains, a flexible N-terminal domain that cannot be detected in the current study and a stable C-terminal domain that forms a stable dimer.     

Crystal structures of cobalamin-independent methionine synthase (MetE) from Streptococcus mutans: a dynamic zinc-inversion model

Cobalamin-independent methionine synthase (MetE) catalyzes the direct transfer of a methyl group from methyltetrahydrofolate to l-homocysteine to form methionine. Previous studies have shown that the MetE active site coordinates a zinc atom, which is thought to act as a Lewis acid and plays a role in the activation of thiol. Extended X-ray absorption fine structure studies and mutagenesis experiments identified the zinc-binding site in MetE from Escherichia coli. Further structural investigations of MetE from Thermotoga maritima lead to the proposition of two models: “induced fit” and “dynamic equilibrium”, to account for the catalytic mechanisms of MetE. Here, we present crystal structures of oxidized and zinc-replete MetE from Streptococcus mutans at the physiological pH. The structures reveal that zinc is mobile in the active center and has the possibility to invert even in the absence of homocysteine. These structures provide evidence for the dynamic equilibrium model.