Comparison of the π-stacking properties of purine versus pyrimidine residues. Some generalizations regarding selectivity

Interaction between a fluoroquinolone derivative KG022 and RNAs: Effect of base pairs 3′ adjacent to the bulged residues

RNA-targeted small molecules are a promising modality in drug discovery. Recently, we found that a fluoroquinolone derivative, KG022, can bind to RNAs with bulged C or G. To clarify the RNA specificity of KG022, we analyzed the effect of the base pair located at the 3′side of the bulged residue. It was found that KG022 prefers G-C and A-U base pairs at the 3′side. Solution structures of the complexes of KG022 with the four RNA molecules with bulged C or G and G-C or A-U base pairs at the 3′side of the bulged residue were determined to find that the fluoroquinolone moiety is located between two purine bases, and this may be the mechanism of the specificity. This work provides an important example of the specificity of RNA-targeted small molecules.

Recognition of single-stranded nucleic acids by small-molecule splicing modulators

Risdiplam is the first approved small-molecule splicing modulator for the treatment of spinal muscular atrophy (SMA). Previous studies demonstrated that risdiplam analogues have two separate binding sites in exon 7 of the SMN2 pre-mRNA: (i) the 5'-splice site and (ii) an upstream purine (GA)-rich binding site. Importantly, the sequence of this GA-rich binding site significantly enhanced the potency of risdiplam analogues. In this report, we unambiguously determined that a known risdiplam analogue, SMN-C2, binds to single-stranded GA-rich RNA in a sequence-specific manner. The minimum required binding sequence for SMN-C2 was identified as GAAGGAAGG. We performed all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method, which captured spontaneous binding of a risdiplam analogue to the target nucleic acids. We uncovered, for the first time, a ligand-binding pocket formed by two sequential GAAG loop-like structures. The simulation findings were highly consistent with experimental data obtained from saturation transfer difference (STD) NMR and structure-affinity-relationship studies of the risdiplam analogues. Together, these studies illuminate us to understand the molecular basis of single-stranded purine-rich RNA recognition by small-molecule splicing modulators with an unprecedented binding mode.

Simulated NIR spectra as sensitive markers of the structure and interactions in nucleobases

Near-infrared (near-IR; NIR) spectroscopy is continuously advancing in biophysical and biochemical fields of investigation. For instance, recent progresses in NIR hyperspectral imaging of biological systems may be noted. However, interpretation of NIR bands for biological samples is difficult and creates a considerable barrier in exploring the full potential of NIR spectroscopy in bioscience. For this reason, we carried out a systematic study of NIR spectra of adenine, cytosine, guanine, and thymine in polycrystalline state. Interpretation of NIR spectra of these nucleobases was supported by anharmonic vibrational analysis using Deperturbed Vibrational Second-Order Perturbation Theory (DVPT2). A number of molecular models of nucleobases was applied to study the effect of the inter-molecular interactions on the NIR spectra. The accuracy of simulated NIR spectra appears to depend on the intra-layer interactions; in contrast, the inter-layer interactions are less influential. The best results were achieved by combining the simulated spectra of monomers and dimers. It is of particular note that in-plane deformation bands are far more populated than out-of-plane ones and the importance of ring modes is relatively small. This trend is in contrast to that observed in mid-IR region. As shown, the local, short-range chemical neighborhood of nucleobase molecules influence their NIR spectra more considerably. This suggests that NIR spectra are more sensitive probe of the nucleobase pairing than mid-IR ones. The obtained results allow, for the first time, to construct a frequency correlation table for NIR spectra of purines and pyrimidines.

Characterization of the one-electron oxidized Cu(II)-salen complexes with a side chain aromatic ring: the effect of the indole ring on the Cu(II)-phenoxyl radical species

To gain insights into the role of the proximal indole ring in the redox-active metal center as seen in galactose oxidase, we prepared the Cu(II)-salen-type complexes having a pendent indol-3-ylmethyl (1), methyl (2) or benzyl (3) group substituted on the ethylenediamine moiety and investigated the structures and redox properties by various physicochemical methods and theoretical calculations. Neutral complexes 1, 2, and 3 showed no significant difference in the UV-Vis-NIR and EPR spectra. One-electron oxidation of 1, 2, and 3 by addition of 1 equiv. of thianthrenyl radical gave [1]SbCl ₆ , [2]SbCl ₆ , and [3]SbCl ₆ , respectively, which could be assigned to relatively localized phenoxyl radical species. The cyclic and differential pulse voltammograms of [1]SbCl ₆ showed two redox waves with a large separation between the first and second redox potentials compared with the separations observed for [2]SbCl ₆ and [3]SbCl ₆ . This suggests that [1]SbCl ₆ is more stabilized than [2]SbCl ₆ and [3]SbCl ₆ . The NIR band of [1]SbCl ₆ showed a larger blue shift than that of [2]SbCl ₆ and [3]SbCl ₆ . The EPR spectrum of [2]SbCl ₆ exhibited an intense signal at the g value of 2 due to partial disproportionation to form the EPR active two-electron oxidized complex [2] ²⁺ , while the EPR intensity of [1]SbCl ₆ was much weaker than that of [2]SbCl ₆ . These results indicate that the pendent indole moiety stabilizes the Cu(II)-phenoxyl radical in [1]SbCl ₆ most probably by stacking with the phenoxyl moiety, which is further supported by DFT calculations.

Differences in Integron Cassette Excision Dynamics Shape a Trade-Off between Evolvability and Genetic Capacitance

Integrons ensure a rapid and "on demand" response to environmental stresses driving bacterial adaptation. They are able to capture, store, and reorder functional gene cassettes due to site-specific recombination catalyzed by their integrase. Integrons can be either sedentary and chromosomally located or mobile when they are associated with transposons and plasmids. They are respectively called sedentary chromosomal integrons (SCIs) and mobile integrons (MIs). MIs are key players in the dissemination of antibiotic resistance genes. Here, we used in silico and in vivo approaches to study cassette excision dynamics in MIs and SCIs. We show that the orientation of cassette arrays relative to replication influences attC site folding and cassette excision by placing the recombinogenic strands of attC sites on either the leading or lagging strand template. We also demonstrate that stability of attC sites and their propensity to form recombinogenic structures also regulate cassette excision. We observe that cassette excision dynamics driven by these factors differ between MIs and SCIs. Cassettes with high excision rates are more commonly found on MIs, which favors their dissemination relative to SCIs. This is especially true for SCIs carried in the Vibrio genus, where maintenance of large cassette arrays and vertical transmission are crucial to serve as a reservoir of adaptive functions. These results expand the repertoire of known processes regulating integron recombination that were previously established and demonstrate that, in terms of cassette dynamics, a subtle trade-off between evolvability and genetic capacitance has been established in bacteria.IMPORTANCE The integron system confers upon bacteria a rapid adaptation capability in changing environments. Specifically, integrons are involved in the continuous emergence of bacteria resistant to almost all antibiotic treatments. The international situation is critical, and in 2050, the annual number of deaths caused by multiresistant bacteria could reach 10 million, exceeding the incidence of deaths related to cancer. It is crucial to increase our understanding of antibiotic resistance dissemination and therefore integron recombination dynamics to find new approaches to cope with the worldwide problem of multiresistance. Here, we studied the dynamics of recombination and dissemination of gene encoding cassettes carried on integrons. By combining in silico and in vivo analyses, we show that cassette excision is highly regulated by replication and by the intrinsic properties of cassette recombination sites. We also demonstrated differences in the dynamics of cassette recombination between mobile and sedentary chromosomal integrons (MIs and SCIs). For MIs, a high cassette recombination rate is favored and timed to conditions when generating diversity (upon which selection can act) allows for a rapid response to environmental conditions and stresses. In contrast, for SCIs, cassette excisions are less frequent, limiting cassette loss and ensuring a large pool of cassettes. We therefore confirm a role of SCIs as reservoirs of adaptive functions and demonstrate that the remarkable adaptive success of integron recombination system is due to its intricate regulation.

Efficiency of integron cassette insertion in correct orientation is ensured by the interplay of the three unpaired features of attC recombination sites

The integron is a bacterial recombination system that allows acquisition, stockpiling and expression of cassettes carrying protein-coding sequences, and is responsible for the emergence and rise of multiresistance in Gram-negative bacteria. The functionality of this system depends on the insertion of promoterless cassettes in correct orientation, allowing their expression from the promoter located upstream of the cassette array. Correct orientation is ensured by strand selectivity of integron integrases for the bottom strand of cassette recombination sites (attC), recombined in form of folded single-stranded hairpins. Here, we investigated the basis of such strand selectivity by comparing recombination of wild-type and mutated attC sites with different lengths, sequences and structures. We show that all three unpaired structural features that distinguish the bottom and top strands contribute to strand selectivity. The localization of Extra-Helical Bases (EHBs) directly favors integrase binding to the bottom strand. The Unpaired Central Spacer (UCS) and the Variable Terminal Structure (VTS) influence strand selectivity indirectly, probably through the stabilization of the bottom strand and the resulting synapse due to the nucleotide skew between the two strands. These results underscore the importance of the single-stranded nature of the attC site that allows such tight control over integron cassette orientation.

Understanding Stacking Interactions between an Aromatic Ring and Nucleobases in Aqueous Solution: Experimental and Theoretical Study

Stacking interactions between aromatic compounds and nucleobases are crucial in recognition of nucleotides and nucleic acids, but a comprehensive understanding of the strength and selectivity of these interactions in aqueous solution has been elusive. To this end, model complexes have been designed and analyzed by experiment and theory. For the first time, stacking free energies between five nucleobases and anthracene were determined experimentally from thermodynamic double mutant cycles. Three different experimental methods were proposed and evaluated. The dye prefers to bind nucleobases in the order (kcal/mol): G (1.3) > T (0.9) > U (0.8) > C (0.5) > A (0.3). The respective trend of interaction free energies extracted from DFT calculations correlates to that obtained experimentally. Analysis of the data suggests that stacking interactions dominate over hydrophobic effects in an aqueous solution and can be predicted with DFT calculations.

Solution properties of metal ion complexes formed with the antiviral and cytostatic nucleotide analogue 9-[2-(phosphonomethoxy)ethyl]-2-amino-6-dimethylaminopurine (PME2A6DMAP)

The acidity constants of protonated 9-[2-(phosphonomethoxy)ethyl]-2-amino-6-dimethylaminopurine (H3(PME2A6DMAP)+) are considered, and the stability constants of the M(H;PME2A6DMAP)+ and M(PME2A6DMAP) complexes (M2+ = Mg2+, Ca2+, Sr2+, Ba2+, Mn2+, Co2+, Ni2+, Cu2+, Zn2+, or Cd2+) were measured by potentiometric pH titrations in aqueous solution (25 °C; I = 0.1 mol/L, NaNO3). In the M(H;PME2A6DMAP)+ species, H+ and M2+ (mainly outersphere) are at the phosphonate group; this is relevant for phosphoryl-diester bridges in nucleic acids because, in the present system, there is no indication for a M2+–purine binding. This contrasts, for example, with the complexes formed by 9-[2-(phosphonomethoxy)ethyl]adenine, M(H;PMEA)+, where M2+ is mainly situated at the adenine residue. Application of log [Formula: see text] vs. [Formula: see text] plots for simple phosph(on)ate ligands, R–PO32− (R being a residue that does not affect M2+ binding), proves that all M(PME2A6DMAP) complexes have larger stabilities than what would be expected for a M2+–phosphonate coordination. Comparisons with M(PME–R) complexes, where R is a noncoordinating residue of the (phosphonomethoxy)ethane chain, allow one to conclude that the increased stability is due to the formation of five-membered chelates involving the ether–oxygen of the –CH2–O–CH2–PO32− residue: the percentages of formation of these M(PME2A6DMAP)cl/O chelates, which occur in intramolecular equilibria, vary between 20% (Sr2+, Ba2+) and 50% (Zn2+, Cd2+), up to a maximum of 67% (Cu2+). Any M2+ interaction with N3 or N7 of the purine moiety, as in the parent M(PMEA) complexes, is suppressed by the (C2)NH2 and (C6)N(CH3)2 substituents. This observation, together with the previously determined stacking properties, offers an explanation why PME2A6DMAP2– has remarkable therapeutic effects.