Evaluation of dispersion-corrected density functional theory (B3LYP-DCP) for compounds of biochemical interest

Understanding the surrounding effects on Raman optical activity signatures of a chiral cage system: Cryptophane-PP-111

Cryptophane molecules are cage-like structures consisting in two hemispheres, each made of three benzene rings. These hemispheres are bound together with three O(CH2)nOlinkers of various lengths giving rise to a plethora of cryptophane derivatives. Moreover, they are able to encapsulate neutral guests: CH2Cl2, CHCl3, …; and charged species: Cs+, Tl+, …. Finally, they exhibit chiroptical properties thanks to the anti arrangement of the linkers between the hemispheres. This work focuses on the Raman optical activity (ROA) signatures of Cryptophane-111 (n=1 for each linker). More specifically, we aim at simulating accurately its ROA spectra with and without a xenon atom inside its cavity. Experimental data (Buffeteau et al., 2017) have already demonstrated the effect of the encapsulation in the low-wavenumbers region. To generate the initial structures, we rely on the novel Conformer-Rotamer Ensemble Sampling Tool (CREST) program, developed by S. Grimme and co-workers. This is required due to the flexibility provided by the linkers. The CREST algorithm seems promising and has already been used to sample the potential energy surface (PES) of target systems before the simulation of their vibrational spectroscopies (Eikås et al., 2022). We observe large similarities between the two sets of conformers (one with and one without Xe encapsulated), demonstrating the robustness of the CREST algorithm. For corresponding structures, the presence of xenon pushed the two hemispheres slightly further apart. After optimization at the DFT level, only one unique conformer has a Boltzmann population ratio greater than 1%, pointing out the relative rigidity of the cage. Based on this unique conformer, our simulations are in good agreement with the experimental data. Regarding xenon encapsulation, the (experimental and theoretical) ROA signatures at low wavenumbers are impacted: slight shifts in wavenumbers are observed as well as a decrease in relative ROA intensity for bands around 150 cm-1. The wavenumber shifts were very well reproduced by our simulations, but the experimental decrease in the ROA intensity was unfortunately not reproduced.

Drug–DNA interaction, a joint DFT-D3/MD study on safranal as an anticancer and DNA nanostructure model

In this research, using a combination of quantum mechanics and molecular dynamic (MD) simulations, the interaction of safranal (2,6,6-trimethylcyclohexa-1,3-dien-1-carboxaldehyde) as an anti-cancer drug and Dickerson B-DNA was studied. MD simulations were executed for 35 ns in water. Binding energy analysis in three definite parts of the B-DNA and comparison between different contributions of the binding energy shows that the van der Waals energy part of the interaction is impressive among the standard molecular mechanic energy terms. On the basis of Gibbs energies, it is confirmed that the most important interactions in the safranal complex are related to the A–T and C–G rich regions, which is in agreement with the experimental data. Quantum theory of atoms in molecules and natural bond orbital analyses were applied. A diminution in the electronic chemical potential of the safranal–DNA complex in comparison with the isolated DNA, 0.026 and 0.022 au for the S1 region and 0.012 and 0.017 au for the S2 region, was obtained in the gas phase and water, respectively, which increases the complex stability. An enhancement in the electrophilicity character, during the complexation process, shows the electron charge flux between the safranal and DNA, especially in water. The strengths of the CH⋯O bonds at the center of safranal–DNA interaction were also evaluated. A mean value of 0.06 au for the electron density of the bond critical point of the H⋯O in the complex confirms the H-bond formation during the complexation.

Molecular dynamic simulation and DFT study on the Drug-DNA interaction; Crocetin as an anti-cancer and DNA nanostructure model

In this research, the interaction of Crocetin as an anti-cancer drug and a Dickerson DNA has been investigated. 25 ns molecular dynamic simulations of Crocetin and DNA composed of 12 base pairs and a sequence of d(CGCGAATTCGCG)₂ were done in water. Three definite parts of the B-DNA were considered in analyzing the best interactive site from the thermodynamic point of view. Binding energy analysis showed that van der Waals interaction is the most important part related to the reciprocal O and H atoms of the Crocetin and DNA. Stabilizing interactions, obtained by ΔG calculations, showed that maximum and minimum interactions are related to the S1 and S3 regions, respectively. This means that the most probable van der Waals interaction site of the Dickerson B-DNA and Crocetin is located in the minor groove of DNA. Two sharp peaks at 2.55 and 1.75 Å in radial distribution functions of the PO⋯HO and NH⋯OC parts are related to new hydrogen bonds between the Crocetin and DNA in the complex which can be considered as the driving force of the anti-cancer mechanism of the Crocetin. Average values of 0.3 au and zero for the electron densities of the H⋯O bonds for DNA and complex, obtained by Quantum theory of atoms in molecules (QTAIM), means that the origin of DNA instability after complexation may be related to the H-bond denaturation by Crocetin. Finally, the evaluation of the dispersion interactions using the dispersion functional, -148.76 kcal.mol^-1, confirmed the importance of the dispersion interaction in drug-DNA complex.

Local NH–π interactions involving aromatic residues of proteins: influence of backbone conformation and ππ* excitation on the π H-bond strength, as revealed from studies of isolated model peptides

Gas phase conformer-selective IR spectroscopy combined and relevant quantum chemistry methods document the NH–π interactions in Phe residues.

Isolated neutral peptides

This chapter examines the structural characterisation of isolated neutral amino-acids and peptides. After a presentation of the experimental and theoretical state-of-the-art in the field, a review of the major structures and shaping interactions is presented. Special focus is made on conformationally-resolved studies which enable one to go beyond simple structural characterisation; probing flexibility and excited-state photophysics are given as examples of promising future directions.

An evaluation of substituent effects on aromatic edge-to-face interactions and CF-π versus CH-π interactions using an imino torsion balance model

A selection of imines derived from phenyl t-butyl ketones and substituted 2-phenylethylamines or phenylalanine exhibit slow rotation around the aryl–imino bond at ambient temperature, resulting in a large non-equivalence of the ortho hydrogens in the 1H NMR spectra. This facilitates assessment of aryl substituent effects on the face tilted-T CH–π interaction between a phenyl ring (A) on the imino carbon proximate to the terminal phenyl ring (B). Analysis of the marked temperature dependence of the chemical shift of the interacting ortho hydrogen affords estimates of the opposing enthalpic and entropic factors involved in the rapid equilibrium between the closed edge-to-face conformation and alternative open conformations devoid of a CH–π interaction while in solution. Above ca. 80 °C the entropy term (TΔS) cancels out the enthalpy (ΔH) favouring the closed conformation and open conformations are preferred. Accordingly, commonly reported binding free energies may not be a good measure of the energetic strength of intramolecular aromatic interactions. Investigation of an ortho fluoro substituted compound indicates that a CF–π interaction is at least 1.0 kcal mol−1 weaker in enthalpy than the CH–π interaction. Several X-ray crystal structures depicting an intramolecular edge-to-face interaction are presented.

Absorption spectra of PTCDI: a combined UV-Vis and TD-DFT study

Absorption spectra of PTCDI (3,4,9,10-perylene-tetracarboxylic-diimide) have been investigated in chloroform, N,N'-dimethylformamide (DMF) and dimethylsulfoxide (DMSO). While no signature of assembled PTCDI molecules is observed in chloroform solution, distinct bands assigned to their aggregation have been identified in DMF and DMSO solutions. PTCDI monomers show very similar absorption patterns in chloroform and DMSO solutions. Experimental data, including the vibronic structure of the absorption spectra were explained based on the Franck-Condon approximation and quantum chemical results obtained at PBE0-DCP/6-31+G(d,p) level of theory. Geometry optimization of the first excited state leads to a nice agreement between the calculated adiabatic transition energies and experimental data.

Transmetallation versus β-hydride elimination: the role of 1,4-benzoquinone in chelation-controlled arylation reactions with arylboronic acids

The formation of an atypical, saturated, diarylated, Heck/Suzuki, domino product produced under oxidative Heck reaction conditions, employing arylboronic acids and a chelating vinyl ether, has been investigated by DFT calculations. The calculations highlight the crucial role of 1,4-benzoquinone (BQ) in the reaction. In addition to its role as an oxidant of palladium, which is necessary to complete the catalytic cycle, this electron-deficient alkene opens up a low-energy reaction pathway from the post-insertion σ-alkyl complex. The association of BQ lowers the free-energy barrier for transmetallation of the σ-alkyl complex to create a pathway that is energetically lower than the oxidative Heck reaction pathway. Furthermore, the calculations showed that the reaction is made viable by BQ-mediated reductive elimination and leads to the saturated diarylated product.

Measurement of the interaction of aqueous copper(II) with a model amyloid-β protein fragment — Interference from buffers

For the formation of a complex of Cu2+ with the amyloid-β (Aβ) proxy N-α-dihydrourocanylhistamine (L) in unbuffered aqueous solution (pH ∼ 5.7, 25.0 °C), UV spectrophotometric measurements give a stability constant of 3.8 × 105 L mol–1. This stability constant is within the lower limit of the range of stability constants reported in the literature for complexes of Aβ with Cu2+ — as expected, in view of the smaller number of coordination sites in L. Computer modeling indicates that the Cu2+–L complex is CuL(H2O)22+, with terdentate L bound to Cu2+ via two Nπ atoms and the O atom of the peptide link. Attempts to make stability constant measurements for Cu2+ with L in aqueous solution buffered with Tris/TrisH+/ClO4– to pH near 7.2 were unsuccessful because the Tris base when in large excess over CuL2+ promoted its dissociation to Cu2+ + L by scavenging free Cu2+ as Cu(Tris)(TrisH–1)+, or when in roughly equimolar concentrations formed a ternary adduct, CuL(Tris)2+. The interactions of Cu2+ with Tris buffer were re-examined spectrophotometrically and with the aid of computations that show that the most stable Cu2+–Tris complexes are the syn- and anti-isomers of Cu(Tris)22+, but in the experimental pH ranges these are present as Cu(Tris)(TrisH–1)+. Since Cu2+(aq) is strongly complexed by almost any base capable of forming a buffer system with near-physiological pH, stability constants reported for Cu2+–Aβ complexes in any buffer solution should be regarded with skepticism unless interactions of the buffer with Cu2+ and with CuAβ2+ are taken quantitatively into account. Moreover, in vivo, biological buffers will reduce the physiological importance of Aβ–Cu2+ complexes by competing with Aβ for Cu2+.

Charge delocalization and enhanced acidity in tricationic superelectrophiles

This Article presents the results from studies related to the chemistry of tricationic superelectrophiles. A series of triaryl methanols were ionized in Brønsted superacids, and the corresponding tricationic intermediates were formed. The trications are found to participate in two types of reactions; both are characteristic of highly charged organic cations. One set of reactions occurs through charge migration. A second set of reactions occurs through deprotonation of an unusually acidic site on the tricationic species. One of the tricationic intermediates has been directly observed by low temperature NMR spectroscopy. These highly charged ions and their reactions have also been studied using density functional theory calculations. As a result of charge migration, electron density at a carbocation site is found to increase with progression from monocationic to pentacationic structures.

The CH/π hydrogen bond in chemistry. Conformation, supramolecules, optical resolution and interactions involving carbohydrates

The CH/π hydrogen bond is an attractive molecular force occurring between a soft acid and a soft base. Contribution from the dispersion energy is important in typical cases where aliphatic or aromatic CH groups are involved. Coulombic energy is of minor importance as compared to the other weak hydrogen bonds. The hydrogen bond nature of this force, however, has been confirmed by AIM analyses. The dual characteristic of the CH/π hydrogen bond is the basis for ubiquitous existence of this force in various fields of chemistry. A salient feature is that the CH/π hydrogen bond works cooperatively. Another significant point is that it works in nonpolar as well as polar, protic solvents such as water. The interaction energy depends on the nature of the molecular fragments, CH as well as π-groups: the stronger the proton donating ability of the CH group, the larger the stabilizing effect. This Perspective focuses on the consequence of this molecular force in the conformation of organic compounds and supramolecular chemistry. Implication of the CH/π hydrogen bond extends to the specificity of molecular recognition or selectivity in organic reactions, polymer science, surface phenomena and interactions involving proteins. Many problems, unsettled to date, will become clearer in the light of the CH/π paradigm.