Minimum MD simulation length required to achieve reliable results in free energy perturbation calculations: Case study of relative binding free energies of fructose‐1,6‐bisphosphatase inhibitors

In an attempt to establish the criteria for the length of simulation to achieve the desired convergence of free energy calculations, two studies were carried out on chosen complexes of FBPase‐AMP mimics. Calculations were performed for varied length of simulations and for different starting configurations using both conventional‐ and QM/MM‐FEP methods. The results demonstrate that for small perturbations, 1248 ps simulation time could be regarded a reasonable yardstick to achieve convergence of the results. As the simulation time is extended, the errors associated with free energy calculations also gradually tapers off. Moreover, when starting the simulation from different initial configurations of the systems, the results are not changed significantly, when performed for 1248 ps. This study carried on FBPase‐AMP mimics corroborates well with our previous successful demonstration of requirement of simulation time for solvation studies, both by conventional and ab initio FEP. The establishment of aforementioned criteria of simulation length serves a useful benchmark in drug design efforts using FEP methodologies, to draw a meaningful and unequivocal conclusion. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Calculation of relative binding free energy differences for fructose 1,6-bisphosphatase inhibitors using the thermodynamic cycle perturbation approach.

An iterative, computer-assisted, drug design strategy that combines molecular design, molecular mechanics, molecular dynamics (MD), and free energy perturbation (FEP) calculations with compound synthesis, biochemical testing of inhibitors, and crystallographic structure determination of protein-inhibitor complexes was successfully used to predict the rank order of a series of nucleoside monophosphate analogues as fructose 1,6-bisphosphatase (FBPase) inhibitors. The X-ray structure of FBPase complexed with 5-aminoimidazole-4-carboxamide-1-beta-D-ribofuranosyl 5'-monophosphate (ZMP) provided structural information used for subsequent analogue design and free energy calculations. The FEP protocol was validated by calculating the free energy differences for the mutation of ZMP (1) to AMP (2). The calculated results showed a net gain of 1.7 kcal/mol, which agreed with the experimental result of 1.3 kcal/mol. FEP calculations were performed for 18 other AMP analogues. Inhibition constants were determined for over half of these analogues, usually after completion of the calculation, and were consistent with the predictions. Solvation free energy differences between AMP and various AMP analogues proved to be an important factor in binding free energies, suggesting that increased desolvation costs associated with the addition of polar groups to an inhibitor must be overcome by stronger ligand-protein interactions if the structural modification is to enhance inhibitor potency. The results indicate that FEP calculations predict relative binding affinities with high accuracy and provide valuable insight into the factors that influence inhibitor binding and therefore should greatly aid efforts to optimize initial lead compounds and reduce the time required for the discovery of new drug candidates.     

Structure-guided design of AMP mimics that inhibit fructose-1,6-bisphosphatase with high affinity and specificity

AMP binding sites are commonly used by nature for allosteric regulation of enzymes controlling the production and metabolism of carbohydrates and lipids. Since many of these enzymes represent potential drug targets for metabolic diseases, efforts were initiated to discover AMP mimics that bind to AMP-binding sites with high affinity and high enzyme specificity. Herein we report the structure-guided design of potent fructose 1,6-bisphosphatase (FBPase) inhibitors that interact with the AMP binding site on FBPase despite their structural dissimilarity to AMP. Molecular modeling, free-energy perturbation calculations, X-ray crystallography, and enzyme kinetic data guided our redesign of AMP, which began by replacing the 5'-phosphate with a phosphonic acid attached to C8 of the adenine base via a 3-atom spacer. Additional binding affinity was gained by replacing the ribose with an alkyl group that formed van der Waals interactions with a hydrophobic region within the AMP binding site and by replacing the purine nitrogens N1 and N3 with carbons to minimize desolvation energy expenditures. The resulting benzimidazole phosphonic acid, 16, inhibited human FBPase (IC50 = 90 nM) 11-fold more potently than AMP and exhibited high specificity for the AMP binding site on FBPase. 16 also inhibited FBPase in primary rat hepatocytes and correspondingly resulted in concentration-dependent inhibition of the gluconeogenesis pathway. Accordingly, these results suggest that the AMP site of FBPase may represent a potential drug target for reducing the excessive glucose produced by the gluconeogenesis pathway in patients with type 2 diabetes.     

白色念珠菌N-肉豆蔻酰基转移酶中药物结合位点的MCSS分析

采用多拷贝同时搜寻方法(MCSS)分析得到了CaNMT活性位点的疏水区域、氢键结合位点和负电性区域. MCSS计算结果显示, CaNMT活性位点有两个疏水性比较强的区域: 一个由Tyr107, Tyr109, Val108, Phe117, Phe123, Ala127, Phe176和Leu337等残基组成; 另一个由Phe115, Phe240和Phe339组成. CaNMT活性位点发现有两个氢键作用区域, 其中Tyr119, His227, Asn392和Leu451是与已有抑制剂的氢键结合位点, Tyr107, Asn175, Thr211和Asp412是新发现的氢键结合位点, 而且在NMT家族中高度稳定, 它们对设计新结构类型的CaNMT抑制剂具有重要作用. Leu451是负电性兼氢键作用位点, 是抑制剂设计时所必需考虑的位点.     

Improved convergence of binding affinities with free energy perturbation: Application to nonpeptide ligands with pp60src SH2 domain

Free Energy Perturbations (FEP) in the context of Monte Carlo (MC) simulations were conducted to predict the relative free energies of binding for a series of human Src SH2 domain ligands. Two procedures for disappearing atoms during a single-topology FEP are investigated and dramatic differences in free energy convergence behavior are seen. Comparison of these two protocols suggests that the coupling of the removal of angular constraints with the disappearance of an atom may significantly slow free energy convergence. The series of ligands under investigation here cover a range of modifications at the 3-position of 4-({[4-(cyclohexyl methoxy)benzyl]amino}carbonyl) phenyl phosphate. Unlike any other compound in this study, the 3-amide analog can form two hydrogen bonds within the region of the perturbation, one to a backbone amide hydrogen and one to a highly coordinated water molecule. Agreement with experimental trends in binding affinity is seen, although the computed relative free energy of binding of the amido compound is underestimated. These results are reconciled by examination of the hydration energies of model systems, which predict primary amides as too hydrophilic.     

Computational insights into the selectivity mechanism of APP-IP over matrix metalloproteinases

In this work, selectivity mechanism of APP-IP inhibitor (β-amyloid precursor protein-derived inhibitory peptide) over matrix metalloproteinases (MMPs including MMP-2, MMP-7, MMP-9 and MMP-14) was investigated by molecular modeling methods. Among MMPs, MMP-2 is the most favorable one for APP-IP interacting based on our calculations. The predicted binding affinities can give a good explanation of the activity difference of inhibitor APP-IP. In Comparison with MMP-2/APP-IP complex, the side chain of Tyr214^MMP-7 makes the binding pocket so shallow that the whole side chain of Tyr3^APP-IP can not be fully embraced, thus unfavorable for the N-terminal of APP-IP binding to MMP-7. The poor selectivity of APP-IP toward MMP-9 is mainly related with the decrease of interaction between the APP-IP C-terminal and MMP-9 due to the bulky side chains of Pro193 and Gln199, which is in agreement with experiment. The mutations at residues P193A and Q199G of MMP-9 alternate the binding pattern of the C-terminal of APP-IP by forming two new hydrogen bonds and hydrophobic interactions with MMP-9. The mutants favor the binding affinity of MMP-9 largely. For MMP-14/APP-IP, the large steric effect of Phe204^MMP-14 and the weak contributions of the polar residues Asn231^MMP-14 and Thr190^MMP-14 could explain why MMP-14 is non-selective for APP-IP interacting. Here, the molecular modeling methods were successfully employed to explore the selective inhibitor of MMPs, and our work gives valuable information for future rational design of selective peptide inhibitors toward individual MMP.     

Effective synthesis of cyclic peptide yunnanin C and analogues via Ser/Thr ligation (STL)-mediated peptide cyclization

Cyclic peptide yunnanin C isolated from the root of Stellaria yunnanensis was efficiently synthesized in which the linear peptide was prepared by Boc-SPPS and the cyclization was realized by serine/threonine ligation (STL)-mediated cyclization. In addition, nine yunnanin C analogues, including mutations of Tyr7Gly, Tyr7Val, Tyr7Pro, Tyr7Phe, Ser1Thr, Pro2Val, Gly5Pro, Phe6Ala and Ile4Ala, were prepared in the same fashion. Here, we demonstrated that STL-mediated peptide cyclization could be an effective approach to construct cyclic peptides. Except that proline at the C-terminus could retard the cyclization process, cyclization of yunnanin C analogues with various C-terminal amino acids proceeded with fast cyclization rate (<4 h) and only trace amount of dimers (<5%) at a working concentration of 5 mM.     

嗜热蛋白酶PH1704别构中心量子化学计算与突变体动力学简

通过量子化学计算,确定嗜热菌Pyrococcus horikoshii OT3的PH1704蛋白酶别构位点的关键残基为Arg113,Tyr120和Asn129. 其中,Arg113及Asn129与别构抑制剂结合,参与别构调控. Tyr120残基位于亚基交界面附近,并与亲核残基Cys100之间以氢键相连,可通过影响亚基聚合来影响酶的亲核催化. DJ-1超家族的4种构建蛋白的结构显示,120位点位于亚基交界面处,影响亚基的聚合,进而影响蛋白酶的活力,并间接参与别构调控. 分子生物学实验显示,突变体R113T/Y120P/N129D的kcat/km(L·μmol-1·min-1)值是野生型kcat/km值的6倍,h系数由野生型的0.86转变为1.3,负协同效应消失. 113和129位点处阴离子别构剂脱离,从而破坏113,120和129位点间的封闭环结构,使AC交界面α7螺旋(124~129,524~529)间聚合度增强;120位点残基由Tyr转变为Pro,与Cys100间氢键断裂,亲核进攻的阻力减小,从而使酶活力提高,别构负调控消失.     

Use of quantum mechanics/molecular mechanics-based FEP method for calculating relative binding affinities of FBPase inhibitors for type-2 diabetes

A quantum mechanics (QM)/molecular mechanics (MM)-based free energy perturbation (FEP) method, developed recently, provides most accurate estimation of binding affinities. The validity of the method was evaluated for a large set of diverse inhibitors for fructose 1,6-bisphosphatase (FBPase), a target enzyme for type-II diabetes mellitus. The validation set comprises of 22 important structurally different mutations. The calculated relative binding free energies using the QM/MM-based FEP method reproduce the experimental values with exceptional precision of less than ±0.5 kcal/mol. The CPU requirements for QM/MM-based FEP are about fivefold greater than conventional FEP methods, but it is superior in accuracy of predictions. In addition, the QM/MM-based FEP method eliminates the need for time-consuming development of MM force field parameters, which are frequently required for novel inhibitors described by MM. Future automation of the method and parallelization of the code for 128/256/512 cluster computers is expected to enhance the speed and increase its use for drug design and lead optimization. The present application of QM/MM-based FEP method for structurally diverse set of analogs serves to enhance the scope of FEP method and demonstrate the utility of QM/MM-based FEP method for its potential in drug discovery.     

A cytochrome c variant resistant to heme degradation by hydrogen peroxide

BACKGROUND: Cytochrome c has peroxidase-like activity and can catalyze the oxidation of a variety of organic substrates, including aromatic, organosulfur and lipid compounds. Like peroxidases, cytochrome c is inactivated by hydrogen peroxide. During this inactivation the heme prosthetic group is destroyed. RESULTS: Variants of the iso-1-cytochrome c were constructed by site-directed mutagenesis and were found to be more stable in the presence of hydrogen peroxide than the wild type. No heme destruction was detected in a triple variant (Tyr67-->Phe/Asn52-->Ile/Cys102-->Thr) with the catalytic hydrogen peroxide concentration of 1 mM, even following the loss of catalytic activity, whereas both double variants Tyr67-->Phe/Cys102-->Thr and Asn52-->Ile/Cys102-->Thr showed a greater rate of peroxide-induced heme destruction than observed with the wild-type protein. CONCLUSIONS: Heme destruction and catalytic inactivation are two independent processes. An internal water molecule (Wat166) is shown to be important in the heme destruction process. The absence of a protein radical in the resistant variant suggests that the protein radical is necessary in the heme destruction process, but presumably is not involved in the reactions leading up to the protein inactivation.     

嗜热蛋白酶PH1704别构中心量子化学计算与突变体动力学

通过量子化学计算, 确定嗜热菌Pyrococcus horikoshii OT3的PH1704蛋白酶别构位点的关键残基为Arg113, Tyr120和Asn129. 其中, Arg113及Asn129与别构抑制剂结合, 参与别构调控. Tyr120残基位于亚基交界面附近, 并与亲核残基Cys100之间以氢键相连, 可通过影响亚基聚合来影响酶的亲核催化. DJ-1超家族的4种构建蛋白的结构显示, 120位点位于亚基交界面处, 影响亚基的聚合, 进而影响蛋白酶的活力, 并间接参与别构调控. 分子生物学实验显示, 突变体R113T/Y120P/N129D的k_cat/k_m(L·μmol^-1·min^-1)值是野生型k_cat/k_m值的6倍, h系数由野生型的0.86转变为1.3, 负协同效应消失. 113和129位点处阴离子别构剂脱离, 从而破坏113, 120和129位点间的封闭环结构, 使AC交界面α7螺旋(124~129, 524~529)间聚合度增强; 120位点残基由Tyr转变为Pro, 与Cys100间氢键断裂, 亲核进攻的阻力减小, 从而使酶活力提高, 别构负调控消失.     

Carbonic anhydrase inhibitors: X-ray and molecular modeling study for the interaction of a fluorescent antitumor sulfonamide with isozyme II and IX

The X-ray crystal structure of the fluorescent antitumor sulfonamide carbonic anhydrase (CA, EC, 4.2.1.1) inhibitor (4-sulfamoylphenylethyl)thioureido fluorescein (1) in complex with the cytosolic isoform hCA II is reported, together with a modeling study of the adduct of 1 with the tumor-associated isoform hCA IX. Its binding to hCA II is similar to that of other benzesulfonamides, with the ionized sulfonamide coordinated to the Zn2+ ion within the enzyme active site, and also participating in a network of hydrogen bonds with residues Thr199 and Glu106. The scaffold of 1 did not establish polar interactions within the enzyme active site but made hydrophobic contacts (<4.5 A) with Gln92, Val121, Phe131, Val135, Leu198, Thr199, Thr200, and Pro202. The substituted 3-carboxy-amino-phenyl functionality was at van der Waals distance from Phe131, Gly132, and Val135. The bulky tricyclic fluorescein moiety was located at the rim of the active site, on the protein surface, and strongly interacted with the alpha-helix formed by residues Asp130-Val135. All these interactions were preserved in the hCA IX-1 adduct, but the carbonyl moiety of the fluorescein tail of 1 participates in a strong hydrogen bond with the guanidine moiety of Arg130, an amino acid characteristic of the hCA IX active site. This may account for the roughly 2 times higher affinity of 1 for hCA IX over hCA II and may explain why in vivo the compound specifically accumulates only in hypoxic tumors overexpressing CA IX and not in the normal tissues. The compound is in clinical studies as an imaging tool for acute hypoxic tumors.     

谷氨酰胺结合蛋白的分子动力学模拟和自由能计算

谷氨酰胺结合蛋白(Glutamine-binding protein, GlnBp)是大肠杆菌透性酶系统中一个细胞外液底物专一性结合蛋白, 对于细胞外液中谷氨酰胺(Gln)的运输和传递至关重要. 本文运用分子动力学(Molecular dynamics, MD)模拟采样, 考察了GlnBp关键残基与底物Gln之间的相互作用和GlnBp两条铰链的功能差别; 并采用MM-PBSA方法计算了GlnBp与底物Gln的结合自由能. 结果表明: Ph13, Phe50, Thr118和Ile69与底物Gln的范德华相互作用和Arg75, Thr70, Asp157, Gly68, Lys115, Ala67, His156与底物Gln的静电相互作用是结合Gln的主要推动力; 复合物的铰链区85～89柔性大, 对构象开合提供了结构基础; 而铰链区181～185柔性小, 其作用更多是在功能上把底物Gln限制在口袋中; 自由能预测值与实验值吻合. 本研究很好地解释了GlnBp结构与功能的关系, 为进一步了解GlnBp的开合及转运Gln的机制提供了重要的结构信息.     

Roles of the respective loops at complementarity determining region on the antigen-antibody recognition

For the rational design of antibody, it is important to clarify the characteristics of the interaction between antigen and antibody. In this study, we evaluated a contribution of the respective complementarity determining region (CDR) loops on the antibody recognition of antigen by performing molecular dynamics simulations for 20 kinds of antigen-antibody complexes. Ser and Tyr showed high appearance rates at CDR loops and the sum of averaged appearance rates of Ser and Tyr was about 20⿿30% at all the loops. For example, Ser and Tyr occupied 23.9% at the light chain first loop (L1) and 23.6% at the heavy chain third loop (H3). The direct hydrogen bonds between antigen and antibody were not equally distributed over heavy and light chains. That is, about 70% of the hydrogen bonds were observed at CDRs of the heavy chain and also the direct hydrogen bond with the shortest distance mainly existed at the loops of the heavy chain for all the complexes. It was revealed from the comparison in contribution to the binding free energy among CDR loops that the heavy chain (especially at H2 and H3) had significant influence on the binding between antigen and antibody because three CDR loops of the heavy chain showed the lowest binding free energy (οG_bind) in 19 complexes out of 20. Tyr in heavy chain (especially in H2 and H3) largely contributed to οG_bind whereas Ser hardly contributed to οG_bind even if the number of the direct hydrogen bond with Ser was the fourth largest and also the appearance rate at CDR was the highest among 20 kinds of amino acid residues. The contributions ofTrp and Phe, which bear aromatic ring in the side chain, were often observed in the heavy chain although the energetic contribution of these residues was not so high as Tyr. The present computational analysis suggests that Tyr plays an outstanding role for the antigen-antibody interaction and the CDR loops of the heavy chain is critically important for antibody recognition of antigen.     

Stabilizing and destabilizing effects of phenylalanine --> F5-phenylalanine mutations on the folding of a small protein

We report a systematic evaluation of phenylalanine-to-pentafluorophenylalanine (Phe --> F5-Phe) mutants for the 35-residue chicken villin headpiece subdomain (c-VHP), the hydrophobic core of which features a cluster of three Phe side chains (residues 6, 10, and 17). Phe --> F5-Phe mutations are interesting because aryl-perfluoroaryl interactions of optimal geometry are intrinsically more favorable than aryl-aryl interactions and because perfluoroaryl units are more hydrophobic than are analogous aryl units. One mutant, Phe-10 --> F5-Phe, provides enhanced tertiary structural stability relative to the native sequence. The other six mutants analyzed caused a decrease in stability.     

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药物分子设计提供指导.     

牛视紫红质蛋白质中视黄醛的活性位点

视紫红质蛋白是一个跨膜蛋白, 视黄醛(RET)在该蛋白中的活性结合位点涉及到视觉过程机理, 与一些眼科疾病病理有关. 基于牛视紫红质蛋白1U19的蛋白质晶体结构数据, 采用密度泛函理论的B3LYP方法计算RET-Lys296残基与视黄醛分子周围半径为0.6 nm的空间范围30个氨基酸残基相互作用和结合能. 数值显示1U19蛋白中的残基Glu113、Glu181和Glu122是质子化的RET-Lys296残基的活性结合位点, 结合能分别为-333.38、-205.67和-194.56 kJ·mol-1. 这些氨基酸残基带有一个负电荷, 与质子化的RET-Lys296残基发生强烈的离子静电相互作用. 另外几个残基Ala292、Cys187、Phe293、Pro291以及Trp265等与质子化RET-Lys296残基也有相互吸引作用. 当RET-Lys296残基非质子化, 上述相互作用消失, 促使视黄醛分子与视蛋白分离. 研究发现残基Glu113和Glu181周围各自有一个结晶水分子通过双氢键形式起着稳定作用.     

Pseudoreceptor model for ryanodine derivatives at calcium release channels

This paper describes the generation of a pseudoreceptor model for ryanodine receptor (RyR) modulating ryanoids in rabbit skeletal muscle. For this purpose, the molecular modelling software PrGen was applied to correlate experimentally determined and calculated free energies of binding for a set of 15 ryanodine derivatives. The final model indicates a narrow cleft with hydrogen bond donor and acceptor capacities (represented by an Asn) as most crucial for binding the pyrrole carboxylate substituent at C3 of ryanodine. In addition, hydrophobic residues flank the aromatic pyrrole ring (Tyr, Phe, and Ile). Two of those residues (Tyr and Ile) interact with the 2-isopropyl moiety, which seems to contribute to binding. Opposite to the pyrrole locus, a second hydrophobic region (represented by a Leu) restricts ryanodine derivatives in their longitudinal axis and leads to the discrimination of equatorial and axial positioned methyl groups and of polar substituents at C9. Finally, a charged glutamate residue generates strong hydrogen bonding and electrostatic interactions with the hydroxyl groups at C10 and C15. For this binding-site model – composed of six amino acid residues – a correlation for the training set ligands of R = 0.99 (Q² = 0.975) and a root mean square (rms) deviation of 0.568 kcal/mol for the prediction of the binding energies of four test set ligands was obtained. Based on this pseudoreceptor model the putative topology of the real binding site of ryanoids will be discussed.     

疏水蛋白在云母表面的吸附

采用分子动力学模拟方法研究了疏水蛋白(HFBI)在亲水云母表面的吸附过程.通过6组平行的分子动力学模拟得到了2种不同的稳定吸附结构,即通过N端和通过亲水的α螺旋与表面吸附,得到了5种吸附残基.进一步用自适应偏置力方法计算了所有吸附残基与表面的结合自由能.结果表明,残基Lys是吸附过程的关键残基,即当HFBI通过含有Lys残基的α螺旋与云母表面作用时,其吸附构象最稳定.静电作用是吸附过程的主要驱动力.在该吸附结构中,HFBI的疏水面暴露在溶液中,有效降低了云母表面的润湿性.     

Amino acid-amino acid interactions studied by charge-transfer chromatography

The interaction between amino acids was studied by charge-transfer reversed-phase thin-layer chromatography. The dependence of the lipophilicity of Trp on the concentration of other amino acids in the eluent was considered to be linearly related to the strength of interaction. Arg, Asn, Glu, Met, Phe and Thr interacted with Trp; Ala, Gly and Ser showed no interaction. Stepwise regression analysis indicated that the pK value of the amino acid side-chain and the lipophilicity of the amino acid had the greatest impact on the interaction, suggesting the simultaneous presence of weak hydrophilic and hydrophobic bonding forces between amino acids. Sodium acetate in the eluent increased the interactive strength between Phe and Trp; acetic acid and sodium chloride did not influence the interaction significantly. No significant difference was found between the effects of L- and D-Asn.