Structural basis for potent and broad inhibition of HIV-1 RT by thiophene[3,2-d]pyrimidine non-nucleoside inhibitors
Abstract
Rapid generation of drug-resistant mutations in HIV-1 reverse transcriptase (RT), a prime target for anti-HIV therapy, poses a major impediment to effective anti-HIV treatment. Our previous efforts have led to the development of two novel non-nucleoside reverse transcriptase inhibitors (NNRTIs) with piperidine-substituted thiophene[3,2-d]pyrimidine scaffolds, compounds K-5a2 and 25a, which demonstrate highly potent anti-HIV-1 activities and improved resistance profiles compared with etravirine and rilpivirine, respectively. Here, we have determined the crystal structures of HIV-1 wild-type (WT) RT and seven RT variants bearing prevalent drug-resistant mutations in complex with K-5a2 or 25a at ~2 Å resolution. These high-resolution structures illustrate the molecular details of the extensive hydrophobic interactions and the network of main chain hydrogen bonds formed between the NNRTIs and the RT inhibitor-binding pocket, and provide valuable insights into the favorable structural features that can be employed for designing NNRTIs that are broadly active against drug-resistant HIV-1 variants.
https://doi.org/10.7554/eLife.36340.001Introduction
HIV-1 reverse transcriptase (RT) (hereinafter referred to as RT) plays an essential role in the viral life cycle by reverse transcribing the single-stranded RNA genome to a double-stranded DNA copy (Deeks et al., 2015; Engelman and Cherepanov, 2012). For this reason, it has been an important target of anti-HIV therapies (Esté and Cihlar, 2010; Gubernick et al., 2016). There are two main types of RT inhibitors: nucleoside RT inhibitors (NRTIs), which act as chain terminators and compete with incoming nucleotides in the polymerase active site (Ren et al., 1998; Sarafianos et al., 1999; Tu et al., 2010; Yarchoan et al., 1988), and non-nucleoside RT inhibitors (NNRTIs), which inhibit the activity of RT noncompetitively (Merluzzi et al., 1990; Spence et al., 1995). NNRTIs are a group of structurally diverse compounds that bind to the non-nucleoside inhibitor-binding pocket (NNIBP) located ~10 Å from the polymerase active site (Ding et al., 1995; Kohlstaedt et al., 1992; Ren et al., 1995). NNIBP is a hydrophobic pocket that emerges only when NNRTIs bind and induce conformational rearrangements of the residues defining the pocket (Esnouf et al., 1995; Hsiou et al., 1996). NNRTIs are key components in highly active antiretroviral therapy (HAART) due to their high specificity, desirable pharmacokinetics and generally good tolerance (Moore and Chaisson, 1999; Pomerantz and Horn, 2003).
Despite the success of NNRTIs in suppressing HIV-1 replication and reducing viral loads, their effectiveness is compromised by the emergence of drug-resistant mutations in RT (Wainberg et al., 2011). Earlier NNRTIs, including nevirapine (NVP), delavirdine (DLV) and efavirenz (EFV), have low genetic barriers for resistance and are extremely susceptible to mutations in the NNIBP of RT (Das and Arnold, 2013a, 2013b). K103N, Y181C and Y188L are among the most prevalent NNRTI-resistant mutations identified in RT (de Béthune, 2010; Wensing et al., 2017). Y181C and Y188L mutations introduce steric hindrances between NNRTIs and the pocket, and/or eliminate critical π-π stacking interactions between side chains of the two tyrosine residues and the aromatic rings in NNRTIs (Hsiou et al., 1998; Ren et al., 2001). As to the K103N mutation, it was long believed that it prevented the entry of NNRTIs by stabilizing the closed conformation of NNIBP (Hsiou et al., 2001). However, a more recent study indicates that the resistance is more likely caused by the electrostatic difference between Asn103 and Lys103 (Lai et al., 2016). In the light of this new piece of data, the K103N mutation seems to utilize the same mechanism as Y181C and Y188L do to confer resistance to NNRTIs: by altering the shape or surface property of the NNIBP.
Next-generation NNRTIs are designed with conformational flexibility and positional adaptability and are able to target the NNIBPs of an array of drug-resistant RT mutants (Das et al., 2008, 2004). Etravirine (ETR, also known as TMC125) and rilpivirine (RPV, also known as TMC278) are two U.S. Food and Drug Administration (FDA)-approved second-generation NNRTIs belonging to the diarylpyrimidine (DAPY) family (Figure 1). Both drugs show potent antiviral activities against wild-type (WT) HIV-1 and many HIV-1 variants displaying significant resistance to first-generation NNRTIs (Janssen et al., 2005; Ludovici et al., 2001). However, some existing resistance-associated RT mutations, such as K101P and Y181I, can still cause substantial decreases in susceptibility to ETR and RPV (Azijn et al., 2010; Giacobbi and Sluis-Cremer, 2017; Smith et al., 2016). Besides, new resistant mutations can arise from prolonged use of ETR and RPV, which undermine their anti-HIV-1 activities (Wensing et al., 2017). In patients who failed ETR- or RPV-based therapies, E138K/Q/R are among the most frequently occurred mutations identified in RT (Xu et al., 2013). Therefore, it is imperative to develop new NNRTIs with improved drug-resistance profiles.
Our previous efforts have led to the design and synthesis of two piperidine-substituted thiophene[3,2-d]pyrimidine NNRTIs using ETR as a lead compound (Figure 1) (Kang et al., 2017, 2016). Compound K-5a2 features a thiophene[3,2-d]pyrimidine central ring, and replaces the cyanophenyl right wing of ETR with a more extended piperidine-linked benzenesulfonamide group, while keeping the 4-cyano-2,6-dimethylphenyl structure in the left wing of ETR. Compound 25a shares the same central ring and right wing structures with K-5a2, but grafts the 4-cyanovinyl-2,6-dimethylphenyl structure of RPV onto its left wing. Compared with ETR, compound K-5a2 displays much lower cytotoxicity and increased anti-HIV-1 potency against WT virus and virus strains with a variety of NNRTI-resistant mutations, except K103N and K103N/Y181C (Kang et al., 2016). The further optimized compound 25a is exceptionally potent in inhibiting WT HIV-1 and exhibits significantly better anti-HIV-1 activities than ETR against all of the tested NNRTI-resistant HIV-1 strains in cellular assays (Kang et al., 2017).
In this study, we demonstrated that 25a is superior to RPV in inhibiting RT bearing a wide range of resistance mutations, including K101P, Y181I and K103N/Y181I, against which RPV loses considerable potency, and determined the crystal structures of WT and mutant RTs in complex with either K-5a2 or 25a. These structures illustrate the detailed interactions between RT and the two inhibitors, and explain why K-5a2 and 25a are resilient to NNRTI-resistant mutations in the NNIBP. Additionally, comparison of the binding modes of K-5a2 and 25a with those of ETR and RPV suggests the possible mechanisms for the susceptibilities of ETR and RPV to E138K and K101P mutations. Our results outline the structural features of NNRTIs that can be employed for future drug design to overcome prevalent NNRTI-resistant mutations.
Results
Structure determination
The complexes of 25a or K-5a2 bound to WT RT or RT variants with drug-resistant mutations were prepared by soaking either NNRTI into the RT crystals. The structures were determined by molecular replacement using the structure of WT RT/RPV complex (PDB ID: 4G1Q) as the search template and were subsequently refined to 1.9–2.23 Å resolution (Supplementary file 1). Overall, the structure of RT in the complexes has the same ‘open-cleft’ conformation as observed in prior RT/NNRTIs structures (Figure 2A and C) (Das et al., 2008, 2004; Ding et al., 1995; Ren et al., 1995). The Fo-Fc electron-density maps unambiguously defined the binding positions and conformations of both inhibitors in the NNIBP (Figure 2B and D and Figure 2—figure supplement 1).
Interactions between piperidine-substituted thiophene[3,2-d]pyrimidine NNRTIs and RT
The RT-bound 25a and K-5a2 adopt a horseshoe conformation, which is similar to that seen with NNRTIs in the DAPY family (Das et al., 2008, 2004). Both inhibitors exhibit remarkable structural complementarity to the NNIBP with substantial extensions into the three channels (tunnel, entrance and groove) characterizing the pocket (Figure 3A and B and Figure 3—figure supplement 1A and B). The left wing structures of 25a and K-5a2 form hydrophobic interactions with Pro95 and Leu234, and project into the tunnel lined by Tyr181, Tyr188, Phe227, and Trp229, forming π-π interactions with these residues. The entrance channel gated by Glu138 in the p51 subunit and Lys101 in the p66 subunit is an underexplored region in the NNIBP. By substituting the central pyrimidine ring of DAPY NNRTIs to a thiophene[3,2-d]pyrimidine heterocyclic structure, 25a and K-5a2 are able to establish nonpolar interactions with the alkyl chain of Glu138, while retaining the favorable hydrophobic contacts with Val179 and Leu100 manifested in the complexes of RT with ETR or RPV (Das et al., 2008, 2004). The piperidine-linked aryl structure of the right wing arches into the groove surrounded by Lys103, Val106, Pro225, Phe227, Pro236, and Tyr318, developing numerous van der Waals contacts with their lipophilic side chains, and directs the terminal sulfonamide group to the solvent-exposed surface of RT.
In addition, binding of 25a and K-5a2 to the NNIBP is stabilized by an extensive hydrogen-bonding network between the inhibitors and the main chains of several key residues around the pocket (Figure 3C and D and Figure 3—figure supplement 1C and D): (i) the surface-positioned sulfonamide group is double-hydrogen bonded with the carbonyl oxygen of Lys104 and the backbone nitrogen of Val106; (ii) the piperidine nitrogen forms hydrogen bonds with the main chains of Lys103 and Pro236 through a bridging water molecule; (iii) the amine group linking the central thiophene pyrimidine and the piperidine ring interacts with the carbonyl oxygen of Lys101, forming a conserved hydrogen bond observed in a number of second-generation NNRTIs/RT complexes (Das et al., 2008; Lansdon et al., 2010); (iv) additionally, the nitrogen and sulfur atoms in the central thiophene pyrimidine ring are involved in two water-mediated hydrogen bonds with the backbone nitrogen of Lys101 and the carbonyl oxygens of Glu138, respectively. These extensive interactions between the two piperidine-substituted thiophene[3,2-d]pyrimidine NNRTIs and RT lock the enzyme in an open-cleft conformation and inhibit its polymerization activity.
The above interactions between RT and the two NNRTIs generally agree with the results from molecular docking (Kang et al., 2017, 2016). Nevertheless, a close inspection of the inhibitors observed in the crystal structures and those predicted by molecular docking reveals a few notable differences in their binding modes. First, the thiophene pyrimidine nitrogen in the inhibitors is not directly hydrogen bonded with the backbone nitrogen of Lys101, as predicted by inhibitor docking, but through a bridging water molecule instead. Second, the crystal structures define a water-mediated hydrogen bond between the carboxyl oxygen of Glu138 and the sulfur group in the inhibitors central ring, which is absent in the predicted binding modes. Due to their free movement and transient involvement in the binding process, it is difficult to predict the role of solvent molecules in the interactions between enzymes and inhibitors using ligand-docking programs. These water-mediated interactions, however, can be critical for the enzyme–inhibitor complex formation and thus can provide important insights in understanding the resistance mechanisms of RT mutants.
Inhibition of HIV-1 RT by piperidine-substituted thiophene[3,2-d]pyrimidine NNRTIs
Our previous MT-4 cell-based antiviral activity evaluations showed that K-5a2 displays ~3-fold greater efficacy than ETR against the WT HIV-1 strain, and higher or similar efficacy against virus variants bearing four prevalent single-residue mutations (L100I, K103N, E138K and Y181C) in RT (hereinafter referred to as L100I RT, K103N RT, E138K RT and Y181C RT, respectively). However, K-5a2 is less effective than ETR in inhibiting HIV-1 strains containing K103N/Y181C RT (Kang et al., 2016). The compound 25a, resulting from further optimization of K-5a2, overcomes the limitations of K-5a2 and exhibits significantly better inhibitory effects on all tested HIV-1 strains (Kang et al., 2017). To better compare the anti-HIV-1 potency of 25a with that of existing NNRTIs, we measured the EC50 values of RPV towards WT HIV-1 and mutant RT-bearing HIV-1 variants using the same method. 25a holds advantages over RPV in most of the tested drug-resistant HIV-1 strains while retaining similar antiviral potency against HIV-1 strains containing WT RT (Table 1). It is noteworthy that the particularly challenging K103N/Y181C double-mutation only causes a ~4.6-fold change in susceptibility to 25a, whereas it reduces the anti-HIV-1 efficacy of RPV by more than 10-fold. The superiority of 25a over RPV in targeting K103N/Y181C RT was further validated in the in vitro RT inhibition assay using purified recombinant RT variants, where K103N/Y181C mutation confers lower level of resistance to 25a (7.2-fold change in the IC50 value) than to RPV (15-fold change) (Figure 4A and B and Table 2).
To further evaluate the resistance profile of 25a, we compared the RT inhibitory activities of 25a and RPV against two additional clinically relevant RT mutants, Y188L RT and V106A/F227L RT. While 25a and RPV exhibit similar inhibitory potency against WT RT, 25a is more resilient to Y188L and V106A/F227L mutations (0.70- and 1.7-fold change, respectively) than RPV (2.2- and 4.0-fold change, respectively) (Figure 4A and B and Table 2). To assess whether mutations against which RPV loses considerable potency would be susceptible to 25a, we tested the inhibitory activities of 25a and RPV against K101P RT, Y181I RT and K103N/Y181I RT, which were shown to cause substantial reduction in susceptibility to RPV (Azijn et al., 2010; Giacobbi and Sluis-Cremer, 2017; Smith et al., 2016). As expected, all three mutations dramatically lower the anti-RT potency of RPV and cause 20-fold, 90-fold and 1805-fold change in the IC50 values, respectively. In contrast, there is considerably less resistance to 25a for all three RT mutants (1.3-, 8.8- and 96-fold, respectively) (Figure 4C and D and Table 2).
The longer right wing of 25a enables its interactions with NNIBP residues that are not contacted by RPV, such as Pro225 and Pro236. To gauge the likely impact of these mutations on 25a efficacy, we measured the RT-inhibiting potency of 25a against RT containing P225H or P236L substitutions, two clinically identified mutations shown no significant reduction in susceptibility to RPV (Basson et al., 2015). Like that of RPV, the potency of 25a was not negatively affected by either P225H or P236L mutations (0.58- and 0.60-fold change, respectively) (Figure 4E and F and Table 2), indicating that 25a has a relatively high genetic barrier to the development of novel drug-resistant mutations.
It is worth mentioning that 25a has steeper dose-response curve slopes than RPV in the inhibition of all above RT variants (Table 2). This characteristic can help 25a achieve greater inhibition of RT activity at higher than IC50 concentrations, which are usually more clinically relevant (Shen et al., 2008). Taken together, by comparing the inhibitory potency of 25a and RPV in a wide range of RT mutants, we have shown that 25a has an improved resistance profile over RPV and is able to effectively inhibit the RT mutants causing high-level resistance to RPV.
Structural basis for improved resistance profile of piperidine-substituted thiophene[3,2-d]pyrimidine NNRTIs
To shed light on the mechanism underlying the outstanding resistance profile of the two piperidine-substituted thiophene[3,2-d]pyrimidine NNRTIs, we determined the crystal structures of K103N RT, E138K RT, and Y188L RT complexed with compound K-5a2, as well as K103N RT, E138K RT, K103N/Y181C RT, V106A/F227L RT, K101P RT, and Y181I RT complexed with 25a. The attempt to obtain the crystal structure of 25a in complex with K103N/Y181I RT proved unsuccessful, possibly due its suboptimal anti-RT potency towards K103N/Y181I RT, although it has displayed marked improvement over RPV in inhibiting this specific mutant (Table 2). Superposition of these mutant RT/NNRTI complexes structures onto their respective WT RT/NNRTI complexes structures shows no major deviation in the conformations of the enzyme and inhibitors (Figure 5 and Figure 5—figure supplement 1). Root-mean-square deviations (RMSDs) for the structural alignments between WT RT/NNRTI complexes and mutant RT/NNRTI complexes range from 0.094 to 0.283 Å for the overall Cα atoms, and from 0.095 to 1.108 Å for the Cα atoms of the NNIBP regions (residues 98–110, 178–190, 226–240 of the p66 subunit, plus residues 137–139 of the p51 subunit) (Table 3).
Examination of the interactions between the RT mutants and the two NNRTIs reveals that all the hydrogen bonds depicted in Figure 3D are preserved, although there are some variations in the bond lengths. To analyze the extent of interactions between the inhibitors and different RT mutants, we measured the buried surface areas between the inhibitors and the whole NNIBP as well as a selection of key residues in the NNIBP of each RT variants (Table 4). In the structures of K103N RT in complex with 25a or K-5a2, the Lys to Asn substitution in RT shortens the aliphatic side chain and reduces the contact interface between residue 103 and the inhibitors, but 25a and K-5a2 are able to establish more contacts with Phe227 and Pro236 by varying their multiple torsion angles (Table 5) to counterbalance the loss (Figure 5A and Figure 5—figure supplement 1A). Similarly, in the structure of Y188L RT/K-5a2 complex, the cyano-dimethylphenyl group in K-5a2 is diverted away from Leu188 to avoid steric clashes, leading to declines in the buried areas between the inhibitor and Leu188 and Phe227. However, this mutation-caused damage is alleviated by its enhanced interactions with Lys103, Val106 and Pro236 (Figure 5—figure supplement 1C, Table 4). In the case of E138K RT, since the mutation does not disrupt the hydrophobic interactions between the inhibitors’ central thiophene ring and residue 138, 25a and K-5a2 maintain almost the same binding poses in the NNIBP and similar buried areas with each of the residues along the pocket as in their complexes with WT RT (Figure 5B and Figure 5—figure supplement 1B).
In regard to RT carrying the more disruptive K103N/Y181C double-mutation, Y181C mutation abolishes the favorable π–π stacking interactions between the Tyr181 side chain and the dimethylphenyl ring of 25a, and greatly reduce the binding interface between 25a and Cys181. Moreover, the dramatic changes in the NNIBP result in a decrease of the buried interface between 25a and residue 103. Nonetheless, the markedly weakened interactions between 25a and both mutated residues are remedied by the increase in the contact areas between 25a and several other residues in the NNIBP, including Tyr183, Phe227 and Pro236. (Figure 5C). In the V106A/F227L RT/25a complex structure, the much smaller side chain of Ala106 buries significantly less surface area with the inhibitor. Furthermore, the double-mutation causes more dramatic changes in the conformations of the NNIBP and the bound 25a. In particular, the cyanovinyl group of 25a is flipped so that it can maintain similar extent of interactions with the mutated Leu227. This torsional change, however, diverts the inhibitor away from the tunnel lined by Tyr181, Tyr183 and Tyr188 and diminishes the contact areas between 25a and all three tyrosine residues. To compensate for the loss, the right wing of 25a shifts closer to Lys101. Such movement shortens the distance between the linker amine group of 25a and the carboxyl oxygen of Lys101 from 3.05 to 2.76 Å and strengthens the hydrogen bond between them (Figure 5D). This hydrogen bond is conserved in the binding of many NNRTIs, including ETR and RPV, and contributes greatly to the binding affinities of NNRTIs (Das et al., 2008; Lansdon et al., 2010).
With respect to K101P RT, the deprotonation of Pro101 main-chain nitrogen attenuates its water-mediated hydrogen bond with the pyrimidine nitrogen in 25a. Nonetheless, the Lys to Pro substitution places its cyclic side chain in the vicinity of the central thiophene pyrimidine ring of 25a and leads to enhanced hydrophobic interactions (Figure 5E). As to the Y181I RT/25a complex structure, the mutation not only removes the π-π stacking interactions in the left wing of 25a, but also introduces steric hinderance with the linker oxygen and thiophene sulfur group in 25a, which pushes the inhibitor slightly away from the tunnel. This unfavorable change in the NNIBP is mitigated by the enlarged contact areas between 25a and the more closely placed side-chains of Lys101 and Tyr183 (Figure 5F).
It is noteworthy that in both K103N/Y181C RT/25a and Y181I RT/25a complex structures, Tyr183, a residue in the conserved YMDD motif at the polymerase active site, moves 2.3–2.5 Å towards the inhibitor, enhancing its hydrophobic interactions with the cyanovinyl group of 25a. The recruitment of Tyr183 by the cyanovinyl group of 25a is reminiscent of that observed in the structure of K103N/Y181C RT in complex with RPV, whose left wing has the same 4-cyanovinyl-2,6-dimethylphenyl structure (Das et al., 2008). Interestingly, in the structure of K103N RT/K-5a2 complex, Tyr183 also undergoes a significant conformational change and is placed closer to the inhibitor, although it is still outside the contact radius of the cyano group in K-5a2. Superposition of the structures of WT RT/K-5a2, K103N RT/K-5a2, WT RT/25a, and K103N/Y181C RT/25a complexes reveals a gradual rotation of Tyr183 from the ‘down’ position in WT RT/K-5a2 complex to the ‘up’ position in K103N/Y181C RT/25a complex (Figure 5—figure supplement 1D). This stepwise movement of Tyr183 is likely triggered by three factors: (i) inhibitor repositioning because of the K103N mutation, (ii) loss of aromatic interactions due to the Y181C or Y181I mutation, and (iii) presence of a cyanovinyl group in the inhibitor left wing. Tyr183 makes the most significant contribution to the NNRTI-binding in the circumstance that all of the above factors are present. The ability of K-5a2 and 25a to recruit Tyr183 is particularly significant for their function of inhibiting the polymerase activity, because Tyr183 is completely conserved among all HIV-1 sequences and makes direct contacts with the nucleic acid substrate (Das et al., 2012; Sarafianos et al., 2001). The repositioning of Tyr183 towards the NNIBP removes this important interaction and destabilizes the binding of nucleic acid.
Comparison of the binding modes for K-5a2, 25a and DAPY NNRTIs
By adopting the typical horseshoe conformation, K-5a2 and 25a substantially overlap the binding sites occupied by ETR and RPV (Figure 6). The thiophene substituent in the central ring of K-5a2 and 25a extends further into the entrance channel and is proximal to Glu138 located at the opening. The positions of their left wing structures are adjusted to small conformational changes of Tyr181, Tyr183 and Tyr188 to maximize the contacts with the pocket residues in this region. The piperidine ring in the right wing of K-5a2 and 25a slants away from the cyanophenyl plane in ETR and RPV through a ~ 90° tilt of τ4 (from 16°/10° to –71°/–79°) (Table 5) (Das et al., 2008; Lansdon et al., 2010), leading to the displacement of Tyr318 ~1 Å away from the inhibitors. The main difference in RT conformation is the uplift of the loop proceeding β9 and that connecting β10–β11. Upon the binding of K-5a2 and 25a, Pro225 and Pro236, two residues sitting at the groove channel opening, are pushed apart to accommodate the benzenesulfonamide group, which protrudes to the solvent-exposed surface of the enzyme.
There is no structure available for E138K RT and K101P RT in complex with either ETR or RPV; however, the structures of these two RT mutants in complex with K-5a2 and 25a provide a structural basis to understand why the mutations can render ETR and RPV less effective. In the structure of WT RT/ETR complex, the amino substituent of the central pyrimidine ring forms a salt bridge with the carboxyl side chain of Glu138. Transplanting Lys138 from the structure of E138K RT/K-5a2 complex into this structure reveals a severe charge-charge repulsion between the amino group in ETR and the side chain of Lys138, which would destabilize the binding of the inhibitor (Figure 6—figure supplement 1A). In the case of PRV, Glu138 contributes to the RT-RPV interactions by bonding Lys101 and placing it in the vicinity of the central pyrimidine ring for several critical interactions: (i) hydrophobic interactions between the pyrimidine ring of RPV and the Cε atom of Lys101, (ii) the hydrogen bond between the pyrimidine nitrogen atom and the main-chain nitrogen of Lys101, and (iii) the hydrogen bond between a linker nitrogen atom of RPV and the carbonyl oxygen of Lys101. Modeling Lys138 from the structure of E138K RT/25a complex indicates that Lys at residue 138 of the p51 subunit would shove Lys101 away from RPV. This conformational change in RT would not only disrupt the hydrophobic interactions but also weaken the two key hydrogen bonds due to the deviation of the Lys101 backbone (Figure 6—figure supplement 1B).
The K101P mutation considerably reduces susceptibility to RPV. Superposition of the WT RT/RPV and K101P RT/25a complex structures reveals that the mutation would remove the hydrophobic interaction between the RPV pyrimidine ring and the long aliphatic side chain of Lys101, and possibly introduces steric clashes between the polar groups in the pyrimidine ring and the non-polar side chain of Pro101. More importantly, Lys to Pro substitution will abrogate the critical hydrogen bond between the pyrimidine nitrogen of RPV and the backbone nitrogen in RT due to deprotonated form of Pro101 main-chain nitrogen and hence its inability to act as a hydrogen bond donor (Figure 6—figure supplement 1C).
The improved resistance profiles of 25a over RPV against other RT mutants, especially Y181I RT, Y188L RT and K103N/Y181C RT, are likely due to the bulkier substituents in its right wing and more extensive hydrogen bond interactions with NNIBP residues. Therefore, the π-π interactions provided by Tyr181 and Tyr188 likely make a much less contribution to the binding affinity of 25a than to that of RPV. Moreover, the higher degree of strategic flexibility of 25a (because of more rotatable bonds) can make it more effective to preserve or even enhance its interactions with other NNIBP residues when Y181I and Y188L mutations displace the left wing structure of 25a.
Discussion
Emergence of drug-resistant mutations in HIV-1 RT remains a major challenge for the design and development of NNRTIs. Using ETR as a lead compound, our previous efforts led to the design of two piperidine-substituted thiophene[3,2-d]pyrimidine derivatives, K-5a2 and 25a, with single-digit nanomolar EC50 values against HIV-1 strains containing either WT RT or RT variants bearing various resistance-associated mutations. In the current study, we have shown that 25a is more effective than RPV against a broad set of RT mutants and determined the crystal structures of both WT RT and a number of RT mutants in complex with either K-5a2 or 25a. These high-resolution structures enable unambiguous determination of the binding modes of K-5a2 and 25a, and accurate illustration of the detailed interactions between RT and these highly potent NNRTIs. By virtue of their structural flexibility, K-5a2 and 25a are able to adapt to the conformational changes of RT induced by mutations in the NNIBP and optimize their complementarity with the mutated pocket by varying their multiple torsion angles. As such, the buried areas between the inhibitors and RT are similar across WT RT and various RT mutants, suggesting K-5a2 and 25a can occupy the NNIBP of RT mutants as effectively as they bind to the pocket of WT RT (Table 4). Energy calculation of the K-5a2 and 25a shows that the NNRTIs in different RT-bound conformations are almost isoenergetic (Table 5), indicating the conformational changes of K-5a2 and 25a induced by NNIBP residue mutations do not bring in significant strain energy penalties. Although in cellular environment, both RT and the bound inhibitors are in constant motion, and the interactions between them are undergoing repeated breaking and re-forming, the binding mode captured in the crystal structure should represent the averaged state of the complex or a state with high likelihood.
Aside from structural flexibility, hydrogen bonding with the main chains of NNIBP residues was previously suggested as another method to design NNRTIs that can overcome the effects of drug-resistant mutations in RT (Zhan et al., 2009). Compared with ETR and RPV, K-5a2 and 25a form considerably more hydrogen bonds between their polar groups (the thiophene sulfur, piperidine nitrogen and solvent-exposed sulfonamide) and the main chains of residues throughout the binding pocket. This extensive network of main-chain hydrogen bonds contributes substantially to the free energy of RT-inhibitor binding and is less susceptible to side chain mutations in the pocket. Furthermore, the more extended right wing structures of K-5a2 and 25a contact with a larger set of NNIBP residues than DAPY NNRTIs. Although it potentially makes K-5a2 and 25a susceptible to mutations of the residues not contacted by ETR or RPV, our results have shown that mutations of Pro225 and Pro236, whose side chains interact with the right wing of 25a, but not that of RPV, do not cause resistance to 25a. Additional RT-25a interactions that are not present in the RT/RPV complex include the hydrogen bonds between the sulfonamide group of 25a and Lys104 and Val106. The mutations of these two residues are unlikely to cause loss of potency of 25a because the two hydrogen bonds are established through the main chains of Lys104 and Val106. Even if the side-chain substitutions deviate the main chains, a minor change in the torsion angle τ1, which was shown to span a wide range without significant energetic penalty (Table 5), would readily place the sulfonamide group of 25a in an optimal position for hydrogen bonds formation.
In conclusion, our study depicts the binding poses of two newly developed NNRTIs, compounds K-5a2 and 25a, in their complexes with WT and mutant RTs, and exemplifies how broadly active NNRTIs retain satisfactory activities against RT containing drug-resistant mutations by taking advantage of the plasticity of both the inhibitors and the NNIBP of RT. Our findings provide a reliable model to analyze the structural effects of drug-resistant mutations in RT, and will contribute to structure-based design of novel NNRTIs that can effectively target multiple variants of RT.
Materials and methods
Cloning, protein preparation and crystallization
Request a detailed protocolAn engineered HIV-1 RT construct, RT52A (Bauman et al., 2008; Das et al., 2008), here referred to as WT RT, was used as the template for site-directed mutagenesis to introduce E138K mutation in the p51 subunit, K101P, K103N, Y181I, Y188L, K103N/Y181C, K103N/Y181I and V106A/F227L mutations in the p66 subunit. WT and mutant RTs were expressed and purified as described previously (Das et al., 2008; Frey et al., 2015). Briefly, the p51 subunit with an N-terminal 6xHis tag followed by a human rhinovirus (HRV) 3C protease cleavage site and un-tagged p66 subunit were co-expressed in E. coli BL21 star (DE3) (Thermo Fisher Scientific, Waltham, MA). Cells were grown at 37°C and induced at 17°C for 16 hr. WT and mutant RTs were purified on a HisTrap affinity column and a HiTrap Heparin affinity column (GE Healthcare), sequentially. The N-terminal 6xHis tag was removed by HRV 3C protease, and the un-tagged RT was purified on a Superdex 200 gel filtration column (GE Healthcare) in buffer containing 10 mM Tris (pH 8.0), 75 mM NaCl and 2 mM Tris(2-carboxyethyl)phosphine (TCEP).
Crystallization of WT and mutant RTs were set up using the sitting drop vapor diffusion method at 4°C, with 2 µl of protein solution added to 2 µl of well buffer containing 50 mM MES or imidazole buffer (pH 6.0–6.6), 10% (v/v) polyethylene glycol (PEG) 8000, 100 mM ammonium sulfate, 15 mM magnesium sulfate, and 10 mM spermine. Crystals were grown for 2 weeks, and RT/NNRTI complexes were prepared by soaking RT crystals in buffer containing 0.5 mM K-5a2 or 25a, 50 mM MES or imidazole buffer (pH 6.0), 12% (v/v) polyethylene glycol (PEG) 8000, 100 mM ammonium sulfate, 15 mM magnesium sulfate, 10 mM spermine, 25% ethylene glycol, and 10% DMSO for 2 days. Soaked crystals were harvested and flash-frozen in liquid nitrogen.
Data collection and structure determination
Request a detailed protocolX-ray diffraction data were collected at the Advanced Photon Source at Argonne National Laboratory on beamline 24ID-E at a wavelength of 0.97918 Å. Data sets were integrated and scaled with XDS software package (Kabsch, 2010). Structures of RT/K-5a2 and RT/25a complexes were determined by molecular replacement in Phaser (McCoy et al., 2007) using the structure of WT RT/RPV complex (PDB ID: 4G1Q) as a search template. One RT molecule was present in the asymmetric unit. The ligand restraints and 3D structures of K-5a2 and 25a were generated in eLBOW (Moriarty et al., 2009) using SMILES strings as inputs. Models of K-5a2 and 25a were built into the structures based on the unbiased Fo – Fc difference Fourier electron density map calculated in the absence of an NNRTI. Models were manually rebuilt in Coot (Emsley et al., 2010), and refined in PHENIX (Adams et al., 2010). The quality of the final models was analyzed with MolProbity (Chen et al., 2010). Data collection and refinement statistics are summarized in Supplementary file 1. All figures were generated using PyMOL, UCSF Chimera (Pettersen et al., 2004) or UCSF ChimeraX (Goddard et al., 2018).
Cell lines
Request a detailed protocolMT-4 cells were obtained from the NIH AIDS Reagent Program and authenticated by the supplier. All cells are tested negative for mycoplasma, bacteria, and fungi.
T cell-based anti-HIV-1 activity assays
Request a detailed protocolThe anti-HIV-1 activities of rilpivirine (RPV) against WT HIV-1 (IIIB strain) as well as seven mutant RT-carrying HIV-1 variants (L100I, K103N, E138K, Y181C and K103N/Y181C) were evaluated in MT-4 cells using MTT method as described previously (Kang et al., 2017, 2016; Pannecouque et al., 2008). Briefly, stock solutions (10 × final concentration) of RPV were added in 25 µl to two series of triplicate wells in order to allow simultaneous evaluation of their effects on mock- and HIV-1-infected cells. Using a Biomek 3000 robot (Beckman Instruments, Fullerton, CA), nine five-fold serial dilutions of RPV (final 200 µl volume per well) were made directly in flat-bottomed 96-well microtiter trays, including untreated control HIV-1 and mock-infected cell samples for each sample. Stock of WT HIV-1 or mutant HIV-1 strains (50 µl at 100–300-fold of 50% cell culture infectious dose (CCID50)) or equal amount of culture medium was added to either HIV-1-infected or mock-infected wells of the microtiter tray. Mock-infected cells were used to evaluate the cytotoxicity of the compounds. Exponentially growing MT-4 cells were centrifuged for 5 min at 220 × g and the supernatant was discarded. The MT-4 cells were resuspended at 6 × 105 cells/ml, and 50 µl aliquots were transferred to the microtiter tray wells. At five days after infection, the viability of mock- and HIV-1-infected cells was determined spectrophotometrically in an Infinite M1000 microplate reader (Tecan, Zürich, Switzerland). All data were calculated using the median optical density (OD) value of triplicate wells. The 50% effective antiviral concentration (EC50) was defined as the concentration of the test compound affording 50% protection from viral cytopathogenicity. The 50% cytotoxic concentration (CC50) was defined as the compound concentration that reduced the absorbance (OD540 nm) of mock-infected MT-4 cells by 50%. The results are presented as mean ± SD (n = 3).
Reverse transcriptase inhibition assays
Request a detailed protocolThe HIV-1 RT inhibition assay was performed using a PicoGreen-based EnzChek Reverse Transcriptase Assay kit (Thermo Fisher Scientific) according to manufacturer’s protocol with minor modifications. Briefly, 58 µl of Recombinant WT or mutant RT (final concentration in reaction is 20 nM) in buffer containing 50 mM Tris (pH 8.0), 50 mM KCl, 6 mM MgCl2, and 10 mM DTT was incubated with 2 µl 25a or RPV (Sigma-Aldrich) with a concentration gradient comprising eleven three-fold serial dilutions of each inhibitor, or equal amount of DMSO at 25°C for 1 hr. 30 µl of pre-annealed poly(rA)•d(T)16 in buffer containing 50 mM Tris (pH 8.0), 50 mM KCl, 6 mM MgCl2, 10 mM DTT, and 100 mM dTTP was added to the RT-inhibitor mixture to start DNA polymerization reaction. After 30 min of incubation at 25°C, 10 µl of 150 mM EDTA was added to stop the reaction. 100 µl of 2x PicoGreen reagent was then added to each reaction and product formation was quantified using a TriStar LB 941 microplate reader (Berthold Technologies) with excitation/emission = 485/535 nm. The activity of WT or each mutant RT in the presence of inhibitors was normalized against the DMSO control. The IC50 and curve slope values were calculated by fitting the data into inhibition dose-response curves with variable slopes using GraphPad Prism version 7.0a. The experiment was repeated three times independently.
Accession numbers
Request a detailed protocolThe atomic coordinates and structure factors have been deposited in the Protein Data Bank under the accession codes 6C0J, 6C0K, 6C0L, 6CGF, 6C0N, 6C0O, 6C0P, 6C0R, 6DUF, 6DUG, and 6DUH.
Data availability
Diffraction data and atomic coordinates have been deposited in the Protein Data Bank under the accession codes 6C0J, 6C0K, 6C0L, 6CGF,6C0N, 6C0O, 6C0P, 6C0R, 6DUF, 6DUG, and 6DUH.
-
Crystal structure of HIV-1 reverse transcriptase in complex with nonnucleoside inhibitor K-5a2Publicly available at the RCSB Protein Data Bank (accession no. 6C0J).
-
Crystal structure of HIV-1 K103N mutant reverse transcriptase in complex with non-nucleoside inhibitor K-5a2Publicly available at the RCSB Protein Data Bank (accession no. 6C0K).
-
Crystal structure of HIV-1 E138K mutant reverse transcriptase in complex with non-nucleoside inhibitor K-5a2Publicly available at the RCSB Protein Data Bank (accession no. 6C0L).
-
Crystal structure of HIV-1 Y188L mutant reverse transcriptase in complex with non-nucleoside inhibitor K-5a2Publicly available at the RCSB Protein Data Bank (accession no. 6CGF).
-
Crystal structure of HIV-1 reverse transcriptase in complex with non-nucleoside inhibitor 25aPublicly available at the RCSB Protein Data Bank (accession no. 6C0N).
-
Crystal structure of HIV-1 K103N mutant reverse transcriptase in complex with non-nucleoside inhibitor 25aPublicly available at the RCSB Protein Data Bank (accession no. 6C0O).
-
Crystal structure of HIV-1 E138K mutant reverse transcriptase in complex with non-nucleoside inhibitor 25aPublicly available at the RCSB Protein Data Bank (accession no. 6C0P).
-
Crystal structure of HIV-1 K103N/Y181C mutant reverse transcriptase in complex with non-nucleoside inhibitor 25aPublicly available at the RCSB Protein Data Bank (accession no. 6C0R).
-
Crystal structure of HIV-1 reverse transcriptase V106A/F227L mutant in complex with non-nucleoside inhibitor 25aPublicly available at the RCSB Protein Data Bank (accession no. 6DUF).
-
Crystal structure of HIV-1 reverse transcriptase K101P mutant in complex with non-nucleoside inhibitor 25aPublicly available at the RCSB Protein Data Bank (accession no. 6DUG).
-
Crystal structure of HIV-1 reverse transcriptase Y181I mutant in complex with non-nucleoside inhibitor 25aPublicly available at the RCSB Protein Data Bank (accession no. 6DUH).
-
Crystal structure of HIV-1 reverse transcriptase (RT) in complex with Rilpivirine (TMC278, Edurant), a non-nucleoside rt-inhibiting drugPublicly available at the RCSB Protein Data Bank (accession no. 4G1Q).
-
HIV-1 Reverse Transcriptase in Complex with TMC125Publicly available at the RCSB Protein Data Bank (accession no. 3MEC).
References
-
PHENIX: a comprehensive Python-based system for macromolecular structure solutionActa Crystallographica Section D Biological Crystallography 66:213–221.https://doi.org/10.1107/S0907444909052925
-
TMC278, a Next-Generation nonnucleoside reverse transcriptase inhibitor (NNRTI), Active against Wild-Type and NNRTI-Resistant HIV-1Antimicrobial Agents and Chemotherapy 54:718–727.https://doi.org/10.1128/AAC.00986-09
-
Impact of drug Resistance-Associated amino acid changes in HIV-1 subtype C on susceptibility to newer nonnucleoside reverse transcriptase inhibitorsAntimicrobial Agents and Chemotherapy 59:960–971.https://doi.org/10.1128/AAC.04215-14
-
Crystal engineering of HIV-1 reverse transcriptase for structure-based drug designNucleic Acids Research 36:5083–5092.https://doi.org/10.1093/nar/gkn464
-
MolProbity: all-atom structure validation for macromolecular crystallographyActa Crystallographica Section D Biological Crystallography 66:12–21.https://doi.org/10.1107/S0907444909042073
-
HIV-1 reverse transcriptase and antiviral drug resistance. part 1Current Opinion in Virology 3:111–118.https://doi.org/10.1016/j.coviro.2013.03.012
-
HIV-1 reverse transcriptase and antiviral drug resistance. part 2Current Opinion in Virology 3:119–128.https://doi.org/10.1016/j.coviro.2013.03.014
-
Roles of conformational and positional adaptability in structure-based design of TMC125-R165335 (etravirine) and related non-nucleoside reverse transcriptase inhibitors that are highly potent and effective against wild-type and drug-resistant HIV-1 variantsJournal of Medicinal Chemistry 47:2550–2560.https://doi.org/10.1021/jm030558s
-
HIV-1 reverse transcriptase complex with DNA and nevirapine reveals non-nucleoside inhibition mechanismNature Structural & Molecular Biology 19:253–259.https://doi.org/10.1038/nsmb.2223
-
Structure of HIV-1 RT/TIBO R 86183 complex reveals similarity in the binding of diverse nonnucleoside inhibitorsNature Structural & Molecular Biology 2:407–415.https://doi.org/10.1038/nsb0595-407
-
Features and development of cootActa Crystallographica. Section D, Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
The structural biology of HIV-1: mechanistic and therapeutic insightsNature Reviews Microbiology 10:279–290.https://doi.org/10.1038/nrmicro2747
-
Mechanism of inhibition of HIV-1 reverse transcriptase by non-nucleoside inhibitorsNature Structural & Molecular Biology 2:303–308.https://doi.org/10.1038/nsb0495-303
-
Structure-based evaluation of non-nucleoside inhibitors with improved potency and solubility that target HIV reverse transcriptase variantsJournal of Medicinal Chemistry 58:2737–2745.https://doi.org/10.1021/jm501908a
-
The HIV therapy marketNature Reviews Drug Discovery 15:451–452.https://doi.org/10.1038/nrd.2016.69
-
The Lys103Asn mutation of HIV-1 RT: a novel mechanism of drug resistanceJournal of Molecular Biology 309:437–445.https://doi.org/10.1006/jmbi.2001.4648
-
XDSActa Crystallographica. Section D, Biological Crystallography 66:125–132.https://doi.org/10.1107/S0907444909047337
-
Crystal structures of HIV-1 reverse transcriptase with etravirine (TMC125) and rilpivirine (TMC278): implications for drug designJournal of Medicinal Chemistry 53:4295–4299.https://doi.org/10.1021/jm1002233
-
Evolution of anti-HIV drug candidates. part 3: Diarylpyrimidine (DAPY) analoguesBioorganic & Medicinal Chemistry Letters 11:2235–2239.https://doi.org/10.1016/S0960-894X(01)00412-7
-
Phaser crystallographic softwareJournal of Applied Crystallography 40:658–674.https://doi.org/10.1107/S0021889807021206
-
electronic ligand builder and optimization workbench (eLBOW): a tool for ligand coordinate and restraint generationActa Crystallographica Section D Biological Crystallography 65:1074–1080.https://doi.org/10.1107/S0907444909029436
-
UCSF chimera--a visualization system for exploratory research and analysisJournal of Computational Chemistry 25:1605–1612.https://doi.org/10.1002/jcc.20084
-
Twenty years of therapy for HIV-1 infectionNature Medicine 9:867–873.https://doi.org/10.1038/nm0703-867
-
High resolution structures of HIV-1 RT from four RT–inhibitor complexesNature Structural & Molecular Biology 2:293–302.https://doi.org/10.1038/nsb0495-293
-
Rilpivirine and doravirine have complementary efficacies against NNRTI-Resistant HIV-1 mutantsJAIDS Journal of Acquired Immune Deficiency Syndromes 72:1–491.https://doi.org/10.1097/QAI.0000000000001031
-
Structural basis of HIV-1 resistance to AZT by excisionNature Structural & Molecular Biology 17:1202–1209.https://doi.org/10.1038/nsmb.1908
-
Development of antiretroviral drug resistanceNew England Journal of Medicine 365:637–646.https://doi.org/10.1056/NEJMra1004180
-
2017 update of the drug resistance mutations in HIV-1Topics in Medicine 24:132–133.
-
Effect of mutations at position E138 in HIV-1 reverse transcriptase and their interactions with the M184I mutation on defining patterns of resistance to nonnucleoside reverse transcriptase inhibitors rilpivirine and etravirineAntimicrobial Agents and Chemotherapy 57:3100–3109.https://doi.org/10.1128/AAC.00348-13
-
Design strategies of novel NNRTIs to overcome drug resistanceCurrent Medicinal Chemistry 16:3903–3917.https://doi.org/10.2174/092986709789178019
Article and author information
Author details
Funding
Howard Hughes Medical Institute (Investigator Program)
- Thomas A Steitz
National Institute of General Medical Sciences (GM022778)
- Thomas A Steitz
National Natural Science Foundation of China (81273354)
- Xinyong Liu
Shandong Province (Key Research and Development Project: 2017CXGC1401)
- Xinyong Liu
Shandong Province (Major Project of Science and Technology: 2015ZDJS04001)
- Xinyong Liu
Shandong University (Young Scholars Program: 2016WLJH32)
- Peng Zhan
National Natural Science Foundation of China (Key Project for International Cooperation: 81420108027)
- Xinyong Liu
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Dr. Stefan Sarafianos and Dr. Karen Anderson for the HIV-1 RT52A construct, Dr. Chang Liu and Dr. Bin Liu for advice and critical reading of the manuscript. We thank the staff of the Advanced Photon Source beamline 24-ID-E for technical assistance with data collection and the Richards Center at Yale University for computational support. This work was supported by Howard Hughes Medical Institute and National Institutes of Health (GM022778 to T.A.S.); National Natural Science Foundation of China (81273354 to X.L.); Key Project of National Natural Science Foundation of China for International Cooperation (81420108027 to X.L.); Key research and development project of Shandong Province (2017CXGC1401 to X.L.); Major Project of Science and Technology of Shandong Province (2015ZDJS04001 to X.L.); Young Scholars Program of Shandong University (2016WLJH32 to P.Z.).
Copyright
© 2018, Yang et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,269
- views
-
- 394
- downloads
-
- 60
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Microbiology and Infectious Disease
Flaviviruses, including Zika virus (ZIKV), are a significant global health concern, yet no licensed antivirals exist to treat disease. The small membrane (M) protein plays well-defined roles during viral egress and remains within virion membranes following release and maturation. However, it is unclear whether M plays a functional role in this setting. Here, we show that M forms oligomeric membrane-permeabilising channels in vitro, with increased activity at acidic pH and sensitivity to the prototypic channel-blocker, rimantadine. Accordingly, rimantadine blocked an early stage of ZIKV cell culture infection. Structure-based channel models, comprising hexameric arrangements of two trans-membrane domain protomers were shown to comprise more stable assemblages than other oligomers using molecular dynamics simulations. Models contained a predicted lumenal rimantadine-binding site, as well as a second druggable target region on the membrane-exposed periphery. In silico screening enriched for repurposed drugs/compounds predicted to bind to either one site or the other. Hits displayed superior potency in vitro and in cell culture compared with rimantadine, with efficacy demonstrably linked to virion-resident channels. Finally, rimantadine effectively blocked ZIKV viraemia in preclinical models, supporting that M constitutes a physiologically relevant target. This could be explored by repurposing rimantadine, or development of new M-targeted therapies.
-
- Microbiology and Infectious Disease
Peptidoglycan (PG) is a giant macromolecule that completely surrounds bacterial cells and prevents lysis in hypo-osmotic environments. This net-like macromolecule is made of glycan strands linked to each other by two types of transpeptidases that form either 4→3 (PBPs) or 3→3 (LDTs) cross-links. Previously, we devised a heavy isotope-based PG full labeling method coupled to mass spectrometry to determine the mode of insertion of new subunits into the expanding PG network (Atze et al., 2022). We showed that PG polymerization operates according to different modes for the formation of the septum and of the lateral cell walls, as well as for bacterial growth in the presence or absence of β-lactams in engineered strains that can exclusively rely on LDTs for PG cross-linking when drugs are present. Here, we apply our method to the resolution of the kinetics of the reactions leading to the covalent tethering of the Braun lipoprotein (Lpp) to PG and the subsequent hydrolysis of that same covalent link. We find that Lpp and disaccharide-peptide subunits are independently incorporated into the expanding lateral cell walls. Newly synthesized septum PG appears to contain small amounts of tethered Lpp. LDTs did mediate intense shuffling of Lpp between PG stems leading to a dynamic equilibrium between the PG-tethered and free forms of Lpp.