In silico characterization of binding mode of CCR8 inhibitor: homology modeling,docking and membrane based MD simulation study

Human CC-chemokine receptor 8 (CCR8) is a crucial drug target in asthma that belongs to G-protein-coupled receptor superfamily, which is characterized by seven transmembrane helices. To date, there is no X-ray crystal structure available for CCR8; this hampers active research on the target. Molecular basis of interaction mechanism of antagonist with CCR8 remains unclear. In order to provide binding site information and stable binding mode, we performed modeling, docking and molecular dynamics (MD) simulation of CCR8. Docking study of biaryl-ether-piperidine derivative (13C) was performed inside predefined CCR8 binding site to get the representative conformation of 13C. Further, MD simulations of receptor and complex (13C-CCR8) inside dipalmitoylphosphatidylcholine lipid bilayers were performed to explore the effect of lipids. Results analyses showed that the Gln91, Tyr94, Cys106, Val109, Tyr113, Cys183, Tyr184, Ser185, Lys195, Thr198, Asn199, Met202, Phe254, and Glu286 were conserved in both docking and MD simulations. This indicated possible role of these residues in CCR8 antagonism. However, experimental mutational studies on these identified residues could be effective to confirm their importance in CCR8 antagonism. Furthermore, calculated Coulombic interactions represented the crucial roles of Glu286, Lys195, and Tyr113 in CCR8 antagonism. Important residues identified in this study overlap with the previous non-peptide agonist (LMD-009) binding site. Though, the non-peptide agonist and currently studied inhibitor (13C) share common substructure, but they differ in their effects on CCR8. So, to get more insight into their agonist and antagonist effects, further side-by-side experimental studies on both agonist (LMD-009) and antagonist (13C) are suggested.     

Probing ligand binding modes of human cytochrome P450 2J2 by homology modeling, molecular dynamics simulation, and flexible molecular docking

Cytochrome P450 (P450) 2J2 catalyzes epoxidation of arachidonic acid to eicosatrienoic acids, which are related to a variety of diseases such as coronary artery disease, hypertension, and carcinogenesis. Recent experimental data also suggest that P450 2J2 could be a novel biomarker and a potential target for cancer therapy. However, the active site topology and substrate specificity of this enzyme remain unclear. In this study, a three-dimensional model of human P450 2J2 was first constructed on the basis of the crystal structure of human P450 2C9 in complex with a substrate using homology modeling method, and refined by molecular dynamics simulation. Flexible docking approaches were then employed to dock four ligands into the active site of P450 2J2 in order to probe the ligand-binding modes. By analyzing the results, active site architecture and certain key residues responsible for substrate specificity were identified on the enzyme, which might be very helpful for understanding the enzyme's biological role and providing insights for designing novel inhibitors of P450 2J2.     

A comparative analysis of binding sites between human PKD1 and PKC1 based on homology modelling,molecular dynamics simulation and docking studies

Human protein kinase C1 (PKC1) and protein kinase D1 (PKD1) are two closely related enzymes, which have emerged as key regulators of many important cellular processes. In this study, 3D models of human PKC1 and PKD1 were constructed based on homology modelling and molecular dynamics simulations. A novel 2,6-naphthyridine is a potent and selective inhibitor for human PKD1 and not for PKC1, which was docked into them and positioned in their active sites with different orientations. By comparison of active site architectures between human PKC1 and PKD1, the possible reasons affecting their inhibitor binding were proposed. In addition, some residues are identified as critical residues for inhibitor binding.     

Structural insights into human microsomal epoxide hydrolase by combined homology modeling,molecular dynamics simulations,and molecular docking calculations

A new homology model of human microsomal epoxide hydrolase was derived based on multiple templates. The model obtained was fully evaluated, including MD simulations and ensemble‐based docking, showing that the quality of the structure is better than that of only previously known model. Particularly, a catalytic triad was clearly identified, in agreement with the experimental information available. Analysis of intermediates in the enzymatic mechanism led to the identification of key residues for substrate binding, stereoselectivity, and intermediate stabilization during the reaction. In particular, we have confirmed the role of the oxyanion hole and the conserved motif (HGXP) in epoxide hydrolases, in excellent agreement with known experimental and computational data on similar systems. The model obtained is the first one that fully agrees with all the experimental observations on the system. Proteins 2017; 85:720–730. © 2016 Wiley Periodicals, Inc.     

CCR5 antagonist TD-0680 uses a novel mechanism for enhanced potency against HIV-1 entry, cell-mediated infection, and a resistant variant

Regardless of the route of transmission, R5-tropic HIV-1 predominates early in infection, rendering C-C chemokine receptor type 5 (CCR5) antagonists as attractive agents not only for antiretroviral therapy but also for prevention. Here, we report the specificity, potency, and underlying mechanism of action of a novel small molecule CCR5 antagonist, TD-0680. TD-0680 displayed the greatest potency against a diverse group of R5-tropic HIV-1 and SIV strains when compared with its prodrug, TD-0232, the Food and Drug Administration-approved CCR5 antagonist Maraviroc, and TAK-779, with EC(50) values in the subnanomolar range (0.09-2.29 nm). Importantly, TD-0680 was equally potent at blocking envelope-mediated cell-cell fusion and cell-mediated viral transmission as well as the replication of a TAK-779/Maraviroc-resistant HIV-1 variant. Interestingly, TD-0232 and TD-0680 functioned differently despite binding to a similar transmembrane pocket of CCR5. Site-directed mutagenesis, drug combination, and antibody blocking assays identified a novel mechanism of action of TD-0680. In addition to binding to the transmembrane pocket, the unique exo configuration of this molecule protrudes and sterically blocks access to the extracellular loop 2 (ECL2) region of CCR5, thereby interrupting the interaction between virus and its co-receptor more effectively. This mechanism of action was supported by the observations of similar TD-0680 potency against CD4-dependent and -independent SIV strains and by molecular docking analysis using a CCR5 model. TD-0680, therefore, merits development as an anti-HIV-1 agent for therapeutic purposes and/or as a topical microbicide for the prevention of sexual transmission of R5-tropic HIV-1.     

Structure prediction of GPCRs using piecewise homologs and application to the human CCR5 chemokine receptor: validation through agonist and antagonist docking

This article describes the construction and validation of a three-dimensional model of the human CC chemokine receptor 5 (CCR5) receptor using multiple homology modeling. A new methodology is presented where we built each secondary structural model of the protein separately from distantly related homologs of known structure. The reliability of our approach for G-protein coupled receptors was assessed through the building of the human C-X-C chemokine receptor type 4 (CXCR4) receptor of known crystal structure. The models are refined using molecular dynamics simulations and energy minimizations using CHARMM, a classical force field for proteins. Finally, docking models of both the natural agonists and the antagonists of the receptors CCR5 and CXCR4 are proposed. This study explores the possible binding process of ligands to the receptor cavity of chemokine receptors at molecular and atomic levels. We proposed few crucial residues in receptors binding to agonist/antagonist for further validation through experimental analysis. In particular, our study provides better understanding of the blockage mechanism of the chemokine receptors CCR5 and CXCR4, and may help the identification of new lead compounds for drug development in HIV infection, inflammatory diseases, and cancer metastasis.     

Structure prediction and molecular dynamics simulations of a G-protein coupled receptor: human CCR2 receptor

CC chemokine receptor type-2 (CCR2) is a member of G-protein coupled receptors superfamily, expressed on the cell surface of monocytes and macrophages. It binds to the monocyte chemoattractant protein-1, a CC chemokine, produced at the sites of inflammation and infection. A homology model of human CCR2 receptor based on the recently available C-X-C chemokine recepor-4 crystal structure has been reported. Ligand information was used as an essential element in the homology modeling process. Six known CCR2 antagonists were docked into the model using simple and induced fit docking procedure. Docked complexes were then subjected to visual inspection to check their suitability to explain the experimental data obtained from site directed mutagenesis and structure-activity relationship studies. The homology model was refined, validated, and assessed for its performance in docking-based virtual screening on a set of CCR2 antagonists and decoys. The docked complexes of CCR2 with the known antagonists, TAK779, a dual CCR2/CCR5 antagonist, and Teijin-comp1, a CCR2 specific antagonist were subjected to molecular dynamics (MD) simulations, which further validated the binding modes of these antagonists. B-factor analysis of 20?ns MD simulations demonstrated that Cys190 is helpful in providing structural rigidity to the extracellular loop (EL2). Residues important for CCR2 antagonism were recognized using free energy decomposition studies. The acidic residue Glu291 from TM7, a conserved residue in chemokine receptors, is favorable for the binding of Teijin-comp1 with CCR2 by ΔG of ?11.4?kcal/mol. Its contribution arises more from the side chains than the backbone atoms. In addition, Tyr193 from EL2 contributes ?0.9?kcal/mol towards the binding of the CCR2 specific antagonist with the receptor. Here, the homology modeling and subsequent molecular modeling studies proved successful in probing the structure of human CCR2 chemokine receptor for the structure-based virtual screening and predicting the binding modes of CCR2 antagonists.     

Binding site analysis of full-length alpha1a adrenergic receptor using homology modeling and molecular docking

The recent availability of crystal structure of bovine rhodopsin offers new opportunities in order to approach the construction of G protein coupled receptors. This study focuses the attention on the modeling of full-length alpha(1a) adrenergic receptor (alpha(1a)-AR) due to its biological role and significant implications in pharmacological treatment of benign prostate hyperplasia. This work could be considered made up by two main steps: (a) the construction of full structure of alpha(1a)-AR, through homology modeling methods; (b) the automated docking of an endogenous agonist, norepinephrine, and of an antagonist, WB-4101, using BioDock program. The obtained results highlight the key residues involved in binding sites of both agonists and antagonists, confirming the mutagenesis data and giving new suggestions for the rational design of selective ligands.     

Screening of promising molecules against MurG as drug target in multi-drug-resistant-Acinetobacter baumannii - insights from comparative protein modeling,molecular docking and molecular dynamics simulation

Abstract

The UDP-N-acetylglucosamine-N-acetylmuramyl-(pentapeptide) pyrophosphoryl-undecaprenol N-acetylglucosamine transferase (MurG) is located in plasma membrane which plays a crucial role for peptidoglycan biosynthesis in Gram-negative bacteria. Recently, this protein is considered as an important and unique drug target in Acinetobacter baumannii since it plays a key role during the synthesis of peptidoglycan as well as which is not found in Homo sapiens. In this study, initially we performed comparative protein modeling approach to predict the three-dimensional model of MurG based on crystal structure of UDP-N-acetylglucosamine-N-acetylmuramyl-(pentapeptide) pyrophosphoryl-undecaprenol N-acetylglucosamine transferase (PDB ID: 1F0K) from E.coli K12. MurG model has two important functional domains located in N and C- terminus which are separated by a deep cleft. Active site residues are located between two domains and they are Gly20, Arg170, Gly200, Ser201, Gln227, Phe254, Leu275, Thr276, and Glu279 which play essential role for the function of MurG. In order to inhibit the function of MurG, we employed the High Throughput Virtual Screening (HTVS) and docking techniques to identify the promising molecules which will further subjected into screening for computing their drug like and pharmacokinetic properties. From the HTVS, we identified 5279 molecules, among these, 12 were passed the drug-like and pharmacokinetic screening analysis. Based on the interaction analysis in terms of binding affinity, inhibition constant and intermolecular interactions, we selected four molecules for further MD simulation to understand the structural stability of protein-ligand complexes. All the analysis of MD simulation suggested that ZINC09186673 and ZINC09956120 are identified as most promising putative inhibitors for MurG protein in A. baumannii.

Communicated by Ramaswamy H. Sarma     

Prediction of HIV-1 entry inhibitors neomycin-arginine conjugates interaction with the CD4-gp120 binding site by molecular modeling and multistep docking procedure

Developing of multi-target HIV-1 entry inhibitors represents an important avenue of drug therapy. Two such inhibitors are hexa-arginine-neomycin-conjugate (NeoR6) and nona-d-arginine-neomycin-conjugate (Neo-r9). Our findings that NeoR6-resistant mutations appear in the gp120 constant regions; and NeoR6 is not CCR5 antagonist, but inhibits CXCR4 and CCR5 HIV-1 using isolates, led us to suggest that NeoR6 may inhibit HIV-1 entry by interfering with the CD4-gp120 binding. To support this notion, we constructed a homology model of unliganded HIV-1_IIIB gp120 and docked NeoR6 and Neo-r9 to it, using a multistep docking procedure: geometric-electrostatic docking by MolFit; flexible ligand docking by Autodock3 and final refinement of the obtained complexes by Discover3. Binding free energies were calculated by MM-PBSA methodology. The model predicts competitive inhibition of CD4-gp120 binding by NeoR6 and Neo-r9. We determined plausible binding sites between constructed CD4-bound gp120 trimer and homology modeled membranal CXCR4, and tested NeoR6 and Neo-r9 interfering with this interaction. These models support our notion that another mechanism of anti-HIV-1 activity of NeoR6 is inhibition of gp120-CXCR4 binding. These structural models and interaction of NeoR6 and Neo-r9 with gp120 and CXCR4 provide a powerful approach for structural based drug design for selective targeting of HIV-1 entry and/or for inhibition of other retroviruses with similar mechanism of entry.     

Restricted variable residues in the C-terminal segment of HIV-1 V3 loop regulate the molecular anatomy of CCR5 utilization

The V3 loop of the HIV-1 envelope glycoprotein (Env) is the major determinant for coreceptor utilization, but the structural basis for this specificity remains to be defined. By characterizing a set of naturally occurring R5 Env variants, we demonstrate that Asp324 in the conserved IIGDIR motif of the V3 loop (CTRPN(300)NNTRKSIHIGP(311)GRAFYTTGEIIGD(324)IRQAHC) C-terminal segment regulates the molecular anatomy of CCR5 utilization. Whereas gp120 subunits with Asp or Asn at position 324 were fusogenic with coreceptor chimeras containing either the N-terminal domain or the body of CCR5, substitution of charged (Glu, Lys) or small hydrophobic (Gly, Ala) residues resulted in complete loss of fusogenic activity with the N terminus and markedly reduced utilization of the body of CCR5, although their ability to use wild-type CCR5 was unchanged. This phenotypic conversion was confirmed in both gain and loss of function experiments using Env from multiple subtypes. Alignment of sequences of R5 V3 loops (n=599) from the HIV database revealed that the mutation of Asp324 in the conserved IIGDIR motif is restricted to Asn324, with proportions of 71.5% and 28%, respectively. Infection of primary CD4(+)T cells demonstrated that Env bearing Asp324 was less sensitive to RANTES, suggesting that Asp or Asn in this position may be crucial for viral fitness. The CD4-dependent gp120 binding to CCR5 was decreased when Asp324 was replaced with a charged or hydrophobic residue, but unchanged when replaced with Asn. Molecular modeling analyses predicted that Asp/Asn324 forms a critical H-bond with Asn300. These findings indicate that Asp or Asn at position 324 of the V3 stem stabilizes the conformation of V3 loop and hence influences the intensities of interaction between CD4-activated gp120 and CCR5 which results in viral entry.     

Targeting Spare CC Chemokine Receptor 5 (CCR5) as a Principle to Inhibit HIV-1 Entry

CCR5 binds the chemokines CCL3, CCL4, and CCL5 and is the major coreceptor for HIV-1 entry into target cells. Chemokines are supposed to form a natural barrier against human immunodeficiency virus, type 1 (HIV-1) infection. However, we showed that their antiviral activity is limited by CCR5 adopting low-chemokine affinity conformations at the cell surface. Here, we investigated whether a pool of CCR5 that is not stabilized by chemokines could represent a target for inhibiting HIV infection. We exploited the characteristics of the chemokine analog PSC-RANTES (N-α-(n-nonanoyl)-des-Ser(1)-[l-thioprolyl(2), l-cyclohexylglycyl(3)]-RANTES(4-68)), which displays potent anti-HIV-1 activity. We show that native chemokines fail to prevent high-affinity binding of PSC-RANTES, analog-mediated calcium release (in desensitization assays), and analog-mediated CCR5 internalization. These results indicate that a pool of spare CCR5 may bind PSC-RANTES but not native chemokines. Improved recognition of CCR5 by PSC-RANTES may explain why the analog promotes higher amounts of β-arrestin 2·CCR5 complexes, thereby increasing CCR5 down-regulation and HIV-1 inhibition. Together, these results highlight that spare CCR5, which might permit HIV-1 to escape from chemokines, should be targeted for efficient viral blockade.     

Computational insights into the binding pattern of mitochondrial calcium uniporter inhibitor through homology modeling,molecular dynamics simulation,binding free energy prediction and density functional theory calculation

Abstract

The mitochondrial calcium uniporter (MCU) is the critical protein of the inner mitochondrial membrane that is the primary mediator for calcium uptake into the mitochondrial matrix. Herein we built the optimal homology model of human MCU which was refined through all-atom molecular dynamics simulation. Then, the binding mode of known inhibitor was predicted through molecular docking method, along with molecular dynamics simulation and binding free energy calculation to verify the docking result and stability of the protein-inhibitor complex. Finally, density functional theory (DFT) calculation enhanced our understanding of the molecular interaction of MCU inhibitor. Our research would provide a deeper insight into the interactions between human MCU and its inhibitor, which boosts to develop novel therapy against MCU related disease.

Communicated by Ramaswamy H. Sarma     

Role of V3 independent domains on a dualtropic human immunodeficiency virus type 1 (HIV-1) envelope gp120 in CCR5 coreceptor utilization and viral infectivity

The molecular mechanism of human immunodeficiency virus type 1 (HIV-1) entry into cells involves specific interactions between the viral envelope glycoprotein gp120 and two target cell proteins, CD4 and either CCR5 or CXCR4 chemokine receptors. In order to delineate the functional role of HIV-1 gp120 subdomains of dualtropic strains in CCR5 coreceptor usage, we used a panel of chimeric viruses in which the V1/V2 and V3 domains of gp120 from the dualtropic HIV-1(KMT) isolate were introduced either alone or in combination into the T-tropic HIV-1(NL4-3) background. These chimeric constructs were employed in cell-cell fusion and cell-free virus infectivity assays using cell lines expressing CD4 and the CCR5 chemokine receptor. In both assays, the V3 domain of HIV-1(KMT) but not the V1/V2 domain proved to be the principal determinant of CCR5 coreceptor usage. However, in the cell-free viral infectivity assay although a chimeric virus with a combined V1/V2 and V3 domains of HIV-1(KMT) efficiently fused with coreceptor expressing cells, yet its infectivity was markedly diminished in CCR5 as well as CXCR4 expressing cells. Restoring a comparable level of infection of such chimeric virus required the C3-V5 domain from HIV-1(KMT) to be introduced. Our present findings confirmed that the V3 domain is the major determinant of fusion activity and cellular tropism, and demonstrated a dispensable role for the V1/V2 domain. In addition the C3-V5 domain appeared to play an important role in viral infectivity when the corresponding V1/V2 and V3 domains are present.     

Small-angle X-ray scattering based structure,modeling and molecular dynamics analyses of a family 5 glycoside hydrolase first endo-mannanase named as RfGH5_7 from Ruminococcus flavefaciens

Abstract

Endo-β-1,4-mannanase named as RfGH5_7 from Ruminococcus flavefaciens cloned, expressed and purified earlier was structurally characterized in present study. The RaptorX modeled structure of RfGH5_7 showed a (β/α)₈ Triose-phosphate Isomerase (TIM) barrel fold. The Ramachandran plot assessment of RfGH5_7 showed that all amino acids fall in allowed region except one, Asn22 in the disallowed region. The superposition of RfGH5_7 modeled structure with its nearest homologues revealed that Glu154 acts as proton donor while Glu249 acts as nucleophile. Secondary structure of RfGH5_7 through Circular Dichroism (CD) analysis revealed 33.5% α-helices, 17% β-strands and 49.5% random coils. Molecular Dynamic (MD) simulation showed Root Mean Square Deviation (RMSD), 0.67?nm and radius of gyration (Rg) between 1.9?nm and 1.85?nm. The binding interaction of mannotetraose on the surface of RfGH5_7 structure displayed polar interactions with His219, Tyr221, Trp278, Ser279 and Gly282 residues. Small-angle X-ray scattering (SAXS) analysis displayed the intact and monodispersed nature of the enzyme RfGH5_7. The radius of gyration (Rg) by Guinier analysis for globular shape was found to be 2.29?±?0.09?nm and for rod-shape it was 0.95?±?0.02?nm. Kratky plot confirmed that RfGH5_7 structure is compact and folded in solution. The ab initio derived dummy model of RfGH5_7 displayed single domain structure of yellow humped fish like shape. The RfGH5_7 modeled structure was well fitted with ab initio derived model from SAXS data.

Communicated by Ramaswamy H. Sarma     

An allosteric rheostat in HIV-1 gp120 reduces CCR5 stoichiometry required for membrane fusion and overcomes diverse entry limitations

Binding of the human immunodeficiency virus (HIV-1) envelope glycoprotein gp120 to the CCR5 co-receptor reduces constraints on the metastable transmembrane subunit gp41, thereby enabling gp41 refolding, fusion of viral and cellular membranes, and infection. We previously isolated adapted HIV-1_JRCSF variants that more efficiently use mutant CCR5s, including CCR5(Δ18) lacking the important tyrosine sulfate-containing amino terminus. Effects of mutant CCR5 concentrations on HIV-1 infectivities were highly cooperative, implying that several may be required. However, because wild-type CCR5 efficiently mediates infections at trace concentrations that were difficult to measure accurately, analyses of its cooperativity were not feasible. New HIV-1_JRCSF variants efficiently use CCR5(HHMH), a chimera containing murine extracellular loop 2. The adapted virus induces large syncytia in cells containing either wild-type or mutant CCR5s and has multiple gp120 mutations that occurred independently in CCR5(Δ18)-adapted virus. Accordingly, these variants interchangeably use CCR5(HHMH) or CCR5(Δ18). Additional analyses strongly support a novel energetic model for allosteric proteins, implying that the adaptive mutations reduce quaternary constraints holding gp41, thus lowering the activation energy barrier for membrane fusion without affecting bonds to specific CCR5 sites. In accordance with this mechanism, highly adapted HIV-1s require only one associated CCR5(HHMH), whereas poorly adapted viruses require several. However, because they are allosteric ensembles, complexes with additional co-receptors fuse more rapidly and efficiently than minimal ones. Similarly, wild-type HIV-1_JRCSF is highly adapted to wild-type CCR5 and minimally requires one. The adaptive mutations cause resistances to diverse entry inhibitors and cluster appropriately in the gp120 trimer interface overlying gp41. We conclude that membrane fusion complexes are allosteric machines with an ensemble of compositions, and that HIV-1 adapts to entry limitations by gp120 mutations that reduce its allosteric hold on gp41. These results provide an important foundation for understanding the mechanisms that control membrane fusion and HIV-1's facile adaptability.     

In silico study of 1-(4-Phenylpiperazin-1-yl)-2-(1H-pyrazol-1-yl) ethanones derivatives as CCR1 antagonist: Homology modeling,docking and 3D-QSAR approach

C–C chemokine receptor type 1 (CCR1) is a chemokine receptor with seven transmembrane helices and it belongs to the G-Protein Coupled receptor (GPCR) family. It plays an important role in rheumatoid arthritis, organ transplant rejection, Alzheimer’s disease and also causes inflammation. Because of its role in disease processes, CCR1 is considered to be an important drug target. In the present study, we have performed three dimensional Quantitative Structure activity relationship (3D-QSAR) studies on a series of 1-(4-Phenylpiperazin-1-yl)-2-(1H-pyrazol-1-yl) ethanone derivatives targeting CCR1. Homology modeling of CCR1 was performed based on a template structure (4EA3) which has a high sequence identity and resolution. The highest active molecule was docked into this model. Ligand-based and Receptor-based quantitative structure–activity relationship (QSAR) study was performed and CoMFA models with reasonable statistics was developed for both ligand-based (q² = 0.606; r² = 0.968) and receptor-guided (q² = 0.640; r² = 0.932) alignment methods. Contour map analyses identified favorable regions for high affinity binding. The docking results highlighted the important active site residues. Tyr113 was found to interact with the ligand through hydrogen bonding. This residue has been considered responsible for anchoring ligands inside the active site. Our results could also be helpful to understand the inhibitory mechanism of 1-(4-Phenylpiperazin-1-yl)-2-(1H-pyrazol-1-yl) ethanone derivatives thereby to design more effective ligands in the future.