In-silico assessment of protein-protein electron transfer. a case study: cytochrome c peroxidase--cytochrome c

The fast development of software and hardware is notably helping in closing the gap between macroscopic and microscopic data. Using a novel theoretical strategy combining molecular dynamics simulations, conformational clustering, ab-initio quantum mechanics and electronic coupling calculations, we show how computational methodologies are mature enough to provide accurate atomistic details into the mechanism of electron transfer (ET) processes in complex protein systems, known to be a significant challenge. We performed a quantitative study of the ET between Cytochrome c Peroxidase and its redox partner Cytochrome c. Our results confirm the ET mechanism as hole transfer (HT) through residues Ala194, Ala193, Gly192 and Trp191 of CcP. Furthermore, our findings indicate the fine evolution of the enzyme to approach an elevated turnover rate of 5.47×10⁶ s⁻¹ for the ET between Cytc and CcP through establishment of a localized bridge state in Trp191.

We have developed a protocol capable of describing long-range electron transfer mechanisms at an atomic detailed level. We demonstrate the maturity of the computational techniques in obtaining a quantitative view of the Cytochrome c Peroxidase/Cytochrome c electron transfer process, known to be a significant challenge. In excellent agreement with experimental data, our results allow for the description of the electron transfer pathway, its mechanism and the electron transfer rate at a quantitative level. The overall protocol is free of parameterization and can be applied to any complex electron transfer process. Furthermore, the results reveal the fine enzyme evolution of this protein-protein complex to optimize its electron transfer rate by a localized bridge state.

DNA base pair stacks with high electric conductance: a systematic structural search

We report a computational search for DNA π-stack structures exhibiting high electric conductance in the hopping regime, based on the INDO/S calculations of electronic coupling and the method of data analysis called k-means clustering. Using homogeneous poly(G)-poly(C) and poly(A)-poly(T) stacks as the simplest structural models, we identify the configurations of neighboring G:C and A:T pairs that allow strong electronic coupling and, therefore, molecular electric conductance much larger than the values reported for the corresponding reference systems in the literature. A computational approach for modeling the impact of thermal fluctuations on the averaged dimer structure was also proposed and applied to the [(G:C),(G:C)] and [(A:T),(A:T)] duplexes. The results of this work may provide guidance for the construction of DNA devices and DNA-based elements of nanoscale molecular circuits. Several factors that cause changes of step parameters favorable to the formation of the predicted stack conformation with high electric conductance of DNA molecules are also discussed; favorable geometries may enhance the conductivity by factors as large as 15.

Electronic coupling for charge transfer in donor-bridge-acceptor systems. Performance of the two-state FCD model

Electronic coupling is a key parameter that determines the rate of electron transfer reactions and electrical conductivity of molecular wires. To examine the performance of a two-state approach based on the orthogonal transformation of adiabatic states to diabatic states, we compare the effective donor-acceptor coupling V(DA) computed with three different approaches in model donor-bridge-acceptor (D-B-A) systems. It is found that V(DA) derived with the two-state method accounts properly for both the direct and superexchange interactions. The approach becomes, however, less accurate with the increasing energy difference of the donor and acceptor states. We suggest a simple diagnostic to identify the situation when the estimated coupling might be inaccurate and consider how to improve the performance of the two-state scheme in such a case.

An in silico design for a DNA nanomechanical switch

We have calculated how the charge transfer properties of DNA change in response to the application of an external stretching force. Since charge transfer occurs through the DNA π-stack, any disruption to this stacking causes dramatic changes in the transport properties of the biomolecule, as our calculations demonstrate. We therefore propose that the mechanical response of DNA to an applied stretching force might be used in the design of a nanomechanical switch.

Can charge transfer in DNA significantly be modulated by varying the pi stack conformation?

DNA is an ideal target for single-molecule manipulations. The conformation of stacked basepairs in DNA depends sensitively on various factors such as temperature, the kind of solvent and counterions, changes in the backbone, applied forces, etc. This raises the question of whether the rate of charge transfer (CT) through the stack can be considerably enhanced by tuning of "observed" DNA conformations. Using a stochastic approach to account for the effects of thermal fluctuations, we study how the efficiency of CT in poly(dA)-poly(dT) and poly(dG)-poly(dC) sequences will change by variation of the pi stack structure. The CT process is shown to be not very sensitive to the torsional angle (twist) while affected more strongly by altering translation modes (shift and slide). We conclude that the design of basepair stacks with significantly improved electrical conductivity, as compared to poly(dG)-poly(dC), appears to be quite elusive. Specific changes of the pi stack structure can increase the efficiency of CT by a factor of approximately 3 (e.g., in systems with negative shift); much stronger effects can hardly be expected. This result, formally derived for molecular ensembles, should also be applicable for single-molecule systems because of the strong effects of dynamical disordering.

Stabilization of radical anion states of nucleobases in DNA

Trapping of an electron by DNA leads to the formation of radical anion states of pyrimidine bases. Because these states play an important role in biological and chemical processes, their computational treatment is of particular interest. We show that simple electrostatic and quantum chemical models can accurately reproduce the adiabatic electron affinities (EAs) of short DNA stacks recently derived from high-level ab initio calculations (M. Kobylecka, J. Leszczynski, and J. Rak, J. Am. Chem. Soc., 2008, 130, 15683). The electrostatic interaction of an excess electron localized on cytosine or thymine with intra- and inter-strand adjacent nucleobases is found to strongly affect the energy of the radical anions. This interaction is the main origin of the dependence of EA of nucleobases on the nature of neighboring base pairs. In particular, the states XT(-)Y and XC(-)Y, where X and Y = C, T, are, by ca. 0.7 eV, more stable than radical anions GT(-)G and GC(-)G. We find that second-neighbor effects can also significantly modulate EAs, although being smaller than the effects of adjacent bases. The strongest destabilizing effect is found for 5'-GC and 3'-GC, while the 5'-AT base pair stabilizes the radical anion states. Using a combined QM/MD approach, we consider how structural fluctuations of DNA influence the stability of the radical anion states. Despite large dispersions of the stabilization energies due to conformational dynamics of DNA, there are only few thermally accessible structures where GT(-)G and GC(-)G are energetically more favorable than the corresponding pyrimidine triplets. Although stabilization energies calculated for stacks of regular structure are in qualitative agreement with the QM/MD results, structural fluctuations of pi stacks should be taken into account for more accurate description of the excess electron trapped by DNA. The results obtained in this study suggest that simple electrostatic models, in combination with MD simulations, can be very helpful to explore the long time scale behavior of radical anions in DNA.