Critical regions for the sweetness of brazzein

Brazzein is a small, heat-stable, intensely sweet protein consisting of 54 amino acid residues. Based on the wild-type brazzein, 25 brazzein mutants have been produced to identify critical regions important for sweetness. To assess their sweetness, psychophysical experiments were carried out with 14 human subjects. First, the results suggest that residues 29-33 and 39-43, plus residue 36 between these stretches, as well as the C-terminus are involved in the sweetness of brazzein. Second, charge plays an important role in the interaction between brazzein and the sweet taste receptor.     

Temperature‐dependent conformational change affecting Tyr11 and sweetness loops of brazzein

The sweet protein brazzein, a member of the Csβα fold family, contains four disulfide bonds that lend a high degree of thermal and pH stability to its structure. Nevertheless, a variable temperature study has revealed that the protein undergoes a local, reversible conformational change between 37 and 3°C with a midpoint about 27°C that changes the orientations and side‐chain hydrogen bond partners of Tyr8 and Tyr11. To test the functional significance of this effect, we used NMR saturation transfer to investigate the interaction between brazzein and the amino terminal domain of the sweet receptor subunit T1R2; the results showed a stronger interaction at 7°C than at 37°C. Thus the low temperature conformation, which alters the orientations of two loops known to be critical for the sweetness of brazzein, may represent the bound state of brazzein in the complex with the human sweet receptor. Proteins 2013; © 2012 Wiley Periodicals, Inc.     

Monkey electrophysiological and human psychophysical responses to mutants of the sweet protein brazzein: delineating brazzein sweetness

Responses to brazzein, 25 brazzein mutants and two forms of monellin were studied in two types of experiments: electrophysiological recordings from chorda tympani S fibers of the rhesus monkey, Macaca mulatta, and psychophysical experiments. We found that different mutations at position 29 (changing Asp29 to Ala, Lys or Asn) made the molecule significantly sweeter than brazzein, while mutations at positions 30 or 33 (Lys30Asp or Arg33Ala) removed all sweetness. The same pattern occurred again at the beta-turn region, where Glu41Lys gave the highest sweetness score among the mutants tested, whereas a mutation two residues distant (Arg43Ala) abolished the sweetness. The effects of charge and side chain size were examined at two locations, namely positions 29 and 36. The findings indicate that charge is important for eliciting sweetness, whereas the length of the side-chain plays a lesser role. We also found that the N- and C-termini are important for the sweetness of brazzein. The close correlation (r = 0.78) between the results of the above two methods corroborates our hypothesis that S fibers convey sweet taste in primates.     

Synthesis and characterization of the sweet protein brazzein

The sweet protein brazzein isolated from the fruit of the African plant, Pentadiplandra brazzeana Baillon is 2000-500 times sweeter than sucrose, and consists of 54 amino acid residues with four intramolecular disulfide bonds. Brazzein was prepared by the fluoren-9-yl-methoxycarbonyl solid-phase method, and was identical to natural brazzein by high performance liquid chromatography, mass spectroscopy, peptide mapping, and taste evaluation. The D enantiomer of brazzein was also synthesized, and was shown to be the mirror image of brazzein. The D enantiomer (ent-brazzein) was devoid of any sweetness and was essentially tasteless. © 1996 John Wiley & Sons, Inc.     

Probing the sweet determinants of brazzein: wild-type brazzein and a tasteless variant, brazzein-ins(R18a-I18b), exhibit different pH-dependent NMR chemical shifts

Brazzein is a small, intensely sweet protein. As a probe of the functional properties of its solvent-exposed loop, two residues (Arg-Ile) were inserted between Leu18 and Ala19 of brazzein. Psychophysical testing demonstrated that this mutant is totally tasteless. NMR chemical shift mapping of differences between this mutant and brazzein indicated that residues affected by the insertion are localized to the mutated loop, the region of the single alpha-helix, and around the Cys16-Cys37 disulfide bond. Residues unaffected by this mutation included those near the C-terminus and in the loop connecting the alpha-helix and the second beta-strand. In particular, several residues of brazzein previously shown to be essential for its sweetness (His31, Arg33, Glu41, Arg43, Asp50, and Tyr54) exhibited negligible chemical shift changes. Moreover, the pH dependence of the chemical shifts of His31, Glu41, Asp50, and Tyr54 were unaltered by the insertion. The insertion led to large chemical shift and pKa perturbation of Glu36, a residue shown previously to be important for brazzein's sweetness. These results serve to refine the known sweetness determinants of brazzein and lend further support to the idea that the protein interacts with a sweet-taste receptor through a multi-site interaction mechanism, as has been postulated for brazzein and other sweet proteins (monellin and thaumatin).     

Sweetness determinant sites of brazzein, a small, heat-stable, sweet-tasting protein

Brazzein, originally isolated from the fruit of the African plant Pentadiplandra brazzeana Baillon, is the smallest, most heat-stable and pH-stable member of the set of proteins known to have intrinsic sweetness. These properties make brazzein an ideal system for investigating the chemical and structural requirements of a sweet-tasting protein. We have used the three-dimensional structure of the protein (J. E. Caldwell et al. (1998) Nat. Struct. Biol. 5, 427-431) as a guide in designing 15 synthetic genes in expression constructs aimed at delineating the sweetness determinants of brazzein. Protein was produced heterologously in Escherichia coli, isolated, and purified as described in the companion paper (Assadi-Porter, F. M., Aceti, D., Cheng, H., and Markley, J. L., this issue). Analysis by one-dimensional (1)H NMR spectroscopy indicated that all but one of these variants had folded properly under the conditions used. A taste panel compared the gustatory properties of solutions of these proteins to those of sucrose and brazzein isolated from fruit. Of the 14 mutations in the des-pGlu1-brazzein background, four exhibited almost no sweetness, six had significantly reduced sweetness, two had taste properties equivalent to des-pGlu1-brazzein (two times as sweet as the major form of brazzein isolated from fruit which contains pGlu1), and two were about twice as sweet as des-pGlu1-brazzein. Overall, the results suggest that two regions of the protein are critical for the sweetness of brazzein: a region that includes the N- and C-termini of the protein, which are located close to one another, and a region that includes the flexible loop around Arg43.     

Hydrogen bonds in rubredoxins from mesophilic and hyperthermophilic organisms

The extent and strength of the hydrogen bond networks in rubredoxins from the hyperthermophile Pyrococcus furiosus (PfRd), and its mesophilic analogue Clostridium pasteurianum (CpRd), are examined and compared using NMR spectroscopy. NMR parameters examined in this study include through-hydrogen bond (h3)J(NC)(') scalar couplings and (1)H, (13)C, and (15)N chemical shifts, as well as covalent (1)J(NH) and (1)J(NC)(') scalar couplings. These parameters have allowed the characterization in solution of 12 hydrogen bonds in each protein. Despite a 83% sequence homology and a low RMSD for the backbone heavy atoms (0.648 A) in the crystalline state, subtle, but definite, changes have been identified in the detailed hydrogen-bonding patterns. CpRd shows an increased number of hydrogen bonds in the triple-stranded beta-sheet and an additional hydrogen bond in the multiple-turn segment including residues 14-32. On the other hand, PfRd exhibits an overall strengthening of N-H...O=C hydrogen bonds in the loops involved at the metal binding site as well as evidence for an additional NH...S(Cys) hydrogen bond involving the alanine residue 44. These data, as well as temperature dependence of the NMR parameters, suggest that the particular NMR hydrogen bond pattern found in the hyperthermophile rubredoxin leads to an increased stabilization at the metal binding pocket. It seems to result from a subtle redistribution of hydrogen-bonding interactions between the triple-stranded beta-sheet and the actual metal binding site.     

Expression of the sweet protein brazzein in maize for production of a new commercial sweetener

The availability of foods low in sugar content yet high in flavour is critically important to millions of individuals conscious of carbohydrate intake for diabetic or dietetic purposes. Brazzein is a sweet protein occurring naturally in a tropical plant that is impractical to produce economically on a large scale, thus limiting its availability for food products. We report here the use of a maize expression system for the production of this naturally sweet protein. High expression of brazzein was obtained, with accumulation of up to 4% total soluble protein in maize seed. Purified corn brazzein possessed a sweetness intensity of up to 1200 times that of sucrose on a per weight basis. In addition, application tests demonstrated that brazzein-containing maize germ flour could be used directly in food applications, providing product sweetness. These results demonstrate that high-intensity sweet protein engineered into food products can give sweetener attributes useful in the food industry.     

Improvement of hydrogen bond geometry in protein NMR structures by residual dipolar couplings – an assessment of the interrelation of NMR restraints

We have examined how the hydrogen bond geometry in three different proteins is affected when structural restraints based on measurements of residual dipolar couplings are included in the structure calculations. The study shows, that including restraints based solely on (1)H(N)-(15)N residual dipolar couplings has pronounced impact on the backbone rmsd and Ramachandran plot but does not improve the hydrogen bond geometry. In the case of chymotrypsin inhibitor 2 the addition of (13)CO-(13)C(alpha) and (15)N-(13)CO one bond dipolar couplings as restraints in the structure calculations improved the hydrogen bond geometry to a quality comparable to that obtained in the 1.8 A resolution X-ray structure of this protein. A systematic restraint study was performed, in which four types of restraints, residual dipolar couplings, hydrogen bonds, TALOS angles and NOEs, were allowed in two states. This study revealed the importance of using several types of residual dipolar couplings to get good hydrogen bond geometry. The study also showed that using a small set of NOEs derived only from the amide protons, together with a full set of residual dipolar couplings resulted in structures of very high quality. When reducing the NOE set, it is mainly the side-chain to side-chain NOEs that are removed. Despite of this the effect on the side-chain packing is very small when a reduced NOE set is used, which implies that the over all fold of a protein structure is mainly determined by correct folding of the backbone.     

Pressure alters electronic orbital overlap in hydrogen bonds

Pressure-induced changes in ^3h J _NC scalar couplings through hydrogen bonds were investigated in the immunoglobulin binding domain of streptococcal protein G. ¹H, ¹⁵N and ¹³C triple-resonance NMR spectroscopy coupled with the on-line high pressure cell technique was used to monitor ^3h J _NC scalar couplings at 30 and 2000 bar in uniformly labeled ¹⁵N and ¹³C protein isotopes. Both increased and decreased ^3h J _NC scalar couplings were observed at high pressure. No correlation with secondary structure was apparent. The difference in coupling constants as well as pressure-induced chemical shift data suggests a compaction of the helix ends and an increase of the helix pitch at its center in response to pressure. Our data provides the first direct evidence that the electronic orbital overlap in protein backbone hydrogen bonds is altered by pressure.     

Hydrogen bonds between short polar side chains and peptide backbone: prevalence in proteins and effects on helix-forming propensities

A survey of 322 proteins showed that the short polar (SP) side chains of four residues, Thr, Ser, Asp, and Asn, have a very strong tendency to form hydrogen bonds with neighboring backbone amides. Specifically, 32% of Thr, 29% of Ser, 26% of Asp, and 19% of Asn engage in such hydrogen bonds. When an SP residue caps the N terminal of a helix, the contribution to helix stability by a hydrogen bond with the amide of the N3 or N2 residue is well established. When an SP residue is in the middle of a helix, the side chain is unlikely to form hydrogen bonds with neighboring backbone amides for steric and geometric reasons. In essence the SP side chain competes with the backbone carbonyl for the same hydrogen-bonding partner (i.e., the backbone amide) and thus SP residues tend to break backbone carbonyl-amide hydrogen bonds. The proposition that this is the origin for the low propensities of SP residues in the middle of alpha helices (relative to those of nonpolar residues) was tested. The combined effects of restricting side-chain rotamer conformations (documented by Creamer and Rose, Proc Acad Sci USA, 1992;89:5937-5941; Proteins, 1994;19:85-97) and excluding side- chain to backbone hydrogen bonds by the helix were quantitatively analyzed. These were found to correlate strongly with four experimentally determined scales of helix-forming propensities. The correlation coefficients ranged from 0.72 to 0.87, which are comparable to those found for nonpolar residues (for which only the loss of side-chain conformational entropy needs to be considered).     

Effects of backbone contacts 3' to the abasic site on the cleavage and the product binding by human apurinic/apyrimidinic endonuclease (APE1)

The mammalian apurinic/apyrimidinic (AP) endonuclease (APE1) is a multifunctional protein that plays essential roles in DNA repair and gene regulation. We decomposed the APEs into 12 blocks of highly conserved sequence and structure (molegos). This analysis suggested that residues in molegos common to all APEs, but not to the less specific nuclease, DNase I, would dictate enhanced binding to damaged DNA. To test this hypothesis, alanine was substituted for N226 and N229, which form hydrogen bonds to the DNA backbone 3' of the AP sites in crystal structures of the APE1/DNA complex. While the cleavage rate at AP sites of both N226A and N229A mutants increased, their ability to bind to damaged DNA decreased. The ability of a double mutant (N226A/N229A) to bind damaged DNA was further decreased, while the V(max) was almost identical to that of the wild-type APE1. A double mutant at N226 and R177, a residue that binds to the same phosphate as N229, had a significantly decreased activity and substrate binding. As the affinity for product DNA was decreased in all the mutants, the enhanced reaction rate of the single mutants could be due to alleviation of product inhibition of the enzyme. We conclude that hydrogen bonds to phosphate groups 3' to the cleavage site is essential for APE1's binding to the product DNA, which may be necessary for efficient functioning of the base excision repair pathway. The results indicate that the molego analysis can aid in the redesign of proteins with altered binding affinity and activity.     

Large increase in brazzein expression achieved by changing the plasmid /strain combination of the NICE system in Lactococcus lactis

Aims:  To evaluate brazzein production in Lactococcus lactis using the nisin-controlled expression (NICE) system. The approach is through analysis of different plasmid/strain combinations.
Methods and Results:  Two plasmid/strain combinations of the NICE system were used in brazzein expression: L. lactis NZ9000 harbouring plasmid pNZ8148, and L. lactis IL1403 harbouring plasmid pMSP3545. The former combination proved superior, with a >800-fold increase in His-tagged brazzein expression (to 1·65 mg l⁻¹ of fermentation broth), comparable to expression levels in Escherichia coli . Improved expression resulted in a minor increase in secretion to the medium with the use of the Usp45 signal peptide. The yield of wild-type brazzein corresponded to that of His-tagged brazzein. Wild-type brazzein was partially soluble and low-intensity sweetness was detected.
Conclusions:  The plasmid/strain combination of the NICE system has a significant impact on the expression of brazzein where a >800-fold increase was achieved. The greatly increased expression of brazzein resulted in minor improvement in secretion and low-intensity sweetness.
Significance and Impact of the Study:  The choice of the plasmid/strain combination of the NICE system was shown to be of extreme importance in brazzein expression.     

Local unfolding of Cu, Zn superoxide dismutase monomer determines the morphology of fibrillar aggregates

Aggregation of Cu, Zn superoxide dismutase (SOD1) is often found in amyotrophic lateral sclerosis patients. The fibrillar aggregates formed by wild type and various disease-associated mutants have recently been found to have distinct cores and morphologies. Previous computational and experimental studies of wild-type SOD1 suggest that the apo-monomer, highly aggregation prone, displays substantial local unfolding dynamics. The residual folded structure of locally unfolded apoSOD1 corresponds to peptide segments forming the aggregation core as identified by a combination of proteolysis and mass spectroscopy. Therefore, we hypothesize that the destabilization of apoSOD1 caused by various mutations leads to distinct local unfolding dynamics. The partially unfolded structure, exposing the hydrophobic core and backbone hydrogen bond donors and acceptors, is prone to aggregate. The peptide segments in the residual folded structures form the "building block" for aggregation, which in turn determines the morphology of the aggregates. To test this hypothesis, we apply a multiscale simulation approach to study the aggregation of three typical SOD1 variants: wild type, G37R, and I149T. Each of these SOD1 variants has distinct peptide segments forming the core structure and features different aggregate morphologies. We perform atomistic molecular dynamics simulations to study the conformational dynamics of apoSOD1 monomer and coarse-grained molecular dynamics simulations to study the aggregation of partially unfolded SOD1 monomers. Our computational studies of monomer local unfolding and the aggregation of different SOD1 variants are consistent with experiments, supporting the hypothesis of the formation of aggregation "building blocks" via apo-monomer local unfolding as the mechanism of SOD1 fibrillar aggregation.     

Quantification of H/D Isotope Effects on Protein Hydrogen-bonds by h3JNC′ and 1JNC′ Couplings and Peptide Group 15N and 13C′ Chemical Shifts

The effect of hydrogen/deuterium exchange on protein hydrogen bond coupling constants (h3)J(NC') has been investigated in the small globular protein ubiquitin. The couplings across deuterated or protonated hydrogen bonds were measured by a long-range quantitative HA(CACO)NCO experiment. The analysis is combined with a determination of the H(N)/D(N) isotope effect on the amide group (1)J(NC') couplings and the (15)N and (13)C' chemical shifts. On average, H-bond deuteration exchange weakens (h3)J(NC') and strengthens (1)J(NC') couplings. A correlation is found between the size of the (15)N isotope shift, the (15)N chemical shift, and the (h3)J(NC') coupling constants. The data are consistent with a reduction of donor-acceptor overlap as expected from the classical Ubbelohde effect and the common understanding that H(N)/D(N) exchange leads to a shortening of the N-hydron bond length.     

Temperature-dependence of protein hydrogen bond properties as studied by high-resolution NMR

The temperature-dependence of a large number of NMR parameters describing hydrogen bond properties in the protein ubiquitin was followed over a range from 5 to 65 degrees C. The parameters comprise hydrogen bond (H-bond) scalar couplings, h3JNC', chemical shifts, amide proton exchange rates, 15N relaxation parameters as well as covalent 1JNC' and 1JNH couplings. A global weakening of the h3JNC' coupling with increasing temperature is accompanied by a global upfield shift of the amide protons and a decrease of the sequential 1JNC' couplings. If interpreted as a linear increase of the N...O distance, the change in h3JNC' corresponds to an average linear thermal expansion coefficient for the NH-->O hydrogen bonds of 1.7 x 10(-4)/K, which is in good agreement with overall volume expansion coefficients observed for proteins. A residue-specific analysis reveals that not all hydrogen bonds are affected to the same extent by the thermal expansion. The end of beta-sheet beta1/beta5 at hydrogen bond E64-->Q2 appears as the most thermolabile, whereas the adjacent hydrogen bond I3-->L15 connecting beta-strands beta1 and beta2 is even stabilized slightly at higher temperatures. Additional evidence for the stabilization of the beta1/beta2 beta-hairpin at higher temperatures is found in reduced hydrogen exchange rates for strand end residue V17. This reduction corresponds to a stabilizing change in free energy of 9.7 kJ/mol for the beta1/beta2 hairpin. The result can be linked to the finding that the beta1/beta2 hairpin behaves as an autonomously folding unit in the A-state of ubiquitin under changed solvent conditions. For several amide groups the temperature-dependencies of the amide exchange rates and H-bond scalar couplings are uncorrelated. Therefore, amide exchange rates are not a sole function of the hydrogen bond "strength" as given by the electronic overlap of donors and acceptors, but are clearly dependent on other blocking mechanisms.     

Peptaibol zervamicin IIb structure and dynamics refinement from transhydrogen bond J couplings

Zervamicin IIB (Zrv-IIB) is a channel-forming peptaibol antibiotic of fungal origin. The measured transhydrogen bond (3h)J(NC') couplings in methanol solution heaving average value of -0.41 Hz indicate that the stability of the Zrv-IIB helix in this milieu is comparable to the stability of helices in globular proteins. The N-terminus of the peptide forms an alpha-helix, whereas 3(10)-helical hydrogen bonds stabilize the C-terminus. However, two weak transhydrogen bond peaks are observed in a long-range HNCO spectrum for HN Aib(12). Energy calculations using the Empirical Conformation Energy Program for Peptides (ECEPP)/2 force field and the implicit solvent model show that the middle of the peptide helix accommodates a bifurcated hydrogen bond that is simultaneously formed between HN Aib(12) and CO Leu(8) and CO Aib(9). Several lowered (3h)J(NC') on a polar face of the helix correlate with the conformational exchange process observed earlier and imply dynamic distortions of a hydrogen bond pattern with the predominant population of a properly folded helical structure. The refined structure of Zrv-IIB on the basis of the observed hydrogen bond pattern has a small ( approximately 20 degrees ) angle of helix bending that is virtually identical to the angle of bending in dodecylphosphocholine (DPC) micelles, indicating the stability of a hinge region in different environments. NMR parameters ((1)HN chemical shifts and transpeptide bond (1)J(NC') couplings) sensitive to hydrogen bonding along with the solvent accessible surface area of carbonyl oxygens indicate a large polar patch on the convex side of the helix formed by three exposed backbone carbonyls of Aib(7), Aib(9), and Hyp(10) and polar side chains of Hyp(10), Gln(11), and Hyp(13). The unique structural features, high helix stability and the enhanced polar patch, set apart Zrv-IIB from other peptaibols (for example, alamethicin) and possibly underlie its biological and physiological properties.     

Ligand-induced strain in hydrogen bonds of the c-Src SH3 domain detected by NMR

Changes in the molecular conformation of proteins can result from a variety of perturbations, and can play crucial roles in the regulation of biological activity. A new solution NMR method has been applied to monitor ligand-induced changes in hydrogen bond geometry in the chicken c-Src SH3 domain. The structural response of this domain to ligand binding has been investigated by measuring trans-hydrogen bond (15)N-(13)C' scalar couplings in the free state and when bound to the high affinity class I ligand RLP2, containing residues RALPPLPRY. A comparison between hydrogen bonds in high resolution X-ray structures of this domain and those observed via (h3)J(NC') couplings in solution shows remarkable agreement. Two backbone-to-side-chain hydrogen bonds are observed in solution, and each appears to play a role in stabilization of loop structure. Reproducible ligand-induced changes in trans-hydrogen bond scalar couplings are observed across the domain that translate into changes in hydrogen bond length ranging between 0.02 to 0.12 A. The observed changes can be rationalized by an induced fit mechanism in which hydrogen bonds across the protein participate in a compensatory response to forces imparted at the protein-ligand interface. Upon ligand binding, mutual intercalation of the two Leu-Pro segments of the ligand between three aromatic side-chains protruding from the SH3 surface wedges apart secondary structural elements within the SH3 domain. This disruption is transmitted in a domino-like effect across the domain through networks of hydrogen bonded peptide planes. The unprecedented resolution obtained demonstrates the ability to characterize subtle structural rearrangements within a protein upon perturbation, and represents a new step in the endeavor to understand how hydrogen bonds contribute to the stabilization and function of biological macromolecules.     

Helix bending in alamethicin: molecular dynamics simulations and amide hydrogen exchange in methanol

Molecular dynamics simulations of alamethicin in methanol were carried out with either a regular alpha-helical conformation or the x-ray crystal structure as starting structures. The structures rapidly converged to a well-defined hydrogen-bonding pattern with mixed alpha-helical and 3(10)-helical hydrogen bonds, consistent with NMR structural characterization, and did not unfold throughout the 1-ns simulation, despite some sizable backbone fluctuations involving reversible breaking of helical hydrogen bonds. Bending of the helical structure around residues Aib10-Aib13 was associated with reversible flips of the peptide bonds involving G11 (Aib10-G11 or G11-L12 peptide bonds), yielding discrete structural states in which the Aib10 carbonyl or (rarely) the G11 carbonyl was oriented away from the peptide helix. These peptide bond reversals could be accommodated without greatly perturbing the adjacent helical structure, and intramolecular hydrogen bonding was generally maintained in bent states through the formation of new (non-alpha or 3[10]) hydrogen bonds with good geometries: G11 NH-V9 CO (inverse gamma turn), Aib13 NH-Aib8 CO (pi-helix) and, rarely, L12 NH- Q7 NH (pi-helix). These observations may reconcile potentially conflicting NMR structural information for alamethicin in methanol, in which evidence for conformational flexibility in the peptide sequence before P14 (G11-Aib13) contrasts with the stability of backbone amide NH groups to exchange with solvent. Similar reversible reorientation of the Thr11-Gly12 peptide bond of melittin is also observed in dynamics simulations in methanol (R. B. Sessions, N. Gibbs, and C. E. Dempsey, submitted). This phenomenon may have some role in the orientation of the peptide carbonyl in solvating the channel lumen in membrane ion channel states of these peptides.