Double Electron-electron Resonance Distance Distributions Research Articles

Copper coordination states affect the flexibility of copper Metallochaperone Atox1: Insights from molecular dynamics simulations.

Copper is an essential element in nature but in excess, it is toxic to the living cell. The human metallochaperone Atox1 participates in copper homeostasis and is responsible for copper transmission. In a previous multiscale simulation study, we noticed a change in the coordination state of the Cu(I) ion, from 4 bound cysteine residues to 3, in agreement with earlier studies. Here, we perform and analyze classical molecular dynamic simulations of various coordination states: 2, 3, and 4. The main observation is an increase in protein flexibility as a result of a decrease in the coordination state. In addition, we identified several populated conformations that correlate well with double electron-electron resonance distance distributions or an X-ray structure of Cu(I)-bound Atox1. We suggest that the increased flexibility might benefit the process of ion transmission between interacting proteins. Further experiments can scrutinize this hypothesis and shed additional light on the mechanism of action of Atox1.

DEER-PREdict: Software for efficient calculation of spin-labeling EPR and NMR data from conformational ensembles.

Owing to their plasticity, intrinsically disordered and multidomain proteins require descriptions based on multiple conformations, thus calling for techniques and analysis tools that are capable of dealing with conformational ensembles rather than a single protein structure. Here, we introduce DEER-PREdict, a software program to predict Double Electron-Electron Resonance distance distributions as well as Paramagnetic Relaxation Enhancement rates from ensembles of protein conformations. DEER-PREdict uses an established rotamer library approach to describe the paramagnetic probes which are bound covalently to the protein.DEER-PREdict has been designed to operate efficiently on large conformational ensembles, such as those generated by molecular dynamics simulation, to facilitate the validation or refinement of molecular models as well as the interpretation of experimental data. The performance and accuracy of the software is demonstrated with experimentally characterized protein systems: HIV-1 protease, T4 Lysozyme and Acyl-CoA-binding protein. DEER-PREdict is open source (GPLv3) and available at github.com/KULL-Centre/DEERpredict and as a Python PyPI package pypi.org/project/DEERPREdict.

Characterization of the ExoU activation mechanism using EPR and integrative modeling

ExoU, a type III secreted phospholipase effector of Pseudomonas aeruginosa, serves as a prototype to model large, dynamic, membrane-associated proteins. ExoU is synergistically activated by interactions with membrane lipids and ubiquitin. To dissect the activation mechanism, structural homology was used to identify an unstructured loop of approximately 20 residues in the ExoU amino acid sequence. Mutational analyses indicate the importance of specific loop amino acid residues in mediating catalytic activity. Engineered disulfide cross-links show that loop movement is required for activation. Site directed spin labeling EPR and DEER (double electron–electron resonance) studies of apo and holo states demonstrate local conformational changes at specific sites within the loop and a conformational shift of the loop during activation. These data are consistent with the formation of a substrate-binding pocket providing access to the catalytic site. DEER distance distributions were used as constraints in RosettaDEER to construct ensemble models of the loop in both apo and holo states, significantly extending the range for modeling a conformationally dynamic loop.

Cu2+-based distance measurements by pulsed EPR provide distance constraints for DNA backbone conformations in solution.

Electron paramagnetic resonance (EPR) has become an important tool to probe conformational changes in nucleic acids. An array of EPR labels for nucleic acids are available, but they often come at the cost of long tethers, are dependent on the presence of a particular nucleotide or can be placed only at the termini. Site directed incorporation of Cu2+-chelated to a ligand, 2,2′dipicolylamine (DPA) is potentially an attractive strategy for site-specific, nucleotide independent Cu2+-labelling in DNA. To fully understand the potential of this label, we undertook a systematic and detailed analysis of the Cu2+-DPA motif using EPR and molecular dynamics (MD) simulations. We used continuous wave EPR experiments to characterize Cu2+ binding to DPA as well as optimize Cu2+ loading conditions. We performed double electron-electron resonance (DEER) experiments at two frequencies to elucidate orientational selectivity effects. Furthermore, comparison of DEER and MD simulated distance distributions reveal a remarkable agreement in the most probable distances. The results illustrate the efficacy of the Cu2+-DPA in reporting on DNA backbone conformations for sufficiently long base pair separations. This labelling strategy can serve as an important tool for probing conformational changes in DNA upon interaction with other macromolecules.

Conformational heterogeneity of the allosteric drug and metabolite (ADaM) site in AMP-activated protein kinase (AMPK)

AMP-activated protein kinase (AMPK) is a master regulator of energy homeostasis and a promising drug target for managing metabolic diseases such as type 2 diabetes. Many pharmacological AMPK activators, and possibly unidentified physiological metabolites, bind to the allosteric drug and metabolite (ADaM) site at the interface between the kinase domain (KD) in the α-subunit and the carbohydrate-binding module (CBM) in the β-subunit. Here, using double electron-electron resonance (DEER) spectroscopy, we demonstrate that the CBM-KD interaction is partially dissociated and the interface highly disordered in the absence of pharmacological ADaM site activators as inferred from a low depth of modulation and broad DEER distance distributions. ADaM site ligands such as 991, and to a lesser degree phosphorylation, stabilize the KD-CBM association and strikingly reduce conformational heterogeneity in the ADaM site. Our findings that the ADaM site, formed by the KD-CBM interaction, can be modulated by diverse ligands and by phosphorylation suggest that it may function as a hub for integrating regulatory signals.

Two-Dimensional Distance Correlation Maps from Pulsed Triple Electron Resonance (TRIER) on Proteins with Three Paramagnetic Centers

Recently, we introduced the pulsed triple electron resonance (TRIER) experiment, which correlates dipolar frequencies of molecules with three electron spins. These correlation patterns contain valuable information: in combination with double electron–electron resonance (DEER) they allow to interpret distance distributions of biological systems that exist in more than one conformation. Together with an improved sequence and new data processing, we were now for the first time able to obtain two-dimensional distance correlation maps of the previously investigated model compounds as well as of spin-labeled proteins. For this we applied two-dimensional approximate Pake transformation to TRIER data. This enabled us to get distance correlation plots from two triple-labeled protein samples that were in good agreement with DEER data and simulations. With such information it should then be possible to assign peaks in DEER distance distributions for macromolecules that can exist in more than one conformation.

Probing Conformational Motions Underlying Anesthetic Drug Actions in a Ligand-Gated Ion Channel

General anesthetics and ethanol exert their functional effects on pentameric ligand-gated ion channels (pLGICs) by binding to distinct allosteric sites that are far from the agonist binding sites. While recent high-resolution structures of prokaryotic and eukaryotic pLGICs, solved in the presence of these drugs have illuminated their binding sites, the structural mechanisms underlying how these drugs enhance or inhibit pLGIC function are largely unknown. Recently, a single residue variant in the bacterial pLGIC, GLIC, was identified that is potentiated by ethanol (PNAS 108:12149-12154, 2011). Taking advantage of this GLIC ethanol-sensitive mutant, we are using site directed spin labeling electron paramagnetic resonance (SDSL EPR) spectroscopy to study motions underlying ethanol potentiation of GLIC currents. Previous studies have shown that agonist-mediated channel opening is associated with an inward tilting of the extracellular binding domain (ECD). Here, we examined if, and how, the ECD moves in response to ethanol binding in the transmembrane channel domain. Individual cysteine residues were introduced in the ECD of an ethanol-sensitive GLIC construct (F14’A) at the top of β-strands (β1, 2, 5, 6, 8), labeled with the spin probe MTSL and ethanol-induced motions were measured using double electron-electron resonance (DEER) spectroscopy and compared to agonist (proton)-induced gating motions. DEER distance distributions for spin labels attached to N19C (β1), S93C (β5) and Q101C (β6) revealed that the attached probes moved closer together at pH 3.0, indicating that the inner β-sheet of GLIC ECD tilts inwards upon agonist activation. Interestingly, DEER distance distributions for N19 (β1) and S93 (β5) showed no significant changes in the presence of ethanol, indicating that structural mechanisms underlying ethanol allosteric drug modulation of pLGIC function are different than agonist-mediated motions that lead to channel activation.

Guiding Double Electron Electron Resonance Experiments using Information Theory and Molecular Dynamics Simulations

Flexible proteins present a unique set of challenges for refinement with either computational or experimental techniques: because these proteins occupy a large conformational space, standard molecular dynamics (MD) simulations can be prohibitively expensive, and high-resolution experimental techniques like x-ray crystallography intentionally capture only small set of low energy states. A new approach is therefore needed to capture the conformational ensembles of flexible proteins. Double electron-electron resonance (DEER) spectroscopy, an experimental technique that measures distance distributions between pairs of spin-labeled residues on a protein, provides quantitative information about less populated but still potentially important conformational states. However, because DEER experiments are costly and time consuming, only a small subset of pairs can be measured, often chosen arbitrarily. If instead of selecting an arbitrary or random subset, we select a subset of pairs such that they provide the maximum amount of information about the conformational ensemble, DEER could become an even more powerful tool for refinement. We have therefore developed 1) a new information-theoretic metric that ranks pairs of residues based on how well they refine a conformational ensemble and 2) a new methodology which identifies a set of highly informative pairs that perform well under this metric. Specifically, we incorporate DEER distance distributions into restrained-ensemble MD simulations, utilize the data from these simulations to identify a set of maximally-informative and minimally-redundant pairs, measure the pairs with DEER, and demonstrate the use of these pairs in refining the conformational ensemble of a flexible protein.

Bayesian Statistical Methods in the Analysis of DEER Data

Double Electron-Electron Resonance (DEER), also known as Pulsed Electron Double Resonance (PELDOR), is an increasingly important pulsed Electron Paramagnetic Resonance (EPR) technique for long distance (20-80 A) measurements in biomacromolecules. Notably, DEER reports a distribution of distances, P(r), corresponding to protein conformational ensembles in frozen glassy solutions. Due to the ill-conditioned, inverse nature of DEER data analysis, numerical regularization is required to obtain P(r). Further, a robust method for quantifying uncertainty in DEER distance distributions does not currently exist, complicating interpretation. Here, we compare methods of regularization parameter selection and demonstrate Bayesian statistical methods to quantify uncertainty in DEER distance distributions arising from measurement noise. We generate corresponding Bayesian credible intervals for visualization and interpretation.

Radial Tilting of the Extracellular Domain of GLIC Revealed by EPR Spectroscopy

Signaling in the brain depends on rapid opening and closing of pentameric ligand-gated ion channels (pLGICs). Defects in these proteins are linked to a variety of diseases including myasthenia, epilepsy and sleep-disorders. In addition, these proteins are therapeutic targets for various drugs including anti-convulsants, anxiolytics, anesthetics and sedatives. While recent crystal structures of prokaryotic and eukaryotic pLGICs solved in detergent micelles have shed light on the molecular architecture of these proteins, describing their conformational dynamics in physiological lipids is essential for understanding their function. Crystal structures of GLIC and the glutamate-activated chloride channel predict that agonist-mediated channel opening is associated with a quaternary twisting and a radial un-blooming (inward tilting) of the extracellular binding domain (ECD). Using site directed spin labeling electron paramagnetic resonance (SDSL EPR) spectroscopy and functional GLIC channels reconstituted in liposomes, we are testing these gating predictions. Residues located at the top of β-strands (β1, 2, 5, 6, 8) were mutated to cysteine, labeled with the spin probe MTSL and proton-induced gating motions were measured using continuous wave (CW) and double electron-electron resonance (DEER) spectroscopy. CW EPR spectroscopy data showed the spin label attached to I92C, S93C, V94C (β5), T99C and Q101C (β6) moved into a slightly more constrained/immobile environment at pH 3.0 as compared to pH 7.6. DEER distance distributions for K48 (β2), S93 (β5), T99 and Q101 (β6) revealed the probes moved closer together at pH 3.0, suggesting that the inner β-sheet of the extracellular domain of GLIC tilt inwards upon agonist activation.

Ensemble-Biased Metadynamics: A Molecular Simulation Method to Sample Experimental Distributions

We introduce an enhanced-sampling method for molecular dynamics (MD) simulations referred to as ensemble-biased metadynamics (EBMetaD). The method biases a conventional MD simulation to sample a molecular ensemble that is consistent with one or more probability distributions known a priori, e.g., experimental intramolecular distance distributions obtained by double electron-electron resonance or other spectroscopic techniques. To this end, EBMetaD adds an adaptive biasing potential throughout the simulation that discourages sampling of configurations inconsistent with the target probability distributions. The bias introduced is the minimum necessary to fulfill the target distributions, i.e., EBMetaD satisfies the maximum-entropy principle. Unlike other methods, EBMetaD does not require multiple simulation replicas or the introduction of Lagrange multipliers, and is therefore computationally efficient and straightforward in practice. We demonstrate the performance and accuracy of the method for a model system as well as for spin-labeled T4 lysozyme in explicit water, and show how EBMetaD reproduces three double electron-electron resonance distance distributions concurrently within a few tens of nanoseconds of simulation time. EBMetaD is integrated in the open-source PLUMED plug-in (www.plumed-code.org), and can be therefore readily used with multiple MD engines.

Allosteric Regulation of the Cyclic Nucleotide-Binding Domain in HCN Channels

Hyperpolarization-activated cyclic nucleotide-gated (HCN) ion channels play an important role in regulating pacemaking activity in the heart and brain. They are regulated by the binding of cyclic nucleotides to a conserved, intracellular cyclic nucleotide-binding domain (CNBD). Binding of cyclic nucleotides increases the rate and extent of activation of the channels and shifts channel activation to less hyperpolarized voltages. In intact channels, cAMP and cGMP are full agonists and cCMP is a partial agonist. We use double electron-electron resonance (DEER) to study the conformational change associated with the binding of these three different cyclic nucleotide species to the CNBD. We find that conformational changes associated with ligand binding vary depending on the cyclic nucleotide bound and the location in the CNBD. We use the restrained ensemble MD (re-MD) simulations method to generate structural models integrating the complete set of experimentally measured data from DEER distance distribution histograms that describe the separation between pairs of spin labels attached to the CNBD. Our results indicate that, in the B-helix, binding of cGMP or cCMP has an effect intermediate to that of cAMP. However, in the C-helix cAMP and cCMP exhibit similar effects, but cGMP produces an intermediate effect between the apo and bound cAMP/cCMP DEER distance distributions. These data provide an interesting lens for studying allostery in HCN channels and indicate that the mechanism of protein allostery in the CNBD varies for different cyclic nucleotides.

Simulation vs. Reality: A Comparison of In Silico Distance Predictions with DEER and FRET Measurements

Site specific incorporation of molecular probes such as fluorescent- and nitroxide spin-labels into biomolecules, and subsequent analysis by Förster resonance energy transfer (FRET) and double electron-electron resonance (DEER) can elucidate the distance and distance-changes between the probes. However, the probes have an intrinsic conformational flexibility due to the linker by which they are conjugated to the biomolecule. This property minimizes the influence of the label side chain on the structure of the target molecule, but complicates the direct correlation of the experimental inter-label distances with the macromolecular structure or changes thereof. Simulation methods that account for the conformational flexibility and orientation of the probe(s) can be helpful in overcoming this problem. We performed distance measurements using FRET and DEER and explored different simulation techniques to predict inter-label distances using the Rpo4/7 stalk module of the M. jannaschii RNA polymerase. This is a suitable model system because it is rigid and a high-resolution X-ray structure is available. The conformations of the fluorescent labels and nitroxide spin labels on Rpo4/7 were modeled using in vacuo molecular dynamics simulations (MD) and a stochastic Monte Carlo sampling approach. For the nitroxide probes we also performed MD simulations with explicit water and carried out a rotamer library analysis. Our results show that the Monte Carlo simulations are in better agreement with experiments than the MD simulations and the rotamer library approach results in plausible distance predictions. Because the latter is the least computationally demanding of the methods we have explored, and is readily available to many researchers, it prevails as the method of choice for the interpretation of DEER distance distributions.