Comparison of some dispersion-corrected and traditional functionals with CCSD(T) and MP2 ab initio methods: dispersion, induction, and basis set superposition error

Novel conjugated microporous polymers for efficient tetracycline adsorption: insights from theoretical investigations

This paper presents a detailed theoretical understanding of the noncovalent interactions between antibiotics tetracycline and conjugated microporous polymer (CMP), which is important to understand the recent experimental finding of efficient removal of antibiotics by CMP materials. We show that the co-work of π-π and H-π interactions determines the final equilibrium structures, when a tetracycline molecule spontaneously adsorbs to the surface or within the pores of the CMP network at physisorption distances. The binding energies for tetracycline/CMP systems are calculated to be -0.31 ∼ -1.15 eV, demonstrating the reliability of the adsorption. The electronic structures of CMP nanostructures remain basically undamaged upon the tetracycline adsorption. The replacement of benzothiadiazole unit with S and N heteroatoms to the phenyl moiety in the linker effectively enhanced the molecular polarity of CMP molecule and increases the interaction area between tetracycline and CMP network, consequently enhancing the average binding energies notably. Our calculations provide useful theoretical guidance for design of novel carbon-based porous adsorbents with good adsorption performance to remove residual tetracycline and other antibiotics in water.

Probing the Directionality of S···O/N Chalcogen Bond and Its Interplay with Weak C-H···O/N/S Hydrogen Bond Using Molecular Electrostatic Potential

The directionality of the chalcogen bond (Ch-bond) formed by S and its interplay with other weak interactions have important chemical and biological implications. Here, dimers made of CH₃-S-X and O/N containing nucleophiles are studied and found to be stabilized by coexisting S···O/N and C-H···O/N interactions. Based on experimentally accessible electron density and molecular electrostatic potentials (MESPs), we showed that reciprocity between S···O/N and C-H···O/N interactions in the stability of cumulative molecular interaction (ΔE) was dependent on the strength of the σ-hole on S (V_s,max). Direct correlation between ΔE of dimers with V_s,max of S supports the electrostatic nature of the Ch-bond. Such interplay of the Ch-bond is necessary for its directionality in complex nucleophiles (carbonyl groups) with multiple electron-rich centers, which is explained using MESP. A correlation between the MESP minima in the π-region and the strength of the S-π interaction explains the directional selectivity of the Ch-bond.

Design of novel conjugated microporous polymers for efficient adsorptive desulfurization of small aromatic sulfur molecules

This paper presents a computational study of the adsorptive desulfurization of small aromatic sulfur compounds by conjugated microporous polymers (CMPs). The density-functional tight-binding method augmented with an R^-6 dispersion correction is employed to investigate the physisorption binding mechanism and electronic properties of the CMP-aromatic sulfur complexes. We show that the widely extended π conjugation in the CMP skeletons is favorable for the non-covalent adsorption of aromatic thiophene and dibenzothiophene via π-π, H-π, and S-π interactions. The average binding energies are calculated to be -6.2 ∼ -15.2 kcal/mol for CMP- thiophene/dibenzothiophene systems. For the dibenzothiophene molecule with larger size and more extended conjugation, it binds more than twice stronger to CMP than the thiophene molecule. We show that the replacement of quinoline unit to the phenylene group in the network linker effectively enhances the average binding capacities by around 0.8-1.8 kcal/mol. Our calculations theoretically demonstrate that CMPs materials are kind of promising candidates for the adsorptive desulfurization of small aromatic sulfur compounds. This paper provides useful theoretical guidance for design of novel carbon-based adsorbents for adsorptive desulfurization.

Stereoselective Bromoboration of Acetylene with Boron Tribromide: Preparation and Cross-Coupling Reactions of (Z)-Bromovinylboronates

The mechanism of acetylene bromoboration in neat boron tribromide was studied carefully by means of experiment and theory. Besides the syn-addition mechanism through a four-center transition state, radical and polar anti-addition mechanisms are postulated, both triggered by HBr, which is evidenced also to take part in the Z/E isomerization of the product. The proposed mechanism is well supported by ab initio calculations at the MP2/6-31+G* level with Ahlrichs' SVP all-electron basis for Br. Implicit solvation in CH₂Cl₂ has been included using the PCM and/or SMD continuum solvent models. Comparative case studies have been performed involving the B3LYP/6-31+G* with Ahlrichs' SVP for Br and MP2/Def2TZVPP levels. The mechanistic studies resulted in development of a procedure for stereoselective bromoboration of acetylene yielding E/Z mixtures of dibromo(bromovinyl)borane with the Z-isomer as a major product (up to 85%). Transformation to the corresponding pinacol and neopentyl glycol boronates and stereoselective decomposition of their E-isomer provided pure (Z)-(2-bromovinyl)boronates in 57-60% overall yield. Their reactivity in a Negishi cross-coupling reaction was tested. An example of the one-pot reaction sequence of Negishi and Suzuki-Miyaura cross-couplings for synthesis of combretastatin A4 is also presented.

Halogen Bonding: A Halogen-Centered Noncovalent Interaction Yet to Be Understood

In addition to the underlying basic concepts and early recognition of halogen bonding, this paper reviews the conflicting views that consistently appear in the area of noncovalent interactions and the ability of covalently bonded halogen atoms in molecules to participate in noncovalent interactions that contribute to packing in the solid-state. It may be relatively straightforward to identify Type-II halogen bonding between atoms using the conceptual framework of σ-hole theory, especially when the interaction is linear and is formed between the axial positive region (σ-hole) on the halogen in one monomer and a negative site on a second interacting monomer. A σ-hole is an electron density deficient region on the halogen atom X opposite to the R–X covalent bond, where R is the remainder part of the molecule. However, it is not trivial to do so when secondary interactions are involved as the directionality of the interaction is significantly affected. We show, by providing some specific examples, that halogen bonds do not always follow the strict Type-II topology, and the occurrence of Type-I and -III halogen-centered contacts in crystals is very difficult to predict. In many instances, Type-I halogen-centered contacts appear simultaneously with Type-II halogen bonds. We employed the Independent Gradient Model, a recently proposed electron density approach for probing strong and weak interactions in molecular domains, to show that this is a very useful tool in unraveling the chemistry of halogen-assisted noncovalent interactions, especially in the weak bonding regime. Wherever possible, we have attempted to connect some of these results with those reported previously. Though useful for studying interactions of reasonable strength, IUPAC’s proposed “less than the sum of the van der Waals radii” criterion should not always be assumed as a necessary and sufficient feature to reveal weakly bound interactions, since in many crystals the attractive interaction happens to occur between the midpoint of a bond, or the junction region, and a positive or negative site.

Is the Fluorine in Molecules Dispersive? Is Molecular Electrostatic Potential a Valid Property to Explore Fluorine-Centered Non-Covalent Interactions?

Can two sites of positive electrostatic potential localized on the outer surfaces of two halogen atoms (and especially fluorine) in different molecular domains attract each other to form a non-covalent engagement? The answer, perhaps counterintuitive, is yes as shown here using the electronic structures and binding energies of the interactions for a series of 22 binary complexes formed between identical or different atomic domains in similar or related halogen-substituted molecules containing fluorine. These were obtained using various computational approaches, including density functional and ab initio first-principles theories with M06-2X, RHF, MP2 and CCSD(T). The physical chemistry of non-covalent bonding interactions in these complexes was explored using both Quantum Theory of Atoms in Molecules and Symmetry Adapted Perturbation Theories. The surface reactivity of the 17 monomers was examined using the Molecular Electrostatic Surface Potential approach. We have demonstrated inter alia that the dispersion term, the significance of which is not always appreciated, which emerges either from an energy decomposition analysis, or from a correlated calculation, plays a structure-determining role, although other contributions arising from electrostatic, exchange-repulsion and polarization effects are also important. The 0.0010 a.u. isodensity envelope, often used for mapping the electrostatic potential is found to provide incorrect information about the complete nature of the surface reactive sites on some of the isolated monomers, and can lead to a misinterpretation of the results obtained.

Ligation state of nickel during C-O bond activation with monodentate phosphines

The oxidative addition of phenolic electrophiles at Ni(0) in the presence of monodentate phosphine ligands was studied with both dispersion-free and dispersion-containing DFT methods. With the popular bulky ligand PCy3, consideration of dispersion has a striking effect on the predicted ligation state of nickel during oxidative addition of aryl sulfamates. Dispersioncontaining methods such as M06L indicate a clear preference for a bis-phosphine ligated transition state (TS), while dispersion free methods like B3LYP strongly favor a monophosphine ligated TS. This discrepancy in predicted ligation state is also found with small phosphines (PMe3) in combination with some aryl electrophiles (carbamates, acetates, pivalates, chlorides), but a bis-PMe3-ligated TS is predicted regardless of dispersion for other electrophiles (sulfamates, mesylates, tosylates). DFT calculations that include dispersion also offer a possible explanation for the observed poor efficacy of P t Bu3 as a ligand in Ni-catalyzed cross-coupling reactions.

Modeling the archetype cysteine protease reaction using dispersion corrected density functional methods in ONIOM-type hybrid QM/MM calculations; the proteolytic reaction of papain

A proteolytic reaction of papain with a simple peptide model substrate N-methylacetamide has been studied. Our aim was twofold: (i) we proposed a plausible reaction mechanism with the aid of potential energy surface scans and second geometrical derivatives calculated at the stationary points, and (ii) we investigated the applicability of the dispersion corrected density functional methods in comparison with the popular hybrid generalized gradient approximations (GGA) method (B3LYP) without such a correction in the QM/MM calculations for this particular problem. In the resting state of papain the ion pair and neutral forms of the Cys-His catalytic dyad have approximately the same energy and they are separated by only a small barrier. Zero point vibrational energy correction shifted this equilibrium slightly to the neutral form. On the other hand, the electrostatic solvation free energy corrections, calculated using the Poisson-Boltzmann method for the structures sampled from molecular dynamics simulation trajectories, resulted in a more stable ion-pair form. All methods we applied predicted at least a two elementary step acylation process via a zwitterionic tetrahedral intermediate. Using dispersion corrected DFT methods the thioester S-C bond formation and the proton transfer from histidine occur in the same elementary step, although not synchronously. The proton transfer lags behind (or at least does not precede) the S-C bond formation. The predicted transition state corresponds mainly to the S-C bond formation while the proton is still on the histidine Nδ atom. In contrast, the B3LYP method using larger basis sets predicts a transition state in which the S-C bond is almost fully formed and the transition state can be mainly featured by the Nδ(histidine) to N(amid) proton transfer. Considerably lower activation energy was predicted (especially by the B3LYP method) for the next amide bond breaking elementary step of acyl-enzyme formation. Deacylation appeared to be a single elementary step process in all the methods we applied.

Computational modelling of singlet excitation energy transfer: a DFT/TD-DFT study of the ground and excited state properties of a syn bimane dimer system using non-empirically tuned range-separated functionals

Range-tuned DFT/TD-DFT improves predicted properties for the sequence of events leading to excitation energy transfer (EET) in bimanes.

Computational study on single molecular spectroscopy of tyrosin-glycine, tryptophane-glycine and glycine-tryptophane

Quantum chemistry calculations play a fundamental role in revealing the molecular structures observed in gas-phase spectroscopic measurements. The supersonic jet cooling widely used in single molecular spectroscopy experiment is a non-equilibrium process and often causes confusion on the theoretical and experimental comparison. A computational approach is proposed here to account for the effect of the non-equilibrium cooling on the experimental spectra and applied to the cases of tyrosin-glycine (YG), tryptophane-glycine (WG) and glycine-tryptophane (GW). The low energy conformers of YG, WG and GW are obtained through thorough conformational searches. The structural features and equilibrium distributions of conformations and the energy barriers for conformer conversions are then determined. Three classes of transition energy barriers, high, medium and low, are found for the conversions among conformers with distinctly different, similar and the same structural types, respectively. The final conformation populations are determined by assuming an initial temperature of about 450 K and allowing for only the conformation conversion with a low energy barrier to occur during the rapid cooling process. The results provide a natural explanation for the numbers of YG, WG and GW conformations observed experimentally. The theoretical conformation assignments are also in good agreement with the experimental IR data.

Testing the CP-correction procedure with different DFT methods on H-bonding complexes of κ-carrabiose with water molecules

Interaction of water molecules with κ-carrabiose disaccharide, within three H-bonding complexes, was investigated. Particular interest was focused on the way with which the BSSE correction has to be performed. Two strategies were used, either performing BSSE correction during or after optimization. For this aim, several DFT-functionals (hybrid GGA and hybrid meta-GGA) and 6-31 + G* basis set were considered. The results demonstrated the uselessness of including of BSSE-CP correction during optimization for all complexes. From a structural point of view, a proper H-bonding description was obtained using the PBE0 functional for all complexes. The basis set effect on the BSSE using B3LYP functional was also investigated. The reliability of B3LYP/6-31 + G** and B3LYP/6-31++G** models for the complexes involving one or two water molecules was reported while the use of B3LYP/6-311 + G** or B3LYP/6-311++G** levels was shown to be more appropriate for larger complexes equivalent to that involving three water molecules. CP-corrected interaction energies were demonstrated to be closer to CBS-4 M interaction energies than the uncorrected ones. Graphical abstract Functional and basis set effects on BSSE.

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.

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.

3c/4e -type long-bonding competes with ω-bonding in noble-gas hydrides HNgY (Ng = He, Ne, Ar, Kr, Xe, Rn; Y = F, Cl, Br, I): a NBO/NRT perspective

Novel resonance bonding for the HNgY molecule is demonstrated based on natural resonance theory analyses. Ng/Y affects the ω-bonding vs. long-bonding propensity in each of the HNgY molecules.

Theoretical study of the physisorption of organic molecules on conjugated microporous polymers: the critical role of skeleton structures on binding strength

By enlarging the size (area) in the linker or adding substituent group in the node, the binding energy between CMPs and adsorbates will increase significantly, which demonstrates the critical role of skeleton structures on binding strength.

Tuning the van der Waals Interaction of Graphene with Molecules via Doping

We use scanning tunneling microscopy to visualize and thermal desorption spectroscopy to quantitatively measure that the binding of naphthalene molecules to graphene, a case of pure van der Waals interaction, strengthens with n and weakens with p doping of graphene. Density-functional theory calculations that include the van der Waals interaction in a seamless, ab initio way accurately reproduce the observed trend in binding energies. Based on a model calculation, we propose that the van der Waals interaction is modified by changing the spatial extent of graphene's π orbitals via doping.

Exploring Strong Interactions in Proteins with Quantum Chemistry and Examples of Their Applications in Drug Design

OBJECTIVES

Three strong interactions between amino acid side chains (salt bridge, cation-π, and amide bridge) are studied that are stronger than (or comparable to) the common hydrogen bond interactions, and play important roles in protein-protein interactions.

METHODS

Quantum chemical methods MP2 and CCSD(T) are used in calculations of interaction energies and structural optimizations.

RESULTS

The energies of three types of amino acid side chain interactions in gaseous phase and in aqueous solutions are calculated using high level quantum chemical methods and basis sets. Typical examples of amino acid salt bridge, cation-π, and amide bridge interactions are analyzed, including the inhibitor design targeting neuraminidase (NA) enzyme of influenza A virus, and the ligand binding interactions in the HCV p7 ion channel. The inhibition mechanism of the M2 proton channel in the influenza A virus is analyzed based on strong amino acid interactions.

CONCLUSION

(1) The salt bridge interactions between acidic amino acids (Glu- and Asp-) and alkaline amino acids (Arg+, Lys+ and His+) are the strongest residue-residue interactions. However, this type of interaction may be weakened by solvation effects and broken by lower pH conditions. (2) The cation- interactions between protonated amino acids (Arg+, Lys+ and His+) and aromatic amino acids (Phe, Tyr, Trp and His) are 2.5 to 5-fold stronger than common hydrogen bond interactions and are less affected by the solvation environment. (3) The amide bridge interactions between the two amide-containing amino acids (Asn and Gln) are three times stronger than hydrogen bond interactions, which are less influenced by the pH of the solution. (4) Ten of the twenty natural amino acids are involved in salt bridge, or cation-, or amide bridge interactions that often play important roles in protein-protein, protein-peptide, protein-ligand, and protein-DNA interactions.

Binding energies of tyrosine kinase inhibitors: Error assessment of computational methods for imatinib and nilotinib binding

The binding energies of imatinib and nilotinib to tyrosine kinase have been determined by quantum mechanical (QM) computations, and compared with literature binding energy studies using molecular mechanics (MM). The potential errors in the computational methods include these critical factors. Errors in X-ray structures such as structural distortions and steric clashes give unrealistically high van der Waals energies, and erroneous binding energies.MM optimization gives a very different configuration to the QM optimization for nilotinib, whereas the imatinib ion gives similar configurations. Solvation energies are a major component of the overall binding energy. The QM based solvent model (PCM/SMD) gives different values from those used in the implicit PBSA solvent MM models. A major error in inhibitor—kinase binding lies in the non-polar solvation terms. Solvent transfer free energies and the required empirical solvent accessible surface area factors for nilotinib and imatinib ion to give the transfer free energies have been reverse calculated. These values differ from those used in the MM PBSA studies.An intertwined desolvation—conformational binding selectivity process is a balance of thermodynamic desolvation and intramolecular conformational kinetic control.The configurational entropies (TΔS) are minor error sources.

A theoretical benchmark study of the spectroscopic constants of the very heavy rare gas dimers

The binding energy of the superheavy dimer Uuo₂ is considerably larger than that of its lighter homologues, despite a 40% reduction due to spin-other orbit interaction.

In silico studies on the origin of selective uptake of carbon dioxide with cucurbit[7]uril amorphous material

The efficient capture and storage of flue gases is of current interest due to environmental problems. The DFT calculation demonstrates the origin of the physisorption of flue gases (CO₂, N₂and CH₄) on amorphous solid cucurbit[7]uril.

The interactions of phenylalanines in β-sheet-like structures from molecular orbital calculations using density functional theory (DFT), MP2, and CCSD(T) methods

We present density functional theory calculations designed to evaluate the importance of π-stacking interactions to the stability of in-register Phe residues within parallel β-sheets, such as amyloids. We have used a model of a parallel H-bonded tetramer of acetylPheNH2 as a model and both functionals that were specifically designed to incorporate dispersion effects (DFs), as well as, several traditional functionals which have not been so designed. None of the functionals finds a global minimum for the π-stacked conformation, although two of the DFs find this to be a local minimum. The stacked phenyls taken from the optimized geometries calculated for each functional have been evaluated using MP2 and CCSD(T) calculations for comparison. The results suggest that π-stacking does not make an important contribution to the stability of this system and (by implication) to amyloid formation.

Capping parallel β-sheets of acetyl(Ala)6NH2 with an acetyl(Ala)5ProNH2 can arrest the growth of the sheet, suggesting a potential for curtailing amyloid growth. An ONIOM and density functional theory study

We present ONIOM calculations using B3LYP/d95(d,p) as the high level and AM1 as the medium level on parallel β-sheets containing four strands of Ac-AAAAAA-NH2 capped with either Ac-AAPAAA-NH2 or Ac-AAAPAA-NH2. Because Pro can form H-bonds from only one side of the peptide linkage (that containing the C═O H-bond acceptor), only one of the two Pro-containing strands can favorably add to the sheet on each side. Surprisingly, when the sheet is capped with AAPAAA-NH2 at one edge, the interaction between the cap and sheet is slightly more stabilizing than that of another all Ala strand. Breaking down the interaction enthalpies into H-bonding and distortion energies shows the favorable interaction to be due to lower distortion energies in both the strand and the four-stranded sheet. Because another strand would be inhibited for attachment to the other side of the capping (Pro-containing) strand, we suggest the possible use of Pro residues in peptides designed to arrest the growth of many amyloids.

Multiexciton Generation in IV-VI Nanocrystals: The Role of Carrier Effective Mass, Band Mixing, and Phonon Emission

We study the role of the effective mass, band mixing, and phonon emission on multiexciton generation in IV-VI nanocrystals. A four-band k · p effective mass model, which allows for an independent variation of these parameters, is adopted to describe the electronic structure of the nanocrystals. Multiexciton generation efficiencies are calculated using a Green's function formalism, providing results that are numerically similar to impact excitation. We find that multiexciton generation efficiencies are maximized when the effective mass of the electron and hole are small and similar. Contact with recent experimental results for multiexciton generation in PbS and PbSe is made.

Comparison of some dispersion-corrected and traditional functionals as applied to peptides and conformations of cyclohexane derivatives

We compare the energetic and structural properties of fully optimized α-helical and antiparallel β-sheet polyalanines and the energetic differences between axial and equatorial conformations of three cyclohexane derivatives (methyl, fluoro, and chloro) as calculated using several functionals designed to treat dispersion (B97-D, ωB97x-D, M06, M06L, and M06-2X) with other traditional functionals not specifically parametrized to treat dispersion (B3LYP, X3LYP, and PBE1PBE) and with experimental results. Those functionals developed to treat dispersion significantly overestimate interaction enthalpies of folding for the α-helix and predict unreasonable structures that contain Ramachandran φ and ψ and C = O...N H-bonding angles that are out of the bounds of databases compiled the β-sheets. These structures are consistent with overestimation of the interaction energies. For the cyclohexanes, these functionals overestimate the stabilities of the axial conformation, especially when used with smaller basis sets. Their performance improves when the basis set is improved from D95∗∗ to aug-cc-pVTZ (which would not be possible with systems as large as the peptides).

The folding of acetyl(Ala)28NH2 and acetyl(Ala)40NH2 extended strand peptides into antiparallel β-sheets. A density functional theory study of β-sheets with β-turns

We report ONIOM calculations using B3LYP/D95** and AM1 on β-sheet formation from acetyl(Ala)(N)NH(2) (N = 28 or 40). The sheets contain from one to four β-turns for N = 28 and up to six for N = 40. We have obtained four types of geometrically optimized structures. All contain only β-turns. They differ from each other in the types of β-turns formed. The unsolvated sheets containing two turns are most stable. Aqueous solvation (using the SM5.2 and CPCM methods) reduces the stabilities of the folded structures compared to the extended strands.