The medium reorganization energy for the charge transfer reactions in proteins

A low static dielectric permittivity of proteins causes the low reorganization energies for the charge transfer reactions inside them. This reorganization energy does not depend on the pre-existing intraprotein electric field. The charge transferred inside the protein interacts with its aqueous surroundings; for many globular proteins, the effect of this surroundings on the reorganization energy is comparable with the effect of reorganization of the protein itself while for the charge transfer in the middle of membrane the aqueous phase plays a minor role. Reorganization energy depends strongly on the system considered, and hence there is no sense to speak on the "protein reorganization energy" as some permanent characteristic parameter. We employed a simple algorithm for calculation of the medium reorganization energy using the numerical solution of the Poisson-Boltzmann equation. Namely, the reaction field energy was computed in two versions - all media having optical dielectric permittivity, and all the media with the static one; the difference of these two quantities gives the reorganization energy. We have calculated reorganization energies for electron transfer in cytochrome c, various ammine-ruthenated cytochromes c, azurin, ferredoxin, cytochrome c oxidase, complex of methylamine dehydrogenase with amicyanin, and for proton transfer in α-chymotrypsin. It is shown that calculation of the medium reorganization energy can be a useful tool in analysis of the mechanisms of the charge transfer reactions in proteins.     

Conformational changes and molecular mobility in plasticized proteins

Most biopolymers exist in a plasticized state, whether it is naturally with water or unnaturally with glycerol or other suitable polyol, to make a flexible material. We have found that the extent to which a biopolymer can be plasticized is dependent on its molecular and higher order structures outside of simply molecular weight. Lactalbumin, ovalbumin, corn zein, wheat gluten, and feather keratin were plasticized with glycerol from very low to very high amounts. The conformation of the proteins was monitored with Fourier transform-infrared (FT-IR) spectroscopy and X-ray powder diffraction (XRD) and correlated with the tensile modulus. Protein conformational changes were pronounced for polar proteins with a low amount of cysteine. FT-IR showed that the conformational changes resulted in ordering of the protein at low to moderate plasticization levels. For proteins with little resistance to conformational changes, additional small-scale ordering occurred around the glass transition, as observed in XRD. Accurate comparison of plasticized proteins was dependent on knowing whether or not the protein was glassy or rubbery at room temperature as no differences arose in the glassy state. The transition from glassy to rubbery behavior with plasticization level can be found from modulus, FT-IR, and XRD data.     

Effects of medium polarization and pre-existing field on activation energy of enzymatic charge-transfer reactions

The highly organized spatial structure of proteins' polar groups results in the existence of a permanent intraprotein electric field and in protein's weak dielectric response, i.e. its low dielectric constant. The first factor affects equilibrium free energy gap of a charge-transfer reaction, the second (medium polarization effect) influences both equilibrium and non-equilibrium (reorganization) energies, decreasing the latter substantially. In the framework of the rigorous 'fixed-charge-density' formalism, the medium polarization component of the reaction activation energy has been calculated, both for the activation energy of the elementary act proper, and the effective activation energy accounting for the charges' transfer from water into a low-dielectric structureless medium. In all typical cases of reactions, the energy spent for charge transfer from water into structureless 'protein' is larger than the gain in activation energy due to the protein's low reorganization energy. Therefore, the low dielectric constant of proteins is not sufficient to ensure their high catalytic activity, and an additional effect of the pre-existing intraprotein electric field, compensating for an excessive charging energy, is necessary. Only a combined action of low reorganization energy and pre-existing electric field provides proteins with their high catalytic activity. The dependence of activation energy on the globule geometry has been analyzed. It is shown that, for each reaction, an optimum set of geometric parameters exists. For five hydrolytic enzymes, the optimum globule radii have been calculated using the experimental geometry of their active sites. The calculated radii agree satisfactorily with the real sizes of these macromolecules, both by absolute and by relative values.     

Effective activation energy of enzymatic and nonenzymatic reactions. Evolution-imposed requirements to enzyme structure

The difference of the activation energies in a protein globule and water has been treated in terms of the theory of an elementary act of charge transfer reaction with regards to the energy spent on the transfer of charged reactants from water into the protein. The protein was treated as a structureless dielectric with a given optical and static dielectric constants surrounded by the aqueous phase. Reactions of different types (charge exchange between reactants, charge separation, neutralization, etc.) have been analyzed both under prevalence of purely electrostatic effects and under considerable nonelectrostatic contributions to the activation energies. It is shown that for all one-electron and most multi-electron reactions involving two reaction centres the energy spent for charged reactant transfer from water into protein is greater than the concomitant activation energy gain. The same effect takes place in a number of cases for multi-centre processes as well. To overcome the entropy hindrances, the reactants and catalysts must combine into multiparticle complexes, i.e. form microscopic regions of low dielectric constant. This results in increased effective activation energy as compared to reactions in water. It has been hypothesized that in order to make up for this loss the evolution has selected the proteins which are characterized by considerable intraglobular permanent electric fields. The presence in proteins of high concentrations of strongly polar peptide groups renders them advantageous in this respect over other polymers that are less polar.     

Conformational change that accompanies pepsinogen activation observed in real time by fluorescence energy transfer

Pig pepsinogen has been reacted with N-carboxymethylisatoic anhydride to form N-carboxymethyl-anthraniloyl-(CMA-) pepsinogen, derivatized at Lysp18, Lysp23, Lysp27, Lysp30, and Lys320. Conformational change associated with activation was detected by following energy transfer from tryptophan residues of the pepsin moiety, excited at 295 nm, to CMA groups, monitored by emission above 415 nm. Efficiency of this energy transfer is a measure of conformational change. For this zymogen derivative the change in efficiency occurs with a first order rate constant of 0.041 s-1 at pH 2.4, 22 degrees, which equals the rate at which, following acidification, alkali-stable potential activity becomes alkali-labile. For the native zymogen the rate of this conversion had been shown to be identical to the rate of cleavage of the scissile bond of pepsinogen. Therefore, the correspondence in this derivative of the rates of conversion to alkali lability and change in energy transfer demonstrates that a conformational change accompanies the peptide bond cleavage of activation.     

Electrostatic effects in proteins: comparison of dielectric and charge models.

Two approaches for calculating electrostatic effects in proteins are compared and ana analysis is presented of the dependence of calculated properties on the model used to define the charge distribution. Changes in electrostatic free energy have been calculated using a screened Coulomb potential (SCP) with a distance-dependent effective dielectric permittivity to model bulk solvent effects and a finite difference approach to solve the Poisson-Boltzmann (FDPB) equation. The properties calculated include shifts in dissociation constants of ionizable groups, the effect of annihilating surface charges on the binding of metals, and shifts in redox potentials due to changes in the charge of ionizable groups. In the proteins considered the charged sites are separated by 3.5-12 A. It is shown that for the systems studied in this distance range the SCP yields calculated values which are at least as accurate as those obtained from solution of the FDPB equation. In addition, in the distance range 3-5 A the SCP gives substantially better results than the FDPB equation. Possible sources of this difference between the two methods are discussed. Shifts in binding constants and redox potentials were calculated with several standard charge sets, and the resulting values show a variation of 20-40% between the 'best' and 'worst' cases. From this study it is concluded that in most applications, changes in electrostatic free energies can be calculated economically and reliably using an SCP approach with a single functional form of the screening function.     

Conformational mobility and enzymatic activity of Ca2+-ATPase from sarcoplasmic reticulum

Purified preparations of Ca2+-dependent ATPase were lipid-deleted and incorporated into egg lecithin (EL) and dipalmitoyl lecithin (DPL) liposomes. The temperature dependences of the catalytic activity and of molecular mobility of the spin label (N-1-hydroxyl-2,2,6,6-tetramethyl-4-piperidyl) maleimide linked to a highly reactive SH-group in the vicinity of the active center (15-16 A) and of the fatty acid spin probe (6-doxylpalmitate) located in the protein-lipid moiety were compared. The molecular mobility of the spin label was measured by the saturation transfer method; that of the spin probe was estimated from the maximal splitting value. It was found that the catalytic activity of DPL is correlated with the molecular mobility of the hydrophobic part of ATPase, while that of EL with the segment flexibility in the vicinity of the active center.     

Electronic coupling and charge transfer of DNA: energy control

The influence of the energetic gap on the effective distance-decay rate of electronic coupling (beta(eff)) in DNA is investigated in the context of the superexchange mechanism. The DNA double helix is described by a tight-binding electronic Hamiltonian model, in which all orbitals have the same energy and interact with one another through an exponentially decaying function of distance. Our numerical results concerning the beta(eff) values obtained for two different DNA molecules are analyzed within the theoretical framework of the "continuous-medium approximation," previously developed by Lopez-Castillo et al. (J.-M. Lopez-Castillo, A. Filali-Mouhim, I.L. Plante, and J.-P. Jay-Gerin. J. Phys. Chem. 99 : 6864-6875, 1995). We find that the intervening DNA bridge between the donor and acceptor sites is defined by a unique dimensionless control parameter gamma/E, where E is the energy of the orbitals of this medium with respect to those of the redox site orbitals (energetic gap) and gamma is the electronic band width of the bridge considered as a continuous medium. In the narrow-band regime, our "through-space" coupling model predicts beta(eff) values that are in good order of magnitude agreement with those calculated by other theoretical approaches as well as with those obtained from experiment. Moreover, under equivalent energetic conditions, the DNA-mediated transfers of holes and electrons differ considerably. This difference depends upon the sign of the parameter gamma/E.     

Conformational regulation of charge recombination reactions in a photosynthetic bacterial reaction center

In bright light the photosynthetic reaction center (RC) of Rhodobacter sphaeroides stabilizes the P(+)(870).Q(-)(A) charge-separated state and thereby minimizes the potentially harmful effects of light saturation. Using X-ray diffraction we report a conformational change that occurs within the cytoplasmic domain of this RC in response to prolonged illumination with bright light. Our observations suggest a novel structural mechanism for the regulation of electron transfer reactions in photosynthesis.     

Isotope effects in the study of enzymatic phosphoryl transfer reactions

CASK, a member of the membrane-associated guanylate kinase (MAGUK) superfamily, binds to the carboxyl-terminus of beta-neurexins on the intracellular side of the presynaptic membrane. The guanylate kinase-like (GUK) domains of MAGUKs lack kinase activities, but might be important for mediating specific protein-protein interaction. By a yeast two-hybrid approach, we identified an interaction between the GUK domain of CASK and the C2B domain of rabphilin3a, a presynaptic protein involved in synaptic vesicle exocytosis. The interaction was confirmed by in vitro GST pull-down and co-immunoprecipitation assays. It was proposed that presynaptic vesicles might be guided to the vicinity of points of exocytosis defined by beta-neurexins via the interaction between rabphilin3a-CASK-beta-neurexins.     

Stereospecificity for nicotinamide nucleotides in enzymatic and chemical hydride transfer reactions

The pyridine nucleotide (NAD and NADP)-linked enzymes are a large class of enzymes constituting approximately 17% of all classified enzymes. When these enzymes catalyze their reactions, the hydride transfer between the substrate and the reaction site (i.e., C-4 of the nicotinamide/dihydronicotinamide ring) of the coenzyme takes place in a stereospecific manner. Thus, in the reaction of oxidation of the reduced coenzyme, one group of enzymes catalyzes the extraction of only the hydrogen having the R configuration at the No. 4 carbon, while the other group catalyzes the removal of only that with the S configuration. Because this aspect of enzyme stereospecificity provides essential information for a given enzyme's reaction mechanism, active site structure, and evolutionary relationship with other enzymes, intensive effort has been made to establish the stereospecificities of as many enzymes as possible. This review presents the compilation of the stereospecificities of these enzymes. Some empirical rules, which are useful but not definitive, in predicting a given enzyme's stereospecificity are also described. In addition, the stereospecificity in enzymatic reactions is compared to the stereo-preference in chemical oxidoreduction of the coenzyme. In order to elucidate the mechanism for the enzyme stereospecificity, the conformations of the coenzyme in free-state and enzyme-bound state are extensively discussed here.     

Conformational change in the activation of lipase: an analysis in terms of low-frequency normal modes.

The interfacial activation of Rhizomucor miehei lipase (RmL) involves the motion of an alpha-helical region (residues 82-96) which acts as a "lid" over the active site of the enzyme, undergoing a displacement from a "closed" to an "open" conformation upon binding of substrate. Normal mode analyses performed in both low and high dielectric media reveal that low-frequency vibrational modes contribute significantly to the conformational transition between the closed and open conformations. In these modes, the lid displacement is coupled to local motions of active site loops as well as global breathing motions. Atomic fluctuations of the first hinge of the lid (residues 83-84) are substantially larger in the low dielectric medium than in the high dielectric medium. Our results also suggest that electrostatic interactions of Arg86 play an important role in terms of both the intrinsic stability of the lid and its displacement, through enhancement of hinge mobility in a high dielectric medium. Additional calculations demonstrate that the observed patterns of atomic fluctuations are an intrinsic feature of the protein structure and not dependent on the nature of specific energy minima.     

Enthalpy array analysis of enzymatic and binding reactions

Enthalpy arrays enable label-free, solution-based calorimetric detection of molecular interactions in a 96-detector array format. The combination of the small size of the detectors and the ability to perform measurements in parallel results in a significant reduction of sample volume and measurement time compared with conventional calorimetry. We have made significant improvements in the technology by reducing the temperature noise of the detectors and improving the fabrication materials and methods. In combination with an automated measurement system, the advances in device performance and data analysis have allowed us to develop basic enzyme assays for substrate specificity and inhibitor activity. We have also performed a full titration of 18-crown-6 with barium chloride. These results point to future applications for enthalpy array technology, including fragment-based screening, secondary assays, and thermodynamic characterization of leads in drug discovery.     

Electron transfer and energy dependent reactions in glutaraldehyde-fixed chloroplasts

The effects of GA fixation on electron transfers in photosystemsI and II in photosynthesis and energy dependent reactions ofchloroplasts, such as changes in light scattering, H⁺ uptakeand 515-nm absorbance, were investigated. Fixation of chloroplastswith GA resulted in a lowering of the DCIP and MV photoreductions.DCIP photoreduction activity in fixed chloroplasts was not restoredin the presence of DPC, an electron donor to photosystem II,but was significantly stimulated by DPC when chloroplasts werefixed after aging. The results suggest that the inhibitory effectof GA fixation on photosystem II differs in its mechanism fromthose of treatments such as heating, Tris-washing and aging.The oxidation-reduction reaction of P700 was depressed by GAfixation. Energy dependent reactions in fixed chloroplasts were more markedlydepressed than were electron transfers. Fixed chloroplasts showeda slight conformational response in the presence of PMS. Analysis of the emission spectrum and the induction of chlorophylla fluorescence in fixed chloroplasts suggested that the twopigment systems were partially disordered and that the correspondingprimary photochemical processes were inhibited. (Received November 21, 1972; )     

Perspectives of data analysis of enzyme inhibition and activation,Part 2: Parametrical classification of types of enzymatic reactions

A parametrical classification of types of enzymatic reactions that counts 15 types—7 inhibited reactions, 7 activated ones, and 1 type of initial (uninhibited and nonactivated) reaction—is considered. The system permits taking into account the number of parameters of enzymatic reactions and symmetricity of a course of their change. © 2009 Wiley Periodicals, Inc. J Biochem Mol Toxicol 23:101–107, 2009; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/jbt.20272     

Excitation energy transfer and charge separation in photosystem II membranes revisited

We have performed time-resolved fluorescence measurements on photosystem II (PSII) containing membranes (BBY particles) from spinach with open reaction centers. The decay kinetics can be fitted with two main decay components with an average decay time of 150 ps. Comparison with recent kinetic exciton annihilation data on the major light-harvesting complex of PSII (LHCII) suggests that excitation diffusion within the antenna contributes significantly to the overall charge separation time in PSII, which disagrees with previously proposed trap-limited models. To establish to which extent excitation diffusion contributes to the overall charge separation time, we propose a simple coarse-grained method, based on the supramolecular organization of PSII and LHCII in grana membranes, to model the energy migration and charge separation processes in PSII simultaneously in a transparent way. All simulations have in common that the charge separation is fast and nearly irreversible, corresponding to a significant drop in free energy upon primary charge separation, and that in PSII membranes energy migration imposes a larger kinetic barrier for the overall process than primary charge separation.     

Effects of polyelectrolyte chain stiffness, charge mobility, and charge sequences on binding to proteins and micelles

The binding affinities of polyanions for bovine serum albumin in NaCl solutions from I = 0.01-0.6 M, were evaluated on the basis of the pH at the point of incipient binding, converting each such pH(c) value into a critical protein charge Zc. Analogous values of critical charge for mixed micelles were obtained as the cationic surfactant mole fraction Yc. The data were well fitted as Yc or Zc = KI a, and values of K and a were considered as a function of normalized polymer charge densities (tau), charge mobility, and chain stiffness. Binding increased with chain flexibility and charge mobility, as expected from simulations and theory. Complex effects of tau were related to intrapolyanion repulsions within micelle-bound loops (seen in the simulations) or negative protein domain-polyanion repulsions. The linearity of Zc with radicalI at I < 0.3 M was explained by using protein electrostatic images, showing that Zc at I < 0.3 M depends on a single positive "patch"; the appearance of multiple positive domains I > 0.3 M (lower pH(c)) disrupts this simple behavior.     

Conformational analysis of lipid-associating proteins in a lipid environment

Two major types of helical structures have been identified in lipid-associating proteins, being either amphipathic or transmembrane domains. A conformational analysis was carried out to characterize some of the properties of these helices. These calculations were performed both on isolated helices and in a lipid environment. According to the results of this analysis, the orientation of the line joining the hydrophobic and hydrophilic centers of the helix seems to determine the orientation of the helix at the lipid/water interface. The calculation of this parameter should be useful to discriminate between an amphipathic helix, parallel to the interface and a transmembrane helix orientated perpendicularly. The membrane-spanning helices are completely immersed in the phospholipid bilayer and their length corresponds to about the thickness of the hydrophobic core of the DPPC bilayer. The energy of interaction, expressed per phospholipid is significantly higher for the transmembrane compared to the amphipathic helices. For the membrane-spanning helices the mean energy of interaction is higher than the interaction energy between two phospholipids, while it is lower for most amphipathic helices. This might account for the stability of these protein-anchoring domains. This computer modeling approach should usefully complement the statistical analysis carried out on these helices, based on their hydrophobicity and hydrophobic moment. It represents a more refined analysis of the domains identified by the prediction techniques and stress the functional character of lipid-associating domains in membrane proteins as well as in soluble plasma lipoproteins.