Three-dimensional structure of aspartyl protease from human immunodeficiency virus HIV-1

Structure-based antivirals for emerging and neglected RNA viruses: an emerging field for medicinal chemistry in academia

"Many [neglected viruses] predominantly hit developing countries in tropical and subtropical regions of the world (40% of the world's population are now at risk of contracting dengue fever), but developed countries are by no means immune to their impact."

Contribution of the 80s loop of HIV-1 protease to the multidrug-resistance mechanism: crystallographic study of MDR769 HIV-1 protease variants

The flexible flaps and the 80s loops (Pro79-Ile84) of HIV-1 protease are crucial in inhibitor binding. Previously, it was reported that the crystal structure of multidrug-resistant 769 (MDR769) HIV-1 protease shows a wide-open conformation of the flaps owing to conformational rigidity acquired by the accumulation of mutations. In the current study, the effect of mutations on the conformation of the 80s loop of MDR769 HIV-1 protease variants is reported. Alternate conformations of Pro81 (proline switch) with a root-mean-square deviation of 3-4.8 Å in the C(α) atoms of the I10V mutant and a side chain with a `flipped-out' conformation in the A82F mutant cause distortion in the S1/S1' binding pockets that affects inhibitor binding. The A82S and A82T mutants show local changes in the electrostatics of inhibitor binding owing to the mutation from nonpolar to polar residues. In summary, the crystallographic studies of four variants of MDR769 HIV-1 protease presented in this article provide new insights towards understanding the drug-resistance mechanism as well as a basis for design of future protease inhibitors with enhanced potency.

Computational modeling of substrate specificity and catalysis of the β-secretase (BACE1) enzyme

In this combined MD simulation and DFT study, interactions of the wild-type (WT) amyloid precursor protein (APP) and its Swedish variant (SW), Lys670 → Asn and Met671 → Leu, with the beta-secretase (BACE1) enzyme and their cleavage mechanisms have been investigated. BACE1 catalyzes the rate-limiting step in the generation of 40-42 amino acid long Alzheimer amyloid beta (Aβ) peptides. All key structural parameters such as position of the flap, volume of the active site, electrostatic binding energy, structures, and positions of the inserts A, D, and F and 10s loop obtained from the MD simulations show that, in comparison to the WT-substrate, BACE1 exhibits greater affinity for the SW-substrate and orients it in a more reactive conformation. The enzyme-substrate models derived from the MD simulations were further utilized to investigate the general acid/base mechanism used by BACE1 to hydrolytically cleave these substrates. This mechanism proceeds through the following two steps: (1) formation of the gem-diol intermediate and (2) cleavage of the peptide bond. For the WT-substrate, the overall barrier of 22.4 kcal/mol for formation of the gem-diol intermediate is 3.3 kcal/mol higher than for the SW-substrate (19.1 kcal/mol). This process is found to be the rate-limiting in the entire mechanism. The computed barrier is in agreement with the measured barrier of ca. 18.00 kcal/mol for the WT-substrate and supports the experimental observation that the cleavage of the SW-substrate is 60 times more efficient than the WT-substrate.

Biomolecularmodeling and simulation: a field coming of age

We assess the progress in biomolecular modeling and simulation, focusing on structure prediction and dynamics, by presenting the field’s history, metrics for its rise in popularity, early expressed expectations, and current significant applications. The increases in computational power combined with improvements in algorithms and force fields have led to considerable success, especially in protein folding, specificity of ligand/biomolecule interactions, and interpretation of complex experimental phenomena (e.g. NMR relaxation, protein-folding kinetics and multiple conformational states) through the generation of structural hypotheses and pathway mechanisms. Although far from a general automated tool, structure prediction is notable for proteins and RNA that preceded the experiment, especially by knowledge-based approaches. Thus, despite early unrealistic expectations and the realization that computer technology alone will not quickly bridge the gap between experimental and theoretical time frames, ongoing improvements to enhance the accuracy and scope of modeling and simulation are propelling the field onto a productive trajectory to become full partner with experiment and a field on its own right.

Effect of structural stress on the flexibility and adaptability of HIV-1 protease

Resistance remains a major issue with regards to HIV-1 protease, despite the availability of numerous HIV-1 protease inhibitors and copious amounts of structural and binding data. In an effort to improve our understanding of how HIV-1 protease is able to "outsmart" new drugs, we have investigated the flexibility of HIV-1 protease and in particular how it adapts to different structural stresses. Our analysis has highlighted the effects of space group on the variability between structures of HIV-1 protease and suggests that consideration of multiple structures and appropriate consideration of different conformations of the Ile50 residue is necessary in any structural analysis. Calculation of the root-mean-square deviation on a per-residue basis has been used to identify 'natural variation', while mutational and ligand analyses have been carried out to identify the effect on structure as a result of specific stresses. It was observed that mutations readily cause changes to occur at sites both close to and distant from a mutation site, with changes more likely to occur at residues that are sites of other major mutations. It is also revealed that HIV-1 protease adaption is dependent on the type and the structure of any bound ligand. Identification of the specific changes that occur due to these stresses will aid in the understanding of resistance and also aid in the design of new drugs.

Molecular Basis for Drug Resistance in HIV-1 Protease

HIV-1 protease is one of the major antiviral targets in the treatment of patients infected with HIV-1. The nine FDA approved HIV-1 protease inhibitors were developed with extensive use of structure-based drug design, thus the atomic details of how the inhibitors bind are well characterized. From this structural understanding the molecular basis for drug resistance in HIV-1 protease can be elucidated. Selected mutations in response to therapy and diversity between clades in HIV-1 protease have altered the shape of the active site, potentially altered the dynamics and even altered the sequence of the cleavage sites in the Gag polyprotein. All of these interdependent changes act in synergy to confer drug resistance while simultaneously maintaining the fitness of the virus. New strategies, such as incorporation of the substrate envelope constraint to design robust inhibitors that incorporate details of HIV-1 protease’s function and decrease the probability of drug resistance, are necessary to continue to effectively target this key protein in HIV-1 life cycle.

Substrate recognition in HIV-1 protease: a computational study

HIV-1 protease is a crucial enzyme for the life cycle of the human immunodeficiency virus, the retrovirus that triggers AIDS. It is well documented that HIV-1 protease mediates the cleavage of Gag, Gag-Pol, and Nef precursor polyproteins and is highly selective concerning the set of 12 different amino acid sequences that cleaves. However, the governing principles and physical parameters, which determine substrate recognition and specificity, remain poorly understood despite the many speculative proposals that abound in the literature. In fact, it has been difficult so far to circumvent the fact that protease's substrates share little sequence identity and lack an obvious consensus binding motif. We have used microsecond time scale MD simulations to quantitatively show that some sequences of the polyprotein Gag-Pol that are not cleaved (nonsubstrates) have in fact a higher affinity to the active site of HIV-1 protease than a substrate; i.e., recognition is not governed by affinity to the active site. On the basis of a detailed analysis of the results and experimental data, we propose that the recognition is based on the geometric specificity of PR:Gag and PR:Gag-Pol multiprotein complex, that selects which residues lie in the specific position that makes them accessible to the active site for cleavage.

Aryloxy phosphoramidate triesters: a technology for delivering monophosphorylated nucleosides and sugars into cells

Prodrug technologies aimed at delivering nucleoside monophosphates into cells (protides) have proved to be effective in improving the therapeutic potential of antiviral and anticancer nucleosides. In these cases, the nucleoside monophosphates are delivered into the cell, where they may then be further converted (phosphorylated) to their active species. Herein, we describe one of these technologies developed in our laboratories, known as the phosphoramidate protide method. In this approach, the charges of the phosphate group are fully masked to provide efficient passive cell-membrane penetration. Upon entering the cell, the masking groups are enzymatically cleaved to release the phosphorylated biomolecule. The application of this technology to various therapeutic nucleosides has resulted in improved antiviral and anticancer activities, and in some cases it has transformed inactive nucleosides to active ones. Additionally, the phosphoramidate technology has also been applied to numerous antiviral nucleoside phosphonates, and has resulted in at least three phosphoramidate-based nucleotides progressing to clinical investigations. Furthermore, the phosphoramidate technology has been recently applied to sugars (mainly glucosamine) in order to improve their therapeutic potential. The development of the phosphoramidate technology, mechanism of action and the application of the technology to various monophosphorylated nucleosides and sugars will be reviewed.

The early years of retroviral protease crystal structures

Soon after its discovery, the attempts to develop anti-AIDS therapeutics focused on the retroviral protease (PR)-an enzyme used by lentiviruses to process the precursor polypeptide into mature viral proteins. An urgent need for the three-dimensional structure of PR to guide rational drug design prompted efforts to produce milligram quantities of this enzyme. However, only minute amounts of PR were present in the HIV-1 and HIV-2 viruses, and initial attempts to express this protein in bacteria were not successful. This review describes X-ray crystallographic studies of the retroviral proteases carried out at NCI-Frederick in the late 1980s and early 1990s and puts into perspective the crucial role that the total protein chemical synthesis played in unraveling the structure, mechanism of action, and inhibition of HIV-1 PR. Notably, the first fully correct structure of HIV-1 PR and the first cocrystal structure of its complex with an inhibitor (a substrate-derived, reduced isostere hexapeptide MVT-101) were determined using chemically synthesized protein. Most importantly, these sets of coordinates were made freely available to the research community and were used worldwide to solve X-ray structures of HIV-1 PR complexes with an array of inhibitors and set in motion a variety of theoretical studies. Publication of the structure of chemically synthesized HIV-1 PR complexed with MVT-101 preceded only by six years the approval of the first PR inhibitor as an anti-AIDS drug.

Insights into drug resistance of mutations D30N and I50V to HIV-1 protease inhibitor TMC-114: free energy calculation and molecular dynamic simulation

The single mutations D30N and I50V are considered as the key residue mutations of the HIV-1 protease drug resistance to inhibitors in clinical use. In this work, molecular dynamics (MD) simulations combined with the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method have been performed to investigate the drug-resistant mechanisms of D30N and I50V to an inhibitor TMC-114. The analyses of absolute binding free energies using the separate trajectory approach suggests that the decrease in the van der Waals energy and electrostatic energy in the gas phase results in the drug resistance of D30N to TMC-114, while for I50V, the decrease in the electrostatic energy mainly drive its drug resistance to TMC-114. Detailed binding free energies between TMC-114 and individual protein residues are computed by using a per-residue basis decomposition method, which provides insights into the inhibitor-protein binding mechanism and also explains the drug-resistant mechanisms of mutations D30N and I50V to TMC-114. The study shows that the loss of the hydrogen bond between TMC-114 and the side chain of Asn30' is the main driving force of the resistance of D30N to TMC-114, and in the case of I50V, the increase in the polar solvation energies between TMC-114 and two residues Val50' and Asp30' definitively drives the resistance of I50V to TMC-114. We expect that this work can provide some helpful insights into the nature of mutational effect and aid the future design of better inhibitors.

Web application for studying the free energy of binding and protonation states of protein-ligand complexes based on HINT

A public web server performing computational titration at the active site in a protein-ligand complex has been implemented. This calculation is based on the Hydropathic interaction noncovalent force field. From 3D coordinate data for the protein, ligand and bridging waters (if available), the server predicts the best combination of protonation states for each ionizable residue and/or ligand functional group as well as the Gibbs free energy of binding for the ionization-optimized protein-ligand complex. The 3D structure for the modified molecules is available as output. In addition, a graph depicting how this energy changes with acidity, i.e., as a function of added protons, can be obtained. This data may prove to be of use in preparing models for virtual screening and molecular docking. A few illustrative examples are presented. In beta secretase (2va7) computational titration flipped the amide groups of Gln12 and Asn37 and protonated a ligand amine yielding an improvement of 6.37 kcal mol(-1) in the protein-ligand binding score. Protonation of Glu139 in mutant HIV-1 reverse transcriptase (2opq) allows a water bridge between the protein and inhibitor that increases the protein-ligand interaction score by 0.16 kcal mol(-1). In human sialidase NEU2 complexed with an isobutyl ether mimetic inhibitor (2f11) computational titration suggested that protonating Glu218, deprotonating Arg237, flipping the amide bond on Tyr334, and optimizing the positions of several other polar protons would increase the protein-ligand interaction score by 0.71 kcal mol(-1).

Using specificity to strategically target proteases

Proteases are a family of naturally occurring enzymes in the body whose dysregulation has been implicated in numerous diseases and cancers. Their ability to selectively and catalytically turnover substrate adds both signal amplification and functionality as parameters for the detection of disease. This review will focus on the development of activity-based methodologies to characterize proteases, and in particular, the use of positional scanning, synthetic combinatorial libraries (PS-SCL's), and substrate activity screening (SAS) assays. The use of these approaches to better understand a protease's natural substrate will be discussed as well as the technologies that emerged.

Anti-HIV drugs for cancer therapeutics: back to the future?

The use of anti-HIV drugs as cancer treatments is not new. Azidothymidine was studied as an antineoplastic in the 1990s, but despite promising in vitro data, clinical trials showed little antitumour activity. HIV protease inhibitors were developed in the early 1990s, and their subsequent incorporation into highly active antiretroviral therapy (HAART) has profoundly changed the natural history of HIV infection. The potential antitumour properties of these drugs have been investigated because of their success in treating HIV-related Kaposi's sarcoma. HAART's effects on Kaposi's sarcoma did not always correlate with immune reconstitution, and activity against other solid and haematological malignancies has been established. Inhibition of tumour-cell invasion and angiogenesis were properties first ascribed to inhibition of HIV protease; however, they have pleiotropic antitumour effects, including inhibition of inflammatory cytokine production, proteasome activity, cell proliferation and survival, and induction of apoptosis. HIV protease inhibitors are thus a new class of anticancer drugs with multiple effects, and other anti-HIV drugs might hold similar promise.

Targeting the open-flap conformation of HIV-1 protease with pyrrolidine-based inhibitors

HIV protease is a well-established drug target in antiviral chemotherapy. Immense research efforts have been made to discover effective inhibitors, thus making the enzyme one of the most studied and best characterized proteins. Although the protease exhibits high flexibility, all approved drugs target virtually the same protein conformation. The development of viral cross-resistance demands the generation of inhibitors with novel scaffolds and deviating modes of binding. Herein we report the design and the short, high-yielding stereoselective synthesis of a series of chiral, symmetric pyrrolidine-based inhibitors targeting the open-flap conformation of the protease. The obtained co-crystal structure with one derivative provides a valuable starting point for further inhibitor design.

Peptidomimetic therapeutic agents targeting the protease enzyme of the human immunodeficiency virus and hepatitis C virus

During the past two decades, great strides have been made in the design of peptidomimetic drugs for the treatment of viral infections, despite the stigma of poor drug-like properties, low oral absorption, and high clearance associated with such compounds. This Account summarizes the progress made toward overcoming such liabilities and highlights the drug discovery efforts that have focused specifically on human immunodeficiency virus (HIV) and hepatitis C virus (HCV) protease inhibitors. The arsenal against the incurable disease AIDS, which is caused by HIV infection, includes peptidomimetic compounds that target the virally encoded aspartic protease enzyme. This enzyme is essential to the production of mature HIV particles and plays a key role in maintaining infectivity. However, because of the rapid genomic evolution of viruses, an inevitable consequence in the treatment of all viral infections is the emergence of resistance to the drugs. Therefore, the incomplete suppression of HIV in treatment-experienced AIDS patients will continue to drive the search for more effective therapeutic agents that exhibit efficacy against the mutants raised by the earlier generation of protease inhibitors. Currently, a number of substrate-based peptidomimetic agents that target the virally encoded HCV NS3/4A protease are in clinical development. Mechanistically, these inhibitors can be generally divided into activated carbonyls that are transition-state mimics or compounds that tap into the feedback mode of enzyme-product inhibition. In the HCV field, there is justified optimism that a number of these compounds will soon reach commercialization as therapeutic agents for the treatment of HCV infections. Structural research has guided the successful design of both HIV and HCV protease inhibitors. X-ray crystallography, NMR, and computational studies have provided valuable insight in to the free-state preorganization of peptidomimetic ligands and their enzyme-bound conformation. Researchers have designed a variety of novel bioisosteric replacements of amino acids and short peptides that contain all of the required pharmacophore moieties and play a key role in inducing conformational changes to the overall molecule. The knowledge gained from these studies will undoubtedly guide the future design of therapeutic agents and further contribute to the success of this field.

Active-site mutations in the South african human immunodeficiency virus type 1 subtype C protease have a significant impact on clinical inhibitor binding: kinetic and thermodynamic study

Human immunodeficiency virus (HIV) infections in sub-Saharan Africa represent about 56% of global infections. Study of active-site mutations (the V82A single mutation and the V82F I84V double mutation) in the less-studied South African HIV type 1 subtype C (C-SA) protease indicated that neither mutation had a significant impact on the proteolytic functioning of the protease. However, the binding affinities of, and inhibition by, saquinavir, ritonavir, indinavir, and nelfinavir were weaker for each variant than for the wild-type protease, with the double mutant exhibiting the most dramatic change. Therefore, our results show that the C-SA V82F I84V double mutation decreased the binding affinities of protease inhibitors to levels significantly lower than that required for effective inhibition.

HIV-1 capsid assembly inhibitor (CAI) peptide: structural preferences and delivery into human embryonic lung cells and lymphocytes

The Human immunodeficiency virus 1 derived capsid assembly inhibitor peptide (HIV-1 CAI-peptide) is a promising lead candidate for anti-HIV drug development. Its drawback, however, is that it cannot permeate cells directly. Here we report the transport of the pharmacologically active CAI-peptide into human lymphocytes and Human Embryonic Lung cells (HEL) using the BioShuttle platform. Generally, the transfer of pharmacologically active substances across membranes, demonstrated by confocal laser scanning microscopy (CLSM), could lead to a loss of function by changing the molecule's structure. Molecular dynamics (MD) simulations and circular dichroism (CD) studies suggest that the CAI-peptide has an intrinsic capacity to form a helical structure, which seems to be critical for the pharmacological effect as revealed by intensive docking calculations and comparison with control peptides. This coupling of the CAI-peptide to a BioShuttle-molecule additionally improved its solubility. Under the conditions described, the HIV-1 CAI peptide was transported into living cells and could be localized in the vicinity of the mitochondria.

Overview of protein structural and functional folds

This overview provides an illustrated, comprehensive survey of some commonly observed protein‐fold families and structural motifs, chosen for their functional significance. It opens with descriptions and definitions of the various elements of protein structure and associated terminology. Following is an introduction into web‐based structural bioinformatics that includes surveys of interactive web servers for protein fold or domain annotation, protein‐structure databases, protein‐structure‐classification databases, structural alignments of proteins, and molecular graphics programs available for personal computers. The rest of the overview describes selected families of protein folds in terms of their secondary, tertiary, and quaternary structural arrangements, including ribbon‐diagram examples, tables of representative structures with references, and brief explanations pointing out their respective biological and functional significance.

Combination of non-natural D-amino acid derivatives and allophenylnorstatine-dimethylthioproline scaffold in HIV protease inhibitors have high efficacy in mutant HIV

Several non-natural D-amino acid derivatives were introduced as P2/P3 residues in allophenylnorstatine-containing (Apns; (2S,3S)-3-amino-2-hydroxy-4-phenylbutyric acid) HIV protease inhibitors. The synthetic analogues exhibited potent inhibitory activity against HIV-1 protease enzyme and HIV-1 replication in MT-4 cells. Structure-activity relationships revealed that D-cysteine or serine derivatives contributed to highly potent anti-HIV activities. Interestingly, anti-HIV activity of all the D-amino acid-introduced inhibitors was remarkably enhanced in their anti-HIV activities against a Nelfinavir-resistant clone, which has a D30N mutation in the protease, over that of the wild-type strain. HIV inhibitory activity of several analogues was moderately affected by an inclusion of alpha1-acid glycoprotein in the test medium.

Complexity in modeling and understanding protonation states: computational titration of HIV-1-protease-inhibitor complexes

The computational-titration (CT) algorithm based on the 'natural' Hydropathic INTeractions (HINT) force field is described. The HINT software model is an empirical, non-Newtonian force field derived from experimentally measured partition coefficients for solvent transfer between octanol and H(2)O (log P(o/w)). The CT algorithm allows the identification, modeling, and optimization of multiple protonation states of residues and ligand functional groups at the protein-ligand active site. The importance of taking into account pH and ionization states of residues, which strongly affect the process of ligand binding, for correctly predicting binding free energies is discussed. The application of the CT protocol to a set of six cyclic inhibitors in their complexes with HIV-1 protease is presented, and the advance of HINT as a virtual-screening tool is outlined.