Tomáš Kubař

Nonadiabatic QM/MM Simulations of Fast Charge Transfer in Escherichia coli DNA Photolyase

In this report, we study the photoactivation process in Escherichia coli DNA photolyase, involving long-range electron transport along a conserved chain of Trp residues between the protein surface and the flavin adenine dinucleotide (FAD) cofactor. Fully coupled nonadiabatic (Ehrenfest) quantum mechanics/molecular mechanics (QM/MM) simulations allow us to follow the time evolution of charge distributions over the natural time scale of multiple charge transfer events and conduct rigorous statistical analysis. Charge transfer rates in excellent agreement with experimental data are obtained without the need for any system-specific parametrization. The simulations are shown to provide a more detailed picture of electron transfer than a classical analysis of Marcus parameters. The protein and solvent both strongly influence the localization and transport properties of a positive charge, but the directionality of the process is mainly caused by solvent polarization. The time scales of charge movement, delocalization, protein relaxation and solvent reorganization overlap and lead to nonequilibrium reaction conditions. All these contributions are explicitly considered and fully resolved in the model used and provide an intricate picture of multistep biochemical electron transfer in a flexible, heterogeneous environment.

Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems

We present an extension to the recent 3OB parametrization of the Density Functional Tight Binding Model DFTB31,2 for biological and organic systems. Parameters for the halogens F, Cl, Br, and I have been developed for use in covalently bound systems and benchmarked on a test set of 106 molecules (the ‘OrgX’ set), using bonding distances, bonding angles, atomization energies, and vibrational frequencies to assess the performance of the parameters. Additional testing has been done with the X40 set of 40 supramolecular systems containing halogens,3 adding a simple correction for the halogen bonds that are strongly overbound in DFTB3. Furthermore, parameters for Ca, K, and Na as counterions in biological systems have been created. To benchmark geometries as well as ligand binding energies a test set ‘BioMe’ of 210 molecules has been created that cover coordination to various functional groups frequently occurring in biological systems. The new DFTB3/3OB parameter set outperforms DFT calculations with a double-ζ basis set in terms of energies and can reproduce DFT geometries, with some minor deviations in bond distances and angles due to the use of a minimal basis set.

Solvent Fluctuations Drive the Hole Transfer in DNA: A Mixed Quantum−Classical Study

We studied the hole transfer across adenine bridges in double-stranded DNA by means of a multiscale approach, propagating the hole in the framework of time-dependent DFT coupled to classical molecular dynamics simulation using a QM/MM scheme. The hole transfer in DNA is codetermined by the large fluctuations of site energies on the order of 0.4 eV, induced by the solvent degrees of freedom. These fluctuations lead to charge-transfer active conformations with large transfer efficiency, which are characterized by a favorable alignment of site energies along the DNA strand. This reduces the barrier for the hole transfer dramatically. Consequently, we find that a charge hopping mechanism is operative already for short bridges with fewer than four adenines, in contrast to the charge-transfer models assuming static DNA structures, where only tunneling occurs. The solvent fluctuations introduce a significant correlation between neighboring sites, enhancing the charge-transfer rate, while the fluctuation of electronic couplings has only a minor impact on the charge-transfer characteristics. Our results emphasize the importance of an accurate description of solvent effects as well as proper sampling, and it is suggested that charge transfer in DNA is gated by the dynamics of solvent.

Semiempirical Quantum Mechanical Methods for Noncovalent Interactions for Chemical and Biochemical Applications

Semiempirical (SE) methods can be derived from either Hartree–Fock or density functional theory by applying systematic approximations, leading to efficient computational schemes that are several orders of magnitude faster than ab initio calculations. Such numerical efficiency, in combination with modern computational facilities and linear scaling algorithms, allows application of SE methods to very large molecular systems with extensive conformational sampling. To reliably model the structure, dynamics, and reactivity of biological and other soft matter systems, however, good accuracy for the description of noncovalent interactions is required. In this review, we analyze popular SE approaches in terms of their ability to model noncovalent interactions, especially in the context of describing biomolecules, water solution, and organic materials. We discuss the most significant errors and proposed correction schemes, and we review their performance using standard test sets of molecular systems for quantum chemical methods and several recent applications. The general goal is to highlight both the value and limitations of SE methods and stimulate further developments that allow them to effectively complement ab initio methods in the analysis of complex molecular systems.

Combined density functional theory and Landauer approach for hole transfer in DNA along classical molecular dynamics trajectories

We investigate in detail the charge transport characteristics of DNA wires with various sequences and lengths in the presence of solvent. Our approach combines large-scale quantum/classical molecular dynamics (MD) simulations with transport calculations based on Landauer theory. The quantum mechanical transmission function of the wire is calculated along MD trajectories and thus encodes the influence of dynamical disorder arising from the environment (water, backbone, counterions) and from the internal base dynamics. We show that the correlated fluctuations of the base pair dynamics are crucial in determining the transport properties of the wire and that the effect of fluctuations can be quite different for sequences with low and high static disorders (differences in base ionization potentials). As a result, in structures with high static disorder as is the case of the studied Dickerson dodecamer, the weight of high-transmissive structures increases due to dynamical fluctuations and so does the calculated average transmission. Our analysis further supports the basic intuition of charge-transfer active conformations as proposed by Barton et al. [J. Am. Chem. Soc. 126, 11471 (2004)]. However, not DNA conformations with good stacking contacts leading to large interbase hopping values are necessarily the most important, but rather those where the average fluctuation of ionization potentials along the base stack is small. The reason behind this is that the ensemble of conformations leads to average electronic couplings, which are large enough for sufficient transmission. On the other hand, the alignment of onsite energies is the critical parameter which gates the charge transport.

A hybrid approach to simulation of electron transfer in complex molecular systems

Electron transfer (ET) reactions in biomolecular systems represent an important class of processes at the interface of physics, chemistry and biology. The theoretical description of these reactions constitutes a huge challenge because extensive systems require a quantum-mechanical treatment and a broad range of time scales are involved. Thus, only small model systems may be investigated with the modern density functional theory techniques combined with non-adiabatic dynamics algorithms. On the other hand, model calculations based on Marcuss seminal theory describe the ET involving several assumptions that may not always be met. We review a multi-scale method that combines a non-adiabatic propagation scheme and a linear scaling quantum-chemical method with a molecular mechanics force field in such a way that an unbiased description of the dynamics of excess electron is achieved and the number of degrees of freedom is reduced effectively at the same time. ET reactions taking nanoseconds in systems with hundreds of quantum atoms can be simulated, bridging the gap between non-adiabatic ab initio simulations and model approaches such as the Marcus theory. A major recent application is hole transfer in DNA, which represents an archetypal ET reaction in a polarizable medium. Ongoing work focuses on hole transfer in proteins, peptides and organic semi-conductors.

Structural stability versus conformational sampling in biomolecular systems: Why is the charge transfer efficiency in G4-DNA better than in double-stranded DNA?

The electrical conduction properties of G4-DNA are investigated using a hybrid approach, which combines electronic structure calculations, molecular dynamics (MD) simulations, and the formulation of an effective tight-binding model Hamiltonian. Charge transport is studied by computing transmission functions along the MD trajectories. Though G4-DNA is structurally more stable than double-stranded DNA (dsDNA), our results strongly suggest that the potential improvement of the electrical transport properties in the former is not necessarily related to an increased stability, but rather to the fact that G4 is able to explore in its conformational space a larger number of charge-transfer active conformations. This in turn is a result of the non-negligible interstrand matrix elements, which allow for additional charge transport pathways. The higher structural stability of G4 can however play an important role once the molecules are contacted by electrodes. In this case, G4 may experience weaker structural distortions than dsDNA and thus preserve to a higher degree its conduction properties.

Parametrization of the SCC-DFTB Method for Halogens

Parametrization of the approximative DFT method SCC-DFTB for halogen elements is presented. The new parameter set is intended to describe halogenated organic as well as inorganic molecules, and it is compatible with the established parametrization of SCC-DFTB for carbon, hydrogen, oxygen, and nitrogen. The performance of the parameter set is tested on a representative set of molecules and discussed.

Solvent driving force ensures fast formation of a persistent and well-separated radical pair in plant cryptochrome.

The photoreceptor protein cryptochrome is thought to host, upon light absorption, a radical pair that is sensitive to very weak magnetic fields, endowing migratory birds with a magnetic compass sense. The molecular mechanism that leads to formation of a stabilized, magnetic field sensitive radical pair has despite various theoretical and experimental efforts not been unambiguously identified yet. We challenge this unambiguity through a unique quantum mechanical molecular dynamics approach where we perform electron transfer dynamics simulations taking into account the motion of the protein upon the electron transfer. This approach allows us to follow the time evolution of the electron transfer in an unbiased fashion and to reveal the molecular driving force that ensures fast electron transfer in cryptochrome guaranteeing formation of a persistent radical pair suitable for magnetoreception. We argue that this unraveled molecular mechanism is a general principle inherent to all proteins of the cryptochrome/photolyase family and that cryptochromes are, therefore, tailored to potentially function as efficient chemical magnetoreceptors.

Charge transfer in model peptides: obtaining Marcus parameters from molecular simulation

Charge transfer within and between biomolecules remains a highly active field of biophysics. Due to the complexities of real systems, model compounds are a useful alternative to study the mechanistic fundamentals of charge transfer. In recent years, such model experiments have been underpinned by molecular simulation methods as well. In this work, we study electron hole transfer in helical model peptides by means of molecular dynamics simulations. A theoretical framework to extract Marcus parameters of charge transfer from simulations is presented. We find that the peptides form stable helical structures with sequence dependent small deviations from ideal PPII helices. We identify direct exposure of charged side chains to solvent as a cause of high reorganization energies, significantly larger than typical for electron transfer in proteins. This, together with small direct couplings, makes long-range superexchange electron transport in this system very slow. In good agreement with experiment, direct transfer between the terminal amino acid side chains can be dicounted in favor of a two-step hopping process if appropriate bridging groups exist.