Showing papers on "Cooperative binding published in 2018"

Small-Molecule Positive Allosteric Modulators of the β2-Adrenoceptor Isolated from DNA-Encoded Libraries.

Abstract: Conventional drug discovery efforts at the β2-adrenoceptor (β2AR) have led to the development of ligands that bind almost exclusively to the receptor’s hormone-binding orthosteric site. However, targeting the largely unexplored and evolutionarily unique allosteric sites has potential for developing more specific drugs with fewer side effects than orthosteric ligands. Using our recently developed approach for screening G protein–coupled receptors (GPCRs) with DNA-encoded small-molecule libraries, we have discovered and characterized the first β2AR small-molecule positive allosteric modulators (PAMs)—compound (Cmpd)-6 [(R)-N-(4-amino-1-(4-(tert-butyl)phenyl)-4-oxobutan-2-yl)-5-(N-isopropyl-N-methylsulfamoyl)-2-((4-methoxyphenyl)thio)benzamide] and its analogs. We used purified human β2ARs, occupied by a high-affinity agonist, for the affinity-based screening of over 500 million distinct library compounds, which yielded Cmpd-6. It exhibits a low micro-molar affinity for the agonist-occupied β2AR and displays positive cooperativity with orthosteric agonists, thereby enhancing their binding to the receptor and ability to stabilize its active state. Cmpd-6 is cooperative with G protein and β-arrestin1 (a.k.a. arrestin2) to stabilize high-affinity, agonist-bound active states of the β2AR and potentiates downstream cAMP production and receptor recruitment of β-arrestin2 (a.k.a. arrestin3). Cmpd-6 is specific for the β2AR compared with the closely related β1AR. Structure–activity studies of select Cmpd-6 analogs defined the chemical groups that are critical for its biologic activity. We thus introduce the first small-molecule PAMs for the β2AR, which may serve as a lead molecule for the development of novel therapeutics. The approach described in this work establishes a broadly applicable proof-of-concept strategy for affinity-based discovery of small-molecule allosteric compounds targeting unique conformational states of GPCRs.

A Type III-B Cmr effector complex catalyzes the synthesis of cyclic oligoadenylate second messengers by cooperative substrate binding.

Abstract: Recently, Type III-A CRISPR-Cas systems were found to catalyze the synthesis of cyclic oligoadenylates (cOAs), a second messenger that specifically activates Csm6, a Cas accessory RNase and confers antiviral defense in bacteria To test if III-B CRISPR-Cas systems could mediate a similar CRISPR signaling pathway, the Sulfolobus islandicus Cmr-α ribonucleoprotein complex (Cmr-α-RNP) was purified from the native host and tested for cOA synthesis We found that the system showed a robust production of cyclic tetra-adenylate (c-A4), and that c-A4 functions as a second messenger to activate the III-B-associated RNase Csx1 by binding to its CRISPR-associated Rossmann Fold domain Investigation of the kinetics of cOA synthesis revealed that Cmr-α-RNP displayed positively cooperative binding to the adenosine triphosphate (ATP) substrate Furthermore, mutagenesis of conserved domains in Cmr2α confirmed that, while Palm 2 hosts the active site of cOA synthesis, Palm 1 domain serves as the primary site in the enzyme-substrate interaction Together, our data suggest that the two Palm domains cooperatively interact with ATP molecules to achieve a robust cOA synthesis by the III-B CRISPR-Cas system

Structural basis for GPR40 allosteric agonism and incretin stimulation.

Abstract: Activation of free fatty acid receptor 1 (GPR40) by synthetic partial and full agonists occur via distinct allosteric sites. A crystal structure of GPR40-TAK-875 complex revealed the allosteric site for the partial agonist. Here we report the 2.76-A crystal structure of human GPR40 in complex with a synthetic full agonist, compound 1, bound to the second allosteric site. Unlike TAK-875, which acts as a Gαq-coupled partial agonist, compound 1 is a dual Gαq and Gαs-coupled full agonist. compound 1 binds in the lipid-rich region of the receptor near intracellular loop 2 (ICL2), in which the stabilization of ICL2 by the ligand is likely the primary mechanism for the enhanced G protein activities. The endogenous free fatty acid (FFA), γ-linolenic acid, can be computationally modeled in this site. Both γ-linolenic acid and compound 1 exhibit positive cooperativity with TAK-875, suggesting that this site could also serve as a FFA binding site.

Switching from Negative-Cooperativity to No-Cooperativity in the Binding of Ion-Pair Dimers by a Bis(calix[4]pyrrole) Macrocycle

Abstract: We report the synthesis of a macrocyclic receptor containing two di-meso-phenylcalix[4]pyrrole units linked by two triazole spacers. The 1,4-substitution of the 1,2,3-triazole spacers conveys different binding affinities to the two heteroditopic binding sites. These features make the receptor an ideal candidate to investigate allosteric cooperativity in the binding of ion-pair dimers. We probed the interaction of tetraalkylammonium salts (TBA·Cl, TBA·OCN, and MTOA·Cl) with the tetra-heterotopic macrocyclic receptor in chloroform solution using 1H NMR spectroscopic titration experiments. The results obtained show that, at millimolar concentration, the addition of 2 equiv of the salt to the receptor’s solution induced the quantitative pairwise binding of the ion-pairs. The 2:1 (ion-pair:receptor) complexes feature different binding geometries and binding cooperativities depending on the nature of the alkylammonium cation. The binding geometries assigned to the complexes of the ion-pair dimers in solution ar...

Functional roles of magnesium binding to extracellular signal-regulated kinase 2 explored by molecular dynamics simulations and principal component analysis.

Abstract: Molecular dynamics (MD) simulations coupled with principal component (PC) analysis were carried out to study functional roles of Mg2+ binding to extracellular signal-regulated kinase 2 (ERK2). The results suggest that Mg2+ binding heavily decreases eigenvalue of the first principal component and totally inhibits motion strength of ERK2, which favors stabilization of ERK2 structure. Binding free energy predictions indicate that Mg2+ binding produces an important effect on binding ability of adenosine triphosphate (ATP) to ERK2 and strengthens the ATP binding. The calculations of residue-based free energy decomposition show that lack of Mg2+ weakens interactions between the hydrophobic rings of ATP and five residues I29, V37, A50, L105, and L154. Hydrogen bond analyses also prove that Mg2+ binding increases occupancies of hydrogen bonds formed between ATP and residues K52, Q103, D104, and M106. We expect that this study can provide a significant theoretical hint for designs of anticancer drugs targeting ERK2.

Characterization of the binding interactions between EvaGreen dye and dsDNA.

Abstract: Understanding the dsDNA·EG binding interaction is important because the EvaGreen (EG) dye is increasingly used in real-time quantitative polymerase chain reaction, high resolution melting analysis, and routine quantification of DNA. In this work, a binding isotherm for the interactions of EG with duplex DNA (poly-dA17·poly-dT17) has been determined from the absorption and fluorescence spectra of the EG and dsDNA·EG complex. The isotherm has a sigmoidal shape and can be modeled with the Hill equation, indicating positive cooperativity for the binding interaction. A Scatchard plot of the binding data yields a concave-down curve in agreement with the Hill analysis of the binding isotherm for a positive cooperative binding interaction. Analysis of the Scatchard plot with the modified McGhee and von Hippel model for a finite one-dimensional homogeneous lattice and nonspecific binding of ligands to duplex DNA yields the intrinsic binding constant, the number of lattice sites occluded by a bound ligand, and the cooperativity parameter of 3.6 × 105 M−1, 4.0, and 8.1, respectively. The occluded site size of 4 indicates that moieties of the EG intercalate into the adjacent base pairs of the duplex DNA with a gap of 1 intercalation site between EG binding sites, as expected for a bifunctional molecule. Interestingly, at high [EG]/[base pair], the intercalation is disrupted. A model is proposed based on the fluorescence spectrum where the formation of anti-parallel stacked chains of EGs bound externally to the duplex DNA occur at these high ratios.

Probing the cooperative mechanism of the μ-δ opioid receptor heterodimer by multiscale simulation.

Abstract: Accumulated experimental evidence indicated that G-protein coupled receptors (GPCRs) could form biologically relevant oligomers and hetero-oligomers possess different functional properties from monomers and homo-oligomers, for example, unique pharmacology. However, the urgent lack of crystal structures of the GPCR oligomers results in very limited knowledge about their structural and functional mechanisms. In this work, we utilized a multiscale simulation strategy coupled with principal component analysis, correlation analysis and a protein structure network to study the hetero-dimerization of the μ-OR and δ-OR. We probed the cooperative mechanism involved in their activations, the allosteric communication pathways, the impact of the interface and differences from the μ-OR homodimer. The result indicates that TM1-TM2-H8 is a stable interface, but some residues of TM7 also participate in the dimer interface. Similar to the homodimer, the hetero-dimerization of the two inactive receptors would enhance the constitutive activation of one subunit but weaken that of the other subunit, both presenting a negative cooperativity. However, in contrast to the homodimer, the hetero-dimerization of the active protomer with the inactive one would weaken the constitutive activation of the inactive unit but maintain the activity of the active subunit. In addition, the hetero-dimerization and the activation of one subunit could significantly alter the types and the numbers of residues participating in the allosteric pathway from the ligand-binding pocket to the G-protein region and the pathway between two subunits. Some important residues were identified, which play important roles in modulating activations and cooperativity between two subunits. The observations from this work indicate that the negative cooperativity should be a common feature for the homodimers and the heterodimers, but the cooperative results would be significantly different between them, depending on the activated extent of one subunit.

Conservation of coactivator engagement mechanism enables small-molecule allosteric modulators.

Abstract: Transcriptional coactivators are a molecular recognition marvel because a single domain within these proteins, the activator binding domain or ABD, interacts with multiple compositionally diverse transcriptional activators. Also remarkable is the structural diversity among ABDs, which range from conformationally dynamic helical motifs to those with a stable core such as a β-barrel. A significant objective is to define conserved properties of ABDs that allow them to interact with disparate activator sequences. The ABD of the coactivator Med25 (activator interaction domain or AcID) is unique in that it contains secondary structural elements that are on both ends of the spectrum: helices and loops that display significant conformational mobility and a seven-stranded β-barrel core that is structurally rigid. Using biophysical approaches, we build a mechanistic model of how AcID forms binary and ternary complexes with three distinct activators; despite its static core, Med25 forms short-lived, conformationally mobile, and structurally distinct complexes with each of the cognate partners. Further, ternary complex formation is facilitated by allosteric communication between binding surfaces on opposing faces of the β-barrel. The model emerging suggests that the conformational shifts and cooperative binding is mediated by a flexible substructure comprised of two dynamic helices and flanking loops, indicating a conserved mechanistic model of activator engagement across ABDs. Targeting a region of this substructure with a small-molecule covalent cochaperone modulates ternary complex formation. Our data support a general strategy for the identification of allosteric small-molecule modulators of ABDs, which are key targets for mechanistic studies as well as therapeutic applications.

Molecular dynamics simulations of apo and holo forms of fatty acid binding protein 5 and cellular retinoic acid binding protein II reveal highly mobile protein, retinoic acid ligand, and water molecules

Abstract: Structural and dynamic properties from a series of 300 ns molecular dynamics, MD, simulations of two intracellular lipid binding proteins, iLBPs, (Fatty Acid Binding Protein 5, FABP5, and Cellular Retinoic Acid Binding Protein II, CRABP-II) in both the apo form and when bound with retinoic acid reveal a high degree of protein and ligand flexibility. The ratio of FABP5 to CRABP-II in a cell may determine whether it undergoes natural apoptosis or unrestricted cell growth in the presence of retinoic acid. As a result, FABP5 is a promising target for cancer therapy. The MD simulations presented here reveal distinct differences in the two proteins and provide insight into the binding mechanism. CRABP-II is a much larger, more flexible protein that closes upon ligand binding, where FABP5 transitions to an open state in the holo form. The traditional understanding obtained from crystal structures of the gap between two β-sheets of the β-barrel common to iLBPs and the α-helix cap that forms the portal to the binding pocket is insufficient for describing protein conformation (open vs. closed) or ligand entry and exit. When the high degree of mobility between multiple conformations of both the ligand and protein are examined via MD simulation, a new mode of ligand motion that improves understanding of binding dynamics is revealed.

Switching of the folding-energy landscape governs the allosteric activation of protein kinase A

Abstract: Protein kinases are dynamic molecular switches that sample multiple conformational states. The regulatory subunit of PKA harbors two cAMP-binding domains [cyclic nucleotide-binding (CNB) domains] that oscillate between inactive and active conformations dependent on cAMP binding. The cooperative binding of cAMP to the CNB domains activates an allosteric interaction network that enables PKA to progress from the inactive to active conformation, unleashing the activity of the catalytic subunit. Despite its importance in the regulation of many biological processes, the molecular mechanism responsible for the observed cooperativity during the activation of PKA remains unclear. Here, we use optical tweezers to probe the folding cooperativity and energetics of domain communication between the cAMP-binding domains in the apo state and bound to the catalytic subunit. Our study provides direct evidence of a switch in the folding-energy landscape of the two CNB domains from energetically independent in the apo state to highly cooperative and energetically coupled in the presence of the catalytic subunit. Moreover, we show that destabilizing mutational effects in one CNB domain efficiently propagate to the other and decrease the folding cooperativity between them. Taken together, our results provide a thermodynamic foundation for the conformational plasticity that enables protein kinases to adapt and respond to signaling molecules.

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

Abstract: 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.

Quantifying the binding stoichiometry and affinity of histo-blood group antigen oligosaccharides for human noroviruses.

Abstract: Human noroviruses (HuNoVs) are a major cause of acute gastroenteritis. Many HuNoVs recognize histo-blood group antigens (HBGAs) as cellular receptors or attachment factors for infection. It was recently proposed that HuNoV recognition of HBGAs involves a cooperative, multistep binding mechanism that exploits both known and previously unknown glycan binding sites. In this study, binding measurements, implemented using electrospray ionization mass spectrometry (ESI-MS) were performed on homodimers of the protruding domain (P dimers) of the capsid protein of three HuNoV strains [Saga (GII.4), Vietnam 026 (GII.10) and VA387 (GII.4)] with the ethyl glycoside of the B trisaccharide (α-d-Gal-(1→3)-[α-l-Fuc-(1→2)]-β-d-Gal-OC2H5) and free B type 1 tetrasaccharide (α-d-Gal-(1→3)-[α-l-Fuc-(1→2)]-β-d-Gal-(1→3)-d-GlcNAc) in an effort to confirm the existence of new HBGA binding sites. After correcting the mass spectra for nonspecific interactions that form in ESI droplets as they evaporate to dryness, all three P dimers were found to bind a maximum of two B trisaccharides at the highest concentrations investigated. The apparent affinities measured for stepwise binding of B trisaccharide suggest positive cooperativity. Similar results were obtained for B type 1 tetrasaccharide binding to Saga P dimer. Based on these results, it is proposed that HuNoV P dimers possess only two HBGA binding sites. It is also shown that nonspecific binding corrections applied to mass spectra acquired using energetic ion source conditions that promote in-source dissociation can lead to apparent HuNoV-HBGA oligosaccharide binding stoichiometries and affinities that are artificially high. Finally, evidence that high concentrations of oligosaccharide can induce conformational changes in HuNoV P dimers is presented.

Quantitative profiling of BATF family proteins/JUNB/IRF hetero-trimers using Spec-seq

Abstract: BATF family transcription factors (BATF, BATF2 and BATF3) form hetero-trimers with JUNB and either IRF4 or IRF8 to regulate cell fate in T cells and dendritic cells in vivo. While each combination of the hetero-trimer has a distinct role, some degree of cross-compensation was observed. The basis for the differential actions of IRF4 and IRF8 with BATF factors and JUNB is still unknown. We propose that the differences in function between these hetero-trimers may be caused by differences in their DNA binding preferences. While all three BATF family transcription factors have similar binding preferences when binding as a hetero-dimer with JUNB, the cooperative binding of IRF4 or IRF8 to the hetero-dimer/DNA complex could change the preferences. We used Spec-seq, which allows for the efficient and accurate determination of relative affinity to a large collection of sequences in parallel, to find differences between cooperative DNA binding of IRF4, IRF8 and BATF family members. We found that without IRF binding, all three hetero-dimer pairs exhibit nearly the same binding preferences to both expected wildtype binding sites TRE (TGA(C/G)TCA) and CRE (TGACGTCA). IRF4 and IRF8 show the very similar DNA binding preferences when binding with any of the three hetero-dimers. No major change of binding preferences was found in the half-sites between different hetero-trimers. IRF proteins bind with substantially lower affinity with either a single nucleotide spacer between IRF and BATF binding site or with an alternative mode of binding in the opposite orientation. In addition, the preference to CRE binding site was reduced with either IRF binding in all BATF–JUNB combinations. The specificities of BATF, BATF2 and BATF3 are all very similar as are their interactions with IRF4 and IRF8. IRF proteins binding adjacent to BATF sites increases affinity substantially compared to sequences with spacings between the sites, indicating cooperative binding through protein–protein interactions. The preference for the type of BATF binding site, TRE or CRE, is also altered when IRF proteins bind. These in vitro preferences aid in the understanding of in vivo binding activities.

Crystallographic snapshots of ligand binding to hexameric purine nucleoside phosphorylase and kinetic studies give insight into the mechanism of catalysis

Abstract: Purine nucleoside phosphorylase (PNP) catalyses the cleavage of the glycosidic bond of purine nucleosides using phosphate instead of water as a second substrate. PNP from Escherichia coli is a homohexamer, build as a trimer of dimers, and each subunit can be in two conformations, open or closed. This conformational change is induced by the presence of phosphate substrate, and very likely a required step for the catalysis. Closing one active site strongly affects the others, by a yet unclear mechanism and order of events. Kinetic and ligand binding studies show strong negative cooperativity between subunits. Here, for the first time, we managed to monitor the sequence of nucleoside binding to individual subunits in the crystal structures of the wild-type enzyme, showing that first the closed sites, not the open ones, are occupied by the nucleoside. However, two mutations within the active site, Asp204Ala/Arg217Ala, are enough not only to significantly reduce the effectiveness of the enzyme, but also reverse the sequence of the nucleoside binding. In the mutant the open sites, neighbours in a dimer of those in the closed conformation, are occupied as first. This demonstrates how important for the effective catalysis of Escherichia coli PNP is proper subunit cooperation.

Structural Basis of Sequential Allosteric Transitions in Tetrameric d-Lactate Dehydrogenases from Three Gram-Negative Bacteria

Abstract: d-Lactate dehydrogenases (d-LDHs) from Fusobacterium nucleatum (FnLDH) and Escherichia coli (EcLDH) exhibit positive cooperativity in substrate binding, and the Pseudomonas aeruginosa enzyme (PaLDH) shows negatively cooperative substrate binding. The apo and ternary complex structures of FnLDH and PaLDH have been determined together with the apo-EcLDH structure. The three enzymes consistently form homotetrameric structures with three symmetric axes, the P-, Q-, and R-axes, unlike Lactobacillus d-LDHs, P-axis-related dimeric enzymes, although apo-FnLDH and EcLDH form asymmetric and distorted quaternary structures. The tetrameric structure allows apo-FnLDH and EcLDH to form wide intersubunit contact surfaces between the opened catalytic domains of the two Q-axis-related subunits in coordination with their asymmetric and distorted quaternary structures. These contact surfaces comprise intersubunit hydrogen bonds and hydrophobic interactions and likely prevent the domain closure motion during initial substrate binding. In contrast, apo-PaLDH possesses a highly symmetrical quaternary structure and partially closed catalytic domains that are favorable for initial substrate binding and forms virtually no intersubunit contact surface between the catalytic domains, which present their negatively charged surfaces to each other at the subunit interface. Complex FnLDH and PaLDH possess highly symmetrical quaternary structures with closed forms of the catalytic domains, which are separate from each other at the subunit interface. Structure-based mutations successfully converted the three enzymes to their dimeric forms, which exhibited no significant cooperativity in substrate binding. These observations indicate that the three enzymes undergo typical sequential allosteric transitions to exhibit their distinctive allosteric functions through the tetrameric structures.

Detection of cooperatively bound transcription factor pairs using ChIP-seq peak intensities and expectation maximization.

Abstract: Transcription factors (TFs) often work cooperatively, where the binding of one TF to DNA enhances the binding affinity of a second TF to a nearby location. Such cooperative binding is important for activating gene expression from promoters and enhancers in both prokaryotic and eukaryotic cells. Existing methods to detect cooperative binding of a TF pair rely on analyzing the sequence that is bound. We propose a method that uses, instead, only ChIP-seq peak intensities and an expectation maximization (CPI-EM) algorithm. We validate our method using ChIP-seq data from cells where one of a pair of TFs under consideration has been genetically knocked out. Our algorithm relies on our observation that cooperative TF-TF binding is correlated with weak binding of one of the TFs, which we demonstrate in a variety of cell types, including E. coli, S. cerevisiae and M. musculus cells. We show that this method performs significantly better than a predictor based only on the ChIP-seq peak distance of the TFs under consideration. This suggests that peak intensities contain information that can help detect the cooperative binding of a TF pair. CPI-EM also outperforms an existing sequence-based algorithm in detecting cooperative binding. The CPI-EM algorithm is available at https://github.com/vishakad/cpi-em.

Crystal structure and pH-dependent allosteric regulation of human β-ureidopropionase, an enzyme involved in anticancer drug metabolism.

Abstract: β-Ureidopropionase catalyzes the third step of the reductive pyrimidine catabolic pathway responsible for breakdown of uracil, thymine and pyrimidine-based antimetabolites such as 5-fluorouracil. Nitrilase-like β-ureidopropionases use a tetrad of conserved residues (Cys233, Lys196, Glu119 and Glu207) for catalysis and occur in a variety of oligomeric states. Positive cooperativity towards the substrate N-carbamoyl-β-alanine and an oligomerization-dependent mechanism of substrate activation and product inhibition have been reported for the enzymes from some species but not others. Here, the activity of recombinant human β-ureidopropionase is shown to be similarly regulated by substrate and product, but in a pH-dependent manner. Existing as homodimer at pH 9, the enzyme increasingly associates to octamers and larger oligomers with decreasing pH. Only at physiological pH it is responsive to effector binding, with N-carbamoyl-β-alanine causing association to more active higher molecular mass species, and β-alanine dissociation to inactive dimers. The parallel between the pH and ligand-induced effects suggests protonation state changes to play a crucial role in the allosteric regulation mechanism. Disruption of dimer-dimer interfaces by site-directed mutagenesis generated dimeric, inactive enzyme variants. The crystal structure of the T299C variant refined to 2.08 A resolution revealed high structural conservation between human and fruit fly β-ureidopropionase, and supports the hypothesis that enzyme activation by oligomer assembly involves ordering of loop regions forming the entrance to the active site at the dimer-dimer interface, effectively positioning the catalytically important Glu207 in the active site.

Cu(II) EPR Reveals Two Distinct Binding Sites and Oligomerization of Innate Immune Protein Calgranulin C

Abstract: S100A12 or Calgranulin C is a homodimeric antimicrobial protein of the S100 family of EF-hand calcium-modulated proteins. S100A12 is involved in many diseases such as inflammation, tumor invasion, cancer and neurological disorders such as Alzheimer’s disease. The binding of transition metal ions to the protein is important as the sequestering of the metal ion induces conformational changes in the protein, inhibiting the growth of various pathogenic microorganisms. In this work, we probe the Cu2+ binding properties of Calgranulin C. We demonstrate that the two Cu2+ binding sites in Calgranulin C show different coordination environments in solution. Continuous wave-electron spin resonance (CW-ESR) spectra of Cu2+-bound protein clearly show two distinct components at higher Cu2+:protein ratios, which is indicative of the two different binding environments for the Cu2+ ions. The g|| and A|| values are also different for the two components, indicating that the number of directly coordinated nitrogen in each site differs. Furthermore, we perform CW-ESR titrations to obtain the binding affinity of the Ca2+-loaded protein to Cu2+ ions. We observe a positive cooperativity in binding of the two Cu2+ ions. To further probe the Cu2+ coordination, we also perform electron spin echo envelope modulation (ESEEM) experiment. We perform ESEEM at two different fields where one Cu2+ binding site dominates the other. At both sites we see distinct signatures of Cu2+–histidine coordination. However, we clearly see that the ESEEM spectra corresponding to the two Cu2+ binding sites are significantly different. There is clear change in the intensity of the double quantum peak with respect to the nuclear quadrupole interaction peak at the two different fields. Furthermore, ESEEM along with hyperfine sublevel correlation show that only one of the two Cu2+ binding sites has backbone coordination, confirming our previous observation. Finally, we perform double electron–electron resonance spectroscopy to probe if the difference in binding environment is due to the Cu2+ binding to different sites in the protein. We obtain a distance distribution with a sharp peak at ~ 3 nm and a broad peak at ~ 4 nm. The shorter distance agrees with the Cu2+–Cu2+ distance expected for a dimer from the crystal structure. The longer distance is consistent with the Cu2+–Cu2+ distance when oligomerization occurs.

Characterizing protein-DNA binding event subtypes in ChIP-exo data

Abstract: Regulatory proteins associate with the genome either by directly binding cognate DNA motifs or via protein-protein interactions with other regulators. Each genomic recruitment mechanism may be associated with distinct motifs, and may also result in distinct characteristic patterns in high-resolution protein-DNA binding assays. For example, the ChIP-exo protocol precisely characterizes protein-DNA crosslinking patterns by combining chromatin immunoprecipitation (ChIP) with 5′ → 3′ exonuclease digestion. Since different regulatory complexes will result in different protein-DNA crosslinking signatures, analysis of ChIP-exo sequencing tag patterns should enable detection of multiple protein-DNA binding modes for a given regulatory protein. However, current ChIP-exo analysis methods either treat all binding events as being of a uniform type, or rely on the presence of DNA motifs to cluster binding events into subtypes. To systematically detect multiple protein-DNA interaction modes in a single ChIP-exo experiment, we introduce the ChIP-exo mixture model (ChExMix). ChExMix probabilistically models the genomic locations and subtype membership of protein-DNA binding events using both ChIP-exo tag enrichment patterns and DNA sequence information, thus offering a principled and robust approach to characterizing binding subtypes in ChIP-exo data. We demonstrate that ChExMix achieves accurate detection and classification of binding event subtypes using in silico mixed ChIP-exo data. We further demonstrate the unique analysis abilities of ChExMix using a collection of ChIP-exo experiments that profile the binding of key transcription factors in MCF-7 cells. In these data, ChExMix detects cooperative binding interactions between FoxA1, ERα, and CTCF, thus demonstrating that ChExMix can effectively stratify ChIP-exo binding events into biologically meaningful subtypes. Availability: ChExMix is available from https://github.com/seqcode/chexmix

A Model for the Binding of Fluorescently Labeled Anti-Human CD4 Monoclonal Antibodies to CD4 Receptors on Human Lymphocytes

Abstract: The CD4 glycoprotein is a component of the T cell receptor complex which plays an important role in the human immune response. This manuscript describes the measurement and modeling of the binding of fluorescently labeled anti-human CD4 monoclonal antibodies (mAb; SK3 clone) to CD4 receptors on the surface of human peripheral blood mononuclear cells (PBMC). CD4 mAb fluorescein isothiocyanate (FITC) and CD4 mAb allophycoerythrin (APC) conjugates were obtained from commercial sources. Four binding conditions were performed, each with the same PBMC sample and different CD4 mAb conjugate. Each binding condition consisted of the PBMC sample incubated for 30 min in labeling solutions containing progressively larger concentrations of the CD4 mAb-label conjugate. After the incubation period, the cells were re-suspended in PBS-based buffer and analyzed using a flow cytometer to measure the mean fluorescence intensity (MFI) of the labeled cell populations. A model was developed to estimate the equilibrium concentration of bound CD4 mAb-label conjugates to CD4 receptors on PBMC. A set of parameters was obtained from the best fit of the model to the measured MFI data and the known number of CD4 receptors on PBMC surface. Divalent and monovalent binding had to be invoked for the APC and FITC CD4 mAb conjugates, respectively. This suggests that the mAb binding depends on the size of the label, which has significant implications for quantitative flow cytometry. The study supports the National Institute of Standards and Technology program to develop quantitative flow cytometry measurements.

Dynamic ensemble of HIV-1 RRE stem IIB reveals non-native conformations that disrupt the Rev binding site

Abstract: The HIV-1 Rev response element (RRE) RNA element mediates the nuclear export of intron containing viral RNAs by forming an oligomeric complex with the viral protein Rev. Stem IIB and nearby stem II three-way junction nucleate oligomerization through cooperative binding of two Rev molecules. Conformational flexibility at this RRE region has been shown to be important for Rev binding. However, the nature of the flexibility has remained elusive. Here, using NMR relaxation dispersion, including a new strategy for directly observing transient conformational states in large RNAs, we find that stem IIB alone or when part of the larger RREII three-way junction robustly exists in dynamic equilibrium with non-native excited state (ES) conformations that have a combined population of ~ 20 %. The ESs disrupt the Rev binding site by changing local secondary structure and their stabilization via point substitution mutations decreases the binding affinity to the Rev arginine-rich motif (ARM) by 15 to 80 fold. The ensemble clarifies the conformational flexibility observed in stem IIB, reveals long-range conformational coupling between stem IIB and the three-way junction that may play roles in cooperative Rev binding, and also identifies non-native RRE conformational states as new targets for the development of anti-HIV therapeutics.

Active regeneration unites high- and low-temperature features in cooperative self-assembly.

Abstract: Cytoskeletal filaments are capable of self-assembly in the absence of externally supplied chemical energy, but the rapid turnover rates essential for their biological function require a constant flux of adenosine triphosphate (ATP) or guanosine triphosphate (GTP) hydrolysis. The same is true for two-dimensional protein assemblies employed in the formation of vesicles from cellular membranes, which rely on ATP-hydrolyzing enzymes to rapidly disassemble upon completion of the process. Recent observations suggest that the nucleolus, p granules, and other three-dimensional membraneless organelles may also demand dissipation of chemical energy to maintain their fluidity. Cooperative binding plays a crucial role in the dynamics of these higher-dimensional structures, but is absent from classic models of one-dimensional cytoskeletal assembly. In this paper, we present a thermodynamically consistent model of active regeneration with cooperative assembly, and compute the maximum turnover rate and minimum disassembly time as a function of the chemical driving force and the binding energy. We find that these driven structures resemble different equilibrium states above and below the nucleation barrier. In particular, we show that the maximal acceleration under large binding energies unites infinite-temperature local fluctuations with low-temperature nucleation kinetics.

Negative Cooperative Binding of Thymidine, Ordered Substrate Binding, and Product Release of Human Mitochondrial Thymidine Kinase 2 Explain Its Complex Kinetic Properties and Physiological Functions

Abstract: Mitochondrial thymidine kinase 2 (TK2) catalyzes the phosphorylation of thymidine (dT) and deoxycytidine (dC) and is essential for mitochondrial function in post-mitotic tissues. The phosphorylatio ...