Solute effects on spin labels at an aqueous-exposed site in the flap region of HIV-1 protease

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.

Effect of Type I Antifreeze Proteins on the Freezing and Melting Processes of Cryoprotective Solutions Studied by Site-Directed Spin Labeling Technique

Antifreeze proteins (AFPs) protect organisms living in subzero environments from freezing injury, which render them potential applications for cryopreservation of living cells, organs, and tissues. Cryoprotective agents (CPAs), such as glycerol and propylene glycol, have been used as ingredients to treat cellular tissues and organs to prevent ice crystal’s formation at low temperatures. To assess AFP’s function in CPA solutions, we have the applied site-directed spin labeling technique to a Type I AFP. A two-step process to prevent bulk freezing of the CPA solutions was observed by the cryo-photo microscopy, i.e., (1) thermodynamic freezing point depression by the CPAs; and (2) inhibition to the growth of seed ice crystals by the AFP. Electron paramagnetic resonance (EPR) experiments were also carried out from room temperature to 97 K, and vice versa. The EPR results indicate that the spin labeled AFP bound to ice surfaces, and inhibit the growths of ice through the bulk freezing processes in the CPA solutions. The ice-surface bound AFP in the frozen matrices could also prevent the formation of large ice crystals during the melting processes of the solutions. Our study illustrates that AFPs can play an active role in CPA solutions for cryopreservation applications.

Mechanisms of antifreeze proteins investigated via the site-directed spin labeling technique

The site-directed spin labeling (SDSL) technique was used to examine the antifreeze mechanisms of type-I antifreeze proteins (AFPs). The effects on the growth of seed ice crystals by the spin-label groups attached to different side chains of the AFPs were observed, and the states of water molecules surrounding the spin-label groups were probed via analyses of variable-temperature (VT) dependent electron paramagnetic resonance (EPR) spectra. The first set of experiments revealed the antifreeze activities of the spin-labeled AFPs at the microscopic level, while the second set of experiments displayed those at the molecular level. The experimental results confirmed the putative ice-binding surface (IBS) of type-I AFPs. The VT EPR spectra indicate that type-I AFPs can inhibit the nucleation of seed ice crystals down to ~ - 20 °C in their aqueous solutions. Thus, the present authors believe that AFPs protect organisms from freezing damage in two ways: (1) inhibiting the nucleation of seed ice crystals, and (2) hindering the growth of seed ice crystals once they have formed. The first mechanism should play a more significant role in protecting against freezing damage among organisms living in cold environments. The VT EPR spectra also revealed that liquid-like water molecules existed around the spin-labeled non-ice-binding side chains of the AFPs frozen within the ice matrices, and ice surrounding the spin-label groups melted at subzero temperatures during the heating process. This manuscript concludes with the proposed model of antifreeze mechanisms of AFPs based on the experimental results.

Effector Roles of Putidaredoxin on Cytochrome P450cam Conformational States

In this study, the effector role of Pdx (putidaredoxin) on cytochrome P450cam conformation is refined by attaching two different spin labels, MTSL or BSL (bifunctional spin-label) onto the F or G helices and using DEER (double electron-electron resonance) to measure the distance between labels. Recent EPR and crystallographic studies have observed that oxidized Pdx induces substrate-bound P450cam to change from the closed to the open state. However, this change was not observed by DEER in the reduced Pdx complex with carbon-monoxide-bound P450cam (Fe(2+)CO). In addition, recent NMR studies have failed to observe a change in P450cam conformation upon binding Pdx. Hence, resolving these issues is important for a full understanding the effector role of Pdx. Here we show that oxidized Pdx induces camphor-bound P450cam to shift from the closed to the open conformation when labeled on either the F or G helices with MTSL. BSL at these sites can either narrow the distance distribution widths dramatically or alter the extent of the conformational change. In addition, we report DEER spectra on a mixed oxidation state containing oxidized Pdx and ferrous CO-bound P450cam, showing that P450cam remains closed. This indicates that CO binding to the heme prevents P450cam from opening, overriding the influence exerted by Pdx binding. Finally, we report the open form P450cam crystal structure with substrate bound, which suggests that crystal packing effects may prevent conformational conversion. Using multiple labeling approaches, DEER provides a unique perspective to resolve how the conformation of P450cam depends on Pdx and ligand states.

Pulsed EPR characterization of HIV-1 protease conformational sampling and inhibitor-induced population shifts

The conformational landscape of HIV-1 protease (PR) can be experimentally characterized by pulsed-EPR double electron-electron resonance (DEER). For this characterization, nitroxide spin labels are attached to an engineered cysteine residue in the flap region of HIV-1 PR. DEER distance measurements from spin-labels contained within each flap of the homodimer provide a detailed description of the conformational sampling of apo-enzyme as well as induced conformational shifts as a function of inhibitor binding. The distance distribution profiles are further interpreted in terms of a conformational ensemble scheme that consists of four unique states termed "curled/tucked", "closed", "semi-open" and "wide-open" conformations. Reported here are the DEER results for a drug-resistant variant clinical isolate sequence, V6, in the presence of FDA approved protease inhibitors (PIs) as well as a non-hydrolyzable substrate mimic, CaP2. Results are interpreted in the context of the current understanding of the relationship between conformational sampling, drug resistance, and kinetic efficiency of HIV-1PR as derived from previous DEER and kinetic data for a series of HIV-1PR constructs that contain drug-pressure selected mutations or natural polymorphisms. Specifically, these collective results support the notion that inhibitor-induced closure of the flaps correlates with inhibitor efficiency and drug resistance. This body of work also suggests DEER as a tool for studying conformational sampling in flexible enzymes as it relates to function.

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

An understanding of macromolecular conformational equilibrium in biological systems is oftentimes essential to understand function, dysfunction, and disease. For the past few years, our lab has been utilizing site-directed spin labeling (SDSL), coupled with electron paramagnetic resonance (EPR) spectroscopy, to characterize the conformational ensemble and ligand-induced conformational shifts of HIV-1 protease (HIV-1PR). The biomedical importance of characterizing the fractional occupancy of states within the conformational ensemble critically impacts our hypothesis of a conformational selection mechanism of drug-resistance evolution in HIV-1PR. The purpose of the following chapter is to give a timeline perspective of our SDSL EPR approach to characterizing conformational sampling of HIV-1PR. We provide detailed instructions for the procedure utilized in analyzing distance profiles for HIV-1PR obtained from pulsed electron-electron double resonance (PELDOR). Specifically, we employ a version of PELDOR known as double electron-electron resonance (DEER). Data are processed with the software package "DeerAnalysis" (http://www.epr.ethz.ch/software), which implements Tikhonov regularization (TKR), to generate a distance profile from electron spin-echo amplitude modulations. We assign meaning to resultant distance profiles based upon a conformational sampling model, which is described herein. The TKR distance profiles are reconstructed with a linear combination of Gaussian functions, which is then statistically analyzed. In general, DEER has proven powerful for observing structural ensembles in proteins and, more recently, nucleic acids. Our goal is to present our advances in order to aid readers in similar applications.

MESMER: minimal ensemble solutions to multiple experimental restraints

MOTIVATION

Macromolecular structures and interactions are intrinsically heterogeneous, temporally adopting a range of configurations that can confound the analysis of data from bulk experiments. To obtain quantitative insights into heterogeneous systems, an ensemble-based approach can be employed, in which predicted data computed from a collection of models is compared to the observed experimental results. By simultaneously fitting orthogonal structural data (e.g. small-angle X-ray scattering, nuclear magnetic resonance residual dipolar couplings, dipolar electron-electron resonance spectra), the range and population of accessible macromolecule structures can be probed.

RESULTS

We have developed MESMER, software that enables the user to identify ensembles that can recapitulate experimental data by refining thousands of component collections selected from an input pool of potential structures. The MESMER suite includes a powerful graphical user interface (GUI) to streamline usage of the command-line tools, calculate data from structure libraries and perform analyses of conformational and structural heterogeneity. To allow for incorporation of other data types, modular Python plugins enable users to compute and fit data from nearly any type of quantitative experimental data.

RESULTS

Conformational heterogeneity in three macromolecular systems was analyzed with MESMER, demonstrating the utility of the streamlined, user-friendly software.

AVAILABILITY AND IMPLEMENTATION

https://code.google.com/p/mesmer/

An Integrated Spin-Labeling/Computational-Modeling Approach for Mapping Global Structures of Nucleic Acids

The technique of site-directed spin labeling (SDSL) provides unique information on biomolecules by monitoring the behavior of a stable radical tag (i.e., spin label) using electron paramagnetic resonance (EPR) spectroscopy. In this chapter, we describe an approach in which SDSL is integrated with computational modeling to map conformations of nucleic acids. This approach builds upon a SDSL tool kit previously developed and validated, which includes three components: (i) a nucleotide-independent nitroxide probe, designated as R5, which can be efficiently attached at defined sites within arbitrary nucleic acid sequences; (ii) inter-R5 distances in the nanometer range, measured via pulsed EPR; and (iii) an efficient program, called NASNOX, that computes inter-R5 distances on given nucleic acid structures. Following a general framework of data mining, our approach uses multiple sets of measured inter-R5 distances to retrieve "correct" all-atom models from a large ensemble of models. The pool of models can be generated independently without relying on the inter-R5 distances, thus allowing a large degree of flexibility in integrating the SDSL-measured distances with a modeling approach best suited for the specific system under investigation. As such, the integrative experimental/computational approach described here represents a hybrid method for determining all-atom models based on experimentally-derived distance measurements.

Experimental mapping of DNA duplex shape enabled by global lineshape analyses of a nucleotide-independent nitroxide probe

Sequence-dependent variation in structure and dynamics of a DNA duplex, collectively referred to as ‘DNA shape’, critically impacts interactions between DNA and proteins. Here, a method based on the technique of site-directed spin labeling was developed to experimentally map shapes of two DNA duplexes that contain response elements of the p53 tumor suppressor. An R5a nitroxide spin label, which was covalently attached at a specific phosphate group, was scanned consecutively through the DNA duplex. X-band continuous-wave electron paramagnetic resonance spectroscopy was used to monitor rotational motions of R5a, which report on DNA structure and dynamics at the labeling site. An approach based on Pearson's coefficient analysis was developed to collectively examine the degree of similarity among the ensemble of R5a spectra. The resulting Pearson's coefficients were used to generate maps representing variation of R5a mobility along the DNA duplex. The R5a mobility maps were found to correlate with maps of certain DNA helical parameters, and were capable of revealing similarity and deviation in the shape of the two closely related DNA duplexes. Collectively, the R5a probe and the Pearson's coefficient-based lineshape analysis scheme yielded a generalizable method for examining sequence-dependent DNA shapes.

Effects of PRE and POST therapy drug-pressure selected mutations on HIV-1 protease conformational sampling

Conformational sampling of pre- and post-therapy subtype B HIV-1 protease sequences derived from a pediatric subject infected via maternal transmission with HIV-1 were characterized by double electron-electron resonance spectroscopy. The conformational ensemble of the PRE construct resembles native-like inhibitor bound states. In contrast, the POST construct, which contains accumulated drug-pressure selected mutations, has a predominantly semi-open conformational ensemble, with increased populations of open-like states. The single point mutant L63P, which is contained in PRE and POST, has decreased dynamics, particularly in the flap region, and also displays a closed-like conformation of inhibitor-bound states. These findings support our hypothesis that secondary mutations accumulate in HIV-1 protease to shift conformational sampling to stabilize open-like conformations, while maintaining the predominant semi-open conformation for activity.

The role of select subtype polymorphisms on HIV-1 protease conformational sampling and dynamics

HIV-1 protease is an essential enzyme for viral particle maturation and is a target in the fight against HIV-1 infection worldwide. Several natural polymorphisms are also associated with drug resistance. Here, we utilized both pulsed electron double resonance, also called double electron-electron resonance, and NMR (15)N relaxation measurements to characterize equilibrium conformational sampling and backbone dynamics of an HIV-1 protease construct containing four specific natural polymorphisms commonly found in subtypes A, F, and CRF_01 A/E. Results show enhanced backbone dynamics, particularly in the flap region, and the persistence of a novel conformational ensemble that we hypothesize is an alternative flap orientation of a curled open state or an asymmetric configuration when interacting with inhibitors.

Pulsed EPR distance measurements in soluble proteins by site-directed spin labeling (SDSL)

The resurgence of pulsed electron paramagnetic resonance (EPR) in structural biology centers on recent improvements in distance measurements using the double electron-electron resonance (DEER) technique. This unit focuses on EPR-based distance measurements by site-directed spin labeling (SDSL) of engineered cysteine residues in soluble proteins, with HIV-1 protease used as a model. To elucidate conformational changes in proteins, experimental protocols were optimized and existing data analysis programs were employed to derive distance-distribution profiles. Experimental considerations, sample preparation, and error analysis for artifact suppression are also outlined herein.

Conformational dynamics and distribution of nitroxide spin labels

Long-range distance measurements based on paramagnetic relaxation enhancement (PRE) in NMR, quantification of surface water dynamics near biomacromolecules by Overhauser dynamic nuclear polarization (DNP) and sensitivity enhancement by solid-state DNP all depend on introducing paramagnetic species into an otherwise diamagnetic NMR sample. The species can be introduced by site-directed spin labeling, which offers precise control for positioning the label in the sequence of a biopolymer. However, internal flexibility of the spin label gives rise to dynamic processes that potentially influence PRE and DNP behavior and leads to a spatial distribution of the electron spin even in solid samples. Internal dynamics of spin labels and their static conformational distributions have been studied mainly by electron paramagnetic resonance spectroscopy and molecular dynamics simulations, with a large body of results for the most widely applied methanethiosulfonate spin label MTSL. These results are critically discussed in a unifying picture based on rotameric states of the group that carries the spin label. Deficiencies in our current understanding of dynamics and conformations of spin labeled groups and of their influence on NMR observables are highlighted and directions for further research suggested.

Elucidating a relationship between conformational sampling and drug resistance in HIV-1 protease

Enzyme targets in rapidly replicating systems, such as retroviruses, commonly respond to drug-selective pressure with mutations arising in the active site pocket that limit inhibitor effectiveness by introducing steric hindrance or by eliminating essential molecular interactions. However, these primary mutations are disposed to compromising pathogenic fitness. Emerging secondary mutations, which are often found outside of the binding cavity, may or can restore fitness while maintaining drug resistance. The accumulated drug pressure selected mutations could have an indirect effect in the development of resistance, such as altering protein flexibility or the dynamics of protein-ligand interactions. Here, we show that accumulation of mutations in a drug-resistant HIV-1 protease (HIV-1 PR) variant, D30N/M36I/A71V, changes the fractional occupancy of the equilibrium conformational sampling ensemble. Correlations are made among populations of the conformational states, namely, closed-like, semiopen, and open-like, with inhibition constants, as well as kinetic parameters. Mutations that stabilize a closed-like conformation correlate with enzymes of lowered activity and with higher affinity for inhibitors, which is corroborated by a further increase in the fractional occupancy of the closed state upon addition of inhibitor or substrate-mimic. Cross-resistance is found to correlate with combinations of mutations that increase the population of the open-like conformations at the expense of the closed-like state while retaining native-like occupancy of the semiopen population. These correlations suggest that at least three states are required in the conformational sampling model to establish the emergence of drug resistance in HIV-1 PR. More importantly, these results shed light on a possible mechanism whereby mutations combine to impart drug resistance while maintaining catalytic activity.

Nonlinear scaling of surface water diffusion with bulk water viscosity of crowded solutions

The translational hydration dynamics within 0.5-1.5 nm of the surface of a DPPC liposome, a model biomacromolecular surface, is analyzed by the recently developed Overhauser dynamic nuclear polarization (ODNP) technique. We find that dramatic changes to the bulk solvent cause only weak changes in the surface hydration dynamics. Specifically, both a >10-fold increase in bulk viscosity and the restriction of diffusion by confinement on a multiple nm length-scale change the local translational diffusion coefficient of the surface water surrounding the lipid bilayer by <2.5-fold. By contrast, previous ODNP studies have shown that changes to the biomacromolecular surface induced by folding, binding, or aggregation can cause local hydration dynamics to vary by factors of up to 30. We suggest that the surface topology and chemistry at the ≤1.5 nm scale, rather than the characteristics of the solvent, nearly exclusively determine the macromolecule's surface hydration dynamics.

Simulation and Modeling of Crowding Effects on the Thermodynamic and Kinetic Properties of Proteins with Atomic Details

Recent experimental studies of protein folding and binding under crowded solutions suggest that crowding agents exert subtle influences on the thermodynamic and kinetic properties of the proteins. While some of the crowding effects can be understood qualitatively from simple models of the proteins, quantitative rationalization of these effects requires an atomistic representation of the protein molecules in modeling their interactions with crowders. A computational approach, known as postprocessing, has opened the door for atomistic modeling of crowding effects. This review summarizes the applications of the postprocessing approach for studying crowding effects on the thermodynamics and kinetics of protein folding, conformational transition, and binding. The integration of atomistic modeling with experiments in crowded solutions promises new insight into biochemical processes in cellular environments.

Inhibitor-induced conformational shifts and ligand-exchange dynamics for HIV-1 protease measured by pulsed EPR and NMR spectroscopy

Double electron-electron resonance (DEER) spectroscopy was utilized to investigate shifts in conformational sampling induced by nine FDA-approved protease inhibitors (PIs) and a nonhydrolyzable substrate mimic for human immunodeficiency virus type 1 protease (HIV-1 PR) subtype B, subtype C, and CRF_01 A/E. The ligand-bound subtype C protease has broader DEER distance profiles, but trends for inhibitor-induced conformational shifts are comparable to those previously reported for subtype B. Ritonavir, one of the strong-binding inhibitors for subtypes B and C, induces less of the closed conformation in CRF_01 A/E. (1)H-(15)N heteronuclear single-quantum coherence (HSQC) spectra were acquired for each protease construct titrated with the same set of inhibitors. NMR (1)H-(15)N HSQC titration data show that inhibitor residence time in the protein binding pocket, inferred from resonance exchange broadening, shifting or splitting correlates with the degree of ligand-induced flap closure measured by DEER spectroscopy. These parallel results show that the ligand-induced conformational shifts resulting from protein-ligand interactions characterized by DEER spectroscopy of HIV-1 PR obtained at the cryogenic temperature are consistent with more physiological solution protein-ligand interactions observed by solution NMR spectroscopy.

Conformational selection underlies recognition of a molybdoenzyme by its dedicated chaperone

Molecular recognition is central to all biological processes. Understanding the key role played by dedicated chaperones in metalloprotein folding and assembly requires the knowledge of their conformational ensembles. In this study, the NarJ chaperone dedicated to the assembly of the membrane-bound respiratory nitrate reductase complex NarGHI, a molybdenum-iron containing metalloprotein, was taken as a model of dedicated chaperone. The combination of two techniques ie site-directed spin labeling followed by EPR spectroscopy and ion mobility mass spectrometry, was used to get information about the structure and conformational dynamics of the NarJ chaperone upon binding the N-terminus of the NarG metalloprotein partner. By the study of singly spin-labeled proteins, the E119 residue present in a conserved elongated hydrophobic groove of NarJ was shown to be part of the interaction site. Moreover, doubly spin-labeled proteins studied by pulsed double electron-electron resonance (DEER) spectroscopy revealed a large and composite distribution of inter-label distances that evolves into a single preexisting one upon complex formation. Additionally, ion mobility mass spectrometry experiments fully support these findings by revealing the existence of several conformers in equilibrium through the distinction of different drift time curves and the selection of one of them upon complex formation. Taken together our work provides a detailed view of the structural flexibility of a dedicated chaperone and suggests that the exquisite recognition and binding of the N-terminus of the metalloprotein is governed by a conformational selection mechanism.

Differential Flap Dynamics in Wild-type and a Drug Resistant Variant of HIV-1 Protease Revealed by Molecular Dynamics and NMR Relaxation

In the rapidly evolving disease of HIV drug resistance readily emerges, nullifying the effectiveness of therapy. Drug resistance has been extensively studied in HIV-1 protease where resistance occurs when the balance between enzyme inhibition and substrate recognition and turn-over is perturbed to favor catalytic activity. Mutations which confer drug resistance can impact the dynamics and structure of both the bound and unbound forms of the enzyme. Flap+ is a multi-drug-resistant variant of HIV-1 protease with a combination of mutations at the edge of the active site, within the active site, and in the flaps (L10I, G48V, I54V, V82A). The impact of these mutations on the dynamics in the unliganded form in comparison with the wild-type protease was elucidated with Molecular Dynamic simulations and NMR relaxation experiments. The comparative analyses from both methods concur in showing that the enzyme's dynamics are impacted by the drug resistance mutations in Flap+ protease. These alterations in the enzyme dynamics, particularly within the flaps, likely modulate the balance between substrate turn-over and drug binding, thereby conferring drug resistance.

DEER distance measurements on proteins

Distance distributions between paramagnetic centers in the range of 1.8 to 6 nm in membrane proteins and up to 10 nm in deuterated soluble proteins can be measured by the DEER technique. The number of paramagnetic centers and their relative orientation can be characterized. DEER does not require crystallization and is not limited with respect to the size of the protein or protein complex. Diamagnetic proteins are accessible by site-directed spin labeling. To characterize structure or structural changes, experimental protocols were optimized and techniques for artifact suppression were introduced. Data analysis programs were developed, and it was realized that interpretation of the distance distributions must take into account the conformational distribution of spin labels. First methods have appeared for deriving structural models from a small number of distance constraints. The present scope and limitations of the technique are illustrated.

Characterization of the disordered-to-α-helical transition of IA₃ by SDSL-EPR spectroscopy

Electron paramagnetic resonance (EPR) spectroscopy coupled with site-directed spin labeling (SDSL) is a valuable tool for characterizing the mobility and conformational changes of proteins but has seldom been applied to intrinsically disordered proteins (IDPs). Here, IA₃ is used as a model system demonstrating SDSL-EPR characterization of conformational changes in small IDP systems. IA₃ has 68 amino acids, is unstructured in solution, and becomes α-helical upon addition of the secondary structural stabilizer 2,2,2-trifluoroethanol (TFE). Two single cysteine substitutions, one in the N-terminus (S14C) and one in the C-terminus (N58C), were generated and labeled with three different nitroxide spin labels. The resultant EPR line shapes of each of the labels were compared and each reported changes in mobility upon addition of TFE. Specifically, the spectral line shape parameters h((+1))/h(₀), the local tumbling volume (V(L)), and the percent change of the h(₋₁) intensity were utilized to quantitatively monitor TFE-induced conformational changes. The values of h((+1)/)h(₀) as a function of TFE titration varied in a sigmoidal manner and were fit to a two-state Boltzmann model that provided values for the midpoint of the transition, thus, reporting on the global conformational change of IA₃. The other parameters provide site-specific information and show that S14C-SL undergoes a conformational change resulting in more restricted motion than N58C-SL, which is consistent with previously published results obtained by studies using NMR and circular dichroism spectroscopy indicating a higher degree of α-helical propensity of the N-terminal segment of IA₃. Overall, the results provide a framework for data analyzes that can be used to study induced unstructured-to-helical conformations in IDPs by SDSL.

Conformational distributions at the N-peptide/boxB RNA interface studied using site-directed spin labeling

In bacteriophage λ, interactions between a 22-amino acid peptide (called the N-peptide) and a stem-loop RNA element (called boxB) play a critical role in transcription anti-termination. The N-peptide/boxB complex has been extensively studied, and serves as a paradigm for understanding mechanisms of protein/RNA recognition. Particularly, ultrafast spectroscopy techniques have been applied to monitor picosecond fluorescence decay behaviors of 2-aminopurines embedded at various positions of the boxB RNA. The studies have led to a model in which the bound N-peptide exists in dynamic equilibrium between two states, with peptide C-terminal fragment either stacking on (i.e., the stacked state) or peeling away from (i.e., the unstacked state) the RNA loop. The function of the N-peptide/boxB complex seems to correlate with the fraction of the stacked state. Here, the N-peptide/boxB system is studied using the site-directed spin labeling technique, in which X-band electron paramagnetic resonance spectroscopy is applied to monitor nanosecond rotational behaviors of stable nitroxide radicals covalently attached to different positions of the N-peptide. The data reveal that in the nanosecond regime the C-terminal fragment of bound N-peptide adopts multiple discrete conformations within the complex. The characteristics of these conformations are consistent with the proposed stacked and unstacked states, and their distributions vary upon mutations within the N-peptide. These results suggest that the dynamic two-state model remains valid in the nanosecond regime, and represents a unique mode of function in the N-peptide/boxB RNA complex. It also demonstrates a connection between picosecond and nanosecond dynamics in a biological complex.

Method to Predict Crowding Effects by Postprocessing Molecular Dynamics Trajectories: Application to the Flap Dynamics of HIV-1 Protease

The internal dynamics of proteins inside of cells may be affected by the crowded intracellular environments. Here, we test a novel approach to simulations of crowding, in which simulations in the absence of crowders are postprocessed to predict crowding effects, against the direct approach of simulations in the presence of crowders. The effects of crowding on the flap dynamics of HIV-1 protease predicted by the postprocessing approach are found to agree well with those calculated by the direct approach. The postprocessing approach presents distinct advantages over the direct approach in terms of accuracy and speed and is expected to have broad impact on atomistic simulations of macromolecular crowding.

Monitoring inhibitor-induced conformational population shifts in HIV-1 protease by pulsed EPR spectroscopy

Double electron-electron resonance (DEER), a pulsed electron paramagnetic resonance (EPR) spectroscopy technique, was utilized to characterize conformational population shifts in HIV-1 protease (HIV-1PR) upon interaction with various inhibitors. Distances between spin-labeled sites in the flap region of HIV-1PR were determined, and detailed analyses provide population percentages for the ensemble flap conformations upon interaction with inhibitor or substrate. Comparisons are made between the percentage of the closed conformer seen with DEER and enzymatic inhibition constants, thermodynamic dissociation constants, and the number of hydrogen bonds identified in crystallographic complexes.