Solution characterization of [methyl-(13)C]methionine HIV-1 reverse transcriptase by NMR spectroscopy

Biophysical Characterization of p51 and p66 Monomers of HIV-1 Reverse Transcriptase with Their Inhibitors

Human immunodeficiency virus (HIV)-1 reverse transcriptase (HIV-1 RT) is responsible for the transcription of viral RNA genomes into DNA genomes and has become an important target for the treatment of acquired immune deficiency syndrome (AIDS). This study used biophysical techniques to characterize the HIV-1 RT structure, monomer forms, and the non-nucleoside reverse transcriptase inhibitors (NNRTIs) bound forms. Inactive p66W401A and p51W401A were selected as models to study the HIV-1 RT monomer structures. Nuclear magnetic resonance (NMR) spectroscopy revealed that the unliganded forms of p66W401A protein and p51W401A protein had similar conformation to each other in solution. The complexes of p66W401A or p51W401A with inhibitors showed similar conformations to p66 in the RT heterodimer bound to the NNRTIs. Furthermore, the results of paramagnetic relaxation enhancement (PRE)-assisted NMR revealed that the unliganded forms of the p66W401A and p51W401A conformations were different from the unliganded heterodimer, characterized by a greater distance between the fingers and thumb subdomains. Small-angle X-ray scattering (SAXS) experiments confirmed that p66W401A and p51W401A can bind with inhibitors, similar to the p66/p51 heterodimer. The findings of this study increase the structural knowledge base of HIV-1 RT monomers, which may be helpful in the future design of potent viral inhibitors.

Targeting the Structural Maturation Pathway of HIV-1 Reverse Transcriptase

Formation of active HIV-1 reverse transcriptase (RT) proceeds via a structural maturation process that involves subdomain rearrangements and formation of an asymmetric p66/p66' homodimer. These studies were undertaken to evaluate whether the information about this maturation process can be used to identify small molecule ligands that retard or interfere with the steps involved. We utilized the isolated polymerase domain, p51, rather than p66, since the initial subdomain rearrangements are largely limited to this domain. Target sites at subdomain interfaces were identified and computational analysis used to obtain an initial set of ligands for screening. Chromatographic evaluations of the p51 homodimer/monomer ratio support the feasibility of this approach. Ligands that bind near the interfaces and a ligand that binds directly to a region of the fingers subdomain involved in subunit interface formation were identified, and the interactions were further characterized by NMR spectroscopy and X-ray crystallography. Although these ligands were found to reduce dimer formation, further efforts will be required to obtain ligands with higher binding affinity. In contrast with previous ligand identification studies performed on the RT heterodimer, subunit interface surfaces are solvent-accessible in the p51 and p66 monomers, making these constructs preferable for identification of ligands that directly interfere with dimerization.

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.

Methyl TROSY spectroscopy: A versatile NMR approach to study challenging biological systems

A major goal in structural biology is to unravel how molecular machines function in detail. To that end, solution-state NMR spectroscopy is ideally suited as it is able to study biological assemblies in a near natural environment. Based on methyl TROSY methods, it is now possible to record high-quality data on complexes that are far over 100 kDa in molecular weight. In this review, we discuss the theoretical background of methyl TROSY spectroscopy, the information that can be extracted from methyl TROSY spectra and approaches that can be used to assign methyl resonances in large complexes. In addition, we touch upon insights that have been obtained for a number of challenging biological systems, including the 20S proteasome, the RNA exosome, molecular chaperones and G-protein-coupled receptors. We anticipate that methyl TROSY methods will be increasingly important in modern structural biology approaches, where information regarding static structures is complemented with insights into conformational changes and dynamic intermolecular interactions.

Mass Spectrometric Characterization of HIV-1 Reverse Transcriptase Interactions with Non-nucleoside Reverse Transcriptase Inhibitors

Non-nucleoside reverse transcriptase inhibitors (NNRTIs) of human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) have been developed for the treatment of acquired immunodeficiency syndrome. HIV-1 RT binding to NNRTIs has been characterized by various biophysical techniques. However, these techniques are often hampered by the low water solubility of the inhibitors, such as the current promising diarylpyrimidine-based inhibitors rilpivirine and etravirine. Hence, a conventional and rapid method that requires small sample amounts is desirable for studying NNRTIs with low water solubility. Here we successfully applied a recently developed mass spectrometric technique under non-denaturing conditions to characterize the interactions between the heterodimeric HIV-1 RT enzyme and NNRTIs with different inhibitory activities. Our data demonstrate that mass spectrometry serves as a semi-quantitative indicator of NNRTI binding affinity for HIV-1 RT using low and small amounts of samples, offering a new high-throughput screening tool for identifying novel RT inhibitors as anti-HIV drugs.

Conformational Plasticity of the NNRTI-Binding Pocket in HIV-1 Reverse Transcriptase: A Fluorine Nuclear Magnetic Resonance Study

HIV-1 reverse transcriptase (RT) is a major drug target in the treatment of HIV-1 infection. RT inhibitors currently in use include non-nucleoside, allosteric RT inhibitors (NNRTIs), which bind to a hydrophobic pocket, distinct from the enzyme's active site. We investigated RT-NNRTI interactions by solution (19)F nuclear magnetic resonance (NMR), using singly (19)F-labeled RT proteins. Comparison of (19)F chemical shifts of fluorinated RT and drug-resistant variants revealed that the fluorine resonance is a sensitive probe for identifying mutation-induced changes in the enzyme. Our data show that in the unliganded enzyme, the NNRTI-binding pocket is highly plastic and not locked into a single conformation. Upon inhibitor binding, the binding pocket becomes rigidified. In the inhibitor-bound state, the (19)F signal of RT is similar to that of drug-resistant mutant enzymes, distinct from what is observed for the free state. Our results demonstrate the power of (19)F NMR spectroscopy to characterize conformational properties using selectively (19)F-labeled protein.

NMR characterization of HIV-1 reverse transcriptase binding to various non-nucleoside reverse transcriptase inhibitors with different activities

Human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) is an important target for antiviral therapy against acquired immunodeficiency syndrome. However, the efficiency of available drugs is impaired most typically by drug-resistance mutations in this enzyme. In this study, we applied a nuclear magnetic resonance (NMR) spectroscopic technique to the characterization of the binding of HIV-1 RT to various non-nucleoside reverse transcriptase inhibitors (NNRTIs) with different activities, i.e., nevirapine, delavirdine, efavirenz, dapivirine, etravirine, and rilpivirine. ¹H-¹³C heteronuclear single-quantum coherence (HSQC) spectral data of HIV-1 RT, in which the methionine methyl groups of the p66 subunit were selectively labeled with ¹³C, were collected in the presence and absence of these NNRTIs. We found that the methyl ¹³C chemical shifts of the M230 resonance of HIV-1 RT bound to these drugs exhibited a high correlation with their anti-HIV-1 RT activities. This methionine residue is located in proximity to the NNRTI-binding pocket but not directly involved in drug interactions and serves as a conformational probe, indicating that the open conformation of HIV-1 RT was more populated with NNRTIs with higher inhibitory activities. Thus, the NMR approach offers a useful tool to screen for novel NNRTIs in developing anti-HIV drugs.

The p66 immature precursor of HIV-1 reverse transcriptase

In contrast to the wealth of structural data available for the mature p66/p51 heterodimeric human immunodeficiency virus type 1 reverse transcriptase (RT), the structure of the homodimeric p66 precursor remains unknown. In all X-ray structures of mature RT, free or complexed, the processing site in the p66 subunit, for generating the p51 subunit, is sequestered into a β-strand within the folded ribonuclease H (RNH) domain and is not readily accessible to proteolysis, rendering it difficult to propose a simple and straightforward mechanism of the maturation step. Here, we investigated, by solution NMR, the conformation of the RT p66 homodimer. Our data demonstrate that the RNH and Thumb domains in the p66 homodimer are folded and possess conformations very similar to those in mature RT. This finding suggests that maturation models which invoke a complete or predominantly unfolded RNH domain are unlikely. The present study lays the foundation for further in-depth mechanistic investigations at the atomic level.

Protein-mediated antagonism between HIV reverse transcriptase ligands nevirapine and MgATP

Nonnucleoside reverse transcriptase inhibitors (NNRTIs) play a central role in the treatment of AIDS, but their mechanisms of action are incompletely understood. The interaction of the NNRTI nevirapine (NVP) with HIV-1 reverse transcriptase (RT) is characterized by a preference for the open conformation of the fingers/thumb subdomains, and a reported variation of three orders of magnitude between the binding affinity of NVP for RT in the presence or absence of primer/template DNA. To investigate the relationship between conformation and ligand binding, we evaluated the use of methionine NMR probes positioned near the tip of the fingers or thumb subdomains. Such probes would be expected to be sensitive to changes in the local environment depending on the fractions of open and closed RT. Comparisons of the NMR spectra of three conservative mutations, I63M, L74M, and L289M, indicated that M63 showed the greatest shift sensitivity to the addition of NVP. The exchange kinetics of the M63 resonance are fast on the chemical shift timescale, but become slow in the presence of NVP due to the slow binding of RT with the inhibitor. The simplest model consistent with this behavior involves a rapid open/closed equilibrium coupled with a slow interaction of the inhibitor with the open conformation. Studies of RT in the presence of both NVP and MgATP indicate a strong negative cooperativity. Binding of MgATP reduces the fraction of RT bound to NVP, as indicated by the intensity of the NVP-perturbed M230 resonance, and enhances the dissociation rate constant of the NVP, resulting in an increase of the open/closed interconversion rate, so that the M63 resonance moves into the fast/intermediate-exchange regime. Protein-mediated interactions appear to explain most of the affinity variation of NVP for RT.

Biophysical Insights into the Inhibitory Mechanism of Non-Nucleoside HIV-1 Reverse Transcriptase Inhibitors

HIV-1 reverse transcriptase (RT) plays a central role in HIV infection. Current United States Federal Drug Administration (USFDA)-approved antiretroviral therapies can include one of five approved non-nucleoside RT inhibitors (NNRTIs), which are potent inhibitors of RT activity. Despite their crucial clinical role in treating and preventing HIV-1 infection, their mechanism of action remains elusive. In this review, we introduce RT and highlight major advances from experimental and computational biophysical experiments toward an understanding of RT function and the inhibitory mechanism(s) of NNRTIs.

Detecting allosteric sites of HIV-1 reverse transcriptase by X-ray crystallographic fragment screening

HIV-1 reverse transcriptase (RT) undergoes a series of conformational changes during viral replication and is a central target for antiretroviral therapy. The intrinsic flexibility of RT can provide novel allosteric sites for inhibition. Crystals of RT that diffract X-rays to better than 2 Å resolution facilitated the probing of RT for new druggable sites using fragment screening by X-ray crystallography. A total of 775 fragments were grouped into 143 cocktails, which were soaked into crystals of RT in complex with the non-nucleoside drug rilpivirine (TMC278). Seven new sites were discovered, including the Incoming Nucleotide Binding, Knuckles, NNRTI Adjacent, and 399 sites, located in the polymerase region of RT, and the 428, RNase H Primer Grip Adjacent, and 507 sites, located in the RNase H region. Three of these sites (Knuckles, NNRTI Adjacent, and Incoming Nucleotide Binding) are inhibitory and provide opportunities for discovery of new anti-AIDS drugs.

Metal and ligand binding to the HIV-RNase H active site are remotely monitored by Ile556

HIV-1 reverse transcriptase (RT) contains a C-terminal ribonuclease H (RH) domain on its p66 subunit that can be expressed as a stable, although inactive protein. Recent studies of several RH enzymes demonstrate that substrate binding plays a major role in the creation of the active site. In the absence of substrate, the C-terminal helix E of the RT RNase H domain is dynamic, characterized by severe exchange broadening of its backbone amide resonances, so that the solution characterization of this region of the protein has been limited. Nuclear magnetic resonance studies of 13C-labeled RH as a function of experimental conditions reveal that the δ1 methyl resonance of Ile556, located in a short, random coil segment following helix E, experiences a large 13C shift corresponding to a conformational change of Ile556 that results from packing of helix E against the central β-sheet. This shift provides a useful basis for monitoring the effects of various ligands on active site formation. Additionally, we report that the RNase H complexes formed with one or both divalent ions can be individually observed and characterized using diamagnetic Zn2+ as a substitute for Mg2+. Ordering of helix E results specifically from the interaction with the lower affinity binding to the A divalent ion site.

Fast methionine-based solution structure determination of calcium-calmodulin complexes

Here we present a novel NMR method for the structure determination of calcium-calmodulin (Ca(2+)-CaM)-peptide complexes from a limited set of experimental restraints. A comparison of solved CaM-peptide structures reveals invariability in CaM's backbone conformation and a structural plasticity in CaM's domain orientation enabled by a flexible linker. Knowing this, the collection and analysis of an extensive set of NOESY spectra is redundant. Although RDCs can define CaM domain orientation in the complex, they lack the translational information required to position the domains on the bound peptide and highlight the necessity of intermolecular NOEs. Here we employ a specific isotope labeling strategy in which the role of methionine in CaM-peptide interactions is exploited to collect these critical NOEs. By (1)H, (13)C-labeling the methyl groups of deuterated methionine against a (2)H, (12)C background, we can acquire a (13)C-edited NOESY characterized by simplified, easily analyzable spectra. Together with measured CaM backbone H(N)-N RDCs and intrapeptide NOE-based distances, these intermolecular NOEs provide restraints for a low temperature torsion-angle dynamics and simulated annealing protocol used to calculate the complex structure. We have applied our method to a CaM complex previously solved through X-ray crystallography: Ca(2+)-CaM bound to the CaM kinase I peptide (PDB code: 1MXE). The resulting structure has a backbone RMSD of 1.6 Å to that previously published. We have also used this test complex to investigate the importance of homologous model selection on the calculated outcome. In addition to having application for fast complex structure determination, this method can be used to determine the structures of difficult complexes characterized by chemical shift overlap and broad signals for which the traditional method based on the use of fully (13)C, (15)N-labeled CaM fails.

Efavirenz binding site in HIV-1 reverse transcriptase monomers

Efavirenz (EFV) is a potent nonnucleoside reverse transcriptase inhibitor (NNRTI) used in the treatment of AIDS. NNRTIs bind in a hydrophobic pocket located in the p66 subunit of reverse transcriptase (RT), which is not present in crystal structures of RT without an inhibitor. Recent studies showed that monomeric forms of the p66 and p51 subunits bind efavirenz with micromolar affinity. The effect of efavirenz on the solution conformations of p66 and p51 monomers was studied by hydrogen-deuterium exchange mass spectrometry (HXMS) and Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS). HXMS data reveal that five peptides, four of which contain efavirenz contact residues seen in the crystal structure of the RT-EFV complex, exhibit a reduced level of exchange in monomer-EFV complexes. Moreover, peptide 232-246 undergoes slow cooperative unfolding-refolding in the bound monomers, but at a rate much slower than that observed in the p66 subunit of the RT heterodimer [Seckler, J. M., Howard, K. J., Barkley, M. D., and Wintrode, P. L. (2009) Biochemistry 48, 7646-7655]. These results suggest that the efavirenz binding site on p66 and p51 monomers is similar to the NNRTI binding pocket in the p66 subunit of RT. Nanoelectrospray ionization FT-ICR mass spectra indicate that the intact monomers each have (at least) two different conformations. In the presence of efavirenz, the mass spectra change significantly and suggest that p51 adopts a single, more compact conformation, whereas p66 undergoes facile, electrospray-induced cleavage. The population shift is consistent with a selected-fit binding mechanism.

Conformational dependence of 13C shielding and coupling constants for methionine methyl groups

Methionine residues fulfill a broad range of roles in protein function related to conformational plasticity, ligand binding, and sensing/mediating the effects of oxidative stress. A high degree of internal mobility, intrinsic detection sensitivity of the methyl group, and low copy number have made methionine labeling a popular approach for NMR investigation of selectively labeled protein macromolecules. However, selective labeling approaches are subject to more limited information content. In order to optimize the information available from such studies, we have performed DFT calculations on model systems to evaluate the conformational dependence of (3)J (CSCC), (3)J (CSCH), and the isotropic shielding, sigma(iso). Results have been compared with experimental data reported in the literature, as well as data obtained on [methyl-(13)C]methionine and on model compounds. These studies indicate that relative to oxygen, the presence of the sulfur atom in the coupling pathway results in a significantly smaller coupling constant, (3)J (CSCC)/(3)J (COCC) approximately 0.7. It is further demonstrated that the (3)J (CSCH) coupling constant depends primarily on the subtended CSCH dihedral angle, and secondarily on the CSCC dihedral angle. Comparison of theoretical shielding calculations with the experimental shift range of the methyl group for methionine residues in proteins supports the conclusion that the intra-residue conformationally-dependent shift perturbation is the dominant determinant of delta(13)Cepsilon. Analysis of calmodulin data based on these calculations indicates that several residues adopt non-standard rotamers characterized by very large approximately 100 degrees chi(3) values. The utility of the delta(13)Cepsilon as a basis for estimating the gauche/trans ratio for chi(3) is evaluated, and physical and technical factors that limit the accuracy of both the NMR and crystallographic analyses are discussed.

Homodimerization of the p51 subunit of HIV-1 reverse transcriptase

The dimerization of HIV reverse transcriptase (RT), required to obtain the active form of the enzyme, is influenced by mutations, non-nucleoside reverse transcriptase inhibitors (NNRTIs), nucleotide substrates, Mg ions, temperature, and specifically designed dimerization inhibitors. In this study, we have utilized nuclear magnetic resonance (NMR) spectroscopy of the [methyl-(13)C]methionine-labeled enzyme and small-angle X-ray scattering (SAXS) to investigate how several of these factors influence the dimerization behavior of the p51 subunit. The (1)H-(13)C HSQC spectrum of p51 obtained at micromolar concentrations indicates that a significant fraction of the p51 adopts a "p66-like" conformation. SAXS data obtained for p51 samples were used to determine the fractions of monomer and dimer in the sample and to evaluate the conformation of the fingers/thumb subdomain. All of the p51 monomer observed was found to adopt the compact, "p51C" conformation observed for the p51 subunit in the RT heterodimer. The NMR and SAXS data indicate that the p51 homodimer adopts a structure that is similar to the p66/p51 heterodimer, with one p51C subunit and a second p51 subunit in an extended, "p51E" conformation that resembles the p66 subunit of the heterodimer. The fractional dimer concentration and the fingers/thumb orientation are found to depend strongly on the experimental conditions and exhibit a qualitative dependence on nevirapine and ionic strength (KCl) that is similar to the behavior reported for the heterodimer and the p66 homodimer. The L289K mutation interferes with p51 homodimer formation as it does with formation of the heterodimer, despite its location far from the dimer interface. This effect is readily interpreted in terms of a conformational selection model, in which p51(L289K) has a much greater preference for the compact, p51C conformation. A reduced level of dimer formation then results from the reduced ratio of the p51E(L289K) to p51C(L289K) monomers.