Search for: All records

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. 

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. 

The HIV-1 capsid protein (CA) assembles into a conical shell during viral maturation, encasing and protecting the viral RNA genome. The C-terminal domain (CTD) of the two-domain capsid protein dimerizes, and this dimer connects individual chains in the mature capsid lattice. Previous NMR studies have shown that different dimer arrangements can be formed by isolated capsid protein chains and in assembled capsid lattices; however, the dynamics and functional relevance of these alternate dimers are unknown. To explore the conformational landscape of the CA-CTD dimer, we carried out atomistic molecular dynamics simulations using the weighted ensemble path sampling strategy, generating an ensemble of conformations. Focusing on the two dimer forms previously observed via solution NMR, we refined the conformational ensemble to highlight two metastable states using a Markov state model. Experimentally, we measured the interconversion rates between the two alternate dimers using¹⁹F NMR, and these rates showed good agreement with the interconversion rates derived from the simulations. After identifying the key interactions that distinguish the dimer states, the alternate dimer was further experimentally verified through disulfide crosslinking. Our results demonstrate the advantages of pairing weighted ensemble path sampling with¹⁹F NMR to gain atomistic insights into the hidden dimer state of the HIV-1 capsid protein.