NMR Evidence for the Displacement of a Conserved Interior Water Molecule in HIV Protease by a Non-Peptide Cyclic Urea-Based Inhibitor

Comprehending the Structure, Dynamics, and Mechanism of Action of Drug-Resistant HIV Protease

Since the emergence of the Human Immunodeficiency Virus (HIV) in the 1980s, strategies to combat HIV-AIDS are continuously evolving. Among the many tested targets to tackle this virus, its protease enzyme (PR) was proven to be an attractive option that brought about numerous research publications and ten FDA-approved drugs to inhibit the PR activity. However, the drug-induced mutations in the enzyme made these small molecule inhibitors ineffective with prolonged usage. The research on HIV PR, therefore, remains a thrust area even today. Through this review, we reiterate the importance of understanding the various structural and functional components of HIV PR in redesigning the structure-based small molecule inhibitors. We also discuss at length the currently available FDA-approved drugs and how these drug molecules induced mutations in the enzyme structure. We then recapitulate the reported mechanisms on how these drug-resistant variants remain sufficiently active to cleave the natural substrates. We end with the future scope covering the recently proposed strategies that show promise to deal with the mutations.

Outliers in SAR and QSAR: 3. Importance of considering the role of water molecules in protein-ligand interactions and quantitative structure-activity relationship studies

It is frequently mentioned that QSARs have not generally lived up to expectations, especially in cases where high predictability is expected yet failed to deliver satisfactory results. Even though outliers can provide an increased opportunity in drug discovery research, outliers in SAR and QSAR can contort predictions and affect the accuracy if proper attention is not given. The percentages of outliers in QSARs have not changed appreciably over the last decade. In our previous studies, we suggested two possible sources of outliers in SAR and QSAR. In this paper, we suggest an additional possible source of outliers in QSAR. We presented several literature examples that show one or more water molecules that play a critical role in protein-ligand binding interactions as observed in their crystal structures. These examples illustrate that failing to account for the effects of water molecules in protein-ligand interactions could mislead interpretation and possibly yield outliers in SAR and QSAR. Examples include cases where QSAR, considering the role of water molecules in protein-ligand crystal structures, provided deeper insight into the understanding and interpretation of the developed QSAR.

Advances in the Treatment of Explicit Water Molecules in Docking and Binding Free Energy Calculations

Background:

The inclusion of direct effects mediated by water during the ligandreceptor recognition is a hot-topic of modern computational chemistry applied to drug discovery and development. Docking or virtual screening with explicit hydration is still debatable, despite the successful cases that have been presented in the last years. Indeed, how to select the water molecules that will be included in the docking process or how the included waters should be treated remain open questions.

Objective:

In this review, we will discuss some of the most recent methods that can be used in computational drug discovery and drug development when the effect of a single water, or of a small network of interacting waters, needs to be explicitly considered.

Results:

Here, we analyse the software to aid the selection, or to predict the position, of water molecules that are going to be explicitly considered in later docking studies. We also present software and protocols able to efficiently treat flexible water molecules during docking, including examples of applications. Finally, we discuss methods based on molecular dynamics simulations that can be used to integrate docking studies or to reliably and efficiently compute binding energies of ligands in presence of interfacial or bridging water molecules.

Conclusions:

Software applications aiding the design of new drugs that exploit water molecules, either as displaceable residues or as bridges to the receptor, are constantly being developed. Although further validation is needed, workflows that explicitly consider water will probably become a standard for computational drug discovery soon.

Filling Blank Spots on the Map: Identification of Ligand Binding Modes and Interacting Water Molecules for Brd4-BD1 by WaterLOGSY Titrations

Fragment-based drug discovery and continuous improvement of existing protein inhibitors rely on the knowledge of exactly how and how strongly a range of small molecules bind to their respective protein targets. By increasing the (perdeuterated) protein concentration, WaterLOGSY titration experiments give access to ligand binding modes even in the case of weak binders as well as to the location of protein-bound water in the surroundings of the ligand. On the basis of these findings, specific chemical modifications of the ligand could be shown to yield significantly enhanced binding affinities.

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Computational challenges of structure-based approaches applied to HIV

Here, we review some of the opportunities and challenges that we face in computational modeling of HIV therapeutic targets and structural biology, both in terms of methodology development and structure-based drug design (SBDD). Computational methods have provided fundamental support to HIV research since the initial structural studies, helping to unravel details of HIV biology. Computational models have proved to be a powerful tool to analyze and understand the impact of mutations and to overcome their structural and functional influence in drug resistance. With the availability of structural data, in silico experiments have been instrumental in exploiting and improving interactions between drugs and viral targets, such as HIV protease, reverse transcriptase, and integrase. Issues such as viral target dynamics and mutational variability, as well as the role of water and estimates of binding free energy in characterizing ligand interactions, are areas of active computational research. Ever-increasing computational resources and theoretical and algorithmic advances have played a significant role in progress to date, and we envision a continually expanding role for computational methods in our understanding of HIV biology and SBDD in the future.

Hydration of protein-RNA recognition sites

We investigate the role of water molecules in 89 protein–RNA complexes taken from the Protein Data Bank. Those with tRNA and single-stranded RNA are less hydrated than with duplex or ribosomal proteins. Protein–RNA interfaces are hydrated less than protein–DNA interfaces, but more than protein–protein interfaces. Majority of the waters at protein–RNA interfaces makes multiple H-bonds; however, a fraction do not make any. Those making H-bonds have preferences for the polar groups of RNA than its partner protein. The spatial distribution of waters makes interfaces with ribosomal proteins and single-stranded RNA relatively ‘dry’ than interfaces with tRNA and duplex RNA. In contrast to protein–DNA interfaces, mainly due to the presence of the 2′OH, the ribose in protein–RNA interfaces is hydrated more than the phosphate or the bases. The minor groove in protein–RNA interfaces is hydrated more than the major groove, while in protein–DNA interfaces it is reverse. The strands make the highest number of water-mediated H-bonds per unit interface area followed by the helices and the non-regular structures. The preserved waters at protein–RNA interfaces make higher number of H-bonds than the other waters. Preserved waters contribute toward the affinity in protein–RNA recognition and should be carefully treated while engineering protein–RNA interfaces.

Accounting for water molecules in drug design

IMPORTANCE OF THE FIELD

Water molecules often appear around ligands in protein crystal structures. Reliable prediction of the effects of water on ligand binding remains a challenge. Solvation effects are crucial for lead optimization where a 100-fold difference in binding affinity is significant but correspond to only ∼3 kcal/mol in binding free energy. Well-known examples, such as nonpeptidic urea inhibitors of HIV protease, prove that careful examination of water molecules and their energetics can contribute significantly to a successful drug design campaign.

AREAS COVERED IN THIS REVIEW

In this review, we examine methods to account for the effect of water in ligand binding at two stages of drug discovery: lead identification via docking calculations and lead optimization. We provide a survey of the models and techniques available to account for water in drug design.

WHAT THE READER WILL GAIN

The reader will become aware of common practices and pitfalls in dealing with water molecules in structure-based drug design.

TAKE HOME MESSAGE

Although solvation effects are not fully understood, some pragmatic recommendations at the end of the article provide guidance for modelers in this area as well as new practitioners.

Water at biomolecular binding interfaces

Water molecules are often found at the binding interface of biomolecular complexes mediating the interaction between polar groups via hydrogen bonds, or simply filling space providing van der Waals interactions. Recent studies have demonstrated the importance of taking such water molecules into account in docking and binding affinity prediction. Here, we review the recent experimental and theoretical work aimed at quantifying the influence of interfacial water on the thermodynamic properties of binding. We highlight especially our recent results obtained by inhomogeneous fluid solvation theory in several systems and the prediction of the thermodynamic consequences of displacement of the bound water molecule by ligand modification. Finally, we discuss possible directions for further progress in this field.

Molecular basis for substrate recognition and drug resistance from 1.1 to 1.6 angstroms resolution crystal structures of HIV-1 protease mutants with substrate analogs

HIV-1 protease (PR) and two drug-resistant variants--PR with the V82A mutation (PR(V82A)) and PR with the I84V mutation (PR(I84V))--were studied using reduced peptide analogs of five natural cleavage sites (CA-p2, p2-NC, p6pol-PR, p1-p6 and NC-p1) to understand the structural and kinetic changes. The common drug-resistant mutations V82A and I84V alter residues forming the substrate-binding site. Eight crystal structures were refined at resolutions of 1.10-1.60 A. Differences in the PR-analog interactions depended on the peptide sequence and were consistent with the relative inhibition. Analog p6(pol)-PR formed more hydrogen bonds of P2 Asn with PR and fewer van der Waals contacts at P1' Pro compared with those formed by CA-p2 or p2-NC in PR complexes. The P3 Gly in p1-p6 provided fewer van der Waals contacts and hydrogen bonds at P2-P3 and more water-mediated interactions. PR(I84V) showed reduced van der Waals interactions with inhibitor compared with PR, which was consistent with kinetic data. The structures suggest that the binding affinity for mutants is modulated by the conformational flexibility of the substrate analogs. The complexes of PR(V82A) showed smaller shifts of the main chain atoms of Ala82 relative to PR, but more movement of the peptide analog, compared to complexes with clinical inhibitors. PR(V82A) was able to compensate for the loss of interaction with inhibitor caused by mutation, in agreement with kinetic data, but substrate analogs have more flexibility than the drugs to accommodate the structural changes caused by mutation. Hence, these structures help to explain how HIV can develop drug resistance while retaining the ability of PR to hydrolyze natural substrates.

Structural biology of HIV

The human immunodeficiency virus (HIV) genome encodes a total of three structural proteins, two envelope proteins, three enzymes, and six accessory proteins. Studies over the past ten years have provided high-resolution three-dimensional structural information for all of the viral enzymes, structural proteins and envelope proteins, as well as for three of the accessory proteins. In some cases it has been possible to solve the structures of the intact, native proteins, but in most cases structural data were obtained for isolated protein domains, peptidic fragments, or mutants. Peptide complexes with two regulatory RNA fragments and a protein complex with an RNA recognition/encapsidation element have also been structurally characterized. This article summarizes the high-resolution structural information that is currently available for HIV proteins and reviews current structure-function and structure-biological relationships.

Molecular recognition of cyclic urea HIV-1 protease inhibitors

As long as the threat of human immunodeficiency virus (HIV) protease drug resistance still exists, there will be a need for more potent antiretroviral agents. We have therefore determined the crystal structures of HIV-1 protease in complex with six cyclic urea inhibitors: XK216, XK263, DMP323, DMP450, XV638, and SD146, in an attempt to identify 1) the key interactions responsible for their high potency and 2) new interactions that might improve their therapeutic benefit. The structures reveal that the preorganized, C2 symmetric scaffolds of the inhibitors are anchored in the active site of the protease by six hydrogen bonds and that their P1 and P2 substituents participate in extensive van der Waals interactions and hydrogen bonds. Because all of our inhibitors possess benzyl groups at P1 and P1', their relative binding affinities are modulated by the extent of their P2 interactions, e.g. XK216, the least potent inhibitor (Ki (inhibition constant) = 4.70 nM), possesses the smallest P2 and the lowest number of P2-S2 interactions; whereas SD146, the most potent inhibitor (Ki = 0.02 nM), contains a benzimidazolylbenzamide at P2 and participates in fourteen hydrogen bonds and approximately 200 van der Waals interactions. This analysis identifies the strongest interactions between the protease and the inhibitors, suggests ways to improve potency by building into the S2 subsite, and reveals how conformational changes and unique features of the viral protease increase the binding affinity of HIV protease inhibitors.

Solvent magnetization artifacts in high-field NMR studies of macromolecular hydration

With the use of high magnetic fields and improved quality factor ratings of the probeheads in modern NMR spectrometers, radiation damping becomes more and more important. In addition, the demagnetizing field effect from protonated solvents gains significance with the increase of the magnetic field strength. During a typical NMR pulse sequence the magnetic fields caused by these effects become time-dependent, which makes the system nonlinear and may, for example, measurably influence the precession frequencies of all nuclei in the sample. Since radiation damping can affect signals that are several kilohertz away from the solvent resonance, the amplitude, phase, and frequency of the desired signals can be disturbed so as to give rise to spectral artifacts. In particular when difference methods are used to obtain the final spectrum, the data sets may be severely deteriorated by such artifacts. This paper investigates effects from the demagnetizing field and from radiation damping with a selection of pulse sequences in use for studies of macromolecular hydration, and strategies are described for the detection and elimination of the ensuing artifacts.

Improved cyclic urea inhibitors of the HIV-1 protease: synthesis, potency, resistance profile, human pharmacokinetics and X-ray crystal structure of DMP 450

BACKGROUND

Effective HIV protease inhibitors must combine potency towards wild-type and mutant variants of HIV with oral bioavailability such that drug levels in relevant tissues continuously exceed that required for inhibition of virus replication. Computer-aided design led to the discovery of cyclic urea inhibitors of the HIV protease. We set out to improve the physical properties and oral bioavailability of these compounds.

RESULTS

We have synthesized DMP 450 (bis-methanesulfonic acid salt), a water-soluble cyclic urea compound and a potent inhibitor of HIV replication in cell culture that also inhibits variants of HIV with single amino acid substitutions in the protease. DMP 450 is highly selective for HIV protease, consistent with displacement of the retrovirus-specific structural water molecule. Single doses of 10 mg kg-1 DMP 450 result in plasma levels in man in excess of that required to inhibit wild-type and several mutant HIVs. A plasmid-based, in vivo assay model suggests that maintenance of plasma levels of DMP 450 near the antiviral IC90 suppresses HIV protease activity in the animal. We did identify mutants that are resistant to DMP 450, however; multiple mutations within the protease gene caused a significant reduction in the antiviral response.

CONCLUSIONS

DMP 450 is a significant advance within the cyclic urea class of HIV protease inhibitors due to its exceptional oral bioavailability. The data presented here suggest that an optimal cyclic urea will provide clinical benefit in treating AIDS if it combines favorable pharmacokinetics with potent activity against not only single mutants of HIV, but also multiply-mutant variants.

Three-dimensional solution structure of the HIV-1 protease complexed with DMP323, a novel cyclic urea-type inhibitor, determined by nuclear magnetic resonance spectroscopy

The three-dimensional solution structure of the HIV-1 protease homodimer, MW 22.2 kDa, complexed to a potent, cyclic urea-based inhibitor, DMP323, is reported. This is the first solution structure of an HIV protease/inhibitor complex that has been elucidated. Multidimensional heteronuclear NMR spectra were used to assemble more than 4,200 distance and angle constraints. Using the constraints, together with a hybrid distance geometry/simulated annealing protocol, an ensemble of 28 NMR structures was calculated having no distance or angle violations greater than 0.3 A or 5 degrees, respectively. Neglecting residues in disordered loops, the RMS deviation (RMSD) for backbone atoms in the family of structures was 0.60 A relative to the average structure. The individual NMR structures had excellent covalent geometry and stereochemistry, as did the restrained minimized average structure. The latter structure is similar to the 1.8-A X-ray structure of the protease/DMP323 complex (Chang CH et al., 1995, Protein Science, submitted); the pairwise backbone RMSD calculated for the two structures is 1.22 A. As expected, the mismatch between the structures is greatest in the loops that are disordered and/or flexible. The flexibility of residues 37-42 and 50-51 may be important in facilitating substrate binding and product release, because these residues make up the respective hinges and tips of the protease flaps. Flexibility of residues 4-8 may play a role in protease regulation by facilitating autolysis.

Flexibility and function in HIV-1 protease

HIV protease is a homodimeric protein whose activity is essential to viral function. We have investigated the molecular dynamics of the HIV protease, thought to be important for proteinase function, bound to high affinity inhibitors using NMR techniques. Analysis of 15N spin relaxation parameters, of all but 13 backbone amide sites, reveals the presence of significant internal motions of the protein backbone. In particular, the flaps that cover the proteins active site of the protein have terminal loops that undergo large amplitude motions on the ps to ns time scale, while the tips of the flaps undergo a conformational exchange on the microsecond time scale. This enforces the idea that the flaps of the proteinase are flexible structures that facilitate function by permitting substrate access to and product release from the active site of the enzyme.

Demonstration of positionally disordered water within a protein hydrophobic cavity by NMR

The presence and location of water of hydration (that is, bound water) in the solution structure of human interleukin-1 beta (hIL-1 beta) was investigated with water-selective two-dimensional heteronuclear magnetic resonance spectroscopy. It is shown here that in addition to water at the surface of the protein and ordered internal water molecules involved in bridging hydrogen bonds, positionally disordered water is present within a large, naturally occurring hydrophobic cavity located at the center of the molecule. These water molecules of hydration have residency times in the range of 1 to 2 nanoseconds to 100 to 200 microseconds and can be readily detected by nuclear magnetic resonance (NMR). Thus, large hydrophobic cavities in proteins may not be truly empty, as analysis of crystal structures appears to show, but may contain mobile water molecules that are crystallographically invisible but detectable by NMR.