Computational study of bindings of olive leaf extract (OLE) to HIV-1 fusion protein gp41

Recent experimental study found that OLE (olive leaf extract) has anti-HIV activity by blocking the HIV virus entry to host cells [Lee-Huang, S., Zhang, L., Huang, P.L., Chang, Y. and Huang, P.L. (2003) Anti-HIV activity of olive leaf extract (OLE) and modulation of host cell gene expression by HIV-1 infection and OLE treatment. Biochem. Biophys. Res. Commun. 307, 1029; Lee-Huang, S., Huang, P.L., Zhang, D., Lee, J.W., Bao, J., Sun, Y., Chang, Y.-Tae, Zhang, J.Z.H. and Huang, P.L. (2007) Discovery of small-molecule HIV-1 fusion and integrase inhibitors oleuropein and hydroxytyrosol. Biochem. Biophys. Res. Commun. 354, 872-878, 879-884]. As part of a joint experimental and theoretical effort, we report here computational study to help identify and characterize the binding complexes of several main compounds of OLE (olive leaf extract) to HIV-1 envelop protein gp41. A number of possible binding modes are found by docking oleuropein and its metabolites, aglycone, elenolic acid and hydroxytyrosol, onto the hydrophobic pocket on gp41. Detailed OLE-gp41 binding interactions and free energies of binding are obtained through molecular dynamics simulation and MM-PBSA calculation. Specific molecular interactions in our predicted OLE/gp41 complexes are identified and hydroxytyrosol is identified to be the main moiety for binding to gp41. This computational study complements the corresponding experimental investigation and helps establish a good starting point for further refinement of OLE-based gp41 inhibitors.     

Binding free energies for nicotine analogs inhibiting cytochrome P450 2A6 by a combined use of molecular dynamics simulations and QM/MM-PBSA calculations

Molecular dynamics (MD) simulations and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations have been performed to explore the dynamic behaviors of cytochrome P450 2A6 (CYP2A6) binding with nicotine analogs (that are typical inhibitors) and to calculate their binding free energies in combination with Poisson–Boltzmann surface area (PBSA) calculations. The combined MD simulations and QM/MM-PBSA calculations reveal that the most important structural parameters affecting the CYP2A6-inhibitor binding affinity are two crucial internuclear distances, that is, the distance between the heme iron atom of CYP2A6 and the coordinating atom of the inhibitor, and the hydrogen-bonding distance between the N297 side chain of CYP2A6 and the pyridine nitrogen of the inhibitor. The combined MD simulations and QM/MM-PBSA calculations have led to dynamic CYP2A6-inhibitor binding structures that are consistent with the observed dynamic behaviors and structural features of CYP2A6-inhibitor binding, and led to the binding free energies that are in good agreement with the experimentally-derived binding free energies. The agreement between the calculated binding free energies and the experimentally-derived binding free energies suggests that the combined MD and QM/MM-PBSA approach may be used as a valuable tool to accurately predict the CYP2A6-inhibitor binding affinities in future computational design of new, potent and selective CYP2A6 inhibitors.     

Insights into the functional role of protonation states in the HIV-1 protease-BEA369 complex: molecular dynamics simulations and free energy calculations

The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method combined with molecular dynamics (MD) simulations were used to investigate the functional role of protonation in human immunodeficiency virus type 1 (HIV-1) protease complexed with the inhibitor BEA369. Our results demonstrate that protonation of two aspartic acids (Asp25/Asp25′) has a strong influence on the dynamics behavior of the complex, the binding free energy of BEA369, and inhibitor–residue interactions. Relative binding free energies calculated using the MM-PBSA method show that protonation of Asp25 results in the strongest binding of BEA369 to HIV-1 protease. Inhibitor–residue interactions computed by the theory of free energy decomposition also indicate that protonation of Asp25 has the most favorable effect on binding of BEA369. In addition, hydrogen-bond analysis based on the trajectories of the MD simulations shows that protonation of Asp25 strongly influences the water-mediated link of a conserved water molecule, Wat301. We expect that the results of this study will contribute significantly to binding calculations for BEA369, and to the design of high affinity inhibitors.     

Study on the binding mode of a pyrrolotriazin derivative with JAK2 by docking and MD simulation

The pyrrolotriazin derivative 2-(4-(4-((7-(3-(N-methylmethylsulfonamido)phenyl)pyrrolo [2,1-f][1,2,4]triazin-2-yl)amino)phenyl)piperidin-1-yl)acetamide (PPA) is a potential Janus kinase 2 (JAK2) inhibitor. The binding mode between PPA and JAK2 was investigated by using a combined method of docking, molecular dynamics (MD) simulation and binding free-energy calculation. The docking calculations preliminarily indicated that there were two possible binding modes 1 and 2; MD simulations and binding free-energy calculations identified that binding mode 1 was more stable and favourable, with the lower MM-PBSA binding free energy of ?34.00?±?0.17?kcal/mol. Moreover, some valuable binding information is revealed as follows: the inhibitor PPA is suitably located at the ATP-binding site of JAK2 and the hydrophobic interaction plays an essential role. PPA not only interacts with residues Leu855, Val863, Ala880, Tyr931, Leu932 and Leu983 via hydrophobic interaction but also interacts with Ser936 and Asp994 by hydrogen bonds. These two factors are advantageous for PPA to strongly bind to JAK2. These results help to understand the action mechanisms and designing new compounds with a higher affinity to JAK2.     

Models for the binary complex of bacteriophage T4 gp59 helicase loading protein: gp32 single-stranded DNA-BINDING protein and ternary complex with pseudo-Y junction DNA

Bacteriophage T4 gp59 helicase assembly protein (gp59) is required for loading of gp41 replicative helicase onto DNA protected by gp32 single-stranded DNA-binding protein. The gp59 protein recognizes branched DNA structures found at replication and recombination sites. Binding of gp32 protein (full-length and deletion constructs) to gp59 protein measured by isothermal titration calorimetry demonstrates that the gp32 protein C-terminal A-domain is essential for protein-protein interaction in the absence of DNA. Sedimentation velocity experiments with gp59 protein and gp32ΔB protein (an N-terminal B-domain deletion) show that these proteins are monomers but form a 1:1 complex with a dissociation constant comparable with that determined by isothermal titration calorimetry. Small angle x-ray scattering (SAXS) studies indicate that the gp59 protein is a prolate monomer, consistent with the crystal structure and hydrodynamic properties determined from sedimentation velocity experiments. SAXS experiments also demonstrate that gp32ΔB protein is a prolate monomer with an elongated A-domain protruding from the core. Fitting structures of gp59 protein and the gp32 core into the SAXS-derived molecular envelope supports a model for the gp59 protein-gp32ΔB protein complex. Our earlier work demonstrated that gp59 protein attracts full-length gp32 protein to pseudo-Y junctions. A model of the gp59 protein-DNA complex, modified to accommodate new SAXS data for the binary complex together with mutational analysis of gp59 protein, is presented in the accompanying article (Dolezal, D., Jones, C. E., Lai, X., Brister, J. R., Mueser, T. C., Nossal, N. G., and Hinton, D. M. (2012) J. Biol. Chem. 287, 18596-18607).     

The HSP90 binding mode of a radicicol-like E-oxime determined by docking,binding free energy estimations,and NMR N chemical shifts

We determine the binding mode of a macrocyclic radicicol-like oxime to yeast HSP90 by combining computer simulations and experimental measurements. We sample the macrocyclic scaffold of the unbound ligand by parallel tempering simulations and dock the most populated conformations to yeast HSP90. Docking poses are then evaluated by the use of binding free energy estimations with the linear interaction energy method. Comparison of QM/MM-calculated NMR chemical shifts with experimental shift data for a selective subset of backbone ¹⁵N provides an additional evaluation criteria. As a final test we check the binding modes against available structure–activity-relationships. We find that the most likely binding mode of the oxime to yeast HSP90 is very similar to the known structure of the radicicol–HSP90 complex.     

Molecular dynamics simulation of interleukin-2 and its complex and determination of the binding free energy

Interleukin-2 (IL-2) protein belongs to the signal modulator cytokine's family and therefore it is prevalent for immunological responses. It has been identified as a centrally important potential drug target for the inhibition of protein-protein interactions; so as to suppress the immunological responses associated with autoimmune, inflammatory and immunological diseases, and cancer. In the present work, we have performed two independent 100?ns of molecular dynamics (MD) simulations on the apo IL-2 protein and its ligand-bound complex (with a potent inhibitor FRG), to study the effect of inhibitor binding on the dynamics and stability of the protein. The calculation of binding free energy via post-processing end state method of Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA) and Molecular Mechanics Generalised Born Surface Area (MM-GBSA) has inferred a good correlation in accordance with the already reported experimental data, demonstrating that the free energy of binding calculated by the two methods has no significant difference. The investigation of individual components of free energy revealed that the association of IL-2 protein with FRG ligand is primarily driven by the van der Waals energy contribution that represents the non-polar/hydrophobic energy contribution as dominant in this case of ligand binding.     

Validation of an automated procedure for the prediction of relative free energies of binding on a set of aldose reductase inhibitors

Among the available methods for predicting free energies of binding of ligands to a protein, the molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) and molecular mechanics generalized Born surface area (MM-GBSA) approaches have been validated for a relatively limited number of targets and compounds in the training set. Here, we report the results of an extensive study on a series of 28 inhibitors of aldose reductase with experimentally determined crystal structures and inhibitory activities, in which we evaluate the ability of MM-PBSA and MM-GBSA methods in predicting binding free energies using a number of different simulation conditions. While none of the methods proved able to predict absolute free energies of binding in quantitative agreement with the experimental values, calculated and experimental free energies of binding were significantly correlated. Comparing the predicted and experimental ΔG of binding, MM-PBSA proved to perform better than MM-GBSA, and within the MM-PBSA methods, the PBSA of Amber performed similarly to Delphi. In particular, significant relationships between experimental and computed free energies of binding were obtained using Amber PBSA and structures minimized with a distance-dependent dielectric function. Importantly, while free energy predictions are usually made on large collections of equilibrated structures sampled during molecular dynamics in water, we have found that a single minimized structure is a reasonable approximation if relative free energies of binding are to be calculated. This finding is particularly relevant, considering that the generation of equilibrated MD ensembles and the subsequent free energy analysis on multiple snapshots is computationally intensive, while the generation and analysis of a single minimized structure of a protein–ligand complex is relatively fast, and therefore suited for high-throughput virtual screening studies. At this aim, we have developed an automated workflow that integrates all the necessary steps required to generate structures and calculate free energies of binding. The procedure is relatively fast and able to screen automatically and iteratively molecules contained in databases and libraries of compounds. Taken altogether, our results suggest that the workflow can be a valuable tool for ligand identification and optimization, being able to automatically and efficiently refine docking poses, which sometimes may not be accurate, and rank the compounds based on more accurate scoring functions.     

Molecular docking study investigating the possible mode of binding of C.I. Acid Red 73 with DNA

C.I. Acid Red 73 is a reactive azo dye with a variable potential carcinogenicity. The mechanism mediating interactions that occur between the dye and DNA have not been completely understood thus far. In this study, molecular docking techniques were applied to describe the most probable mode of DNA binding as well as the sequence selectivity of the C.I. Acid Red 73 dye. These docking experiments revealed that the dye is capable of interacting with the minor groove of the DNA on the basis of its curved shape, which fits well with the topology of double-stranded DNA. In addition, the dye can bind selectively to the minor groove of the DNA by applying CGT sequence selectivity. Further, the minor groove can be recognized although DNA targets present intercalation gaps. However, intercalative binding can also occur when the DNA target possesses an appropriate intercalation gap. Compared with the other eight DNA sequences that were studied, the DNA dodecamer d(CGCGATATCGCG)₂ (PDB ID: 1DNE) presents a very favorable target for the binding of C.I. Acid Red 73 to the minor groove, with the lowest binding free energy −9.19 kcal/mol. Results reported from this study are expected to provide useful information for research involving further simulations of molecular dynamics and toxicology investigations of the dye.     

SpltMNPV日本分离株gp41的克隆表达及gp41和ph的进化分析

从斜纹夜蛾核型多角体病毒(Spodopteralituramulticapsidnucleopolyhedrovirus,SpltMNPV)日本分离株(C3)基因组中克隆了gp41基因。该基因编码区含993bp核苷酸,编码分子量为36.9kDa的多肽。将该基因克隆至原核表达载体pET28a,经IPTG诱导后在大肠杆菌BL21(DE3)中获得了表达。应用CLUSTAL程序分析表明,SpltMNPV日本株(C3)gp41的核苷酸序列和氨基酸序列与SpltMNPV中国G2株相似性最高,均达99.9%。用MEGA分别构建了基于gp41和ph的聚类分析图和分子进化树,发现它们具有相似的拓扑结构。将这两个基因序列结合在一起构建进化树,该树的结构与基于gp41的进化树相似。突变率分析显示gp41的突变率高于ph,这意味着在杆状病毒进化过程中,gp41和ph面临不同的选择压力。     

Density functional theory calculations on entire proteins for free energies of binding: Application to a model polar binding site

In drug optimization calculations, the molecular mechanics Poisson‐Boltzmann surface area (MM‐PBSA) method can be used to compute free energies of binding of ligands to proteins. The method involves the evaluation of the energy of configurations in an implicit solvent model. One source of errors is the force field used, which can potentially lead to large errors due to the restrictions in accuracy imposed by its empirical nature. To assess the effect of the force field on the calculation of binding energies, in this article we use large‐scale density functional theory (DFT) calculations as an alternative method to evaluate the energies of the configurations in a “QM‐PBSA” approach. Our DFT calculations are performed with a near‐complete basis set and a minimal parameter implicit solvent model, within the self‐consistent calculation, using the ONETEP program on protein–ligand complexes containing more than 2600 atoms. We apply this approach to the T4‐lysozyme double mutant L99A/M102Q protein, which is a well‐studied model of a polar binding site, using a set of eight small aromatic ligands. We observe that there is very good correlation between the MM and QM binding energies in vacuum but less so in the solvent. The relative binding free energies from DFT are more accurate than the ones from the MM calculations, and give markedly better agreement with experiment for six of the eight ligands. Furthermore, in contrast to MM‐PBSA, QM‐PBSA is able to correctly predict a nonbinder. Proteins 2014; 82:3335–3346. © 2014 Wiley Periodicals, Inc.     

Prediction of the binding mode between BMS-378806 and HIV-1 gp120 by docking and molecular dynamics simulation

BMS-378806 is a newly discovered small molecule that effectively blocks the binding of CD4 with gp120. The binding mode of this kind of inhibitor remains unknown. In this paper, AutoDock 3.0 in conjunction with molecular dynamics simulation, accommodating the receptor's flexibility, was used to explore the binding mode between BMS-378806 and gp120. Two structures, Mode I and Mode II, with the lowest docking energy were selected as different representative binding modes. The analysis of the results from the molecular dynamics simulation indicated that the binding of BMS-348806 in Mode II is more stable. The average structure of Mode II was analyzed and compared with the experimental data. The conclusion was that BMS-378806 inserts the azaindole ring deeply into the PHE43 cavity and makes contact with a number of residues in the cavity, on the cavity and near the cavity. This study benefits the understanding of the mechanism of this kind of inhibitor and may provide useful information for rational drug design.     

Studies on the binding of 5-N-methylated quindoline derivative to human telomeric G-quadruplex

Quindoline derivatives as telomeric quadruplex ligands have shown good biological activity for telomerase inhibition. In the present study, we used spectroscopic and calorimetric methods to investigate the interactions between a quindoline derivative (5-methyl-11-(2-morpholinoethylamino)-10-H-indolo-[3,2-b]quinolin-5-ium iodide, compound 1) and human telomeric G-quadruplex. The thermodynamic studies using isothermal titration calorimetry (ITC) indicated that their binding process was temperature-dependent and enthalpy–entropy co-driven. The significant negative heat capacity was obtained experimentally from the temperature dependence of enthalpy changes, which was consistent with that from theoretical calculation, and all suggesting significant hydrophobic contribution to the molecular recognition process. Based on the results from UV–vis, ITC, steady-state and time-resolved fluorescence, their binding mode was determined as two ligand molecules stacking on the quartets on both ends of the quadruplex. These results shed light on rational design and development of quindoline derivatives as G-quadruplex binding ligands.     

Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes

Absolute binding free energy calculations and free energy decompositions are presented for the protein-protein complexes H-Ras/C-Raf1 and H-Ras/RalGDS. Ras is a central switch in the regulation of cell proliferation and differentiation. In our study, we investigate the capability of the molecular mechanics (MM)-generalized Born surface area (GBSA) approach to estimate absolute binding free energies for the protein-protein complexes. Averaging gas-phase energies, solvation free energies, and entropic contributions over snapshots extracted from trajectories of the unbound proteins and the complexes, calculated binding free energies (Ras-Raf: -15.0(+/-6.3)kcal mol(-1); Ras-RalGDS: -19.5(+/-5.9)kcal mol(-1)) are in fair agreement with experimentally determined values (-9.6 kcal mol(-1); -8.4 kcal mol(-1)), if appropriate ionic strength is taken into account. Structural determinants of the binding affinity of Ras-Raf and Ras-RalGDS are identified by means of free energy decomposition. For the first time, computationally inexpensive generalized Born (GB) calculations are applied in this context to partition solvation free energies along with gas-phase energies between residues of both binding partners. For selected residues, in addition, entropic contributions are estimated by classical statistical mechanics. Comparison of the decomposition results with experimentally determined binding free energy differences for alanine mutants of interface residues yielded correlations with r(2)=0.55 and 0.46 for Ras-Raf and Ras-RalGDS, respectively. Extension of the decomposition reveals residues as far apart as 25A from the binding epitope that can contribute significantly to binding free energy. These "hotspots" are found to show large atomic fluctuations in the unbound proteins, indicating that they reside in structurally less stable regions. Furthermore, hotspot residues experience a significantly larger-than-average decrease in local fluctuations upon complex formation. Finally, by calculating a pair-wise decomposition of interactions, interaction pathways originating in the binding epitope of Raf are found that protrude through the protein structure towards the loop L1. This explains the finding of a conformational change in this region upon complex formation with Ras, and it may trigger a larger structural change in Raf, which is considered to be necessary for activation of the effector by Ras.     

人类免疫缺陷病毒1型gp41结构与功能的研究进展

人类免疫缺陷病毒1型(HIV-1)通过其包膜糖蛋白(Env)介导侵入靶细胞.Env由受体特异性结合单位gp120和膜融合单位gp41组成.HIV-1的gp41分为3个功能区:膜外区、跨膜区和膜内区.膜外区是病毒感染时膜融合的主要结构基础;跨膜区通过疏水残基使Env锚定在脂质膜上;膜内区则表现多重功能,参与病毒的感染、复...     

Insights into the binding mode of sulphamates and sulphamides to hCA II: crystallographic studies and binding free energy calculations

Sulphamate and sulphamide derivatives have been largely investigated as carbonic anhydrase inhibitors (CAIs) by means of different experimental techniques. However, the structural determinants responsible for their different binding mode to the enzyme active site were not clearly defined so far. In this paper, we report the X-ray crystal structure of hCA II in complex with a sulphamate inhibitor incorporating a nitroimidazole moiety. The comparison with the structure of hCA II in complex with its sulphamide analogue revealed that the two inhibitors adopt a completely different binding mode within the hCA II active site. Starting from these results, we performed a theoretical study on sulphamate and sulphamide derivatives, demonstrating that electrostatic interactions with residues within the enzyme active site play a key role in determining their binding conformation. These findings open new perspectives in the design of effective CAIs using the sulphamate and sulphamide zinc binding groups as lead compounds.     

Identification of fragments targeting an alternative pocket on HIV-1 gp41 by NMR screening and similarity searching

The HIV-1 envelope glycoprotein gp41 fusion intermediate is a promising drug target for inhibiting viral entry. However, drug development has been impeded by challenges inherent in mediating the underlying protein–protein interaction. Here we report on the identification of fragments that bind to a C-terminal sub-pocket adjacent to the well-known hydrophobic pocket on the NHR coiled coil. Using a specifically designed assay and ligand-based NMR screening of a fragment library, we identified a thioenylaminopyrazole compound with a dissociation constant of ～500 μM. Interaction with the C-terminal sub-pocket was confirmed by paramagnetic relaxation enhancement NMR experiments, which also yielded the binding mode. Shape-based similarity searching detected additional phenylpyrazole and phenyltriazole fragments within the library, enriching the hit rate over random screening, and revealing molecular features required for activity. Discovery of the novel scaffolds and binding mechanism suggests avenues for extending the interaction surface and improving the potency of a hydrophobic pocket binding inhibitor.     

Protein-inhibitor flexible docking by a multicanonical sampling: native complex structure with the lowest free energy and a free-energy barrier distinguishing the native complex from the others

Flexible docking between a protein (lysozyme) and an inhibitor (tri-N-acetyl-D-glucosamine, tri-NAG) was carried out by an enhanced conformational sampling method, multicanonical molecular dynamics simulation. We used a flexible all-atom model to express lysozyme, tri-NAG, and water molecules surrounding the two bio-molecules. The advantages of this sampling method are as follows: the conformation of system is widely sampled without trapping at energy minima, a thermally equilibrated conformational ensemble at an arbitrary temperature can be reconstructed from the simulation trajectory, and the thermodynamic weight can be assigned to each sampled conformation. During the simulation, exchanges between the binding and free (i.e., unbinding) states of the protein and the inhibitor were repeatedly observed. The conformational ensemble reconstructed at 300 K involved various conformational clusters. The main outcome of the current study is that the most populated conformational cluster (i.e., the cluster of the lowest free energy) was assigned to the native complex structure (i.e., the X-ray complex structure). The simulation also produced non-native complex structures, where the protein and the inhibitor bound with different modes from that of the native complex structure, as well as the unbinding structures. A free-energy barrier (i.e., activation free energy) was clearly detected between the native complex structures and the other structures. The thermal fluctuations of tri-NAG in the lowest free-energy complex correlated well with the X-ray B-factors of tri-NAG in the X-ray complex structure. The existence of the free-energy barrier ensures that the lowest free-energy structure can be discriminated naturally from the other structures. In other words, the multicanonical molecular dynamics simulation can predict the native complex structure without any empirical objective function. The current study also manifested that the flexible all-atom model and the physico-chemically defined atomic-level force field can reproduce the native complex structure. A drawback of the current method is that it requires a time consuming computation due to the exhaustive conformational sampling. We discussed a possibility for combining the current method with conventional docking methods.     

The free energies for mutating S27 and W79 to alanine in streptavidin and its biotin complex: the relative size of polar and nonpolar free energies on biotin binding.

We have carried out calculations on the relative free energy of binding of biotin and its S27A and W79A mutants to streptavidin. Consistent with earlier suggestions by Miyamoto and Kollman from free energy component analysis and recent experiments by Stayton and coworkers, the reduction in binding strength by the W79A mutant is significantly larger than that of the S27A mutant. Proteins 1999;36:471-473.