Modeling of minimal systems based on ATP-Zn coordination for chemically fueled self-assembly

Following nature's example, there is currently strong interest in using adenosine 5'-triphosphate (ATP) as a fuel for the self-assembly of functional materials with transient/non-equilibrium behaviours. These hold great promise for applications, e.g. in catalysis and drug delivery. In a recent seminal work [Maiti et al., Nat. Chem., 2016, 8, 725], binding of ATP to the metallosurfactant zinc hexadecyl-1,4,7-triazacyclononane ([ZnC₁₆ TACN]²⁺) was exploited to produce ATP-fueled transient vesicles. Crucial to the complex formation is the ability of ATP to bind to the metal ion. As a first step to unveil the key elements underlying this process, we investigate the interaction of ATP with Zn²⁺ and with methyl-1,4,7-triazacyclononane ([ZnCH₃ TACN]²⁺), using all-atom molecular dynamics simulations. The free energy landscape of the complex formation is sampled using well-tempered metadynamics with three collective variables, corresponding to the coordination numbers of Zn²⁺ with the oxygen atoms of the three phosphate groups. We find that the structure of the ternary complex is controlled by direct triphosphate coordination to zinc, with a minor role played by the interactions between ATP and CH₃ TACN which, however, may be important for the build-up of supramolecular assemblies.

Short-Range Cooperative Slow-down of Water Solvation Dynamics Around SO4 2--Mg2+ Ion Pairs

The presence of ions affects the structure and dynamics of water on a multitude of length and time scales. In this context, pairs of Mg²⁺ and SO₄ ^2- ions in water constitute a prototypical system for which conflicting pictures of hydration geometries and dynamics have been reported. Key issues are the molecular pair and solvation shell geometries, the spatial range of electric interactions, and their impact on solvation dynamics. Here, we introduce asymmetric SO₄ ^2- stretching vibrations as new and most specific local probes of solvation dynamics that allow to access ion hydration dynamics at the dilute concentration (0.2 M) of a native electrolyte environment. Highly sensitive heterodyne 2D-IR spectroscopy in the fingerprint region of the SO₄ ^2- ions around 1100 cm^-1 reveals a specific slow-down of solvation dynamics for hydrated MgSO₄ and for Na₂SO₄ in the presence of Mg²⁺ ions, which manifests as a retardation of spectral diffusion compared to aqueous Na₂SO₄ solutions in the absence of Mg²⁺ ions. Extensive molecular dynamics and density functional theory QM/MM simulations provide a microscopic view of the observed ultrafast dephasing and hydration dynamics. They suggest a molecular picture where the slow-down of hydration dynamics arises from the structural peculiarities of solvent-shared SO₄ ^2--Mg²⁺ ion pairs.

Sodium site occupancy and phosphate speciation in natrophosphate are invariant to changes in NaF and Na3PO4 concentration

Analysis of multimodal characterization of Natrophosphate suggests that the crystalline structure is preserved across a range of synthesis conditions.

The mutual interactions of RNA, counterions and water - quantifying the electrostatics at the phosphate-water interface

The structure and dynamics of polyanionic biomolecules, like RNA, are decisively determined by their electric interactions with the water molecules and the counterions in the environment. The solvation dynamics of the biomolecules involves a subtle balance of non-covalent and many-body interactions with structural fluctuations due to thermal motion occurring in a femto- to subnanosecond time range. This complex fluctuating many particle scenario is crucial in defining the properties of biological interfaces with far reaching significance for the folding of RNA structures and for facilitating RNA-protein interactions. Given the inherent complexity, suited model systems, carefully calibrated and benchmarked by experiments, are required to quantify the relevant interactions of RNA with the aqueous environment. In this feature article we summarize our recent progress in the understanding of the electrostatics at the biological interface of double stranded RNA (dsRNA) and transfer RNA (tRNA). Dimethyl phosphate (DMP) is introduced as a viable and rigorously accessible model system allowing the interaction strength with water molecules and counterions, their relevant fluctuation timescales and the spatial reach of interactions to be established. We find strong (up to ≈90 MV cm-1) interfacial electric fields with fluctuations extending up to ≈20 THz and demonstrate how the asymmetric stretching vibration νAS(PO2)- of the polarizable phosphate group can serve as the most sensitive probe for interfacial interactions, establishing a rigorous link between simulations and experiment. The approach allows for the direct interfacial observation of interactions of biologically relevant Mg2+ counterions with phosphate groups in contact pair geometries via the rise of a new absorption band imposed by exchange repulsion interactions at short interatomic distances. The systematic extension to RNA provides microscopic insights into the changes of the hydration structure that accompany the temperature induced melting of the dsRNA double helix and quantify the ionic interactions in the folded tRNA. The results show that pairs of negatively charged phosphate groups and Mg2+ ions represent a key structural feature of RNA embedded in water. They highlight the importance of binding motifs made of contact pairs in the electrostatic stabilization of RNA structures that have a strong impact on the surface potential and enable the fine tuning of the local electrostatic properties which are expected to be relevant for mediating the interactions between biomolecules.

Unravelling the dissolution mechanism of polyphosphate glasses by 31P NMR spectroscopy: ligand competition and reactivity of intermediate complexes

Phosphate glass dissolution can be tailored via compositional and subsequent structural changes, which is of interest for biomedical applications such as therapeutic ion delivery. Here, solid-state ³¹P nuclear magnetic resonance characterisation of 45P₂O₅-xCaO - (55 -x)Na₂O glasses was correlated with dissolution studies using time-dependent liquid ³¹P NMR spectroscopy and quantitative chemical analysis. Glasses dissolved congruently in aqueous media, and the first dissolution stage was the hydration of phosphate chains. In deionised water and Tris buffer (pH₀ 7.4 or 7.9), trimetaphosphate rings and orthophosphates were the predominant species in solution, indicating relatively fast degradation. By contrast, long phosphate chains were identified in EDTA (pH₀ 10.0). Besides pH differences, coordination of phosphate species by metal cations appears to play a catalytic role in the hydrolysis mechanism via turning phosphorus atoms into suitable electrophiles for the subsequent nucleophilic attack by water. Hydrolysis rates were proportional to phosphate complex stability, with stronger complexes for chains than for rings. A competition between solvent and phosphate species for the metal ion occurred in the order EDTA > Tris > deionised water.

Magnesium Contact Ions Stabilize the Tertiary Structure of Transfer RNA: Electrostatics Mapped by Two-Dimensional Infrared Spectra and Theoretical Simulations

Ions interacting with hydrated RNA play a central role in defining its secondary and tertiary structure. While spatial arrangements of ions, water molecules, and phosphate groups have been inferred from X-ray studies, the role of electrostatic and other noncovalent interactions in stabilizing compact folded RNA structures is not fully understood at the molecular level. Here, we demonstrate that contact ion pairs of magnesium (Mg²⁺) and phosphate groups embedded in local water shells stabilize the tertiary equilibrium structure of transfer RNA (tRNA). Employing dialyzed tRNA^Phe from yeast and tRNA from Escherichia coli, we follow the population of Mg²⁺ sites close to phosphate groups of the ribose-phosphodiester backbone step by step, combining linear and nonlinear infrared spectroscopy of phosphate vibrations with molecular dynamics simulations and ab initio vibrational frequency calculations. The formation of up to six Mg²⁺/phosphate contact pairs per tRNA and local field-induced reorientations of water molecules balance the phosphate–phosphate repulsion in nonhelical parts of tRNA, thus stabilizing the folded structure electrostatically. Such geometries display limited sub-picosecond fluctuations in the arrangement of water molecules and ion residence times longer than 1 μs. At higher Mg²⁺ excess, the number of contact ion pairs per tRNA saturates around 6 and weakly interacting ions prevail. Our results suggest a predominance of contact ion pairs over long-range coupling of the ion atmosphere and the biomolecule in defining and stabilizing the tertiary structure of tRNA.