Systematic dissection of sequence features affecting binding specificity of a pioneer factor reveals binding synergy between FOXA1 and AP-1.

Despite the unique ability of pioneer factors (PFs) to target nucleosomal sites in closed chromatin, they only bind a small fraction of their genomic motifs. The underlying mechanism of this selectivity is not well understood. Here, we design a high-throughput assay called chromatin immunoprecipitation with integrated synthetic oligonucleotides (ChIP-ISO) to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1, in human A549 cells. Combining ChIP-ISO with in vitro and neural network analyses, we find that (1) FOXA1 binding is strongly affected by co-binding transcription factors (TFs) AP-1 and CEBPB; (2) FOXA1 and AP-1 show binding cooperativity in vitro; (3) FOXA1's binding is determined more by local sequences than chromatin context, including eu-/heterochromatin; and (4) AP-1 is partially responsible for differential binding of FOXA1 in different cell types. Our study presents a framework for elucidating genetic rules underlying PF binding specificity and reveals a mechanism for context-specific regulation of its binding.

The Transcription Factor ThPOK Orchestrates Stochastic Interchromosomal Interactions Required for IFNB1 Virus-Inducible Gene Expression.

Virus infection induces stochastic activation of the interferon-ß gene. Three previously identified Alu-like DNA elements called NRCs (NF-κB reception centers) function by capturing and delivering NF-κB to the IFNB1 enhancer via stochastic interchromosomal interactions. We show that the transcription factor ThPOK binds cooperatively with NF-κB to NRCs and mediates their physical proximity with the IFNB1 gene via its ability to oligomerize when bound to DNA. ThPOK knockdown significantly decreased the frequency of interchromosomal interactions, NF-κB DNA binding to the IFNB1 enhancer, and virus-induced IFNB1 gene activation. We also demonstrate that cooperative DNA binding between ThPOK and NF-κB on the same face of the double DNA helix is required for interchromosomal interactions and distinguishes NRCs from various other Alu elements bearing κB sites. These studies show how DNA binding cooperativity of stereospecifically aligned transcription factors provides the necessary ultrasensitivity for the all-or-none stochastic cell responses to virus infection.

Molecular mechanism of glycolytic flux control intrinsic to human phosphoglycerate kinase.

Glycolysis plays a fundamental role in energy production and metabolic homeostasis. The intracellular [adenosine triphosphate]/[adenosine diphosphate] ([ATP]/[ADP]) ratio controls glycolytic flux; however, the regulatory mechanism underlying reactions catalyzed by individual glycolytic enzymes enabling flux adaptation remains incompletely understood. Phosphoglycerate kinase (PGK) catalyzes the reversible phosphotransfer reaction, which directly produces ATP in a near-equilibrium step of glycolysis. Despite extensive studies on the transcriptional regulation of PGK expression, the mechanism in response to changes in the [ATP]/[ADP] ratio remains obscure. Here, we report a protein-level regulation of human PGK (hPGK) by utilizing the switching ligand-binding cooperativities between adenine nucleotides and 3-phosphoglycerate (3PG). This was revealed by nuclear magnetic resonance (NMR) spectroscopy at physiological salt concentrations. MgADP and 3PG bind to hPGK with negative cooperativity, whereas MgAMPPNP (a nonhydrolyzable ATP analog) and 3PG bind to hPGK with positive cooperativity. These opposite cooperativities enable a shift between different ligand-bound states depending on the intracellular [ATP]/[ADP] ratio. Based on these findings, we present an atomic-scale description of the reaction scheme for hPGK under physiological conditions. Our results indicate that hPGK intrinsically modulates its function via ligand-binding cooperativities that are finely tuned to respond to changes in the [ATP]/[ADP] ratio. The alteration of ligand-binding cooperativities could be one of the self-regulatory mechanisms for enzymes in bidirectional pathways, which enables rapid adaptation to changes in the intracellular environment.

Structure and dynamics of UBA5-UFM1 complex formation showing new insights in the UBA5 activation mechanism.

Ubiquitin fold modifier 1 (UFM1) is an ubiquitin-like protein (Ubl) involved especially in endoplasmic stress response. Activation occurs via a three-step mechanism like other Ubls. Data obtained reveal that UFM1 regulates the oligomeric state of ubiquitin activating enzyme 5 (UBA5) to initiate the activation step. Mixtures of homodimers and heterotrimers are observed in solution at the equilibrium state, demonstrating that the UBA5-UFM1 complex undergoes several concentration dependent oligomeric translational states to form a final functional complex. The oligomerization state of unbound UBA5 is also concentration dependent and shifts from the monomeric to the dimeric state. Data describing different oligomeric states are complemented with binding studies that reveal a negative cooperativity for the complex formation and thereby provide more detailed insights into the complex formation mechanism.

Measuring the affinity of protein-protein interactions on a single-molecule level by mass photometry.

Measurements of biomolecular interactions are crucial to understand the mechanisms of the biological processes they facilitate. Bulk-based methods such as ITC and SPR provide important information on binding affinities, stoichiometry, and kinetics of interactions. However, the ensemble averaging approaches are not able to probe the intrinsic heterogeneity often displayed by biological systems. Interactions that involve cooperativity or result in the formation of multicomponent complexes pose additional experimental challenges. Single-molecule techniques have previously been applied to solve these problems. However, single-molecule experiments are often technically demanding and require labeling or immobilization of the molecules under study. A recently developed single-molecule method, mass photometry (MP), overcomes these limitations. Here we applied MP to measure the affinities of biomolecular interactions. We have demonstrated how MP allows the user to study multivalent complexes and quantify the affinities of different binding sites in a single measurement. Results obtained from this single-molecule technique have been validated by ITC and BLI. The quality and information content of the MP data, combined with simple and fast measurements and low sample consumption makes MP a new preferred method for measuring strong protein-protein interactions.

Cooperative Binding and Stepwise Encapsulation of Drug Molecules by Sulfonylcalixarene-Based Metal-Organic Supercontainers.

The cooperative binding behavior of a face-directed octahedral metal-organic supercontainer featuring one endo cavity and six exo cavities was thoroughly examined in chloroform solution through ultraviolet-visible (UV-Vis) titration technique using two representative drug molecules as the guests. The titration curves and their nonlinear fit to Hill equation strongly suggest the efficient encapsulation of the guest molecules by the synthetic host, which exhibit interesting cooperative and stepwise binding behavior. Based on the control experiments using tetranuclear complex as a reference, it is clear that two equivalents of the guest molecules are initially encapsulated inside the endo cavity, followed by the trapping of six additional equivalents of the drug molecules through six exo cavities (1 eq. per exo cavity), and the remaining guests are entrapped by the external pockets. The results provide an in-depth understanding of the cooperative binding behavior of metal-organic supercontainers, which opens up new opportunities for designing synthetic receptors for truly biomimetic functional applications.

SEDPHAT--a platform for global ITC analysis and global multi-method analysis of molecular interactions.

Isothermal titration calorimetry experiments can provide significantly more detailed information about molecular interactions when combined in global analysis. For example, global analysis can improve the precision of binding affinity and enthalpy, and of possible linkage parameters, even for simple bimolecular interactions, and greatly facilitate the study of multi-site and multi-component systems with competition or cooperativity. A pre-requisite for global analysis is the departure from the traditional binding model, including an 'n'-value describing unphysical, non-integral numbers of sites. Instead, concentration correction factors can be introduced to account for either errors in the concentration determination or for the presence of inactive fractions of material. SEDPHAT is a computer program that embeds these ideas and provides a graphical user interface for the seamless combination of biophysical experiments to be globally modeled with a large number of different binding models. It offers statistical tools for the rigorous determination of parameter errors, correlations, as well as advanced statistical functions for global ITC (gITC) and global multi-method analysis (GMMA). SEDPHAT will also take full advantage of error bars of individual titration data points determined with the unbiased integration software NITPIC. The present communication reviews principles and strategies of global analysis for ITC and its extension to GMMA in SEDPHAT. We will also introduce a new graphical tool for aiding experimental design by surveying the concentration space and generating simulated data sets, which can be subsequently statistically examined for their information content. This procedure can replace the 'c'-value as an experimental design parameter, which ceases to be helpful for multi-site systems and in the context of gITC.

Resolving the contributions of two cooperative mechanisms to the DNA binding of AGT.

The O(6)-alkylguanine DNA alkyltransferase (AGT) is a DNA repair enzyme that binds DNA with moderate cooperativity. This cooperativity is important for its search for alkylated bases. A structural model of the cooperative complex of AGT with DNA predicts short-range interactions between nearest protein neighbors and long-range interactions between proteins separated in the array. DNA substrates ranging from 11bp to 30bp allowed us to use differences in binding stoichiometry to resolve short- and long-range protein contributions to the stability of AGT complexes. We found that the short-range component of ΔG°(coop) was nearly independent of DNA length and protein packing density. In contrast the long-range component oscillated with DNA length, with a period equal to the occluded binding site size (4bp). The amplitude of the long-range component decayed from ∼-4 kcal/mole of interaction to ∼-1.2 kcal/mol of interaction as the size of cooperative unit increased from 4 to 7 proteins, suggesting a mechanism to limit the size of cooperative clusters. These features allow us to make testable predictions about AGT distributions and interactions with chromatin structures in vivo.

The αC-ß4 loop controls the allosteric cooperativity between nucleotide and substrate in the catalytic subunit of protein kinase A.

Allosteric cooperativity between ATP and substrates is a prominent characteristic of the cAMP-dependent catalytic subunit of protein kinase A (PKA-C). This long-range synergistic action is involved in substrate recognition and fidelity, and it may also regulate PKA's association with regulatory subunits and other binding partners. To date, a complete understanding of this intramolecular mechanism is still lacking. Here, we integrated NMR(Nuclear Magnetic Resonance)-restrained molecular dynamics simulations and a Markov State Model to characterize the free energy landscape and conformational transitions of PKA-C. We found that the apoenzyme populates a broad free energy basin featuring a conformational ensemble of the active state of PKA-C (ground state) and other basins with lower populations (excited states). The first excited state corresponds to a previously characterized inactive state of PKA-C with the αC helix swinging outward. The second excited state displays a disrupted hydrophobic packing around the regulatory (R) spine, with a flipped configuration of the F100 and F102 residues at the αC-ß4 loop. We validated the second excited state by analyzing the F100A mutant of PKA-C, assessing its structural response to ATP and substrate binding. While PKA-CF100A preserves its catalytic efficiency with Kemptide, this mutation rearranges the αC-ß4 loop conformation, interrupting the coupling of the two lobes and abolishing the allosteric binding cooperativity. The highly conserved αC-ß4 loop emerges as a pivotal element to control the synergistic binding of nucleotide and substrate, explaining how mutations or insertions near or within this motif affect the function and drug sensitivity in homologous kinases.

Aberrant homeodomain-DNA cooperative dimerization underlies distinct developmental defects in two dominant CRX retinopathy models.

Paired-class homeodomain transcription factors (HD TFs) play essential roles in vertebrate development, and their mutations are linked to human diseases. One unique feature of paired-class HD is cooperative dimerization on specific palindrome DNA sequences. Yet, the functional significance of HD cooperative dimerization in animal development and its dysregulation in diseases remain elusive. Using the retinal TF Cone-rod Homeobox (CRX) as a model, we have studied how blindness-causing mutations in the paired HD, p.E80A and p.K88N, alter CRX's cooperative dimerization, lead to gene misexpression and photoreceptor developmental deficits in dominant manners. CRXE80A maintains binding at monomeric WT CRX motifs but is deficient in cooperative binding at dimeric motifs. CRXE80A's cooperativity defect impacts the exponential increase of photoreceptor gene expression in terminal differentiation and produces immature, non-functional photoreceptors in the CrxE80A retinas. CRXK88N is highly cooperative and localizes to ectopic genomic sites with strong enrichment of dimeric HD motifs. CRXK88N's altered biochemical properties disrupt CRX's ability to direct dynamic chromatin remodeling during development to activate photoreceptor differentiation programs and silence progenitor programs. Our study here provides in vitro and in vivo molecular evidence that paired-class HD cooperative dimerization regulates neuronal development and dysregulation of cooperative binding contributes to severe dominant blinding retinopathies.

The αC-ß4 loop controls the allosteric cooperativity between nucleotide and substrate in the catalytic subunit of protein kinase A.

Allosteric cooperativity between ATP and substrates is a prominent characteristic of the cAMP-dependent catalytic (C) subunit of protein kinase A (PKA). Not only this long-range synergistic action is involved in substrate recognition and fidelity, but it is likely to regulate PKA association with regulatory subunits and other binding partners. To date, a complete understanding of the molecular determinants for this intramolecular mechanism is still lacking. Here, we used an integrated NMR-restrained molecular dynamics simulations and a Markov Model to characterize the free energy landscape and conformational transitions of the catalytic subunit of protein kinase A (PKA-C). We found that the apo-enzyme populates a broad free energy basin featuring a conformational ensemble of the active state of PKA-C (ground state) and other basins with lower populations (excited states). The first excited state corresponds to a previously characterized inactive state of PKA-C with the αC helix swinging outward. The second excited state displays a disrupted hydrophobic packing around the regulatory (R) spine, with a flipped configuration of the F100 and F102 residues at the tip of the αC-ß4 loop. To experimentally validate the second excited state, we mutated F100 into alanine and used NMR spectroscopy to characterize the binding thermodynamics and structural response of ATP and a prototypical peptide substrate. While the activity of PKA-CF100A toward a prototypical peptide substrate is unaltered and the enzyme retains its affinity for ATP and substrate, this mutation rearranges the αC-ß4 loop conformation interrupting the allosteric coupling between nucleotide and substrate. The highly conserved αC-ß4 loop emerges as a pivotal element able to modulate the synergistic binding between nucleotide and substrate and may affect PKA signalosome. These results may explain how insertion mutations within this motif affect drug sensitivity in other homologous kinases.

Allosteric binding cooperativity in a kinetic context.

Allosteric modulators are of prime interest in drug discovery. These drugs regulate the binding and function of endogenous ligands, with some advantages over orthosteric ligands. A typical pharmacological parameter in allosteric modulation is binding cooperativity. This property can yield unexpected but illuminating results when decomposed into its kinetic parameters. Using two reference models (the allosteric ternary complex receptor model and a heterodimer receptor model), a relationship has been derived for the cooperativity rate constant parameters. This relationship allows many combinations of the cooperativity kinetic parameters for a single binding cooperativity value obtained under equilibrium conditions. This assessment may help understand striking experimental results involving allosteric modulation and suggest further investigations in the field.

Crystallographic and biophysical analyses of Pseudomonas aeruginosa ketopantoate reductase: Implications of ligand induced conformational changes in cofactor recognition.

Ketopantoate reductases (KPRs) catalyse NADPH-dependent reduction of ketopantoate to pantoate, the rate-limiting step of pantothenate biosynthetic pathway. In our recent study, we showed KPRs are under dynamic evolutionary selection and highlighted the possible role of ordered substrate binding kinetics for cofactor selection. To further delineate this at molecular level, here, we perform X-ray crystallographic and biophysical analyses of KPR in presence of non-canonical cofactor NAD+. In our structure, NAD+ was found to be highly dynamic in catalytic pocket of KPR, which could attain stable conformation only in presence of ketopantoate. Further, isothermal calorimetric (ITC) titrations showed that affinity of KPR for ketopantoate is higher in presence of NADP+ than in presence of NAD+ and lowest in absence of redox cofactors. In sum, our results clearly depict two modes of redox cofactor selections in KPRs, firstly by specific salt bridge interactions with unique phosphate moiety of NADP+ and secondly via ordered sequential heterotrophic cooperative binding of substrate ketopantoate.

Interplay of receptor-ligand binding and lipid domain formation during cell adhesion.

Cell adhesion involved in biological processes such as cell migration, immune responses, and cancer metastasis, is mediated by the specific binding of receptor and ligand proteins. Some of these proteins exhibit affinity for nanoscale lipid clusters in cell membranes. A key question is how these nanoscale lipid clusters influence and react to the receptor-ligand binding during cell adhesion. In this article, we review recent computational studies that shed new light on the interplay of the receptor-ligand binding and the formation of lipid domains in adhering membranes. These studies indicate that the receptor-ligand binding promotes coalescence of lipid clusters into mesoscale domains, which, in turn, enhances both the affinity and cooperativity of the receptor-ligand binding in cell-cell adhesion with mobile ligands. In contrast, in the case of cell-extracellular matrix adhesion with immobile ligands, the receptor-ligand binding and the lipid cluster coalescence can be correlated or anti-correlated, depending strongly on the ligand distribution. These findings deepen our understanding of correlations between cell adhesion and membrane heterogeneities.

Rely on Each Other: DNA Binding Cooperativity Shapes p53 Functions in Tumor Suppression and Cancer Therapy.

p53 is a tumor suppressor that is mutated in half of all cancers. The high clinical relevance has made p53 a model transcription factor for delineating general mechanisms of transcriptional regulation. p53 forms tetramers that bind DNA in a highly cooperative manner. The DNA binding cooperativity of p53 has been studied by structural and molecular biologists as well as clinical oncologists. These experiments have revealed the structural basis for cooperative DNA binding and its impact on sequence specificity and target gene spectrum. Cooperativity was found to be critical for the control of p53-mediated cell fate decisions and tumor suppression. Importantly, an estimated number of 34,000 cancer patients per year world-wide have mutations of the amino acids mediating cooperativity, and knock-in mouse models have confirmed such mutations to be tumorigenic. While p53 cancer mutations are classically subdivided into "contact" and "structural" mutations, "cooperativity" mutations form a mechanistically distinct third class that affect the quaternary structure but leave DNA contacting residues and the three-dimensional folding of the DNA-binding domain intact. In this review we discuss the concept of DNA binding cooperativity and highlight the unique nature of cooperativity mutations and their clinical implications for cancer therapy.

Conformational flexibility of the conserved hydrophobic pocket of HIV-1 gp41. Implications for the discovery of small-molecule fusion inhibitors.

During HIV-1 infection, the envelope glycoprotein subunit gp41 folds into a six-helix bundle structure (6HB) formed by the interaction between its N-terminal (NHR) and C-terminal (CHR) heptad-repeats, promoting viral and cell membranes fusion. A highly preserved, hydrophobic pocket (HP) on the NHR surface is crucial in 6HB formation and, therefore, HP-binding compounds constitute promising therapeutics against HIV-1. Here, we investigated the conformational and dynamic properties of the HP using a rationally designed single-chain protein (named covNHR) that mimics the gp41 NHR structure. We found that the fluorescent dye 8-anilino-naphtalene-1-sulfonic acid (ANS) binds specifically to the HP, suggesting that ANS derivatives may constitute lead compounds to inhibit 6HB formation. ANS shows different binding modes to the HP, depending on the occupancy of other NHR pockets. Moreover, in presence of a CHR peptide bound to the N-terminal pockets in gp41, two ANS molecules can occupy the HP showing cooperative behavior. This binding mode was assessed using molecular docking and molecular dynamics simulations. The results show that the HP is conformationally flexible and connected allosterically to other NHR regions, which strongly influence the binding of potential ligands. These findings could guide the development of small-molecule HIV-1 inhibitors targeting the HP.

Is Disrupted Nucleotide-Substrate Cooperativity a Common Trait for Cushing's Syndrome Driving Mutations of Protein Kinase A?

Somatic mutations in the PRKACA gene encoding the catalytic α subunit of protein kinase A (PKA-C) are responsible for cortisol-producing adrenocortical adenomas. These benign neoplasms contribute to the development of Cushing's syndrome. The majority of these mutations occur at the interface between the two lobes of PKA-C and interfere with the enzyme's ability to recognize substrates and regulatory (R) subunits, leading to aberrant phosphorylation patterns and activation. Rarely, patients with similar phenotypes carry an allosteric mutation, E31V, located at the C-terminal end of the αA-helix and adjacent to the αC-helix, but structurally distinct from the PKA-C/R subunit interface mutations. Using a combination of solution NMR, thermodynamics, kinetic assays, and molecular dynamics simulations, we show that the E31V allosteric mutation disrupts central communication nodes between the N- and C- lobes of the enzyme as well as nucleotide-substrate binding cooperativity, a hallmark for kinases' substrate fidelity and regulation. For both orthosteric (L205R and W196R) and allosteric (E31V) Cushing's syndrome mutants, the loss of binding cooperativity is proportional to the density of the intramolecular allosteric network. This structure-activity relationship suggests a possible common mechanism for Cushing's syndrome driving mutations in which decreased nucleotide/substrate binding cooperativity is linked to loss in substrate fidelity and dysfunctional regulation.

Insights into the molecular mechanism of positive cooperativity between partial agonist MK-8666 and full allosteric agonist AP8 of hGPR40 by Gaussian accelerated molecular dynamics (GaMD) simulations.

Activation of human free fatty acid receptor 1 (FFAR1, also called hGPR40) enhances insulin secretion in a glucose-dependent manner. Hence, the development of selective agonist targeting hGPR40 has been proposed as a therapeutic strategy of type 2 diabetes mellitus. Some agonists targeting hGPR40 were reported. The radioligand-binding studies and the crystal structures reveal that there are multiple sites on GPR40, and there exists positive binding cooperativity between the partial agonist MK-8666 and full allosteric agonist (AgoPAM) AP8. In this work, we carried out long-time Gaussian accelerated molecular dynamics (GaMD) simulations on hGPR40 to shed light on the mechanism of the cooperativity between the two agonists at different sites. Our results reveal that the induced-fit conformational coupling is bidirectional between the two sites. The movements and rotations of TM3, TM4, TM5 and TM6 due to their inherent flexibility are crucial in coupling the conformational changes of the two agonists binding sites. These helices adopt similar conformational states upon alternative ligand or both ligands binding. The Leu1384.57, Leu1865.42 and Leu1905.46 play roles in coordinating the rearrangements of residues in the two pockets, which makes the movements of residues in the two sites like gear movements. These results provide detailed information at the atomic level about the conformational coupling between different sites of GPR40, and also provide the structural information for further design of new agonists of GPR40. In addition, these results suggest that it is necessary by considering the effect of other site bound in structure-based ligands discovery.

Evaluation of structurally different brominated flame retardants interacting with the transthyretin and their toxicity on HepG2 cells.

Brominated flame retardants (BFRs) are found at quantifiable levels in both humans and wildlife and may potentially cause a health risk. For BFRs and their derivatives, limited information regarding the relationship among the structure, binding affinity to the target protein and toxicity is currently available. In the present work, representative BFRs with different hydroxyl- or bromo-substituents, namely 2, 2', 4, 4'-tetrabromodiphenyl ether (BDE-47), 3-hydroxy-2, 2', 4, 4'-tetrabromodiphenyl ether (3-OH-BDE-47) and tetrabromobisphenol A (TBBPA), were selected to investigate the interactions with transthyretin (TTR) by electrospray ionization mass spectrometry (ESI-MS) and cytotoxicity on HepG2 cells. It was noted that BDE-47 had a weak binding affinity to TTR, while 3-OH-BDE-47 and TBBPA had a stronger binding affinity than BDE-47 and thyroxine (T4). Hence, 3-OH-BDE-47 and TBBPA could affect the binding of TTR with its native ligand T4 by competitive binding to TTR, even at equal concentrations, which might be associated with BFR toxicity of endocrine disruption. Negative cooperativity was found for 3-OH-BDE-47 and TBBPA binding to TTR, similar to T4 with a well-established negatively cooperative binding mechanism. The tendency of toxic effects on HepG2 cells for these three BFRs was, 3-OH-BDE-47 > TBBPA > BDE-47, and this order was in good agreement with the binding ability explored by ESI-MS experiments and molecular docking simulation. The observations obtained by this study demonstrate that the binding properties of these BFRs to TTR and their cytotoxicity are correlated with structure differentials and functional substituents.

Revealing the Positive Binding Cooperativity Mechanism between the Orthosteric and the Allosteric Antagonists of CCR2 by Metadynamics and Gaussian Accelerated Molecular Dynamics Simulations.

CC chemokine receptor 2 (CCR2) and its endogenous CC chemokine ligands are associated with numerous inflammatory, neurodegenerative diseases, and cancer. CCR2 is becoming an attractive target in the treatment of autoimmune disease and neurodegenerative diseases. The orthosteric antagonist BMS-681 and allosteric antagonist CCR2-RA-[R] of CCR2 show positive binding cooperativity. We performed well-tempered metadynamics simulations and Gaussian accelerated MD simulations to reveal the influence of the orthosteric antagonist on the unbinding of allosteric antagonist of CCR2. We revealed different unbinding pathways of CCR2-RA-[R] in binary complex CCR2-VT5 and ternary complex CCR2-73R-VT5. The different unbinding pathways of CCR2-RA-[R] are due to the conformational dynamics of TM6. We obtained the significant conformational differences of the intracellular side of TM6 upon CCR2 binding to different ligands by GaMD simulation. The conformational dynamics of TM6 are consistent with the unbinding pathway analysis. GaMD simulations indicate that BMS-681 binding restricts the bend of intracellular side of TM6 by stabilizing the extracellular sides of TM6 and TM7. The charged residues Arg2065.43 of TM5 and Glu2917.39 of TM7 play key roles in stabling TM7 and TM6. TM6 and TM7 are crucial components in the orthosteric and allosteric binding sites. Our results illustrate the conformational details about the effect of the orthosteric antagonist on the allosteric antagonist of CCR2. The conformational dynamics of CCR2 upon binding to different ligands can provide a rational basis for development of allosteric ligands of CCR2.