Mesoscale Modeling of Nucleosome-Binding Antibody PL2-6: Mono- versus Bivalent Chromatin Complexes

Interactions of chromatin with bivalent immunoglobin nucleosome-binding antibodies and their monovalent (papain-derived) antigen-binding fragment analogs are useful probes for examining chromatin conformational states. To help interpret antibody-chromatin interactions and explore how antibodies might compete for interactions with chromatin components, we incorporate coarse-grained PL2-6 antibody modeling into our mesoscale chromatin model. We analyze interactions and fiber structures for the antibody-chromatin complexes in open and condensed chromatin, with and without H1 linker histone (LH). Despite minimal and transient interactions at physiological salt, we capture significant differences in antibody-chromatin complex configurations in open fibers, with more intense interactions between the bivalent antibody and chromatin compared to monovalent antigen-binding fragments. For these open chromatin fiber morphologies, antibody binding to histone tails is increased and compaction is greater for bivalent compared to monovalent and antibody-free systems. Differences between monovalent and bivalent binding result from antibody competition with internal chromatin fiber components (nucleosome core and linker DNA) for histone tail (H3, H4, H2A, H2B) interactions. This antibody competition for tail contacts reduces tail-core and tail-linker interactions and increases tail-antibody interactions. Such internal structural changes in open fibers resemble mechanisms of LH condensation, driven by charge screening and entropy changes. For condensed fibers at physiological salt, the three systems are much more similar overall, but some subtle tail interaction differences can be noted. Adding LH results in less-dramatic changes for all systems, except that the bivalent complex at physiological salt shows cooperative effects between LH and the antibodies in condensing chromatin fibers. Such dynamic interactions that depend on the internal structure and complex-stabilizing interactions within the chromatin fiber have implications for gene regulation and other chromatin complexes such as with LH, remodeling proteins, and small molecular chaperones that bind and modulate chromatin structure.     

Dynamic condensation of linker histone C-terminal domain regulates chromatin structure

The basic and intrinsically disordered C-terminal domain (CTD) of the linker histone (LH) is essential for chromatin compaction. However, its conformation upon nucleosome binding and its impact on chromatin organization remain unknown. Our mesoscale chromatin model with a flexible LH CTD captures a dynamic, salt-dependent condensation mechanism driven by charge neutralization between the LH and linker DNA. Namely, at low salt concentration, CTD condenses, but LH only interacts with the nucleosome and one linker DNA, resulting in a semi-open nucleosome configuration; at higher salt, LH interacts with the nucleosome and two linker DNAs, promoting stem formation and chromatin compaction. CTD charge reduction unfolds the domain and decondenses chromatin, a mechanism in consonance with reduced counterion screening in vitro and phosphorylated LH in vivo. Divalent ions counteract this decondensation effect by maintaining nucleosome stems and expelling the CTDs to the fiber exterior. Additionally, we explain that the CTD folding depends on the chromatin fiber size, and we show that the asymmetric structure of the LH globular head is responsible for the uneven interaction observed between the LH and the linker DNAs. All these mechanisms may impact epigenetic regulation and higher levels of chromatin folding.     

Non-specific binding of Na+ and Mg2+ to RNA determined by force spectroscopy methods

RNA duplex stability depends strongly on ionic conditions, and inside cells RNAs are exposed to both monovalent and multivalent ions. Despite recent advances, we do not have general methods to quantitatively account for the effects of monovalent and multivalent ions on RNA stability, and the thermodynamic parameters for secondary structure prediction have only been derived at 1M [Na(+)]. Here, by mechanically unfolding and folding a 20 bp RNA hairpin using optical tweezers, we study the RNA thermodynamics and kinetics at different monovalent and mixed monovalent/Mg(2+) salt conditions. We measure the unfolding and folding rupture forces and apply Kramers theory to extract accurate information about the hairpin free energy landscape under tension at a wide range of ionic conditions. We obtain non-specific corrections for the free energy of formation of the RNA hairpin and measure how the distance of the transition state to the folded state changes with force and ionic strength. We experimentally validate the Tightly Bound Ion model and obtain values for the persistence length of ssRNA. Finally, we test the approximate rule by which the non-specific binding affinity of divalent cations at a given concentration is equivalent to that of monovalent cations taken at 100-fold concentration for small molecular constructs.     

Existence of efficient divalent metal ion-catalyzed and inefficient divalent metal ion-independent channels in reactions catalyzed by a hammerhead ribozyme

The hammerhead ribozyme is generally accepted as a well characterized metalloenzyme. However, the precise nature of the interactions of the RNA with metal ions remains to be fully defined. Examination of metal ion-catalyzed hammerhead reactions at limited concentrations of metal ions is useful for evaluation of the role of metal ions, as demonstrated in this study. At concentrations of Mn²⁺ ions from 0.3 to 3 mM, addition of the ribozyme to the reaction mixture under single-turnover conditions enhances the reaction with the product reaching a fixed maximum level. Further addition of the ribozyme inhibits the reaction, demonstrating that a certain number of divalent metal ions is required for proper folding and also for catalysis. At extremely high concentrations, monovalent ions, such as Na⁺ ions, can also serve as cofactors in hammerhead ribozyme-catalyzed reactions. However, the catalytic efficiency of monovalent ions is extremely low and, thus, high concentrations are required. Furthermore, addition of monovalent ions to divalent metal ion-catalyzed hammerhead reactions inhibits the divalent metal ion-catalyzed reactions, suggesting that the more desirable divalent metal ion–ribozyme complexes are converted to less desirable monovalent metal ion–ribozyme complexes via removal of divalent metal ions, which serve as a structural support in the ribozyme complex. Even though two channels appear to exist, namely an efficient divalent metal ion-catalyzed channel and an inefficient monovalent metal ion-catalyzed channel, it is clear that, under physiological conditions, hammerhead ribozymes are metalloenzymes that act via the significantly more efficient divalent metal ion-dependent channel. Moreover, the observed kinetic data are consistent with Lilley’s and DeRose’s two-phase folding model that was based on ground state structure analyses.     

Ciprofloxacin affects conformational equilibria of DNA gyrase A in the presence of magnesium ions

The conformational equilibria of the A subunit of DNA gyrase (GyrA), of its 59 kDa N-terminal fragment (GyrA59) and of the quinolone-resistant Ser-Trp83 mutant (GyrATrp83), were investigated in the presence of mono- and divalent metal ions and ciprofloxacin, a clinically useful antibacterial quinolone. The stability of the proteins was estimated from temperature denaturation, monitoring unfolding with circular dichroism spectroscopy. Two transitions were observed in GyrA and GyrATrp83, which likely reflect unfolding of the N and C-terminal protein domains. Accordingly, one thermal transition is observed for GyrA59.The melting profile of the GyrA subunit is dramatically affected by monovalent and divalent metal ions, both transitions being shifted to lower temperature upon increasing salt concentration. This effect is much more pronounced with divalent ions (Mg(2+)) and cannot be accounted for by changes in ionic strength only. The presence of ciprofloxacin shifts the melting transitions of the wild-type subunit to higher temperatures when physiological concentrations of Mg(2+) are present. In contrast, both the mutant protein and the 59 kDa fragment do not show evidence for quinolone-driven changes. These data suggest that ciprofloxacin binds to the wild-type subunit in an interaction that involves Ser83 of GyrA and that both C and N-terminal domains may be required for effective drug-protein interactions. The bell-shaped dependence of the binding process upon Mg(2+) concentration, with a maximum centred at 3-4 mM [Mg(2+)], is consistent with a metal-ion mediated GyrA-quinolone-interaction. Affinity chromatography data fully support these findings and additionally confirm the requirement for a free carboxylate to elicit binding of the quinolone to GyrA.We infer that the Mg(2+)-GyrA interaction at physiological metal ion concentration could bear biological relevance, conferring more conformational flexibility to the active enzyme. The results obtained in the presence of ciprofloxacin additionally suggest that the Mg(2+)-mediated quinolone binding to the enzyme might be involved in the mechanism of action of this family of drugs.     

Kinetics and thermodynamics of thermal denaturation in acyl carrier protein.

The denaturation of Escherichia coli acyl carrier protein (ACP) in buffers containing both monovalent and divalent cations was followed by variable-temperature NMR and differential scanning calorimetry. Both high concentrations of monovalent salts (Na+) and moderate concentrations of divalent salts (Ca2+) raise the denaturation temperature, but calorimetry indicates that a significant increase in the enthalpy of denaturation is obtained only with the addition of a divalent salt. NMR experiments in both low ionic strength monovalent buffers and low ionic strength monovalent buffers containing calcium ions show exchange between native and denatured forms to be slow on the NMR time scale. However, in high ionic strength monovalent buffers, where the temperature of denaturation is elevated as it is in the presence of Ca2+, the transition is fast on the NMR time scale. These results suggest that monovalent and divalent cations may act to stabilize ACP in different ways. Monovalent ions may nonspecifically balance the intrinsic negative charge of this protein in a way that is similar for native, denatured, and intermediate forms. Divalent cations provide stability by binding to specific sites present only in the native state.     

Daunomycin-induced unfolding and aggregation of chromatin.

Using equilibrium dialysis and sedimentation velocity analysis, we have characterized the binding of the anti-tumor drug daunomycin to chicken erythrocyte chromatin before and after depletion of linker histones and to its constitutive DNA under several ionic strengths (5, 25, and 75 mM NaCl). The equilibrium dialysis experiments reveal that the drug binds cooperatively to both the chromatin fractions and to the DNA counterpart within the range of ionic strength used in this study. A significant decrease in the binding affinity was observed at 75 mM NaCl. At any given salt concentration, daunomycin exhibits higher binding affinity for DNA than for linker histone-depleted chromatin or chromatin (in decreasing order). Binding of daunomycin to DNA does not significantly affect the sedimentation coefficient of the molecule. This is in contrast to binding to chromatin and to its linker histone-depleted counterpart. In these instances, preferential binding of the drug to the linker DNA regions induces an unfolding of the chromatin fiber that is followed by aggregation, presumably because of histone-DNA interfiber interactions.     

Competitive Binding of Mg2+ and Na+ Ions to Nucleic Acids: From Helices to Tertiary Structures

Nucleic acids generally reside in cellular aqueous solutions with mixed divalent/monovalent ions, and the competitive binding of divalent and monovalent ions is critical to the structures of nucleic acids because of their polyanionic nature. In this work, we first proposed a general and effective method for simulating a nucleic acid in mixed divalent/monovalent ion solutions with desired bulk ion concentrations via molecular dynamics (MD) simulations and investigated the competitive binding of Mg²⁺/Na⁺ ions to various nucleic acids by all-atom MD simulations. The extensive MD-based examinations show that single MD simulations conducted using the proposed method can yield desired bulk divalent/monovalent ion concentrations for various nucleic acids, including RNA tertiary structures. Our comprehensive analyses show that the global binding of Mg²⁺/Na⁺ to a nucleic acid is mainly dependent on its structure compactness, as well as Mg²⁺/Na⁺ concentrations, rather than the specific structure of the nucleic acid. Specifically, the relative global binding of Mg²⁺ over Na⁺ is stronger for a nucleic acid with higher effective surface charge density and higher relative Mg²⁺/Na⁺ concentrations. Furthermore, the local binding of Mg²⁺/Na⁺ to a phosphate of a nucleic acid mainly depends on the local phosphate density in addition to Mg²⁺/Na⁺ concentrations.     

Modeling Studies of Chromatin Fiber Structure as a Function of DNA Linker Length

Chromatin fibers encountered in various species and tissues are characterized by different nucleosome repeat lengths (NRLs) of the linker DNA connecting the nucleosomes. While single cellular organisms and rapidly growing cells with high protein production have short NRL ranging from 160 to 189 bp, mature cells usually have longer NRLs ranging between 190 and 220 bp. Recently, various experimental studies have examined the effect of NRL on the internal organization of chromatin fiber. Here, we investigate by mesoscale modeling of oligonucleosomes the folding patterns for different NRL, with and without linker histone (LH), under typical monovalent salt conditions using both one-start solenoid and two-start zigzag starting configurations. We find that short to medium NRL chromatin fibers (173 to 209 bp) with LH condense into zigzag structures and that solenoid-like features are viable only for longer NRLs (226 bp). We suggest that medium NRLs are more advantageous for packing and various levels of chromatin compaction throughout the cell cycle than their shortest and longest brethren; the former (short NRLs) fold into narrow fibers, while the latter (long NRLs) arrays do not easily lead to high packing ratios due to possible linker DNA bending. Moreover, we show that the LH has a small effect on the condensation of short-NRL arrays but has an important condensation effect on medium-NRL arrays, which have linker lengths similar to the LH lengths. Finally, we suggest that the medium-NRL species, with densely packed fiber arrangements, may be advantageous for epigenetic control because their histone tail modifications can have a greater effect compared to other fibers due to their more extensive nucleosome interaction network.     

Role of divalent ions in folding of tRNA.

The native structure of tRNA is not achieved in low salt (4.5 mM Na+, 25 degrees C), but can be restored by addition of divalent ions. We have explored the structure of the central region in Escherichia coli tRNAfMet by absorption and emission spectroscopy of 4-thiouracil, and the structure of the anticodon loop in yeast tRNAPhe by fluorescence of the 'Y' base, versus the number of manganese ions bound to tRNA, which was derived from electron spin resonance. The fluorescence of the reduced 8-13 photoproduct (in which 4-thiouracil at position 8 is crosslinked to cytosine at position 13) was also analysed. In low salt (e.g. 4.5 mM Na+), the region of 4-thiouracil is affected strongly as the first eight Mn2+ bind to tRNA, whereas the fluorescence of the 'Y' base is affected only after four Mn2+ are bound. Considering the structural similarities of the two tRNAs, this suggests that the reorganisation brought about by divalent ions starts in the central region, the anticodon loop being affected later. The binding of divalent ions to each region starts together with its restructuration. Monovalent ions can substitute for divalent ions in this process, a 15 mM sodium concentration being equivalent to the binding of the first five Mn2+. If divalent ions are then added, even the first ones distribute themselves between both the central and the anticodon region. Alternatively, the renaturation may be achieved by monovalent ions only, implying that no sites exist whose occupancy by divalent ions is crucial for the native structure. These observations suggest that the role and means of divalent ion binding to tRNA are largely explainable in terms of a simple maganese-phosphate binding supplemented by electrostatic interaction with distant phosphates.     

DNA like-charge attraction and overcharging by divalent counterions in the presence of divalent co-ions

Strongly correlated electrostatics of DNA systems has drawn the interest of many groups, especially the condensation and overcharging of DNA by multivalent counterions. By adding counterions of different valencies and shapes, one can enhance or reduce DNA overcharging. In this paper, we focus on the effect of multivalent co-ions, specifically divalent co-ions such as SO\(_{4}^{2-}\). A computational experiment of DNA condensation using Monte Carlo simulation in grand canonical ensemble is carried out where the DNA system is in equilibrium with a bulk solution containing a mixture of salt of different valency of co-ions. Compared to systems with purely monovalent co-ions, the influence of divalent co-ions shows up in multiple aspects. Divalent co-ions lead to an increase of monovalent salt in the DNA condensate. Because monovalent salts mostly participate in linear screening of electrostatic interactions in the system, more monovalent salt molecules enter the condensate leads to screening out of short-range DNA–DNA like charge attraction and weaker DNA condensation free energy. The overcharging of DNA by multivalent counterions is also reduced in the presence of divalent co-ions. Strong repulsions between DNA and divalent co-ions and among divalent co-ions themselves lead to a depletion of negative ions near the DNA surface as compared to the case without divalent co-ions. At large distances, the DNA–DNA repulsive interaction is stronger in the presence of divalent co-ions, suggesting that divalent co-ions’ role is not only that of simple stronger linear screening.     

Effects of various ionic media on extraction of soluble nuclear proteins and on nuclear ultrastructure

Isolated liver nuclei were extracted 3 times at pH 7.2 with solutions containing either (1) monovalent cations, (2) both mono- and divalent cations, or (3) sucrose solutions containing only divalent cations. The extracted proteins were analysed by two-dimensional acrylamide gel electrophoresis and the ultrastructural alterations of the treated nuclei were examined by electron microscopy. The solutions containing Na⁺ or K⁺ monovalent and Ca²⁺ and Mg²⁺ divalent ions extracted the same amount (18–22 %) of the nuclear proteins. The two-dimensional gel electrophoretic patterns of these extracts were nearly identical and the structures of the nuclear components were well preserved even after 3 times repeated extractions. The solution containing only Na⁺ extracted less protein (14–15 %) than the solutions containing both mono- and divalent cations. Extraction with isotonic NaCl solution altered the nuclear and nucleolar morphology; unlike the other solutions employed, this solution extracted some DNA and histones. The isotonic sucrose solution containing only divalent cations extracted less protein than the other solutions (9–11 %) and produced marked condensation of the chromatin. These analytical and electron microscopic studies showed that mono- and divalent cations play a role in structural organization of chromatin.     

Differential stability of DNA crossovers in solution mediated by divalent cations

The assembly of DNA duplexes into higher-order structures plays a major role in many vital cellular functions such as recombination, chromatin packaging and gene regulation. However, little is currently known about the molecular structure and stability of direct DNA–DNA interactions that are required for such functions. In nature, DNA helices minimize electrostatic repulsion between double helices in several ways. Within crystals, B-DNA forms either right-handed crossovers by groove–backbone interaction or left-handed crossovers by groove–groove juxtaposition. We evaluated the stability of such crossovers at various ionic concentrations using large-scale atomistic molecular dynamics simulations. Our results show that right-handed DNA crossovers are thermodynamically stable in solution in the presence of divalent cations. Attractive forces at short-range stabilize such crossover structures with inter-axial separation of helices less than 20 Å. Right-handed crossovers, however, dissociate swiftly in the presence of monovalent ions only. Surprisingly, left-handed crossovers, assembled by sequence-independent juxtaposition of the helices, appear unstable even at the highest concentration of Mg²⁺studied here. Our study provides new molecular insights into chiral association of DNA duplexes and highlights the unique role divalent cations play in differential stabilization of crossover structures. These results may serve as a rational basis to understand the role DNA crossovers play in biological processes.     

An energy dependent divalent cation-hydrogen ion binding in the outer membranes of rat liver mitochondria and its dependence on monovalent cation-hydrogen ion binding

Intact mitochondria are able to bind monovalent and divalent metal cations and to release protons in an energy-independent exchange process. Directly accessible binding sites exist in the outer membrane. They seem to be identical for monovalent and divalent metal ions. The inner membrane-matrix-fraction possesses exchange sites after ultrasonic disruption only for monovalent cations, but not for divalent cations.     

Mechanism and Role of Divalent Cation Binding of Bacteriorhodopsin

Several observations have already suggested that the carboxyl groups are involved in the association of divalent cations with bacteriorhodopsin (Chang et al., 1985). Here we show that at least part of the protons released from deionized purple membrane (`blue membrane') samples when salt is added are from carboxyl groups. We find that the apparent pK of magnesium binding to purple membrane in the presence of 0.5 mM buffer is 5.85. We suggest this is the pK of the carboxyl groups shifted from their usual pK because of the proton concentrating effect of the large negative surface potential of the purple membrane. Divalent cations may interact with negatively charged sites on the surface of purple membrane through the surface potential and/or through binding either by individual ligands or by conformation-dependent chelation. We find that divalent cations can be released from purple membrane by raising the temperature. Moreover, purple membrane binds only about half as many divalent cations after bleaching. Neither of these operations is expected to decrease the surface potential and thus these experiments suggest that some specific conformation in purple membrane is essential for the binding of a substantial fraction of the divalent cations. Divalent cations in purple membrane can be replaced by monovalent, (Na⁺ and K⁺), or trivalent, (La⁺⁺⁺) cations. Flash photolysis measurements show that the amplitude of the photointermediate, O, is affected by the replacement of the divalent cations by other ions, especially by La⁺⁺⁺. The kinetics of the M photointermediate and light-induced H⁺ uptake are not affected by Na⁺ and K⁺, but they are drastically lengthened by La⁺⁺⁺ substitution, especially at alkaline pHs. We suggest that the surface charge density and thus the surface potential is controlled by divalent cation binding. Removal of the cations (to make deionized blue membrane) or replacement of them (e.g. La⁺⁺⁺-purple membrane) changes the surface potential and hence the proton concentration near the membrane surface. An increase in local proton concentration could cause the protonation of critical carboxyl groups, for example the counter-ion to the protonated Schiff's base, causing the red shift associated with the formation of both deionized and acid blue membrane. Similar explanations based on regulation of the surface proton concentration can explain many other effects associated with the association of different cations with bacteriorhodopsin.     

Structural basis of filamin A functions

Filamin A (FLNa) can effect orthogonal branching of F-actin and bind many cellular constituents. FLNa dimeric subunits have N-terminal spectrin family F-actin binding domains (ABDs) and an elongated flexible segment of 24 immunoglobulin (Ig) repeats. We generated a library of FLNa fragments to examine their F-actin binding to define the structural properties of FLNa that enable its various functions. We find that Ig repeats 9–15 contain an F-actin–binding domain necessary for high avidity F-actin binding. Ig repeats 16–24, where most FLNa-binding partners interact, do not bind F-actin, and thus F-actin does not compete with Ig repeat 23 ligand, FilGAP. Ig repeats 16–24 have a compact structure that suggests their unfolding may accommodate pre-stress–mediated stiffening of F-actin networks, partner binding, mechanosensing, and mechanoprotection properties of FLNa. Our results also establish the orientation of FLNa dimers in F-actin branching. Dimerization, mediated by FLNa Ig repeat 24, accounts for rigid high-angle FLNa/F-actin branching resistant to bending by thermal forces, and high avidity F-actin binding and cross-linking.