How well can new-generation density functional methods describe stacking interactions in biological systems?

QSAR Studies, Molecular Docking, Molecular Dynamics, Synthesis, and Biological Evaluation of Novel Quinolinone-Based Thiosemicarbazones against Mycobacterium tuberculosis

In this study, a series of novel quinolinone-based thiosemicarbazones were designed in silico and their activities tested in vitro against Mycobacterium tuberculosis (M. tuberculosis). Quantitative structure-activity relationship (QSAR) studies were performed using quinolinone and thiosemicarbazide as pharmacophoric nuclei; the best model showed statistical parameters of R² = 0.83; F = 47.96; s = 0.31, and was validated by several different methods. The van der Waals volume, electron density, and electronegativity model results suggested a pivotal role in antituberculosis (anti-TB) activity. Subsequently, from this model a new series of quinolinone-thiosemicarbazone 11a-e was designed and docked against two tuberculosis protein targets: enoyl-acyl carrier protein reductase (InhA) and decaprenylphosphoryl-β-D-ribose-2'-oxidase (DprE1). Molecular dynamics simulation over 200 ns showed a binding energy of -71.3 to -12.7 Kcal/mol, suggesting likely inhibition. In vitro antimycobacterial activity of quinolinone-thiosemicarbazone for 11a-e was evaluated against M. bovis, M. tuberculosis H37Rv, and six different strains of drug-resistant M. tuberculosis. All compounds exhibited good to excellent activity against all the families of M. tuberculosis. Several of the here synthesized compounds were more effective than the standard drugs (isoniazid, oxafloxacin), 11d and 11e being the most active products. The results suggest that these compounds may contribute as lead compounds in the research of new potential antimycobacterial agents.

Antiviral activities of natural compounds and ionic liquids to inhibit the Mpro of SARS-CoV-2: a computational approach

The recalcitrant spread of the COVID-19 pandemic produced by the novel coronavirus SARS-CoV-2 is one of the most destructive occurrences in history. Despite the availability of several effective vaccinations and their widespread use, this line of immunization often faces questions about its long-term efficacy. Since coronaviruses rapidly change, and multiple SARS-CoV-2 variants have emerged around the world. Therefore, finding a new target-based medication became a priority to prevent and control COVID-19 infections. The main protease (Mpro) is a salient enzyme in coronaviruses that plays a vital role in viral replication, making it a fascinating therapeutic target for SARS-CoV-2. We screened 0.2 million natural products against the Mpro of SARS-CoV-2 using the Universal Natural Product Database (UNPD). As well, we studied the role of ionic liquids (ILs) on the structural stabilization of Mpro. Cholinium-based ILs are biocompatible and used for a variety of biomedical applications. Molecular docking was employed for the initial screening of natural products and ILs against Mpro. To predict the drug-likeness features of lead compounds, we calculated the ADMET properties. We performed MD simulations for the selected complexes based on the docking outcomes. Using MM/PBSA approaches, we conclude that compounds NP-Hit2 (-25.6 kcal mol-1) and NP-Hit3 (-25.3 kcal mol-1) show stronger binding affinity with Mpro. The hotspot residues of Thr25, Leu27, His41, Met49, Cys145, Met165, and Gln189 strongly interacted with the natural compounds. Furthermore, naproxenate, ketoprofenate, and geranate, cholinium-based ILs strongly interact with Mpro and these ILs have antimicrobial properties. Our findings will aid in the development of effective Mpro inhibitors.

Benchmark of the Extension of Frozen-Density Embedding Theory to Nonvariational Correlated Methods: The Embedded-MP2 Case

The extension of the frozen-density embedding theory for nonvariational methods [J. Chem. Theory Comput. 2020, 16, 6880] was utilized to evaluate intermolecular interaction energies for complexes in the Zhao-Truhlar basis set. In the applied method (FDET-MP2-FAT-LDA), the same auxiliary system is used to evaluate the correlation energy by means of the second-order Møller-Plesset perturbation theory (MP2), as in our previous work [J. Chem. Phys. 2019, 150, 121101]. Local density approximation is used for E_xcT^nad[ρ_A,ρ_B] in both cases. Additionally, the contribution to the energy due to the neglected correlation potential was evaluated and analyzed. The domain of applicability of the local density approximation for E_xcT^nad[ρ_A,ρ_B] was determined based on deviations from the interaction energies from the conventional MP2 calculations. The local density approximation for E_xcT^nad[ρ_A,ρ_B] performs well for hydrogen- or dipole-bound complexes. The relative errors in the interaction energy lie within 3-30%. While for charge-transfer complexes, this approximation fails consistently, and for other types of complexes, the performance of this approximation is not systematic. The sources of error are discussed in detail.

Multicomponent solid forms of the uric acid reabsorption inhibitor lesinurad and cocrystal polymorphs with urea: DFT simulation and solubility study

Lesinurad (systematic name: 2-{[5-bromo-4-(4-cyclopropylnaphthalen-1-yl)-4H-1,2,4-triazol-3-yl]sulfanyl}acetic acid, C₁₇H₁₄BrN₃O₂S) is a selective uric acid reabsorption inhibitor related to gout, which exhibits poor aqueous solubility. High-throughput solid-form screening was performed to screen for new solid forms with improved pharmaceutically relevant properties. During polymorph screening, we obtained two solvates with methanol (CH₃OH) and ethanol (C₂H₅OH). Binary systems with caffeine (systematic name: 3,7-dihydro-1,3,7-trimethyl-1H-purine-2,6-dione, C₈H₁₀N₄O₂) and nicotinamide (C₆H₆N₂O), polymorphs with urea (CH₄N₂O) and eutectics with similar drugs, like allopurinol and febuxostat, were prepared using the crystal engineering approach. All these novel solid forms were confirmed by XRD, DSC and FT-IR. The crystal structures were solved by single-crystal and powder X-ray diffraction. The crystal structures indicate that the lesinurad molecule is highly flexible and the triazole moiety, along with the rotatable thioacetic acid (side chain) and cyclopropane ring, is almost perpendicular to the planar naphthalene moiety. The carboxylic acid-triazole heterosynthon in the drug is interrupted by the presence of methanol and ethanol molecules in their crystal structures and forms intermolecular macrocyclic rings. The caffeine cocrystal maintains the consistency of the acid-triazole heterosynthons as in the drug and, in addition, they are bound by several auxiliary interactions. In the binary system of nicotinamide and urea, the acid-triazole heterosynthon is replaced by an acid-amide synthon. Among the urea cocrystal polymorphs, Form I (P-1, 1:1) consists of an acid-amide (urea) heterodimer, whereas in Form II (P2₁/c, 2:2), both acid-amide heterosynthons and urea-urea dimers co-exist. Density functional theory (DFT) calculations further support the experimentally observed synthon hierarchies in the cocrystals. Aqueous solubility experiments of lesinurad and its binary solids in pH 5 acetate buffer medium indicate the apparent solubility order lesinurad-urea Form I (43-fold) > lesinurad-caffeine (20-fold) > lesinurad-allopurinol (12-fold) ≃ lesinurad-nicotinamide (11-fold) > lesinurad, and this order is correlated with the crystal structures.

Reassessment of the Mechanisms of Thermal C-H Bond Activation of Methane by Cationic Magnesium Oxides: A Critical Evaluation of the Suitability of Different Density Functionals

The mechanisms of the thermal reactions of the two iconic magnesium oxide cations MgO^.+ and Mg₂ O₂ ^.+ with methane have been re-evaluated at the CCSD(T)/CBS//CCSD/def2-TZVP level of theory. For the reaction of MgO^.+ with CH₄ , only the classical hydrogen-atom transfer (HAT) was found; in contrast, for the Mg₂ O₂ ^.+ /CH₄ couple, both HAT and proton-coupled electron-transfer (PCET) exist as mechanistic variants. In order to evaluate the suitability of density functional theory (DFT) methods, the reactions were computed by using 27 density functionals. The results obtained demonstrate that the various DFT methods often deliver rather different results for both geometric and energetic features. As to the prediction of the apparent barriers, pure functionals give the largest mean absolute errors. BMK, ωB97XD, and the double-hybrid functional mPW2PLYP were confirmed to come closest to the results provided by CCSD(T)/CBS. Thus, mechanistic conclusions based on a single DFT method should be viewed with great caution. In summary, this study may assist in the selection of a suitable quantum chemical method to unravel the mechanistic details of C-H bond activation by charged metal oxides.

Theoretical Insights into the Role of Water Molecules in the Guanidinium-Based Protein Denaturation Process in Specific to Aromatic Amino Acids

Noncovalent interactions between the guanidinium cation (Gdm⁺) and aromatic amino acids (AAs) in the water molecules have been studied using quantum chemical calculation and molecular dynamics (MD) simulations. Our studies show that there are two different modes of interactions between Gdm⁺ and AAs with and without water molecules. It is observed that nonhydrated Gdm⁺ interacts with AAs through N-H···π interactions, whereas hydrated clusters of Gdm⁺ are stabilized by stacking interactions with the help of the water-mediated hydrogen bond. Thus, different hydration patterns have significant effects on the predominant cation···π interactions in AAs-Gdm⁺ complexes. Findings from MD simulation elicit that the interaction pattern of Gdm⁺ with AAs varies as Phe < Tyr < Trp. Both the QM and MD calculations show a similar trend in the interaction of AAs with Gdm⁺. Moreover, the interaction of AAs with Gdm⁺ depends on the spatial orientation of AAs in the protein and the concomitant local structure, that is, the AAs present in the unstructured region of protein such as coils and bends exhibit higher binding for Gdm⁺ when compared to the AAs present in the structured region of the protein such as the α-helix and the β-sheet. Our study clearly reveals that H-bonded water molecules and the hydration pattern of Gdm⁺ as well as the positional presence of these AAs in the protein structure context play determining roles in the denaturation of protein by the Gdm⁺ cation.

Quantum-Chemical Insights into the Self-Assembly of Carbon-Based Supramolecular Complexes

Understanding how molecular systems self-assemble to form well-organized superstructures governed by noncovalent interactions is essential in the field of supramolecular chemistry. In the nanoscience context, the self-assembly of different carbon-based nanoforms (fullerenes, carbon nanotubes and graphene) with, in general, electron-donor molecular systems, has received increasing attention as a means of generating potential candidates for technological applications. In these carbon-based systems, a deep characterization of the supramolecular organization is crucial to establish an intimate relation between supramolecular structure and functionality. Detailed structural information on the self-assembly of these carbon-based nanoforms is however not always accessible from experimental techniques. In this regard, quantum chemistry has demonstrated to be key to gain a deep insight into the supramolecular organization of molecular systems of high interest. In this review, we intend to highlight the fundamental role that quantum-chemical calculations can play to understand the supramolecular self-assembly of carbon-based nanoforms through a limited selection of supramolecular assemblies involving fullerene, fullerene fragments, nanotubes and graphene with several electron-rich π-conjugated systems.

Role of 6-Mercaptopurine in the potential therapeutic targets DNA base pairs and G-quadruplex DNA: insights from quantum chemical and molecular dynamics simulations

The theoretical studies on DNA with the anticancer drug 6-Mercaptopurine (6-MP) are investigated using theoretical methods to shed light on drug designing. Among the DNA base pairs considered, 6-MP is stacked with GC with the highest interaction energy of -46.19 kcal/mol. Structural parameters revealed that structure of the DNA base pairs is deviated from the planarity of the equilibrium position due to the formation of hydrogen bonds and stacking interactions with 6-MP. These deviations are verified through the systematic comparison between X-H bond contraction and elongation and the associated blue shift and red shift values by both NBO analysis and vibrational analysis. Bent's rule is verified for the C-H bond contraction in the 6-MP interacted base pairs. The AIM results disclose that the higher values of electron density (ρ) and Laplacian of electron density (∇²ρ) indicate the increased overlap between the orbitals that represent the strong interaction and positive values of the total electron density show the closed-shell interaction. The relative sensitivity of the chemical shift values for the DNA base pairs with 6-MP is investigated to confirm the hydrogen bond strength. Molecular dynamics simulation studies of G-quadruplex DNA d(TGGGGT)₄ with 6-MP revealed that the incorporation of 6-MP appears to cause local distortions and destabilize the G-quadruplex DNA.

Ab initio and density functional theory (DFT) studies on triflic acid with water and protonated water clusters

The structure, stability and infrared spectral signatures of triflic acid (TA) with water clusters (W_n) and protonated water clusters (TAH⁺W_n, n = 1 - 6) were computed using DFT and MP2 methods. Our calculations show that a minimum of three water molecules are necessary to stabilize the dissociated zwitterionic form of TA. It can be seen from the results that there is no significant movement of protons in smaller (n = 1 and 2) and linear (n = 1 - 6) types of water clusters. Further, the geometries of TAW_n clusters first form a neutral pair (NP) to contact ion pair (CIP), then form a solvent separated ion pair (SSIP) in a water hexamer. These findings reveal that proton transfer may take place through NP to CIP and then CIP to SSIP. The calculated binding energies (BEs) of ion pair clusters is always higher than that of NP clusters (i.e., more stable than the NP). Existing excess proton linear chain clusters transfer a proton to adjacent water molecules via a Grotthuss mechanism, whereas the same isomers in the branched motifs do not conduct protons. Examination of geometrical parameters and infrared frequencies reveals hydronium ion (H₃O⁺ also called Eigen cation) formation in both TAW_n and protonated TAW_n clusters. The stability of Eigen water clusters is three times higher than that of other non-Eigen water clusters. Our study shows clearly that formation of ion pairs in TAW_n and TAH⁺W_n clusters greatly favors proton transfer to neighboring water molecules and also enhances the stability of these complexes.

Onset of the Electronic Absorption Spectra of Isolated and π-Stacked Oligomers of 5,6-Dihydroxyindole: An Ab Initio Study of the Building Blocks of Eumelanin

Eumelanin is a naturally occurring skin pigment which is responsible for developing a suntan. The complex structure of eumelanin consists of π-stacked oligomers of various indole derivatives, such as the monomeric building block 5,6-dihydroxyindole (DHI). In this work, we present an ab initio wave-function study of the absorption behavior of DHI oligomers and of doubly and triply π-stacked species of these oligomers. We have simulated the onset of the electronic absorption spectra by employing the MP2 and the linear-response CC2 methods. Our results demonstrate the effect of an increasing degree of oligomerization of DHI and of an increasing degree of π-stacking of DHI oligomers on the onset of the absorption spectra and on the degree of red-shift toward the visible region of the spectrum. We find that π-stacking of DHI and its oligomers substantially red-shifts the onset of the absorption spectra. Our results also suggest that the optical properties of biological eumelanin cannot be simulated by considering the DHI building blocks alone, but instead the building blocks indole-semiquinone and indole-quinone have to be considered as well. This study contributes to advancing the understanding of the complex photophysics of the eumelanin biopolymer.

Microscopic investigations of site and functional selectivity of triazole for CO2 capture and catalytic applications

Ab initio and DFT studies on CO₂ interacting with different tautomers and isomers of triazole (TZ) are carried out to understand the adsorption mechanism, site selectivity and their mutual preferential attracting sites.

Equilibrium Geometries of Noncovalently Bound Intermolecular Complexes Derived from Subsystem Formulation of Density Functional Theory

The subsystem formulation of density functional theory is used to obtain equilibrium geometries and interaction energies for a representative set of noncovalently bound intermolecular complexes. The results are compared with literature benchmark data. The range of applicability of two considered approximations to the exchange-correlation- and nonadditive kinetic energy components of the total energy is determined. Local density approximation, which does not involve any empirical parameters, leads to excellent intermolecular equilibrium distances for hydrogen-bonded complexes (maximal error 0.13 Å for NH3-NH3). It is a method of choice for a wide class of weak intermolecular complexes including also dipole-bound and the ones formed by rare gas atoms or saturated hydrocarbons. The range of applicability of the chosen generalized gradient approximation, which was shown in our previous works to lead to good interaction energies in such complexes, where π-electrons are involved in the interaction, remains limited to this group because it improves neither binding energies nor equilibrium geometries in the wide class of complexes for which local density approximation is adequate. An efficient energy minimization procedure, in which optimization of the geometry and the electron density of each subsystem is made simultaneously, is proposed and tested.

Density Functionals for Noncovalent Interaction Energies of Biological Importance

Forty density functionals and one wavefunction method are assessed against a recently published database of accurate noncovalent interaction energies of biological importance. The comparison shows that two newly developed density functional theory (DFT) methods, PWB6K and M05-2X, give the best performance for this benchmark database of 22 noncovalent complexes, including both hydrogen-bonding and dispersion-dominated complexes. In contrast, the more popular B3LYP and PBEh functionals fail to describe the interactions in the dispersion-dominated complexes. The local spin density approximation and BHandH functionals give good performance for dispersion-dominated interactions at the expense of a large error for hydrogen bonding. PWB6K and M05-2X constitute a new generation of DFT methods based on simultaneously optimized exchange and correlation functionals that include kinetic energy density in both the exchange and correlation functional, and the present study confirms that they have greatly improved performance for noncovalent interactions as compared to previous DFT methods. We interpret this as being due to an improved treatment of medium-range correlation effects by the exchange-correlation functional. We recommend the PWB6K and M05-2X methods for investigating large biological systems and soft materials.

Conformational Stability and Spin States of Cobalt(II) Acetylacetonate: CASPT2 and DFT Study

Electronic structure and conformation of bis(acetylacetonate) cobalt(II), Co(acac)2, a prototypical mediator in controlled radical polymerization of olefins, is reinvestigated. The ab initio multiconfigurational CASSCF/CASPT2 method is used to resolve the doubts stemming from density functional theory results. We assign the quartet ground state for a single molecule and point at tetrahedral conformation as the preferred one. Several density functionals are tested against the ab initio calculations, and their performance is assessed. The strength of intermolecular interactions in the crystal structure composed of square-planar Co(acac)2 molecules ( Burgess , J. ; et al. Acta Crystallogr. 2000 , C56 , 649 - 650 ) is estimated to be sufficient for their planarization (suggested by Matyjaszewski , K. ; et al. Chem.-Eur. J. 2007 , 13 , 2480 - 2492 ).

Structure, Spectroscopy, and Bonding within the Zn(q+)-Imidazole(n) (q = 0, 1, 2; n = 1-4) Clusters and Implications for Zeolitic Imidazolate Frameworks and Zn-Enzymes

Using density functional theory (DFT) with dispersion correction and ab initio post Hartree-Fock methods, we treat the bonding, the structure, the stability, and the spectroscopy of the complexes between Zn(q+) and imidazole (Im), Zn(q+)Imn (where q = 0, 1 and 2; n = 1-4). These entities are subunits of zeolitic imidazolate frameworks (ZIFs) and Zn-enzymes, which possess relevant roles in industrial and biological domains, respectively. We also investigate the Imn (n = 2-4) clusters for comparison. For each species, we determine several new structures that were not found previously. Our calculations show a competition between atomic metal solvation, by either σ-type interactions or π-stacking type interaction, and proton transfer through hydrogen bonding (H-bonding) in charged species. This results in several geometrical environments around the metal. These are connected with structural properties and the functional role of Zn cation within ZIFs and Zn-enzymes. Moreover, we show that the Zn(2+)Imn subunits do not absorb in the visible domain, which may be related to the photostability of ZIFs. Our findings are important for the development of new applications of ZIFs and metalloenzymes.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

Substituent Effects in π-Stacking of Histidine on Functionalized-SWNT and Graphene

Adsorptions of histidine on the functionalized (10,0) single-walled carbon nanotube (SWNT) and graphene were investigated using density function theory methods, M05-2x and DFT-D. The results show that the binding of the histidine ring to the functionalized SWNT is weaker than that to the pristine SWNT for both singlet and triplet complexes, regardless of the electron-donating (-OH, -NH₂) or electron-withdrawing (-COOH) character and their attached sites. The present decreased binding is opposite to the well-known enhanced binding in the substituted benzene dimers. Since the atoms of the histidine are distant from the substituent atoms by over 6Å, there would be no direct interaction between histidine and the substituent as in the case of the substituted benzene systems. The decreased binding can be mainly driven by the aromaticity of the functionalized SWNT. The nucleus-independent chemical shift (NICS) index analysis for the functionalized SWNTs in deed shows that local aromaticity of SWNT is decreased because of the electron redistribution induced by functional groups, and the π-π stacking between the histidine ring and functionalized-SWNT is therefore decreased as compared to the pristine SWNT. However, the above trend does not remain for the binding between the histidine and graphene. The binding of the histidine to the functionalized graphene with -OH and -NH₂ is just slightly weaker than that to the pristine graphene, while its binding to COOH-SWNT becomes a little bit stronger.

Subsystem-DFT potential-energy curves for weakly interacting systems

Subsystem DFT for weakly interacting systems: good at reference geometries, but not necessarily good for potential-energy curves.

Interaction of collagen like peptides with gold nanosurfaces: a molecular dynamics investigation

In this study, an attempt has been made to understand the interaction between collagen like peptides (CPs) with a gold nanosurface (AuNS) using a classical molecular dynamics simulation.

Binding to gold nanoclusters alters the hydrogen bonding interactions and electronic properties of canonical and size-expanded DNA base pairs

Structural perturbations, in terms of size expansion and metal binding, lead to exciting electronic properties which can be exploited in designing novel nano-electronic devices.

Characterization of Znq+–imidazole (q = 0, 1, 2) organometallic complexes: DFT methods vs. standard and explicitly correlated post-Hartree–Fock methods

Benchmarking DFts for the characterization of the Zn^q+–imidazole (q= 0, 1, 2) complexes.

Role of size and shape selectivity in interaction between gold nanoclusters and imidazole: a theoretical study

We present a theoretical study on the structure, stability, spectra and electronic properties of imidazole (Im) adsorbed on gold nanoclusters (Aun, n = 2, 4, 6, 8, 10, and 20). These computations were performed using various density functional theories with and without inclusion of Grimme's (D3) dispersion correction. For small clusters, we also carried out wavefunction-based ab initio (MP2 and SCS-MP2) computations for comparison. Vibrational, atoms in molecules (AIM) and natural bond orbital (NBO) analyses clearly reveal the occurrence of charge transfer (CT) through covalent (N1-Au) and noncovalent interactions that play important roles in the stability of the Im@Aun complexes with anchor assisted H-bonds (Cα-H · Au). Therefore, gold clusters can act as H-bond acceptors with biomolecules for development of new materials and applications. Our study establishes also the ability and reliability of PBE0 and M05-2X functionals compared to B3LYP and PBE for an accurate description of covalent and noncovalent interactions between Im and gold clusters since they lead to close agreement with MP2. Finally, we show that the Au8 cluster may be viewed as large enough to mimic the 3D gold surface.

Experimental and theoretical approaches to the influence of the addition of pyrene to a series of Pd and Ni NHC-based complexes: catalytic consequences

A series of Ni and Pd complexes with three different N-heterocyclic carbene (NHC)-based ligands (imidazolylidene, benzimidazolylidene and pyrene-imidazolylidene) has been prepared and fully characterized. The influence of the addition of pyrene to solutions containing these complexes is studied by means of NMR and UV/Vis spectroscopies and by cyclic voltammetry. The addition of pyrene to the pyrene-NHC-containing Pd and Ni complexes gives rise to the formation of adducts by π-π stacking interactions between pyrene and the pyrene group of the NHC ligand. This interaction causes a modification of the electronic properties of the metal, as demonstrated by cyclic voltammetric studies of the Ni-NHC complexes. Theoretical calculations support this type of π-interactions, and justify the higher interactions observed with the pyrene-NHC containing complexes. The catalytic activities of the complexes were tested in the Suzuki-Miyaura C-C coupling and in the α-arylation of ketones. The addition of pyrene as an external π-stacking additive does not affect the activities of the complexes in the Suzuki-Miyaura coupling, but this observation may be justified due to the fact that the process is heterogeneously catalyzed, as indicated by the mercury-drop test. The addition of pyrene to the catalytic α-arylation of ketones results in a decrease in the activity of the reactions catalyzed by the pyrene-imidazolylidene palladium complex, whereas the other two catalysts do not modify their activity in the presence of this π-stacking additive.

Density functional theory with van der waals corrections study of the adsorption of alkyl, alkylthiol, alkoxyl, and amino-alkyl chains on the H:Si(111) surface

Surface modification of silicon with organic monolayers tethered to the surface by different linkers is an important process in realizing future miniaturized electronic and sensor devices. Understanding the roles played by the nature of the linking group and the chain length on the adsorption structures and stabilities of these assemblies is vital to advance this technology. This paper presents a density functional theory (DFT) study of the hydrogen passivated Si(111) surface modified with alkyl chains of the general formula H:Si-(CH2)n-CH2 and H:Si-X-(CH2)n-CH3, where X = NH, O, S and n = (0, 1, 3, 5, 7, 9, 11), at half coverage. For (X)-hexane and (X)-dodecane functionalization, we also examined various coverages up to full monolayer grafting in order to validate the result of half covered surface and the linker effect on the coverage. We find that it is necessary to take into account the van der Waals interaction between the alkyl chains. The strongest binding is for the oxygen linker, followed by S, N, and C, irrespective of chain length. The result revealed that the sequence of the stability is independent of coverage; however, linkers other than carbon can shift the optimum coverage considerably and allow further packing density. For all linkers apart from sulfur, structural properties, in particular, surface-linker-chain angles, saturate to a single value once n > 3. For sulfur, we identify three regimes, namely, n = 0-3, n = 5-7, and n = 9-11, each with its own characteristic adsorption structures. Where possible, our computational results are shown to be consistent with the available experimental data and show how the fundamental structural properties of modified Si surfaces can be controlled by the choice of linking group and chain length.

Modulation of structural, energetic and electronic properties of DNA and size-expanded DNA bases upon binding to gold clusters

A strong modulation in electronic properties, indicating that such complexes have the potential to serve as scaffolds for building nano electronic devices.