Long-range electrostatic interactions can influence the folding, stability, and cooperativity of dihydrofolate reductase

DNA Framework-Engineered Long-Range Electrostatic Interactions for DNA Hybridization Reactions

Long-range electrostatic interactions beyond biomolecular interaction interfaces have not been extensively studied due to the limitation in engineering electric double layers in physiological fluids. Here we find that long-range electrostatic interactions play an essential role in kinetic modulation of DNA hybridizations. Protein and gold nanoparticles with different charges are encapsulated in tetrahedral frameworks to exert diverse electrostatic effects on site-specifically tethered single DNA strands. Using this strategy, we have successfully modulated the hybridization kinetics in both bulk solution and single molecule level. Experimental and theoretical studies reveal that long-range Coulomb interactions are the key factor for hybridization rates. This work validates the important role of long-range electrostatic forces in nucleic acid-biomacromolecule complexes, which may encourage new strategies of gene regulation, antisense therapy, and nucleic acid detection.

One Plus One Makes Three: Triangular Coupling of Correlated Amino Acid Mutations

Correlated mutations have played a pivotal role in the recent success in protein fold prediction. Understanding nonadditive effects of mutations is crucial for altering protein structure, as mutations of multiple residues may change protein stability or binding affinity in a manner unforeseen by the investigation of single mutants. While the couplings between amino acids can be inferred from homologous protein sequences, the physical mechanisms underlying these correlations remain elusive. In this work we demonstrate that calculations based on the first-principles of statistical mechanics are capable of capturing the effects of nonadditivities in protein mutations. The identified thermodynamic couplings cover the short-range as well as previously unknown long-range correlations. We further explore a set of mutations in staphyloccocal nuclease to unravel an intricate interaction pathway underlying the correlations between amino acid mutations.

The expression and functional analysis of the sialyl-T antigen in prostate cancer

Aberrant glycosylation is a featured characteristic of cancer and plays a role in cancer pathology; thus an understanding of the compositions and functions of glycans is critical for discovering diagnostic biomarkers and therapeutic targets for cancer. In this study, we used MALDI-TOF-MS analysis to determine the O-glycan profiles of prostate cancer cells metastasized to bone (PC-3), brain (DU145), lymph node (LNCaP), and vertebra (VCaP) in comparison to immortalized RWPE-1 cells derived from normal prostatic tissue. Prostate cancer (CaP) cells exhibited an elevation of simple/short O-glycans, with a reduction of complex O-glycans, increased O-glycan sialylation and decreased fucosylation. Core 1 sialylation was increased dramatically in all CaP cells, and especially in PC-3 cells. The expression of Neu5Acα2-3Galβ1-3GalNAc- (sialyl-3T antigen) which is the product of α2,3-sialyltransferase-I (ST3Gal-I) was substantially increased. We therefore focused on exploring the possible function of ST3Gal-I in PC-3 cells. ST3Gal-I silencing studies showed that ST3Gal-I was associated with PC-3 cell proliferation, migration and apoptosis. Further in vivo studies demonstrated that down regulation of ST3Gal-I reduced the tumor size in xenograft mouse model, indicating that sialyl-3T can serve as a biomarker for metastatic prostate cancer prognosis, and that ST3Gal-I could be a target for therapeutic intervention in cancer treatment.

How a multimeric macromolecule is affected by divalent salts? Experimental and simulation study

Salts exist in any cell and living organism in contact with biological macromolecules. How these salts affect biomolecules such as enzyme is important from both basic sciences and practical technologies. It was observed that divalent salts can change structure and function of protein at higher concentrations. Here, we investigated the effect of divalent salt on the behavior of a multimeric enzyme. We treated glucose oxidase as dimer-active enzyme in different CaCl₂ concentration and seen that the enzyme become inactive at high concentration of salt. These experimental results are in agreement with recently published researches. To find a possible mechanism, a series of molecular dynamics simulation of the enzyme were performed at different salt concentration. According to the MD simulation, the conformational changes at the active site and FAD-binding site support the hypothesis of enzyme inactivation at high CaCl₂ concentration. MD simulations also showed that enzyme has an unstable conformation at higher salt concentration which is in agreement with our experimental data. Detailed structural properties of the enzyme have been analyzed under different conditions. To the best of our knowledge, this is the first study that bears detailed structural mechanism about the salt effects on multimeric macromolecules.

Quality matters: extension of clusters of residues with good hydrophobic contacts stabilize (hyper)thermophilic proteins

Identifying determinant(s) of protein thermostability is key for rational and data-driven protein engineering. By analyzing more than 130 pairs of mesophilic/(hyper)thermophilic proteins, we identified the quality (residue-wise energy) of hydrophobic interactions as a key factor for protein thermostability. This distinguishes our study from previous ones that investigated predominantly structural determinants. Considering this key factor, we successfully discriminated between pairs of mesophilic/(hyper)thermophilic proteins (discrimination accuracy: ∼80%) and searched for structural weak spots in E. coli dihydrofolate reductase (classification accuracy: 70%).

From systems to structure: bridging networks and mechanism

There is a wide gap between the generation of large-scale biological data sets and more-detailed, structural and mechanistic studies. However, recent studies that explicitly combine data from systems and structural biological approaches are having a profound effect on our ability to predict how mutations and small molecules affect atomic-level mechanisms, disrupt systems-level networks, and ultimately lead to changes in organismal fitness. In fact, we argue that a shared framework for analysis of nonadditive genetic and thermodynamic responses to perturbations will accelerate the integration of reductionist and global approaches. A stronger bridge between these two areas will allow for a deeper and more-complete understanding of complex biological phenomenon and ultimately provide needed breakthroughs in biomedical research.

Urea denatured state ensembles contain extensive secondary structure that is increased in hydrophobic proteins

The goal of this article is to gain a better understanding of the denatured state ensemble (DSE) of proteins through an experimental and computational study of their denaturation by urea. Proteins unfold to different extents in urea and the most hydrophobic proteins have the most compact DSE and contain almost as much secondary structure as folded proteins. Proteins that unfold to the greatest extent near pH 7 still contain substantial amounts of secondary structure. At low pH, the DSE expands due to charge-charge interactions and when the net charge per residue is high, most of the secondary structure is disrupted. The proteins in the DSE appear to contain substantial amounts of polyproline II conformation at high urea concentrations. In all cases considered, including staph nuclease, the extent of unfolding by urea can be accounted for using the data and approach developed in the laboratory of Wayne Bolen (Auton et al., Proc Natl Acad Sci 2007; 104:15317-15323).

Interaction of the aromatics Tyr-72/Trp-288 in the interface of the extracellular and transmembrane domains is essential for proton gating of acid-sensing ion channels

Acid-sensing ion channels are proton-activated ion channels expressed in the nervous system. They belong to the family of ENaC/Degenerins whose members share a conserved structure but are activated by widely diverse stimuli. We show that interaction of two aromatic residues, Tyr-72, located immediately after the first transmembrane segment, and Trp-288, located at the tip of a loop of the extracellular domain directed toward the first transmembrane segment, is essential for proton activation of the acid-sensing ion channels. The subdomain containing Trp-288 is a module tethered to the rest of the extracellular domain by short linkers and intrasubunit interactions between residues in the putative "proton sensor." Mutations in these two areas shift the apparent affinity of protons toward a more acidic range and change the kinetics of activation and desensitization. These results are consisting with displacement of the module relative to the rest of the extracellular domain to allow interaction of Trp-288 with Tyr-72 during gating. We propose that such interaction may provide functional coupling between the extracellular domain and the pore domain.

Evolutional Design of a Hyperactive Cysteine- and Methionine-free Mutant of Escherichia coli Dihydrofolate Reductase

We developed a strategy for finding out the adapted variants of enzymes, and we applied it to an enzyme, dihydrofolate reductase (DHFR), in terms of its catalytic activity so that we successfully obtained several hyperactive cysteine- and methionine-free variants of DHFR in which all five methionyl and two cysteinyl residues were replaced by other amino acid residues. Among them, a variant (M1A/M16N/M20L/M42Y/C85A/M92F/C152S), named as ANLYF, has an approximately seven times higher k(cat) value than wild type DHFR. Enzyme kinetics and crystal structures of the variant were investigated for elucidating the mechanism of the hyperactivity. Steady-state and transient binding kinetics of the variant indicated that the kinetic scheme of the catalytic cycle of ANLYF was essentially the same as that of wild type, showing that the hyperactivity was brought about by an increase of the dissociation rate constants of tetrahydrofolate from the enzyme-NADPH-tetrahydrofolate ternary complex. The crystal structure of the variant, solved and refined to an R factor of 0.205 at 1.9-angstroms resolution, indicated that an increased structural flexibility of the variant and an increased size of the N-(p-aminobenzoyl)-L-glutamate binding cleft induced the increase of the dissociation constant. This was consistent with a large compressibility (volume fluctuation) of the variant. A comparison of folding kinetics between wild type and the variant showed that the folding of these two enzymes was similar to each other, suggesting that the activity enhancement of the enzyme can be attained without drastic changes of the folding mechanism.

Local stability identification and the role of key acidic amino acid residues in staphylococcal nuclease unfolding

Staphylococcal nuclease is a single domain protein with 149 amino acids. It has no disulfide bonds, which makes it a simple model for the study of protein folding. In this study, 20 mutants of this protein were generated each with a single base substitution of glycine for negatively charged glutamic acid or aspartic acid. Using differential scanning microcalorimetry in thermal denaturation experiments, we identified two mutants, E75G and E129G, having approximately 43% and 44%, respectively, lower DeltaH(cal) values than the wild-type protein. Furthermore, two mutants, E75Q and E129Q, were created and the results imply that substitution of the Gly residue has little influence on destabilization of the secondary structure that leads to the large perturbation of the tertiary protein structure stability. Two local stable areas formed by the charge-charge interactions around E75 and E129 with particular positively charged amino acids are thus identified as being significant in maintenance of the three-dimensional structure of the protein.

Complex interactions at the helix-helix interface stabilize the glycophorin A transmembrane dimer

To explore the residue interactions in the glycophorin A dimerization motif, an alanine scan double mutant analysis at the helix-helix interface was carried out. These data reveal a combination of additive and coupled effects. The majority of the double mutants are found to be equally or slightly more stable than would be predicted by the sum of the energetic cost of the single-point mutants. The proximity of the mutated sites is not related to the presence of coupling between those sites. Previous studies reveal that a single face of the glycophorin A monomer contains a specific glycine-containing motif (GxxxG) that is thought to be a driving force for the association of transmembrane helices. Double mutant cycles suggest that the relationship of the GxxxG motif to the remainder of the helix-helix interface is complex. Sequences containing mutations that abolish the GxxxG motif retain an ability to dimerize, while a sequence containing a GxxxG motif appears unable to form dimers. The energetic effects of weakly coupled and additive double mutants can be explained by changes in van der Waals interactions at the dimer interface. These results emphasize the fact that the sequence context of the dimer interface modulates the strength of the glycophorin A GxxxG-mediated transmembrane dimerization reaction.

Conferring thermostability to mesophilic proteins through optimized electrostatic surfaces

Recently, there have been several experimental reports of proteins displaying appreciable stability gains through mutation of one or two amino acid residues. Here, we employ a simple theoretical model to quickly screen mutant structures for increased thermostability through optimization of the protein's electrostatic surface. Our results are able to reproduce the experimental observation that elimination of like-charge repulsions and creation of opposite-charge attractions on the protein surface is an efficient method to confer thermostability to a mesophilic protein. Using Poisson-Boltzmann electrostatics, we calculate relative protein stabilities for the exhaustive surface mutagenesis of the cold shock, RNase T1, and CheY proteins. Comparison with 25 experimentally characterized cold shock protein mutants reveals an average correlation of 0.86. The model is also quantitatively accurate when reproducing the experimental D49A and D49H mutant stabilities of RNase T1. This work represents the first comprehensive in silico screening of mutant candidates likely to confer thermostability to mesophilic proteins through optimization of surface electrostatics. Systematic single mutant, followed by double mutant, screening yields a limited number of mutant structures displaying significant stability gains suitable for subsequent experimental characterization.

Thermodynamic analysis of the contributions of the copper ion and the disulfide bridge to azurin stability: synergism among multiple depletions

The stabilizing potential of the copper ion and the disulfide bridge in azurin has been explored with the aim of inspecting the ways in which these two factors influence one another. Specifically, whether copper and disulfide contributions to protein stability are additive has been examined. To this aim, the thermal unfolding of a copper-depleted mutant lacking the disulfide bridge between Cys3 and Cys26 (apo C3A/C26A azurin) was studied by differential scanning calorimetry. A comparison of the unfolding parameters of holo and apo C3A/C26A azurin with the apo C3A/C26A protein has shown that the effects of simultaneous copper and disulfide depletion are additive only at two temperatures: T=15 degrees C and T=67 degrees C. Within this range the presence of the copper ion and the disulfide bridge has a positive synergistic effect on azurin stability. These findings might have implications for the rational use of the stabilizing potential of copper and disulfides in copper protein engineering.

The effect of salts on the activity and stability of Escherichia coli and Haloferax volcanii dihydrofolate reductases

The extremely halophilic Archae require near-saturating concentrations of salt in the external environment and in their cytoplasm, potassium being the predominant intracellular cation. The proteins of these organisms have evolved to function in concentrations of salt that inactivate or precipitate homologous proteins from non-halophilic species. It has been proposed that haloadaptation is primarily due to clustering of acidic residues on the surface of the protein, and that these clusters bind networks of hydrated ions. The dihydrofolate reductases from Escherichia coli (ecDHFR) and two DHFR isozymes from Haloferax volcanii (hvDHFR1 and hvDHFR2) have been used as a model system to compare the effect of salts on a mesophilic and halophilic enzyme. The KCl-dependence of the activity and substrate affinity was investigated. ecDHFR is largely inactivated above 1M KCl, with no major effect on substrate affinity. hvDHFR1 and hvDHFR2 unfold at KCl concentrations below approximately 0.5M. Above approximately 1M, the KCl dependence of the hvDHFR activities can be attributed to the effect of salt on substrate affinity. The abilities of NaCl, KCl, and CsCl to enhance the stability to urea denaturation were determined, and similar efficacies of stabilization were observed for all three DHFR variants. The DeltaG degrees (H(2)O) values increased linearly with increasing KCl and CsCl concentrations. The increase of DeltaG degrees (H(2)O) as a function of the smallest cation, NaCl, is slightly curved, suggesting a minor stabilization from cation binding or screening of electrostatic repulsion. At their respective physiological ionic strengths, the DHFR variants exhibit similar stabilities. Salts stabilize ecDHFR and the hvDHFRs by a common mechanism, not a halophile-specific mechanism, such as the binding of hydrated salt networks. The primary mode of salt stabilization of the mesophilic and halophilic DHFRs appears to be through preferential hydration and the Hofmeister effect of salt on the activity and entropy of the aqueous solvent. In support of this conclusion, all three DHFRs are similarly stabilized by the non-ionic cosolute, sucrose.

Effect of ligand binding on the flexibility of dihydrofolate reductase as revealed by compressibility

The partial specific volume, v, and adiabatic compressibility, beta(s), of Escherichia coli dihydrofolate reductase were measured at 30 degrees C in the presence of various ligands (folate, dihydrofolate, tetrahydrofolate, NADPH, NADP, methotrexate, and KCl). Binding of these ligands (binary and ternary complexes) brought about large changes of v (0.734-0.754 cm(3) g(-1)) and beta(s) (6. 6x10(-6)-9.8x10(-6) bar(-1)), keeping a linear relationship between the two parameters. The values of v and beta(s) increased with an increase in internal cavity, V(cav), and a decrease in accessible surface area, ASA, which were calculated from the X-ray crystal structures of the complexes. A large variation of V(cav) relative to ASA by ligand binding suggested that the cavity is a dominant factor and the effect of hydration might be small for the ligand-induced changes of v and beta(s). The beta(s) values of the binary and ternary complexes suggested a characteristic conformational flexibility of the kinetic intermediates in the enzyme reaction coordinate. Comparison of beta(s) with the cavity distribution in the crystal structures revealed that the flexibility of the intermediates was mainly determined by the total cavity volume with minor contributions of the number, position, and size of cavities. These results demonstrate that the compressibility is a useful measure of the conformational flexibility of the intermediates in the enzyme reaction and that the combined study of compressibility and X-ray crystallography gives new insight into the protein dynamics through the behavior of the cavities.

Asymmetric and independent contribution of the second transmembrane segment 12' residues to diliganded gating of acetylcholine receptor channels: a single-channel study with choline as the agonist

Mutagenesis studies have suggested that the second transmembrane segment (M2) plays a critical role during acetylcholine receptor liganded gating. An adequate description of the relationship between gating and structure of the M2 domain, however, has been hampered by the fact that many M2 mutations increase the opening rate constant to levels that, in the presence of acetylcholine, are unresolvably fast. Here, we show that the use of saturating concentrations of choline, a low-efficacy agonist, is a convenient tool to circumvent this problem. In the presence of 20 mM choline: (a) single-channel currents occur in clusters; (b) fast blockade by choline itself reduces the single-channel conductance by approximately 50%, yet the excess open-channel noise is only moderate; (c) the kinetics of gating are fitted best by a single-step, C <--> O model; and (d) opening and closing rate constants are within a well resolvable range. Application of this method to a series of recombinant adult mouse muscle M2 12' mutants revealed that: (a) the five homologous M2 12' positions make independent and asymmetric contributions to diliganded gating, the delta subunit being the most sensitive to mutation; (b) mutations at delta12' increase the diliganded gating equilibrium constant in a manner that is consistent with the sensitivity of the transition state to mutation being approximately 30% like that of the open state and approximately 70% like that of the closed state; (c) the relationship between delta12' amino acid residue volume, hydrophobicity or alpha-helical tendency, and the gating equilibrium constant of the corresponding mutants is not straightforward; however, (d) rate and equilibrium constants for the mutant series are linearly correlated (on log-log plots), which suggests that the conformational rearrangements upon mutation are mostly local and that the position of the transition state along the gating reaction coordinate is unaffected by these mutations.

Evolutionarily conserved pathways of energetic connectivity in protein families

For mapping energetic interactions in proteins, a technique was developed that uses evolutionary data for a protein family to measure statistical interactions between amino acid positions. For the PDZ domain family, this analysis predicted a set of energetically coupled positions for a binding site residue that includes unexpected long-range interactions. Mutational studies confirm these predictions, demonstrating that the statistical energy function is a good indicator of thermodynamic coupling in proteins. Sets of interacting residues form connected pathways through the protein fold that may be the basis for efficient energy conduction within proteins.

Mutational analysis of the BPTI folding pathway: II. Effects of aromatic-->leucine substitutions on folding kinetics and thermodynamics

The rates of the individual steps in the disulfide-coupled folding and unfolding of eight BPTI variants, each containing a single aromatic to leucine amino acid replacement, were measured. From this analysis, the contributions of the four phenylalanine and four tyrosine residues to the stabilities of the native protein and the disulfide-bonded folding intermediates were determined. While the substitutions were found to destabilize the native protein by 2 to 7 kcal/mol, they had significantly smaller effects on the intermediates that represent the earlier stages of folding, even when the site of the substitution was located within the ordered regions of the intermediates. These results suggest that stabilizing interactions contribute less to conformational stability in the context of a partially folded intermediate than in a fully folded native protein, perhaps because of decreased cooperativity among the individual interactions. The kinetic analysis also provides new information about the transition states associated with the slowest steps in folding and unfolding, supporting previous suggestions that these transition states are extensively unfolded. Although the substitutions caused large changes in the distribution of folding intermediates and in the rates of some steps in the folding pathway, the kinetically-preferred pathway for all of the variants involved intramolecular disulfide rearrangements, as observed previously for the wild-type protein. These results suggest that the predominance of the rearrangement mechanism reflects conformational constraints present relatively early in the folding pathway.

Circularly permuted dihydrofolate reductase possesses all the properties of the molten globule state, but can resume functional tertiary structure by interaction with its ligands

It is obvious that functional activity of a protein molecule is closely related to its structure. On the other hand, the understanding of structure-function relationship still remains one of the intriguing problems of molecular biology. There is widespread belief that mutagenesis presents a real way to solve this problem. Following this assumption, we have investigated the effect of circular permutation in dihydrofolate reductase from E. coli on protein structure and functioning. It has been shown that in the absence of ligands two circularly permuted variants of dihydrofolate reductase possess all the properties of the molten globule state. However, after addition of ligands they gain the native-like structural properties and specific activity. This means that the in vitro folding of permuted dihydrofolate reductase is terminated at the stage of the molten globule formation. Interaction of permuted protein with ligands leads to the structural adjustment and formation of active protein molecules.

Probing the determinants of protein stability: comparison of class A beta-lactamases

Five class A beta-lactamases produced by various mesophilic bacterial species have been compared. Although closely related in primary and overall structures, these enzymes exhibit very different stabilities. In order to investigate the factors responsible for these differences, several features deduced from the amino acid composition and three-dimensional structures were studied for the five proteins. This analysis revealed that higher stability appeared to correlate with increased numbers of intramolecular hydrogen bonds and of salt bridges. By contrast, the global hydrophobicity of the protein seemed to play a relatively minor role. A strongly unfavourable balance between charged residues and the presence of a cis-peptide bond preceding a non-proline residue might also contribute to the particularly low stability of two of the enzymes.

Receptor binding properties of four-helix-bundle growth factors deduced from electrostatic analysis

Hormones of the hematopoietin class mediate signal transduction by binding to specific transmembrane receptors. Structural data show that the human growth hormone (hGH) forms a complex with a homodimeric receptor and that hGH is a member of a class of hematopoietins possessing an antiparallel 4-alpha-helix bundle fold. Mutagenesis experiments suggest that electrostatic interactions may have an important influence on hormone-receptor recognition. In order to examine the specificity of hormone-receptor complexation, an analysis was made of the electrostatic potentials of hGH, interleukin-2 (IL-2), interleukin-4 (IL-4), granulocyte colony-stimulating factor (G-CSF), granulocyte-macrophage colony-stimulating factor (GM-CSF), and the hGH and IL-4 receptors. The binding surfaces of hGH and its receptor, and of IL-4 and its receptor, show complementary electrostatic potentials. The potentials of the hGH and its receptor display approximately 2-fold rotational symmetry because the receptor subunits are identical. In contrast, the potentials of GM-CSF and IL-2 lack such symmetry, consistent with their known high affinity for hetero-oligomeric receptors. Analysis of the electrostatic potentials supports a recently proposed hetero-oligomeric model for a high-affinity IL-4 receptor and suggests a possible new receptor binding mode for G-CSF; it also provides valuable information for guiding structural and mutagenesis studies of signal-transducing proteins and their receptors.

19F NMR spectroscopy of [6-19F]tryptophan-labeled Escherichia coli dihydrofolate reductase: equilibrium folding and ligand binding studies

Escherichia coli dihydrofolate reductase contains five tryptophan residues distributed throughout its structure. In order to examine the regions of the protein surrounding these tryptophan residues, we have incorporated 6-fluorotryptophan into the protein. To assign the five resonances observed in the 19F NMR spectrum, five site-directed mutants of the enzyme were made, each with one tryptophan replaced by a phenylalanine. The 19F NMR spectra of the apoprotein, two binary complexes (with NADPH or methotrexate), and one ternary complex (with NADPH and methotrexate) were obtained. The chemical shifts of two of the tryptophan resonances (at positions 22 and 74) are particularly sensitive to ligand binding, while the remaining three (at positions 30, 47, and 133) change, but by less. Since several of the tryptophans are distant from the binding site, these results suggest that 19F NMR can detect ligand-induced changes that are propagated throughout the structure. In the apoprotein, the resonances of the tryptophans at positions 22 and 30 are broadened. In the binary complex with NADPH, the resonances of tryptophans 30 and 74 are broadened while that of tryptophan 22 almost disappears. The line broadening of the tryptophan 22 resonance may reflect motion in that part of the protein, since it is near a region that is disordered in the crystal structure of the apoprotein and its NADP+ complex. In contrast, in the ternary complex this region has a defined structure, and all resonances are of equal intensity and line width. The 19F NMR spectra of the apoprotein and the three ligand complexes were also examined as a function of urea concentration.(ABSTRACT TRUNCATED AT 250 WORDS)

Investigating the role of conserved residue Asp134 in Escherichia coli ribonuclease HI by site-directed random mutagenesis

The role of the conserved Asp134 residue in Escherichia coli ribonuclease HI, which is located at the center of the alpha V helix and lies close to the active site, was analyzed by means of site-directed random mutagenesis. Mutant rnhA genes encoding proteins with ribonuclease H activities were screened by their ability to suppress the ribonuclease-H-dependent, temperature-sensitive growth phenotype of E. coli strain MIC3001. Based on the DNA sequences, nine mutant proteins were predicted to have ribonuclease H activity in vivo. All of these mutant proteins were purified to homogeneity and examined for enzymic activity and protein stability. Among them, only the mutant proteins [D134H]RNase H and [D134N]RNase H were shown to have considerable ribonuclease H activities. Determination of the kinetic parameters revealed that replacement of Asp134 by amino acid residues other than asparagine and histidine dramatically decreased the enzymic activity without seriously affecting the substrate binding. Determination of the CD spectra indicated that none of the mutations seriously affected secondary and tertiary structure. The protein stability was determined from the thermal denaturation curves. All mutant proteins were more stable than the wild-type protein. Such stabilization effects would be a result of a reduction in the negative charge repulsion between Asp134 and the active-site residues, and/or an enhancement of the stability of the alpha V helix. These results strongly suggest that Asp134 does not contribute to the maintenance of the molecular architecture but the carboxyl oxygen at its delta 1 position impacts catalysis.

Comparison of kinetics of formation of helices and hydrophobic core during the folding of staphylococcal nuclease from acid

Our previous kinetic study of the acid and base-induced folding/unfolding transitions of staphylococcal nuclease (SNase) has monitored Trp-140 fluorescence. Trp-140 is located near the flexible COOH terminus and whether or not its fluorescence reflects the overall conformation of the protein has yet to be established. Here we show that the fluorescence intensity of Try-140 correlated closely with the thermal stability (i.e., the calorimetric enthalpy, delta Hcal, of unfolding) of the protein in the pH range 7 to 2.5, confirming that it is a good measure of the overall protein structure. Circular dichroism (CD) at 222 nm, which reflects the helical content of the protein molecule, was used to follow the same folding/unfolding transition in order to compare kinetics of the helix formation and of the appearance of the hydrophobic core. In addition to the three kinetic phases reported earlier with the fluorescence detection, there were a rapid reaction (completed within the 25 ms mixing time of the instrument), which comprised 15% of the signal, and a very slow reaction (time constant > 300 s), which comprised 19% of the signal. With the fluorescence detection for the folding from acid, only 5% of the signal occurred in the rapid phase and there was no reaction slower than 300 s. By comparing kinetics of folding at pH 7 by the CD and fluorescence detection methods, we concluded that: (a) Roughly 15% of the helix content of SNase accumulated before significant changes in the hydrophobic environment (< 5%) of Trp-140 could be detected. The rapid appearance of CD signal reminiscent of helix formation within 25ms would be consistent with the framework model of protein folding. Note, however, that, 15% of the 22% helix content of the protein amounts to an equivalent of fewer than 5 amino acid residues. (b) For the time-resolved signal between 2 ms and 300s, kinetics measured by both properties are consistent with the sequential model, D4 = D3 = D2 = D1 = No (the four Ds are the four substates of the denatured protein and No is the native state). The major folding step by both signals is the D1 to No transition, which gave approximately a 50% change in fluorescence and CD and had a time constant of 160 ms at 25 degree C, pH7.0. (c) The slow phase with the CD signal (>300 s), which is insensitive to Trp-140 fluorescence, has been assigned to be the cis/trans isomerization of Pro-1 17 by other studies. (d) Kinetics in the unfolding direction are consistent with the above interpretation.

Contribution of long-range electrostatic interactions to the stabilization of the catalytic transition state of the serine protease subtilisin BPN'

The possible role of long-range electrostatic interactions on the catalytic activity of the serine protease subtilisin BPN' is investigated using protein engineering techniques. Charged residues on the surface of the enzyme some 13-15 A from the active site were mutated to either neutral or oppositely charged residues. The effect of these mutations on the stability of a complex formed between subtilisin BPN' and Z-Ala-Ala-Pro-Phe-trifluoromethyl ketone, a transition-state inhibitor of the enzyme, was measured. The values of Ki for the complex between the trifluoromethyl ketone and wild-type and mutant subtilisins were used to study the possible contribution of long-range electrostatics in stabilizing the charge distribution in the complex and thus, by analogy, on the transition state of hydrolysis for subtilisin BPN'. Measurement of kon, koff, and Ki for the inhibition of wild-type and mutant subtilisins showed that charged mutations distant from the active site can affect koff and Ki but have little effect on kon. The experimental results show that there is a small, 0.10-0.46 kcal mol-1, but significant contribution to the binding energy from distant surface charges, at low ionic strength. The experimental results were compared to theoretical results, calculated using the DelPhi program for different charge distributions in the complex. The experimental results were found to be most consistent with a complex in which an ion pair is formed between the protonated active site histidine and the ionized oxyanion. Both experimental and theoretical results suggest that long-range electrostatic interactions do play a role in stabilizing the transition-state complex formed between enzyme and inhibitor.(ABSTRACT TRUNCATED AT 250 WORDS)

pH-dependent thermostabilization of Escherichia coli ribonuclease HI by histidine to alanine substitutions

Thermal stabilities of mutant ribonuclease HI proteins from Escherichia coli, in which each of five histidine residues was replaced with alanine, were examined at various pHs. Increases in the Tm values were observed at pH 3.0 for four of the mutant proteins, in which each of the four histidine residues exposed to the solvent was mutated, as compared to the Tm of the wild-type protein. The thermostabilization of three of the mutant proteins was dependent on pH, and only observed at low pH. The thermostabilizing effects of the His-->Ala substitutions were cumulative. The temperature of the midpoint of the transition in the thermal unfolding curves, Tm, of the most stable mutant enzyme, in which His 62, His 83, His 124, and His 127 were replaced by Ala, was 5.5 degrees C higher than that of the wild-type enzyme at pH 3.0. The stability of the wild-type protein decreased as the pH was lowered below pH 4, a condition favoring the protonation of carboxyl groups, probably due to unfavorable electrostatic interactions introduced by the increase in positive charges on the protein. Since imidazole groups are positively charged at pH 3.0, it seems likely that thermal stabilization at pH 3.0 by a His-->Ala substitution would be the result of a reduction in such unfavorable electrostatic interactions. These results suggest that amino acid substitutions that cause a decrease in the number of positive charges on the surface of a protein can be used as a general strategy to enhance protein stability at pH values below pH 4.

External ions and direct membrane effects of enkephalins on Amoeba proteus

We have shown that amphiphilic hormones like vasopressins and endorphins increase water permeability and activate the contractile vacuole (CV) in Amoeba by direct action on the plasma membrane. Using our standard CV assay, the effects of nine opioids, morphine, naloxone and 7 enkephalin derivatives, have been compared in normal, ion-containing growth medium (Chalkley's) and in glass-distilled water. While the absence of external ions does not affect the activity of molecules with net positive charges, opioids with no net charge are devoid of action in glass-distilled water. This shows that, in addition to amphiphilicity, which permits insertion into the lipid core of the membrane, electrostatic interactions with ions and with negative charges of phospholipids membrane and glycocalyx, the thick glycoproteic cell coat of Amoeba, are important in the direct action of these compounds.

Contributions of the polar, uncharged amino acids to the stability of staphylococcal nuclease: evidence for mutational effects on the free energy of the denatured state

In order to quantitate the contributions of the polar, uncharged amino acids to the stability of the native state of staphylococcal nuclease, each of the 13 alanines, 9 glycines, 9 threonines, 6 prolines, 6 glutamines, 6 asparagines, and 3 serines was substituted, either with both alanine and glycine or with 1 of these 2 amino acids plus valine. For each mutant, the stability to reversible denaturation (delta GH2O) was quantitated by determining the Kapp for this reaction as a function of guanidine hydrochloride concentration. In addition, the parameter mGuHCl (= d(delta G)/d[GuHCl]) was calculated from the data. To identify the local structural features responsible for the relatively large and variable changes in delta GH2O and mGuHCl observed for the same type of substitution at different locations in nuclease, statistical correlations were sought between delta GH2O, mGuHCl, and a number of descriptors of the local structure. As with substitutions of the large hydrophobic amino acids [Shortle, D., Stites, W. E., & Meeker, A. K. (1990) Biochemistry 29, 8033-8041], mutation of polar, uncharged residues to Gly leads to a change in stability that, on average, correlates well with the degree to which the wild-type residue is buried. This correlation is especially significant for threonine, an amino acid with both polar and hydrophobic character, but is not demonstrated for the more typically hydrophobic residue alanine. As reported in the previous study of alanine/glycine substitutions of hydrophobic residues, a significant correlation between changes in stability and changes in the value of mGuHCl is again observed, strengthening the conclusion that the putative structural changes in the denatured state which lead to increases or decreases in mGuHCl are responsible for a significant fraction of the stability loss for an average mutant. The existence of this correlation is consistent with the denatured state of wild-type staphylococcal nuclease having evolved to a relatively high free energy via optimization of a balance between a maximal exposure of hydrophobic surface and a minimal gain in chain entropy. On average, mutations are less stable in proportion to the extent of which they perturb this balance. A new and puzzling correlation is reported between the extent of buriedness of a residue in the wild-type native state versus the difference in mGuHCl between the Ala mutation and the Gly mutation at that position.

Mutational studies of protein structures and their stabilities

The fundamental relationship between structure and function has served to guide investigations into the workings of living systems at all levels - from the whole organism to individual cells on down to individual molecules. When X-ray crystallography began to reveal the three-dimensional structures of proteins like myoglobin, lysozyme and RNase A, protein chemists were well prepared to draw inferences about functional mechanisms from the precise positioning of amino acid residues they could see. The close proximity between an amino acid side chain and a chemical group on a bound ligand strongly suggests a functional role for that side chain in binding affinity and specificity. Likewise, the nearly universal finding of large clusters of hydrophobic side chains buried in the core of proteins strongly supports a major functional role of hydrophobic interactions in protein folding and stability. Even though eminently plausible hypotheses like these, grounded in the most fundamental principles of chemistry and the logic of structure–function relationships, become widely accepted and make their way into textbooks, protein chemists have felt compelled to search for ways to test them and put them on a more quantitative basis.

Folding kinetics of T4 lysozyme and nine mutants at 12 degrees C

The kinetics of unfolding and refolding of T4 lysozyme and nine of its mutants have been investigated as a function of guanidinium chloride concentration at 12 degrees C. All show simple two-state, first-order kinetics. Two types of mutants were studied: proline-alanine interchanges and substitutions at position 3 with side chains of varying hydrophobicity. Crystal structures are available for seven of the ten proteins. The effect of mutations on the folding kinetics is more pronounced and complex than on equilibrium thermodynamics. The proteins fall into two broad kinetic classes with one class rather close to the wild type. P86A is a mutant with marked changes in kinetics but only a very small change in stability. Since the 86 position is in the middle of an alpha-helix, the indications are that the helix containing an A residue is more stable in the transition state than one containing a P residue. The other mutants are more complicated, with the refolding and unfolding rates unequally affected by the mutations. On the basis of comparisons with other investigations, we conclude that the rate-determining step in the presence of guanidinium chloride is not the same as in aqueous solution and that it most likely precedes it. The indications are that we are studying the formation of a transition intermediate which is destabilized by the denaturant and which resembles the A intermediate of the framework or molten globule models for protein folding.

Cumulative site-directed charge-change replacements in bacteriophage T4 lysozyme suggest that long-range electrostatic interactions contribute little to protein stability

Bacteriophage T4 lysozyme is a basic molecule with an isoelectric point above 9.0, and an excess of nine positive charges at neutral pH. It might be expected that it would be energetically costly to bring these out-of-balance charges from the extended, unfolded, form of the protein into the compact folded state. To determine the contribution of such long-range electrostatic interactions to the stability of the protein, five positively charged surface residues, Lys16, Arg119, Lys135, Lys147 and Arg154, were individually replaced with glutamic acid. Eight selected double, triple and quadruple mutants were also constructed so as to sequentially reduce the out-of-balance formal charge on the molecule from +9 to +1 units. Each of the five single variant proteins was crystallized and high-resolution X-ray analysis confirmed that each mutant structure was, in general, very similar to the wild-type. In the case of R154E, however, the Arg154 to Glu replacement caused a rearrangement in which Asp127 replaced Glu128 as the capping residue of a nearby alpha-helix. The thermal stabilities of all 13 variant proteins were found to be fairly similar, ranging from 0.5 kcal/mol more stable than wild-type to 1.7 kcal/mol less stable than wild-type. In the case of the five single charge-change variants, for which the structures were determined, the changes in stability can be rationalized in terms of changes in local interactions at the site of the replacement. There is no evidence that the reduction in the out-of-balance charge on the molecule increases the stability of the folded relative to the unfolded form, either at pH 2.8 or at pH 5.3. This indicates that long-range electrostatic interactions between the substituted amino acid residues and other charged groups on the surface of the molecule are weak or non-existent. Furthermore, the relative stabilities of the multiple charge replacement mutant proteins were found to be almost exactly equal to the sums of the relative stabilities of the constituent single mutant proteins. This also clearly indicates that the electrostatic interactions between the replaced charges are negligibly small. The activities of the charge-change mutant lysozymes, as measured by the rate of hydrolysis of cell wall suspensions, are essentially equal to that of the wild-type lysozyme, but on a lysoplate assay the mutant enzymes appear to have higher activity.(ABSTRACT TRUNCATED AT 400 WORDS)

Transient intermediates in the folding of dihydrofolate reductase as detected by far-ultraviolet circular dichroism spectroscopy

The kinetics of the reversible folding and unfolding of Escherichia coli dihydrofolate reductase have been studied by stopped-flow circular dichroism in the peptide region at pH 7.8 and 15 degrees C. The reactions were induced by concentration jumps of a denaturant, urea. The method can detect various intermediates transiently populated in the reactions although the equilibrium unfolding of the protein is apparently approximated by a two-state reaction. The results can be summarized as follows. (1) From transient circular dichroism spectra measured as soon as the refolding is started, a substantial amount of secondary structure is formed in the burst phase, i.e., within the dead time of stopped-flow mixing (18 ms). (2) The kinetics from this burst-phase intermediate to the native state are multiphasic, consisting of five phases designated as tau 1, tau 2, tau 3, tau 4, and tau 5 in increasing order of the reaction rate. Measurements of the kinetics at various wavelengths have provided kinetic difference circular dichroism spectra for the individual phases. (3) The tau 5 phase shows a kinetic difference spectrum consistent with an exciton contribution of two aromatic residues in the peptide CD region. The absence of the tau 5 phase in a mutant protein, in which Trp 74 is replaced by leucine, suggests that Trp 74 is involved in the exciton pair and that the tau 5 phase reflects the formation of a hydrophobic cluster around Trp 74. From the similarity of the kinetic difference spectrum to the difference between the native spectra of the mutant and wild-type proteins, it appears that Trp 47 is the partner in the exciton pair and that the structure formed in the tau 5 phase persists during the later stages of folding. (4) The later stages of folding show kinetic difference spectra that can be interpreted by rearrangement of secondary structure, particularly the central beta sheet of the protein. The pairwise similarities in the spectrum between the tau 3 and tau 4 phases, and between the tau 1 and tau 2 phases, also suggest the presence of two parallel folding channels for refolding. (5) The unfolding kinetics show three to four phases and are interpreted in terms of the presence of multiple native species. The total ellipticity change in kinetic unfolding reaction, however, agrees with the ellipticity difference between the native and unfolding states, indicating the absence of the burst phase in unfolding.(ABSTRACT TRUNCATED AT 400 WORDS)

Effects of point mutations in a hinge region on the stability, folding, and enzymatic activity of Escherichia coli dihydrofolate reductase

The role of a hinge region in the folding, stability, and activity of Escherichia coli dihydrofolate reductase was investigated with three site-directed mutants at valine-88, the central residue of the hinge. The three mutants, V88A and V88I and a valine-88 deletion, were created to perturb the packing of hydrophobic residues in the interior of a loose turn formed by residues 85-91. Deleting the valine-88 residue destabilized the protein by 2.93 +/- 0.6 kcal/mol as determined by equilibrium unfolding transitions in urea monitored by circular dichroism at 20 degrees C. Substitution of alanine for valine-88 stabilized the protein by -0.20 +/- 0.02 kcal/mol, and the isoleucine substitution was mildly destabilizing by 1.73 +/- 0.2 kcal/mol. Although there was no clear correlation between side-chain volume and stability, these results suggest that side-chain interactions in the interior of the turn influence the folding and stability of dihydrofolate reductase. The specific activity of the valine deletion mutant was approximately twice that of the wild-type protein while the specific activities of the V88A and V88I proteins were only slightly greater than the wild type. The full time courses of the reactions catalyzed by the mutants were almost identical with that for the wild type, indicating no major changes in the kinetic mechanism. Additionally, the rate constants associated with interconversion between various forms of the apoenzyme were identical for the mutant and wild-type enzymes. The rate constants for refolding transitions were examined by dilution of urea-inactivated protein.(ABSTRACT TRUNCATED AT 250 WORDS)

Isolation and in vitro characterization of temperature-sensitive mutants of the bacteriophage f1 gene V protein

In vivo selections were used to isolate 43 temperature-sensitive gene V mutants of the bacteriophage f1 from a collection of mutants constructed by saturation mutagenesis of the gene. The sites of temperature-sensitive substitutions are found in both the beta-sheets and the turns of the protein, and some sites are exposed to the solvent while others are not. Thirteen of the variant proteins were purified and characterized to evaluate their free energy changes upon unfolding and their affinities for single-stranded DNA, and eight were tested for their tendencies to aggregate at 42 degrees C. Each of the three temperature-sensitive mutants at buried sites and six of ten at surface sites had free energy changes of unfolding substantially lower (less stabilizing) than the wild-type at 25 degrees C. A seventh mutant at a surface site had a substantially altered unfolding transition and its free energy of unfolding was not estimated. The affinities of the mutant proteins for single-stranded DNA varied considerably, but two mutants at a surface site, Lys69, had much weaker binding to single-stranded DNA than any of the other mutants, while two mutants at another surface site, Glu30, had the highest DNA-binding affinities. The wild-type gene V protein is stable at 42 degrees C, but six of the eight mutants tested aggregated within a few minutes and the remaining two aggregated within 30 minutes at this temperature. Overall, each of the temperature-sensitive proteins tested had a tendency to aggregate at 42 degrees C, and most also had either a low free energy of unfolding (at 25 degrees C), or weak DNA binding. We suggest that any of these properties can lead to a temperature-sensitive gene V phenotype.