Structural analysis of antiviral agents that interact with the capsid of human rhinoviruses

X-Ray diffraction data have been obtained for nine related antiviral agents ("WIN compounds") while bound to human rhinovirus 14 (HRV14). These compounds can inhibit both viral attachment to host cells and uncoating. To calculate interpretable electron density maps it was necessary to account for (1) the low (approximately 60%) occupancies of these compounds in the crystal, (2) the large (up to 7.9 A) conformational changes induced at the attachment site, and (3) the incomplete diffraction data. Application of a density difference map technique, which exploits the 20-fold noncrystallographic redundancy in HRV14, resulted in clear images of the HRV14:WIN complexes. A real-space refinement procedure was used to fit atomic models to these maps. The binding site of WIN compounds in HRV14 is a hydrophobic pocket composed mainly from residues that form the beta-barrel of VP1. Among rhinoviruses, the residues associated with the binding pocket are far more conserved than external residues and are mostly contained within regular secondary structural elements. Molecular dynamics simulations of three HRV14:WIN complexes suggest that portions of the WIN compounds and viral protein near the entrance of the binding pocket are more flexible than portions deeper within the beta-barrel.     

Human rhinovirus 14 complexed with antiviral compound R 61837.

The binding of the antirhinoviral agent R 61837 to human rhinovirus 14 has been examined by X-ray crystallographic methods. The compound R 61837 binds in the same pocket (underneath the canyon floor) as the "WIN" antirhinoviral agents. It does not penetrate as far into the pocket but causes similar conformational changes in the virus capsid. The movement of residues 1217 to 1221 of viral protein 1 (in the "FMDV loop") is more pronounced for R 61837 than for WIN compounds. Although both R 61837 and WIN antiviral agents partially fill the same hydrophobic pocket, atomic binding interactions differ, showing that considerable diversity in the nature of antiviral agents is possible.     

Pocket factors are unlikely to play a major role in the life cycle of human rhinovirus

Human rhinovirus 14 (HRV14) is a member of the rhinovirus genus, which belongs to the picornavirus family, which includes clinically and economically important members, such as poliovirus, foot-and-mouth disease virus, and endomyocarditis virus. Capsid stability plays an important role in the viral infection process, in that it needs to be stable enough to move from cell to cell and yet be able to release its genetic material upon the appropriate environmental cues from the host cell. It has been suggested that certain host cell molecules, "pocket factors," bind to the WIN drug-binding cavity beneath the canyon floor and provide transient stability to a number of the picornaviruses. To directly test this hypothesis, HRV14 was mutated in (V1188M, C1199W, and V1188M/C1199W) and around (S1223G) the drug-binding pocket. Infectivity, limited proteolysis, and matrix-assisted laser desorption ionization analyses indicate that filling the drug-binding pocket with bulky side chains is not deleterious to the viral life cycle and lends some stabilization to the capsid. In contrast, studies with the S1223G mutant suggest that this mutation at least partially overcomes WIN drug-mediated inhibition of cell attachment and capsid breathing. Finally, HRV16, which is inherently more stable than HRV14 in a number of respects, was found to "breathe" only at 37 degrees C and did not tolerate stabilizing mutations in the drug-binding cavity. These results suggest that it is the drug-binding cavity itself and not the putative pocket factor that is crucial for the capsid dynamics, which is, in turn, necessary for infection.     

Structural and virological studies of the stages of virus replication that are affected by antirhinovirus compounds

Pleconaril is a broad-spectrum antirhinovirus and antienterovirus compound that binds into a hydrophobic pocket within viral protein 1, stabilizing the capsid and resulting in the inhibition of cell attachment and RNA uncoating. When crystals of human rhinovirus 16 (HRV16) and HRV14 are incubated with pleconaril, drug occupancy in the binding pocket is lower than when pleconaril is introduced during assembly prior to crystallization. This effect is far more marked in HRV16 than in HRV14 and is more marked with pleconaril than with other compounds. These observations are consistent with virus yield inhibition studies and radiolabeled drug binding studies showing that the antiviral effect of pleconaril against HRV16 is greater on the infectivity of progeny virions than the parent input viruses. These data suggest that drug integration into the binding pocket during assembly, or at some other late stage in virus replication, may contribute to the antiviral activity of capsid binding compounds.     

Molecular dynamics simulations of human rhinovirus and an antiviral compound

The human rhinovirus 14 (HRV14) protomer, with or without the antiviral compound WIN 52084s, was simulated using molecular dynamics and rotational symmetry boundary conditions to model the effect of the entire icosahedral capsid. The protein asymmetrical unit, comprising four capsid proteins (VP1, VP2, VP3, and VP4) and two calcium ions, was solvated both on the exterior and the interior to fill the inside of the capsid. The stability of the simulations of this large system (~800 residues and 6,650 water molecules) is comparable to more conventional globular protein simulations. The influence of the antiviral compound on compressibility and positional fluctuations is reported. The compressibility, estimated from the density fluctuations in the region of the binding pocket, was found to be greater with WIN 52084s bound than without the drug, substantiating previous computations on reduced viral systems. An increase in compressibility correlates with an entropically more favorable system. In contrast to the increase in density fluctuations and compressibility, the positional fluctuations decreased dramatically for the external loops of VP1 and the N-terminus of VP3 when WIN 52084s is bound. Most of these VP1 and VP3 loops are found near the fivefold axis, a region whose mobility was not considered in reduced systems, but can be observed with this simulation of the full viral protomer. Altered loop flexibility is consistent with changes in proteolytic sensitivity observed experimentally. Moreover, decreased flexibility in these intraprotomeric loops is noteworthy since the externalization of VP4, part of VP1, and RNA during the uncoating process is thought to involve areas near the fivefold axis. Both the decrease in positional fluctuations at the fivefold axis and the increase in compressibility near the WIN pocket are discussed in relationship to the antiviral activity of stabilizing the virus against uncoating.     

WIN 52035-Dependent Human Rhinovirus 16: Assembly Deficiency Caused by Mutations near the Canyon Surface

Three drug-dependent mutants of human rhinovirus 16 (HRV16) were characterized by sequence analyses of spontaneous mutant isolates and were genetically reconstructed from a parental cDNA plasmid. These mutants formed plaques in the presence but not in the absence of the selecting antiviral drug, WIN 52035, which binds to the capsid of wild-type virus and inhibits its attachment to the host cell. The drug-dependent phenotype of each mutant was caused by a single amino acid substitution in the VP1 coat protein. The three independent mutations conferring drug dependence are M1103T, T1208A, and V1210A. Single-step growth experiments involving rescue of one of the three mutants (V1210A) by delayed drug addition suggested (i) that the drug dependence lesion is at the stage of virus assembly and (ii) that one or more components of the viral assembly pool decay in the absence of drug. RNA accumulation and infectivity were unaffected by the absence of drug in all three mutants, suggesting that the labile assembly component is coat protein.     

The drug-binding pocket of the human multidrug resistance P-glycoprotein is accessible to the aqueous medium

P-Glycoprotein (P-gp) is an ATP-dependent drug pump that transports a broad range of compounds out of the cell. Cross-linking studies have shown that the drug-binding pocket is at the interface between the transmembrane (TM) domains and can simultaneously bind two different drug substrates. Here, we determined whether cysteine residues within the drug-binding pocket were accessible to the aqueous medium. Cysteine mutants were tested for their reactivity with the charged thiol-reactive compounds sodium (2-sulfonatoethyl)methanethiosulfonate (MTSES) and [2-(trimethylammonium)ethyl)]methanethiosulfonate (MTSET). Residue Ile-306(TM5) is close to the verapamil-binding site. It was changed to cysteine, reacted with MTSES or MTSET, and assayed for verapamil-stimulated ATPase activity. Reaction of mutant I306C(TM5) with either compound reduced its affinity for verapamil. We confirmed that the reduced affinity for verapamil was indeed due to introduction of a charge at position 306 by demonstrating that similar effects were observed when Ile-306 was replaced with arginine or glutamic acid. Mutant I306R showed a 50-fold reduction in affinity for verapamil and very little change in the affinity for rhodamine B or colchicine. MTSES or MTSET modification also affected the cross-linking pattern between pairs of cysteines in the drug-binding pocket. For example, both MTSES and MTSET inhibited cross-linking between I306C(TM5) and I868C(TM10). Inhibition was enhanced by ATP hydrolysis. By contrast, cross-linking of cysteine residues located outside the drug-binding pocket (such as G300C(TM5)/F770C(TM8)) was not affected by MTSES or MTSET. These results indicate that the drug-binding pocket is accessible to water.     

Insight into Pleiotropic Drug Resistance ATP-binding Cassette Pump Drug Transport through Mutagenesis of Cdr1p Transmembrane Domains

The fungal ATP-binding cassette (ABC) transporter Cdr1 protein (Cdr1p), responsible for clinically significant drug resistance, is composed of two transmembrane domains (TMDs) and two nucleotide binding domains (NBDs). We have probed the nature of the drug binding pocket by performing systematic mutagenesis of the primary sequences of the 12 transmembrane segments (TMSs) found in the TMDs. All mutated proteins were expressed equally well and localized properly at the plasma membrane in the heterologous host Saccharomyces cerevisiae, but some variants differed significantly in efflux activity, substrate specificity, and coupled ATPase activity. Replacement of the majority of the amino acid residues with alanine or glycine yielded neutral mutations, but about 42% of the variants lost resistance to drug efflux substrates completely or selectively. A predicted three-dimensional homology model shows that all the TMSs, apart from TMS4 and TMS10, interact directly with the drug-binding cavity in both the open and closed Cdr1p conformations. However, TMS4 and TMS10 mutations can also induce total or selective drug susceptibility. Functional data and homology modeling assisted identification of critical amino acids within a drug-binding cavity that, upon mutation, abolished resistance to all drugs tested singly or in combinations. The open and closed Cdr1p models enabled the identification of amino acid residues that bordered a drug-binding cavity dominated by hydrophobic residues. The disposition of TMD residues with differential effects on drug binding and transport are consistent with a large polyspecific drug binding pocket in this yeast multidrug transporter.     

Distinct effects of two HIV-1 capsid assembly inhibitor families that bind the same site within the N-terminal domain of the viral CA protein

The emergence of resistance to existing classes of antiretroviral drugs necessitates finding new HIV-1 targets for drug discovery. The viral capsid (CA) protein represents one such potential new target. CA is sufficient to form mature HIV-1 capsids in vitro, and extensive structure-function and mutational analyses of CA have shown that the proper assembly, morphology, and stability of the mature capsid core are essential for the infectivity of HIV-1 virions. Here we describe the development of an in vitro capsid assembly assay based on the association of CA-NC subunits on immobilized oligonucleotides. This assay was used to screen a compound library, yielding several different families of compounds that inhibited capsid assembly. Optimization of two chemical series, termed the benzodiazepines (BD) and the benzimidazoles (BM), resulted in compounds with potent antiviral activity against wild-type and drug-resistant HIV-1. Nuclear magnetic resonance (NMR) spectroscopic and X-ray crystallographic analyses showed that both series of inhibitors bound to the N-terminal domain of CA. These inhibitors induce the formation of a pocket that overlaps with the binding site for the previously reported CAP inhibitors but is expanded significantly by these new, more potent CA inhibitors. Virus release and electron microscopic (EM) studies showed that the BD compounds prevented virion release, whereas the BM compounds inhibited the formation of the mature capsid. Passage of virus in the presence of the inhibitors selected for resistance mutations that mapped to highly conserved residues surrounding the inhibitor binding pocket, but also to the C-terminal domain of CA. The resistance mutations selected by the two series differed, consistent with differences in their interactions within the pocket, and most also impaired virus replicative capacity. Resistance mutations had two modes of action, either directly impacting inhibitor binding affinity or apparently increasing the overall stability of the viral capsid without affecting inhibitor binding. These studies demonstrate that CA is a viable antiviral target and demonstrate that inhibitors that bind within the same site on CA can have distinct binding modes and mechanisms of action.     

Simultaneous binding of two different drugs in the binding pocket of the human multidrug resistance P-glycoprotein

The human multidrug resistance P-glycoprotein (P-gp, ABCB1) transports a wide variety of structurally diverse compounds out of the cell. The drug-binding pocket of P-gp is located in the transmembrane domains. Although occupation of the drug-binding pocket by one molecule is sufficient to activate the ATPase activity of P-gp, the drug-binding pocket may be large enough to accommodate two different substrates at the same time. In this study, we used cysteine-scanning mutagenesis to test whether P-gp could simultaneously interact with the thiol-reactive drug substrate, Tris-(2-maleimidoethyl)amine (TMEA) and a second drug substrate. TMEA is a cross-linker substrate of P-gp that allowed us to test for stimulation of cross-linking by a second substrate such as calcein-acetoxymethyl ester, colchicine, demecolcine, cyclosporin A, rhodamine B, progesterone, and verapamil. We report that verapamil induced TMEA cross-linking of mutant F343C(TM6)/V982C(TM12). By contrast, no cross-linked product was detected in mutants F343C(TM6), V982C(TM12), or F343C(TM6)/V982C(TM12) in the presence of TMEA alone. The verapamil-stimulated ATPase activity of mutant F343C(TM6)/V982C(TM12) in the presence of TMEA decreased with increased cross-linking of the mutant protein. These results show that binding of verapamil must induce changes in the drug-binding pocket (induced-fit mechanism) resulting in exposure of residues F343C(TM6)/V982C(TM12) to TMEA. The results also indicate that the common drug-binding pocket in P-gp is large enough to accommodate both verapamil and TMEA simultaneously and suggests that the substrates must occupy different regions in the common drug-binding pocket.     

Three-dimensional structures of drug-resistant mutants of human rhinovirus 14

Mutants of human rhinovirus 14 were isolated and characterized by searching for resistance to compounds that inhibit viral uncoating. The portions of the RNA that code for amino acids that surround the antiviral compound binding site were sequenced. X-ray analysis of two of these mutants, 1188 Val----Leu and 1199 Cys----Tyr, shows that these were single-site substitutions which would sterically hinder drug binding. Differences in the resistance of mutant viruses to various antiviral compounds may be rationalized in terms of the three-dimensional structures of these mutants. Predictions of the structures of mutant rhinovirus 14 with the substitutions 1188 Val----Leu, 1199 Cys----Tyr and 1199 Cys----Trp in VP1 were made using a molecular dynamics technique. The predicted structure of the 1199 Cys----Tyr mutant was consistent with the electron density map, while the 1188 Val----Leu prediction was not. Large (up to 1.4 A) conformational differences between native rhinovirus 14 and the 1199 Cys----Tyr mutant occurred in main-chain atoms near the mutation site. These changes, as well as the orientation of the 1199 tyrosine side-chain, were correctly predicted by the molecular dynamics calculation. The structure of the predicted 1199 Cys----Trp mutation is consistent with the drug-resistant properties of this virus.     

Nonconserved lysine residues attenuate the biological function of the low-risk human papillomavirus E7 protein

The mucosotrophic human papillomaviruses (HPVs) are classified as high-risk (HR) or low-risk (LR) genotypes based on their neoplastic properties. We have demonstrated previously that the E7 protein destabilizes p130, a pRb-related pocket protein, thereby promoting S-phase reentry in postmitotic, differentiated keratinocytes of squamous epithelia, and that HR HPV E7 does so more efficiently than LR HPV E7. The E7 proteins of LR HPV-11 and -6b uniquely possess lysine residues following a casein kinase II phosphorylation motif which is critical for the biological function of E7. We now show that mutations of these lysine residues elevated the efficiency of S-phase reentry, independent of their charge. An 11E7 K39,42R mutation moderately increased the association with and the destabilization of p130. Unexpectedly, polyubiquitination on these lysine residues did not attenuate E7 activity, as their mutation caused elevated proteasomal degradation and decreased protein stability. In this regard, the biologically more potent HR HPV E7 proteins were also less stable than the LR HPV E7 proteins. We infer that these lysine residues impede functional protein-protein interactions. A G22D mutation of 11E7 at the pocket protein binding motif possessed augmented efficiency in promoting S-phase reentry and strongly enhanced association with p130 and pRb. The combined effects of these two classes of 11E7 mutations exhibited an efficiency of S-phase reentry comparable to that of HR HPV E7. Thus, these nonconserved residues are primarily responsible for the differential abilities of LR and HR HPV E7 proteins to promote unscheduled DNA replication in organotypic raft cultures.     

Distribution of drug resistance mutations in type 3 poliovirus identifies three regions involved in uncoating functions.

We have previously described the use of an uncoating inhibitor, WIN 51711, to select drug-resistant mutants of the Sabin strain of poliovirus type 3. Two-thirds of the mutants proved to be dependent on the drug for plaque formation because of extreme thermolability (A. G. Mosser and R. R. Rueckert, J. Virol. 67:1246-1254, 1993). Here we report the responsible mutations; all were traced to single amino acid substitutions. Mutations conferring dependence and thermolability occurred in all four capsid proteins (VP1 to VP4), but all were clustered near residue 53 of VP4 at the inner capsid surface. Amino acid substitutions of the remaining non-drug-dependent mutants were mapped to three distinct loci: (i) on or near the inner capsid surface, at VP4 residue 46 or VP1 residue 129, in the vicinity of the drug dependence substitutions; (ii) at residues 192, 194, and 260 in the lining of the VP1 beta barrel, which is the drug-binding site; and (iii) at VP1 residue 105 on the edge of the canyon surrounding the fivefold axis of symmetry, the putative receptor-binding site. All of the mutations increased the eclipse rate of cell-attached virus. Such mutants help identify parts of the capsid that play a role in viral uncoating functions.     

Viral cell recognition and entry.

Rhinovirus infection is initiated by the recognition of a specific cell-surface receptor. The major group of rhinovirus serotypes attach to intercellular adhesion molecule-1 (ICAM-1). The attachment process initiates a series of conformational changes resulting in the loss of genomic RNA from the virion. X-ray crystallography and sequence comparisons suggested that a deep crevice or canyon is the site on the virus recognized by the cellular receptor molecule. This has now been verified by electron microscopy of human rhinovirus 14 (HRV14) and HRV16 complexed with a soluble component of ICAM-1. A hydrophobic pocket underneath the canyon is the site of binding of various hydrophobic drug compounds that can inhibit attachment and uncoating. This pocket is also associated with an unidentified, possibly cellular in origin, "pocket factor." The pocket factor binding site overlaps the binding site of the receptor. It is suggested that competition between the pocket factor and receptor regulates the conformational changes required for the initiation of the entry of the genomic RNA into the cell.     

Do drug substrates enter the common drug-binding pocket of P-glycoprotein through "gates"?

Overexpression of P-glycoprotein (P-gp; ABCB1) can cause multidrug resistance during cancer and AIDS chemotherapy because of its ability to transport a broad range of structurally unrelated compounds from the cell. P-gp is a member of the ABC family of proteins. It is a single polypeptide containing four domains—two transmembrane (TM) domains each of which contains six TM segments and two nucleotide-binding domains. Chemical modification and cross-linking studies of cysteine mutants of P-gp indicate that the common drug-binding pocket is at the interface between the TM domains. It has been postulated that drug substrates enter the lipid bilayer, are extracted by P-gp and transported to the extracellular medium. It is not clear how drug substrates enter the drug-binding pocket. Here, we propose that drug-substrates diffuse from the lipid bilayer into the drug-binding pocket through “gates” formed by TM segments at either end of the drug-binding pocket.     

Synthesis and Activity of Piperazine-containing Antirhinoviral Agents and Crystal Structure of SDZ 880-061 Bound to Human Rhinovirus 14

A series of antipicornaviral agents containing piperazinyl moieties was synthesized with the objective of obtaining a compound with a broad spectrum of antirhinovirus activity, high potency (≤0.003 μg/ml), and low cytotoxicity (≥30 μg/ml). Five compounds of this series were evaluated in detail for efficacy against various HRV serotypes. The agent SDZ 880-061, containing the benzothiazine moiety SDZ 108-075, which is particularly active against HRV14, and the thiazolyl acetic acid ester group of SDZ 89 – 124, which is potent against HRV1B, indeed has a relatively broad antiviral spectrum. SDZ 880-061 inhibited 85% of 89 HRV serotypes tested at a concentration of ≤3 μg/ml. The 3.0 Å resolution X-ray structure of SDZ 880-061 bound to HRV14 has revealed the binding characteristics of this potent compound. It binds in the same pocket as other capsid-binding antiviral agents characterized to date, leaving the innermost portion of the pocket vacant. The binding causes similar, although less extensive, alterations of the HRV14 VP1 backbone conformation (residues 100 to 110, 151 to 159, and 213 to 224) compared to other antiviral agents analyzed structurally. Although the contacts between SDZ 880-061 and HRV14 are mostly of hydrophobic character, the inhibitor has three relatively short polar interactions with residues of VP1 that represent potential hydrogen bonds. The amount of solvent-accessible surface area of SDZ 880-061 buried in the complex (613 Ų) is within the range of that observed in protein–protein interfaces. The observed influence of time of addition or removal of SDZ 880-061 on virus yield and on the infectious-center formation indicates that the compound primarily interferes with HRV14 cellular attachment. Since it is assumed that uncoating requires virion instability and/or flexibility, the finding that SDZ 880-061 has only a marginal effect on uncoating may be due to the fact that it does not completely fill the hydrophobic pocket.     

Two groups of rhinoviruses revealed by a panel of antiviral compounds present sequence divergence and differential pathogenicity.

A variety of chemically different compounds inhibit the replication of several serotypes of rhinoviruses (common-cold viruses). We noticed that one of these antiviral compounds, WIN 51711, had an antiviral spectrum clearly distinctive from a consensus spectrum or other capsid-binding compounds, although all of them were shown to share the same binding site. A systematic evaluation of all known rhinovirus capsid-binding compounds against all serotyped rhinoviruses was therefore initiated. Multivariate analysis of the results revealed the existence of two groups of rhinoviruses, which we will call antiviral groups A and B. The differential sensitivity of members of these groups to antiviral compounds suggests the existence of a dimorphic binding site. The antiviral groups turned out to be a reflection of a divergence of rhinovirus serotypes on a much broader level. Similarities in antiviral spectra were highly correlated with sequence similarities, not only of amino acids lining the antiviral compound-binding-site, but also of amino acids of the whole VP1 protein. Furthermore, analysis of epidemiological data indicated that group B rhinoviruses produced more than twice as many clinical infections per serotype than group A rhinoviruses did. Rhinoviruses belonging to the minor receptor group were without exception all computed to lie in the same region of antiviral group B.     

Transmembrane segment 1 of human P-glycoprotein contributes to the drug-binding pocket

P-glycoprotein (P-gp; ABCB1) actively transports a broad range of structurally unrelated compounds out of the cell. An important step in the transport cycle is coupling of drug binding with ATP hydrolysis. Drug substrates such as verapamil bind in a common drug-binding pocket at the interface between the TM (transmembrane) domains of P-gp and stimulate ATPase activity. In the present study, we used cysteine-scanning mutagenesis and reaction with an MTS (methanethiosulphonate) thiol-reactive analogue of verapamil (MTS-verapamil) to test whether the first TM segment [TM1 (TM segment 1)] forms part of the drug-binding pocket. One mutant, L65C, showed elevated ATPase activity (10.7-fold higher than an untreated control) after removal of unchanged MTS-verapamil. The elevated ATPase activity was due to covalent attachment of MTS-verapamil to Cys65 because treatment with dithiothreitol returned the ATPase activity to basal levels. Verapamil covalently attached to Cys65 appears to occupy the drug-binding pocket because verapamil protected mutant L65C from modification by MTS-verapamil. The ATPase activity of the MTS-verapamil-modified mutant L65C could not be further stimulated with verapamil, calcein acetoxymethyl ester or demecolcine. The ATPase activity could be inhibited by cyclosporin A but not by trans-(E)-flupentixol. These results suggest that TM1 contributes to the drug-binding pocket.     

Binding of an antiviral agent to a sensitive and a resistant human rhinovirus. Computer simulation studies with sampling of amino acid side-chain conformations. II. Calculation of free-energy differences by thermodynamic integration.

Thermodynamic-cycle perturbation theory and molecular dynamics simulations were used to calculate the difference in the free energy of binding of the antiviral compound WIN53338 to the wild-type human rhinovirus 14 and to a drug-resistant mutant of the virus in which valine 188 of the viral protein 1 is mutated to leucine. Because of the difficulty of achieving adequate sampling of all of the rotational isomers of amino acid side-chains in molecular dynamics simulations, an explicit treatment of the effects of the existence of multiple rotational isomers of residue 188 on the calculated free energies was used. The rotamers of residue 188 were first mapped by steric and energetic techniques as described in the accompanying article. Thermodynamic integration was then carried out during simulations of the virus, both with and without the antiviral compound bound, by mutating residue 188 while restraining its side-chain to one conformation. The contributions of the other rotamers of residue 188 to the free-energy changes for this mutation were then added to those calculated by thermodynamic integration as correction factors. Binding of WIN53338 to the wild-type virus was calculated to be favored over binding to the mutant virus by 1.7(+/- 3.0) kcal/mol. This is consistent with experimental data which, if differences in activity are assumed to be due to differences in binding, indicate that the binding affinity of WIN53338 for the wild-type virus is at least 0.15 to 1.7 kcal/mol greater than for the mutant virus. Thermodynamic integration was also performed in the conventional manner without restraints and was found to give less accurate results.