Probing pH-dependent functional elements in proteins: modification of carboxylic acid pairs in Trichoderma reesei cellobiohydrolase Cel6A

Two carboxylic acid side chains can, depending on their geometry and environment, share a proton in a hydrogen bond and form a carboxyl-carboxylate pair. In the Trichoderma reesei cellobiohydrolase Cel6A structure, five carboxyl-carboxylate pairs are observed. One of these pairs (D175-D221) is involved in catalysis, and three other pairs are found in, or close to the two surface loops covering the active site tunnel of the catalytic domain. To stabilize Cel6A at alkaline pH values, where deprotonation of the carboxylic acids leads to repulsion of their side chains, we designed two mutant enzymes. In the first mutant, one carboxyl-carboxylate pair (E107-E399) was replaced by a corresponding amide-carboxylate pair (Q107-E399), and in the second mutant, all three carboxyl-carboxylate pairs (E107-E399, D170-E184, and D366-D419) were mutated in a similar manner. The unfolding studies using both intrinsic tryptophan fluorescence and far-ultraviolet circular dichroism spectroscopy at different pH values demonstrate that the unfolding temperature (T(m)) of both mutants has changed, resulting in destabilization of the mutant enzymes at acidic pH and stabilization at alkaline pH. The effect of stabilization seems additive, as a Cel6A triple mutant is the most stable enzyme variant. This increased stability is also reflected in the 2- or 4-fold increased half-life of the two mutants at alkaline pH, while the catalytic rate on cellotetraose (at t = 0) has not changed. Increased operational stability at alkaline pH was also observed on insoluble cellulosic substrates. Local conformational changes are suggested to take place in the active site loops of Cel6A wild-type enzyme at elevated pHs (pH 7), affecting to the end-product spectrum on insoluble cellulose. The triple mutant does not show such pH-dependent behavior. Overall, our results demonstrate that carboxyl-carboxylate pair engineering is a useful tool to alter pH-dependent protein behavior.     

Nearest-neighbor effects and structural preferences in dipeptides are a function of the electronic properties of amino acid side-chains

The electronic properties of amino acid side-chains are emerging as an important factor in the preference for secondary structure in proteins. These properties have not been fully characterized, nor has their role in the behavior of peptides been explored in any detail. The present studies sought to evaluate several possibilities: 1) that hydrophilicity can be expressed solely in electronic terms, 2) that substituent effects of side-chains extend across the peptide bond, and (3) nearest-neighbor effects in dipeptides correlate with secondary structural preferences. Quantum mechanics (QM) calculations were used to define the electronic properties of individual amino acids and dipeptides. It was found that the hydrophilicity of an amino acid side-chain can be accurately represented as a function of the electron densities of its component atoms. In addition, the nature of an amino acid in the second position of a dipeptide affects the electronic properties (Mulliken populations and electron densities) of the main-chain atoms of the first residue. Certain electronic features of the dipeptides strongly correlated with propensity for secondary structure. Specifically, Mulliken population data at the Calpha atom and N atom predicted preference for alpha-helices versus coil and strand conformations, respectively. Analysis of dipeptides arrayed in either helical or extended structures revealed lengthening of main-chain bonds in the alpha-helical conformations. A thorough characterization of the electronic properties of amino acids and short peptide segments may provide a better understanding of the forces that determine secondary structure in proteins.     

Simplified methods for pKa and acid pH-dependent stability estimation in proteins: removing dielectric and counterion boundaries

Much computational research aimed at understanding ionizable group interactions in proteins has focused on numerical solutions of the Poisson-Boltzmann (PB) equation, incorporating protein exclusion zones for solvent and counterions in a continuum model. Poor agreement with measured pKas and pH-dependent stabilities for a (protein, solvent) relative dielectric boundary of (4,80) has lead to the adoption of an intermediate (20,80) boundary. It is now shown that a simple Debye-Huckel (DH) calculation, removing both the low dielectric and counterion exclusion regions associated with protein, is equally effective in general pKa calculations. However, a broad-based discrepancy to measured pH-dependent stabilities is maintained in the absence of ionizable group interactions in the unfolded state. A simple model is introduced for these interactions, with a significantly improved match to experiment that suggests a potential utility in predicting and analyzing the acid pH-dependence of protein stability. The methods are applied to the relative pH-dependent stabilities of the pore-forming domains of colicins A and N. The results relate generally to the well-known preponderance of surface ionizable groups with solvent-mediated interactions. Although numerical PB solutions do not currently have a significant advantage for overall pKa estimations, development based on consideration of microscopic solvation energetics in tandem with the continuum model could combine the large deltapKas of a subset of ionizable groups with the overall robustness of the DH model.     

The role of a conserved tyrosine residue in high-potential iron sulfur proteins.

Conserved tyrosine-12 of Ectothiorhodospira halophila high-potential iron sulphur protein (HiPIP) iso-I was substituted with phenylalanine (Y12F), histidine (Y12H), tryptophan (Y12W), isoleucine (Y12I), and alanine (Y12A). Variants Y12A and Y12I were expressed to reasonable levels in cells grown at lower temperatures, but decomposed during purification. Variants Y12F, Y12H, and Y12W were substantially destabilized with respect to the recombinant wild-type HiPIP (rcWT) as determined by differential scanning calorimetry over a pH range of 7.0-11.0. Characterization of the Y12F variant by NMR indicates that the principal structural differences between this variant and the rcWT HiPIP result from the loss of the two hydrogen bonds of the Tyr-12 hydroxyl group with Asn-14 O delta 1 and Lys-59 NH, respectively. The effect of the loss of the latter interaction is propagated through the Lys-59/Val-58 peptide bond, thereby perturbing Gly-46. The delta delta GDapp of Y12F of 2.3 kcal/mol with respect to rcWT HiPIP (25 degrees C, pH 7.0) is entirely consistent with the contribution of these two hydrogen bonds to the stability of the latter. CD measurements show that Tyr-12 influences several electronic transitions within the cluster. The midpoint reduction potentials of variants Y12F, Y12H, and Y12W were 17, 19, and 22 mV (20 mM MOPS, 0.2 M sodium chloride, pH 6.98, 25 degrees C), respectively, higher than that of rcWT HiPIP. The current results indicate that, although conserved Tyr-12 modulates the properties of the cluster, its principle function is to stabilize the HiPIP through hydrogen bonds involving its hydroxyl group and electrostatic interactions involving its aromatic ring.     

Deuterium isotope effects and fractionation factors of hydrogen-bonded A:T base pairs of DNA

Deuterium isotope effects and fractionation factors of N1...H3–N3 hydrogen bonded Watson–Crick A:T base pairs of two DNA dodecamers are presented here. Specifically, two-bond deuterium isotope effects on the chemical shifts of ¹³C2 and ¹³C4, ²¹³C2 and ²¹³C4, and equilibrium deuterium/protium fractionation factors of H3, , were measured and seen to correlate with the chemical shift of the corresponding imino proton, _H3. Downfield-shifted imino protons associated with larger values of ²¹³C2 and ²¹³C4 and smaller values, which together suggested that the effective H3–N3 vibrational potentials were more anharmonic in the stronger hydrogen bonds of these DNA molecules. We anticipate that ²¹³C2, ²¹³C4 and values can be useful gauges of hydrogen bond strength of A:T base pairs.     

Evaluation of C-H cdots,three dots,centered O hydrogen bonds in native and misfolded proteins

Non-traditional C-H cdots, three dots, centered Y hydrogen bonds, in which a carbon atom acts as the hydrogen donor and an electronegative atom Y (Y=N, O or S) acts as the acceptor, have been reported in proteins, but their importance in protein structures is not well established. Here, we present the results of three computational tests that examine the significance of C-H cdots, three dots, centered Y bonds involving the C(alpha) in proteins. First, we compared the number of C(alpha)-H cdots, three dots, centered Y bonds in native structures with two sets of compact, energy-minimized decoy structures. The decoy structures contain about as many C(alpha)-H cdots, three dots, centered Y bonds as the native structures, indicating that the constraints of chain connectivity and compactness can lead to incidental formation of C(alpha)-H cdots, three dots, centered Y bonds. Secondly, we examined whether short C(alpha)-H cdots, three dots, centered Y bonds have a tendency to be linear, as is expected for a cohesive hydrogen-bonding interaction. The native structures do show this trend, but so does one of the decoy sets, suggesting that this criterion is also not sufficient to indicate a stabilizing interaction. Finally, we examined the preference for C(alpha)-H cdots, three dots, centered Y bond donors to be near to strong hydrogen bond acceptors. In the native proteins, the alpha protons attract strong acceptors like oxygen atoms more than weak acceptors. In contrast, hydrogen bond donors in the decoy structures do not distinguish between strong and weak acceptors. Thus, any individual C(alpha)-H cdots, three dots, centered Y bond may be fortuitous and occur due to the polypeptide connectivity and compactness. Taken collectively, however, C(alpha)-H cdots, three dots, centered Y bonds provide a weakly cohesive force that stabilizes proteins.     

Hydrogen bonding markedly reduces the pK of buried carboxyl groups in proteins

The ionizable groups in proteins with the lowest pKs are the carboxyl groups of aspartic acid side-chains. One of the lowest, pK=0.6, is observed for Asp76 in ribonuclease T1. This low pK appeared to result from hydrogen bonds to a water molecule and to the side-chains of Asn9, Tyr11, and Thr91. The results here confirm this by showing that the pK of Asp76 increases to 1.7 in N9A, to 4.0 in Y11F, to 4.2 in T91V, to 4.4 in N9A+Y11F, to 4.9 in N9A+T91V, to 5.9 in Y11F+T91V, and to 6.4 in the triple mutant: N9A+Y11F+T91V. In ribonuclease Sa, the lowest pK=2.4 for Asp33. This pK increases to 3.9 in T56A, which removes the hydrogen bond to Asp33, and to 4.4 in T56V, which removes the hydrogen bond and replaces the -OH group with a -CH(3) group. It is clear that hydrogen bonds are able to markedly lower the pK values of carboxyl groups in proteins. These same hydrogen bonds make large contributions to the conformational stability of the proteins. At pH 7, the stability of D76A ribonuclease T1 is 3.8 kcal mol(-1) less than wild-type, and the stability of D33A ribonuclease Sa is 4.1 kcal mol(-1) less than wild-type. There is a good correlation between the changes in the pK values and the changes in stability. The results suggest that the pK values for these buried carboxyl groups would be greater than 8 in the absence of hydrogen bonds, and that the hydrogen bonds and other interactions of the carboxyl groups contribute over 8 kcal mol(-1) to the stability.     

Interhelical hydrogen bonds and spatial motifs in membrane proteins: polar clamps and serine zippers

Polar and ionizable amino acid residues are frequently found in the transmembrane (TM) regions of membrane proteins. In this study, we show that they help to form extensive hydrogen bond connections between TM helices. We find that almost all TM helices have interhelical hydrogen bonding. In addition, we find that a pair of contacting TM helices is packed tighter when there are interhelical hydrogen bonds between them. We further describe several spatial motifs in the TM regions, including "Polar Clamp" and "Serine Zipper," where conserved Ser residues coincide with tightly packed locations in the TM region. With the examples of halorhodopsin, calcium-transporting ATPase, and bovine cytochrome c oxidase, we discuss the roles of hydrogen bonds in stabilizing helical bundles in polytopic membrane proteins and in protein functions.     

Life at the border: Adaptation of proteins to anisotropic membrane environment

This review discusses main features of transmembrane (TM) proteins which distinguish them from water‐soluble proteins and allow their adaptation to the anisotropic membrane environment. We overview the structural limitations on membrane protein architecture, spatial arrangement of proteins in membranes and their intrinsic hydrophobic thickness, co‐translational and post‐translational folding and insertion into lipid bilayers, topogenesis, high propensity to form oligomers, and large‐scale conformational transitions during membrane insertion and transport function. Special attention is paid to the polarity of TM protein surfaces described by profiles of dipolarity/polarizability and hydrogen‐bonding capacity parameters that match polarity of the lipid environment. Analysis of distributions of Trp resides on surfaces of TM proteins from different biological membranes indicates that interfacial membrane regions with preferential accumulation of Trp indole rings correspond to the outer part of the lipid acyl chain region—between double bonds and carbonyl groups of lipids. These “midpolar” regions are not always symmetric in proteins from natural membranes. We also examined the hydrophobic effect that drives insertion of proteins into lipid bilayer and different free energy contributions to TM protein stability, including attractive van der Waals forces and hydrogen bonds, side‐chain conformational entropy, the hydrophobic mismatch, membrane deformations, and specific protein–lipid binding.     

A mechanistic analysis of the increase in the thermal stability of proteins in aqueous carboxylic acid salt solutions

The stability of proteins is known to be affected significantly in the presence of high concentration of salts and is highly pH dependent. Extensive studies have been carried out on the stability of proteins in the presence of simple electrolytes and evaluated in terms of preferential interactions and increase in the surface tension of the medium. We have carried out an in-depth study of the effects of a series of carboxylic acid salts: ethylene diamine tetra acetate, butane tetra carboxylate, propane tricarballylate, citrate, succinate, tartarate, malonate, and gluconate on the thermal stability of five different proteins that vary in their physico-chemical properties: RNase A, cytochrome c, trypsin inhibitor, myoglobin, and lysozyme. Surface tension measurements of aqueous solutions of the salts indicate an increase in the surface tension of the medium that is very strongly correlated with the increase in the thermal stability of proteins. There is also a linear correlation of the increase in thermal stability with the number of carboxylic groups in the salt. Thermal stability has been found to increase by as much as 22 C at 1 M concentration of salt. Such a high thermal stability at identical concentrations has not been reported before. The differences in the heat capacities of denaturation, deltaCp for RNase A, deduced from the transition curves obtained in the presence of varying concentrations of GdmCl and that of carboxylic acid salts as a function of pH, indicate that the nature of the solvent medium and its interactions with the two end states of the protein control the thermodynamics of protein denaturation. Among the physico-chemical properties of proteins, there seems to be an interplay of the hydrophobic and electrostatic interactions that lead to an overall stabilizing effect. Increase in surface free energy of the solvent medium upon addition of the carboxylic acid salts appears to be the dominant factor in governing the thermal stability of proteins.     

Amino acid pairing at the N- and C-termini of helical segments in proteins

A systematic survey was carried out in an unbiased sample of 815 protein chains with a maximum of 20% homology selected from the Protein Data Bank, whose structures were solved at a resolution higher than 1.6 A and with a R-factor lower than 25%. A set of 5556 subsequences with alpha-helix or 3(10)-helix motifs was extracted from the protein chains considered. Global and local propensities were then calculated for all possible amino acid pairs of the type (i, i + 1), (i, i + 2), (i, i + 3), and (i, i + 4), starting at the relevant helical positions N1, N2, N3, C3, C2, C1, and N-int (interior positions), and also at the first nonhelical positions in both termini of the helices, namely, N-cap and C-cap. The statistical analysis of the propensity values has shown that pairing is significantly dependent on the type of the amino acids and on the position of the pair. A few sequences of three and four amino acids were selected and their high prevalence in helices is outlined in this work. The Glu-Lys-Tyr-Pro sequence shows a peculiar distribution in proteins, which may suggest a relevant structural role in alpha-helices when Pro is located at the C-cap position. A bioinformatics tool was developed, which updates automatically and periodically the results and makes them available in a web site.     

Hydrogen exchange properties of proteins in native and denatured states monitored by mass spectrometry and NMR.

The extent of deuterium labeling of hen lysozyme, its three-disulfide derivative, and the homologous alpha-lactalbumins, has been measured by both mass spectrometry and NMR. Different conformational states of the proteins were produced by varying the solution conditions. Alternate protein conformers were found to contain different numbers of 2H atoms. Furthermore, measurement in the gas phase of the mass spectrometer or directly in solution by NMR gave consistent results. The unique ability of mass spectrometry to distinguish distributions of 2H atoms in protein molecules is exemplified using samples prepared to contain different populations of 2H-labeled protein. A comparison of the peak widths of bovine alpha-lactalbumin in alternate solution conformations but containing the same average number of 2H atoms showed dramatic differences due to different 2H distributions in the two protein conformers. Measurement of 2H distributions by ESI-MS enabled characterization of conformational averaging and structural heterogeneity. In addition, a time course for hydrogen exchange was examined and the variation in distributions of 2H atom compared with simulations for different hydrogen exchange models. The results clearly show that exchange from the native state of bovine alpha-lactalbumin at 15 degrees C is dominated by local unfolding events.     

Stability of secondary structural elements in a solvent-free environment. II: The β-pleated sheets

The stability of single β-strands and multistrand β-pleated sheets as elements of secondary structure is examined in the absence of intermolecular interactions. Such experimental conditions (e.g., complete removal of solvent molecules and counterions) are achieved by placing the peptide ions in the gas phase. The metastable multiply- charged peptide ions produced by electrospray ionization undergo unimolecular dissociation. Intercharge repulsion within the precursor ions gives rise to the elevated kinetic energy of fragment ions, which is measured using Mass-analyzed Ion Kinetic Energy (MIKE) spectrometry. Intercharge distances calculated based on these measurements are compared to the numbers derived from molecular mechanics calculations with charge site assignments based on relative proton affinities. Evidence is presented suggesting that single β-strands form collapsed structures in the absence of solvents, while multistrand β-pleated sheets are likely to retain “native-like” secondary structures under the same conditions. These results indicate that intramolecular hydrogen bonds are the major factor determining the three-dimensional arrangements of polypeptides in the gas phase, compensating both long- and short-range electrostatic repulsions. This is in good agreement with our earlier findings (Proteins 27:165–170, 1997) concerning stability of helical conformation of melittin in the absence of solvent. Proteins Suppl. 2:22–27, 1998. © 1998 Wiley-Liss, Inc.     

Characterization of disease-associated single amino acid polymorphisms in terms of sequence and structure properties.

In the present work, we use structural information to characterize a set of disease-associated single amino acid polymorphisms exhaustively. The analysis of different properties, such as substitution matrix elements, secondary structure, accessibility, free energies of transfer from water to octanol, amino acid volume, etc., suggests that many disease-causing mutations are associated with extreme changes in the value of parameters relating to protein stability. Overall, our results indicate that, while knowledge of protein structure clearly helps in understanding these mutations, a finer understanding can come only from a quantitative knowledge of protein stability and of the protein environment in the cell. Interestingly, use of evolutionary information from multiple sequence alignments can be used to increase our knowledge of disease-associated mutations.     

Role of methoxypolyethylene glycol on the hydration, activity, conformation and dynamic properties of a lipase in a dry film

A combined approach based on the use of ATR-FT/IR and steady-state fluorescence spectroscopy allowed to shed light on the effects of the additive methoxypolyethylene glycol (MePEG) on the hydration, conformation and dynamic properties of lipase from Burkholderia cepacia dehydrated to form a film. Spectroscopic data show that the additive has little effect on the structure of the protein; however, H/D exchange kinetic and fluorescence anisotropy suggest a more flexible enzyme molecule when in the presence of MePEG. By infrared spectroscopy, we estimated that, after conditioning the films at water activity of 1, the water content in the lipase dehydrated with MePEG is 5.4- and 4.7-fold higher than in the absence of the additive and the additive alone, respectively. Additionally, our infrared data suggest that MePEG acts by hindering intermolecular protein-protein interactions and contributing to increase the accessibility and flexibility of the lipase in the dehydrated solid film. These factors also explain the enhancement of the enzyme catalytic activity (i.e., up to 3.7-fold in neat organic solvent) when in the presence of MePEG. The method and results presented might better address the use of additives for the preparation of enzymes employed in non-aqueous media or of proteins used in a dry form in different fields of biotechnology.     

Sensitivity enhanced NMR spectroscopy by quenching scalar coupling mediated relaxation: Application to the direct observation of hydrogen bonds in 13C/15N-labeled proteins

In this paper, we demonstrate that the sensitivity of triple-resonance NMR experiments can be enhanced significantly through quenching scalar coupling mediated relaxation by using composite-pulse decoupling (CPD) or an adiabatic decoupling sequence on aliphatic, in particular alpha-carbons in ¹³C/¹⁵N-labeled proteins. The CPD-HNCO experiment renders 50% sensitivity enhancement over the conventional CT-HNCO experiment performed on a 12 kDa FK506 binding protein, when a total of 266 ms of amide nitrogen–carbonyl carbon defocusing and refocusing periods is employed. This is a typical time period for the direct detection of hydrogen bonds in proteins via trans-hydrogen bond ^3h J _NC couplings. The experimental data fit theoretical analysis well. The significant enhancement in sensitivity makes the experiment more applicable to larger-sized proteins without resorting to perdeuteration.     

Receptins: a novel term for an expanding spectrum of natural and engineered microbial proteins with binding properties for mammalian proteins

A new term ‘receptin’, derived from recipere (lat.), is proposed to denote microbial binding proteins that interact with mammalian target proteins. An example of such a ‘receptin’ is staphyloccocal protein A which binds to the Fc part of many mammalian immunoglobulins. Several other types of ‘receptins’ are listed. This term may easily be distinguished from the similar term ‘receptor’, describing a binding site on a cell surface, mostly eukaryotic, where a secondary effect is induced inside the cell upon binding to a ligand. A receptin, however, does not necessarily have to induce a secondary event. Receptins include so called MSCRAMMs, adhesins, and also engineered receptins, affibodies, and engineered ligands. It denotes any protein of microbial origin, cell‐bound or soluble, which can bind to a mammalian protein. It fulfills the need for an umbrella terminology for a large group of binding structures. In contrast, the term ‘lectin’ represents a group of proteins with affinity for carbohydrate structures. The new term ‘receptin’ includes a number of key microbial proteins involved in host–parasite interactions and in virulence. Some receptins are promising vaccine candidates. Copyright © 1999 John Wiley & Sons, Ltd.     

One- and two-dimensional NMR spectral analysis of the consequences of single amino acid replacements in proteins

The traditional approach of using homologous sequences to elucidate the role of specific amino acid residues in protein structure and function becomes more meaningful as the number of differences is minimized, with the limit being alteration of a single residue. For small proteins in solution, NMR spectroscopy offers a means of obtaining detailed information about each residue and its response to a given change in the protein sequence. Extraction of this information has been aided by recent progress in spectrometer technology (higher magnetic fields, more sensitive signal detection, more sophisticated computers) and experimental strategies (new NMR pulse sequences including multiple-quantum and two-dimensional NMR methods). The set of avian ovomucoid third domains, which consists of the third domain proper plus a short leader (connecting peptide) and has a maximum of 56 amino acid residues, offers an attractive system for developing experimental methods for investigating sequence-structure and structure-function relationships in proteins. Our NMR results provide examples of sequence effects on pKa' values, average conformation, and internal motion of amino acid side chains.     

The role of position a in determining the stability and oligomerization state of alpha-helical coiled coils: 20 amino acid stability coefficients in the hydrophobic core of proteins

We describe here a systematic investigation into the role of position a in the hydrophobic core of a model coiled-coil protein in determining coiled-coil stability and oligomerization state. We employed a model coiled coil that allowed the formation of an extended three-stranded trimeric oligomerization state for some of the analogs; however, due to the presence of a Cys-Gly-Gly linker, unfolding occurred from the same two-stranded monomeric oligomerization state for all of the analogs. Denaturation from a two-stranded state allowed us to measure the relative contribution of 20 different amino acid side chains to coiled-coil stability from chemical denaturation profiles. In addition, the relative hydrophobicity of the substituted amino acid side chains was assessed by reversed-phase high-performance liquid chromatography and found to correlate very highly (R = 0.95) with coiled-coil stability. We also determined the effect of position a in specifying the oligomerization state using ultracentrifugation as well as high-performance size-exclusion chromatography. We found that nine of the analogs populated one oligomerization state exclusively at peptide concentrations of 50 microM under benign buffer conditions. The Leu-, Tyr-, Gln-, and His-substituted analogs were found to be exclusively three-stranded trimers, while the Asn-, Lys-, Orn-, Arg-, and Trp-substituted analogs formed exclusively two-stranded monomers. Modeling results for the Leu-substituted analog showed that a three-stranded oligomerization state is preferred due to increased side-chain burial, while a two-stranded oligomerization state was observed for the Trp analog due to unfavorable cavity formation in the three-stranded state.