Hybrid Quantum Mechanical/ Molecular Mechanical Methods for Studying Energy Transduction in Biomolecular Machines

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems

Electrostatic interactions contribute to the control of intramolecular thiol-disulfide isomerization in a protein

The roles of structural factors and of electrostatic interactions with the environment on the outcome of thiol–disulfide exchange reactions were investigated in a mutated immunoglobulin domain (I27*) under mechanical stress. An extensive ensemble of molecular dynamics trajectories was generated by means of QM/MM simulations for a total sampling of 5.7 μs. A significant number of thiol–disulfide exchanges were observed, and the Cys32 thiolate preferred to attack Cys55 over Cys24, in agreement with previous experimental and computational studies. The structural features as well as electronic structures of the thiol–disulfide system along the reaction were analyzed, as were the electrostatic interactions with the environment. The previous findings of better accessibility of Cys55 were confirmed. Additionally, the reaction was found to be directed by the electrostatic interactions of the involved sulfur atoms with the molecular environment. The relationships of atomic charges, which stem from the electrostatic interactions, lead to the kinetic preference of the attack on Cys55. Further, QM/MM metadynamics simulations of thiol–disulfide exchange in a small model system with varied artificial external electric potentials revealed changes in reaction kinetics of the same magnitude as in I27*. Therefore, the electrostatic interactions are confirmed to play a role in the regioselectivity of the thiol–disulfide exchange reactions in the protein

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green's functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives

Fragment orbital based description of charge transfer in peptides including backbone orbitals

Charge transfer in peptides and proteins can occur on different pathways, depending on the energetic landscape as well as the coupling between the involved orbitals. Since details of the mechanism and pathways are difficult to access experimentally, different modeling strategies have been successfully applied to study these processes in the past. These can be based on a simple empirical pathway model, efficient tight binding type atomic orbital Hamiltonians or ab initio and density functional calculations. An interesting strategy, which allows an efficient calculations of charge transfer parameters, is based on a fragmentation of the system into functional units. While this works well for systems like DNA, where the charge transfer pathway is naturally divided into distinct molecular fragments, this is less obvious for charge transfer along peptide and protein backbones. In this work, we develop and access a strategy for an effective fragmentation approach, which allows one to compute electronic couplings for large systems along nanosecond time scale molecular dynamics trajectories. The new methodology is applied to a solvated peptide, for which charge transfer properties have been studied recently using an empirical pathway model. As could be expected, dynamical effects turn out to be important, which emphasizes the importance of using effective quantum approaches which allow for sufficient sampling. However, the computed rates are orders of magnitude smaller than experimentally determined, which indicates the shortcomings of present modeling approaches

Dispersion interactions govern the strong thermal stability of a protein

Rubredoxin from the hyperthermophile Pyrococcus furiosus (Pf Rd) is an extremely thermostable protein, which makes it an attractive subject of protein folding and stability studies. A fundamental question arises as to what the reason for such extreme stability is and how it can be elucidated from a complex set of interatomic interactions, We addressed this issue first theoretically through a computational analysis of the hydrophobic core of the protein and its mutants, including the interactions taking place inside the core, Here we show that a single mutation of one of phenylalanine\u27s residues inside the protein\u27s hydrophobic core results in a dramatic decrease in its thermal stability. The calculated unfolding Gibbs energy as well as the stabilization energy differences between a few core residues follows the same trend as the melting temperature of protein variants determined experimentally by microcalorimetry measurements. NMR spectroscopy experiments have shown that the only part of the protein affected by mutation is the reasonably rearranged hydrophobic core. It is hence concluded that stabilization energies, which are dominated by London dispersion, represent the main source of stability of this protein. © 2007 Wiley-VCH Verlag GmbH & Co. KGaA,