Abstract
The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.
Introduction
It is widely accepted that high-level quantum mechanical (QM) methods provide the most detailed and accurate description of molecular structures, dynamics, and functions. However, for many biochemical systems that are often too complex, and/or biochemical processes that are too long, classical approaches are more commonly employed due to their efficiency and reasonable accuracy. To model biochemical systems classically, both long-range polar and short-range non-polar interactions are important for accurate and transferrable models (Perutz, 1978; Davis and McCammon, 1990; Honig and Nicholls, 1995). Among the classical simulation methods, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach (Srinivasan et al., 1998; Kollman et al., 2000; Gohlke and Case, 2004; Yang et al., 2011; Miller et al., 2012; Wang C. H. et al., 2016) has emerged as an efficient and reliable method to model molecular recognition, such as for protein-ligand binding interactions.
The development of binding free energy calculation methods has been a central focus in molecular simulations. Indeed, theoretically rigorous, but computationally expensive, free energy perturbation, and thermodynamic integration methods were both proposed much earlier than the MMPBSA method (Zwanzig, 1954; Bennett, 1976; Straatsma and McCammon, 1991). For simple and small systems where they can be reliably applied, these “exact” methods have been shown to be more accurate than the MMPBSA method. However, their applications to typically large and complex biomolecular recognition problems are quite limited. This is due to low efficiency and slow convergence. There are several key approximations utilized in the MMPBSA method that allow it to be used as an efficient and reasonable approximation for free energy simulations. The PBSA model is used so that the solvation contribution to the free energy is approximated by using a continuum solvent model. In addition, the method separately approximates enthalpy and entropy contributions to the free energy, with varying degrees of success as discussed below. There are other approximate methods that are similar to MMPBSA, such as the Mining Minima (M2) method and the Linear Interaction Energy (LIE) method. These are all categorized as end-point methods that focus on the end states of processes to compute the free energy changes (Head et al., 1997; Luo et al., 1999, 2001; Luo and Gilson, 2000; Aqvist and Marelius, 2001; Chen et al., 2004; Chang et al., 2007; Moghaddam et al., 2011; Mikulskis et al., 2012; Muddana and Gilson, 2012; Muddana et al., 2014).
In this review, we focus on recent developments and applications of the MMPBSA method that have been reported since 2014, which was roughly the time when the last major MMPBSA review was published (Genheden and Ryde, 2015). In the following, we first review the improvements made to the MMPBSA method over the last few years. This is followed by recent applications of the method in various fields, with a focus on biomedical applications.
Improvements of MMPBSA
Overview of MMPBSA
The MMPBSA method is most often applied to the calculation of binding free energies (ΔGbind) of small molecule ligands bound to large biomolecule receptors, although large inter-biomolecular recognitions are also often reported as reviewed below. The binding free energy of the bound ligand-receptor complex in an aqueous solvent (ΔGbind, aq) can be approximated as (Srinivasan et al., 1998; Kollman et al., 2000):
where ΔEMM, ΔGbind, solv, and − TΔS represent the gas-phase molecular mechanical energy change, the solvation free energy change, and the conformational entropy change upon binding, respectively. All of these changes are computed via ensemble averaging over a large set of sampled conformations. ΔEMM includes three terms calculated using molecular mechanics (MM): the covalent energy change (ΔEcovalent), the electrostatic energy change (ΔEelectrostatic), and the van der Waals energy change (ΔEvdW). ΔEcovalent consists of changes in the bond terms (ΔEbond), the angle terms (ΔEangle), and the torsion terms (ΔEtorsion). The solvation free energy change (ΔGbind, solv) is usually separated into polar and non-polar contributions (ΔGpolar and ΔGnon-polar). The entropy term is the most difficult to compute, and it is often approximated with a normal mode method using a few selected snapshots.
To compute all the ensemble averages in Equations (1–4) for a binding affinity calculation, the MMPBSA method usually begins with a molecular dynamics (MD) simulation of the complex using the single-trajectory approach, or three separate MD simulations of the complex, receptor, and ligand, respectively, in the multi-trajectory approach. A snapshot of a structure is taken at various time points during the production portion of the MD simulations. These snapshots are then used to calculate average values and uncertainties of various quantities of interest. The MD simulations are almost always conducted in an explicit solvent model to obtain the most accurate snapshots possible before carrying out any calculations that make use of them. It is important to obtain many different conformations, or snapshots, over a suitable timeframe for use in later statistical analysis as it is often not trivial to observe converged averaging, even in the relatively easier single-trajectory approach (Wang C. H. et al., 2016).
The polar solvation term
For the calculation of the solvation free energy change (ΔGbind, solv), the explicit solvent is removed and replaced with an implicit continuum solvent to greatly speed up the calculation time. The polar solvation term (ΔGpolar) in Equation (4) is calculated using a finite-difference solution or a Generalized Born (GB) pairwise approximation of the Poisson-Boltzmann equation (PBE) (Warwicker and Watson, 1982; Bashford and Karplus, 1990; Jeancharles et al., 1991; Gilson, 1995; Edinger et al., 1997; Luo et al., 1997, 2002; Lu and Luo, 2003; Tan et al., 2006; Cai et al., 2009, 2010, 2011; Wang et al., 2009, 2010, 2012, 2013; Ye et al., 2009, 2010; Wang and Luo, 2010; Hsieh and Luo, 2011; Botello-Smith et al., 2012; Liu et al., 2013; Wang C. H. et al., 2017). In addition to modeling binding interactions, PBE methods have also been applied to the prediction of pKa values for ionizable groups in biomolecules (Luo et al., 1998; Georgescu et al., 2002; Nielsen and McCammon, 2003; Warwicker, 2004; Tang et al., 2007), solvation free energies (Nicholls et al., 2008; Shivakumar et al., 2009), and protein folding and design (Hsieh and Luo, 2004; Wen and Luo, 2004; Wen et al., 2004; Lwin and Luo, 2005, 2006; Marshall et al., 2005; Lwin et al., 2006; Korman et al., 2008; Tan and Luo, 2008, 2009).
The PBE is based on the more fundamental Poisson equation:
where is a predefined dielectric distribution function for the solvated molecular system, is the potential distribution function, and is the fixed atomic charge density. To model the electrostatic interaction due to the additional presence of salt in the aqueous solution, the electrostatic potential can be described by the PBE as:
with
Here, is a predefined ion-exclusion function with a value of 0 within the Stern layer and the molecular interior and a value of 1 outside the Stern layer. The salt-related term is a function of the potential, the valence, zi, of ion type i, and the bulk concentration, ci, at a given temperature T. When both the ionic strength and electric field are weak, the PBE can be linearized for easier numerical solutions:
where . Here εν denotes the solvent dielectric constant, and I represents the ionic strength of the solution.
Over the past few years, a few new algorithm developments were reported for the numerical solution of the PBE (Xie, 2014; Fisicaro et al., 2016; Xie and Jiang, 2016). To deal with the singularity and nonlinearity of the PBE, Xie proposed a new decomposition and minimization scheme, together with a new proof on the existence and uniqueness of the PBE solution. A new PBE finite element solver was developed based on these solution decomposition and minimization techniques (Xie, 2014). Fisicaro et al. presented a preconditioned conjugate gradient technique to solve the generalized Poisson problem, and the linear regime of the PBE, in some 10 iterations. In combination with a self-consistent procedure, this technique was able to solve the non-linear Poisson–Boltzmann problem in a formulation including ionic steric effects (Fisicaro et al., 2016). Later Xie et al. incorporated nonlocal dielectric effects into the classic PBE for a protein in ionic solvent to derive a nonlocal modified Poisson–Boltzmann equation (NMPBE) and developed a finite element algorithm with a related package for solving the NMPBE (Xie and Jiang, 2016). Their results demonstrate the potential for the NMPBE to be a better predictor of electrostatic solvation and binding free energies compared to the standard PBE. It is worth noting that there has been a community wide push to explore alternative hardware for biomolecular simulations, such as the graphics processing units (GPU), which have a parallel architecture and are suited for high-performance computation with dense data parallelism (Colmenares et al., 2014a,b; Qi R. et al., 2017). A finite difference scheme with the successive over-relaxation method was implemented on the CUDA-based GPUs in the DelPhi package, which achieved a speedup of ~10 times in the linear and non-linear cases (Colmenares et al., 2014b). More recently, Qi et al. implemented and analyzed commonly used linear PBE solvers on CUDA GPUs for biomolecular simulations, including both standard and preconditioned conjugate gradient (CG) solvers with several alternative preconditioners (Qi R. et al., 2017). After extensive testing, the optimal GPU performance was observed using the Jacobi-preconditioned CG solver with a significant speedup that was up to 50 times faster than the standard CG solver on CPU. These progressive efforts on efficient numerical PBE solvers show great potential for accelerating MMPBSA computation.
Since the prior review (Genheden and Ryde, 2015), the numerical procedure and related factors for the widely used finite-difference method were also investigated for their impact on the MMPBSA method (Wang C. H. et al., 2016). This study showed that the impact of grid spacing on the quality of MMPBSA calculations is small in protein-ligand binding calculations; the agreement with experiment changed by a negligible amount when the grid spacing was changed from 0.50 to 0.25 Å. This indicated that the widely adopted default value of 0.50 Å used by the community was sufficient. The impact of different atomic radius sets and different molecular surface definitions was also analyzed, and weak influences were found on the agreement with experiment (Wang C. H. et al., 2016). This is probably due to the use of high protein dielectrics for the often-charged ligands and/or active sites as discussed below.
The effect of the solute dielectric constant was also investigated. A higher solute dielectric constant (using 2 or 4 instead of 1) was found to perform better in the virtual screening of ligands for tyrosine kinases (Sun et al., 2014a). Our own analysis of six groups of receptors reached a similar conclusion; the binding affinities using high dielectric constants (4 and 20) agreed better with experiment. The difference between calculations using dielectric constants of 4 and 20 was not very apparent except for the case of a highly charged binding pocket in one receptor (Wang C. H. et al., 2016). Aside from the study of higher solute dielectric constants, a residue-dependent dielectric model was also developed for use in an alanine scanning protocol with the MMPBSA method (Simoes et al., 2017). An attempt to modify the solute dielectric environment by incorporating structurally important, explicit water molecules in protein-ligand pockets for MMPBSA calculations was also reported, and it was found to improve the modeling of binding affinities for a series of JNK3 kinase inhibitors (Zhu Y. L. et al., 2014).
A hybrid QM/MM solute was also used in MMPBSA applications for predicting the binding affinities of FabI inhibitors (Su et al., 2015). The study suggested that the prediction results are sensitive to radii sets, GB methods, QM Hamiltonians, sampling protocols, and simulation length. The finding here appears to contradict our prior discussion. This is because the solute dielectric constant is set to 1 (vaccum) due to the explicit consideration of polarization in QM/MM studies. In the study of (Wang C. H. et al., 2016), a high solute dielectric constant was used to mimic polarization implicitly. In general, the use of high dielectric constants washes out the sensitivity to radii and surface definitions.
The non-polar solvation term
The non-polar (non-electrostatic) solvation free energy contribution (ΔGnon-polar) to the solvation free energy change (ΔGbind, solv) arises from the solute cavity formation within the solvent and van der Waals interactions between the solute and the solvent around the cavity (Weeks et al., 1971; Smith and Tanford, 1973; Pratt and Chandler, 1977, 1980; Widom, 1982; Kang et al., 1987; Floris and Tomasi, 1989; Floris et al., 1991; Ashbaugh et al., 1999; Hummer, 1999; Lum et al., 1999; Gallicchio et al., 2000, 2002; Levy et al., 2003; Zacharias, 2003; Gallicchio and Levy, 2004; Su and Gallicchio, 2004; Dzubiella et al., 2006; Wagoner and Baker, 2006; Tan et al., 2007). Until recent times, the non-polar solvation free energy has been simply estimated to be proportional to the solvent accessible surface area (SASA) of the solute:
This is referred to as the classical approach below. The surface tension γ and the correction term b are usually set to be constant for all solute molecules; for example, these are 0.00542 kcal/mol-Å2 and 0.92 kcal/mol, respectively, in the AMBER package (Case et al., 2016).
In more modern approaches (Tan et al., 2007), the cavity formation free energy and van der Waals (dispersion) free energy are modeled as separate terms because they scale differently vs. solute size. One way to correlate the cavity formation free energy is to use the volume (SAV) enclosed by the solvent accessible surface. A solvent accessible volume integration, or a solvent accessible surface integration, can be used to compute the dispersion term (ΔGdispersion), so the total non-polar solvation free energy can be estimated as:
These scaling factors apparently depend on the choices of atomic and solvent probe radii, for example, they are set to 0.0378 kcal/mol-Å2 and −0.569 kcal/mol in the AMBER package (Case et al., 2016). It is interesting to note that the coefficient is much higher compared to the classical model since the dispersion term in solvation is always a negative term. The overall non-polar solvation free energies are similar between the classical and modern models, at least for small molecules where all atoms are exposed. The performance of both the classical and modern non-polar solvation models was analyzed, and it was found that the modern approach reduced the root-mean-square deviations of computed binding affinities (both relative and absolute) from experimental values while the correlations changed little from those computed using the classical approach (Wang C. H. et al., 2016).
The entropy term
The configurational entropy (S) in Equation (1) is often approximated by normal mode or quasi-harmonic analysis. Many proposed methods exist to calculate entropy (Kassem et al., 2015), but it is notoriously difficult to obtain a converged quantity. Thus, further approximations are often utilized; for example, usually only residues within a small sphere (radius of 8–12 Å) centered at the ligand and a limited number of snapshots (<100) are used for normal mode analysis. Furthermore, the single trajectory approach (i.e., that for the complex) is most often used in MMPBSA studies that do not consider any binding-induced structural changes. In such an approach, configurational entropy computed by the normal mode analysis is often omitted completely in the ranking of relative binding affinities as its inclusion often does not improve the agreement with experiment (Yang et al., 2011).
A few new ideas to approximate entropy were reported recently. For example, a new method termed BEERT (Binding Entropy Estimation of Rotation and Translation) was proposed to approximate configurational entropy changes in terms of the reduction in translational and rotational freedom of the ligand upon protein-ligand binding, starting from the flexible molecule approach (Ben-Shalornit et al., 2017). An interaction entropy (IE) method was also proposed to investigate the entropy change upon binding (Duan L. et al., 2016; Duan et al., 2017). The interaction entropy (IE) contribution to binding free energy was defined as
where is the fluctuation of the protein-ligand interaction energy for both electrostatic and van der Waals interactions. The ensemble average of can be extracted from MD simulations, avoiding normal mode calculations. Both developers of BEERT and IE claim that these methods are highly efficient.
Extension to membrane proteins
The development of implicit membrane models based on existing continuum solvent approaches has advanced, making it possible to extend the MMPBSA method to biological membrane systems. The presence of an implicit membrane adds a complication to the numerical solution of the underlying PBE that is brought about by dielectric inhomogeneities that appear on the boundary surfaces of the computation grid. This issue can be alleviated by employing the periodic boundary condition, which is a common practice for electrostatic computations in MD simulations. The conjugate gradient and successive over-relaxation methods can be adjusted to take into account periodic calculations, but the convergence rate using either method is quite low. This limits their application to MMPBSA calculations which require that a large number of conformations be processed. To improve the convergence rate for use in biomolecular applications, the Incomplete Cholesky preconditioning method and the geometric multigrid method were both extended to incorporate periodicity (Botello-Smith and Luo, 2015). Applications to protein-ligand binding utilizing the newly developed membrane MMPBSA method were also reported (Greene et al., 2016; Xiao et al., 2017).
Extension to the high-speed screening of ligands
The MMPBSA method has also been utilized as a rescoring method for docking applications. In protein-protein docking, it was found that MMPBSA rescoring was more capable at distinguishing correct complex structures from decoys than ZDOCK scoring in a test set of 46 protein-protein complexes for certain combinations of force fields, solvation models, and dielectric constants (Chen et al., 2016). A new binding affinity estimator, PBSA_E, that is based on MMPBSA inputs, was also optimized for protein-ligand binding. The optimization made use of a training set consisting of high-quality experimental data that was gathered from 145 complex structures. When the predicted binding affinities using PBSA_E were compared with the predicted binding affinities using other popular scoring functions such as GlideXP, GlideSP, and SYBYL_F, the PBSA_E method exhibited improved accuracy in terms of both achieving higher correlations with measured binding affinities and lower root-mean-square deviations (Liu X. et al., 2016). A study on MMPBSA in hierarchical virtual screening (HVS) was reported. This study demonstrated the predictability and validity of using the MMPBSA method for lead discovery as it identified novel inhibitors of the p38 MAP kinase by employing a physics-based scoring function combined with a knowledge-based structural filter (Cao et al., 2014).
A study on the performance of MMPBSA in docking applications found that it depends on the choice of the solvent models among many factors that were analyzed. Specifically, the authors found that the choice of solvent models plays a minor role for one-protein-family/one-ligand cases which represent the unbiased protein–ligand complex sampling. However, for the total dataset with biased sampling, where some proteins and their homologs have an overabundant presence in the dataset, they found that numerical PB solvent methods do not perform as well as GB solvent methods. In addition, they showed that numerical PB methods were more sensitive to whether MD simulations were used for averaging. Such methods may be currently more suitable for individual protein binding free energy rankings where MD simulations can be easily conducted. This study also demonstrated that the numerical noise from screening applications that utilize only one or a few structures should be addressed in future developments of the dielectric model for numerical PB methods (Sun et al., 2014b).
Extension to other high-performance analyses
MMPBSA was also adapted for high-performance mutational analysis. Single Amino Acid Mutation based change in Binding free Energy (SAAMBE) was developed to predict changes in binding free energy that are brought about by point mutations. SAAMBE utilized 3D structures of protein-protein complexes in a sequence- and structure-based approach. The method was centered around two components: a MMPBSA-based component, and a set of statistical terms obtained from the physico-chemical properties of protein complexes. A better agreement with experimentally determined binding free energy changes over a set of 1,300 mutations in 43 proteins indicated a significant improvement for predictions made using SAAMBE (Li M. H. et al., 2014; Petukh et al., 2015).
Additionally, MMPBSA was adapted into a novel structure-based multiscale approach to identify the key specificity determining residues (SDRs) of PDZ domains that appeared in explicit solvent MD simulations on PDZ-peptide complexes. SDRs were then used together with knowledge-based scoring functions in a proteome-wide search to locate their interaction partners (Tiwari and Mohanty, 2014).
New toolkits
Over the past few years, several new toolkits were released to facilitate the use of MMPBSA calculations. A free energy workflow tool, FEW, was developed for AMBER to assist in the setup of molecular dynamics simulations in explicit membrane environments. FEW also assists in the setup and execution of effective binding free energy calculations for a single-trajectory implicit solvent/implicit membrane MMPBSA approach that involves multiple ligands binding to the same membrane protein (Homeyer and Gohlke, 2015). g_mmpbsa was developed for GROMACS, and it implemented the MMPBSA approach using subroutines written in-house or sourced from the GROMACS and APBS packages (Kumari et al., 2014). GMXPBSA is another user-friendly suite of Bash/Perl scripts for streamlining MMPBSA calculations for GROMACS users (Paissoni et al., 2014). An easy-to-use pipeline tool named Calculation of Free Energy (CaFE) was published to facilitate both MMPBSA and LIE calculations. CaFe is capable of handling numerous static structure and molecular dynamics trajectory file formats generated by different molecular simulation packages, and it also supports various force field parameters (Liu and Hou, 2016).
Applications of MMPBSA
Protein-ligand binding interactions
Since the last comprehensive review of MMPBSA applications in 2015 (Genheden and Ryde, 2015), MMPBSA has continued to see its role in pharmaceutical research and development increase. MMPBSA was widely applied in the development of anticancer compounds where kinases were identified as the most promising targets. Recent inhibitor design includes the epidermal growth factor receptor kinase domain (Li et al., 2015; Moonrin et al., 2015; Zhao et al., 2017), anaplastic lymphoma kinase (Kong et al., 2015), cyclin-dependent kinases (Li X. L. et al., 2014; Czelen, 2016; Dong et al., 2016a,b, 2017; Arba et al., 2017), extracellular signal-regulated kinase 2 (Chen, 2017), casein kinase 2 (Wang X. W. et al., 2014), sphingosine kinases (Fang et al., 2016), Src/Abl tyrosine kinases (Fong, 2015; Ma et al., 2015), RET tyrosine kinase (Gao et al., 2015), PRK1 (Slynko et al., 2014), Akt kinase (Lu et al., 2015), phosphatidylinositol 3 kinase (Bian et al., 2014), and Myt1 kinase (Wichapong et al., 2014). As an illustration of the method's performance, Slynko et al. observed a correlation coefficient of 0.78 between the experimental pIC50 of 26 PRK1 inhibitors and computational MMPBSA binding affinities. They also utilized a quantitative structure activity relationship (QSAR) model to improve the correlation to 0.88, combining affinities from MMPBSA, QM/MMGBSA, and Glide scoring (Slynko et al., 2014). There were other attractive targets including indoleamine 2,3-dioxygenase 1 (Zou et al., 2017), translationally controlled tumor protein (Kumar R. et al., 2017), estrogen receptor (Anbarasu and Jayanthi, 2017), MutT homolog 1 (Zhou et al., 2016), survivin (Sarvagalla et al., 2016), CD44 (Nguyen et al., 2016), calmodulin (Gonzalez-Andrade et al., 2016), androgen receptor (Liu H. L. et al., 2016), human topoisomerase I (Guruge et al., 2016), Mcl-1 (Zhao et al., 2015), vascular endothelial growth factor receptor-2 (Wu et al., 2015), tubulin (Liao et al., 2014a,b; Santoshi and Naik, 2014; Suri and Naik, 2015; Suri et al., 2015), the Hsp70 protein family (Bhattacharjee et al., 2015; Schneider et al., 2016), the Hsp90 protein family (Arba et al., 2015), glucose 6-phosphate dehydrogenase (Obiol-Pardo et al., 2014; Zhao et al., 2014), lysozyme (Zhan et al., 2015), p53 (Verma S. et al., 2016), wheat germ agglutinin (Parasuraman et al., 2014), bromodomains (Muvva et al., 2014), matrix metalloproteinases (Zhou et al., 2014), protein arginine methyltransferases (Hong et al., 2014; Yan et al., 2014), human arsenic methyltransferase (Abro and Azam, 2016), Atox1 proteins (Wang X. L. et al., 2014), tyrosyl-DNA phosphodiesterase 2 (Kossmann et al., 2016), and urokinase-type plasminogen activator (Sa et al., 2014).
Applications in the development of antibacterial, antiviral, and antiparasitic drugs are also common, including inhibitor designs to targets such as enoyl reductase (Kamsri et al., 2014; Yang et al., 2017), succinate-ubiquinone oxidoreductase (Zhu X. L. et al., 2014; Xiong et al., 2015), thymidylate kinase (Biswas et al., 2017), peptidoglycan recognition proteins (Sahoo et al., 2014a), Mycobacterium tuberculosis pantothenate synthetase (Ntie-Kang et al., 2014), glutamine synthetase (Moreira et al., 2016), FtsZ protein (Zhang H. et al., 2015), and the cytochrome bc1 complex (Zhu et al., 2015). Other efforts focused on HIV viruses, including HIV-1 protease (Li D. et al., 2014; Tzoupis et al., 2014; Chen J. Z. et al., 2015; Meher and Wang, 2015; Chen, 2016; Hu et al., 2016; Sroczynski et al., 2016; Xanthopoulos et al., 2016), gp120 (Wang J. H. et al., 2015), gp41 (Song et al., 2014), reverse transcriptase (Bernardo and Silva, 2014), HIV integrase (Quevedo et al., 2014; Han et al., 2016), and a comparative analysis of inhibitors for HIV-1 and HIV-2 proteases (Chen et al., 2014a). Wright et al. studied 9 inhibitor-bound HIV-1 proteases and compared the absolute binding free energies computed via MMPBSA and MMGBSA with and without normal mode-derived configurational entropy. They found only the values of MMPBSA using the normal mode method were close to experimental values (Wright et al., 2014). Influenza was also an important research topic in antiviral compounds, with targets such as neuraminidase protein (Wang and Chen, 2014; Tran et al., 2015; Chintakrindi et al., 2016; Yang et al., 2016), non-structural proteins (Ai et al., 2014), H7N9 neuraminidase (Phanich et al., 2016), and the M2 proton channel (Homeyer et al., 2016). Other antiviral targets were also examined, such as RNA-dependent RNA polymerase (Yu et al., 2014; Wang J. H. et al., 2016), NS3/4A hepatitis C virus protease (Xue et al., 2014; Fu and Wei, 2015; Meeprasert et al., 2016), human furin (Omotuyi, 2015), and the capsomere of virus-like particles (Li Y. Y. et al., 2014). The method was additionally applied to the identification of novel protease targets for antiviral compounds (Pethe et al., 2017). Studies of antimalarials were also reported. For example, the reduced effectiveness of pyrimethamine, and its relationship to mutation in Plasmodium falciparum dihydrofolate reductase, was studied (Mokmak et al., 2014; Abbat et al., 2015). Falcipain-2, a papain family cysteine protease, was also examined as an antimalarial drug target (Omotuyi, 2014).
A third common group of applications are for studies of targets in neural disorders. Targets such as dopamine D2 receptor (Salmas et al., 2017), monoamine oxidase enzymes (Marsavelski and Vianello, 2017), and opioid G protein-coupled receptors (Leonis et al., 2014) were analyzed as antipsychotic compounds. Here, Leonis et al. employed the recently reported crystal structure of the human κ-opioid receptor (κ-OR) to explain the binding mechanism with its antagonist JDTic and agonist SalA. Both JDTic and SalA show that they are capable of forming a favorable complex in the MMPBSA analysis, which was later confirmed by experiment (Leonis et al., 2014). For Parkinson's disease, ligand-binding to the SUR1 receptor (Santos et al., 2016) and the adenosine receptor (Zhang L. H. et al., 2014) was analyzed. Many targets were studied for Alzheimer's disease including β-secretase (Koukoulitsa et al., 2016), angiotensin-converting enzyme (Bhavaraju et al., 2016), acetylcholinesterase and butyrylcholinesterase (Kurt et al., 2017), and inhibitors directly targeting amyloid aggregation (Berhanu and Masunov, 2015). Ataxin-2 protein was studied for the treatment of spinocerebellar ataxia (Sinha et al., 2017), superoxide dismutase for amyotrophic lateral sclerosis (Zhuang et al., 2016), and Niemann-Pick type Cl and C2 proteins (Poongavanam et al., 2016) that occur in rare neurodegenerative diseases.
MMPBSA was also widely applied in studies of many other major diseases. For blood disorders, targets included human serum albumin (Roy et al., 2015; Kragh-Hansen et al., 2016; Yu et al., 2016a,b), collagen-binding alpha 2 beta 1 integrin (Zhang and Sun, 2014), and human alpha-thrombin (Duan L. L. et al., 2016). In the latter study, Duan et al. investigated the binding affinity between human alpha-thrombin and ligand L86 by employing a nonpolarizable AMBER force field and the polarized protein-specific charge (PPC) force field. They found that the PPC binding affinity was closer to the experimental value reported by Nantermet et al. (2003). For immune disorders, studies on Janus kinases (Zhang W. et al., 2016; Wang J. L. et al., 2017), interleukin 10 cytokine (Ni et al., 2017), and receptor-related orphan receptor-gamma-t (Wang F. F. et al., 2015), were reported. For inflammatory disorders, targets such as COX-2 (Chaudhary and Aparoy, 2017), interleukin 6 (Verma R. et al., 2016), toll-like receptors (Shen et al., 2016), human leukocyte antigen (Kongkaew et al., 2015), chymase enzyme (Verma et al., 2017), tumor necrosis factor (Ivanisenko et al., 2014), and Nalp3 (Sahoo et al., 2014b) were studied. For diabetes, analyses of targets included phosphorylase kinase (Begum et al., 2015), glycogen synthase kinase (Arfeen et al., 2015), protein kinase C beta II (Grewal and Sobhia, 2014), and dipeptidyl peptidase-4 enzyme (Gu et al., 2014; Sneha and Doss, 2016). In a search for compounds for the treatment of chronic kidney disease, Vitamin D receptor and cytochrome P450 were analyzed (Nagamani et al., 2016). Malonyl-CoA decarboxylase (Ling et al., 2016), sirtuins (Karaman and Sippl, 2015), adipocyte fatty-acid binding protein (Chen et al., 2014b), 11β-hydroxysteroid dehydrogenase type 1 (Qian H. Y. et al., 2016), and protein tyrosine phosphatase 1B (Kocakaya, 2014) were reported for treatments of other metabolic diseases. For other diseases and disorders, targets such as pyrroline-5-carboxylate reductase for cutis laxa (Sang et al., 2017), renin complexes for hypertension (Tzoupis et al., 2015), and the type 1 receptor of TGF beta for wound healing (Gesteira et al., 2017) were examined.
There were a wide range of applications to targets outside the scope of pharmaceutical research. Although these studies are not related to drug discovery, accurate modeling of protein-ligand binding is also important in a wide variety of other contexts. For example, Starovoytov et al. calculated MMPBSA binding affinities of BPA-A, BPA-C, and BPA-D bound to human estrogen-related receptor γ in their toxicology research. Their analysis showed that BPA-A was the strongest binder to the receptor (Starovoytov et al., 2014). Other similar studies include globin, for its catalytic mechanism in the hydrolysis of substituted phenyl hexanoates (Ercan et al., 2014), the role of conserved residues in substrate binding to Brassica rapa auxin amidohydrolase (Smolko et al., 2016), the effect of hydrophobic interactions in substrate binding to recombinant enzyme carboxylesterase (Shao et al., 2014), the substrate-enzyme interactions of endo-1,4-β-xylanase (Zhan et al., 2014), azoreductase protein for the biodegradation of azo dyes (Dehghanian et al., 2016; Haghshenas et al., 2016), cytochrome P450 2A6 for nicotine addiction (Lu et al., 2014), Cel48F for producing bioethanol via fiber degradation (Qian M. D. et al., 2016), rubisco for biofuel production (Siqueira et al., 2016), streptavidin-biotin complex in biochemical sensing (Liu F. J. et al., 2016), sorotidine 5-monophosphate decarboxylase for its impressive rate enhancement (Jamshidi et al., 2014), tyrosyl-tRNA synthetases for genetic encoding of unnatural amino acids (Ren et al., 2015), T7 RNA polymerase for generating RNA labels (Borkotoky et al., 2016), AF9 in the YEATS family for the recognition of H3K9ac (Wang Q. et al., 2016), cysteine protease 1 precursor from Zea for the hydrolysate of corn gluten meal (Liu et al., 2014), folate receptor alpha for producing milk with high folate concentration (Sahoo et al., 2014c), and both acid amido synthetase (Wang X. et al., 2015) and brassinosteroid (Lei et al., 2015) for plant growth.
Protein-protein binding interactions
There are many applications of MMPBSA that involve the calculation of protein-protein or protein-peptide binding affinities. For anti-cancer applications, we saw studies on complexes of an ankyrin repeat with integrin-linked kinase binding to PINCH1 (Gautam et al., 2014, 2015), cyclin-dependent kinase 8 with cyclin C (Xu et al., 2014), Hsp90 with Cdc37-derived peptides (Wang L. et al., 2015), and p53 with MDMX (Shi et al., 2015). As an illustration, Wang et al. designed an eleven-residue peptide that was able to bind to Hsp90 with a predicted affinity of 6.9 mM, which was comparable to the experimental value of 3.0 mM (Wang L. et al., 2015).
For anti-viral applications, research efforts were focused on analyzing interactions for ankyrin and domain III of the envelope protein of dengue virus II (Chong et al., 2015; Dubey et al., 2017), a modular capsomere of a murine polyomavirus (MPV) VLP designed to protect against influenza (Zhang L. et al., 2014), antibody recognition of immunoglobulin 2D1 for influenza virus H1N1 (Leong et al., 2015), differential structural dynamics and antigenicity of two HA-specific CTL epitopes binding to HLA-A*0201 for influenza virus H5N1 (Sun and Liu, 2015), a complex of histone deacetylase 6 and ubiquitin-specific protease 5 for influenza virus A (Passos et al., 2016), and for a complex of an antibody and epitope variant of HIV-1 p24 capsid protein (Karim et al., 2015). For antibacterial applications, there were studies of salamander PGRP1 with its splice variant adPGRP1a for the innate immune system (Qi Z. T. et al., 2017), and an adsorption mechanism of human beta-defensin-3 on bacterial membranes (Lee et al., 2016). Based on their MMPBSA analysis, Lee et al. found that the binding affinity of human beta-defensin-3 bound to a gram-positive membrane is over 3 times higher than when it is bound to a gram-negative membrane. Interestingly, this difference is mostly derived from electrostatic interactions, consistent with a net charge that is three times larger for gram-positive membranes compared to gram-negative membranes (Lee et al., 2016).
Studies of protein-protein interactions were also found for other diseases and disorders including serum paraoxonase 1 with high-density lipoprotein for antiatherosclerotic activity (Patra et al., 2014), vascular endothelial growth factor A with binding domains of anti-angiogenic agents for retinal neovascular degenerative diseases (Platania et al., 2015), angiotensin-converting enzyme with Angiotensin-II for hypertension (Guan et al., 2016), titin with T-cap/telethonin for dilated cardiomyopathy (Kumar D. T. et al., 2017), CC chemokine ligand 5 (CCL5) with human neutrophil peptide-1 (HNP1) for chronic inflammatory diseases (Wichapong et al., 2016), and ANKS6 with ANKS3 for polycystic kidney disease (Kan et al., 2016). In the study of CCL5 binding with human HNP1, Wichapong et al. reported a correlation coefficient of r = 0.66 between experiment and MMPBSA results for different species of CCL5. In addition, they confirmed that the entropy change upon binding was negligible in this study, so the entropy term could be ignored when a relative binding free energy was considered (Wichapong et al., 2016).
In addition, there were research efforts that focused on basic mechanisms such as protein stability and conformational dynamics (Bhavaraju and Hansmann, 2015; Getov et al., 2016). Many other studies were also reported, including an analysis of full length amylin oligomer aggregation (Berhanu and Masunov, 2014), effects of single point mutations on amyloid formation (Bhavaraju and Hansmann, 2015; Getov et al., 2016; Petukh et al., 2016), interactions between methylated histone H3 and effector domains of the PHD family in pursuit of a molecular mechanism of epigenetics (Grauffel et al., 2015), the binding mechanism of actin-depolymerizing factor 1 and G-actin (Du et al., 2016), the interaction between phosphotyrosine binding domains and peptides for neuronal development, immune responses, tissue homeostasis, and cell growth (Sain et al., 2016), the cognate transducer complex srII-htrII for the downstream signaling mechanism of sensory rhodopsin (Sahoo and Fujiwara, 2017), the binding of CBP to c-Myb for understanding the exact function of CBP and its interaction with c-Myb (Odoux et al., 2016), the catalytic stability of the tetrameric complex of cystathionine gamma-lyase (El-Sayed et al., 2015), the study of the mechanism of a molecular chaperone using acid-stress chaperone HdeA and its substrate protein (Zhou et al., 2017), and the binding of thiopeptide to a ribosomal subunit in order to understand the structure–activity relationship of thiostrepton derivatives (Wolf et al., 2014).
There are also several non-biomedical applications such as a study on the interactions between two adjacent gamma tubulins within the gamma-tubulin ring complex for growing yeast (Suri et al., 2014), between beta-sheet regions in corneous beta proteins of sauropsids to explain its stability and polymerization into filaments (Calvaresi et al., 2016), and the recognition of Avr protein by eukaryotic transcription factor xa5 of rice to understand the gene-for-gene mechanism that governs the direct interaction of R-Avr protein (Dehury et al., 2015). In the study by Suri et al. computational alanine scanning was employed to determine hotspot contributions in the interaction between gamma tubulins. Their analysis showed that most hotspot mutations reduce affinity by 1.16 kcal/mol, while for very crucial amino acids, such as Asp252 and Arg341, the affinity was decreased by more than 10 kcal/mol, which correlated well with experiment (Suri et al., 2014).
Complexes involving nucleic acids
For DNA-protein interactions, there were mechanistic studies involving proteins such as HU, one of the major nucleoid-associated proteins for stabilizing DNA bending (Kim et al., 2014), highly conserved chromatin protein Cren7 for cellular processes such as transcription, replication, and repair (Chen L. et al., 2015), MutS, which recognizes mismatched DNA in DNA repair using ATP (Ishida and Matsumoto, 2016), copper nucleases for predicting electrostatic interactions with B-DNA (Liu C. M. et al., 2016), and metallopeptides binding to the Drew-Dickerson dodecamer (Galindo-Murillo and Cheatham, 2014). For RNA-protein interactions, the binding of toll-like receptor 3 and 22 was studied in relation to fish viral diseases (Sahoo et al., 2015), and the RNA-recognition motif of RNA-binding protein was analyzed for anti-cancer actitivies (Chang et al., 2016). Two studies on DNA-DNA interactions were also reported, including the calculation of binding affinities of short double stranded oligonucleotides (Yesudas et al., 2015), of three-quartet intramolecular human telomeric DNA G-quadruplexes (Islam et al., 2016), and of strand-strand interactions in a human telomeric tetrameric quadruplex (Chaubey et al., 2015). Yesudas et al. studied oligonucleotides (9–20 mers DNA) with the “3-trajectory” approach, without counter ions, by using GBSA, PBSA, and 3D-RISM-KH methods to calculate binding free energies. They showed that the 3D-RISM-KH method was in better agreement with the experimental data for larger oligonucleotides while the GBSA method performed better for smaller oligonucleotides (Yesudas et al., 2015).
There were also several studies involving small ligands binding to nucleic acids. Interactions of anticancer compounds with DNA, such as cisplatin and oxaliplatin (Jalili et al., 2016), distamycin with the DNA minor groove (Jalili and Maddah, 2017), and plant alkaloid chelerythrine with the human telomere sequence (Ghosh et al., 2015) were reported. Aminobenzimidazole binding to an internal ribosome entry site was studied for its anti-hepatitis C virus effect (Henriksen et al., 2014). The interactions of histone-derived antimicrobial peptides buforin II and DesHDAP1 with DNA were also investigated (Sim et al., 2017). Basic mechanistic studies were carried out, including cationic porphyrin-anthraquinone hybrids binding to DNA duplexes (Arba and Tjahjono, 2015) and G-quadruplexes (Arba et al., 2016), as well as ligands binding to riboswitches upon mutation (Hu et al., 2017). The interactions of reactive metabolites of anticancer compounds with DNA were also analyzed (Tumbi et al., 2014). In the binding of ligands to riboswitches, Hu et al. analyzed the relative binding free energies for a guanine riboswitch (GR) and a GUA complex relative to three complexes: 6GU (3.4 kcal/mol), 2BP (5.48 kcal/mol), and XAN (6.19 kcal/mol). These values were in good agreement with experimental observations of 3.21, 4.12, and 5.47 kcal/mol, respectively (Hu et al., 2017).
Guest-host and nano systems
For guest-host systems, several studies were reported. Two octa acid hosts complexed with six guest molecules were analyzed by Bhakat and Soderhjelm for resolving the problem of trapped water in binding cavities. They performed a well-tempered funnel metadynamics (WT-FM) and MMPBSA analyses for the two octa-acid hosts, OAH (without methyl groups) and OAMe (with methyl groups), using both GAFF and OPLS force fields. Their analyses showed that the binding affinities of WT-FM are basically similar to experimental values while MMPBSA results have errors in the range of 5–10 kcal/mol due to its approximate nature (Bhakat and Soderhjelm, 2017). A binding interaction was analyzed to show that electrostatic interactions have the largest contribution to the stability of the cucurbituril-pseudorotaxane complex (Malhis et al., 2015). Complex stability was analyzed for multiple systems, such as the 1:1 and 1:2 inclusion complexes formed by nor-Seco-cucurbit[10]uril and 1-adamantanmethylammonium in water (El-Barghouthi et al., 2015), cyclodextrin-Ibuprofen complexes (Wang R. M. et al., 2015), E-selectin-oligosaccharide complexes (Barra et al., 2017), and the naringenin-2,6-dimethyl β-cyclodextrin inclusion complex (Sangpheak et al., 2014). Finally, enantiomeric discrimination of chiral organic salts by chiral aza-15-crown-5 ether with C1 symmetry was reported (Kocakaya et al., 2015).
For nano systems, a mechanism explaining how C-60 can block potassium ion channels was proposed (Calvaresi et al., 2015). The authors showed that a new binding site for C-60 exists in the channel cavity at the intracellular entrance of the selectivity filter. The escape barrier from the binding site is ~21 kcal/mol as calculated via the umbrella sampling method, in good agreement with the MMPBSA result. MMPBSA was also used in the theoretical design of the cyclic lipopeptide nanotube as a molecular channel in the lipid bilayer (Izadyar et al., 2015; Khavani et al., 2015, 2017), in the study of the enzyme immobilization mechanism of alpha-chymotrypsin onto carbon nanotubes in organic media (Zhang L. Y. et al., 2015), and the mechanism of carbon nanotube activation of subtilisin Carlsberg in polar and non-polar organic media (Zhang L. Y. et al., 2016).
Monomer stability
For single biomolecule stability research, we saw reports on the structures and energies of the alternate frame folding calbindin-D9k protein (Tong et al., 2015), the effect of copper ions in the stability and structural change of human growth protein (Tazikeh-Lemeski, 2016), and the role of potassium in stabilizing the human telomeric intra-molecular G-quadruplex structure (Wang Z. et al., 2015; Wang and Liu, 2017). In their study of calbindin-D9k, Tong et al. found two transition states and an intermediate state with a first rate-controlling barrier of 4.7 kcal/mol and a second barrier of 1.7 kcal/mol using MMPBSA, both of which are in good agreement with experiment (Tong et al., 2015).
Current limitations and future directions
MMPBSA methods are widely applied to calculate binding affinities at a reasonable computational cost. These computational analyses have provided a large volume of valuable predictive results in a wide variety of studies. Even though MMPBSA is known to be less accurate than some of the more computationally expensive methods, like the free energy perturbation and thermodynamic integration methods, the qualitative agreement is often good enough to aid collaborative efforts involving both computational and experimental researchers. Developers are also actively working to improve MMPBSA methods for higher accuracy and efficiency by introducing better solute and solvent models, by porting the expensive energy computation (mostly involving the solvation terms) to faster GPU platforms, and by improving entropy estimations. Efforts are also needed to extend the MMPBSA method for various screening purposes that involve a large number of ligands and/or mutations to achieve a higher overall level of accuracy and efficiency. It is apparent that applications of the MMPBSA method have grown considerably in many different areas of biomolecular study. Most of these applications involve protein-ligand binding affinity calculations due to their utility in drug discovery research efforts. There are also many applications in the study of biomacromolecular complexes. It is also noteworthy that a few guest-host and nano systems utilized MMPBSA calculations, indicating a wider development space for this method in the future.
Statements
Author contributions
CW: Did the literature research and drafted the review; DG: Extensively revised the review; LX: Wrote the membrane protein extension review; RQ: Wrote the solvation term review; RL: Designed the overall structure of the review and was involved in drafting the review. All people polished the final version.
Acknowledgments
This work was supported by National Institute of Health/NIGMS (GM093040 & GM079383 to RL).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The reviewer WC and handling Editor declared their shared affiliation.
References
1
AbbatS.JainV.BharatamP. V. (2015). Origins of the specificity of inhibitor P218 toward wild-type and mutant PfDHFR: a molecular dynamics analysis. J. Biomol. Struct. Dyn.33, 1913–1928. 10.1080/07391102.2014.979231
2
AbroA.AzamS. S. (2016). Binding free energy based analysis of arsenic (+3 oxidation state) methyltransferase with S-adenosylmethionine. J. Mol. Liq.220, 375–382. 10.1016/j.molliq.2016.04.109
3
AiH. X.ZhangL.ChangA. K.WeiH. Y.CheY. C.LiuH. S. (2014). Virtual screening of potential inhibitors from TCM for the CPSF30 binding site on the NS1A protein of influenza A virus. J. Mol. Model.20:10. 10.1007/s00894-014-2142-7
4
AnbarasuK.JayanthiS. (2017). Designing and optimization of novel human LMTK3 inhibitors against breast cancer - a computational approach. J. Recept. Signal Transduct.37, 51–59. 10.3109/10799893.2016.1155069
5
AqvistJ.MareliusJ. (2001). The linear interaction energy method for predicting ligand binding free energies. Comb. Chem. High Throughput Screen.4, 613–626. 10.2174/1386207013330661
6
ArbaM.TjahjonoD. H. (2015). The binding modes of cationic porphyrin-anthraquinone hybrids to DNA duplexes: in silico study. J. Biomol. Struct. Dyn.33, 657–665. 10.1080/07391102.2014.887480
7
ArbaM.IhsanS.RamadhanL. A. N.TjahjonoD. H. (2017). In silico study of porphyrin-anthraquinone hybrids as CDK2 inhibitor. Comput. Biol. Chem.67, 9–14. 10.1016/j.compbiolchem.2016.12.005
8
ArbaM.KartasasmitaR. E.TjahjonoD. H. (2015). Molecular docking and molecular dynamics simulation of the interaction of cationic imidazolium porphyrin-anthraquinone and Hsp90, in Proceedings of the 3rd International Conference on Computation for Science and Technology, ed TjahjonoD. H. (Paris: Atlantis Press), 1–5. 10.2991/iccst-15.2015.1
9
ArbaM.KartasasmitaR. E.TjahjonoD. H. (2016). Molecular docking and dynamics simulations on the interaction of cationic porphyrin-anthraquinone hybrids with DNA G-quadruplexes. J. Biomol. Struct. Dyn.34, 427–438. 10.1080/07391102.2015.1033015
10
ArfeenM.PatelR.KhanT.BharatamP. V. (2015). Molecular dynamics simulation studies of GSK-3 beta ATP competitive inhibitors: understanding the factors contributing to selectivity. J. Biomol. Struct. Dyn.33, 2578–2593. 10.1080/07391102.2015.1063457
11
AshbaughH. S.KalerE. W.PaulaitisM. E. (1999). A “universal” surface area correlation for molecular hydrophobic phenomena. J. Am. Chem. Soc.121, 9243–9244. 10.1021/ja992119h
12
BarraP. A.RibeiroA. J. M.RamosM. J.JimenezV. A.AldereteJ. B.FernandesP. A. (2017). Binding free energy calculations on E-selectin complexes with sLe(x) oligosaccharide analogs. Chem. Biol. Drug Des.89, 114–123. 10.1111/cbdd.12837
13
BashfordD.KarplusM. (1990). Pkas of ionizable groups in proteins - atomic detail from a continuum electrostatic model. Biochemistry29, 10219–10225. 10.1021/bi00496a010
14
BegumJ.SkamnakiV. T.MoffattC.BischlerN.SarrouJ.SkaltsounisA. L.et al. (2015). An evaluation of indirubin analogues as phosphorylase kinase inhibitors. J. Mol. Graph. Model.61, 231–242. 10.1016/j.jmgm.2015.07.010
15
BennettC. H. (1976). Efficient estimation of free-energy differences from monte-carlo data. J. Comput. Phys.22, 245–268. 10.1016/0021-9991(76)90078-4
16
Ben-ShalornitI. Y.Pfeiffer-MarekS.BaringhausK. H.GohlketH. (2017). Efficient approximation of ligand rotational and translational entropy changes upon binding for use in MM-PBSA calculations. J. Chem. Inf. Model.57, 170–189. 10.1021/acs.jcim.6b00373
17
BerhanuW. M.MasunovA. E. (2014). Full length amylin oligomer aggregation: insights from molecular dynamics simulations and implications for design of aggregation inhibitors. J. Biomol. Struct. Dyn.32, 1651–1669. 10.1080/07391102.2013.832635
18
BerhanuW. M.MasunovA. E. (2015). Atomistic mechanism of polyphenol amyloid aggregation inhibitors: molecular dynamics study of Curcumin, Exifone, and Myricetin interaction with the segment of tau peptide oligomer. J. Biomol. Struct. Dyn.33, 1399–1411. 10.1080/07391102.2014.951689
19
BernardoC. E. P.SilvaP. J. (2014). Computational development of rubromycin-based lead compounds for HIV-1 reverse transcriptase inhibition. PeerJ2:18. 10.7717/peerj.470
20
BhakatS.SoderhjelmP. (2017). Resolving the problem of trapped water in binding cavities: prediction of host-guest binding free energies in the SAMPL5 challenge by funnel metadynamics. J. Comput. Aided Mol. Des.31, 119–132. 10.1007/s10822-016-9948-6
21
BhattacharjeeR.DeviA.MishraS. (2015). Molecular docking and molecular dynamics studies reveal structural basis of inhibition and selectivity of inhibitors EGCG and OSU-03012 toward glucose regulated protein-78 (GRP78) overexpressed in glioblastoma. J. Mol. Model.21:17. 10.1007/s00894-015-2801-3
22
BhavarajuM.HansmannU. H. E. (2015). Effect of single point mutations in a form of systemic amyloidosis. Protein Sci.24, 1451–1462. 10.1002/pro.2730
23
BhavarajuM.PhillipsM.BowmanD.Aceves-HernandezJ. M.HansmannU. H. E. (2016). Binding of ACE-inhibitors to in vitro and patient-derived amyloid-beta fibril models. J. Chem. Phys.144:13. 10.1063/1.4938261
24
BianX. L.DongW. Q.ZhaoY.SunR.KongW. J.LiY. P. (2014). Definition of the binding mode of phosphoinositide 3-kinase alpha-selective inhibitor A-66S through molecular dynamics simulation. J. Mol. Model.20:10. 10.1007/s00894-014-2166-z
25
BiswasA.ShuklaA.VijayanR. S. K.JeyakanthanJ.SekarK. (2017). Crystal structures of an archaeal thymidylate kinase from Sulfolobus tokodaii provide insights into the role of a conserved active site Arginine residue. J. Struct. Biol.197, 236–249. 10.1016/j.jsb.2016.12.001
26
BorkotokyS.MeenaC. K.MuraliA. (2016). Interaction analysis of T7 RNA polymerase with heparin and its low molecular weight derivatives - an in silico approach. Bioinform. Biol. Insights10, 155–166. 10.4137/BBI.S40427
27
Botello-SmithW. M.LuoR. (2015). Applications of MMPBSA to membrane proteins i: efficient numerical solutions of periodic poisson-boltzmann equation. J. Chem. Inf. Model.55, 2187–2199. 10.1021/acs.jcim.5b00341
28
Botello-SmithW. M.LiuX.CaiQ.LiZ.ZhaoH.LuoR. (2012). Numerical poisson-boltzmann model for continuum membrane systems. Chem. Phys. Lett.555, 274–281. 10.1016/j.cplett.2012.10.081
29
CaiQ.HsiehM.-J.WangJ.LuoR. (2010). Performance of nonlinear finite-difference poisson-boltzmann solvers. J. Chem. Theory Comput.6, 203–211. 10.1021/ct900381r
30
CaiQ.WangJ.ZhaoH.-K.LuoR. (2009). On removal of charge singularity in Poisson-Boltzmann equation. J. Chem. Phys.130:145101. 10.1063/1.3099708
31
CaiQ.YeX.WangJ.LuoR. (2011). On-the-fly numerical surface integration for finite-difference poisson-boltzmann methods. J. Chem. Theory Comput.7, 3608–3619. 10.1021/ct200389p
32
CalvaresiM.EckhartL.AlibardiL. (2016). The molecular organization of the beta-sheet region in Corneous beta proteins (beta-keratins) of sauropsids explains its stability and polymerization into filaments. J. Struct. Biol.194, 282–291. 10.1016/j.jsb.2016.03.004
33
CalvaresiM.FuriniS.DomeneC.BottoniA.ZerbettoF. (2015). Blocking the passage: C-60 geometrically clogs K+ channels. ACS Nano9, 4827–4834. 10.1021/nn506164s
34
CaoR.HuangN.WangY. L. (2014). Evaluation and application of MD-PB/SA in structure-based hierarchical virtual screening. J. Chem. Inf. Model.54, 1987–1996. 10.1021/ci5003203
35
CaseD. A.BetzR. M.Botello-SmithW.CeruttiD. S.CheathamT. E.DardenT. A.et al. (2016). “Amber 2016.”San Francisco, CA: University of California.
36
ChangC. E. A.ChenW.GilsonM. K. (2007). Ligand configurational entropy and protein binding. Proc. Natl. Acad. Sci. U.S.A.104, 1534–1539. 10.1073/pnas.0610494104
37
ChangS.ZhangD. W.XuL.WanH.HouT. J.KongR. (2016). Exploring the molecular basis of RNA recognition by the dimeric RNA-binding protein via molecular simulation methods. RNA Biol.13, 1133–1143. 10.1080/15476286.2016.1223007
38
ChaubeyA. K.DubeyK. D.OjhaR. P. (2015). MD simulation of LNA-modified human telomeric G-quadruplexes: a free energy calculation. Med. Chem. Res.24, 753–763. 10.1007/s00044-014-1182-y
39
ChaudharyN.AparoyP. (2017). Deciphering the mechanism behind the varied binding activities of COXIBs through Molecular Dynamic Simulations, MM-PBSA binding energy calculations and per-residue energy decomposition studies. J. Biomol. Struct. Dyn.35, 868–882. 10.1080/07391102.2016.1165736
40
ChenF.LiuH.SunH. Y.PanP. C.LiY. Y.LiD.et al. (2016). Assessing the performance of the MM/PBSA and MM/GBSA methods. 6. Capability to predict protein-protein binding free energies and re-rank binding poses generated by protein-protein docking. Phys. Chem. Chem. Phys.18, 22129–22139. 10.1039/C6CP03670H
41
ChenJ. Z. (2016). Drug resistance mechanisms of three mutations V32I, I47V and V82I in HIV-1 protease toward inhibitors probed by molecular dynamics simulations and binding free energy predictions. RSC Adv.6, 58573–58585. 10.1039/C6RA09201B
42
ChenJ. Z. (2017). Clarifying binding difference of ATP and ADP to extracellular signal-regulated kinase 2 by using molecular dynamics simulations. Chem. Biol. Drug Des.89, 548–558. 10.1111/cbdd.12877
43
ChenJ. Z.LiangZ. Q.WangW.YiC. H.ZhangS. L.ZhangQ. G. (2014a). Revealing origin of decrease in potency of darunavir and amprenavir against HIV-2 relative to HIV-1 protease by molecular dynamics simulations. Sci. Rep.4:6872. 10.1038/srep06872
44
ChenJ. Z.WangJ. N.ZhuW. L. (2014b). Binding modes of three inhibitors 8CA, F8A and I4A to A-FABP studied based on molecular dynamics simulation. PLoS ONE9:e99862. 10.1371/journal.pone.0099862
45
ChenJ. Z.WangX. Y.ZhuT.ZhangQ. G.ZhangJ. Z. H. (2015). A comparative insight into amprenavir resistance of mutations V32I, G48V, I50V, I54V, and I84V in HIV-1 protease based on thermodynamic integration and MM-PBSA methods. J. Chem. Inf. Model.55, 1903–1913. 10.1021/acs.jcim.5b00173
46
ChenL.ZhengQ. C.ZhangH. X. (2015). Insights into the effects of mutations on Cren7-DNA binding using molecular dynamics simulations and free energy calculations. Phys. Chem. Chem. Phys.17, 5704–5711. 10.1039/C4CP05413J
47
ChenW.ChangC. E.GilsonM. K. (2004). Calculation of cyclodextrin binding affinities: energy, entropy, and implications for drug design. Biophys. J.87, 3035–3049. 10.1529/biophysj.104.049494
48
ChintakrindiA. S.MartisE. A.GohilD. J.KothariS. T.ChowdharyA. S.CoutinhoE. C.et al. (2016). A computational model for docking of noncompetitive neuraminidase inhibitors and probing their binding interactions with neuraminidase of influenza virus H5N1. Curr. Comput. Aided Drug Des.12, 272–281. 10.2174/1573409912666160713111242
49
ChongW. L.ZainS. M.Abd RahmanN.OthmanR.OthmanS. B.NimmanpipugP.et al. (2015). Exploration of residue binding energy of potential ankyrin for Dengue virus II from MD simulations, in Proceedings of the 3rd International Conference on Computation for Science and Technology, ed TjahjonoD. H. (Paris: Atlantis Press), 100–103.
50
ColmenaresJ.GaliziaA.OrtizJ.ClematisA.RocchiaW. (2014a). A combined MPI-CUDA parallel solution of linear and nonlinear Poisson-Boltzmann equation. Biomed Res. Int.2014:560987. 10.1155/2014/560987
51
ColmenaresJ.OrtizJ.RocchiaW. (2014b). GPU linear and non-linear Poisson-Boltzmann solver module for DelPhi. Bioinformatics30, 569–570. 10.1093/bioinformatics/btt699
52
CzelenP. (2016). Molecular dynamics study on inhibition mechanism of CDK-2 and GSK-3 beta by CHEMBL272026 molecule. Struct. Chem.27, 1807–1818. 10.1007/s11224-016-0803-0
53
DavisM. E.McCammonJ. A. (1990). Electrostatics in biomolecular structure and dynamics. Chem. Rev.90, 509–521. 10.1021/cr00101a005
54
DehghanianF.HaghshenasH.KayM.TavakolH. (2016). A molecular dynamics investigation on the interaction properties of AzrC and its cofactor. J. Iran. Chem. Soc.13, 2143–2153. 10.1007/s13738-016-0932-9
55
DehuryB.MaharanaJ.SahooB. R.SahuJ.SenP.ModiM. K.et al. (2015). Molecular recognition of avirulence protein (avrxa5) by eukaryotic transcription factor xa5 of rice (Oryza sativa L.): insights from molecular dynamics simulations. J. Mol. Graph. Model.57, 49–61. 10.1016/j.jmgm.2015.01.005
56
DongK. K.WangX.YangX. Y.ZhuX. L. (2016a). Binding mechanism of CDK5 with roscovitine derivatives based on molecular dynamics simulations and MM/PBSA methods. J. Mol. Graph. Model.68, 57–67. 10.1016/j.jmgm.2016.06.007
57
DongK. K.YangX. Y.ZhaoT. T.ZhuX. L. (2016b). Exploring the selectivity of Tetrahydropyrido 1,2-a isoindolone derivatives to GSK3 beta and CDK5 by computational methods. Chin. J. Inorg. Chem.32, 1919–1930. 10.11862/CJIC.2016.263
58
DongK. K.YangX. Y.ZhaoT. T.ZhuX. L. (2017). An insight into the inhibitory selectivity of 4-(Pyrazol-4-yl)-pyrimidines to CDK4 over CDK2. Mol. Simul.43, 599–609. 10.1080/08927022.2017.1279283
59
DuJ.WangX.DongC. H.YangJ. M.YaoX. J. (2016). Computational study of the binding mechanism of actin-depolymerizing factor 1 with actin in Arabidopsis thaliana. PLoS ONE11:e159053. 10.1371/journal.pone.0159053
60
DuanL.LiuX.ZhangJ. Z. H. (2016). Interaction entropy: a new paradigm for highly efficient and reliable computation of protein–ligand binding free energy. J. Am. Chem. Soc.138, 5722–5728. 10.1021/jacs.6b02682
61
DuanL. L.FengG. Q.ZhangQ. G. (2016). Large-scale molecular dynamics simulation: effect of polarization on thrombin-ligand binding energy. Sci. Rep.6:31488. 10.1038/srep31488
62
DuanL. L.FengG. Q.WangX. W.WangL. Z.ZhangQ. G. (2017). Effect of electrostatic polarization and bridging water on CDK2-ligand binding affinities calculated using a highly efficient interaction entropy method. Phys. Chem. Chem. Phys.19, 10140–10152. 10.1039/C7CP00841D
63
DubeyK. D.TiwariG.OjhaR. P. (2017). Targeting domain-III hinging of dengue envelope (DENV-2) protein by MD simulations, docking and free energy calculations. J. Mol. Model.23:102. 10.1007/s00894-017-3259-2
64
DzubiellaJ.SwansonJ. M. J.McCammonJ. A. (2006). Coupling nonpolar and polar solvation free energies in implicit solvent models. J. Chem. Phys.124:084905. 10.1063/1.2171192
65
EdingerS. R.CortisC.ShenkinP. S.FriesnerR. A. (1997). Solvation free energies of peptides: comparison of approximate continuum solvation models with accurate solution of the Poisson-Boltzmann equation. J. Phys. Chem. B101, 1190–1197. 10.1021/jp962156k
66
El-BarghouthiM. I.Abdel-HalimH. M.Haj-IbrahimF. J.BodoorK.AssafK. I. (2015). Molecular dynamics of nor-Seco-cucurbit 10 uril complexes. J. Incl. Phenom. Macrocycl. Chem.82, 323–333. 10.1007/s10847-015-0488-9
67
El-SayedA. S. A.Abdel-AzeimS.IbrahimH. M.YassinM. A.Abdel-GhanyS. E.EsenerS.et al. (2015). Biochemical stability and molecular dynamic characterization of Aspergillus fumigatus cystathionine gamma-lyase in response to various reaction effectors. Enzyme Microb. Technol.81, 31–46. 10.1016/j.enzmictec.2015.08.004
68
ErcanS.ArslanN.KocakayaS. O.PirincciogluN.WilliamsA. (2014). Experimental and theoretical study of the mechanism of hydrolysis of substituted phenyl hexanoates catalysed by globin in the presence of surfactant. J. Mol. Model.20:2096. 10.1007/s00894-014-2096-9
69
FangL.WangX. J.XiM. Y.LiuT. Q.YinD. L. (2016). Assessing the ligand selectivity of sphingosine kinases using molecular dynamics and MM-PBSA binding free energy calculations. Mol. Biosyst.12, 1174–1182. 10.1039/C6MB00067C
70
FisicaroG.GenoveseL.AndreussiO.MarzariN.GoedeckerS. (2016). A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments. J. Chem. Phys.144:014103. 10.1063/1.4939125
71
FlorisF. M.TomasiJ.AhuirJ. L. P. (1991). Dispersion and repulsion contributions to the solvation energy - refinements to a simple computational model in the continuum approximation. J. Comput. Chem.12, 784–791. 10.1002/jcc.540120703
72
FlorisF.TomasiJ. (1989). Evaluation of the dispersion contribution to the solvation energy - a simple computational model in the continuum approximation. J. Comput. Chem.10, 616–627. 10.1002/jcc.540100504
73
FongC. W. (2015). Binding energies of tyrosine kinase inhibitors: error assessment of computational methods for imatinib and nilotinib binding. Comput. Biol. Chem.58, 40–54. 10.1016/j.compbiolchem.2015.05.002
74
FuJ. J.WeiJ. (2015). Molecular dynamics study on drug resistance mechanism of HCV NS3/4A protease inhibitor: BI201335. Mol. Simul.41, 674–682. 10.1080/08927022.2014.917298
75
Galindo-MurilloR.CheathamT. E. (2014). DNA binding dynamics and energetics of cobalt, nickel, and copper metallopeptides. ChemMedChem9, 1252–1259. 10.1002/cmdc.201402020
76
GallicchioE.LevyR. M. (2004). AGBNP: An analytic implicit solvent model suitable for molecular dynamics simulations and high-resolution modeling. J. Comput. Chem.25, 479–499. 10.1002/jcc.10400
77
GallicchioE.KuboM. M.LevyR. M. (2000). Enthalpy-entropy and cavity decomposition of alkane hydration free energies: numerical results and implications for theories of hydrophobic solvation. J. Phys. Chem. B104, 6271–6285. 10.1021/jp0006274
78
GallicchioE.ZhangL. Y.LevyR. M. (2002). The SGB/NP hydration free energy model based on the surface generalized born solvent reaction field and novel nonpolar hydration free energy estimators. J. Comput. Chem.23, 517–529. 10.1002/jcc.10045
79
GaoC. X.GrotliM.ErikssonL. A. (2015). Characterization of interactions and pharmacophore development for DFG-out inhibitors to RET tyrosine kinase. J. Mol. Model.21:167. 10.1007/s00894-015-2708-z
80
GautamV.ChongW. L.WisitponchaiT.NimmanpipugP.ZainS. M.RahmanN. A.et al. (2014). GPU-enabled molecular dynamics simulations of ankyrin kinase complex, in 3rd International Conference on Fundamental and Applied Sciences, eds DassS. C.GuanB. H.BhatA. H.FayeI.SoleimaniH.YahyaN. (Melville, WA: American Institute of Physics), 112–115.
81
GautamV.SabriN. H.ChongW. L.ZainS. M.RahmanN. A.LeeV. S.et al. (2015). Computational alanine scanning mutagenesis: characterizing the hotspots of ILK-ankyrin repeat and PINCH1 complex, in Proceedings of the 3rd International Conference on Computation for Science and Technology, ed TjahjonoD. H. (Paris: Atlantis Press), 92–94.
82
GenhedenS.RydeU. (2015). The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin. Drug Discov.10, 449–461. 10.1517/17460441.2015.1032936
83
GeorgescuR. E.AlexovE. G.GunnerM. R. (2002). Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins. Biophys. J.83, 1731–1748. 10.1016/S0006-3495(02)73940-4
84
GesteiraT. F.Coulson-ThomasV. J.YuanY.ZhangJ. H.NaderH. B.KaoW. W. Y. (2017). Lumican peptides: rational design targeting ALK5/TGFBRI. Sci. Rep.7:42057. 10.1038/srep42057
85
GetovI.PetukhM.AlexovE. (2016). SAAFEC: predicting the effect of single point mutations on protein folding free energy using a knowledge-modified MM/PBSA approach. Int. J. Mol. Sci.17:512. 10.3390/ijms17040512
86
GhoshS.JanaJ.KarR. K.ChatterjeeS.DasguptaD. (2015). Plant alkaloid chelerythrine induced aggregation of human telomere sequence-A unique mode of association between a small molecule and a quadruplex. Biochemistry54, 974–986. 10.1021/bi501117x
87
GilsonM. K. (1995). Theory of electrostatic interactions in macromolecules. Curr. Opin. Struct. Biol.5, 216–223. 10.1016/0959-440X(95)80079-4
88
GohlkeH.CaseD. A. (2004). Converging free energy estimates: MM-PB(GB)SA studies on the protein-protein complex Ras-Raf. J. Comput. Chem.25, 238–250. 10.1002/jcc.10379
89
Gonzalez-AndradeM.Rodriguez-SotresR.Madariaga-MazonA.Rivera-ChavezJ.MataR.Sosa-PeinadoA.et al. (2016). Insights into molecular interactions between CaM and its inhibitors from molecular dynamics simulations and experimental data. J. Biomol. Struct. Dyn.34, 78–91. 10.1080/07391102.2015.1022225
90
GrauffelC.StoteR. H.DejaegereA. (2015). Molecular dynamics for computational proteomics of methylated histone H3. Biochim. Biophys. Acta1850, 1026–1040. 10.1016/j.bbagen.2014.09.015
91
GreeneD.Botello-SmithW. M.FollmerA.XiaoL.LambrosE.LuoR. (2016). Modeling membrane protein-ligand binding interactions: the human purinergic platelet receptor. J. Phys. Chem. B120, 12293–12304. 10.1021/acs.jpcb.6b09535
92
GrewalB. K.SobhiaM. E. (2014). Scaffold hopping for identification of novel PKC beta II inhibitors based on ligand and structural approaches, virtual screening and molecular dynamics study. Comb. Chem. High Throughput Screen.17, 2–11. 10.2174/1386207311301010008
93
GuY. L.WangW.ZhuX. L.DongK. K. (2014). Molecular dynamic simulations reveal the mechanism of binding between xanthine inhibitors and DPP-4. J. Mol. Model.20:2075. 10.1007/s00894-014-2075-1
94
GuanS. S.HanW. W.ZhangH.WangS.ShanY. M. (2016). Insight into the interactive residues between two domains of human somatic Angiotensin-converting enzyme and Angiotensin II by MM-PBSA calculation and steered molecular dynamics simulation. J. Biomol. Struct. Dyn.34, 15–28. 10.1080/07391102.2015.1007167
95
GurugeA. G.UdawatteC.WeerasingheS. (2016). An in silico approach of coumarin-derived inhibitors for human DNA topoisomerase I. Aust. J. Chem.69, 1005–1015. 10.1071/CH16232
96
HaghshenasH.KayM.DehghanianF.TavakolH. (2016). Molecular dynamics study of biodegradation of azo dyes via their interactions with AzrC azoreductase. J. Biomol. Struct. Dyn.34, 453–462. 10.1080/07391102.2015.1039585
97
HanD.SuM.TanJ. J.LiC. H.ZhangX. Y.WangC. X. (2016). Structure-activity relationship and binding mode studies for a series of diketo-acids as HIV integrase inhibitors by 3D-QSAR, molecular docking and molecular dynamics simulations. RSC Adv.6, 27594–27606. 10.1039/C6RA00713A
98
HeadM. S.GivenJ. A.GilsonM. K. (1997). ”Mining minima”: direct computation of conformational free energy. J. Phys. Chem. A101, 1609–1618. 10.1021/jp963817g
99
HenriksenN. M.HayatshahiH. S.DavisD. R.CheathamT. E. (2014). Structural and energetic analysis of 2-Aminobenzimidazole inhibitors in complex with the Hepatitis C virus IRES RNA using molecular dynamics simulations. J. Chem. Inf. Model.54, 1758–1772. 10.1021/ci500132c
100
HomeyerN.GohlkeH. (2015). Extension of the free energy workflow FEW towards implicit solvent/implicit membrane MM-PBSA calculations. Biochim. Biophys. Acta1850, 972–982. 10.1016/j.bbagen.2014.10.013
101
HomeyerN.IoannidisH.KolarovF.GauglitzG.ZikosC.KolocourisA.et al. (2016). Interpreting thermodynamic profiles of aminoadamantane compounds inhibiting the M2 proton channel of Influenza A by free energy calculations. J. Chem. Inf. Model.56, 110–126. 10.1021/acs.jcim.5b00467
102
HongW.LiJ. Y.LaughtonC. A.YapL. F.PatersonI. C.WangH. (2014). Investigating the binding preferences of small molecule inhibitors of human protein arginine methyltransferase 1 using molecular modelling. J. Mol. Graph. Model.51, 193–202. 10.1016/j.jmgm.2014.05.010
103
HonigB.NichollsA. (1995). Classical electrostatics in biology and chemistry. Science268, 1144–1149. 10.1126/science.7761829
104
HsiehM. J.LuoR. (2004). Physical scoring function based on AMBER force field and Poisson-Boltzmann implicit solvent for protein structure prediction. Proteins56, 475–486. 10.1002/prot.20133
105
HsiehM. J.LuoR. (2011). Exploring a coarse-grained distributive strategy for finite-difference Poisson-Boltzmann calculations. J. Mol. Model.17, 1985–1996. 10.1007/s00894-010-0904-4
106
HuG. D.MaA. J.WangJ. H. (2017). Ligand selectivity mechanism and conformational changes in guanine riboswitch by molecular dynamics simulations and free energy calculations. J. Chem. Inf. Model.57, 918–928. 10.1021/acs.jcim.7b00139
107
HuG. D.MaA. J.DouX. H.ZhaoL. L.WangJ. H. (2016). Computational studies of a mechanism for binding and drug resistance in the wild type and four mutations of HIV-1 protease with a GRL-0519 inhibitor. Int. J. Mol. Sci.17:819. 10.3390/ijms17060819
108
HummerG. (1999). Hydrophobic force field as a molecular alternative to surface-area models. J. Am. Chem. Soc.121, 6299–6305. 10.1021/ja984414s
109
IshidaH.MatsumotoA. (2016). Mechanism for verification of mismatched and homoduplex DNAs by nucleotides-bound MutS analyzed by molecular dynamics simulations. Proteins84, 1287–1303. 10.1002/prot.25077
110
IslamB.StadlbauerP.NeidleS.HaiderS.SponerJ. (2016). Can we execute reliable MM-PBSA free energy computations of relative stabilities of different guanine quadruplex folds?J. Phys. Chem. B120, 2899–2912. 10.1021/acs.jpcb.6b01059
111
IvanisenkoN. V.TregubchakT. V.SaikO. V.IvanisenkoV. A.ShchelkunovS. N. (2014). Exploring interaction of TNF and orthopoxviral CrmB protein by surface plasmon resonance and free energy calculation. Protein Pept. Lett.21, 1273–1281.
112
IzadyarM.KhavaniM.HousaindokhtM. R. (2015). A combined molecular dynamic and quantum mechanic study of the solvent and guest molecule effect on the stability and length of heterocyclic peptide nanotubes. Phys. Chem. Chem. Phys.17, 11382–11391. 10.1039/C5CP00973A
113
JaliliS.MaddahM. (2017). Molecular dynamics simulation of the sliding of distamycin anticancer drug along DNA: interactions and sequence selectivity. J. Iran. Chem. Soc.14, 531–540. 10.1007/s13738-016-1001-0
114
JaliliS.MaddahM.SchofieldJ. (2016). Molecular dynamics simulation and free energy analysis of the interaction of platinum-based anti-cancer drugs with DNA. J. Theor. Comput. Chem.15:1650054. 10.1142/S0219633616500541
115
JamshidiS.Rafii-TabarH.JaliliS. (2014). Investigation into mechanism of orotidine 5 '-monophosphate decarboxylase enzyme by MM-PBSA/MM-GBSA and molecular docking. Mol. Simul.40, 469–476. 10.1080/08927022.2013.819579
116
JeancharlesA.NichollsA.SharpK.HonigB.TempczykA.HendricksonT. F.et al. (1991). Electrostatic contributions to solvation energies - comparison of free-energy perturbation and continuum calculations. J. Am. Chem. Soc.113, 1454–1455. 10.1021/ja00004a079
117
KamsriP.KoohatammakunN.SrisupanA.MeewongP.PunkvangA.SaparpakornP.et al. (2014). Rational design of InhA inhibitors in the class of diphenyl ether derivatives as potential anti-tubercular agents using molecular dynamics simulations. SAR QSAR Environ. Res.25, 473–488. 10.1080/1062936X.2014.898690
118
KanW.FangF. Q.ChenL.WangR. G.DengQ. G. (2016). Influence of the R823W mutation on the interaction of the ANKS6-ANKS3: insights from molecular dynamics simulation and free energy analysis. J. Biomol. Struct. Dyn.34, 1113–1122. 10.1080/07391102.2015.1071281
119
KangY. K.NemethyG.ScheragaH. A. (1987). Free-energies of hydration of solute molecules. 1. Improvement of the hydration shell-model by exact computations of overlapping volumes. J. Phys. Chem.91, 4105–4109. 10.1021/j100299a032
120
KaramanB.SipplW. (2015). Docking and binding free energy calculations of sirtuin inhibitors. Eur. J. Med. Chem.93, 584–598. 10.1016/j.ejmech.2015.02.045
121
KarimH. A. A.TayapiwatanaC.NimmanpipugP.ZainS. M.RahmanN. A.LeeV. S. (2015). Molecular dynamics simulation on designed antibodies of HIV-1 capsid protein (p24), in Proceedings of the 3rd International Conference on Computation for Science and Technology, ed TjahjonoD. H. (Paris: Atlantis Press), 85–88.
122
KassemS.AhmedM.El-SheikhS.BarakatK. H. (2015). Entropy in bimolecular simulations: a comprehensive review of atomic fluctuations-based methods. J. Mol. Graph. Model.62, 105–117. 10.1016/j.jmgm.2015.09.010
123
KhavaniM.IzadyarM.HousaindokhtM. R. (2015). Theoretical design of the cyclic lipopeptide nanotube as a molecular channel in the lipid bilayer, molecular dynamics and quantum mechanics approach. Phys. Chem. Chem. Phys.17, 25536–25549. 10.1039/C5CP03136B
124
KhavaniM.IzadyarM.HousaindokhtM. R. (2017). Glucose derivatives substitution and cyclic peptide diameter effects on the stability of the self-assembled cyclic peptide nanotubes; a joint QM/MD study. J. Mol. Graph. Model.71, 28–39. 10.1016/j.jmgm.2016.10.019
125
KimD. H.ImH.JeeJ. G.JangS. B.YoonH. J.KwonA. R.et al. (2014). beta-Arm flexibility of HU from Staphylococcus aureus dictates the DNA-binding and recognition mechanism. Acta Crystallogr. Sec. D Biol. Crystallogr.70, 3273–3289. 10.1107/S1399004714023931
126
KocakayaS. O. (2014). The molecular modeling of novel inhibitors of protein tyrosine phosphatase 1B based on catechol by MD and MM-GB (PB)/SA calculations. Bull. Korean Chem. Soc.35, 1769–1776. 10.5012/bkcs.2014.35.6.1769
127
KocakayaS. O.TurgutY.PirincciogluN. (2015). Enantiomeric discrimination of chiral organic salts by chiral aza-15-crown-5 ether with C-1 symmetry: experimental and theoretical approaches. J. Mol. Model.21:13. 10.1007/s00894-015-2604-6
128
KollmanP. A.MassovaI.ReyesC.KuhnB.HuoS. H.ChongL.et al. (2000). Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. Acc. Chem. Res.33, 889–897. 10.1021/ar000033j
129
KongX. T.PanP. C.LiD.TianS.LiY. Y.HouT. J. (2015). Importance of protein flexibility in ranking inhibitor affinities: modeling the binding mechanisms of piperidine carboxamides as Type I1/2 ALK inhibitors. Phys. Chem. Chem. Phys.17, 6098–6113. 10.1039/C4CP05440G
130
KongkaewS.YotmaneeP.RungrotmongkolT.KaiyawetN.MeeprasertA.KaburakiT.et al. (2015). Molecular dynamics simulation reveals the selective binding of human leukocyte antigen alleles associated with BehCet's disease. PLoS ONE10:e135575. 10.1371/journal.pone.0135575
131
KormanT. P.TanY.-H.WongJ.LuoR.TsaiS.-C. (2008). Inhibition kinetics and emodin cocrystal structure of a type II polyketide ketoreductase. Biochemistry47, 1837–1847. 10.1021/bi7016427
132
KossmannB. R.AbdelmalakM.LopezS.TenderG.YanC. L.PommierY.et al. (2016). Discovery of selective inhibitors of tyrosyl-DNA phosphodiesterase 2 by targeting the enzyme DNA-binding cleft. Bioorg. Med. Chem. Lett.26, 3232–3236. 10.1016/j.bmcl.2016.05.065
133
KoukoulitsaC.Villalonga-BarberC.CsonkaR.AlexiX.LeonisG.DellisD.et al. (2016). Biological and computational evaluation of resveratrol inhibitors against Alzheimer's disease. J. Enzyme Inhib. Med. Chem.31, 67–77. 10.3109/14756366.2014.1003928
134
Kragh-HansenU.MinchiottiL.ColettaA.BienkK.GallianoM.SchiottB.et al. (2016). Mutants and molecular dockings reveal that the primary L-thyroxine binding site in human serum albumin is not the one which can cause familial dysalbuminemic hyperthyroxinemia. Biochim. Biophys. Acta1860, 648–660. 10.1016/j.bbagen.2016.01.001
135
KumarD. T.DossC. G. P.SnehaP.TayubiI. A.SivaR.ChakrabortyC.et al. (2017). Influence of V54M mutation in giant muscle protein titin: a computational screening and molecular dynamics approach. J. Biomol. Struct. Dyn.35, 917–928. 10.1080/07391102.2016.1166456
136
KumarR.MauryaR.SaranS. (2017). Identification of novel inhibitors of the translationally controlled tumor protein (TCTP): insights from molecular dynamics. Mol. Biosyst.13, 510–524. 10.1039/C6MB00850J
137
KumariR.KumarR.LynnA.Open Source Drug DiscoveryC. (2014). g_mmpbsa-A GROMACS tool for high-throughput MM-PBSA calculations. J. Chem. Inf. Model.54, 1951–1962. 10.1021/ci500020m
138
KurtB. Z.GaziogluI.DagA.SalmasR. E.KayikG.DurdagiS.et al. (2017). Synthesis, anticholinesterase activity and molecular modeling study of novel carbamate-substituted thymol/carvacrol derivatives. Bioorg. Med. Chem.25, 1352–1363. 10.1016/j.bmc.2016.12.037
139
LeeJ.JungS. W.ChoA. E. (2016). Molecular insights into the adsorption mechanism of human beta-Defensin-3 on bacterial membranes. Langmuir32, 1782–1790. 10.1021/acs.langmuir.5b04113
140
LeiB. L.LiuJ. Y.YaoX. J. (2015). Unveiling the molecular mechanism of brassinosteroids: insights from structure-based molecular modeling studies. Steroids104, 111–117. 10.1016/j.steroids.2015.09.002
141
LeongP.AmaroR. E.LiW. W. (2015). Molecular dynamics analysis of antibody recognition and escape by human H1N1 influenza hemagglutinin. Biophys. J.108, 2704–2712. 10.1016/j.bpj.2015.04.025
142
LeonisG.AvramopoulosA.SalmasR. E.DurdagiS.YurtseverM.PapadopoulosM. G. (2014). Elucidation of conformational states, dynamics, and mechanism of binding in human kappa-opioid receptor complexes. J. Chem. Inf. Model.54, 2294–2308. 10.1021/ci5002873
143
LevyR. M.ZhangL. Y.GallicchioE.FeltsA. K. (2003). On the nonpolar hydration free energy of proteins: surface area and continuum solvent models for the solute-solvent interaction energy. J. Am. Chem. Soc.125, 9523–9530. 10.1021/ja029833a
144
LiD.ZhangY.ZhaoR. N.FanS.HanJ. G. (2014). Investigation on the mechanism for the binding and drug resistance of wild type and mutations of G86 residue in HIV-1 protease complexed with Darunavir by molecular dynamic simulation and free energy calculation. J. Mol. Model.20:2122. 10.1007/s00894-014-2122-y
145
LiM. H.PetukhM.AlexovE.PanchenkoA. R. (2014). Predicting the impact of missense mutations on protein-protein binding affinity. J. Chem. Theory Comput.10, 1770–1780. 10.1021/ct401022c
146
LiX. L.WangX. W.TianZ. B.ZhaoH. L.LiangD.LiW. S.et al. (2014). Structural basis of valmerins as dual inhibitors of GSK3 beta/CDK5. J. Mol. Model.20:2407. 10.1007/s00894-014-2407-1
147
LiY. Y.LiuX. D.DongX. Y.ZhangL.SunY. (2014). Biomimetic design of affinity peptide ligand for capsomere of virus-like particle. Langmuir30, 8500–8508. 10.1021/la5017438
148
LiY.LiX.DongZ. G. (2015). Statistical analysis of EGFR structures' performance in virtual screening. J. Comput. Aided Mol. Des.29, 1045–1055. 10.1007/s10822-015-9877-9
149
LiaoS. Y.MoG. Q.ChenJ. C.ZhengK. C. (2014a). Docking and molecular dynamics studies of the binding between Peloruside A and tubulin. J. Enzyme Inhib. Med. Chem.29, 702–709. 10.3109/14756366.2013.845816
150
LiaoS. Y.MoG. Q.ChenJ. C.ZhengK. C. (2014b). Exploration of the binding mode between (-)-zampanolide and tubulin using docking and molecular dynamics simulation. J. Mol. Model.20:2070. 10.1007/s00894-014-2070-6
151
LingB. P.LiuY. X.LiX. P.WangZ. G.BiS. W. (2016). Identification of the active site of human mitochondrial malonyl-coenzyme a decarboxylase: a combined computational study. Proteins84, 792–802. 10.1002/prot.25029
152
LiuC. M.ZhuY. Y.TangM. S. (2016). Theoretical studies on binding modes of copper-based nucleases with DNA. J. Mol. Graph. Model.64, 11–29. 10.1016/j.jmgm.2015.12.003
153
LiuF. J.ZhangJ. Z. H.MeiY. (2016). The origin of the cooperativity in the streptavidin-biotin system: a computational investigation through molecular dynamics simulations. Sci. Rep.6:27190. 10.1038/srep27190
154
LiuH. L.HanR.LiJ. Z.LiuH. X.ZhengL. F. (2016). Molecular mechanism of R-bicalutamide switching from androgen receptor antagonist to agonist induced by amino acid mutations using molecular dynamics simulations and free energy calculation. J. Comput. Aided Mol. Des.30, 1189–1200. 10.1007/s10822-016-9992-2
155
LiuH. M.ChenL. C.LiQ.ZhengM. Z.LiuJ. S. (2014). Computational study on substrate specificity of a novel cysteine Protease 1 precursor from Zea mays. Int. J. Mol. Sci.15, 10459–10478. 10.3390/ijms150610459
156
LiuH.HouT. J. (2016). CaFE: a tool for binding affinity prediction using end-point free energy methods. Bioinformatics32, 2216–2218. 10.1093/bioinformatics/btw215
157
LiuX.LiuJ. F.ZhuT.ZhangL. J.HeX.ZhangJ. Z. H. (2016). PBSA_E: a PBSA-based free energy estimator for protein-ligand binding affinity. J. Chem. Inf. Model.56, 854–861. 10.1021/acs.jcim.6b00001
158
LiuX.WangC.WangJ.LiZ.ZhaoH.LuoR. (2013). Exploring a charge-central strategy in the solution of Poisson's equation for biomolecular applications. Phys. Chem. Chem. Phys. 15, 129–141. 10.1039/C2CP41894K
159
LuH. T.HuangX. Q.AbdulhameedM. D. M.ZhanC. G. (2014). Binding free energies for nicotine analogs inhibiting cytochrome P450 2A6 by a combined use of molecular dynamics simulations and QM/MM-PBSA calculations. Bioorg. Med. Chem.22, 2149–2156. 10.1016/j.bmc.2014.02.037
160
LuJ.ZhangZ. Y.NiZ.ShenH. J.TuZ. G.LiuH. Q.et al. (2015). The non-additive contribution of hydroxyl substituents to Akt kinase-apigenin affinity. Mol. Simul.41, 653–662. 10.1080/08927022.2014.913099
161
LuQ.LuoR. (2003). A Poisson-Boltzmann dynamics method with nonperiodic boundary condition. J. Chem. Phys.119, 11035–11047. 10.1063/1.1622376
162
LumK.ChandlerD.WeeksJ. D. (1999). Hydrophobicity at small and large length scales. J. Phys. Chem. B103, 4570–4577. 10.1021/jp984327m
163
LuoR.GilsonM. K. (2000). Synthetic adenine receptors: direct calculation of binding affinity and entropy. J. Am. Chem. Soc.122, 2934–2937. 10.1021/ja994034m
164
LuoR.DavidL.GilsonM. K. (2002). Accelerated Poisson-Boltzmann calculations for static and dynamic systems. J. Comput. Chem.23, 1244–1253. 10.1002/jcc.10120
165
LuoR.DavidL.HungH.DevaneyJ.GilsonM. K. (1999). Strength of solvent-exposed salt-bridges. J. Phys. Chem. B103, 727–736. 10.1021/jp982715i
166
LuoR.GilsonH. S. R.PotterM. J.GilsonM. K. (2001). The physical basis of nucleic acid base stacking in water. Biophys. J.80, 140–148. 10.1016/S0006-3495(01)76001-8
167
LuoR.HeadM. S.MoultJ.GilsonM. K. (1998). pK(a) shifts in small molecules and HIV protease: electrostatics and conformation. J. Am. Chem. Soc.120, 6138–6146. 10.1021/ja974307i
168
LuoR.MoultJ.GilsonM. K. (1997). Dielectric screening treatment of electrostatic solvation. J. Phys. Chem. B101, 11226–11236. 10.1021/jp9724838
169
LwinT. Z.LuoR. (2005). Overcoming entropic barrier with coupled sampling at dual resolutions. J. Chem. Phys.123:194904. 10.1063/1.2102871
170
LwinT. Z.LuoR. (2006). Force field influences in β-hairpin folding simulations. Protein Sci.15, 2642–2655. 10.1110/ps.062438006
171
LwinT. Z.ZhouR. H.LuoR. (2006). Is Poisson-Boltzmann theory insufficient for protein folding simulations?J. Chem. Phys. 124. 10.1063/1.2161202
172
MaS. J.ZengG. H.FangD. Q.WangJ. P.WuW. J.XieW. G.et al. (2015). Studies of N-9-arenthenyl purines as novel DFG-in and DFG-out dual Src/Abl inhibitors using 3D-QSAR, docking and molecular dynamics simulations. Mol. Biosyst.11, 394–406. 10.1039/C4MB00350K
173
MalhisL. D.BodoorK.AssafK. I.Al-SakhenN. A.El-BarghouthiM. I. (2015). Molecular dynamics simulation of a cucurbituril based molecular switch triggered by pH changes. Comput. Theor. Chem.1066, 104–112. 10.1016/j.comptc.2015.05.010
174
MarsavelskiA.VianelloR. (2017). What a difference a methyl group makes: the selectivity of monoamine oxidase B towards histamine and N-Methylhistamine. Chemistry23, 2915–2925. 10.1002/chem.201605430
175
MarshallS. A.VizcarraC. L.MayoS. L. (2005). One- and two-body decomposable Poisson-Boltzmann methods for protein design calculations. Protein Sci.14, 1293–1304. 10.1110/ps.041259105
176
MeeprasertA.HannongbuaS.KungwanN.RungrotmongkolT. (2016). Effect of D168V mutation in NS3/4A HCV protease on susceptibilities of faldaprevir and danoprevir. Mol. Biosyst.12, 3666–3673. 10.1039/C6MB00610H
177
MeherB. R.WangY. X. (2015). Exploring the drug resistance of V32I and M46L mutant HIV-1 protease to inhibitor TMC114: flap dynamics and binding mechanism. J. Mol. Graph. Model.56, 60–73. 10.1016/j.jmgm.2014.11.003
178
MikulskisP.GenhedenS.RydbergP.SandbergL.OlsenL.RydeU. (2012). Binding affinities in the SAMPL3 trypsin and host-guest blind tests estimated with the MM/PBSA and LIE methods. J. Comput. Aided Mol. Des.26, 527–541. 10.1007/s10822-011-9524-z
179
MillerB. R.McgeeT. D.SwailsJ. M.HomeyerN.GohlkeH.RoitbergA. E. (2012). MMPBSA.py: an efficient program for end-state free energy calculations. J. Chem. Theory Comput.8, 3314–3321. 10.1021/ct300418h
180
MoghaddamS.YangC.RekharskyM.KoY. H.KimK.InoueY.et al. (2011). New ultrahigh affinity host-guest complexes of cucurbit[7]uril with bicyclo[2.2.2]octane and adamantane guests: thermodynamic analysis and evaluation of M2 affinity calculations. J. Am. Chem. Soc.133, 3570–3581. 10.1021/ja109904u
181
MokmakW.ChunsrivirotS.HannongbuaS.YuthavongY.TongsimaS.KamchonwongpaisanS. (2014). Molecular dynamics of interactions between rigid and flexible antifolates and dihydrofolate reductase from pyrimethamine-sensitive and pyrimethamine-resistant Plasmodium falciparum. Chem. Biol. Drug Des.84, 450–461. 10.1111/cbdd.12334
182
MoonrinN.SongtaweeN.RattanabunyongS.ChunsrivirotS.MokmakW.TongsimaS.et al. (2015). Understanding the molecular basis of EGFR kinase domain/MIG-6 peptide recognition complex using computational analyses. BMC Bioinformatics16:103. 10.1186/s12859-015-0528-x
183
MoreiraC.RamosM. J.FernandesP. A. (2016). Glutamine synthetase drugability beyond its active site: exploring oligomerization interfaces and pockets. Molecules21:e1028. 10.3390/molecules21081028
184
MuddanaH. S.GilsonM. K. (2012). Calculation of host-guest binding affinities using a quantum-mechanical energy model. J. Chem. Theory Comput.8, 2023–2033. 10.1021/ct3002738
185
MuddanaH. S.YinJ.SapraN. V.FenleyA. T.GilsonM. K. (2014). Blind prediction of SAMPL4 cucurbit[7]uril binding affinities with the mining minima method. J. Comput. Aided Mol. Des.28, 463–474. 10.1007/s10822-014-9726-2
186
MuvvaC.SingamE. R. A.RamanS. S.SubramanianV. (2014). Structure-based virtual screening of novel, high-affinity BRD4 inhibitors. Mol. Biosyst.10, 2384–2397. 10.1039/C4MB00243A
187
NagamaniS.MuthusamyK.MarshalJ. J. (2016). E-pharmacophore filtering and molecular dynamics simulation studies in the discovery of potent drug-like molecules for chronic kidney disease. J. Biomol. Struct. Dyn.34, 2233–2250. 10.1080/07391102.2015.1111168
188
NantermetP. G.BarrowJ. C.NewtonC. L.PellicoreJ. M.YoungM. B.LewisS. D.et al. (2003). Design and synthesis of potent and selective macrocyclic thrombin inhibitors. Bioorg. Med. Chem. Lett.13, 2781–2784. 10.1016/S0960-894X(03)00506-7
189
NguyenT. T.TranD. P.HuyP. D. Q.HoangZ.CarloniP.PhamP. V.et al. (2016). Ligand binding to anti-cancer target CD44 investigated by molecular simulations. J. Mol. Model.22:165. 10.1007/s00894-016-3029-6
190
NiG. Y.WangY. J.CumminsS.WaltonS.MounseyK.LiuX. S.et al. (2017). Inhibitory mechanism of peptides with a repeating hydrophobic and hydrophilic residue pattern on interleukin-10. Hum. Vaccin. Immunother.13, 518–527. 10.1080/21645515.2016.1238537
191
NichollsA.MobleyD. L.GuthrieJ. P.ChoderaJ. D.BaylyC. I.CooperM. D.et al. (2008). Predicting small-molecule solvation free energies: an informal blind test for computational chemistry. J. Med. Chem.51, 769–779. 10.1021/jm070549+
192
NielsenJ. E.McCammonJ. A. (2003). On the evaluation and optimization of protein X-ray structures for pKa calculations. Protein Sci.12, 313–326. 10.1110/ps.0229903
193
Ntie-KangF.KannanS.WichapongK.OwonoL. C. O.SipplW.MegnassanE. (2014). Binding of pyrazole-based inhibitors to Mycobacterium tuberculosis pantothenate synthetase: docking and MM-GB(PB)SA analysis. Mol. Biosyst.10, 223–239. 10.1039/C3MB70449A
194
Obiol-PardoC.Alcarraz-VizanG.Diaz-MoralliS.CascanteM.Rubio-MartinezJ. (2014). Design of an interface peptide as new inhibitor of human glucose-6-phosphate dehydrogenase. J. Mol. Graph. Model.49, 110–117. 10.1016/j.jmgm.2014.02.004
195
OdouxA.JindalD.TamasT. C.LimB. W. H.PollardD.XuW. (2016). Experimental and molecular dynamics studies showed that CBP KIX mutation affects the stability of CBP:c-Myb complex. Comput. Biol. Chem.62, 47–59. 10.1016/j.compbiolchem.2016.03.004
196
OmotuyiI. O. (2015). Ebola virus envelope glycoprotein derived peptide in human Furin-bound state: computational studies. J. Biomol. Struct. Dyn.33, 461–470. 10.1080/07391102.2014.981207
197
OmotuyiO. I. (2014). Methyl-methoxylpyrrolinone and flavinium nucleus binding signatures on falcipain-2 active site. J. Mol. Model.20:2386. 10.1007/s00894-014-2386-2
198
PaissoniC.SpiliotopoulosD.MuscoG.SpitaleriA. (2014). GMXPBSA 2.0: a GROMACS tool to perform MM/PBSA and computational alanine scanning. Comput. Phys. Commun.185, 2920–2929. 10.1016/j.cpc.2014.06.019
199
ParasuramanP.MuruganV.SelvinJ. F. A.GromihaM. M.FukuiK.VelurajaK. (2014). Insights into the binding specificity of wild type and mutated wheat germ agglutinin towards Neu5Ac alpha(2-3)Gal: a study by in silico mutations and molecular dynamics simulations. J. Mol. Recognit.27, 482–492. 10.1002/jmr.2369
200
PassosC. D.Simoes-PiresC. A.CarruptP. A.NurissoA. (2016). Molecular dynamics of zinc-finger ubiquitin binding domains: a comparative study of histone deacetylase 6 and ubiquitin-specific protease 5. J. Biomol. Struct. Dyn.34, 2581–2598. 10.1080/07391102.2015.1124051
201
PatraM. C.RathS. N.PradhanS. K.MaharanaJ.DeS. (2014). Molecular dynamics simulation of human serum paraoxonase 1 in DPPC bilayer reveals a critical role of transmembrane helix H1 for HDL association. Eur. Biophys. J. Biophys. Lett.43, 35–51. 10.1007/s00249-013-0937-6
202
PerutzM. F. (1978). Electrostatic effects in proteins. Science201, 1187–1191. 10.1126/science.694508
203
PetheM. A.RubensteinA. B.KhareS. D. (2017). Large-scale structure-based prediction and identification of novel protease substrates using computational protein design. J. Mol. Biol.429, 220–236. 10.1016/j.jmb.2016.11.031
204
PetukhM.DaiL. G.AlexovE. (2016). SAAMBE: Webserver to predict the charge of binding free energy caused by amino acids mutations. Int. J. Mol. Sci.17:547. 10.3390/ijms17040547
205
PetukhM.LiM. H.AlexovE. (2015). Predicting binding free energy change caused by point mutations with knowledge-modified MM/PBSA method. PLoS Comput. Biol.11:e4276. 10.1371/journal.pcbi.1004276
206
PhanichJ.RungrotmongkolT.KungwanN.HannongbuaS. (2016). Role of R292K mutation in influenza H7N9 neuraminidase toward oseltamivir susceptibility: MD and MM/PB(GB)SA study. J. Comput. Aided Mol. Des.30, 917–926. 10.1007/s10822-016-9981-5
207
PlataniaC. B. M.Di PaolaL.LeggioG. M.RomanoG. L.DragoF.SalomoneS.et al. (2015). Molecular features of interaction between VEGFA and anti-angiogenic drugs used in retinal diseases: a computational approach. Front. Pharmacol.6:248. 10.3389/fphar.2015.00248
208
PoongavanamV.KongstedJ.WustnerD. (2016). Computational analysis of sterol ligand specificity of the niemann pick C2 protein. Biochemistry55, 5165–5179. 10.1021/acs.biochem.6b00217
209
PrattL. R.ChandlerD. (1977). Theory of hydrophobic effect. J. Chem. Phys.67, 3683–3704. 10.1063/1.435308
210
PrattL. R.ChandlerD. (1980). Effects of solute-solvent attractive forces on hydrophobic correlations. J. Chem. Phys.73, 3434–3441. 10.1063/1.440541
211
QiR.Botello-SmithW. M.LuoR. (2017). Acceleration of linear finite-difference Poisson-Boltzmann methods on graphics processing units. J. Chem. Theory Comput.13, 3378–3387. 10.1021/acs.jctc.7b00336
212
QiZ. T.MengF. C.ZhangQ. H.WangZ. S.QiaoG.XuW.et al. (2017). Structural insights into ligand binding of PGRP1 splice variants in Chinese giant salamander (Andrias davidianus) from molecular dynamics and free energy calculations. J. Mol. Model.23:315. 10.1007/s00894-017-3315-y
213
QianH. Y.ChenJ. J.PanY. L.ChenJ. Z. (2016). Molecular modeling studies of 11 beta-Hydroxysteroid dehydrogenase type 1 inhibitors through receptor-based 3D-QSAR and molecular dynamics simulations. Molecules21:e1222. 10.3390/molecules21091222
214
QianM. D.GuanS. S.ShanY. M.ZhangH.WangS. (2016). Structural and molecular basis of cellulase Cel48F by computational modeling: insight into catalytic and product release mechanism. J. Struct. Biol.194, 347–356. 10.1016/j.jsb.2016.03.012
215
QuevedoM. A.RiboneS. R.BrinonM. C.DehaenW. (2014). Development of a receptor model for efficient in silico screening of HIV-1 integrase inhibitors. J. Mol. Graph. Model.52, 82–90. 10.1016/j.jmgm.2014.06.007
216
RenW.TruongT. M.AiH. W. (2015). Study of the binding energies between unnatural amino acids and engineered orthogonal tyrosyl-tRNA synthetases. Sci. Rep.5:12632. 10.1038/srep12632
217
RoyD.KumarV.JamesJ.ShihabudeenM. S.KulshresthaS.GoelV.et al. (2015). Evidence that chemical Chaperone 4-Phenylbutyric acid binds to human serum albumin at fatty acid binding sites. PLoS ONE10:e133012. 10.1371/journal.pone.0133012
218
SaR. J.FangL.HuangM. D.LiQ. H.WeiY. Q.WuK. C. (2014). Evaluation of interactions between urokinase plasminogen and inhibitors using molecular dynamic simulation and free-energy calculation. J. Phys. Chem. A118, 9113–9119. 10.1021/jp5064319
219
SahooB. R.FujiwaraT. (2017). Conformational states of HAMP domains interacting with sensory rhodopsin membrane systems: an integrated all-atom and coarse-grained molecular dynamics simulation approach. Mol. Biosyst.13, 193–207. 10.1039/C6MB00730A
220
SahooB. R.DikhitM. R.BhoiG. K.MaharanaJ.LenkaS. K.DubeyP. K.et al. (2015). Understanding the distinguishable structural and functional features in zebrafish TLR3 and TLR22, and their binding modes with fish dsRNA viruses: an exploratory structural model analysis. Amino Acids47, 381–400. 10.1007/s00726-014-1872-2
221
SahooB. R.DubeyP. K.GoyalS.BhoiG. K.LenkaS. K.MaharanaJ.et al. (2014a). Exploration of the binding modes of buffalo PGRP1 receptor complexed with meso-diaminopimelic acid and lysine-type peptidoglycans by molecular dynamics simulation and free energy calculation. Chem. Biol. Interact.220, 255–268. 10.1016/j.cbi.2014.06.028
222
SahooB. R.MaharanaJ.BhoiG. K.LenkaS. K.PatraM. C.DikhitM. R.et al. (2014b). A conformational analysis of mouse Nalp3 domain structures by molecular dynamics simulations, and binding site analysis. Mol. Biosyst.10, 1104–1116. 10.1039/C3MB70600A
223
SahooB. R.MaharanaJ.PatraM. C.BhoiG. K.LenkaS. K.DubeyP. K.et al. (2014c). Structural and dynamic investigation of bovine folate receptor alpha (FOLR1), and role of ultra-high temperature processing on conformational and thermodynamic characteristics of FOLR1-folate complex. Colloids Surf. B Biointerfaces121, 307–318. 10.1016/j.colsurfb.2014.05.028
224
SainN.TiwariG.MohantyD. (2016). Understanding the molecular basis of substrate binding specificity of PTB domains. Sci. Rep.6:31418. 10.1038/srep31418
225
SalmasR. E.SeemanP.AksoydanB.SteinM.YurtseverM.DurdagiS. (2017). Biological Insights of the dopaminergic stabilizer ACR16 at the binding pocket of dopamine D2 receptor. ACS Chem. Neurosci.8, 826–836. 10.1021/acschemneuro.6b00396
226
SangP.XieY. H.LiL. H.YeY. J.HuW.WangJ.et al. (2017). Effect of the R119G mutation on human P5CR structure and its interactions with NAD: insights derived from molecular dynamics simulation and free energy analysis. Comput. Biol. Chem.67, 141–149. 10.1016/j.compbiolchem.2016.12.015
227
SangpheakW.KhuntaweeW.WolschannP.PongsawasdiP.RungrotmongkolT. (2014). Enhanced stability of a naringenin/2,6-dimethyl beta-cyclodextrin inclusion complex: molecular dynamics and free energy calculations based on MM- and QM-PBSA/GBSA. J. Mol. Graph. Model.50, 10–15. 10.1016/j.jmgm.2014.03.001
228
SantosG.Giraldez-AlvarezL. D.Avila-RodriguezM.CapaniF.GalembeckE.NetoA. G.et al. (2016). SUR1 receptor interaction with hesperidin and linarin predicts possible mechanisms of action of valeriana officinalis in Parkinson. Front. Aging Neurosci.8:97. 10.3389/fnagi.2016.00097
229
SantoshiS.NaikP. K. (2014). Molecular insight of isotypes specific beta-tubulin interaction of tubulin heterodimer with noscapinoids. J. Comput. Aided Mol. Des.28, 751–763. 10.1007/s10822-014-9756-9
230
SarvagallaS.CheungC. H. A.TsaiJ. Y.HsiehH. P.CoumarM. S. (2016). Disruption of protein-protein interactions: hot spot detection, structure-based virtual screening and in vitro testing for the anti-cancer drug target survivin. RSC Adv.6, 31947–31959. 10.1039/C5RA22927H
231
SchneiderM.RosamM.GlaserM.PatronovA.ShahH.BackK. C.et al. (2016). BiPPred: combined sequence- and structure-based prediction of peptide binding to the Hsp70 chaperone BiP. Proteins84, 1390–1407. 10.1002/prot.25084
232
ShaoH.XuL.YanY. J. (2014). Biochemical characterization of a carboxylesterase from the Archaeon Pyrobaculum sp. 1860 and a rational explanation of its substrate specificity and thermostability. Int. J. Mol. Sci.15, 16885–16910. 10.3390/ijms150916885
233
ShenP. C.LiW. W.WangY.HeX.HeL. Q. (2016). Binding mode of chitin and TLR2 via molecular docking and dynamics simulation. Mol. Simul.42, 936–941. 10.1080/08927022.2015.1124102
234
ShiS. H.ZhangS. L.ZhangQ. G. (2015). Probing difference in binding modes of inhibitors to MDMX by molecular dynamics simulations and different free energy methods. PLoS ONE10:e0141409. 10.1371/journal.pone.0141409
235
ShivakumarD.DengY. Q.RouxB. (2009). Computations of absolute solvation free energies of small molecules using explicit and implicit solvent model. J. Chem. Theory Comput.5, 919–930. 10.1021/ct800445x
236
SimS.WangP.BeyerB. N.CutronaK. J.RadhakrishnanM. L.ElmoreD. E. (2017). Investigating the nucleic acid interactions of histone-derived antimicrobial peptides. FEBS Lett.591, 706–717. 10.1002/1873-3468.12574
237
SimoesI. C. M.CostaI. P. D.CoimbraJ. T. S.RamosM. J.FernandesP. A. (2017). New parameters for higher accuracy in the computation of binding free energy differences upon alanine scanning mutagenesis on protein-protein interfaces. J. Chem. Inf. Model.57, 60–72. 10.1021/acs.jcim.6b00378
238
SinhaS.GoyalS.SomvanshiP.GroverA. (2017). Mechanistic insights into the binding of class IIa HDAC Inhibitors toward spinocerebellar ataxia type-2: A 3D-QSAR and pharmacophore modeling approach. Front. Neurosci.10:606. 10.3389/fnins.2016.00606
239
SiqueiraA. S.LimaA. R. J.Dall'agnolL. T.De AzevedoJ. S. N.VianezJ.GoncalvesE. C. (2016). Comparative modeling and molecular dynamics suggest high carboxylase activity of the Cyanobium sp CACIAM14 RbcL protein. J. Mol. Model.22:68. 10.1007/s00894-016-2943-y
240
SlynkoI.ScharfeM.RumpfT.EibJ.MetzgerE.SchuleR.et al. (2014). Virtual screening of PRK1 inhibitors: ensemble docking, rescoring using binding free energy calculation and QSAR model development. J. Chem. Inf. Model.54, 138–150. 10.1021/ci400628q
241
SmithR.TanfordC. (1973). Hydrophobicity of long-chain alkyl carboxylic-acids, as measured by their distribution between heptane and aqueous-solutions. Proc. Natl. Acad. Sci. U.S.A.70, 289–293. 10.1073/pnas.70.2.289
242
SmolkoA.SupljikaF.MartincicJ.Jajcanin-JozicN.Grabar-BranilovicM.TomicS.et al. (2016). The role of conserved Cys residues in Brassica rapa auxin amidohydrolase: Cys139 is crucial for the enzyme activity and Cys320 regulates enzyme stability. Phys. Chem. Chem. Phys.18, 8890–8900. 10.1039/C5CP06301A
243
SnehaP.DossC. G. P. (2016). Gliptins in managing diabetes - reviewing computational strategy. Life Sci.166, 108–120. 10.1016/j.lfs.2016.10.009
244
SongK. Z.BaoJ.SunY. M.ZhangJ. Z. H. (2014). Binding of N-substituted pyrrole derivatives to HIV-1 gp41. J. Theor. Comput. Chem.13:1450018. 10.1142/S0219633614500187
245
SrinivasanJ.CheathamT. E.CieplakP.KollmanP. A.CaseD. A. (1998). Continuum solvent studies of the stability of DNA, RNA, and phosphoramidate - DNA helices. J. Am. Chem. Soc.120, 9401–9409. 10.1021/ja981844+
246
SroczynskiD.MalinowskiZ.SzczesniakA. K.PakulskaW. (2016). New 1(2H)-phthalazinone derivatives as potent nonpeptidic HIV-1 protease inhibitors: molecular docking studies, molecular dynamics simulation, oral bioavailability and ADME prediction. Mol. Simul.42, 628–641. 10.1080/08927022.2015.1067808
247
StarovoytovO. N.LiuY. L.TanL. X.YangS. Z. (2014). Effects of the hydroxyl group on phenyl based Ligand/ERR gamma protein binding. Chem. Res. Toxicol.27, 1371–1379. 10.1021/tx500082r
248
StraatsmaT. P.McCammonJ. A. (1991). Multiconfiguration thermodynamic integration. J. Chem. Phys.95, 1175–1188. 10.1063/1.461148
249
SuP. C.TsaiC. C.MehboobS.HevenerK. E.JohnsonM. E. (2015). Comparison of radii sets, entropy, QM methods, and sampling on MM-PBSA, MM-GBSA, and QM/MM-GBSA ligand binding energies of F-tularensis enoyl-ACP reductase (FabI). J. Comput. Chem.36, 1859–1873. 10.1002/jcc.24011
250
SuY.GallicchioE. (2004). The non-polar solvent potential of mean force for the dimerization of alanine dipeptide: the role of solute-solvent van der Waals interactions. Biophys. Chem.109, 251–260. 10.1016/j.bpc.2003.11.007
251
SunH. Y.LiY. Y.ShenM. Y.TianS.XuL.PanP. C.et al. (2014a). Assessing the performance of MM/PBSA and MM/GBSA methods. 5. Improved docking performance using high solute dielectric constant MM/GBSA and MM/PBSA rescoring. Phys. Chem. Chem. Phys.16, 22035–22045. 10.1039/C4CP03179B
252
SunH. Y.LiY. Y.TianS.XuL.HouT. J. (2014b). Assessing the performance of MM/PBSA and MM/GBSA methods. 4. Accuracies of MM/PBSA and MM/GBSA methodologies evaluated by various simulation protocols using PDBbind data set. Phys. Chem. Chem. Phys.16, 16719–16729. 10.1039/C4CP01388C
253
SunY. P.LiuQ. (2015). Differential structural dynamics and antigenicity of two similar influenza H5N1 virus HA-specific HLA-A*0201-restricted CLT epitopes. RSC Adv.5, 2318–2327. 10.1039/C4RA08874C
254
SuriC.NaikP. K. (2015). Combined molecular dynamics and continuum solvent approaches (MM-PBSA/GBSA) to predict noscapinoid binding to gamma-tubulin dimer. SAR QSAR Environ. Res.26, 507–519. 10.1080/1062936X.2015.1070200
255
SuriC.HendricksonT. W.JoshiH. C.NaikP. K. (2014). Molecular insight into gamma-gamma tubulin lateral interactions within the gamma-tubulin ring complex (gamma-TuRC). J. Comput. Aided Mol. Des.28, 961–972. 10.1007/s10822-014-9779-2
256
SuriC.JoshiH. C.NaikP. K. (2015). Molecular modeling reveals binding interface of gamma-tubulin with GCP4 and interactions with noscapinoids. Proteins83, 827–843. 10.1002/prot.24773
257
TanC.TanY. H.LuoR. (2007). Implicit nonpolar solvent models. J. Phys. Chem. B111, 12263–12274. 10.1021/jp073399n
258
TanC.YangL.LuoR. (2006). How well does Poisson-Boltzmann implicit solvent agree with explicit solvent? A quantitative analysis. J. Phys. Chem. B110, 18680–18687. 10.1021/jp063479b
259
TanY.LuoR. (2009). Structural and functional implications of p53 missense cancer mutations. BMC Biophys.2:5. 10.1186/1757-5036-2-5
260
TanY.-H.LuoR. (2008). Protein stability prediction: a Poisson-Boltzmann approach. J. Phys. Chem. B112, 1875–1883. 10.1021/jp709660v
261
TangC. L.AlexovE.PyleA. M.HonigB. (2007). Calculation of pK(a)s in RNA: on the structural origins and functional roles of protonated nucleotides. J. Mol. Biol.366, 1475–1496. 10.1016/j.jmb.2006.12.001
262
Tazikeh-LemeskiE. (2016). Binding free energy and the structural changes determination in hGH protein with different concentrations of copper ions (A molecular dynamics simulation study). J. Theor. Comput. Chem.15:1650045. 10.1142/S0219633616500450
263
TiwariG.MohantyD. (2014). Structure-based multiscale approach for identification of interaction partners of PDZ domains. J. Chem. Inf. Model.54, 1143–1156. 10.1021/ci400627y
264
TongM. Q.WangQ.WangY.ChenG. J. (2015). Structures and energies of the transition between two conformations of the alternate frame folding calbindin-D-9k protein: a theoretical study. RSC Adv.5, 65798–65810. 10.1039/C5RA11234F
265
TranN.VanT.NguyenH.LeL. (2015). Identification of Novel Compounds against an R294K Substitution of Influenza A (H7N9) virus using ensemble based drug virtual screening. Int. J. Med. Sci.12, 163–176. 10.7150/ijms.10826
266
TumbiK. M.NandekarP. P.ShaikhN.KesharwaniS. S.SangamwarA. T. (2014). Molecular dynamics simulation studies for DNA sequence recognition by reactive metabolites of anticancer compounds. J. Mol. Recogn.27, 138–150. 10.1002/jmr.2342
267
TzoupisH.LeonisG.AvramopoulosA.MavromoustakosT.PapadopoulosM. G. (2014). Systematic molecular dynamics, MM-PBSA, and Ab initio approaches to the saquinavir resistance mechanism in HIV-1 PR due to 11 double and multiple mutations. J. Phys. Chem. B118, 9538–9552. 10.1021/jp502687q
268
TzoupisH.LeonisG.AvramopoulosA.ReisH.CzyznikowskaZ.ZervaS.et al. (2015). Elucidation of the binding mechanism of renin using a wide array of computational techniques and biological assays. J. Mol. Graph. Model.62, 138–149. 10.1016/j.jmgm.2015.09.015
269
VermaR.YadavM.PradhanD.BhuyanR.AggarwalS.NayekA.et al. (2016). Probing binding mechanism of interleukin-6 and olokizumab: in silico design of potential lead antibodies for autoimmune and inflammatory diseases. J. Recept. Signal Trans.36, 601–616. 10.3109/10799893.2016.1147584
270
VermaS.GroverS.TyagiC.GoyalS.JamalS.SinghA.et al. (2016). Hydrophobic interactions are a key to MDM2 inhibition by polyphenols as revealed by molecular dynamics simulations and MM/PBSA free energy calculations. PLoS ONE11:14. 10.1371/journal.pone.0149014
271
VermaS.SinghA.KumariA.GoyalS.JamalS.SinhaS.et al. (2017). Dissecting the role of mutations in chymase inhibition: free energy and decomposition analysis. Gene609, 68–79. 10.1016/j.gene.2017.01.030
272
WagonerJ. A.BakerN. A. (2006). Assessing implicit models for nonpolar mean solvation forces: the importance of dispersion and volume terms. Proc. Natl. Acad. Sci. U.S.A.103, 8331–8336. 10.1073/pnas.0600118103
273
WangC. H.NguyenP. H.PhamK.HuynhD.LeT. B. N.WangH. L.et al. (2016). Calculating protein-ligand binding affinities with MMPBSA: method and error analysis. J. Comput. Chem.37, 2436–2446. 10.1002/jcc.24467
274
WangC. H.XiaoL.LuoR. (2017). Numerical interpretation of molecular surface field in dielectric modeling of solvation. J. Comput. Chem.38, 1057–1070. 10.1002/jcc.24782
275
WangC.WangJ.CaiQ.LiZ. L.ZhaoH.LuoR. (2013). Exploring high accuracy Poisson-Boltzmann methods for biomolecular simulations. Comput. Theor. Chem.1024, 34–44. 10.1016/j.comptc.2013.09.021
276
WangF. F.YangW.ShiY. H.LeG. W. (2015). 3D-QSAR, molecular docking and molecular dynamics studies of a series of ROR gamma t inhibitors. J. Biomol. Struct. Dyn.33, 1929–1940. 10.1080/07391102.2014.980321
277
WangJ. H.LiY.YangY. F.ZhangJ. X.DuJ.ZhangS. W.et al. (2015). Profiling the interaction mechanism of indole-based derivatives targeting the HIV-1 gp120 receptor. RSC Adv.5, 78278–78298. 10.1039/C5RA04299B
278
WangJ. H.YangY. F.LiY.WangY. H. (2016). Computational study exploring the interaction mechanism of benzimidazole derivatives as potent cattle bovine viral diarrhea virus inhibitors. J. Agric. Food Chem.64, 5941–5950. 10.1021/acs.jafc.6b01067
279
WangJ. L.ChengL. P.WangT. C.DengW.WuF. H. (2017). Molecular modeling study of CP-690550 derivatives as JAK3 kinase inhibitors through combined 3D-QSAR, molecular docking, and dynamics simulation techniques. J. Mol. Graph. Model.72, 178–186. 10.1016/j.jmgm.2016.12.020
280
WangJ.LuoR. (2010). Assessment of linear finite-difference Poisson-Boltzmann solvers. J. Comput. Chem.31, 1689–1698. 10.1002/jcc.21456
281
WangJ.CaiQ.LiZ.-L.ZhaoH.-K.LuoR. (2009). Achieving energy conservation in Poisson-Boltzmann molecular dynamics: accuracy and precision with finite-difference algorithms. Chem. Phys. Lett.468, 112–118. 10.1016/j.cplett.2008.12.049
282
WangJ.CaiQ.XiangY.LuoR. (2012). Reducing Grid Dependence in finite-difference Poisson-Boltzmann calculations. J. Chem. Theor. Comput.8, 2741–2751. 10.1021/ct300341d
283
WangJ.TanC.ChancoE.LuoR. (2010). Quantitative analysis of Poisson-Boltzmann implicit solvent in molecular dynamics. Phys. Chem. Chem. Phys.12, 1194–1202. 10.1039/B917775B
284
WangL.BaoQ. C.XuX. L.JiangF.GuK.JiangZ. Y.et al. (2015). Discovery and identification of Cdc37-derived peptides targeting the Hsp90-Cdc37 protein-protein interaction. RSC Adv.5, 96138–96145. 10.1039/C5RA20408A
285
WangQ.ZhengQ. C.ZhangH. X. (2016). Exploring the mechanism how AF9 recognizes and binds H3K9ac by molecular dynamics simulations and free energy calculations. Biopolymers105, 779–786. 10.1002/bip.22896
286
WangR. M.ZhouH.SiuS. W. I.GanY.WangY. T.OuyangD. F. (2015). Comparison of three molecular simulation approaches for cyclodextrin-ibuprofen complexation. J. Nanomater.8:193049. 10.1155/2015/193049
287
WangX. L.LiC. Q.WangY.ChenG. J. (2014). Interaction of classical platinum agents with the monomeric and dimeric Atox1 proteins: a molecular dynamics simulation study. Int. J. Mol. Sci.15, 75–99. 10.3390/ijms15010075
288
WangX. W.PanP. C.LiY. Y.LiD.HouT. J. (2014). Exploring the prominent performance of CX-4945 derivatives as protein kinase CK2 inhibitors by a combined computational study. Mol. Biosyst.10, 1196–1210. 10.1039/C4MB00013G
289
WangX.DuJ.YaoX. J. (2015). Structural and dynamic basis of acid amido synthetase GH3.1: an investigation of substrate selectivity and major active site access channels. Mol. Biosyst.11, 809–818. 10.1039/C4MB00608A
290
WangY. T.ChenY. C. (2014). Insights from QM/MM modeling the 3D structure of the 2009 H1N1 influenza A virus neuraminidase and its binding interactions with antiviral drugs. Mol. Inform.33, 240–249. 10.1002/minf.201300117
291
WangZ. G.LiuJ. P. (2017). Effects of the central potassium ions on the G-quadruplex and stabilizer binding. J. Mol. Graph. Model.72, 168–177. 10.1016/j.jmgm.2017.01.006
292
WangZ.ChenR.HouL.LiJ.LiuJ. P. (2015). Molecular dynamics and principal components of potassium binding with human telomeric intra-molecular G-quadruplex. Protein Cell6, 423–433. 10.1007/s13238-015-0155-3
293
WarwickerJ. (2004). Improved pK(a) calculations through flexibility based sampling of a water-dominated interaction scheme. Protein Sci.13, 2793–2805. 10.1110/ps.04785604
294
WarwickerJ.WatsonH. C. (1982). Calculation of the electric-potential in the active-site cleft due to alpha-helix dipoles. J. Mol. Biol.157, 671–679. 10.1016/0022-2836(82)90505-8
295
WeeksJ. D.ChandlerD.AndersenH. C. (1971). Role of repulsive forces in determining equilibrium structure of simple liquids. J. Chem. Phys.54:5237. 10.1063/1.1674820
296
WenE. Z.LuoR. (2004). Interplay of secondary structures and side-chain contacts in the denatured state of BBA1. J. Chem. Phys.121, 2412–2421. 10.1063/1.1768151
297
WenE. Z.HsiehM. J.KollmanP. A.LuoR. (2004). Enhanced ab initio protein folding simulations in Poisson-Boltzmann molecular dynamics with self-guiding forces. J. Mol. Graph. Model.22, 415–424. 10.1016/j.jmgm.2003.12.008
298
WichapongK.AlardJ. E.Ortega-GomezA.WeberC.HackengT. M.SoehnleinO.et al. (2016). Structure-based design of peptidic inhibitors of the interaction between CC chemokine ligand 5 (CCL5) and human neutrophil peptides 1 (HNP1). J. Med. Chem.59, 4289–4301. 10.1021/acs.jmedchem.5b01952
299
WichapongK.RoheA.PlatzerC.SlynkoI.ErdmannF.SchmidtM.et al. (2014). Application of docking and QM/MM-GBSA rescoring to screen for novel Myt1 kinase inhibitors. J. Chem. Inf. Model.54, 881–893. 10.1021/ci4007326
300
WidomB. (1982). Potential-distribution theory and the statistical-mechanics of fluids. J. Phys. Chem.86, 869–872. 10.1021/j100395a005
301
WolfA.SchoofS.BaumannS.ArndtH. D.KirschnerK. N. (2014). Structure-activity relationships of thiostrepton derivatives: implications for rational drug design. J. Comput. Aided Mol. Des.28, 1205–1215. 10.1007/s10822-014-9797-0
302
WrightD. W.HallB. A.KenwayO. A.JhaS.CoveneyP. V. (2014). Computing clinically relevant binding free energies of HIV-1 protease inhibitors. J. Chem. Theory Comput.10, 1228–1241. 10.1021/ct4007037
303
WuX. Y.WanS. H.WangG. F.JinH.LiZ. H.TianY. X.et al. (2015). Molecular dynamics simulation and free energy calculation studies of kinase inhibitors binding to active and inactive conformations of VEGFR-2. J. Mol. Graph. Model.56, 103–112. 10.1016/j.jmgm.2014.12.006
304
XanthopoulosD.KritsiE.SupuranC. T.PapadopoulosM. G.LeonisG.ZoumpoulakisP. (2016). Discovery of HIV type1 aspartic protease hit compounds through combined computational approaches. ChemMedChem11, 1646–1652. 10.1002/cmdc.201600220
305
XiaoL.DiaoJ.GreeneD. A.WangJ.LuoR. (2017). A continuum Poisson-Boltzmann model for membrane channel proteins. J. Chem. Theory Comput.13, 3398–3412. 10.1021/acs.jctc.7b00382
306
XieD. X. (2014). New solution decomposition and minimization schemes for Poisson-Boltzmann equation in calculation of biomolecular electrostatics. J. Comput. Phys.275, 294–309. 10.1016/j.jcp.2014.07.012
307
XieD. X.JiangY. (2016). A nonlocal modified Poisson-Boltzmann equation and finite element solver for computing electrostatics of biomolecules. J. Comput. Phys.322, 1–20. 10.1016/j.jcp.2016.06.028
308
XiongL.ZhuX. L.ShenY. Q.WishwaW.LiK.YangG. F. (2015). Discovery of N-benzoxazol-5-yl-pyrazole-4-carboxamides as nanomolar SQR inhibitors. Eur. J. Med. Chem.95, 424–434. 10.1016/j.ejmech.2015.03.060
309
XuW.Amire-BrahimiB.XieX. J.HuangL. Y.JiJ. Y. (2014). All-atomic molecular dynamic studies of human CDK8: insight into the A-loop, point mutations and binding with its partner CycC. Comput. Biol. Chem.51, 1–11. 10.1016/j.compbiolchem.2014.03.003
310
XueW. W.YangY.WangX. T.LiuH. X.YaoX. J. (2014). Computational study on the inhibitor binding mode and allosteric regulation mechanism in hepatitis C virus NS3/4A protein. PLoS ONE9:e87077. 10.1371/journal.pone.0087077
311
YanL. L.YanC. L.QianK.SuH. R.Kofsky-WoffordS. A.LeeW. C.et al. (2014). Diamidine compounds for selective inhibition of protein arginine methyltransferase 1. J. Med. Chem.57, 2611–2622. 10.1021/jm401884z
312
YangT. Y.WuJ. C.YanC. L.WangY. F.LuoR.GonzalesM. B.et al. (2011). Virtual screening using molecular simulations. Proteins Struct. Funct. Bioinform.79, 1940–1951. 10.1002/prot.23018
313
YangX. Y.LuJ. R.YingM.MuJ. B.LiP. C.LiuY. (2017). Docking and molecular dynamics studies on triclosan derivatives binding to FabI. J. Mol. Model.23:13. 10.1007/s00894-016-3192-9
314
YangZ. W.WuF.YuanX. H.ZhangL.ZhangS. L. (2016). Novel binding patterns between ganoderic acids and neuraminidase: insights from docking, molecular dynamics and MM/PBSA studies. J. Mol. Graph. Model.65, 27–34. 10.1016/j.jmgm.2016.02.006
315
YeX.CaiQ.YangW.LuoR. (2009). Roles of boundary conditions in DNA simulations: analysis of ion distributions with the finite-difference Poisson-Boltzmann method. Biophys. J.97, 554–562. 10.1016/j.bpj.2009.05.012
316
YeX.WangJ.LuoR. (2010). A revised density function for molecular surface calculation in continuum solvent models. J. Chem. Theory Comput.6, 1157–1169. 10.1021/ct900318u
317
YesudasJ. P.BlinovN.DewS. K.KovalenkoA. (2015). Calculation of binding free energy of short double stranded oligonucleotides using MM/3D-RISM-KH approach. J. Mol. Liq.201, 68–76. 10.1016/j.molliq.2014.11.017
318
YuH. J.FangY.LuX.LiuY. J.ZhangH. B. (2014). Combined 3D-QSAR, molecular docking, molecular dynamics simulation, and binding free energy calculation studies on the 5-hydroxy-2h-pyridazin-3-one derivatives as HCV NS5B polymerase inhibitors. Chem. Biol. Drug Des.83, 89–105. 10.1111/cbdd.12203
319
YuH. X.LiD. X.GuoB. S.LiuL.PanQ.LiuB. L.et al. (2016a). Effect of Solvent Water Molecules on Human Serum Albumin Complex Docked Paclitaxel by MM-PBSA Method. Singapore: World Scientific Publ Co Pte Ltd.
320
YuH. X.LiD. X.XuF.PanQ.ChaiP.LiuB. L.et al. (2016b). The binding affinity of human serum albumin and paclitaxel through MMPBSA based on docked complex. Mol. Simul.42, 1460–1467. 10.1080/08927022.2016.1198479
321
ZachariasM. (2003). Continuum solvent modeling of nonpolar solvation: improvement by separating surface area dependent cavity and dispersion contributions. J. Phys. Chem. A107, 3000–3004. 10.1021/jp027598c
322
ZhanD. L.YuL.JinH. Y.GuanS. S.HanW. W. (2014). Molecular modeling and MM-PBSA free energy analysis of endo-1,4-beta-xylanase from Ruminococcus albus 8. Int. J. Mol. Sci.15, 17284–17303. 10.3390/ijms151017284
323
ZhanM. Z.GuoM.JiangY. K.WangX. M. (2015). Characterization of the interaction between gallic acid and lysozyme by molecular dynamics simulation and optical spectroscopy. Int. J. Mol. Sci.16, 14786–14807. 10.3390/ijms160714786
324
ZhangH.LuJ.-R.MuJ.-B.LiuJ.-B.YangX.-Y.WangM.-J.et al. (2015). Molecular dynamics simulation and antibacterial mechanism of 3MBA derivatives as FtsZ protein inhibitors. Acta Phys. Chim. Sin.31, 566–575. 10.3866/PKU.WHXB201501061
325
ZhangL.SunY. (2014). Biomimetic design of platelet adhesion inhibitors to block integrin alpha 2 beta 1-collagen interactions: I. Construction of an affinity binding model. Langmuir30, 4725–4733. 10.1021/la404599s
326
ZhangL.TangR. H.BaiS.ConnorsN. K.LuaL. H. L.ChuanY. P.et al. (2014). Energetic changes caused by antigenic module insertion in a virus-like particle revealed by experiment and molecular dynamics simulations. PLoS ONE9:e107313. 10.1371/journal.pone.0107313
327
ZhangL. H.LiuT. J.WangX.WangJ. N.LiG. H.LiY.et al. (2014). Insight into the binding mode and the structural features of the pyrimidine derivatives as human A(2A) adenosine receptor antagonists. Biosystems115, 13–22. 10.1016/j.biosystems.2013.04.003
328
ZhangL. Y.LiY. Z.YuanY.JiangY. Y.GuoY. Z.LiM. L.et al. (2016). Molecular mechanism of carbon nanotube to activate Subtilisin Carlsberg in polar and non-polar organic media. Sci. Rep.6:36838. 10.1038/srep36838
329
ZhangL. Y.XiaoX. C.YuanY.GuoY. Z.LiM. L.PuX. M. (2015). Probing immobilization mechanism of alpha-chymotrypsin onto carbon nanotube in organic media by molecular dynamics simulation. Sci. Rep.5:9297. 10.1038/srep09297
330
ZhangW.LiJ. Z.HuangZ. X.WangH. Y.LuoH.WangX.et al. (2016). Computer-aided identification of potential TYK2 inhibitors from drug database. J. Mol. Struct.1122, 309–317. 10.1016/j.molstruc.2016.05.099
331
ZhaoF. L.YangG. H.XiangS.GaoD. D.ZengC. (2017). In silico analysis of the effect of mutation on epidermal growth factor receptor in non-small-cell lung carcinoma: from mutational analysis to drug designing. J. Biomol. Struct. Dyn.35, 427–434. 10.1080/07391102.2016.1146165
332
ZhaoR. N.FanS.HanJ. G.LiuG. (2015). Molecular dynamics study of segment peptides of Bax, Bim, and Mcl-1 BH3 domain of the apoptosis-regulating proteins bound to the anti-apoptotic Mcl-1 protein. J. Biomol. Struct. Dyn.33, 1067–1081. 10.1080/07391102.2014.929028
333
ZhaoZ. B.LiuY.YaoY. (2014). Computational determination of binding structures and free energies of glucose 6-phosphate dehydrogenase with novel steroid inhibitors. J. Mol. Graph. Model.51, 168–172. 10.1016/j.jmgm.2014.05.009
334
ZhouD.-D.YuY.-Q.WuH.LiY.-N.QiaoJ.-J. (2017). Simulation study on the mechanism of molecular chaperone HdeA and SurA. Prog. Biochem. Biophys.44, 242–252. 10.16476/j.pibb.2016.0347
335
ZhouS. L.WangM.TongZ. F.WangJ. Y. (2016). The recognition mechanism of crizotinib on MTH1: influence of chirality on the bioactivity. Mol. Phys.114, 2364–2372. 10.1080/00268976.2016.1145750
336
ZhouZ. G.YaoQ. Z.LeiD.ZhangQ. Q.ZhangJ. (2014). Investigations on the mechanisms of interactions between matrix metalloproteinase 9 and its flavonoid inhibitors using a combination of molecular docking, hybrid quantum mechanical/molecular mechanical calculations, and molecular dynamics simulations. Can. J. Chem.92, 821–830. 10.1139/cjc-2014-0180
337
ZhuX. L.XiongL.LiH.SongX. Y.LiuJ. J.YangG. F. (2014). Computational and experimental insight into the molecular mechanism of carboxamide inhibitors of succinate-ubquinone oxidoreductase. ChemMedChem9, 1512–1521. 10.1002/cmdc.201300456
338
ZhuX. L.ZhangM. M.LiuJ. J.GeJ. M.YangG. F. (2015). Ametoctradin is a Potent Q(o) Site inhibitor of the mitochondria! respiration complex III. J. Agric. Food Chem.63, 3377–3386. 10.1021/acs.jafc.5b00228
339
ZhuY. L.BerozaP.ArtisD. R. (2014). Including explicit water molecules as part of the protein structure in MM/PBSA calculations. J. Chem. Inf. Model.54, 462–469. 10.1021/ci4001794
340
ZhuangX. Y.ZhaoB.LiuS.SongF. R.CuiF. C.LiuZ. Q.et al. (2016). Noncovalent interactions between superoxide dismutase and flavonoids studied by native mass spectrometry combined with molecular simulations. Anal. Chem.88, 11720–11726. 10.1021/acs.analchem.6b03359
341
ZouY.WangF.WangY.GuoW. J.ZhangY. H.XuQ.et al. (2017). Systematic study of imidazoles inhibiting IDO1 via the integration of molecular mechanics and quantum mechanics calculations. Eur. J. Med. Chem.131, 152–170. 10.1016/j.ejmech.2017.03.021
342
ZwanzigR. W. (1954). High-temperature equation of state by a perturbation method. I. Nonpolar gases. J. Chem. Phys.22, 1420–1426. 10.1063/1.1740409
Summary
Keywords
molecular recognition, binding affinity, free energy simulation, MMPBSA, continuum solvation model
Citation
Wang C, Greene D, Xiao L, Qi R and Luo R (2018) Recent Developments and Applications of the MMPBSA Method. Front. Mol. Biosci. 4:87. doi: 10.3389/fmolb.2017.00087
Received
18 September 2017
Accepted
30 November 2017
Published
10 January 2018
Volume
4 - 2017
Edited by
Chia-en Chang, University of California, Riverside, United States
Reviewed by
Wei Chen, University of California, Riverside, United States; Haifeng Chen, Shanghai Jiao Tong University, China
Updates
Copyright
© 2018 Wang, Greene, Xiao, Qi and Luo.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ray Luo rluo@uci.edu
This article was submitted to Molecular Recognition, a section of the journal Frontiers in Molecular Biosciences
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.