O to bR transition in bacteriorhodopsin occurs through a proton hole mechanism.

Extensive classical and quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations are used to establish the structural features of the O state in bacteriorhodopsin (bR) and its conversion back to the bR ground state. The computed free energy surface is consistent with available experimental data for the kinetics and thermodynamics of the O to bR transition. The simulation results highlight the importance of the proton release group (PRG, consisting of Glu194/204) and the conserved arginine 82 in modulating the hydration level of the protein cavity. In particular, in the O state, deprotonation of the PRG and downward rotation of Arg82 lead to elevated hydration level and a continuous water network that connects the PRG to the protonated Asp85. Proton exchange through this water network is shown by ∼0.1-µs semiempirical QM/MM free energy simulations to occur through the generation and propagation of a proton hole, which is relayed by Asp212 and stabilized by Arg82. This mechanism provides an explanation for the observation that the D85S mutant of bacteriorhodopsin pumps chloride ions. The electrostatics-hydration coupling mechanism and the involvement of all titration states of water are likely applicable to many biomolecules involved in bioenergetic transduction.

Tracking the Electron Transfer Cascade in European Robin Cryptochrome 4 Mutants.

The primary step in the mechanism by which migratory birds sense the Earth's magnetic field is thought to be the light-induced formation of long-lived magnetically sensitive radical pairs within cryptochrome flavoproteins located in the birds' retinas. Blue-light absorption by the non-covalently bound flavin chromophore triggers sequential electron transfers along a chain of four tryptophan residues toward the photoexcited flavin. The recently demonstrated ability to express cryptochrome 4a from the night-migratory European robin (Erithacus rubecula), ErCry4a, and to replace each of the tryptophan residues by a redox-inactive phenylalanine offers the prospect of exploring the roles of the four tryptophans. Here, we use ultrafast transient absorption spectroscopy to compare wild type ErCry4a and four mutants having a phenylalanine at different positions in the chain. We find that each of the three tryptophan residues closest to the flavin adds a distinct relaxation component (time constants: 0.5, 30, and 150 ps) in the transient absorption data. The dynamics of the mutant containing a phenylalanine at the fourth position, furthest from the flavin, are very similar to those of wild type ErCry4a, except for a reduced concentration of long-lived radical pairs. The experimental results are evaluated and discussed in the framework of real-time quantum mechanical/molecular mechanical electron transfer simulations based on the density functional-based tight binding approach. This comparison between simulation results and experimental measurements provides a detailed microscopic insight into the sequential electron transfers along the tryptophan chain. Our results offer a route to the study of spin transport and dynamical spin correlations in flavoprotein radical pairs.

Extended-sampling QM/MM simulation of biochemical reactions involving P-N bonds.

The description of the phosphate group and its reactions with nitrogen species appears to be challenging using semi-empirical quantum chemical methods, and this holds for DFTB3 too. A new parameterization of DFTB3, consisting of a new P-N repulsive function, has been developed to improve its performance for reactions in which a P-N bond is replaced by a P-O bond or vice versa. Extended-sampling QM/MM simulations using the new parameterization of DFTB3 represent biochemical phosphorylation and hydrolysis reactions involving P-N bonds accurately. The parameter set is benchmarked on a reaction modeling the autophosphorylation of histidine, and is applied to study the complex mechanism of the acidic hydrolysis of an anticancer drug, as well as to the autophosphorylation of a genuine histidine kinase protein.

Mechanism of proton-coupled electron transfer described with QM/MM implementation of coupled-perturbed density-functional tight-binding.

Coupled-perturbed equations for degenerate orbitals were implemented for third order density-functional tight binding, which allowed the use of Mulliken charges as reaction coordinates. The method was applied to proton-coupled electron transfer (PCET) reactions in a model system and thoroughly tested for QM and QM/MM setups (i.e., coupled quantum and molecular mechanics). The performed enhanced sampling simulations were stable, and the obtained potentials of the mean force were able to address the thermodynamic and kinetic features of the reactions by showing the expected topography and energy barriers. Hence, this method has the potential to distinguish between concerted and sequential mechanisms and could next be applied to proton-coupled electron transfer reactions in more complex systems like proteins.

Unravelling the mechanism of glucose binding in a protein-based fluorescence probe: molecular dynamics simulation with a tailor-made charge model.

Fluorophores linked to the glucose/galactose-binding protein (GGBP) are a promising class of glucose sensors with potential application in medical devices for diabetes patients. Several different fluorophores at different positions in the protein were tested experimentally so far, but a deeper molecular understanding of their function is still missing. In this work, we use molecular dynamics simulations to investigate the mechanism of glucose binding in the GGBP-Badan triple mutant and make a comparison to the GGBP wild-type protein. The aim is to achieve a detailed molecular understanding of changes in the glucose binding site due to the mutations and their effect on glucose binding. Free simulations give an insight into the changes of the hydrogen-bonding network in the active site and into the mechanisms of glucose binding. Additionally, metadynamics simulations for wild type and mutant unravel the energetics of binding/unbinding in these proteins. Computed free energies for the opening of the binding pocket for the wild-type and the mutant agree well with the experimental data. Further, the simulations also give an insight into the changes of the chromophore conformations upon glucose binding, which can help to understand fluorescence changes. Therefore, the molecular details unravelled in this work may support effective optimisation strategies for the construction of more efficient glucose sensors.

Reduction pathway of glutaredoxin 1 investigated with QM/MM molecular dynamics using a neural network correction.

Glutaredoxins are small enzymes that catalyze the oxidation and reduction of protein disulfide bonds by the thiol-disulfide exchange mechanism. They have either one or two cysteines in their active site, resulting in different catalytic reaction cycles that have been investigated in many experimental studies. However, the exact mechanisms are not yet fully known, and to our knowledge, no theoretical studies have been performed to elucidate the underlying mechanism. In this study, we investigated a proposed mechanism for the reduction of the disulfide bond in the protein HMA4n by a mutated monothiol Homo sapiens glutaredoxin and the co-substrate glutathione. The catalytic cycle involves three successive thiol-disulfide exchanges that occur between the molecules. To estimate the regioselectivity of the different attacks, classical molecular dynamics simulations were performed and the trajectories analyzed regarding the sulfur-sulfur distances and the attack angles between the sulfurs. The free energy profile of each reaction was obtained with hybrid quantum mechanical/molecular mechanical metadynamics simulations. Since this required extensive phase space sampling, the semi-empirical density functional tight-binding method was used to describe the reactive cysteines. For an accurate description, we used specific reaction parameters fitted to B3LYP energies of the thiol-disulfide exchange and a machine learned energy correction that was trained on coupled-cluster single double perturbative triple [CCSD(T)] energies of thiol-disulfide exchanges. Our calculations show the same regiospecificity as observed in the experiment, and the obtained barrier heights are about 12 and 20 kcal/mol for the different reaction steps, which confirms the proposed pathway.

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.

What accounts for the different functions in photolyases and cryptochromes: a computational study of proton transfers to FAD.

Photolyases (PL) and cryptochromes (CRY) are light-sensitive flavoproteins, respectively, involved in DNA repair and signal transduction. Their activation is triggered by an electron transfer process, which partially or fully reduces the photo-activated FAD cofactor. The full reduction additionally requires a proton transfer to the isoalloxazine ring. In plant CRY, an efficient proton transfer takes place within several µs, enabled by a conserved aspartate working as a proton donor, whereas in E. coli PL a proton transfer occurs in the 4 s timescale without any obvious proton donor, indicating the presence of a long-range proton transfer pathway. Unexpectedly, the insertion of an aspartate as a proton donor in a suitable position for proton transfer in E. coli PL does not initiate a transfer process similar to plant CRY, but even prevents the formation of a protonated FAD. In the present work, thanks to a combination of classical molecular dynamics and state-of-the-art DFTB3/MM simulations, we identify a proton transfer pathway from bulk to FAD in E. coli PL associated with a free energy profile in agreement with the experimental kinetics data. The free energy profiles of the proton transfer between aspartate and FAD show an inversion of the driving force between plant CRY and E. coli PL mutants. In the latter, the proton transfer from the aspartate is faster than in plant CRY but also thermodynamically disfavoured, in agreement with the experimental data. Our results further illustrate the fine tuning of the electrostatic FAD environment and the adaptability of the FAD pocket to ensure the divergent functions of the members of the PL-CRY family.

Molecular Insights into Variable Electron Transfer in Amphibian Cryptochrome.

Cryptochrome proteins are activated by the absorption of blue light, leading to the formation of radical pairs through electron transfer in the active site. Recent experimental studies have shown that once some of the amino acid residues in the active site of Xenopus laevis cryptochrome DASH are mutated, radical-pair formation is still observed. In this study, we computationally investigate electron-transfer pathways in the X. laevis cryptochrome DASH by extensively equilibrating a previously established homology model using molecular dynamics simulations and then mutating key amino acids involved in the electron transfer. The electron-transfer pathways are then probed by using tight-binding density-functional theory. We report the alternative electron-transfer pathways resolved at the molecular level and, through comparison of amino acid sequences for cryptochromes from different species, we demonstrate that one of these alternative electron-transfer pathways could be general for all cryptochrome DASH proteins.

Orthogonal 19 F-Labeling for Solid-State NMR Spectroscopy Reveals the Conformation and Orientation of Short Peptaibols in Membranes.

Peptaibols are promising drug candidates in view of their interference with cellular membranes. Knowledge of their lipid interactions and membrane-bound structure is needed to understand their activity and should be, in principle, accessible by solid-state NMR spectroscopy. However, their unusual amino acid composition and noncanonical conformations make it very challenging to find suitable labels for NMR spectroscopy. Particularly in the case of short sequences, new strategies are required to maximize the structural information that can be obtained from each label. Herein, l-3-(trifluoromethyl)bicyclopent[1.1.1]-1-ylglycine, (R)- and (S)-trifluoromethylalanine, and 15 N-backbone labels, each probing a different direction in the molecule, have been combined to elucidate the conformation and membrane alignment of harzianin HK-VI. For the short sequence of 11 amino acids, 12 orientational constraints have been obtained by using 19 F and 15 N NMR spectroscopy. This strategy revealed a ß-bend ribbon structure, which becomes realigned in the membrane from a surface-parallel state towards a membrane-spanning state, with increasing positive spontaneous curvature of the lipids.

"siRNA traffic lights": arabino-configured 2'-anchors for fluorescent dyes are key for dual color readout in cell imaging.

Two fluorescent dyes covalently attached in diagonal interstrand orientation to siRNA undergo energy transfer and thereby enable a dual color fluorescence readout (red/green) for hybridization. Three different structural variations were carried out and compared by their optical properties, including (i) the base surrogate approach with an acyclic linker as a substitute of the 2-deoxyriboside between the phosphodiester bridges, (ii) the 2'-modification of conventional ribofuranosides and (iii) the arabino-configured 2'-modification. The double stranded siRNA with the latter type of modification delivered the best energy transfer efficiency, which was explained by molecular dynamics simulations that showed that the two dyes are more flexible at the arabino-configured sugars compared to the completely stacked situation at the ribo-configured ones. Single molecule fluorescence lifetime measurements indicate their application in fluorescence cell imaging, which reveals a red/green fluorescence contrast in particular for the arabino-configured 2'-modification by the two dyes, which is key for tracking of siRNA transport into HeLa cells.

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.

On the mechanism of spontaneous thiol-disulfide exchange in proteins.

The thiol-disulfide exchange reaction in model systems and small peptides was investigated by means of a combined QM/MM metadynamics scheme. The free energy landscapes of these systems were generated, providing the structures of reactants and products with atomic detail, as well as the heights of free energy barriers (or, activation energies) opposing the spontaneous exchange. A QM/MM scheme with purely classical water turned out to be an efficient and accurate compromise solution. The calculations yielded the expected symmetric trisulfide transition state at S-S distances of 2.7 Å, interestingly, with a slight deviation from linearity at an S-S-S angle of 165°. The structure of the transition state as well as the free energy barrier were very similar for the intramolecular thiol-disulfide reactions in model peptides. While CXC disulfide bonds were found sterically unfavorable, CXXC were favored over longer-range disulfide bonds along the peptide backbone, in line with the high abundance of CXXC motifs in redox proteins.

Coupled-perturbed DFTB-QM/MM metadynamics: Application to proton-coupled electron transfer.

We present a new concept of free energy calculations of chemical reactions by means of extended sampling molecular dynamics simulations. Biasing potentials are applied on partial atomic charges, which may be combined with atomic coordinates either in a single collective variable or in multi-dimensional biasing simulations. The necessary additional gradients are obtained by solving coupled-perturbed equations within the approximative density-functional tight-binding method. The new computational scheme was implemented in a combination of Gromacs and Plumed. As a prospective application, proton-coupled electron transfer in a model molecular system is studied. Two collective variables are introduced naturally, one for the proton transfer and the other for the electron transfer. The results are in qualitative agreement with the extended free simulations performed for reference. Free energy minima as well as the mechanism of the process are identified correctly, while the topology of the transition region and the height of the energy barrier are only reproduced qualitatively. The application also illustrates possible difficulties with the new methodology. These may be inefficient sampling of spatial coordinates when atomic charges are biased exclusively and a decreased stability of the simulations. Still, the new approach represents a viable alternative for free energy calculations of a certain class of chemical reactions, for instance a proton-coupled electron transfer in proteins.

Structural Behavior of the Peptaibol Harzianin HK VI in a DMPC Bilayer: Insights from MD Simulations.

Microsecond molecular dynamics simulations of harzianin HK VI (HZ) interacting with a dimyristoylphosphatidylcholine bilayer were performed at the condition of low peptide-to-lipid ratio. Two orientations of HZ molecule in the bilayer were found and characterized. In the orientation perpendicular to the bilayer surface, HZ induces a local thinning of the bilayer. When inserted into the bilayer parallel to its surface, HZ is located nearly completely within the hydrophobic region of the bilayer. A combination of solid-state NMR and circular dichroism experiments found the latter orientation to be dominant. An extended sampling simulation provided qualitative results and showed the same orientation to be a global minimum of free energy. The secondary structure of HZ was characterized, and it was found to be located in the 310-helical family. The specific challenges of computer simulation of nonpolar peptides are discussed briefly.

Cation solvation with quantum chemical effects modeled by a size-consistent multi-partitioning quantum mechanics/molecular mechanics method.

In the condensed phase, quantum chemical properties such as many-body effects and intermolecular charge fluctuations are critical determinants of the solvation structure and dynamics. Thus, a quantum mechanical (QM) molecular description is required for both solute and solvent to incorporate these properties. However, it is challenging to conduct molecular dynamics (MD) simulations for condensed systems of sufficient scale when adapting QM potentials. To overcome this problem, we recently developed the size-consistent multi-partitioning (SCMP) quantum mechanics/molecular mechanics (QM/MM) method and realized stable and accurate MD simulations, using the QM potential to a benchmark system. In the present study, as the first application of the SCMP method, we have investigated the structures and dynamics of Na+, K+, and Ca2+ solutions based on nanosecond-scale sampling, a sampling 100-times longer than that of conventional QM-based samplings. Furthermore, we have evaluated two dynamic properties, the diffusion coefficient and difference spectra, with high statistical certainty. Furthermore the calculation of these properties has not previously been possible within the conventional QM/MM framework. Based on our analysis, we have quantitatively evaluated the quantum chemical solvation effects, which show distinct differences between the cations.

The effect of the environment on the methyl transfer reaction mechanism between trimethylsulfonium and phenolate.

Methyl transfer reactions play an important role in biology and are catalyzed by various enzymes. Here, the influence of the molecular environment on the reaction mechanism was studied using advanced ab initio methods, implicit solvation models and QM/MM molecular dynamics simulations. Various conceptual DFT and electronic structure descriptors identified different processes along the reaction coordinate e.g. electron transfer. The results show that the polarity of the solvent increases the energy required for the electron transfer and that this spontaneous process is located in the transition state region identified by the (mean) reaction force analysis and takes place through the bonds which are broken and formed. The inclusion of entropic contributions and hydrogen bond interactions in QM/MM molecular dynamics simulations with a validated DFTB3 Hamiltonian yields activation barriers in good agreement with the experimental values in contrast to the values obtained using two implicit solvation models.

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.

New QM/MM implementation of the DFTB3 method in the gromacs package.

The approximate density-functional tight-binding theory method DFTB3 has been implemented in the quantum mechanics/molecular mechanics (QM/MM) framework of the Gromacs molecular simulation package. We show that the efficient smooth particle-mesh Ewald implementation of Gromacs extends to the calculation of QM/MM electrostatic interactions. Further, we make use of the various free-energy functionalities provided by Gromacs and the PLUMED plugin. We exploit the versatility and performance of the current framework in three typical applications of QM/MM methods to solve biophysical problems: (i) ultrafast proton transfer in malonaldehyde, (ii) conformation of the alanine dipeptide, and (iii) electron-induced repair of a DNA lesion. Also discussed is the further development of the framework, regarding mostly the options for parallelization.

Force distribution analysis of mechanochemically reactive dimethylcyclobutene.

Internal molecular forces can guide chemical reactions, yet are not straightforwardly accessible within a quantum mechanical description of the reacting molecules. Here, we present a force-matching force distribution analysis (FM-FDA) to analyze internal forces in molecules. We simulated the ring opening of trans-3,4-dimethylcyclobutene (tDCB) with on-the-fly semiempirical molecular dynamics. The self-consistent density functional tight binding (SCC-DFTB) method accurately described the force-dependent ring-opening kinetics of tDCB, showing quantitative agreement with both experimental and computational data at higher levels. Mechanical force was applied in two different ways, namely, externally by a constant pulling force and internally by embedding tDCB within a strained macrocycle-containing stiff stilbene. We analyzed the distribution of tDCB internal forces in the two different cases by FM-FDA and found that external force gave rise to a symmetric force distribution in the cyclobutene ring, which also scaled linearly with the external force, indicating that the force distribution was uniquely determined by the symmetric architecture of tDCB. In contrast, internal forces due to stiff stilbene resulted in an asymmetric force distribution within tDCB, which indicated a different geometry of force application and supported the important role of linkers in the mechanochemical reactivity of tDCB. In addition, three coordinates were identified through which the distributed forces contributed most to rate acceleration. These coordinates are mostly parallel to the coordinate connecting the two CH3 termini of tDCB. Our results confirm previous observations that the linker outside of the reactive moiety, such as a stretched polymer or a macrocycle, affects its mechanochemical reactivity. We expect FM-FDA to be of wide use to understand and quantitatively predict mechanochemical reactivity, including the challenging cases of systems within strained macrocycles.