Effects of Extending the Computational Model on DNA–Protein T-shaped Interactions: The Case of Adenine–Histidine Dimers

Radiationless deactivation pathways versus H-atom elimination from the N-H bond photodissociation in PhNH2-(Py)n (n = 1,2) complexes

Minimum energy structures of the ground and lowest excited states of aniline (PhNH₂) solvated by pyridine (Py) show that the clusters formed are stabilized by hydrogen bonds in which only one or both hydrogen atoms of the NH₂ group take part. Two different N-H bonds photodissociation in PhNH₂-(Py)_n (n = 1,2) complexes, free and hydrogen bonded have been studied by analyzing excited state potential energy surfaces. In the first one, only N-H bonds engaged in hydrogen bonding in these complexes are considered. RICC2 calculations of potential energy (PE) profiles indicate that all photochemical reaction paths along N-H stretching occur mainly via the proton-coupled electron transfer (PCET) mechanism. The repulsive charge transfer ¹ππ*(CT) state dominates the PE profiles, leading to low-lying ¹ππ*(CT)/S₀ conical intersections and thus provide channels for ultrafast radiationless deactivation of the electronic excitation or stabilization to biradical complexes. The second photoreaction consists of a direct dissociation along the free N-H bond of the NH₂ group. It has been shown that this process is played by excited singlet states of ¹πσ* character having repulsive potential energy profiles with respect to the stretching of N-H bond, which dissociates over an exit barrier about 0.5 eV giving rise to the formation of a ¹πσ*/S₀ conical intersection. This may cause an internal conversion to the ground state or may lead to H-atom elimination. This photophysical process is the same in both planar and T-shaped conformers of the PhNH₂-Py monomer complex. Our findings reveal that there is no single dominating path in the photodissociation of N-H bonds in PhNH₂-(Py)_n complexes, but rather a variety of paths involving H-atom elimination and several quenching mechanisms.

RNA-protein interactions in an unstructured context

Despite their importance, our understanding of noncovalent RNA-protein interactions is incomplete. This especially concerns the binding between RNA and unstructured protein regions, a widespread class of such interactions. Here, we review the recent experimental and computational work on RNA-protein interactions in an unstructured context with a particular focus on how such interactions may be shaped by the intrinsic interaction affinities between individual nucleobases and protein side chains. Specifically, we articulate the claim that the universal genetic code reflects the binding specificity between nucleobases and protein side chains and that, in turn, the code may be seen as the Rosetta stone for understanding RNA-protein interactions in general.

Absolute binding-free energies between standard RNA/DNA nucleobases and amino-acid sidechain analogs in different environments

Despite the great importance of nucleic acid-protein interactions in the cell, our understanding of their physico-chemical basis remains incomplete. In order to address this challenge, we have for the first time determined potentials of mean force and the associated absolute binding free energies between all standard RNA/DNA nucleobases and amino-acid sidechain analogs in high- and low-dielectric environments using molecular dynamics simulations and umbrella sampling. A comparison against a limited set of available experimental values for analogous systems attests to the quality of the computational approach and the force field used. Overall, our analysis provides a microscopic picture behind nucleobase/sidechain interaction preferences and creates a unified framework for understanding and sculpting nucleic acid-protein interactions in different contexts. Here, we use this framework to demonstrate a strong relationship between nucleobase density profiles of mRNAs and nucleobase affinity profiles of their cognate proteins and critically analyze a recent hypothesis that the two may be capable of direct, complementary interactions.

DNA-protein π-interactions in nature: abundance, structure, composition and strength of contacts between aromatic amino acids and DNA nucleobases or deoxyribose sugar

Four hundred twenty-eight high-resolution DNA-protein complexes were chosen for a bioinformatics study. Although 164 crystal structures (38% of those searched) contained no interactions, 574 discrete π-contacts between the aromatic amino acids and the DNA nucleobases or deoxyribose were identified using strict criteria, including visual inspection. The abundance and structure of the interactions were determined by unequivocally classifying the contacts as either π-π stacking, π-π T-shaped or sugar-π contacts. Three hundred forty-four nucleobase-amino acid π-π contacts (60% of all interactions identified) were identified in 175 of the crystal structures searched. Unprecedented in the literature, 230 DNA-protein sugar-π contacts (40% of all interactions identified) were identified in 137 crystal structures, which involve C-H···π and/or lone-pair···π interactions, contain any amino acid and can be classified according to sugar atoms involved. Both π-π and sugar-π interactions display a range of relative monomer orientations and therefore interaction energies (up to -50 (-70) kJ mol(-1) for neutral (charged) interactions as determined using quantum chemical calculations). In general, DNA-protein π-interactions are more prevalent than perhaps currently accepted and the role of such interactions in many biological processes may yet to be uncovered.

Aromatic interactions modulate the 5'-base selectivity of the DNA-binding autoantibody ED-10

We present detailed computational analyses of the binding of four dinucleotides to a highly sequence-selective single-stranded DNA (ssDNA) binding antibody (ED-10) and selected point mutants. Anti-DNA antibodies are central to the pathogenesis of systemic lupus erythematosus (SLE), and a more complete understanding of the mode of binding of DNA and other ligands will be necessary to elucidate the role of anti-DNA antibodies in the kidney inflammation associated with SLE. Classical molecular mechanics based molecular dynamics simulations and density functional theory (DFT) computations were applied to pinpoint the origin of selectivity for the 5'-nucleotide. In particular, the strength of interactions between each nucleotide and the surrounding residues were computed using MMGBSA as well as DFT applied to a cluster model of the binding site. The results agree qualitatively with experimental binding free energies, and indicate that π-stacking, CH/π, NH/π, and hydrogen-bonding interactions all contribute to 5'-base selectivity in ED-10. Most importantly, the selectivity for dTdC over dAdC arises primarily from differences in the strength of π-stacking and XH/π interactions with the surrounding aromatic residues; hydrogen bonds play little role. These data suggest that a key Tyr residue, which is not present in other anti-DNA antibodies, plays a key role in the 5'-base selectivity, while we predict that the mutation of a single Trp residue can tune the selectivity for dTdC over dAdC.

Examination of tyrosine/adenine stacking interactions in protein complexes

The π-stacking interactions between tyrosine amino acid side chains and adenine-bearing ligands are examined. Crystalline protein structures from the protein data bank (PDB) exhibiting face-to-face tyrosine/adenine arrangements were used to construct 20 unique 4-methylphenol/N9-methyladenine (p-cresol/9MeA) model systems. Full geometry optimization of the 20 crystal structures with the M06-2X density functional theory method identified 11 unique low-energy conformations. CCSD(T) complete basis set (CBS) limit interaction energies were estimated for all of the structures to determine the magnitude of the interaction between the two ring systems. CCSD(T) computations with double-ζ basis sets (e.g., 6-31G*(0.25) and aug-cc-pVDZ) indicate that the MP2 method overbinds by as much as 3.07 kcal mol(-1) for the crystal structures and 3.90 kcal mol(-1) for the optimized structures. In the 20 crystal structures, the estimated CCSD(T) CBS limit interaction energy ranges from -4.00 to -6.83 kcal mol(-1), with an average interaction energy of -5.47 kcal mol(-1), values remarkably similar to the corresponding data for phenylalanine/adenine stacking interactions. Geometry optimization significantly increases the interaction energies of the p-cresol/9MeA model systems. The average estimated CCSD(T) CBS limit interaction energy of the 11 optimized structures is 3.23 kcal mol(-1) larger than that for the 20 crystal structures.