A solution NMR study of the binding kinetics and the internal dynamics of an HIV-1 protease-substrate complex

The structural, dynamic, and thermodynamic basis of darunavir resistance of a heavily mutated HIV-1 protease using molecular dynamics simulation

The human immunodeficiency virus type 1 protease (HIV-1 PR) is an important enzyme in the life cycle of the HIV virus. It cleaves inactive pre-proteins of the virus and changes them into active proteins. Darunavir (DRV) suppresses the wild-type HIV-1 PR (WT-Pr) activity but cannot inhibit some mutant resistant forms (MUT-Pr). Increasing knowledge about the resistance mechanism can be helpful for designing more effective inhibitors. In this study, the mechanism of resistance of a highly MUT-Pr strain against DRV was investigated. For this purpose, complexes of DRV with WT-Pr (WT-Pr-D) and MUT-Pr (MUT-Pr-D) were studied by all-atom molecular dynamics simulation in order to extract the dynamic and energetic properties. Our data revealed that mutations increased the flap-tip flexibility due to the reduction of the flap-flap hydrophobic interactions. So, the protease’s conformation changed from a closed state to a semi-open state that can facilitate the disjunction of DRV from the active site. On the other hand, energy analysis limited to the final basins of the energy landscape indicated that the entropy of binding of DRV to MUT-Pr was more favorable than that of WT-Pr. However, the enthalpy penalty overcomes it and makes binding more unfavorable relative to the WT-Pr. The unfavorable interaction of DRV with R8, I50, I84, D25′, and A28′ residues in MUT-Pr-D relative to WT-Pr-D is the reason for this enthalpy penalty. Thus, mutations drive resistance to DRV. The hydrogen bond analysis showed that compared with WT-Pr, the hydrogen bonds between DRV and the active-site residues of MUT-Pr were disrupted.

Drug Resistance Mechanism of M46I-Mutation-Induced Saquinavir Resistance in HIV-1 Protease Using Molecular Dynamics Simulation and Binding Energy Calculation

Drug-resistance-associated mutation in essential proteins of the viral life cycle is a major concern in anti-retroviral therapy. M46I, a non-active site mutation in HIV-1 protease has been clinically associated with saquinavir resistance in HIV patients. A 100 ns molecular dynamics (MD) simulation and MM-PBSA calculations were performed to study the molecular mechanism of M46I-mutation-based saquinavir resistance. In order to acquire deeper insight into the drug-resistance mechanism, the flap curling, closed/semi-open/open conformations, and active site compactness were studied. The M46I mutation significantly affects the energetics and conformational stability of HIV-1 protease in terms of RMSD, RMSF, Rg, SASA, and hydrogen formation potential. This mutation significantly decreased van der Waals interaction and binding free energy (∆G) in the M46I–saquinavir complex and induced inward flap curling and a wider opening of the flaps for most of the MD simulation period. The predominant open conformation was reduced, but inward flap curling/active site compactness was increased in the presence of saquinavir in M46I HIV-1 protease. In conclusion, the M46I mutation induced structural dynamics changes that weaken the protease grip on saquinavir without distorting the active site of the protein. The produced information may be utilized for the discovery of inhibitor(s) against drug-resistant HIV-1 protease.

Ritonavir and xk263 Binding-Unbinding with HIV-1 Protease: Pathways, Energy and Comparison

Understanding non-covalent biomolecular recognition, which includes drug-protein bound states and their binding/unbinding processes, is of fundamental importance in chemistry, biology, and medicine. Fully revealing the factors that govern the binding/unbinding processes can further assist in designing drugs with desired binding kinetics. HIV protease (HIVp) plays an integral role in the HIV life cycle, so it is a prime target for drug therapy. HIVp has flexible flaps, and the binding pocket can be accessible by a ligand via various pathways. Comparing ligand association and dissociation pathways can help elucidate the ligand-protein interactions such as key residues directly involved in the interaction or specific protein conformations that determine the binding of a ligand under certain pathway(s). Here, we investigated the ligand unbinding process for a slow binder, ritonavir, and a fast binder, xk263, by using unbiased all-atom accelerated molecular dynamics (aMD) simulation with a re-seeding approach and an explicit solvent model. Using ritonavir-HIVp and xk263-HIVp ligand-protein systems as cases, we sampled multiple unbinding pathways for each ligand and observed that the two ligands preferred the same unbinding route. However, ritonavir required a greater HIVp motion to dissociate as compared with xk263, which can leave the binding pocket with little conformational change of HIVp. We also observed that ritonavir unbinding pathways involved residues which are associated with drug resistance and are distal from catalytic site. Analyzing HIVp conformations sampled during both ligand-protein binding and unbinding processes revealed significantly more overlapping HIVp conformations for ritonavir-HIVp rather than xk263-HIVp. However, many HIVp conformations are unique in xk263-HIVp unbinding processes. The findings are consistent with previous findings that xk263 prefers an induced-fit model for binding and unbinding, whereas ritonavir favors a conformation selection model. This study deepens our understanding of the dynamic process of ligand unbinding and provides insights into ligand-protein recognition mechanisms and drug discovery.

Molecular Mechanics Study of Flow and Surface Influence in Ligand-Protein Association

Ligand–protein association is the first and critical step for many biological and chemical processes. This study investigated the molecular association processes under different environments. In biology, cells have different compartments where ligand–protein binding may occur on a membrane. In experiments involving ligand–protein binding, such as the surface plasmon resonance and continuous flow biosynthesis, a substrate flow and surface are required in experimental settings. As compared with a simple binding condition, which includes only the ligand, protein, and solvent, the association rate and processes may be affected by additional ligand transporting forces and other intermolecular interactions between the ligand and environmental objects. We evaluated these environmental factors by using a ligand xk263 binding to HIV protease (HIVp) with atomistic details. Using Brownian dynamics simulations, we modeled xk263 and HIVp association time and probability when a system has xk263 diffusion flux and a non-polar self-assembled monolayer surface. We also examined different protein orientations and accessible surfaces for xk263. To allow xk263 to access to the dimer interface of immobilized HIVp, we simulated the system by placing the protein 20Å above the surface because immobilizing HIVp on a surface prevented xk263 from contacting with the interface. The non-specific interactions increased the binding probability while the association time remained unchanged. When the xk263 diffusion flux increased, the effective xk263 concentration around HIVp, xk263–HIVp association time and binding probability decreased non-linearly regardless of interacting with the self-assembled monolayer surface or not. The work sheds light on the effects of the solvent flow and surface environment on ligand–protein associations and provides a perspective on experimental design.

Increasing the Sampling Efficiency of Protein Conformational Change by Combining a Modified Replica Exchange Molecular Dynamics and Normal Mode Analysis

Understanding conformational change at an atomic level is significant when determining a protein functional mechanism. Replica exchange molecular dynamics (REMD) is a widely used enhanced sampling method to explore protein conformational space. However, REMD with an explicit solvent model requires huge computational resources, immensely limiting its application. In this study, a variation of parallel tempering metadynamics (PTMetaD) with the omission of solvent-solvent interactions in exchange attempts and the use of low-frequency modes calculated by normal-mode analysis (NMA) as collective variables (CVs), namely ossPTMetaD, is proposed with the aim to accelerate MD simulations simultaneously in temperature and geometrical spaces. For testing the performance of ossPTMetaD, five protein systems with diverse biological functions and motion patterns were selected, including large-scale domain motion (AdK), flap movement (HIV-1 protease and BACE1), and DFG-motif flip in kinases (p38α and c-Abl). The simulation results showed that ossPTMetaD requires much fewer numbers of replicas than temperature REMD (T-REMD) with a reduction of ∼70% to achieve a similar exchange ratio. Although it does not obey the detailed balance condition, ossPTMetaD provides consistent results with T-REMD and experimental data. The high accessibility of the large conformational change of protein systems by ossPTMetaD, especially in simulating the very challenging DFG-motif flip of protein kinases, demonstrated its high efficiency and robustness in the characterization of the large-scale protein conformational change pathway and associated free energy profile.

Inhibition of the activity of HIV-1 protease through antibody binding and mutations probed by molecular dynamics simulations

HIV-1 protease is an essential enzyme in the life cycle of the HIV-1 virus. The conformational dynamics of the flap region of the protease is critical for the ligand binding mechanism, as well as for the catalytic activity. The monoclonal antibody F11.2.32 raised against HIV-1 protease inhibits its activity on binding. We have studied the conformational dynamics of protease in its free, inhibitor ritonavir and antibody bound forms using molecular dynamics simulations. We find that upon Ab binding to the epitope region (residues 36-46) of protease, the overall flexibility of the protease is decreased including the flap region and the active site, which is similar to the decrease in flexibility observed by inhibitor binding to the protease. This suggests an allosteric mechanism to inhibit protease activity. Further, the protease mutants G40E and G40R are known to have decreased activity and were also subjected to MD simulations. We find that the loss of flexibility in the mutants is similar to that observed in the protease bound to the Ab/inhibitor. These insights highlight the role played by dynamics in the function of the protease and how control of flexibility through Ab binding and site specific mutations can inhibit protease activity.

Identification of tunnels as in potato hydrolases

Enzymes with an active center hidden in the middle of the molecule in a tunnel-like cavity constitute an interesting object of analysis due to the highly specialized environment for the course of the catalytic reaction. Identifying the tunnel is a challenge in itself. Moreover, the structural conditioning for the course of the reaction provides information on the diversity of the environment, which must necessarily meet the conditions of high specificity. The use of a fuzzy oil drop model to identify residues constituting the walls of the tunnel located in the center of the protein seems highly justified. The fuzzy oil drop model, which assumes the highest concentration of hydrophobicity in the center of the molecule, in these enzymes shows a significant hydrophobicity deficit resulting from the absence of any residues in the central part of the molecule. Comparison of the expected distribution in consistent with the 3D Gaussian distribution where the observed distribution resulting from the interaction of residues in the protein shows significant differences precisely in the positions of residues located near the center of the molecule. The inside characteristics of the tunnel are the background for the enzymatic reaction. This environment additionally constitutes an external force field, which creates favorable conditions for carrying out the catalytic process. The use of fuzzy oil drop model has been verified using the potato (solanum tuberosum) epoxide hydrolase I. This forms the preliminary basis for testing the fuzzy oil drop model. The data presented here provides an impetus for a large scale analysis of all proteins containing tunnels in enzyme structures available in the Protein Data Bank (PDB).

The Impact of HIV-1 Drug Escape on the Global Treatment Landscape

The rising prevalence of HIV drug resistance (HIVDR) could threaten gains made in combating the HIV epidemic and compromise the 90-90-90 target proposed by United Nations Programme on HIV/AIDS (UNAIDS) to have achieved virological suppression in 90% of all persons receiving antiretroviral therapy (ART) by the year 2020. HIVDR has implications for the persistence of HIV, the selection of current and future ART drug regimens, and strategies of vaccine and cure development. Focusing on drug classes that are in clinical use, this Review critically summarizes what is known about the mechanisms the virus utilizes to escape drug control. Armed with this knowledge, strategies to limit the expansion of HIVDR are proposed.

Enzymatic activity of human immunodeficiency virus type 1 protease in crowded solutions

Cells are crowded with various macromolecules and metabolites, which affect biochemical reactions in many ways, from the diffusion of substrates to catalytic activities of enzymes. We herein investigated the proteolytic activity of the human immunodeficiency virus type 1 protease (HIV-1 PR) under non-crowded and crowded conditions. The latter environment was mimicked with various (poly)ethylene glycol molecules as crowding agents. We found that these crowding agents affect the kinetic parameters of the HIV-1 PR catalyzed reaction by increasing the Michaelis-Menten constant and decreasing the maximum velocity. The influence of crowding was concentration dependent. We explain this effect by the dynamics of the HIV-1 PR flexible flaps that cover the peptide substrate binding site and are crucial for enzyme activity, and by a possibly slower substrate-enzyme association time in the crowded conditions.

Drug Resistance Mutation L76V Alters Nonpolar Interactions at the Flap-Core Interface of HIV-1 Protease

Four HIV-1 protease (PR) inhibitors, clinical inhibitors lopinavir and tipranavir, and two investigational compounds 4 and 5, were studied for their effect on the structure and activity of PR with drug-resistant mutation L76V (PR_L76V). Compound 5 exhibited the best K _i value of 1.9 nM for PR_L76V, whereas the other three inhibitors had K _i values of 4.5-7.6 nM, 2-3 orders of magnitude worse than for wild-type enzymes. Crystal structures showed only minor differences in interactions of inhibitors with PR_L76V compared to wild-type complexes. The shorter side chain of Val76 in the mutant lost hydrophobic interactions with Lys45 and Ile47 in the flap, and with Asp30 and Thr74 in the protein core, consistent with decreased stability. Inhibitors forming additional polar interactions with the flaps or dimer interface of PR_L76V were unable to compensate for the decrease in internal hydrophobic contacts. These structures provide insights for inhibitor design.

Sequence, Structural Analysis and Metrics to Define the Unique Dynamic Features of the Flap Regions Among Aspartic Proteases

Aspartic proteases are a class of hydrolytic enzymes that have been implicated in a number of diseases such as HIV, malaria, cancer and Alzheimer's. The flap region of aspartic proteases is a characteristic unique structural feature of these enzymes; and found to have a profound impact on protein overall structure, function and dynamics. Flap dynamics also plays a crucial role in drug binding and drug resistance. Therefore, understanding the structure and dynamic behavior of this flap regions is crucial in the design of potent and selective inhibitors against aspartic proteases. Defining metrics that can describe the flap motion/dynamics has been a challenging topic in literature. This review is the first attempt to compile comprehensive information on sequence, structure, motion and metrics used to assess the dynamics of the flap region of different aspartic proteases in "one pot". We believe that this review would be of critical importance to the researchers from different scientific domains.

Binding kinetics and substrate selectivity in HIV-1 protease-Gag interactions probed at atomic resolution by chemical exchange NMR

The conversion of immature noninfectious HIV-1 particles to infectious virions is dependent upon the sequential cleavage of the precursor group-specific antigen (Gag) polyprotein by HIV-1 protease. The precise mechanism whereby protease recognizes distinct Gag cleavage sites, located in the intrinsically disordered linkers connecting the globular domains of Gag, remains unclear. Here, we probe the dynamics of the interaction of large fragments of Gag and various variants of protease (including a drug resistant construct) using Carr-Purcell-Meiboom-Gill relaxation dispersion and chemical exchange saturation transfer NMR experiments. We show that the conformational dynamics within the flaps of HIV-1 protease that form the lid over the catalytic cleft play a significant role in substrate specificity and ordered Gag processing. Rapid interconversion between closed and open protease flap conformations facilitates the formation of a transient, sparsely populated productive complex between protease and Gag substrates. Flap closure traps the Gag cleavage sites within the catalytic cleft of protease. Modulation of flap opening through protease-Gag interactions fine-tunes the lifetime of the productive complex and hence the likelihood of Gag proteolysis. A productive complex can also be formed in the presence of a noncognate substrate but is short-lived owing to lack of optimal complementarity between the active site cleft of protease and the substrate, resulting in rapid flap opening and substrate release, thereby allowing protease to differentiate between cognate and noncognate substrates.

Protein-ligand (un)binding kinetics as a new paradigm for drug discovery at the crossroad between experiments and modelling

In the last three decades, protein and nucleic acid structure determination and comprehension of the mechanisms, leading to their physiological and pathological functions, have become a cornerstone of biomedical sciences. A deep understanding of the principles governing the fates of cells and tissue at the molecular level has been gained over the years, offering a solid basis for the rational design of drugs aimed at the pharmacological treatment of numerous diseases. Historically, affinity indicators (i.e. Kd and IC50/EC50) have been assumed to be valid indicators of the in vivo efficacy of a drug. However, recent studies pointed out that the kinetics of the drug-receptor binding process could be as important or even more important than affinity in determining the drug efficacy. This eventually led to a growing interest in the characterisation and prediction of the rate constants of protein-ligand association and dissociation. For instance, a drug with a longer residence time can kinetically select a given receptor over another, even if the affinity for both receptors is comparable, thus increasing its therapeutic index. Therefore, understanding the molecular features underlying binding and unbinding processes is of central interest towards the rational control of drug binding kinetics. In this review, we report the theoretical framework behind protein-ligand association and highlight the latest advances in the experimental and computational approaches exploited to investigate the binding kinetics.

Mechanism of the Association Pathways for a Pair of Fast and Slow Binding Ligands of HIV-1 Protease

Equilibrium constants, together with kinetic rate constants of binding, are key factors in the efficacy and safety of drug compounds, informing drug design. However, the association pathways of protein-ligand binding, which contribute to their kinetic behaviors, are little understood. In this work, we used unbiased all-atom molecular dynamics (MD) simulations with an explicit solvent model to study the association processes of protein-ligand binding. Using the HIV protease (HIVp)-xk263 and HIVp-ritonavir protein-ligand systems as cases, we observed that ligand association is a multistep process involving diffusion, localization, and conformational rearrangements of the protein, ligand, and water molecules. Moreover, these two ligands preferred different routes of binding, which reflect two well-known binding mechanisms: induced-fit and conformation selection models. Our study shows that xk263 has a stronger capacity for desolvating surrounding water molecules, thereby inducing a semiopen conformation of the HIVp flaps (induced-fit model). In contrast, the slow dehydration characteristic of ritonavir allows for gradual association with the binding pocket of HIVp when the protein's flap conformation is fully open (conformation selection model). By studying the mechanism of ligand association and understanding the role of solvent molecules during the binding event, we can obtain a different perspective on the mechanism of macromolecule recognition, providing insights into drug discovery.

Structural insights into HIV-1 protease flap opening processes and key intermediates

The study provided an integrated view of the transition pathway of the flap opening of HIV-1 protease using MD simulation.

Enzyme Tunnels and Gates As Relevant Targets in Drug Design

Many enzymes contain tunnels and gates that are essential to their function. Gates reversibly switch between open and closed conformations and thereby control the traffic of small molecules-substrates, products, ions, and solvent molecules-into and out of the enzyme's structure via molecular tunnels. Many transient tunnels and gates undoubtedly remain to be identified, and their functional roles and utility as potential drug targets have received comparatively little attention. Here, we describe a set of general concepts relating to the structural properties, function, and classification of these interesting structural features. In addition, we highlight the potential of enzyme tunnels and gates as targets for the binding of small molecules. The different types of binding that are possible and the potential pharmacological benefits of such targeting are discussed. Twelve examples of ligands bound to the tunnels and/or gates of clinically relevant enzymes are used to illustrate the different binding modes and to explain some new strategies for drug design. Such strategies could potentially help to overcome some of the problems facing medicinal chemists and lead to the discovery of more effective drugs.

Transient HIV-1 Gag-protease interactions revealed by paramagnetic NMR suggest origins of compensatory drug resistance mutations

Cleavage of the group-specific antigen (Gag) polyprotein by HIV-1 protease represents the critical first step in the conversion of immature noninfectious viral particles to mature infectious virions. Selective pressure exerted by HIV-1 protease inhibitors, a mainstay of current anti-HIV-1 therapies, results in the accumulation of drug resistance mutations in both protease and Gag. Surprisingly, a large number of these mutations (known as secondary or compensatory mutations) occur outside the active site of protease or the cleavage sites of Gag (located within intrinsically disordered linkers connecting the globular domains of Gag to one another), suggesting that transient encounter complexes involving the globular domains of Gag may play a role in guiding and facilitating access of the protease to the Gag cleavage sites. Here, using large fragments of Gag, as well as catalytically inactive and active variants of protease, we probe the nature of such rare encounter complexes using intermolecular paramagnetic relaxation enhancement, a highly sensitive technique for detecting sparsely populated states. We show that Gag-protease encounter complexes are primarily mediated by interactions between protease and the globular domains of Gag and that the sites of transient interactions are correlated with surface exposed regions that exhibit a high propensity to mutate in the presence of HIV-1 protease inhibitors.

Pulsed EPR characterization of HIV-1 protease conformational sampling and inhibitor-induced population shifts

The conformational landscape of HIV-1 protease (PR) can be experimentally characterized by pulsed-EPR double electron-electron resonance (DEER). For this characterization, nitroxide spin labels are attached to an engineered cysteine residue in the flap region of HIV-1 PR. DEER distance measurements from spin-labels contained within each flap of the homodimer provide a detailed description of the conformational sampling of apo-enzyme as well as induced conformational shifts as a function of inhibitor binding. The distance distribution profiles are further interpreted in terms of a conformational ensemble scheme that consists of four unique states termed "curled/tucked", "closed", "semi-open" and "wide-open" conformations. Reported here are the DEER results for a drug-resistant variant clinical isolate sequence, V6, in the presence of FDA approved protease inhibitors (PIs) as well as a non-hydrolyzable substrate mimic, CaP2. Results are interpreted in the context of the current understanding of the relationship between conformational sampling, drug resistance, and kinetic efficiency of HIV-1PR as derived from previous DEER and kinetic data for a series of HIV-1PR constructs that contain drug-pressure selected mutations or natural polymorphisms. Specifically, these collective results support the notion that inhibitor-induced closure of the flaps correlates with inhibitor efficiency and drug resistance. This body of work also suggests DEER as a tool for studying conformational sampling in flexible enzymes as it relates to function.

Enhanced conformational sampling technique provides an energy landscape view of large-scale protein conformational transitions

A novel in silico approach (NMA–ITS) is introduced to rapidly and effectively sample the configuration space and give quantitative data for exploring the conformational changes of proteins.

Conformational Motions of HIV-1 Protease Identified Using Reversible Digitally Filtered Molecular Dynamics

HIV-1 protease performs a vital step in the propagation of the HIV virus and is therefore an important drug target in the treatment of AIDS. It consists of a homodimer, with access to the active site limited by two protein flaps. NMR studies have identified two time scales of motions that occur in these flaps, and it is thought that the slower of these is responsible for a conformational change that makes the protein ligand-accessible. This motion occurs on a time scale outside that achievable using traditional molecular dynamics simulations. Reversible Digitally Filtered Molecular Dynamics (RDFMD) is a method that amplifies low frequency motions associated with conformational change and has recently been applied to, among others, E. coli dihydrofolate reductase, inducing a conformational change between known crystal structures. In this paper, the conformational motions of HIV-1 protease produced during MD and RDFMD simulations are presented, including movement between the known semiopen and closed conformations, and the opening and closing of the protein flaps.

Spin labeling and Double Electron-Electron Resonance (DEER) to Deconstruct Conformational Ensembles of HIV Protease

An understanding of macromolecular conformational equilibrium in biological systems is oftentimes essential to understand function, dysfunction, and disease. For the past few years, our lab has been utilizing site-directed spin labeling (SDSL), coupled with electron paramagnetic resonance (EPR) spectroscopy, to characterize the conformational ensemble and ligand-induced conformational shifts of HIV-1 protease (HIV-1PR). The biomedical importance of characterizing the fractional occupancy of states within the conformational ensemble critically impacts our hypothesis of a conformational selection mechanism of drug-resistance evolution in HIV-1PR. The purpose of the following chapter is to give a timeline perspective of our SDSL EPR approach to characterizing conformational sampling of HIV-1PR. We provide detailed instructions for the procedure utilized in analyzing distance profiles for HIV-1PR obtained from pulsed electron-electron double resonance (PELDOR). Specifically, we employ a version of PELDOR known as double electron-electron resonance (DEER). Data are processed with the software package "DeerAnalysis" (http://www.epr.ethz.ch/software), which implements Tikhonov regularization (TKR), to generate a distance profile from electron spin-echo amplitude modulations. We assign meaning to resultant distance profiles based upon a conformational sampling model, which is described herein. The TKR distance profiles are reconstructed with a linear combination of Gaussian functions, which is then statistically analyzed. In general, DEER has proven powerful for observing structural ensembles in proteins and, more recently, nucleic acids. Our goal is to present our advances in order to aid readers in similar applications.

Effects of drug-resistant mutations on the dynamic properties of HIV-1 protease and inhibition by Amprenavir and Darunavir

Molecular dynamics simulations are performed to investigate the dynamic properties of wild-type HIV-1 protease and its two multi-drug-resistant variants (Flap + (L10I/G48V/I54V/V82A) and Act (V82T/I84V)) as well as their binding with APV and DRV inhibitors. The hydrophobic interactions between flap and 80 s (80’s) loop residues (mainly I50-I84’ and I50’-I84) play an important role in maintaining the closed conformation of HIV-1 protease. The double mutation in Act variant weakens the hydrophobic interactions, leading to the transition from closed to semi-open conformation of apo Act. APV or DRV binds with HIV-1 protease via both hydrophobic and hydrogen bonding interactions. The hydrophobic interactions from the inhibitor is aimed to the residues of I50 (I50’), I84 (I84’), and V82 (V82’) which create hydrophobic core clusters to further stabilize the closed conformation of flaps, and the hydrogen bonding interactions are mainly focused with the active site of HIV-1 protease. The combined change in the two kinds of protease-inhibitor interactions is correlated with the observed resistance mutations. The present study sheds light on the microscopic mechanism underlying the mutation effects on the dynamics of HIV-1 protease and the inhibition by APV and DRV, providing useful information to the design of more potent and effective HIV-1 protease inhibitors.

Resistance mechanism of human immunodeficiency virus type-1 protease to inhibitors: A molecular dynamic approach

Human immunodeficiency virus type 1 (HIV-1) protease inhibitors comprise an important class of drugs used in HIV treatments. However, mutations of protease genes accelerated by low fidelity of reverse transcriptase yield drug resistant mutants of reduced affinities for the inhibitors. This problem is considered to be a serious barrier against HIV treatment for the foreseeable future. In this study, molecular dynamic simulation method was used to examine the combinational and additive effects of all known mutations involved in drug resistance against FDA approved inhibitors. Results showed that drug resistant mutations are not randomly distributed along the protease sequence; instead, they are localized on flexible or hot points of the protein chain. Substitution of more hydrophobic residues in flexible points of protease chains tends to increase the folding, lower the flexibility and decrease the active site area of the protease. The reduced affinities of HIV-1 protease for inhibitors seemed to be due to substantial decrease in the size of the active site and flap mobility. A correlation was found between the binding energy of inhibitors and their affinities for each mutant suggesting the distortion of the active site geometry in drug resistance by preventing effective fitting of inhibitors into the enzymes' active site. To overcome the problem of drug resistance of HIV-1 protease, designing inhibitors of variable functional groups and configurations is proposed.

Design, synthesis and evaluation of a potent substrate analog inhibitor identified by scanning Ala/Phe mutagenesis, mimicking substrate co-evolution, against multidrug-resistant HIV-1 protease

Multidrug-resistant (MDR) clinical isolate-769, human immunodeficiency virus type-1 (HIV-1) protease (PDB ID: 1TW7), was shown to exhibit wide-open flaps and an expanded active site cavity, causing loss of contacts with protease inhibitors. In the current study, the expanded active site cavity of MDR769 HIV-1 protease was screened with a series of peptide-inhibitors that were designed to mimic the natural substrate cleavage site, capsid/p2. Scanning Ala/Phe chemical mutagenesis approach was incorporated into the design of the peptide series to mimic the substrate co-evolution. Among the peptides synthesized and evaluated, a lead peptide (6a) with potent activity (IC50: 4.4nM) was identified against the MDR769 HIV-1 protease. Isothermal titration calorimetry data showed favorable binding profile for 6a against both wild type and MDR769 HIV-1 protease variants. Nuclear magnetic resonance spectrum of (15)N-labeled MDR769 HIV-1 protease in complex with 6a showed some major perturbations in chemical shift, supporting the peptide induced conformational changes in protease. Modeling analysis revealed multiple contacts between 6a and MDR769 HIV-1 protease. The lead peptide-inhibitor, 6a, with high potency and good binding profile can be used as the basis for developing potent small molecule inhibitors against MDR variants of HIV.

Resurrection of an Urbilaterian U1A/U2B″/SNF protein

The U1A/U2B″/SNF family of proteins found in the U1 and U2 spliceosomal small nuclear ribonucleoproteins is highly conserved. In spite of the high degree of sequence and structural conservation, modern members of this protein family have unique RNA binding properties. These differences have necessarily resulted from evolutionary processes, and therefore, we reconstructed the protein phylogeny in order to understand how and when divergence occurred and how protein function has been modulated. Contrary to the conventional understanding of an ancient human U1A/U2B″ gene duplication, we show that the last common ancestor of bilaterians contained a single ancestral protein (URB). The gene for URB was synthesized, the protein was overexpressed and purified, and we assessed RNA binding to modern snRNA sequences. We find that URB binds human and Drosophila U1 snRNA SLII and U2 snRNA SLIV with higher affinity than do modern homologs, suggesting that both Drosophila SNF and human U1A/U2B″ have evolved into weaker binders of one RNA or both RNAs.

Elucidating a relationship between conformational sampling and drug resistance in HIV-1 protease

Enzyme targets in rapidly replicating systems, such as retroviruses, commonly respond to drug-selective pressure with mutations arising in the active site pocket that limit inhibitor effectiveness by introducing steric hindrance or by eliminating essential molecular interactions. However, these primary mutations are disposed to compromising pathogenic fitness. Emerging secondary mutations, which are often found outside of the binding cavity, may or can restore fitness while maintaining drug resistance. The accumulated drug pressure selected mutations could have an indirect effect in the development of resistance, such as altering protein flexibility or the dynamics of protein-ligand interactions. Here, we show that accumulation of mutations in a drug-resistant HIV-1 protease (HIV-1 PR) variant, D30N/M36I/A71V, changes the fractional occupancy of the equilibrium conformational sampling ensemble. Correlations are made among populations of the conformational states, namely, closed-like, semiopen, and open-like, with inhibition constants, as well as kinetic parameters. Mutations that stabilize a closed-like conformation correlate with enzymes of lowered activity and with higher affinity for inhibitors, which is corroborated by a further increase in the fractional occupancy of the closed state upon addition of inhibitor or substrate-mimic. Cross-resistance is found to correlate with combinations of mutations that increase the population of the open-like conformations at the expense of the closed-like state while retaining native-like occupancy of the semiopen population. These correlations suggest that at least three states are required in the conformational sampling model to establish the emergence of drug resistance in HIV-1 PR. More importantly, these results shed light on a possible mechanism whereby mutations combine to impart drug resistance while maintaining catalytic activity.

Small molecule regulation of protein conformation by binding in the Flap of HIV protease

The fragment indole-6-carboxylic acid (1F1), previously identified as a flap site binder in a fragment-based screen against HIV protease (PR), has been cocrystallized with pepstatin-inhibited PR and with apo-PR. Another fragment, 3-indolepropionic acid (1F1-N), predicted by AutoDock calculations and confirmed in a novel inhibition of nucleation crystallization assay, exploits the same interactions in the flap site in two crystal structures. Both 1F1 and 1F1-N bind to the closed form of apo-PR and to pepstatin:PR. In solution, 1F1 and 1F1-N raise the Tm of apo-PR by 3.5-5 °C as assayed by differential scanning fluorimetry (DSF) and show equivalent low-micromolar binding constants to both apo-PR and pepstatin:PR, assayed by backscattering interferometry (BSI). The observed signal intensities in BSI are greater for each fragment upon binding to apo-PR than to pepstatin-bound PR, consistent with greater conformational change in the former binding event. Together, these data indicate that fragment binding in the flap site favors a closed conformation of HIV PR.

Inhibitor-induced conformational shifts and ligand-exchange dynamics for HIV-1 protease measured by pulsed EPR and NMR spectroscopy

Double electron-electron resonance (DEER) spectroscopy was utilized to investigate shifts in conformational sampling induced by nine FDA-approved protease inhibitors (PIs) and a nonhydrolyzable substrate mimic for human immunodeficiency virus type 1 protease (HIV-1 PR) subtype B, subtype C, and CRF_01 A/E. The ligand-bound subtype C protease has broader DEER distance profiles, but trends for inhibitor-induced conformational shifts are comparable to those previously reported for subtype B. Ritonavir, one of the strong-binding inhibitors for subtypes B and C, induces less of the closed conformation in CRF_01 A/E. (1)H-(15)N heteronuclear single-quantum coherence (HSQC) spectra were acquired for each protease construct titrated with the same set of inhibitors. NMR (1)H-(15)N HSQC titration data show that inhibitor residence time in the protein binding pocket, inferred from resonance exchange broadening, shifting or splitting correlates with the degree of ligand-induced flap closure measured by DEER spectroscopy. These parallel results show that the ligand-induced conformational shifts resulting from protein-ligand interactions characterized by DEER spectroscopy of HIV-1 PR obtained at the cryogenic temperature are consistent with more physiological solution protein-ligand interactions observed by solution NMR spectroscopy.

Kinetic characterization of the critical step in HIV-1 protease maturation

HIV maturation requires multiple cleavage of long polyprotein chains into functional proteins that include the viral protease itself. Initial cleavage by the protease dimer occurs from within these precursors, and yet only a single protease monomer is embedded in each polyprotein chain. Self-activation has been proposed to start from a partially dimerized protease formed from monomers of different chains binding its own N termini by self-association to the active site, but a complete structural understanding of this critical step in HIV maturation is missing. Here, we captured the critical self-association of immature HIV-1 protease to its extended amino-terminal recognition motif using large-scale molecular dynamics simulations, thus confirming the postulated intramolecular mechanism in atomic detail. We show that self-association to a catalytically viable state requires structural cooperativity of the flexible β-hairpin "flap" regions of the enzyme and that the major transition pathway is first via self-association in the semiopen/open enzyme states, followed by enzyme conformational transition into a catalytically viable closed state. Furthermore, partial N-terminal threading can play a role in self-association, whereas wide opening of the flaps in concert with self-association is not observed. We estimate the association rate constant (k(on)) to be on the order of ∼1 × 10(4) s(-1), suggesting that N-terminal self-association is not the rate-limiting step in the process. The shown mechanism also provides an interesting example of molecular conformational transitions along the association pathway.

NMR study of xenotropic murine leukemia virus-related virus protease in a complex with amprenavir

Xenotropic murine leukemia virus-related virus (XMRV) is a virus created through recombination of two murine leukemia proviruses under artificial conditions during the passage of human prostate cancer cells in athymic nude mice. The homodimeric protease (PR) of XMRV plays a critical role in the production of functional viral proteins and is a prerequisite for viral replication. We synthesized XMRV PR using the wheat germ cell-free expression system and carried out structural analysis of XMRV PR in a complex with an inhibitor, amprenavir (APV), by means of NMR. Five different combinatorially (15)N-labeled samples were prepared and backbone resonance assignments were made by applying Otting's method, with which the amino acid types of the [(1)H, (15)N] HSQC resonances were automatically identified using the five samples (Wu et al., 2006) [14]. A titration experiment involving APV revealed that one APV molecule binds to one XMRV PR dimer. For many residues, two distinct resonances were observed, which is thought to be due to the structural heterogeneity between the two protomers in the APV:XMRV PR=1:2 complex. PR residues at the interface with APV have been identified on the basis of chemical shift perturbation and identification of the intermolecular NOEs by means of filtered NOE experiments. Interestingly, chemical shift heterogeneity between the two protomers of XMRV PR has been observed not only at the interface with APV but also in regions apart from the interface. This indicates that the structural heterogeneity induced by the asymmetry of the binding of APV to the XMRV PR dimer is transmitted to distant regions. This is in contrast to the case of the APV:HIV-1 PR complex, in which the structural heterogeneity is only localized at the interface. Long-range transmission of the structural change identified for the XMRV PR complex might be utilized for the discovery of a new type of drug.

Identifying the molecular mechanics and binding dynamics characteristics of potent inhibitors to HIV-1 protease

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the primary inhibition targets for chemotherapy of AIDS because of its critical role in the replication cycle of the HIV. In this work, a combinatory coarse-grained and atomistic simulation method was developed for dissecting molecular mechanisms and binding process of inhibitors to the active site of HIV-1 PR, in which 35 typical inhibitors were trialed. We found that the molecular size and stiffness of the inhibitors and the binding energy between the inhibitors and PR play important roles in regulating the binding process. Comparatively, the smaller and more flexible inhibitors have larger binding energy and higher binding rates; they even bind into PR without opening the flaps. In contrast, the larger and stiffer inhibitors have lower binding energy and lower binding rate, and their binding is subjected to the opening and gating of the PR flaps. Furthermore, the components of binding free energy were quantified and analyzed by their dependence on the molecular size, structures, and hydrogen bond networks of inhibitors. Our results also deduce significant dynamics descriptors for determining the quantitative structure and property relationship in potent drug ligands for HIV-1 PR inhibition.

Inhibition of HIV-1 protease: the rigidity perspective

MOTIVATION

HIV-1 protease is a key drug target due to its role in the life cycle of the HIV-1 virus. Rigidity analysis using the software First is a computationally inexpensive method for inferring functional information from protein crystal structures. We evaluate the rigidity of 206 high-resolution (2 Å or better) X-ray crystal structures of HIV-1 protease and compare the effects of different inhibitors binding to the enzyme.

RESULTS

Inhibitor binding has little effect on the overall rigidity of the protein homodimer, including the rigidity of the active site. The principal effect of inhibitor binding on rigidity is to constrain the flexibility of the β-hairpin flaps, which move to allow access to the active site of the enzyme. We show that commercially available antiviral drugs which target HIV-1 protease can be divided into two classes, those which significantly affect flap rigidity and those which do not. The non-peptidic inhibitor tipranavir is distinctive in its consistently strong effect on flap rigidity.

CONTACT

jack.heal@warwick.ac.uk; r.roemer@warwick.ac.uk

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Intramolecular regulation of the sequence-specific mRNA interferase activity of MazF fused to a MazE fragment with a linker cleavable by specific proteases

The genomes of human immunodeficiency virus type 1 (HIV-1) and hepatitis C virus (HCV) consist of single-stranded RNA encoding polyproteins, which are processed to individual functional proteins by virus-encoded specific proteases. These proteases have been used as targets for drug development. Here, instead of targeting these proteases to inhibit viral infection, we utilized the protease activity to activate a toxic protein to prevent viral infection. We engineered the MazE-MazF antitoxin-toxin system of Escherichia coli to fuse a C-terminal 41-residue fragment of antitoxin MazE to the N-terminal end of toxin MazF with a linker having a specific protease cleavage site for either HIV PR (HIV-1 protease), NS3 protease (HCV protease), or factor Xa. These fusion proteins formed a stable dimer (instead of the MazF(2)-MazE(2)-MazF(2) heterohexamer in nature) to inactivate the ACA (sequence)-specific mRNA interferase activity of MazF. When the fusion proteins were incubated with the corresponding proteases, the MazE fragment was cleaved from the fusion proteins, releasing active MazF, which then acted as an ACA-specific mRNA interferase cleaving single-stranded MS2 phage RNA. The intramolecular regulation of MazF toxicity by proteases as demonstrated may provide a novel approach for preventive and therapeutic treatments of infection by HIV-1, HCV, and other single-stranded RNA viruses.

Identifying and studying ubiquitin receptors by NMR

Ubiquitin and ubiquitin chains are recognized by a large and growing family of receptor proteins. NMR spectroscopy provides a powerful means to evaluate whether and how a protein binds to ubiquitin. It can be used to measure binding affinities, to map interaction surfaces, and to solve the three-dimensional structure of ubiquitin:receptor complexes. Herein, we describe three NMR techniques of varying complexity that are valuable tools to characterize protein:protein complexes. These include heteronuclear correlation experiments, paramagnetic relaxation enhancement (PRE) experiments via spin labeling, and techniques designed to obtain intermolecular dipole-dipole interactions by nuclear Overhauser effects (NOEs).

The evolution of catalytic function in the HIV-1 protease

The evolution of species is a complex phenomenon based on the optimization of a multidimensional function referred to as fitness. At the level of biomolecular evolution, the fitness function can be reduced to include physiochemical properties relevant to the biological function of a particular molecule. In this work, questions involving the physical-chemical mechanisms underlying the evolution of HIV-1 protease are addressed through molecular simulation and subsequent analysis of thermodynamic properties related to the activity of the enzyme. Specifically, the impact of 40 single amino acid mutations on the binding affinity toward the matrix/capsid (MA/CA) substrate and corresponding transition state intermediate has been characterized using a molecular mechanics Poisson-Boltzmann surface area approach. We demonstrate that this approach is capable of extracting statistically significant information relevant to experimentally determined catalytic activity. Further, no correlation was observed between the effect of mutations on substrate and transition state binding, suggesting independent evolutionary pathways toward optimizing substrate specificity and catalytic activity. In addition, a detailed analysis of calculated binding affinity data suggests that ground-state destabilization (reduced binding affinity for the substrate) could be a contributing factor in the evolutionary optimization of HIV-1 protease. A numerical model is developed to demonstrate that ground-state destabilization is a valid mechanism for activity optimization given the high concentrations of substrate experienced by the functional enzyme in vivo.

Dynamic flaps in HIV-1 protease adopt unique ordering at different stages in the catalytic cycle

The flexibility of HIV-1 protease flaps is known to be essential for the enzymatic activity. Here we attempt to capture a multitude of conformations of the free and substrate-bound HIV-1 protease that differ drastically in their flap arrangements. The substrate binding process suggests the opening of active site gate in conjunction with a reversal of flap tip ordering, from the native semiopen state. The reversed-flap, open-gated enzyme readily transforms to a closed conformation after proper placement of the substrate into the binding cleft. After substrate processing, the closed state protease which possessed opposite flap ordering relative to the semiopen state, encounters another flap reversal via a second open conformation that facilitates the evolution of native semiopen state of correct flap ordering. The complicated transitional pathway, comprising of many high and low energy states, is explored by combining standard and activated molecular dynamics (MD) simulation techniques. The study not only complements the existing findings from X-ray, NMR, EPR, and MD studies but also provides a wealth of detailed information that could help the structure-based drug design process.

Explicit solvent dynamics and energetics of HIV-1 protease flap opening and closing

An accurate description of the conformational dynamics of the β-hairpin flaps of HIV-1 protease is of central importance in elucidating the functional recognition of the enzyme by ligands. Using all-atom molecular dynamics simulations in explicit solvent, with a total of 461 trajectories of ∼50 ns each, we report the closed, semiopen, open, and wide-open flap conformation of the free wild-type protease. The free energy of flap opening and closing from the semiopen state is 0.9 ± 0.2 and 2.4 ± 0.4 kcal/mol, respectively. The mean relaxation time of opening is ∼8 ns, in good agreement with NMR data. The explicit solvent simulations quantitatively confirm the hypothesis that the semiopen state is the dominant population in the free protease whilst fast flap tip fluctuations lead frequently to an open state. More pronounced flap rearrangements lead to a rare wide-open state with the catalytic site completely exposed to the solvent. The structures of the different flap conformations provided herein are of general interest for improved drug design of HIV-1 protease, in particular, the wide-open conformation could be favored by the large Gag and GagPol polyprotein chains. Strategies that take into account multiple flap-gating mechanisms may lead to more effective inhibitors.

Crucial roles of the subnanosecond local dynamics of the flap tips in the global conformational changes of HIV-1 protease

To understand the underlying mechanisms of the open and closed conformational change of HIV-1 protease (HIV-1 PR) at multiple time scales, we performed serial fully unrestrained, extremely long time molecular dynamics simulations with an explicit solvent model. Spontaneous semiopen to closed conformational transition and inhibitor-collision-induced opening of the flaps were simulated in a real time scale. We found that the rapid, local subnanosecond fluctuations of the flap tips might be the mechanisms triggering the global open and close conformational transitions at the 100-ns time scale. The subnanosecond fluctuation is induced by the Phi-Psi rotations of the residues at the flap tips, mainly Psi of Gly49 and Phi of Ile50, disturbing the interactions between the two tips and then their stability. We further showed that the water molecule W301 is helpful for the stability of the PR-inhibitor complex by acting as a collision buffer for the dynamic interaction between flap tips and the inhibitor. These results might help gain a better insight into the dynamics of HIV-1 PR, especially the local dynamics of the flap tips, which may provide important guidelines for design of novel potent inhibitors.

Solution conformation and dynamics of the HIV-1 integrase core domain

The human immunodeficiency virus type 1 (HIV-1) integrase (IN) is a critical enzyme involved in infection. It catalyzes two reactions to integrate the viral cDNA into the host genome, 3' processing and strand transfer, but the dynamic behavior of the active site during catalysis of these two processes remains poorly characterized. NMR spectroscopy can reveal important structural details about enzyme mechanisms, but to date the IN catalytic core domain has proven resistant to such an analysis. Here, we present the first NMR studies of a soluble variant of the catalytic core domain. The NMR chemical shifts are found to corroborate structures observed in crystals, and confirm prior studies suggesting that the alpha4 helix extends toward the active site. We also observe a dramatic improvement in NMR spectra with increasing MgCl(2) concentration. This improvement suggests a structural transition not only near the active site residues but also throughout the entire molecule as IN binds Mg(2+). In particular, the stability of the core domain is linked to the conformation of its C-terminal helix, which has implications for relative domain orientation in the full-length enzyme. (15)N relaxation experiments further show that, although conformationally flexible, the catalytic loop of IN is not fully disordered in the absence of DNA. Indeed, automated chemical shift-based modeling of the active site loop reveals several stable clusters that show striking similarity to a recent crystal structure of prototype foamy virus IN bound to DNA.

Highly conserved glycine 86 and arginine 87 residues contribute differently to the structure and activity of the mature HIV-1 protease

The structural and functional role of conserved residue G86 in HIV-1 protease (PR) was investigated by NMR and crystallographic analyses of substitution mutations of glycine to alanine and serine (PR(G86A) and PR(G86S)). While PR(G86S) had undetectable catalytic activity, PR(G86A) exhibited approximately 6000-fold lower catalytic activity than PR. (1)H-(15)N NMR correlation spectra revealed that PR(G86A) and PR(G86S) are dimeric, exhibiting dimer dissociation constants (K(d)) of approximately 0.5 and approximately 3.2 muM, respectively, which are significantly lower than that seen for PR with R87K mutation (K(d) > 1 mM). Thus, the G86 mutants, despite being partially dimeric under the assay conditions, are defective in catalyzing substrate hydrolysis. NMR spectra revealed no changes in the chemical shifts even in the presence of excess substrate, indicating very poor binding of the substrate. Both NMR chemical shift data and crystal structures of PR(G86A) and PR(G86S) in the presence of active-site inhibitors indicated high structural similarity to previously described PR/inhibitor complexes, except for specific perturbations within the active site loop and around the mutation site. The crystal structures in the presence of the inhibitor showed that the region around residue 86 was connected to the active site by a conserved network of hydrogen bonds, and the two regions moved further apart in the mutants. Overall, in contrast to the role of R87 in contributing significantly to the dimer stability of PR, G86 is likely to play an important role in maintaining the correct geometry of the active site loop in the PR dimer for substrate binding and hydrolysis. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Revealing the dimer dissociation and existence of a folded monomer of the mature HIV-2 protease

Purification and in vitro protein-folding schemes were developed to produce monodisperse samples of the mature wild-type HIV-2 protease (PR2), enabling a comprehensive set of biochemical and biophysical studies to assess the dissociation of the dimeric protease. An E37K substitution in PR2 significantly retards autoproteolytic cleavage during expression. Furthermore, it permits convenient measurement of the dimer dissociation of PR2(E37K) (elevated K(d) approximately 20 nM) by enzyme kinetics. Differential scanning calorimetry reveals a T(m) of 60.5 for PR2 as compared with 65.7 degrees C for HIV-1 protease (PR1). Consistent with weaker binding of the clinical inhibitor darunavir (DRV) to PR2, the T(m) of PR2 increases by 14.8 degrees C in the presence of DRV as compared with 22.4 degrees C for PR1. Dimer interface mutations, such as a T26A substitution in the active site (PR2(T26A)) or a deletion of the C-terminal residues 96-99 (PR2(1-95)), drastically increase the K(d) (>10(5)-fold). PR2(T26A) and PR2(1-95) consist predominantly of folded monomers, as determined by nuclear magnetic resonance (NMR) and size-exclusion chromatography coupled with multiangle light scattering and refractive index measurements (SMR), whereas wild-type PR2 and its active-site mutant PR2(D25N) are folded dimers. Addition of twofold excess active-site inhibitor promotes dimerization of PR2(T26A) but not of PR2(1-95), indicating that subunit interactions involving the C-terminal residues are crucial for dimer formation. Use of SMR and NMR with PR2 facilitates probing for potential inhibitors that restrict protein folding and/or dimerization and, thus, may provide insights for the future design of inhibitors to circumvent drug resistance.

Fragment-based screen against HIV protease

We have employed a fragment-based screen against wild-type (NL4-3) HIV protease (PR) using the Active Sight fragment library and X-ray crystallography. The experiments reveal two new binding sites for small molecules. PR was co-crystallized with fragments, or crystals were soaked in fragment solutions, using five crystal forms, and 378 data sets were collected to 2.3-1.3 A resolution. Fragment binding induces a distinct conformation and specific crystal form of TL-3 inhibited PR during co-crystallization. One fragment, 2-methylcyclohexanol, binds in the 'exo site' adjacent to the Gly(16)Gly(17)Gln(18)loop where the amide of Gly(17)is a specific hydrogen bond donor, and hydrophobic contacts occur with the side chains of Lys(14)and Leu(63). Another fragment, indole-6-carboxylic acid, binds on the 'outside/top of the flap' via hydrophobic contacts with Trp(42), Pro(44), Met(46), and Lys(55), a hydrogen bond with Val(56), and a salt-bridge with Arg(57). 2-acetyl-benzothiophene also binds at this site. This study is the first fragment-based crystallographic screen against HIV PR, and the first time that fragments were screened against an inhibitor-bound drug target to search for compounds that both bind to novel sites and stabilize the inhibited conformation of the target.

Interactions of different inhibitors with active-site aspartyl residues of HIV-1 protease and possible relevance to pepsin

The importance of the active site region aspartyl residues 25 and 29 of the mature HIV-1 protease (PR) for the binding of five clinical and three experimental protease inhibitors [symmetric cyclic urea inhibitor DMP323, nonhydrolyzable substrate analog (RPB) and the generic aspartic protease inhibitor acetyl-pepstatin (Ac-PEP)] was assessed by differential scanning calorimetry. DeltaT(m) values, defined as the difference in T(m) for a given protein in the presence and absence of inhibitor, for PR with DRV, ATV, SQV, RTV, APV, DMP323, RPB, and Ac-PEP are 22.4, 20.8, 19.3, 15.6, 14.3, 14.7, 8.7, and 6.5 degrees C, respectively. Binding of APV and Ac-PEP is most sensitive to the D25N mutation, as shown by DeltaT(m) ratios [DeltaT(m)(PR)/DeltaT(m)(PR(D25N))] of 35.8 and 16.3, respectively, whereas binding of DMP323 and RPB (DeltaT(m) ratios of 1-2) is least affected. Binding of the substrate-like inhibitors RPB and Ac-PEP is nearly abolished (DeltaT(m)(PR)/DeltaT(m)(PR(D29N)) > or = 44) by the D29N mutation, whereas this mutation only moderately affects binding of the smaller inhibitors (DeltaT(m) ratios of 1.4-2.2). Of the nine FDA-approved clinical HIV-1 protease inhibitors screened, APV, RTV, and DRV competitively inhibit porcine pepsin with K(i) values of 0.3, 0.6, and 2.14 microM, respectively. DSC results were consistent with this relatively weak binding of APV (DeltaT(m) 2.7 degrees C) compared with the tight binding of Ac-PEP (DeltaT(m) > or = 17 degrees C). Comparison of superimposed structures of the PR/APV complex with those of PR/Ac-PEP and pepsin/pepstatin A complexes suggests a role for Asp215, Asp32, and Ser219 in pepsin, equivalent to Asp25, Asp25', and Asp29 in PR in the binding and stabilization of the pepsin/APV complex.

NMR studies on the interaction of sugars with the C-terminal domain of an R-type lectin from the earthworm Lumbricus terrestris

The R-type lectin EW29, isolated from the earthworm Lumbricus terrestris, consists of two homologous domains (14,500 Da) showing 27% identity with each other. The C-terminal domain (Ch; C-half) of EW29 (EW29Ch) has two sugar-binding sites in subdomains alpha and gamma, and the protein uses these sugar-binding sites for its function as a single-domain-type hemagglutinin. In order to determine the sugar-binding ability and specificity for each of the two sugar-binding sites in EW29Ch, ligand-induced chemical-shift changes in EW29Ch were monitored using (1)H-(15)N HSQC spectra as a function of increasing concentrations of lactose, melibiose, D-galactose, methyl alpha-D-galactopyranoside and methyl beta-D-galactopyranoside. Shift perturbation patterns for well-resolved resonances confirmed that all of these sugars associated independently with the two sugar-binding sites of EW29Ch. NMR titration experiments showed that the sugar-binding site in subdomain alpha had a slow or intermediate exchange regime on the chemical-shift timescale (K(d) = 10(-2) to 10(-1) mM), whereas that in subdomain gamma had a fast exchange regime for these sugars (K(d) = 2-6 mM). Thus, our results suggest that the two sugar-binding sites of EW29Ch in the same molecule retain its hemagglutinating activity, but this activity is 10-fold lower than that of the whole protein because EW29Ch has two sugar-binding sites in the same molecule, one of which has a weak binding mode.

Solute effects on spin labels at an aqueous-exposed site in the flap region of HIV-1 protease

The effects of solutes on spin-label mobility and protein conformation have been investigated with X-band continuous-wave and pulsed electron paramagnetic resonance (EPR) spectroscopy for spin labels attached to an aqueous-exposed site in the beta-hairpin flap region of HIV-1 protease. Specifically, we examined the effects of glycerol, sucrose, PEG3000, and Ficoll400 for four commonly used nitroxide spin labels and found that the largest perturbations to the EPR line shapes occur for solutions containing PEG3000 and glycerol. From comparisons of the spectral line shapes and distance distribution profiles of spin-labeled HIV-1 protease with and without inhibitor, it was concluded that solutes such as glycerol and PEG3000 alter the line shapes of the spin label in the beta-hairpin flaps of HIV-1 PR by modulation of spin-label mobility through changes in preferential interactions with the solutes. It is noteworthy that the high osmolality of the 40% glycerol solution did not alter the conformation of the flaps as determined from pulsed EPR distance measurements.