An official website of the United States government
Abstract Electron paramagnetic resonance (EPR) has become a powerful probe of conformational heterogeneity and dynamics of biomolecules. In this Review, we discuss different computational modeling techniques that enrich the interpretation of EPR measurements of dynamics or distance restraints. A variety of spin labels are surveyed to provide a background for the discussion of modeling tools. Molecular dynamics (MD) simulations of models containing spin labels provide dynamical properties of biomolecules and their labels. These simulations can be used to predict EPR spectra, sample stable conformations and sample rotameric preferences of label sidechains. For molecular motions longer than milliseconds, enhanced sampling strategies and de novo prediction software incorporating or validated by EPR measurements are able to efficiently refine or predict protein conformations, respectively. To sample large‐amplitude conformational transition, a coarse‐grained or an atomistic weighted ensemble (WE) strategy can be guided with EPR insights. Looking forward, we anticipate an integrative strategy for efficient sampling of alternate conformations by de novo predictions, followed by validations by systematic EPR measurements and MD simulations. Continuous pathways between alternate states can be further sampled by WE‐MD including all intermediate states.
more »
« less
Direct observation of negative cooperativity in a detoxification enzyme at the atomic level by Electron Paramagnetic Resonance spectroscopy and simulation
Abstract The catalytic activity of human glutathione S‐transferase A1‐1 (hGSTA1‐1), a homodimeric detoxification enzyme, is dependent on the conformational dynamics of a key C‐terminal helix α9 in each monomer. However, the structural details of how the two monomers interact upon binding of substrates is not well understood and the structure of the ligand‐free state of the hGSTA1‐1 homodimer has not been resolved. Here, we used a combination of electron paramagnetic resonance (EPR) distance measurements and weighted ensemble (WE) simulations to characterize the conformational ensemble of the ligand‐free state at the atomic level. EPR measurements reveal a broad distance distribution between a pair of Cu(II) labels in the ligand‐free state that gradually shifts and narrows as a function of increasing ligand concentration. These shifts suggest changes in the relative positioning of the two α9 helices upon ligand binding. WE simulations generated unbiased pathways for the seconds‐timescale transition between alternate states of the enzyme, leading to the generation of atomically detailed structures of the ligand‐free state. Notably, the simulations provide direct observations of negative cooperativity between the monomers of hGSTA1‐1, which involve the mutually exclusive docking of α9 in each monomer as a lid over the active site. We identify key interactions between residues that lead to this negative cooperativity. Negative cooperativity may be essential for interaction of hGSTA1‐1 with a wide variety of toxic substrates and their subsequent neutralization. More broadly, this work demonstrates the power of integrating EPR distances with WE rare‐events sampling strategy to gain mechanistic information on protein function at the atomic level.
more »
« less
- Award ID(s):
- 2112871
- PAR ID:
- 10466021
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Protein Science
- Volume:
- 32
- Issue:
- 10
- ISSN:
- 0961-8368
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)We present a quantitative study comparing the binding of 4-methoxypyridine, MeOPy, ligand to Co( ii )octaethylporphyrin, CoOEP, at the phenyloctane/HOPG interface and in toluene solution. Scanning tunneling microscopy (STM) was used to study the ligand binding to the porphyrin receptors adsorbed on graphite. Electronic spectroscopy was employed for examining this process in fluid solution. The on surface coordination reaction was completely reversible and followed a simple Langmuir adsorption isotherm. Ligand affinities (or Δ G ) for the binding processes in the two different chemical environments were determined from the respective equilibrium constants. The free energy value of −13.0 ± 0.3 kJ mol −1 for the ligation reaction of MeOPy to CoOEP at the solution/HOPG interface is less negative than the Δ G for cobalt porphyrin complexed to the ligand in solution, −16.8 ± 0.2 kJ mol −1 . This result indicates that the MeOPy–CoOEP complex is more stable in solution than on the surface. Additional thermodynamic values for the formation of the surface ligated species (Δ H c = −50 kJ mol −1 and Δ S c = −120 J mol −1 ) were extracted from temperature dependent STM measurements. Density functional computational methods were also employed to explore the energetics of both the solution and surface reactions. At high concentrations of MeOPy the monolayer was observed to be stripped from the surface. Computational results indicate that this is not because of a reduction in adsorption energy of the MeOPy–CoOEP complex. Nearest neighbor analysis of the MeOPy–CoOEP in the STM images revealed positive cooperative ligand binding behavior. Our studies bring new insights to the general principles of affinity and cooperativity in the ligand–receptor interactions at the solution/solid interface. Future applications of STM will pave the way for new strategies designing highly functional multisite receptor systems for sensing, catalysis, and pharmacological applications.more » « less
-
Abstract 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.more » « less
-
null (Ed.)Abstract We show how equilibrium binding curves of receptor homodimers can be expressed as rational polynomial functions of the equilibrium binding curves of the constituent monomers, without approximation and without assuming independence of receptor monomers. Using a distinguished spanning tree construction for reduced graph powers, the method properly accounts for thermodynamic constraints and allosteric interactions between receptor monomers (i.e. conformational coupling). The method is completely general; it begins with an arbitrary undirected graph representing the topology of a monomer state-transition diagram and ends with an algebraic expression for the equilibrium binding curve of a receptor oligomer composed of two or more identical and indistinguishable monomers. Several specific examples are analysed, including guanine nucleotide-binding protein-coupled receptor dimers and tetramers composed of multiple ‘ternary complex’ monomers.more » « less
-
Thiamine pyrophosphate (TPP) riboswitches regulate thiamine metabolism by inhibiting the translation of enzymes essential to thiamine synthesis pathways upon binding to thiamine pyrophosphate in cells across all domains of life. Recent work on the Arabidopsis thaliana TPP riboswitch suggests a multistep TPP binding process involving multiple riboswitch configurational ensembles and Mg 2+ dependence underlies the mechanism of TPP recognition and subsequent transition to the expression-inhibiting state of the aptamer domain followed by changes in the expression platform. However, details of the relationship between TPP riboswitch conformational changes and interactions with TPP and Mg 2+ in the aptamer domain constituting this mechanism are unknown. Therefore, we integrated single-molecule multiparameter fluorescence and force spectroscopy with atomistic molecular dynamics simulations and found that conformational transitions within the aptamer domain's sensor helices associated with TPP and Mg 2+ ligand binding occurred between at least five different ensembles on timescales ranging from µs to ms. These dynamics are orders of magnitude faster than the 10 sec-timescale folding kinetics associated with expression-state switching in the switch helix. Together, our results show that a TPP and Mg 2+ dependent mechanism determines dynamic configurational state ensemble switching of the aptamer domain's sensor helices that regulate the switch helix's stability, which ultimately may lead to the expression-inhibiting state of the riboswitch. Additionally, we propose that two pathways exist for ligand recognition and that this mechanism underlies a kinetic rheostat-like behavior of the Arabidopsis thaliana TPP riboswitch.more » « less
