A multi-iron enzyme installs copper-binding oxazolone/thioamide pairs on a nontypeable Haemophilus influenzae virulence factor.

The multinuclear nonheme iron-dependent oxidases (MNIOs) are a rapidly growing family of enzymes involved in the biosynthesis of ribosomally synthesized, posttranslationally modified peptide natural products (RiPPs). Recently, a secreted virulence factor from nontypeable Haemophilus influenzae (NTHi) was found to be expressed from an operon, which we designate the hvf operon, that also encodes an MNIO. Here, we show by Mössbauer spectroscopy that the MNIO HvfB contains a triiron cofactor. We demonstrate that HvfB works together with HvfC [a RiPP recognition element (RRE)-containing partner protein] to perform six posttranslational modifications of cysteine residues on the virulence factor precursor peptide HvfA. Structural characterization by tandem mass spectrometry and NMR shows that these six cysteine residues are converted to oxazolone and thioamide pairs, similar to those found in the RiPP methanobactin. Like methanobactin, the mature virulence factor, which we name oxazolin, uses these modified residues to coordinate Cu(I) ions. Considering the necessity of oxazolin for host cell invasion by NTHi, these findings point to a key role for copper during NTHi infection. Furthermore, oxazolin and its biosynthetic pathway represent a potential therapeutic target for NTHi.

Insights on the G protein-coupled receptor helix 8 solution structure and orientation using a neurotensin receptor 1 peptide.

G-protein coupled receptors (GPCRs) are the largest class of membrane proteins encoded in the human genome with high pharmaceutical relevance and implications to human health. These receptors share a prevalent architecture of seven transmembrane helices followed by an intracellular, amphipathic helix 8 (H8) and a disordered C-terminal tail (Ctail). Technological advancements have led to over 1000 receptor structures in the last two decades, yet frequently H8 and the Ctail are conformationally heterogeneous or altogether absent. Here we synthesize a peptide comprising the neurotensin receptor 1 (NTS1) H8 and Ctail (H8-Ctail) to investigate its structural stability, conformational dynamics, and orientation in the presence of detergent and phospholipid micelles, which mimic the membrane. Circular dichroism (CD) and nuclear magnetic resonance (NMR) measurements confirm that zwitterionic 1,2-diheptanoyl-sn-glycero-3-phosphocholine is a potent stabilizer of H8 structure, whereas the commonly-used branched detergent lauryl maltose neopentyl glycol (LMNG) is unable to completely stabilize the helix - even at amounts four orders of magnitude greater than its critical micellar concentration. We then used NMR spectroscopy to assign the backbone chemical shifts. A series of temperature and lipid titrations were used to define the H8 boundaries as F376-R392 from chemical shift perturbations, changes in resonance intensity, and chemical-shift-derived phi/psi angles. Finally, the H8 azimuthal and tilt angles, defining the helix orientation relative of the membrane normal were measured using paramagnetic relaxation enhancement NMR. Taken together, our studies reveal the H8-Ctail region is sensitive to membrane physicochemical properties and is capable of more adaptive behavior than previously suggested by static structural techniques.

Characterizing the monomer-dimer equilibrium of UbcH8/Ube2L6: A combined SAXS and NMR study.

Interferon-stimulated gene-15 (ISG15) is an interferon-induced protein with two ubiquitin-like (Ubl) domains linked by a short peptide chain, and the conjugated protein of the ISGylation system. Similar to ubiquitin and other Ubls, ISG15 is ligated to its target proteins with a series of E1, E2, and E3 enzymes known as Uba7, Ube2L6/UbcH8, and HERC5, respectively. Ube2L6/UbcH8 plays a literal central role in ISGylation, underscoring it as an important drug target for boosting innate antiviral immunity. Depending on the type of conjugated protein and the ultimate target protein, E2 enzymes have been shown to function as monomers, dimers, or both. UbcH8 has been crystalized in both monomeric and dimeric forms, but the functional state is unclear. Here, we used a combined approach of small-angle X-ray scattering (SAXS) and nuclear magnetic resonance (NMR) spectroscopy to characterize UbcH8's oligomeric state in solution. SAXS revealed a dimeric UbcH8 structure that could be dissociated when fused with an N-terminal glutathione S-transferase molecule. NMR spectroscopy validated the presence of a concentration-dependent monomer-dimer equilibrium and suggested a backside dimerization interface. Chemical shift perturbation and peak intensity analysis further suggest dimer-induced conformational dynamics at ISG15 and E3 interfaces - providing hypotheses for the protein's functional mechanisms. Our study highlights the power of combining NMR and SAXS techniques in providing structural information about proteins in solution.

Chemically induced partial unfolding of the multifunctional Apurinic/apyrimidinic endonuclease 1.

Targeting of the multifunctional enzyme apurinic/apyrimidinic endonuclease I/redox factor 1 (APE1) has produced small molecule inhibitors of both its endonuclease and redox activities. While one of the small molecules, the redox inhibitor APX3330, completed a Phase I clinical trial for solid tumors and a Phase II clinical trial for Diabetic Retinopathy/Diabetic Macular Edema, the mechanism of action for this drug has yet to be fully understood. Here, we demonstrate through HSQC NMR studies that APX3330 induces chemical shift perturbations (CSPs) of both surface and internal residues in a concentration-dependent manner, with a cluster of surface residues defining a small pocket on the opposite face from the endonuclease active site of APE1. Furthermore, APX3330 induces partial unfolding of APE1 as evidenced by a time-dependent loss of chemical shifts for approximately 35% of the residues within APE1 in the HSQC NMR spectrum. Notably, regions that are partially unfolded include adjacent strands within one of two beta sheets that comprise the core of APE1. One of the strands comprises residues near the N-terminal region and a second strand is contributed by the C-terminal region of APE1, which serves as a mitochondrial targeting sequence. These terminal regions converge within the pocket defined by the CSPs. In the presence of a duplex DNA substrate mimic, removal of excess APX3330 resulted in refolding of APE1. Our results are consistent with a reversible mechanism of partial unfolding of APE1 induced by the small molecule inhibitor, APX3330, defining a novel mechanism of inhibition.

Stabilization of pre-existing neurotensin receptor conformational states by ß-arrestin-1 and the biased allosteric modulator ML314.

The neurotensin receptor 1 (NTS1) is a G protein-coupled receptor (GPCR) with promise as a drug target for the treatment of pain, schizophrenia, obesity, addiction, and various cancers. A detailed picture of the NTS1 structural landscape has been established by X-ray crystallography and cryo-EM and yet, the molecular determinants for why a receptor couples to G protein versus arrestin transducers remain poorly defined. We used 13CεH3-methionine NMR spectroscopy to show that binding of phosphatidylinositol-4,5-bisphosphate (PIP2) to the receptor's intracellular surface allosterically tunes the timescale of motions at the orthosteric pocket and conserved activation motifs - without dramatically altering the structural ensemble. ß-arrestin-1 further remodels the receptor ensemble by reducing conformational exchange kinetics for a subset of resonances, whereas G protein coupling has little to no effect on exchange rates. A ß-arrestin biased allosteric modulator transforms the NTS1:G protein complex into a concatenation of substates, without triggering transducer dissociation, suggesting that it may function by stabilizing signaling incompetent G protein conformations such as the non-canonical state. Together, our work demonstrates the importance of kinetic information to a complete picture of the GPCR activation landscape.

Ligands selectively tune the local and global motions of neurotensin receptor 1 (NTS1).

Nuclear magnetic resonance (NMR) studies have revealed that fast methyl sidechain dynamics can report on entropically-driven allostery. Yet, NMR applications have been largely limited to the super-microsecond motional regimes of G protein-coupled receptors (GPCRs). We use 13Cε-methionine chemical shift-based global order parameters to test if ligands affect the fast dynamics of a thermostabilized GPCR, neurotensin receptor 1 (NTS1). We establish that the NTS1 solution ensemble includes substates with lifetimes on several, discrete timescales. The longest-lived states reflect those captured in agonist- and inverse agonist-bound crystal structures, separated by large energy barriers. We observe that the rapid fluctuations of individual methionine residues, superimposed on these long-lived states, respond collectively with the degree of fast, global dynamics correlating with ligand pharmacology. This approach lends confidence to interpreting spectra in terms of local structure and methyl dihedral angle geometry. The results suggest a role for sub-microsecond dynamics and conformational entropy in GPCR ligand discrimination.

A method for selective 19 F-labeling absent of probe sequestration (SLAPS).

Fluorine (19 F) offers several distinct advantages for biomolecular nuclear magnetic resonance spectroscopy such as no background signal, 100% natural abundance, high sensitivity, and a large chemical shift range. Exogenous cysteine-reactive 19 F-probes have proven especially indispensable for characterizing large, challenging systems that are less amenable to other isotopic labeling strategies such as G protein-coupled receptors. As fluorine linewidths are inherently broad, limiting reactions with offsite cysteines is critical for spectral simplification and accurate deconvolution of component peaks-especially when analyzing systems with intermediate to slow timescale conformational exchange. Here, we uncovered noncovalent probe sequestration by detergent proteomicelles as a second source of offsite labeling when using the popular 19 F-probe BTFMA (2-bromo-N-(4-[trifluoromethyl]phenyl)acetamide). The chemical shift and relaxation rates of these unreacted 19 F-BTFMA molecules are insufficient to distinguish them from protein-conjugates, but they can be easily identified using mass spectrometry. We present a simple four-step protocol for Selective Labeling Absent of Probe Sequestration (SLAPS): physically disrupt cell membranes in the absence of detergent, incubate membranes with cysteine-reactive 19 F-BTFMA, remove excess unreacted 19 F-BTFMA molecules via ultracentrifugation, and finally solubilize in the detergent of choice. Our approach builds upon the in-membrane chemical modification method with the addition of one crucial step: removal of unreacted 19 F-probes by ultracentrifugation prior to detergent solubilization. SLAPS is broadly applicable to other lipophilic cysteine-reactive probes and membrane protein classes solubilized in detergent micelles or lipid mimetics.

Effect of Ligands and Transducers on the Neurotensin Receptor 1 Conformational Ensemble.

Using a discrete, intracellular 19F nuclear magnetic resonance (NMR) probe on transmembrane helix 6 of the neurotensin receptor 1 (NTS1), we aim to understand how ligands and transducers modulate the receptor's structural ensemble in a solution. For apo NTS1, 19F NMR spectra reveal an ensemble of at least three conformational substates (one inactive and two active-like) in equilibrium that exchange on the millisecond to second timescale. Dynamic NMR experiments reveal that these substates follow a linear three-site exchange process that is both thermodynamically and kinetically remodeled by orthosteric ligands. As previously observed in other G protein-coupled receptors (GPCRs), the full agonist is insufficient to completely stabilize the active-like state. The inactive substate is abolished upon coupling to ß-arrestin-1 (ßArr1) or the C-terminal helix of Gαq, which comprises ≳60% of the GPCR/G protein interface surface area. Whereas ßArr1 exclusively selects for pre-existing active-like substates, the Gαq peptide induces a new substate. Both transducer molecules promote substantial line broadening of active-like states, suggesting contributions from additional microsecond to millisecond exchange processes. Together, our study suggests that (i) the NTS1 allosteric activation mechanism may be alternatively dominated by induced fit or conformational selection depending on the coupled transducer, and (ii) the available static structures do not represent the entire conformational ensemble observed in a solution.

Expanding the genotype and phenotype spectrum of SYT1-associated neurodevelopmental disorder.

PURPOSE: Synaptotagmin-1 (SYT1) is a critical mediator of neurotransmitter release in the central nervous system. Previously reported missense SYT1 variants in the C2B domain are associated with severe intellectual disability, movement disorders, behavioral disturbances, and electroencephalogram abnormalities. In this study, we expand the genotypes and phenotypes and identify discriminating features of this disorder. METHODS: We describe 22 individuals with 15 de novo missense SYT1 variants. The evidence for pathogenicity is discussed, including the American College of Medical Genetics and Genomics/Association for Molecular Pathology criteria, known structure-function relationships, and molecular dynamics simulations. Quantitative behavioral data for 14 cases were compared with other monogenic neurodevelopmental disorders. RESULTS: Four variants were located in the C2A domain with the remainder in the C2B domain. We classified 6 variants as pathogenic, 4 as likely pathogenic, and 5 as variants of uncertain significance. Prevalent clinical phenotypes included delayed developmental milestones, abnormal eye physiology, movement disorders, and sleep disturbances. Discriminating behavioral characteristics were severity of motor and communication impairment, presence of motor stereotypies, and mood instability. CONCLUSION: Neurodevelopmental disorder-associated SYT1 variants extend beyond previously reported regions, and the phenotypic spectrum encompasses a broader range of severities than initially reported. This study guides the diagnosis and molecular understanding of this rare neurodevelopmental disorder and highlights a key role for SYT1 function in emotional regulation, motor control, and emergent cognitive function.

Molecular basis for the evolved instability of a human G-protein coupled receptor.

Membrane proteins are prone to misfolding and degradation. This is particularly true for mammalian forms of the gonadotropin-releasing hormone receptor (GnRHR). Although they function at the plasma membrane, mammalian GnRHRs accumulate within the secretory pathway. Their apparent instability is believed to have evolved through selection for attenuated GnRHR activity. Nevertheless, the molecular basis of this adaptation remains unclear. We show that adaptation coincides with a C-terminal truncation that compromises the translocon-mediated membrane integration of its seventh transmembrane domain (TM7). We also identify a series of polar residues in mammalian GnRHRs that compromise the membrane integration of TM2 and TM6. Reverting a lipid-exposed polar residue in TM6 to an ancestral hydrophobic residue restores expression with no impact on function. Evolutionary trends suggest variations in the polarity of this residue track with reproductive phenotypes. Our findings suggest that the marginal energetics of cotranslational folding can be exploited to tune membrane protein fitness.

TRACT revisited: an algebraic solution for determining overall rotational correlation times from cross-correlated relaxation rates.

Accurate rotational correlation times ([Formula: see text]) are critical for quantitative analysis of fast timescale NMR dynamics. As molecular weights increase, the classic derivation of [Formula: see text] using transverse and longitudinal relaxation rates becomes increasingly unsuitable due to the non-trivial contribution of remote dipole-dipole interactions to longitudinal relaxation. Derivations using cross-correlated relaxation experiments, such as TRACT, overcome these limitations but are erroneously calculated in 65% of the citing literature. Herein, we developed an algebraic solutions to the Goldman relationship that facilitate rapid, point-by-point calculations for straightforward identification of appropriate spectral regions where global tumbling is likely to be dominant. The rigid-body approximation of the Goldman relationship has been previously shown to underestimate TRACT-based rotational correlation time estimates. This motivated us to develop a second algebraic solution that employs a simplified model-free spectral density function including an order parameter term that could, in principle, be set to an average backbone S2 ≈ 0.9 to further improve the accuracy of [Formula: see text] estimation. These solutions enabled us to explore the boundaries of the Goldman relationship as a function of the H-N internuclear distance ([Formula: see text]), difference of the two principal components of the axially-symmetric 15N CSA tensor ([Formula: see text]), and angle of the CSA tensor relative to the N-H bond vector ([Formula: see text]). We hope our algebraic solutions and analytical strategies will increase the accuracy and application of the TRACT experiment.

Biased antagonism of CXCR4 avoids antagonist tolerance.

Repeated dosing of drugs targeting G protein-coupled receptors can stimulate antagonist tolerance, which reduces their efficacy; thus, strategies to avoid tolerance are needed. The efficacy of AMD3100, a competitive antagonist of the chemokine receptor CXCR4 that mobilizes leukemic blasts from the bone marrow into the blood to sensitize them to chemotherapy, is reduced after prolonged treatment. Tolerance to AMD3100 increases the abundance of CXCR4 on the surface of leukemic blasts, which promotes their rehoming to the bone marrow. AMD3100 inhibits both G protein signaling by CXCR4 and ß-arrestin1/2-dependent receptor endocytosis. We demonstrated that biased antagonists of G protein-dependent chemotaxis but not ß-arrestin1/2 recruitment and subsequent receptor endocytosis avoided tolerance. The peptide antagonist X4-2-6, which is derived from transmembrane helix 2 and extracellular loop 1 of CXCR4, limited chemotaxis and signaling but did not promote CXCR4 accumulation on the cell surface or cause tolerance. The activity of X4-2-6 was due to its distinct mechanism of inhibition of CXCR4. The peptide formed a ternary complex with the receptor and its ligand, the chemokine CXCL12. Within this complex, X4-2-6 released the portion of CXCL12 critical for receptor-mediated activation of G proteins but enabled the rest of the chemokine to recruit ß-arrestins to the receptor. In contrast, AMD3100 displaced all components of the chemokine responsible for CXCR4 activation. We further identified a small molecule with similar biased antagonist properties to those of X4-2-6, which may provide a viable alternative to patients when antagonist tolerance prevents drugs from reaching efficacy.

Bicelle size modulates the rate of bacteriorhodopsin folding.

The conformational equilibria of integral membrane proteins have proven extremely difficult to characterize within native lipid bilayers. To circumvent technical issues, investigations of the structure and stability of α-helical membrane proteins are often carried out in mixed micelle or bicelle solvents that mimic the membrane and facilitate measurements of reversible folding. Under these conditions, the energetics of membrane protein folding are typically proportional to the mole fraction of an anionic detergent in the micelle. However, investigations of the folding and unfolding of bacteriorhodopsin (bR) surprisingly revealed that the folding rate is also highly sensitive to the bulk molar concentration of lipids and detergents. We show here that this rate enhancement coincides with changes in bicelle size and suggest this effect arises through restriction of the conformational search space during folding. In conjunction with previous mutagenic studies, these results provide additional evidence that a topological search limits the rate of bR folding. Furthermore, this finding provides insights into the manner by which micellar and bicellar environments influence the conformational stability of polytopic membrane proteins.

15N detection harnesses the slow relaxation property of nitrogen: Delivering enhanced resolution for intrinsically disordered proteins.

Studies over the past decade have highlighted the functional significance of intrinsically disordered proteins (IDPs). Due to conformational heterogeneity and inherent dynamics, structural studies of IDPs have relied mostly on NMR spectroscopy, despite IDPs having characteristics that make them challenging to study using traditional 1H-detected biomolecular NMR techniques. Here, we develop a suite of 3D 15N-detected experiments that take advantage of the slower transverse relaxation property of 15N nuclei, the associated narrower linewidth, and the greater chemical shift dispersion compared with those of 1H and 13C resonances. The six 3D experiments described here start with aliphatic 1H magnetization to take advantage of its higher initial polarization, and are broadly applicable for backbone assignment of proteins that are disordered, dynamic, or have unfavorable amide proton exchange rates. Using these experiments, backbone resonance assignments were completed for the unstructured regulatory domain (residues 131-294) of the human transcription factor nuclear factor of activated T cells (NFATC2), which includes 28 proline residues located in functionally important serine-proline (SP) repeats. The complete assignment of the NFATC2 regulatory domain enabled us to study phosphorylation of NFAT by kinase PKA and phosphorylation-dependent binding of chaperone protein 14-3-3 to NFAT, providing mechanistic insight on how 14-3-3 regulates NFAT nuclear translocation.

Recent developments in solution nuclear magnetic resonance (NMR)-based molecular biology.

Visualizing post-translational modifications, conformations, and interaction surfaces of protein structures at atomic resolution underpins the development of novel therapeutics to combat disease. As computational resources expand, in silico calculations coupled with experimentally derived structures and functional assays have led to an explosion in structure-based drug design (SBDD) with several compounds in clinical trials. It is increasingly clear that "hidden" transition-state structures along activation trajectories can be harnessed to develop novel classes of allosteric inhibitors. The goal of this mini-review is to empower the clinical researcher with a general knowledge of the strengths and weaknesses of nuclear magnetic resonance (NMR) spectroscopy in molecular medicine. Although NMR can determine protein structures at atomic resolution, its unrivaled strength lies in sensing subtle changes in a nuclei's chemical environment as a result of intrinsic conformational dynamics, solution conditions, and binding interactions. These can be recorded at atomic resolution, without explicit structure determination, and then incorporated with static structures or molecular dynamics simulations to produce a complete biological picture.

Structural basis for chemokine recognition by a G protein-coupled receptor and implications for receptor activation.

Chemokines orchestrate cell migration for development, immune surveillance, and disease by binding to cell surface heterotrimeric guanine nucleotide-binding protein (G protein)-coupled receptors (GPCRs). The array of interactions between the nearly 50 chemokines and their 20 GPCR targets generates an extensive signaling network to which promiscuity and biased agonism add further complexity. The receptor CXCR4 recognizes both monomeric and dimeric forms of the chemokine CXCL12, which is a distinct example of ligand bias in the chemokine family. We demonstrated that a constitutively monomeric CXCL12 variant reproduced the G protein-dependent and ß-arrestin-dependent responses that are associated with normal CXCR4 signaling and lead to cell migration. In addition, monomeric CXCL12 made specific contacts with CXCR4 that are not present in the structure of the receptor in complex with a dimeric form of CXCL12, a biased agonist that stimulates only G protein-dependent signaling. We produced an experimentally validated model of an agonist-bound chemokine receptor that merged a nuclear magnetic resonance-based structure of monomeric CXCL12 bound to the amino terminus of CXCR4 with a crystal structure of the transmembrane domains of CXCR4. The large CXCL12:CXCR4 protein-protein interface revealed by this structure identified previously uncharacterized functional interactions that fall outside of the classical "two-site model" for chemokine-receptor recognition. Our model suggests a mechanistic hypothesis for how interactions on the extracellular face of the receptor may stimulate the conformational changes required for chemokine receptor-mediated signal transduction.

New paradigms in chemokine receptor signal transduction: Moving beyond the two-site model.

Chemokine receptor (CKR) signaling forms the basis of essential immune cellular functions, and dysregulated CKR signaling underpins numerous disease processes of the immune system and beyond. CKRs, which belong to the seven transmembrane domain receptor (7TMR) superfamily, initiate signaling upon binding of endogenous, secreted chemokine ligands. Chemokine-CKR interactions are traditionally described by a two-step/two-site mechanism, in which the CKR N-terminus recognizes the chemokine globular core (i.e. site 1 interaction), followed by activation when the unstructured chemokine N-terminus is inserted into the receptor TM bundle (i.e. site 2 interaction). Several recent studies challenge the structural independence of sites 1 and 2 by demonstrating physical and allosteric links between these supposedly separate sites. Others contest the functional independence of these sites, identifying nuanced roles for site 1 and other interactions in CKR activation. These developments emerge within a rapidly changing landscape in which CKR signaling is influenced by receptor PTMs, chemokine and CKR dimerization, and endogenous non-chemokine ligands. Simultaneous advances in the structural and functional characterization of 7TMR biased signaling have altered how we understand promiscuous chemokine-CKR interactions. In this review, we explore new paradigms in CKR signal transduction by considering studies that depict a more intricate architecture governing the consequences of chemokine-CKR interactions.

Structural analysis of a novel small molecule ligand bound to the CXCL12 chemokine.

CXCL12 binds to CXCR4, promoting both chemotaxis of lymphocytes and metastasis of cancer cells. We previously identified small molecule ligands that bind CXCL12 and block CXCR4-mediated chemotaxis. We now report a 1.9 Å resolution X-ray structure of CXCL12 bound by such a molecule at a site normally bound by sY21 of CXCR4. The complex structure reveals binding hot spots for future inhibitor design and suggests a new approach to targeting CXCL12-CXCR4 signaling in drug discovery.

Discovery and characterization of a disulfide-locked C(2)-symmetric defensin peptide.

We report the discovery of HD5-CD, an unprecedented C2-symmetric ß-barrel-like covalent dimer of the cysteine-rich host-defense peptide human defensin 5 (HD5). Dimerization results from intermonomer disulfide exchange between the canonical α-defensin Cys(II)-Cys(IV) (Cys(5)-Cys(20)) bonds located at the hydrophobic interface. This disulfide-locked dimeric assembly provides a new element of structural diversity for cysteine-rich peptides as well as increased protease resistance, broad-spectrum antimicrobial activity, and enhanced potency against the opportunistic human pathogen Acinetobacter baumannii.

Chemokine cooperativity is caused by competitive glycosaminoglycan binding.

Chemokines comprise a family of secreted proteins that activate G protein-coupled chemokine receptors and thereby control the migration of leukocytes during inflammation or immune surveillance. The positional information required for such migratory behavior is governed by the binding of chemokines to membrane-tethered glycosaminoglycans (GAGs), which establishes a chemokine concentration gradient. An often observed but incompletely understood behavior of chemokines is the ability of unrelated chemokines to enhance the potency with which another chemokine subtype can activate its cognate receptor. This phenomenon has been demonstrated to occur between many chemokine combinations and across several model systems and has been dubbed chemokine cooperativity. In this study, we have used GAG binding-deficient chemokine mutants and cell-based functional (migration) assays to demonstrate that chemokine cooperativity is caused by competitive binding of chemokines to GAGs. This mechanistic explanation of chemokine cooperativity provides insight into chemokine gradient formation in the context of inflammation, in which multiple chemokines are secreted simultaneously.