The properties of bio-energy transport and influence of structure nonuniformity and temperature of systems on energy transport along polypeptide chains

A new theory of bio-energy transport along protein molecules in living systems, where the energy is released by hydrolysis of adenosine triphosphate (ATP), is proposed based on some physical and biological reasons. In the new theory, the Davydov's Hamiltonian and wave function of the systems are simultaneously modified and extended. A new interaction has been added into Davydov's Hamiltonian. The wave function of the excitation state of single particles for the excitons in the Davydov model is replaced by a new wave function of two-quanta quasicoherent state. In such a case, the bio-energy is transported by the new soliton, which differs from the Davydov's soliton. The soliton is formed through self- trapping of two excitons interacting amino acid residues. The exciton is generated by vibrations of amide-I (CO stretching) arising from the energy of hydrolysis of ATP. The properties of the new soliton are extensively studied by analytical method and its lifetime is calculated using the nonlinear quantum perturbation theory and a wide ranges of parameter values relevant to protein molecules. The lifetime of the new soliton at the biological temperature 300 K is enough large and belongs to the order of 10?1? s, or τ/τ?≥700, in which the soliton can transports over several hundreds amino acid residues. These studied results show clearly that the new soliton is thermally stable and has so larger lifetime that it can play an important role in biological processes. Thus the new model is a candidate of the bio-energy transport mechanism in protein molecules. In the meanwhile, the influences of structure nonuniformity in protein molecules and temperature of the systems on the states and properties of the soliton transport of bio-energy are numerically simulated and studied by the fourth-order Runge-Kutta method. The structure nonuniformity arises from the disorder distributions of masses of amino acid residues, side groups and impurities, which results also in the fluctuations of the spring constant of protein molecules, dipole-dipole interaction between the neighboring amides, exciton-phonon (vibration of amino acids)interaction, chain-chain interaction among the three channels and ground state energy of the systems. We investigated the behaviors and states of the new solitons in a single protein molecular chain and α-Helix protein molecules with three channels under influences of the structure nonuniformity. We prove first that the bio-energy is transported by a soliton, which can move without dispersion, retaining its shape, velocity and energy in a uniform and periodic protein molecule. When the structure nonuniformity exists, although the fluctuations of the spring constant, dipole-dipole interaction constant, exciton-phonon coupling constant and ground state energy and the nonuniformity distributions of masses of amino acid residues can change the states and properties of motion of new soliton, they are still quite stable and very robust against these structure nonuniformities, i.e., even there are a larger structure nonuniformity in the protein molecules, the new solitons cannot be still dispersed. If the effects of thermal perturbation of medium on the soliton in nonuniform proteins is considered again, the new soliton can transport also over a larger spacing of 400 amino acids and has a longer time period of 300 ps, it is still thermally stable up to 320 K under the influence of the above structure nonuniformities. However, the new soliton disperses in the case of a higher temperature of 325 K and in more large structure nonuniformity. Thus, we determine that the new soliton's lifetime and critical temperature are 300 ps and 320 K, respectively. These results are also consistent with analytical data obtained via quantum perturbed theory. For α-Helix protein molecules with three channels, the results obtained show that the structure nonuniformity and quantum fluctuation can change the states and features of the new solitons, for example, the amplitudes, energies and velocities of the new soliton are decreased, but the solitons have been not destroyed, they can still transport steadily along the molecular chains retaining energy and momentum. When the quantum fluctuations are larger, such as, structure disorders and quantum fluctuations of 0.67<α(K)<2, ΔW=±8%Wˉ, ΔJ=±1%Jˉ, Δ(χ?+χ?)=±3%(χˉ?+χˉ?) and ΔL=±1%Lˉ and Δ??=?|β(n)|, ?=0.1 meV, |β(n)|<0.5, the new soliton is still stable. Therefore, the new solitons are quite robust against these nonuniform effects. However, they will be dispersed or disrupted in cases of very large structure nonuniformity. When the influence of temperature on solitons is considered, we find that the new solitons can transport steadily over 333 amino acid residues in the case of a long time period of 120 ps, in which the soliton can retain its shape and energy to travel forward along protein molecules after their mutual collision at the biological temperature of 300 K. However, the soliton disperses in cases of higher temperatures 325 K under action of a larger structure disorder. Thus, its critical temperature is about 320 K. When the effects of structure nonuniformity and temperature are considered simultaneously, then the new soliton has still high thermal stability and can transport also along the protein molecular chains retaining its amplitude, energy and velocity, they will disperses in the larger fluctuations, for example, 0.67 Mˉ相似文献   

Nonlinear mechanism for weak photon emission from biosystems

The nonlinear mechanism for the origin of the weak biophoton emission from biological systems is suggested. The mechanism is based on the properties of solitons that provide energy transfer and charge transport in metabolic processes. Such soliton states are formed in alpha-helical proteins. Account of the electron-phonon interaction in macromolecules results in the self-trapping of electrons in a localized soliton-like state, known as Davydov's solitons. The important role of the helical symmetry of macromolecules is elucidated for the formation, stability and dynamical properties of solitons. It is shown that the soliton with the lowest energy has an inner structure with the many-hump envelope. The total probability of the excitation in the helix is characterized by interspine oscillations with the frequency of oscillations, proportional to the soliton velocity. The radiative life-time of a soliton is calculated and shown to exceed the life-time of an excitation on an isolated peptide group by several orders of magnitude.     

The nature of phonons and solitary waves in alpha-helical proteins.

A parametric study of the Davydov model of energy transduction in alpha-helical proteins is described. Previous investigations have shown that the Davydov model predicts that nonlinear interactions between phonons and amide-I excitations can stabilize the latter and produce a long-lived combined excitation (the so-called Davydov soliton), which propagates along the helix. The dynamics of this solitary wave are approximately those of solitons described using the nonlinear Schr?dinger equation. The present study extends these previous investigations by analyzing the effect of helix length and nonlinear coupling efficiency on the phonon spectrum in short and medium length alpha-helical segments. The phonon energy accompanying amide-I excitation shows periodic variation in time with fluctuations that follow three different time scales. The phonon spectrum is highly dependent upon chain length but a majority of the energy remains localized in normal mode vibrations even in the long chain alpha-helices. Variation of the phonon-exciton coupling coefficient changes the amplitudes but not the frequencies of the phonon spectrum. The computed spectra contain frequencies ranging from 200 GHz to 6 THz, and as the chain length is increased, the long period oscillations increase in amplitude. The most important prediction of this study, however, is that the dynamics predicted by the numerical calculations have more in common with dynamics described by using the Frohlich polaron model than by using the Davydov soliton. Accordingly, the relevance of the Davydov soliton model was applied to energy transduction in alpha-helical proteins is questionable.(ABSTRACT TRUNCATED AT 250 WORDS)     

电离辐射对蛋白质分子作用的量子孤波分析

本文运用量子力学方法对蛋白质分子中孤波传播的非线性动力学特征进行了探讨。研究表明：电离辐射产生的自由基对蛋白质分子的伤害将会对携带能量、信息的孤立子波传播产生较为显著的影响。     

Catalytic role of proton transfers in the formation of a tetrahedral adduct in a serine carboxyl peptidase

Quantum mechanical/molecular mechanical molecular dynamics and 2D free energy simulations are performed to study the formation of a tetrahedral adduct by an inhibitor N-acetyl-isoleucyl-prolyl-phenylalaninal (AcIPF) in a serine-carboxyl peptidase (kumamolisin-As) and elucidate the role of proton transfers during the nucleophilic attack by the Ser278 catalytic residue. It is shown that although the serine-carboxyl peptidases have a fold resembling that of subtilisin, the proton transfer processes during the nucleophilic attack by the Ser residue are likely to be more complex for these enzymes compared to the case in classical serine proteases. The computer simulations demonstrate that both general base and acid catalysts are required for the formation and stabilization of the tetrahedral adduct. The 2D free energy maps further demonstrate that the proton transfer from Ser278 to Glu78 (the general base catalyst) is synchronous with the nucleophilic attack, whereas the proton transfer from Asp164 (the general acid catalyst) to the inhibitor is not. The dynamics of the protons at the active site in different stages of the nucleophilic attack as well as the motions of the corresponding functional groups are also studied. It is found that the side chain of Glu78 is generally rather flexible, consistent with its possible multifunctional role during catalysis. The effects of proton shuffling from Asp82 to Glu78 and from Glu32 to Asp82 are examined, and the results indicate that such proton shuffling may not play an important role in the stabilization of the tetrahedral intermediate analogue.     

Ab initio molecular dynamics study of proton transfer in a polyglycine analog of the ion channel gramicidin A.

Proton transfer in biological systems is thought to often proceed through hydrogen-bonded chains of water molecules. The ion channel, gramicidin A (gA), houses within its helical structure just such a chain. Using the density functional theory based ab initio molecular dynamics Car-Parrinello method, the structure and dynamics of proton diffusion through a polyglycine analog of the gA ion channel has been investigated. In the channel, a proton, which is initially present as hydronium (H3O+), rapidly forms a strong hydrogen bond with a nearest neighbor water, yielding a transient H5O2+ complex. As in bulk water, strong hydrogen bonding of this complex to a second neighbor solvation shell is required for proton transfer to occur. Within gA, this second neighbor shell included not only a channel water molecule but also a carbonyl of the channel backbone. The present calculations suggest a transport mechanism in which a priori carbonyl solvation is a requirement for proton transfer.     

A QM/MM study of proton transport pathways in a [NiFe] hydrogenase

A theoretical QM/MM study of the [NiFe] hydrogenase from Desulfovibrio fructosovorans has been performed to investigate possible routes of proton transfer between the active site and the protein surface. We obtained the minimum energy paths, with a modified version of the nudged elastic band method, for a set of proposed pathways. The calculations were carried out for the crystallographic structure and for several structures of the protein obtained from a molecular dynamics simulation. The results show one of the studied pathways to be preferred for transport from the active site to the surface, but the preference is not so strong when transport occurs in the opposite direction.     

Initial stage of coupling electron and proton transport in the vicinity of the secondary quinone in photosynthesis of bacteria and plants

The coupling of electron and proton transport in the vicinity of the secondary quinone QB in the reaction center of bacteria and photosystem II of higher plants was investigated. The energy levels and wave functions of the proton in the system QB--histidine L 190 were calculated. It was shown that the proton of histidine forms a hydrogen bond with the doubly reduced quinone QB2-. A new scheme of proton transport through histidine L 190 and its coupling with electron transport was proposed.     

Transport properties of membrane vesicles fromAcholeplasma laidlawii

Membrane vesicles obtained from Acholeplasma laidlawii accumulate glucose as well as maltose and fructose against their concentration gradient in the absence of exogenous energy sources. Glucose uptake by membrane vesicles is inhibited by anaerobiosis and by electron transfer inhibitors, such as rotenone and amytal, but not by 2-heptyl-4-hydroxyquinoline N-oxide, antimycin A, cyanide and azide. Rotenone, cyanide and amytal also produce a rapid efflux of glucose from the membrane vesicles. Arsenate, oligomycin and N,N'-dicyclohexylcarbodimide do not inhibit glucose transport. Transport of glucose is markedly inhibited by proton conductors such as CCCP and pentachlorophenol. It is concluded that glucose transport can be driven by a high-energy state of the membrane or by the membrane potential.     

Dynamics of protein conformational fluctuation in enzyme catalysis with special attention to proton transfers in serine proteinases

We have provided a quantum mechanical model for proteinase-catalyzed peptide, amide and ester hydrolysis. The model rests on electron and atom transfer theory, but incorporates the dynamics of conformational nuclear modes as a new element. The model is applied to acylation, but can straightaway be extended to deacylation, and is substantiated by recent structural and kinetic data for proteinase enzyme catalysis. The role of the conformational modes is found to be two-fold. First, the crystallographic distances for the proton transfers involved are far too large for direct transfer. His-57 mobility, handled stochastically, to bring the donor and acceptor groups within suitable reach, is therefore a crucial element of the theory. Secondly, the charge alignment in the Asp-102/His-57/tetrahedral intermediate system implies that the curvature of the potential surface along the conformational coordinates in this state is much lower than in the initial enzyme-substrate and final acyl states. A consequence of this is that the activation energy liberated after the first proton transfer is not dissipated, but stored in the conformational system and used in the second proton transfer step.     

Hydration of CO2 by carbonic anhydrase: intramolecular proton transfer between Zn2+-bound H2O and histidine 64 in human carbonic anhydrase II

The energy barrier for the intramolecular proton transfer between zinc-bound water and His 64 in the active site of human carbonic anhydrase II (HCA II) has been studied at the partial retention of diatomic differential overlap (PRDDO) level. The most important stabilizing factor for the intramolecular proton transfer is the zinc ion, which lowers the pKa of zinc-bound water and electrostatically repels the proton. The energy barrier of 127.5 kcal/mol for proton transfer between a water dimer is completely removed in the presence of the zinc ion. The zinc ligands, which donate electrons to the zinc ion, raise the barrier slightly to 34 kcal/mol for a 4-coordinated zinc complex including three imidazole ligands from His 94, His 96, and His 119 and to 54 kcal/mol for the 5-coordinated zinc complex including the fifth water ligand. A few model calculations indicate that these energy barriers are expected to be reduced to within experimental range (approximately 10 kcal/mol) when large basis set, correlation energies, and molecular dynamics are considered. The proton-transfer group, which functions as proton receiver in the intramolecular proton transfer, helps to attract the proton; and the partially ordered active site water molecules are important for proton relay function.     

Green fluorescent protein variants as ratiometric dual emission pH sensors. 2. Excited-state dynamics

In the preceding paper [Hanson, G. T., McAnaney, T. B., Park, E. S., Rendell, M. E. P., Yarbrough, D. K., Chu, S., Xi, L., Boxer, S. G., Montrose, M. H., and Remington, S. J. (2002) Biochemistry 41, 15477-15488], novel mutants of the green fluorescent protein (GFP) that exhibit dual steady-state emission properties were characterized structurally and discussed as potential intracellular pH probes. In this work, the excited-state dynamics of one of these new dual emission GFP variants, deGFP4 (C48S/S65T/H148C/T203C), is studied by ultrafast fluorescence upconversion spectroscopy. Following excitation of the high-energy absorption band centered at 398 nm and assigned to the neutral form of the chromophore, time-resolved emission was monitored from the excited state of both the neutral and intermediate anionic chromophores at both high and low pH and upon deuteration of exchangeable protons. The time-resolved emission dynamics and isotope effect appear to be very different from those of wild-type GFP [Chattoraj, M., King, B. A., Bublitz, G. U., and Boxer, S. G. (1996) Proc. Natl. Acad. Sci. U.S.A. 93, 8362-8367]; however, due to overlapping emission bands, the apparent difference can be analyzed quantitatively within the same framework used to describe GFP excited-state dynamics. The results indicate that the pH-sensitive steady-state emission characteristics of deGFP4 are a result of a pH-dependent modulation of the rate of excited-state proton transfer. At high pH, a rapid interconversion from the excited state of the higher energy neutral chromophore to the lower energy intermediate anionic chromophore is achieved by proton transfer. At low pH, excited-state proton transfer is slowed to the point where it is no longer rate limiting.     

Charge transfer in green fluorescent protein.

Charge transfer reactions that contribute to the photoreactions of the wild type green fluorescent protein (GFP) do not occur in the isolated p-hydroxybenzylidene-imidazolidinone chromophore, demonstrating the role of the protein environment. The high quantum efficiency of the fluorescence photocycle that includes excited state proton transfer and the suppression of non-radiative pathways by the protein environment have been correlated with structural dynamics in the chromophore environment. A low quantum efficiency competing phototransformation reaction of GFP is accompanied by both proton and electron transfer, and closely mimics the charge redistribution that is occurring in the fluorescence photocycle. The protein response to this destabilising event has been demonstrated by cryo-trapping of early products in the reaction pathway and is found to be strong even at 100 K, including displacements of chromophore, protein, solvent and a photogenerated CO2 molecule derived from the decarboxylated Glu 222 side chain. We discuss the ramifications of the observation of strong conformational perturbations below the protein dynamical transition at approximately 200 K, in view of low temperature work on other light sensitive proteins such as myoglobin and bacteriorhodopsin. The proton and electron transfer in the phototransformation pathway mimics the proton and charge transfer which occurs during the fluorescence cycle, which leads to common structural responses in both photoreactions as shown by ultrafast spectroscopy. We review and discuss literature on light-induced and thermal charge transfer events, focusing on recent findings addressing conformational dynamics and implications for thermodynamic properties.     

Biophysical aspects of intra-protein proton transfer

The passage of proton trough proteins is common to all membranal energy conserving enzymes. While the routes differ among the various proteins, the mechanism of proton propagation is based on the same chemical-physical principles. The proton progresses through a sequence of dissociation association steps where the protein and water molecules function as a solvent that lowers the energy penalty associated with the generation of ions in the protein. The propagation of the proton in the protein is a random walk, between the temporary proton binding sites that make the conducting path, that is biased by the intra-protein electrostatic potential. Kinetic measurements of proton transfer reactions, in the sub-ns up to micros time frame, allow to monitor the dynamics of the partial reactions of an overall proton transfer through a protein.     

Light harvesting in photosystem I: modeling based on the 2.5-A structure of photosystem I from Synechococcus elongatus

The structure of photosystem I from the thermophilic cyanobacterium Synechococcus elongatus has been recently resolved by x-ray crystallography to 2.5-A resolution. Besides the reaction center, photosystem I consists also of a core antenna containing 90 chlorophyll and 22 carotenoid molecules. It is their function to harvest solar energy and to transfer this energy to the reaction center (RC) where the excitation energy is converted into a charge separated state. Methods of steady-state optical spectroscopy such as absorption, linear, and circular dichroism have been applied to obtain information on the spectral properties of the complex, whereas transient absorption and fluorescence studies reported in the literature provide information on the dynamics of the excitation energy transfer. On the basis of the structure, the spectral properties and the energy transfer kinetics are simultaneously modeled by application of excitonic coupling theory to reveal relationships between structure and function. A spectral assignment of the 96 chlorophylls is suggested that allows us to reproduce both optical spectra and transfer and emission spectra and lifetimes of the photosystem I complex from S. elongatus. The model calculation allowed to study the influence of the following parameters on the excited state dynamics: the orientation factor, the heterogeneous site energies, the modifications arising from excitonic coupling (redistribution of oscillator strength, energetic splitting, reorientation of transition dipoles), and presence or absence of the linker cluster chlorophylls between antenna and reaction center. For the F?rster radius and the intrinsic primary charge separation rate, the following values have been obtained: R(0) = 7.8 nm and k(CS) = 0.9 ps(-1). Variations of these parameters indicate that the excited state dynamics is neither pure trap limited, nor pure transfer (to-the-trap) limited but seems to be rather balanced.     

Insight into the phosphodiesterase mechanism from combined QM/MM free energy simulations

Molecular dynamics simulations employing a combined quantum mechanical and molecular mechanical potential have been carried out to elucidate the reaction mechanism of the hydrolysis of a cyclic nucleotide cAMP substrate by phosphodiesterase 4B (PDE4B). PDE4B is a member of the PDE superfamily of enzymes that play crucial roles in cellular signal transduction. We have determined a two-dimensional potential of mean force (PMF) for the coupled phosphoryl bond cleavage and proton transfer through a general acid catalysis mechanism in PDE4B. The results indicate that the ring-opening process takes place through an S(N)2 reaction mechanism, followed by a proton transfer to stabilize the leaving group. The computed free energy of activation for the PDE4B-catalyzed cAMP hydrolysis is about 13 kcal·mol(-1) and an overall reaction free energy is about -17 kcal·mol(-1), both in accord with experimental results. In comparison with the uncatalyzed reaction in water, the enzyme PDE4B provides a strong stabilization of the transition state, lowering the free energy barrier by 14 kcal·mol(-1). We found that the proton transfer from the general acid residue His234 to the O3' oxyanion of the ribosyl leaving group lags behind the nucleophilic attack, resulting in a shallow minimum on the free energy surface. A key contributing factor to transition state stabilization is the elongation of the distance between the divalent metal ions Zn(2+) and Mg(2+) in the active site as the reaction proceeds from the Michaelis complex to the transition state.     

Cysteine biosynthetic enzymes are the pieces of a metabolic energy pump

Understanding the mechanisms of free energy transfer in metabolism is fundamental to understanding how the chemical forces that sustain the molecular organization of the cell are distributed. Recent studies of molecular motors (1-3) and ATP-driven proton transport (4-6) describe how chemical potential is transferred at the molecular level. These systems catalyze energy transfer through structural change and appear to be dedicated exclusively to their coupling tasks (7, 8). Here we report the discovery of a new class of energy-transfer system. It is a biosynthetic pump composed of cysteine biosynthesis enzymes, ATP sulfurylase and O-acetylserine sulfhydrylase, each with its own catalytic function and from whose interactions emerge new function: the hydrolysis of ATP. The hydrolysis is kinetically and energetically linked to the chemistry catalyzed by ATP sulfurylase, the first enzyme in the cysteine biosynthetic pathway, in such a way that each molecule of ATP hydrolyzed, each stroke of the pump, produces 1 equivalent of that enzyme's product. These findings integrate cysteine metabolism and broaden our understanding of the ways in which higher order allostery is used to effect free energy transfer.