Assessment of Density Functionals for Intramolecular Dispersion-Rich Interactions

Neutral Peptides in the Gas Phase: Conformation and Aggregation Issues

Combined IR and UV laser spectroscopic techniques in molecular beams merged with theoretical approaches have proven to be an ideal tool to elucidate intrinsic structural properties on a molecular level. It offers the possibility to analyze structural changes, in a controlled molecular environment, when successively adding aggregation partners. By this, it further makes these techniques a valuable starting point for a bottom-up approach in understanding the forces shaping larger molecular systems. This bottom-up approach was successfully applied to neutral amino acids starting around the 1990s. Ever since, experimental and theoretical methods developed further, and investigations could be extended to larger peptide systems. Against this background, the review gives an introduction to secondary structures and experimental methods as well as a summary on theoretical approaches. Vibrational frequencies being characteristic probes of molecular structure and interactions are especially addressed. Archetypal biologically relevant secondary structures investigated by molecular beam spectroscopy are described, and the influences of specific peptide residues on conformational preferences as well as the competition between secondary structures are discussed. Important influences like microsolvation or aggregation behavior are presented. Beyond the linear α-peptides, the main results of structural analysis on cyclic systems as well as on β- and γ-peptides are summarized. Overall, this contribution addresses current aspects of molecular beam spectroscopy on peptides and related species and provides molecular level insights into manifold issues of chemical and biochemical relevance.

Binding of Sulfoxaflor to Aplysia californica-AChBP: Computational Insights from Multiscale Approaches

Structural features and binding properties of sulfoxaflor (SFX) with Ac-AChBP, the surrogate of the insect nAChR ligand binding domain (LBD), are reported herein using various complementary molecular modeling approaches (QM, molecular docking, molecular dynamics, and QM/QM'). The different SFX stereoisomers show distinct behaviors in terms of binding and interactions with Ac-AChBP. Molecular docking and Molecular Dynamics (MD) simulations highlight the specific intermolecular contacts involved in the binding of the different SFX isomers and the relative contribution of the SFX functional groups. QM/QM' calculations provide further insights and a significant refinement of the geometric and energetic contributions of the various residues leading to a preference for the SS and RR stereoisomers. Notable differences in terms of binding interactions are pointed out for the four stereoisomers. The results point out the induced fit of the Ac-AChBP binding site according to the SFX stereoisomer. In this process, the water molecules-mediated contacts play a key role, their energetic contribution being among the most important for the various stereoisomers. In all cases, the interaction with Trp147 is the major binding component, through CH···π and π···π interactions. This study provides a rationale for the binding of SFX to insect nAChR, in particular with respect to the new class of sulfoximine-based insect nAChR competitive modulators, and points out the requirements of various levels of theory for an accurate description of ligand-receptor interactions.

Stereochemistry-dependent hydrogen bonds stabilise stacked conformations in jet-cooled cyclic dipeptides: (LD) vs. (LL) cyclo tyrosine-tyrosine

Tyrosine-containing cyclic dipeptides based on a diketopiperazine (DKP) ring are studied under jet-cooled conditions using resonance-enhanced multi-photon ionisation (REMPI), conformer-selective IR-UV double resonance vibrational spectroscopy and quantum chemical calculations. The conformational landscape of the dipeptide containing natural L tyrosine (Tyr), namely c-LTyr-LTyr strongly differs from that of its diastereomer c-LTyr-DTyr. A similar family of conformers exists in both systems, with one aromatic ring folded on the dipeptide DKP ring and the other one extended. Weak NHπ and CHπ interactions are observed, which are slightly different in c-LTyr-LTyr and c-LTyr-DTyr. These structures are identical to those of LL and LD cyclo diphenylalanine, which only differ from c-Tyr-Tyr by the absence of hydroxyl on the benzene rings. While this is the only conformation observed for c-LTyr-DTyr, c-LTyr-LTyr exhibits an additional form stabilised by the interaction of the two hydroxyls, in which the two aromatic rings are in a stacked geometry. Stereochemical effects are still visible in the radical cation, for which one structure is observed for c-LTyr-DTyr, while the spectrum of the c-LTyr-LTyr radical cation is explained in terms of two co-existing structures.

The influence of solvent on conformational properties of peptides with Aib residue-a DFT study

The conformational propensities of the Aib residue on the example of two model peptides Ac-Aib-NHMe (1) and Ac-Aib-NMe2 (2), were studied by B3LYP and M06-2X functionals, in the gas phase and in the polar solvents. To verify the reliability of selected functionals, we also performed MP2 calculations for the tested molecules in vacuum. Polarizable continuum models (PCM and SMD) were used to estimate the solvent effect. Ramachandran maps were calculated to find all energy minima. Noncovalent intramolecular interactions due to hydrogen-bonds and dipole attractions between carbonyl groups are responsible for the relative stabilities of the conformers. In order to verify the theoretical results, the available conformations of similar X-ray structures from the Cambridge Crystallographic Data Center (CCDC) were analyzed. The results of the calculations show that both derivatives with the Aib residue in the gas phase prefer structures stabilized by intramolecular N-H⋯O hydrogen bonds, i.e., C5 and C7 conformations, while polar solvent promotes helical conformation with φ, ψ values equal to +/-60°, +/-40°. In addition, in the case of molecule 2, the helical conformation is the only one available in the polar environment. This result is fully consistent with the X-ray data. Graphical abstract Effect of solvent on the Ramachandran maps of the model peptides with Aib residue.

Use of Crystal Structure Informatics for Defining the Conformational Space Needed for Predicting Crystal Structures of Pharmaceutical Molecules

Determining the range of conformations that a flexible pharmaceutical-like molecule could plausibly adopt in a crystal structure is a key to successful crystal structure prediction (CSP) studies. We aim to use conformational information from the crystal structures in the Cambridge Structural Database (CSD) to facilitate this task. The conformations produced by the CSD Conformer Generator are reduced in number by considering the underlying rotamer distributions, an analysis of changes in molecular shape, and a minimal number of molecular ab initio calculations. This method is tested for five pharmaceutical-like molecules where an extensive CSP study has already been performed. The CSD informatics-derived set of crystal structure searches generates almost all the low-energy crystal structures previously found, including all experimental structures. The workflow effectively combines information on individual torsion angles and then eliminates the combinations that are too high in energy to be found in the solid state, reducing the resources needed to cover the solid-state conformational space of a molecule. This provides insights into how the low-energy solid-state and isolated-molecule conformations are related to the properties of the individual flexible torsion angles.

Exploring short intramolecular interactions in alkylaromatic substrates

From proteins and peptides to semiconducting polymers, aliphatic chains on aromatic groups are recurring motifs in macromolecules from very diverse application fields. Fields in which molecular folding and packing determine the macroscopic physical properties that make such advanced materials appealing in the first place. Within each macromolecule, the intrinsic structure of each unit defines how it interacts with its neighbours, ultimately opening up or denying certain backbone conformations. This eventually also determines how macromolecules interact with each other. This account deals specifically with the conformational problem of many common alkylaromatic units, examining the features of an intramolecular interaction involving a side chain with as few as three methylene groups. A set of 23 model compounds featuring an intramolecular interaction between an aliphatic X-H (X = C, N, O, and S) bond and an aromatic ring was considered. Quantitative computational analysis was made possible, thanks to complete basis set extrapolated CCSD(T) calculations and NCI topological analysis, the latter of which revealed an elaborate network of dispersive and steric interactions leading to somewhat unintuitive and unexpected results, such as the higher energetic stability of certain twisted conformational isomers over those with extended side chains. Vicinal covalent effects from polarizing groups and various heteroatoms, along with the occurrence of non-dispersive phenomena, were also investigated. The conclusions drawn from the investigation include a comprehensive set of guidelines intended to aid in the prediction of the most stable conformation for this class of building blocks. Our findings affect a variety of different research fields, including the tailoring of functional materials for organic electronics and photovoltaics, with insights into a rational treatment of conformational disorder, and the study of protein- and peptide-folding preferences, putting an emphasis on peculiar interactions between the backbone and aromatic residues.

Can Cyclen Bind Alkali Metal Azides? A DFT Study as a Precursor to Synthesis

Can cyclen (1,4,7,10-tetraazacyclododecane) bind alkali metal azides? This question is addressed by studying the geometric and electronic structures of the alkali metal azide-cyclen [M(cyclen)N3] complexes using density functional theory (DFT). The effects of adding a second cyclen ring to form the sandwich alkali metal azide-cyclen [M(cyclen)2N3] complexes are also investigated. N3(-) is found to bind to a M(+) (cyclen) template to give both end-on and side-on structures. In the end-on structures, the terminal nitrogen atom of the azide group (N1) bonds to the metal as well as to a hydrogen atom of the cyclen ring through a hydrogen bond in an end-on configuration to the cyclen ring. In the side-on structures, the N3 unit is bonded (in a side-on configuration to the cyclen ring) to the metal through the terminal nitrogen atom of the azide group (N1), and through the other terminal nitrogen atom (N3) of the azide group by a hydrogen bond to a hydrogen atom of the cyclen ring. For all the alkali metals, the N3-side-on structure is lowest in energy. Addition of a second cyclen unit to [M(cyclen)N3] to form the sandwich compounds [M(cyclen)2N3] causes the bond strength between the metal and the N3 unit to decrease. It is hoped that this computational study will be a precursor to the synthesis and experimental study of these new macrocyclic compounds; structural parameters and infrared spectra were computed, which will assist future experimental work.

Multi-conformer molecules in solutions: an NMR-based DFT/MP2 conformational study of two glucopyranosides of a vitamin E model compound

Overall geometries of both glucosyl derivatives of PMC were found on the basis of their NMR spectra in CDCl₃and relatedδ_H,C/ⁿJ_HHIEF-PCM(UFF,CHCl₃)/DFT calculational results.

An Atomic Counterpoise Method for Estimating Inter- and Intramolecular Basis Set Superposition Errors

An atomic counterpoise method is proposed to calculate estimates of inter- and intramolecular basis set superposition errors. The method estimates the basis set superposition error as a sum of atomic contributions and can be applied for both independent particle and electron correlation models. It is shown that the atomic counterpoise method provides results very similar to the molecular counterpoise method for intermolecular basis set superposition errors at both the HF and MP2 levels of theory with a sequence of increasingly larger basis sets. The advantage of the atomic counterpoise method is that it can be applied with equal ease to estimate intramolecular basis set superposition errors, for which few other methods exist. The atomic counterpoise method is computationally quite efficient, requiring typically double the amount of computer time as required for calculating the uncorrected energy.

Molecular recognition of thiaclopride by Aplysia californica AChBP: new insights from a computational investigation

The binding of thiaclopride (THI), a neonicotinoid insecticide, with Aplysia californica acetylcholine binding protein (Ac-AChBP), the surrogate of the extracellular domain of insects nicotinic acetylcholine receptors, has been studied with a QM/QM' hybrid methodology using the ONIOM approach (M06-2X/6-311G(d):PM6). The contributions of Ac-AChBP key residues for THI binding are accurately quantified from a structural and energetic point of view. The importance of water mediated hydrogen-bond (H-bond) interactions involving two water molecules and Tyr55 and Ser189 residues in the vicinity of the THI nitrile group, is specially highlighted. A larger stabilization energy is obtained with the THI-Ac-AChBP complex compared to imidacloprid (IMI), the forerunner of neonicotinoid insecticides. Pairwise interaction energy calculations rationalize this result with, in particular, a significantly more important contribution of the pivotal aromatic residues Trp147 and Tyr188 with THI through CH···π/CH···O and π-π stacking interactions, respectively. These trends are confirmed through a complementary non-covalent interaction (NCI) analysis of selected THI-Ac-AChBP amino acid pairs.

Structure and conformational analysis of the anti-HIV reverse transcriptase inhibitor AZT using MP2 and DFT methods. Differences with the natural nucleoside thymidine. Simulation of the 1st phosphorylation step with ATP

A comprehensive quantum-chemical investigation of the conformational landscape of the HIV-1 reverse transcriptase inhibitor AZT (3'-azido-3'-deoxythymidine) nucleoside analogue was carried out. The whole conformational parameters (χ, γ, β, δ, ϕ, P, νmax) were analysed as well as the NBO charges. The search located at least 55 stable structures, 9 of which were by MP2 within a 1 kcal mol(-1) electronic energy range of the global minimum. Most conformers were anti or high-anti around the glycoside bond and with North sugar ring puckering angles. The distribution of all the conformers according to the ranges of stability of the characteristic torsional angles was established. The results obtained were in accordance with those found in related anti-HIV nucleoside analogues. The best conformer in the anti form corresponded to the calculated values by MP2 of χ = -126.9°, β = 176.4° and γ = 49.1°. An analysis of the lowest vibrations in conformer C1 was carried out. The first hydration shell was simulated and the structural differences with the natural nucleoside deoxythymidine (dT) were determined. The first phosphorylation step was simulated by interacting ATP with the best hydrated clusters of AZT and dT. The Na cations act as a bridge between the phosphate moieties of ATP making it easy for -P3O3 to receive the H5' proton from AZT or dT. A proton-transfer mechanism is proposed through the water molecules. When the number of the water molecules surrounding AZT is lower than 8, the first phosphorylation step of AZT can be carried out. However, the appropriate orientation of the O5'-H in dT avoids this limitation and it can be performed with large numbers of water molecules.

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.

Analysis of the conformational profiles of fenamates shows route towards novel, higher accuracy, force-fields for pharmaceuticals

The conformational barriers of the fenamates which lead to conformational polymorphism can be represented by a novel, physically motivated, model intramolecular potential suitable for extension to other pharmaceuticals.

Conformational analysis of a polyconjugated protein-binding ligand by joint quantum chemistry and polarizable molecular mechanics. Addressing the issues of anisotropy, conjugation, polarization, and multipole transferability

We investigate the conformational properties of a potent inhibitor of neuropilin-1, a protein involved in cancer processes and macular degeneration. This inhibitor consists of four aromatic/conjugated fragments: a benzimidazole, a methylbenzene, a carboxythiourea, and a benzene-linker dioxane, and these fragments are all linked together by conjugated bonds. The calculations use the SIBFA polarizable molecular mechanics procedure. Prior to docking simulations, it is essential to ensure that variations in the ligand conformational energy upon rotations around its six main-chain torsional bonds are correctly represented (as compared to high-level ab initio quantum chemistry, QC). This is done in two successive calibration stages and one validation stage. In the latter, the minima identified following independent stepwise variations of each of the six main-chain torsion angles are used as starting points for energy minimization of all the torsion angles simultaneously. Single-point QC calculations of the minimized structures are then done to compare their relative energies ΔE conf to the SIBFA ones. We compare three different methods of deriving the multipoles and polarizabilities of the central, most critical moiety of the inhibitor: carboxythiourea (CTU). The representation that gives the best agreement with QC is the one that includes the effects of the mutual polarization energy E pol between the amide and thioamide moieties. This again highlights the critical role of this contribution. The implications and perspectives of these findings are discussed.

Absorption and fluorescence signatures of 1,2,3-triazole based regioisomers: challenging compounds for TD-DFT

Photophysical signatures of 1,2,3-triazole regioisomers: a surprisingly challenging problem for TD-DFT.

Mapping of the interaction sites of galanthamine: a quantitative analysis through pairwise potentials and quantum chemistry

A quantitative analysis of the interaction sites of the anti-Alzheimer drug galanthamine with molecular probes (water and benzene molecules) representative of its surroundings in the binding site of acetylcholinesterase (AChE) has been realized through pairwise potentials calculations and quantum chemistry. This strategy allows a full and accurate exploration of the galanthamine potential energy surface of interaction. Significantly different results are obtained according to the distances of approaches between the various molecular fragments and the conformation of the galanthamine N-methyl substituent. The geometry of the most relevant complexes has then been fully optimized through MPWB1K/6-31 + G(d,p) calculations, final energies being recomputed at the LMP2/aug-cc-pVTZ(-f) level of theory. Unexpectedly, galanthamine is found to interact mainly from its hydrogen-bond donor groups. Among those, CH groups in the vicinity of the ammonium group are prominent. The trends obtained provide rationales to the predilection of the equatorial orientation of the galanthamine N-methyl substituent for binding to AChE. The analysis of the interaction energies pointed out the independence between the various interaction sites and the rigid character of galanthamine. The comparison between the cluster calculations and the crystallographic observations in galanthamine-AChE co-crystals allows the validation of the theoretical methodology. In particular, the positions of several water molecules appearing as strongly conserved in galanthamine-AChE co-crystals are predicted by the calculations. Moreover, the experimental position and orientation of lateral chains of functionally important aminoacid residues are in close agreement with the ones predicted theoretically. Our study provides relevant information for a rational drug design of galanthamine based AChE inhibitors.

Modelling zwitterions in solution: 3-fluoro-γ-aminobutyric acid (3F-GABA)

The conformations and relative stabilities of folded and extended 3-fluoro-γ-aminobutyric acid (3F-GABA) conformers were studied using explicit solvation models. Geometry optimisations in the gas phase with one or two explicit water molecules favour folded and neutral structures containing intramolecular NH···O-C hydrogen bonds. With three or five explicit water molecules zwitterionic minima are obtained, with folded structures being preferred over extended conformers. The stability of folded versus extended zwitterionic conformers increases on going from a PCM continuum solvation model to the microsolvated complexes, though extended structures become less disfavoured with the inclusion of more water molecules. Full explicit solvation was studied with a hybrid quantum-mechanical/molecular-mechanical (QM/MM) scheme and molecular dynamics simulations, including more than 6000 TIP3P water molecules. According to free energies obtained from thermodynamic integration at the PM3/MM level and corrected for B3LYP/MM total energies, the fully extended conformer is more stable than folded ones by about -4.5 kJ mol(-1). B3LYP-computed (3)J(F,H) NMR spin-spin coupling constants, averaged over PM3/MM-MD trajectories, agree best with experiment for this fully extended form, in accordance with the original NMR analysis. The seeming discrepancy between static PCM calculations and experiment noted previously is now resolved. That the inexpensive semiempirical PM3 method performs so well for this archetypical zwitterion is encouraging for further QM/MM studies of biomolecular systems.

New insights on the molecular features and electrophysiological properties of dinotefuran, imidacloprid and acetamiprid neonicotinoid insecticides

Structural features and hydrogen-bond interactions of dinotefuran (DIN), imidacoloprid (IMI) and acetamiprid (ACE) have been investigated experimentally through analyses of new crystal structures and observations in structural databases, as well as by Density Functional Theory quantum chemical calculations. Several conformations are observed experimentally in the solid state, highlighting the large flexibility of these compounds. This feature is confirmed by the theoretical calculations in the gas phase, the numerous and different energetic minima of the three neonicotinoids being located within a 10kJ/mol range. Comparisons of the observed and simulated data sheds light on the hydrogen-bond (HB) strength of the functional group at the tip of the electronegative fragment of each pharmacophore (NO(2) for DIN and IMI and CN for ACE). This effect originates in the 'push-pull' nature of these fragments and the related extensive electron delocalization. Molecular electrostatic potential calculations provide a ranking of the two fragments of the three neonicotinoid in terms of HB strength. Thus, the NO(2) group of DIN is the strongest HB acceptor of the electronegative fragment, closely followed by the cyano group of ACE. These two groups are significantly more potent than the NO(2) group of IMI. With respect to the other fragments of the three neonicotinoids, the nitrogen atom of the pyridine of IMI and ACE are stronger HB acceptors than the oxygen atom of the furanyl moiety of DIN. Finally, compared to electrophysiological studies obtained from cockroach synaptic and extrasynaptic receptors, DIN appears more effective than IMI and ACE because it strongly increases dose-dependently the ganglionic depolarisation and the currents amplitudes. These data suggest that DIN, IMI and ACE belong to two subgroups which act differently as agonists of insect nicotinic receptors.

The Accuracy of Density Functional Theory in the Description of Cation-π and π-Hydrogen Bond Interactions

Cation-π and π-hydrogen bond interactions are ubiquitous in protein folding, molecular recognition, and ligand-receptor associations. As such systems are routinely studied at the DFT level, it becomes essential to understand the underlying accuracy of the plethora of density functionals currently available for the description of these interactions. For that purpose, we carried out theoretical calculations on two small model systems (benzene-Na(+) and benzene-H2O) that represent a paradigm for those intermolecular interactions and systematically tested 46 density functionals against the results of high-level post-HF methods, ranging from MP2 to extrapolated CCSD(T)/CBS. A total of 13 basis sets were also tested to examine the convergence of the interaction energy with basis set size. The convergence was surprisingly fast, with deviations below 0.2 kcal/mol for double-ζ polarized basis sets with diffuse functions. Concerning functional benchmarking, the Truhlar group functionals were particularly well suited for the description of the π-hydrogen bond interactions. In the case of cation-π interactions, there was not a clear correlation between accuracy and functional sophistication. Despite the large number of functionals predicting interaction energies within chemical accuracy (five for π-hydrogen bond and 20 for cation-π interactions), not a single functional has shown chemical accuracy in both cases. Moreover, if we calculate the average error for these two interactions, only two density functionals resulted in an average error below 1.0 kcal/mol (M06 and HCTH, with average errors of 0.6 and 0.8 kcal/mol). The obtained results serve as a guide for future computer simulations on this kind of system.

Gas-phase folding of a two-residue model peptide chain: on the importance of an interplay between experiment and theory

In order to assess the ability of theory to describe properly the dispersive interactions that are ubiquitous in peptide and protein systems, an isolated short peptide chain has been studied using both gas-phase laser spectroscopy and quantum chemistry. The experimentally observed coexistence of an extended form and a folded form in the supersonic expansion was found to result from comparable Gibbs free energies for the two species under the high-temperature conditions (< or = 320 K) resulting from the laser desorption technique used to vaporize the molecules. These data have been compared to results obtained using a series of quantum chemistry methods, including DFT, DFT-D, and post-Hartree-Fock methods, which give rise to a wide range of relative stabilities predicted for these two forms. The experimental observation was best reproduced by an empirically dispersion-corrected functional (B97-D) and a hybrid functional with a significant Hartree-Fock exchange term (M06-2X). In contrast, the popular post-Hartree-Fock method MP2, which is often used for benchmarking these systems, had to be discarded because of a very large basis-set superposition error. The applicability of the atomic counterpoise correction (ACP) is also discussed. This work also introduces the mandatory theoretical examination of experimental abundances. DeltaH(0 K) predictions are clearly not sufficient for discussion of folding, as the conformation inversion temperature is crucial to the conformation determination and requires taking into account thermodynamical corrections (DeltaG) in order to computationally isolate the most stable conformation.