An irregular beta-bulge common to a group of bacterial RNases is an important determinant of stability and function in barnase

β-Bulges: extensive structural analyses of β-sheets irregularities

β-Sheets are quite frequent in protein structures and are stabilized by regular main-chain hydrogen bond patterns. Irregularities in β-sheets, named β-bulges, are distorted regions between two consecutive hydrogen bonds. They disrupt the classical alternation of side chain direction and can alter the directionality of β-strands. They are implicated in protein-protein interactions and are introduced to avoid β-strand aggregation. Five different types of β-bulges are defined. Previous studies on β-bulges were performed on a limited number of protein structures or one specific family. These studies evoked a potential conservation during evolution. In this work, we analyze the β-bulge distribution and conservation in terms of local backbone conformations and amino acid composition. Our dataset consists of 66 times more β-bulges than the last systematic study (Chan et al. Protein Science 1993, 2:1574-1590). Novel amino acid preferences are underlined and local structure conformations are highlighted by the use of a structural alphabet. We observed that β-bulges are preferably localized at the N- and C-termini of β-strands, but contrary to the earlier studies, no significant conservation of β-bulges was observed among structural homologues. Displacement of β-bulges along the sequence was also investigated by Molecular Dynamics simulations.

Identifying the hotspots on the top faces of WD40-repeat proteins from their primary sequences by β-bulges and DHSW tetrads

The analysis of 36 available crystal structures of WD40 repeat proteins reveals widespread existence of a beta-bulge formed at the beginning of strand a and the end of strand b, termed as WD_b–a bulge: among a total of 259 WD40 blades, there are 243 such β-bulges. The R₁ positions in these WD_b–a bulges have fair distributions of Arg, His, Ile, Leu, Lys, Met, Phe, Trp, Tyr and Val residues. These residues protrude on the top face of the WD40 proteins and can serve as hotspots for protein-protein interactions. An analysis of 29 protein complexes formed by 17 WD proteins reveals that these R₁ residues, along with two other residues (R₁-2 and D-1), are indeed widely involved in protein-protein interactions. Interestingly, these WD_b–a bulges can be easily identified by the 4-amino acid sequences of (V, L, I), R₁, R₂, (V, L, I), along with some other significant amino acids. Thus, the hotspots of WD40 proteins on the top face can be readily predicted based on the primary sequences of the proteins. The literature-reported mutagenesis studies for Met30, MDV1, Tup11, COP1 and SPA1, which crystal structures are not available, can be readily understood based on the feature-based method. Applying the method, the twelve potential hotspots on the top face of Tup11 from S. japonicas have been identified. Our ITC measurements confirm seven of them, Tyr382, Arg284, Tyr426, Tyr508, Leu559, Lys575 and Ile601, are essential for recognizing Fep1. The ITC measurements further convinced that the feature-based method provides accurate prediction of hotspots on the top face.

Protein characterization of a candidate mechanism SNP for Crohn's disease: the macrophage stimulating protein R689C substitution

High throughput genome wide associations studies (GWAS) are now identifying a large number of genome loci related to risk of common human disease. Each such locus presents a challenge in identifying the relevant underlying mechanism. Here we report the experimental characterization of a proposed causal single nucleotide polymorphism (SNP) in a locus related to risk of Crohn's disease and ulcerative colitis. The SNP lies in the MST1 gene encoding Macrophage Stimulating Protein (MSP), and results in an R689C amino acid substitution within the β-chain of MSP (MSPβ). MSP binding to the RON receptor tyrosine kinase activates signaling pathways involved in the inflammatory response. We have purified wild-type and mutant MSPβ proteins and compared biochemical and biophysical properties that might impact the MSP/RON signaling pathway. Surface plasmon resonance (SPR) binding studies showed that MSPβ R689C affinity to RON is approximately 10-fold lower than that of the wild-type MSPβ and differential scanning fluorimetry (DSF) showed that the thermal stability of the mutant MSPβ was slightly lower than that of wild-type MSPβ, by 1.6 K. The substitution was found not to impair the specific Arg483-Val484 peptide bond cleavage by matriptase-1, required for MSP activation, and mass spectrometry of tryptic fragments of the mutated protein showed that the free thiol introduced by the R689C mutation did not form an aberrant disulfide bond. Together, the studies indicate that the missense SNP impairs MSP function by reducing its affinity to RON and perhaps through a secondary effect on in vivo concentration arising from reduced thermodynamic stability, resulting in down-regulation of the MSP/RON signaling pathway.

Propagation of Dynamic Changes in Barnase Upon Binding of Barstar: An NMR and Computational Study

NMR spectroscopy and computer simulations were used to examine changes in chemical shifts and in dynamics of the ribonuclease barnase that result upon binding to its natural inhibitor barstar. Although the spatial structures of free and bound barnase are very similar, binding results in changes of the dynamics of both fast side-chains, as revealed by (2)H relaxation measurements, and NMR chemical shifts in an extended beta-sheet that is located far from the binding interface. Both side-chain dynamics and chemical shifts are sensitive to variations in the ensemble populations of the inter-converting molecular states, which can escape direct structural observation. Molecular dynamics simulations of free barnase and barnase in complex with barstar, as well as a normal mode analysis of barnase using a Gaussian network model, reveal relatively rigid domains that are separated by the extended beta-sheet mentioned above. The observed changes in NMR parameters upon ligation can thus be rationalized in terms of changes in inter-domain dynamics and in populations of exchanging states, without measurable structural changes. This provides an alternative model for the propagation of a molecular response to ligand binding across a protein that is based exclusively on changes in dynamics.

Mechanism of allosteric regulation of transglutaminase 2 by GTP

Allosteric regulation is a fundamental mechanism of biological control. Here, we investigated the allosteric mechanism by which GTP inhibits cross-linking activity of transglutaminase 2 (TG2), a multifunctional protein, with postulated roles in receptor signaling, extracellular matrix assembly, and apoptosis. Our findings indicate that at least two components are involved in functionally coupling the allosteric site and active center of TG2, namely (i) GTP binding to mask a conformationally destabilizing switch residue, Arg-579, and to facilitate interdomain interactions that promote adoption of a compact, catalytically inactive conformation and (ii) stabilization of the inactive conformation by an uncommon H bond between a cysteine (Cys-277, an active center residue) and a tyrosine (Tyr-516, a residue located on a loop of the beta-barrel 1 domain that harbors the GTP-binding site). Although not essential for GTP-mediated inhibition of cross-linking, this H bond enhances the rate of formation of the inactive conformer.

Cooperative fluctuations point to the dimerization interface of p53 core domain

Elastic network models are used for investigation of the p53 core domain functional dynamics. Global modes of motion indicate high positive correlations for residue fluctuations across the A-B interface, which are not observed at the B-C interface. Major hinge formation is observed at the A-B interface upon dimerization indicating stability of the A-B dimer. These findings imply A-B as the native dimerization interface, whereas B-C is the crystal interface. The A-B dimer exhibits an opening-closing motion about DNA, supporting the previously suggested clamp-like model of nonspecific DNA binding followed by diffusion. Monomer A has limited positive correlations with DNA, while monomer B exhibits high positive correlations with DNA in the functionally significant slow modes. Thus, monomer B might seem to maintain the stability of the dimer-DNA complex by forming the relatively fixed arm of the dimer clamp, whereas the other arm of the clamp, monomer A, might allow sliding via continuous association/dissociation mechanisms.

Reconstruction of the Conserved β-Bulge in Mammalian Defensins Using d-Amino Acids

Defensins are cationic antimicrobial mini-proteins that play important roles in the innate immune defense against microbial infection. Six invariant Cys residues in each defensin form three structurally indispensable intramolecular disulfide bridges. The only other residue invariant in all known mammalian defensins is a Gly. Structural studies indicate that the invariant Gly residue is located in an atypical, classic-type beta-bulge with the backbone torsion angles (Phi, Psi) disallowed for L-amino acids but permissible for D-enantiomers. We replaced the invariant Gly17 residue in human neutrophil alpha-defensin 2 (HNP2) by L-Ala or one of the D-amino acids Ala, Glu, Phe, Arg, Thr, Val, or Tyr. Although L-Ala17-HNP2 could not be folded, resulting in massive aggregation, all of the D-amino acid-substituted analogs folded with high efficiency. The high resolution x-ray crystal structures of dimeric D-Ala17-HNP2 were determined in three different crystal forms, showing a well preserved beta-bulge identical to those found in other defensins. The seven D-analogs of HNP2 exhibited highly variable bactericidal activity against Gram-positive and Gram-negative test strains, consistent with the premise that interplay between charge and hydrophobicity dictates how amphiphilic defensins kill. Further, the bactericidal activity of these d-amino acid analogs of HNP2 correlated well with their ability to induce leakage from large unilamellar vesicles, supporting membrane permeabilization as the lethal event in microbial killing by HNP2. Our findings identify a conformational prerequisite in the beta-bulge of defensins essential for correct folding and native structure, thereby explaining the molecular basis of the Gly-Xaa-Cys motif conserved in all mammalian defensins.

Contribution of single tryptophan residues to the fluorescence and stability of ribonuclease Sa

Ribonuclease Sa (RNase Sa) contains no tryptophan (Trp) residues. We have added single Trp residues to RNase Sa at sites where Trp is found in four other microbial ribonucleases, yielding the following variants of RNase Sa: Y52W, Y55W, T76W, and Y81W. We have determined crystal structures of T76W and Y81W at 1.1 and 1.0 A resolution, respectively. We have studied the fluorescence properties and stabilities of the four variants and compared them to wild-type RNase Sa and the other ribonucleases on which they were based. Our results should help others in selecting sites for adding Trp residues to proteins. The most interesting findings are: 1), Y52W is 2.9 kcal/mol less stable than RNase Sa and the fluorescence intensity emission maximum is blue-shifted to 309 nm. Only a Trp in azurin is blue-shifted to a greater extent (308 nm). This blue shift is considerably greater than observed for Trp71 in barnase, the Trp on which Y52W is based. 2), Y55W is 2.1 kcal/mol less stable than RNase Sa and the tryptophan fluorescence is almost completely quenched. In contrast, Trp59 in RNase T1, on which Y55W is based, has a 10-fold greater fluorescence emission intensity. 3), T76W is 0.7 kcal/mol more stable than RNase Sa, indicating that the Trp side chain has more favorable interactions with the protein than the threonine side chain. The fluorescence properties of folded Y76W are similar to those of the unfolded protein, showing that the tryptophan side chain in the folded protein is largely exposed to solvent. This is confirmed by the crystal structure of the T76W which shows that the side chain of the Trp is only approximately 7% buried. 4), Y81W is 0.4 kcal/mol less stable than RNase Sa. Based on the crystal structure of Y81W, the side chain of the Trp is 87% buried. Although all of the Trp side chains in the variants contribute to the unusual positive circular dichroism band observed near 235 nm for RNase Sa, the contribution is greatest for Y81W.

Essential domain motions in barnase revealed by MD simulations

The wealth of data accumulated on the bacterial ribonuclease barnase is complemented by molecular dynamics trajectories starting from four different experimental structures and covering a total of >10 ns. Using principal component analysis, the simulations are interpreted in view of dynamic domains and hinges promoting relative motions of these domains. Two domains with residues 7-22 and 52-108 for the first domain and residues 25-51 for the second domain were consistently observed. Hinge regions consist primarily of Tyr24, Ser50, Ile51, and Gly52. Earlier mutation studies have demonstrated that the residues of the hinge regions play essential roles for the stability and activity of barnase. The domain motions are correlated to inter-domain interactions involving functionally important active site residues, such as Lys27 and Glu73. A model is presented that combines the observation of dynamic domains and their motions with the extensive mutation data from the literature. Enthalpic energy contributions originating from specific inter-domain interactions as well as entropic energy contributions due to the domain motions are discussed in the frame of this model and compared with destabilization energies measured for corresponding mutants.

The role of a beta-bulge in the folding of the beta-hairpin structure in ubiquitin

It is known that the peptide corresponding to the N-terminal beta-hairpin of ubiquitin, U(1-17), can populate the monomeric beta-hairpin conformation in aqueous solution. In this study, we show that the Gly-10 that forms the bulge of the beta-turn in this hairpin is very important to the stability of the hairpin. The deletion of this residue to desG10(1-16) unfolds the structure of the peptide in water. Even under denaturing conditions, this bulge appears to be important in maintaining the residual structure of ubiquitin, which involves tertiary interactions within the sequence 1 to 34 in the denatured state. We surmise that this residual structure functions as one of the nucleation centers in the folding process and is important in stabilizing the transition state. In accordance with this idea, deleting Gly-10 slows down the refolding and unfolding rate by about one half.

cis-Inhibition of proteasomal degradation by viral repeats: impact of length and amino acid composition

The Gly-Ala repeat (GAr) of the Epstein-Barr virus nuclear antigen 1 is a cis acting inhibitor of ubiquitin-proteasome proteolysis. We have investigated the capacity of various repeats to inhibit the turnover of the proteasomal substrate IkappaBalpha. Inhibition of TNFalpha-induced degradation was achieved by insertion of octamers containing three alanines or valines, interspersed by no more then three consecutive glycines. The inhibitory activity was abolished by increasing the length of the spacer, by eliminating the spacers, or by substitution of a single hydrophobic residue with a polar or charged residue. A serine containing octamer was inactive but inhibition was partially restored by insertion of three consecutive repeats. These findings suggest a model where inhibition requires the interaction of at least three alanine residues of the GAr in a beta-strand conformation with adjacent hydrophobic binding pockets of a putative receptor.

Extreme functional sensitivity to conservative amino acid changes on enzyme exteriors

Mutagenesis studies and alignments of homologous sequences have demonstrated that protein function typically is compatible with a variety of amino-acid residues at most exterior non-active-site positions. These observations have led to the current view that functional constraints on sequence are minimal at these positions. Here, it is shown that this inference assumes that the set of acceptable residues at each position is independent of the overall sequence context. Two approaches are used to test this assumption. First, highly conservative replacements of exterior residues, none of which would cause significant functional disruption alone, are combined until roughly one in five have been changed. This is found to cause complete loss of function in vivo for two unrelated monomeric enzymes: barnase (a bacterial RNase) and TEM-1 beta-lactamase. Second, a set of hybrid sequences is constructed from the 50 %-identical TEM-1 and Proteus mirabilis beta-lactamases. These hybrids match the TEM-1 sequence except for a region at the C-terminal end, where they are random composites of the two parents. All of these hybrids are biologically inactive. In both experiments, complete loss of activity demonstrates the importance of sequence context in determining whether substitutions are functionally acceptable. Contrary to the prevalent view, then, enzyme function places severe constraints on residue identities at positions showing evolutionary variability, and at exterior non-active-site positions, in particular. Homologues sharing less than about two-thirds sequence identity should probably be viewed as distinct designs with their own sets of optimising features.