Spin labeling and Double Electron-Electron Resonance (DEER) to Deconstruct Conformational Ensembles of HIV Protease

Natural Polymorphisms D60E and I62V Stabilize a Closed Conformation in HIV-1 Protease in the Absence of an Inhibitor or Substrate

HIV infection remains a global health issue plagued by drug resistance and virological failure. Natural polymorphisms (NPs) contained within several African and Brazilian protease (PR) variants have been shown to induce a conformational landscape of more closed conformations compared to the sequence of subtype B prevalent in North America and Western Europe. Here we demonstrate through experimental pulsed EPR distance measurements and molecular dynamic (MD) simulations that the two common NPs D60E and I62V found within subtypes F and H can induce a closed conformation when introduced into HIV-1PR subtype B. Specifically, D60E alters the conformation in subtype B through the formation of a salt bridge with residue K43 contained within the nexus between the flap and hinge region of the HIV-1 PR fold. On the other hand, I62V modulates the packing of the hydrophobic cluster of the cantilever and fulcrum, also resulting in a more closed conformation.

A guideline for the distance measurement plans of site-directed spin labels for structural prediction of nucleic acids

CONTEXT AND RESULTS

Site-directed spin labeling (SDSL) combined with electron paramagnetic resonance spectroscopy methods has been successfully used to predict the structures of nucleic acids. These methods measure the distances between spin labels yielding distance equations that are solved using numerical algorithms to provide one or several structural predictions. In this work, the minimum number of SDSL distance measurements and distance measurement types required to predict a unique nucleic acid structure were investigated. Our results indicate that at least six distance measurements should be obtained given that the distance measurements do not connect one SDSL on one arm with more than three SDSLs on the other arm. Moreover, there may be a preference for 1-to-1 SLs distance measurements rather than 1-to-many SLs as the latter was linked to undefined structures discussed in this study.

METHODS

Pairs of double-helical arms of nucleic acid were simulated using the finite element software Pro/ENGINEER (PTC Inc., Boston, MA). In each simulation, a specific SDSL distance measurement plan was adopted and the resulting structure was tested for movability. Immovable structures indicate that this plan will potentially result in a unique structural prediction of the nucleic acid. All the possible plans for SDSL distance measurements were investigated either by direct measurement or by extrapolation.

Molecular-level interplay between intrinsically disordered clients and Hsp90

Proteostasis is maintained by a network of molecular chaperones, a prominent member of which is the 90-kilodalton heat shock protein Hsp90. The chaperone function of Hsp90 has been extensively reviewed previously, emphasizing its ATPase activity and remodeling of folded client proteins. Experimental evidence implicating Hsp90 in neurodegenerative diseases has bolstered interest in the noncanonical chaperoning of intrinsically disordered protein (IDPs), however the interplay between Hsp90 and its disordered clients remains poorly understood. In this review we describe recent advances that have contributed to our understanding of the intricate mechanisms characterizing Hsp90-mediated chaperoning of the IDPs tau and α-synuclein and survey emerging insights into the modulation of the chaperone-client interplay in the context of neurodegeneration.

Simple Cryoprotectant-Free Method to Advance Pulsed Dipolar ESR Spectroscopy for Capturing Protein Conformational Ensembles

Double electron-electron resonance (DEER) is a powerful technique for studying protein conformations. To preserve the room-temperature ensemble, proteins are usually shock-frozen in liquid nitrogen prior to DEER measurements. The use of cryoprotectant additives is, therefore, necessary to ensure the formation of a vitrified state. Here, we present a simple modification of the freezing process using a flexible fused silica microcapillary, which increases the freezing rates and thus enables DEER measurement without the use of cryoprotectants. The Bid protein, which is highly sensitive to cryoprotectant additives, is used as a model. We show that DEER with the simple modification can successfully reveal the cold denaturation of Bid, which was not possible with the conventional DEER preparations. The DEER result reveals the nature of Bid folding. Our method advances DEER for capturing the chemically and thermally induced conformational changes of a protein in a cryoprotectant-free medium.

Deciphering Complex Mechanisms of Resistance and Loss of Potency through Coupled Molecular Dynamics and Machine Learning

Drug resistance threatens many critical therapeutics through mutations in the drug target. The molecular mechanisms by which combinations of mutations, especially those remote from the active site, alter drug binding to confer resistance are poorly understood and thus difficult to counteract. A machine learning strategy was developed that coupled parallel molecular dynamics simulations with experimental potency to identify specific conserved mechanisms underlying resistance. Physical features were extracted from the simulations, analyzed, and integrated into one consistent and interpretable elastic network model. To rigorously test this strategy, HIV-1 protease variants with diverse mutations were used, with potencies ranging from picomolar to micromolar to the drug darunavir. Feature reduction resulted in a model with four specific features that predicts for both the training and test sets inhibitor binding free energy within 1 kcal/mol of the experimental value over this entire range of potency. These predictive features are physically interpretable, as they vary specifically with affinity and diagonally transverse across the protease homodimer. This physics-based strategy of parallel molecular dynamics and machine learning captures mechanisms by which complex combinations of mutations confer resistance and identify critical features that serve as bellwethers of affinity, which will be critical in future drug design.

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.

Structure and regulation of the BsYetJ calcium channel in lipid nanodiscs

BsYetJ is a bacterial homolog of transmembrane BAX inhibitor-1 motif-containing 6 (TMBIM6) membrane protein that plays a key role in the control of calcium homeostasis. However, the BsYetJ (or TMBIM6) structure embedded in a lipid bilayer is uncharacterized, let alone the molecular mechanism of the calcium transport activity. Herein, we report structures of BsYetJ in lipid nanodiscs identified by double electron-electron resonance spectroscopy. Our results reveal that BsYetJ in lipid nanodiscs is structurally different from those crystallized in detergents. We show that BsYetJ conformation is pH-sensitive in apo state (lacking calcium), whereas in a calcium-containing solution it is stuck in an intermediate, inert to pH changes. Only when the transmembrane calcium gradient is established can the calcium-release activity of holo-BsYetJ occur and be mediated by pH-dependent conformational changes, suggesting a dual gating mechanism. Conformational substates involved in the process and a key residue D171 relevant to the gating of calcium are identified. Our study suggests that BsYetJ/TMBIM6 is a pH-dependent, voltage-gated calcium channel.

Viewing rare conformations of the β2 adrenergic receptor with pressure-resolved DEER spectroscopy

The β₂ adrenergic receptor (β₂AR) is an archetypal G protein coupled receptor (GPCR). One structural signature of GPCR activation is a large-scale movement (ca. 6 to 14 Å) of transmembrane helix 6 (TM6) to a conformation which binds and activates a cognate G protein. The β₂AR exhibits a low level of agonist-independent G protein activation. The structural origin of this basal activity and its suppression by inverse agonists is unknown but could involve a unique receptor conformation that promotes G protein activation. Alternatively, a conformational selection model proposes that a minor population of the canonical active receptor conformation exists in equilibrium with inactive forms, thus giving rise to basal activity of the ligand-free receptor. Previous spin-labeling and fluorescence resonance energy transfer experiments designed to monitor the positional distribution of TM6 did not detect the presence of the active conformation of ligand-free β₂AR. Here we employ spin-labeling and pressure-resolved double electron-electron resonance spectroscopy to reveal the presence of a minor population of unliganded receptor, with the signature outward TM6 displacement, in equilibrium with inactive conformations. Binding of inverse agonists suppresses this population. These results provide direct structural evidence in favor of a conformational selection model for basal activity in β₂AR and provide a mechanism for inverse agonism. In addition, they emphasize 1) the importance of minor populations in GPCR catalytic function; 2) the use of spin-labeling and variable-pressure electron paramagnetic resonance to reveal them in a membrane protein; and 3) the quantitative evaluation of their thermodynamic properties relative to the inactive forms, including free energy, partial molar volume, and compressibility.

Darunavir-Resistant HIV-1 Protease Constructs Uphold a Conformational Selection Hypothesis for Drug Resistance

Multidrug resistance continues to be a barrier to the effectiveness of highly active antiretroviral therapy in the treatment of human immunodeficiency virus 1 (HIV-1) infection. Darunavir (DRV) is a highly potent protease inhibitor (PI) that is oftentimes effective when drug resistance has emerged against first-generation inhibitors. Resistance to darunavir does evolve and requires 10–20 amino acid substitutions. The conformational landscapes of six highly characterized HIV-1 protease (PR) constructs that harbor up to 19 DRV-associated mutations were characterized by distance measurements with pulsed electron double resonance (PELDOR) paramagnetic resonance spectroscopy, namely double electron–electron resonance (DEER). The results show that the accumulated substitutions alter the conformational landscape compared to PI-naïve protease where the semi-open conformation is destabilized as the dominant population with open-like states becoming prevalent in many cases. A linear correlation is found between values of the DRV inhibition parameter K_i and the open-like to closed-state population ratio determined from DEER. The nearly 50% decrease in occupancy of the semi-open conformation is associated with reduced enzymatic activity, characterized previously in the literature.

Conformational landscape of non-B variants of HIV-1 protease: A pulsed EPR study

HIV infection is a global health epidemic with current FDA-approved HIV-1 Protease inhibitors (PIs) designed against subtype B protease, yet they are used in HIV treatment world-wide regardless of patient HIV classification. In this study, double electron-electron resonance (DEER) electron paramagnetic resonance (EPR) spectroscopy was utilized to gain insights in how natural polymorphisms in several African and Brazilian protease (PR) variants affect the conformational landscape both in the absence and presence of inhibitors. Findings show that Subtypes F and H HIV-1 PR adopt a primarily closed conformation in the unbound state with two secondary mutations, D60E and I62V, postulated to be responsible for the increased probability for closed conformation. In contrast, subtype D, CRF_AG, and CRF_BF HIV-1 PR adopt a primarily semi-open conformation, as observed for PI-naïve-subtype B when unbound by substrate or inhibitor. The impact that inhibitor binding has on shifting the conformational land scape of these variants is also characterized, where analysis provides classification of inhibitor induced shifts away from the semi-open state into weak, moderate and strong effects. The findings are compared to those for prior studies of inhibitor induced conformational shifts in PI-naïve Subtype B, C and CRF_AE.

Experimental Validation of the ALLNOX Program for Studying Protein-Nucleic Acid Complexes

Measurement of distances between spectroscopic labels (e.g., spin labels, fluorophores) attached to specific sites of biomolecules is an important method for studying biomolecular complexes. ALLNOX (Addition of Labels and Linkers) has been developed as a program to model interlabel distances based on an input macromolecule structure. Here, we report validation of ALLNOX using measured distances between nitroxide spin labels attached to specific sites of a protein-DNA complex. The results demonstrate that ALLNOX predicts average interspin distances that matched with values measured with pairs of labels attached at the protein and/or DNA. This establishes a solid foundation for using spin labeling in conjunction with ALLNOX to investigate complexes without high-resolution structures. With its high degree of flexibility for the label or the target biomolecule, ALLNOX provides a useful tool for investigating the structure-function relationship in a large variety of biological molecules.

Physical methods for studying flavoprotein photoreceptors

Molecular mechanisms of dark-to-light state transitions in flavoprotein photoreceptors have been the subject of intense investigation. Blue-light sensing flavoproteins fall into three general classes that share aspects of their activation processes: LOV domains, BLUF proteins, and cryptochromes. In all cases, light-induced changes in flavin redox, protonation, and bonding states result in hydrogen-bond and conformational rearrangements important for regulation of downstream targets. Physical characterization of these flavoprotein states can provide valuable insights into biological function, but clear conclusions are often challenging to draw owing to complexities of data collection and interpretation. In this chapter, we briefly review the three classes of flavoprotein photoreceptors and provide methods for their recombinant production, reconstitution with flavin cofactor, and characterization. We then relate best practices and special considerations for the application of several types of spectroscopies, redox potential measurements, and X-ray scattering experiments to photosensitive flavoproteins. The methods presented are generally accessible to most laboratories.

Optimal Tikhonov regularization for DEER spectroscopy

Tikhonov regularization is the most commonly used method for extracting distance distributions from experimental double electron-electron resonance (DEER) spectroscopy data. This method requires the selection of a regularization parameter, α, and a regularization operator, L. We analyze the performance of a large set of α selection methods and several regularization operators, using a test set of over half a million synthetic noisy DEER traces. These are generated from distance distributions obtained from in silico double labeling of a protein crystal structure of T4 lysozyme with the spin label MTSSL. We compare the methods and operators based on their ability to recover the model distance distributions from the noisy time traces. The results indicate that several α selection methods perform quite well, among them the Akaike information criterion and the generalized cross validation method with either the first- or second-derivative operator. They perform significantly better than currently utilized L-curve methods.

Identification of a ubiquitin-binding interface using Rosetta and DEER

ExoU is a type III-secreted cytotoxin expressing A₂ phospholipase activity when injected into eukaryotic target cells by the bacterium Pseudomonas aeruginosa The enzymatic activity of ExoU is undetectable in vitro unless ubiquitin, a required cofactor, is added to the reaction. The role of ubiquitin in facilitating ExoU enzymatic activity is poorly understood but of significance for designing inhibitors to prevent tissue injury during infections with strains of P. aeruginosa producing this toxin. Most ubiquitin-binding proteins, including ExoU, demonstrate a low (micromolar) affinity for monoubiquitin (monoUb). Additionally, ExoU is a large and dynamic protein, limiting the applicability of traditional structural techniques such as NMR and X-ray crystallography to define this protein-protein interaction. Recent advancements in computational methods, however, have allowed high-resolution protein modeling using sparse data. In this study, we combine double electron-electron resonance (DEER) spectroscopy and Rosetta modeling to identify potential binding interfaces of ExoU and monoUb. The lowest-energy scoring model was tested using biochemical, biophysical, and biological techniques. To verify the binding interface, Rosetta was used to design a panel of mutations to modulate binding, including one variant with enhanced binding affinity. Our analyses show the utility of computational modeling when combined with sensitive biological assays and biophysical approaches that are exquisitely suited for large dynamic proteins.

Effects of Hinge-region Natural Polymorphisms on Human Immunodeficiency Virus-Type 1 Protease Structure, Dynamics, and Drug Pressure Evolution

Multidrug resistance to current Food and Drug Administration-approved HIV-1 protease (PR) inhibitors drives the need to understand the fundamental mechanisms of how drug pressure-selected mutations, which are oftentimes natural polymorphisms, elicit their effect on enzyme function and resistance. Here, the impacts of the hinge-region natural polymorphism at residue 35, glutamate to aspartate (E35D), alone and in conjunction with residue 57, arginine to lysine (R57K), are characterized with the goal of understanding how altered salt bridge interactions between the hinge and flap regions are associated with changes in structure, motional dynamics, conformational sampling, kinetic parameters, and inhibitor affinity. The combined results reveal that the single E35D substitution leads to diminished salt bridge interactions between residues 35 and 57 and gives rise to the stabilization of open-like conformational states with overall increased backbone dynamics. In HIV-1 PR constructs where sites 35 and 57 are both mutated (e.g. E35D and R57K), x-ray structures reveal an altered network of interactions that replace the salt bridge thus stabilizing the structural integrity between the flap and hinge regions. Despite the altered conformational sampling and dynamics when the salt bridge is disrupted, enzyme kinetic parameters and inhibition constants are similar to those obtained for subtype B PR. Results demonstrate that these hinge-region natural polymorphisms, which may arise as drug pressure secondary mutations, alter protein dynamics and the conformational landscape, which are important thermodynamic parameters to consider for development of inhibitors that target for non-subtype B PR.

Conformations of Human Telomeric G-Quadruplex Studied Using a Nucleotide-Independent Nitroxide Label

Guanine-rich oligonucleotides can form a unique G-quadruplex (GQ) structure with stacking units of four guanine bases organized in a plane through Hoogsteen bonding. GQ structures have been detected in vivo and shown to exert their roles in maintaining genome integrity and regulating gene expression. Understanding GQ conformation is important for understanding its inherent biological role and for devising strategies to control and manipulate functions based on targeting GQ. Although a number of biophysical methods have been used to investigate structure and dynamics of GQs, our understanding is far from complete. As such, this work explores the use of the site-directed spin labeling technique, complemented by molecular dynamics simulations, for investigating GQ conformations. A nucleotide-independent nitroxide label (R5), which has been previously applied for probing conformations of noncoding RNA and DNA duplexes, is attached to multiple sites in a 22-nucleotide DNA strand derived from the human telomeric sequence (hTel-22) that is known to form GQ. The R5 labels are shown to minimally impact GQ folding, and inter-R5 distances measured using double electron-electron resonance spectroscopy are shown to adequately distinguish the different topological conformations of hTel-22 and report variations in their occupancies in response to changes of the environment variables such as salt, crowding agent, and small molecule ligand. The work demonstrates that the R5 label is able to probe GQ conformation and establishes the base for using R5 to study more complex sequences, such as those that may potentially form multimeric GQs in long telomeric repeats.