Transcription factors ERα and Sox2 have differing multiphasic DNA and RNA binding mechanisms

Many transcription factors (TFs) have been shown to bind RNA, leading to open questions regarding the mechanism(s) of this RNA binding and its role in regulating TF activities. Here we use biophysical assays to interrogate the k o n , k o f f , and K d for DNA and RNA binding of two model human transcription factors, ERα and Sox2. Unexpectedly, we found that both proteins exhibited multiphasic nucleic acid binding kinetics. We propose that Sox2 RNA and DNA multiphasic binding kinetics could be explained by a conventional model for sequential Sox2 monomer association and dissociation. In contrast, ERα nucleic acid binding exhibited biphasic dissociation paired with novel triphasic association behavior, where two apparent binding transitions are separated by a 10-20 min "lag" phase depending on protein concentration. We considered several conventional models for the observed kinetic behavior, none of which adequately explained all the ERα nucleic acid binding data. Instead, simulations with a model incorporating sequential ERα monomer association, ERα nucleic acid complex isomerization, and product "feedback" on isomerization rate recapitulated the general kinetic trends for both ERα DNA and RNA binding. Collectively, our findings reveal that Sox2 and ERα bind RNA and DNA with previously unappreciated multiphasic binding kinetics, and that their reaction mechanisms differ with ERα binding nucleic acids via a novel reaction mechanism.

Guardians of the Genome: How the Single-Stranded DNA-Binding Proteins RPA and CST Facilitate Telomere Replication

Telomeres act as the protective caps of eukaryotic linear chromosomes; thus, proper telomere maintenance is crucial for genome stability. Successful telomere replication is a cornerstone of telomere length regulation, but this process can be fraught due to the many intrinsic challenges telomeres pose to the replication machinery. In addition to the famous "end replication" problem due to the discontinuous nature of lagging strand synthesis, telomeres require various telomere-specific steps for maintaining the proper 3' overhang length. Bulk telomere replication also encounters its own difficulties as telomeres are prone to various forms of replication roadblocks. These roadblocks can result in an increase in replication stress that can cause replication forks to slow, stall, or become reversed. Ultimately, this leads to excess single-stranded DNA (ssDNA) that needs to be managed and protected for replication to continue and to prevent DNA damage and genome instability. RPA and CST are single-stranded DNA-binding protein complexes that play key roles in performing this task and help stabilize stalled forks for continued replication. The interplay between RPA and CST, their functions at telomeres during replication, and their specialized features for helping overcome replication stress at telomeres are the focus of this review.

The RNA-Binding Domain of hnRNP U Extends beyond the RGG/RG Motifs

Heterogeneous nuclear ribonucleoprotein U (hnRNP U) is a ubiquitously expressed protein that regulates chromatin architecture through its interactions with numerous DNA, protein, and RNA partners. The RNA-binding domain (RBD) of hnRNP U was previously mapped to an RGG/RG motif within its disordered C-terminal region, but little is understood about its binding mode and potential for selective RNA recognition. Analysis of publicly available hnRNP U enhanced UV cross-linking and immunoprecipitation (eCLIP) data identified high-confidence binding sites within human RNAs. We synthesized a set of diverse RNAs encompassing 11 of these identified cross-link sites for biochemical characterization using a combination of fluorescence anisotropy and electrophoretic mobility shift assays. These in vitro binding experiments with a rationally designed set of RNAs and hnRNP U domains revealed that the RGG/RG motif is a small part of a more expansive RBD that encompasses most of the disordered C-terminal region. This RBD contains a second, previously experimentally uncharacterized RGG/RG motif with RNA-binding properties comparable to those of the canonical RGG/RG motif. These RGG/RG motifs serve redundant functions, with neither serving as the primary RBD. While in isolation, each RGG/RG motif has modest affinity for RNA, together they significantly enhance the association of hnRNP U with RNA, enabling the binding of most of the designed RNA set with low to midnanomolar binding affinities. Identification and characterization of the complete hnRNP U RBD highlight the perils of a reductionist approach to defining biochemical activities in this system and pave the way for a detailed investigation of its RNA-binding specificity.

The RNA-binding domain of hnRNP U extends beyond the RGG/RG motifs

Heterogeneous nuclear ribonucleoprotein U (hnRNP U) is a ubiquitously expressed protein that regulates chromatin architecture through its interactions with numerous DNA, protein, and RNA partners. The RNA-binding domain (RBD) of hnRNP U was previously mapped to an RGG/RG element within its disordered C-terminal region, but little is understood about its binding mode and potential for selective RNA recognition. Analysis of publicly available hnRNP U enhanced UV crosslinking and immunoprecipitation (eCLIP) data identified high-confidence binding sites within human RNAs. We synthesized a set of diverse RNAs encompassing eleven of these identified crosslink sites for biochemical characterization using a combination of fluorescence anisotropy and electrophoretic mobility shift assays. These in vitro binding experiments with a rationally designed set of RNAs and hnRNP U domains revealed that the RGG/RG element is a small part of a more expansive RBD that encompasses most of the disordered C-terminal region. This RBD contains a second, previously experimentally uncharacterized RGG/RG element with RNA-binding properties comparable to the canonical RGG/RG element. These RGG/RG elements serve redundant functions, with neither serving as the primary RBD. While in isolation each RGG/RG element has modest affinity for RNA, together they significantly enhance the association of hnRNP U with RNA, enabling binding of most of the designed RNA set with low to mid-nanomolar binding affinities. Identification and characterization of the complete hnRNP U RBD highlights the perils of a reductionist approach to defining biochemical activities in this system and paves the way for a detailed investigation of its RNA-binding specificity.

The DNA binding high mobility group box protein family functionally binds RNA

Nucleic acid binding proteins regulate transcription, splicing, RNA stability, RNA localization, and translation, together tailoring gene expression in response to stimuli. Upon discovery, these proteins are typically classified as either DNA or RNA binding as defined by their in vivo functions; however, recent evidence suggests dual DNA and RNA binding by many of these proteins. High mobility group box (HMGB) proteins have a DNA binding HMGB domain, act as transcription factors and chromatin remodeling proteins, and are increasingly understood to interact with RNA as means to regulate gene expression. Herein, multiple layers of evidence that the HMGB family are dual DNA and RNA binding proteins is comprehensively reviewed. For example, HMGB proteins directly interact with RNA in vitro and in vivo, are localized to RNP granules involved in RNA processing, and their protein interactors are enriched in RNA binding proteins involved in RNA metabolism. Importantly, in cell-based systems, HMGB-RNA interactions facilitate protein-protein interactions, impact splicing outcomes, and modify HMGB protein genomic or cellular localization. Misregulation of these HMGB-RNA interactions are also likely involved in human disease. This review brings to light that as a family, HMGB proteins are likely to bind RNA which is essential to HMGB protein biology. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

Quantification of transcriptome changes to investigate the role of glucocorticoid receptor-RNA binding during dexamethasone treatment

OBJECTIVES

The glucocorticoid receptor (GR) is a well-studied, ligand-activated transcription factor and a common target of anti-inflammatory treatments. Recently, several studies have drawn attention the effects of binding of GR to RNA rather than DNA and the potential implications of this activity for GR function. The objective of our study was to further characterize the relationship between GR function and RNA binding by measuring changes in the glucocorticoid-driven transcriptome in the presence of a GR mutant that exhibited reduced RNA affinity.

DATA DESCRIPTION

GR was activated in three cell lines containing GR constructs (GR-HaloTag). One of the cell lines contained a wild-type GR-HaloTag. Another contained GR-HaloTag with a mutation that reduced RNA affinity and slightly reduced DNA affinity. The third cell line contained GR-HaloTag with a mutation that only slightly reduced DNA affinity. All three cell lines were treated with dexamethasone, a GR agonist. RNA-seq samples were collected every hour for 3 h. Moreover, transcriptome quantification was accomplished via labeling of RNAs transcribed in the final hour of dexamethasone treatment using 4-thiouridine. These labeled RNAs were then purified and sequenced. This data set is the first of its kind for GR and contains valuable insights into the function of RNA binding by GR.

RPA-like single-stranded DNA-binding protein complexes including CST serve as specialized processivity factors for polymerases

Telomeres and other single-stranded regions of the genome require specialized management to maintain stability and for proper progression of DNA metabolism pathways. Human Replication Protein A and CTC1-STN1-TEN1 are structurally similar heterotrimeric protein complexes that have essential ssDNA-binding roles in DNA replication, repair, and telomeres. Yeast and ciliates have related ssDNA-binding proteins with strikingly conserved structural features to these human heterotrimeric protein complexes. Recent breakthrough structures have extended our understanding of these commonalities by illuminating a common mechanism used by these proteins to act as processivity factors for their partner polymerases through their ability to manage ssDNA.

RNA binding by the glucocorticoid receptor attenuates dexamethasone-induced gene activation

The glucocorticoid receptor (GR) is a ligand-activated transcription factor that regulates a suite of genes through direct binding of GR to specific DNA promoter elements. GR also interacts with RNA, but the function of this RNA-binding activity remains elusive. Current models speculate that RNA could repress the transcriptional activity of GR. To investigate the function of the GR-RNA interaction on GR's transcriptional activity, we generated cells that stably express a mutant of GR with reduced RNA binding affinity and treated the cells with the GR agonist dexamethasone. Changes in the dexamethasone-driven transcriptome were quantified using 4-thiouridine labeling of RNAs followed by high-throughput sequencing. We find that while many genes are unaffected, GR-RNA binding is repressive for specific subsets of genes in both dexamethasone-dependent and independent contexts. Genes that are dexamethasone-dependent are activated directly by chromatin-bound GR, suggesting a competition-based repression mechanism in which increasing local concentrations of RNA may compete with DNA for binding to GR at sites of transcription. Unexpectedly, genes that are dexamethasone-independent instead display a localization to specific chromosomal regions, which points to changes in chromatin accessibility or architecture. These results show that RNA binding plays a fundamental role in regulating GR function and highlights potential functions for transcription factor-RNA interactions.

RPA engages telomeric G-quadruplexes more effectively than CST

G-quadruplexes (G4s) are a set of stable secondary structures that form within guanine-rich regions of single-stranded nucleic acids that pose challenges for DNA maintenance. The G-rich DNA sequence at telomeres has a propensity to form G4s of various topologies. The human protein complexes Replication Protein A (RPA) and CTC1-STN1-TEN1 (CST) are implicated in managing G4s at telomeres, leading to DNA unfolding and allowing telomere replication to proceed. Here, we use fluorescence anisotropy equilibrium binding measurements to determine the ability of these proteins to bind various telomeric G4s. We find that the ability of CST to specifically bind G-rich ssDNA is substantially inhibited by the presence of G4s. In contrast, RPA tightly binds telomeric G4s, showing negligible changes in affinity for G4 structure compared to linear ssDNAs. Using a mutagenesis strategy, we found that RPA DNA-binding domains work together for G4 binding, and simultaneous disruption of these domains reduces the affinity of RPA for G4 ssDNA. The relative inability of CST to disrupt G4s, combined with the greater cellular abundance of RPA, suggests that RPA could act as a primary protein complex responsible for resolving G4s at telomeres.

Replication Protein A Utilizes Differential Engagement of Its DNA-Binding Domains to Bind Biologically Relevant ssDNAs in Diverse Binding Modes

Replication protein A (RPA) is a ubiquitous ssDNA-binding protein that functions in many DNA processing pathways to maintain genome integrity. Recent studies suggest that RPA forms a highly dynamic complex with ssDNA that can engage with DNA in many modes that are orchestrated by the differential engagement of the four DNA-binding domains (DBDs) in RPA. To understand how these modes influence RPA interaction with biologically relevant ligands, we performed a comprehensive and systematic evaluation of RPA's binding to a diverse set of ssDNA ligands that varied in sequence, length, and structure. These equilibrium binding data show that WT RPA binds structured ssDNA ligands differently from its engagement with minimal ssDNAs. Next, we investigated each DBD's contributions to RPA's binding modes through mutation of conserved, functionally important aromatic residues. Mutations in DBD-A and -B have a much larger effect on binding when ssDNA is embedded into DNA secondary structures compared to their association with unstructured minimal ssDNA. As our data support a complex interplay of binding modes, it is critical to define the trimer core DBDs' role in binding these biologically relevant ligands. We found that DBD-C is important for engaging DNA with diverse binding modes, including, unexpectedly, at short ssDNAs. Thus, RPA uses its constituent DBDs to bind biologically diverse ligands in unanticipated ways. These findings lead to a better understanding of how RPA carries out its functions at diverse locations of the genome and suggest a mechanism through which dynamic recognition can impact differential downstream outcomes.

An Extended DNA Binding Domain of the Estrogen Receptor Alpha Directly Interacts with RNAs in Vitro

Estrogen receptor alpha (ERα) is a ligand-responsive transcription factor critical for sex determination and development. Recent reports challenge the canonical view of ERα function by suggesting an activity beyond binding dsDNA at estrogen-responsive promotor elements: association with RNAs in vivo. Whether these interactions are direct or indirect remains unknown, which limits the ability to understand the extent, specificity, and biological role of ERα-RNA binding. Here we demonstrate that an extended DNA-binding domain of ERα directly binds a wide range of RNAs in vitro with structural specificity. ERα binds RNAs that adopt a range of hairpin-derived structures independent of sequence, while interacting poorly with single- and double-stranded RNA. RNA affinities are only 4-fold weaker than consensus dsDNA and significantly tighter than nonconsensus dsDNA sequences. Moreover, RNA binding is competitive with DNA binding. Together, these data show that ERα utilizes an extended DNA-binding domain to achieve a high-affinity/low-specificity mode for interacting with RNA.

The DNA-Binding High-Mobility Group Box Domain of Sox Family Proteins Directly Interacts with RNA In Vitro

There is a growing body of evidence that a substantial number of protein domains identified as DNA-binding also interact with RNA to regulate biological processes. Several recent studies have revealed that the Sox2 transcription factor binds RNA through its high-mobility group box (HMGB) domain in vitro and in vivo. A high degree of conservation of this domain among members of the Sox family of transcription factors suggests that RNA-binding activity may be a general feature of these proteins. To address this hypothesis, we examined a subset of HMGB domains from human Sox family of proteins for their ability to bind both DNA and RNA in vitro. We observed selective, high-affinity interactions between Sox family HMGB domains and various model RNA elements, including a four-way junction RNA, a hairpin RNA with an internal bulge, G-quadruplex RNA, and a fragment of long noncoding RNA ES2, which is known to directly interact with Sox2. Importantly, the HMGB domains bind these RNA ligands significantly tighter than nonconsensus dsDNA and in some cases with affinities rivaling those of their consensus dsDNA sequences. These data suggest that RNA binding is a conserved feature of the Sox family of transcription factors with the potential to modulate unappreciated biological functions.

Cyp33 binds AU-rich RNA motifs via an extended interface that competitively disrupts the gene repressive Cyp33-MLL1 interaction in vitro

Cyp33 is an essential human cyclophilin prolyl isomerase that plays myriad roles in splicing and chromatin remodeling. In addition to a canonical cyclophilin (Cyp) domain, Cyp33 contains an RNA-recognition motif (RRM) domain, and RNA-binding triggers proline isomerase activity. One prominent role for Cyp33 is through a direct interaction with the mixed lineage leukemia protein 1 (MLL1, also known as KMT2A) complex, which is a histone methyltransferase that serves as a global regulator of human transcription. MLL activity is regulated by Cyp33, which isomerizes a key proline in the linker between the PHD3 and Bromo domains of MLL1, acting as a switch between gene activation and repression. The direct interaction between MLL1 and Cyp33 is critical, as deletion of the MLL1-PHD3 domain responsible for this interaction results in oncogenesis. The Cyp33 RRM is central to these activities, as it binds both the PHD3 domain and RNA. To better understand how RNA binding drives the action of Cyp33, we performed RNA-SELEX against full-length Cyp33 accompanied by deep sequencing. We have identified an enriched Cyp33 binding motif (AAUAAUAA) broadly represented in the cellular RNA pool as well as tightly binding RNA aptamers with affinities comparable and competitive with the Cyp33 MLL1-PHD3 interaction. RNA binding extends beyond the canonical RRM domain, but not to the Cyp domain, suggesting an indirect mechanism of interaction. NMR chemical shift mapping confirms an overlapping, but not identical, interface on Cyp33 for RNA and PHD3 binding. This finding suggests RNA can disrupt the gene repressive Cyp33-MLL1 complex providing another layer of regulation for chromatin remodeling by MLL1.

hnRNPK recognition of the B motif of Xist and other biological RNAs

Heterogeneous nuclear ribonuclear protein K (hnRNPK) is an abundant RNA-binding protein crucial for a wide variety of biological processes. While its binding preference for multi-cytosine-patch (C-patch) containing RNA is well documented, examination of binding to known cellular targets that contain C-patches reveals an unexpected breadth of binding affinities. Analysis of in-cell crosslinking data reinforces the notion that simple C-patch preference is not fully predictive of hnRNPK localization within transcripts. The individual RNA-binding domains of hnRNPK work together to interact with RNA tightly, with the KH3 domain being neither necessary nor sufficient for binding. Rather, the RG/RGG domain is implicated in providing essential contributions to RNA-binding, but not DNA-binding, affinity. hnRNPK is essential for X chromosome inactivation, where it interacts with Xist RNA specifically through the Xist B-repeat region. We use this interaction with an RNA motif derived from this B-repeat region to determine the RNA-structure dependence of C-patch recognition. While the location preferences of hnRNPK for C-patches are conformationally restricted within the hairpin, these structural constraints are relieved in the absence of RNA secondary structure. Together, these results illustrate how this multi-domain protein's ability to accommodate and yet discriminate between diverse cellular RNAs allows for its broad cellular functions.

The structure of human CST reveals a decameric assembly bound to telomeric DNA

The CTC1-STN1-TEN1 (CST) complex is essential for telomere maintenance and resolution of stalled replication forks genome-wide. Here, we report the 3.0-angstrom cryo-electron microscopy structure of human CST bound to telomeric single-stranded DNA (ssDNA), which assembles as a decameric supercomplex. The atomic model of the 134-kilodalton CTC1 subunit, built almost entirely de novo, reveals the overall architecture of CST and the DNA-binding anchor site. The carboxyl-terminal domain of STN1 interacts with CTC1 at two separate docking sites, allowing allosteric mediation of CST decamer assembly. Furthermore, ssDNA appears to staple two monomers to nucleate decamer assembly. CTC1 has stronger structural similarity to Replication Protein A than the expected similarity to yeast Cdc13. The decameric structure suggests that CST can organize ssDNA analogously to the nucleosome's organization of double-stranded DNA.

Spen links RNA-mediated endogenous retrovirus silencing and X chromosome inactivation

The Xist lncRNA mediates X chromosome inactivation (XCI). Here we show that Spen, an Xist-binding repressor protein essential for XCI , binds to ancient retroviral RNA, performing a surveillance role to recruit chromatin silencing machinery to these parasitic loci. Spen loss activates a subset of endogenous retroviral (ERV) elements in mouse embryonic stem cells, with gain of chromatin accessibility, active histone modifications, and ERV RNA transcription. Spen binds directly to ERV RNAs that show structural similarity to the A-repeat of Xist, a region critical for Xist-mediated gene silencing. ERV RNA and Xist A-repeat bind the RRM domains of Spen in a competitive manner. Insertion of an ERV into an A-repeat deficient Xist rescues binding of Xist RNA to Spen and results in strictly local gene silencing in cis. These results suggest that Xist may coopt transposable element RNA-protein interactions to repurpose powerful antiviral chromatin silencing machinery for sex chromosome dosage compensation.

The Sox2 transcription factor binds RNA

Certain transcription factors are proposed to form functional interactions with RNA to facilitate proper regulation of gene expression. Sox2, a transcription factor critical for maintenance of pluripotency and neurogenesis, has been found associated with several lncRNAs, although it is unknown whether these interactions are direct or via other proteins. Here we demonstrate that human Sox2 interacts directly with one of these lncRNAs with high affinity through its HMG DNA-binding domain in vitro. These interactions are primarily with double-stranded RNA in a non-sequence specific fashion, mediated by a similar but not identical interaction surface. We further determined that Sox2 directly binds RNA in mouse embryonic stem cells by UV-cross-linked immunoprecipitation of Sox2 and more than a thousand Sox2-RNA interactions in vivo were identified using fRIP-seq. Together, these data reveal that Sox2 employs a high-affinity/low-specificity paradigm for RNA binding in vitro and in vivo.

Some transcription factors have been proposed to functionally interact with RNA to facilitate proper regulation of gene expression. Here the authors demonstrate that human Sox2 interact directly and with high affinity to RNAs through its HMG DNA-binding domain.

Nonspecific Binding of RNA to PARP1 and PARP2 Does Not Lead to Catalytic Activation

Poly(ADP-ribose) polymerases 1 and 2 (PARP1 and PARP2, respectively), upon binding damaged DNA, become activated to add long chains of poly(ADP-ribose) (PAR) to themselves and other nuclear proteins. This activation is an essential part of the DNA damage response. The PAR modifications recruit the DNA repair machinery to sites of DNA damage and result in base excision and single-strand break repair, homologous recombination, nucleotide excision repair, and alternative nonhomologous end joining. More recently, both PARP1 and PARP2 have been shown to bind to or be activated by RNA, a property that could interfere with the function of PARP1 and PARP2 in the response to DNA damage or lead to necrosis by depletion of cellular NAD⁺. We have quantitatively evaluated the in vitro binding of a variety of RNAs to PARP1 and PARP2 and queried the ability of these RNAs to switch on enzymatic activity. We find that while both proteins bind RNAs without specificity toward sequence or structure, their interaction with RNA does not lead to auto-PARylation. Thus, although PARP1 and PARP2 are promiscuous with respect to activation by DNA, they both demonstrate exquisite selectivity against activation by RNA.

The glucocorticoid receptor DNA-binding domain recognizes RNA hairpin structures with high affinity

The glucocorticoid receptor (GR) binds the noncoding RNA Gas5 via its DNA-binding domain (DBD) with functional implications in pro-apoptosis signaling. Here, we report a comprehensive in vitro binding study where we have determined that GR-DBD is a robust structure-specific RNA-binding domain. GR-DBD binds to a diverse range of RNA hairpin motifs, both synthetic and biologically derived, with apparent mid-nanomolar affinity while discriminating against uniform dsRNA. As opposed to dimeric recognition of dsDNA, GR-DBD binds to RNA as a monomer and confers high affinity primarily through electrostatic contacts. GR-DBD adopts a discrete RNA-bound state, as assessed by NMR, distinct from both free and DNA-bound. NMR and alanine mutagenesis suggest a heightened involvement of the C-terminal α-helix of the GR-DBD in RNA-binding. RNA competes for binding with dsDNA and occurs in a similar affinity range as dimer binding to the canonical DNA element. Given the prevalence of RNA hairpins within the transcriptome, our findings strongly suggest that many RNAs have potential to impact GR biology.

Discrimination against RNA Backbones by a ssDNA Binding Protein

Pot1 is the shelterin component responsible for the protection of the single-stranded DNA (ssDNA) overhang at telomeres in nearly all eukaryotic organisms. The C-terminal domain of the DNA-binding domain, Pot1pC, exhibits non-specific ssDNA recognition, achieved through thermodynamically equivalent alternative binding conformations. Given this flexibility, it is unclear how specificity for ssDNA over RNA, an activity required for biological function, is achieved. Examination of the ribose-position specificity of Pot1pC shows that ssDNA specificity is additive but not uniformly distributed across the ligand. High-resolution structures of several Pot1pC complexes with RNA-DNA chimeric ligands reveal Pot1pC discriminates against RNA by utilizing non-compensatory binding modes that feature significant rearrangement of the binding interface. These alternative conformations, accessed through both ligand and protein flexibility, recover much, but not all, of the binding energy, leading to the observed reduction in affinities. These findings suggest that intermolecular interfaces are remarkably sophisticated in their tuning of specificity toward flexible ligands.

Diversity in peptide recognition by the SH2 domain of SH2B1

SH2B1 is a multidomain protein that serves as a key adaptor to regulate numerous cellular events, such as insulin, leptin, and growth hormone signaling pathways. Many of these protein-protein interactions are mediated by the SH2 domain of SH2B1, which recognizes ligands containing a phosphorylated tyrosine (pY), including peptides derived from janus kinase 2, insulin receptor, and insulin receptor substrate-1 and -2. Specificity for the SH2 domain of SH2B1 is conferred in these ligands either by a hydrophobic or an acidic side chain at the +3 position C-terminal to the pY. This specificity for chemically disparate species suggests that SH2B1 relies on distinct thermodynamic or structural mechanisms to bind to peptides. Using binding and structural strategies, we have identified unique thermodynamic signatures for each peptide binding mode, and several SH2B1 residues, including K575 and R578, that play distinct roles in peptide binding. The high-resolution structure of the SH2 domain of SH2B1 further reveals conformationally plastic protein loops that may contribute to the ability of the protein to recognize dissimilar ligands. Together, numerous hydrophobic and electrostatic interactions, in addition to backbone conformational flexibility, permit the recognition of diverse peptides by SH2B1. An understanding of this expanded peptide recognition will allow for the identification of novel physiologically relevant SH2B1/peptide interactions, which can contribute to the design of obesity and diabetes pharmaceuticals to target the ligand-binding interface of SH2B1 with high specificity.

Human CST Prefers G-Rich but Not Necessarily Telomeric Sequences

The human CST (CTC1-STN1-TEN1) heterotrimeric complex plays roles in both telomere maintenance and DNA replication through its ability to interact with single-stranded DNA (ssDNA) of a variety of sequences. The precise sequence specificity required to execute these functions is unknown. Telomere-binding proteins have been shown to specifically recognize key telomeric sequence motifs within ssDNA while accommodating nonspecifically recognized sequences through conformationally plastic interfaces. To better understand the role CST plays in these processes, we have produced a highly purified heterotrimer and elucidated the sequence requirements for CST recognition of ssDNA in vitro. CST discriminates against random sequence and binds a minimal ssDNA comprised of three repeats of telomeric sequence. Replacement of individual nucleotides with their complement reveals that guanines are specifically recognized in a largely additive fashion and that specificity is distributed uniformly throughout the ligand. Unexpectedly, adenosines are also well tolerated at these sites, but cytosines are disfavored. Furthermore, sequences unrelated to the telomere repeat, yet still G-rich, bind CST well. Thus, CST is not inherently telomere-specific, but rather is a G-rich sequence binder. This biochemical activity is reminiscent of the yeast t-RPA and Tetrahymena thermophila CST complexes and is consistent with roles at G-rich sites throughout the genome.

Multimodal Recognition of Diverse Peptides by the C-Terminal SH2 Domain of Phospholipase C-γ1 Protein

SH2 domains recognize phosphotyrosine (pY)-containing peptide ligands and play key roles in the regulation of receptor tyrosine kinase pathways. Each SH2 domain has individualized specificity, encoded in the amino acids neighboring the pY, for defined targets that convey their distinct functions. The C-terminal SH2 domain (PLCC) of the phospholipase C-γ1 full-length protein (PLCγ1) typically binds peptides containing small and hydrophobic amino acids adjacent to the pY, including a peptide derived from platelet-derived growth factor receptor B (PDGFRB) and an intraprotein recognition site (Y783 of PLCγ1) involved in the regulation of the protein's lipase activity. Remarkably, PLCC also recognizes unexpected peptides containing amino acids with polar or bulky side chains that deviate from this pattern. This versatility in recognition specificity may allow PLCγ1 to participate in diverse, previously unrecognized, signaling pathways in response to binding chemically dissimilar partners. We have used structural approaches, including nuclear magnetic resonance and X-ray crystallography, to elucidate the mechanisms of noncognate peptide binding to PLCC by ligands derived from receptor tyrosine kinase ErbB2 and from the insulin receptor. The high-resolution peptide-bound structures reveal that PLCC has a relatively static backbone but contains a chemically rich protein surface comprised of a combination of hydrophobic pockets and amino acids with charged side chains. We demonstrate that this expansive and chemically diverse PLCC interface, in addition to peptide conformational plasticity, permits PLCC to recognize specific noncognate peptide ligands with multimodal specificity.

NMR Chemical Shift Mapping of SH2 Peptide Interactions

Heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) experiments offer a rapid and high resolution approach to gaining binding and conformational insights into a protein-peptide interaction. By tracking ¹H and ¹⁵N chemical shift changes over the course of a peptide titration into isotopically labeled protein, amide NH pairs of amino acids whose chemical environment changes upon peptide binding can be identified. When mapped onto a structure of the protein, this approach can identify the peptide-binding interface or regions undergoing conformation changes within a protein upon ligand binding. Monitoring NMR chemical shift changes can also serve as a screening technique to identify novel interaction partners for a protein or to determine the binding affinity of a weak protein-peptide interaction. Here, we describe the application of NMR chemical shift mapping to the study of peptide binding to the C-terminal SH2 domain of PLCγ1.

Calorimetric Measurement of SH2 Domain Ligand Affinities

Isothermal titration calorimetry (ITC) has emerged as a leading approach in the characterization of protein/ligand interactions. This technique measures the heat change of a system upon binding of a ligand to a biomolecule, and thereby requires no immobilization, intrinsic fluorescence, or labeling of any kind of either species. If properly designed, a single experiment can not only measure the binding affinity, but also determine additional binding and thermodynamic parameters, including the enthalpy, entropy, and the stoichiometry of the interaction. Here, we describe the protocol for the collection of calorimetric data for the binding of peptides to SH2 protein domains.

Tying up the Ends: Plasticity in the Recognition of Single-Stranded DNA at Telomeres

Telomeres terminate nearly exclusively in single-stranded DNA (ssDNA) overhangs comprised of the G-rich 3' end. This overhang varies widely in length from species to species, ranging from just a few bases to several hundred nucleotides. These overhangs are not merely a remnant of DNA replication but rather are the result of complex further processing. Proper management of the telomeric overhang is required both to deter the action of the DNA damage machinery and to present the ends properly to the replicative enzyme telomerase. This Current Topic addresses the biochemical and structural features used by the proteins that manage these variable telomeric overhangs. The Pot1 protein tightly binds the single-stranded overhang, preventing DNA damage sensors from binding. Pot1 also orchestrates the access of telomerase to that same substrate. The remarkable plasticity of the binding interface exhibited by the Schizosaccharomyces pombe Pot1 provides mechanistic insight into how these roles may be accomplished, and disease-associated mutations clustered around the DNA-binding interface in the hPOT1 highlight the importance of this function. The budding yeast Cdc13-Stn1-Ten1, a telomeric RPA complex closely associated with telomere function, also interacts with ssDNA in a fashion that allows degenerate sequences to be recognized. A related human complex composed of hCTC1, hSTN1, and hTEN1 has recently emerged with links to both telomere maintenance and general DNA replication and also exhibits mutations associated with telomere pathologies. Overall, these sequence-specific ssDNA binders exhibit a range of recognition properties that allow them to perform their unique biological functions.

Telomere Replication Stress Induced by POT1 Inactivation Accelerates Tumorigenesis

Genome sequencing studies have revealed a number of cancer-associated mutations in the telomerebinding factor POT1. Here, we show that when combined with p53 deficiency, depletion of murine POT1a in common lymphoid progenitor cells fosters genetic instability, accelerates the onset, and increases the severity of T cell lymphomas. In parallel, we examined human and mouse cells carrying POT1 mutations found in cutaneous T cell lymphoma (CTCL) patients. Inhibition of POT1 activates ATRdependent DNA damage signaling and induces telomere fragility, replication fork stalling, and telomere elongation. Our data suggest that these phenotypes are linked to impaired CST (CTC1-STN1-TEN1) function at telomeres. Lastly, we show that proliferation of cancer cells lacking POT1 is enabled by the attenuation of the ATR kinase pathway. These results uncover a role for defective telomere replication during tumorigenesis.

Pinzaru et al. define a role for POT1 inactivation in the onset of thymic lymphomas. Inhibition of POT1 causes replication defects at telomeres resulting in telomere fragility, replication fork stalling, and genome instability. These results suggest a role of defective telemore replication during tumorigenesis

Less is more: structures of difficult targets with minimal constraints

By merging recent experimental and computational methodology advances, resolution-adapted structural recombination Rosetta has emerged as a powerful strategy for solving the structure of traditionally challenging targets. In this issue of Structure, Sgourakis and colleagues solve the structure of one such target, the immunoevasin protein m04, using this approach.

The telomeric protein Pot1 from Schizosaccharomyces pombe binds ssDNA in two modes with differing 3' end availability

Telomere protection and length regulation are important processes for aging, cancer and several other diseases. At the heart of these processes lies the single-stranded DNA (ssDNA)-binding protein Pot1, a component of the telomere maintenance complex shelterin, which is present in species ranging from fission yeast to humans. Pot1 contains a dual OB-fold DNA-binding domain (DBD) that fully confers its high affinity for telomeric ssDNA. Studies of S. pombe Pot1-DBD and its individual OB-fold domains revealed a complex non-additive behavior of the two OB-folds in the context of the complete Pot1 protein. This behavior includes the use of multiple distinct binding modes and an ability to form higher order complexes. Here we use NMR and biochemical techniques to investigate the structural features of the complete Pot1-DBD. These experiments reveal one binding mode characterized by only subtle alternations to the individual OB-fold subdomain structures, resulting in an inaccessible 3′ end of the ssDNA. The second binding mode, which has equivalent affinity, interacts differently with the 3′ end, rendering it available for interaction with other proteins. These findings suggest a structural switch that contributes to telomere end-protection and length regulation.

Single-stranded DNA-binding proteins: multiple domains for multiple functions

The recognition of single-stranded DNA (ssDNA) is integral to myriad cellular functions. In eukaryotes, ssDNA is present stably at the ends of chromosomes and at some promoter elements. Furthermore, it is formed transiently by several cellular processes including telomere synthesis, transcription, and DNA replication, recombination, and repair. To coordinate these diverse activities, a variety of proteins have evolved to bind ssDNA in a manner specific to their function. Here, we review the recognition of ssDNA through the analysis of high-resolution structures of proteins in complex with ssDNA. This functionally diverse set of proteins arises from a limited set of structural motifs that can be modified and arranged to achieve distinct activities, including a range of ligand specificities. We also investigate the ways in which these domains interact in the context of large multidomain proteins/complexes. These comparisons reveal the structural features that define the range of functions exhibited by these proteins.

The tenacious recognition of yeast telomere sequence by Cdc13 is fully exerted by a single OB-fold domain

Cdc13, the telomere end-binding protein from Saccharomyces cerevisiae, is a multidomain protein that specifically binds telomeric single-stranded DNA (ssDNA) with exquisitely high affinity to coordinate telomere maintenance. Recent structural and genetic data have led to the proposal that Cdc13 is the paralog of RPA70 within a telomere-specific RPA complex. Our understanding of Cdc13 structure and biochemistry has been largely restricted to studies of individual domains, precluding analysis of how each domain influences the activity of the others. To better facilitate a comparison to RPA70, we evaluated the ssDNA binding of full-length S. cerevisiae Cdc13 to its minimal substrate, Tel11. We found that, unlike RPA70 and the other known telomere end-binding proteins, the core Cdc13 ssDNA-binding activity is wholly contained within a single tight-binding oligosaccharide/oligonucleotide/oligopeptide binding (OB)-fold. Because two OB-folds are implicated in dimerization, we also evaluated the relationship between dimerization and ssDNA-binding activity and found that the two activities are independent. We also find that Cdc13 binding exhibits positive cooperativity that is independent of dimerization. This study reveals that, while Cdc13 and RPA70 share similar domain topologies, the corresponding domains have evolved different and specialized functions.

Nonspecific recognition is achieved in Pot1pC through the use of multiple binding modes

Pot1 is the protein responsible for binding to and protecting the 3' single-stranded DNA (ssDNA) overhang at most eukaryotic telomeres. Here, we present the crystal structure of one of the two oligonucleotide/oligosaccharide-binding folds (Pot1pC) that make up the ssDNA-binding domain in S. pombe Pot1. Comparison with the homologous human domain reveals unexpected structural divergence in the mode of ligand binding that explains the differing ligand requirements between species. Despite the presence of apparently base-specific hydrogen bonds, Pot1pC is able to bind a wide range of ssDNA sequences with thermodynamic equivalence. To address how Pot1pC binds ssDNA with little to no specificity, multiple structures of Pot1pC bound to noncognate ssDNA ligands were solved. These structures reveal that this promiscuity is implemented through new binding modes that thermodynamically compensate for base-substitutions through alternate stacking interactions and new H-bonding networks.

A small molecule inhibitor of Pot1 binding to telomeric DNA

Chromosome ends are complex structures, consisting of repetitive DNA sequence terminating in an ssDNA overhang with many associated proteins. Because alteration of the regulation of these ends is a hallmark of cancer, telomeres and telomere maintenance have been prime drug targets. The universally conserved ssDNA overhang is sequence-specifically bound and regulated by Pot1 (protection of telomeres 1), and perturbation of Pot1 function has deleterious effects for proliferating cells. The specificity of the Pot1/ssDNA interaction and the key involvement of this protein in telomere maintenance have suggested directed inhibition of Pot1/ssDNA binding as an efficient means of disrupting telomere function. To explore this idea, we developed a high-throughput time-resolved fluorescence resonance energy transfer (TR-FRET) screen for inhibitors of Pot1/ssDNA interaction. We conducted this screen with the DNA-binding subdomain of Schizosaccharomyces pombe Pot1 (Pot1pN), which confers the vast majority of Pot1 sequence-specificity and is highly similar to the first domain of human Pot1 (hPOT1). Screening a library of ∼20 000 compounds yielded a single inhibitor, which we found interacted tightly with sub-micromolar affinity. Furthermore, this compound, subsequently identified as the bis-azo dye Congo red (CR), was able to competitively inhibit hPOT1 binding to telomeric DNA. Isothermal titration calorimetry and NMR chemical shift analysis suggest that CR interacts specifically with the ssDNA-binding cleft of Pot1, and that alteration of this surface disrupts CR binding. The identification of a specific inhibitor of ssDNA interaction establishes a new pathway for targeted telomere disruption.

Telomerase and telomere-associated proteins: structural insights into mechanism and evolution

Recent advances in our structural understanding of telomerase and telomere-associated proteins have contributed significantly to elucidating the molecular mechanisms of telomere maintenance. The structures of telomerase TERT domains have provided valuable insights into how experimentally identified conserved motifs contribute to the telomerase reverse transcriptase reaction. Additionally, structures of telomere-associated proteins in a variety of organisms have revealed that, across evolution, telomere-maintenance mechanisms employ common structural elements. For example, the single-stranded 3' overhang of telomeric DNA is specifically and tightly bound by an OB-fold in nearly all species, including ciliates (TEBP and Pot1a), fission yeast (SpPot1), budding yeast (Cdc13), and humans (hPOT1). Structures of the yeast Cdc13, Stn1, and Ten1 proteins demonstrated that telomere maintenance is regulated by a complex that bears significant similarity to the RPA heterotrimer. Similarly, proteins that specifically bind double-stranded telomeric DNA in divergent species use homeodomains to execute their functions (human TRF1 and TRF2 and budding yeast ScRap1). Likewise, the conserved protein Rap1, which is found in budding yeast, fission yeast, and humans, contains a structural motif that is known to be critical for protein-protein interaction. In addition to revealing the common underlying themes of telomere maintenance, structures have also elucidated the specific mechanisms by which many of these proteins function, including identifying a telomere-specific domain in Stn1 and how the human TRF proteins avoid heterodimerization. In this review, we summarize the high-resolution structures of telomerase and telomere-associated proteins and discuss the emergent common structural themes among these proteins. We also address how these high-resolution structures complement biochemical and cellular studies to enhance our understanding of telomere maintenance and function.

The role of adjuvant in mediating antigen structure and stability

The purpose of this study was to probe the fate of a model antigen, a cysteine-free mutant of bacteriophage T4 lysozyme, to the level of fine structural detail, as a consequence of its interaction with an aluminum (Al)-containing adjuvant. Fluorescence spectroscopy and differential scanning calorimetry were used to compare the thermal stability of the protein in solution versus adsorbed onto an Al-containing adjuvant. Differences in accessible hydrophobic surface areas were investigated using an extrinsic fluorescence probe, 8-Anilino-1-naphthalenesulfonic acid (ANS). As has been observed with other model antigens, the apparent thermal stability of the protein decreased following adsorption onto the adjuvant. ANS spectra suggested that adsorption onto the adjuvant caused an increase in exposure of hydrophobic regions of the protein. Electrostatic interactions drove the adsorption, and disruption of these interactions with high ionic strength buffers facilitated the collection of two-dimensional (15) N heteronuclear single quantum coherence nuclear magnetic resonance data of protein released from the adjuvant. Although the altered stability of the adsorbed protein suggested changes to the protein's structure, the fine structure of the desorbed protein was nearly identical to the protein's structure in the adjuvant-free formulation. Thus, the adjuvant-induced changes to the protein that were responsible for the reduced thermal stability were not observed upon desorption.

Schizosaccharomyces pombe protection of telomeres 1 utilizes alternate binding modes to accommodate different telomeric sequences

The ends of eukaryotic chromosomes consist of long tracts of repetitive GT-rich DNA with variable sequence homogeneity between and within organisms. Telomeres terminate in a conserved 3'-ssDNA overhang that, regardless of sequence variability, is specifically and tightly bound by proteins of the telomere-end protection family. The high affinity ssDNA-binding activity of S. pombe Pot1 protein (SpPot1) is conferred by a DNA-binding domain consisting of two subdomains, Pot1pN and Pot1pC. Previous work has shown that Pot1pN binds a single repeat of the core telomere sequence (GGTTAC) with exquisite specificity, while Pot1pC binds an extended sequence of nine nucleotides (GGTTACGGT) with modest specificity requirements. We find that full-length SpPot1 binds the composite 15mer, (GGTTAC)(2)GGT, and a shorter two-repeat 12mer, (GGTTAC)(2), with equally high affinity (<3 pM), but with substantially different kinetic and thermodynamic properties. The binding mode of the SpPot1/15mer complex is more stable than that of the 12mer complex, with a 2-fold longer half-life and increased tolerance to nucleotide and amino acid substitutions. Our data suggest that SpPot1 protection of heterogeneous telomeres is mediated through 5'-sequence recognition and the use of alternate binding modes to maintain high affinity interaction with the G-strand, while simultaneously discriminating against the complementary strand.

Sequence-specific binding to telomeric DNA is not a conserved property of the Cdc13 DNA binding domain

In the budding yeast Saccharomyces cerevisiae, chromosome end protection is provided by a heterotrimeric complex composed of Cdc13 in association with the RPA-like proteins Stn1 and Ten1. We report here that the high affinity and specificity of the S. cerevisiae Cdc13 DNA binding domain for single-stranded telomeric DNA are not widely shared by other fungal Cdc13 proteins, suggesting that restriction of this complex to telomeres may be limited to the Saccharomyces clade. We propose that the evolutionarily conserved task of Stn1 and Ten1 (and their associated large subunit) is a genome-wide role in DNA replication rather than a telomere-dedicated activity.

Effect of benzyl alcohol on recombinant human interleukin-1 receptor antagonist structure and hydrogen-deuterium exchange

Benzyl alcohol, a preservative commonly added to multidose therapeutic protein formulations, can accelerate aggregation of recombinant human interleukin-1 receptor antagonist (rhIL-1ra). To investigate the interactions between benzyl alcohol and rhIL-1ra, we used nuclear magnetic resonance to observe the effect of benzyl alcohol on the chemical shifts of amide resonances of rhIL-1ra and to measure hydrogen-deuterium exchange rates of individual rhIL-1ra residues. Addition of 0.9% benzyl alcohol caused significant chemical shifts of amide resonances for residues 90-97, suggesting that these solvent-exposed residues participate in the binding of benzyl alcohol. In contrast, little perturbation of exchange rates was observed in the presence of either sucrose or benzyl alcohol.

Electron-tunneling pathways in cytochrome C

Distant Fe(2+)-Ru(3+) electronic couplings have been extracted from intramolecular electrontransfer rates in Ru(histidine(x)) (where X = 33, 39, 62, and 72) derivatives of cytochrome c. The couplings increase according to 62 (0.0060) < 72 (0.057) < 33 (0.097) < 39 (0.11 per wave numbers); however, this order is out of line with the histidine to heme edge-edge distances [62 (14.8) > 39 (12.3) > 33 (11.1) > 72 (8.4 angstroms)]. The rates (and the couplings) correlate with the lengths of sigma-tunneling pathways comprised of covalent bonds, hydrogen bonds, and through-space jumps from the histidines to the heme group. Space jumps greatly decrease couplings: One from Pro(71) to Met(80) extends the sigma-tunneling length of the His(72) pathway by roughly 10 covalent-bond units.

Investigating the role of the Est3 protein in yeast telomere replication

The Est3 subunit of yeast telomerase, which adopts a predicted OB-fold, is essential for telomere replication. To assess the possible contributions that Est3 might make to enzyme catalysis, we compared telomerase activity from wild type and est3-Delta strains of Saccharomyces castellii, which revealed that loss of the Est3 subunit results in a 2- to 3-fold decline in nucleotide addition. This effect was not primer-specific, based on assessment of a panel of primers that spanned the template of the S. castellii telomerase RNA. Furthermore, using nuclear magnetic resonance chemical shift perturbation, no chemical shift change was observed at any site in the protein upon addition of single-stranded DNA, arguing against a role for Est3 in recognition of telomeric substrates by telomerase. Addition of exogenous Est3 protein, including mutant Est3 proteins that are severely impaired for telomere replication in vivo, fully restored activity in est3-Delta telomerase reactions. Thus, Est3 performs an in vivo regulatory function in telomere replication, which is distinct from any potential contribution that Est3 might make to telomerase activity.

Nonadditivity in the recognition of single-stranded DNA by the schizosaccharomyces pombe protection of telomeres 1 DNA-binding domain, Pot1-DBD

The Schizosaccharomyces pombe protection of telomeres 1 (SpPot1) protein recognizes the 3' single-stranded ends of telomeres and provides essential protective and regulatory functions. The ssDNA-binding activity of SpPot1 is conferred by its ssDNA-binding domain, Pot1-DBD (residues 1-389), which can be further separated into two distinct domains, Pot1pN (residues 1-187) and Pot1pC (residues 188-389). Here we show that Pot1pC, like Pot1pN, can function independently of Pot1-DBD and binds specifically to a minimal nonameric oligonucleotide, d(GGTTACGGT), with a K(D) of 400 +/- 70 nM (specifically recognized nucleotides in bold). NMR chemical shift perturbation analysis indicates that the overall structures of the isolated Pot1pN and Pot1pC domains remain intact in Pot1-DBD. Furthermore, alanine scanning reveals modest differences in the ssDNA-binding contacts provided by isolated Pot1pN and within Pot1-DBD. Although the global character of both Pot1pN and Pot1pC is maintained in Pot1-DBD, chemical shift perturbation analysis highlights localized structural differences within the G1/G2 and T3/T4 binding pockets of Pot1pN in Pot1-DBD, which correlate with its distinct ssDNA-binding activity. Furthermore, we find evidence for a putative interdomain interface on Pot1pN that mediates interactions with Pot1pC that ultimately result in the altered ssDNA-binding activity of Pot1-DBD. Together, these data provide insight into the mechanisms underlying the activity and regulation of SpPot1 at the telomere.

Inhibition of yeast telomerase action by the telomeric ssDNA-binding protein, Cdc13p

Appropriate control of the chromosome end-replicating enzyme telomerase is crucial for maintaining telomere length and genomic stability. The essential telomeric DNA-binding protein Cdc13p both positively and negatively regulates telomere length in budding yeast. Here we test the effect of purified Cdc13p on telomerase action in vitro. We show that the full-length protein and its DNA-binding domain (DBD) inhibit primer extension by telomerase. This inhibition occurs by competitive blocking of telomerase access to DNA. To further understand the requirements for productive telomerase 3′-end access when Cdc13p or the DBD is bound to a telomerase substrate, we constrained protein binding at various distances from the 3′-end on two sets of increasingly longer oligonucleotides. We find that Cdc13p inhibits the action of telomerase through three distinct biochemical modes, including inhibiting telomerase even when a significant tail is available, representing a novel ‘action at a distance’ inhibitory activity. Thus, while yeast Cdc13p exhibits the same general activity as human POT1, providing an off switch for telomerase when bound near the 3′-end, there are significant mechanistic differences in the ways telomere end-binding proteins inhibit telomerase action.

High concentration formulations of recombinant human interleukin-1 receptor antagonist: I. Physical characterization

At relatively high protein concentrations (i.e., up to 100 mg/mL), recombinant human interleukin-1 receptor antagonist (rhIL-1ra) was found to exist in a monomer-dimer equilibrium controlled by solution ionic strength. Sedimentation equilibrium at 25 degrees C was used to measure the increase in the dimer dissociation constant (K(d)) as a function of ionic strength. K(d) increased from 2.0 to 12.6 mM as the solution ionic strength was increased from 0.011 to 0.184 molal. These K(d) values were used with both static light scattering and membrane osmometry data collected over a protein concentration range of 1-100 mg/mL to determine second osmotic virial coefficients. Expanding the second osmotic virial coefficient model to account for separate monomer-monomer (B(22)), monomer-dimer (B(23)), and dimer-dimer (B(33)) interactions reveals net monomer-dimer interactions are attractive, whereas the others are repulsive. Lastly, isothermal titration calorimetry dilution experiments showed that rhIL-1ra dimerization is enthalpically driven (DeltaH(dimerization) << 0), which is consistent with intermolecular cation-pi interactions previously proposed as the monomer-monomer contact sites in dimers.

Accurate structural correlations from maximum likelihood superpositions

The cores of globular proteins are densely packed, resulting in complicated networks of structural interactions. These interactions in turn give rise to dynamic structural correlations over a wide range of time scales. Accurate analysis of these complex correlations is crucial for understanding biomolecular mechanisms and for relating structure to function. Here we report a highly accurate technique for inferring the major modes of structural correlation in macromolecules using likelihood-based statistical analysis of sets of structures. This method is generally applicable to any ensemble of related molecules, including families of nuclear magnetic resonance (NMR) models, different crystal forms of a protein, and structural alignments of homologous proteins, as well as molecular dynamics trajectories. Dominant modes of structural correlation are determined using principal components analysis (PCA) of the maximum likelihood estimate of the correlation matrix. The correlations we identify are inherently independent of the statistical uncertainty and dynamic heterogeneity associated with the structural coordinates. We additionally present an easily interpretable method (“PCA plots”) for displaying these positional correlations by color-coding them onto a macromolecular structure. Maximum likelihood PCA of structural superpositions, and the structural PCA plots that illustrate the results, will facilitate the accurate determination of dynamic structural correlations analyzed in diverse fields of structural biology.

Biological macromolecules comprise extensive networks of interconnected atoms. These complex coupled networks result in correlated structural dynamics, where atoms and residues move and evolve together as concerted conformational changes. The availability of a wealth of macromolecular structures necessitates the use of robust strategies for analyzing the correlated modes of motion found in molecular ensembles. Current strategies use a combination of least-squares superpositions and statistical analysis of the structural covariance matrix. However, the least-squares treatment implicitly requires that atoms are uncorrelated and that each atom has the same positional uncertainty, two assumptions which are violated in structural ensembles. For example, the atoms in the proteins are connected by chemical bonds, covalent and non-covalent, resulting in strong correlations. Furthermore, different atoms have different variances, because some atoms are known with less precision or have greater mobility. Using maximum likelihood (ML) analysis, we have developed a technique that is markedly more accurate than the classical least-squares approach by accounting for both correlations and heterogeneous variances. The improved ability to accurately analyze the major modes of dynamic structural correlations will benefit a diverse range of biological disciplines, including nuclear magnetic resonance (NMR) spectroscopy, crystallography, molecular dynamics, and molecular evolution.

Deciphering the mechanism of thermodynamic accommodation of telomeric oligonucleotide sequences by the Schizosaccharomyces pombe protection of telomeres 1 (Pot1pN) protein

Linear chromosomes terminate in specialized nucleoprotein structures called telomeres, which are required for genomic stability and cellular proliferation. Telomeres end in an unusual 3' single-strand overhang that requires a special capping mechanism to prevent inappropriate recognition by the DNA damage machinery. In Schizosaccharomyces pombe, this protective function is mediated by the Pot1 protein, which binds specifically and with high affinity to telomeric ssDNA. We have characterized the thermodynamics and accommodation of both cognate and noncognate telomeric single-stranded DNA (ssDNA) sequences by Pot1pN, an autonomous ssDNA-binding domain (residues 1-187) found in full-length S. pombe Pot1. Direct calorimetric measurements of cognate telomeric ssDNA binding to Pot1pN show favorable enthalpy, unfavorable entropy, and a negative heat-capacity change. Thermodynamic analysis of the binding of noncognate telomeric ssDNA to Pot1pN resulted in unexpected changes in free energy, enthalpy, and entropy. Chemical-shift perturbation and structural analysis of these bound noncognate sequences show that these thermodynamic changes result from the structural rearrangement of both Pot1pN and the bound oligonucleotide. These data suggest that the ssDNA-binding interface is highly dynamic and, in addition to the conformation observed in the crystal structure of the Pot1pN/d(GGTTAC) complex, capable of adopting alternative thermodynamically equivalent conformations.

Probing the mechanism of recognition of ssDNA by the Cdc13-DBD

The Saccharomyces cerevisiae protein Cdc13 tightly and specifically binds the conserved G-rich single-stranded overhang at telomeres and plays an essential role in telomere end-protection and length regulation. The 200 residue DNA-binding domain of Cdc13 (Cdc13-DBD) binds an 11mer single-stranded representative of the yeast telomeric sequence [Tel11, d(GTGTGGGTGTG)] with a 3 pM affinity and specificity for three bases (underlined) at the 5′ end. The structure of the Cdc13-DBD bound to Tel11 revealed a large, predominantly aromatic protein interface with several unusual features. The DNA adopts an irregular, extended structure, and the binding interface includes a long (∼30 amino acids) structured loop between strands β2-β3 (L_2–3) of an OB-fold. To investigate the mechanism of ssDNA binding, we studied the free and bound states of Cdc13-DBD using NMR spectroscopy. Chemical shift changes indicate that the basic topology of the domain, including L_2–3, is essentially intact in the free state. Changes in slow and intermediate time scale dynamics, however, occur in L_2–3, while conformational changes distant from the DNA interface suggest an induced fit mechanism for binding in the ‘hot spot’ for binding affinity and specificity. These data point to an overall binding mechanism well adapted to the heterogeneous nature of yeast telomeres.

Characterization of protein-protein interactions critical for poliovirus replication: analysis of 3AB and VPg binding to the RNA-dependent RNA polymerase

Two critical interactions within the poliovirus RNA replication complex are those of the RNA-dependent RNA polymerase 3D with the viral proteins 3AB and VPg. 3AB is a membrane-binding protein responsible for the localization of the polymerase to the membranous vesicles at which replication occurs. VPg (a peptide comprising the 3B region of 3AB) is the 22-residue soluble product of 3AB cleavage and serves as the protein primer for RNA replication. The detailed interactions of these proteins with the RNA-dependent RNA polymerase 3D were analyzed to elucidate the precise roles of 3AB and VPg in the viral RNA replication complex. Using a membrane-based pull-down assay, we have identified a binding "hot-spot" spanning residues 100 to 104 in the 3B (VPg) region of 3AB which plays a critical role in mediating the interaction of 3AB with the polymerase. Isothermal titration calorimetry shows that the interaction of VPg with 3D is enthalpically driven, with a dissociation constant of 11 microM. Mutational analyses of VPg indicate that a subset of the residues important for 3AB-3D binding are also important for VPg-3D binding. Two residues in particular, P14 and R17, were shown to be absolutely critical for the binding interaction. This work provides the direct characterization of two binding interactions critical for the replication of this important class of viruses and identifies a conserved polymerase binding sequence responsible for targeting the polymerase.

Empirical Bayes hierarchical models for regularizing maximum likelihood estimation in the matrix Gaussian Procrustes problem

Procrustes analysis involves finding the optimal superposition of two or more "forms" via rotations, translations, and scalings. Procrustes problems arise in a wide range of scientific disciplines, especially when the geometrical shapes of objects are compared, contrasted, and analyzed. Classically, the optimal transformations are found by minimizing the sum of the squared distances between corresponding points in the forms. Despite its widespread use, the ordinary unweighted least-squares (LS) criterion can give erroneous solutions when the errors have heterogeneous variances (heteroscedasticity) or the errors are correlated, both common occurrences with real data. In contrast, maximum likelihood (ML) estimation can provide accurate and consistent statistical estimates in the presence of both heteroscedasticity and correlation. Here we provide a complete solution to the nonisotropic ML Procrustes problem assuming a matrix Gaussian distribution with factored covariances. Our analysis generalizes, simplifies, and extends results from previous discussions of the ML Procrustes problem. An iterative algorithm is presented for the simultaneous, numerical determination of the ML solutions.