Coupling of local folding to site-specific binding of proteins to DNA

Unveiling the multifaceted interactions of antitumor drug mitoxantrone with ct-DNA through biophysical and in silico studies

Mitoxantrone, an anthraquinone derivative, is a widely used anticancer drug with its well-known ability to engage in complex interactions with DNA. Although known for its intercalating ability, the enigma surrounding its binding modes with DNA persists. The existing corpus of literature primarily focuses on mitoxantrone-DNA interactions with short DNA sequences, thereby yielding insights into its interactive nature is limited to this specific sequence. This study aims to elucidate the diverse modes with which mitoxantrone interacts with calf thymus DNA using a combination of spectroscopy, calorimetry and in silico studies. The findings from spectroscopic, calorimetric and molecular dynamic results in correlation with existing literature, unveil a fascinating narrative: mitoxantrone intercalates at lower concentrations but promotes condensation at higher concentrations. Although intercalation with side chains positioned in the minor/major groove is the major binding mode in GC-rich sequences, molecular modelling studies hint at an alternative binding mode in AT-rich sequences where it exclusively displays pure electrostatic interaction. These findings underscore the pivotal role of both drug structure and base sequence in dictating binding mode and affinity. Such insights not only deepen the understanding of structure-activity relationships but also hold promise for guiding future drug design strategies.

Insights into Ligand-Mediated Activation of an Oligomeric Ring-Shaped Gene-Regulatory Protein from Solution- and Solid-State NMR

The 91 kDa oligomeric ring-shaped ligand binding protein TRAP (trp RNA binding attenuation protein) regulates the expression of a series of genes involved in tryptophan (Trp) biosynthesis in bacilli. When cellular Trp levels rise, the free amino acid binds to sites buried in the interfaces between each of the 11 (or 12, depending on the species) protomers in the ring. Crystal structures of Trp-bound TRAP show the Trp ligands are sequestered from solvent by a pair of loops from adjacent protomers that bury the bound ligand via polar contacts to several threonine residues. Binding of the Trp ligands occurs cooperatively, such that successive binding events occur with higher apparent affinity but the structural basis for this cooperativity is poorly understood. We used solution methyl-TROSY NMR relaxation experiments focused on threonine and isoleucine sidechains, as well as magic angle spinning solid-state NMR 13C-13C and 15N-13C chemical shift correlation spectra on uniformly labeled samples recorded at 800 and 1200 MHz, to characterize the structure and dynamics of the protein. Methyl 13C relaxation dispersion experiments on ligand-free apo TRAP revealed concerted exchange dynamics on the µs-ms time scale, consistent with transient sampling of conformations that could allow ligand binding. Cross-correlated relaxation experiments revealed widespread disorder on fast timescales. Chemical shifts for methyl-bearing side chains in apo- and Trp-bound TRAP revealed subtle changes in the distribution of sampled sidechain rotameric states. These observations reveal a pathway and mechanism for induced conformational changes to generate homotropic Trp-Trp binding cooperativity.

Early Events in G-quadruplex Folding Captured by Time-Resolved Small-Angle X-Ray Scattering

Time-resolved small-angle X-ray experiments (TR-SAXS) are reported here that capture and quantify a previously unknown rapid collapse of the unfolded oligonucleotide as an early step in G4 folding of hybrid 1 and hybrid 2 telomeric G-quadruplex structures. The rapid collapse, initiated by a pH jump, is characterized by an exponential decrease in the radius of gyration from 20.6 to 12.6 Å. The collapse is monophasic and is complete in less than 600 ms. Additional hand-mixing pH-jump kinetic studies show that slower kinetic steps follow the collapse. The folded and unfolded states at equilibrium were further characterized by SAXS studies and other biophysical tools, to show that G4 unfolding was complete at alkaline pH, but not in LiCl solution as is often claimed. The SAXS Ensemble Optimization Method (EOM) analysis reveals models of the unfolded state as a dynamic ensemble of flexible oligonucleotide chains with a variety of transient hairpin structures. These results suggest a G4 folding pathway in which a rapid collapse, analogous to molten globule formation seen in proteins, is followed by a confined conformational search within the collapsed particle to form the native contacts ultimately found in the stable folded form.

Crystal structure of archaeal IF5A-DHS complex reveals insights into the hypusination mechanism

The translation factor IF5A is highly conserved in Eukarya and Archaea and undergoes a unique post-translational hypusine modification by the deoxyhypusine synthase (DHS) enzyme. DHS transfers the butylamine moiety from spermidine to IF5A using NAD as a cofactor, forming a deoxyhypusine intermediate. IF5A is a key player in protein synthesis, preventing ribosome stalling in proline-rich sequences during translation elongation and facilitating translation elongation and termination. Additionally, human eIF5A participates in various essential cellular processes and contributes to cancer metastasis, with inhibiting hypusination showing anti-proliferative effects. The hypusination pathway of IF5A is therefore an attractive new therapeutic target. We elucidated the 2.0 Å X-ray crystal structure of the archaeal DHS-IF5A complex, revealing hetero-octameric architecture and providing a detailed view of the complex active site including the hypusination loop. This structure, along with biophysical data and molecular dynamics simulations, provides new insights into the catalytic mechanism of the hypusination reaction.

Structural Adaptation of Secondary p53 Binding Sites on MDM2 and MDMX

The thermodynamics of secondary p53 binding sites on MDM2 and MDMX were evaluated using p53 peptides containing residues 16-29, 17-35, and 1-73. All the peptides had large, negative heat capacity (ΔCp), consistent with the burial of p53 residues F19, W23, and L26 in the primary binding sites of MDM2 and MDMX. MDMX has a higher affinity and more negative ΔCp than MDM2 for p5317-35, which is due to MDMX stabilization and not additional interactions with the secondary binding site. ΔCp measurements show binding to the secondary site is inhibited by the disordered tails of MDM2 for WT p53 but not a more helical mutant where proline 27 is changed to alanine. This result is supported by all-atom molecular dynamics simulations showing that p53 residues 30-35 turn away from the disordered tails of MDM2 in P27A17-35 and make direct contact with this region in p5317-35. Molecular dynamics simulations also suggest that an intramolecular methionine-aromatic motif found in both MDM2 and MDMX structurally adapts to support multiple p53 binding modes with the secondary site. ΔCp measurements also show that tighter binding of the P27A mutant to MDM2 and MDMX is due to increased helicity, which reduces the energetic penalty associated with coupled folding and binding. Our results will facilitate the design of selective p53 inhibitors for MDM2 and MDMX.

Uncovering the Association Mechanism between Two Intrinsically Flexible Proteins

The understanding of protein-protein interaction mechanisms is key to the atomistic description of cell signaling pathways and for the development of new drugs. In this context, the mechanism of intrinsically disordered proteins folding upon binding has attracted attention. The VirB9 C-terminal domain (VirB9Ct) and the VirB7 N-terminal motif (VirB7Nt) associate with VirB10 to form the outer membrane core complex of the Type IV Secretion System injectisome. Despite forming a stable and rigid complex, VirB7Nt behaves as a random coil, while VirB9Ct is intrinsically dynamic in the free state. Here we combined NMR, stopped-flow fluorescence, and computer simulations using structure-based models to characterize the VirB9Ct-VirB7Nt coupled folding and binding mechanism. Qualitative data analysis suggested that VirB9Ct preferentially binds to VirB7Nt by way of a conformational selection mechanism at lower temperatures. However, at higher temperatures, energy barriers between different VirB9Ct conformations are more easily surpassed. Under these conditions the formation of non-native initial encounter complexes may provide alternative pathways toward the native complex conformation. These observations highlight the intimate relationship between folding and binding, calling attention to the fact that the two molecular partners must search for the most favored intramolecular and intermolecular interactions on a rugged and funnelled conformational energy landscape, along which multiple intermediates may lead to the final native state.

The molecular basis for cellular function of intrinsically disordered protein regions

Intrinsically disordered protein regions exist in a collection of dynamic interconverting conformations that lack a stable 3D structure. These regions are structurally heterogeneous, ubiquitous and found across all kingdoms of life. Despite the absence of a defined 3D structure, disordered regions are essential for cellular processes ranging from transcriptional control and cell signalling to subcellular organization. Through their conformational malleability and adaptability, disordered regions extend the repertoire of macromolecular interactions and are readily tunable by their structural and chemical context, making them ideal responders to regulatory cues. Recent work has led to major advances in understanding the link between protein sequence and conformational behaviour in disordered regions, yet the link between sequence and molecular function is less well defined. Here we consider the biochemical and biophysical foundations that underlie how and why disordered regions can engage in productive cellular functions, provide examples of emerging concepts and discuss how protein disorder contributes to intracellular information processing and regulation of cellular function.

Dissection of integrated readout reveals the structural thermodynamics of DNA selection by transcription factors

Nucleobases such as inosine have been extensively utilized to map direct contacts by proteins in the DNA groove. Their deployment as targeted probes of dynamics and hydration, which are dominant thermodynamic drivers of affinity and specificity, has been limited by a paucity of suitable experimental models. We report a joint crystallographic, thermodynamic, and computational study of the bidentate complex of the arginine side chain with a Watson-Crick guanine (Arg×GC), a highly specific configuration adopted by major transcription factors throughout the eukaryotic branches in the Tree of Life. Using the ETS-family factor PU.1 as a high-resolution structural framework, inosine substitution for guanine resulted in a sharp dissection of conformational dynamics and hydration and elucidated their role in the DNA specificity of PU.1. Our work suggests an under-exploited utility of modified nucleobases in untangling the structural thermodynamics of interactions, such as the Arg×GC motif, where direct and indirect readout are tightly integrated.

Conformational entropy in molecular recognition of intrinsically disordered proteins

Broad conformational ensembles make intrinsically disordered proteins or regions entropically intriguing. Although methodologically challenging and understudied, emerging studies into their changes in conformational entropy (ΔS°conf) upon complex formation have provided both quantitative and qualitative insight. Recent work based on thermodynamics from isothermal titration calorimetry and NMR spectroscopy uncovers an expanded repertoire of regulatory mechanisms, where ΔS°conf plays roles in partner selection, state behavior, functional buffering, allosteric regulation, and drug design. We highlight these mechanisms to display the large entropic reservoir of IDPs for the regulation of molecular communication. We call upon the field to make efforts to contribute to this insight as more studies are needed for forwarding mechanistic decoding of intrinsically disordered proteins and their complexes.

Functional unfoldomics: Roles of intrinsic disorder in protein (multi)functionality

Intrinsically disordered proteins (IDPs), which are functional proteins without stable tertiary structure, and hybrid proteins containing ordered domains and intrinsically disordered regions (IDRs) constitute prominent parts of all proteomes collectively known as unfoldomes. IDPs/IDRs exist as highly dynamic structural ensembles of rapidly interconverting conformations and are characterized by the exceptional structural heterogeneity, where their different parts are (dis)ordered to different degree, and their overall structure represents a complex mosaic of foldons, inducible foldons, inducible morphing foldons, non-foldons, semifoldons, and even unfoldons. Despite their lack of unique 3D structures, IDPs/IDRs play crucial roles in the control of various biological processes and the regulation of different cellular pathways and are commonly involved in recognition and signaling, indicating that the disorder-based functional repertoire is complementary to the functions of ordered proteins. Furthermore, IDPs/IDRs are frequently multifunctional, and this multifunctionality is defined by their structural flexibility and heterogeneity. Intrinsic disorder phenomenon is at the roots of the structure-function continuum model, where the structure continuum is defined by the presence of differently (dis)ordered regions, and the function continuum arises from the ability of all these differently (dis)ordered parts to have different functions. In their everyday life, IDPs/IDRs utilize a broad spectrum of interaction mechanisms thereby acting as interaction specialists. They are crucial for the biogenesis of numerous proteinaceous membrane-less organelles driven by the liquid-liquid phase separation. This review introduces functional unfoldomics by representing some aspects of the intrinsic disorder-based functionality.

Quantifying Amide-Aromatic Interactions at Molecular and Atomic Levels: Experimentally Determined Enthalpic and Entropic Contributions to Interactions of Amide sp2O, N, C and sp3C Unified Atoms with Naphthalene sp2C Atoms in Water

In addition to amide hydrogen bonds and the hydrophobic effect, interactions involving π-bonded sp2 atoms of amides, aromatics, and other groups occur in protein self-assembly processes including folding, oligomerization, and condensate formation. These interactions also occur in aqueous solutions of amide and aromatic compounds, where they can be quantified. Previous analysis of thermodynamic coefficients quantifying net-favorable interactions of amide compounds with other amides and aromatics revealed that interactions of amide sp2O with amide sp2N unified atoms (presumably C═O···H-N hydrogen bonds) and amide/aromatic sp2C (lone pair π, n-π*) are particularly favorable. Sp3C-sp3C (hydrophobic), sp3C-sp2C (hydrophobic, CH-π), sp2C-sp2C (hydrophobic, π-π), and sp3C-sp2N interactions are favorable, sp2C-sp2N interactions are neutral, while sp2O-sp2O and sp2N-sp2N self-interactions and sp2O-sp3C interactions are unfavorable. Here, from determinations of favorable effects of 14 amides on naphthalene solubility at 10, 25, and 45 °C, we dissect amide-aromatic interaction free energies into enthalpic and entropic contributions and find these vary systematically with amide composition. Analysis of these results yields enthalpic and entropic contributions to intrinsic strengths of interactions of amide sp2O, sp2N, sp2C, and sp3C unified atoms with aromatic sp2C atoms. For each interaction, enthalpic and entropic contributions have the same sign and are much larger in magnitude than the interaction free energy itself. The amide sp2O-aromatic sp2C interaction is enthalpy-driven and entropically unfavorable, consistent with direct chemical interaction (e.g., lone pair-π), while amide sp3C- and sp2C-aromatic sp2C interactions are entropy-driven and enthalpically unfavorable, consistent with hydrophobic effects. These findings are relevant for interactions involving π-bonded sp2 atoms in protein processes.

Asymmetric dimerization in a transcription factor superfamily is promoted by allosteric interactions with DNA

Transcription factors, such as nuclear receptors achieve precise transcriptional regulation by means of a tight and reciprocal communication with DNA, where cooperativity gained by receptor dimerization is added to binding site sequence specificity to expand the range of DNA target gene sequences. To unravel the evolutionary steps in the emergence of DNA selection by steroid receptors (SRs) from monomeric to dimeric palindromic binding sites, we carried out crystallographic, biophysical and phylogenetic studies, focusing on the estrogen-related receptors (ERRs, NR3B) that represent closest relatives of SRs. Our results, showing the structure of the ERR DNA-binding domain bound to a palindromic response element (RE), unveil the molecular mechanisms of ERR dimerization which are imprinted in the protein itself with DNA acting as an allosteric driver by allowing the formation of a novel extended asymmetric dimerization region (KR-box). Phylogenetic analyses suggest that this dimerization asymmetry is an ancestral feature necessary for establishing a strong overall dimerization interface, which was progressively modified in other SRs in the course of evolution.

Quantifying Amide-Aromatic Interactions at Molecular and Atomic Levels: Experimentally-determined Enthalpic and Entropic Contributions to Interactions of Amide sp 2 O, N, C and sp 3 C Unified Atoms with Naphthalene sp 2 C Atoms in Water

In addition to amide hydrogen bonds and the hydrophobic effect, interactions involving π-bonded sp 2 atoms of amides, aromatics and other groups occur in protein self-assembly processes including folding, oligomerization and condensate formation. These interactions also occur in aqueous solutions of amide and aromatic compounds, where they can be quantified. Previous analysis of thermodynamic coefficients quantifying net-favorable interactions of amide compounds with other amides and aromatics revealed that interactions of amide sp 2 O with amide sp 2 N unified atoms (presumably C=O···H-N hydrogen bonds) and amide/aromatic sp 2 C (lone pair-π, n-π * ) are particularly favorable. Sp 3 C-sp 3 C (hydrophobic), sp 3 C-sp 2 C (hydrophobic, CH-π), sp 2 C-sp 2 C (hydrophobic, π-π) and sp 3 C-sp 2 N interactions are favorable, sp 2 C-sp 2 N interactions are neutral, while sp 2 O-sp 2 O and sp 2 N-sp 2 N self-interactions and sp 2 O-sp 3 C interactions are unfavorable. Here, from determinations of favorable effects of fourteen amides on naphthalene solubility at 10, 25 and 45 °C, we dissect amide-aromatic interaction free energies into enthalpic and entropic contributions and find these vary systematically with amide composition. Analysis of these results yields enthalpic and entropic contributions to intrinsic strengths of interactions of amide sp 2 O, sp 2 N, sp 2 C and sp 3 C unified atoms with aromatic sp 2 C atoms. For each interaction, enthalpic and entropic contributions have the same sign and are much larger in magnitude than the interaction free energy itself. The amide sp 2 O-aromatic sp 2 C interaction is enthalpy-driven and entropically unfavorable, consistent with direct chemical interaction (e.g. lone pair-π) while amide sp 3 C- and sp 2 C-aromatic sp 2 C interactions are entropy-driven and enthalpically unfavorable, consistent with hydrophobic effects. These findings are relevant for interactions involving π-bonded sp 2 atoms in protein processes.

Table of Contents Graphic

Analysis of Aptamer-Small Molecule Binding Interactions Using Isothermal Titration Calorimetry

Isothermal titration calorimetry (ITC) is a technique where the heat given off, or absorbed, during a binding event is measured and used to determine the binding thermodynamics and affinity associated with binding. This protocol focuses on ITC applications for studying aptamer interactions with small molecule ligands where ITC has the advantage of being a label-free solution-based technique. The limitation of ITC using a relatively large amount of material compared to other analytical techniques is not applicable here as large amounts of nucleic acids, especially DNA, can be readily obtained. In this chapter we describe how to use ITC methods to measure the thermodynamics and affinity of binding using the interaction of quinine with a DNA cocaine-binding aptamer as an example.

The RGG motif proteins: Interactions, functions, and regulations

Proteins with motifs rich in arginines and glycines were discovered decades ago and are functionally involved in a staggering range of essential processes in the cell. Versatile, specific, yet adaptable molecular interactions enabled by the unique combination of arginine and glycine, combined with multiplicity of molecular recognition conferred by repeated di-, tri-, and multiple peptide motifs, allow RGG motif proteins to interact with a broad range of proteins and nucleic acids. Furthermore, posttranslational modifications at the arginines in the motif extend the RGG protein's capacity for a fine-tuned regulation. In this review, we focus on the biochemical properties of the RGG motif, its molecular interactions with RNAs and proteins, and roles of the posttranslational modification in modulating their interactions. We discuss current knowledge of the RGG motif proteins involved in mRNA transport and translation, highlight our merging understanding of their molecular functions in translational regulation and summarize areas of research in the future critical in understanding this important family of proteins. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications Translation > Mechanisms.

Silica adsorption tag derived from the silica polycondensation protein glassin for the immobilization of soluble proteins

Glassin is a water-soluble protein from the siliceous skeleton of a marine sponge that adsorbs tightly to silica at pH 7.0-9.0 and accelerates silica particle formation from silicic acid. Glassin comprises three distinct domains: a His and Asp-rich (HD) domain, a Pro-rich (P) domain, and a His and Thr-rich (HT) domain. Here, we investigated the roles of these three domains in silica adsorption by using glutathione S-transferase (GST) fused with glassin or with each domain. GST fused with the HD domain exhibited tight adsorption, equivalent to that of GST fused with the full-length glassin sequence at values above 7.0. The apparent K_d values for the binding of full-length glassin and HD to fumed silica at pH 7.0 were 20.8 and 22.7 nM, respectively, indicating that this domain greatly contributes to the silica adsorption ability of glassin. In addition, no internal cleavage was observed in the HD domain, whereas GST fused with the full-length glassin sequence exhibited internal cleavage. The HD domain adsorbed on silica did not dissociate even at pH 6.0. Given these findings, we concluded that the HD domain has potential as a tag for the immobilization of soluble proteins.

Computer simulation of molecular recognition in biomolecular system: from in silico screening to generalized ensembles

Prediction of ligand-receptor complex structure is important in both the basic science and the industry such as drug discovery. We report various computation molecular docking methods: fundamental in silico (virtual) screening, ensemble docking, enhanced sampling (generalized ensemble) methods, and other methods to improve the accuracy of the complex structure. We explain not only the merits of these methods but also their limits of application and discuss some interaction terms which are not considered in the in silico methods. In silico screening and ensemble docking are useful when one focuses on obtaining the native complex structure (the most thermodynamically stable complex). Generalized ensemble method provides a free-energy landscape, which shows the distribution of the most stable complex structure and semi-stable ones in a conformational space. Also, barriers separating those stable structures are identified. A researcher should select one of the methods according to the research aim and depending on complexity of the molecular system to be studied.

Intrinsic protein disorder uncouples affinity from binding specificity

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins often function by molecular recognition, in which they undergo induced folding. Based on prior generalizations, the idea prevails in the IDP field that due to the entropic penalty of induced folding, the major functional advantage associated with this binding mode is "uncoupling" specificity from binding strength. Nevertheless, both weaker binding and high specificity of IDPs/IDRs rest on limited experimental observations, making these assumptions more speculations than evidence-supported facts. The issue is also complicated by the rather vague concept of specificity that lacks an exact measure, such as the Kd for binding strength. We addressed these issues by creating and analyzing a comprehensive dataset of well-characterized ID/globular protein complexes, for which both the atomic structure of the complex and free energy (ΔG, Kd ) of interaction is known. Through this analysis, we provide evidence that the affinity distributions of IDP/globular and globular/globular complexes show different trends, whereas specificity does not connote to weaker binding strength of IDPs/IDRs. Furthermore, protein disorder extends the spectrum in the direction of very weak interactions, which may have important regulatory consequences and suggest that, in a biological sense, strict correlation of specificity and binding strength are uncoupled by structural disorder.

Binding-Induced Diversity of a Human Telomeric G-Quadruplex Stability Phase Space

The structural polymorphism of G-quadruplex nucleic acids is an important factor in their recognition by proteins and small-molecule ligands. However, it is not clear why the binding of several ligands alters G-quadruplex topology. We addressed this question by following the (un)folding and binding of the human telomeric fragment 5′-(GGGTTA)₃GGGT-3′ (22GT) by calorimetry (DSC, ITC) and spectroscopy (CD). A thermodynamic analysis of the obtained data led to a detailed description of the topological phase space of stability (phase diagram) of 22GT and shows how it changes in the presence of a specific bisquinolinium ligand (360A). Various 1:1 and 2:1 ligand–quadruplex complexes were observed. With increasing temperature, the 1:1 complexes transformed into 2:1 complexes, which is attributed to the preferential binding of the ligand to the folding intermediates. Overall, the dissection of the thermodynamic parameters in combination with molecular modelling clarified the driving forces of the topological quadruplex transformations in a wide range of ligand concentrations and temperatures.

Fly casting with ligand sliding and orientational selection supporting complex formation of a GPCR and a middle sized flexible molecule

A GA-guided multidimensional virtual-system coupled molecular dynamics (GA-mD-VcMD) simulation was conducted to elucidate binding mechanisms of a middle-sized flexible molecule, bosentan, to a GPCR protein, human endothelin receptor type B (hETB). GA-mD-VcMD is a generalized ensemble method that produces a free-energy landscape of the ligand-receptor binding by searching large-scale motions accompanied with stable maintenance of the fragile cell-membrane structure. All molecular components (bosentan, hETB, membrane, and solvent) were represented with an all-atom model. Then sampling was conducted from conformations where bosentan was distant from the binding site in the hETB binding pocket. The deepest basin in the resultant free-energy landscape was assigned to native-like complex conformation. The following binding mechanism was inferred. First, bosentan fluctuating randomly in solution is captured using a tip region of the flexible N-terminal tail of hETB via nonspecific attractive interactions (fly casting). Bosentan then slides occasionally from the tip to the root of the N-terminal tail (ligand–sliding). During this sliding, bosentan passes the gate of the binding pocket from outside to inside of the pocket with an accompanying rapid reduction of the molecular orientational variety of bosentan (orientational selection). Last, in the pocket, ligand–receptor attractive native contacts are formed. Eventually, the native-like complex is completed. The bosentan-captured conformations by the tip-region and root-region of the N-terminal tail correspond to two basins in the free-energy landscape. The ligand-sliding corresponds to overcoming of a free-energy barrier between the basins.

Structural and dynamic studies of DNA recognition by NF-κB p50 RHR homodimer using methyl NMR spectroscopy

Protein dynamics involving higher-energy sparsely populated conformational substates are frequently critical for protein function. This study describes the dynamics of the homodimer (p50)2 of the p50 Rel homology region (RHR) of the transcription factor NF-κB, using 13C relaxation dispersion experiments with specifically (13C, 1H)-labeled methyl groups of Ile (δ), Leu and Val. Free (p50)2 is highly dynamic in solution, showing μs-ms relaxation dispersion consistent with exchange between the ground state and higher energy substates. These fluctuations propagate from the DNA-binding loops through the core of the domain. The motions are damped in the presence of κB DNA, but the NMR spectra of the DNA complexes reveal multiple local conformations of the p50 RHR homodimer bound to certain κB DNA sequences. Varying the length and sequence of κB DNA revealed two factors that promote a single bound conformation for the complex: the length of the κB site in the duplex and a symmetrical sequence of guanine nucleotides at both ends of the recognition motif. The dynamic nature of the DNA-binding loops, together with the multiple bound conformations of p50 RHR with certain κB sites, is consistent with variations in the transcriptional activity of the p50 homodimer with different κB sequences.

The free energy folding penalty accompanying binding of intrinsically disordered α-helical motifs

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes and preform critical roles in many cellular processes, most often through the association with globular proteins. Despite lacking a stable three-dimensional structure by themselves, they may acquire a defined conformation upon binding globular targets. The most common type of secondary structure acquired by these binding motifs entails formation of an α-helix. It has been hypothesized that such disorder-to-order transitions are associated with a significant free energy penalty due to IDP folding, which reduces the overall IDP-target affinity. However, the exact magnitude of IDP folding penalty in α-helical binding motifs has not been systematically estimated. Here, we report the folding penalty contributions for 30 IDPs undergoing folding-upon-binding and find that the average IDP folding penalty is +2.0 kcal/mol and ranges from 0.7 to 3.5 kcal/mol. We observe that the folding penalty scales approximately linearly with the change in IDP helicity upon binding, which provides a simple empirical way to estimate folding penalty. We analyze to what extent do pre-structuring and target-bound IDP dynamics (fuzziness) reduce the folding penalty and find that these effects combined, on average, reduce the folding cost by around half. Taken together, the presented analysis provides a quantitative basis for understanding the role of folding penalty in IDP-target interactions and introduces a method estimate this quantity. Estimation and reduction of IDP folding penalty may prove useful in the rational design of helix-stabilized inhibitors of IDP-target interactions. STATEMENT: The α-helical binding motifs are ubiquitous among the intrinsically disordered proteins (IDPs). Upon binding their targets, they undergo a disorder-to-order transition, which is accompanied by a significant folding penalty whose magnitude is generally not known. Here, we use recently developed statistical-thermodynamic model to estimate the folding penalties for 30 IDPs and clarify the roles of IDP pre-folding and bound-state dynamics in reducing the folding penalty.

An Outlook on the Complexity of Protein Morphogenesis in Health and Disease

The study of the mechanisms whereby proteins achieve their native functionally competent conformation has been a key issue in molecular biosciences over the last 6 decades. Nevertheless, there are several debated issues and open problems concerning some aspects of this fundamental problem. By considering the emerging complexity of the so-called “native state,” we attempt hereby to propose a personal account on some of the key topics in the field, ranging from the relationships between misfolding and diseases to the significance of protein disorder. Finally, we briefly describe the recent and exciting advances in predicting protein structures from their amino acid sequence.

Initial heat analysis in dissociation isothermal titration calorimetry: An analytical tool for thermodynamic dissection of biomolecular condensates

Multi-domain proteins or intrinsically disordered proteins (IDPs) often undergo liquid-liquid phase separation (LLPS) and form membraneless organelles or protein condensates. Such compartmentalization is considered critical in many cellular processes dynamically modulated by various external signals. However, molecular mechanisms underlying potential regulatory functions of the protein condensates remain obscure due to a limited understanding of the driving forces for their assembly. Here we propose isothermal titration calorimetry (ITC) as an efficient analytical tool to dissociate condensates and measure the corresponding dissociation heat. Subsequent analysis of the initial dissociation heat as a function of total protein concentration allows simple and accurate determination of the thermodynamic parameters for cooperative condensate formations including the dissociation (or condensation) enthalpy and the critical protein concentration. By performing systematic simulations, we further demonstrate that the initial heat analysis is sufficiently robust to quantitatively dissect protein condensates with a broad range of thermodynamic properties. Therefore, our proposed method analyzing the initial heat measured in dissociation ITC provides opportunities to further scrutinize the thermodynamic quantities as functions of solution variables to explore the molecular driving forces of LLPS.

Geometric algebra generation of molecular surfaces

Geometric algebra is a powerful framework that unifies mathematics and physics. Since its revival in the 1960s, it has attracted great attention and has been exploited in fields like physics, computer science and engineering. This work introduces a geometric algebra method for the molecular surface generation that uses the Clifford-Fourier transform (CFT) which is a generalization of the classical Fourier transform. Notably, the classical Fourier transform and CFT differ in the derivative property in [Formula: see text] for k even. This distinction is due to the non-commutativity of geometric product of pseudoscalars with multivectors and has significant consequences in applications. We use the CFT in [Formula: see text] to benefit from the derivative property in solving partial differential equations (PDEs). The CFT is used to solve the mode decomposition process in PDE transform. Two different initial cases are proposed to make the initial shapes in the present method. The proposed method is applied first to small molecules and proteins. To validate the method, the molecular surfaces generated are compared to surfaces of other definitions. Applications are considered to protein electrostatic surface potentials and solvation free energy. This work opens the door for further applications of geometric algebra and CFT in biological sciences.

Following the tracks: how transcription factor binding dynamics control transcription

Transcription, the process of copying genetic information from DNA to messenger RNA, is regulated by sequence-specific DNA binding proteins known as transcription factors (TFs). Recent advances in single-molecule tracking (SMT) technologies have enabled visualization of individual TF molecules as they diffuse and interact with the DNA in the context of living cells. These SMT studies have uncovered multiple populations of DNA binding events characterized by their distinctive DNA residence times. In this perspective, we review recent insights into how these residence times relate to specific and non-specific DNA binding, as well as the contribution of TF domains on the DNA binding dynamics. We discuss different models that aim to link transient DNA binding by TFs to bursts of transcription and present an outlook for how future advances in microscopy development may broaden our understanding of the dynamics of the molecular steps that underlie transcription activation.

Metal Ion Binding Induces Local Protein Unfolding and Destabilizes Human Carbonic Anhydrase II

Human carbonic anhydrase II (HCA) is a robust metalloprotein and an excellent biological model system to study the thermodynamics of metal ion coordination. Apo-HCA binds one zinc ion or two copper ions. We studied these binding processes at five temperatures (15-35 °C) using isothermal titration calorimetry, yielding thermodynamic parameters corrected for pH and buffer effects. We then sought to identify binding-induced structural changes. Our data suggest that binding at the active site organizes 6-8 residues; however, copper binding near the N-terminus results in a net unfolding of 6-7 residues. This surprising destabilization was confirmed using circular dichroism and protein stability measurements. Metal binding induced unfolding may represent an important regulatory mechanism, but it may be easily missed by NMR and X-ray crystallography. Thus, in addition to highlighting a highly novel binding-induced unfolding event, we demonstrate the value of calorimetry for studying the structural implications of metal binding.

The Non-continuum Nature of Eukaryotic Transcriptional Regulation

Eukaryotic transcription factors are versatile mediators of specificity in gene regulation. This versatility is achieved through mutual specification by context-specific DNA binding on the one hand, and identity-specific protein-protein partnerships on the other. This interactivity, known as combinatorial control, enables a repertoire of complex transcriptional outputs that are qualitatively disjoint, or non-continuum, with respect to binding affinity. This feature contrasts starkly with prokaryotic gene regulators, whose activities in general vary quantitatively in step with binding affinity. Biophysical studies on prokaryotic model systems and more recent investigations on transcription factors highlight an important role for folded state dynamics and molecular hydration in protein/DNA recognition. Analysis of molecular models of combinatorial control and recent literature in low-affinity gene regulation suggest that transcription factors harbor unique conformational dynamics that are inaccessible or unused by prokaryotic DNA-binding proteins. Thus, understanding the intrinsic dynamics involved in DNA binding and co-regulator recruitment appears to be a key to understanding how transcription factors mediate non-continuum outcomes in eukaryotic gene expression, and how such capability might have evolved from ancient, structurally conserved counterparts.

The Honey Bee Gene Bee Antiviral Protein-1 Is a Taxonomically Restricted Antiviral Immune Gene

Insects have evolved a wide range of strategies to combat invading pathogens, including viruses. Genes that encode proteins involved in immune responses often evolve under positive selection due to their co-evolution with pathogens. Insect antiviral defense includes the RNA interference (RNAi) mechanism, which is triggered by recognition of non-self, virally produced, double-stranded RNAs. Indeed, insect RNAi genes (e.g., dicer and argonaute-2) are under high selective pressure. Honey bees (Apis mellifera) are eusocial insects that respond to viral infections via both sequence specific RNAi and a non-sequence specific dsRNA triggered pathway, which is less well-characterized. A transcriptome-level study of virus-infected and/or dsRNA-treated honey bees revealed increased expression of a novel antiviral gene, GenBank: MF116383, and in vivo experiments confirmed its antiviral function. Due to in silico annotation and sequence similarity, MF116383 was originally annotated as a probable cyclin-dependent serine/threonine-protein kinase. In this study, we confirmed that MF116383 limits virus infection, and carried out further bioinformatic and phylogenetic analyses to better characterize this important gene-which we renamed bee antiviral protein-1 (bap1). Phylogenetic analysis revealed that bap1 is taxonomically restricted to Hymenoptera and Blatella germanica (the German cockroach) and that the majority of bap1 amino acids are evolving under neutral selection. This is in-line with the results from structural prediction tools that indicate Bap1 is a highly disordered protein, which likely has relaxed structural constraints. Assessment of honey bee gene expression using a weighted gene correlation network analysis revealed that bap1 expression was highly correlated with several immune genes-most notably argonaute-2. The coexpression of bap1 and argonaute-2 was confirmed in an independent dataset that accounted for the effect of virus abundance. Together, these data demonstrate that bap1 is a taxonomically restricted, rapidly evolving antiviral immune gene. Future work will determine the role of bap1 in limiting replication of other viruses and examine the signal cascade responsible for regulating the expression of bap1 and other honey bee antiviral defense genes, including coexpressed ago-2, and determine whether the virus limiting function of bap1 acts in parallel or in tandem with RNAi.

Quantification of Conformational Entropy Unravels Effect of Disordered Flanking Region in Coupled Folding and Binding

Intrinsic disorder (ID) constitutes a new dimension to the protein structure-function relationship. The ability to undergo conformational changes upon binding is a key property of intrinsically disordered proteins and remains challenging to study using conventional methods. A 1994 paper by R. S. Spolar and M. T. Record presented a thermodynamic approach for estimating changes in conformational entropy based on heat capacity changes, allowing quantification of residues folding upon binding. Here, we adapt the method for studies of intrinsically disordered proteins. We integrate additional data to provide a broader experimental foundation for the underlying relations and, based on >500 protein-protein complexes involving disordered proteins, reassess a key relation between polar and nonpolar surface area changes, previously determined using globular protein folding. We demonstrate the improved suitability of the adapted method to studies of the folded αα-hub domain RST from radical-induced cell death 1, whose interactome is characterized by ID. From extensive thermodynamic data, quantifying the conformational entropy changes upon binding, and comparison to the NMR structure, the adapted method improves accuracy for ID-based studies. Furthermore, we apply the method, in conjunction with NMR, to reveal hitherto undetected effects of interaction-motif context. Thus, inclusion of the disordered context of the DREB2A RST-binding motif induces structuring of the binding motif, resulting in major enthalpy-entropy compensation in the interaction interface. This study, also evaluating additional interactions, demonstrates the strength of the ID-adapted Spolar-Record thermodynamic approach for dissection of structural features of ID-based interactions, easily overlooked in traditional studies, and for translation of these into mechanistic knowledge.

Thermodynamic Models for Assembly of Intrinsically Disordered Protein Hubs with Multiple Interaction Partners

Prevalent in diverse protein interactomes, intrinsically disordered proteins or regions (IDPs or IDRs) often drive assembly of higher-order macromolecular complexes, using multiple target-binding motifs. Such IDP hubs are suggested to process various cellular signals and coordinate relevant biological processes. However, the mechanism of assembly and functional regulation of IDP hubs remains elusive due to the challenges in dissecting their intricate protein-protein interaction networks. Here we present basic thermodynamic models for the assembly of simple IDP hubs with multiple target proteins, constructing partition functions from fundamental binding parameters. We combined these basic functions to develop advanced thermodynamic models to analyze the assembly of the Nup153 hubs interacting with multiple karyopherin β1 (Kap) molecules, critical components of nucleocytoplasmic transport. The thermodynamic analysis revealed a complex organization of the Kap binding sites within the C-terminal IDR of Nup153 including a high-affinity 1:1 interaction site and a series of low-affinity sites for binding of multiple Kaps with negative cooperativity. The negative cooperativity arises from the overlapping nature of the low-affinity sites where Kap occupies multiple dipeptide motifs. The quantitative dissection culminated in construction of the Nup153 hub ensemble, elucidating how distribution among various Kap-bound states is modulated by Kap concentration and competing nuclear proteins. In particular, the Kap occupancy of the IDR can be fine-tuned by varying the location of competition within the overlapping sites, suggesting coupling of specific nuclear processes to distinct transport activities. In general, our results demonstrate the feasibility and a potential mechanism for manifold regulation of IDP functions by diverse cellular signals.

Testing mechanisms of DNA sliding by architectural DNA-binding proteins: dynamics of single wild-type and mutant protein molecules in vitro and in vivo

Architectural DNA-binding proteins (ADBPs) are abundant constituents of eukaryotic or bacterial chromosomes that bind DNA promiscuously and function in diverse DNA reactions. They generate large conformational changes in DNA upon binding yet can slide along DNA when searching for functional binding sites. Here we investigate the mechanism by which ADBPs diffuse on DNA by single-molecule analyses of mutant proteins rationally chosen to distinguish between rotation-coupled diffusion and DNA surface sliding after transient unbinding from the groove(s). The properties of yeast Nhp6A mutant proteins, combined with molecular dynamics simulations, suggest Nhp6A switches between two binding modes: a static state, in which the HMGB domain is bound within the minor groove with the DNA highly bent, and a mobile state, where the protein is traveling along the DNA surface by means of its flexible N-terminal basic arm. The behaviors of Fis mutants, a bacterial nucleoid-associated helix-turn-helix dimer, are best explained by mobile proteins unbinding from the major groove and diffusing along the DNA surface. Nhp6A, Fis, and bacterial HU are all near exclusively associated with the chromosome, as packaged within the bacterial nucleoid, and can be modeled by three diffusion modes where HU exhibits the fastest and Fis the slowest diffusion.

What's in the BAGs? Intrinsic disorder angle of the multifunctionality of the members of a family of chaperone regulators

In humans, the family of Bcl-2 associated athanogene (BAG) proteins includes six members characterized by exceptional multifunctionality and engagement in the pathogenesis of various diseases. All of them are capable of interacting with a multitude of often unrelated binding partners. Such binding promiscuity and related functional and pathological multifacetedness cannot be explained or understood within the frames of the classical "one protein-one structure-one function" model, which also fails to explain the presence of multiple isoforms generated for BAG proteins by alternative splicing or alternative translation initiation and their extensive posttranslational modifications. However, all these mysteries can be solved by taking into account the intrinsic disorder phenomenon. In fact, high binding promiscuity and potential to participate in a broad spectrum of interactions with multiple binding partners, as well as a capability to be multifunctional and multipathogenic, are some of the characteristic features of intrinsically disordered proteins and intrinsically disordered protein regions. Such functional proteins or protein regions lacking unique tertiary structures constitute a cornerstone of the protein structure-function continuum concept. The aim of this paper is to provide an overview of the functional roles of human BAG proteins from the perspective of protein intrinsic disorder which will provide a means for understanding their binding promiscuity, multifunctionality, and relation to the pathogenesis of various diseases.

Substitutions at Nonconserved Rheostat Positions Modulate Function by Rewiring Long-Range, Dynamic Interactions

Amino acid substitutions at nonconserved protein positions can have noncanonical and "long-distance" outcomes on protein function. Such outcomes might arise from changes in the internal protein communication network, which is often accompanied by changes in structural flexibility. To test this, we calculated flexibilities and dynamic coupling for positions in the linker region of the lactose repressor protein. This region contains nonconserved positions for which substitutions alter DNA-binding affinity. We first chose to study 11 substitutions at position 52. In computations, substitutions showed long-range effects on flexibilities of DNA-binding positions, and the degree of flexibility change correlated with experimentally measured changes in DNA binding. Substitutions also altered dynamic coupling to DNA-binding positions in a manner that captured other experimentally determined functional changes. Next, we broadened calculations to consider the dynamic coupling between 17 linker positions and the DNA-binding domain. Experimentally, these linker positions exhibited a wide range of substitution outcomes: Four conserved positions tolerated hardly any substitutions ("toggle"), ten nonconserved positions showed progressive changes from a range of substitutions ("rheostat"), and three nonconserved positions tolerated almost all substitutions ("neutral"). In computations with wild-type lactose repressor protein, the dynamic couplings between the DNA-binding domain and these linker positions showed varied degrees of asymmetry that correlated with the observed toggle/rheostat/neutral substitution outcomes. Thus, we propose that long-range and noncanonical substitutions outcomes at nonconserved positions arise from rewiring long-range communication among functionally important positions. Such calculations might enable predictions for substitution outcomes at a range of nonconserved positions.

The importance of hydrophobic interactions in the structure of transcription systems

Hydrophobic forces play a crucial role in both the stability of B DNA and its interactions with proteins. In the present study, we postulate that the hydrophobic effect is an essential component in establishing specificity in the interaction transcription factor proteins with their consensus DNA sequence partners. The PDB coordinates of more than 50 transcription systems have been used to analyze the hydrophobic attraction of proteins towards their DNA consensus. This analysis includes computing the hydrophobic energy of the interacting molecules by means of their hydrophobic moments. Hydrophobic moments have successfully been used in previous studies involving self-assembly protein systems. In the present case, in spite of some variability, we found specificity in transcription factors when interacting with their respective consensus DNA sequences. By applying our model of biological membrane pattern for hydrophobic interactions, we postulate that hydrophobic forces constitute the necessary intermediate interaction between the unspecific electrostatic attraction for DNA phosphate groups and the very short-range interaction promoting hydrogen bonds. We conclude that hydrophobic interactions serve as the intermediate force guiding transcriptions factors towards the proper hydrogen bonds to their DNAs.

Halogen Atoms in the Protein-Ligand System. Structural and Thermodynamic Studies of the Binding of Bromobenzotriazoles by the Catalytic Subunit of Human Protein Kinase CK2

Binding of a family of brominated benzotriazoles to the catalytic subunit of human protein kinase CK2 (hCK2α) was used as a model system to assess the contribution of halogen bonding to protein–ligand interaction. CK2 is a constitutively active pleiotropic serine/threonine protein kinase that belongs to the CMGC group of eukaryotic protein kinases (EPKs). Due to the addiction of some cancer cells, CK2 is an attractive and well-characterized drug target. Halogenated benzotriazoles act as ATP-competitive inhibitors with unexpectedly good selectivity for CK2 over other EPKs. We have characterized the interaction of bromobenzotriazoles with hCK2α by X-ray crystallography, low-volume differential scanning fluorimetry, and isothermal titration calorimetry. Properties of free ligands in solution were additionally characterized by volumetric and RT-HPLC measurements. Thermodynamic data indicate that the affinity increases with bromo substitution, with greater contributions from 5- and 6-substituents than 4- and 7-substituents. Except for 4,7-disubstituted compounds, the bromobenzotriazoles adopt a canonical pose with the triazole close to lysine 68, which precludes halogen bonding. More highly substituted benzotriazoles adopt many additional noncanonical poses, presumably driven by a large hydrophobic contribution to binding. Some noncanonical ligand orientations allow the formation of halogen bonds with the hinge region. Consistent with a predominantly hydrophobic interaction, the isobaric heat capacity decreases upon ligand binding, the more so the higher the substitution.

Probing coupled conformational transitions of intrinsically disordered proteins in their interactions with target proteins

Intrinsically disordered proteins or regions (IDPs or IDRs) are abundant in the eukaryotic proteome and critical in regulation of dynamic cellular processes. Intensive structural investigations have proposed the molecular mechanisms of the interaction between IDRs and their binding partners. Here we extract the distinct thermodynamic features of coupled conformational transitions of IDRs founding the interaction mechanisms. We also present simulation tools to facilitate a design of the calorimetric experiments probing and quantifying the conformational transitions of IDRs. The suggested thermodynamic approach will further advance our understanding of distribution among multiple states of IDRs in their interactions with target molecules.

Conformational Plasticity and DNA-Binding Specificity of the Eukaryotic Transcription Factor Pax5

The eukaryotic transcription factor Pax5 has a DNA-binding Paired domain composed of two independent helical bundle subdomains joined by a flexible linker. Previously, we showed distinct biophysical properties of the N-terminal (NTD) and C-terminal (CTD) subdomains, with implications for how these two regions cooperate to distinguish nonspecific and cognate DNA sites [Perez-Borrajero, C., et al. (2016) J. Mol. Biol. 428, 2372-2391]. In this study, we combined experimental methods and molecular dynamics (MD) simulations to dissect the mechanisms underlying the functional differences between the Pax5 subdomains. Both subdomains showed a similar dependence of DNA-binding affinity on ionic strength. However, due to a greater contribution of non-ionic interactions, the NTD bound its cognate DNA half-site with an affinity approximately 10-fold higher than that of the CTD with its half-site. These interactions involve base-mediated contacts as evidenced by nuclear magnetic resonance spectroscopy-monitored chemical shift perturbations. Isothermal titration calorimetry revealed that favorable enthalpic and compensating unfavorable entropic changes were substantially larger for DNA binding by the NTD than by the CTD. Complementary MD simulations indicated that the DNA recognition helix H3 of the NTD is particularly flexible in the absence of DNA and undergoes the largest changes in conformational dynamics upon binding. Overall, these data suggest that the differences observed for the subdomains of Pax5 are due to the coupling of DNA binding with dampening of motions in the NTD required for specific base contacts. Thus, the conformational plasticity of the Pax5 Paired domain underpins the differing roles of its subdomains in association with nonspecific versus cognate DNA sites.

Mechanisms of Macromolecular Interactions Mediated by Protein Intrinsic Disorder

Intrinsically disordered proteins or regions (IDPs or IDRs) are widespread in the eukaryotic proteome. Although lacking stable three-dimensional structures in the free forms, IDRs perform critical functions in various cellular processes. Accordingly, mutations and altered expression of IDRs are associated with many pathological conditions. Hence, it is of great importance to understand at the molecular level how IDRs interact with their binding partners. In particular, discovering the unique interaction features of IDRs originating from their dynamic nature may reveal uncharted regulatory mechanisms of specific biological processes. Here we discuss the mechanisms of the macromolecular interactions mediated by IDRs and present the relevant cellular processes including transcription, cell cycle progression, signaling, and nucleocytoplasmic transport. Of special interest is the multivalent binding nature of IDRs driving assembly of multicomponent macromolecular complexes. Integrating the previous theoretical and experimental investigations, we suggest that such IDR-driven multiprotein complexes can function as versatile allosteric switches to process diverse cellular signals. Finally, we discuss the future challenges and potential medical applications of the IDR research.

Functions of short lifetime biological structures at large: the case of intrinsically disordered proteins

Although for more than a century a protein function was intimately associated with the presence of unique structure in a protein molecule, recent years witnessed a skyrocket rise of the appreciation of protein intrinsic disorder concept that emphasizes the importance of the biologically active proteins without ordered structures. In different proteins, the depth and breadth of disorder penetrance are different, generating an amusing spatiotemporal heterogeneity of intrinsically disordered proteins (IDPs) and intrinsically disordered protein region regions (IDPRs), which are typically described as highly dynamic ensembles of rapidly interconverting conformations (or a multitude of short lifetime structures). IDPs/IDPRs constitute a substantial part of protein kingdom and have unique functions complementary to functional repertoires of ordered proteins. They are recognized as interaction specialists and global controllers that play crucial roles in regulation of functions of their binding partners and in controlling large biological networks. IDPs/IDPRs are characterized by immense binding promiscuity and are able to use a broad spectrum of binding modes, often resulting in the formation of short lifetime complexes. In their turn, functions of IDPs and IDPRs are controlled by various means, such as numerous posttranslational modifications and alternative splicing. Some of the functions of IDPs/IDPRs are briefly considered in this review to shed some light on the biological roles of short-lived structures at large.

Experimentally determined strengths of favorable and unfavorable interactions of amide atoms involved in protein self-assembly in water

Folding and other protein self-assembly processes are driven by favorable interactions between O, N, and C unified atoms of the polypeptide backbone and side chains. These processes are perturbed by solutes that interact with these atoms differently than water does. Amide NH···O=C hydrogen bonding and various π-system interactions have been better characterized structurally or by simulations than experimentally in water, and unfavorable interactions are relatively uncharacterized. To address this situation, we previously quantified interactions of alkyl ureas with amide and aromatic compounds, relative to interactions with water. Analysis yielded strengths of interaction of each alkylurea with unit areas of different hybridization states of unified O, N, and C atoms of amide and aromatic compounds. Here, by osmometry, we quantify interactions of 10 pairs of amides selected to complete this dataset. An analysis yields intrinsic strengths of six favorable and four unfavorable atom-atom interactions, expressed per unit area of each atom and relative to interactions with water. The most favorable interactions are sp2O-sp2C (lone pair-π, presumably n-π*), sp2C-sp2C (π-π and/or hydrophobic), sp2O-sp2N (hydrogen bonding) and sp3C-sp2C (CH-π and/or hydrophobic). Interactions of sp3C with itself (hydrophobic) and with sp2N are modestly favorable, while sp2N interactions with sp2N and with amide/aromatic sp2C are modestly unfavorable. Amide sp2O-sp2O interactions and sp2O-sp3C interactions are more unfavorable, indicating the preference of amide sp2O to interact with water. These intrinsic interaction strengths are used to predict interactions of amides with proteins and chemical effects of amides (including urea, N-ethylpyrrolidone [NEP], and polyvinylpyrrolidone [PVP]) on protein stability.

Origin of heat capacity increment in DNA folding: The hydration effect

BACKGROUND

Understanding DNA folding thermodynamics is crucial for prediction of DNA thermal stability. It is now well established that DNA folding is accompanied by a decrease of the heat capacity ∆c_{p, F}, however its molecular origin is not understood. In analogy to protein folding it has been assumed that this is due to dehydration of DNA constituents, however no evidence exists to support this conclusion.

METHODS

Here we analyze partial molar heat capacity of nucleic bases and nucleosides in aqueous solutions obtained from calorimetric experiments and calculate the hydration heat capacity contribution ∆c_p^hyd.

RESULTS

We present hydration heat capacity contributions of DNA constituents and show that they correlate with the solvent accessible surface area. The average contribution for nucleic base dehydration is +0.56 J mol^-1 K^-1 Å^-2 and can be used to estimate the ∆c_{p, F} contribution for DNA folding.

CONCLUSIONS

We show that dehydration is one of the major sources contributing to the observed ∆c_{p, F} increment in DNA folding. Other possible sources contributing to the overall ∆c_{p, F} should be significant but appear to compensate each other to high degree. The calculated ∆c_p^hyd for duplexes and noncanonical DNA structures agree excellently with the overall experimental ∆c_{p, F} values. By contrast, empirical parametrizations developed for proteins result in poor ∆c_p^hyd predictions and should not be applied to DNA folding.

GENERAL SIGNIFICANCE

Heat capacity is one of the main thermodynamic quantities that strongly affects thermal stability of macromolecules. At the molecular level the heat capacity in DNA folding stems from removal of water from nucleobases.

Isothermal titration calorimetry: practical approaches and current applications in soft matter

Isothermal Titration Calorimetry (ITC) elucidates the thermodynamic profile (ΔH, ΔS, ΔG, Ka, and stoichiometry) of binding and dissociation reactions in solution. While ITC has primarily been used to investigate the thermodynamics of interactions between biological macromolecules and small molecules, it has become increasingly common for measuring binding interactions between synthetic polymers and small molecules, ions, or nanoparticles. This tutorial review describes applications of ITC in studying synthetic macromolecules and provides experimental guidelines for performing ITC experiments. We also highlight specific examples of using ITC to study soft matter, then discuss the limitations and the future of ITC in this field.

Dynamics of Ionic Interactions at Protein-Nucleic Acid Interfaces

Molecular association of proteins with nucleic acids is required for many biological processes essential to life. Electrostatic interactions via ion pairs (salt bridges) of nucleic acid phosphates and protein side chains are crucial for proteins to bind to DNA or RNA. Counterions around the macromolecules are also key constituents for the thermodynamics of protein–nucleic acid association. Until recently, there had been only a limited amount of experiment-based information about how ions and ionic moieties behave in biological macromolecular processes. In the past decade, there has been significant progress in quantitative experimental research on ionic interactions with nucleic acids and their complexes with proteins. The highly negatively charged surfaces of DNA and RNA electrostatically attract and condense cations, creating a zone called the ion atmosphere. Recent experimental studies were able to examine and validate theoretical models on ions and their mobility and interactions with macromolecules. The ionic interactions are highly dynamic. The counterions rapidly diffuse within the ion atmosphere. Some of the ions are released from the ion atmosphere when proteins bind to nucleic acids, balancing the charge via intermolecular ion pairs of positively charged side chains and negatively charged backbone phosphates. Previously, the release of counterions had been implicated indirectly by the salt-concentration dependence of the association constant.

Recently, direct detection of counterion release by NMR spectroscopy has become possible and enabled more accurate and quantitative analysis of the counterion release and its entropic impact on the thermodynamics of protein–nucleic acid association. Recent studies also revealed the dynamic nature of ion pairs of protein side chains and nucleic acid phosphates. These ion pairs undergo transitions between two major states. In one of the major states, the cation and the anion are in direct contact and form hydrogen bonds. In the other major state, the cation and the anion are separated by water. Transitions between these states rapidly occur on a picosecond to nanosecond time scale. When proteins interact with nucleic acids, interfacial arginine (Arg) and lysine (Lys) side chains exhibit considerably different behaviors. Arg side chains show a higher propensity to form rigid contacts with nucleotide bases, whereas Lys side chains tend to be more mobile at the molecular interfaces. The dynamic ionic interactions may facilitate adaptive molecular recognition and play both thermodynamic and kinetic roles in protein–nucleic acid interactions.

Structural basis of a novel repressor, SghR, controlling Agrobacterium infection by cross-talking to plants

Agrobacterium tumefaciens infects various plants and causes crown gall diseases involving temporal expression of virulence factors. SghA is a newly identified virulence factor enzymatically releasing salicylic acid from its glucoside conjugate and controlling plant tumor development. Here, we report the structural basis of SghR, a LacI-type transcription factor highly conserved in Rhizobiaceae family, regulating the expression of SghA and involved in tumorigenesis. We identified and characterized the binding site of SghR on the promoter region of sghA and then determined the crystal structures of apo-SghR, SghR complexed with its operator DNA, and ligand sucrose, respectively. These results provide detailed insights into how SghR recognizes its cognate DNA and shed a mechanistic light on how sucrose attenuates the affinity of SghR with DNA to modulate the expression of SghA. Given the important role of SghR in mediating the signaling cross-talk during Agrobacterium infection, our results pave the way for structure-based inducer analog design, which has potential applications for agricultural industry.

Disorder for Dummies: Functional Mutagenesis of Transient Helical Segments in Disordered Proteins

Most cytosolic eukaryotic proteins contain a mixture of ordered and disordered regions. Disordered regions facilitate cell signaling by concentrating sites for posttranslational modifications and protein-protein interactions into arrays of short linear motifs that can be reorganized by RNA splicing. The evolution of disordered regions looks different from their ordered counterparts. In some cases, selection is focused on maintaining protein binding interfaces and PTM sites, but sequence heterogeneity is common. In other cases, simple properties like charge, length, or end-to-end distance are maintained. Many disordered protein binding sites contain some transient secondary structure that may resemble the structure of the bound state. α-Helical secondary structure is common and a wide range of fractional helicity is observed in different disordered regions. Here we provide a simple protocol to identify transient helical segments and design mutants that can change their structure and function.