Ligand Photo-Isomerization Triggers Conformational Changes in iGluR2 Ligand Binding Domain

Neurological glutamate receptors bind a variety of artificial ligands, both agonistic and antagonistic, in addition to glutamate. Studying their small molecule binding properties increases our understanding of the central nervous system and a variety of associated pathologies. The large, oligomeric multidomain membrane protein contains a large and flexible ligand binding domains which undergoes large conformational changes upon binding different ligands. A recent application of glutamate receptors is their activation or inhibition via photo-switchable ligands, making them key systems in the emerging field of optochemical genetics. In this work, we present a theoretical study on the binding mode and complex stability of a novel photo-switchable ligand, ATA-3, which reversibly binds to glutamate receptors ligand binding domains (LBDs). We propose two possible binding modes for this ligand based on flexible ligand docking calculations and show one of them to be analogues to the binding mode of a similar ligand, 2-BnTetAMPA. In long MD simulations, it was observed that transitions between both binding poses involve breaking and reforming the T686-E402 protein hydrogen bond. Simulating the ligand photo-isomerization process shows that the two possible configurations of the ligand azo-group have markedly different complex stabilities and equilibrium binding modes. A strong but slow protein response is observed after ligand configuration changes. This provides a microscopic foundation for the observed difference in ligand activity upon light-switching

Biological relevance of charge transfer branching pathways in photolyases

In PhrA, a class III CPD photolyase, two branching tryptophan charge transfer pathways have been characterized in the mechanism of FAD photoreduction. To provide a molecular explanation of the charge transfer abilities of both pathways, we performed simulations where the protein motion and the positive charge are simultaneously propagated. Our computational approach reveals that one pathway drives a very fast charge transfer whereas the other pathway provides a very good thermodynamic stabilization of the positive charge. During the simulations, the positive charge firstly moves on the fast triad, allowing the stabilization of reduced FAD. After one nanosecond, we observe a nearly equal probability to find the charge at ending tryptophan of either pathway. Our results highlight the role of the protein environment, which manages the association of a kinetic and a thermodynamic pathways to trigger a fast and efficient FAD photoreduction

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

Molecular Dynamics Investigation of gluazo, a Photo-Switchable Ligand for the Glutamate Receptor GluK2

Photochromic ligands (PCLs), defined as photoswitchable molecules that are able to endow native receptors with a sensitivity towards light, have become a promising photo-pharmacological tool for various applications in biology. In general, PCLs consist of a ligand of the target receptor covalently linked to an azobenzene, which can be reversibly switched between two configurations upon light illumination. Gluazo, as a PCL that targets excitatory amino acid receptors, in its dark-adapted trans iso-form was characterized to be a partial agonist of the kainate glutamate receptor GluK2. Application of UV light leads to the formation of the cis form, with remarkedly reduced affinity towards GluK2. The mechanism of the change of ligand affinity induced by the photoisomerization was unresolved. The presented computational study explains how the somerization of such a PCL affects the structural changes in the target receptor that lead to its activation

HAB79: A new molecular dataset for benchmarking DFT and DFTB electronic couplings against high-level ab initio calculations

A new molecular dataset called HAB79 is introduced to provide ab initio reference values for electronic couplings (transfer integrals) and to benchmark density functional theory (DFT) and density functional tight-binding (DFTB) calculations. The HAB79 dataset is composed of 79 planar heterocyclic polyaromatic hydrocarbon molecules frequently encountered in organic (opto)electronics, arranged to 921 structurally diverse dimer configurations. We show that CASSCF/NEVPT2 with a minimal active space provides a robust reference method that can be applied to the relatively large molecules of the dataset. Electronic couplings are largest for cofacial dimers, in particular, sulfur-containing polyaromatic hydrocarbons, with values in excess of 0.5 eV, followed by parallel displaced cofacial dimers. V-shaped dimer motifs, often encountered in the herringbone layers of organic crystals, exhibit medium-sized couplings, whereas T-shaped dimers have the lowest couplings. DFT values obtained from the projector operator-based diabatization (POD) method are initially benchmarked against the smaller databases HAB11 (HAB7-) and found to systematically improve when climbing Jacob’s ladder, giving mean relative unsigned errors (MRUEs) of 27.7% (26.3%) for the generalized gradient approximation (GGA) functional BLYP, 20.7% (15.8%) for hybrid functional B3LYP, and 5.2% (7.5%) for the long-range corrected hybrid functional omega-B97X. Cost-effective POD in combination with a GGA functional and very efficient DFTB calculations on the dimers of the HAB79 database give a good linear correlation with the CASSCF/NEVPT2 reference data, which, after scaling with a multiplicative constant, gives reasonably small MRUEs of 17.9% and 40.1%, respectively, bearing in mind that couplings in HAB79 vary over 4 orders of magnitude. The ab initio reference data reported here are expected to be useful for benchmarking other DFT or semi-empirical approaches for electronic coupling calculations

Complex-band structure: a method to determine the off-resonant electron transport in oligomers

We validate that off-resonant electron transport across {\it ultra-short} oligomer molecular junctions is characterised by a conductance which decays exponentially with length, and we discuss a method to determine the damping factor via the energy spectrum of a periodic structure as a function of complex wavevector. An exact mapping to the complex wavevector is demonstrated by first-principle-based calculations of: a) the conductance of molecular junctions of phenyl-ethynylene wires covalently bonded to graphitic ribbons as a function of the bridge length, and b) the complex-band structure of poly-phenyl-ethynylene.Comment: version to appear in Chem Phys Lett; 8 pages, 4 figures; minor changes to the 06/08/03 submission (nomenclature and added concluding remark