Adenine versus guanine quartets in aqueous solution: dispersion-corrected DFT study on the differences in π-stacking and hydrogen-bonding behavior

Effects of Chalcogen Atoms on Excited-State Double-Proton Transfer Behavior for 3,6-bis(4,5-Dihydroxyoxazo-2-yl)benzene-1,2-diol Derivatives: A Computational Investigation

The impact of the chalcogen atomic electronegativity (O, S, and Se atoms) of new organic molecules on excited-state dynamical reactions is self-evident. Inspired by this kind of distinguished photochemical characteristic, in this work, we performed a computational investigation of chalcogen-substituted 3,6-bis(4,5-dihydroxyoxazo-2-yl)benzene-1,2-diol (BDYBD) derivatives (i.e., BDYBD-O, BDYBD-S, and BDYBD-Se). In this paper, we pay close attention to characteristic BDYBD derivatives that contain intramolecular double hydrogen bonds (O1-H2···N3 and O4-H5···N6). The main goal of this study was to explore how changes in atomic electronegativity affect the way hydrogen bonds interact and how excited molecules affect transfer protons. We go into further detail in the main text of the paper. By fixing our attention to geometrical variations and infrared (IR) vibrational spectra between the S0 and S1 states, exploring hydrogen bonding behaviors using the core-valence bifurcation (CVB) index, and simulating hydrogen bonding energy (EHB) via the atom in molecule (AIM) method, we clarified the photo-induced strengthened dual hydrogen bonding interactions that facilitate the excited-state dual-proton transfer (ESDPT) behavior of BDYBD derivatives. The reorganization of charge stemming from photoexcitation further verifies the tendencies of ESDPT reactions. We relied on constructing potential energy surfaces (PESs) by adopting a restrictive optimization approach, and herein, we finally clarify the gradual ESDPT mechanism of BDYBD derivatives. Particularly, we confirm that the variation in chalcogen atomic electronegativity has a regulatory effect on the ESDPT behavior of BDYBD derivatives; that is, the lower the atomic electronegativity, the more favorable it is for the gradual ESDPT reaction.

Strengthened cooperativity of DNA-based cyclic hydrogen-bonded rosettes by subtle functionalization

Cooperative effects cause extra stabilization of hydrogen-bonded supramolecular systems. In this work we have designed hydrogen-bonded rosettes derived from a guanine-cytosine Janus-type motif with the aim of finding a monomer that enhances the synergy of supramolecular systems. For this, relativistic dispersion-corrected density functional theory computations have been performed. Our proposal involves a monomer with three hydrogen-bonds pointing in the same direction, which translates into shorter bonds, stronger donor-acceptor interactions, and more attractive electrostatic interactions, thus giving rise to rosettes with strengthened cooperativity. This newly designed rosette has triple the cooperativity found for the naturally occurring guanine quadruplex.

Understanding the Interactions of Guanine Quadruplexes with Peptides as Novel Strategies for Diagnosis or Tuning Biological Functions

Guanine quadruplexes (G4s) are nucleic acid structures exhibiting a complex structural behavior and exerting crucial biological functions in both cells and viruses. The specific interactions of peptides with G4s, as well as an understanding of the factors driving the specific recognition are important for the rational design of both therapeutic and diagnostic agents. In this review, we examine the most important studies dealing with the interactions between G4s and peptides, highlighting the strengths and limitations of current analytic approaches. We also show how the combined use of high-level molecular simulation techniques and experimental spectroscopy is the best avenue to design specifically tuned and selective peptides, thus leading to the control of important biological functions.

Structural Insights into the Antiparallel G-Quadruplex in the Presence of K+ and Mg2+ Ions

G-Quadruplex (GQ) is a secondary structural unit of DNA, formed at the telomere region of the chromosome with a high guanine content. It is reported that the GQs can hinder many biological processes. Thus, research thrives to explore the structural stability of GQs. Studies based on circular dichroism (CD) and nuclear magnetic resonance (NMR) experiments established the vital role of cations such as K⁺ and Mg²⁺ in the stability of antiparallel G-quadruplexes (AGQs). However, there is a need to understand how stability in AGQ is attained in the presence of cations. Here, we employed molecular dynamics (MD) simulations, steered MD (SMD) simulations, and QM/MM calculations to understand the biophysical and electronic bases of the stability imparted to AGQs via cation binding. Our results showed that Mg²⁺ prefers to bind in the plane with the guanine tetrad, whereas K⁺ binds in between the AGQ tetrads. Thus, three Mg²⁺ cations or two K⁺ ions are needed to stabilize an AGQ molecule, where each and every tetrad binds to Mg²⁺ more robustly with a higher binding affinity. SMD revealed that the traversal of K⁺ through the AGQ central channel required less force than that of Mg²⁺, illustrating the presence of more strong interactions between Mg²⁺ and AGQ tetrads compared to K⁺. The stabilization in the AGQ tetrads due to cation binding is reassessed by employing ab initio simulations. Mixed QM/MM calculations confirmed that Mg²⁺ binds strongly to AGQ compared to K⁺, and it induces higher interactions between the guanine tetrads. However, K⁺ binding to AGQ induces a higher stabilization energy than Mg²⁺ binding to AGQ tetrads. Despite the higher binding energy, Mg²⁺ binding imparts lower stabilization to AGQ due to its unfavorable fermionic quantum energy.

Reading and erasing of the phosphonium analogue of trimethyllysine by epigenetic proteins

N ε-Methylation of lysine residues in histones plays an essential role in the regulation of eukaryotic transcription. The 'highest' methylation mark, N ε-trimethyllysine, is specifically recognised by N ε-trimethyllysine binding 'reader' domains, and undergoes demethylation, as catalysed by 2-oxoglutarate dependent JmjC oxygenases. We report studies on the recognition of the closest positively charged N ε-trimethyllysine analogue, i.e. its trimethylphosphonium derivative (KPme3), by N ε-trimethyllysine histone binding proteins and Nε-trimethyllysine demethylases. Calorimetric and computational studies with histone binding proteins reveal that H3KP4me3 binds more tightly than the natural H3K4me3 substrate, though the relative differences in binding affinity vary. Studies with JmjC demethylases show that some, but not all, of them can accept the phosphonium analogue of their natural substrates and that the methylation state selectivity can be changed by substitution of nitrogen for phosphorus. The combined results reveal that very subtle changes, e.g. substitution of nitrogen for phosphorus, can substantially affect interactions between ligand and reader domains / demethylases, knowledge that we hope will inspire the development of highly selective small molecules modulating their activity.

Through-Space Stabilization of an Imidazolium Cation by Aromatic Rings

Imidazole-based compounds are widely found in natural products, synthetic molecules, and biomolecules. Noncovalent interactions between the imidazole ring and other functional groups play an important role in determining the function of diverse molecules. However, there is a limited understanding of the underlying noncovalent interactions between imidazoles and aromatic systems. In this work, we report physical-organic chemistry studies on 2-(2,6-diarylphenyl)-1H-imidazoles and their protonated forms to investigate the noncovalent interactions between the central imidazole ring and two flanking aromatic rings possessing substituents at the para/meta position. Hammett analysis revealed that pK_a values and proton affinities correlate well with Hammett σ values of para-substituents at the flanking rings. Additional quantitative Kohn-Sham molecular orbital and energy decomposition analyses reveal that through-space π-π interactions and NH-π interactions contribute to the intramolecular stabilization of the imidazolium cation. The results are important because they clearly demonstrate that the imidazolium cation forms energetically favorable noncovalent interactions with aromatic rings via the through-space effect, a knowledge that can be used in rational drug and catalyst design.

Probing Noncovalent Interactions in [3,3]Metaparacyclophanes

Arene-arene interactions are fundamentally important in molecular recognition. To precisely probe arene-arene interactions in cyclophanes, we designed and synthesized (2,6-phenol)paracyclophanes and (2,6-aniline)paracyclophanes that possess two aromatic rings in close proximity. Fine-tuning the aromatic character of one aromatic ring by fluorine substituents enables investigations on the intramolecular interactions between the electron-rich phenol and aniline with tetra-H- and tetra-F-substituted benzene. pK_a measurements revealed that the tetra-F-template increases the acidity of the phenol (ΔpK_a = 0.55). X-ray crystallography and computational analyses demonstrated that all [3,3]metaparacyclophanes adopt cofacial parallel conformations, implying the presence of π-π stacking interactions. Advanced quantum chemical analyses furthermore revealed that both electrostatic interactions and orbital interactions provide the key contribution to the structure and stability of [3,3]metaparacyclophanes.

B-DNA Structure and Stability: The Role of Nucleotide Composition and Order

We have quantum chemically analyzed the influence of nucleotide composition and sequence (that is, order) on the stability of double-stranded B-DNA triplets in aqueous solution. To this end, we have investigated the structure and bonding of all 32 possible DNA duplexes with Watson-Crick base pairing, using dispersion-corrected DFT at the BLYP-D3(BJ)/TZ2P level and COSMO for simulating aqueous solvation. We find enhanced stabilities for duplexes possessing a higher GC base pair content. Our activation strain analyses unexpectedly identify the loss of stacking interactions within individual strands as a destabilizing factor in the duplex formation, in addition to the better-known effects of partial desolvation. Furthermore, we show that the sequence-dependent differences in the interaction energy for duplexes of the same overall base pair composition result from the so-called "diagonal interactions" or "cross terms". Whether cross terms are stabilizing or destabilizing depends on the nature of the electrostatic interaction between polar functional groups in the pertinent nucleobases.

Probing the Lewis Acidity of Boronic Acids through Interactions with Arene Substituents

Boronic acids are Lewis acids that exist in equilibrium with boronate forms in aqueous solution. Here we experimentally and computationally investigated the Lewis acidity of 2,6‐diarylphenylboronic acids; specially designed phenylboronic acids that possess two flanking aromatic rings with tunable aromatic character. Hammett analysis of 2,6‐diarylphenylboronic acids reveals that their Lewis acidity remains unchanged upon the introduction of EWG/EDG at the distant para position of the flanking aromatic rings. Structural and computational studies demonstrate that polar‐π interactions and solvation effects contribute to the stabilization of boronic acids and boronate forms by aromatic rings. Our physical‐organic chemistry work highlights that boronic acids and boronates can be stabilized by aromatic systems, leading to an important molecular knowledge for rational design and development of boronic acid‐based catalysts and inhibitors of biomedically important proteins.

How Divalent Cations Interact with the Internal Channel Site of Guanine Quadruplexes

The formation of guanine quadruplexes (GQ) in DNA is crucial in telomere homeostasis and regulation of gene expression. Pollution metals can interfere with these DNA superstructures upon coordination. In this work, we study the affinity of the internal GQ channel site towards alkaline earth metal (Mg²⁺, Ca²⁺, Sr²⁺, and Ba²⁺), and (post‐)transition metal (Zn²⁺, Cd²⁺, Hg²⁺, and Pb²⁺) cations using density functional theory computations. We find that divalent cations generally bind to the GQ cavity with a higher affinity than conventional monovalent cations (e. g. K⁺). Importantly, we establish the nature of the cation‐GQ interaction and highlight the relationship between ionic and nuclear charge, and the electrostatic and covalent interactions. The covalent interaction strength plays an important role in the cation affinity and can be traced back to the relative stabilization of cations’ unoccupied atomic orbitals. Overall, our findings contribute to a deeper understanding of how pollution metals could induce genomic instability.

New Insights on the Interaction of Phenanthroline Based Ligands and Metal Complexes and Polyoxometalates with Duplex DNA and G-Quadruplexes

This work provides new insights from our team regarding advances in targeting canonical and non-canonical nucleic acid structures. This modality of medical treatment is used as a form of molecular medicine specifically against the growth of cancer cells. Nevertheless, because of increasing concerns about bacterial antibiotic resistance, this medical strategy is also being explored in this field. Up to three strategies for the use of DNA as target have been studied in our research lines during the last few years: (1) the intercalation of phenanthroline derivatives with duplex DNA; (2) the interaction of metal complexes containing phenanthroline with G-quadruplexes; and (3) the activity of Mo polyoxometalates and other Mo-oxo species as artificial phosphoesterases to catalyze the hydrolysis of phosphoester bonds in DNA. We demonstrate some promising computational results concerning the favorable interaction of these small molecules with DNA that could correspond to cytotoxic effects against tumoral cells and microorganisms. Therefore, our results open the door for the pharmaceutical and medical applications of the compounds we propose.

A unified computational view of DNA duplex, triplex, quadruplex and their donor-acceptor interactions

DNA can assume various structures as a result of interactions at atomic and molecular levels (e.g., hydrogen bonds, π–π stacking interactions, and electrostatic potentials), so understanding of the consequences of these interactions could guide development of ways to produce elaborate programmable DNA for applications in bio- and nanotechnology. We conducted advanced ab initio calculations to investigate nucleobase model structures by componentizing their donor-acceptor interactions. By unifying computational conditions, we compared the independent interactions of DNA duplexes, triplexes, and quadruplexes, which led us to evaluate a stability trend among Watson–Crick and Hoogsteen base pairing, stacking, and even ion binding. For a realistic solution-like environment, the influence of water molecules was carefully considered, and the potassium-ion preference of G-quadruplex was first analyzed at an ab initio level by considering both base-base and ion-water interactions. We devised new structure factors including hydrogen bond length, glycosidic vector angle, and twist angle, which were highly effective for comparison between computationally-predicted and experimentally-determined structures; we clarified the function of phosphate backbone during nucleobase ordering. The simulated tendency of net interaction energies agreed well with that of real world, and this agreement validates the potential of ab initio study to guide programming of complicated DNA constructs.

Intercalation Ability of Novel Monofunctional Platinum Anticancer Drugs: A Key Step in Their Biological Action

Phenanthriplatin (PtPPH) is a monovalent platinum(II)-based complex with a large cytotoxicity against cancer cells. Although the aqua-activated drug has been assumed to be the precursor for DNA damage, it is still under debate whether the way in which that metallodrug attacks to DNA is dominated by a direct binding to a guanine base or rather by an intercalated intermediate product. Aiming to capture the mechanism of action of PtPPH, the present contribution used theoretical tools to systematically assess the sequence of all possible mechanisms on drug activation and reactivity, for example, hydrolysis, intercalation, and covalent damage to DNA. Ab initio quantum mechanical (QM) methods, hybrid QM/QM′ schemes, and independent gradient model approaches are implemented in an unbiased protocol. The performed simulations show that the cascade of reactions is articulated in three well-defined stages: (i) an early and fast intercalation of the complex between the DNA bases, (ii) a subsequent hydrolysis reaction that leads to the aqua-activated form, and (iii) a final formation of the covalent bond between PtPPH and DNA at a guanine site. The permanent damage to DNA is consequently driven by that latter bond to DNA but with a simultaneous π–π intercalation of the phenanthridine into nucleobases. The impact of the DNA sequence and the lateral backbone was also discussed to provide a more complete picture of the forces that anchor the drug into the double helix.

Do Sulfonamides Interact with Aromatic Rings?

Aromatic rings form energetically favorable interactions with many polar groups in chemical and biological systems. Recent molecular studies have shown that sulfonamides can chelate metal ions and form hydrogen bonds, however, it is presently not established whether the polar sulfonamide functionality also interacts with aromatic rings. Here, synthetic, spectroscopic, structural, and quantum chemical analyses on 2,6-diarylbenzenesulfonamides are reported, in which two flanking aromatic rings are positioned close to the central sulfonamide moiety. Fine-tuning the aromatic character by substituents on the flanking rings leads to linear trends in acidity and proton affinity of sulfonamides. This physical-organic chemistry study demonstrates that aromatic rings have a capacity to stabilize sulfonamides via through-space NH-π interactions. These results have implications in rational drug design targeting electron-rich aromatic rings in proteins.

Learning to Model G-Quadruplexes: Current Methods and Perspectives

G-quadruplexes have raised considerable interest during the past years for the development of therapies against cancer. These noncanonical structures of DNA may be found in telomeres and/or oncogene promoters, and it has been observed that the stabilization of such G-quadruplexes may disturb tumor cell growth. Nevertheless, the mechanisms leading to folding and stabilization of these G-quadruplexes are still not well established, and they are the focus of much current work in this field. In seminal works, stabilization was observed to be produced by cations. However, subsequent studies showed that different kinds of small molecules, from planar and nonplanar organic molecules to square-planar and octahedral metal complexes, may also lead to the stabilization of G-quadruplexes. Thus, the comprehension and rationalization of the interaction of these small molecules with G-quadruplexes are also important topics of current interest in medical applications. To shed light on the questions arising from the literature on the formation of G-quadruplexes, their stabilization, and their interaction with small molecules, synergies between experimental studies and computational works are needed. In this review, we mainly focus on in silico approaches and provide a broad compilation of different leading studies carried out to date by different computational methods. We divide these methods into twomain categories: (a) classical methods, which allow for long-timescale molecular dynamics simulations and the corresponding analysis of dynamical information, and (b) quantum methods (semiempirical, quantum mechanics/molecular mechanics, and density functional theory methods), which allow for the explicit simulation of the electronic structure of the system but, in general, are not capable of being used in long-timescale molecular dynamics simulations and, therefore, give a more static picture of the relevant processes.

Understanding alkali metal cation affinities of multi-layer guanine quadruplex DNA

To gain better understanding of the stabilizing interactions between metal ions and DNA quadruplexes, dispersion-corrected density functional theory (DFT-D) based calculations were performed on double-, triple- and four-layer guanine tetrads interacting with alkali metal cations. All computations were performed in aqueous solution that mimics artificial supramolecular conditions where guanine bases assemble into stacked quartets as well as biological environments in which telomeric quadruplexes are formed. To facilitate the computations on these significant larger systems, optimization of the DFT description was performed first by evaluating the performance of partial reduced basis sets. Analysis of the stabilizing interactions between alkali cations and the DNA bases in double and triple-layer guanine quadruplex DNA reproduced the experimental affinity trend of the order Li⁺< Rb⁺ < Na⁺ < K⁺. The desolvation and the size of alkali metal cations are thought to be responsible for the order of affinity. Nevertheless, for the alkali metal cation species individually, the magnitude of the bond energy stays equal for binding as first, second or third cation in double, triple and four-layer guanine quadruplexes, respectively. This is the result of an interplay between a decreasingly stabilizing interaction energy and increasingly stabilizing solvation effects, along the consecutive binding events. This diminished interaction energy is the result of destabilizing electrostatic repulsion between the hosted alkali metal cations. This work emphasizes the stabilizing effect of aqueous solvent on large highly charged biomolecules.

Stabilizing solvent effects and electrostatic repulsion are responsible for the constant alkali metal cation affinity in multi-layer guanine quadruplexes.

Large Non-planar Conjugated Molecule with Strong Intermolecular Interactions Achieved with Homoleptic Zn(II) Complex of Di(5-quinolylethynyl)-tetraphenylazadipyrromethene

Zinc(II) complexes of tetraphenylazadipyrromethenes are potential non-planar n-type conjugated materials. To tune the properties, we installed 5-quinolylethynyl groups at the pyrrolic positions. Compared to the complex with 1-napthylethynyl, we found evidence for stronger intermolecular interactions in the new complex, including much higher overlap integrals in crystals. X-ray analysis revealed unconventional C-H···N hydrogen bonding between two quinolyls of neighboring molecules, pointing to a new strategy for the development of non-planar molecular semiconductors with stronger intermolecular interactions.

Tuning the Binding Strength of Even and Uneven Hydrogen-Bonded Arrays with Remote Substituents

We investigated the tunability of hydrogen bond strength by altering the charge accumulation around the frontier atoms with remote substituents. For pyridine···H2O with NH2 and CN substituted at different positions on pyridine, we find that the electron-withdrawing CN group decreases the negative charge accumulation around the frontier atom N, resulting in weakening of the hydrogen bond, whereas the electron-donating NH2 group increases the charge accumulation around N, resulting in strengthening of the hydrogen bond. By applying these design principles on DDAA-AADD, DADA-ADAD, DAA-ADD, and ADA-DAD hydrogen-bonded dimers, we find that the effect of the substituent is delocalized over the whole molecular system. As a consequence, systems with an equal number of hydrogen bond donor (D) and acceptor (A) atoms are not tunable in a predictable way because of cancellation of counteracting strengthening and weakening effects. Furthermore, we show that the position of the substituent and long-range electrostatics can play an important role as well. Overall, the design principles presented in this work are suitable for monomers with an unequal number of donor and acceptor atoms and can be exploited to tune the binding strength of supramolecular building blocks.

Critical Role of Protons for Emission Quenching of Indoline Dyes in Solution and on Semiconductor Surfaces

By combining time-correlated single photon counting (TCSPC) measurements, density functional theory (DFT), and time-dependent DFT (TD-DFT) calculations, we herein investigate the role of protons, in solutions and on semiconductor surfaces, for the emission quenching of indoline dyes. We show that the rhodanine acceptor moieties, and in particular the carbonyl oxygens, undergo protonation, leading to nonradiative excited-state deactivation. The presence of the carboxylic acid anchoring group, close to the rhodanine moiety, further facilitates the emission quenching, by establishing stable H-bond complexes with carboxylic acid quenchers, with high association constants, in both ground and excited states. This complexation favors the proton transfer process, at a low quencher concentration, in two ways: bringing close to the rhodanine unit the quencher and assisting the proton release from the acid by a partial-concerted proton donation from the close-by carboxylic group to the deprotonated acid. Esterification of the carboxylic group, indeed, inhibits the ground-state complex formation with carboxylic acids and thus the quenching at a low quencher concentration. However, the rhodanine moiety in the ester form can still be the source of emission quenching through dynamic quenching mechanism with higher concentrations of protic solvents or carboxylic acids. Investigating this quenching process on mesoporous ZrO₂, for solar cell applications, also reveals the sensitivity of the adsorbed excited rhodanine dyes toward adsorbed protons on surfaces. This has been confirmed by using an organic base to remove surface protons and utilizing cynao-acrylic dye as a reference dye. Our study highlights the impact of selecting such acceptor group in the structural design of organic dyes for solar cell applications and the overlooked role of protons to quench the excited state for such chemical structures.

Effect of Alkali Metal Cations on Length and Strength of Hydrogen Bonds in DNA Base Pairs

For many years, non-covalently bonded complexes of nucleobases have attracted considerable interest. However, there is a lack of information about the nature of hydrogen bonding between nucleobases when the bonding is affected by metal coordination to one of the nucleobases, and how the individual hydrogen bonds and aromaticity of nucleobases respond to the presence of the metal cation. Here we report a DFT computational study of nucleobase pairs interacting with alkali metal cations. The metal cations contribute to the stabilization of the base pairs to varying degrees depending on their position. The energy decomposition analysis revealed that the nature of bonding between nucleobases does not change much upon metal coordination. The effect of the cations on individual hydrogen bonds were described by changes in VDD charges on frontier atoms, H-bond length, bond energy from NBO analysis, and the delocalization index from QTAIM calculations. The aromaticity changes were determined by a HOMA index.

Substituted adenine quartets: interplay between substituent effect, hydrogen bonding, and aromaticity

Adenine, one of the components of DNA/RNA helices, has the ability to form self-organizing structures with cyclic hydrogen bonds (A₄), similar to guanine quartets. Here, we report a computational investigation of the effect of substituents (X = NO₂, Cl, F, H, Me, and NH₂) on the electronic structure of 9H-adenine and its quartets (A₄-N1, A₄-N3, and A₄-N7). DFT calculations were used to show the relationships between the electronic nature of the substituents, strength of H-bonds in the quartets, and aromaticity of five- and six-membered rings of adenine. We demonstrated how the remote substituent X modifies the proton-donating properties of the NH₂ group involved in the H-bonds within quartets and how the position of the substituent and its electronic nature affect the stability of the quartets. We also showed the possible changes in electronic properties of the substituent and aromaticity of adenine rings caused by tetramer formation. The results indicate that the observed relationships depend on the A₄ type. Moreover, the same substituent can both strengthen and weaken intermolecular interactions, depending on the substitution position.

Substituent effects on hydrogen bonds in adenine quartets and aromaticity of adenine rings depend on the quartet type. A₄-N3 and A₄-N7 quartets are more responsive to the electronic nature of substituents than A₄-N1.

Mechanism of biomolecular recognition of trimethyllysine by the fluorinated aromatic cage of KDM5A PHD3 finger

The understanding of biomolecular recognition of posttranslationally modified histone proteins is centrally important to the histone code hypothesis. Despite extensive binding and structural studies on the readout of histones, the molecular language by which posttranslational modifications on histone proteins are read remains poorly understood. Here we report physical-organic chemistry studies on the recognition of the positively charged trimethyllysine by the electron-rich aromatic cage containing PHD3 finger of KDM5A. The aromatic character of two tryptophan residues that solely constitute the aromatic cage of KDM5A was fine-tuned by the incorporation of fluorine substituents. Our thermodynamic analyses reveal that the wild-type and fluorinated KDM5A PHD3 fingers associate equally well with trimethyllysine. This work demonstrates that the biomolecular recognition of trimethyllysine by fluorinated aromatic cages is associated with weaker cation-π interactions that are compensated by the energetically more favourable trimethyllysine-mediated release of high-energy water molecules that occupy the aromatic cage.

Through-Space Polar-π Interactions in 2,6-Diarylthiophenols

Molecular recognition between polar groups and aromatic molecules is fundamentally important to rational drug design. Although it has been well established that many polar functionalities interact with electron-rich aromatic residues through energetically favorable polar-π interactions, there is a limited understanding of the association between thiols and aromatic systems. Herein we report physical-organic chemistry studies on 2,6-diarylthiophenols that possess the central thiophenol ring and two flanking aromatic rings with tunable electronic properties caused by substituents at distant para position. Hammett analysis revealed that pK_a values and proton affinities correlate well with Hammett sigma values of substituents. Additional energy decomposition analysis supported the conclusion that both through-space SH-π interactions and S^- -π interactions contribute to intramolecular stabilization of 2,6-diarylthiophenols.

Comparison of Molecular Recognition of Trimethyllysine and Trimethylthialysine by Epigenetic Reader Proteins

Gaining a fundamental insight into the biomolecular recognition of posttranslationally modified histones by epigenetic reader proteins is of crucial importance to understanding the regulation of the activity of human genes. Here, we seek to establish whether trimethylthialysine, a simple trimethyllysine analogue generated through cysteine alkylation, is a good trimethyllysine mimic for studies on molecular recognition by reader proteins. Histone peptides bearing trimethylthialysine and trimethyllysine were examined for binding with five human reader proteins employing a combination of thermodynamic analyses, molecular dynamics simulations and quantum chemical analyses. Collectively, our experimental and computational findings reveal that trimethylthialysine and trimethyllysine exhibit very similar binding characteristics for the association with human reader proteins, thereby justifying the use of trimethylthialysine for studies aimed at dissecting the origin of biomolecular recognition in epigenetic processes that play important roles in human health and disease.

Origin of Hydrocarbons Stability from a Computational Perspective: A Case Study of Ortho-Xylene Isomers

It is shown herein that intuitive and text-book steric-clash based interpretation of the higher energy "in-in" xylene isomer (as arising solely from the repulsive CH⋅⋅⋅HC contact) with respect to the corresponding global-minimum "out-out" configuration (where the clashing C-H bonds are tilted out) is misleading. It is demonstrated that the two hydrogen atoms engaged in the CH⋅⋅⋅HC contact in "in-in" are involved in attractive interaction so they cannot explain the lower stability of this isomer. We have proven, based on the arsenal of modern bonding descriptors (EDDB, HOMA, NICS, FALDI, ETS-NOCV, DAFH, FAMSEC, IQA), that in order to understand the relative stability of "in-in" versus "out-out" xylenes isomers one must consider the changes in the electronic structure encompassing the entire molecules as arising from the cooperative action of hyperconjugation, aromaticity and unintuitive London dispersion plus charge delocalization based intra-molecular CH⋅⋅⋅HC interactions.

Human DNA Telomeres in Presence of Oxidative Lesions: The Crucial Role of Electrostatic Interactions on the Stability of Guanine Quadruplexes

By using all atom molecular dynamics simulations, we studied the behavior of human DNA telomere sequences in guanine quadruplex (G4) conformation and in the presence of oxidative lesions, namely abasic sites. In particular, we evidenced that while removing one guanine base induces a significant alteration and destabilization of the involved leaflet, human telomere oligomers tend, in most cases, to maintain at least a partial quadruplex structure, eventually by replacing the empty site with undamaged guanines of different leaflets. This study shows that (i) the disruption of the quadruplex leaflets induces the release of at least one of the potassium cations embedded in the quadruplex channel and that (ii) the electrostatic interactions of the DNA sequence with the aforementioned cations are fundamental to the maintenance of the global quadruplex structure.

The Nature of Hydrogen Bonds: A Delineation of the Role of Different Energy Components on Hydrogen Bond Strengths and Lengths

Hydrogen bonds are a complex interplay between different energy components, and their nature is still subject of an ongoing debate. In this minireview, we therefore provide an overview of the different perspectives on hydrogen bonding. This will be done by discussing the following individual energy components: 1) electrostatic interactions, 2) charge-transfer interactions, 3) π-resonance assistance, 4) steric repulsion, 5) cooperative effects, 6) dispersion interactions and 7) secondary electrostatic interactions. We demonstrate how these energetic factors are essential in a correct description of the hydrogen bond, and discuss several examples of systems whose energetic and geometrical features are not captured by easy-to-use predictive models.

Oxidation of organic diselenides and ditellurides by H2O2 for bioinspired catalyst design

The reactivity of diselenides and ditellurides of general formula (RX)2 (X = Se, Te; R = H, CH3, Ph) toward hydrogen peroxide was studied through a computational approach based on accurate Density Functional Theory (DFT) calculations. The aliphatic and aromatic dichalcogenides have been chosen in light of their activity in glutathione peroxidase (GPx)-like catalytic cycles and their promising features as efficient antioxidant compounds. The reaction products, the energetics and the mechanistic details of these oxidations are discussed. Analogous disulfides are included in our analysis for completeness. We find that the barrier for oxidation of dichalcogenides decreases from disulfides to diselenides to ditellurides. On the other hand, variation of the substituents at the chalcogen nucleus has relatively little effect on the reactivity.

Secondary Electrostatic Interaction Model Revised: Prediction Comes Mainly from Measuring Charge Accumulation in Hydrogen-Bonded Monomers

The secondary electrostatic interaction (SEI) model is often used to predict and explain relative hydrogen bond strengths of self-assembled systems. The SEI model oversimplifies the hydrogen-bonding mechanisms by viewing them as interacting point charges, but nevertheless experimental binding strengths are often in line with the model’s predictions. To understand how this rudimentary model can be predictive, we computationally studied two tautomeric quadruple hydrogen-bonded systems, DDAA-AADD and DADA-ADAD. Our results reveal that when the proton donors D (which are electron-donating) and the proton acceptors A (which are electron-withdrawing) are grouped together as in DDAA, there is a larger accumulation of charge around the frontier atoms than when the proton donor and acceptor groups are alternating as in DADA. This accumulation of charge makes the proton donors more positive and the proton acceptors more negative, which enhances both the electrostatic and covalent interactions in the DDAA dimer. The SEI model is thus predictive because it provides a measure for the charge accumulation in hydrogen-bonded monomers. Our findings can be understood from simple physical organic chemistry principles and provide supramolecular chemists with meaningful understanding for tuning hydrogen bond strengths and thus for controlling the properties of self-assembled systems.

Designing Self-Assembled Rosettes: Why Ammeline is a Superior Building Block to Melamine

In supramolecular chemistry, the rational design of self-assembled systems remains a challenge. Herein, hydrogen-bonded rosettes of melamine and ammeline have been theoretically examined by using dispersion-corrected density functional theory (DFT-D). Our bonding analyses, based on quantitative Kohn-Sham molecular orbital theory and corresponding energy decomposition analyses (EDA), show that ammeline is a much better building block than melamine for the fabrication of cyclic complexes based on hydrogen bonds. This superior capacity is explained by both stronger hydrogen bonding and the occurrence of a strong synergy.

Supramolecular Self-Sorting Networks using Hydrogen-Bonding Motifs

A current objective in supramolecular chemistry is to mimic the transitions between complex self-sorted systems that represent a hallmark of regulatory function in nature. In this work, a self-sorting network, comprising linear hydrogen motifs, was created. Selecting six hydrogen-bonding motifs capable of both high-fidelity and promiscuous molecular recognition gave rise to a complex self-sorting system, which included motifs capable of both narcissistic and social self-sorting. Examination of the interactions between individual components, experimentally and computationally, provided a rationale for the product distribution during each phase of a cascade. This reasoning holds through up to five sequential additions of six building blocks, resulting in the construction of a biomimetic network in which the presence or absence of different components provides multiple unique pathways to distinct self-sorted configurations.

RNA versus DNA G-Quadruplex: The Origin of Increased Stability

DNA quadruplexes have been the subject of investigation because of their biological relevance and because of their potential application in supramolecular chemistry. Similarly, RNA quadruplexes are now gaining increasing attention. Although DNA and RNA quadruplexes are structurally very similar, the latter show higher stability. In this study we report dispersion‐corrected density functional theory (DFT‐D) quantum chemical calculations that were undertaken to understand the difference in stabilities of RNA and DNA quadruplexes. The smallest meaningful model of a stack of quartets, interacting with alkali metal cations, was simulated in an aqueous environment. The energy decomposition analysis allows for in‐depth examination of the interaction energies, emphasising the role of noncovalent interactions and better electrostatics in determining RNA‐GQs higher stabilities, particularly pinpointing the role of the extra 2′‐OH groups. Furthermore, our computations present new insights on why the cation is required for self‐assembly: unexpectedly the cation is not necessary to relieve the repulsion between the oxygen atoms in the central cavity, but it is needed to overcome the entropic penalty.

Comparison of the DFT-SAPT and Canonical EDA Schemes for the Energy Decomposition of Various Types of Noncovalent Interactions

Interaction energies computed with density functional theory can be divided into physically meaningful components by symmetry-adapted perturbation theory (DFT-SAPT) or the canonical energy decomposition analysis (EDA). In this work, the decomposition results obtained by these schemes were compared for more than 200 hydrogen-, halogen-, and pnicogen-bonded, dispersion-bound, and mixed complexes to investigate their similarity in the evaluation of the nature of noncovalent interactions. BLYP functional with D3(BJ) correction was used for the EDA scheme, whereas asymptotically corrected PBE0 functional for DFT-SAPT provided some of the best combinations for description of noncovalent interactions. Both schemes provide similar results concerning total interaction energies and insight into the individual energy components. For most complexes, the dominant energetic term was identified equally by both decomposition schemes. Because the canonical EDA is computationally less demanding than the DFT-SAPT, the former can be especially used in cases where the systems investigated are very large.

The influence of substituents and the environment on the NMR shielding constants of supramolecular complexes based on A-T and A-U base pairs

In the present study, we have theoretically analyzed supramolecular complexes based on the Watson-Crick A-T and A-U base pairs using dispersion-corrected density functional theory (DFT). Hydrogen atoms H8 and/or H6 in the natural adenine and thymine/uracil bases were replaced, respectively, by substituents X8, Y6 = NH^-, NH₂, NH₃⁺ (N series), O^-, OH, OH₂⁺ (O series), F, Cl or Br (halogen series). We examined the effect of the substituents on the hydrogen-bond lengths, strength and bonding mechanism, and the NMR shielding constants of the C2-adenine and C2-thymine/uracil atoms in the base pairs. The general belief in the literature that there is a direct connection between changes in the hydrogen-bond strength and the C2-adenine shielding constant is conclusively rejected by our computations.

Nature of Intramolecular Resonance Assisted Hydrogen Bonding in Malonaldehyde and Its Saturated Analogue

The nature of resonance-assisted hydrogen bonds (RAHB) is still subject of an ongoing debate. We therefore analyzed the σ and π charge redistributions associated with the formation of intramolecular hydrogen bonds in malonaldehyde (MA) and its saturated analogue 3-hydroxypropanal (3-OH) and addressed the question whether there is a resonance assistance phenomenon in the sense of a synergistic interplay between the σ and π electron systems. Our quantum chemical calculations at the BP86/TZ2P level of theory show that the π charge flow is indeed in line with the Lewis structure as proposed by the RAHB model. This typical rearrangement of charge is only present in the unsaturated system, and not in its saturated analogue. Resonance in the π electron system assists the intramolecular hydrogen bond by reducing the hydrogen bond distance, and by providing an additional stabilizing component to the net bonding energy. The σ orbital interaction plays an important role in the enhanced hydrogen bond strength in MA as well. However, there is no resonance assistance in the sense of an interplay between σ charge transfer and π polarization; σ and π contribute independently from each other.

Symmetry-adapted perturbation theory potential for the adenine dimer

A new ab initio intermolecular interaction potential for the adenine dimer has been developed.

Extended lead(ii) architectures engineered via tetrel bonding interactions

The cooperative action of multiple non-covalent interactions allowed for the assembly of a new series of Pb^II coordination compounds.

Interstrand cross-linking implies contrasting structural consequences for DNA: insights from molecular dynamics

Oxidatively-generated interstrand cross-links rank among the most deleterious DNA lesions. They originate from abasic sites, whose aldehyde group can form a covalent adduct after condensation with the exocyclic amino group of purines, sometimes with remarkably high yields. We use explicit solvent molecular dynamics simulations to unravel the structures and mechanical properties of two DNA sequences containing an interstrand cross-link. Our simulations palliate the absence of experimental structural and stiffness information for such DNA lesions and provide an unprecedented insight into the DNA embedding of lesions that represent a major challenge for DNA replication, transcription and gene regulation by preventing strand separation. Our results based on quantum chemical calculations also suggest that the embedding of the ICL within the duplex can tune the reaction profile, and hence can be responsible for the high difference in yields of formation.

Insights into the self-assembly steps of cyanuric acid toward rosette motifs: a DFT study

The nature of non-covalent interactions in self-assembling systems is a topic that has aroused great attention in literature. In this field, the 1,3,5-triazinane-2,4,6-trione or cyanuric acid (CA) is one of the most widely used molecules to formulate self-assembled materials or monolayers. In the present work, a variety of molecular aggregates of CA are examined using three different DFT functionals (B3LYP, B3LYP-D3, and ω-B97XD) in the framework of the quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analysis. Herein, a step by step aggregation path is proposed and the origin of cooperative effects is also examined. It is shown that a greater cooperativity is not always associated with a greater binding energy, and the greatest cooperative effect occurs with highly directional hydrogen bonds. The intramolecular charge transfers play a key role in this effect. Graphical abstract The noncovalent interactions in cyanuric acid supramolecules were analyzed. The calculations provide insights into the self-assembly steps from dimers to rosette-like motif. The complexes with collinear hydrogen bonds show positive cooperativity, while in the arrangement with double hydrogen bonds the cooperative effect is essentially zero.

Hydrogen-Bond Strength of CC and GG Pairs Determined by Steric Repulsion: Electrostatics and Charge Transfer Overruled

Theoretical and experimental studies have elucidated the bonding mechanism in hydrogen bonds as an electrostatic interaction, which also exhibits considerable stabilization by charge transfer, polarization, and dispersion interactions. Therefore, these components have been used to rationalize the differences in strength of hydrogen‐bonded systems. A completely new viewpoint is presented, in which the Pauli (steric) repulsion controls the mechanism of hydrogen bonding. Quantum chemical computations on the mismatched DNA base pairs CC and GG (C=cytosine, G=guanine) show that the enhanced stabilization and shorter distance of GG is determined entirely by the difference in the Pauli repulsion, which is significantly less repulsive for GG than for CC. This is the first time that evidence is presented for the Pauli repulsion as decisive factor in relative hydrogen‐bond strengths and lengths.

DFT study of water adsorption on lignite molecule surface

High moisture content is a main characteristic of low-rank coal, such as lignite. Numerous oxygen containing functional groups in lignite make it represent some special properties, and these functional groups affect the adsorption mechanisms of water molecules on lignite surface. This study reports some typical water · · · lignite conformations, along with a detailed analysis of the geometry, electrostatic potential distribution, reduced density gradient of interaction, and interaction energy decomposition. The results show that water molecules tend to aggregate around functional groups, and hydrogen bonds play a dominant role in the interaction. The adsorption energy of water cluster on lignite surface is larger than that of isolated water molecule, a good linear relationship between the interaction distance and adsorption energy of layers has been found. Since water is a polar molecule, the local minima and maxima of electrostatic potential in conformations increase along with more water adsorbing on lignite surface. Reduced density gradient analysis shows that H-bonds, van der Waals interaction, and a little steric make up the interaction between water cluster and lignite molecule. In these studied conformations which mainly are H-bond complexes, electrostatic and exchange repulsion play a dominant role, whereas polarization and dispersion make relatively small contribution to the interaction. Attractive and repulsive interaction both affect the stability of water · · · lignite conformations.