The G₄C₂ hexanucleotide repeat expansion mutation (HREM) in C9ORF72, represents the most common mutation associated with amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD). Three main disease mechanisms have been proposed to date: C9ORF72 haploinsufficiency, RNA toxicity and accumulation of dipeptide repeat proteins. Pure GC content of the HREM potentially enables the formation of various non-B DNA structures such as G-quadruplexes and i-motifs. These structures are proposed to act as promoters and regulatory elements affecting replication, transcription and translation of the surrounding region. G-quadruplexes have already been shown on the G-rich sense DNA and RNA strands (G₄C₂)_n, the structure of the anti-sense (G₂C₄)_n strand remains unresolved. Similar C-rich sequences may, under acidic conditions, form i-motifs consisting of two parallel duplexes in a head to tail orientation held together by hemi-protonated C⁺-C pairs. We show that d(G₂C₄)_n repeats do form i-motif and protonated hairpins even under near-physiological conditions. Rather than forming a DNA duplex, i-motifs persist even in the presence of the sense strand. This preferential formation of G-quadruplex and i-motif/hairpin structures over duplex DNA, may explain HREM replicational and transcriptional instability. Furthermore, i-motifs/hairpins can represent a novel pharmacological target for C9ORF72 associated ALS and FTLD.

Adenosine tripolyphosphate (ATP) is a small polyvalent anion that has recently been shown to interact with proteins and have a major impact on assembly processes involved in biomolecular condensate formation and protein aggregation. However, the nature of non-specific protein-ATP interactions and their effects on protein solubility are largely unknown. Here, the binding of ATP to the globular model protein is characterized in detail using X-ray crystallography and nuclear magnetic resonance (NMR). Using NMR, we identified six ATP binding sites on the lysozyme surface, with one known high-affinity nucleic acid binding site and five non-specific previously unknown sites with millimolar affinities that also bind tripolyphosphate (TPP). ATP binding occurs primarily through the polyphosphate moiety, which was confirmed by the X-ray structure of the lysozyme-ATP complex. Importantly, ATP binds preferentially to arginine over lysine in non-specific binding sites. ATP and TPP have similar effects on solution-phase protein-protein interactions. At low salt concentrations, ion binding to lysozyme causes precipitation, while at higher salt concentrations, redissolution occurs. The addition of an equimolar concentration of magnesium to ATP does not alter ATP binding affinities but prevents lysozyme precipitation. These findings have important implications for both protein crystallization and cell biology. Crystallization occurs readily in ATP solutions outside the well-established crystallization window. In the context of cell biology, the findings suggest that ATP binds non-specifically to folded proteins in physiological conditions. Based on the nature of the binding sites identified by NMR, we propose several mechanisms for how ATP binding can prevent the aggregation of natively folded proteins.

Integral equation approaches, based on the Ornstein-Zernike equation, provide a fast way to calculate phase diagrams and thermodynamic properties of systems as opposed to time-consuming and computationally expensive computer simulations. However, when employing integral equations it is necessary to introduce simplifications. The Ornstein-Zernike equation merely relates two unknown functions h(r) and c(r), and another relation (closer) between these two functions is needed. The later function cannot be obtained in a closed form and it is always in some approximations. Various approximations exist with each of its own advantages and disadvantages. In this work we extensively tested hyper-netted chain, Percus-Yevick, Kovalenko-Hirata, and Rogers-Young closure on an interaction model with core-softened potential. Convergence domain was established for each method. We calculated pair distribution functions, pressure, and excess energy. Results were compared with Monte Carlo simulation results and literature data from molecular dynamics simulations.