Binding modes of CCR5-targetting HIV entry inhibitors: partial and full antagonists

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.

Binding Interactions of Ergotamine and Dihydroergotamine to 5-Hydroxytryptamine Receptor 1B (5-HT1b) Using Molecular Dynamics Simulations and Dynamic Network Analysis

Ergotamine (ERG) and dihydroergotamine (DHE), common migraine drugs, have small structural differences but lead to clinically important distinctions in their pharmacological profiles. For example, DHE is less potent than ERG by about 10-fold at the 5-hydroxytrptamine receptor 1B (5-HT_1B). Although the high-resolution crystal structures of the 5-HT_1B receptor with both ligands have been solved, the high similarity between these two complex structures does not sufficiently explain their activity differences and the activation mechanism of the receptor. Hence, an examination of the dynamic motion of both drugs with the receptor is required. In this study, we ran a total of 6.0 μs molecular dynamics simulations on each system. Our simulation data show the subtle variations between the two systems in terms of the ligand-receptor interactions and receptor secondary structures. More importantly, the ligand and protein root-mean-square fluctuations (RMSFs) for the two systems were distinct, with ERG having a trend of lower RMSF values, indicating it to be bound tighter to 5-HT_1B with less fluctuations. The molecular mechanism-general born surface area (MM-GBSA) binding energies illustrate this further, proving ERG has an overall stronger MM-GBSA binding energy. Analysis of several different microswitches has shown that the 5-HT_1B-ERG complex is in a more active conformation state than 5-HT_1B-DHE, which is further supported by the dynamic network model, with reference to mutagenesis data with the critical nodes and the first three low-energy modes from the normal mode analysis. We also identify Trp327^6.48 and Phe331^6.52 as key residues involved in the active state 5-HT_1B for both ligands. Using the detailed dynamic information from our analysis, we made predictions for possible modifications to DHE and ERG that yielded five derivatives that might have more favorable binding energies and reduced structural fluctuations.

Pharmacophore modeling, atom based 3D-QSAR and Molecular docking approaches to screen C-X-C chemokine receptor type 4 antagonists as microbicides for human immunodeficiency virus-1

The C-X-C chemokine receptor type 4 receptor CXCR4 which acts as a co-receptor for human immunodeficiency virus-1, expressed in the later stages of infection is considered as an attractive and new target for drug design. Microbicides acting as co-receptor blockers are highly significant as these drugs block HIV lifecycle at early stage itself. The urgent need for a safe and effective microbicide urges to explore new CXCR4 antagonists which could be developed as microbicides. The pharmacophore based 3D-QSAR models and docking models were developed using PHASE and GLIDE modules of Schrodinger software. Three-dimensional quantitative structure-activity relationships (3D-QSAR) studies and pharmacophore modelling was carried out on a dataset of 114 CXCR4 antagonists with the intention of exploring entry inhibitors with better therapeutic potential. A training set of 43 compounds was used to create 3D-QSAR models and they were validated using a test set of 28 compounds. CXCR4 antagonists with good inhibitory activity could be designed and structurally modified based upon the QSAR model developed with necessary pharmacophore features. The results revealed that the common pharmacophore hypothesis ADHPR.1 was used for 3D-QSAR model development and the most active compound, CXCR4 antagonist no.44 which is a imidazopyridine-tetrahydro-8-quinolinamine derivative interacted with the CXCR4 receptor residue ASP 97 by the formation of a hydrogen bond. Also, the docking studies were carried out for the dataset for analyzing the binding conformation of CXCR4 and 114 antagonists. The results obtained from the 3D-QSAR studies and docking simulation can be used for designing new and potent CXCR4 antagonists. The compound identified from this study can be taken up further for validation by in vitro/in vivo studies.

Investigation of Inhibition Mechanism of Chemokine Receptor CCR5 by Micro-second Molecular Dynamics Simulations

Chemokine receptor 5 (CCR5) belongs to G protein coupled receptors (GPCRs) and plays an important role in treatment of human immunodeficiency virus (HIV) infection since HIV uses CCR5 protein as a co-receptor. Recently, the crystal structure of CCR5-bound complex with an approved anti-retroviral drug (maroviroc) was resolved. During the crystallization procedure, amino acid residues (i.e., Cys224, Arg225, Asn226 and Glu227) at the third intra-cellular loop were replaced by the rubredoxin for stability reasons. In the current study, we aimed to understand the impact of the incorporated rubredoxin on the conformations of TM domains of the target protein. For this reason, rubredoxin was deleted from the crystal structure and the missing amino acids were engineered. The resultant structure was subjected to long (μs) molecular dynamics (MD) simulations to shed light into the inhibitory mechanism. The derived model structure displayed a significant deviation in the cytoplasmic domain of TM5 and IC3 in the absence of rubredoxin. The principal component analyses (PCA) and MD trajectory analyses revealed important structural and dynamical differences at apo and holo forms of the CCR5.

From base pair to bedside: molecular simulation and the translation of genomics to personalized medicine

Despite the promises made that genomic sequencing would transform therapy by introducing a new era of personalized medicine, relatively few tangible breakthroughs have been made. This has led to the recognition that complex interactions at multiple spatial, temporal, and organizational levels may often combine to produce disease. Understanding this complexity requires that existing and future models are used and interpreted within a framework that incorporates knowledge derived from investigations at multiple levels of biological function. It also requires a computational infrastructure capable of dealing with the vast quantities of data generated by genomic approaches. In this review, we discuss the use of molecular modeling to generate quantitative and qualitative insights at the smallest scales of the systems biology hierarchy, how it can play an important role in the development of a systems understanding of disease and in the application of such knowledge to help discover new therapies and target existing ones on a personal level.

Structural insights from binding poses of CCR2 and CCR5 with clinically important antagonists: a combined in silico study

Chemokine receptors are G protein-coupled receptors that contain seven transmembrane domains. In particular, CCR2 and CCR5 and their ligands have been implicated in the pathophysiology of a number of diseases, including rheumatoid arthritis and multiple sclerosis. Based on their roles in disease, they have been attractive targets for the pharmaceutical industry, and furthermore, targeting both CCR2 and CCR5 can be a useful strategy. Owing to the importance of these receptors, information regarding the binding site is of prime importance. Structural studies have been hampered due to the lack of X-ray crystal structures, and templates with close homologs for comparative modeling. Most of the previous models were based on the bovine rhodopsin and β2-adrenergic receptor. In this study, based on a closer homolog with higher resolution (CXCR4, PDB code: 3ODU 2.5 Å), we constructed three-dimensional models. The main aim of this study was to provide relevant information on binding sites of these receptors. Molecular dynamics simulation was done to refine the homology models and PROCHECK results indicated that the models were reasonable. Here, binding poses were checked with some established inhibitors of high pharmaceutical importance against the modeled receptors. Analysis of interaction modes gave an integrated interpretation with detailed structural information. The binding poses confirmed that the acidic residues Glu291 (CCR2) and Glu283 (CCR5) are important, and we also found some additional residues. Comparisons of binding sites of CCR2/CCR5 were done sequentially and also by docking a potent dual antagonist. Our results can be a starting point for further structure-based drug design.

[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.

Theoretical studies on the interactions and interferences of HIV-1 glycoprotein gp120 and its coreceptor CCR5

The interaction between the HIV gp120 protein and coreceptor CCR5 or CXCR4 of the host cell is critical in mediating the HIV entry process. A model for the CCR5-gp120 complex has been developed. In the model, the N-terminus of CCR5 binds to three discontinuous domains of gp120, including the fourth conserved (C4) region, β19/β20 connecting loop, and V3 loop. The second extra-cellular loop (ECL2) of CCR5 also interacts with the crown part of the gp120 V3 loop. The bindings of the three CCR5 antagonists, maraviroc, aplaviroc, and vicriviroc, to the trans-membrane domain of CCR5 have been modeled. The bindings are found to affect the conformation of the ECL2 domain, which in turn drives the N-terminus of CCR5 to an altered state. Aplaviroc is more hydrophilic than maraviroc and vicriviroc, and its binding is more interfered by solvent, resulting in a quite different effect to the structure of CCR5 compared with those of the other two molecules. The above results are in accord with experimental observations and provide a structural basis for further design of CCR5 antagonists.

Pharmacotherapy of HIV-1 Infection: Focus on CCR5 Antagonist Maraviroc

Sustained inhibition of HIV-1, the goal of antiretroviral therapy, is often impeded by the emergence of viral drug resistance. For patients infected with HIV-1 resistant to conventional drugs from the viral reverse transcriptase and protease inhibitor classes, the recently approved entry and integration inhibitors effectively suppress HIV-1 and offer additional therapeutic options. Entry inhibitors are particularly attractive because, unlike conventional antiretrovirals, they target HIV-1 extracellularly, thereby sparing cells from both viral- and drug-induced toxicities. The fusion inhibitor enfuvirtide and the CCR5 antagonist maraviroc are the first entry inhibitors licensed for patients with drug-resistant HIV-1, with maraviroc restricted to those infected with CCR5-tropic HIV-1 (R5 HIV-1) only. Vicriviroc (another CCR5 antagonist) is in Phase III clinical trials, whereas the CCR5 antibodies PRO 140 and HGS 004 are in early stages of clinical development. Potent antiviral synergy between maraviroc and CCR5 antibodies, coupled with distinct patterns of resistance, suggest their combinations might be particularly effective in patients. In addition, given that oral administration of maraviroc achieves high drug levels in cervicovaginal fluid, combinations of maraviroc and other CCR5 inhibitors could be effective in preventing HIV-1 transmission. Moreover, since CCR5 antagonists prevent rejection of transplanted organs, maraviroc could both suppress HIV-1 and prolong organ survival for the growing number of HIV-1 patients with kidney or liver failure necessitating organ transplantation. Thus, maraviroc offers an important treatment option for patients with drug-resistant R5 HIV-1, who presently account for >50% of drug-resistance cases.

Quantitative structure-activity relationships studies of CCR5 inhibitors and toxicity of aromatic compounds using gene expression programming

Quantitative structure-activity relationship (QSAR) study of chemokine receptor 5 (CCR5) binding affinity of substituted 1-(3,3-diphenylpropyl)-piperidinyl amides and ureas and toxicity of aromatic compounds have been performed. The gene expression programming (GEP) was used to select variables and produce nonlinear QSAR models simultaneously using the selected variables. In our GEP implementation, a simple and convenient method was proposed to infer the K-expression from the number of arguments of the function in a gene, without building the expression tree. The results were compared to those obtained by artificial neural network (ANN) and support vector machine (SVM). It has been demonstrated that the GEP is a useful tool for QSAR modeling.

Closing the door on flaviviruses: entry as a target for antiviral drug design

With the emergence and rapid spread of West Nile virus in the United States since 1999, and the 50-100 million infections per year caused by dengue virus globally, the threat of flaviviruses as re-emerging human pathogens has become a reality. To support the efforts that are currently being pursued to develop effective vaccines against these viruses, researchers are also actively pursuing the development of small molecule compounds that target various aspects of the virus life cycle. Recent advances in the structural characterization of the flaviviruses have provided a strong foundation towards these efforts. These studies have provided the pseudo-atomic structures of virions from several members of the genus as well as atomic resolution structures of several viral proteins. Most importantly, these studies have highlighted specific structural rearrangements that occur within the virion that are necessary for the virus to complete its life cycle. These rearrangements occur when the virus must transition from immature, to mature, to fusion-active states and rely heavily on the conformational flexibility of the envelope (E) protein that forms the outer glycoprotein shell of the virus. Analysis of these conformational changes can suggest promising targets for structure-based antiviral design. For instance, by targeting the flexibility of the E protein, it might be possible to inhibit required rearrangements of this protein and trap the virus in a specific state. This would interfere with a productive flaviviral infection. This review presents a structural perspective of the flavivirus life cycle and focuses on the role of the E protein as an opportune target for structure-based antiviral drug design.