Reaction Path Averaging: Characterizing the Structural Response of the DNA Double Helix to Electron Transfer

A polarizable environment, prominently the solvent, responds to electronic changes in biomolecules rapidly. The knowledge of conformational relaxation of the biomolecule itself, however, may be scarce or missing. In this work, we describe in detail the structural changes in DNA undergoing electron transfer between two adjacent nucleobases. We employ an approach based on averaging of tens to hundreds of thousands of nonequilibrium trajectories generated with molecular dynamics simulation, and a reduction of dimensionality suitable for DNA. We show that the conformational response of the DNA proceeds along a single collective coordinate that represents the relative orientation of two consecutive base pairs, namely, a combination of helical parameters shift and tilt. The structure of DNA relaxes on time scales reaching nanoseconds, contributing marginally to the relaxation of energies, which is dominated by the modes of motion of the aqueous solvent. The concept of reaction path averaging (RPA), conveniently exploited in this context, makes it possible to filter out any undesirable noise from the nonequilibrium data, and is applicable to any chemical process in general.Comment: 45 pages, 20 figures, published, added Supplementary informatio

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

Disulfide bond reduction and exchange in C4 domain of von Willebrand factor undermines platelet binding

Background The von Willebrand factor (VWF) is a key player in regulating hemostasis through adhesion of platelets to sites of vascular injury. It is a large, multi-domain, mechano-sensitive protein that is stabilized by a net of disulfide bridges. Binding to platelet integrin is achieved by the VWF-C4 domain, which exhibits a fixed fold, even under conditions of severe mechanical stress, but only if critical internal disulfide bonds are closed. Objective To determine the oxidation state of disulfide bridges in the C4 domain of VWF and implications for VWF’s platelet binding function. Methods We combined classical molecular dynamics and quantum mechanical simulations, mass spectrometry, site-directed mutagenesis, and platelet binding assays. Results We show that 2 disulfide bonds in the VWF-C4 domain, namely the 2 major force-bearing ones, are partially reduced in human blood. Reduction leads to pronounced conformational changes within C4 that considerably affect the accessibility of the integrin-binding motif, and thereby impair integrin-mediated platelet binding. We also reveal that reduced species in the C4 domain undergo specific thiol/disulfide exchanges with the remaining disulfide bridges, in a process in which mechanical force may increase the proximity of specific reactant cysteines, further trapping C4 in a state of low integrin-binding propensity. We identify a multitude of redox states in all 6 VWF-C domains, suggesting disulfide bond reduction and swapping to be a general theme. Conclusions Our data suggests a mechanism in which disulfide bonds dynamically swap cysteine partners and control the interaction of VWF with integrin and potentially other partners, thereby critically influencing its hemostatic function

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