Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “electrostatics”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55Linked to original sources

Powering DNA repair through substrate electrostatic interactions.

The reaction catalyzed by the DNA repair enzyme uracil DNA glycosylase (UDG) proceeds through an unprecedented stepwise mechanism involving a positively charged oxacarbenium ion sugar and uracil anion leaving group. Here we use a novel approach to evaluate the catalytic contribution of electrostatic interactions between four essential phosphodiester groups of the DNA substrate and the cationic transition state. Our strategy was to substitute each of these phosphate groups with an uncharged (R)- or (S)-methylphosphonate linkage (MeP). We then compared the damaging effects of these methylphosphonate substitutions on catalysis with their damaging effects on binding of a cationic 1-azadeoxyribose (1-aza-dR(+)) oxacarbenium ion analogue to the UDG-uracil anion binary complex. A plot of log k(cat)/K(m) for the series of MeP-substituted substrates against log K(D) for binding of the 1-aza-dR(+) inhibitors gives a linear correlation of unit slope, confirming that the electronic features of the transition state resemble that of the 1-aza-dR(+), and that the anionic backbone of DNA is used in transition state stabilization. We estimate that all of the combined phosphodiester interactions with the substrate contribute 6-8 kcal/mol toward lowering the activation barrier, a stabilization that is significant compared to the 16 kcal/mol catalytic power of UDG. However, unlike groups of the enzyme that selectively stabilize the charged transition state by an estimated 7 kcal/mol, these phosphodiester groups also interact strongly in the ground state. To our knowledge, these results provide the first experimental evidence for electrostatic stabilization of a charged enzymatic transition state and intermediate using the anionic backbone of DNA.

Catalysis↗

Dissecting structural and electrostatic interactions of charged groups in alpha-sarcin. An NMR study of some mutants involving the catalytic residues.

The cytotoxic ribonuclease alpha-sarcin is the best characterized member of the ribotoxin family. Ribotoxins share a common structural core, catalytic residues, and active site topology with members of the broader family of nontoxic microbial extracellular RNases. They are, however, much more specific in their biological action. To shed light on the highly specific alpha-sarcin activity, we have evaluated the structural and electrostatic interactions of its charged groups, by combining the structural and pK(a) characterization by NMR of several variants with theoretical calculations based on the Tanford-Kirkwood and Poisson-Boltzmann models. The NMR data reveal that the global conformation of wild-type alpha-sarcin is preserved in the H50Q, E96Q, H137Q, and H50/137Q variants, and that His137 is involved in an H-bond that is crucial in maintaining the active site structure and in reinforcing the stability of the enzyme. The loss of this H-bond in the H137Q and H50/137Q variants modifies the local structure of the active site. The pK(a) values of active site groups H50, E96, and H137 in the four variants have been determined by two-dimensional NMR. The catalytic dyad of E96 and H137 is not sensitive to charge replacements, since their pK(a) values vary less than +/-0.3 pH unit with respect to those of the wild type. On the contrary, the pK(a) of His50 undergoes drastic changes when compared to its value in the intact protein. These amount to an increase of 0.5 pH unit or a decrease of 1.1 pH units depending on whether a positive or negative charge is substituted at the active site. The main determinants of the pK(a) values of most of the charged groups in alpha-sarcin have been established by considering the NMR results in conjunction with those derived from theoretical pK(a) calculations. With regard to the active site residues, the H50 pK(a) is chiefly influenced by electrostatic interactions with E96 and H137, whereas the effect of the low dielectric constant and the interaction with R121 appear to be the main determinants of the altered pK(a) value of E96 and H137. Charge-charge interactions and an increased level of burial perturb the pK(a) values of the active site residues of alpha-sarcin, which can account for its reduced ribonucleolytic activity and its high specificity.

Amino Acid Substitution↗

Critical balance of electrostatic and hydrophobic interactions is required for beta 2-microglobulin amyloid fibril growth and stability.

Investigation of factors that modulate amyloid formation of proteins is important to understand and mitigate amyloid-related diseases. To understand the role of electrostatic interactions and the effect of ionic cosolutes, especially anions, on amyloid formation, we have investigated the effect of salts such as NaCl, NaI, NaClO(4), and Na(2)SO(4) on the amyloid fibril growth of beta(2)-microglobulin, the protein involved in dialysis-related amyloidosis. Under acidic conditions, these salts exhibit characteristic optimal concentrations where the fibril growth is favored. The presence of salts leads to an increase in hydrophobicity of the protein as reported by 8-anilinonaphthalene-1-sulfonic acid, indicating that the anion interaction leads to the necessary electrostatic and hydrophobic balance critical for amyloid formation. However, high concentrations of salts tilt the balance to high hydrophobicity, leading to partitioning of the protein to amorphous aggregates. Such amorphous aggregates are not competent for fibril growth. The order of anions based on the lowest concentration at which fibril formation is favored is SO(4)(2)(-) > ClO(4)(-) > I(-) > Cl(-), consistent with the order of their electroselectivity series, suggesting that preferential anion binding, rather than general ionic strength effect, plays an important role in the amyloid fibril growth. Anion binding is also found to stabilize the amyloid fibrils under acidic condition. Interestingly, sulfate promotes amyloid growth of beta(2)-microglobulin at pH between 5 and 6, closer to its isoelectric point. Considering the earlier studies on the role of glycosaminoglycans and proteoglycans (i.e., sulfated polyanions) on amyloid formation, our study suggests that preferential interaction of sulfate ions with amyloidogenic proteins may have biological significance.

Acids↗

Apparent Debye-Huckel electrostatic effects in the folding of a simple, single domain protein.

We have monitored the effects of salts and denaturants on the folding of the simple, two-state protein FynSH3. As predicted by Debye-Huckel limiting law, both the stability and (log) folding rate of FynSH3 increase nearly perfectly linearly (r(2)> 0.99) with the square root of ionic strength upon increasing concentrations of the relatively nonchaotropic salt sodium chloride. The stability of FynSH3 is also linear in square root ionic strength when the relatively nonchaotropic salts sodium bromide, potassium bromide, and potassium chloride are employed. Comparison of the kinetic and equilibrium effects of sodium chloride suggests that the electrostatic interactions formed in the folding transition state are approximately 50% as destabilizing as those formed in the native state, presumably reflecting the more compact nature of the latter. In contrast, the relationship between concentration and folding kinetics is more complex when the highly chaotropic salt guanidine hydrochloride (GuHCl) is employed. At moderate to high GuHCl concentrations the net effect of the linear, presumably chaotrope-induced deceleration and the presumed, square root-dependent ionic strength-induced acceleration is well approximated as linear, thus accounting for the observation of "chevron behavior" (log folding rate linear in denaturant concentration) typically reported for the folding of single domain proteins. At very low GuHCl concentrations, however, significant kinetic rollover is observed. This rollover is reasonably well fitted as a sum of a linear, presumably chaotropic effect and a square root-dependent, presumably electrostatic effect. These results thus not only provide insight into the nature of the folding transition state but also suggest that caution is in order when extrapolating GuHCl-based chevrons to estimate folding rates in the absence of denaturant and in interpreting deviations from chevron linearity as evidence for non-two-state kinetics.

Bromides↗

Role of electrostatics and salt bridges in stabilizing the compound I radical in ascorbate peroxidase.

Cytochrome c (CcP) and ascorbate peroxidase (APX) are heme peroxidases which have very similar active site structures yet differ substantially in the properties of compound I, the intermediate formed upon reaction with peroxides. Although both peroxidases have a tryptophan in the proximal heme pocket, Trp191 in CcP and Trp179 in APX, only Trp191 in CcP forms a stable cation radical while APX forms the more traditional porphyrin pi-cation radical. Previous work [Barrows, T. P., et al. (2004)Biochemistry 43, 8826-8834] has shown that converting three methionine residues in the cytochrome c peroxidase (CcP) proximal heme pocket to the corresponding residues in APX dramatically decreased the stability of the Trp191 radical in CcP compound I. On the basis of these results, we reasoned that replacing the analogous residues at positions 160, 203, and 204 in APX with methionine should stabilize a Trp179 radical in APX compound I. Steady- and transient-state kinetics of this mutant (designated APX3M) show a significant destabilization of the native porphyrin pi-radical, while electron paramagnetic resonance (EPR) studies show an increase in the intensity of the signal at g = 2.006 with characteristics consistent with formation of a Trp radical. This hypothesis was tested by replacing Trp179 with Phe in the APX3M background. The EPR spectrum of this mutant was very similar to that of the CcP W191G mutant which is known to form a tyrosine radical. Previously published theoretical studies [Guallar, V., et al. (2003) Proc. Natl. Acad. Sci. U.S.A. 100, 6998-7002] suggest that electrostatic shielding of the heme propionates also plays a role in the stability of the porphyrin radical. Arg172 in APX hydrogen bonds with one of the heme propionates. Replacing Arg172 with an asparagine residue in the APX3M background generates a mutant which no longer forms the full complement of the compound I porphyrin pi-radical. These results suggest that the electrostatics of the proximal pocket and the shielding of propionate groups by salt bridges are critical factors controlling the location of a stable compound I radical in heme peroxidases.

Arginine↗

Role of Ca2+ in the electrostatic stability and the functional activity of the globular domain of human C1q.

C1q is the recognition subunit of the classical pathway of the complement system and a major connecting link between classical pathway-driven innate immunity and IgG- or IgM-mediated acquired immunity. The basic structural subunit of C1q is composed of an N-terminal triple-helical collagen-like region and a C-terminal heterotrimeric globular head domain (gC1q) that is made up of individual A, B, and C chains. Recent crystallographic studies have revealed that the gC1q domain, which is the main target-binding region of C1q, has a compact and spherical heterotrimeric assembly, held together by both electrostatic and nonpolar interactions, with quasi-3-fold symmetry. A characteristic feature of the gC1q domain is the presence of a exposed Ca(2+) located near the apex. We have investigated, using theoretical and experimental approaches, the role of Ca(2+) in the electrostatic stability and target-binding properties of the native C1q as well as recombinant monomeric forms of the C-terminal regions of the A, B, and C chains. Here, we report that Ca(2+) primarily influences the target recognition properties of C1q toward IgG, IgM, C-reactive protein, and pentraxin 3. At pH 7.4, the loss of Ca(2+) leads to changes in the direction of electric moment from coaxial (where the putative C-reactive protein-binding site is located) to perpendicular to the molecular axis (toward the most likely IgG-binding site), which appears important for target recognition by C1q and subsequent complement activation.

Calcium↗

Electrostatic environment of hemes in proteins: pK(a)s of hydroxyl ligands.

The pK(a)s of ferric aquo-heme and aquo-heme electrochemical midpoints (E(m)s) at pH 7 in sperm whale myoglobin, Aplysia myoblogin, hemoglobin I, heme oxygenase 1, horseradish peroxidase and cytochrome c oxidase were calculated with Multi-Conformation Continuum Electrostatics (MCCE). The pK(a)s span 3.3 pH units from 7.6 in heme oxygenase 1 to 10.9 in peroxidase, and the E(m)s range from -250 mV in peroxidase to 125 mV in Aplysia myoglobin. Proteins with higher in situ ferric aquo-heme pK(a)s tend to have lower E(m)s. Both changes arise from the protein stabilizing a positively charged heme. However, compared with values in solution, the protein shifts the aquo-heme E(m)s more than the pK(a)s. Thus, the protein has a larger effective dielectric constant for the protonation reaction, showing that electron and proton transfers are coupled to different conformational changes that are captured in the MCCE analysis. The calculations reveal a breakdown in the classical continuum electrostatic analysis of pairwise interactions. Comparisons with DFT calculations show that Coulomb's law overestimates the large unfavorable interactions between the ferric water-heme and positively charged groups facing the heme plane by as much as 60%. If interactions with Cu(B) in cytochrome c oxidase and Arg 38 in horseradish peroxidase are not corrected, the pK(a) calculations are in error by as much as 6 pH units. With DFT corrected interactions calculated pK(a)s and E(m)s differ from measured values by less than 1 pH unit or 35 mV, respectively. The in situ aquo-heme pK(a) is important for the function of cytochrome c oxidase since it helps to control the stoichiometry of proton uptake coupled to electron transfer [Song, Michonova-Alexova, and Gunner (2006) Biochemistry 45, 7959-7975].

Animals↗

Electrostatic contributions to residue-specific protonation equilibria and proton binding capacitance for a small protein.

Charge-charge interactions in proteins are important in a host of biological processes. Here we use 13C NMR chemical shift data for individual aspartate and glutamate side chain carboxylate groups to accurately detect site-specific protonation equilibria in a variant of the B1 domain of protein G (PGB1-QDD). Carbon chemical shifts are dominated by changes in the electron distribution within the side chain and therefore excellent reporters of the charge state of individual groups, and the data are of high precision. We demonstrate that it is possible to detect local charge interactions within this small protein domain that stretch and skew the chemical shift titration curves away from "ideal" behavior and introduce a framework for the analysis of such convoluted data to study local charge-charge interactions and electrostatic coupling. It is found that, due to changes in electrostatic potential, the proton binding affinity, Ka, of each carboxyl group changes throughout the titration process and results in a linearly pH dependent pKa value. This result could be readily explained by calculations of direct charge-charge interactions based on Coulomb's law. In addition, the slope of pKa versus pH was dependent on screening by salt, and this dependence allowed the selective study of charge-charge interactions. For PGB1-QDD, it was established that mainly differences in self-energy, and not direct charge-charge interactions, are responsible for shifted pKa values within the protein environment.

Hydrogen-Ion Concentration↗

Characterization of four covalently-linked yeast cytochrome c/cytochrome c peroxidase complexes: Evidence for electrostatic interaction between bound cytochrome c molecules.

Four covalent complexes between recombinant yeast cytochrome c and cytochrome c peroxidase (rCcP) were synthesized via disulfide bond formation using specifically designed protein mutants (Papa, H. S., and Poulos, T. L. (1995) Biochemistry 34, 6573-6580). One of the complexes, designated V5C/K79C, has cysteine residues replacing valine-5 in rCcP and lysine-79 in cytochrome c with disulfide bond formation between these residues linking the two proteins. The V5C/K79C complex has the covalently bound cytochrome c located on the back-side of cytochrome c peroxidase, approximately 180 degrees from the primary cytochrome c-binding site as defined by the crystallographic structure of the 1:1 noncovalent complex (Pelletier, H., and Kraut J. (1992) Science 258, 1748-1755). Three other complexes have the covalently bound cytochrome c located approximately 90 degrees from the primary binding site and are designated K12C/K79C, N78C/K79C, and K264C/K79C, respectively. Steady-state kinetic studies were used to investigate the catalytic properties of the covalent complexes at both 10 and 100 mM ionic strength at pH 7.5. All four covalent complexes have catalytic activities similar to those of rCcP (within a factor of 2). A comprehensive study of the ionic strength dependence of the steady-state kinetic properties of the V5C/K79C complex provides evidence for significant electrostatic repulsion between the two cytochromes bound in the 2:1 complex at low ionic strength and shows that the electrostatic repulsion decreases as the ionic strength of the buffer increases.

Cytochrome-c Peroxidase↗

Shift of the special pair redox potential: electrostatic energy computations of mutants of the reaction center from Rhodobacter sphaeroides.

Shifts of the special pair redox potential of the photosynthetic reaction center of Rhodobacter sphaeroides are considered for several point mutations [Lin. X., Murchison, H. A., Nagarijan, V., Parson, W. W., Allen, J. P., & Williams, J. C. (1994) Proc. Natl. Acad. Sci. U.S.A. 91, 10265-10269] in the neighborhood of the special pair. The shifts are calculated from electrostatic energies by solving Poisson's equation for energy-minimized structures of the reaction center. Different conditions for the evaluation of the electrostatic energy are probed. To test the influence of the hydrogen bonding at the acetyl groups of the special pair, the orientation and torsion potential of the acetyl groups are varied. The calculated shifts of the midpoint potential of double and triple mutants can approximately be obtained from the corresponding shifts of the single point mutations. The calculated shifts agree with the measured values for all single and double mutants considered. However, a clear decision between different acetyl group conformations was only possible for the mutants HF(L168) and HF(L168) + LH(L131) where the calculated shifts of the redox potential agree with experiments only if the acetyl oxygen atom at DM points toward the Mg2+ ion of DL. This is corroborated by computations of the interaction energy of the acetyl group at DM, which adopts a lower value in the wild-type reaction center if its oxygen atom is bonded to the Mg2+ ion of DL.

Bacteriochlorophylls↗

Electrostatic environment of the tryptophylquinone cofactor in methylamine dehydrogenase: evidence from resonance Raman spectroscopy of model compounds.

Methylamine dehydrogenase (MADH) utilizes its endogenous tryptophan tryptophylquinone (TTQ) as a cofactor in enzymatic catalysis, with the C6 carbonyl of the quinone implicated as the site of attack by substrates and other nucleophiles. Resonance Raman (RR) spectroscopy provides an ideal method for investigating the state of this carbonyl group whose C==O stretch is distinct from other vibrational modes of the cofactor and is readily identified by its shift to lower energy in H218O. In a series of indole 6,7-quinone models for TTQ, the in-phase stretching vibration of the two C==O groups occurs at 1650 cm-1 in nonpolar solvents and shifts to 1638 cm-1 in H2O. The absorption maximum undergoes an analogous shift from 400 to 425 mm. The spectral properties of the indole quinones in H2O approach the corresponding values in Thiobacillus versutus MADH (C==O stretch at 1612 cm-1, lamdamax at 440mm) and are indicative of strongly hydrogen bonding of the C==O and NH groups of the cofactor in the native enzyme. Addition of monovalent cations [NH4+,Cs+, and (CH3)3NH+] to MADH causes further increases in the lamdamax and decreases in the frequency of the C==O stretch[1590 cm-1 with (CH3)3NH+]. This implies a strong electrostatic interaction between monovalent cations and a carbonyl oxygen (most likely at C6) in TTQ. The fact that these cations behave as competitive inhibitors of the methylamine substrate suggests that methylamine binds to the same location in the enzyme prior to its covalent reaction with the cofactor. Addition of monovalent cations to the one-electron-reduced semiquinone form MADH results in RR spectral shifts for a number of vibrational modes of the cofactor. Thus, the ability of monovalent cations to promote and stabilize the formation of the semiquinone intermediate is also due to their direct electrostatic interaction with the TTQ cofactor.

Catalysis↗

A symmetry-driven search for electrostatic interaction partners in charybdotoxin and a voltage-gated K+ channel.

A structural model of charybdotoxin bound to a Shaker K+ channel has emerged from mechanistic and mutagenic analysis of toxin-channel interactions. We test this model by predicting through-space electrostatic interactions between specific pairs of channel-toxin residues. Dissociation constants of channel-toxin variants, determined by radiolabeled toxin binding to Shaker-transfected COS membrane fragments, were used to identify pairs of residues located closely enough to interact electrostatically. The results further refine the structural model of the bound complex and produce a more detailed view of the vestibule of the Shaker channel.

Animals↗

The role of electrostatic interactions in the assembly of the factor X activating complex on both activated platelets and negatively-charged phospholipid vesicles.

Factor X was activated by factor IXa on the surface of either activated platelets or artificial lipid vesicles in the presence of different NaCl concentrations. The Vmax of reactions using platelets was optimal at physiologic [NaCl] both in the absence and in the presence of factor VIIIa. In contrast, the Vmax of reactions using vesicles decreased with increasing [NaCl] in the absence of factor VIIIa, and increased with increasing [NaCl] when cofactor was present. In the absence of factor VIIIa, the EC50FIXa' although stable to changes in [NaCl] in platelet-supported reactions, was found to increase significantly as [NaCl] increased in vesicle-supported reactions and correlated with the decreased Vmax. Thus, in contrast to platelet-supported reactions, enzyme interaction with negatively-charged vesicles was highly dependent upon electrostatic interactions. In the presence of factor VIIIa, the EC50FIXa of vesicle-supported reactions decreased with increasing [NaCl], indicating that interactions between FIXa and FVIIIa can increase enzyme affinity when fewer ionic interactions are favored. The EC50FVIIIa was insensitive to changes in [NaCl] on both surfaces. The K(m)app derived from platelet-supported titrations of factor X was lowest just above physiological [NaCl], whereas on vesicles K(m)app was minimal at the lowest [NaCl] tested. Thus, the direct interaction of factor X and factor IXa with the artificial lipid surface is highly dependent upon ionic interactions with the negatively-charged polar heads of phospholipids. However, the interaction of factor IXa and factor X with the activated platelet surface must rely both on electrostatic interactions with lipid and on other interactions provided by surface proteins.

Animals↗

Electrostatic effects in the kinetics of coenzyme binding to isozymes of alcohol dehydrogenase from horse liver.

The kinetic mechanism for the binding of NAD+ and NADH to the EE and SS isozymes of alcohol dehydrogenase (LADH) was studied between pH 7 and pH 10 by monitoring the quenching of tryptophan fluorescence. A consistent interpretation of all data was only possible by introducing a two-step binding mechanism. The first binding step is related to docking of the adenosine part of the coenzymes and the subsequent isomerization to the binding of the nicotinamide part. At high NADH concentrations an additional slow isomerization was identified as a conformational transition of the protein. A pH dependence for NADH binding is observed which is restricted to changes in the binding kinetics of the adenosine moiety going from pH 7 to pH 10, a tendency which is similar also for NAD+. This is attributed to pH-dependent variations in electrostatic attractions acting as a steering force of the docking process. The nicotinamide docking of NADH is equally fast for both isozymes and pH-independent over the measured range, whereas this docking equilibrium for NAD+ is pH-dependent for EE- and SS-LADH alike and the rate of association comparable. Presumably, a GluEE-366-LysSS substitution results in a stronger binding and faster association of both oxidized and reduced cofactor to the SS isozyme. A structural proof is presented for coenzyme-competitive binding of a sulfate ion, resulting in electrostatic shielding.

Alcohol Dehydrogenase↗

Use of 1H NMR spectroscopy and computer simulations To analyze histidine pKa changes in a protein tyrosine phosphatase: experimental and theoretical determination of electrostatic properties in a small protein.

In bovine low Mr protein tyrosine phosphatase, the pKa values of His-66 and His-72 are 8.3 and 9.2, respectively. These unusually high values were hypothesized to be caused by electrostatic interactions with several nearby negatively charged groups. To test this, mutant enzymes were made in which one or more carboxylate side chains were removed or introduced near the histidines. Michaelis kinetic parameters, measured using p-nitrophenyl phosphate as a substrate, indicated that all mutant enzymes retained approximately 50% or more of the activity of wild-type enzyme. The effect that each mutation had on the pKa of the nearby histidine was monitored by 1H NMR spectroscopy using the MLEV-17 pulse sequence to filter out the broad interfering amide resonances in the spectra. Independently, computer simulations of the pKas were obtained using the finite difference method to solve the linear Poisson-Boltzmann equation. The proximity of a charged residue to the titrating histidine imidazole largely determined the extent of the pKa perturbation. The change in pKa for His-72 in the mutant enzymes was -1.69 units for D42A, -2.36 units for E23A, -2.99 units for E23A/D42A, and unchanged for E139A and Q143E. Thus, the pKa of His-72 in the double mutant E23A/D42A decreased to nearly that of a free histidine imidazole group. The His-66 pKa change was -1.25 units for E139A and was not significant for the other mutants. His-66, Glu-139, and Gln-143 are at the protein surface and much more exposed to the higher solvent dielectric compared to His-72, Glu-23, and Asp-42. These structural characteristics explain the smaller decrease in the observed pKa of His-66 for the E139A mutant compared to the decrease in the pKa of His-72 when a single nearby carboxylate was removed. These observations were adequately predicted by theoretical electrostatic calculations using the Poisson-Boltzmann equation as a model for a solute molecule of low dielectric in solution of high dielectric.

Animals↗

Electrostatic channeling of substrates between enzyme active sites: comparison of simulation and experiment.

Recent simulation work has indicated that channeling of charged substrates between the active sites of bifunctional enzymes or bienzyme complexes can be significantly enhanced by favorable interactions with the electrostatic field of the enzymes. The results of such simulations are expressed in terms of transfer efficiencies, which describe the probability that substrate leaving the active site of the first enzyme will reach the active site of the second enzyme before escaping out into bulk solution. The experimental indicators of channeling, on the other hand, are factors such as a decrease in the transient (lag) time for appearance of the final product of the coupled enzyme reaction or a decrease in the susceptibility of the overall reaction rate to the presence of competing enzymes or competitive inhibitors. The work reported here aims to establish a connection between the transfer efficiencies obtained from simulation, with the above-mentioned experimental observables. This is accomplished by extending previously reported analytical approaches to combine the simulated transfer efficiency with the Michaelis-Menten kinetic parameters Km and Vmax of the enzymes involved; expressions are derived to allow both transient times and steady state rates to be calculated. These results are applied to the two systems that have been studied both theoretically and experimentally. In the first case, that of the bifunctional enzyme dihydrofolate reductase-thymidylate synthase (DHFR-TS), the experimentally observed decrease in transient times is found to be consistent with a transfer efficiency of >/=80%. In the second case, that of a citrate synthase-malate dehydrogenase fusion protein, a transfer efficiency of 73% is consistent with the experimental transient time measurements. Separate and independent analysis of the effects of adding the competing enzyme aspartate aminotransferase gives a transfer efficiency of 69%, in excellent agreement with the transient time results. The transfer efficiencies thus obtained from experimental results are in both cases in good agreement with those obtained from simulations that include electrostatic interactions. One important discrepancy between simulation and experiment, is however, found in the reported effects of adding a competitive inhibitor in the DHFR-TS system: qualitatively different results are expected from the theoretical analysis. A possible reason for this apparent contradiction is discussed.

Aspartate Aminotransferases↗

Electrostatic effects in DNA bending by GCN4 mutants.

DNA architecture has been shown to be important for cellular processes such as activation of transcription, recombination, and replication. Many proteins reconfigure the shape of duplex DNA upon binding. Previous experiments have shown that some members of the eukaryotic bZIP family of DNA binding proteins appear to bend DNA, while others do not. We are exploring the role of electrostatic effects in DNA bending by bZIP proteins. The yeast bZIP transcription factor GCN4 does not induce DNA bending in vitro. Previously we substituted basic residues for three neutral amino acids in GCN4 to produce a GCN4 derivative that bends DNA by approximately 15 degrees. This result is consistent with a model of induced DNA bending wherein excess positive charge in proximity to one face of the double helix neutralizes local phosphate diester anions resulting in a laterally-asymmetric charge distribution along the DNA. Such an unbalanced charge distribution can result in collapse of the DNA toward the neutralized surface. We now present a more comprehensive analysis of electrostatic effects in DNA bending by GCN4 derivatives. It is shown that the direction and extent of DNA bending by these derivatives are a linear function of the charges of the amino acids adjacent to the basic domain of the protein. This relation holds over the charge range +6 (16 degrees bend toward the minor groove) to -6 (25 degrees bend toward the major groove).

Amino Acid Sequence↗

Electrostatic stabilization in methionine aminopeptidase from hyperthermophile Pyrococcus furiosus.

The thermostability of methionine aminopeptidase from a hyperthermophile P. furiosus (PfMAP) was extremely high: the denaturation temperature was 106.2 degreesC at pH 10.2. To explore the contribution of electrostatic interaction to the superior thermostability of PfMAP, the thermostability of PfMAP was examined by differential scanning calorimetry (DSC) in various salt concentrations in the acidic region far from the isoelectric point of PfMAP. (1) In 20 mM glycine buffer, the DSC curve of PfMAP exhibited a single peak. Transition temperatures (Tm) were lowered with decreasing pH from 4 to 3. The heat denaturation of PfMAP was not reversible. (2) Denaturation enthalpy (DeltaH) measured at different pHs linearly correlated with Tm up to 102 degreesC, suggesting that the denaturation heat capacity (DeltaCp) for PfMAP is constant up to 100 degreesC. DeltaCp was estimated to be 0.82 J K-1 g-1. (3) In the presence of 10-100 mM KCl at pH 3.2, two peaks appeared on the DSC curves. The first peak shifted to lower temperatures with increasing concentration of KCl and, oppositely, the second one to higher temperatures. It was found that the first and second peaks originated from the heat denaturation of the native form of PfMAP and the melting of the non-native associated form having molten globule-like structure, respectively, judged from the CD spectra and ultracentrifugation analyses. This indicates the following: first, the attractive electrostatic interaction is an important factor in stabilizing the native form of PfMAP; second, the presence of KCl stimulates the formation of the molten globule-like state of PfMAP and stabilizes it. (4) In a comparison of the sequence and crystal structure of PfMAP, which has been recently determined (1xgs.pdb), with those of MAP from Escherichia coli (EcMAP), it was predicted that the extra four short-range ion pairs less than 3 A involved in PfMAP are crucial candidates as determinants for the superior thermostability of PfMAP.

Aminopeptidases↗