首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
We examined CH/π hydrogen bonds in protein/ligand complexes involving at least one proline residue using the ab initio fragment molecular orbital (FMO) method and the program CHPI. FMO calculations were carried out at the Hartree–Fock (HF)/6‐31G*, HF/6‐31G**, second‐order Møller–Plesset perturbation (MP2)/6‐31G*, and MP2/6‐31G** levels for three Src homology 3 (SH3) domains and five proline‐recognition domains (PRDs) complexed with their corresponding ligand peptides. PRDs use a conserved set of aromatic residues to recognize proline‐rich sequences of specific ligands. Many CH/π hydrogen bonds were identified in these complexes. CH/π hydrogen bonds occurred, in particular, in the central part of the proline‐rich motifs. Our results suggest that CH/π hydrogen bonds are important in the recognition of SH3 and PRDs by their ligand peptides and play a vital role in the signal transduction system. Combined use of the FMO method and CHPI analysis is a valuable tool for the study of protein/protein and protein/ligand interactions and may be useful in rational drug design. © 2011 Wiley Periodicals, Inc. J Comput Chem 2011  相似文献   

2.
All electron calculations were performed on the photosynthetic reaction center of Blastochloris viridis, using the fragment molecular orbital (FMO) method. The protein complex of 20,581 atoms and 77,754 electrons was divided into 1398 fragments, and the two‐body expansion of FMO/6‐31G* was applied to calculate the ground state. The excited electronic states of the embedded electron transfer system were separately calculated by the configuration interaction singles approach with the multilayer FMO method. Despite the structural symmetry of the system, asymmetric excitation energies were observed, especially on the bacteriopheophytin molecules. The asymmetry was attributed to electrostatic interaction with the surrounding proteins, in which the cytoplasmic side plays a major role. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010  相似文献   

3.
The CCSD(T) level interaction energies of CH/pi complexes at the basis set limit were estimated. The estimated interaction energies of the benzene complexes with CH(4), CH(3)CH(3), CH(2)CH(2), CHCH, CH(3)NH(2), CH(3)OH, CH(3)OCH(3), CH(3)F, CH(3)Cl, CH(3)ClNH(2), CH(3)ClOH, CH(2)Cl(2), CH(2)FCl, CH(2)F(2), CHCl(3), and CH(3)F(3) are -1.45, -1.82, -2.06, -2.83, -1.94, -1.98, -2.06, -2.31, -2.99, -3.57, -3.71, -4.54, -3.88, -3.22, -5.64, and -4.18 kcal/mol, respectively. Dispersion is the major source of attraction, even if substituents are attached to the carbon atom of the C-H bond. The dispersion interaction between benzene and chlorine atoms, which is not the CH/pi interaction, is the cause of the very large interaction energy of the CHCl(3) complex. Activated CH/pi interaction (acetylene and substituted methanes with two or three electron-withdrawing groups) is not very weak. The nature of the activated CH/pi interaction may be similar to the hydrogen bond. On the other hand, the nature of other typical (nonactivated) CH/pi interactions is completely different from that of the hydrogen bond. The typical CH/pi interaction is significantly weaker than the hydrogen bond. Dispersion interaction is mainly responsible for the attraction in the CH/pi interaction, whereas electrostatic interaction is the major source of attraction in the hydrogen bond. The orientation dependence of the interaction energy of the typical CH/pi interaction energy is very small, whereas the hydrogen bond has strong directionality. The weak directionality suggests that the hydrogen atom of the interacting C-H bond is not essential for the attraction and that the typical CH/pi interaction does not play critical roles in determining the molecular orientation in molecular assemblies.  相似文献   

4.
The structures of 2-phenylethanethiol (PET, PhCH(2)CH(2)SH) and its 1:1 water clusters have been studied using resonant two-photon ionization spectroscopy including band contour analysis and UV-UV holeburning, combined with extensive ab initio calculations on ground and excited states. The most populated conformer, labeled Ggpi, has a gauche arrangement about the SCCC and HSCC bonds that permits a stabilizing SH...dpi type of hydrogen bond. The other observed conformer, Ag, is anti with respect to the SCCC bond. In the dominant 1:1 water cluster, a water molecule binds to the Ggpi conformer via an OH...S hydrogen bond and two significant CH...O interactions. There is also evidence for water binding to conformer Ag with a similar arrangement, and for a second Ggpi cluster where water inserts between the SH and the aromatic ring. The additional interactions to the water molecules result in net D(e) binding energies approximately double those resulting from a single thiol-water hydrogen bond. The (1)(pi,pi(*)) excited state lifetimes in the bare molecules are very short because of internal conversion to a dissociative (1)(n,pi(*)) state related to the thiol. In the dominant Gw(1) cluster, the lifetime is significantly increased from <1 to approximately 4 ns. Hydrogen bonding to the thiol, which raises the energy of the dissociative (1)(n,pi(*)) state, accounts for this behavior.  相似文献   

5.
Peptides and proteins may contain post-translationally modified phosphorylated amino acid residues, in particular phosphorylated serine (pSer), threonine (pThr) and tyrosine (pTyr). Following earlier work by Lehmann et al., the [M-H]- anions of peptides containing pSer and pThr functionality show loss of the elements of H3PO4. This process, illustrated for Ser (and using a model system), is CH3CONH-C(CH2OPO3H2)CONHCH(3) --> [CH3CONHC(==CH2)CONHCH3 (-OPO3H2)] (a) --> [CH3CONHC(==CH2)CONHCH3-H]- + H3PO4, a process endothermic by 83 kJ mol(-1) at the MP2/6-31++G(d,p)//HF/6-31++G(d,p) level of theory. In addition, intermediate (a) may decompose to yield CH3CONHC(==CH2)CONHCH3 + H2PO4 - in a process exothermic by 3 kJ mol(-1). The barrier to the transition state for these two processes is 49 kJ mol(-1). Characteristic cleavages of pSer and pThr are more energetically favourable than the negative ion backbone cleavages of peptides described previously. In contrast, loss of HPO3 from [M-H]- is characteristic of pTyr. The cleavage [NH2CH(CH2-C6H4-OPO3H-)CO2H] --> [NH2C(CH2-C6H4-O-)CO2H (HPO3)] (b) --> NH2CH(CH2-C6H4-O-)CO2H + HPO3 is endothermic by 318 kJ mol(-1) at the HF/6-31+G(d)//AM1 level of theory. In addition, intermediate (b) also yields NH2CH(CH2-C6H4-OH)CO2H + PO3 - (reaction endothermic by 137 kJ mol(-1)). The two negative ion cleavages of pTyr have a barrier to the transition state of 198 kJ mol(-1) (at the HF/6-31+G(d)//AM1 level of theory) comparable with those already reported for negative ion backbone cleavages.  相似文献   

6.
在DFT-B3LYP/6-311++G**水平上分别求得(CH3)2S…ClOH卤键复合物和(CH3)2S…HOCl氢键复合物势能面上的稳定构型. 频率分析表明, 与单体HOCl相比, 在两种复合物中, 10Cl—11O和12H—11O键伸缩振动频率发生显著的红移. 经MP2/6-311++G**水平计算的含基组重叠误差(BSSE)校正的气相中相互作用能分别为-11.69和-24.16 kJ·mol-1. 自然键轨道理论(NBO)分析表明, 在(CH3)2S…ClOH卤键复合物中, 引起10Cl—11O键变长的因素包括两种电荷转移: (i) 孤对电子LP(1S)1→σ*(10Cl—11O); (ii) 孤对电子LP(1S)2→σ*(10Cl—11O), 其中孤对电子LP(1S)2→σ*(10Cl—11O)转移占主要作用, 总的结果是使σ*(10Cl—11O)的自然布居数增加0.14035e, 同时11O原子的再杂化使其与10Cl成键时s成分增加, 即具有与电荷转移作用同样的“拉长效应”; 在(CH3)2S…HOCl氢键复合物中也存在类似的电荷转移, 但是11O原子的再杂化不同于前者. 自然键共振理论(NRT)进行键序分析表明, 在卤键复合物和氢键复合物中, 10Cl—11O和12H—11O键的键序都减小. 通过分子中原子理论(AIM)分析了复合物中卤键和氢键的电子密度拓扑性质.  相似文献   

7.
Several methods including molecular mechanics, molecular dynamics, ONIOM that combines quantum chemistry with molecular mechanics and standard quantum chemistry are used to study the configuration and electron structures of an adduct of the DMA segment d(ATACATG*G*TACATA)-d(TATGTACCATGTAT) with cis-[Pt(NH3)(2-Picoline)]2+. The investigation shows that the configuration optimized by ONIOM is similar to that determined by NMR. Strong chemical bonds between Pt of the complex and two N7s of neighboring guanines in the DNA duplex and hydrogen bond between the NH3of the complex and O6 of a nearby guanine have a large impact on the configuration of the adduct. Chemical bonds, the aforementioned hydrogen bond, and the interaction between a methyl of the complex and a methyl of the base in close proximity are critical for the complex to specifically recognize DNA.  相似文献   

8.
A quantum chemical model is introduced to predict the H‐bond donor strength of monofunctional organic compounds from their ground‐state electronic properties. The model covers ? OH, ? NH, and ? CH as H‐bond donor sites and was calibrated with experimental values for the Abraham H‐bond donor strength parameter A using the ab initio and density functional theory levels HF/6‐31G** and B3LYP/6‐31G**. Starting with the Morokuma analysis of hydrogen bonding, the electrostatic (ES), polarizability (PL), and charge transfer (CT) components were quantified employing local molecular parameters. With hydrogen net atomic charges calculated from both natural population analysis and the ES potential scheme, the ES term turned out to provide only marginal contributions to the Abraham parameter A, except for weak hydrogen bonds associated with acidic ? CH sites. Accordingly, A is governed by PL and CT contributions. The PL component was characterized through a new measure of the local molecular hardness at hydrogen, η(H), which in turn was quantified through empirically defined site‐specific effective donor and acceptor energies, EEocc and EEvac. The latter parameter was also used to address the CT contribution to A. With an initial training set of 77 compounds, HF/6‐31G** yielded a squared correlation coefficient, r2, of 0.91. Essentially identical statistics were achieved for a separate test set of 429 compounds and for the recalibrated model when using all 506 compounds. B3LYP/6‐31G** yielded slightly inferior statistics. The discussion includes subset statistics for compounds containing ? OH, ? NH, and active ? CH sites and a nonlinear model extension with slightly improved statistics (r2 = 0.92). © 2008 Wiley Periodicals, Inc. J Comput Chem 2009  相似文献   

9.
The hydrogen-bonding ability of five-membered heteroaromatic molecules containing one chalcogen and two heteroatoms with nitrogen in addition to chalcogen, respectively, have been analyzed using density functional and molecular orbital methods through adduct formation with water. The stabilization energies for all the adducts are established at B3LYP/6-31+G* and MP2/6-31+G* levels after correcting for the basis set superposition error by using the counterpoise method and also corrected for zero-point vibrational energies. A natural bond orbital analysis at B3LYP/6-31+G* level and natural energy decomposition analysis at HF/6-31+G* using MP2/6-31+G* geometries have been carried out to understand the nature of hydrogen-bonding interaction in monohydrated heterocyclic adducts. Nucleus-independent chemical shift have been evaluated to understand the correlation between hydrogen bond formation and aromaticity.  相似文献   

10.
Recently reported high-level ab initio calculations and gas phase spectroscopic measurements show that the nature of CH/pi interactions is considerably different from conventional hydrogen bonds, although the CH/pi interactions were often regarded as the weakest class of hydrogen bonds. The major source of attraction in the CH/pi interaction is the dispersion interaction and the electrostatic contribution is small, while the electrostatic interaction is mainly responsible for the attraction in the conventional hydrogen bonds. The nature of the "typical" CH/pi interactions is similar to that of van der Waals interactions, if some exceptional "activated" CH/pi interactions of highly acidic C-H bonds are excluded. Shifts of C-H vibrational frequencies and electronic spectra also support the similarity. The hydrogen bond is important in controlling structures of molecular assemblies, since the hydrogen bond is sufficiently strong and directional due to the large electrostatic contribution. On the other hand, the directionality of the "typical" CH/pi interaction is very weak. Although the "typical" CH/pi interaction is often regarded as an important interaction in controlling the structures of molecular assemblies as in the cases of conventional hydrogen bonds, the importance of the "typical" CH/pi interactions is questionable.  相似文献   

11.
The gas-phase electron transmission (ET) and dissociative electron attachment (DEA) spectra are reported for the series of (bromoalkyl)benzenes C6H5(CH2)nBr (n = 0-3), where the bromine atom is directly bonded to a benzene ring or separated from it by 1-3 CH2 groups, and the dihalo derivative 1-Br-4-Cl-benzene. The relative DEA cross sections (essentially due to the Br- fragment) are reported, and the absolute cross sections are also evaluated. HF/6-31G and B3LYP/6-31G* calculations are employed to evaluate the virtual orbital energies (VOEs) for the optimized geometries of the neutral state molecules. The pi* VOEs, scaled with empirical equations, satisfactorily reproduce the corresponding experimental vertical electron attachment energies (VAEs). According to the calculated localization properties, the LUMO (as well as the singly occupied MO of the lowest lying anion state) of C6H5(CH2)3Br is largely localized on both the benzene ring and the C-Br bond, despite only a small pi*/sigma*C-Br interaction and in contrast to the chlorine analogue where the LUMO is predicted to possess essentially ring pi character. This would imply a less important role of intramolecular electron transfer in the bromo derivative for production of the halogen negative fragment through dissociation of the first resonant state. The VAEs calculated as the anion/neutral energy difference with the 6-31+G* basis set which includes diffuse functions are relatively close to the experimental values but do not parallel their sequence. In addition the SOMO of some compounds is not described as a valence MO with large pi* character but as a diffuse sigma* MO.  相似文献   

12.
在DFT-B3LYP及MP2/6-311++G**水平上分别求得CH3SH…HOCl氢键复合物和CH3SH…ClOH卤键复合物势能面上的稳定构型. 频率分析表明, 与单体HOCl相比, 在两种复合物中, Cl(9)—O(7)和H(8)—O(7)键伸缩振动频率发生显著的红移. 经MP2/6-311++G**水平计算的含基组重叠误差(BSSE)校正的气相中相互作用能分别为-19.23和-6.85 kJ&#8226;mol-1. 自然键轨道理论(NBO)分析表明, 在CH3SH…ClOH卤键复合物中, 引起Cl(9)—O(7)键变长的因素包括2种电荷转移: (i)孤对电子LP[S(1)]1→σ*[Cl(9)—O(7)]; (ii)孤对电子LP[S(1)]2→σ*[Cl(9)—O(7)], 其中孤对电子LP[S(1)]2→σ*[Cl(9)—O(7)]转移占主要作用, 总的结果是使σ*[Cl(10)—O(11)]的自然布居数增加, 同时O(7)和Cl(9)原子s成分均增加的杂化重优具有与电荷转移作用相同的“拉长效应”; 在CH3SH…HOCl氢键复合物中也存在类似的电荷转移, 但是O(7)原子的再杂化效应不同于前者. 自然键共振理论(NRT)进行键序分析表明, 在氢键复合物和卤键复合物中, H(8)—O(7)和Cl(9)—O(7)键的键序都减小. 通过分子中原子理论(AIM)分析了复合物中氢键和卤键的电子密度拓扑性质.  相似文献   

13.
The molecular structures of the diphosphines P(2)[CH(SiH(3))(2)](4), P(2)[C(SiH(3))(3)](4), P(2)[SiH(CH(3))(2)](4), and P(2)[Si(CH(3))(3)](4) and the corresponding radicals P[CH(SiH(3))(2)](2), P[C(SiH(3))(3)](2), P[SiH(CH(3))(2)](2), and P[Si(CH(3))(3)](2) were predicted by theoretical quantum chemical calculations at the HF/3-21G*, B3LYP/3-21G*, and MP2/6-31+G* levels. The conformational analyses of all structures found the gauche conformers of the diphosphines with C(2) symmetry to be the most stable. The most stable conformers of the phosphido radicals were also found to possess C(2) symmetry. The structural changes upon dissociation allow the release of some of the energy stored in the substituents and therefore contribute to the decrease of the P-P bond dissociation energy. The P-P bond dissociation enthalpies at 298 K in the compounds studied were calculated to vary from -11.4 kJ mol(-1) (P(2)[C(SiH(3))(3)](4)) to 179.0 kJ mol(-1) (P(2)[SiH(CH(3))(2)](4)) at the B3LYP/3-21G* level. The MP2/6-31+G* calculations predict them to be in the range of 52.8-207.9 kJ mol(-1). All the values are corrected for basis set superposition error. The P-P bond energy defined by applying a mechanical analogy of the flexible substituents connected by a spring shows less variation, between 191.3 and 222.6 kJ mol(-1) at the B3LYP/3-21G level and between 225.6 and 290.4 kJ mol(-1) at the MP2/6-31+G* level. Its average value can be used to estimate bond dissociation energies from the energetics of structural relaxation.  相似文献   

14.
在DFT-B3LYP/6-311++G**水平上分别求得(CH3)2S…HOO和(CH3)2O…HOO开壳型氢键复合物势能面上的稳定构型. 频率分析表明, 与单体HOO自由基相比, 复合物中H10-O11键伸缩振动频率发生显著的红移, 红移值分别为424.21和374.22 cm-1. 在MP2/6-311++G**水平计算得到, 含基组重叠误差(BSSE)校正和零点振动能(ZPVE)校正的相互作用能分别为-24.68和-31.01 kJ·mol-1. 自然键轨道(NBO)理论分析表明, 在(CH3)2S…HOO复合物中, 引起H10-O11键变长的因素包括两种电荷转移: (1) LP(S1)1→σ*(H10-O11); (2) LP(S1)2→σ*(H10-O11), 其中LP(S1)2→σ*(H10-O11)占主要作用, 总的结果是使σ*(H10-O11)的自然布居数增加了37.27 me; 在(CH3)2O…HOO中也有相似的电荷转移的超共轭作用. AIM理论分析表明, S1…H10间和O1…H10间都存在键鞍点, ▽2ρ(r)分别为0.06196和0.03745, 说明这种相互作用介于共价键和离子键之间, 偏于静电作用.  相似文献   

15.
The analytic energy gradients for the combined fragment molecular orbital and polarizable continuum model (FMO/PCM) method are derived and implemented. Applications of FMO/PCM geometry optimization to polyalanine show that the structures obtained with the FMO/PCM method are very close to those obtained with the corresponding full ab initio PCM methods. FMO/PCM (RHF/6‐31G* level) is used to optimize the solution structure of the 304‐atom Trp‐cage miniprotein and the result is in agreement with NMR experiments. The key factors determining the relative stability of the α‐helix, β‐turn and the extended form in solution are elucidated for polyalanine. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010  相似文献   

16.
片断的UHF运算不能保证每个片断轨道具有确切的电子占据数,故Kost定域化是必需的.当片断产生于多键断裂时,在确保目标轨道单占据性的同时,Kost定域破坏了轨道基组的对称性.为此,在Kost定域化后,必须对单占据轨道作2×2对称化旋转后,再作有条件的RHF运算.以乙烯基片断CHCH(波二烯分子中的一个片断)为例,详述了对称化的方法、原理和计算程序.以C-H片断为例,细述三单键片断轨道基组对称化的特殊性.介绍C-HR参考键长选择的判据,探讨键长与选择Gaussian基组大小的关系.  相似文献   

17.
Ab initio calculations at the Hartree-Fock, MP2 and MP4 levels were performed to find structures of the equilibrium and transition states and the reaction energies and energies of activation of several competing reaction pathways of O (3P)+CH3SH. A 6-31G* basis set was used in all calculations. The mechanism of hydrogen atom abstraction from the S-H group methanethiol was found to be very competitive with the oxygen atom addition to the sulfur atom.  相似文献   

18.
The CH3S* + O2 reaction system is considered an important process in atmospheric chemistry and in combustion as a pathway for the exothermic conversion of methane-thiyl radical, CH3S*. Several density functional and ab initio computational methods are used in this study to determine thermochemical parameters, reaction paths, and kinetic barriers in the CH3S* + O2 reaction system. The data are also used to evaluate feasibility of the DFT methods for higher molecular weight oxy-sulfur hydrocarbons, where sulfur presents added complexity from its many valence states. The methods include: B3LYP/6-311++G(d,p), B3LYP/6-311++G(3df,2p), CCSD(T)/6-311G(d,p)//MP2/6-31G(d,p), B3P86/6-311G(2d,2p)//B3P86/6-31G(d), B3PW91/6-311++G(3df,2p), G3MP2, and CBS-QB3. The well depth for the CH3S* + 3O2 reaction to the syn-CH3SOO* adduct is found to be 9.7 kcal/mol. Low barrier exit channels from the syn-CH3SOO* adduct include: CH2S + HO2, (TS6, E(a) is 12.5 kcal/mol), CH3 + SO2 via CH3SO2 (TS2', E(a) is 17.8) and CH3SO + O (TS17, E(a) is 24.7) where the activation energy is relative to the syn-CH3SOO* stabilized adduct. The transition state (TS5) for formation of the CH3SOO adduct from CH3S* + O2 and the reverse dissociation of CH3SOO to CH3S* + O2 is relatively tight compared to typical association and simple bond dissociation reactions; this is a result of the very weak interaction. Reverse reaction is the dominant dissociation path due to enthalpy and entropy considerations. The rate constants from the chemical activation reaction and from the stabilized adduct to these products are estimated as functions of temperature and pressure. Our forward rate constant and CH3S loss profile are in agreement with the experiments under similar conditions. Of the methods above, the G3MP2 and CBS-QB3 composite methods are recommended for thermochemical determinations on these carbon-sulfur-oxygen systems, when they are feasible.  相似文献   

19.
吡啶-BH~3相互作用复合物的理论研究   总被引:2,自引:2,他引:2  
对吡啶-BH~3复合物分别用MP2/6-31+G^*和B3LYP/6-31+G^*进行理论计算以预测该复合物的构型及解离能,得到四种构型,在MP2优化构型基础上作CCSD/6-31+G^*单点能量计算以验证MP2与B3LYP结果的可靠性,然后用B3LYP作振动频率分析,计算了各构型的垂直电离势,最后用更大基组作单点能量计算和自然键轨道(NBO)分析。结果表明,N-B直接相连的构型最稳定,其解离能为141.50kJ/mol,MP2和B3LYP对N-H接近的构型结果相关较大,另外两种构型稳定性介于二者之间,解离能分别为15.18kJ/mol,14.06kJ/mol(MP2/6-31+G^*)。  相似文献   

20.
The title compound, C31H37NO4S [systematic name: (R)-tert-butyl-2-[(tert-butoxycarbonyl)amino]-3-(tritylsulfanyl)propanoate] is an L-cysteine derivative with three functions: NH2, COOH and SH, blocked by protecting groups tert-butoxycarbonyl, tert-butyl and trityl, respectively. The main chain of the molecule adopts the extended, nearly all-trans C5 conformation with the intramolecular N-H...O=C hydrogen bond. The urethane group is not involved in any intermolecular hydrogen bonding. Only weak intermolecular hydrogen bonds and hydrophobic contacts are observed in the crystal structure. These are C-H...O hydrogen bonds and CH/pi interactions with donor...acceptor distances, C...O ca. 3.5 A and C...C ca. 3.7 A, respectively. The first type of interaction links phenyl H-atoms and carbonyl groups. The second type of interaction is formed between a methyl group of the tert-butyl fragment and a trityl phenyl ring. The resulting molecular conformation in the crystal is very close to an ab initio minimum energy conformer of the isolated molecule. The extended C5 conformation of the main peptide chain is the same and there is slight discrepancy in the disposition of trityl phenyl rings. Their small dislocation creates the possibility of forming the entire network above of extensive, specific, weak intermolecular interactions; these constrain the molecule and permit it to retain the minimum energy C5 conformation of its main chain in the solid state. In contrast, in n-hexane solution, where such specific interactions cannot occur, only a small population of the molecules adopts the extended C5 conformation.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号