research papers
Protonated nucleobases are not fully ionized in their chloride salt crystals and form metastable base pairs further stabilized by the surrounding anions
aBiological and Chemical Research Center, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa 02-089, Poland
*Correspondence e-mail: pdomin@chem.uw.edu.pl
This paper presents experimental charge-density studies of cytosinium chloride, adeninium chloride hemihydrate and guaninium dichloride crystals based on ultra-high-resolution X-ray diffraction data and extensive theoretical calculations. The results confirm that the cohesive energies of the studied systems are dominated by contributions from intermolecular electrostatic interactions, as expected for ionic crystals. Electrostatic interaction energies (Ees) usually constitute 95% of the total interaction energy. The Ees energies in this study were several times larger in absolute value when compared, for example, with dimers of neutral nucleobases. However, they were not as large as some theoretical calculations have predicted. This was because the molecules appeared not to be fully ionized in the studied crystals. Apart from charge transfer from chlorine to the protonated nucleobases, small but visible charge redistribution within the nucleobase cations was observed. Some dimers of singly protonated bases in the studied crystals, namely a cytosinium–cytosinium trans sugar/sugar edge pair and an adeninium–adeninium trans Hoogsteen/Hoogsteen edge pair, exhibited attractive interactions (negative values of Ees) or unusually low repulsion despite identical molecular charges. The pairs are metastable as a result of strong hydrogen bonding between bases which overcompensates the overall cation–cation repulsion, the latter being weakened due to charge transfer and molecular charge-density polarization.
Keywords: charge density; multipole refinement; quantum crystallography; electrostatic potential; electrostatic interaction energy; like-charged ions; nucleobases; cytosine; adenine; guanine; base pairs; intermolecular interactions; spin density; momentum density; hydrogen bonding; pharmaceutical solids; nucleic acid structures.
1. Introduction
Nucleobases are naturally occurring chemical compounds that, together with phosphoric acid and sugars, constitute et al., 2002). Because of their crucial roles, nucleobases have been of considerable scientific interest since their discovery (Avery et al., 1944; Watson & Crick, 1953a,b). Numerous reports concerning the noncovalent interactions of nucleobases (in general, hydrogen-bonding and π–π stacking interactions between bases, and cation–anion interactions with surrounding molecules) create an ever-expanding area of research in the scientific community (Cech, 1993a,b,c; Cantor et al., 1980). It is also now well established that nucleobases in are prone to forming not only canonical Watson–Crick pairs but also a variety of non-canonical pairs (Leontis & Westhof, 2001; Leontis et al., 2002).
the main information-carrying molecules of the cell which, by directing protein synthesis, determine the inherited characteristics of every living organism. However, they are also required for numerous other important functions within the cell, such as catalysing biological reactions or sensing and transmitting responses to cellular signals (AlbertsIn the last two decades, the significance of the protonation (and deprotonation) of nucleobases for the biological functions of et al., 2015, 2014; Gehring et al., 1993; Berger et al., 1995; Chawla et al., 2011). It was found that many nucleobases have elevated pKa in and can be protonated at physiological pH (Acharya et al., 2004; Goh et al., 2013; Muth et al., 2000; Pechlaner et al., 2015; Siegfried et al., 2010; Tang et al., 2007; Ward et al., 2014; Wilcox & Bevilacqua, 2013; Wolter et al., 2017). For example, a near-neutral pKa of 7.6 ± 0.2 was observed for a single adenosine in the ribosomal peptidyl transferase centre, which indicates a shift of ∼4.0 pKa units from the unperturbed pKa (Muth et al., 2000).
has been recognized (HalderProtonated nucleobases play an important role in defining three-dimensional structures, mediating conformational changes, controlling proper folding and stability, and facilitating catalysis (Bevilacqua et al., 2004; Asensio et al., 1998; Wilcox et al., 2011). Changes in the protonation state, and resulting changes in local charge, have a significant influence on the propensity of nucleobases to form specific noncovalent interactions.
In the light of the above observations, we found it interesting to investigate the interactions of protonated bases with themselves and with anions surrounding them from a charge-density point of view. We focused our investigations on how nucleobases assemble themselves to form small-molecule crystals. In the present work, we show the results of our research on cytosinium chloride (CC), adeninium chloride hemihydrate (ACH) and guaninium dichloride (GDC) salts in the crystalline state for which it was possible to determine crystal charge densities experimentally. The crystal structures of CC, ACH and GDC were reported previously (Matković-Čalogović & Sanković, 1999; Broomhead, 1948; Cochran, 1951; Kistenmacher & Shigematsu, 1974; Cunane & Taylor, 1993; Mandel, 1977; Nieger, 2006). In addition, Cunane & Taylor (1993) published high-resolution X-ray diffraction data (sinθ/λ = 1.32 Å−1 at 123 K) for a single crystal of ACH. This first charge-density study mainly presented the strategy and procedure of the studied system as well as its molecular electrostatic potential features, which was regarded as a great success at the time. For the other two, CC and GDC, only structural parameters have been reported until now.
The present paper is the first in a series of two. The series is dedicated to the qualitative and quantitative study of the interactions of the chosen protonated nucleobases with themselves, chlorine anions or water molecules by X-ray diffraction complemented with quantum mechanics calculations. Our aim is to find the relationship between molecular and crystal architecture, and analyse the nature of the intermolecular interactions. Here, we show detailed structural descriptions of the molecular motifs found in the studied crystals, to give a basis for further charge-density, electrostatic potential and interaction energy considerations. Then we focus on net molecular charges, chloride to protonated nucleobase charge transfer, molecular, dimer and crystal electrostatic potentials, and charge redistribution within nucleobases. We finish with an analysis of electrostatic interaction energies for dimers and whole crystals in order to relate, in a quantitative manner, the observed molecular electrostatic properties with the energies. Among the dimers analysed, in the present paper we focus on protonated nucleobase pairs. The next paper of the series will give a more detailed analysis of the interactions present in a variety of molecular dimers identified in the studied crystals. A comprehensive topological analysis of atom–atom intermolecular bonding and of the nature of electrostatic intermolecular interactions with respect to molecular and atomic multipole moments and charge-density penetration phenomena will be presented there. The analysis reveals how short-range (at single-atom, functional-group and whole-molecule levels) and long-range (at the level of second-nearest neighbours and further away) interactions contribute to the stability of the studied ionic crystals. The series aims to give a wide perspective on the electrostatic interactions which contribute to supramolecular assembly and thus on their crucial role in crystal engineering and, perhaps, structural biology.
2. Experimental and computational methods
For a detailed description of the methods used, please see the supporting information.
2.1. Experimental charge-density models
2.1.1. Crystallization, data collection and processing
The nucleobases were each added to a mixture of distilled water with a few drops of 38% HCl, and the mixtures were stirred and heated until the compounds had dissolved completely. The solutions were left to evaporate at room temperature in the case of cytosinium chloride (CC), and at 310 K in the cases of adeninium chloride hemihydrate (ACH) and guaninium dichloride (GDC). Crystals were obtained after one month.
For CC and ACH, single-crystal X-ray measurements were performed at 90 K on an Agilent Technologies SuperNova four-circle diffractometer equipped with a low-temperature nitrogen gas-flow device (Oxford Cryosystems Cryostream Plus). For guaninium dichloride GDC, X-ray measurements were carried out on a Bruker APEXII ULTRA single-crystal diffractometer with a TXS rotating anode (Mo Kα) equipped with a CCD-type area detector, multilayer optics and an Oxford Cryostream low-temperature attachment set to 100 K.
For CC and ACH, the determination of unit-cell parameters, integration of reflection intensities and data reduction, including multi-scan absorption corrections, were performed using CrysAlis PRO (Version 1.171.36.32; Agilent, 2013). Finally, reflection merging was carried out with the SORTAV program (Blessing, 1987, 1989, 1995, 1997). For GDC, the determination of unit-cell parameters, integration of reflection intensities and data reduction were performed with the APEX2 suite of programs (integration was carried out with SAINT, Version 8.27B; Bruker, 2013), and the multi scan absorption correction, scaling and merging of reflection data were carried out with the SORTAV program.
2.1.2. Structure solution and refinement
Using OLEX2 (Dolomanov et al., 2009), the structures were solved with SHELXS (Sheldrick, 2008, 2015) with and refined with olex2.refine using the independent-atom model (IAM).
Multipole refinements were performed in MoPro (Guillot et al., 2001; Jelsch et al., 2005) with the use of the Stewart–Hansen–Coppens multipolar model (Stewart et al., 1975; Hansen & Coppens, 1978). Refinements were performed against amplitudes (F) with the I ≥ 2σ(I) cut-off. The initial atomic coordinates, x, y and z, for all atoms, the anisotropic atomic displacement parameters (Uij) for non-hydrogen atoms and the isotropic atomic displacement parameters for hydrogen atoms were taken from the IAM Local Cartesian coordinate systems and the initial multipolar and contraction–expansion parameters for nucleobases and water molecules were defined by LSDB (Volkov, Li et al., 2004) combined with UBDB2011 (Jarzembska & Dominiak, 2012). Each atom was assigned core and spherical valence scattering factors derived from the Su & Coppens (1998) atomic wavefunctions for neutral-atom configurations, except for the chlorine atoms. For chlorine atoms two possibilities were investigated: the Cl−1 ion scattering radial function and ion or the Cl0 neutral scattering radial function and neutral The coefficients were taken from Kissel et al. (1995). The contraction–expansion parameters, κ and κ′, for all hydrogen atoms were kept fixed at the UBDB2011 values during The κ′ parameters of the chlorine atoms were restrained to 1.0. The values of the Uij parameters for the hydrogen atoms were estimated from the SHADE 3.0 server (Madsen, 2006). The X—H distances were restrained to the average values obtained from neutron diffraction studies (Allen & Bruno, 2010) with a restraint σ of 0.004 Å. As atom Cl1 of ACH undergoes noticeable anharmonic motion, Gram–Charlier (GC) coefficients (Kuhs, 1983; Johnson, 1969; Scheringer, 1985) up to the third order were used to model it, while the physical reliability of the anharmonic model was confirmed by the probability density function computed with MolIso (Hübschle & Luger, 2006).
The outcomes of the multipole refinements were verified by examining the R factors, goodness of fit (GOF) and residual densities (Table 1). The evaluation was additionally supported by JNK2RDA (Meindl & Henn, 2008) and XDRKplot implemented in the WinGX package (Farrugia, 2012) (see Figs. S1–S4 in the supporting information). PLATON diagrams (Spek, 2009) with the atom-labelling schemes are shown in Fig. 1. The files can be retrieved from the Cambridge Structural Database (CSD) (Groom et al., 2016) (deposition numbers 1539172–1539174). Sets of raw diffraction frames and the associated data are accessible online under the following DOIs: https://doi.org/10.18150/repod.7313736, https://doi.org/10.18150/repod.0481200 and https://doi.org/10.18150/repod.8020589 (Repository for Open Data, ICM, University of Warsaw, Poland).
|
2.2. UBDB charge-density models
UBDB models of charge densities for CC, ACH and GDC crystals were built on experimental geometries using LSDB (Volkov, Li et al., 2004) and UBDB2011 (Jarzembska & Dominiak, 2012). Chlorine atoms were treated as spherical anions with a Cl−1 ion scattering radial function and ion Core and spherical valence factors for each atom were taken from Su and Coppens' atomic wavefunctions (Su & Coppens, 1998). For the majority of calculations, the net charge of each molecule was scaled separately to its formal value: +1 e for cytosinium and adeninium, +2 e for guaninium, 0 e for water molecules and −1 e for chloride anions. For alternative simulations (as explained further in the text), net molecular charges were scaled to experimentally observed values.
2.3. Periodic quantum mechanical calculations
Periodic quantum mechanical calculations were applied to compute theoretical crystal charge densities, theoretical structure factors and crystal cohesive energies (Ecoh). The CRYSTAL14 package (Dovesi, Orlando et al., 2014; Dovesi, Saunders et al., 2014) was used at the DFT-B3LYP/pVDZ level of theory (Becke, 1988; Perdew, 1986; Lee et al., 1988; Dunning, 1989). The computations were done in two versions, for experimental and optimized geometries. During optimization, the cell parameters were fixed while the atom positions were allowed to vary. The Grimme dispersion correction (Grimme, 2006) was applied for all calculations, whereas a correction for basis set superposition error (BSSE) (Civalleri et al., 2008) was applied for the computation of cohesive energy following the supermolecular approach (Maschio et al., 2011). The irreducible was sampled using 170 k points (the shrinking factor of the net was set to 8). See Table S1 for further details. Static theoretical structure factors were computed up to sinθ/λ = 1.30 Å−1 using the XFAC option in CRYSTAL14 (Flack, 1984; Le Page & Gabe, 1979). Multipole refinements against theoretical structure factors were performed using MoPro (Guillot et al., 2001; Jelsch et al., 2005) with a strategy analogous to that of the experimental structure factors. The scale factors and the x, y and z parameters were not refined. Basic of the are given in Table 2 and more information can be found in the supporting information (see text and Figs. S5–S12 and S14).
|
2.4. Hirshfeld surface and QTAIM analyses
Hirshfeld surface analyses (Spackman & McKinnon, 2002; McKinnon et al., 2007) were performed for experimental geometries using Crystal Explorer (Version 3.3; Wolff et al., 2012). Wavefunction calculations were done for each molecule separately, applying their formal charge, at the HF/6-31G** level (Hehre et al., 1972) using GAUSSIAN09 (Frisch et al., 2016).
A QTAIM (quantum theory of atoms in molecules; Bader, 1990) analysis was carried out on all electron-density models of CC, ACH and GDC crystals. For models in the multipole representation (experimental, UBDB and the one fitted to the theoretical structure factors), integrated charges were computed using the WinXPRO program (Stash & Tsirelson, 2002, 2007). For exact (not approximated by the multipolar model) theoretical crystal electron densities, integrated charges were calculated for experimental and optimized geometries using TOPOND14 (Gatti et al., 1994; Tsirelson, 2002).
2.5. Electrostatic potentials
Electrostatic potentials for single molecules, selected dimers and entire crystals were computed from charge-density models in multipole representation with the use of the XDPROP module of XD2016 (Volkov et al., 2016) and visualized by MoleCoolQt (Hübschle & Dittrich, 2011).
2.6. Intermolecular interaction energies and electrostatic contributions to them
To compute intermolecular interaction energies for isolated dimers with geometries as found in the studied crystals, the density functional theory-based symmetry-adapted perturbation theory (DFT-SAPT) (Jansen & Hesselmann, 2001; Williams & Chabalowski, 2001; Jeziorski et al., 1994) method was applied. Formal charges were assigned to particular molecules. The total intermolecular interaction energy in the DFT-SAPT method is given as the sum of the first- (E1) and second-order (E2) perturbation energy terms and the energy term, specifically electrostatic (Epol1), induction (Eind2) and dispersion (Edisp2) energy terms, together with exchange-repulsion (Eexch1, Eexch-ind2 and Eexch -disp2) terms:
The correction was applied in all cases to estimate the polarization effect beyond the second order. The DFT-SAPT calculations applied Kohn–Sham (KS) orbitals which were determined using the PBE0AC exchange-correlation potential (Heßelmann & Jansen, 2002; Gross et al., 1996). For neutral and cationic molecules, the asymptotic behaviour of the functional was corrected. The values of the highest occupied (HOMO) and ionization energies were calculated using the PBE0 functional. The anionic systems (chloride anions) were left without asymptotic correction (Lee & Burke, 2010). All DFT-SAPT and HOMO energy calculations were done in MOLPRO2012.1 (Werner et al., 2012) with the aug-cc-pVTZ Dunning basis set (Kendall et al., 1992; Dunning, 1989).
For the calculation of electrostatic, polarization, dispersion and repulsion contributions to the lattice energy, the semi-classical density sum (the PIXEL method; Gavezzotti, 2011) was used, which relies on a dimer approximation. Molecular electron densities of molecules bearing their formal charges were obtained using GAUSSIAN09 at the MP2/6-31** level (Hariharan & Pople, 1973). The electron densities were then analysed using the PIXELc module (Gavezzotti, 2003a,b,c,d) of the Coulomb–London–Pauli (CLP) program (Gavezzotti, 2011) which allows the calculation of lattice energies. The total intermolecular interaction energy was defined in PIXEL as follows
where Ees is the electrostatic interaction energy, Epol is the polarization energy, Edisp is the dispersion energy and Erep is the repulsion energy between interacting molecules. Due to the limitations of the program, i.e. a maximum of two molecules are allowed in the the crystal cohesive energy was only computed for CC.
Electrostatic contributions to intermolecular interaction energies and to crystal cohesive energies were also computed from experimental, UBDB and theoretical multipolar models of charge densities with the XDPROP module of the XD2016 package (Volkov et al., 2016). The exact electrostatic energy (Ees) was computed with the use of the EPMM method (Volkov, Koritsanszky & Coppens, 2004). The electrostatic energies from molecular monopole moments (point charges, ECoul) were obtained by extracting monopole–monopole contributions from calculations with the mMM option or, alternatively, computed directly from Coulomb's law, ECoul = q1q2/4π∊0r. To differentiate between results for dimers and results for the whole crystal, energies for the latter will be abbreviated as Ecoh,es and Ecoh,Coul.
It is to be noted that Epol1 from equation (1) refers to the exact electrostatic interaction energy and from now onwards will be abbreviated as Ees.
2.7. Transition-state search
For two selected dimers of the protonated nucleobases, a transition-state search was performed in GAUSSIAN16 (Frisch et al., 2016). All the geometry optimizations, transition-state searches and vibrational frequency analyses were performed with the hybrid B3LYP functional using the 6-311G(d,p) basis set. The B3LYP method has repeatedly been shown to yield results that are at least equal to MP2 calculations (Singleton et al., 1997). To describe dispersion interactions properly, the Becke–Johnson damping function approach DFT-D3(BJ) was used (Grimme et al., 2011).
3. Results and discussion
3.1. Packing and structural details
Crystal structures based on conventional ; Cochran, 1951; Kistenmacher & Shigematsu, 1974; Cunane & Taylor, 1993; Matković-Čalogović & Sanković, 1999; Mandel, 1977; Nieger, 2006). Here, we will repeat some of the previously discussed aspects of the structures and point to some new ones, in order to compare all three samples and give a structural basis for further discussion of the interactions in terms of electron density and energy.
analyses of all three samples have been reported previously (Broomhead, 1948Cytosinium chloride (CC) crystallizes in the monoclinic P21/n with one 1H,3H-cytosinium cation and one chloride anion in the Despite their presumably overall repulsive interaction due to their formal charges, the two cytosinium cations are arranged into a centrosymmetric trans sugar/sugar edge base pair linked by two N1—H1⋯O2 hydrogen bonds [see dimer AA1 in Fig. 2(a), Table 3 for geometric information, and Fig. S15 for nucleobase edge nomenclature conventions (Leontis & Westhof, 2001; Leontis et al., 2002)]. By the term `base pair' we mean dimers of bases in which the bases lie roughly in the same plane and there are at least two hydrogen bonds (classical and/or weak) between the bases. The second possible base pair to be formed by two cytosinium cations, in which two N3—H3⋯O2 hydrogen bonds could develop to create a trans Watson–Crick/Watson–Crick edge base pair, is not observed here in the CC structure. Instead, the Watson–Crick and Hoogsteen edges of the cytosinium cations are involved in direct interactions with the chloride anions.
|
The AA1 dimers form an infinite one-dimensional ribbon, held together by two chloride anions interacting with the Hoogsteen edges of the cytosinium cations through N7—H7B⋯Cl1 and C5—H5⋯Cl1 hydrogen bonds on one side of the chloride anions (dimers AB1) and through an N7—H7B⋯Cl1 bond on the other side of the anions (dimers AB2, Fig. 2a). These ribbons further stack together to form a layer, shown in Fig. 2(b). Inside a single stack, the ribbons are shifted by ca 4.3 Å along the ribbons and towards one of the edges by ca 1.2 Å, with a distance between two ribbons of 3.375 Å. The relative orientation of the ribbons in a stack is probably related to optimizing the distribution of molecules of opposite charge: the centre of a cytosinium cation in one ribbon lies almost exactly above a chloride anion in another ribbon, thus allowing for π⋯Cl–contacts (dimers AB3) (Schottel et al., 2008; Mooibroek et al., 2008; Patrick, 2014; Chifotides & Dunbar, 2013). In addition, repulsive cation–cation interactions between ribbons are, presumably, partially counterbalanced by stacking interactions between two cytosinium rings (dimers AA2).
Two adjacent layers of ribbon stacks interact with each other in such a way that each chloride anion is in close contact through N7—H7A⋯Cl1 and N3—H3⋯Cl1 hydrogen bonds with yet another cytosinium cation (dimer AB4, see Fig. 2c). The N3—H3⋯Cl1 hydrogen bond is the shortest among all of the N—H⋯Cl hydrogen bonds observed in the CC structure and is generally among the shortest N—H⋯Cl bonds found in the CSD (Steiner, 1998; Mascal, 1997). The cations of one stack, and the ribbons to which they belong, are tilted by ca 57° with respect to the cations (and their ribbons) of another stack. In addition to cation–anion interactions between the layers of the stacks, cytosinium cations from two adjacent layers are in direct contact with each other through C6—H6⋯O2 hydrogen bonds (dimer AA3).
Adeninium chloride hemihydrate (ACH) crystallizes in the P2/n, with the consisting of one adeninium cation, one chloride anion and one half of a water molecule located on a twofold rotation axis. The cation adopts the most common 1H,9H-adeninium tautomer (Marian et al., 2005). From the variety of bidentate hydrogen-bonded base pairs which could be formed by adeninium cations, only two are observed in the structure of ACH: a centrosymmetric trans Hoogsteen/Hoogsteen edge pair with two N10—H10A⋯N7 hydrogen bonds (dimer AA1) and a centrosymmetric trans sugar/sugar edge pair with two C2—H2⋯N3 hydrogen bonds (dimer AA2, Fig. 3a). Atom C8 from the Hoogsteen edge, atom N9 from the sugar edge and the whole Watson–Crick edge are exposed to chloride anions and water molecules.
Both adeninium base pairs are arranged together into infinite one-dimensional ribbons (Fig. 3a) where the ribbon edges are framed by chloride ions and water molecules. Each AA1 dimer of cations is embraced by two chloride anions through N10—H10B⋯Cl1 and C8—H8⋯Cl1 hydrogen bonds (dimers AB1 and AB2, respectively), while each AA2 dimer is surrounded by two further chloride ions through N9—H9⋯Cl1 and C2—H2⋯Cl1 hydrogen bonds (dimers AB3 and AB4, respectively). To satisfy the hydrogen donor remaining at the Watson–Crick edge of the adeninium cation, a water molecule is positioned to form an N1—H1⋯O1 hydrogen bond (dimer AW1). Additionally, the oxygen atom donates a proton to the aforementioned chloride anions to form an O1—H1A⋯Cl1 hydrogen bond (dimer BW1). According to the geometric data, the O1—H1A⋯Cl1 and N9—H9⋯Cl1 hydrogen bonds seem to be the strongest hydrogen bonds involving chloride ions in the ACH structure.
Similarly to the CC structure, ribbons in the ACH crystal stack together to form two-dimensional layers (Fig. 3b). Stacked ribbons are shifted along by almost an entire molecule (ca 2.95 Å) and very slightly (ca 0.5 Å) towards one of the edges. The distance between the ribbons is ca 3.19 Å. For ACH, the relative shift of the ribbons is probably associated with minimization of anion–anion and cation–cation repulsion. Chloride anions of adjacent ribbons are placed above and below the empty space between the chloride anions or above and below the water molecule of the central ribbon. The water molecule forms a second O1—H1A⋯Cl1 hydrogen bond (symmetry-related to BW1 mentioned above) with one of the chloride ions of an adjacent ribbon. For adeninium cations from adjacent ribbons, only close contacts between atoms on the sugar edges were observed (dimers AA3).
Again, as in the case of the CC structure, adjacent stacks of ribbons are connected to each other by chloride and water molecules, here a common component in both interacting stacks (Fig. 3b). Thus, no new types of dimer are formed, albeit each chloride ion involved in dimers AB1 and AB2 in one ribbon forms dimers AB3, AB4 and BW1 with a ribbon from an adjacent layer of stacks. The ribbons from two adjacent stacks are tilted with respect to each other by ca 83°. There is no direct contact between cations from neighbouring layers of stacks.
For guaninium dichloride (GDC), the Pnma consists of one half of a guaninium dication and two halves of chloride anions. All atoms of GDC lie in a special position on the crystallographic mirror plane. All heterocyclic nitrogen atoms are protonated in the guaninium dication.
in the orthorhombicNeither of the two possible bidentate hydrogen-bonded base pairs of guaninium dications is present in the structure of GDC. Contrary to the two previously discussed structures, there is no ribbon architecture observed here. Instead, two-dimensional layers of guaninium dications and chloride anions being in contact with each other through hydrogen bonds (Fig. 4a) are present. Each chloride anion in the layer is in direct hydrogen-bonding contact with three guaninium dications. For the Cl1 ion, the N11—H11B⋯Cl1 and N3—H3⋯Cl1 hydrogen bonds are present in dimer AB1, N1—H1⋯Cl1 in dimer AB2 and C8—H8⋯Cl1 in dimer AB3. For the Cl2 ion, the hydrogen bonds are N9—H9⋯Cl2 in dimer AB4, N7—H7⋯Cl2 in dimer AB5 and N11—H11A⋯Cl2 in dimer AB6. Among the N—H⋯Cl hydrogen bonds, N1—H1⋯Cl1 and N7—H7⋯Cl2 seem to be the strongest ones. The remaining close contacts within the single layer are between two guanine dications through N9—H9⋯O10 and C8—H8⋯O10 hydrogen bonds (dimer AA1). The layers of ions stack with each other in such a way that the Cl1 anions are almost exactly below and above atom C8 of the five-membered ring of the guaninium dication (dimers AB7), and the Cl2 anions are almost exactly below and above the six-membered ring (dimers AB8) (Fig. 4b). There are no close atom–atom contacts indicating stacking interactions between the guaninium dications. The distance between the layers is 3.242 Å.
To understand the architecture of these crystals further, it is essential to take a closer look at the overall molecular packing features. The Hirshfeld surface is a useful tool to investigate detailed aspects of molecular packing. First of all it allows identification of all molecules which are in direct contact with the molecule under consideration, regardless of whether they are close enough to be found using the van der Waals radii concept or not. After computing the Hirshfeld surfaces for each molecule separately in the studied structures, it appears that the molecules are surrounded by more close neighbours than have been discussed until now. For example, in the case of the chloride anion in CC (Fig. 5), six cytosinium cations, not four, contribute to its Hirshfeld surface. The dnorm function mapped on the Hirshfeld surface nicely shows several red spots, indicating very close intermolecular interactions. However, it is obvious from the fingerprint plots associated with each Hirshfeld surface (Fig. 5 and Fig. S16) that, besides atom–atom contacts which fall below the van der Waals limit and manifest on the plots as sharp spikes, there are more intermolecular atom–atom contacts present in these structures. Analysis of the percentage contribution of the observed contacts to the Hirshfeld surfaces (Fig. S17) confirms that the previously identified hydrogen bonds of D—H⋯Cl− and D—H⋯N/O types are very common in all the studied structures. Also, D—H⋯(C)π contacts are seen in all three structures. Not only previously identified (C)π⋯Cl− contacts, but also (N)π⋯Cl− contacts appear in the studied crystals, mainly in CC and GDC. In addition to the known CC and ACH (C)π⋯(C)π contacts, ones involving nitrogen atoms also appear in these two structures. What is qualitatively new with respect to atom–atom contacts identified on the basis of van der Waals radii (Table 3) is that the Hirshfeld surface analysis shows that H⋯H contacts between nucleobases are the most common in these crystals. This suggests their importance in stabilization. It should also be noted that, from a chloride anion point of view, the CC and GDC structures are more similar to each other than to ACH. In the latter, intermolecular interactions between chloride ions are facilitated almost exclusively by hydrogen bonds, whereas in the former, (C)π⋯Cl− interactions also seem to be important.
3.2. Molecular charges
The studied crystals are ionic in nature, so for them in particular electrostatic interactions are expected to play a major role in crystal stability. One of the most important factors when analysing intermolecular interactions is molecular charge. According to the experimental charge densities of the studied crystals, the charges on the nucleobase and chloride ions are not equal to the formal ones, i.e. the molecules are not fully ionized, as was already suggested for ACH by Cunane & Taylor (1993). Some degree of electron-density transfer from the chloride ions to the protonated nucleobase cations was observed. In all the studied structures, the chloride anions bear less negative charge, between −0.78 (4) e and −0.88 (5) e as computed directly from the Pval parameters of the multipolar model (Fig. 6). Accordingly, the nucleobase ions are less positively charged. Experimental charge-density integration over atomic basins calculated using QTAIM analysis leads to similar values of molecular charge. Moreover, topological charges computed directly from periodic calculations provide independent support to the finding that the chloride anions are less negative in these crystals; their charges are in the range of −0.70 e to −0.80 e according to theory. Values of charge departing from the formal −1 e for chloride anions in organic salts have recently been reported (Niranjana Devi et al., 2017; Wang et al., 2017; Nelyubina & Lyssenko, 2015; Nelyubina et al., 2015a,b, 2007).
Interestingly, the charge on the water molecule in ACH equals almost 0.0 e according to both experiment and theory.
In order to gain more evidence supporting our findings regarding charge transfer, we decided to compare experimental charge-density models with those in multipolar model representations simulated from periodic quantum mechanical calculations. On the basis of the periodic calculations, we computed static theoretical structure factors and fitted multipolar models to them. In addition, we checked how usage of the Cl0 neutral scattering radial functions instead of Cl−1 ionic functions influences the results.
First of all, it appears that multipolar models fitted in et al., 2011), as visible on residual maps (see the text and Figs. S5, S7, S9, S11 and S14 in the supporting information). Molecular topological charges computed from multipolar models differ by only ca 0.08 e from the target ones and are always too negative for chloride ions (and too positive for nucleobase cations).
to the theoretical charge densities reproduce the molecular topological charges very well, despite the known limitations of multipolar models in describing electron densities in the core region (FischerThe above is true only if Cl−1 ionic scattering radial functions are used. When Cl0 neutral scattering radial functions are used, the charge on the chloride ions in the multipolar model representation is substantially reduced, up to ca −0.42 e, and correspondingly the same happens for the molecular charges of the appropriate nucleobases. However, from R-factor statistics and residual-density analysis (Table 2, and Fig. S5–S12) it can be concluded that, in the case of theoretical charge densities, Cl−1 ionic scattering radial functions are more appropriate. They allow of multipolar models which provide charge densities closer to the target ones, yet they slightly overestimate the absolute values of the molecular charges.
There is no substantial difference in all the above findings when periodic calculations performed for experimental geometries are compared with calculations that include a geometry optimization step. Apparently, the differences in hydrogen-atom positions (the X—H distances differ the most among covalent bonds when these two geometries are compared) and the fact that only the latter calculations are based on systems in equilibrium do not matter here; the experimental geometry is accurate enough.
In every case, topological molecular charges computed from multipolar models are very similar (ca ±0.02 e) to those derived directly from Pval. This proves that multipolar models of particular molecules are well behaved and do not show too much delocalization compared with molecules of another kind. This is important for further intermolecular interaction energy calculations.
In the case of the experimental data, the use of Cl0 neutral scattering radial functions also leads to a reduction in the molecular charge from the formal one, but the reduction is not as large as in the case of the theoretical simulations. For all three structures, the chloride anions bear charges between −0.46 (4) e and −0.58 (5) e as computed directly from the Pval parameters of the multipolar model (Fig. 6). In contrast with the theoretical simulations, (R factors, statistical plots and residual-density analysis; Table 1 and Figs. S1–S4) are slightly better for models with neutral scattering factors than with ionic, but the difference is not large enough to disregard (because of, for example, possible functional or basis-set limitations; Medvedev et al., 2017; Cutini et al., 2016) conclusions resulting from the theoretical simulations or to oppose the principle of Ockham's razor: upon crystal formation, chloride anions most probably retain their ionic nature and they are still more ionic than neutral. Thus, we decided that multipolar models using Cl−1 ionic scattering radial functions are the best models of charge density which can be achieved from the fits to the collected experimental diffraction data of the studied crystals.
To summarize, in our opinion it is reasonably well proven by the data presented here that the chloride ions in the studied crystals are not fully ionized, and accordingly the same holds for the protonated nucleobases. The phenomenon might be general for all possible nucleobase chlorides and some other organic hydrochloride salts.
The theoretical simulations also helped to increase understanding of the possible causes of residual electron density observed around chloride anions in the case of the experimental data. On residual maps from the theoretical simulations there are no such pronounced residual features close to the chloride-ion positions that are visible on the experimental maps (Fig. S13). Thus, the assumption that experimental static electron densities are properly modelled by a multipolar model is justified. The features might be attributed either to anharmonic motion or to experimental errors like, for example, poor absorption corrections. Indeed, of the third-order anharmonic parameters in the case of ACH (the only data set with a resolution high enough to justify the refinement) clears the residual plots a lot (Figs. S1 and S13).
3.3. Electrostatic potential of molecules and crystals
The complementarity of molecular electrostatic potentials is believed to be one of the most important factors influencing how molecules interact with each other (Politzer & Murray, 2015).
In the case of ionic molecules, complementarity of molecular charge, i.e. positive charge complementing negative and vice versa, is an obvious observation. But in the crystals studied here, direct contacts between molecules with like charges, i.e. the protonated nucleobases, are observed. Interestingly for this type of molecular contact, some degree of complementarity can be observed (Fig. 7 top, and Figs. 2 and 3). Moreover, the electrostatic potential at the molecular surface for molecules in a crystal differs noticeably from that of isolated molecules (Fig. 7 top and bottom, and Fig. S17). The values of the electrostatic potentials of the former are, obviously, shifted towards less positive ones due to the fact that molecules in a crystal are less positive than isolated molecules with formal charges assigned to them (see the values of Vav.+ in Fig. 7). However, the electrostatic potential is redistributed within the nucleobase molecules. For example, in the case of the cytosinium cation in a crystal, its electrostatic potential around the N7 amino group is more positive than around the other N—H (and C—H) groups, including even the N3 protonation site. For an isolated cation, the difference is much smaller, if any. In the CC crystal, chloride ions cluster around the N7 amino group (Fig. 2), thus cation–anion interactions may in addition be enhanced by charge-density redistribution in a cationic molecule. Moreover, the redistribution may reduce cation–cation electrostatic repulsions: in dimers AA1 and AA3 the regions of molecules with the least positive electrostatic potentials cluster together, and in the case of dimer AA2 the region of the molecule having the most positive electrostatic potential (C4—N7H2) stacks with the least positive one (C2=O2). Similar effects can be observed for the adeninium cation in ACH, i.e. the charges redistribute in such a way that regions of the nucleobase interacting directly with the chloride have the highest positive potential and that potential is more uniformly spread over N—H and C—H fragments. These latter observations remarkably resemble the results of natural bond analysis of intermolecular N—H⋯Cl and C—H⋯Cl interactions in nicotinohydrazide dihydrochloride crystals (Kruszynski, 2011). Additionally, regions which are in direct contact with another cation have relatively less positive electrostatic potential (Figs. 3 and 7).
It is also enlightening to analyse the electrostatic potential itself (Fig. 8) and the gradient field trajectories of the electrostatic potential (Fig. 9) of single chloride or nucleobase ions, dimers and entire crystals. From the plots of a single cytosinium cation and of the AA1 dimer, it is clearly visible that it is only after formation of the dimer that the electrostatic potential gradient lines in the plane of the nucleobase surround the dimer uniformly. Hence, the dimer electrostatic potential field outside of the dimer surface resembles more the field of an ideal cation (a cation with a uniformly distributed charge inside it). This may explain why this type of dimer is formed in nucleobase chlorides.
3.4. Interaction energies for dimers
Following our analysis of the charge densities and electrostatic potentials of the studied molecules and crystals, we computed the energies of the intermolecular electrostatic interactions between whole molecules forming dimers as extracted from the crystal structures (Table 4). The list of analysed dimers includes not only dimers identified by atom–atom contacts being shorter than the sum of the van der Waals radii, but also from Hirshfeld surface analysis. In addition, interactions between the two closest chloride anions in the CC and ACH structures have been added to the analysis for the sake of comparison. To understand their role in intermolecular interactions in the studied dimers, we compared electrostatic energies with total interaction energies (Etot).
‡Values in parentheses in this column are energies computed from the multipole model fitted to theoretical structure factors computed from the optimized |
As one might expect for ionic structures, dimer interaction energies (Etot) are in general dominated by the net charges of the molecules, i.e. the energies are negative for dimers of molecules having opposite charges (AB-type dimers) and positive for dimers of molecules whose charges have the same sign (AA- or BB-type dimers). The absolute values of the energies are very high, much higher than those of neutral molecules (Jarzembska et al., 2012; Jarzembska, Goral et al., 2013; Jarzembska, Kamiński et al., 2013), reaching values of around ±200 kcal mol−1 (1 kcal mol−1 = 4.184 kJ mol−1) for dimers in the GDC structure containing dication(s). Decomposition of the total interaction energies within the perturbation theory scheme (Table S3) confirms that it is indeed the electrostatic component (Ees, the electrostatic interaction energy between two charge densities of isolated molecules) which contributes the most to the total energy, the other components being – with two exceptions (dimers AA1 in CC and ACH) – at least three times smaller in absolute value.
The Ees energies from the experimental crystal charge densities are systematically lower in absolute value than those from theoretical approaches relying on monomer and dimer gas-phase (not periodic) approximations (UBDB, DFT-SAPT). To understand the main source of this discrepancy, Ees energies can be contrasted with energies computed from Coulomb's law applied to two point molecular charges placed at the distances of the centres of molecular mass; see ECoul in Table 4. Although the point molecular charge is a very crude approximation of the true charge density, it adequately explains the overall trends observed for the Ees values in relation to the molecular charges and to distance. The ECoul energies computed from experimental molecular charges are shifted towards smaller absolute values with respect to the ECoul energies computed from formal charges. The shift is, in general, of the same magnitude as the shift in the Ees energies from experimental crystal charge densities with respect to the Ees energies from UBDB or DFT-SAPT. Thus, the fact that the molecules are not fully ionized in the crystal structures has a noticeable effect on the electrostatic energies computed for particular dimers.
A detailed analysis of why Ees deviates more from ECoul for some dimers than for others, whether it can be explained by atom–atom interactions, and whether topological analysis of intermolecular bonding lines will tell more in that respect will be presented in our next paper.
A closer look at the aforementioned shifts identifies that there is a secondary effect influencing the Ees energies that is unrelated to differences in molecular charge. The effect can be attributed to a charge-density redistribution within the molecules upon crystal formation, a redistribution already visible in the electrostatic potential at the molecular surfaces (see previous subsection). To quantify the effect we did the following simulation. We used the UBDB approach to build molecular electron densities of nucleobases rescaled to the molecular charges observed experimentally. In this way we obtained electron densities which have net charges as do those from the experiment, but they miss the effect of charge redistribution inside the molecule caused by the presence of neighbouring molecules. For the selected dimers we then computed the Ees energies from these charge densities (Table 4, values in parentheses in the UBDB column). The differences between these energies and those computed from the experimental charge densities quantify the effect of charge redistribution on the electrostatic interaction energies. The differences are small, up to around 5 kcal mol−1. Their values follow, to some extent, basic assumptions based on the electrostatic potential at a molecular surface, especially in the case of CC. However, not all dimers gain from charge redistribution. In the case of GDC, the role of charge redistribution is unclear, both from the electrostatic potential and from the energetic point of view. It has to be noted, however, that the multipole model is generally able to describe electronic polarization qualitatively, but has some difficulties with a quantitative description of it (Bąk et al., 2012).
The AA1 dimer in CC (and in ACH to some extent) deserves special attention. Despite the fact that it is a dimer of two cations, the negative electrostatic energy suggests that these two cations may form stable interactions. This is not only because the protonated nucleobases have charges lower than a formal +1 e in the crystalline state, but also because the charge redistributes within the cation. Another characteristic feature of the AA1 dimer in CC, and also in ACH, is the fact that the electrostatic interaction energies have the smallest absolute values when compared with the other dimers of AA type analysed here. In addition, the remaining components of the total interaction energy, like induction, exchange-repulsion and dispersion, contribute similar amounts of energy in absolute terms as the electrostatic component does, thus resembling the behaviour of dimers of neutral molecules (Czyżnikowska et al., 2010). Apparently, protonated nucleobases in highly charged environments may form (meta)stable base-pair dimers and their interactions may not differ significantly from those among neutral (non-protonated) nucleobases. Similarity in interactions between neutral dimers and dimers of like-charged molecules is not a new concept and has already been reported in the literature for dimers of hydrogen oxalate (Macchi et al., 2000), phosphates (Mata et al., 2012, 2013), oxalic acids (Mata et al., 2015) and (Alkorta et al., 2016), and for like-charged halogen-bonded complexes (Quiñonero et al., 2016).
To explore further the possibility of forming (meta)stable dimers with the AA1 dimers in CC and ACH, we have simulated interaction energy curves for the dissociating dimers as an isolated system in the gas phase using the DFT-SAPT approach (Fig. 10). The energies were computed for a series of dimers whose geometry was obtained by gradual extension or contraction of the hydrogen bond donor⋯acceptor distances, starting from the dimers' geometry as observed in the crystals. The interaction energy curves show local minima and local maxima. Thus, the simulation confirms that the studied dimers are indeed capable of forming metastable dimers in the gas phase and the geometry of these dimers very much resembles the geometry present in the crystalline state. The energy barrier is high enough to suggest that these dimers might be kinetically stable and experimentally observable, even in isolation in the gas phase and most probably in solution. The existence of the true local minima and transition states, and the depths of the local energy wells, were further confirmed by geometry optimization and transition-state search in the gas phase using the DFT approach (Fig. 10). Analogous observations made for other kinds of same-charge dimers, not for nucleobases, are present in the literature (Gamrad et al., 2015; Weinhold, 2017, 2018; Alkorta et al., 2016; Grimme & Djukic, 2011; Shokri et al., 2013; Mata et al., 2015).
To understand the influence of charge transfer on the stability of the studied pairs, we have corrected the energies obtained for Coulombic repulsion by subtracting from them the difference between the ECoul energy computed from the formal molecular charges and the ECoul energy computed from the partial molecular charges found for the nucleobases in the crystals (Fig. 10). The charge transfer lowers the absolute values of the interaction energy at the minima and maxima and seems to increase the dissociation barriers. Moreover, the charge transfer as observed in crystals brings the energy of the local minimum very close to the negative values and very little is needed to achieve truly attractive interactions. The effect resembles, to some extent, the influence of solvent discussed in the literature cited above. However, to our knowledge, we are the first to show that intermolecular charge transfer is responsible for the conversion of metastable dimers of same-charge molecules into stable ones.
All the nucleobase dimers analysed here (which will be shown in our next paper), and the AA1 dimers in CC and ACH in particular, show that ionic molecules in ionic crystals are not only passive hostages of Coulombic compression (Braga et al., 1998, 1999, 2002; D'Oria et al., 2012), but they actively respond to it through intermolecular charge transfer, intramolecular electronic polarization and local attractive interactions like hydrogen bonding. The possibility of forming protonated nucleobase dimers in solution prior to crystallization should not be excluded.
3.5. Cohesive energies for crystals
To sum up the analysis of the energies of intermolecular interactions in the crystals studied here, we computed the total crystal cohesive energies (Ecoh) as well as the electrostatic contributions to them (Ecoh,es) (Table 5). Crystal cohesive energies, analogous to dimers, result almost entirely from electrostatic interactions according to the PIXEL method. The other energy components, including polarization, are a couple of times smaller in their absolute values. As observed for the interaction energies for dimers, the electrostatic contribution seems to be a useful predictor for crystal cohesive energies in the ionic crystals included in this study. This could explain why the cohesive energy for GDC is much larger in absolute value than that of CC or ACH.
‡The remaining contributions to Ecoh are: Ecoh,pol (PIXEL) = −26.0 kcal mol−1, Ecoh,disp (PIXEL) = −35.0 kcal mol−1 and Ecoh,rep (PIXEL) = 55.0 kcal mol−1. |
As expected for stable ionic crystals, the negative electrostatic energies between the cations and anions in the studied crystals dominate over the positive electrostatic energies among the cations and anions. Interestingly, a large portion of the electrostatic contribution to the cohesive energy in GDC can be explained simply by invoking the interactions of point molecular charges placed in the crystal at centres of molecular mass (see Ecoh,es and Ecoh,Coul). In the case of CC, the shape of the cation and the charge distribution inside the molecule have almost as much influence on the overall electrostatic interactions in the crystal as molecular charges alone. For ACH, the molecular charges contribute even less than in the case of GDC and CC, because in this structure the charges are diluted by almost-neutral water molecules.
As in the case of electrostatic energies for dimers, the electrostatic contributions to the cohesive energies of the crystals are also very much influenced by the fact that the chlorides and nucleobases are not fully ionized (see values of Ecoh,Coul for the experimental and UBDB models).
The electrostatic contributions to the crystal cohesive energies (Ecoh,es) computed from the experimental densities are very similar to these from the UBDB simulation, in which the molecular charge is the same as in the experiment but charge redistribution within the molecules is missing (Table 4, values in parentheses in the UBDB column). This suggests that electronic polarization, after intermolecular charge transfer is separated out, does not sum to large values when the whole crystal is analysed.
The electrostatic contributions to crystal cohesive energies (Ecoh,es) computed from multipolar models fitted to simulated theoretical structure factors are very similar to those from experiment, except for GDC. Most probably, the discrepancies in the level of charge transfer as established from these two approaches are the main source of disparity in the Ecoh,es energies. It is uncertain which method (theory or experiment) is closer to the true value.
To sum up, two independent methods, i.e. experimental charge densities and periodic quantum mechanical calculations, gave values of crystal cohesive energies for the studied crystals at comparable levels. However, the relative difference of 5–20% might still not be very gratifying. There are many reasons for these discrepancies. Firstly, the Ecoh energies from the periodic calculations were computed with reference to fully ionized molecules, which is not the case for Ecoh,es from experiment. The energy of partial charge transfer from chloride to protonated nucleobase ions is missing from these calculations. Secondly, periodic calculations are likely to be affected by errors related to dispersion and/or correlation effects, whereas experimental results, apart from experimental errors, are limited by multipolar modelling. It is also uncertain how dispersion and exchange-repulsion effects present in the experimental crystal charge densities should be taken into account to compute their contribution to cohesive energies via summation over molecule–molecule interactions. Further studies are necessary to improve both methods and achieve better agreement.
4. Conclusions and outlook
The present work contributes to our understanding of intermolecular interactions among protonated nucleobases and chloride counter-ions surrounding them in the solid state. For single crystals of cytosinium chloride (CC), adeninium chloride hemihydrate (ACH) and guaninium dichloride (GDC), we collected high-resolution X-ray diffraction data to determine crystal charge densities experimentally. On the basis of these, we have extensively analysed the molecular charges, electrostatic potentials and electrostatic interaction energies for molecular dimers and whole crystals. These experimental results were further supported by extensive theoretical calculations, including UBDB charge-density modelling, dimer interaction energy computations based on perturbation theory (PIXEL, DFT-SAPT) and periodic quantum mechanical computations (CRYSTAL) based on the variational principle.
As expected for ionic crystals, our results confirm that crystal cohesive energies and the interaction energies of ionic dimers identified in the crystalline state are dominated by contributions from electrostatic interactions. Electrostatic energies (Ees) usually constitute 95% of the total interaction energy in these systems. However, the electrostatic interactions are not as strong as perturbation theory-based calculations suggest. This is because, according to the experimental data, the molecules included in this study are not fully ionized in the crystalline state. To our knowledge, this is the first experimental evidence, deeply explored in a quantitative manner from the point of view of interaction energy, that protonated nucleobases may bear a charge smaller than the formal one. Secondary effects of intramolecular charge redistribution, besides intermolecular charge transfer, are visible on molecular electrostatic potentials but seem to make minor contributions to the crystal cohesive energies.
Moreover, we show for the first time that singly protonated nucleobases are still capable of forming (meta)stable base pairs while surrounded by negative ions. This is because they interact strongly with each other through hydrogen bonding which overcomes cation–cation repulsion. The nature of this interaction does not differ much from those between neutral molecules. Chloride-to-nucleobase charge transfer and intramolecular charge-density polarization also contribute to this effect, which has not been observed experimentally before and which is discussed in detail in the context of same-charge dimers.
From a biological point of view, our finding is of great importance. Our study suggests that at least protonated adenine, as part of a nucleic acid, might be able to form stable homodimers of the type observed in the studied crystal, namely trans Hoogsteen/Hoogsteen edge base pairs. In the case of protonated cytosine, the particular base pair observed in the studied crystal cannot be formed by as it involves nitrogen atoms, which are bonded to a sugar moiety in But it is still not ruled out that other types of base pair, including heterodimers, could be formed between two protonated bases and have metastable characteristics. are highly charged polymers, bearing large negative charges due to their phosphate groups, often surrounded by many small cations and anions. Thus, they may impose a similar electric field on protonated bases as the chloride anions do in the studied crystals. What is important is that protonation of nucleobases is a well known mechanism by which these molecules change their properties towards specific recognition of particular partners. While analysing structures of and looking for the molecular basis of their specific interactions and functions, the possibility of attractive interactions between protonated nucleobases should not be excluded.
One may further speculate that base pairs which may exist in their protonated and neutral forms and which, while protonated, are metastable may have a unique ability to fine tune interactions between two nucleic acid chains. One can imagine that, in the case of the adenine–adenine trans Hoogsteen/Hoogsteen edge base pair (dimer AA1 in ACH), any changes in the hydrogen bonding at the N1 position will influence the partial charge of the base and thus regulate how stable the pair is. In addition, the N10⋯N7 distance might be very elongated before the pair will leave its metastable zone, even if stabilizing anions are removed from the close environment. Once the energy barrier of dissociation for the protonated pair has been overcome, a large amount of energy will be released (equivalent to the height of the maximum). This energy might be utilized for a conformational change of the interacting nucleic acid chains.
In fact, the adenine–adenine trans Hoogsteen/Hoogsteen edge base pair is the most common of all the types of adenine–adenine pairs present in the non-redundant set of RNA structures (https://ndbserver.rutgers.edu/ndbmodule/services/BPCatalog/bpCatalog.html): around 300 instances of this studied base pair were noted. For example, the studied pair is the only one present in the parallel double helix formed by two polyadenylic acid (poly A) chains (Gleghorn et al., 2016; Ke et al., 2009; Copp et al., 2017; Petrovic & Polavarapu, 2005, 2006; Kim et al., 2011; Safaee et al., 2013). Moreover, the existence of the pair in its protonated form in the double helix, exactly as in the studied crystals, with phosphate groups in the position of chlorides and with a water molecule in the same place, is well proven in acidic media. However, there is quite a debate in the literature as to whether the pair is protonated or not at neutral pH. It is worth stressing that poly A tails are very common in some types of RNA. Also, the studied base pair is very common in crystal structures of small molecules deposited in the CSD. There are ca 70 such structures containing the pair, among which more than half are built from two adenines protonated at the N1 position or N9 derivatives of adenines. Currently, we are doing more comprehensive statistics on the homopairs of nucleobases and their derivatives in their neutral and protonated forms present in the CSD and we will publish the results soon. The results presented here explain why this particular base pair is observed so frequently.
The question of whether the observed protonated base pairs are structure-determining remains open. Further studies of the dimerization of like-charged nucleobases in solution are also necessary. Our experimental results show that the presence of appropriate external electric fields substantially influences how charged nucleobases interact with each other, thus allowing for attractive interactions between molecules of the same net charge. Most probably, the same mechanism may occur in states other than crystalline ones. Nucleobase pairing is definitely RNA/DNA structure-determining, and it would be interesting to investigate whether it still holds true for pairings of protonated nucleobases.
Supporting information
Link https://doi.org/10.18150/repod.7313736
Raw diffraction frames
Link https://doi.org/10.18150/repod.0481200
Raw diffraction frames
Link https://doi.org/10.18150/repod.8020589
Raw diffraction frames
https://doi.org/10.1107/S2052252518006346/lc5097sup1.cif
contains datablocks CC, ACH, GDC. DOI:Structure factors: contains datablock CC. DOI: https://doi.org/10.1107/S2052252518006346/lc5097CCsup2.hkl
Structure factors: contains datablock ACH. DOI: https://doi.org/10.1107/S2052252518006346/lc5097ACHsup3.hkl
Structure factors: contains datablock GDC. DOI: https://doi.org/10.1107/S2052252518006346/lc5097GDCsup4.hkl
Supporting information file. DOI: https://doi.org/10.1107/S2052252518006346/lc5097CCsup5.cml
Supporting information file. DOI: https://doi.org/10.1107/S2052252518006346/lc5097ACHsup6.cml
Supporting information file. DOI: https://doi.org/10.1107/S2052252518006346/lc5097GDCsup7.cml
Experimental details, and additional figures and tables. DOI: https://doi.org/10.1107/S2052252518006346/lc5097sup8.pdf
Data collection: CrysAlis PRO, Agilent Technologies, Version 1.171.37.31 (release 14-01-2014 CrysAlis171 .NET) (compiled Jan 14 2014,18:38:05) for CC; CrysAlis PRO, Agilent Technologies, Version 1.171.36.32 (release 02-08-2013 CrysAlis171 .NET) (compiled Aug 2 2013,16:46:58) for ACH; APEX2 v2013.10-0 (Bruker AXS) for GDC. Cell
CrysAlis PRO, Agilent Technologies, Version 1.171.37.31 (release 14-01-2014 CrysAlis171 .NET) (compiled Jan 14 2014,18:38:05) for CC; CrysAlis PRO, Agilent Technologies, Version 1.171.36.32 (release 02-08-2013 CrysAlis171 .NET) (compiled Aug 2 2013,16:46:58) for ACH; APEX2 v2013.10-0 (Bruker AXS) for GDC. Data reduction: CrysAlis PRO, Agilent Technologies, Version 1.171.37.31 (release 14-01-2014 CrysAlis171 .NET) (compiled Jan 14 2014,18:38:05) for CC; CrysAlis PRO, Agilent Technologies, Version 1.171.36.32 (release 02-08-2013 CrysAlis171 .NET) (compiled Aug 2 2013,16:46:58) for ACH; SAINT V8.34A (Bruker AXS Inc., 2013) for GDC. For all structures, program(s) used to solve structure: SHELXS (Sheldrick, 2008); program(s) used to refine structure: MoPro (J. Appl. Cryst. 2005, 38, 38-54).C4H6N3O·Cl | F(000) = 304 |
Mr = 147.56 | Dx = 1.600 Mg m−3 |
Monoclinic, P21/n | Mo Kα radiation, λ = 0.71073 Å |
a = 8.1481 (1) Å | Cell parameters from 14655 reflections |
b = 6.8774 (1) Å | θ = 3.0–45.5° |
c = 10.9947 (1) Å | µ = 0.54 mm−1 |
β = 95.970 (1)° | T = 90 K |
V = 612.78 (1) Å3 | Block, transparent colorless |
Z = 4 | 0.17 × 0.13 × 0.05 mm |
SuperNova, Single source at offset, Eos diffractometer | 5144 independent reflections |
Radiation source: fine-focus sealed tube | 4590 reflections with > 2.0σ(I) |
Mirror monochromator | Rint = 0.042 |
Detector resolution: 16.0026 pixels mm-1 | θmax = 45.3°, θmin = 3.0° |
ω scans | h = −16→16 |
Absorption correction: multi-scan CrysAlisPro, Agilent Technologies, Version 1.171.37.31 (release 14-01-2014 CrysAlis171 .NET) (compiled Jan 14 2014,18:38:05) Empirical absorption correction using spherical harmonics, implemented in SCALE3 ABSPACK scaling algorithm. | k = −13→13 |
Tmin = 0.839, Tmax = 1.000 | l = −22→22 |
82389 measured reflections |
Refinement on F | Primary atom site location: structure-invariant direct methods |
Least-squares matrix: full | Secondary atom site location: difference Fourier map |
R[F2 > 2σ(F2)] = 0.016 | Hydrogen site location: difference Fourier map |
wR(F2) = 0.024 | Only H-atom coordinates refined |
S = 1.21 | Weighting scheme based on measured s.u.'s |
4590 reflections | (Δ/σ)max = −0.002 |
276 parameters | Δρmax = 0.27 e Å−3 |
6 restraints | Δρmin = −0.26 e Å−3 |
Refinement. Refinement of F1 against reflections. The threshold expression of F2 > 2sigma(F2) is used for calculating R-factors(gt) and is not relevant to the choice of reflections for refinement. R-factors based on F2 are statistically about twice as large as those based on F, and R-factors based on ALL data will be even larger. |
x | y | z | Uiso*/Ueq | ||
CL1 | 0.557537 (11) | 0.124370 (11) | 0.662831 (8) | 0.013327 (12) | |
O2 | −0.01306 (4) | 0.84121 (5) | 0.37276 (3) | 0.01386 (3) | |
N1 | 0.14888 (4) | 0.81248 (5) | 0.55465 (3) | 0.01102 (3) | |
N3 | 0.13971 (4) | 0.56699 (4) | 0.40919 (3) | 0.01013 (3) | |
C6 | 0.26193 (5) | 0.70920 (5) | 0.62854 (4) | 0.01156 (3) | |
C2 | 0.08729 (4) | 0.74804 (5) | 0.44161 (3) | 0.00996 (3) | |
N7 | 0.29604 (4) | 0.28850 (6) | 0.43622 (3) | 0.01401 (4) | |
C5 | 0.31723 (4) | 0.53264 (5) | 0.59493 (4) | 0.01146 (3) | |
C4 | 0.25282 (4) | 0.45840 (5) | 0.47871 (3) | 0.00996 (3) | |
H6 | 0.3068 (19) | 0.778 (2) | 0.712 (3) | 0.03081 | |
H7B | 0.3816 (9) | 0.2050 (12) | 0.4856 (8) | 0.02514 | |
H5 | 0.4092 (9) | 0.4506 (14) | 0.6519 (8) | 0.02559 | |
H3 | 0.1026 (15) | 0.5185 (17) | 0.3236 (6) | 0.02810 | |
H7A | 0.2507 (12) | 0.2467 (14) | 0.3508 (6) | 0.02918 | |
H1 | 0.105 (2) | 0.940 (5) | 0.5828 (17) | 0.02782 |
U11 | U22 | U33 | U12 | U13 | U23 | |
CL1 | 0.01730 (4) | 0.01272 (4) | 0.00966 (4) | 0.00343 (2) | −0.00001 (2) | 0.00259 (2) |
O2 | 0.01618 (11) | 0.01359 (10) | 0.01046 (12) | 0.00583 (10) | −0.00504 (10) | −0.00220 (9) |
N1 | 0.01336 (10) | 0.01056 (10) | 0.00830 (11) | 0.00315 (8) | −0.00287 (9) | −0.00192 (8) |
N3 | 0.01177 (10) | 0.00965 (9) | 0.00844 (10) | 0.00209 (8) | −0.00140 (9) | −0.00172 (8) |
C6 | 0.01347 (11) | 0.01184 (11) | 0.00848 (12) | 0.00254 (10) | −0.00307 (10) | −0.00130 (9) |
C2 | 0.01109 (10) | 0.01003 (10) | 0.00806 (11) | 0.00271 (9) | −0.00230 (10) | −0.00122 (9) |
N7 | 0.01654 (12) | 0.01083 (12) | 0.01447 (14) | 0.00417 (10) | 0.00071 (14) | −0.00190 (9) |
C5 | 0.01250 (11) | 0.01135 (11) | 0.00986 (12) | 0.00271 (9) | −0.00197 (10) | −0.00005 (9) |
C4 | 0.01071 (10) | 0.00924 (10) | 0.00972 (12) | 0.00169 (8) | 0.00009 (10) | −0.00045 (8) |
H6 | 0.03878 | 0.02996 | 0.02071 | 0.00803 | −0.01104 | −0.00967 |
H7B | 0.02773 | 0.02312 | 0.02300 | 0.00818 | −0.00487 | 0.00151 |
H5 | 0.02819 | 0.02543 | 0.02136 | 0.00853 | −0.00594 | 0.00206 |
H3 | 0.03629 | 0.02692 | 0.01871 | 0.00942 | −0.00843 | −0.00787 |
H7A | 0.03588 | 0.02809 | 0.02107 | 0.00837 | −0.00899 | −0.00776 |
H1 | 0.03460 | 0.02277 | 0.02410 | 0.01103 | −0.00636 | −0.00804 |
O2—C2 | 1.2343 (5) | C6—C5 | 1.3596 (5) |
N1—C2 | 1.3649 (5) | C6—H6 | 1.06 (3) |
N1—C6 | 1.3628 (5) | N7—C4 | 1.3198 (6) |
N1—H1 | 1.01 (4) | N7—H7A | 1.015 (5) |
N3—C2 | 1.3754 (5) | N7—H7B | 1.016 (5) |
N3—C4 | 1.3583 (5) | C5—C4 | 1.4246 (5) |
N3—H3 | 1.015 (5) | C5—H5 | 1.084 (5) |
C2—N1—C6 | 122.36 (3) | N1—C2—N3 | 115.80 (3) |
C2—N1—H1 | 117 (1) | C4—N7—H7A | 119.6 (5) |
C6—N1—H1 | 121 (1) | C4—N7—H7B | 120.3 (5) |
C2—N3—C4 | 124.24 (3) | H7A—N7—H7B | 119.8 (7) |
C2—N3—H3 | 117.9 (7) | C4—C5—C6 | 117.53 (3) |
C4—N3—H3 | 117.4 (7) | C4—C5—H5 | 120.0 (5) |
N1—C6—C5 | 121.74 (3) | C6—C5—H5 | 122.5 (5) |
N1—C6—H6 | 115.6 (10) | N3—C4—N7 | 118.55 (3) |
C5—C6—H6 | 123 (1) | N3—C4—C5 | 118.22 (3) |
O2—C2—N1 | 123.18 (3) | N7—C4—C5 | 123.22 (3) |
O2—C2—N3 | 120.98 (3) | ||
O2—C2—N1—C6 | −178.771 (14) | C6—C5—C4—N7 | −179.57 (6) |
O2—C2—N1—H1 | 3 (2) | C2—N1—C6—C5 | −1.65 (6) |
O2—C2—N3—C4 | 178.22 (4) | C2—N1—C6—H6 | 176.1 (6) |
O2—C2—N3—H3 | 6 (1) | C2—N3—C4—N7 | −178.39 (6) |
N1—C2—N3—C4 | −3.85 (5) | C2—N3—C4—C5 | 2.50 (5) |
N1—C2—N3—H3 | −175.9 (4) | N7—C4—N3—H3 | −6 (1) |
N1—C6—C5—C4 | 0.12 (5) | N7—C4—C5—H5 | 2 (1) |
N1—C6—C5—H5 | 178.56 (12) | C5—C4—N3—H3 | 174.6 (4) |
N3—C2—N1—C6 | 3.35 (6) | C5—C4—N7—H7A | −176.6 (3) |
N3—C2—N1—H1 | −175.3 (4) | C5—C4—N7—H7B | −2 (1) |
N3—C4—N7—H7A | 4 (1) | C5—C6—N1—H1 | 176.9 (3) |
N3—C4—N7—H7B | 179.08 (15) | C4—C5—C6—H6 | −177.5 (4) |
N3—C4—C5—C6 | −0.51 (6) | H6—C6—N1—H1 | −5 (2) |
N3—C4—C5—H5 | −178.99 (19) | H6—C6—C5—H5 | 1 (1) |
D—H···A | D—H | H···A | D···A | D—H···A |
N3—H3···CL1i | 1.01 (1) | 2.02 (1) | 3.0232 (4) | 169 (1) |
N7—H7A···CL1i | 1.01 (1) | 2.62 (1) | 3.4577 (4) | 140 (1) |
N7—H7B···CL1 | 1.02 (1) | 2.36 (1) | 3.3042 (4) | 154 (1) |
C5—H5···CL1 | 1.08 (1) | 2.55 (1) | 3.4600 (4) | 142 (1) |
C6—H6···O2ii | 1.06 (3) | 2.33 (2) | 3.1093 (6) | 129 (1) |
N1—H1···O2iii | 1.01 (4) | 1.77 (4) | 2.7782 (5) | 176 (1) |
N1—H1···C2iii | 1.01 (4) | 2.66 (4) | 3.5856 (5) | 154 (1) |
N7—H7B···CL1iv | 1.02 (1) | 2.87 (1) | 3.3079 (5) | 107 (1) |
Symmetry codes: (i) x−1/2, −y+1/2, z−1/2; (ii) x+1/2, −y+3/2, z+1/2; (iii) −x, −y+2, −z+1; (iv) −x+1, −y, −z+1. |
C5H6N5·Cl·0.5(H2O) | F(000) = 372 |
Mr = 180.60 | Dx = 1.622 Mg m−3 |
Monoclinic, P2/n | Mo Kα radiation, λ = 0.71073 Å |
a = 8.6936 (1) Å | Cell parameters from 35568 reflections |
b = 4.8189 (1) Å | θ = 2.5–65.7° |
c = 17.6971 (1) Å | µ = 0.46 mm−1 |
β = 93.526 (1)° | T = 90 K |
V = 739.99 (1) Å3 | Block, transparent colorless |
Z = 4 | 0.34 × 0.20 × 0.07 mm |
SuperNova, Single source at offset, Eos diffractometer | 13139 independent reflections |
Radiation source: SuperNova (Mo) X-ray Source | 11241 reflections with > 2.0σ(I) |
Mirror monochromator | Rint = 0.036 |
Detector resolution: 8.0013 pixels mm-1 | θmax = 66.0°, θmin = 2.3° |
ω scans | h = −22→22 |
Absorption correction: multi-scan CrysAlisPro, Agilent Technologies, Version 1.171.36.32 (release 02-08-2013 CrysAlis171 .NET) (compiled Aug 2 2013,16:46:58) Empirical absorption correction using spherical harmonics, implemented in SCALE3 ABSPACK scaling algorithm. | k = −11→12 |
Tmin = 0.916, Tmax = 1.000 | l = −45→45 |
101278 measured reflections |
Refinement on F | Primary atom site location: structure-invariant direct methods |
Least-squares matrix: full | Secondary atom site location: difference Fourier map |
R[F2 > 2σ(F2)] = 0.022 | Hydrogen site location: difference Fourier map |
wR(F2) = 0.028 | Only H-atom coordinates refined |
S = 1.21 | Weighting scheme based on measured s.u.'s |
11241 reflections | (Δ/σ)max = −0.003 |
358 parameters | Δρmax = 0.29 e Å−3 |
7 restraints | Δρmin = −0.24 e Å−3 |
Refinement. Refinement of F1 against reflections. The threshold expression of F2 > 2sigma(F2) is used for calculating R-factors(gt) and is not relevant to the choice of reflections for refinement. R-factors based on F2 are statistically about twice as large as those based on F, and R-factors based on ALL data will be even larger. |
x | y | z | Uiso*/Ueq | ||
CL1 | 0.493124 (7) | 0.818437 (16) | 0.274354 (3) | 0.013112 (12) | |
O1 | 0.25000 | 0.75833 (8) | 0.75000 | 0.01202 (5) | |
N7 | 0.12201 (2) | 0.19947 (4) | 0.433142 (10) | 0.00963 (3) | |
N1 | 0.27391 (2) | 0.59967 (4) | 0.597488 (10) | 0.00953 (3) | |
C5 | 0.19298 (2) | 0.37546 (4) | 0.485889 (10) | 0.00821 (3) | |
N3 | 0.38286 (2) | 0.75373 (5) | 0.485314 (11) | 0.01020 (3) | |
C4 | 0.29266 (2) | 0.55510 (5) | 0.451440 (10) | 0.00874 (3) | |
C6 | 0.18266 (2) | 0.39733 (4) | 0.564757 (10) | 0.00830 (3) | |
N9 | 0.28228 (2) | 0.48721 (5) | 0.376700 (10) | 0.01119 (3) | |
N10 | 0.09544 (2) | 0.24113 (5) | 0.605898 (10) | 0.01101 (3) | |
C2 | 0.36843 (2) | 0.76663 (5) | 0.558240 (12) | 0.01004 (3) | |
C8 | 0.17931 (3) | 0.27374 (5) | 0.368659 (11) | 0.01130 (4) | |
H2 | 0.4328 (12) | 0.9109 (18) | 0.5949 (7) | 0.02447 | |
H8 | 0.1493 (14) | 0.181 (2) | 0.3139 (4) | 0.02647 | |
H9 | 0.3414 (10) | 0.579 (2) | 0.3361 (5) | 0.02454 | |
H10A | 0.0272 (10) | 0.0980 (16) | 0.5790 (7) | 0.02379 | |
H1A | 0.1633 (8) | 0.8834 (17) | 0.7517 (5) | 0.02411 | |
H10B | 0.0911 (16) | 0.274 (3) | 0.6620 (2) | 0.02323 | |
H1 | 0.2672 (14) | 0.625 (3) | 0.65379 (18) | 0.02315 |
U11 | U22 | U33 | U12 | U13 | U23 | |
CL1 | 0.01463 (2) | 0.01601 (2) | 0.009058 (16) | −0.002779 (15) | 0.003707 (12) | 0.001535 (13) |
O1 | 0.01330 (8) | 0.01247 (10) | 0.01043 (7) | 0 | 0.00189 (5) | 0 |
N7 | 0.01117 (5) | 0.00999 (5) | 0.00783 (4) | −0.00302 (4) | 0.00139 (3) | −0.00062 (4) |
N1 | 0.01148 (5) | 0.00962 (5) | 0.00750 (4) | −0.00132 (4) | 0.00064 (3) | −0.00093 (3) |
C5 | 0.00939 (5) | 0.00855 (5) | 0.00677 (4) | −0.00152 (4) | 0.00124 (3) | −0.00004 (4) |
N3 | 0.01107 (5) | 0.01020 (6) | 0.00940 (5) | −0.00325 (4) | 0.00113 (4) | −0.00026 (4) |
C4 | 0.00954 (5) | 0.00932 (6) | 0.00747 (4) | −0.00206 (4) | 0.00129 (4) | 0.00002 (4) |
C6 | 0.00989 (5) | 0.00825 (5) | 0.00684 (4) | −0.00054 (4) | 0.00118 (3) | 0.00015 (4) |
N9 | 0.01304 (5) | 0.01338 (6) | 0.00739 (4) | −0.00447 (5) | 0.00262 (4) | −0.00017 (4) |
N10 | 0.01410 (5) | 0.01139 (6) | 0.00779 (4) | −0.00281 (4) | 0.00279 (4) | 0.00057 (4) |
C2 | 0.01109 (5) | 0.00952 (6) | 0.00945 (5) | −0.00198 (4) | 0.00028 (4) | −0.00107 (4) |
C8 | 0.01354 (6) | 0.01306 (7) | 0.00745 (5) | −0.00399 (5) | 0.00179 (4) | −0.00131 (4) |
H2 | 0.02848 | 0.02477 | 0.01984 | −0.01006 | −0.00097 | −0.00492 |
H8 | 0.03342 | 0.03087 | 0.01508 | −0.00889 | 0.00114 | −0.00643 |
H9 | 0.02900 | 0.02874 | 0.01652 | −0.00954 | 0.00668 | 0.00114 |
H10A | 0.02697 | 0.02301 | 0.02143 | −0.00961 | 0.00186 | −0.00321 |
H1A | 0.02540 | 0.02379 | 0.02323 | 0.00611 | 0.00233 | −0.00075 |
H10B | 0.02983 | 0.02651 | 0.01373 | −0.00425 | 0.00439 | −0.00109 |
H1 | 0.03036 | 0.02644 | 0.01290 | −0.00635 | 0.00343 | −0.00189 |
O1—H1A | 0.967 (3) | N3—C2 | 1.3058 (3) |
N7—C8 | 1.3223 (3) | C4—N9 | 1.3601 (3) |
N7—C5 | 1.3786 (3) | C6—N10 | 1.3189 (3) |
N1—C6 | 1.3634 (3) | N9—C8 | 1.3653 (3) |
N1—C2 | 1.3694 (3) | N9—H9 | 1.009 (3) |
N1—H1 | 1.009 (3) | N10—H10B | 1.009 (3) |
C5—C4 | 1.3918 (3) | N10—H10A | 1.009 (3) |
C5—C6 | 1.4081 (3) | C2—H2 | 1.083 (3) |
N3—C4 | 1.3537 (3) | C8—H8 | 1.083 (3) |
C8—N7—C5 | 103.954 (16) | N10—C6—C5 | 125.334 (15) |
C6—N1—C2 | 123.799 (17) | C4—N9—C8 | 106.647 (17) |
C6—N1—H1 | 116.0 (8) | C4—N9—H9 | 125.5 (6) |
C2—N1—H1 | 120.1 (7) | C8—N9—H9 | 127.9 (6) |
N7—C5—C4 | 110.539 (16) | C6—N10—H10B | 120.8 (8) |
N7—C5—C6 | 131.612 (15) | C6—N10—H10A | 118.1 (7) |
C4—C5—C6 | 117.841 (15) | H10B—N10—H10A | 121 (1) |
C4—N3—C2 | 112.525 (17) | N1—C2—N3 | 124.782 (17) |
N9—C4—N3 | 127.064 (17) | N1—C2—H2 | 112.2 (7) |
N9—C4—C5 | 105.693 (16) | N3—C2—H2 | 123.0 (7) |
N3—C4—C5 | 127.237 (17) | N9—C8—N7 | 113.166 (18) |
N1—C6—N10 | 120.851 (17) | N9—C8—H8 | 121.6 (7) |
N1—C6—C5 | 113.814 (16) | N7—C8—H8 | 125.2 (6) |
N7—C8—N9—C4 | 0.19 (3) | C5—C6—N10—H10A | −2 (1) |
N7—C8—N9—H9 | 180 (1) | N3—C4—N9—C8 | 179.13 (3) |
N7—C5—C4—N9 | −0.06 (3) | N3—C4—N9—H9 | 0 (1) |
N7—C5—C4—N3 | −179.26 (3) | N3—C4—C5—C6 | −0.19 (3) |
N7—C5—C6—N1 | 179.32 (3) | N3—C2—N1—C6 | 0.11 (3) |
N7—C5—C6—N10 | −0.33 (2) | N3—C2—N1—H1 | −178.6 (6) |
N1—C6—N10—H10B | 3 (2) | C4—N9—C8—H8 | 180 (1) |
N1—C6—N10—H10A | 178.1 (3) | C4—N3—C2—H2 | 179 (1) |
N1—C6—C5—C4 | 0.50 (3) | C4—C5—N7—C8 | 0.17 (3) |
N1—C2—N3—C4 | 0.23 (3) | C4—C5—C6—N10 | −179.16 (3) |
C5—N7—C8—N9 | −0.22 (3) | C6—N1—C2—H2 | −179 (1) |
C5—N7—C8—H8 | −180 (1) | C6—C5—N7—C8 | −178.73 (3) |
C5—C4—N9—C8 | −0.07 (3) | C6—C5—C4—N9 | 179.01 (3) |
C5—C4—N9—H9 | −179 (1) | N9—C4—N3—C2 | −179.22 (4) |
C5—C4—N3—C2 | −0.18 (3) | N10—C6—N1—C2 | 179.19 (3) |
C5—C6—N1—C2 | −0.48 (3) | N10—C6—N1—H1 | −2 (2) |
C5—C6—N1—H1 | 178.3 (3) | H2—C2—N1—H1 | 2 (2) |
C5—C6—N10—H10B | −177.7 (5) | H8—C8—N9—H9 | −1 (2) |
D—H···A | D—H | H···A | D···A | D—H···A |
N10—H10B···CL1i | 1.01 (1) | 2.26 (1) | 3.1760 (2) | 151 (1) |
O1—H1A···CL1ii | 0.97 (1) | 2.12 (1) | 3.0742 (3) | 170 (1) |
N9—H9···CL1 | 1.01 (1) | 2.11 (1) | 3.0978 (2) | 166 (1) |
N1—H1···O1 | 1.01 (1) | 1.83 (1) | 2.8251 (2) | 166 (1) |
N10—H10A···N7iii | 1.01 (1) | 1.94 (1) | 2.8980 (3) | 158 (1) |
N10—H10A···C8iii | 1.01 (1) | 2.74 (1) | 3.4924 (3) | 132 (1) |
N1—H1···O1iv | 1.01 (1) | 1.83 (1) | 2.8251 (2) | 166 (1) |
Symmetry codes: (i) x−1/2, −y+1, z+1/2; (ii) x−1/2, −y+2, z+1/2; (iii) −x, −y, −z+1; (iv) −x+1/2, y, −z+3/2. |
0.5(C5H7N5O)·Cl | Dx = 1.712 Mg m−3 |
Mr = 112.25 | Mo Kα radiation, λ = 0.71073 Å |
Orthorhombic, Pnma | Cell parameters from 1075 reflections |
a = 13.5939 (5) Å | θ = 2.6–41.1° |
b = 6.4841 (2) Å | µ = 0.72 mm−1 |
c = 9.8600 (4) Å | T = 100 K |
V = 869.10 (5) Å3 | Block, opaque colorless |
Z = 8 | 0.16 × 0.16 × 0.06 mm |
F(000) = 456.0 |
Bruker APEX II ULTRA diffractometer | 5605 independent reflections |
Radiation source: Enhance (Mo) X-ray Source | 4538 reflections with > 2.0σ(I) |
Mirror monochromator | Rint = 0.059 |
φ and ω scans | θmax = 54.0°, θmin = 2.6° |
Absorption correction: multi-scan R.H. Blessing, Acta Cryst. (1995), A51, 33-38 | h = −14→14 |
Tmin = 0.863, Tmax = 0.907 | k = −22→22 |
232874 measured reflections | l = −30→30 |
Refinement on F | Primary atom site location: structure-invariant direct methods |
Least-squares matrix: full | Secondary atom site location: difference Fourier map |
R[F2 > 2σ(F2)] = 0.018 | Hydrogen site location: difference Fourier map |
wR(F2) = 0.019 | Only H-atom coordinates refined |
S = 1.05 | Weighting scheme based on measured s.u.'s |
4482 reflections | (Δ/σ)max = 0.004 |
322 parameters | Δρmax = 0.44 e Å−3 |
7 restraints | Δρmin = −0.32 e Å−3 |
Refinement. Refinement of F1 against reflections. The threshold expression of F2 > 2sigma(F2) is used for calculating R-factors(gt) and is not relevant to the choice of reflections for refinement. R-factors based on F2 are statistically about twice as large as those based on F, and R-factors based on ALL data will be even larger. |
x | y | z | Uiso*/Ueq | ||
O10 | 0.64794 (4) | 0.75000 | 0.65461 (7) | 0.02046 (10) | |
CL1 | 0.536194 (6) | 0.75000 | 0.320999 (10) | 0.015803 (12) | |
N7 | 0.84807 (2) | 0.75000 | 0.79937 (4) | 0.01118 (4) | |
H7 | 0.7985 (5) | 0.75000 | 0.8764 (6) | 0.02535 | |
CL2 | 0.674142 (6) | 0.75000 | 0.998482 (8) | 0.011765 (10) | |
N11 | 0.81770 (2) | 0.75000 | 0.25365 (4) | 0.01534 (5) | |
H11A | 0.7536 (4) | 0.75000 | 0.2011 (8) | 0.02798 | |
H11B | 0.8843 (3) | 0.75000 | 0.2072 (9) | 0.02658 | |
N3 | 0.90940 (2) | 0.75000 | 0.45196 (4) | 0.01036 (4) | |
H3 | 0.9690 (4) | 0.75000 | 0.3902 (8) | 0.02473 | |
N9 | 0.98376 (2) | 0.75000 | 0.68086 (4) | 0.01062 (4) | |
H9 | 1.0559 (3) | 0.75000 | 0.6544 (13) | 0.02407 | |
N1 | 0.73677 (2) | 0.75000 | 0.45887 (4) | 0.01132 (4) | |
H1 | 0.6723 (4) | 0.75000 | 0.4073 (9) | 0.02649 | |
C4 | 0.90760 (2) | 0.75000 | 0.59068 (4) | 0.00946 (3) | |
C8 | 0.94524 (2) | 0.75000 | 0.80731 (4) | 0.01193 (4) | |
H8 | 0.9876 (6) | 0.75000 | 0.9003 (5) | 0.02654 | |
C2 | 0.82117 (2) | 0.75000 | 0.38686 (4) | 0.01051 (4) | |
C6 | 0.72806 (2) | 0.75000 | 0.60010 (4) | 0.01181 (4) | |
C5 | 0.82225 (3) | 0.75000 | 0.66424 (4) | 0.01031 (4) |
U11 | U22 | U33 | U12 | U13 | U23 | |
O10 | 0.00495 (14) | 0.0447 (4) | 0.0118 (2) | 0 | 0.00112 (13) | 0 |
CL1 | 0.00951 (3) | 0.02118 (4) | 0.01671 (4) | 0 | −0.00449 (3) | 0 |
N7 | 0.00618 (9) | 0.01790 (14) | 0.00947 (11) | 0 | 0.00022 (8) | 0 |
H7 | 0.01575 | 0.04202 | 0.01827 | 0 | 0.00220 | 0 |
CL2 | 0.00892 (3) | 0.01657 (4) | 0.00981 (3) | 0 | −0.00050 (2) | 0 |
N11 | 0.00934 (11) | 0.02753 (19) | 0.00914 (12) | 0 | −0.00030 (7) | 0 |
H11A | 0.01961 | 0.04275 | 0.02159 | 0 | −0.00577 | 0 |
H11B | 0.01895 | 0.03827 | 0.02253 | 0 | 0.00420 | 0 |
N3 | 0.00577 (9) | 0.01599 (13) | 0.00932 (11) | 0 | 0.00069 (8) | 0 |
H3 | 0.01550 | 0.03713 | 0.02157 | 0 | 0.00434 | 0 |
N9 | 0.00494 (8) | 0.01606 (13) | 0.01086 (11) | 0 | −0.00027 (9) | 0 |
H9 | 0.01097 | 0.03531 | 0.02594 | 0 | 0.00073 | 0 |
N1 | 0.00573 (9) | 0.01859 (14) | 0.00965 (12) | 0 | −0.00031 (9) | 0 |
H1 | 0.01445 | 0.04376 | 0.02126 | 0 | −0.00559 | 0 |
C4 | 0.00513 (8) | 0.01358 (11) | 0.00966 (10) | 0 | 0.00010 (8) | 0 |
C8 | 0.00659 (9) | 0.01884 (14) | 0.01036 (11) | 0 | −0.00086 (8) | 0 |
H8 | 0.01875 | 0.04105 | 0.01982 | 0 | −0.00762 | 0 |
C2 | 0.00656 (9) | 0.01560 (12) | 0.00936 (10) | 0 | 0.00003 (8) | 0 |
C6 | 0.00465 (8) | 0.02059 (14) | 0.01020 (11) | 0 | 0.00021 (8) | 0 |
C5 | 0.00525 (10) | 0.01649 (15) | 0.00919 (13) | 0 | 0.00022 (10) | 0 |
O10—C6 | 1.2146 (7) | N9—C4 | 1.3647 (5) |
N7—C8 | 1.3232 (5) | N9—C8 | 1.3522 (5) |
N7—C5 | 1.3778 (6) | N9—H9 | 1.0150 (11) |
N7—H7 | 1.0151 (11) | N1—C6 | 1.3976 (5) |
N11—C2 | 1.3143 (5) | N1—C2 | 1.3492 (5) |
N11—H11A | 1.0140 (11) | N1—H1 | 1.0140 (11) |
N11—H11B | 1.0141 (11) | C4—C5 | 1.3683 (5) |
N3—C2 | 1.3604 (5) | C8—H8 | 1.0830 (11) |
N3—C4 | 1.3680 (5) | C6—C5 | 1.4282 (5) |
N3—H3 | 1.0140 (11) | ||
C8—N7—C5 | 108.15 (3) | N9—C4—N3 | 129.63 (3) |
C8—N7—H7 | 128.2 (4) | N9—C4—C5 | 107.33 (3) |
C5—N7—H7 | 123.7 (4) | N3—C4—C5 | 123.04 (3) |
C2—N11—H11A | 122.8 (5) | N7—C8—N9 | 109.39 (3) |
C2—N11—H11B | 114.8 (5) | N7—C8—H8 | 125.5 (4) |
H11A—N11—H11B | 122.4 (7) | N9—C8—H8 | 125.1 (4) |
C2—N3—C4 | 117.13 (3) | N1—C2—N3 | 120.10 (3) |
C2—N3—H3 | 114.9 (5) | N1—C2—N11 | 119.70 (3) |
C4—N3—H3 | 127.9 (4) | N3—C2—N11 | 120.21 (3) |
C4—N9—C8 | 107.88 (3) | O10—C6—N1 | 121.13 (4) |
C4—N9—H9 | 124.5 (6) | O10—C6—C5 | 127.45 (4) |
C8—N9—H9 | 127.7 (6) | N1—C6—C5 | 111.42 (3) |
C6—N1—C2 | 126.61 (3) | N7—C5—C6 | 131.04 (3) |
C6—N1—H1 | 115.3 (5) | N7—C5—C4 | 107.25 (3) |
C2—N1—H1 | 118.1 (5) | C6—C5—C4 | 121.70 (3) |
O10—C6—N1—C2 | 180.00 (6) | N3—C2—N1—H1 | 180.00 (15) |
O10—C6—N1—H1 | 0 (1) | N3—C4—N9—C8 | 180.00 (5) |
O10—C6—C5—N7 | 0.00 (7) | N3—C4—N9—H9 | 0 (1) |
O10—C6—C5—C4 | 180.00 (6) | N3—C4—C5—C6 | 0.00 (5) |
N7—C8—N9—C4 | 0.00 (6) | H3—N3—C2—N1 | 180.0 (4) |
N7—C8—N9—H9 | 180 (1) | H3—N3—C4—N9 | 0.0 (9) |
N7—C5—C6—N1 | 180.00 (6) | H3—N3—C4—C5 | 180.0 (3) |
N7—C5—C4—N9 | 0.00 (6) | N9—C4—N3—C2 | 180.00 (6) |
N7—C5—C4—N3 | 180.00 (6) | N9—C4—C5—C6 | 180.00 (6) |
H7—N7—C8—N9 | 180.0 (8) | N9—C8—N7—C5 | 0.00 (5) |
H7—N7—C8—H8 | 0 (1) | H9—N9—C4—C5 | 180 (1) |
H7—N7—C5—C6 | 0.0 (9) | H9—N9—C8—H8 | 0 (1) |
H7—N7—C5—C4 | 180.0 (8) | N1—C6—C5—C4 | 0.00 (5) |
N11—C2—N1—C6 | 180.00 (5) | N1—C2—N3—C4 | 0.00 (5) |
N11—C2—N1—H1 | 0 (1) | H1—N1—C6—C5 | 180.0 (3) |
N11—C2—N3—C4 | 180.00 (4) | C4—N9—C8—H8 | 180.0 (9) |
N11—C2—N3—H3 | 0 (1) | C4—C5—N7—C8 | 0.00 (6) |
H11A—N11—C2—N1 | 0 (1) | C8—N7—C5—C6 | 180.00 (6) |
H11A—N11—C2—N3 | 180.00 (19) | C8—N9—C4—C5 | 0.00 (6) |
H11B—N11—C2—N1 | 180.0 (4) | H8—C8—N7—C5 | 180.0 (9) |
H11B—N11—C2—N3 | 0 (1) | C2—N1—C6—C5 | 0.00 (5) |
N3—C2—N1—C6 | 0.00 (5) | C2—N3—C4—C5 | 0.00 (5) |
D—H···A | D—H | H···A | D···A | D—H···A |
N11—H11A···CL2i | 1.01 (1) | 2.27 (1) | 3.1841 (4) | 149 (1) |
N7—H7···CL2 | 1.01 (1) | 2.08 (1) | 3.0732 (4) | 167 (1) |
N1—H1···CL1 | 1.01 (1) | 2.04 (1) | 3.0467 (4) | 175 (1) |
N11—H11B···CL1ii | 1.01 (1) | 2.08 (1) | 3.0600 (3) | 161 (1) |
N3—H3···CL1ii | 1.01 (1) | 2.27 (1) | 3.1960 (4) | 151 (1) |
N9—H9···O10iii | 1.01 (1) | 2.26 (1) | 2.7591 (7) | 109 (1) |
N9—H9···CL2iii | 1.01 (1) | 2.20 (1) | 3.1345 (4) | 152 (1) |
C8—H8···O10iii | 1.08 (1) | 2.25 (1) | 2.7809 (6) | 108 (1) |
N11—H11A···CL2iv | 1.01 (1) | 2.27 (1) | 3.1841 (4) | 149 (1) |
N7—H7···CL2v | 1.01 (1) | 2.08 (1) | 3.0732 (4) | 167 (1) |
N1—H1···CL1v | 1.01 (1) | 2.04 (1) | 3.0467 (4) | 175 (1) |
N11—H11B···CL1vi | 1.01 (1) | 2.08 (1) | 3.0600 (3) | 161 (1) |
N3—H3···CL1vi | 1.01 (1) | 2.27 (1) | 3.1960 (4) | 151 (1) |
N9—H9···O10vii | 1.01 (1) | 2.26 (1) | 2.7591 (7) | 109 (1) |
N9—H9···CL2vii | 1.01 (1) | 2.20 (1) | 3.1345 (4) | 152 (1) |
C8—H8···O10vii | 1.08 (1) | 2.25 (1) | 2.7809 (6) | 108 (1) |
Symmetry codes: (i) x, y, z−1; (ii) x+1/2, y, −z+1/2; (iii) x+1/2, y, −z+3/2; (iv) x, −y+3/2, z−1; (v) x, −y+3/2, z; (vi) x+1/2, −y+3/2, −z+1/2; (vii) x+1/2, −y+3/2, −z+3/2. |
Acknowledgements
The authors would like to thank Dr Sławomir Domagała, Dr Anna Makal and Dr Damian Trzybiński for valuable discussions and suggestions regarding data collection and reduction, and the anonymous referees for their valuable comments and suggestions which helped us to describe our results in a more precise and clear way.
Funding information
The following funding is acknowledged: Narodowe Centrum Nauki (grant No. NN204129138 to Paulina Maria Dominiak); Department of Chemistry, University of Warsaw (grant No. 120000-501/86-DSM-105000 to Prashant Kumar). The authors also thank the PL-Grid infrastructure (grant Nos. ELESAPT2016 and UBDB2017) for providing computation facilities.
References
Acharya, S., Barman, J., Cheruku, P., Chatterjee, S., Acharya, P., Isaksson, J. & Chattopadhyaya, J. (2004). J. Am. Chem. Soc. 126, 8674–8681. Web of Science CrossRef Google Scholar
Agilent (2013). CrysAlis PRO. Agilent Technologies, Yarnton, Oxfordshire, England. Google Scholar
Alberts, B., Johnson, A., Lewis, J., Raff, M., Roberts, K. & Walter, P. (2002). Molecular Biology of the Cell, 4th ed. New York: Garland Science. Google Scholar
Alkorta, I., Mata, I., Molins, E. & Espinosa, E. (2016). Chem. Eur. J. 22, 9226–9234. Web of Science CrossRef Google Scholar
Allen, F. H. & Bruno, I. J. (2010). Acta Cryst. B66, 380–386. Web of Science CSD CrossRef CAS IUCr Journals Google Scholar
Asensio, J. L., Lane, A. N., Dhesi, J., Bergqvist, S. & Brown, T. (1998). J. Mol. Biol. 275, 811–822. Web of Science CrossRef Google Scholar
Avery, O., Macleod, C. & McCarty, M. (1944). J. Exp. Med. 79, 137–158. CrossRef Google Scholar
Bader, R. F. W. (1990). Atoms in Molecules: A Quantum Theory. Oxford: Clarendon Press. Google Scholar
Bąk, J. M., Czyżnikowska, Ż. & Dominiak, P. M. (2012). Acta Cryst. A68, 705–714. Web of Science CrossRef IUCr Journals Google Scholar
Becke, A. D. (1988). Phys. Rev. A, 38, 3098–3100. CrossRef CAS PubMed Web of Science Google Scholar
Berger, I., Kang, C., Fredian, A., Ratliff, R., Moyzis, R. & Rich, A. (1995). Nat. Struct. Mol. Biol. 2, 416–425. CrossRef Web of Science Google Scholar
Bevilacqua, P. C., Brown, T. S., Nakano, S. & Yajima, R. (2004). Biopolymers, 73, 90–109. Web of Science CrossRef Google Scholar
Blessing, R. A. (1987). Crystallogr. Rev. 1, 3–58. CrossRef Google Scholar
Blessing, R. H. (1989). J. Appl. Cryst. 22, 396–397. CrossRef Web of Science IUCr Journals Google Scholar
Blessing, R. H. (1995). Acta Cryst. A51, 33–38. CrossRef CAS Web of Science IUCr Journals Google Scholar
Blessing, R. H. (1997). J. Appl. Cryst. 30, 421–426. CrossRef CAS Web of Science IUCr Journals Google Scholar
Bondi, A. (1964). J. Phys. Chem. 68, 441–451. CrossRef CAS Web of Science Google Scholar
Braga, D., Bazzi, C., Grepioni, F. & Novoa, J. J. (1999). New J. Chem. 23, 577–579. Web of Science CrossRef Google Scholar
Braga, D., Grepioni, F. & Novoa, J. (1998). Chem. Commun. pp. 1959–1960. Web of Science CrossRef Google Scholar
Braga, D., Maini, L. & Grepioni, F. (2002). Chem. Eur. J. 8, 1804–1812. Web of Science CSD CrossRef PubMed CAS Google Scholar
Broomhead, J. M. (1948). Acta Cryst. 1, 324–329. CrossRef IUCr Journals Web of Science Google Scholar
Bruker (2013). SAINT and SADABS. Bruker AXS Inc., Madison, Wisconsin, USA. Google Scholar
Cantor, C. R., Wollenzien, P. L. & Hearst, J. E. (1980). Nucleic Acids Res. 8, 1855–1872. CrossRef Web of Science Google Scholar
Cech, T. R. (1993a). Biochem. Soc. Trans. 21, 229–234. CrossRef Web of Science Google Scholar
Cech, T. R. (1993b). Gene, 135, 33–36. CrossRef Web of Science Google Scholar
Cech, T. R. (1993c). Nature, 365, 204–205. CrossRef Web of Science Google Scholar
Chawla, M., Sharma, P., Halder, S., Bhattacharyya, D. & Mitra, A. (2011). J. Phys. Chem. B, 115, 1469–1484. Web of Science CrossRef Google Scholar
Chifotides, H. T. & Dunbar, K. R. (2013). Acc. Chem. Res. 46, 894–906. Web of Science CrossRef CAS PubMed Google Scholar
Civalleri, B., Zicovich-Wilson, C. M., Valenzano, L. & Ugliengo, P. (2008). CrystEngComm, 10, 405–410. Web of Science CrossRef CAS Google Scholar
Cochran, W. (1951). Acta Cryst. 4, 81–92. CrossRef IUCr Journals Web of Science Google Scholar
Copp, W., Denisov, A. Y., Xie, X., Noronha, A. M., Liczner, C., Safaee, N., Wilds, C. J. & Gehring, K. (2017). Nucleic Acids Res. 45, 10321–10331. Web of Science CrossRef Google Scholar
Cunane, L. M. & Taylor, M. R. (1993). Acta Cryst. B49, 524–530. CSD CrossRef CAS Web of Science IUCr Journals Google Scholar
Cutini, M., Civalleri, B., Corno, M., Orlando, R., Brandenburg, J. G., Maschio, L. & Ugliengo, P. (2016). J. Chem. Theor. Comput. 12, 3340–3352. Web of Science CrossRef Google Scholar
Czyżnikowska, Z., Góra, R. W., Zaleśny, R., Lipkowski, P., Jarzembska, K. N., Dominiak, P. M. & Leszczynski, J. (2010). J. Phys. Chem. B, 114, 9629–9644. Web of Science PubMed Google Scholar
Dolomanov, O. V., Bourhis, L. J., Gildea, R. J., Howard, J. A. K. & Puschmann, H. (2009). J. Appl. Cryst. 42, 339–341. Web of Science CrossRef CAS IUCr Journals Google Scholar
D'Oria, E., Braga, D. & Novoa, J. J. (2012). CrystEngComm, 14, 792–798. CAS Google Scholar
Dovesi, R., Orlando, R., Erba, A., Zicovich-Wilson, C. M., Civalleri, B., Casassa, S., Maschio, L., Ferrabone, M., De La Pierre, M., D'Arco, P., Noël, Y., Causà, M., Rérat, M. & Kirtman, B. (2014). Int. J. Quantum Chem. 114, 1287–1317. Web of Science CrossRef CAS Google Scholar
Dovesi, R., Saunders, V. R., Roetti, C., Orlando, R., Zicovich-Wilson, C. M., Pascale, F., Civalleri, B., Doll, K., Harrison, N. M., Bush, I. J., DArco, P., Llunell, M., Causa, M. & Noel, Y. (2014). CRYSTAL14 User's Manual. University of Torino, Italy. Google Scholar
Dunning, J. (1989). J. Chem. Phys. 90, 1007–1023. CrossRef Web of Science Google Scholar
Farrugia, L. J. (2012). J. Appl. Cryst. 45, 849–854. Web of Science CrossRef CAS IUCr Journals Google Scholar
Fischer, A., Tiana, D., Scherer, W., Batke, K., Eickerling, G., Svendsen, H., Bindzus, N. & Iversen, B. B. (2011). J. Phys. Chem. A, 115, 13061–13071. Web of Science CrossRef CAS PubMed Google Scholar
Flack, H. D. (1984). J. Appl. Cryst. 17, 361–362. CrossRef Web of Science IUCr Journals Google Scholar
Frisch, M. J. et al. (2016). GAUSSIAN09. Gaussian Inc., Wallingford, Connecticut, USA. Google Scholar
Gamrad, W., Dreier, A., Goddard, R. & Pörschke, K. (2015). Angew. Chem. Int. Ed. 54, 4482–4487. Web of Science CrossRef Google Scholar
Gatti, C., Saunders, V. R. & Roetti, C. (1994). J. Chem. Phys. 101, 10686–10696. CrossRef CAS Web of Science Google Scholar
Gavezzotti, A. (2003a). J. Phys. Chem. B, 107, 2344–2353. Web of Science CrossRef CAS Google Scholar
Gavezzotti, A. (2003b). Z. Kristallogr. 220, 499–510. Google Scholar
Gavezzotti, A. (2003c). CrystEngComm, 5, 429–438. Web of Science CrossRef Google Scholar
Gavezzotti, A. (2003d). CrystEngComm, 5, 439–446. Web of Science CrossRef Google Scholar
Gavezzotti, A. (2011). New J. Chem. 35, 1360–1368. Web of Science CrossRef CAS Google Scholar
Gehring, K., Leroy, J. L. & Guéron, M. (1993). Nature, 363, 561–565. CrossRef CAS PubMed Web of Science Google Scholar
Gleghorn, M. L., Zhao, J., Turner, D. H. & Maquat, L. E. (2016). Nucleic Acids Res. 44, 8417–8424. Web of Science CrossRef Google Scholar
Goh, G. B., Knight, J. L. & Brooks, C. L. III (2013). J. Phys. Chem. Lett. 4, 760–766. Web of Science CrossRef Google Scholar
Grimme, S. (2006). J. Comput. Chem. 27, 1787–1799. Web of Science CrossRef PubMed CAS Google Scholar
Grimme, S. & Djukic, J. P. (2011). Inorg. Chem. 50, 2619–2628. Web of Science CrossRef Google Scholar
Grimme, S., Ehrlich, S. & Goerigk, L. (2011). J. Comput. Chem. 32, 1456–1465. Web of Science CrossRef CAS PubMed Google Scholar
Groom, C. R., Bruno, I. J., Lightfoot, M. P. & Ward, S. C. (2016). Acta Cryst. B72, 171–179. Web of Science CSD CrossRef IUCr Journals Google Scholar
Gross, E. K. U., Dobson, J. F. & Petersilka, M. (1996). Density Functional Theory of Time-Dependent Phenomena, pp. 81–172. Berlin, Heidelberg: Springer. Google Scholar
Guillot, B., Viry, L., Guillot, R., Lecomte, C. & Jelsch, C. (2001). J. Appl. Cryst. 34, 214–223. Web of Science CrossRef CAS IUCr Journals Google Scholar
Halder, A., Bhattacharya, S., Datta, A., Bhattacharyya, D. & Mitra, A. (2015). Phys. Chem. Chem. Phys. 17, 26249–26263. Web of Science CrossRef Google Scholar
Halder, A., Halder, S., Bhattacharyya, D. & Mitra, A. (2014). Phys. Chem. Chem. Phys. 16, 18383–18396. Web of Science CrossRef Google Scholar
Hansen, N. K. & Coppens, P. (1978). Acta Cryst. A34, 909–921. CrossRef CAS IUCr Journals Web of Science Google Scholar
Hariharan, P. C. & Pople, J. A. (1973). Theor. Chim. Acta, 28, 213–222. CrossRef CAS Web of Science Google Scholar
Hehre, W. J., Ditchfield, R. & Pople, J. A. (1972). J. Chem. Phys. 56, 2257–2261. CrossRef CAS Web of Science Google Scholar
Hesselmann, A. & Jansen, G. (2002). Chem. Phys. Lett. 357, 464–470. Google Scholar
Hübschle, C. B. & Dittrich, B. (2011). J. Appl. Cryst. 44, 238–240. Web of Science CrossRef IUCr Journals Google Scholar
Hübschle, C. B. & Luger, P. (2006). J. Appl. Cryst. 39, 901–904. Web of Science CrossRef IUCr Journals Google Scholar
Jansen, G. & Hesselmann, A. (2001). J. Phys. Chem. A, 105, 11156–11157. Web of Science CrossRef Google Scholar
Jarzembska, K. N. & Dominiak, P. M. (2012). Acta Cryst. A68, 139–147. Web of Science CrossRef CAS IUCr Journals Google Scholar
Jarzembska, K. N., Goral, A. M., Gajda, R. & Dominiak, P. M. (2013). Cryst. Growth Des. 13, 239–254. Web of Science CSD CrossRef CAS Google Scholar
Jarzembska, K. N., Kamiński, R., Wenger, E., Lecomte, C. & Dominiak, P. M. (2013). J. Phys. Chem. C, 117, 7764–7775. Web of Science CSD CrossRef CAS Google Scholar
Jarzembska, K. N., Kubsik, M., Kamiński, R., Woźniak, K. & Dominiak, P. M. (2012). Cryst. Growth Des. 12, 2508–2524. Web of Science CSD CrossRef CAS Google Scholar
Jelsch, C., Guillot, B., Lagoutte, A. & Lecomte, C. (2005). J. Appl. Cryst. 38, 38–54. Web of Science CrossRef IUCr Journals Google Scholar
Jeziorski, B., Moszynski, R. & Szalewicz, K. (1994). Chem. Rev. 94, 1887–1930. CrossRef CAS Web of Science Google Scholar
Johnson, C. K. (1969). Acta Cryst. A25, 187–194. CrossRef CAS IUCr Journals Web of Science Google Scholar
Ke, C., Loksztejn, A., Jiang, Y., Kim, M., Humeniuk, M., Rabbi, M. & Marszalek, P. E. (2009) Biophys. J. 96, 2918–2925. Web of Science CrossRef Google Scholar
Kendall, R. A., Dunning, J. T. H. & Harrison, R. J. (1992). J. Chem. Phys. 96, 6796–6806. CrossRef CAS Web of Science Google Scholar
Kim, S., Choi, J. & Majima, T. (2011). J. Phys. Chem. B, 115, 15399–15405. Web of Science CrossRef Google Scholar
Kissel, L., Zhou, B., Roy, S. C., Sen Gupta, S. K. & Pratt, R. H. (1995). Acta Cryst. A51, 271–288. CrossRef CAS Web of Science IUCr Journals Google Scholar
Kistenmacher, T. J. & Shigematsu, T. (1974). Acta Cryst. B30, 166–168. CSD CrossRef CAS IUCr Journals Web of Science Google Scholar
Kruszynski, R. (2011). Acta Cryst. C67, o52–o56. Web of Science CSD CrossRef IUCr Journals Google Scholar
Kuhs, W. F. (1983). Acta Cryst. A39, 148–158. CrossRef CAS Web of Science IUCr Journals Google Scholar
Lee, D. & Burke, K. (2010). Mol. Phys. 108, 2687–2701. Web of Science CrossRef Google Scholar
Lee, C., Yang, W. & Parr, R. G. (1988). Phys. Rev. B, 37, 785–789. CrossRef CAS Web of Science Google Scholar
Leontis, N. B., Stombaugh, J. & Westhof, E. (2002). Nucleic Acids Res. 30, 3497–3531. Web of Science CrossRef PubMed CAS Google Scholar
Leontis, N. B. & Westhof, E. (2001). RNA, 7, 499–512. Web of Science CrossRef PubMed CAS Google Scholar
Le Page, Y. & Gabe, E. J. (1979). J. Appl. Cryst. 12, 464–466. CrossRef CAS IUCr Journals Web of Science Google Scholar
Macchi, P., Iversen, B. B., Sironi, A., Chakoumakos, B. C. & Larsen, F. K. (2000). Angew. Chem. Int. Ed. 39, 2719–2722. Web of Science CrossRef CAS Google Scholar
Madsen, A. O. (2006). J. Appl. Cryst. 36, 757–758. Web of Science CrossRef IUCr Journals Google Scholar
Mandel, N. S. (1977). Acta Cryst. B33, 1079–1082. CSD CrossRef CAS IUCr Journals Web of Science Google Scholar
Marian, C., Nolting, D. & Weinkauf, R. (2005). Phys. Chem. Chem. Phys. 7, 3306–3316. Web of Science CrossRef Google Scholar
Mascal, M. (1997). J. Chem. Soc. Perkin Trans. 2, 12, 1999–2001. CrossRef Web of Science Google Scholar
Maschio, L., Civalleri, B., Ugliengo, P. & Gavezzotti, A. (2011). J. Phys. Chem. A, 115, 11179–11186. Web of Science CrossRef CAS PubMed Google Scholar
Mata, I., Alkorta, I., Molins, E. & Espinosa, E. (2012). ChemPhysChem, 13, 1421–1424. Web of Science CrossRef Google Scholar
Mata, I., Alkorta, I., Molins, E. & Espinosa, E. (2013). Chem. Phys. Lett. 555, 106–109. Web of Science CrossRef Google Scholar
Mata, I., Molins, E., Alkorta, I. & Espinosa, E. (2015). J. Phys. Chem. A, 119, 183–194. Web of Science CrossRef Google Scholar
Matković-Čalogović, D. & Sanković, K. (1999). Acta Cryst. C55, 467–469. Web of Science CSD CrossRef IUCr Journals Google Scholar
McKinnon, J. J., Jayatilaka, D. & Spackman, M. A. (2007). Chem. Commun. pp. 3814–3816. Web of Science CrossRef Google Scholar
Medvedev, M. G., Bushmarinov, I. S., Sun, J., Perdew, J. P. & Lyssenko, K. A. (2017). Science, 355, 49–52. Web of Science CrossRef Google Scholar
Meindl, K. & Henn, J. (2008). Acta Cryst. A64, 404–418. Web of Science CrossRef CAS IUCr Journals Google Scholar
Mooibroek, T. J., Black, C. A., Gamez, P. & Reedijk, J. (2008). Cryst. Growth Des. 8, 1082–1093. Web of Science CrossRef CAS Google Scholar
Muth, G. W., Ortoleva-Donnelly, L. & Strobel, S. A. (2000). Science, 289, 947–950. Web of Science CrossRef Google Scholar
Nelyubina, Y. V., Antipin, M. Yu. & Lyssenko, K. A. (2007). J. Phys. Chem. A, 111, 1091–1095. Web of Science CrossRef Google Scholar
Nelyubina, Y. V., Korlyukov, A. A. & Lyssenko, K. A. (2015a). Chem. Phys. Chem. 16, 676–681. Web of Science CrossRef Google Scholar
Nelyubina, Y. V., Korlyukov, A. A. & Lyssenko, K. A. (2015b). RSC Adv. 5, 75360–75373. Web of Science CrossRef Google Scholar
Nelyubina, Y. V. & Lyssenko, K. A. (2015). Chem. Eur. J. 21, 9733–9741. Web of Science CSD CrossRef CAS PubMed Google Scholar
Nieger, M. (2006). Private communication (refcode CYTOSC02). CCDC, Union Road, Cambridge, England. Google Scholar
Niranjana Devi, R., Jelsch, C., Israel, S., Aubert, E., Anzline, C. & Hosamani, A. A. (2017). Acta Cryst. B73, 10–22. Web of Science CSD CrossRef IUCr Journals Google Scholar
Patrick, G. (2014). Inorg. Chem. Front. 1, 35–43. Google Scholar
Pechlaner, M., Donghi, D., Zelenay, V. & Sigel, R. K. O. (2015). Angew. Chem. 127, 9823–9826. CrossRef Google Scholar
Perdew, J. P. (1986). Phys. Rev. B, 33, 8822–8824. CrossRef CAS Web of Science Google Scholar
Petrovic, A. G. & Polavarapu, P. L. (2005). J. Phys. Chem. B, 109, 23968–23705. Web of Science CrossRef Google Scholar
Petrovic, A. G. & Polavarapu, P. L. (2006). J. Phys. Chem. B, 110, 22826–22833. Web of Science CrossRef Google Scholar
Politzer, P. & Murray, J. S. (2015). Cryst. Growth Des. 15, 3767–3774. Web of Science CrossRef Google Scholar
Politzer, P., Murray, J. S. & Peralta-Inga, Z. (2001). Int. J. Quantum Chem. 85, 676–684. Web of Science CrossRef CAS Google Scholar
Quiñonero, D., Alkorta, I. & Elguero, J. (2016). Phys. Chem. Chem. Phys. 18, 27939–27950. Google Scholar
Safaee, N., Noronha, A. M., Rodionov, D., Kozlov, G., Wilds, C. J., Sheldrick, G. M. & Gehring, K. (2013). Angew. Chem. Int. Ed. 52, 10370–10373. Web of Science CrossRef CAS Google Scholar
Scheringer, C. (1985). Acta Cryst. A41, 79–81. CrossRef CAS Web of Science IUCr Journals Google Scholar
Schottel, B. L., Chifotides, H. T. & Dunbar, K. R. (2008). Chem. Soc. Rev. 37, 68–83. Web of Science CrossRef PubMed CAS Google Scholar
Sheldrick, G. M. (2008). Acta Cryst. A64, 112–122. Web of Science CrossRef CAS IUCr Journals Google Scholar
Sheldrick, G. M. (2015). Acta Cryst. C71, 3–8. Web of Science CrossRef IUCr Journals Google Scholar
Shokri, A., Ramezani, M., Fattahi, A. & Kass, S. R. (2013). J. Phys. Chem. A, 117, 9252–9258. Web of Science CrossRef Google Scholar
Siegfried, N. A., O'Hare, B. & Bevilacqua, P. C. (2010). Biochemistry, 49, 3225–3236. Web of Science CrossRef Google Scholar
Singleton, D. A., Merrigan, S. R., Liu, J. & Houk, K. N. (1997). J. Am. Chem. Soc. 119, 3385–3386. CrossRef Web of Science Google Scholar
Spackman, M. A. & McKinnon, J. J. (2002). CrystEngComm, 4, 378–392. Web of Science CrossRef CAS Google Scholar
Spek, A. L. (2009). Acta Cryst. D65, 148–155. Web of Science CrossRef CAS IUCr Journals Google Scholar
Stash, A. & Tsirelson, V. (2002). J. Appl. Cryst. 35, 371–373. Web of Science CrossRef CAS IUCr Journals Google Scholar
Stash, A. I. & Tsirelson, V. (2007). Crystallogr. Rep. 50, 177–184. Web of Science CrossRef Google Scholar
Steiner, T. (1998). Acta Cryst. B54, 456–463. Web of Science CrossRef CAS IUCr Journals Google Scholar
Stewart, R. F., Bentley, J. & Goodman, B. (1975). J. Chem. Phys. 63, 3786–3793. CrossRef CAS Web of Science Google Scholar
Su, Z. & Coppens, P. (1998). Acta Cryst. A54, 646–652. Web of Science CrossRef CAS IUCr Journals Google Scholar
Tang, C. L., Alexov, E., Pyle, A. M. & Honig, B. (2007). J. Mol. Biol. 366, 1475–1496. Web of Science CrossRef Google Scholar
Tsirelson, V. G. (2002). Acta Cryst. B58, 632–639. Web of Science CrossRef CAS IUCr Journals Google Scholar
Volkov, A., Koritsanszky, T. & Coppens, P. (2004). Chem. Phys. Lett. 391, 170–175. Web of Science CrossRef CAS Google Scholar
Volkov, A., Li, X., Koritsanszky, T. & Coppens, P. (2004). J. Phys. Chem. A, 108, 4283–4300. Web of Science CrossRef CAS Google Scholar
Volkov, A., Macchi, P., Farrugia, L. J., Gatti, C., Mallinson, P., Richter, T. & Koritsanszky, T. (2016). XD2016. University at Buffalo, State University of New York, USA, University of Milan, Italy, University of Glasgow, UK, CNRISTM, Milan, Italy, and Middle Tennessee State University, USA. Google Scholar
Wang, A., Wang, R., Kalf, I., Dreier, A., Lehmann, C. V. & Englert, U. (2017). Cryst. Growth Des. 17, 5, 2357–2364. Google Scholar
Ward, W. L., Plakos, K. & DeRose, V. (2014). Chem. Rev. 114, 4318–4342. Web of Science CrossRef Google Scholar
Watson, J. D. & Crick, F. H. C. (1953a). Nature, 171, 737–738. CrossRef Web of Science Google Scholar
Watson, J. D. & Crick, F. H. C. (1953b). Nature, 171, 964–967. CrossRef Web of Science Google Scholar
Weinhold, F. (2017). Angew. Chem. 129, 14769–14773. CrossRef Google Scholar
Weinhold, F. (2018). Inorg. Chem. 57, 2035–2044. Web of Science CrossRef Google Scholar
Werner, H.-J., Knowles, P. J., Knizia, G., Manby, F. R. & Schütz, M. (2012). WIREs Comput. Mol. Sci. 2, 242–253. Web of Science CrossRef Google Scholar
Wilcox, J. L., Ahluwalia, A. K. & Bevilacqua, P. C. (2011). Acc. Chem. Res. 44, 1270–1279. Web of Science CrossRef Google Scholar
Wilcox, J. L. & Bevilacqua, P. C. (2013). J. Am. Chem. Soc. 135, 7390–7393. Web of Science CrossRef Google Scholar
Williams, H. L. & Chabalowski, C. F. (2001). J. Phys. Chem. A, 105, 646–659. Web of Science CrossRef Google Scholar
Wolff, S., Grimwood, D., McKinnon, J., Turner, M., Jayatilaka, D. & Spackman, M. (2012). CrystalExplorer. University of Western Australia. Google Scholar
Wolter, A. C., Weickhmann, A. K., Nasiri, A. H., Hantke, K., Ohlenschläger, O., Wunderlich, C. H., Kreutz, C., Duchardt-Ferner, E. & Wöhnert, J. (2017). Angew. Chem. Int. Ed. 56, 401–404. Web of Science CrossRef Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.