N−H···π Interactions in Indole···Benzene-h6,d6 and Indole···Benzene-h6,d6 Radical Cation Complexes. Mass Analyzed Threshold Ionization Experiments and Correlated ab Initio Quantum Chemical Calculations

Cumulative Millisecond-Long Sampling for a Comprehensive Energetic Evaluation of Aqueous Ionic Liquid Effects on Amino Acid Interactions

The interactions of amino acid side-chains confer diverse energetic contributions and physical properties to a protein's stability and function. Various computational tools estimate the effect of changing a given amino acid on the protein's stability based on parametrized (free) energy functions. When parametrized for the prediction of protein stability in water, such energy functions can lead to suboptimal results for other solvents, such as ionic liquids (IL), aqueous ionic liquids (aIL), or salt solutions. However, to our knowledge, no comprehensive data are available describing the energetic effects of aIL on intramolecular protein interactions. Here, we present the most comprehensive set of potential of mean force (PMF) profiles of pairwise protein-residue interactions to date, covering 50 relevant interactions in water, the two biotechnologically relevant aIL [BMIM/Cl] and [BMIM/TfO], and [Na/Cl]. These results are based on a cumulated simulation time of >1 ms. aIL and salt ions can weaken, but also strengthen, specific residue interactions by more than 3 kcal mol^-1, depending on the residue pair, residue-residue configuration, participating ions, and concentration, necessitating considering such interactions specifically. These changes originate from a complex interplay of competitive or cooperative noncovalent ion-residue interactions, changes in solvent structural dynamics, or unspecific charge screening effects and occur at the contact distance but also at larger, solvent-separated distances. This data provide explanations at the atomistic and energetic levels for complex IL effects on protein stability and should help improve the prediction accuracies of computational tools that estimate protein stability based on (free) energy functions.

Different self-assemblies and absorption–emission properties of the picrate salts of aromatic amine or heterocycle linked oximes

Different sub-assemblies and fluorescence quenching in picrate salts of an aromatic amine and of three different heterocycle tethered aldoximes are described.

Interplay among Electrostatic, Dispersion, and Steric Interactions: Spectroscopy and Quantum Chemical Calculations of π-Hydrogen Bonded Complexes

π-Hydrogen bonding interactions are ubiquitous in both materials and biology. Despite their relatively weak nature, great progress has been made in their investigation by experimental and theoretical methods, but this becomes significantly more complicated when secondary intermolecular interactions are present. In this study, the effect of successive methyl substitution on the supramolecular structure and interaction energy of indole⋅⋅⋅methylated benzene (ind⋅⋅⋅n-mb, n=1-6) complexes is probed through a combination of supersonic jet experiments and benchmark-quality quantum chemical calculations. It is demonstrated that additional secondary interactions introduce a subtle interplay among electrostatic and dispersion forces, as well as steric repulsion, which fine-tunes the overall structural motif. Resonant two-photon ionization and IR-UV double-resonance spectroscopy techniques are used to probe jet-cooled ind⋅⋅⋅n-mb (n=2, 3, 6) complexes, with redshifting of the N-H IR stretching frequency showing that increasing the degree of methyl substitution increases the strength of the primary N-H⋅⋅⋅π interaction. Ab initio harmonic frequency and binding energy calculations confirm this trend for all six complexes. Electronic spectra of the three dimers are broad and structureless, with quantum chemical calculations revealing that this is likely to be due to multiple tilted conformations of each dimer possessing similar stabilization energies.

A tyrosine-tryptophan dyad and radical-based charge transfer in a ribonucleotide reductase-inspired maquette

In class 1a ribonucleotide reductase (RNR), a substrate-based radical is generated in the α2 subunit by long-distance electron transfer involving an essential tyrosyl radical (Y122O·) in the β2 subunit. The conserved W48 β2 is ∼10 Å from Y122OH; mutations at W48 inactivate RNR. Here, we design a beta hairpin peptide, which contains such an interacting tyrosine–tryptophan dyad. The NMR structure of the peptide establishes that there is no direct hydrogen bond between the phenol and the indole rings. However, electronic coupling between the tyrosine and tryptophan occurs in the peptide. In addition, downshifted ultraviolet resonance Raman (UVRR) frequencies are observed for the radical state, reproducing spectral downshifts observed for β2. The frequency downshifts of the ring and CO bands are consistent with charge transfer from YO· to W or another residue. Such a charge transfer mechanism implies a role for the β2 Y-W dyad in electron transfer.

Tyrosine-tryptophan dyads are known to mediate electron transfer reactions in a range of different proteins. Here, the authors study a beta hairpin peptide, probing the tyrosine-tryptophan interaction and showing no hydrogen bonding but rather charge transfer between the tyrosyl radical and tryptophan'.

Halogen vs hydrogen bonding in thiazoline-2-thione stabilization with σ- and π-electron acceptors adducts: theoretical and experimental study

Molecular charge-transfer complexes (CT) between thiazoline-2-thione (THZ) and different σ- (I2) and π-acceptors (Tetracyanoethylene (TCNE), 2,3-dichloro-5,6-dicyano-1,4-benzoquinone (DDQ), and 2,3,5,6-tetrachloro-1,4-benzoquinone (CHL)) were investigated. UV-Vis absorption spectroscopy and theoretical calculations using both MP2/aug-cc-pVDZ-PP and B3LYP/6-311++G(d,p) level of theory were corroborated to study the nature of the stabilizing forces for THZ-I2, THZ-DDQ, THZ-TCNE, and THZ-CHL. Halogen bonding (XB) was the stabilizing attractive force in THZ-I2 and THZ-CHL whereas; hydrogen bonding (HB) was dominated in both THZ-TCNE, and THZ-DDQ complexes. Formation constant (K), extinction coefficient (ɛ), thermodynamic parameters such as enthalpy change (ΔH), entropy (ΔS), and Gibbs free energy (ΔG) were measured in different solvents.

Mimicking anesthetic-receptor interactions in jets: the propofol-isopropanol cluster

The interaction of the general anesthetic propofol with an individual residue of threonine in the membrane receptors has been modeled in the gas phase by examining the adduct of propofol with the isopropanol side-chain. We determined the structural preferences of the cluster using a combination of mass-resolved laser spectroscopy and quantum mechanical calculations. The first electronic transition of propofol-isopropanol was recorded with vibrational resolution using resonant two-photon ionization (R2PI) and ion dip IR spectroscopy. The spectra obtained were compared with density-functional calculations (DFT) using the M06-2X functional in order to obtain the cluster's structure. Three isomers have been detected. The results suggest that propofol acts as a Brønsted acid, donating a proton to the isopropanol molecule in a conformation that resembles that of propofol-water, but displaced towards the aromatic ring, due to the interaction with the aliphatic side of isopropanol. The higher affinity of propofol for isopropanol compared to water may correlate with the biological role of propofol at the protein binding site. On the other hand, propofol shows a similar affinity for isopropanol and phenol, which could explain the mobility that propofol experiences inside the GABAA cavity.

Oligonucleotides as ‘bio-solvent’ for in situ extraction and functionalisation of carbon nanoparticles

We used ssDNA oligonucleotides as a ‘bio-solvent’ for CNP extraction and in situ functionalisation, developing efficient, eco-friendly, biocompatible fluorescence probes.

Effect of acceptor heteroatoms on π-hydrogen bonding interactions: a study of indole···thiophene heterodimer in a supersonic jet

Resonant two photon ionization (R2PI), IR-UV, and UV-UV double resonance spectroscopic techniques combined with quantum chemistry calculations have been used to determine the structure of indole···thiophene dimer observed in a supersonic jet. With the help of combined experimental and theoretical IR spectra it has been found that the observed dimer has a N-H···π hydrogen bonded slanted T-shaped structure. The present study demonstrates the effect of heteroatoms present in the acceptors on the strength of the π-hydrogen bonding interactions. It was concluded by Sherrill and co-workers from their theoretical study of benzene···pyridine dimer that aromatic rings containing heteroatoms are poorest π-hydrogen bond acceptors [E. G. Hohenstein and C. D. Sherrill, J. Phys. Chem. A 113, 878 (2009)]. But the current spectroscopic investigation exhibits that five membered aromatic heterocycles are favorable π-hydrogen bond acceptors. In this study, it has also been shown that thiophene is a better π-hydrogen bond acceptor than furan. The present work has immense biological significance as indole is the chromophore of tryptophan residue in the proteins and thiophene derivatives have potential therapeutic applications. Thus, understanding the binding motif between indole and thiophene in the heterodimer studied in this work may help in designing efficient drugs.

Thiourea-catalysed ring opening of episulfonium ions with indole derivatives by means of stabilizing non-covalent interactions

Small organic and metal-containing molecules (molecular mass <1,000) can catalyse synthetically useful reactions with the high levels of stereoselectivity typically associated with macromolecular enzymatic catalysts. Whereas enzymes are generally understood to accelerate reactions and impart selectivity as they stabilize specific transition structures through networks of cooperative interactions, enantioselectivity with chiral, small-molecule catalysts is rationalized typically by the steric destabilization of all but one dominant pathway. However, it is increasingly apparent that stabilizing effects also play an important role in small-molecule catalysis, although the mechanistic characterization of such systems is rare. Here, we show that arylpyrrolidino amido thiourea catalysts catalyse the enantioselective nucleophilic ring opening of episulfonium ions by indoles. Evidence is provided for the selective transition-state stabilization of the major pathway by the thiourea catalyst in the rate- and selectivity-determining step. Enantioselectivity is achieved through a network of attractive anion binding, cation-π and hydrogen-bond interactions between the catalyst and the reacting components in the transition-structure assembly.

Substituent effects on non-covalent interactions with aromatic rings: insights from computational chemistry

Non-covalent interactions with aromatic rings pervade modern chemical research. The strength and orientation of these interactions can be tuned and controlled through substituent effects. Computational studies of model complexes have provided a detailed understanding of the origin and nature of these substituent effects, and pinpointed flaws in entrenched models of these interactions in the literature. Here, we provide a brief review of efforts over the last decade to unravel the origin of substituent effects in π-stacking, XH/π, and ion/π interactions through detailed computational studies. We highlight recent progress that has been made, while also uncovering areas where future studies are warranted.

Accurate Interaction Energies for Problematic Dispersion-Bound Complexes: Homogeneous Dimers of NCCN, P2, and PCCP

All intermolecular interactions involve London dispersion forces. The accurate treatment of dispersion is essential for the computation of realistic interaction potentials. In general, the most reliable method for computing intermolecular interactions is coupled-cluster singles and doubles with perturbative triples [CCSD(T)] in conjunction with a sufficiently flexible Gaussian atomic orbital basis set, a combination which is not routinely applicable due to its excessive computational demands (CPU time, memory, storage). Recently, many theoretical methods have been developed that attempt to account for dispersion in a more efficient manner. It is well-known that dispersion interactions are more difficult to compute in some systems than others; for example, π-π dispersion is notoriously difficult, while alkane-alkane dispersion is relatively simple to compute. In this work, numerous theoretical methods are tested for their ability to compute reliable interaction energies in particularly challenging systems, namely, the P2, PCCP, and NCCN dimers. Symmetry-adapted perturbation theory (SAPT) is applied to these dimers to demonstrate their sensitivity to the treatment of dispersion. Due to the small size of these systems, highly accurate CCSD(T) potential energy curves could be estimated at the complete basis set limit. Numerous theoretical methods are tested against the reliable CCSD(T) benchmarks. Methods using a treatment of dispersion that relies on time-dependent density functional theory (TDDFT) response functions are found to be the most reliable.

Aromatic rings in chemical and biological recognition: energetics and structures

This review describes a multidimensional treatment of molecular recognition phenomena involving aromatic rings in chemical and biological systems. It summarizes new results reported since the appearance of an earlier review in 2003 in host-guest chemistry, biological affinity assays and biostructural analysis, data base mining in the Cambridge Structural Database (CSD) and the Protein Data Bank (PDB), and advanced computational studies. Topics addressed are arene-arene, perfluoroarene-arene, S⋅⋅⋅aromatic, cation-π, and anion-π interactions, as well as hydrogen bonding to π systems. The generated knowledge benefits, in particular, structure-based hit-to-lead development and lead optimization both in the pharmaceutical and in the crop protection industry. It equally facilitates the development of new advanced materials and supramolecular systems, and should inspire further utilization of interactions with aromatic rings to control the stereochemical outcome of synthetic transformations.

Density fitting of intramonomer correlation effects in symmetry-adapted perturbation theory

Symmetry-adapted perturbation theory (SAPT) offers insight into the nature of intermolecular interactions. In addition, accurate energies can be obtained from the wave function-based variant of SAPT provided that intramonomer electron correlation effects are included. We apply density-fitting (DF) approximations to the intramonomer correlation corrections in SAPT. The introduction of this approximation leads to an improvement in the computational cost of SAPT by reducing the scaling of certain SAPT terms, reducing the amount of disk I/O, and avoiding the explicit computation of certain types of MO integrals. We have implemented all the intramonomer correlation corrections to SAPT through second-order under the DF approximation. Additionally, leading third-order terms are also implemented. The accuracy of this truncation of SAPT is tested against the S22 test set of Hobza and co-workers [Phys. Chem. Chem. Phys. 8, 1985 (2006)]. When the intramonomer corrections to dispersion are included in SAPT, a mean absolute deviation of 0.3-0.4 kcal mol(-1) is observed for the S22 test set when using an aug-cc-pVDZ basis. The computations on the adenine-thymine complexes in the S22 test set with an aug-cc-pVDZ basis represent the largest SAPT computations to date that include this degree of intramonomer correlation. Computations of this size can now be performed routinely with our newly developed DF-SAPT program.

Sulfur, not too far behind O, N, and C: SH...pi hydrogen bond

We report hydrogen-bonded complexes of H(2)S with indole and 3-methyl indole stabilized by the S-H...pi interaction. It is interesting to discover that although sulfur and its hydrides are known as poor hydrogen-bond donor/acceptors, sulfur is not too far behind oxygen, nitrogen, and carbon in regard to forming the pi-type hydrogen bonds. This report also extends the scope of our earlier studies from sigma-type hydrogen-bonded complexes of sulfur (O-H...S and N-H...S sigma-type hydrogen-bonded complexes) to pi-type hydrogen-bonded complexes of sulfur (S-H...pi pi-type hydrogen-bonded complexes). The experiments were carried out using the supersonic jet expansion technique, and the complexes were probed using laser-induced spectroscopy such as laser-induced fluorescence (LIF), resonant two-photon inonization (R2PI), and fluorescence dip infrared spectroscopy (FDIRS). The FDIR spectroscopy revealed that while there was no shift in the N-H stretch, the S-H stretch was red shifted by about 21 cm(-1). For the H(2)O complexes of indole and 3-methylindole, however, there was a significant red shift in the N-H stretch. These observations suggest that H(2)O forms a NH...O type complex, whereas H(2)S prefers to form a SH...pi type complex. The experimental results were complemented by ab initio calculations and energy decomposition analysis. The binding energies for both the sigma-type and pi-type hydrogen-bonded M.L complexes (M = indole and 3-methylindole; L = H(2)O and H(2)S) were calculated by extrapolating MP2 interaction energies to the complete basis set limit. The calculated M.H(2)S (sigma-type) interaction energy (2.74 kcal/mol) was considerably smaller than that of the M.H(2)S pi-type hydrogen-bonded complex (4.89 kcal/mol), which is exactly opposite of the trend found for the M.H(2)O complexes. This is consistent with the experimental observations. Comparison of the S-H...pi interaction with the other type of X-H...pi (X = C, N, and O) shows that the S-H...pi interaction is the strongest among them. In all of the pi-type HB complexes, the dispersion energy component has significant contribution to the total binding energy.

Polarizable empirical force field for nitrogen-containing heteroaromatic compounds based on the classical Drude oscillator

The polarizable empirical CHARMM force field based on the classical Drude oscillator has been extended to the nitrogen-containing heteroaromatic compounds pyridine, pyrimidine, pyrrole, imidazole, indole, and purine. Initial parameters for the six-membered rings were based on benzene with nonbond parameter optimization focused on the nitrogen atoms and adjacent carbons and attached hydrogens. In the case of five-member rings, parameters were first developed for imidazole and transferred to pyrrole. Optimization of all parameters was performed against an extensive set of quantum mechanical and experimental data. Ab initio data were used for the determination of initial electrostatic parameters, the vibrational analysis, and in the optimization of the relative magnitudes of the Lennard-Jones (LJ) parameters, through computations of the interactions of dimers of model compounds, model compound-water interactions, and interactions of rare gases with model compounds. The absolute values of the LJ parameters were determined targeting experimental heats of vaporization, molecular volumes, heats of sublimation, crystal lattice parameters, and free energies of hydration. Final scaling of the polarizabilities from the gas-phase values by 0.85 was determined by reproduction of the dielectric constants of pyridine and pyrrole. The developed parameter set was extensively validated against additional experimental data such as diffusion constants, heat capacities, and isothermal compressibilities, including data as a function of temperature.

Specific and nonspecific interactions in a molecule with flexible side chain: 2-phenylethanol and its 1:1 complex with argon studied by high-resolution UV spectroscopy

Using high-resolution resonance-enhanced two-photon ionization spectroscopy in combination with genetic-algorithm-based computer-aided rotational fit analysis and ab initio quantum chemistry calculations we determined the conformational structure and transition moment orientation in 2-phenylethanol and its 1:1 clusters with argon. The results clearly demonstrate that the gauche structure of 2-phenylethanol, which is stabilized by the intramolecular pi-hydrogen bond between the folded side chain and the benzene ring, is the most abundant in the cold molecular beam. In this conformer the transition moment is rotated by 18 degrees from the short axis of the aromatic ring. Two distinct 1:1 complexes of 2-phenylethanol with argon in a cis- and trans-configuration with respect to the side chain have been found. Employing the Kraitchman [Am. J. Phys. 21, 17 (1953)] analysis we have found that the structure of the 2-phenylethanol moiety and the orientation of the transition moment do not change after the complexation with argon within the experimental accuracy. From the measured band intensities we conclude that in addition to the dispersion interaction of the argon atom with the aromatic ring a hydrogen-bond-type interaction with the terminal -OH group of the side chain stabilizes the cis-structure of the 1:1 complex of 2-phenylethanol with argon.

The World of Non-Covalent Interactions: 2006

The review focusses on the fundamental importance of non-covalent interactions in nature by illustrating specific examples from chemistry, physics and the biosciences. Laser spectroscopic methods and both ab initio and molecular modelling procedures used for the study of non-covalent interactions in molecular clusters are briefly outlined. The role of structure and geometry, stabilization energy, potential and free energy surfaces for molecular clusters is extensively discussed in the light of the most advanced ab initio computational results for the CCSD(T) method, extrapolated to the CBS limit. The most important types of non-covalent complexes are classified and several small and medium size non-covalent systems, including H-bonded and improper H-bonded complexes, nucleic acid base pairs, and peptides and proteins are discussed with some detail. Finally, we evaluate the interpretation of experimental results in comparison with state of the art theoretical models: this is illustrated for phenol...Ar, the benzene dimer and nucleic acid base pairs. A review with 270 references.

Structure and Binding Energy of Anion−π and Cation−π Complexes: A Comparison of MP2, RI-MP2, DFT, and DF-DFT Methods

Several complexes of benzene with cations, hexafluorobenzene with anions, 1,3,5-trifluorobenzene with cations and anions, and s-triazine with cations and anions have been evaluated and compared at the MP2 and resolution of the identity MP2 (RI-MP2) levels. The RI-MP2 method is considerably faster than the MP2 and the interaction energies and equilibrium distances are almost identical for both methods. A similar result is found when comparing DFT and density fitting DFT (DF-DFT) levels. Therefore RI-MP2 and DF-DFT methods are well suited for the study of ion-pi interactions.

Accurate description of van der Waals complexes by density functional theory including empirical corrections

An empirical method to account for van der Waals interactions in practical calculations with the density functional theory (termed DFT-D) is tested for a wide variety of molecular complexes. As in previous schemes, the dispersive energy is described by damped interatomic potentials of the form C6R(-6). The use of pure, gradient-corrected density functionals (BLYP and PBE), together with the resolution-of-the-identity (RI) approximation for the Coulomb operator, allows very efficient computations for large systems. Opposed to previous work, extended AO basis sets of polarized TZV or QZV quality are employed, which reduces the basis set superposition error to a negligible extend. By using a global scaling factor for the atomic C6 coefficients, the functional dependence of the results could be strongly reduced. The "double counting" of correlation effects for strongly bound complexes is found to be insignificant if steep damping functions are employed. The method is applied to a total of 29 complexes of atoms and small molecules (Ne, CH4, NH3, H2O, CH3F, N2, F2, formic acid, ethene, and ethine) with each other and with benzene, to benzene, naphthalene, pyrene, and coronene dimers, the naphthalene trimer, coronene. H2O and four H-bonded and stacked DNA base pairs (AT and GC). In almost all cases, very good agreement with reliable theoretical or experimental results for binding energies and intermolecular distances is obtained. For stacked aromatic systems and the important base pairs, the DFT-D-BLYP model seems to be even superior to standard MP2 treatments that systematically overbind. The good results obtained suggest the approach as a practical tool to describe the properties of many important van der Waals systems in chemistry. Furthermore, the DFT-D data may either be used to calibrate much simpler (e.g., force-field) potentials or the optimized structures can be used as input for more accurate ab initio calculations of the interaction energies.

True stabilization energies for the optimal planar hydrogen-bonded and stacked structures of guanine...cytosine, adenine...thymine, and their 9- and 1-methyl derivatives: complete basis set calculations at the MP2 and CCSD(T) levels and comparison with experiment

Planar H-bonded and stacked structures of guanine...cytosine (G.C), adenine...thymine (A...T), 9-methylguanine...1-methylcytosine (mG...mC), and 9-methyladenine...1-methylthymine (mA...mT) were optimized at the RI-MP2 level using the TZVPP ([5s3p2d1f/3s2p1d]) basis set. Planar H-bonded structures of G...C, mG...mC, and A...T correspond to the Watson-Crick (WC) arrangement, in contrast to mA...mT for which the Hoogsteen (H) structure is found. Stabilization energies for all structures were determined as the sum of the complete basis set limit of MP2 energies and a (DeltaE(CCSD(T)) - DeltaE(MP2)) correction term evaluated with the cc-pVDZ(0.25,0.15) basis set. The complete basis set limit of MP2 energies was determined by two-point extrapolation using the aug-cc-pVXZ basis sets for X = D and T and X = T and Q. This procedure is required since the convergency of the MP2 interaction energy for the present complexes is rather slow, and it is thus important to include the extrapolation to the complete basis set limit. For the MP2/aug-cc-pVQZ level of theory, stabilization energies for all complexes studied are already very close to the complete basis set limit. The much cheaper D-->T extrapolation provided a complete basis set limit close (by less than 0.7 kcal/mol) to the more accurate T-->Q term, and the D-->T extrapolation can be recommended for evaluation of complete basis set limits of more extended complexes (e.g. larger motifs of DNA). The convergency of the (DeltaE(CCSD(T)) - DeltaE(MP2)) term is known to be faster than that of the MP2 or CCSD(T) correlation energy itself, and the cc-pVDZ(0.25,0.15) basis set provides reasonable values for planar H-bonded as well as stacked structures. Inclusion of the CCSD(T) correction is essential for obtaining reliable relative values for planar H-bonding and stacking interactions; neglecting the CCSD(T) correction results in very considerable errors between 2.5 and 3.4 kcal/mol. Final stabilization energies (kcal/mol) for the base pairs studied are very substantial (A...T WC, 15.4; mA...mT H, 16.3; A...T stacked, 11.6; mA...mT stacked, 13.1; G...C WC, 28.8; mG...mC WC, 28.5; G...C stacked, 16.9; mG...mC stacked, 18.0), much larger than published previously. On the basis of comparison with experimental data, we conclude that our values represent the lower boundary of the true stabilization energies. On the basis of error analysis, we expect the present H-bonding energies to be fairly close to the true values, while stacked energies are still expected to be about 10% too low. The stacking energy for the mG...mC pair is considerably lower than the respective H-bonding energy, but it is larger than the mA...mT H-bonding energy. This conclusion could significantly change the present view on the importance of specific H-bonding interactions and nonspecific stacking interactions in nature, for instance, in DNA. Present stabilization energies for H-bonding and stacking energies represent the most accurate and reliable values and can be considered as new reference data.