Kinetic coevolutionary models predict the temporal emergence of HIV-1 resistance mutations under drug selection pressure

Drug resistance in HIV type 1 (HIV-1) is a pervasive problem that affects the lives of millions of people worldwide. Although records of drug-resistant mutations (DRMs) have been extensively tabulated within public repositories, our understanding of the evolutionary kinetics of DRMs and how they evolve together remains limited. Epistasis, the interaction between a DRM and other residues in HIV-1 protein sequences, is key to the temporal evolution of drug resistance. We use a Potts sequence-covariation statistical-energy model of HIV-1 protein fitness under drug selection pressure, which captures epistatic interactions between all positions, combined with kinetic Monte-Carlo simulations of sequence evolutionary trajectories, to explore the acquisition of DRMs as they arise in an ensemble of drug-naive patient protein sequences. We follow the time course of 52 DRMs in the enzymes protease, RT, and integrase, the primary targets of antiretroviral therapy. The rates at which DRMs emerge are highly correlated with their observed acquisition rates reported in the literature when drug pressure is applied. This result highlights the central role of epistasis in determining the kinetics governing DRM emergence. Whereas rapidly acquired DRMs begin to accumulate as soon as drug pressure is applied, slowly acquired DRMs are contingent on accessory mutations that appear only after prolonged drug pressure. We provide a foundation for using computational methods to determine the temporal evolution of drug resistance using Potts statistical potentials, which can be used to gain mechanistic insights into drug resistance pathways in HIV-1 and other infectious agents.

Design of a SARS-CoV-2 papain-like protease inhibitor with antiviral efficacy in a mouse model

The emergence of SARS-CoV-2 variants and drug-resistant mutants calls for additional oral antivirals. The SARS-CoV-2 papain-like protease (PLpro) is a promising but challenging drug target. We designed and synthesized 85 noncovalent PLpro inhibitors that bind to a recently discovered ubiquitin binding site and the known BL2 groove pocket near the S4 subsite. Leads inhibited PLpro with the inhibitory constant Ki values from 13.2 to 88.2 nanomolar. The co-crystal structures of PLpro with eight leads revealed their interaction modes. The in vivo lead Jun12682 inhibited SARS-CoV-2 and its variants, including nirmatrelvir-resistant strains with EC50 from 0.44 to 2.02 micromolar. Oral treatment with Jun12682 improved survival and reduced lung viral loads and lesions in a SARS-CoV-2 infection mouse model, suggesting that PLpro inhibitors are promising oral SARS-CoV-2 antiviral candidates.

Structures, functions and adaptations of the human LINE-1 ORF2 protein

The LINE-1 (L1) retrotransposon is an ancient genetic parasite that has written around one-third of the human genome through a 'copy and paste' mechanism catalysed by its multifunctional enzyme, open reading frame 2 protein (ORF2p)1. ORF2p reverse transcriptase (RT) and endonuclease activities have been implicated in the pathophysiology of cancer2,3, autoimmunity4,5 and ageing6,7, making ORF2p a potential therapeutic target. However, a lack of structural and mechanistic knowledge has hampered efforts to rationally exploit it. We report structures of the human ORF2p 'core' (residues 238-1061, including the RT domain) by X-ray crystallography and cryo-electron microscopy in several conformational states. Our analyses identified two previously undescribed folded domains, extensive contacts to RNA templates and associated adaptations that contribute to unique aspects of the L1 replication cycle. Computed integrative structural models of full-length ORF2p show a dynamic closed-ring conformation that appears to open during retrotransposition. We characterize ORF2p RT inhibition and reveal its underlying structural basis. Imaging and biochemistry show that non-canonical cytosolic ORF2p RT activity can produce RNA:DNA hybrids, activating innate immune signalling through cGAS/STING and resulting in interferon production6-8. In contrast to retroviral RTs, L1 RT is efficiently primed by short RNAs and hairpins, which probably explains cytosolic priming. Other biochemical activities including processivity, DNA-directed polymerization, non-templated base addition and template switching together allow us to propose a revised L1 insertion model. Finally, our evolutionary analysis demonstrates structural conservation between ORF2p and other RNA- and DNA-dependent polymerases. We therefore provide key mechanistic insights into L1 polymerization and insertion, shed light on the evolutionary history of L1 and enable rational drug development targeting L1.

Design of SARS-CoV-2 papain-like protease inhibitor with antiviral efficacy in a mouse model

The emergence of SARS-CoV-2 variants and drug-resistant mutants calls for additional oral antivirals. The SARS-CoV-2 papain-like protease (PLpro) is a promising but challenging drug target. In this study, we designed and synthesized 85 noncovalent PLpro inhibitors that bind to the newly discovered Val70Ub site and the known BL2 groove pocket. Potent compounds inhibited PLpro with inhibitory constant Ki values from 13.2 to 88.2 nM. The co-crystal structures of PLpro with eight leads revealed their interaction modes. The in vivo lead Jun12682 inhibited SARS-CoV-2 and its variants, including nirmatrelvir-resistant strains with EC50 from 0.44 to 2.02 μM. Oral treatment with Jun12682 significantly improved survival and reduced lung viral loads and lesions in a SARS-CoV-2 infection mouse model, suggesting PLpro inhibitors are promising oral SARS-CoV-2 antiviral candidates.

Discovery of diarylpyrimidine derivatives bearing piperazine sulfonyl as potent HIV-1 nonnucleoside reverse transcriptase inhibitors

HIV-1 reverse transcriptase is one of the most attractive targets for the treatment of AIDS. However, the rapid emergence of drug-resistant strains and unsatisfactory drug-like properties seriously limit the clinical application of HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs). Here we show that a series of piperazine sulfonyl-bearing diarylpyrimidine-based NNRTIs were designed to improve the potency against wild-type and NNRTI-resistant strains by enhancing backbone-binding interactions. Among them, compound 18b1 demonstrates single-digit nanomolar potency against the wild-type and five mutant HIV-1 strains, which is significantly better than the approved drug etravirine. The co-crystal structure analysis and molecular dynamics simulation studies were conducted to explain the broad-spectrum inhibitory activity of 18b1 against reverse transcriptase variants. Besides, compound 18b1 demonstrates improved water solubility, cytochrome P450 liability, and other pharmacokinetic properties compared to the currently approved diarylpyrimidine (DAPY) NNRTIs. Therefore, we consider compound 18b1 a potential lead compound worthy of further study.

Halo Library, a Tool for Rapid Identification of Ligand Binding Sites on Proteins Using Crystallographic Fragment Screening

X-ray crystallographic fragment screening (XCFS) uses fragment-sized molecules (∼60 to 300 Da) to access binding sites on proteins that may be inaccessible to larger drug-like molecules (>300 Da). Previous studies have shown that fragments containing halogen atoms bind more often to proteins than non-halogenated fragments. Here, we designed the Halo Library containing 46 halogenated fragments (including the "universal fragment" 4-bromopyrazole), a majority of which have been reported to bind to or inhibit one or more targets. The library was screened against the crystals of HIV-1 reverse transcriptase with the drug rilpivirine, yielding an overall hit rate of 26%. Two new binding sites were discovered, and several hot spots were identified. This small library may thus provide a convenient tool for rapidly assessing the feasibility of a target for XCFS, mapping hot spots and cryptic sites, as well as finding fragment binders that can be useful for developing drug leads.

Targeting HIV-1 Reverse Transcriptase Using a Fragment-Based Approach

Human immunodeficiency virus type I (HIV-1) is a retrovirus that infects cells of the host's immune system leading to acquired immunodeficiency syndrome and potentially death. Although treatments are available to prevent its progression, HIV-1 remains a major burden on health resources worldwide. Continued emergence of drug-resistance mutations drives the need for novel drugs that can inhibit HIV-1 replication through new pathways. The viral protein reverse transcriptase (RT) plays a fundamental role in the HIV-1 replication cycle, and multiple approved medications target this enzyme. In this study, fragment-based drug discovery was used to optimize a previously identified hit fragment (compound B-1), which bound RT at a novel site. Three series of compounds were synthesized and evaluated for their HIV-1 RT binding and inhibition. These series were designed to investigate different vectors around the initial hit in an attempt to improve inhibitory activity against RT. Our results show that the 4-position of the core scaffold is important for binding of the fragment to RT, and a lead compound with a cyclopropyl substitution was selected and further investigated. Requirements for binding to the NNRTI-binding pocket (NNIBP) and a novel adjacent site were investigated, with lead compound 27-a minimal but efficient NNRTI-offering a starting site for the development of novel dual NNIBP-Adjacent site inhibitors.

Biochemical and structural insights into SARS-CoV-2 polyprotein processing by Mpro

SARS-CoV-2, a human coronavirus, is the causative agent of the COVID-19 pandemic. Its genome is translated into two large polyproteins subsequently cleaved by viral papain-like protease and main protease (Mpro). Polyprotein processing is essential yet incompletely understood. We studied Mpro-mediated processing of the nsp7-11 polyprotein, whose mature products include cofactors of the viral replicase, and identified the order of cleavages. Integrative modeling based on mass spectrometry (including hydrogen-deuterium exchange and cross-linking) and x-ray scattering yielded a nsp7-11 structural ensemble, demonstrating shared secondary structural elements with individual nsps. The pattern of cross-links and HDX footprint of the C145A Mpro and nsp7-11 complex demonstrate preferential binding of the enzyme active site to the polyprotein junction sites and additional transient contacts to help orient the enzyme on its substrate for cleavage. Last, proteolysis assays were used to characterize the effect of inhibitors/binders on Mpro processing/inhibition using the nsp7-11 polyprotein as substrate.

Cryo-EM structure of the HIV-1 Pol polyprotein provides insights into virion maturation

Key proteins of retroviruses and other RNA viruses are translated and subsequently processed from polyprotein precursors by the viral protease (PR). Processing of the HIV Gag-Pol polyprotein yields the HIV structural proteins and enzymes. Structures of the mature enzymes PR, reverse transcriptase (RT), and integrase (IN) aided understanding of catalysis and design of antiretrovirals, but knowledge of the Pol precursor architecture and function before PR cleavage is limited. We developed a system to produce stable HIV-1 Pol and determined its cryo-electron microscopy structure. RT in Pol has a similar arrangement to the mature RT heterodimer, and its dimerization may draw together two PR monomers to activate proteolytic processing. HIV-1 thus may leverage the dimerization interfaces in Pol to regulate assembly and maturation of polyprotein precursors.

Evolution of the SARS-CoV-2 proteome in three dimensions (3D) during the first 6 months of the COVID-19 pandemic

Understanding the molecular evolution of the SARS-CoV-2 virus as it continues to spread in communities around the globe is important for mitigation and future pandemic preparedness. Three-dimensional structures of SARS-CoV-2 proteins and those of other coronavirusess archived in the Protein Data Bank were used to analyze viral proteome evolution during the first 6 months of the COVID-19 pandemic. Analyses of spatial locations, chemical properties, and structural and energetic impacts of the observed amino acid changes in >48 000 viral isolates revealed how each one of 29 viral proteins have undergone amino acid changes. Catalytic residues in active sites and binding residues in protein-protein interfaces showed modest, but significant, numbers of substitutions, highlighting the mutational robustness of the viral proteome. Energetics calculations showed that the impact of substitutions on the thermodynamic stability of the proteome follows a universal bi-Gaussian distribution. Detailed results are presented for potential drug discovery targets and the four structural proteins that comprise the virion, highlighting substitutions with the potential to impact protein structure, enzyme activity, and protein-protein and protein-nucleic acid interfaces. Characterizing the evolution of the virus in three dimensions provides testable insights into viral protein function and should aid in structure-based drug discovery efforts as well as the prospective identification of amino acid substitutions with potential for drug resistance.

Structural basis of HIV inhibition by L-nucleosides: opportunities for drug development and repurposing

Infection with HIV can cripple the immune system and lead to AIDS. Hepatitis B virus (HBV) is a hepadnavirus that causes human liver diseases. Both pathogens are major public health problems affecting millions of people worldwide. The polymerases from both viruses are the most common drug target for viral inhibition, sharing common architecture at their active sites. The L-nucleoside drugs emtricitabine and lamivudine are widely used HIV reverse transcriptase (RT) and HBV polymerase (Pol) inhibitors. Nevertheless, structural details of their binding to RT(Pol)/nucleic acid remained unknown until recently. Here, we discuss the implications of these structures, alongside related complexes with L-dNTPs, for the development of novel L-nucleos(t)ide drugs, and prospects for repurposing them.

HIV-1 gp120 Antagonists Also Inhibit HIV-1 Reverse Transcriptase by Bridging the NNRTI and NRTI Sites

HIV-1 infection is typically treated using ≥2 drugs, including at least one HIV-1 reverse transcriptase (RT) inhibitor. Drugs targeting RT comprise nucleos(t)ide RT inhibitors (NRTIs) and non-nucleoside RT inhibitors (NNRTIs). NRTI-triphosphates bind at the polymerase active site and, following incorporation, inhibit DNA elongation. NNRTIs bind at an allosteric pocket ∼10 Å away from the polymerase active site. This study focuses on compounds ("NBD derivatives") originally developed to bind to HIV-1 gp120, some of which inhibit RT. We have determined crystal structures of three NBD compounds in complex with HIV-1 RT, correlating with RT enzyme inhibition and antiviral activity, to develop structure-activity relationships. Intriguingly, these compounds bridge the dNTP and NNRTI-binding sites and inhibit the polymerase activity of RT in the enzymatic assays (IC50 < 5 μM). Two of the lead compounds, NBD-14189 and NBD-14270, show potent antiviral activity (EC50 < 200 nM), and NBD-14270 shows low cytotoxicity (CC50 > 100 μM).

Crystal Structure of a Retroviral Polyprotein: Prototype Foamy Virus Protease-Reverse Transcriptase (PR-RT)

In most cases, proteolytic processing of the retroviral Pol portion of the Gag-Pol polyprotein precursor produces protease (PR), reverse transcriptase (RT), and integrase (IN). However, foamy viruses (FVs) express Pol separately from Gag and, when Pol is processed, only the IN domain is released. Here, we report a 2.9 Å resolution crystal structure of the mature PR-RT from prototype FV (PFV) that can carry out both proteolytic processing and reverse transcription but is in a configuration not competent for proteolytic or polymerase activity. PFV PR-RT is monomeric and the architecture of PFV PR is similar to one of the subunits of HIV-1 PR, which is a dimer. There is a C-terminal extension of PFV PR (101-145) that consists of two helices which are adjacent to the base of the RT palm subdomain, and anchors PR to RT. The polymerase domain of PFV RT consists of fingers, palm, thumb, and connection subdomains whose spatial arrangements are similar to the p51 subunit of HIV-1 RT. The RNase H and polymerase domains of PFV RT are connected by flexible linkers. Significant spatial and conformational (sub)domain rearrangements are therefore required for nucleic acid binding. The structure of PFV PR-RT provides insights into the conformational maturation of retroviral Pol polyproteins.

Revealing the Structural Plasticity of SARS-CoV-2 nsp7 and nsp8 Using Structural Proteomics

Coronavirus (CoV) nonstructural proteins (nsps) assemble to form the replication-transcription complex (RTC) responsible for viral RNA synthesis. nsp7 and nsp8 are important cofactors of the RTC, as they interact and regulate the activity of RNA-dependent RNA polymerase and other nsps. To date, no structure of the full-length SARS-CoV-2 nsp7:nsp8 complex has been published. The current understanding of this complex is based on structures from truncated constructs, with missing electron densities, or from related CoV species where SARS-CoV-2 nsp7 and nsp8 share upward of 90% sequence identity. Despite available structures solved using crystallography and cryo-EM representing detailed static snapshots of the nsp7:nsp8 complex, it is evident that the complex has a high degree of structural plasticity. However, relatively little is known about the conformational dynamics of the individual proteins and how they complex to interact with other nsps. Here, the solution-based structural proteomic techniques, hydrogen-deuterium exchange mass spectrometry (HDX-MS) and cross-linking mass spectrometry (XL-MS), illuminate the dynamics of SARS-CoV-2 full-length nsp7 and nsp8 proteins and the nsp7:nsp8 protein complex. Results presented from the two techniques are complementary and validate the interaction surfaces identified from the published three-dimensional heterotetrameric crystal structure of the SARS-CoV-2 truncated nsp7:nsp8 complex. Furthermore, mapping of XL-MS data onto higher-order complexes suggests that SARS-CoV-2 nsp7 and nsp8 do not assemble into a hexadecameric structure as implied by the SARS-CoV full-length nsp7:nsp8 crystal structure. Instead, our results suggest that the nsp7:nsp8 heterotetramer can dissociate into a stable dimeric unit that might bind to nsp12 in the RTC without significantly altering nsp7-nsp8 interactions.

Co-crystallization with diabodies: A case study for the introduction of synthetic symmetry

This work presents a method for introducing synthetic symmetry into protein crystallization samples using an antibody fragment termed a diabody (Dab). These Dabs contain two target binding sites, and engineered disulfide bonds have been included to modulate Dab flexibility. The impacts of Dab engineering have been observed through assessment of thermal stability, small-angle X-ray scattering, and high-resolution crystal structures. Complexes between the engineered Dabs and HIV-1 reverse transcriptase (RT) bound to a high-affinity DNA aptamer were also generated to explore the capacity of engineered Dabs to enable the crystallization of bound target proteins. This strategy increased the crystallization hit frequency obtained for RT-aptamer, and the structure of a Dab-RT-aptamer complex was determined to 3.0-Å resolution. Introduction of synthetic symmetry using a Dab could be a broadly applicable strategy, especially when monoclonal antibodies for a target have previously been identified.

2,4,5-Trisubstituted Pyrimidines as Potent HIV-1 NNRTIs: Rational Design, Synthesis, Activity Evaluation, and Crystallographic Studies

There is an urgent unmet medical need for novel human immunodeficiency virus type 1 (HIV-1) inhibitors that are effective against a variety of NNRTI-resistance mutations. We report our research efforts aimed at discovering a novel chemotype of anti-HIV-1 agents with improved potency against a variety of NNRTI-resistance mutations in this paper. Structural modifications of the lead K-5a2 led to the identification of a potent inhibitor 16c. 16c yielded highly potent anti-HIV-1 activities and improved resistance profiles compared with the approved drug etravirine. The co-crystal structure revealed the key role of the water networks surrounding the NNIBP for binding and for resilience against resistance mutations, while suggesting further extension of 16c toward the NNRTI-adjacent site as a lead development strategy. Furthermore, 16c demonstrated favorable pharmacokinetic and safety properties, suggesting the potential of 16c as a promising anti-HIV-1 drug candidate.

Resolving the Dynamic Motions of SARS-CoV-2 nsp7 and nsp8 Proteins Using Structural Proteomics

Coronavirus (CoV) non-structural proteins (nsps) assemble to form the replication-transcription complex (RTC) responsible for viral RNA synthesis. nsp7 and nsp8 are important cofactors of the RTC, as they interact and regulate the activity of RNA-dependent RNA polymerase (RdRp) and other nsps. To date, no structure of full-length SARS-CoV-2 nsp7:nsp8 complex has been published. Current understanding of this complex is based on structures from truncated constructs or with missing electron densities and complexes from related CoV species with which SARS-CoV-2 nsp7 and nsp8 share upwards of 90% sequence identity. Despite available structures being solved using crystallography and cryo-EM representing detailed snapshots of the nsp7:nsp8 complex, it is evident that the complex has a high degree of structural plasticity. However, relatively little is known about the conformational dynamics of the complex and how it assembles to interact with other nsps. Here, the solution-based structural proteomic techniques, hydrogen-deuterium exchange mass spectrometry (HDX-MS) and crosslinking mass spectrometry (XL-MS), illuminate the structural dynamics of the SARS-CoV-2 full-length nsp7:nsp8 complex. The results presented from the two techniques are complementary and validate the interaction surfaces identified from the published three-dimensional heterotetrameric crystal structure of SARS-CoV-2 truncated nsp7:nsp8 complex. Furthermore, mapping of XL-MS data onto higher order complexes suggests that SARS-CoV-2 nsp7 and nsp8 do not assemble into a hexadecameric structure as implied by the SARS-CoV full-length nsp7:nsp8 crystal structure. Instead our results suggest that the nsp7:nsp8 heterotetramer can dissociate into a stable dimeric unit that might bind to nsp12 in the RTC without altering nsp7-nsp8 interactions.

Evolution of the SARS-CoV-2 proteome in three dimensions (3D) during the first six months of the COVID-19 pandemic

Three-dimensional structures of SARS-CoV-2 and other coronaviral proteins archived in the Protein Data Bank were used to analyze viral proteome evolution during the first six months of the COVID-19 pandemic. Analyses of spatial locations, chemical properties, and structural and energetic impacts of the observed amino acid changes in >48,000 viral proteome sequences showed how each one of the 29 viral study proteins have undergone amino acid changes. Structural models computed for every unique sequence variant revealed that most substitutions map to protein surfaces and boundary layers with a minority affecting hydrophobic cores. Conservative changes were observed more frequently in cores versus boundary layers/surfaces. Active sites and protein-protein interfaces showed modest numbers of substitutions. Energetics calculations showed that the impact of substitutions on the thermodynamic stability of the proteome follows a universal bi-Gaussian distribution. Detailed results are presented for six drug discovery targets and four structural proteins comprising the virion, highlighting substitutions with the potential to impact protein structure, enzyme activity, and functional interfaces. Characterizing the evolution of the virus in three dimensions provides testable insights into viral protein function and should aid in structure-based drug discovery efforts as well as the prospective identification of amino acid substitutions with potential for drug resistance.

Integrative structural biology studies of HIV-1 reverse transcriptase binding to a high-affinity DNA aptamer

The high-resolution crystal structure of HIV-1 reverse transcriptase (RT) bound to a 38-mer DNA hairpin aptamer with low pM affinity was previously described. The high-affinity binding aptamer contained 2'-O-methyl modifications and a seven base-pair GC-rich tract and the structure of the RT-aptamer complex revealed specific contacts between RT and the template strand of the aptamer. Similar to all crystal structures of RT bound to nucleic acid template-primers, the aptamer bound RT with a bend in the duplex DNA. To understand the structural basis for the ultra-high-affinity aptamer binding, an integrative structural biology approach was used. Hydrogen-deuterium exchange coupled to liquid chromatography-mass spectrometry (HDX-MS) was used to examine the structural dynamics of RT alone and in the presence of the DNA aptamer. RT was selectively labeled with 15N to unambiguously identify peptides from each subunit. HDX of unliganded RT shows a mostly stable core. The p66 fingers and thumb subdomains, and the RNase H domain are relatively dynamic. HDX indicates that both the aptamer and a scrambled version significantly stabilize regions of RT that are dynamic in the absence of DNA. No substantial differences in RT dynamics are observed between aptamer and scrambled aptamer binding, despite a large difference in binding affinity. Small-angle X-ray scattering and circular dichroism spectroscopy were used to investigate the aptamer conformation in solution and revealed a pre-bent DNA that possesses both A- and B-form helical character. Both the 2'-O-methyl modifications and the GC tract appear to contribute to an energetically favorable conformation for binding to RT that contributes to the aptamer's ultra-high affinity for RT. The X-ray structure of RT with an RNA/DNA version of the aptamer at 2.8 Å resolution revealed a potential role of the hairpin positioning in affinity. Together, the data suggest that both the 2'-O-methyl modifications and the GC tract contribute to an energetically favorable conformation for high-affinity binding to RT.

Non-Nucleoside Reverse Transcriptase Inhibitors Join Forces with Integrase Inhibitors to Combat HIV

In the treatment of acquired immune deficiency syndrome (AIDS), the diarylpyrimidine (DAPY) analogs etravirine (ETR) and rilpivirine (RPV) have been widely effective against human immunodeficiency virus (HIV) variants that are resistant to other non-nucleoside reverse transcriptase inhibitors (NNRTIs). With non-inferior or improved efficacy, better safety profiles, and lower doses or pill burdens than other NNRTIs in the clinic, combination therapies including either of these two drugs have led to higher adherence than other NNRTI-containing treatments. In a separate development, HIV integrase strand transfer inhibitors (INSTIs) have shown efficacy in treating AIDS, including raltegravir (RAL), elvitegravir (EVG), cabotegravir (CAB), bictegravir (BIC), and dolutegravir (DTG). Of these, DTG and BIC perform better against a wide range of resistance mutations than other INSTIs. Nevertheless, drug-resistant combinations of mutations have begun to emerge against all DAPYs and INSTIs, attributable in part to non-adherence. New dual therapies that may promote better adherence combine ETR or RPV with an INSTI and have been safer and non-inferior to more traditional triple-drug treatments. Long-acting dual- and triple-therapies combining ETR or RPV with INSTIs are under study and may further improve adherence. Here, highly resistant emergent mutations and efficacy data on these novel treatments are reviewed. Overall, ETR or RPV, in combination with INSTIs, may be treatments of choice as long-term maintenance therapies that optimize efficacy, adherence, and safety.

Binding interface and impact on protease cleavage for an RNA aptamer to HIV-1 reverse transcriptase

RNA aptamers that bind HIV-1 reverse transcriptase (RT) inhibit RT in enzymatic and viral replication assays. Some aptamers inhibit RT from only a few viral clades, while others show broad-spectrum inhibition. Biophysical determinants of recognition specificity are poorly understood. We investigated the interface between HIV-1 RT and a broad–spectrum UCAA-family aptamer. SAR and hydroxyl radical probing identified aptamer structural elements critical for inhibition and established the role of signature UCAA bulge motif in RT-aptamer interaction. HDX footprinting on RT ± aptamer shows strong contacts with both subunits, especially near the C-terminus of p51. Alanine scanning revealed decreased inhibition by the aptamer for mutants P420A, L422A and K424A. 2D proton nuclear magnetic resonance and SAXS data provided constraints on the solution structure of the aptamer and enable computational modeling of the docked complex with RT. Surprisingly, the aptamer enhanced proteolytic cleavage of precursor p66/p66 by HIV-1 protease, suggesting that it stabilizes the productive conformation to allow maturation. These results illuminate features at the RT-aptamer interface that govern recognition specificity by a broad-spectrum antiviral aptamer, and they open new possibilities for accelerating RT maturation and interfering with viral replication.

Enumeration of Total Bacteria and Coliforms in Milk by Dry Rehydratable Film Methods: Collaborative Study

Eleven laboratories participated in a collaborative study to compare the dry rehydratable film (Petrifilm® SM and Petrifilm® VRB) methods, respectively, to the standard plate count (SPC) and violet red bile agar (VRBA) standard methods for estimation of total bacteria and coliform counts in raw and homogenized pasteurized milk. Each laboratory analyzed 16 samples (8 different samples in blind duplicate) for total count by both the SPC and Petrifilm SM methods. A second set of 16 samples was analyzed by the VRBA and Petrifilm VRB methods. The repeatability standard deviations (the square root of the between-replicates variance) of the SPC, Petrifilm SM, VRBA, and Petrifilm VRB methods were 0.0S104, 0.0444, 0.14606, and 0.13806, respectively; the reproducibility standard deviations were 0.7197, C.06380, 0.15326, and 0.13806, respectively. The difference between the mean Iog10 SPC and the mean logio Petrifilm SM results was 0.027. For the VRBA and Petrifilm VRB methods, the mean log10 difference was 0.013. These results generally indicate the suitability of the dry rehydratable film methods as alternatives to the SPC and VRBA methods for milk samples. The methods have been adopted official first action.

Discovery and Characterization of Fluorine-Substituted Diarylpyrimidine Derivatives as Novel HIV-1 NNRTIs with Highly Improved Resistance Profiles and Low Activity for the hERG Ion Channel

Our previous efforts have led to the development of two potent NNRTIs, K-5a2 and 25a, exhibiting effective anti-HIV-1 potency and resistance profiles compared with etravirine. However, both inhibitors suffered from potent hERG inhibition and short half-life. In this article, with K-5a2 and etravirine as leads, series of novel fluorine-substituted diarylpyrimidine derivatives were designed via molecular hybridization and bioisosterism strategies. The results indicated 24b was the most active inhibitor, exhibiting broad-spectrum activity (EC50 = 3.60-21.5 nM) against resistant strains, significantly lower cytotoxicity (CC50= 155 μM), and reduced hERG inhibition (IC50 > 30 μM). Crystallographic studies confirmed the binding of 24b and the role of the fluorine atom, as well as optimal contacts of a nitrile group with the main-chain carbonyl group of H235. Furthermore, 24b showed longer half-life and favorable safety properties. All the results demonstrated that 24b has significant promise in circumventing drug resistance as an anti-HIV-1 candidate.

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

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

Structural Basis of HIV-1 Inhibition by Nucleotide-Competing Reverse Transcriptase Inhibitor INDOPY-1

HIV-1 reverse transcriptase (RT) is an essential enzyme, targeting half of approved anti-AIDS drugs. While nucleoside RT inhibitors (NRTIs) are DNA chain terminators, the nucleotide-competing RT inhibitor (NcRTI) INDOPY-1 blocks dNTP binding to RT. Lack of structural information hindered INDOPY-1 improvement. Here we report the HIV-1 RT/DNA/INDOPY-1 crystal structure, revealing a unique mode of inhibitor binding at the polymerase active site without involving catalytic metal ions. The structure may enable new strategies for developing NcRTIs.

Epistasis and entrenchment of drug resistance in HIV-1 subtype B

The development of drug resistance in HIV is the result of primary mutations whose effects on viral fitness depend on the entire genetic background, a phenomenon called 'epistasis'. Based on protein sequences derived from drug-experienced patients in the Stanford HIV database, we use a co-evolutionary (Potts) Hamiltonian model to provide direct confirmation of epistasis involving many simultaneous mutations. Building on earlier work, we show that primary mutations leading to drug resistance can become highly favored (or entrenched) by the complex mutation patterns arising in response to drug therapy despite being disfavored in the wild-type background, and provide the first confirmation of entrenchment for all three drug-target proteins: protease, reverse transcriptase, and integrase; a comparative analysis reveals that NNRTI-induced mutations behave differently from the others. We further show that the likelihood of resistance mutations can vary widely in patient populations, and from the population average compared to specific molecular clones.

Aryl and Arylalkyl Substituted 3-Hydroxypyridin-2(1H)-ones: Synthesis and Evaluation as Inhibitors of Influenza A Endonuclease

Seasonal influenza infections are associated with an estimated 250-500 000 deaths annually. Resistance to the antiviral M2 ion-channel inhibitors has largely invalidated their clinical utility. Resistance to neuraminidase inhibitors has also been observed in several influenza A virus (IAV) strains. These data have prompted research on inhibitors that target the cap-snatching endonuclease activity of the polymerase acidic protein (PA). Baloxavir marboxil (Xofluza®), recently approved for clinical use, inhibits cap-snatching endonuclease. Resistance to Xofluza® has been reported in both in vitro systems and in the clinic. An X-ray crystallographic screening campaign of a fragment library targeting IAV endonuclease identified 5-chloro-3-hydroxypyridin-2(1H)-one as a bimetal chelating agent at the active site. We have reported the structure-activity relationships for 3-hydroxypyridin-2(1H)-ones and 3-hydroxyquinolin-2(1H)-ones as endonuclease inhibitors. These studies identified two distinct binding modes associated with inhibition of this enzyme that are influenced by the presence of substituents at the 5- and 6-positions of 3-hydroxypyridin-2(1H)-ones. Herein we report the structure-activity relationships associated with various para-substituted 5-phenyl derivatives of 6-(p-fluorophenyl)-3-hydroxypyridin-2(1H)-ones and the effect of using naphthyl, benzyl, and naphthylmethyl groups as alternatives to the p-fluorophenyl substituent on their activity as endonuclease inhibitors.

Michael G. Rossmann (1930-2019), pioneer in macromolecular and virus crystallography: scientist, mentor and friend

Michael George Rossmann, who made monumental contributions to science, passed away peacefully in West Lafayette, Indiana on 14 May 2019 at the age of 88, following a courageous five-year battle with cancer. Michael was born in Frankfurt, Germany on 30 July 1930. As a young boy, he emigrated to England with his mother just as World War II ignited. Michael was a highly innovative and energetic person, well known for his intensity, persistence and focus in pursuing his research goals. Michael was a towering figure in crystallography as a highly distinguished faculty member at Purdue University for 55 years. Michael made many seminal contributions to crystallography in a career that spanned the entirety of structural biology, beginning in the 1950s at Cambridge where the first protein structures were determined in the laboratories of Max Perutz (hemoglobin, 1960) and John Kendrew (myoglobin, 1958). Michael's work was central in establishing and defining the field of structural biology, which amazingly has described the structures of a vast array of complex biological molecules and assemblies in atomic detail. Knowledge of three-dimensional biological structure has important biomedical significance including understanding the basis of health and disease at the molecular level, and facilitating the discovery of many drugs.

Structure of HIV-1 reverse transcriptase/d4TTP complex: Novel DNA cross-linking site and pH-dependent conformational changes

Stavudine (d4T, 2',3'-didehydro-2',3'-dideoxythymidine) was one of the first chain-terminating nucleoside analogs used to treat HIV infection. We present the first structure of the active, triphosphate form of d4T (d4TTP) bound to a catalytic complex of HIV-1 RT/dsDNA template-primer. We also present a new strategy for disulfide (S-S) chemical cross-linking between N⁶ of a modified adenine at the second overhang base to I63C in the fingers subdomain of RT. The cross-link site is upstream of the duplex-binding region of RT, however, the structure is very similar to published RT structures with cross-linking to Q258C in the thumb, which suggests that cross-linking at either site does not appreciably perturb the RT/DNA structures. RT has a catalytic maximum at pH 7.5. We determined the X-ray structures of the I63C-RT/dsDNA/d4TTP cross-linked complexes at pH 7, 7.5, 8, 8.5, 9, and 9.5. We found small (~0.5 Å), pH-dependent motions of the fingers subdomain that folds in to form the dNTP-binding pocket. We propose that the pH-activity profile of RT relates to this motion of the fingers. Due to side effects of neuropathy and lipodystrophy, use of d4T has been stopped in most countries, however, chemical modification of d4T might lead to the development of a new class of nucleoside analogs targeting RNA and DNA polymerases.

Publisher Correction: Structural Implications of Mutations Conferring Rifampin Resistance in Mycobacterium leprae

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Structural Implications of Mutations Conferring Rifampin Resistance in Mycobacterium leprae

The rpoB gene encodes the β subunit of RNA polymerase holoenzyme in Mycobacterium leprae (M. leprae). Missense mutations in the rpoB gene were identified as etiological factors for rifampin resistance in leprosy. In the present study, we identified mutations corresponding to rifampin resistance in relapsed leprosy cases from three hospitals in southern India which treat leprosy patients. DNA was extracted from skin biopsies of 35 relapse/multidrug therapy non-respondent leprosy cases, and PCR was performed to amplify the 276 bp rifampin resistance-determining region of the rpoB gene. PCR products were sequenced, and mutations were identified in four out of the 35 cases at codon positions D441Y, D441V, S437L and H476R. The structural and functional effects of these mutations were assessed in the context of three-dimensional comparative models of wild-type and mutant M. leprae RNA polymerase holoenzyme (RNAP), based on the recently solved crystal structures of RNAP of Mycobacterium tuberculosis, containing a synthetic nucleic acid scaffold and rifampin. The resistance mutations were observed to alter the hydrogen-bonding and hydrophobic interactions of rifampin and the 5' ribonucleotide of the growing RNA transcript. This study demonstrates that rifampin-resistant strains of M. leprae among leprosy patients in southern India are likely to arise from mutations that affect the drug-binding site and stability of RNAP.

Structural Basis of Mycobacterium tuberculosis Transcription and Transcription Inhibition

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis, which kills 1.8 million annually. Mtb RNA polymerase (RNAP) is the target of the first-line antituberculosis drug rifampin (Rif). We report crystal structures of Mtb RNAP, alone and in complex with Rif, at 3.8-4.4 Å resolution. The results identify an Mtb-specific structural module of Mtb RNAP and establish that Rif functions by a steric-occlusion mechanism that prevents extension of RNA. We also report non-Rif-related compounds-Nα-aroyl-N-aryl-phenylalaninamides (AAPs)-that potently and selectively inhibit Mtb RNAP and Mtb growth, and we report crystal structures of Mtb RNAP in complex with AAPs. AAPs bind to a different site on Mtb RNAP than Rif, exhibit no cross-resistance with Rif, function additively when co-administered with Rif, and suppress resistance emergence when co-administered with Rif.

A New Class of Allosteric HIV-1 Integrase Inhibitors Identified by Crystallographic Fragment Screening of the Catalytic Core Domain

HIV-1 integrase (IN) is essential for virus replication and represents an important multifunctional therapeutic target. Recently discovered quinoline-based allosteric IN inhibitors (ALLINIs) potently impair HIV-1 replication and are currently in clinical trials. ALLINIs exhibit a multimodal mechanism of action by inducing aberrant IN multimerization during virion morphogenesis and by competing with IN for binding to its cognate cellular cofactor LEDGF/p75 during early steps of HIV-1 infection. However, quinoline-based ALLINIs impose a low genetic barrier for the evolution of resistant phenotypes, which highlights a need for discovery of second-generation inhibitors. Using crystallographic screening of a library of 971 fragments against the HIV-1 IN catalytic core domain (CCD) followed by a fragment expansion approach, we have identified thiophenecarboxylic acid derivatives that bind at the CCD-CCD dimer interface at the principal lens epithelium-derived growth factor (LEDGF)/p75 binding pocket. The most active derivative (5) inhibited LEDGF/p75-dependent HIV-1 IN activity in vitro with an IC₅₀ of 72 μm and impaired HIV-1 infection of T cells at an EC₅₀ of 36 μm The identified lead compound, with a relatively small molecular weight (221 Da), provides an optimal building block for developing a new class of inhibitors. Furthermore, although structurally distinct thiophenecarboxylic acid derivatives target a similar pocket at the IN dimer interface as the quinoline-based ALLINIs, the lead compound, 5, inhibited IN mutants that confer resistance to quinoline-based compounds. Collectively, our findings provide a plausible path for structure-based development of second-generation ALLINIs.

Conformational States of HIV-1 Reverse Transcriptase for Nucleotide Incorporation vs Pyrophosphorolysis-Binding of Foscarnet

HIV-1 reverse transcriptase (RT) catalytically incorporates individual nucleotides into a viral DNA strand complementing an RNA or DNA template strand; the polymerase active site of RT adopts multiple conformational and structural states while performing this task. The states associated are dNTP binding at the N site, catalytic incorporation of a nucleotide, release of a pyrophosphate, and translocation of the primer 3'-end to the P site. Structural characterization of each of these states may help in understanding the molecular mechanisms of drug activity and resistance and in developing new RT inhibitors. Using a 38-mer DNA template-primer aptamer as the substrate mimic, we crystallized an RT/dsDNA complex that is catalytically active, yet translocation-incompetent in crystals. The ability of RT to perform dNTP binding and incorporation in crystals permitted obtaining a series of structures: (I) RT/DNA (P-site), (II) RT/DNA/AZTTP ternary, (III) RT/AZT-terminated DNA (N-site), and (IV) RT/AZT-terminated DNA (N-site)/foscarnet complexes. The stable N-site complex permitted the binding of foscarnet as a pyrophosphate mimic. The Mg(2+) ions dissociated after catalytic addition of AZTMP in the pretranslocated structure III, whereas ions A and B had re-entered the active site to bind foscarnet in structure IV. The binding of foscarnet involves chelation with the Mg(2+) (B) ion and interactions with K65 and R72. The analysis of interactions of foscarnet and the recently discovered nucleotide-competing RT inhibitor (NcRTI) α-T-CNP in two different conformational states of the enzyme provides insights for developing new classes of polymerase active site RT inhibitors.

Exploring the role of the α-carboxyphosphonate moiety in the HIV-RT activity of α-carboxy nucleoside phosphonates

As α-carboxy nucleoside phosphonates (α-CNPs) have demonstrated a novel mode of action of HIV-1 reverse transcriptase inhibition, structurally related derivatives were synthesized, namely the malonate 2, the unsaturated and saturated bisphosphonates 3 and 4, respectively and the amide 5. These compounds were evaluated for inhibition of HIV-1 reverse transcriptase in cell-free assays. The importance of the α-carboxy phosphonoacetic acid moiety for achieving reverse transcriptase inhibition, without the need for prior phosphorylation, was confirmed. The malonate derivative 2 was less active by two orders of magnitude than the original α-CNPs, while displaying the same pattern of kinetic behavior; interestingly the activity resides in the “L”-enantiomer of 2, as seen with the earlier series of α-CNPs. A crystal structure with an RT/DNA complex at 2.95 Å resolution revealed the binding of the “L”-enantiomer of 2, at the polymerase active site with a weaker metal ion chelation environment compared to 1a (T-α-CNP) which may explain the lower inhibitory activity of 2.

Rapid experimental SAD phasing and hot-spot identification with halogenated fragments

Through X-ray crystallographic fragment screening, 4-bromopyrazole was discovered to be a 'magic bullet' that is capable of binding at many of the ligand 'hot spots' found in HIV-1 reverse transcriptase (RT). The binding locations can be in pockets that are 'hidden' in the unliganded crystal form, allowing rapid identification of these sites for in silico screening. In addition to hot-spot identification, this ubiquitous yet specific binding provides an avenue for X-ray crystallographic phase determination, which can be a significant bottleneck in the determination of the structures of novel proteins. The anomalous signal from 4-bromopyrazole or 4-iodopyrazole was sufficient to determine the structures of three proteins (HIV-1 RT, influenza A endonuclease and proteinase K) by single-wavelength anomalous dispersion (SAD) from single crystals. Both compounds are inexpensive, readily available, safe and very soluble in DMSO or water, allowing efficient soaking into crystals.

Factors influencing the efficacy of rilpivirine in HIV-1 subtype C in low- and middle-income countries

OBJECTIVES

The use of the NNRTI rilpivirine in low- and middle-income countries (LMICs) is under debate. The main objective of this study was to provide further clinical insights and biochemical evidence on the usefulness of rilpivirine in LMICs.

PATIENTS AND METHODS

Rilpivirine resistance was assessed in 5340 therapy-naive and 13,750 first-generation NNRTI-failed patients from Europe and therapy-naive HIV-1 subtype C (HIV-1C)-infected individuals from India (n = 617) and Ethiopia (n = 127). Rilpivirine inhibition and binding affinity assays were performed using patient-derived HIV-1C reverse transcriptases (RTs).

RESULTS

Primary rilpivirine resistance was rare, but the proportion of patients with >100,000 HIV-1 RNA copies/mL pre-ART was high in patients from India and Ethiopia, limiting the usefulness of rilpivirine as a first-line drug in LMICs. In patients failing first-line NNRTI treatments, cross-resistance patterns suggested that 73% of the patients could benefit from switching to rilpivirine-based therapy. In vitro inhibition assays showed ∼ 2-fold higher rilpivirine IC50 for HIV-1C RT than HIV-1B RT. Pre-steady-state determination of rilpivirine-binding affinities revealed 3.7-fold lower rilpivirine binding to HIV-1C than HIV-1B RT. Structural analysis indicated that naturally occurring polymorphisms close to the NNRTI-binding pocket may reduce rilpivirine binding, leading to lower susceptibility of HIV-1C to rilpivirine.

CONCLUSIONS

Our clinical and biochemical findings indicate that the usefulness of rilpivirine has limitations in HIV-1C-dominated epidemics in LMICs, but the drug could still be beneficial in patients failing first-line therapy if genotypic resistance testing is performed.

Negative-Strand RNA Virus L Proteins: One Machine, Many Activities

Structures of L proteins from La Crosse orthobunyavirus and vesicular stomatitis virus reveal insights into RNA synthesis and distinctive mRNA capping mechanisms of segmented and non-segmented negative-sense single-strand RNA viruses.

Pronounced Inhibition Shift from HIV Reverse Transcriptase to Herpetic DNA Polymerases by Increasing the Flexibility of α-Carboxy Nucleoside Phosphonates

Alpha-carboxynucleoside phosphonates (α-CNPs) are novel viral DNA polymerase inhibitors that do not need metabolic conversion for enzyme inhibition. The prototype contains a cyclopentyl linker between nucleobase and α-carboxyphosphonate and preferentially (50- to 100-fold) inhibits HIV-1 RT compared with herpetic DNA polymerases. A synthesis methodology involving three steps has been developed for the synthesis of a series of novel α-CNPs, including a Rh(II)-catalyzed O-H insertion that connects the carboxyphosphonate group to a linker moiety and an attachment of a nucleobase to the other end of the linker by a Mitsunobu reaction followed by final deprotection. Replacing the cyclopentyl moiety in the prototype α-CNPs by a more flexible entity results in a selectivity shift of ∼ 100-fold in favor of the herpetic DNA polymerases when compared to selectivity for HIV-1 RT. The nature of the kinetic interaction of the acyclic α-CNPs against the herpetic DNA polymerases differs from the nature of the nucleobase-specific kinetic interaction of the cyclopentyl α-CNPs against HIV RT.

Structure of HIV-1 reverse transcriptase bound to a novel 38-mer hairpin template-primer DNA aptamer

The development of a modified DNA aptamer that binds HIV-1 reverse transcriptase (RT) with ultra-high affinity has enabled the X-ray structure determination of an HIV-1 RT-DNA complex to 2.3 Å resolution without the need for an antibody Fab fragment or RT-DNA cross-linking. The 38-mer hairpin-DNA aptamer has a 15 base-pair duplex, a three-deoxythymidine hairpin loop, and a five-nucleotide 5'-overhang. The aptamer binds RT in a template-primer configuration with the 3'-end positioned at the polymerase active site and has 2'-O-methyl modifications at the second and fourth duplex template nucleotides that interact with the p66 fingers and palm subdomains. This structure represents the highest resolution RT-nucleic acid structure to date. The RT-aptamer complex is catalytically active and can serve as a platform for studying fundamental RT mechanisms and for development of anti-HIV inhibitors through fragment screening and other approaches. Additionally, the structure allows for a detailed look at a unique aptamer design and provides the molecular basis for its remarkably high affinity for RT.

Differential isotopic enrichment to facilitate characterization of asymmetric multimeric proteins using hydrogen/deuterium exchange mass spectrometry

Hydrogen/deuterium exchange (HDX) coupled to mass spectrometry has emerged as a powerful tool for analyzing the conformational dynamics of protein-ligand and protein-protein interactions. Recent advances in instrumentation and methodology have expanded the utility of HDX for the analysis of large and complex proteins; however, asymmetric dimers with shared amino acid sequence present a unique challenge for HDX because assignment of peptides with identical sequence to their subunit of origin remains ambiguous. Here we report the use of differential isotopic labeling to facilitate HDX analysis of multimers using HIV-1 reverse transcriptase (RT) as a model. RT is an asymmetric heterodimer of 51 kDa (p51) and 66 kDa (p66) subunits. The first 440 residues of p51 and p66 are identical. In this study differentially labeled RT was reconstituted from isotopically enriched ((15)N-labeled) p51 and unlabeled p66. To enable detection of (15)N-deuterated RT peptides, the software HDX Workbench was modified to follow a 100% (15)N model. Our results demonstrated that (15)N enrichment of p51 did not affect its conformational dynamics compared to unlabeled p51, but (15)N-labeled p51 did show different conformational dynamics than p66 in the RT heterodimer. Differential HDX-MS of isotopically labeled RT in the presence of the non-nucleoside reverse transcriptase inhibitor (NNRTI) efavirenz (EFV) showed subunit-specific perturbation in the rate of HDX consistent with previously published results and the RT-EFV cocrystal structure.