On the relative contribution of ionic interactions over iron-sulfur clusters to ferredoxin stability

Metal centres play an important structural role in maintaining the native conformation of a protein. Here we use biophysical methods to investigate what is the relative contribution of iron-sulfur clusters in respect to ionic interactions in a thermophilic di-cluster ferredoxin model. Changes in protonation affect both the stability and the conformational dynamics of the protein fold. In the pH 5.5-8 interval, the protein has a high melting temperature (T(m) approximately 120 degrees C), which decreases towards pH extremes. Acidification triggers events in two steps: down to the isoelectric point (pH 3.5) the Fe-S clusters remain unchanged, the secondary structure content increases and the single Trp becomes more solvent shielded, denoting a more compact fold. Further acidification down to pH 2 sets off exposure of the hydrophobic core and Fe-S cluster disintegration, yielding a molten globule state. The relative stabilising contribution of the clusters becomes evident when stabilising ionic interactions are switched off as a result of poising the protein at pH 3.5, at an overall null charge: under these conditions, the Fe-S clusters disassemble at T(m)=72 degrees C, whereas the protein unfolds at T(m)=52 degrees C. Overall, this ferredoxin denotes a considerable structural plasticity around its native conformation, a property which appears to depend more on the integrity of its metal clusters rather than on the status of its stabilising electrostatic interactions. The latter however play a relevant role in determining the protein thermal stability.     

The contribution of thymine-thymine interactions to the stability of folded dimeric quadruplexes.

The loop of four thymines in the sodium form of the dimeric folded quadruplex [d(G3T4G3)]2 assumes a well-defined structure in which hydrogen bonding between the thymine bases appears to contribute to the stability and final conformation of the quadruplex. We have investigated the importance of the loop interactions by systematically replacing each thymine in the loop with a cytosine. The quadruplexes formed by d(G3CT3G3), d(G3TCT2G3), d(G3T2CTG3) and d(G3T3CG3) in the presence of 150 mM Na+ were studied by gel mobility, circular dichroism and 1H NMR spectroscopy. The major species formed by d(G3CT3G3), d(G3TCT2G3) and d(G3T3CG3) at 1 mM strand concentration at neutral pH is a dimeric folded quadruplex. d(G3T2CTG3) has anomalous behaviour and associates into a greater percentage of linear four-stranded quadruplex than the other three oligonucleotides at neutral pH and at the same concentration. The linear four-stranded quadruplex has a greater tendency to oligomerize to larger ill-defined structures, as demonstrated by broad 1H NMR resonances. At pH 4, when the cytosine is protonated, there is a greater tendency for each of the oligonucleotides to form some four-stranded linear quadruplex, except for d(G3T2CTG3), which has the reverse tendency. The experimental results are discussed in terms of hydrogen bonding within the thymine loop.     

Estimating the contribution of engineered surface electrostatic interactions to protein stability by using double-mutant cycles

Coulombic interactions between charges on the surface of proteins contribute to stability. It is difficult, however, to estimate their importance by protein engineering methods because mutation of one residue in an ion pair alters the energetics of many interactions in addition to the coulombic energy between the two components. We have estimated the interaction energy between two charged residues, Asp-12 and Arg-16, in an alpha-helix on the surface of a barnase mutant by invoking a double-mutant cycle involving wild-type enzyme (Asp-12, Thr-16), the single mutants Thr----Arg-16 and Asp----Ala-12, and the double mutant Asp----Ala-12, Thr----Arg-16. The changes in free energy of unfolding of the single mutants are not additive because of the coulombic interaction energy. Additivity is restored at high concentrations of salt that shield electrostatic interactions. The geometry of the ion pair in the mutant was assumed to be the same as that in the highly homologous ribonuclease from Bacillus intermedius, binase, which has Asp-12 and Arg-16 in the native enzyme. The ion pair does not form a hydrogen-bonded salt bridge, but the charges are separated by 5-6 A. The mutant barnase containing the ion pair Asp-12/Arg-16 is more stable than wild type by 0.5 kcal/mol, but only a part of the increased stability is attributable to the electrostatic interaction. We present a formal analysis of how double-mutant cycles can be used to measure the energetics of pairwise interactions.     

The contribution of ionic interactions to the conformational stability and function of polygalacturonase from A. niger

Aspergillus niger produces multiple forms of polygalacturonases with molecular masses ranging from 30 to 60 kDa. The high molecular weight polygalacturonase (61 ± 2 kDa) from A. niger possesses a pH optimum of 4.3 and a pI of 3.9. The enzyme exhibited high sensitivity, both in terms of activity and structure, in the pH range of 4.3–7.0. The enzyme was irreversibly inactivated at pH 7.0. The enzyme is predominantly rich in parallel β structure. There is unfolding of the enzyme molecule between 4.3 and 7.0 resulting in irreversible loss of secondary and tertiary structure with the exposure of hydrophobic surfaces. ANS binding measurements, intrinsic fluorescence and acrylamide quenching measurements have confirmed the unfolding and exposure of hydrophobic surfaces. The midpoint of pH transition for both activity and secondary structure is 6.2 ± 0.1. The pH-induced changes of polygalacturonase confirm the role of histidine residues in structure and activity of the enzyme. The irreversible nature of inactivation is due to the unfolding induced exposure of hydrophobic surfaces leading to association/aggregation of the molecule. Size exclusion chromatography measurements have established the association of enzyme at higher pH. Urea induced unfolding measurements at pH 4.3 and 7.0 have confirmed the loss in stability as we approach neutral pH.     

Histidine-aromatic interactions in barnase. Elevation of histidine pKa and contribution to protein stability.

Aromatic side-chains are found in the vicinity of histidine residues in many proteins and protein complexes. We have studied the interaction between a histidine residue (His18) and aromatic residues at position 94 in barnase. Three different techniques have been applied to show that Trp94 interacts more strongly with the protonated form of His18. The aromatic-histidine interaction stabilizes the protonated form of histidine by 0.8 to 1 kcal mol-1 relative to the unprotonated and, thereby, increases its pKa value. This was shown indirectly from the pH dependence of the stability of the wild-type protein and the mutant Trp94----Leu; and directly from the difference in pKa of His18 between wild-type barnase and the same mutant protein, and from double-mutant cycles that measure the total interaction energy of Trp94 with His18 at both low and high pH. When Trp94 is replaced by other aromatic amino acids, the strength of the interaction decreases in the series His-Trp greater than His-Tyr greater than His-Phe. The interaction is not masked by high salt concentrations. The raising of the pKa value of His18 by interaction with Trp94 is shown to be consistent with solution studies with model compounds. The histidine-aromatic interaction could have implications in binding and catalysis for modulation of the histidine pKa value.     

Electrostatic contribution to the binding stability of protein-protein complexes

To investigate roles of electrostatic interactions in protein binding stability, electrostatic calculations were carried out on a set of 64 mutations over six protein-protein complexes. These mutations alter polar interactions across the interface and were selected for putative dominance of electrostatic contributions to the binding stability. Three protocols of implementing the Poisson-Boltzmann model were tested. In vdW4 the dielectric boundary between the protein low dielectric and the solvent high dielectric is defined as the protein van der Waals surface and the protein dielectric constant is set to 4. In SE4 and SE20, the dielectric boundary is defined as the surface of the protein interior inaccessible to a 1.4-A solvent probe, and the protein dielectric constant is set to 4 and 20, respectively. In line with earlier studies on the barnase-barstar complex, the vdW4 results on the large set of mutations showed the closest agreement with experimental data. The agreement between vdW4 and experiment supports the contention of dominant electrostatic contributions for the mutations, but their differences also suggest van der Waals and hydrophobic contributions. The results presented here will serve as a guide for future refinement in electrostatic calculation and inclusion of nonelectrostatic effects.     

The electrostatic contribution to DNA base-stacking interactions.

Base-stacking and phosphate-phosphate interactions in B-DNA are studied using the finite difference Poisson-Boltzmann equation. Interaction energies and dielectric constants are calculated and compared to the predictions of simple dielectric models. No extant simple dielectric model adequately describes phosphate-phosphate interactions. Electrostatic effects contribute negligibly to the sequence and conformational dependence of base-stacking interactions. Electrostatic base-stacking interactions can be adequately modeled using the Hingerty screening function. The repulsive and dispersive Lennard-Jones interactions dominate the dependence of the stacking interactions on roll, tilt, twist, and propellor. The Lennard-Jones stacking energy in ideal B-DNA is found to be essentially independent of sequence.     

The folding of an enzyme. II. Substructure of barnase and the contribution of different interactions to protein stability.

Barnase is described anatomically in terms of its substructures and their mode of packing. The surface area of hydrophobic residues buried on formation and packing of the structural elements has been calculated. Changes in stability have been measured for 64 mutations, 41 constructed in this study, strategically located over the protein. The purpose is to provide: (1) information on the magnitudes of changes in stabilization energy for mutations of residues that are important in maintaining the structure; and (2) probes for the folding pathway to be used in subsequent studies. The majority of mutations delete functional moieties of side-chains or make isosteric changes. The energetics of the interactions are variable and context-dependent. The following general conclusions may be drawn, however, from this study about the classes of interactions that stabilize the protein. (1) Truncation of buried hydrophobic side-chains has, in general, the greatest effect on stability. For fully buried residues, this averages at 1.5 kcal mol-1 per methylene group with a standard deviation of +/- 0.6 kcal mol-1. Truncation of partly exposed leucine, isoleucine or valine residues that are in the range of 50 to 80 A2 of solvent-accessible area (30 to 50% of the total solvent-accessible area on a Gly-X-Gly tripeptide, i.e. those packed against the surface) has a smaller, but relatively constant effect on stability, at 0.81 kcal mol-1 per methylene group with a statistical standard deviation of +/- 0.18 kcal mol-1. (2) There is a very poor correlation between hydrophobic surface area buried and the free energy change for an extensive data set of hydrophobic mutants. The best correlation is found to be between the free energy change and the number of methylene groups within a 6 A radius of the hydrophobic groups deleted. (3) Burial of the hydroxyl group of threonine in a pocket that is intended for a gamma-methyl group of valine costs 2.5 kcal mol-1, in the range expected for the loss of two hydrogen bonds.(ABSTRACT TRUNCATED AT 400 WORDS)     

The solvent contribution to the free energy of protein--ligand interactions.

The intrinsic solvent contribution to the free energy of protein-ligand interactions in solution is shown to be related to a free energy per unit area term, obtained from analysis of the solution to gas phase process, and the change in accessible area on association. Analysis of the free energy data on a per unit area basis for the solution to gas phase process leads to the conclusion that the aliphatic CH₂ group is only slightly intrinsically hydrophobic, $\text{δΔG°/}\text{A}\text{?}{}^2\text{= 6}\text{cal mol}{}^{\mathrm{?1}}\text{A}\text{?}{}^2$, whereas the aromatic compound are actually intrinsically hydrophilic, $\text{δΔG°/}\text{A}\text{?}{}^2\text{= -26}\text{cal mol}{}^{\mathrm{?1}}\text{A}\text{?}{}^2$. This leads to the conclusion that, for the interaction of benzene, naphthalene and anthracene with the binding site of α-chymotrypsin, the ligand-solvent free energy contribution is actually unfavorable. Since the protein-solvent contribution is small or unfavorable, the central conclusion is that the solvent contribution to protein-ligand interactions is small or unfavorable and that it is the protein-ligand non-bonded interactions that provide the driving force for association.     

The contribution of halogen atoms to protein-ligand interactions.

The three-dimensional structure of para-fluoro-D-phenylalanine (PFF) in its complex with the zinc protease carboxypeptidase A (CPA) has been determined at 2.0 A resolution by X-ray crystallographic methods. The structure reveals that the para-fluorobenzyl side chain of the inhibitor is buried in the S'1 hydrophobic pocket of the enzyme. Intriguingly, this ligand molecule inhibits CPA better than its amino acid analogues D-phenylalanine (D-Phe) and D-tyrosine (D-Tyr) by factors of 4 and 5, respectively. Moreover, the para-fluoro derivative is a better inhibitor than para-chloro- or para-bromo-D-phenylalanine by nearly a factor of 50. This result is consistent with binding enhancements realized in other protein complexes involving halogenated ligand molecules, regardless of whether the carbon-halogen group of the ligand makes specific polar interactions or non-specific hydrophobic interactions with its protein host. In the CPA-PFF complex, the fluorine atom of PFF does not make any direct polar contact with the enzyme, and the contact surface area of the protein-ligand interface is only slightly greater, although more hydrophobic, than that of D-Phe and D-Tyr. Therefore, we conclude that the slight binding enhancement measured for PFF relative to D-Phe and D-Tyr arises predominantly from increasing the hydrophobic character of the protein-ligand interface, and not solely from increasing the degree of protein-ligand contact.     

Role of invariant water molecules and water-mediated ionic interactions in D-xylose isomerase from Streptomyces rubiginosus

The enzyme, D-xylose isomerase (D-xylose keto-isomerase; EC 5.3.1.5) is a soluble enzyme that catalyzes the conversion of the aldo-sugar D-xylose to the keto-sugar D-xylulose. A total of 27 subunits of D-xylose isomerase from Streptomyces rubiginosus were analyzed in order to identify the invariant water molecules and their water-mediated ionic interactions. A total of 70 water molecules were found to be invariant. The structural and/or functional roles of these water molecules have been discussed. These invariant water molecules and their ionic interactions may be involved in maintaining the structural stability of the enzyme D-xylose isomerase. Fifty-eight of the 70 invariant water molecules (83%) have at least one interaction with the main chain polar atom.     

Contribution of cation-pi interactions to protein stability

Calculations predict that cation- interactions make an important contribution to protein stability. While there have been some attempts to experimentally measure strengths of cation-pi interactions using peptide model systems, much less experimental data are available for globular proteins. We have attempted to determine the magnitude of cation-pi interactions of Lys with aromatic amino acids in four different proteins (LIVBP, MBP, RBP, and Trx). In each case, Lys was replaced with Gln and Met. In a separate series of experiments, the aromatic amino acid in each cation-pi pair was replaced by Leu. Stabilities of wild-type (WT) and mutant proteins were characterized by both thermal and chemical denaturation. Gln and aromatic --> Leu mutants were consistently less stable than corresponding Met mutants, reflecting the nonisosteric nature of these substitutions. The strength of the cation-pi interaction was assessed by the value of the change in the free energy of unfolding [DeltaDeltaG(degrees) = DeltaG(degrees)(Met) - DeltaG(degrees)(WT)]. This ranged from +1.1 to -1.9 kcal/mol (average value -0.4 kcal/mol) at 298 K and +0.7 to -2.6 kcal/mol (average value -1.1 kcal/mol) at the Tm of each WT. It therefore appears that the strength of cation-pi interactions increases with temperature. In addition, the experimentally measured values are appreciably smaller in magnitude than calculated values with an average difference /DeltaG(degrees)expt - DeltaG(degrees)calc/av of 2.9 kcal/mol. At room temperature, the data indicate that cation-pi interactions are at best weakly stabilizing and in some cases are clearly destabilizing. However, at elevated temperatures, close to typical Tm's, cation-pi interactions are generally stabilizing.     

Contribution of hydrophobic interactions to protein stability

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6 ± 0.3 kcal/mol per -CH₂- group), than to the stability of a large protein, VlsE (1.6 ± 0.3 kcal/mol per -CH₂- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH₂- group on folding contributes, on average, 1.1 ± 0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG_tr values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60 ± 4% and hydrogen bonds contribute 40 ± 4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.     

Salt contribution to RNA tertiary structure folding stability

Accurate quantification of the ionic contribution to RNA folding stability could greatly enhance our ability to understand and predict RNA functions. Recently, motivated by the potential importance of ion correlation and fluctuation in RNA folding, we developed the tightly bound ion (TBI) model. Extensive experimental tests showed that the TBI model can lead to better treatment of multivalent ions than the Poisson-Boltzmann equation. In this study, we use the model to quantify the contribution of salt (Na⁺ and Mg²⁺) to the RNA tertiary structure folding free energy. Folding of the RNA tertiary structure often involves intermediates. We focus on the folding transition from an intermediate state to the native state, and compute the electrostatic folding free energy of the RNA. Based on systematic calculations for a variety of RNA molecules, we derive a set of formulas for the electrostatic free energy for tertiary structural folding as a function of the sequence length and compactness of the RNA and the Na⁺ and Mg²⁺ concentrations. Extensive comparisons with experimental data suggest that our model and the extracted empirical formulas are quite reliable.     

Environmental bacteriophage-host interactions: factors contribution to natural transduction

Over the past two decades the potential for the exchange of bacterial genes in natural environments through transduction (bacteriophage-mediated gene transfer) has been well established. Studies carried out by various laboratories throughout the world have demonstrated that both chromosomal and plasmid DNA can be successfully transduced in natural environments ranging from sewer plants to rivers and lakes. Transduction has been shown to take place in the gills of oysters and the kidneys of mice. Model studies have demonstrated the ability of transduction to maintain genetic material in bacterial gene pools that would otherwise be lost because of negative fitness. Thus, transduction may affect the course of bacterial evolution. Identification of natural transduction has led to the investigation of the dynamics of bacteriophage host interactions in natural aquatic environments and to the exploration of various environmental factors that affect virus-host interactions. Two important environmental factors which affect virus-host interactions are the metabolic state of the host and the exposure of the host to DNA-damaging stresses such as solar UV light. Recent researches on these two areas of virus-host relationships are reviewed.     

Dissecting the contribution of diffusion and interactions to the mobility of nuclear proteins

Quantitative characterization of protein interactions under physiological conditions is vital for systems biology. Fluorescence photobleaching/activation experiments of GFP-tagged proteins are frequently used for this purpose, but robust analysis methods to extract physicochemical parameters from such data are lacking. Here, we implemented a reaction-diffusion model to determine the contributions of protein interaction and diffusion on fluorescence redistribution. The model was validated and applied to five chromatin-interacting proteins probed by photoactivation in living cells. We found that very transient interactions are common for chromatin proteins. Their observed mobility was limited by the amount of free protein available for diffusion but not by the short residence time of the bound proteins. Individual proteins thus locally scan chromatin for binding sites, rather than diffusing globally before rebinding at random nuclear positions. By taking the real cellular geometry and the inhomogeneous distribution of binding sites into account, our model provides a general framework to analyze the mobility of fluorescently tagged factors. Furthermore, it defines the experimental limitations of fluorescence perturbation experiments and highlights the need for complementary methods to measure transient biochemical interactions in living cells.     

Role of cation-pi interactions to the stability of thermophilic proteins

Elucidating the factors responsible for exhibiting extreme thermal stability of thermophilic proteins is very important for an understanding of the mechanism of protein stability, as well as to design stable proteins. In this work, we have analyzed the influence of cation-pi interactions to enhance the stability from mesophilic to thermophilic proteins. The favorable residue pairs forming such a system of interactions have been brought out. We found that the Tyr has a greater number of such interactions with Lys in thermophilic proteins. Specifically, the same Lys would experience a greater number of cation-pi interactions with several Tyr residues in thermophiles. On the other hand, the influence of Phe in making cation-pi interactions is higher in mesophiles than in thermophiles. Further, a network of cation-pi interactions are maintained by Lys in thermophiles, whereas Arg plays a major role in mesophilic proteins. Moreover, atoms that have a substantial positive charge in both Lys and Arg make a more significant contribution for cation-pi interactions than do cationic group atoms.