Computational study on the interaction between CCR5 and HIV-1 entry inhibitor maraviroc: insight from accelerated molecular dynamics simulation and free energy calculation

Structural dynamics of chemokine receptors

Membrane proteins such as G protein-coupled receptors (GPCRs) are involved in awide range of physiological and pathological cellular processes. Binding of extracellular signals to GPCRs, including hormones, neurotransmitters, peptides and proteins, can activate intracellular signaling cascades via G protein interaction. Chemokine receptors are key GPCRs implicated in cancers, immune responses, cell migration and inflammation. Specifically, the CCR5 and CXCR4 chemokine receptors serve as important therapeutic targets against Human Immunodeficiency virus (HIV) entry into human cells. Maraviroc and Vicriviroc, two clinically used HIV entry inhibitors, are antagonists of the CCR5 receptor. These drugs block HIV entry, but ultimately resistance develops, due to emergence of viruses that can utilize the CXCR4 co-receptor. Unfortunately, development of chemokine receptor antagonists as selective drugs of HIV infection has been greatly hindered as their target orthosteric site is conserved among different receptor subtypes. Accordingly, it is important to understand the structural dynamics of these receptors to develop more effective therapeutics. In this chapter, we describe the latest advances in studies of these two key chemokine receptors with respect to their structures, dynamics and function.

Computational study of the structural ensemble of CC chemokine receptor type 5 (CCR5) and its interactions with different ligands

CC Chemokine receptor 5 (CCR5), a member of the Superfamily of G Protein-Coupled Receptors (GPCRs), is an important effector in multiple physiopathological processes such as inflammatory and infectious entities, including central nervous system neuroinflammatory diseases such as Alzheimer's disease, recovery from nervous injuries, and in the HIV-AIDS infective processes. Thus, CCR5 is an attractive target for pharmacological modulation. Since maraviroc was described as a CCR5 ligand that modifies the HIV-AIDS progression, multiple efforts have been developed to describe the functionality of the receptor. In this work, we characterized key structural features of the CCR5 receptor employing extensive atomistic molecular dynamics (MD) in its apo form and in complex with an endogenous agonist, the chemokine CCL5/RANTES, an HIV entry inhibitor, the partial inverse agonist maraviroc, and the experimental antagonists Compound 21 and 34, aiming to elucidate the structural features and mechanistic processes that constitute its functional states, contributing with structural details and a general understanding of this relevant system.

Towards Selective Binding to the GLUT5 Transporter: Synthesis, Molecular Dynamics and In Vitro Evaluation of Novel C-3-Modified 2,5-Anhydro-D-mannitol Analogs

Deregulation and changes in energy metabolism are emergent and important biomarkers of cancer cells. The uptake of hexoses in cancer cells is mediated by a family of facilitative hexose membrane-transporter proteins known as Glucose Transporters (GLUTs). In the clinic, numerous breast cancers do not show elevated glucose metabolism (which is mediated mainly through the GLUT1 transporter) and may use fructose as an alternative energy source. The principal fructose transporter in most cancer cells is GLUT5, and its mRNA was shown to be elevated in human breast cancer. This offers an alternative strategy for early detection using fructose analogs. In order to selectively scout GLUT5 binding-pocket requirements, we designed, synthesized and screened a new class of fructose mimics based upon the 2,5-anhydromannitol scaffold. Several of these compounds display low millimolar IC50 values against the known high-affinity 18F-labeled fructose-based probe 6-deoxy-6-fluoro-D-fructose (6-FDF) in murine EMT6 breast cancer cells. In addition, this work used molecular docking and molecular dynamics simulations (MD) with previously reported GLUT5 structures to gain better insight into hexose–GLUT interactions with selected ligands governing their preference for GLUT5 compared to other GLUTs. The improved inhibition of these compounds, and the refined model for their binding, set the stage for the development of high-affinity molecular imaging probes targeting cancers that express the GLUT5 biomarker.

Conformation Transition of Intracellular Part of Glucagon Receptor in Complex With Agonist Glucagon by Conventional and Accelerated Molecular Dynamics Simulations

The inactive conformations of glucagon receptor (GCGR) are widely reported by crystal structures that support the precision structure for drug discovery of type 2 diabetes. The previous study shows that the intracellular part is open in the glucagon-bound GCGR (glu-GCGR) and closed in the apo-GCGR by accelerated molecular dynamics (aMD) simulations. However, the crystal structure of GCGR in complex with partial agonist shows that the intracellular part is closed in the inactive conformation. To understand the differences between the studies of aMD simulations and crystal structure, the 2,500 ns conventional molecular dynamics (cMD) simulations are performed on the simulated model of glu-GCGR. The result shows that the transmembrane helices (TMH) 6 of glu-GCGR is outward ~4 Å to drive the intracellular part of glu-GCGR open until ~390 ns cMD simulations. The (TMH) 6 of glu-GCGR becomes closed after ~490 ns cMD simulations, which are consistent with the crystal structure of GCGR in complex with the partial agonist. To further elucidate the activation mechanism of GCGR deeply, the simulated models of glu-GCGR, apo-GCGR, and antagonist-bound GCGR (ant-GCGR) are constructed to perform 10 of parallel 300 ns aMD simulations, respectively. The results show that both of glu-GCGR and apo-GCGR can generate the open conformations of the intracellular part. But the glu-GCGR has the higher percentage of open conformations than apo-GCGR. The ant-GCGR is restricted to generate the open conformations of the intracellular part by antagonist MK-0893. It indicates that the glu-GCGR, apo-GCGR, and ant-GCGR can be distinguished by the aMD simulated method. Free energy landscape shows that the open conformations of the intracellular part of GCGR are in intermediate state. Our results show that aMD simulations enhance the space samplings of open conformations of GCGR via adding extra boost potential. It indicates that the aMD simulations are an effective way for drug discovery of GCGR.

Probing the Druggablility on the Interface of the Protein-Protein Interaction and Its Allosteric Regulation Mechanism on the Drug Screening for the CXCR4 Homodimer

Modulating protein–protein interactions (PPIs) with small drug-like molecules targeting it exhibits great promise in modern drug discovery. G protein-coupled receptors (GPCRs) are the largest family of targeted proteins and could form dimers in living biological cells through PPIs. However, compared to drug development of the orthosteric site, there has been lack of investigations on the druggability of the PPI interface for GPCRs and its functional implication on experiments. Thus, in order to address these issues, we constructed a novel computational strategy, which involved in molecular dynamics simulation, virtual screening and protein structure network (PSN), to study one representative GPCR homodimer (CXCR4). One druggable pocket was identified in the PPI interface and one small molecule targeting it was screened, which could strengthen PPI mainly through hydrophobic interaction between the benzene rings of the PPI molecule and TM4 of the receptor. The PSN results further reveals that the PPI molecule could increase the number of the allosteric regulation pathways between the druggable pocket of the dimer interface to the orthostatic site for the subunit A but only play minor role for the other subunit B, leading to the asymmetric change in the volume of the binding pockets for the two subunits (increase for the subunit A and minor change for the subunit B). Consequently, the screening performance of the subunit A to the antagonists is enhanced while the subunit B is unchanged nearly, implying that the PPI molecule may be beneficial to enhance the drug efficacies of the antagonists. In addition, one main regulation pathway with the highest frequency was identified for the subunit A, which consists of Trp195^5.34–Tyr190^ECL2–Val196^5.35–Gln200^5.39–Asp262^6.58–Cys28^N-term, revealing their importance in the allosteric regulation from the PPI molecule. The observations from the work could provide valuable information for the development of the PPI drug-like molecule for GPCRs.

Computational approach towards understanding structural and functional role of cytokinin oxidase/dehydrogenase 2 (CKX2) in enhancing grain yield in rice plant

Cytokinin oxidase/dehydrogenase (CKX) is the only known enzyme associated with irreversible degradation of cytokinins in plants. CKX2 contains flavin adenine dinucleotide (FAD) domain. Earlier studies utilising antisense & hpRNAi suppression techniques in mutant/transgenic rice plants revealed that when CKX2 binds with FAD, CKX2 expression reduces, which in turn causes cytokinin aggregation in inflorescence meristem that subsequently enhances both branches and grain number resulting in increased grain yield. Owing to the non-existence of complete three-dimensional structure of CKX2, insight into the structure and function of CKX2 and its relationship with its cofactor FAD is still a topic of debate. In the present study, computational approach was employed to estimate the three-dimensional structure of CKX2 through comparative modelling approach. Later, CKX2 and FAD interaction study was performed to understand the underlying mechanism involved with reduced expression of CKX2. Molecular dynamic simulation studies of both CKX2 and CKX-FAD complex revealed that after binding with FAD, CKX2 experienced increased pressure and reduced RMSD, potential energy and free energy landscape energy, which in turn lessen anti-correlation between almost all α and β strands and random motion of C-α, subsequently reducing CKX2 expression. In near future, these information can be utilised for increasing rice yield under irrigated field condition by introgression of Gn1a gene through marker assisted back-crossing breeding. AbbreviationsGROMACSGROningen MAchine for Chemical SimulationsNPTConstant Number of Particles, Volume and TemperatureRMSDRoot Mean Square DeviationRMSFRoot Mean Square FluctuationsQTLquantitative trait lociFADflavin adenine dinucleotideNVTConstant Number of Particles, Pressure and TemperatureLINCSLinear Constraint SolverCKX2Cytokinin oxidase/dehydrogenase 2MM/PBSAMolecular Mechanics/Poisson-Boltzmann surface areaSDFStructure Data FileCommunicated by Ramaswamy H. Sarma.

Use multiscale simulation to explore the effects of the homodimerizations between different conformation states on the activation and allosteric pathway for the μ-opioid receptor

Recently, oligomers of G-protein coupled receptors (GPCRs) have been an important topic in the GPCR fields. However, knowledge about their structures and activation mechanisms is very limited due to the absence of crystal structures reported. In this work, we used multiscale simulations to study the effects of homodimerization between different conformation states on their activation, dynamic behaviors, and allosteric communication pathways for μ-OR. The results indicated that the dimerization of one inactive monomer with either one inactive monomer or one active one could enhance its constitutive activation. However, the conformation state of the other protomer (e.g., active or inactive) can influence the activated extent. The dimerization between the two inactive protomers leads to a negative cooperativity for their activation, which should contribute to the asymmetric activation of GPCR dimers observed in some experiments. On the other hand, for the active monomer, its dimerization with one inactive receptor could alleviate its deactivation, whereby negative and positive cooperativities can be observed between the two subunits of the dimer, depending on the different regions. Observations from protein structure network (PSN) analysis indicated that the dimerization of one inactive monomer with one active one would cause a significant drop in the number of main pathways from the ligand binding pocket to the G-protein coupled region for the inactive protomer, while the impact is minor for the active protomer. But, for the active monomer or the inactive one, its dimerization with one inactive monomer would significantly change the types of residues participating in the pathway with the highest frequency.

Deciphering the binding behavior of flavonoids to the cyclin dependent kinase 6/cyclin D complex

Flavonoids, a class of natural compounds with variable phenolic structures, have been found to possess anti-cancer activities by modulating different enzymes and receptors like CDK6. To understand the binding behavior of flavonoids that inhibit the active CDK6, molecular dynamics (MD) simulations were performed on six inhibitors, chrysin (M01), fisetin (M03), galangin (M04), genistein (M05), quercetin (M06) and kaempferol (M07), complexed with CDK6/cyclin D. For all six flavonoids, the 3'-OH and 4'-OH of B-ring were found to be favorable for hydrogen bond formation, but the 3-OH on the C-ring and 5-OH on the A-ring were unfavorable, which were confirmed by the MD simulation results of the test molecule, 3', 4', 7-trihydroxyflavone (M15). The binding efficiencies of flavonoids against the CDK6/cyclin D complex were mainly through the electrostatic (especially the H-bond force) and vdW interactions with residues ILE19, VAL27, ALA41, GLU61, PHE98, GLN103, ASP163 and LEU152. The order of binding affinities of these flavonoids toward the CDK6/cyclin D was M03 > M01 > M07 > M15 > M06 > M05 > M04. It is anticipated that the binding features of flavonoid inhibitors studied in the present work may provide valuable insights for the development of CDK6 inhibitors.

Computational studies on horseshoe shape pocket of human orexin receptor type 2 and boat conformation of suvorexant by molecular dynamics simulations

The FDA approved drug suvorexant binds to the horseshoe shape pocket of OX₂ R with the boat conformation. The horseshoe shape pocket plays an important role on the biological activity of OX₂ R in the cell membrane. To study the binding mechanism between the horseshoe shape pocket of OX₂ R and boat conformation of suvorexant, the crystal structures of wild type and N324A mutant of OX₂ R in complex with antagonist suvorexant are chosen to perform molecular dynamics (MD) simulations, QM/MM, and MMGBSA calculations. By comparison with the wild type of OX₂ R, the results show the 1,2,3-triazole and p-toluamide groups of suvorexant are changed in the N324A mutant of OX₂ R during 200 ns MD simulations. The QM/MM and weak interaction analysis are employed to calculate the non-covalent bonds interaction between suvorexant and key residues in the wild type and N324A mutant of OX₂ R. The MMGBSA calculations indicate the entropy energy is an important influence factor for suvorexant affinity in the distorted horseshoe shape pocket of OX₂ R. Our results not only show the horseshoe shape pocket of OX₂ R is the necessary conformation for the binding of antagonist suvorexant, but also give the important sites and structural features for antagonist design of OX₂ R.

Importance of protein flexibility on molecular recognition: modeling binding mechanisms of aminopyrazine inhibitors to Nek2

NIMA-related kinase 2 (Nek2) plays a significant role in cell cycle regulation, and overexpression of Nek2 has been observed in several types of carcinoma, suggesting it is a potential target for cancer therapy. In this study, we attempted to gain more insight into the binding mechanisms of a series of aminopyrazine inhibitors of Nek2 through multiple molecular modeling techniques, including molecular docking, molecular dynamics (MD) simulations and free energy calculations. The simulation results showed that the induced fit docking and ensemble docking based on multiple protein structures yield better predictions than conventional rigid receptor docking, highlighting the importance of incorporating receptor flexibility into the accurate predictions of the binding poses and binding affinities of Nek2 inhibitors. Additionally, we observed that the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) calculations did not show better performance than the docking scoring to rank the binding affinities of the studied inhibitors, suggesting that MM/GBSA is system-dependent and may not be the best choice for the Nek2 systems. Moreover, the detailed information on protein-ligand binding was characterized by the MM/GBSA free energy decomposition, and a number of derivatives with improved docking scores were designed. It is expected that our study can provide valuable information for the future rational design of novel and potent inhibitors of Nek2.

Binding Kinetics and Pathways of Ligands to GPCRs

Previously, drugs were developed focusing on target affinity and selectivity. However, it is becoming evident that the drug-target residence time, related to the off-rate, is an important parameter for successful drug development. The residence time influences both the on-rate and overall effectiveness of drugs. Furthermore, ligand binding is now appreciated to be a multistep process because metastable and/or intermediate binding sites in the extracellular region have been identified. In this review, we summarize experimental ligand-binding data for G-protein-coupled receptors (GPCRs), and their binding pathways, analyzed by molecular dynamics (MD). The kinetics of drug binding to GPCRs are complex and depend on several factors, including charge distribution on the receptor surface, ligand-receptor interactions in the binding channel and the binding site, or solvation.

Frequency and Env determinants of HIV-1 subtype C strains from antiretroviral therapy-naive subjects that display incomplete inhibition by maraviroc

BACKGROUND

Entry of human immunodeficiency virus type 1 (HIV-1) into cells involves the interaction of the viral gp120 envelope glycoproteins (Env) with cellular CD4 and a secondary coreceptor, which is typically one of the chemokine receptors CCR5 or CXCR4. CCR5-using (R5) HIV-1 strains that display reduced sensitivity to CCR5 antagonists can use antagonist-bound CCR5 for entry. In this study, we investigated whether naturally occurring gp120 alterations in HIV-1 subtype C (C-HIV) variants exist in antiretroviral therapy (ART)-naïve subjects that may influence their sensitivity to the CCR5 antagonist maraviroc (MVC).

RESULTS

Using a longitudinal panel of 244 R5 Envs cloned from 20 ART-naïve subjects with progressive C-HIV infection, we show that 40% of subjects (n = 8) harbored viruses that displayed incomplete inhibition by MVC, as shown by plateau's of reduced maximal percent inhibitions (MPIs). Specifically, when pseudotyped onto luciferase reporter viruses, 16 Envs exhibited MPIs below 98% in NP2-CCR5 cells (range 79.7-97.3%), which were lower still in 293-Affinofile cells that were engineered to express high levels of CCR5 (range 15.8-72.5%). We further show that Envs exhibiting reduced MPIs to MVC utilized MVC-bound CCR5 less efficiently than MVC-free CCR5, which is consistent with the mechanism of resistance to CCR5 antagonists that can occur in patients failing therapy. Mutagenesis studies identified strain-specific mutations in the gp120 V3 loop that contributed to reduced MPIs to MVC.

CONCLUSIONS

The results of our study suggest that some ART-naïve subjects with C-HIV infection harbor HIV-1 with reduced MPIs to MVC, and demonstrate that the gp120 V3 loop region contributes to this phenotype.

Localization of Millisecond Dynamics: Dihedral Entropy from Accelerated MD

Here, we demonstrate a method to capture local dynamics on a time scale 3 orders of magnitude beyond state-of-the-art simulation approaches. We apply accelerated molecular dynamics simulations for conformational sampling and extract reweighted backbone dihedral distributions. Local dynamics are characterized by torsional probabilities, resulting in residue-wise dihedral entropies. Our approach is successfully validated for three different protein systems of increasing size: alanine dipeptide, bovine pancreatic trypsin inhibitor (BPTI), and the major birch pollen allergen Bet v 1a. We demonstrate excellent agreement of flexibility profiles with both large-scale computer simulations and NMR experiments. Thus, our method provides efficient access to local protein dynamics on extended time scales of high biological relevance.

Investigation of allosteric modulation mechanism of metabotropic glutamate receptor 1 by molecular dynamics simulations, free energy and weak interaction analysis

Metabotropic glutamate receptor 1 (mGlu₁), which belongs to class C G protein-coupled receptors (GPCRs), can be coupled with G protein to transfer extracellular signal by dimerization and allosteric regulation. Unraveling the dimer packing and allosteric mechanism can be of great help for understanding specific regulatory mechanism and designing more potential negative allosteric modulator (NAM). Here, we report molecular dynamics simulation studies of the modulation mechanism of FITM on the wild type, T815M and Y805A mutants of mGlu₁ through weak interaction analysis and free energy calculation. The weak interaction analysis demonstrates that van der Waals (vdW) and hydrogen bonding play an important role on the dimer packing between six cholesterol molecules and mGlu₁ as well as the interaction between allosteric sites T815, Y805 and FITM in wild type, T815M and Y805A mutants of mGlu₁. Besides, the results of free energy calculations indicate that secondary binding pocket is mainly formed by the residues Thr748, Cys746, Lys811 and Ser735 except for FITM-bound pocket in crystal structure. Our results can not only reveal the dimer packing and allosteric regulation mechanism, but also can supply useful information for the design of potential NAM of mGlu₁.

Computational insights into different inhibition modes of the κ-opioid receptor with antagonists LY2456302 and JDTic

Residence time calculations were carried out based on binding free energy scanning of the metadynamics simulations on LY2456302–κ-OR and JDTic–κ-OR systems.