Superimposing the 27 crystal protein/inhibitor complexes of β-secretase to calculate the binding affinities by the linear interaction energy method

G-score: A function to solve the puzzle of modeling the protonation states of β-secretase binding pocket

The population density concept has emerged as a proposal for the analysis of molecular dynamics results, the key characteristic of population density is the evaluation of the simultaneous occurrence of a set of relevant parameters for a system. However, despite its statistical strength, selection of the tolerance level for the comparison of different models may appear as arbitrary. This work introduces the G-score, a function which summarizes and categorizes the results of population density analysis. Additionally, it incorporates parameters based on rmsd and dihedral angles, besides the protein-protein and protein-ligand interatomic distances conventionally used, which complement each other to provide a better description of the behavior of the system. These newly-proposed tools were applied to determine the most probable protonation state of the aspartic dyad of BACE1, Asp93 and Asp289, in the presence of three types of transition state inhibitors namely: reduced amides, tertiary carbinamines and hydroxyethylamines. The results show a full agreement between G-score values and population density charts, with the advantage of allowing a quick and direct comparison among all the considered models. We anticipate that the simplicity of calculating the parameters employed in this study will permit the extensive use of population density and the G-score for other molecular systems.

A comparative study of binding affinities for 6,7-dimethoxy-4-pyrrolidylquinazolines as phosphodiesterase 10A inhibitors using the linear interaction energy method

The linear interaction energy (LIE) method was used to estimate the free energies of binding for a set of 27 pyrrolidylquinazoline derivatives as phosphodiesterase 10A inhibitors. Twenty-six X-ray crystal structures of phosphodiesterase 10A and two sampling methods, minimization and Hybrid Monte Carlo, were used to assess the affinity models based on the linear interaction energies. The best model was obtained based on the parameters α=0.16 and β=0.04, which represent non-polar and polar interactions, respectively, with a root mean square error (RMSE) of 0.42kcal/mol (R(2)=0.71) and 0.52kcal/mol (R(2)=0.86) for the training and test sets, respectively. In addition, the applicability domain of the model was investigated. After validation of the models, the best model was subsequently used in a virtual screening process, which resulted in a set of optimized compounds. The models developed in this study could be useful as filter for virtual screening and lead optimization processes for phosphodiesterase 10A drug developments.

Exploring the binding of BACE-1 inhibitors using comparative binding energy analysis (COMBINE)

BACKGROUND

The inhibition of the activity of β-secretase (BACE-1) is a potentially important approach for the treatment of Alzheimer disease. To explore the mechanism of inhibition, we describe the use of 46 X-ray crystallographic BACE-1/inhibitor complexes to derive quantitative structure-activity relationship (QSAR) models. The inhibitors were aligned by superimposing 46 X-ray crystallographic BACE-1/inhibitor complexes, and gCOMBINE software was used to perform COMparative BINding Energy (COMBINE) analysis on these 46 minimized BACE-1/inhibitor complexes. The major advantage of the COMBINE analysis is that it can quantitatively extract key residues involved in binding the ligand and identify the nature of the interactions between the ligand and receptor.

RESULTS

By considering the contributions of the protein residues to the electrostatic and van der Waals intermolecular interaction energies, two predictive and robust COMBINE models were developed: (i) the 3-PC distance-dependent dielectric constant model (built from a single X-ray crystal structure) with a q2 value of 0.74 and an SDEC value of 0.521; and (ii) the 5-PC sigmoidal electrostatic model (built from the actual complexes present in the Brookhaven Protein Data Bank) with a q2 value of 0.79 and an SDEC value of 0.41.

CONCLUSIONS

These QSAR models and the information describing the inhibition provide useful insights into the design of novel inhibitors via the optimization of the interactions between ligands and those key residues of BACE-1.

Role of computational methods in pharmaceutical sciences

Over the past two decades computational methods have eased up the financial and experimental burden of early drug discovery process. The in silico methods have provided support in terms of databases, data mining of large genomes, network analysis, systems biology on the bioinformatics front and structure-activity relationship, similarity analysis, docking, and pharmacophore methods for lead design and optimization. This review highlights some of the applications of bioinformatics and chemoinformatics methods that have enriched the field of drug discovery. In addition, the review also provided insights into the use of free energy perturbation methods for efficiently computing binding energy. These in silico methods are complementary and can be easily integrated into the traditional in vitro and in vivo methods to test pharmacological hypothesis.