Spatial domain organization in the HIV-1 reverse transcriptase p66 homodimer precursor probed by double electron-electron resonance EPR

Interplay between protease and reverse transcriptase dimerization in a model HIV-1 polyprotein

The Gag-Pol polyprotein in human immunodeficiency virus type I (HIV-1) encodes enzymes that are essential for virus replication: protease (PR), reverse transcriptase (RT), and integrase (IN). The mature forms of PR, RT and IN are homodimer, heterodimer and tetramer, respectively. The precise mechanism underlying the formation of dimer or tetramer is not yet understood. Here, to gain insight into the dimerization of PR and RT in the precursor, we prepared a model precursor, PR-RT, incorporating an inactivating mutation at the PR active site, D25A, and including two residues in the p6* region, fused to a SUMO-tag, at the N-terminus of the PR region. We also prepared two mutants of PR-RT containing a dimer dissociation mutation either in the PR region, PR(T26A)-RT, or in the RT region, PR-RT(W401A). Size exclusion chromatography showed both monomer and dimer fractions in PR-RT and PR(T26A)-RT, but only monomer in PR-RT(W401A). SEC experiments of PR-RT in the presence of protease inhibitor, darunavir, significantly enhanced the dimerization. Additionally, SEC results suggest an estimated PR-RT dimer dissociation constant that is higher than that of the mature RT heterodimer, p66/p51, but slightly lower than the premature RT homodimer, p66/p66. Reverse transcriptase assays and RT maturation assays were performed as tools to assess the effects of the PR dimer-interface on these functions. Our results consistently indicate that the RT dimer-interface plays a crucial role in the dimerization in PR-RT, whereas the PR dimer-interface has a lesser role.

Xplor-NIH: Better parameters and protocols for NMR protein structure determination

The present work describes an update to the protein covalent geometry and atomic radii parameters in the Xplor-NIH biomolecular structure determination package. In combination with an improved treatment of selected non-bonded interactions between atoms three bonds apart, such as those involving methyl hydrogens, and a previously developed term that affects the system's gyration volume, the new parameters are tested using structure calculations on 30 proteins with restraints derived from nuclear magnetic resonance data. Using modern structure validation criteria, including several formally adopted by the Protein Data Bank, and a clear measure of structural accuracy, the results show superior performance relative to previous Xplor-NIH implementations. Additionally, the Xplor-NIH structures compare favorably against originally determined NMR models.

Integrating Electron Paramagnetic Resonance Spectroscopy and Computational Modeling to Measure Protein Structure and Dynamics

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.

Direct observation of negative cooperativity in a detoxification enzyme at the atomic level by Electron Paramagnetic Resonance spectroscopy and simulation

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.

Potent targeted activator of cell kill molecules eliminate cells expressing HIV-1

Antiretroviral therapy inhibits HIV-1 replication but is not curative due to establishment of a persistent reservoir after virus integration into the host genome. Reservoir reduction is therefore an important HIV-1 cure strategy. Some HIV-1 nonnucleoside reverse transcriptase inhibitors induce HIV-1 selective cytotoxicity in vitro but require concentrations far exceeding approved dosages. Focusing on this secondary activity, we found bifunctional compounds with HIV-1-infected cell kill potency at clinically achievable concentrations. These targeted activator of cell kill (TACK) molecules bind the reverse transcriptase-p66 domain of monomeric Gag-Pol and act as allosteric modulators to accelerate dimerization, resulting in HIV-1⁺ cell death through premature intracellular viral protease activation. TACK molecules retain potent antiviral activity and selectively eliminate infected CD4⁺ T cells isolated from people living with HIV-1, supporting an immune-independent clearance strategy.

Quantitative analysis of sterol-modulated monomer-dimer equilibrium of the β1-adrenergic receptor by DEER spectroscopy

G protein-coupled receptors (GPCR) activate numerous intracellular signaling pathways. The oligomerization properties of GPCRs, and hence their cellular functions, may be modulated by various components within the cell membrane (such as the presence of cholesterol). Modulation may occur directly via specific interaction with the GPCR or indirectly by affecting the physical properties of the membrane. Here, we use pulsed Q-band double electron-electron resonance (DEER) spectroscopy to probe distances between R1 nitroxide spin labels attached to Cys163 and Cys344 of the β1-adrenergic receptor (β1AR) in n-dodecyl-β-D-maltoside micelles upon titration with two soluble cholesterol analogs, cholesteryl hemisuccinate (CHS) and sodium cholate. The former, like cholesterol, inserts itself into the lipid membrane, parallel to the phospholipid chains; the latter is aligned parallel to the surface of membranes. Global quantitative analysis of DEER echo curves upon titration of spin-labeled β1AR with CHS and sodium cholate reveal the following: CHS binds specifically to the β1AR monomer at a site close to the Cys163-R1 spin label with an equilibrium dissociation constant [Formula: see text] ~1.4 ± 0.4 mM. While no direct binding of sodium cholate to the β1AR receptor was observed by DEER, sodium cholate induces specific β1AR dimerization ([Formula: see text] ~35 ± 6 mM and a Hill coefficient n ~ 2.5 ± 0.4) with intersubunit contacts between transmembrane helices 1 and 2 and helix 8. Analysis of the DEER data obtained upon the addition of CHS to the β1AR dimer in the presence of excess cholate results in dimer dissociation with species occupancies as predicted from the individual KD values.

Human endogenous retrovirus-K (HERV-K) reverse transcriptase (RT) structure and biochemistry reveals remarkable similarities to HIV-1 RT and opportunities for HERV-K-specific inhibition

Human endogenous retroviruses (HERVs) comprise nearly 8% of the human genome and are derived from ancient integrations of retroviruses into the germline. The biology of HERVs is poorly defined, but there is accumulating evidence supporting pathological roles in diverse diseases, such as cancer, autoimmune, and neurodegenerative diseases. Functional proteins are produced by HERV-encoded genes, including reverse transcriptases (RTs), which could be a contributor to the pathology attributed to aberrant HERV-K expression. To facilitate the discovery and development of HERV-K RT potent and selective inhibitors, we expressed active HERV-K RT and determined the crystal structure of a ternary complex of this enzyme with a double-stranded DNA substrate. We demonstrate a range of RT inhibition with antiretroviral nucleotide analogs, while classic nonnucleoside analogs do not inhibit HERV-K RT. Detailed comparisons of HERV-K RT with other known RTs demonstrate similarities to diverse RT families and a striking similarity to the HIV-1 RT asymmetric heterodimer. Our analysis further reveals opportunities for selective HERV-K RT inhibition.

Quantitative Agreement between Conformational Substates of Holo Calcium-Loaded Calmodulin Detected by Double Electron-Electron Resonance EPR and Predicted by Molecular Dynamics Simulations

Calcium-loaded calmodulin (CaM/4Ca²⁺) comprises two domains that undergo rigid body reorientation from a predominantly extended conformation to a compact one upon binding target peptides. A recent replica-exchange molecular dynamics (MD) simulation on holo CaM/4Ca²⁺ suggested the existence of distinct structural clusters (substates) along the path from extended to compact conformers in the absence of substrates. Here, we experimentally demonstrate the existence of CaM/4Ca²⁺ substates trapped in local minima by three freezing/annealing regimes (slow, 40 s; intermediate, 1.5 s; fast, 0.5 ms) using pulsed Q-band double electron-electron resonance (DEER) EPR spectroscopy to measure interdomain distances between nitroxide spin-labels positioned at A17C and A128C in the N- and C-terminal domains, respectively. The DEER echo curves were directly fit to population-optimized P(r) pairwise distance distributions calculated from the coordinates of the MD clusters and compact crystal structure. DEER data on fully deuterated CaM/4Ca²⁺ were acquired at multiple values of the second echo period (10-35 μs) and analyzed globally to eliminate instrumental and overfitting artifacts and ensure accurate populations, peak positions, and widths. The DEER data for all three freezing regimes are quantitatively accounted for within experimental error by 5-6 distinct conformers comprising a predominantly populated extended form (60-75%) and progressively more compact states whose populations decrease as the degree of compactness increases. The shortest interdomain separation is found in the compact crystal structure, which has an occupancy of 4-6%. Thus, CaM/4Ca²⁺ samples high energy local minima comprising a few discrete substates of increasing compactness in a rugged energy landscape.

An optimal acquisition scheme for Q-band EPR distance measurements using Cu2+-based protein labels

Recent advances in site-directed Cu²⁺ labeling of proteins and nucleic acids have added an attractive new methodology to measure the structure-function relationship in biomolecules. Despite the promise, accessing the higher sensitivity of Q-band Double Electron Electron Resonance (DEER) has been challenging for Cu²⁺ labels designed for proteins. Q-band DEER experiments on this label typically require many measurements at different magnetic fields, since the pulses can excite only a few orientations at a given magnetic field. Herein, we analyze such orientational effects through simulations and show that three DEER measurements, at strategically selected magnetic fields, are generally sufficient to acquire an orientational-averaged DEER time trace for this spin label at Q-band. The modeling results are experimentally verified on Cu²⁺ labeled human glutathione S-transferase (hGSTA1-1). The DEER distance distribution measured at the Q-band shows good agreement with the distance distribution sampled by molecular dynamics (MD) simulations and X-band experiments. The concordance of MD sampled distances and experimentally measured distances adds growing evidence that MD simulations can accurately predict distances for the Cu²⁺ labels, which remains a key bottleneck for the commonly used nitroxide label. In all, this minimal collection scheme reduces data collection time by as much as six-fold and is generally applicable to many octahedrally coordinated Cu²⁺ systems. Furthermore, the concepts presented here may be applied to other metals and pulsed EPR experiments.

Time-resolved DEER EPR and solid-state NMR afford kinetic and structural elucidation of substrate binding to Ca2+-ligated calmodulin

Complex formation between calmodulin and target proteins underlies numerous calcium signaling processes in biology, yet structural and mechanistic details, which entail major conformational changes in both calmodulin and its substrates, have been unclear. We show that a combination of time-resolved electron paramagnetic and NMR measurements can elucidate the molecular mechanism, at the quantitative kinetic and structural levels, of the binding pathway of a peptide substrate from skeletal muscle myosin light-chain kinase to calcium-loaded calmodulin. The mechanism involves coupled folding and binding and comprises a bifurcated process, with rapid, direct complex formation when the peptide interacts first with the C-terminal domain of calmodulin or a slower, two-step complex formation when the peptide interacts initially with the N-terminal domain.

Recent advances in rapid mixing and freeze quenching have opened the path for time-resolved electron paramagnetic resonance (EPR)-based double electron-electron resonance (DEER) and solid-state NMR of protein–substrate interactions. DEER, in conjunction with phase memory time filtering to quantitatively extract species populations, permits monitoring time-dependent probability distance distributions between pairs of spin labels, while solid-state NMR provides quantitative residue-specific information on the appearance of structural order and the development of intermolecular contacts between substrate and protein. Here, we demonstrate the power of these combined approaches to unravel the kinetic and structural pathways in the binding of the intrinsically disordered peptide substrate (M13) derived from myosin light-chain kinase to the universal eukaryotic calcium regulator, calmodulin. Global kinetic analysis of the data reveals coupled folding and binding of the peptide associated with large spatial rearrangements of the two domains of calmodulin. The initial binding events involve a bifurcating pathway in which the M13 peptide associates via either its N- or C-terminal regions with the C- or N-terminal domains, respectively, of calmodulin/4Ca²⁺ to yield two extended “encounter” complexes, states A and A*, without conformational ordering of M13. State A is immediately converted to the final compact complex, state C, on a timescale τ ≤ 600 μs. State A*, however, only reaches the final complex via a collapsed intermediate B (τ ∼ 1.5 to 2.5 ms), in which the peptide is only partially ordered and not all intermolecular contacts are formed. State B then undergoes a relatively slow (τ ∼ 7 to 18 ms) conformational rearrangement to state C.

Reconstruction of Coupled Intra- and Interdomain Protein Motion from Nuclear and Electron Magnetic Resonance

Proteins composed of multiple domains allow for structural heterogeneity and interdomain dynamics that may be vital for function. Intradomain structures and dynamics can influence interdomain conformations and vice versa. However, no established structure determination method is currently available that can probe the coupling of these motions. The protein Pin1 contains separate regulatory and catalytic domains that sample "extended" and "compact" states, and ligand binding changes this equilibrium. Ligand binding and interdomain distance have been shown to impact the activity of Pin1, suggesting interdomain allostery. In order to characterize the conformational equilibrium of Pin1, we describe a novel method to model the coupling between intra- and interdomain dynamics at atomic resolution using multistate ensembles. The method uses time-averaged nuclear magnetic resonance (NMR) restraints and double electron-electron resonance (DEER) data that resolve distance distributions. While the intradomain calculation is primarily driven by exact nuclear Overhauser enhancements (eNOEs), J couplings, and residual dipolar couplings (RDCs), the relative domain distribution is driven by paramagnetic relaxation enhancement (PREs), RDCs, interdomain NOEs, and DEER. Our data support a 70:30 population of the compact and extended states in apo Pin1. A multistate ensemble describes these conformations simultaneously, with distinct conformational differences located in the interdomain interface stabilizing the compact or extended states. We also describe correlated conformations between the catalytic site and interdomain interface that may explain allostery driven by interdomain contact.

Large Multidomain Protein NMR: HIV-1 Reverse Transcriptase Precursor in Solution

NMR studies of large proteins, over 100 kDa, in solution are technically challenging and, therefore, of considerable interest in the biophysics field. The challenge arises because the molecular tumbling of a protein in solution considerably slows as molecular mass increases, reducing the ability to detect resonances. In fact, the typical ¹H-¹³C or ¹H-¹⁵N correlation spectrum of a large protein, using a ¹³C- or ¹⁵N-uniformly labeled protein, shows severe line-broadening and signal overlap. Selective isotope labeling of methyl groups is a useful strategy to reduce these issues, however, the reduction in the number of signals that goes hand-in-hand with such a strategy is, in turn, disadvantageous for characterizing the overall features of the protein. When domain motion exists in large proteins, the domain motion differently affects backbone amide signals and methyl groups. Thus, the use of multiple NMR probes, such as ¹H, ¹⁹F, ¹³C, and ¹⁵N, is ideal to gain overall structural or dynamical information for large proteins. We discuss the utility of observing different NMR nuclei when characterizing a large protein, namely, the 66 kDa multi-domain HIV-1 reverse transcriptase that forms a homodimer in solution. Importantly, we present a biophysical approach, complemented by biochemical assays, to understand not only the homodimer, p66/p66, but also the conformational changes that contribute to its maturation to a heterodimer, p66/p51, upon HIV-1 protease cleavage.

Orientation and dynamics of Cu2+ based DNA labels from force field parameterized MD elucidates the relationship between EPR distance constraints and DNA backbone distances

Pulsed electron paramagnetic resonance (EPR) based distance measurements using the recently developed Cu2+-DPA label present a promising strategy for measuring DNA backbone distance constraints. Herein we develop force field parameters for Cu2+-DPA in order to understand the features of this label at an atomic level. We perform molecular dynamics (MD) simulations using the force field parameters of Cu2+-DPA on four different DNA duplexes. The distance between the Cu2+ centers, extracted from the 2 μs MD trajectories, agrees well with the experimental distance for all the duplexes. Further analyses of the trajectory provide insight into the orientation of the Cu2+-DPA inside the duplex that leads to such agreement with experiments. The MD results also illustrate the ability of the Cu2+-DPA to report on the DNA backbone distance constraints. Furthermore, measurement of fluctuations of individual residues showed that the flexibility of Cu2+-DPA in a DNA depends on the position of the label in the duplex, and a 2 μs MD simulation is not sufficient to fully capture the experimental distribution in some cases. Finally, the MD trajectories were utilized to understand the key aspects of the double electron electron resonance (DEER) results. The lack of orientational selectivity effects of the Cu2+-DPA at Q-band frequency is rationalized in terms of fluctuations in the Cu2+ coordination environment and rotameric fluctuations of the label linker. Overall, a combination of EPR and MD simulations based on the Cu2+-DPA labelling strategy can contribute towards understanding changes in DNA backbone conformations during protein-DNA interactions.

Probing the Interaction between HIV-1 Protease and the Homodimeric p66/p66' Reverse Transcriptase Precursor by Double Electron-Electron Resonance EPR Spectroscopy

Following excision from the Gag-Pol polyprotein, HIV-1 reverse transcriptase is released as an asymmetric homodimer comprising two p66 subunits that are structurally dissimilar but identical in amino acid sequence. Subsequent cleavage of the RNase H domain from only one of the subunits, denoted p66', results in the formation of the mature p66/p51 enzyme in which catalytic activity resides in the p66 subunit, and the p51 subunit (derived from p66') provides a supporting structural scaffold. Here, we probe the interaction of the p66/p66' asymmetric reverse transcriptase precursor with HIV-1 protease by pulsed Q-band double electron-electron resonance EPR spectroscopy to measure distances between nitroxide labels introduced at surface-engineered cysteine residues. The data suggest that the flexible, exposed linker between the RNaseH and connection domains in the open state of the p66' subunit binds to the active site of protease in a configuration that is similar to that of extended peptide substrates.

Submillisecond Freezing Permits Cryoprotectant-Free EPR Double Electron-Electron Resonance Spectroscopy

Double electron-electron resonance (DEER) EPR spectroscopy is a powerful method for obtaining distance distributions between pairs of engineered nitroxide spin-labels in proteins and other biological macromolecules. These measurements require the use of cryogenic temperatures (77 K or less) to prolong the phase memory relaxation time (T_m ) sufficiently to enable detection of a DEER echo curve. Generally, a cryoprotectant such as glycerol is added to protein samples to facilitate glass formation and avoid protein clustering (which can result in a large decrease in T_m ) during relatively slow flash freezing in liquid N₂ . However, cryoprotectants are osmolytes and can influence protein folding/unfolding equilibria, as well as species populations in weak multimeric systems. Here we show that submillisecond rapid freezing, achieved by high velocity spraying of the sample onto a rapidly spinning, liquid nitrogen cooled copper disc obviates the requirement for cryoprotectants and permits high quality DEER data to be obtained in absence of glycerol. We demonstrate this approach on five different protein systems: protein A, the metastable drkN SH3 domain, urea-unfolded drkN SH3, HIV-1 reverse transcriptase, and the transmembrane domain of HIV-1 gp41 in lipid bicelles.

Molecular Dynamics Simulations Based on Newly Developed Force Field Parameters for Cu2+ Spin Labels Provide Insights into Double-Histidine-Based Double Electron-Electron Resonance

Electron paramagnetic resonance (EPR) in combination with the recently developed double-histidine (dHis)-based Cu²⁺ spin labeling has provided valuable insights into protein structure and conformational dynamics. To relate sparse distance constraints measured by EPR to protein fluctuations in solution, modeling techniques are needed. In this work, we have developed force field parameters for Cu²⁺-nitrilotriacetic and Cu²⁺-iminodiacetic acid spin labels. We employed molecular dynamics (MD) simulations to capture the atomic-level details of dHis-labeled protein fluctuations. The interspin distances extracted from 200 ns MD trajectories show good agreement with the experimental results. The MD simulations also illustrate the dramatic rigidity of the Cu²⁺ labels compared to the standard nitroxide spin label. Further, the relative orientations between spin-labeled sites were measured to provide insight into the use of double electron-electron resonance (DEER) methods for such labels. The relative mean angles, as well as the standard deviations of the relative angles, agree well in general with the spectral simulations published previously. The fluctuations of relative orientations help rationalize why orientation selectivity effects are minimal at X-band frequencies, but observable at the Q-band for such labels. In summary, the results show that by combining the experimental results with MD simulations precise information about protein conformations as well as flexibility can be obtained.

Evolving understanding of HIV-1 reverse transcriptase structure, function, inhibition, and resistance

The essential role of reverse transcription in the HIV life cycle is illustrated by the fact that half of the ∼30 FDA-approved drugs for HIV treatment target HIV-1 reverse transcriptase (RT). Even though more than 160 structures of RT deposited in the Protein Data Bank (PDB) have revealed the molecular architecture of RT in great detail, some key states of RT function and inhibition remain still unknown. Recent structures of RT initiation complexes, RT poised for RNA hydrolysis, and RT with approved drugs and investigational compounds have provided a deeper understanding of RT function and inhibition, suggesting novel avenues for targeting this central enzyme of HIV.