Molecular structural formulas as one-electron density and hamiltonian operators: the VIF method extended

The valency interaction formula (VIF) method is given a broader and more general interpretation in which these simple molecular structural formulas implicitly include all overlaps between valence atomic orbitals even for interactions not drawn in the VIF picture. This applies for VIF pictures as one-electron Hamiltonian operators as well as VIF pictures as one-electron density operators that constitute a new implementation of the VIF method simpler in its application and more accurate in its results than previous approaches. A procedure for estimating elements of the effective charge density-bond order matrix, Pmunu, from electron configurations in atoms is presented, and it is shown how these lead to loop and line constants in the VIF picture. From these structural formulas, one finds the number of singly, doubly, and unoccupied molecular orbitals, as well as the number of molecular orbitals with energy lower, equal, and higher than -1/2Eh, the negative of the hydrogen atom's ionization energy. The VIF results for water are in qualitative agreement with MP2/6311++G3df3pd, MO energy levels where the simple VIF for water presented in the earlier literature does not agree with computed energy levels. The method presented here gives the simplest accurate VIF pictures for hydrocarbons. It is shown how VIF can be used to predict thermal barriers to chemical reactions. Insertion of singlet carbene into H2 is given as an example. VIF pictures as one-electron density operators describe the ground-state multiplicities of B2, N2, and O2 molecules and as one-electron Hamiltonian operators give the correct electronegativity trend across period two. Previous implementations of VIF do not indicate singly occupied molecular orbitals directly from the pictorial VIF rules for these examples. The direct comparison between structural formulas that represent electron density and those that represent energy is supported by comparison of a simple electronegativity scale, chiD=N/n2, with well-known electronegativity scales of Pauling, Mulliken, and Allen. This scale comes from the method used to calculate Pmumu for sp3 hybridized period-two elements and is comparable to electronegativity because it has the same form as <1/r> for hydrogenic orbitals. It therefore provides a physical basis for the representation of one electron density and Hamiltonian operators by the same VIF picture.     

用分子对接方法研究HIV-1整合酶与病毒DNA的结合模式

用分子对接方法研究了HIV-1整合酶(Integrase, IN)二聚体与3’ 端加工(3’ Processing, 3’-P)前的8 bp及27 bp病毒DNA的相互作用, 并获得IN与27 bp病毒DNA的特异性结合模式. 模拟结果表明, IN有特异性DNA结合区和非特异性DNA结合区; IN二聚体B链的K14, R20, K156, K159, K160, K186, K188, R199和A链的K219, W243, K244, R262, R263是IN结合病毒DNA的关键残基; 并从结构上解释了能使IN发挥活性的病毒DNA的最小长度是15 bp. 通过分析结合能发现, IN与DNA稳定结合的主要因素是非极性相互作用, 而关键残基与病毒DNA相互识别主要依赖于极性相互作用. 模拟结果与实验数据较吻合.     

Molecular mechanism of HIV-1 integrase-vDNA interactions and strand transfer inhibitor action: a molecular modeling perspective

Human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an essential enzyme for splicing a viral DNA (vDNA) replica of its genome into host cell chromosomal DNA (hDNA) and has been recently recognized as a promising therapeutic target for developing anti-AIDS agents. The interaction between HIV-1 IN and vDNA plays an important role in the integration process of the virus. However, a detailed understanding about the mechanism of this interactions as well as the action of the anti-HIV drug raltegravir (RAL, approved by FDA in 2007) targeting HIV-1 IN in the inhibition of the vDNA strand transfer is still absent. In the present work, a molecular modeling study by combining homology modeling, molecular dynamics (MD) simulations with molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), and molecular mechanics Generalized-Born surface area (MM-GBSA) calculations was performed to investigate the molecular mechanism of HIV-1 IN-vDNA interactions and the inhibition action of vDNA strand transfer inhibitor (INSTI) RAL. The structural analysis showed that RAL did not influence the interaction between vDNA and HIV-1 IN, but rather targeted a special conformation of HIV-1 IN to compete with host DNA and block the function of HIV-1 IN by forcing the 3'-OH of the terminal A17 nucleotide away from the three catalytic residues (Asp64, Asp116, and Glu152) and two Mg(2+) ions. Thus, the obtained results could be helpful for understanding of the integration process of the HIV-1 virus and provide some new clues for the rational design and discovery of potential compounds that would specifically block HIV-1 virus replication.     

用分子对接方法预测HIV-1整合酶与金精三羧酸抑制剂的相互作用

用分子对接方法(Docking)研究了HIV-1整合酶与其抑制剂金精三羧酸的结合过程.为弄清金属离子在结合中所起的作用,选择含有一个Mg+2或不含Mg+2的两种不同的整合酶受体分别与金精三羧酸对接.结果表明, Mg+2对稳定配体与受体的结合起了重要作用. 金精三羧酸配体与含有一个金属Mg+2的整合酶受体对接,最优结合自由能为-45.19 kJ/mol. 当Mg+2失去后,整合酶的活性中心构象将发生变化,使金精三羧酸抑制剂与整合酶的结合自由能(-24.35 kJ/mol)明显增加. 预测了未知的HIV-1整合酶与其抑制剂金精三羧酸的复合物结构, 并可对基于结构的抗HIV-1整合酶的药物设计提供重要信息.     

HIV-1整合酶与抑制剂金精三羧酸复合物的分子动力学模拟

用分子对接程序(Autodock)将含有一个Mg²⁺的HIV-1整合酶核心区(以下简称IN-A)与抑制剂小分子金精三羧酸(简称Aurin)进行对接,预测其未知的复合物结构,然后用分子动力学(MD)方法对IN-A与Aurin的对接结果进行了950 ps的模拟.MD模拟结果发现,IN-A与Aurin形成了两个稳定的氢键,Mg²⁺也与Aurin上的氧原子形成了稳定的配键,IN-A与Aurin之间的静电相互作用能和范德华相互作用能的平均值分别为-205.8和-162.7 kJ/mol.根据MD模拟得到的IN-A与Aurin相互作用后的构象变化信息,我们对对接复合物结构进行了修正,给出了更加合理和稳定的复合物预测结构.本工作得到的HIV-1整合酶与抑制剂Aurin的结合模式信息将有助于设计和改造出效果更好的抗HIV-1整合酶的先导化合物.     

Molecular dynamic and free energy studies of primary resistance mutations in HIV-1 protease-ritonavir complexes

To understand the basis of drug resistance of the HIV-1 protease, molecular dynamic (MD) and free energy calculations of the wild-type and three primary resistance mutants, V82F, I84V, and V82F/I84V, of HIV-1 protease complexed with ritonavir were carried out. Analysis of the MD trajectories revealed overall structures of the protein and the hydrogen bonding of the catalytic residues to ritonavir were similar in all four complexes. Substantial differences were also found near the catalytic binding domain, of which the double mutant complex has the greatest impact on conformational changes of the protein and the inhibitor. The tip of the HIV-1 protease flap of the double mutant has the greater degree of opening with respect to that of the others. Additionally, the phenyl ring of Phe82 moves away from the binding pocket S1', and the conformational change of ritonavir subsite P1' consequently affects the cavity size of the protein and the conformational energy of the inhibitor. Calculations of binding free energy using the solvent continuum model were able to reproduce the same trend of the experimental inhibition constant. The results show that the resistance mutants require hydrophobic residues to maintain the interactions in the binding pocket. Changes of the cavity volume correlate well with free energy penalties due to the mutation and are responsible for the loss of drug susceptibility.     

金属离子对HIV-1整合酶与硫氮硫扎平抑制剂作用模式影响的分子模拟研究

HIV-1整合酶（IN）通过依赖金属离子的两步反应将病毒DNA整合入宿主细胞过程中。结合于HIV-1上的金属离子个数的变化直接影响整合酶与抑制剂之间的结合。本工作用同源模建方法搭建了每条单链核心区具有两个Mg2+ 的（2Mg-IN-Core）和具有一个Mg2+ 的HIV-1 IN二聚体核心区模型（1Mg-IN-Core）。分子对接分别得到它们与硫氮硫扎平类化合物能量较低的复合物结构,把对接结果进行了比较。研究发现：当整合酶中结合的Mg2+个数改变时,它与抑制剂的结合模式也会发生很大的变化;抑制剂能够特异的且稳定的与2Mg-IN-Core模型的活性位点结合;同时与ASP64和GLU152螯合的那个Mg2+离子对于硫氮硫扎平抑制剂与整合酶上的结合有很大的影响。2Mg-IN-Core模型与抑制剂的复合物平均结构进行了2000 ps的 分子动力学模拟,分析发现同时与ASP64及ASP116螯合的Mg2+与IN蛋白形成了四个稳定的螯合键;同时与ASP64及GLU152螯合的Mg2+可与IN结合、也可与抑制剂形成稳定的配位键,这个Mg2+对IN与硫氮硫扎平抑制剂之间的结合有较大影响。     

Structure-based QSAR analysis of a set of 4-hydroxy-5,6-dihydropyrones as inhibitors of HIV-1 protease: an application of the receptor-dependent (RD) 4D-QSAR formalism

Receptor-dependent (RD) 4D-QSAR models were constructed for a set of 39 4-hydroxy-5,6-dihydropyrone analogue HIV-1 protease inhibitors. The receptor model used in this QSAR analysis was derived from the HIV-1 protease (PDB ID ) crystal structure. The bound ligand in the active site of the enzyme, also a 4-hydroxy-5,6-dihydropyrone analogue, was used as the reference ligand for docking the data set compounds. The optimized RD 4D-QSAR models are not only statistically significant (r(2) = 0.86, q(2) = 0.80 for four- and greater-term models) but also possess reasonable predictivity based on test set predictions. The proposed "active" conformations of the docked analogues in the active site of the enzyme are consistent in overall molecular shape with those suggested from crystallographic studies. Moreover, the RD 4D-QSAR models also "capture" the existence of specific induced-fit interactions between the enzyme active site and each specific inhibitor. Hydrophobic interactions, steric shape requirements, and hydrogen bonding of the 4-hydroxy-5,6-dihydropyrone analogues with the HIV-1 protease binding site model dominate the RD 4D-QSAR models in a manner again consistent with experimental conclusions. Some possible hypotheses for the development of new lead HIV-1 protease inhibitors can be inferred from the RD 4D-QSAR models.     

A quantum mechanics/molecular mechanics study of the protein-ligand interaction for inhibitors of HIV-1 integrase

Human immunodeficiency virus type-1 integrase (HIV-1 IN) is an essential enzyme for effective viral replication. Diketo acids such as L-731,988 and S-1360 are potent and selective inhibitors of HIV-1 IN. In this study, we used molecular dynamics simulations, within the hybrid quantum mechanics/molecular mechanics (QM/MM) approach, to determine the protein-ligand interaction energy between HIV-1 IN and L-731,988 and 10 of its derivatives and analogues. This hybrid methodology has the advantage that it includes quantum effects such as ligand polarisation upon binding, which can be very important when highly polarisable groups are embedded in anisotropic environments, as for example in metal-containing active sites. Furthermore, an energy decomposition analysis was performed to determine the contributions of individual residues to the enzyme-inhibitor interactions on averaged structures obtained from rather extensive conformational sampling. Analysis of the results reveals first that there is a correlation between protein-ligand interaction energy and experimental strand transfer into human chromosomes and secondly that the Asn-155, Lys-156 and Lys-159 residues and the Mg(2+) ion are crucial to anti-HIV IN activity. These results may explain the available experimental data.     

表皮生长因子受体和抑制剂之间分子对接的研究

研究了表皮生长因子受体(EGFR)和4-苯胺喹唑啉类抑制剂之间的相互作用模式，表皮生长因子受体的三维结构通过同源蛋白模建的方法得到，而抑制剂和靶酶结合复合物结构则通过分子力学和分子动力学结合的方法计算得到。从模拟结果得到的抑制剂和靶酶之间的相互作用模式表明范德华相互作用、疏水相互作用以及氢键相互作用对抑制剂的活性都有重要的影响，抑制剂的苯胺部分位于活性口袋的底部，能够与受体残基的非极性侧链产生很强的范德华和疏水相互作用，抑制剂双环上的取代基团也能和活性口袋外部的部分残基形成一定的范德华和疏水性相互作用，而抑制剂喹唑啉环上的氮原子能和周围的残基形成较强的氢键相互作用，对抑制剂的活性有较大的影响，计算得到抑制剂和靶酶之间的非键相互作用能以及抑制剂和靶酶之间的相互作用信息能够很好地解释抑制剂活性和结构的关系，为全新抑制剂的设计提供了重要的结构信息。     

Molecular Dynamics Simulation of Temperature-dependent Flexibility of Thermophilic Xylose Isomerase

The complex model of Thermus thermophilus xylose isomerase （TtXI） with D-xylose was constructed, and molecular dynamics （MD） simulations were carried out at 300 and 360 K for 10 ns by NAMD2.5. The radius of gyration （Rg）, subunit interactions, and residue flexibility were analyzed. The results show that residues 60-69, 142-148, 169-172, and 332-340 have high flexibility at 300 and 360 K. Residues with higher flexibility at 360 K than that at 300 K can mainly be divided into two groups： one locates in the helix-loophelix region consisting of residues 55-80 in catalytic domain; the other at subunit interfaces. The Rg of catalytic domain at 360 K shows 0.16 A higher than that at 300 K, but Rg of small C-terminal domain has no obvious difference. The results indicate that enhanced Rg of catalytic domain may lead to the intense motion of the active site of TtXI and promote the D-xylose isomization reaction. Eight hydrogen bonds and five ion pairs are reduced at subunit interfaces at 360 K compared with 300 K, that may be the main reason for the decrease in rigidity and increase in activity at high temperature of TtXI. This result also help to explain the cold-adaption phenomenon of TtXI E372G mutant reported previously. Our results reveal the relationship between temperature and structure flexibility of TtXI, and play an important role in understanding the thermostability of thermophile protein with multiple subunits.     

Predicted solution structure of zymogen human coagulation FVII

A model solution structure for the complete tissue factor-free calcium ion-bound human zymogen FVII (residues 1-406) (FVII) has been constructed to study possible conformational changes associated with the activation process and tissue factor (TF) binding. The initial structure for the present model was constructed using the X-ray crystallographic structure of human coagulation FVIIa/TF complex bound with calcium ions (Banner et al., Nature 1996, 380, 41-46). This model was subsequently subjected to lengthy molecular dynamics simulations. The Amber force field in conjunction with the PME electrostatic summation method was employed. The estimated TF free solution structure was then compared with the currently available X-ray crystal structures of FVIIa (with or without TF, variable inhibitor bound) to estimate the restructuring of FVII due to TF binding and activation. The solution structure of the zymogen FVII in the absence of TF is predicted to be an extended domain structure similar to that of the TF-bound X-ray crystal structure. An additional extension of the serine protease (SP) domain of the zymogen above a reference lipid surface by approximately 7 A was in agreement with experiment. Significant Gla-EGF1 and EGF1-EGF2 interdomain motions in the zymogen were observed. Carbohydrate dimers attached to Ser-52 and Ser-60 did not cause restructuring in this domain. Minimal restructuring of the SP domain is found upon inference of the zymogen from the activated form. The catalytic triad residues maintain the H-bonded network while Lys-341 occupies the S1 specific site in the zymogen.     

The Catalytic Mechanism of HIV-1 Integrase for DNA 3'-End Processing Established by QM/MM Calculations

The development of HIV-1 integrase (INT) inhibitors has been hampered by incomplete structural and mechanistic information. Despite the efforts made to overcome these limitations, only one compound has been approved for clinical use so far. In this work, we have used all experimental information available for INT and similar enzymes, to build a model of the holo-integrase:DNA complex that includes an entire central core domain, a ssDNA GCAGT substrate, and two magnesium ions. Subsequently, we used a large array of computational techniques, which included molecular dynamics, thermodynamic integration, and high-level quantum mechanics/molecular mechanics (QM/MM) calculations to study the possible pathways for the mechanism of 3' end processing catalyzed by INT. We found that the only viable mechanism to hydrolyze the DNA substrate is a nucleophilic attack of an active site water molecule to the phosphorus atom of the scissile phosphoester bond, with the attacking water being simultaneously deprotonated by an Mg(2+)-bound hydroxide ion. The unstable leaving oxoanion is protonated by an Mg(2+)-bound water molecule within the same elementary reaction step. This reaction has an activation free energy of 15.4 kcal/mol, well within the limits imposed by the experimental turnover. This work significantly improves the fundamental knowledge on the integrase chemistry. It can also contribute to the discovery of leads against HIV-1 infection as it provides, for the first time, accurate transition states structures that can be successfully used as templates for high-throughput screening of new INT inhibitors.     

HIV-1整合酶与芳香二酮酸类抑制剂相互作用的分子模拟研究

用分子对接方法研究了一系列芳香二酮酸类抑制剂与HIV-1整合酶的识别及相互作用. 结果表明, 抑制剂结合到整合酶Asp64～Leu68, Thr115～Phe121, Gln148～Lys159和Mg2＋所构成的口袋区, 抑制机理与5CITEP相似. 采用分子动力学模拟和MM/PBSA方法计算了芳香二酮酸类抑制剂与整合酶之间的结合自由能, 计算结果与实验值相吻合, 平均绝对偏差为3.6 kJ/mol, 体系范德华相互作用和溶剂化效应的非极性项是利于形成复合物的主要因素. 相关性分析结果表明, 结合自由能值与疏水相互作用有较强的线性相关(R＝0.61), 基于此, 用多元线性回归方法给出了一个能较强预测芳香二酮酸类抑制剂与HIV-1整合酶的结合自由能预测模型, 为后续基于抑制剂结构的抗HIV-1药物分子设计提供指导.     

Insights into the Functions of M-T Hook Structure in HIV Fusion Inhibitor Using Molecular Modeling

HIV-1 membrane fusion plays an important role in the process that HIV-1 entries host cells. As a treatment strategy targeting HIV-1 entry process, fusion inhibitors have been proposed. Nevertheless, development of a short peptide possessing high anti-HIV potency is considered a daunting challenge. He et al. found that two residues, Met626 and Thr627, located the upstream of the C-terminal heptad repeat of the gp41, formed a unique hook-like structure (M-T hook) that can dramatically improve the binding stability and anti-HIV activity of the inhibitors. In this work, we explored the molecular mechanism why M-T hook structure could improve the anti-HIV activity of inhibitors. Firstly, molecular dynamic simulation was used to obtain information on the time evolution between gp41 and ligands. Secondly, based on the simulations, molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) and molecular mechanics Generalized Born surface area (MM-GBSA) methods were used to calculate the binding free energies. The binding free energy of the ligand with M-T hook was considerably higher than the other without M-T. Further studies showed that the hydrophobic interactions made the dominant contribution to the binding free energy. The numbers of Hydrogen bonds between gp41 and the ligand with M-T hook structure were more than the other. These findings should provide insights into the inhibition mechanism of the short peptide fusion inhibitors and be useful for the rational design of novel fusion inhibitors in the future.