Energetics and kinetics of primary charge separation in bacterial photosynthesis

We report the results of molecular dynamics (MD) simulations and formal modeling of the free-energy surfaces and reaction rates of primary charge separation in the reaction center of Rhodobacter sphaeroides. Two simulation protocols were used to produce MD trajectories. Standard force-field potentials were employed in the first protocol. In the second protocol, the special pair was made polarizable to reproduce a high polarizability of its photoexcited state observed by Stark spectroscopy. The charge distribution between covalent and charge-transfer states of the special pair was dynamically adjusted during the simulation run. We found from both protocols that the breadth of electrostatic fluctuations of the protein/water environment far exceeds previous estimates, resulting in about 1.6 eV reorganization energy of electron transfer in the first protocol and 2.5 eV in the second protocol. Most of these electrostatic fluctuations become dynamically frozen on the time scale of primary charge separation, resulting in much smaller solvation contributions to the activation barrier. While water dominates solvation thermodynamics on long observation times, protein emerges as the major thermal bath coupled to electron transfer on the picosecond time of the reaction. Marcus parabolas were obtained for the free-energy surfaces of electron transfer by using the first protocol, while a highly asymmetric surface was obtained in the second protocol. A nonergodic formulation of the diffusion-reaction electron-transfer kinetics has allowed us to reproduce the experimental results for both the temperature dependence of the rate and the nonexponential decay of the population of the photoexcited special pair.     

A secretory phospholipase A2-mediated neuroprotection and anti-apoptosis

Background  

Phospholipase A₂ liberates free fatty acids and lysophospholipids upon hydrolysis of phospholipids and these products are often associated with detrimental effects such as inflammation and cerebral ischemia. The neuroprotective effect of neutral phospholipase from snake venom has been investigated.     

Glassy protein dynamics and gigantic solvent reorganization energy of plastocyanin

We report the results of molecular dynamics simulations of electron-transfer activation parameters of plastocyanin metalloprotein involved as an electron carrier in natural photosynthesis. We have discovered that slow, non-ergodic conformational fluctuations of the protein, coupled to hydrating water, result in a very broad distribution of donor-acceptor energy gaps far exceeding those observed for commonly studied inorganic and organic donor-acceptor complexes. The Stokes shift is not affected by these fluctuations and can be calculated from solvation models in terms of the linear response of the solvent dipolar polarization. The non-ergodic character of large-amplitude protein/water mobility breaks the strong link between the Stokes shift and the reorganization energy characteristic of equilibrium (ergodic) theories of electron transfer. This mechanism might be responsible for fast electronic transitions in natural electron-transfer proteins characterized by low reaction free energy.     

Redox entropy of plastocyanin: developing a microscopic view of mesoscopic polar solvation

We report applications of analytical formalisms and molecular dynamics (MD) simulations to the calculation of redox entropy of plastocyanin metalloprotein in aqueous solution. The goal of our analysis is to establish critical components of the theory required to describe polar solvation at the mesoscopic scale. The analytical techniques include a microscopic formalism based on structure factors of the solvent dipolar orientations and density and continuum dielectric theories. The microscopic theory employs the atomistic structure of the protein with force-field atomic charges and solvent structure factors obtained from separate MD simulations of the homogeneous solvent. The MD simulations provide linear response solvation free energies and reorganization energies of electron transfer in the temperature range of 280-310 K. We found that continuum models universally underestimate solvation entropies, and a more favorable agreement is reported between the microscopic calculations and MD simulations. The analysis of simulations also suggests that difficulties of extending standard formalisms to protein solvation are related to the inhomogeneous structure of the solvation shell at the protein-water interface combining islands of highly structured water around ionized residues along with partial dewetting of hydrophobic patches. Quantitative theories of electrostatic protein hydration need to incorporate realistic density profile of water at the protein-water interface.