Novel binding mode of highly potent HIV-proteinase inhibitors incorporating the (R)-hydroxyethylamine isostere

Unveiling a Novel Solvatomorphism of Anti-inflammatory Flufenamic Acid: X-ray Structure, Quantum Chemical, and In Silico Studies

This paper delves into the polymorphism of 2-[3-(trifluoromethyl)anilino]benzoic acid, commonly referred to as flufenamic acid (FA), a pharmaceutical agent employed in treating inflammatory conditions. The central focus of the study is on a newly unearthed solvatomorphic structure of FA in methanol (FAM), and a thorough comparison is conducted with the commercially available standard structure. Employing a comprehensive approach, including X-ray crystallography, Hirshfeld surface analysis, density functional theory (DFT), molecular docking, and molecular dynamics (MD) simulations, the research aims to unravel the structural and functional implications of solvatomorphism. The X-ray crystal structure analysis brings to light notable differences between the standard FA and solvatomorphic FAM, showcasing variations in intermolecular interactions and crystal packing. Key features such as hydrogen bonding, π···π stacking, and C-H···π interactions are identified as influential factors shaping the stability and conformation of the compounds. Hirshfeld surface analysis further quantifies the nature and contribution of intermolecular interactions, providing a comprehensive perspective on molecular stability. Density functional theory offers valuable electronic structure insights, highlighting disparities in frontier molecular orbitals between FA and FAM. Molecular docking studies against prostaglandin D2 11-ketoreductase explore potential drug interactions, unveiling distinct binding modes and hydrogen bonding patterns that shed light on how the solvatomorphic structure may impact drug-target interactions. In-depth molecular dynamics simulations over 100 ns investigate the stability of the protein-ligand complex, with root mean square deviation and root mean square fluctuation analyses revealing minimal deviations and affirming the stability of FAM within the active site of the target protein.

Target-based drug design strategies to overcome resistance to antiviral agents: opportunities and challenges

Viral infections have a major impact in human health. Ongoing viral transmission and escalating selective pressure have the potential to favor the emergence of vaccine- and antiviral drug-resistant viruses. Target-based approaches for the design of antiviral drugs can play a pivotal role in combating drug-resistant challenges. Drug design computational tools facilitate the discovery of novel drugs. This review provides a comprehensive overview of current drug design strategies employed in the field of antiviral drug resistance, illustrated through the description of a series of successful applications. These strategies include technologies that enhance compound-target affinity while minimizing interactions with mutated binding pockets. Furthermore, emerging approaches such as virtual screening, targeted protein/RNA degradation, and resistance analysis during drug design have been harnessed to curtail the emergence of drug resistance. Additionally, host targeting antiviral drugs offer a promising avenue for circumventing viral mutation. The widespread adoption of these refined drug design strategies will effectively address the prevailing challenge posed by antiviral drug resistance.

Scalable Inhibitors of the Nsp3-Nsp4 Coupling in SARS-CoV-2

The human Betacoronavirus SARS-CoV-2 is a novel pathogen claiming millions of lives and causing a global pandemic that has disrupted international healthcare systems, economies, and communities. The virus is fast mutating and presenting more infectious but less lethal versions. Currently, some small-molecule therapeutics have received FDA emergency use authorization for the treatment of COVID-19, including Lagevrio (molnupiravir) and Paxlovid (nirmaltrevir/ritonavir), which target the RNA-dependent RNA polymerase and the 3CLpro main protease, respectively. Proteins downstream in the viral replication process, specifically the nonstructural proteins (Nsps1-16), are potential drug targets due to their crucial functions. Of these Nsps, Nsp4 is a particularly promising drug target due to its involvement in the SARS-CoV viral replication and double-membrane vesicle formation (mediated via interaction with Nsp3). Given the degree of sequence conservation of these two Nsps across the Betacoronavirus clade, their protein-protein interactions and functions are likely to be conserved as well in SARS-CoV-2. Through AlphaFold2 and its recent advancements, protein structures were generated of Nsp3 and 4 lumenal loops of interest. Then, using a combination of molecular docking suites and an existing library of lead-like compounds, we virtually screened 7 million ligands to identify five putative ligand inhibitors of Nsp4, which could present an alternative pharmaceutical approach against SARS-CoV-2. These ligands exhibit promising lead-like properties (ideal molecular weight and log P profiles), maintain fixed-Nsp4-ligand complexes in molecular dynamics (MD) simulations, and tightly associate with Nsp4 via hydrophobic interactions. Additionally, alternative peptide inhibitors based on Nsp3 were designed and shown in MD simulations to provide a highly stable binding to the Nsp4 protein. Finally, these therapeutics were attached to dendrimer structures to promote their multivalent binding with Nsp4, especially its large flexible luminal loop (Nsp4LLL). The therapeutics tested in this study represent many different approaches for targeting large flexible protein structures, especially those localized to the ER. This study is the first work targeting the membrane rearrangement system of viruses and will serve as a potential avenue for treating viruses with similar replicative function.

Revealing the drug resistance mechanism of saquinavir due to G48V and V82F mutations in subtype CRF01_AE HIV-1 protease: molecular dynamics simulation and binding free energy calculations

Human immunodeficiency virus-1 (HIV-1) protease is one of the important targets in AIDS therapy. The majority of HIV infections are caused due to non-B subtypes in developing countries. The co-occurrence of mutations along with naturally occurring polymorphisms in HIV-1 protease cause resistance to the FDA approved drugs, thereby posing a major challenge in the treatment of antiretroviral therapy. In this work, the resistance mechanism against SQV due to active site mutations G48V and V82F in CRF01_AE (AE) protease was explored. The binding free energy calculations showed that the direct substitution of valine at position 48 introduces a bulkier side chain, directly impairing the interaction with SQV in the binding pocket. Also, the intramolecular hydrogen bonding network of the neighboring residues is altered, indirectly affecting the binding of SQV. Interestingly, the substitution of phenylalanine at position 82 induces conformational changes in the 80's loop and the flap region, thereby favoring the binding of SQV. The V82F mutant structure also maintains similar intramolecular hydrogen bond interactions as observed in AE-WT.Communicated by Ramaswamy H. Sarma.

In silico methods and tools for drug discovery

In the past, conventional drug discovery strategies have been successfully employed to develop new drugs, but the process from lead identification to clinical trials takes more than 12 years and costs approximately $1.8 billion USD on average. Recently, in silico approaches have been attracting considerable interest because of their potential to accelerate drug discovery in terms of time, labor, and costs. Many new drug compounds have been successfully developed using computational methods. In this review, we briefly introduce computational drug discovery strategies and outline up-to-date tools to perform the strategies as well as available knowledge bases for those who develop their own computational models. Finally, we introduce successful examples of anti-bacterial, anti-viral, and anti-cancer drug discoveries that were made using computational methods.

Single-Agent and Fixed-Dose Combination HIV-1 Protease Inhibitor Drugs in Fission Yeast (Schizosaccharomyces pombe)

Successful combination antiretroviral therapies (cART) eliminate active replicating HIV-1, slow down disease progression, and prolong lives. However, cART effectiveness could be compromised by the emergence of viral multidrug resistance, suggesting the need for new drug discoveries. The objective of this study was to further demonstrate the utility of the fission yeast cell-based systems that we developed previously for the discovery and testing of HIV protease (PR) inhibitors (PIs) against wild-type or multi-PI drug resistant _M11PR that we isolated from an infected individual. All thirteen FDA-approved single-agent and fixed-dose combination HIV PI drugs were tested. The effect of these drugs on HIV PR activities was tested in pure compounds or formulation drugs. All FDA-approved PI drugs, except for a prodrug FPV, were able to suppress the wild-type PR-induced cellular and enzymatic activities. Relative drug potencies measured by EC₅₀ in fission yeast were discussed in comparison with those measured in human cells. In contrast, none of the FDA-approved drugs suppressed the multi-PI drug resistant _M11PR activities. Results of this study show that fission yeast is a reliable cell-based system for the discovery and testing of HIV PIs and further demonstrate the need for new PI drugs against viral multi-PI resistance.

In silico validation of novel inhibitors of malarial aspartyl protease, plasmepsin V and antimalarial efficacy prediction

Plasmepsin V (Plm V) is an essential aspartic protease required for survival of the malaria parasite, Plasmodium falciparum (Pf). Plm V is required for cleaving the PEXEL motifs of many Pf proteins and its inhibition leads to a knockout effect, indicating its suitability as potential drug target. To decipher new inhibitors of PfPlm V, molecular docking of four HIV-1 protease inhibitors active against PfPlmV was performed on Glide module of Schrödinger suite that supported saquinavir as a lead drug, and therefore, selected as a control. Saquinavir contains an important hydroxyethylamine (HEA) pharmacophore, which was utilized as backbone coupled with piperazine scaffold to build new library of compounds. Newly designed HEA compounds were screened virtually against Plm V. Molecular docking led to a few hits (1 and 3) with higher docking score over the control drug. Notably, compound 1 showed the highest docking score (-11.90 kcal/mol) and XP Gscore (-11.948 kcal/mol). The Prime MMGBSA binding free energy for compound 1 (-60.88 kcal/mol) and 3 (-50.96 kcal/mol) was higher than saquinavir (-37.51 kcal/mol). The binding free energy for the last frame of molecular dynamic simulation supported compound 1 (-92.88 kcal/mol) as potent inhibitor of PfPlm V over saquinavir (-72.77 kcal/mol), and thus, deserves experimental validations in culture and subsequently in animal models.Communicated by Ramaswamy H. Sarma.

Structural Insights to Human Immunodeficiency Virus (HIV-1) Targets and Their Inhibition

Human immunodeficiency virus (HIV) is a deadly virus that attacks the body's immune system, subsequently leading to AIDS (acquired immunodeficiency syndrome) and ultimately death. Currently, there is no vaccine or effective cure for this infection; however, antiretrovirals that act at various phases of the virus life cycle have been useful to control the viral load in patients. One of the major problems with antiretroviral therapies involves drug resistance. The three-dimensional structure from crystallography studies are instrumental in understanding the structural basis of drug binding to various targets. This chapter provides key insights into different targets and drugs used in the treatment from a structural perspective. Specifically, an insight into the binding characteristics of drugs at the active and allosteric sites of different targets and the importance of targeting allosteric sites for design of new-generation antiretrovirals to overcome complex and resistant forms of the virus has been reviewed.

An Overview of Antiretroviral Agents for Treating HIV Infection in Paediatric Population

Paediatric Acquired ImmunoDeficiency Syndrome (AIDS) is a life-threatening and infectious disease in which the Human Immunodeficiency Virus (HIV) is mainly transmitted through Mother-To- Child Transmission (MTCT) during pregnancy, labour and delivery, or breastfeeding. This review provides an overview of the distinct therapeutic alternatives to abolish the systemic viral replication in paediatric HIV-1 infection. Numerous classes of antiretroviral agents have emerged as therapeutic tools for downregulation of different steps in the HIV replication process. These classes encompass Non- Nucleoside Analogue Reverse Transcriptase Inhibitors (NNRTIs), Nucleoside/Nucleotide Analogue Reverse Transcriptase Inhibitors (NRTIs/NtRTIs), INtegrase Inhibitors (INIs), Protease Inhibitors (PIs), and Entry Inhibitors (EIs). Co-administration of certain antiretroviral drugs with Pharmacokinetic Enhancers (PEs) may boost the effectiveness of the primary therapeutic agent. The combination of multiple antiretroviral drug regimens (Highly Active AntiRetroviral Therapy - HAART) is currently the standard therapeutic approach for HIV infection. So far, the use of HAART offers the best opportunity for prolonged and maximal viral suppression, and preservation of the immune system upon HIV infection. Still, the frequent administration of high doses of multiple drugs, their inefficient ability to reach the viral reservoirs in adequate doses, the development of drug resistance, and the lack of patient compliance compromise the complete HIV elimination. The development of nanotechnology-based drug delivery systems may enable targeted delivery of antiretroviral agents to inaccessible viral reservoir sites at therapeutic concentrations. In addition, the application of Computer-Aided Drug Design (CADD) approaches has provided valuable tools for the development of anti-HIV drug candidates with favourable pharmacodynamics and pharmacokinetic properties.

Protease Inhibitors for the Treatment of HIV/AIDS: Recent Advances and Future Challenges

Acquired Immunodeficiency Syndrome (AIDS) is a chronic disease characterized by multiple life-threatening illnesses caused by a retro-virus, Human Immunodeficiency Virus (HIV). HIV infection slowly destroys the immune system and increases the risk of various other infections and diseases. Although, there is no immediate cure for HIV infection/AIDS, several drugs targeting various cruxes of HIV infection are used to slow down the progress of the disease and to boost the immune system. One of the key therapeutic strategies is Highly Active Antiretroviral Therapy (HAART) or ' AIDS cocktail' in a general sense, which is a customized combination of anti-retroviral drugs designed to combat the HIV infection. Since HAART's inception in 1995, this treatment was found to be effective in improving the life expectancy of HIV patients over two decades. Among various classes of HAART treatment regimen, Protease Inhibitors (PIs) are known to be widely used as a major component and found to be effective in treating HIV infection/AIDS. For the past several years, a variety of protease inhibitors have been reported. This review outlines the drug design strategies of PIs, chemical and pharmacological characteristics of some mechanism-based inhibitors, summarizes the recent developments in small molecule based drug discovery with HIV protease as a drug target. Further discussed are the pharmacology, PI drug resistance on HIV PR, adverse effects of HIV PIs and challenges/impediments in the successful application of HIV PIs as an important class of drugs in HAART regimen for the effective treatment of AIDS.

A Structural View on Medicinal Chemistry Strategies against Drug Resistance

The natural phenomenon of drug resistance is a widespread issue that hampers the performance of drugs in many major clinical indications. Antibacterial and antifungal drugs are affected, as well as compounds for the treatment of cancer, viral infections, or parasitic diseases. Despite the very diverse set of biological targets and organisms involved in the development of drug resistance, the underlying molecular mechanisms have been identified to understand the emergence of resistance and to overcome this detrimental process. Detailed structural information on the root causes for drug resistance is nowadays frequently available, so next-generation drugs can be designed that are anticipated to suffer less from resistance. This knowledge-based approach is essential for fighting the inevitable occurrence of drug resistance.

Nature Inspired Molecular Design: Stereoselective Synthesis of Bicyclic and Polycyclic Ethers for Potent HIV-1 Protease Inhibitors

We have developed a conceptually new generation of non-peptidic HIV-1 protease inhibitors incorporating novel structural templates inspired by nature. This has resulted in protease inhibitors with exceptional potency and excellent pharmacological and drug-resistance profiles. The design of a stereochemically defined bis-tetrahydrofuran (bis-THF) scaffold followed by modifications to promote hydrogen bonding interactions with the backbone atoms of HIV-1 protease led to darunavir, the first clinically approved drug for treatment of drug resistant HIV. Subsequent X-ray crystal structure-based design efforts led us to create a range of exceptionally potent inhibitors incorporating other intriguing molecular templates possessing fused ring polycyclic ethers with multiple stereocenters. These structural templates are critical to inhibitors' exceptional potency and drug-like properties. Herein, we will highlight the synthetic strategies that provided access to these complex scaffolds in a stereoselective and optically active form, enabling our medicinal chemistry and drug development efforts.

Proteases and protease inhibitors in infectious diseases

There are numerous proteases of pathogenic organisms that are currently targeted for therapeutic intervention along with many that are seen as potential drug targets. This review discusses the chemical and biological makeup of some key druggable proteases expressed by the five major classes of disease causing agents, namely bacteria, viruses, fungi, eukaryotes, and prions. While a few of these enzymes including HIV protease and HCV NS3-4A protease have been targeted to a clinically useful level, a number are yet to yield any clinical outcomes in terms of antimicrobial therapy. A significant aspect of this review discusses the chemical and pharmacological characteristics of inhibitors of the various proteases discussed. A total of 25 inhibitors have been considered potent and safe enough to be trialed in humans and are at different levels of clinical application. We assess the mechanism of action and clinical performance of the protease inhibitors against infectious agents with their developmental strategies and look to the next frontiers in the use of protease inhibitors as anti-infective agents.

Application of computational methods for anticancer drug discovery, design, and optimization

Developing a novel drug is a complex, risky, expensive and time-consuming venture. It is estimated that the conventional drug discovery process ending with a new medicine ready for the market can take up to 15 years and more than a billion USD. Fortunately, this scenario has recently changed with the arrival of new approaches. Many novel technologies and methodologies have been developed to increase the efficiency of the drug discovery process, and computational methodologies have become a crucial component of many drug discovery programs. From hit identification to lead optimization, techniques such as ligand- or structure-based virtual screening are widely used in many discovery efforts. It is the case for designing potential anticancer drugs and drug candidates, where these computational approaches have had a major impact over the years and have provided fruitful insights into the field of cancer. In this paper, we review the concept of rational design presenting some of the most representative examples of molecules identified by means of it. Key principles are illustrated through case studies including specifically successful achievements in the field of anticancer drug design to demonstrate that research advances, with the aid of in silico drug design, have the potential to create novel anticancer drugs.

Application of computational methods for anticancer drug discovery, design, and optimization

Developing a novel drug is a complex, risky, expensive and time-consuming venture. It is estimated that the conventional drug discovery process ending with a new medicine ready for the market can take up to 15 years and more than a billion USD. Fortunately, this scenario has recently changed with the arrival of new approaches. Many novel technologies and methodologies have been developed to increase the efficiency of the drug discovery process, and computational methodologies have become a crucial component of many drug discovery programs. From hit identification to lead optimization, techniques such as ligand- or structure-based virtual screening are widely used in many discovery efforts. It is the case for designing potential anticancer drugs and drug candidates, where these computational approaches have had a major impact over the years and have provided fruitful insights into the field of cancer. In this paper, we review the concept of rational design presenting some of the most representative examples of molecules identified by means of it. Key principles are illustrated through case studies including specifically successful achievements in the field of anticancer drug design to demonstrate that research advances, with the aid of in silico drug design, have the potential to create novel anticancer drugs.

Computer simulations of the interaction of human immunodeficiency virus (HIV) aspartic protease with spherical gold nanoparticles: implications in acquired immunodeficiency syndrome (AIDS)

The interaction of gold nanoparticles (AuNP) with human immune-deficiency virus aspartic protease (HIVPR) is modelled using a regime of molecular dynamics simulations. The simulations of the 'docking', first as a rigid-body complex, and eventually through flexible-fit analysis, creates 36 different complexes from four initial orientations of the nanoparticle strategically positioned around the surface of the enzyme. The structural deviations of the enzymes from the initial x-ray crystal structure during each docking simulation are assessed by comparative analysis of secondary structural elements, root mean square deviations, B-factors, interactive bonding energies, dihedral angles, radius of gyration (R g), circular dichroism (CD), volume occupied by C α , electrostatic potentials, solvation energies and hydrophobicities. Normalisation of the data narrows the selection from the initial 36 to one 'final' probable structure. It is concluded that, after computer simulations on each of the 36 initial complexes incorporating the 12 different biophysical techniques, the top five complexes are the same no matter which technique is explored. The significance of the present work is an expansion of an earlier study on the molecular dynamic simulation for the interaction of HIVPR with silver nanoparticles. This work is supported by experimental evidence since the initial 'orientation' of the AgNP with the enzyme is the same as the 'final' AuNP-HIVPR complex generated in the present study. The findings will provide insight into the forces of the binding of the HIVPR to AuNP. It is anticipated that the protocol developed in this study will act as a standard process for the interaction of any nanoparticle with any biomedical target.

Recent Progress in the Development of HIV-1 Protease Inhibitors for the Treatment of HIV/AIDS

HIV-1 protease inhibitors continue to play an important role in the treatment of HIV/AIDS, transforming this deadly ailment into a more manageable chronic infection. Over the years, intensive research has led to a variety of approved protease inhibitors for the treatment of HIV/AIDS. In this review, we outline current drug design and medicinal chemistry efforts toward the development of next-generation protease inhibitors beyond the currently approved drugs.

Unique Flap Conformation in an HIV-1 Protease with High-Level Darunavir Resistance

Darunavir (DRV) is one of the most powerful protease inhibitors (PIs) for treating human immunodeficiency virus type-1 (HIV-1) infection and presents a high genetic barrier to the generation of resistant viruses. However, DRV-resistant HIV-1 infrequently emerges from viruses exhibiting resistance to other protease inhibitors. To address this resistance, researchers have gathered genetic information on DRV resistance. In contrast, few structural insights into the mechanism underlying DRV resistance are available. To elucidate this mechanism, we determined the crystal structure of the ligand-free state of a protease with high-level DRV resistance and six DRV resistance-associated mutations (including I47V and I50V), which we generated by in vitro selection. This crystal structure showed a unique curling conformation at the flap regions that was not found in the previously reported ligand-free protease structures. Molecular dynamics simulations indicated that the curled flap conformation altered the flap dynamics. These results suggest that the preference for a unique flap conformation influences DRV binding. These results provide new structural insights into elucidating the molecular mechanism of DRV resistance and aid to develop PIs effective against DRV-resistant viruses.

Binding Free Energy Calculations of Nine FDA-approved Protease Inhibitors Against HIV-1 Subtype C I36T↑T Containing 100 Amino Acids Per Monomer

In this work, have investigated the binding affinities of nine FDA-approved protease inhibitor drugs against a new HIV-1 subtype C mutated protease, I36T↑T. Without an X-ray crystal structure, homology modelling was used to generate a three-dimensional model of the protease. This and the inhibitor models were employed to generate the inhibitor/I36T↑T complexes, with the relative positions of the inhibitors being superimposed and aligned using the X-ray crystal structures of the inhibitors/HIV-1 subtype B complexes as a reference. Molecular dynamics simulations were carried out on the complexes to calculate the average binding free energies for each inhibitor using the molecular mechanics generalized Born surface area (MM-GBSA) method. When compared to the binding free energies of the HIV-1 subtype B and subtype C proteases (calculated previously by our group using the same method), it was clear that the I36T↑T proteases mutations and insertion had a significant negative effect on the binding energies of the non-pepditic inhibitors nelfinavir, darunavir and tipranavir. On the other hand, ritonavir, amprenavir and indinavir show improved calculated binding energies in comparison with the corresponding data for wild-type C-SA protease. The computational model used in this study can be used to investigate new mutations of the HIV protease and help in establishing effective HIV drug regimes and may also aid in future protease drug design.

HIV protease inhibitors: a review of molecular selectivity and toxicity

Highly active antiretroviral therapy (HAART) is recognized as the most effective treatment method for AIDS, and protease inhibitors play a very important role in HAART. However, poor bioavailability and unbearable toxicity are their common disadvantages. Thus, the development of safer and potentially promising protease inhibitors is eagerly needed. In this review, we introduced the chemical characteristics and associated side effects of HIV protease inhibitors, as well as the possible off-target mechanisms causing the side effects. From the chemical structures of HIV protease inhibitors and their possible off-target molecules, we could obtain hints for optimizing the molecular selectivity of the inhibitors, to provide help in the design of new compounds with enhanced bioavailability and reduced side effects.

The impact of active site mutations of South African HIV PR on drug resistance: Insight from molecular dynamics simulations, binding free energy and per-residue footprints

Molecular dynamics simulations and binding free energy calculations were used to provide an understanding of the impact of active site drug-resistant mutations of the South African HIV protease subtype C (C-SA HIV PR), V82A and V82F/I84V on drug resistance. Unique per-residue interaction energy 'footprints' were developed to map the overall drug-binding profiles for the wild type and mutants. Results confirmed that these mutations altered the overall binding landscape of the amino acid residues not only in the active site region but also in the flaps as well. Four FDA-approved drugs were investigated in this study; these include ritonavir (RTV), saquinavir (SQV), indinavir (IDV), and nelfinavir (NFV). Computational results compared against experimental findings were found to be complementary. Against the V82F/I84V variant, saquinavir, indinavir, and nelfinavir lose remarkable entropic contributions relative to both wild-type and V82A C-SA HIV PRs. The per-residue energy 'footprints' and the analysis of ligand-receptor interactions for the drug complexes with the wild type and mutants have also highlighted the nature of drug interactions. The data presented in this study will prove useful in the design of more potent inhibitors effective against drug-resistant HIV strains.

Defective hydrophobic sliding mechanism and active site expansion in HIV-1 protease drug resistant variant Gly48Thr/Leu89Met: mechanisms for the loss of saquinavir binding potency

HIV drug resistance continues to emerge; consequently, there is an urgent need to develop next generation antiretroviral therapeutics.1 Here we report on the structural and kinetic effects of an HIV protease drug resistant variant with the double mutations Gly48Thr and Leu89Met (PR_G48T/L89M), without the stabilizing mutations Gln7Lys, Leu33Ile, and Leu63Ile. Kinetic analyses reveal that PR_G48T/L89M and PR_WT share nearly identical Michaelis–Menten parameters; however, PR_G48T/L89M exhibits weaker binding for IDV (41-fold), SQV (18-fold), APV (15-fold), and NFV (9-fold) relative to PR_WT. A 1.9 Å resolution crystal structure was solved for PR_G48T/L89M bound with saquinavir (PR_{G48T/L89M-SQV}) and compared to the crystal structure of PR_WT bound with saquinavir (PR_WT-SQV). PR_{G48T/L89M-SQV} has an enlarged active site resulting in the loss of a hydrogen bond in the S3 subsite from Gly48 to P3 of SQV, as well as less favorable hydrophobic packing interactions between P1 Phe of SQV and the S1 subsite. PR_{G48T/L89M-SQV} assumes a more open conformation relative to PR_WT-SQV, as illustrated by the downward displacement of the fulcrum and elbows and weaker interatomic flap interactions. We also show that the Leu89Met mutation disrupts the hydrophobic sliding mechanism by causing a redistribution of van der Waals interactions in the hydrophobic core in PR_{G48T/L89M-SQV}. Our mechanism for PR_{G48T/L89M-SQV} drug resistance proposes that a defective hydrophobic sliding mechanism results in modified conformational dynamics of the protease. As a consequence, the protease is unable to achieve a fully closed conformation that results in an expanded active site and weaker inhibitor binding.

The triple threat of HIV-1 protease inhibitors

Newly released human immunodeficiency virus type 1 (HIV-1) particles obligatorily undergo a maturation process to become infectious. The HIV-1 protease (PR) initiates this step, catalyzing the cleavage of the Gag and Gag-Pro-Pol structural polyproteins. Proper organization of the mature virus core requires that cleavage of these polyprotein substrates proceeds in a highly regulated, specific series of events. The vital role the HIV-1 PR plays in the viral life cycle has made it an extremely attractive target for inhibition and has accordingly fostered the development of a number of highly potent substrate-analog inhibitors. Though the PR inhibitors (PIs) inhibit only the HIV-1 PR, their effects manifest at multiple different stages in the life cycle due to the critical importance of the PR in preparing the virus for these subsequent events. Effectively, PIs masquerade as entry inhibitors, reverse transcription inhibitors, and potentially even inhibitors of post-reverse transcription steps. In this chapter, we review the triple threat of PIs: the intermolecular cooperativity in the form of a cooperative dose-response for inhibition in which the apparent potency increases with increasing inhibition; the pleiotropic effects of HIV-1 PR inhibition on entry, reverse transcription, and post-reverse transcription steps; and their potency as transition state analogs that have the potential for further improvement that could lead to an inability of the virus to evolve resistance in the context of single drug therapy.

Evaluation of novel Saquinavir analogs for resistance mutation compatibility and potential as an HIV-Protease inhibitor drug

A fundamental issue related to therapy of HIV-1 infection is the emergence of viral mutations which severely limits the long term efficiency of the HIV-protease (HIV-PR) inhibitors. Development of new drugs is therefore continuously needed. Chemoinformatics enables to design and discover novel molecules analogous to established drugs using computational tools and databases. Saquinavir, an anti-HIV Protease drug is administered for HIV therapy. In this work chemoinformatics tools were used to design structural analogs of Saquinavir as ligand and molecular dockings at AutoDock were performed to identify potential HIV-PR inhibitors. The analogs S1 and S2 when docked with HIV-PR had binding energies of -4.08 and -3.07 kcal/mol respectively which were similar to that for Saquinavir. The molecular docking studies revealed that the changes at N2 of Saquinavir to obtain newly designed analogs S1 (having N2 benzoyl group at N1) and S2 (having 3-oxo-3phenyl propanyl group at N2) were able to dock with HIV-PR with similar affinity as that of Saquinavir. Docking studies and computationally derived pharmacodynamic and pharmacokinetic properties׳ comparisons at ACD/I-lab establish that analog S2 has more potential to evade the problem of drug resistance mutation against HIV-1 PR subtype-A. S2 can be further developed and tested clinically as a real alternative drug for HIV-1 PR across the clades in future.

A sensitive assay using a native protein substrate for screening HIV-1 maturation inhibitors targeting the protease cleavage site between the matrix and capsid

The matrix/capsid processing site in the HIV-1 Gag precursor is likely the most sensitive target to inhibit HIV-1 replication. We have previously shown that modest incomplete processing at the site leads to a complete loss of virion infectivity. In the study presented here, a sensitive assay based on fluorescence polarization that can monitor cleavage at the MA/CA site in the context of the folded protein substrate is described. The substrate, an MA/CA fusion protein, was labeled with the fluorescein-based FlAsH (fluorescein arsenical hairpin) reagent that binds to a tetracysteine motif (CCGPCC) that was introduced within the N-terminal domain of CA. By limiting the size of CA and increasing the size of MA (with an N-terminal GST fusion), we were able to measure significant differences in polarization values as a function of HIV-1 protease cleavage. The sensitivity of the assay was tested in the presence of increasing amounts of an HIV-1 protease inhibitor, which resulted in a gradual decrease in the fluorescence polarization values demonstrating that the assay is sensitive in discerning changes in protease processing. The high-throughput screening assay validation in 384-well plates showed that the assay is reproducible and robust with an average Z' value of 0.79 and average coefficient of variation values of <3%. The robustness and reproducibility of the assay were further validated using the LOPAC(1280) compound library, demonstrating that the assay provides a sensitive high-throughput screening platform that can be used with large compound libraries for identifying novel maturation inhibitors targeting the MA/CA site of the HIV-1 Gag polyprotein.

Trendspotting in the Protein Data Bank

The Protein Data Bank (PDB) was established in 1971 as a repository for the three dimensional structures of biological macromolecules. Since then, more than 85000 biological macromolecule structures have been determined and made available in the PDB archive. Through analysis of the corpus of data, it is possible to identify trends that can be used to inform us abou the future of structural biology and to plan the best ways to improve the management of the ever-growing amount of PDB data.

Extreme entropy-enthalpy compensation in a drug-resistant variant of HIV-1 protease

The development of HIV-1 protease inhibitors has been the historic paradigm of rational structure-based drug design, where structural and thermodynamic analyses have assisted in the discovery of novel inhibitors. While the total enthalpy and entropy change upon binding determine the affinity, often the thermodynamics are considered in terms of inhibitor properties only. In the current study, profound changes are observed in the binding thermodynamics of a drug-resistant variant compared to wild-type HIV-1 protease, irrespective of the inhibitor bound. This variant (Flap+) has a combination of flap and active site mutations and exhibits extremely large entropy-enthalpy compensation compared to wild-type protease, 5-15 kcal/mol, while losing only 1-3 kcal/mol in total binding free energy for any of six FDA-approved inhibitors. Although entropy-enthalpy compensation has been previously observed for a variety of systems, never have changes of this magnitude been reported. The co-crystal structures of Flap+ protease with four of the inhibitors were determined and compared with complexes of both the wild-type protease and another drug-resistant variant that does not exhibit this energetic compensation. Structural changes conserved across the Flap+ complexes, which are more pronounced for the flaps covering the active site, likely contribute to the thermodynamic compensation. The finding that drug-resistant mutations can profoundly modulate the relative thermodynamic properties of a therapeutic target independent of the inhibitor presents a new challenge for rational drug design.

Protein pharmacophore selection using hydration-site analysis

Virtual screening using pharmacophore models is an efficient method to identify potential lead compounds for target proteins. Pharmacophore models based on protein structures are advantageous because a priori knowledge of active ligands is not required and the models are not biased by the chemical space of previously identified actives. However, in order to capture most potential interactions between all potentially binding ligands and the protein, the size of the pharmacophore model, i.e. number of pharmacophore elements, is typically quite large and therefore reduces the efficiency of pharmacophore based screening. We have developed a new method to select important pharmacophore elements using hydration-site information. The basic premise is that ligand functional groups that replace water molecules in the apo protein contribute strongly to the overall binding affinity of the ligand, due to the additional free energy gained from releasing the water molecule into the bulk solvent. We computed the free energy of water released from the binding site for each hydration site using thermodynamic analysis of molecular dynamics (MD) simulations. Pharmacophores which are colocalized with hydration sites with estimated favorable contributions to the free energy of binding are selected to generate a reduced pharmacophore model. We constructed reduced pharmacophore models for three protein systems and demonstrated good enrichment quality combined with high efficiency. The reduction in pharmacophore model size reduces the required screening time by a factor of 200-500 compared to using all protein pharmacophore elements. We also describe a training process using a small set of known actives to reliably select the optimal set of criteria for pharmacophore selection for each protein system.

Enhancing protein backbone binding--a fruitful concept for combating drug-resistant HIV

The evolution of drug resistance is one of the most fundamental problems in medicine. In HIV/AIDS, the rapid emergence of drug-resistant HIV-1 variants is a major obstacle to current treatments. HIV-1 protease inhibitors are essential components of present antiretroviral therapies. However, with these protease inhibitors, resistance occurs through viral mutations that alter inhibitor binding, resulting in a loss of efficacy. This loss of potency has raised serious questions with regard to effective long-term antiretroviral therapy for HIV/AIDS. In this context, our research has focused on designing inhibitors that form extensive hydrogen-bonding interactions with the enzyme's backbone in the active site. In doing so, we limit the protease's ability to acquire drug resistance as the geometry of the catalytic site must be conserved to maintain functionality. In this Review, we examine the underlying principles of enzyme structure that support our backbone-binding concept as an effective means to combat drug resistance and highlight their application in our recent work on antiviral HIV-1 protease inhibitors.

THE ϕ, ψ SPACE OF BORON ISOSTERES OF AMINO ACIDS: ANAB INITIOSTUDY

Boron isosteres [ CH3- CO - BH - CH(R) - CONHCH3; R = CH3, CH2OH , CH2SH ] of natural amino acids have been designed as inhibitors of serine protease based on boron's ability to mimic the transition state of the normal enzyme catalyzed reaction. The conformation and electronic properties of these isosteres have been determined by ab initio calculations at the HF/6-31G** level of theory. The global minimum energy structures of the boron isosteres (which are all identical) and the natural amino acids are radically different. Interestingly, the boron isosteres have favorable conformations with positive values for the ϕ dihedral, which are "disallowed states" for the natural amino acids. This is possible because of the smaller atomic size of boron compared to nitrogen. The molecular electrostatic potential of the boron isosteres is also very different from the natural amino acids. There is a positive potential around the boron atom, which provokes the attack of the Ser195– OH group in the enzyme, resulting in formation of an irreversible enzyme-inhibitor adduct.

DESIGN OF HYBRID INHIBITORS TO HIV-1 PROTEASE

A series of HIV-1 protease (PR) inhibitors are designed to increase the binding affinity with PR subsites based on the quantum analysis of the contributions of molecular fragments in six FDA-approved PR drugs to the total binding interaction. The binding free energies were estimated by modified linear interaction energy approach [Zoete H, Michielin O, Karplus M, J Comput Aided Mol Des17:861, 2003], in which the binding free energy is written as a linear combination of the electrostatic interaction energy between PR and the ligand, Eelec, the van der Waals interaction energy between PR and the ligand, E vdW , and the difference of the solvation free energies of the complex, the receptor and the isolated ligand, ΔG solv . The parameters of these energy terms were fitted for a training set of 14 HIV-1 protease–inhibitor complexes of known 3D structure with a correlation coefficient of 0.91 and an unsigned mean error of 0.83 kcal/mol.

[Bioinformatics studies on drug resistance against anti-HIV-1 drugs]

More than 20 drugs have been available for anti-HIV-1 treatment in Japan. Combination therapy with these drugs dramatically decreases in morbidity and mortality of AIDS. However, due to high mutation rate of HIV-1, treatment with ineffective drugs toward patients infected with HIV-1 causes accumulation of mutations in the virus, and emergence of drug resistant viruses. Thus, to achieve appropriate application of the drugs toward the respective patients living with HIV-1, methods for predicting the level of drug-resistance using viral sequence information has been developed on the basis of bioinformatics. Furthermore, ultra-deep sequencing by next-generation sequencer whose data analysis is also based on bioinformatics, or in silico structural modeling have been achieved to understand drug resistant mechanisms. In this review, I overview the bioinformatics studies about drug resistance against anti-HIV-1 drugs.

Computational mutation scanning and drug resistance mechanisms of HIV-1 protease inhibitors

The drug resistance of various clinically available HIV-1 protease inhibitors has been studied using a new computational protocol, that is, computational mutation scanning (CMS), leading to valuable insights into the resistance mechanisms and structure-resistance correction of the HIV-1 protease inhibitors associated with a variety of active site and nonactive site mutations. By using the CMS method, the calculated mutation-caused shifts of the binding free energies linearly correlate very well with those derived from the corresponding experimental data, suggesting that the CMS protocol may be used as a generalized approach to predict drug resistance associated with amino acid mutations. Because it is essentially important for understanding the structure-resistance correlation and for structure-based drug design to develop an effective computational protocol for drug resistance prediction, the reasonable and computationally efficient CMS protocol for drug resistance prediction should be valuable for future structure-based design and discovery of antiresistance drugs in various therapeutic areas.

Rationale for more diverse inhibitors in competition with substrates in HIV-1 protease

The structural fluctuations of HIV-1 protease in interaction with its substrates versus inhibitors were analyzed using the anisotropic network model. The directions of fluctuations in the most cooperative functional modes differ mainly around the dynamically key regions, i.e., the hinge axes, which appear to be more flexible in substrate complexes. The flexibility of HIV-1 protease is likely optimized for the substrates' turnover, resulting in substrate complexes being dynamic. In contrast, in an inhibitor complex, the inhibitor should bind and lock down to inactivate the active site. Protease and ligands are not independent. Substrates are also more flexible than inhibitors and have the potential to meet the dynamic distributions that are inherent in the protease. This may suggest a rationale and guidelines for designing inhibitors that can better fit the ensemble of binding sites that are dynamically accessible to the protease.

Amprenavir complexes with HIV-1 protease and its drug-resistant mutants altering hydrophobic clusters

The structural and kinetic effects of amprenavir (APV), a clinical HIV protease (PR) inhibitor, were analyzed with wild-type enzyme and mutants with single substitutions of V32I, I50V, I54V, I54M, I84V and L90M that are common in drug resistance. Crystal structures of the APV complexes at resolutions of 1.02-1.85 Å reveal the structural changes due to the mutations. Substitution of the larger side chains in PR(V32I) , PR(I54M) and PR(L90M) resulted in the formation of new hydrophobic contacts with flap residues, residues 79 and 80, and Asp25, respectively. Mutation to smaller side chains eliminated hydrophobic interactions in the PR(I50V) and PR(I54V) structures. The PR(I84V)-APV complex had lost hydrophobic contacts with APV, the PR(V32I)-APV complex showed increased hydrophobic contacts within the hydrophobic cluster and the PR(I50V) complex had weaker polar and hydrophobic interactions with APV. The observed structural changes in PR(I84V)-APV, PR(V32I)-APV and PR(I50V)-APV were related to their reduced inhibition by APV of six-, 10- and 30-fold, respectively, relative to wild-type PR. The APV complexes were compared with the corresponding saquinavir complexes. The PR dimers had distinct rearrangements of the flaps and 80's loops that adapt to the different P1' groups of the inhibitors, while maintaining contacts within the hydrophobic cluster. These small changes in the loops and weak internal interactions produce the different patterns of resistant mutations for the two drugs.

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.