SiH4与HX(X=F,Cl,Br,I)形成的二氢键复合物的结构特征及本质

本文就SiH4与HX形成的二氢键复合物的结构特征及本质进行了探讨。在MP2/6-311++G(3d,3p)水平优化、频率验证得到复合物的分子结构,通过分子间距离及电子密度等值线图,我们确认SiH4与卤化氢已形成了二氢键复合物。MP2/6-311++G(3d,3p)水平下进行BSSE校正后的结合能为2.703－4.439 KJ/mol。用对称匹配微绕理论(SAPT)对结合能进行分解,分解结果显示,SiH4匟X(X=F,Cl,Br,I)二氢键复合物中静电能对总吸引能的贡献小于28％,并且相对稳定,这就是说SiH4匟X二氢键复合物的本质并非静电作用,而是静电能、诱导能、色散能、交换能对总结合能的贡献都非常重要。     

BeH2与HX(X=F,Cl, Br,I)形成的二氢键复合物的结构特征与本质

对BeH2与HX(X=F, Cl, Br, I)形成的二氢键复合物的结构特征及本质进行了探讨. 在MP2/6-311++G(3d,3p)水平优化、频率验证, 得到复合物的分子结构, 用分子间距离及电子密度拓扑理论确认BeH2与卤化氢已形成了二氢键型复合物. 在MP2/6-311++G(3d, 3p)水平下进行基函数重叠误差(BSSE)校正后的结合能在-14.468 kJ·mol-1到-5.464 kJ·mol-1之间.用对称匹配微扰理论(SAPT)对复合物的结合能进行分解, 结果表明, BeH2…HX二氢键复合物中静电能对总吸引能的贡献都是最主要的, 但交换排斥能、诱导能、色散能对总结合能的贡献也很重要. 从BeH2…HF到BeH2…HI, 诱导能对总吸引能的贡献从37.8%逐渐减小到24.0%. 而色散能对总吸引能的贡献从BeH2…HF体系中的16.0%逐渐增加到BeH2…HI体系中的33.8%.     

Pairwise decomposition of residue interaction energies using semiempirical quantum mechanical methods in studies of protein-ligand interaction

Pairwise decomposition of the interaction energy between molecules is shown to be a powerful tool that can increase our understanding of macromolecular recognition processes. Herein we calculate the pairwise decomposition of the interaction energy between the protein human carbonic anhydrase II (HCAII) and the fluorine-substituted ligand N-(4-sulfamylbenzoyl)benzylamine (SBB) using semiempirical quantum mechanics based methods. We dissect the interaction between the ligand and the protein by dividing the ligand and the protein into subsystems to understand the structure-activity relationships as a result of fluorine substitution. In particular, the off-diagonal elements of the Fock matrix that is composed of the interaction between the ionic core and the valence electrons and the exchange energy between the subsystems or atoms of interest is examined in detail. Our analysis reveals that the fluorine-substituted benzylamine group of SBB does not directly affect the binding energy. Rather, we find that the strength of the interaction between Thr199 of HCAII and the sulfamylbenzoyl group of SBB affects the binding affinity between the protein and the ligand. These observations underline the importance of the sulfonamide group in binding affinity as shown by previous experiments (Maren, T. H.; Wiley: C. E. J. Med. Chem. 1968, 11, 228-232). Moreover, our calculations qualitatively agree with the structural aspects of these protein-ligand complexes as determined by X-ray crystallography.     

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.     

Interactions of CO2 with cluster models of metal–organic frameworks

The interactions between carbon dioxide and cluster models of coordinatively unsaturated metal–organic frameworks (MOFs) were studied using a variety of ab initio methods. Three metal species and three organic linkers in four structures were considered in these models as a representation of the tunable nature of MOFs and the potential multireference character of such systems. Common single-reference methods, such as MP2 and CCSD(T), were compared with multireference methods based on complete active space self-consistent field theory, going as far as multireference configuration interaction with single and double excitations (MRCISD). Special consideration is taken to avoid issues of size inconsistency in the CI results, where an alternate reference is used in the interaction energy definition. The benchmark values are used to judge the adequacy of a selection of density functionals for the current systems. Symmetry-adapted perturbation theory (SAPT) decomposition was performed to elucidate the important effects that comprise the binding interactions. The systems proved to have very limited multireference character, and MP2 values were closer to the CCSD(T) benchmark than the more difficult MRCISD results. Though the SAPT total energies prove to be relatively poor approximations to the benchmark interaction energies, they reveal (in most cases) the correct trends with respect to the choice of the metal. The SAPT energy decompositions indicate that the CO₂ binding is primarily driven by electrostatics, but induction and dispersion also provide sizable, and quite similar, attractive contributions. Importantly, the small diformate model provides a faithful representation of complexes with large aromatic linkers, both in terms of the total interaction energy and the SAPT decomposition.     

Interaction of molecular hydrogen with open transition metal centers for enhanced binding in metal-organic frameworks: a computational study

Molecular hydrogen is known to form stable, "nonclassical" sigma complexes with transition metal centers that are stabilized by donor-acceptor interactions and electrostatics. In this computational study, we establish that strong H2 sorption sites can be obtained in metal-organic frameworks by incorporating open transition metal sites on the organic linkers. Using density functional theory and energy decomposition analysis, we investigate the nature and characteristics of the H2 interaction with models of exposed open metal binding sites {half-sandwich piano-stool shaped complexes of the form (Arene)ML(3- n)(H2)n [M=Cr, Mo, V(-), Mn(+); Arene = C6H5X (X=H, F, Cl, OCH3, NH2, CH3, CF3) or C6H3Y2X (Y=COOH, X=CF3, Cl; L=CO; n=1-3]}. The metal-H2 bond dissociation energy of the studied complexes is calculated to be between 48 and 84 kJ/mol, based on the introduction of arene substituents, changes to the metal core, and of charge-balancing ligands. Thus, design of the binding site controls the H2 binding affinity and could be potentially used to control the magnitude of the H2 interaction energy to achieve reversible sorption characteristics at ambient conditions. Energy decomposition analysis illuminates both the possibilities and present challenges associated with rational materials design.     

The influence of CH…π interaction on hydrogen bonding ability of –CONH<Subscript>2</Subscript> functional group of benzamide

Interplay between CH…π and hydrogen bond interactions of benzamide has been investigated by quantum mechanical calculations. The effect of the substituents on geometrical parameters has also been studied at the B3LYP level with 6-311++G(d,p) basis set. The electron-withdrawing substituents enhance the total interaction energy of the complexes. The results indicated that the cooperativity of interactions leads to extra stability of the ternary complexes. The CH…π interaction and the hydrogen bond energies have been estimated using the electron densities calculated by the atoms in molecules (AIM) method at hydrogen bond critical points. The strength of hydrogen bonding increases in the presence of CH…π interaction in the ternary complexes. The effect of CH…π interaction on the hydrogen bond interaction has also been studied by the natural bond orbital, AIM and the molecular electrostatic potential analyses.     

Structures and interaction energies of stacked graphene-nucleobase complexes

The noncovalent interactions of nucleobases and hydrogen-bonded (Watson-Crick) base-pairs on graphene are investigated with the DFT-D method, i.e., all-electron density functional theory (DFT) in generalized gradient approximation (GGA) combined with an empirical correction for dispersion (van der Waals) interactions. Full geometry optimization is performed for complexes with graphene sheet models of increasing size (up to C(150)H(30)). Large Gaussian basis sets of at least polarized triple-zeta quality are employed. The interaction energies are extrapolated to infinite lateral size of the sheets. Comparisons are made with B2PLYP-D and SCS-MP2 single point energies for coronene and C(54)H(18) substrates. The contributions to the binding (Pauli exchange repulsion, electrostatic and induction, dispersion) are analyzed. At a frozen inter-fragment distance, the interaction energy surface of the rigid C(96)H(24) and base monomers is explored in three dimensions (two lateral and one rotational). Methodologically and also regarding the results of an energy decomposition analysis, the complexes behave like other pi-stacked systems examined previously. The sequence obtained for the interaction energy of bases with graphene (G > A > T > C > U) is the same for all methods and supports recent experimental findings. The absolute values are rather large (about -20 to -25 kcal mol(-1)) but in the expected range for pi-systems of that size. The relatively short equilibrium inter-plane distance (about 3 A) is consistent with atomic force microscopy results of monolayer guanine and adenine on graphite. In graphene ... Watson-Crick pair complexes, the bases lie differently from their isolated energy minima leading to geometrical anti-cooperativity. Together with an electronic contribution of 2 and 6%, this adds up to total binding anti-cooperativities of 7 and 12% for AT and CG, respectively, on C(96)H(24). Hydrogen bonds themselves are merely affected by binding on graphene.     

Rational in silico design of aptamers for organophosphates based on the example of paraoxon

Poisoning by organophosphates (OPs) takes one of the leading places in the total number of exotoxicoses. Detoxication of OPs at the first stage of the poison entering the body could be achieved with the help of DNA- or RNA-aptamers, which are able to bind poisons in the bloodstream. The aim of the research was to develop an approach to rational in silico design of aptamers for OPs based on the example of paraoxon. From the published sequence of an aptamer binding organophosphorus pesticides, its three-dimensional model has been constructed. The most probable binding site for paraoxon was determined by molecular docking and molecular dynamics (MD) methods. Then the nucleotides of the binding site were mutated consequently and the values of free binding energy have been calculated using MD trajectories and MM-PBSA approach. On the basis of the energy values, two sequences that bind paraoxon most efficiently have been selected. The value of free binding energy of paraoxon with peripheral anionic site of acetylcholinesterase (AChE) has been calculated as well. It has been revealed that the aptamers found bind paraoxon more effectively than AChE. The peculiarities of paraoxon interaction with the aptamers nucleotides have been analyzed. The possibility of improving in silico approach for aptamer selection is discussed.     

Elucidation of fatty acid amide hydrolase inhibition by potent alpha-ketoheterocycle derivatives from Monte Carlo simulations

Fatty acid amide hydrolase (FAAH) is a serine hydrolase responsible for the degradation of anandamide, an endogenous cannabinoid agonist, and oleamide, a sleep-inducing lipid. Recently, Boger and co-workers reported a potent, selective, and efficacious class of reversible alpha-ketoheterocycle inhibitors of FAAH that produce analgesia in animal models (J. Med. Chem. 2005, 48, 1849-1856; Bioorg. Med. Chem. Lett. 2005, 15, 1423-1428). Key aspects of the structure-activity data are addressed here through computational analysis of FAAH inhibition using Monte Carlo (MC) simulations in conjunction with free energy perturbation (FEP) calculations. The MC/FEP simulations demonstrate that incorporation of pyridine at the C5 position of the 2-keto-oxazole and 2-keto-1,3,4-oxadiazole derivatives significantly enhances binding affinity by formation of a hydrogen-bonded array between the pyridyl nitrogen and Lys142 and Thr236. The results also attribute the activity boost upon substitution of oxazole by oxadiazole to reduced steric interactions in the active site and a lower torsional energy penalty upon binding.     

Influence of the water molecule on cation-pi interaction: ab initio second order Møller-Plesset perturbation theory (MP2) calculations

The influence of introducing water molecules into a cation-pi complex on the interaction between the cation and the pi system was investigated using the MP2/6-311++G method to explore how a cation-pi complex changes in terms of both its geometry and its binding strength during the hydration. The calculation on the methylammonium-benzene complex showed that the cation-pi interaction is weakened by introducing H(2)O molecules into the system. For example, the optimized interaction distance between the cation and the benzene becomes longer and longer, the transferred charge between them becomes less and less, and the cation-pi binding strength becomes weaker and weaker as the water molecule is introduced one by one. Furthermore, the introduction of the third water molecule leads to a dramatic change in both the complex geometry and the binding energy, resulting in the destruction of the cation-pi interaction. The decomposition on the binding energy shows that the influence is mostly brought out through the electrostatic and induction interactions. This study also demonstrated that the basis set superposition error, thermal energy, and zero-point vibrational energy are significant and needed to be corrected for accurately predicting the binding strength in a hydrated cation-pi complex at the MP2/6-311++G level. Therefore, the results are helpful to better understand the role of water molecules in some biological processes involving cation-pi interactions.     

Estimates of the ab initio limit for sulfur-pi interactions: the H2S-benzene dimer

The interaction between aromatic rings and sulfur atoms in the side chains of amino acids is a factor in the formation and stabilization of alpha-helices in proteins. We studied the H(2)S-benzene dimer as the simplest possible prototype of sulfur-pi interactions. High-quality potential energy curves were obtained using coupled-cluster theory with single, double, and perturbative triple substitutions (CCSD(T)) and a large, augmented quadruple-zeta basis set (aug-cc-pVQZ). The equilibrium intermonomer distance for the hydrogens-down C(2)(v) configuration is 3.8 A with an interaction energy of -2.74 kcal mol(-1). Extrapolating the binding energy to the complete basis set limit gives -2.81 kcal mol(-1). This binding energy is comparable to that of H(2)O-benzene or of the benzene dimer, and the equilibrium distance is in close agreement with experiment. Other orientations of the dimer were also considered at less complete levels of theory. A considerable reduction in binding for the sulfur-down configuration, together with an energy decomposition analysis, indicates that the attraction in H(2)S-benzene is best thought of as arising from a favorable electrostatic interaction between partially positive hydrogens in H(2)S with the negatively charged pi-cloud of the benzene.     

Comprehensive study of threonine adsorption on carbon nanotube: A dispersion complemented density functional theory‐based treatment

The interaction mechanism of threonine (Thr) on the sidewall of (8, 8) single‐walled carbon nanotubes (CNTs) was investigated by density functional tight‐binding method. All the functional groups of Thr were used to interact with the surface of CNT. The structural parameters were analyzed to identify the noncovalent interactions, and the binding energy and strain energy were used to indicate the binding properties. We found that the CH/π interactions play more important roles than NH/π and OH/π interactions in stabilizing the complex structures. Furtherly, the charge transfer properties, density of states (DOS) and partial density of states, and highest occupied molecular orbitals and lowest unoccupied molecular orbitals were also studied to illustrate the adsorbed interactions. The results show that the DOS structure of CNT could be modified by the adsorption of Thr, and, therefore, the conductivity of CNT will be improved by introducing proper amino acids. Our data should be helpful for the design of biocompatible molecules for CNT modification. © 2015 Wiley Periodicals, Inc.     

Towards a force field based on density fitting

Total intermolecular interaction energies are determined with a first version of the Gaussian electrostatic model (GEM-0), a force field based on a density fitting approach using s-type Gaussian functions. The total interaction energy is computed in the spirit of the sum of interacting fragment ab initio (SIBFA) force field by separately evaluating each one of its components: electrostatic (Coulomb), exchange repulsion, polarization, and charge transfer intermolecular interaction energies, in order to reproduce reference constrained space orbital variation (CSOV) energy decomposition calculations at the B3LYP/aug-cc-pVTZ level. The use of an auxiliary basis set restricted to spherical Gaussian functions facilitates the rotation of the fitted densities of rigid fragments and enables a fast and accurate density fitting evaluation of Coulomb and exchange-repulsion energy, the latter using the overlap model introduced by Wheatley and Price [Mol. Phys. 69, 50718 (1990)]. The SIBFA energy scheme for polarization and charge transfer has been implemented using the electric fields and electrostatic potentials generated by the fitted densities. GEM-0 has been tested on ten stationary points of the water dimer potential energy surface and on three water clusters (n = 16,20,64). The results show very good agreement with density functional theory calculations, reproducing the individual CSOV energy contributions for a given interaction as well as the B3LYP total interaction energies with errors below kBT at room temperature. Preliminary results for Coulomb and exchange-repulsion energies of metal cation complexes and coupled cluster singles doubles electron densities are discussed.     

Activation mechanism of a signaling protein at atomic resolution from advanced computations

Advanced computational techniques including transition path sampling and free energy calculations are combined synergistically to reveal the activation mechanism at unprecedented resolution for a small signaling protein, chemotaxis protein Y. In the conventional "Y-T coupling" model for response regulators, phosphorylation induces the displacement of the conserved Thr87 residue through hydrogen-bond formation, which in turn makes it sterically possible for Tyr106 to isomerize from a solvent exposed configuration to a buried rotameric state. More than 160 unbiased activation trajectories show, however, that the rotation of Tyr106 does not rely on the displacement of Thr87 per se. Free energy calculations reveal that the Tyr106 rotation is a low-barrier process in the absence of the Thr87-phosphate hydrogen bond, although the rotation is stabilized by the formation of this interaction. The simulations also find that structural change in the beta4-alpha4 loop does not gate the Tyr106 rotation as suggested previously; rather, the rotation of Tyr106 stabilizes the activated configuration of this loop. The computational strategy used and mechanistic insights obtained have an impact on the study of signaling proteins and allosteric systems in general.     

A quantum chemical study of the interactions of uracil as a constituent of ribonucleic acid (RNA) with thiazolidinedione and rhodanine bioactive molecules: an insight into energetic and structural features

In this study, the biologically active configurations composed of Thiazolidinedione–Uracil (TU) and Rhodanine–Uracil (RU) have been fully investigated from the energetic and structural points of view, employing B3LYP and M062X functionals in combination with the different basis sets. Dispersion corrections to the interaction energy using M062X–GD3 and double hybrid density functionals (B2PLYP–GD2, B2PLYP–GD3 and mPW2PLYP–GD2) are also taking into account. The basis set superposition error-corrected interaction energy for hydrogen bonded configurations ranges from ??5.27 to ??13.53 and ??5.25 to ??12.93 kcal/mol for TU and RU complexes respectively as calculated at M062X/6–311++G(df,pd) level. The charge transfer process within all of the TU and RU configurations were analyzed using Natural Bond Orbital (NBO) calculations. The nature of the interactions is analyzed with NBO and Atoms in Molecules (AIM) analysis at M062X/6–311++G(df,pd) and energy decomposition analysis at BP86–D3/TZ2P(ZORA)//M062X/6–311++G(df,pd) level of theory. The results confirm that the nature of the interactions is nearly electrostatic, with a contribution of about 51–56% of the total interaction energy. The orbital interactions (ΔE_orb) for the considered TU and RU complexes have a contribution of about 24–38% of the total interaction energy. Based on the AIM and NBO results, the interactions were defined as electrostatic H-bonds with partially covalent character. In addition, correlation between interaction energies and vibrational frequency changes was investigated.     

Electron density and energy decomposition analysis in hydrogen-bonded complexes of azabenzenes with water, acetamide, and thioacetamide

Ab initio and density functional theoretical studies on hydrogen-bonded complexes of azabenzenes with water, acetamide, and thioacetamide have been carried out to explore the controversy involved in the relative order of their stability in a systematic way. The interaction energies of these complexes have been analyzed using the Morokuma energy decomposition method, and the nature of the various hydrogen bonds formed has been investigated through topological aspects using Bader's atom in a molecule (AIM) theory. Morokuma energy decomposition analysis reveals that the major contributions to the energetics are from the polarization (PL) and charge transfer (CT) energies. From the calculated topological results, excellent linear correlation is shown to exist between the hydrogen-bond length, electron density [rho(r)], and its Laplacian [nabla(2)rho(r)] at the bond critical points for all the complexes considered.     

硝仿肼离子对相互作用的密度泛函理论研究

用密度泛函理论(DFT)方法,在B3LYP/6-31+G^**水平下,求得硝仿肼离子对体系势能面上2种全优化构型.经基组叠加误差(BSSE)和零点能校正,求得离子对最大相互作用能为-420.03kJ/mol,肼和硝仿离子化所需能量可由该值得到完全补偿.离子对间键的主要贡献为库仑作用,但键鞍点上的电子密度值表明共价作用也有显著的贡献.基于统计热力学求得相关体系的热力学性质,298.2K时由自由离子形成最稳定离子对的最大焓变和最大自由能变化分别为-419.72和-376.61kJ/mol     

乙酰胆碱酯酶AChE与1,7-二氮杂咔唑抑制剂的作用机理的分子动力学模拟

通过分子对接、分子动力学(MD)模拟以及成键自由能分析方法,从原子水平上模拟研究了3种1,7-二氮杂咔唑衍生物(分别记为M1、M2和M3)与ACh E的结合模式及相互作用机理,分析和讨论了研究体系的静电相互作用和范德华相互作用(vd W)。用MM-PBSA方法计算的3种抑制剂与ACh E之间的结合自由能与抑制剂的实验生物活性数据(IC50值)相对应。分析结果表明,残基S286与抑制剂之间形成的氢键作用有利于抑制剂与ACh E之间的结合。范德华相互作用,尤其是抑制剂与关键残基W279和Y334的作用,对抑制剂与ACh E之间的结合自由能有较大的贡献,在区分抑制剂M1(或M2)和M3的生物活性上发挥着重要的作用。