Role of nitrogen atom in aromatic stacking

Systematic Study of Heteroarene Stacking Using a Congeneric Set of Molecular Glues for Procaspase-6

Heteroaromatic stacking interactions are important in drug binding, supramolecular chemistry, and materials science, making protein-ligand model systems of these interactions of considerable interest. Here we studied 30 congeneric ligands that each present a distinct heteroarene for stacking between tyrosine residues at the dimer interface of procaspase-6. Complex X-ray crystal structures of 10 analogs showed that stacking geometries were well conserved, while high-accuracy computations showed that heteroarene stacking energy was well correlated with predicted overall ligand binding energies. Empirically determined KD values in this system thus provide a useful measure of heteroarene stacking with tyrosine. Stacking energies are discussed in the context of torsional strain, the number and positioning of heteroatoms, tautomeric state, and coaxial orientation of heteroarene in the stack. Overall, this study provides an extensive data set of empirical and high-level computed binding energies in a versatile new protein-ligand system amenable to studies of other intermolecular interactions.

Atomistic View of the Energy Transfer in a Fluorophore-Functionalized Gold Nanocluster

Understanding the dynamics of Förster resonance energy transfer (FRET) in fluorophore-functionalized nanomaterials is critical for developing and utilizing such materials in biomedical imaging and optical sensing applications. However, structural dynamics of noncovalently bound systems have a significant effect on the FRET properties affecting their applications in solutions. Here, we study the dynamics of the FRET in atomistic detail by disclosing the structural dynamics of the noncovalently bound azadioxotriangulenium dye (KU) and atomically precise gold nanocluster (Au₂₅(p-MBA)₁₈, p-MBA = para-mercaptobenzoic acid) with a combination of experimental and computational methods. Two distinct subpopulations involved in the energy transfer process between the KU dye and the Au₂₅(p-MBA)₁₈ nanoclusters were resolved by time-resolved fluorescence experiments. Molecular dynamics simulations revealed that KU is bound to the surface of Au₂₅(p-MBA)₁₈ by interacting with the p-MBA ligands as a monomer and as a π-π stacked dimer where the center-to-center distance of the monomers to Au₂₅(p-MBA)₁₈ is separated by ∼0.2 nm, thus explaining the experimental observations. The ratio of the observed energy transfer rates was in reasonably good agreement with the well-known 1/R⁶ distance dependence for FRET. This work discloses the structural dynamics of the noncovalently bound nanocluster-based system in water solution, providing new insight into the dynamics and energy transfer mechanism of the fluorophore-functionalized gold nanocluster at an atomistic level.

Effects of heteroatoms on aromatic pi-pi interactions: benzene-pyridine and pyridine dimer

Heteroatoms are found in many noncovalent complexes which are of biological importance. The effect of heteroatoms on pi-pi interactions is assessed via highly accurate quantum chemical computations for the two simplest cases of interactions between aromatic molecules containing heteroatoms, namely, benzene-pyridine and pyridine dimer. Benchmark quality estimated coupled-cluster through perturbative triples [CCSD(T)] binding energies are computed near the complete basis set limit. Comparisons to the benzene dimer are made to determine the contributions from heteroatoms. The presence of a heteroatom reduces the spatial extent of the pi-electron cloud and polarizability of pyridine as compared to benzene. As a result, the magnitude of the dispersion, exchange, and induction interactions in benzene-pyridine and pyridine dimer is generally reduced as compared to those for the benzene dimer. Benzene-pyridine and pyridine dimer bind more strongly than the benzene dimer in several configurations, and in contrast to the benzene dimer, parallel-displaced configurations can be significantly preferred over T-shaped configurations. Hydrogens para to a heteroatom are more effective "pi-hydrogen bond" donors, but aromatic rings with heteroatoms are worse "pi-hydrogen bond" acceptors.

Probing phenylalanine/adenine pi-stacking interactions in protein complexes with explicitly correlated and CCSD(T) computations

To examine the effects of pi-stacking interactions between aromatic amino acid side chains and adenine bearing ligands in crystalline protein structures, 26 toluene/(N9-methyl)adenine model configurations have been constructed from protein/ligand crystal structures. Full geometry optimizations with the MP2 method cause the 26 crystal structures to collapse to six unique structures. The complete basis set (CBS) limit of the CCSD(T) interaction energies has been determined for all 32 structures by combining explicitly correlated MP2-R12 computations with a correction for higher-order correlation effects from CCSD(T) calculations. The CCSD(T) CBS limit interaction energies of the 26 crystal structures range from -3.19 to -6.77 kcal mol (-1) and average -5.01 kcal mol (-1). The CCSD(T) CBS limit interaction energies of the optimized complexes increase by roughly 1.5 kcal mol (-1) on average to -6.54 kcal mol (-1) (ranging from -5.93 to -7.05 kcal mol (-1)). Corrections for higher-order correlation effects are extremely important for both sets of structures and are responsible for the modest increase in the interaction energy after optimization. The MP2 method overbinds the crystal structures by 2.31 kcal mol (-1) on average compared to 4.50 kcal mol (-1) for the optimized structures.

Effect of base stacking on the acid-base properties of the adenine cation radical [A*+] in solution: ESR and DFT studies

In this study, the acid-base properties of the adenine cation radical are investigated by means of experiment and theory. Adenine cation radical (A*(+)) is produced by one-electron oxidation of dAdo and of the stacked DNA-oligomer (dA)6 by Cl2*(-) in aqueous glass (7.5 M LiCl in H2O and in D2O) and investigated by ESR spectroscopy. Theoretical calculations and deuterium substitution at C8-H and N6-H in dAdo aid in our assignments of structure. We find the pKa value of A*(+) in this system to be ca. 8 at 150 K in seeming contradiction to the accepted value of < or = 1 at ambient temperature. However, upon thermal annealing to > or = 160 K, complete deprotonation of A*(+) occurs in dAdo in these glassy systems even at pH ca. 3. A*(+) found in (dA)6 at 150 K also deprotonates on thermal annealing. The stability of A*(+) at 150 K in these systems is attributed to charge delocalization between stacked bases. Theoretical calculations at various levels (DFT B3LYP/6-31G*, MPWB95, and HF-MP2) predict binding energies for the adenine stacked dimer cation radical of 12 to 16 kcal/mol. Further DFT B3LYP/6-31G* calculations predict that, in aqueous solution, monomeric A*(+) should deprotonate spontaneously (a predicted pKa of ca. -0.3 for A*(+)). However, the charge resonance stabilized dimer AA*(+) is predicted to result in a significant barrier to deprotonation and a calculated pKa of ca. 7 for the AA*(+) dimer which is 7 pH units higher than the monomer. These theoretical and experimental results suggest that A*(+) isolated in solution and A*(+) in adenine stacks have highly differing acid-base properties resulting from the stabilization induced by hole delocalization within adenine stacks.

Intriguing relations of interaction energy components in stacked nucleic acids

Major components of the interaction energy that define several approximate levels starting from second order Möller-Plesset theory were studied for 58 stacked nucleic acid dimers. They included typical B-DNA and A-DNA structures, and selected published geometries. A survey of the various terms yields an unexpected correlation between the Pauli exchange and dispersion or correlation terms, which holds for each class of similar planar geometries and for various basis sets. The geometries that exhibit these correlations span a specific range of molecular overlaps when compared to a model benzene-pyridine stacked dimer. Also, the relationship between electrostatic interactions and MP2 stabilization energies reported earlier is confirmed and a prediction interval of practical relevance is estimated.