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 1H-13C or 1H-15N correlation spectrum of a large protein, using a 13C- or 15N-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 1H, 19F, 13C, and 15N, 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.

Distinctive phosphoinositide- and Ca2+-binding properties of normal and cognitive performance-linked variant forms of KIBRA C2 domain.

Kidney- and brain-expressed protein (KIBRA), a multifunctional scaffold protein with around 20 known binding partners, is involved in memory and cognition, organ size control via the Hippo pathway, cell polarity, and membrane trafficking. KIBRA includes tandem N-terminal WW domains, a C2 domain, and motifs for binding atypical PKC and PDZ domains. A naturally occurring human KIBRA variant involving residue changes at positions 734 (Met-to-Ile) and 735 (Ser-to-Ala) within the C2 domain affects cognitive performance. We have elucidated 3D structures and calcium- and phosphoinositide-binding properties of human KIBRA C2 domain. Both WT and variant C2 adopt a canonical type I topology C2 domain fold. Neither Ca2+ nor any other metal ion was bound to WT or variant KIBRA C2 in crystal structures, and Ca2+ titration produced no significant reproducible changes in NMR spectra. NMR and X-ray diffraction data indicate that KIBRA C2 binds phosphoinositides via an atypical site involving ß-strands 5, 2, 1, and 8. Molecular dynamics simulations indicate that KIBRA C2 interacts with membranes via primary and secondary sites on the same domain face as the experimentally identified phosphoinositide-binding site. Our results indicate that KIBRA C2 domain association with membranes is calcium-independent and involves distinctive C2 domain-membrane relative orientations.

19F NMR relaxation studies of fluorosubstituted tryptophans.

We present 19F longitudinal and transverse relaxation studies for four differently fluorosubstituted L-tryptophans, which carry single F atoms in the indole ring, both in the context of the free amino acid and when located in the cyclophilin A protein. For the free 4F-, 5F-, 6F-, 7F-L-Trp, satisfactory agreement between experimentally measured and calculated relaxation rates was obtained, suggesting that the parameters used for calculating the rates for the indole frame are sufficiently accurate. We also measured and calculated relaxation rates for four differently 19F-tryptophan labeled cyclophilin A proteins, transferring the parameters from the free amino acid to the protein-bound moiety. Our results suggest that 19F relaxation data of the large and rigid indole ring in Trp are only moderately affected by protein motions and provide critical reference points for evaluating fluorine NMR relaxation in the future, especially in fluorotryptophan labeled proteins.

NMR and MD studies combined to elucidate inhibitor and water interactions of HIV-1 protease and their modulations with resistance mutations.

Over the last two decades, both the sensitivity of NMR and the time scale of molecular dynamics (MD) simulation have increased tremendously and have advanced the field of protein dynamics. HIV-1 protease has been extensively studied using these two methods, and has presented a framework for cross-evaluation of structural ensembles and internal dynamics by integrating the two methods. Here, we review studies from our laboratories over the last several years, to understand the mechanistic basis of protease drug-resistance mutations and inhibitor responses, using NMR and crystal structure-based parallel MD simulations. Our studies demonstrate that NMR relaxation experiments, together with crystal structures and MD simulations, significantly contributed to the current understanding of structural/dynamic changes due to HIV-1 protease drug resistance mutations.

An NMR strategy to detect conformational differences in a protein complexed with highly analogous inhibitors in solution.

This manuscript presents an NMR strategy to investigate conformational differences in protein-inhibitor complexes, when the inhibitors tightly bind to a protein at sub-nanomolar dissociation constants and are highly analogous to each other. Using HIV-1 protease (PR), we previously evaluated amide chemical shift differences, ΔCSPs, of PR bound to darunavir (DRV) compared to PR bound to several DRV analogue inhibitors, to investigate subtle but significant long-distance conformation changes caused by the inhibitor's chemical moiety variation [Khan, S. N., Persons, J. D. Paulsen, J. L., Guerrero, M., Schiffer, C. A., Kurt-Yilmaz, N., and Ishima, R., Biochemistry, (2018), 57, 1652-1662]. However, ΔCSPs are not ideal for investigating subtle PR-inhibitor interface differences because intrinsic differences in the electron shielding of the inhibitors affect protein ΔCSPs. NMR relaxation is also not suitable as it is not sensitive enough to detect small conformational differences in rigid regions among similar PR-inhibitor complexes. Thus, to gain insight into conformational differences at the inhibitor-protein interface, we recorded 15N-half filtered NOESY spectra of PR bound to two highly analogous inhibitors and assessed NOEs between PR amide protons and inhibitor protons, between PR amide protons and hydroxyl side chains, and between PR amide protons and water protons. We also verified the PR amide-water NOEs using 2D water-NOE/ROE experiments. Differences in water-amide proton NOE peaks, possibly due to amide-protein hydrogen bonds, were observed between subunit A and subunit B, and between the DRV-bound form and an analogous inhibitor-bound form, which may contribute to remote conformational changes.

Probing Structural Changes among Analogous Inhibitor-Bound Forms of HIV-1 Protease and a Drug-Resistant Mutant in Solution by Nuclear Magnetic Resonance.

In the era of state-of-the-art inhibitor design and high-resolution structural studies, detection of significant but small protein structural differences in the inhibitor-bound forms is critical to further developing the inhibitor. Here, we probed differences in HIV-1 protease (PR) conformation among darunavir and four analogous inhibitor-bound forms and compared them with a drug-resistant mutant using nuclear magnetic resonance chemical shifts. Changes in amide chemical shifts of wild-type (WT) PR among these inhibitor-bound forms, ΔCSP, were subtle but detectable and extended >10 Å from the inhibitor-binding site, asymmetrically between the two subunits of PR. Molecular dynamics simulations revealed differential local hydrogen bonding as the molecular basis of this remote asymmetric change. Inhibitor-bound forms of the drug-resistant mutant also showed a similar long-range ΔCSP pattern. Differences in ΔCSP values of the WT and the mutant (ΔΔCSPs) were observed at the inhibitor-binding site and in the surrounding region. Comparing chemical shift changes among highly analogous inhibitors and ΔΔCSPs effectively eliminated local environmental effects stemming from different chemical groups and enabled exploitation of these sensitive parameters to detect subtle protein conformational changes and to elucidate asymmetric and remote conformational effects upon inhibitor interaction.

The HIV-1 p66 homodimeric RT exhibits different conformations in the binding-competent and -incompetent NNRTI site.

Non-nucleoside inhibitors of human immunodeficiency virus type 1 reverse transcriptase (RT), NNRTIs, which bind to the p66/p51 heterodimeric RT, also interact with the p66/p66 homodimer, whose structure is unknown. 19 F nuclear magnetic resonance of a single 4-trifluoromethylphenylalanine (tfmF) residue, incorporated into the NNRTI binding pocket of the p66/p66 homodimer at position 181, was used to investigate NNRTI binding. In the NNRTI-bound homodimer complex, two different 19 F signals are observed, with the resonance frequencies matching those of the NNRTI-bound p66/p51 heterodimer spectra, in which the individual p66-subunit or p51-subunit were labeled with tfmF at positions 181. These data suggest that the NNRTI-bound p66/p66 homodimer conformation, particularly around residue 181, is very similar to that in the p66/p51 heterodimer, explaining why NNRTI binding to p66/p66 enhances dimer formation.

A 2-Hydroxyisoquinoline-1,3-Dione Active-Site RNase H Inhibitor Binds in Multiple Modes to HIV-1 Reverse Transcriptase.

The RNase H (RNH) function of HIV-1 reverse transcriptase (RT) plays an essential part in the viral life cycle. We report the characterization of YLC2-155, a 2-hydroxyisoquinoline-1,3-dione (HID)-based active-site RNH inhibitor. YLC2-155 inhibits both polymerase (50% inhibitory concentration [IC50] = 2.6 µM) and RNH functions (IC50 = 0.65 µM) of RT but is more effective against RNH. X-ray crystallography, nuclear magnetic resonance (NMR) analysis, and molecular modeling were used to show that YLC2-155 binds at the RNH-active site in multiple conformations.

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.

Structural integrity of the ribonuclease H domain in HIV-1 reverse transcriptase.

The mature form of reverse transcriptase (RT) is a heterodimer comprising the intact 66-kDa subunit (p66) and a smaller 51-kDa subunit (p51) that is generated by removal of most of the RNase H (RNH) domain from a p66 subunit by proteolytic cleavage between residues 440 and 441. Viral infectivity is eliminated by mutations such as F440A and E438N in the proteolytic cleavage sequence, while normal processing and virus infectivity are restored by a compensatory mutation, T477A, that is located more than 10 Å away from the processing site. The molecular basis for this compensatory effect has remained unclear. We therefore investigated structural characteristics of RNH mutants using computational and experimental approaches. Our Nuclear Magnetic Resonance and Differential Scanning Fluorimetry results show that both F440A and E438N mutations disrupt RNH folding. Addition of the T477A mutation restores correct folding of the RNH domain despite the presence of the F440A or E438N mutations. Molecular dynamics simulations suggest that the T477A mutation affects the processing site by altering relative orientations of secondary structure elements. Predictions of sequence tolerance suggest that phenylalanine and tyrosine are structurally preferred at residues 440 and 441, respectively, which are the P1 and P1' substrate residues known to require bulky side chains for substrate specificity. Interestingly, our study demonstrates that the processing site residues, which are critical for protease substrate specificity and must be exposed to the solvent for efficient processing, also function to maintain proper RNH folding in the p66/p51 heterodimer.

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.

A probe to monitor performance of ¹5N longitudinal relaxation experiments for proteins in solution.

The magnitude of the ¹5N longitudinal relaxation rate typically decreases as magnetic field strength increases in globular proteins in solution. Thus, it is important to test the performance of ¹5N longitudinal relaxation experiments at high field strength. Herein, a tool to investigate systematic errors in ¹5N longitudinal relaxation rate, R1, is introduced. The tool, a difference in R1 values between the two components of the ¹H-coupled ¹5N magnetizations, R1(1)-R1(2), conveniently detects inefficiencies in cancellation of cross correlation between ¹H-¹5N dipolar coupling and ¹5N chemical shift anisotropy. Experiments, in varying conditions, and simulations of a two-spin system indicate that insufficient cancellation of the cross correlation is due to (1) ¹H pulse imperfection and (2) ¹H off-resonance effect, and (3) is further amplified by residual ¹5N transverse magnetization that is caused by the ¹5N off-resonance effect. Results also show that this problem can be easily and practically remedied by discarding the initial decay points when recording ¹5N longitudinal relaxation in proteins.

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.

Assessment of Novel Mesothelin-Specific Human Antibody Domain VH-Fc Fusion Proteins-Based PET Agents.

Mesothelin (MSLN) is a tumor-associated antigen found in a variety of cancers and is a target for imaging and therapeutic applications in MSLN-expressing tumors. We have developed high affinity anti-MSLN human VH domain antibodies, providing alternative targeting vectors to conventional IgG antibodies that are associated with long-circulating half-lives and poor penetration of tumors, limiting antitumor activity in clinical trials. Based on two newly identified anti-MSLN VH binders (3C9, 2A10), we generated VH-Fc fusion proteins and modified them for zirconium-89 radiolabeling to create anti-MSLN VH-Fc PET tracers. The focus of this study was to assess the ability of PET-imaging to compare the in vivo performance of anti-MSLN VH-Fc fusion proteins (2A10, 3C9) targeting different epitopes of MSLN vs IgG1 (m912; a clinical benchmark antibody with an overlapped epitope as 2A10) for PET imaging in a mouse model of colorectal cancer (CRC). The anti-MSLN VH-Fc fusion proteins were successfully modified and radiolabeled with zirconium-89. The resulting MSLN-targeted PET-imaging agents demonstrated specific uptake in the MSLN-expressing HCT116 tumors. The in vivo performance of the MSLN-targeted PET-imaging agents utilizing VH-Fc showed more rapid and greater accumulation and deeper penetration within the tumor than the full-length IgG1 m912-based PET-imaging agent. Furthermore, PET imaging allowed us to compare the pharmacokinetics of epitope-specific VH domain-based PET tracers. Overall, these data are encouraging for the incorporation of PET imaging to assess modified VH domain structures to develop novel anti-MSLN VH domain-based therapeutics in MSLN-positive cancers as well as their companion PET imaging agents.

Preclinical assessment of a novel human antibody VH domain targeting mesothelin as an antibody-drug conjugate.

Mesothelin (MSLN) has been a validated tumor-associated antigen target for several solid tumors for over a decade, making it an attractive option for therapeutic interventions. Novel antibodies with high affinity and better therapeutic properties are needed. In the current study, we have isolated and characterized a novel heavy chain variable (VH) domain 3C9 from a large-size human immunoglobulin VH domain library. 3C9 exhibited high affinity (KD [dissociation constant] <3 nM) and binding specificity in a membrane proteome array (MPA). In a mouse xenograft model, 3C9 fused to human IgG1 Fc was detected at tumor sites as early as 8 h post-infusion and remained at the site for over 10 days. Furthermore, 3C9 fused to a human Fc domain drug conjugate effectively inhibited MSLN-positive tumor growth in a mouse xenograft model. The X-ray crystal structure of full-length MSLN in complex with 3C9 reveals interaction of the 3C9 domains with two distinctive residue patches on the MSLN surface. This newly discovered VH antibody domain has a high potential as a therapeutic candidate for MSLN-expressing cancers.

Quantitative comparison of errors in 15N transverse relaxation rates measured using various CPMG phasing schemes.

Nitrogen-15 Carr-Purcell-Meiboom-Gill (CPMG) transverse relaxation experiment are widely used to characterize protein backbone dynamics and chemical exchange parameters. Although an accurate value of the transverse relaxation rate, R(2), is needed for accurate characterization of dynamics, the uncertainty in the R(2) value depends on the experimental settings and the details of the data analysis itself. Here, we present an analysis of the impact of CPMG pulse phase alternation on the accuracy of the (15)N CPMG R(2). Our simulations show that R(2) can be obtained accurately for a relatively wide spectral width, either using the conventional phase cycle or using phase alternation when the r.f. pulse power is accurately calibrated. However, when the r.f. pulse is miscalibrated, the conventional CPMG experiment exhibits more significant uncertainties in R(2) caused by the off-resonance effect than does the phase alternation experiment. Our experiments show that this effect becomes manifest under the circumstance that the systematic error exceeds that arising from experimental noise. Furthermore, our results provide the means to estimate practical parameter settings that yield accurate values of (15)N transverse relaxation rates in the both CPMG experiments.

Recent developments in (15)N NMR relaxation studies that probe protein backbone dynamics.

Nuclear Magnetic Resonance (NMR) relaxation is a powerful technique that provides information about internal dynamics associated with configurational energetics in proteins, as well as site-specific information involved in conformational equilibria. In particular, (15)N relaxation is a useful probe to characterize overall and internal backbone dynamics of proteins because the relaxation mainly reflects reorientational motion of the N-H bond vector. Over the past 20 years, experiments and protocols for analysis of (15)N R (1), R 2, and the heteronuclear (15)N-{(1)H} NOE data have been well established. The development of these methods has kept pace with the increase in the available static-magnetic field strength, providing dynamic parameters optimized from data fitting at multiple field strengths. Using these methodological advances, correlation times for global tumbling and order parameters and correlation times for internal motions of many proteins have been determined. More recently, transverse relaxation dispersion experiments have extended the range of NMR relaxation studies to the milli- to microsecond time scale, and have provided quantitative information about functional conformational exchange in proteins. Here, we present an overview of recent advances in (15)N relaxation experiments to characterize protein backbone dynamics.

Relative domain orientation of the L289K HIV-1 reverse transcriptase monomer.

HIV-1 reverse transcriptase (RT) is a heterodimer comprised p66 and p51 subunits (p66/p51). Several single amino acid substitutions in RT, including L289K, decrease p66/p51 dimer affinity, and reduce enzymatic functioning. Here, small-angle X-ray scattering (SAXS) with proton paramagnetic relaxation enhancement (PRE), 19 F site-specific NMR, and size exclusion chromatography (SEC) were performed for the p66 monomer with the L289K mutation, p66L289K . NMR and SAXS experiments clearly elucidated that the thumb and RNH domains in the monomer do not rigidly interact with each other but are spatially close to the RNH domain. Based on this structural model of the monomer, p66L289K and p51 were predicted to form a heterodimer while p66 and p51L289K not. We tested this hypothesis by SEC analysis of p66 and p51 containing L289K in different combinations and clearly demonstrated that L289K substitution in the p51 subunit, but not in the p66 subunit, reduces p66/p51 formation. Based on the derived monomer model and the importance of the inter-subunit RNH-thumb domain interaction in p66/p51, validated by SEC, the mechanism of p66 homodimer formation was discussed.

A highly-specific fully-human antibody and CAR-T cells targeting CD66e/CEACAM5 are cytotoxic for CD66e-expressing cancer cells in vitro and in vivo.

Neuro-endocrine prostate cancer (NEPC) accounts for about 20% of lethal metastatic castration-resistant prostate cancer (CRPC). NEPC has the most aggressive biologic behavior of all prostate cancers and is associated with poor patient outcome. Effective treatment for NEPC is not available because NEPC exhibit distinct cell-surface expression profiles compared to other types of prostate cancer. Recently, the carcinoembryonic antigen-related cell adhesion molecule 5 (CEACAM5) (known as CEA or CD66e) was suggested to be a specific surface protein marker for NEPC. Therefore, we identified a new, fully-human anti-CEACAM5 monoclonal antibody, 1G9, which bound to the most proximal membrane domains, A3 and B3, of CEACAM5 with high affinity and specificity. It shows no off-target binding to other CEACAM family members, membrane distal domains of CEACAM5, or 5800 human membrane proteins. IgG1 1G9 exhibited CEACAM5-specific ADCC activity toward CEACAM5-positive prostate cancer cells in vitro and in vivo. Chimeric antigen receptor T cells (CAR-T) based on scFv 1G9 induced specific and strong antitumor activity in a mouse model of prostate cancer. Our results suggest that IgG1 and CAR-T cells based on 1G9 are promising candidate therapeutics for CEACAM5-positive NEPC and other cancers.

Inhibitors of HIV-1 Nef-Mediated Activation of the Myeloid Src-Family Kinase Hck Block HIV-1 Replication in Macrophages and Disrupt MHC-I Downregulation.

HIV-1 Nef is an attractive target for antiretroviral drug discovery because of its role in promoting HIV-1 infectivity, replication, and host immune system avoidance. Here, we applied a screening strategy in which recombinant HIV-1 Nef protein was coupled to activation of the Src-family tyrosine kinase Hck, which enhances the HIV-1 life cycle in macrophages. Nef stimulates recombinant Hck activity in vitro, providing a robust assay for chemical library screening. High-throughput screening of more than 730 000 compounds using the Nef·Hck assay identified six unique hit compounds that bound directly to recombinant Nef by surface plasmon resonance (SPR) in vitro and inhibited HIV-1 replication in primary macrophages in the 0.04 to 5 µM range without cytotoxicity. Eighty-four analogs were synthesized around an isothiazolone scaffold from this series, many of which bound to recombinant Nef and inhibited HIV-1 infectivity in the low to submicromolar range. Compounds in this series restored MHC-I to the surface of HIV-infected primary cells and disrupted a recombinant protein complex of Nef with the C-terminal tail of MHC-I and the µ1 subunit of the AP-1 endocytic trafficking protein. Nef inhibitors in this class have the potential to block HIV-1 replication in myeloid cells and trigger recognition of HIV-infected cells by the adaptive immune system in vivo.