Molecular mechanics PBSA ligand binding energy and interaction of Efavirenz derivatives with HIV-1 reverse transcriptase

In order to evaluate the properties of several HIV-1 reverse transcripase(RT) inhibitors, Efavirenz (SUSTIVA) and a set of its derivatives (benzoxazinones) have been placed into the nonnucleoside analogue binding site of the enzyme by molecular docking. The resulting geometries were used for a molecular dynamics simulation and binding energy calculations. The enzyme-inhibitor binding energies were estimated from experimental inhibitory activities (IC90). The correlation of the predicted and experimental binding energies were satisfactory acceptable as indicated by r2 = 0.865. Based on MD simulations, the obtained results indicate that the tight association of the ligand to the HIV-1 RT binding pocket was based on hydrogen bonding between Efavirenz's N1 and the oxygen of the backbone of Lys 101, with an estimated average distance of 1.88 A. Moreover, electrostatic interaction was mainly contributed by two amino acid residues in the binding site; Lys 101 and His 235. MD simulations open the possibility to study the reaction of the flexible enzyme to those substances as well as the overall affinity.     

Impact of ligand and protein desolvation on ligand binding to the S1 pocket of thrombin

In the present study, we investigate the impact of a tightly bound water molecule on ligand binding in the S1 pocket of thrombin. The S1 pocket contains a deeply buried deprotonated aspartate residue (Asp189) that is, due to its charged state, well hydrated in the uncomplexed state. We systematically studied the importance of this water molecule by evaluating a series of ligands that contains pyridine-type P1 side chains that could potentially alter the binding properties of this water molecule. All of the pyridine derivatives retain the original hydration state albeit sometimes with a slight perturbance. In order to prevent a direct H-bond formation with Asp189, and to create a permanent positive charge on the P1 side chain that is positioned adjacent to the Asp189 carboxylate anion, we methylated the pyridine nitrogen. This methylation resulted in displacement of water but was accompanied by a loss in binding affinity. Quantum chemical calculations of the ligand solvation free energy showed that the positively charged methylpyridinium derivatives suffer a large penalty of desolvation upon binding. Consequently, they have a substantially less favorable enthalpy of binding. In addition to the ligand desolvation penalty, the hydration shell around Asp189 has to be overcome, which is achieved in nearly all pyridinium derivatives. Only for the ortho derivative is a partial population of a water next to Asp189 found. Possibly, the gain of electrostatic interactions between the charged P1 side chain and Asp189 helps to compensate for the desolvation penalty. In all uncharged pyridine derivatives, the solvation shell remains next to Asp189, partly mediating interactions between ligand and protein. In the case of the para-pyridine derivative, a strongly disordered cluster of water sites is observed between ligand and Asp189.     

Optimizing the affinity and specificity of ligand binding with the inclusion of solvation effect

Solvation effect is an important factor for protein–ligand binding in aqueous water. Previous scoring function of protein–ligand interactions rarely incorporates the solvation model into the quantification of protein–ligand interactions, mainly due to the immense computational cost, especially in the structure‐based virtual screening, and nontransferable application of independently optimized atomic solvation parameters. In order to overcome these barriers, we effectively combine knowledge‐based atom–pair potentials and the atomic solvation energy of charge‐independent implicit solvent model in the optimization of binding affinity and specificity. The resulting scoring functions with optimized atomic solvation parameters is named as specificity and affinity with solvation effect (SPA‐SE). The performance of SPA‐SE is evaluated and compared to 20 other scoring functions, as well as SPA. The comparative results show that SPA‐SE outperforms all other scoring functions in binding affinity prediction and “native” pose identification. Our optimization validates that solvation effect is an important regulator to the stability and specificity of protein–ligand binding. The development strategy of SPA‐SE sets an example for other scoring function to account for the solvation effect in biomolecular recognitions. Proteins 2015; 83:1632–1642. © 2015 Wiley Periodicals, Inc.     

Stacking-interactions in the control-gear binding of kynurenic acid with NR2A- and GluR2-subunits of glutamate ionotropic receptors

Kynurenine products in tryptophan metabolism are of crucial importance in modulation of neurodegenerative processes in the CNS. Kynurenic acid (KYNA): the endogenous antagonist of ionotropic glutamate receptors, displays specific affinity towards glycine site ofNMDA-receptor NR1 subunit. Mechanisms for the selective interaction of KYNA and its derivatives with other glutamate receptor subtypes are studied insufficiently. Ab initio quantum chemical calculations for KYNA-imidazole dimer, as a model for ligand interaction with His88 fragment of NR2A-subunit, along with KYNA-phenol dimer, as a model for ligand interaction with Tyr61 fragment of GluR2-subunit, were carried out in order to investigate stacking-interaction role of KYNA binding by NR2A subunit of NMDA-receptor and GluR2 subunit of AMPA-receptor. Stacking-interaction energy of KYNA-H88 for the assumed ligand orientation in the binding site is 3.0-5.0 kcal/mol and 102. kcal/mol for the optimized dimer KYNA-imidazole geometry. Stacking-interaction energy of KYNA-Tyr61 for the assumed ligand orientation in the binding site is 6.7-8.5 kcal/mol. The obtained values are comparable with the energies of hydrogen bonds. Thus, stacking-interaction should be taken into account while studing ligand glutamate receptor binding mechanisms. Stacking-interaction is evidently important for the initial ligand orientation inside the receptor binding site after which the delicate tuning of hydrogen bonding pattern is realized. On the other hand, the specific affinity of KYNA derivatives to the receptor subunits could be explained by ligand-aromatic receptor aminoacid stacking-interaction geometry difference.     

Receptor binding kinetics and cellular responses of six N-formyl peptide agonists in human neutrophils

The goal of this study was to elucidate the relationships between early ligand binding/receptor processing events and cellular responses for the N-formyl peptide receptor system on human neutrophils as a model of a GPCR system in a physiologically relevant context. Binding kinetics of N-formyl-methionyl-leucyl-phenylalanyl-phenylalanyl-lysine-fluorescein and N-formyl-valyl-leucyl-phenylalanyl-lysine-fluorescein to the N-formyl peptide receptor on human neutrophils were characterized and combined with previously published binding data for four other ligands. Binding was best fit by an interconverting two-receptor state model that included a low affinity receptor state that converted to a high affinity state. Response behaviors elicited at 37 degrees C by the six different agonists for the N-formyl peptide receptor were measured. Dose response curves for oxidant production, actin polymerization, and G-protein activation were obtained for each ligand; whereas all ligands showed equal efficacy for all three responses, the ED(50) values varied as much as 7000-fold. The level of agonism and rank order of potencies of ligands for actin and oxidant responses were the same as for the G-protein activation assay, suggesting that the differences in abilities of ligands to mediate responses were determined upstream of G-protein activation at the level of ligand-receptor interactions. The rate constants governing ligand binding and receptor affinity conversion were ligand-dependent. Analysis of the forward and reverse rate constants governing binding to the proposed signaling receptor state showed that it was of a similar energy for all six ligands, suggesting the hypothesis that ligand efficacy is dictated by the energy state of this ligand-receptor complex. However, the interconverting two-receptor state model was not sufficient to predict response potency, suggesting the presence of receptor states not discriminated by the binding data.