Tingjun Hou

Assessing the performance of the MM/PBSA and MM/GBSA methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations

The Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) and the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) methods calculate binding free energies for macromolecules by combining molecular mechanics calculations and continuum solvation models. To systematically evaluate the performance of these methods, we report here an extensive study of 59 ligands interacting with six different proteins. First, we explored the effects of the length of the molecular dynamics (MD) simulation, ranging from 400 to 4800 ps, and the solute dielectric constant (1, 2, or 4) on the binding free energies predicted by MM/PBSA. The following three important conclusions could be observed: (1) MD simulation length has an obvious impact on the predictions, and longer MD simulation is not always necessary to achieve better predictions. (2) The predictions are quite sensitive to the solute dielectric constant, and this parameter should be carefully determined according to the characteristics of the protein/ligand binding interface. (3) Conformational entropy often show large fluctuations in MD trajectories, and a large number of snapshots are necessary to achieve stable predictions. Next, we evaluated the accuracy of the binding free energies calculated by three Generalized Born (GB) models. We found that the GB model developed by Onufriev and Case was the most successful model in ranking the binding affinities of the studied inhibitors. Finally, we evaluated the performance of MM/GBSA and MM/PBSA in predicting binding free energies. Our results showed that MM/PBSA performed better in calculating absolute, but not necessarily relative, binding free energies than MM/GBSA. Considering its computational efficiency, MM/GBSA can serve as a powerful tool in drug design, where correct ranking of inhibitors is often emphasized.

Assessing the performance of the molecular mechanics/Poisson Boltzmann surface area and molecular mechanics/generalized Born surface area methods. II. The accuracy of ranking poses generated from docking

In molecular docking, it is challenging to develop a scoring function that is accurate to conduct high‐throughput screenings. Most scoring functions implemented in popular docking software packages were developed with many approximations for computational efficiency, which sacrifices the accuracy of prediction. With advanced technology and powerful computational hardware nowadays, it is feasible to use rigorous scoring functions, such as molecular mechanics/Poisson Boltzmann surface area (MM/PBSA) and molecular mechanics/generalized Born surface area (MM/GBSA) in molecular docking studies. Here, we systematically investigated the performance of MM/PBSA and MM/GBSA to identify the correct binding conformations and predict the binding free energies for 98 protein–ligand complexes. Comparison studies showed that MM/GBSA (69.4%) outperformed MM/PBSA (45.5%) and many popular scoring functions to identify the correct binding conformations. Moreover, we found that molecular dynamics simulations are necessary for some systems to identify the correct binding conformations. Based on our results, we proposed the guideline for MM/GBSA to predict the binding conformations. We then tested the performance of MM/GBSA and MM/PBSA to reproduce the binding free energies of the 98 protein–ligand complexes. The best prediction of MM/GBSA model with internal dielectric constant 2.0, produced a Spearmans correlation coefficient of 0.66, which is better than MM/PBSA (0.49) and almost all scoring functions used in molecular docking. In summary, MM/GBSA performs well for both binding pose predictions and binding free‐energy estimations and is efficient to re‐score the top‐hit poses produced by other less‐accurate scoring functions.

Assessing the performance of MM/PBSA and MM/GBSA methods. 3. The impact of force fields and ligand charge models.

Here, we systematically investigated how the force fields and the partial charge models for ligands affect the ranking performance of the binding free energies predicted by the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) approaches. A total of 46 small molecules targeted to five different protein receptors were employed to test the following issues: (1) the impact of five AMBER force fields (ff99, ff99SB, ff99SB-ILDN, ff03, and ff12SB) on the performance of MM/GBSA, (2) the influence of the time scale of molecular dynamics (MD) simulations on the performance of MM/GBSA with different force fields, (3) the impact of five AMBER force fields on the performance of MM/PBSA, and (4) the impact of four different charge models (RESP, ESP, AM1-BCC, and Gasteiger) for small molecules on the performance of MM/PBSA or MM/GBSA. Based on our simulation results, the following important conclusions can be obtained: (1) for short time-scale MD simulations (1 ns or less), the ff03 force field gives the best predictions by both MM/GBSA and MM/PBSA; (2) for middle time-scale MD simulations (2-4 ns), MM/GBSA based on the ff99 force field yields the best predictions, while MM/PBSA based on the ff99SB force field does the best; however, longer MD simulations, for example, 5 ns or more, may not be quite necessary; (3) for most cases, MM/PBSA with the Tans parameters shows better ranking capability than MM/GBSA (GB(OBC1)); (4) the RESP charges show the best performance for both MM/PBSA and MM/GBSA, and the AM1-BCC and ESP charges can also give fairly satisfactory predictions. Our results provide useful guidance for the practical applications of the MM/GBSA and MM/PBSA approaches.

Recent Advances in Free Energy Calculations with a Combination of Molecular Mechanics and Continuum Models

Recently, the combination of state-of-the-art molecular mechanical force fields with continuum solvation models enables us to make relatively accurate predictions of both structures and free energies for macromolecules from molecular dynamics trajectories. The first part of this review is focused on the history and basic theory of free energy calculations based on physically effective energy functions. The second part illustrates the applications of free energy calculations on many biological systems, including proteins, DNA, RNA, protein-ligand, protein-protein, protein-nucleic acid complexes, etc. Finally, the prospective and possible strategies to improve the techniques of MM-PBSA and MM-GBSA is discussed. Molecular mechanics (MM), though its theoretical background is not as solid as quantum mechanics (QM), has broad applications in studying biological systems for its simplicity and efficiency. The harmonic function form (Equation 1), which is widely-used in many popular molecular mechanical force fields, describes molecular energy using bond stretching, bond angle bending, torsional angle twisting as well as non-bonded electrostatic and van der Waals interactions. It is well-known that most biological procedures take place in aqueous solutions. Therefore, solvation effect cannot be neglected in studying the structures and interactions of biological systems. There are two basic ways to take the solvent effect into account: with either an explicit water model or an implicit water model. For the first approach, a biological system is usually immersed in a periodical water box; the second approach does not apply water molecules explicitly. Instead, solvation energy and force due to the solvent effect are calculated using some formulas or by solving some equations numerically or analytically. The widely used implicit solvation models include the generalized Born surface area (GBSA) (1-5) and the Poisson-Bolzmann surface area (PBSA) (6-15).

Characterization of domain-peptide interaction interface: a case study on the amphiphysin-1 SH3 domain.

Many important protein-protein interactions are mediated by peptide recognition modular domains, such as the Src homology 3 (SH3), SH2, PDZ, and WW domains. Characterizing the interaction interface of domain-peptide complexes and predicting binding specificity for modular domains are critical for deciphering protein-protein interaction networks. Here, we propose the use of an energetic decomposition analysis to characterize domain-peptide interactions and the molecular interaction energy components (MIECs), including van der Waals, electrostatic, and desolvation energy between residue pairs on the binding interface. We show a proof-of-concept study on the amphiphysin-1 SH3 domain interacting with its peptide ligands. The structures of the human amphiphysin-1 SH3 domain complexed with 884 peptides were first modeled using virtual mutagenesis and optimized by molecular mechanics (MM) minimization. Next, the MIECs between domain and peptide residues were computed using the MM/generalized Born decomposition analysis. We conducted two types of statistical analyses on the MIECs to demonstrate their usefulness for predicting binding affinities of peptides and for classifying peptides into binder and non-binder categories. First, combining partial least squares analysis and genetic algorithm, we fitted linear regression models between the MIECs and the peptide binding affinities on the training data set. These models were then used to predict binding affinities for peptides in the test data set; the predicted values have a correlation coefficient of 0.81 and an unsigned mean error of 0.39 compared with the experimentally measured ones. The partial least squares-genetic algorithm analysis on the MIECs revealed the critical interactions for the binding specificity of the amphiphysin-1 SH3 domain. Next, a support vector machine (SVM) was employed to build classification models based on the MIECs of peptides in the training set. A rigorous training-validation procedure was used to assess the performances of different kernel functions in SVM and different combinations of the MIECs. The best SVM classifier gave satisfactory predictions for the test set, indicated by average prediction accuracy rates of 78% and 91% for the binding and non-binding peptides, respectively. We also showed that the performance of our approach on both binding affinity prediction and binder/non-binder classification was superior to the performances of the conventional MM/Poisson-Boltzmann solvent-accessible surface area and MM/generalized Born solvent-accessible surface area calculations. Our study demonstrates that the analysis of the MIECs between peptides and the SH3 domain can successfully characterize the binding interface, and it provides a framework to derive integrated prediction models for different domain-peptide systems.

Characterization of Domain–Peptide Interaction Interface: Prediction of SH3 Domain-Mediated Protein–Protein Interaction Network in Yeast by Generic Structure-Based Models

Determination of the binding specificity of SH3 domain, a peptide recognition module (PRM), is important to understand their biological functions and reconstruct the SH3-mediated protein-protein interaction network. In the present study, the SH3-peptide interactions for both class I and II SH3 domains were characterized by the intermolecular residue-residue interaction network. We developed generic MIEC-SVM models to infer SH3 domain-peptide recognition specificity that achieved satisfactory prediction accuracy. By investigating the domain-peptide recognition mechanisms at the residue level, we found that the class-I and class-II binding peptides have different binding modes even though they occupy the same binding site of SH3. Furthermore, we predicted the potential binding partners of SH3 domains in the yeast proteome and constructed the SH3-mediated protein-protein interaction network. Comparison with the experimentally determined interactions confirmed the effectiveness of our approach. This study showed that our sophisticated computational approach not only provides a powerful platform to decipher protein recognition code at the molecular level but also allows identification of peptide-mediated protein interactions at a proteomic scale. We believe that such an approach is general to be applicable to other domain-peptide interactions.

Recent Advances in Computational Prediction of Drug Absorption and Permeability in Drug Discovery

Approximately 40%-60% of developing drugs failed during the clinical trials because of ADME/Tox deficiencies. Virtual screening should not be restricted to optimize binding affinity and improve selectivity; and the pharmacokinetic properties should also be included as important filters in virtual screening. Here, the current development in theoretical models to predict drug absorption-related properties, such as intestinal absorption, Caco-2 permeability, and blood-brain partitioning are reviewed. The important physicochemical properties used in the prediction of drug absorption, and the relevance of predictive models in the evaluation of passive drug absorption are discussed. Recent developments in the prediction of drug absorption, especially with the application of new machine learning methods and newly developed software are also discussed. Future directions for research are outlined.

Recent Development and Application of Virtual Screening in Drug Discovery: An Overview

Virtual screening, especially the structure-based virtual screening, has emerged as a reliable, cost-effective and time-saving technique for the discovery of lead compounds. Here, the basic ideas and computational tools for virtual screening have been briefly introduced, and emphasis is placed on aspects of recent development of docking-based virtual screening, scoring functions in molecular docking and ADME/Tox-based virtual screening in the past three years (2000 to 2003). Moreover, successful examples are provided to further demonstrate the effectiveness of virtual screening in drug discovery.

Computational analysis and prediction of the binding motif and protein interacting partners of the Abl SH3 domain

Protein-protein interactions, particularly weak and transient ones, are often mediated by peptide recognition domains, such as Src Homology 2 and 3 (SH2 and SH3) domains, which bind to specific sequence and structural motifs. It is important but challenging to determine the binding specificity of these domains accurately and to predict their physiological interacting partners. In this study, the interactions between 35 peptide ligands (15 binders and 20 non-binders) and the Abl SH3 domain were analyzed using molecular dynamics simulation and the Molecular Mechanics/Poisson-Boltzmann Solvent Area method. The calculated binding free energies correlated well with the rank order of the binding peptides and clearly distinguished binders from non-binders. Free energy component analysis revealed that the van der Waals interactions dictate the binding strength of peptides, whereas the binding specificity is determined by the electrostatic interaction and the polar contribution of desolvation. The binding motif of the Abl SH3 domain was then determined by a virtual mutagenesis method, which mutates the residue at each position of the template peptide relative to all other 19 amino acids and calculates the binding free energy difference between the template and the mutated peptides using the Molecular Mechanics/Poisson-Boltzmann Solvent Area method. A single position mutation free energy profile was thus established and used as a scoring matrix to search peptides recognized by the Abl SH3 domain in the human genome. Our approach successfully picked ten out of 13 experimentally determined binding partners of the Abl SH3 domain among the top 600 candidates from the 218,540 decapeptides with the PXXP motif in the SWISS-PROT database. We expect that this physical-principle based method can be applied to other protein domains as well.

Adme evaluation in drug discovery. 7. Prediction of oral absorption by correlation and classification

A critically evaluated database of human intestinal absorption for 648 chemical compounds is reported in this study, among which 579 are believed to be transported by passive diffusion. The correlation analysis between the intestinal absorption and several important molecular properties demonstrated that no single molecular property could be used as a good discriminator to efficiently distinguish the poorly absorbed compounds from those that are well absorbed. The theoretical correlation models for a training set of 455 compounds were proposed by using the genetic function approximation technique. The best prediction model contains four molecular descriptors: topological polar surface area, the predicted distribution coefficient at pH = 6.5, the number of violations of the Lipinskis rule-of-five, and the square of the number of hydrogen-bond donors. The model was able to predict the fractional absorption with an r = 0.84 and a prediction error (absolute mean error) of 11.2% for the training set. Moreover, it achieves an r = 0.90 and a prediction error of 7.8% for a 98-compound test set. The recursive partitioning technique was applied to find the simple hierarchical rules to classify the compounds into poor (%FA < or = 30%) and good (%FA > 30%) intestinal absorption classes. The high quality of the classification model was validated by the satisfactory predictions on the training set (correctly identifying 95.9% of the compounds in the poor-absorption class and 96.1% of the compounds in the good-absorption class) and on the test set (correctly identifying 100% of the compounds in the poor-absorption class and 96.8% of the compounds in the good-absorption class). We expect that, in the future, the rules for the prediction of carrier-mediated transporting and first pass metabolism can be integrated into the current hierarchical rules, and the classification model may become more powerful in the prediction of intestinal absorption or even human bioavailability. The databases of human intestinal absorption reported here are available for download from the supporting Web site: http://modem.ucsd.edu/adme.