DNA double-helix-mediated long-range electron transfer

A theoretical analysis based upon large-scale self-consistent Hartree-Fock calculations at a semiempirical quantum theory level (CNDO/S) is performed to investigate long-range electron transfer in a donor-DNA-acceptor molecule, where the donor and acceptor moieties are tethered to the DNA. The π-stacked base pairs are found to dominate the long-range electronic coupling. Despite the π-electron mediated coupling, the exponential distance decay constant of the electron transfer rate is ∼ 1.2–1.6 Å⁻¹, values typical of electron transfer proteins. The calculated long-range electron transfer rate of the order of 10⁶ s⁻¹ for a metal-to-metal distance of 21 Å is found to be in agreement with kinetic measurements by Meade and Kayyem. Based on the current analysis, the π-electrons dominate the long-range electronic coupling interactions in DNA, but they do not lead to one-dimensional molecular wire-like properties. © 1996 John Wiley & Sons, Inc.     

Evidence for a near-resonant charge transfer mechanism for double-stranded peptide nucleic acid

We present evidence for a near-resonant mechanism of charge transfer in short peptide nucleic acid (PNA) duplexes obtained through electrochemical, STM break junction (STM-BJ), and computational studies. A seven base pair (7-bp) PNA duplex with the sequence (TA)(3)-(XY)-(TA)(3) was studied, in which XY is a complementary nucleobase pair. The experiments showed that the heterogeneous charge transfer rate constant (k(0)) and the single-molecule conductance (σ) correlate with the oxidation potential of the purine base in the XY base pair. The electrochemical measurements showed that the enhancement of k(0) is independent, within experimental error, of which of the two PNA strands contains the purine base of the XY base pair. 7-bp PNA duplexes with one or two GC base pairs had similar measured k(0) and conductance values. While a simple superexchange model, previously used to rationalize charge transfer in single stranded PNA (Paul et al. J. Am. Chem. Soc. 2009, 131, 6498-6507), describes some of the experimental observations, the model does not explain the absence of an enhancement in the experimental k(0) and σ upon increasing the G content in the duplexes from one to two. Moreover, the superexchange model is not consistent with other studies (Paul et al. J. Phys. Chem. B 2010, 114, 14140), that showed a hopping charge transport mechanism is likely important for PNA duplexes longer than seven base pairs. A quantitative computational analysis shows that a near-resonant charge transfer regime, wherein a mix of superexchange and hopping mechanisms are expected to coexist, can rationalize all of the experimental results.     

NCIPLOT: a program for plotting non-covalent interaction regions

Non-covalent interactions hold the key to understanding many chemical, biological, and technological problems. Describing these non-covalent interactions accurately, including their positions in real space, constitutes a first step in the process of decoupling the complex balance of forces that define non-covalent interactions. Because of the size of macromolecules, the most common approach has been to assign van der Waals interactions (vdW), steric clashes (SC), and hydrogen bonds (HBs) based on pairwise distances between atoms according to their van der Waals radii. We recently developed an alternative perspective, derived from the electronic density: the Non-Covalent Interactions (NCI) index [J. Am. Chem. Soc. 2010, 132, 6498]. This index has the dual advantages of being generally transferable to diverse chemical applications and being very fast to compute, since it can be calculated from promolecular densities. Thus, NCI analysis is applicable to large systems, including proteins and DNA, where analysis of non-covalent interactions is of great potential value. Here, we describe the NCI computational algorithms and their implementation for the analysis and visualization of weak interactions, using both self-consistent fully quantum-mechanical, as well as promolecular, densities. A wide range of options for tuning the range of interactions to be plotted is also presented. To demonstrate the capabilities of our approach, several examples are given from organic, inorganic, solid state, and macromolecular chemistry, including cases where NCI analysis gives insight into unconventional chemical bonding. The NCI code and its manual are available for download at http://www.chem.duke.edu/~yang/software.htm.     

Conformationally averaged score functions for electronic propagation in proteins

We explore the influence of conformational dynamics on protein-mediated electron donor-acceptor interactions. We introduce a thermally averaged score function to characterize electronic propagation from redox cofactors into the protein and solvent. The score function is explored for myoglobin at the extended-Hückel level, and the results are compared with those of simpler models. The conformationally averaged quantum results are consistent with the empirical analysis of the Pathways model. Notably, subtle effects of quantum interference among multiple coupling pathways that arise in static structures are largely averaged out when protein thermal motion is included. Propagation through bulk water near the single-protein interface decays rapidly with distance.     

Persistence of structure over fluctuations in biological electron-transfer reactions

In the soft-wet environment of biomolecular electron transfer, it is possible that structural fluctuations could wash out medium-specific electronic effects on electron tunneling rates. We show that beyond a transition distance (2-3 A in water and 6-7 A in proteins), fluctuation contributions to the mean-squared donor-to-acceptor tunneling matrix element are likely to dominate over the average matrix element. Even though fluctuations dominate the tunneling mechanism at larger distances, we find that the protein fold is "remembered" by the electronic coupling, and structure remains a key determinant of electron transfer kinetics.     

Hepatitis C virus NS5B polymerase: QM/MM calculations show the important role of the internal energy in ligand binding

The inter- and intramolecular interactions that determine the experimentally observed binding mode of the ligand (2Z)-2-(benzoylamino)-3-[4-(2-bromophenoxy)phenyl]-2-propenoate in complex with hepatitis C virus NS5B polymerase have been studied using QM/MM calculations. DFT-based QM/MM optimizations were performed on a number of ligand conformers in the protein-ligand complex. Using these initial poses, our aim is 2-fold. First, we identify the minimum energy pose. Second, we dissect the energetic contributions to this pose using QM/MM methods. The study reveals the critical importance of internal energy for the proper energy ranking of the docked poses. Using this protocol, we successfully identified three poses that have low RMSD with respect to the crystallographic structure from among the top 20 initially docked poses. We show that the most important energetic component contributing to binding for this particular protein-ligand system is the conformational (i.e., QM internal) energy.     

Computational design, synthesis and biological evaluation of para-quinone-based inhibitors for redox regulation of the dual-specificity phosphatase Cdc25B

Quinoid inhibitors of Cdc25B were designed based on the Linear Combination of Atomic Potentials (LCAP) methodology. In contrast to a published hypothesis, the biological activities and hydrogen peroxide generation in reducing media of three synthetic models did not correlate with the quinone half-wave potential, E(1/2).     

Solution structure of a peptide nucleic acid duplex from NMR data: features and limitations

This paper describes the results of a 1D and 2D NMR spectroscopy study of a palindromic 8-base pair PNA duplex GGCATGCC in H2O and H2O-D2O solutions. The (1)H NMR peaks have been assigned for most of the protons of the six central base pairs, as well as for several amide protons of the backbone. The resulting 36 interbase and base-backbone distance restraints were used together with Watson-Crick restraints to generate the PNA duplex structure in the course of 10 independent simulated annealing runs followed by restrained molecular dynamics (MD) simulations in explicit water. The resulting PNA structures correspond to a P-type helix with helical parameters close to those observed in the crystal structures of PNA. Based on the current limited number of restraints obtained from NMR spectra, alternative structures obtained by MD from starting PNA models based on DNA cannot be ruled out and are also discussed.     

Dynamic docking and electron-transfer between cytochrome b5 and a suite of myoglobin surface-charge mutants. Introduction of a functional-docking algorithm for protein-protein complexes

Horse myoglobin (Mb) provides a convenient "workbench" for probing the effects of electrostatics on binding and reactivity in the dynamic [Mb, cytochrome b(5)] electron-transfer (ET) complex. We have combined mutagenesis and heme neutralization to prepare a suite of six Mb surface-charge variants: the [S92D]Mb and [V67R]Mb mutants introduce additional charges on the "front" face, and incorporation of the heme di-ester into each of these neutralizes the charge on the heme propionates which further increases the positive charge on the "front" face. For this set of mutants, the nominal charge of Mb changes by -1 to +3 units relative to that for native Mb. For each member of this set, we have measured the bimolecular quenching rate constant (k(2)) for the photoinitiated (3)ZnDMb --> Fe(3+)b(5) ET reaction as a function of ionic strength. We find: (i) a dramatic decoupling of binding and reactivity, in which k(2) varies approximately 10(3)-fold within the suite of Mbs without a significant change in binding affinity; (ii) the ET reaction occurs within the "thermodynamic" or "rapid exchange" limit of the "Dynamic Docking" model, in which a large ensemble of weakly bound protein-protein configurations contribute to binding, but only a few are reactive, as shown by the fact that the zero-ionic-strength bimolecular rate constant varies exponentially with the net charge on Mb; (iii) Brownian dynamic docking profiles allow us to visualize the microscopic basis of dynamic docking. To describe these results we present a new theoretical approach which mathematically combines PATHWAY donor/acceptor coupling calculations with Poisson-Boltzmann-based electrostatics estimates of the docking energetics in a Monte Carlo (MC) sampling framework that is thus specially tailored to the intermolecular ET problem. This procedure is extremely efficient because it targets only the functionally active complex geometries by introducing a "reactivity filter" into the computations themselves, rather than as a subsequent step. This efficiency allows us to employ more computationally expensive and accurate methods to describe the relevant intermolecular interaction energies and the protein-mediated donor/acceptor coupling interactions. It is employed here to compute the changes in the bimolecular rate constant for ET between Mb and cyt b(5) upon variations in the myoglobin surface charge, pH, and ionic strength.     

Ab initio based calculations of electron-transfer rates in metalloproteins

A long-standing challenge in electron-transfer theory is to compute accurate rates of long-distance reactions in proteins. We describe an ab initio Hartree-Fock approach to compute electronic-coupling interactions and electron-transfer rates in proteins that allows the favorable comparison with experiment. The method includes the following key features; each is essential for reliable rate computations: (1) summing contributions over multiple tunneling pathways, (2) averaging couplings over thermally accessible protein conformations, (3) describing donor and acceptor electronic structure explicitly, including solvation effects, and averaging coupling over multiple energy-level crossings of the nearly degenerate donor-acceptor ligand-field states, and (4) eliminating basis set artifacts associated with diffuse basis functions. The strong dependence of coupling on donor-acceptor distance and on pathway interferences causes large variations of the computed electron-coupling values with protein geometry, and the strongest coupled conformers dominate the electron-transfer rate. As such, averaging over thermally accessible conformers of the protein and of the redox cofactors is essential. This approach was tested on six ruthenium-modified azurin derivatives using the high temperature nonadiabatic rate expression and compared with simpler pathways, average barrier, and semiempirical INDO models. Results of ab initio Hartree-Fock calculations with a split-valence basis set are in good agreement with the experimental rates. Predicted rates in the longer-distance derivatives are underestimated by 3-8-fold. This analysis indicates that the key ingredients needed for quantitatively reliable protein electron-transfer rate calculations are accessible.