Predicting relative binding affinities of non-peptide HIV protease inhibitors with free energy perturbation calculations

The relative binding free energies in HIV protease of haloperidol thioketal (THK) and three of its derivatives were examined with free energy calculations. THK is a weak inhibitor (IC50 = 15 M) for which two cocrystal structures with HIV type 1 proteases have been solved [Rutenber, E. et al., J. Biol. Chem., 268 (1993) 15343]. A THK derivative with a phenyl group on C2 of the piperidine ring was expected to be a poor inhibitor based on experiments with haloperidol ketal and its 2- phenyl derivative (Caldera, P., personal communication). Our calculations predict that a 5-phenyl THK derivative, suggested based on examination of the crystal structure, will bind significantly better than THK. Although there are large error bars as estimated from hysteresis, the calculations predict that the 5-phenyl substituent is clearly favored over the 2-phenyl derivative as well as the parent compound. The unfavorable free energies of solvation of both phenyl THK derivatives relative to the parent compound contributed to their predicted binding free energies. In a third simulation, the change in binding free energy for 5-benzyl THK relative to THK was calculated. Although this derivative has a lower free energy in the protein, its decreased free energy of solvation increases the predicted G(bind) to the same range as that of the 2-phenyl derivative.     

Surface free energy effect on bacterial retention

Bacterial infection is one of the most frequent and severe complications in the long-term effectiveness of medical implants and devices, greatly increasing treatment cost and inconvenience to the patient. Surface physical and chemical properties are known to influence the extent and form of bacterial infection, although the exact correlation with specific properties is difficult due to the complexity of the system. One approach in the attempt to reduce the bacterial colonisation is to modify the surface energy and chemistry, so as to influence the interactions between the surface and the bacteria that come into contact with it. Five types of coatings were investigated in this study, together with silicone, and polished and non-polished stainless steel 316L. Surfaces were tested for retention of Pseudomonas aeruginosa AK1 after 1 h. A good correlation (>90%) was found between P. aeruginosa AK1 retention and total surface free energy, as well as its polar and dispersive components. The minimum level of P. aeruginosa AK1 retention was found for a range of total surface free energy in the range 20–27 mN/m.     

吲哚并喹啉衍生物与G-四链体相互作用的分子模拟研究

G-四链体是富含鸟嘌呤碱基的DNA序列通过氢键相互作用形成的四链螺旋结构. 通过小分子化合物诱导与稳定端粒G-四链体从而抑制端粒酶活性是一种新的抗癌策略. 为了研究一系列吲哚并喹啉衍生物与端粒G-四链体的相互作用, 探究其相互作用模式, 从而为实现基于G-四链体结构的药物合理设计提供依据, 使用分子对接的方法构建了吲哚并喹啉衍生物与G-四链体复合物结构, 在此基础上进行分子动力学模拟, 并使用线性相互作用能(LIE)方法计算了化合物与G-四链体的结合自由能. 结果表明: 化合物与G-四链体的主要相互作用方式由氢键、静电与π-π堆积作用构成, 侧链末端基团类型和侧链的长短是影响相互作用强弱的重要因素. 通过LIE方法计算的结合自由能与实验结果基本吻合, 相关度达到r²＝0.79. 并且, 基于预测的结合模式, 总结了拥有更高活性的新型吲哚并喹啉衍生物应具有的几个结构特征.     

Theoretical calculations: Can Gibbs free energy for intermolecular complexes be predicted efficiently and accurately?

The theoretical study has been performed to refine the procedure for calculations of Gibbs free energy with a relative accuracy of less than 1 kcal/mol. Three benchmark intermolecular complexes are examined via several quantum-chemical methods, including the second-order Moller-Plesset perturbation (MP2), coupled cluster (CCSD(T)), and density functional (BLYP, B3LYP) theories augmented by Dunnings correlation-consistent basis sets. The effects of electron correlation, basis set size, and anharmonicity are systematically analyzed, and the results are compared with available experimental data. The results of the calculations suggest that experimental accuracy can be reached only by extrapolation of MP2 and CCSD(T) total energies to the complete basis set. The contribution of anharmonicity to the zero point energy and TDeltaSint values is fairly small. The new, economic way to reach chemical accuracy in the calculations of the thermodynamic parameters of intermolecular interactions is proposed. In addition, interaction energy (De) and free energy change (DeltaA) for considered species have been evaluated by Carr-Parrinello molecular dynamics (CPMD) simulations and static BLYP-plane wave calculations. The free energy change along the reaction paths were determined by the thermodynamic integration/\"Blue Moon Ensemble\" technique. Comparison between obtained values, and available experimental and conventional ab initio results has been made. We found that the accuracy of CPMD simulations is affected by several factors, including statistical uncertainty and convergence of constrained forces (TD integration), and the nature of DFT (density functional theory) functional. The results show that CPMD technique is capable of reproducing interaction and free energy with an accuracy of 1 kcal/mol and 2-3 kcal/mol respectively.     

MMP-2和Hydroxamate类抑制剂绝对自由能的计算

采用基于线性响应近似的自由能计算方法计算了一类hydroxamate抑制剂和MMP-2的绝对结合自由能。计算中,催化锌离子和MMP-2以及配体之间采用了非键模型。分子动力学模拟结果显示,采用非键模型时,催化Zn离子采用五配位的形式,但配位键的形式和初始结构比较有很大的差别。通过拟合,分别得到了单参数、双参数以及三参数的自由能预测模型,其中,含有常数校正项的三参数模型具有最佳的预测能力,预测自由能和实际自由能之间平均绝对误差仅为2.38kJ/mol。     

Deciphering the Selectivity of Bombyx mori Pheromone Binding Protein for Bombykol over Bombykal: A Theoretical Approach

In this article we report calculations dedicated to estimate the selectivity of the Bombyx mori pheromone binding protein towards the two closely related pheromonal components Bombykol and Bombykal. The selectivity is quantified by the binding free‐energy difference, obtained either by the thermodynamic integration or by the MM‐GBSA approach. In the latter, the selectivity is decomposed on a per‐residue basis, which identifies the residues considered crucial for the selectivity of the protein for Bombykol over Bombykal. A discussion on the role of Bombyx mori pheromone binding protein is provided on the basis of these results.     

Structure and Thermodynamics of Mg:Phosphate Interactions in Water: A Simulation Study

The association of Mg²⁺ and H₂PO₄^? in water can give insights into Mg:phosphate interactions in general, which are very widespread, but for which experimental data is surprisingly sparse. It is studied through molecular dynamics simulations (>100 ns) by using the polarizable AMOEBA force field, and the association free energy is computed for the first time. Explicit consideration of outer‐sphere and two types of inner‐sphere association provides considerable insight into the dynamics and thermodynamics of ion pairing. After careful assessment of the computational approximations, the agreement with experimental values indicates that the methodology can be extended to other inorganic and biological Mg:phosphate interactions in solution.     

Computation of affinity and selectivity: Binding of 2,4-diaminopteridine and 2,4-diaminoquinazoline inhibitors to dihydrofolate reductases

Binding energy calculations for complexes of mutant and wild-type human dihydrofolate reductases with 2,4-diaminopteridine and 2,4-diaminoquinazoline inhibitors are reported. Quantitative insight into binding energetics of these molecules is obtained from calculations based on force field energy evaluation and thermal sampling by molecular dynamics simulations. The calculated affinity of methotrexate for wild-type and mutant enzymes is reasonably well reproduced. Truncation of the methotrexate glutamate tail results in a loss of affinity by several orders of magnitude. No major difference in binding strength is predicted between the pteridines and the quinazolines, while the N-methyl group present in methotrexate appears to confer significantly stronger binding. The recent improvement, which is used here, of our linear interaction energy method for binding affinity prediction, as well as problems with treating charged and flexible ligands are discussed. This approach should be suitable in a drug discovery context for prediction of binding energies of new inhibitors prior to their synthesis, when some information about the binding mode is available.