Protein polarization is critical to stabilizing AF-2 and helix-2' domains in ligand binding to PPAR-gamma

Fragment-Based Deep Learning for Simultaneous Prediction of Polarizabilities and NMR Shieldings of Macromolecules and Their Aggregates

Simultaneous prediction of the molecular response properties, such as polarizability and the NMR shielding constant, at a low computational cost is an unresolved issue. We propose to combine a linear-scaling generalized energy-based fragmentation (GEBF) method and deep learning (DL) with both molecular and atomic information-theoretic approach (ITA) quantities as effective descriptors. In GEBF, the total molecular polarizability can be assembled as a linear combination of the corresponding quantities calculated from a set of small embedded subsystems in GEBF. In the new GEBF-DL(ITA) protocol, one can predict subsystem polarizabilities based on the corresponding molecular wave function (thus electron density and ITA quantities) and DL model rather than calculate them from the computationally intensive coupled-perturbed Hartree-Fock or Kohn-Sham equations and finally obtain the total molecular polarizability via a linear combination equation. As a proof-of-concept application, we predict the molecular polarizabilities of large proteins and protein aggregates. GEBF-DL(ITA) is shown to be as accurate enough as GEBF, with mean absolute percentage error <1%. For the largest protein aggregate (>4000 atoms), GEBF-DL(ITA) gains a speedup ratio of 3 compared with GEBF. It is anticipated that when more advanced electronic structure methods are used, this advantage will be more appealing. Moreover, one can also predict the NMR chemical shieldings of proteins with reasonably good accuracy. Overall, the cost-efficient GEBF-DL(ITA) protocol should be a robust theoretical tool for simultaneously predicting polarizabilities and NMR shieldings of large systems.

Advanced molecular mechanisms of modified DRV compounds in targeting HIV-1 protease mutations and interrupting monomer dimerization

HIV-1 protease (PR) plays a crucial role in the treatment of HIV as a key target. The global issue of emerging drug resistance is escalating, and PR mutations pose a substantial challenge to the effectiveness of inhibitors. HIV-1 PR is an ideal model for studying drug resistance to inhibitors. The inhibitor, darunavir (DRV), exhibits a high genetic barrier to viral resistance, but with mutations of residues in the PR, there is also some resistance to DRV. Inhibitors can impede PR in two ways: one involves binding to the active site of the dimerization protease, and the other involves binding to the PR monomer, thereby preventing dimerization. In this study, we aimed to investigate the inhibitory effect of DRV with a modified inhibitor on PR, comparing the differences between wild-type and mutated PR, using molecular dynamics simulations. The inhibitory effect of the inhibitors on PR monomers was subsequently investigated. And molecular mechanics Poisson-Boltzmann surface area evaluated the binding free energy. The energy contribution of individual residues in the complex was accurately calculated by the alanine scanning binding interaction entropy method. The results showed that these inhibitors had strong inhibitory effects against PR mutations, with GRL-142 exhibiting potent inhibition of both the PR monomer and dimer. Improved inhibitors could strengthen hydrogen bonds and interactions with PR, thereby boosting inhibition efficacy. The binding of the inhibitor and mutation of the PR affected the distance between D25 and I50, preventing their dimerization and the development of drug resistance. This study could accelerate research targeting HIV-1 PR inhibitors and help to further facilitate drug design targeting both mechanisms.

The Importance of Charge Transfer and Solvent Screening in the Interactions of Backbones and Functional Groups in Amino Acid Residues and Nucleotides

Quantum mechanical (QM) calculations at the level of density-functional tight-binding are applied to a protein–DNA complex (PDB: 2o8b) consisting of 3763 atoms, averaging 100 snapshots from molecular dynamics simulations. A detailed comparison of QM and force field (Amber) results is presented. It is shown that, when solvent screening is taken into account, the contributions of the backbones are small, and the binding of nucleotides in the double helix is governed by the base–base interactions. On the other hand, the backbones can make a substantial contribution to the binding of amino acid residues to nucleotides and other residues. The effect of charge transfer on the interactions is also analyzed, revealing that the actual charge of nucleotides and amino acid residues can differ by as much as 6 and 8% from the formal integer charge, respectively. The effect of interactions on topological models (protein -residue networks) is elucidated.

Understanding the avidin-biotin binding based on polarized protein-specific charge

In this study, charge updating schemes based on the local polarized protein-specific charge (LPPC) were introduced to vary the atomic charges of the biotin molecule and the residues in close contact during the simulation of the avidin-biotin complexes. The need of the charge variation of the ligand in response to changes in its surroundings was thoroughly studied. The results show that the calculated binding energy difference between biotin (BTN1) and 2'-iminobiotin (BTN2) and avidin is in excellent agreement with the experimental value, thus verifying the feasibility of updating the atomic charges of ligands during the simulation.

A fixed multi-site interaction charge model for an accurate prediction of the QM/MM interactions

A fixed multi-site interaction charge (FMIC) model was proposed for the accurate prediction of intermolecular electrostatic interactions based on the quantum mechanical linear response of a molecule to an external electric field. In such a model, some additional off-center interaction sites were added for capturing multipole interactions for a given molecule. By multivariate least-square fitting analysis of the calculated QM/MM interactions of a given molecule with the electrostatic environment and the electrostatic potentials of the environment at the pre-defined distributed interaction sites, the FMIC of the molecule was obtained. The model system of CO in myoglobin (Mb) was utilized to demonstrate the derivation of the FMIC. The accuracy of FMIC in predicting the electrostatic interactions between CO and the Mb environment was investigated using 10 000 different Mb-CO configurations generated from the 400 ps QM/MM MD simulation. In comparison to the QM/MM calculations at the B3LYP/aug-cc-pVTZ/ff99SB level, the mean unsigned error (MUE) of the results based on the FMIC model was merely 0.10 kcal mol^-1, and the root mean square error (RMSE) was only 0.13 kcal mol^-1, which are significantly lower than the results predicted by the ESP charge model (MUE = 1.45 kcal mol^-1, and RMSE = 1.7 kcal mol^-1, respectively). The transferability of FMIC was tested by applying the obtained FMIC in the wild type Mb-CO system to the mutants of V68F and H64L Mb-CO systems. The MUEs of the obtained results for 10 000 different configurations are both smaller than 0.2 kcal mol^-1 for the V68F and H64L Mb-CO systems in comparison to the B3LYP/aug-cc-pVTZ/ff99SB calculations, and the RMSEs are also lower than 0.2 kcal mol^-1 for both mutants. The applications of FMIC were extended to model the electrostatic interactions between a hydrogen fluoride molecule and 492 waters in a truncated octahedron box; our study showed that the FMIC could give satisfactory results with a MUE of 0.12 kcal mol^-1 and a RMSE of 0.16 kcal mol^-1 in comparison to the B3LYP/aug-cc-pVDZ/TIP3P calculations for 10 000 different configurations generated using the 10 ns classical MD simulation. Therefore, the FMIC method provides an accurate and efficient tool for predicting intermolecular electrostatic interactions, which can be utilized in the future development of molecular force fields.

Role of Peroxisome Proliferator-Activated Receptor Gamma (PPARγ) in Different Disease States: Recent Updates

Peroxisome proliferator-activated receptor (PPAR), a ligand dependant transcription factor, is a member of the nuclear receptor superfamily. PPAR exists in three isoforms i.e. PPAR alpha (PPARα), PPAR beta (PPARβ), and PPAR gamma (PPARγ). These are multi-functional transcription factors and help in regulating inflammation, type 2 diabetes, lipid concentration in the body, metastasis, and tumor growth or angiogenesis. Activation of PPARγ causes inhibition of growth of cultured human breast, gastric, lung, prostate, and other cancer cells. PPARγ is mainly involved in fatty acid storage, glucose metabolism, and homeostasis and adipogenesis regulation. A large number of natural and synthetic ligands bind to PPARγ and modulate its activity. Ligands such as thiazolidinedione, troglitazone, rosiglitazone, pioglitazone effectively bind to PPARγ; however, most of these were found to display severe side effects such as hepatotoxicity, weight gain, cardiovascular complications and bladder tumor. Now the focus is shifted towards the development of dual-acting or pan PPAR ligands. The current review article describes the functions and role of PPARγ in various disease states. In addition, recently reported PPARγ ligands and pan PPAR ligands were discussed in detail. It is envisaged that the present review article may help in the development of potent PPAR ligands with no or minimal side effects.

An electrostatic energy-based charge model for molecular dynamics simulation

The interactions of the polar chemical bonds such as C=O and N-H with an external electric field were investigated, and a linear relationship between the QM/MM interaction energies and the electric field along the chemical bond is established in the range of weak to intermediate electrical fields. The linear relationship indicates that the electrostatic interactions of a polar group with its surroundings can be described by a simple model of a dipole with constant moment under the action of an electric field. This relationship is employed to develop a general approach to generating an electrostatic energy-based charge (EEC) model for molecules containing single or multiple polar chemical bonds. Benchmark test studies of this model were carried out for (CH₃)₂-CO and N-methyl acetamide in explicit water, and the result shows that the EEC model gives more accurate electrostatic energies than those given by the widely used charge model based on fitting to the electrostatic potential (ESP) in direct comparison to the energies computed by the QM/MM method. The MD simulations of the electric field at the active site of ketosteroid isomerase based on EEC demonstrated that EEC gave a better representation of the electrostatic interaction in the hydrogen-bonding environment than the Amber14SB force field by comparison with experiment. The current study suggests that EEC should be better suited for molecular dynamics study of molecular systems with polar chemical bonds such as biomolecules than the widely used ESP or RESP (restrained ESP) charge models.

Sulfur-containing therapeutics in the treatment of Alzheimer's disease

Sulfur is widely existent in natural products and synthetic organic compounds as organosulfur, which are often associated with a multitude of biological activities. OBenzothiazole, in which benzene ring is fused to the 4,5-positions of the thiazolerganosulfur compounds continue to garner increasing amounts of attention in the field of medicinal chemistry, especially in the development of therapeutic agents for Alzheimer's disease (AD). AD is a fatal neurodegenerative disease and the primary cause of age-related dementia posing severe societal and economic burdens. Unfortunately, there is no cure for AD. A lot of research has been conducted on sulfur-containing compounds in the context of AD due to their innate antioxidant potential and some are currently being evaluated in clinical trials. In this review, we have described emerging trends in the field, particularly the concept of multi-targeting and formulation of disease-modifying strategies. SAR, pharmacological targets, in vitro/vivo ADMET, efficacy in AD animal models, and applications in clinical trials of such sulfur compounds have also been discussed. This article provides a comprehensive review of organosulfur-based AD therapeutic agents and provides insights into their future development.

Entropic effect and residue specific entropic contribution to the cooperativity in streptavidin-biotin binding

Molecular dynamics (MD) simulations were performed employing the polarized protein-specific charge (PPC) to explore the origin of the cooperativity in streptavidin-biotin systems (wild type, two single mutations and one double-mutation). The results of the experiment found that the existence of cooperativity is mainly the result of the entropic effect. In this study, the entropic contribution to the binding free energy was calculated using the recently developed interaction entropy (IE) method, and computational results are in excellent agreement with the experimental observations and are further verified by the calculation of the thermodynamic integration. Comparison of different force fields in terms of predicted binding strength ordering, cooperativity of energy and the stability of hydrogen bonding suggests that the PPC force field combined IE method is a suitable choice. In addition, the IE method enables us to obtain the residue-specific entropic contributions to the streptavidin-biotin binding affinity and identify ten hot-spot residues providing the dominant contribution to the cooperative binding. Importantly, the overall cooperativity obtained from the ten residues also comes mainly from the entropic effect in our study. The calculation of the potential of mean force shows that the unbinding of streptavidin-biotin is a multi-step process, and each step corresponds to the formation and rupture of the hydrogen bond network. And S45A mutation may increase the rigidity of the linker region, making the flap region relatively difficult to open. The present study provides significant molecular insight into the binding cooperativity of the streptavidin-biotin complex.

Computational investigation reveals Picrasidine C as selective PPARα lead: binding pattern, selectivity mechanism and ADME/tox profile

Natural products and their derivatives have been recognized as an important source of therapeutic agents for many years. Previously we isolated a dimeric β-carboline-type alkaloid Picrasidine C from the root of Picrasma quassioides as subtype-selective peroxisome proliferator-activated receptor α (PPARα) agonist. In order to modify this natural product for better affinity and druggability, we investigated a series of properties exhibited by Picrasidine C, such as its binding mode with PPARα, the selectivity mechanism over PPARγ, as well as ADME/Tox profile through computational methods including sequence alignment, molecular docking, pharmacophore modeling and molecular dynamics simulations. The detailed information of binding pattern and affinity for Picrasidine C elucidated here will be valuable for chemical modification. Besides, the steric hindrance of residue Phe363 in PPARγ pocket was speculated as the main isoform selectivity mechanism for Picrasidine C, which would be helpful for the design of selective derivatives. ADME/Tox prediction was conducted to avoid potential undesirable pharmacokinetic properties for reducing the risk of failure. Finally, novel skeletons were derived from lead compound by core hopping method, validated through molecular dynamic simulations and MM-GBSA calculation. In short, the information obtained from computational strategy would be valuable for us to find more potent, safe and selective PPARα agonists during structural optimization.

Modification, Antitumor Activity, and Targeted PPARγ Study of 18β-Glycyrrhetinic Acid, an Important Active Ingredient of Licorice

Licorice is a traditional Chinese medicine, which is often used as sweetener and cosmetic ingredients in food and pharmaceutical industries. Among them, glycyrrhetic acid is one of the most important agents. Studies have shown that glycyrrhetic acid exhibited antitumor activities as PPARγ agonist. However, the limited number of PPARγ glycyrrhetinic agonists and their high toxicity greatly limit the design based on the structure. Therefore, clarifying the binding mode between PPARγ and small molecules, we focused on the introduction of a natural active piperazine skeleton in the position of glycyrrhetinic acid C-3. According to the Combination Principle and the Structure-Based Drug Design, 19 glycyrrhetic acid derivatives were designed and synthesized as potential PPARγ agonists. Compounds 4c and 4q were screened as high-efficiency and low-toxicity lead compounds.

Insight Into the Binding Mechanism of p53/pDIQ-MDMX/MDM2 With the Interaction Entropy Method

The study of the p53-MDMX/MDM2 binding sites is a research hotspot for tumor drug design. The inhibition of p53-targeted MDMX/MDM2 has become an effective approach in anti-tumor drug development. In this paper, a theoretically rigorous and computationally accurate method, namely, the interaction entropy (IE) method, combined with the polarized protein-specific charge (PPC) force field, is used to explore the difference in the binding mechanism between p53-MDMX and p53-MDM2. The interaction of a 12mer peptide inhibitor (pDIQ), which is similar to p53 in structure, with MDMX/MDM2 is also studied. The results demonstrate that p53/pDIQ with MDM2 generates a stronger interaction than with MDMX. Compared to p53, pDIQ has larger binding free energies with MDMX and MDM2. According to the calculated binding free energies, the differences in the binding free energy among the four complexes that are obtained from the combination of PPC and IE are more consistent with the experimental values than with the results from the combination of the non-polarizable AMBER force field and IE. In addition, according to the decomposition of the binding free energy, the van der Waals (vdW) interactions are the main driving force for the binding of the four complexes. They are also the main source of the weaker binding affinity of p53/pDIQ-MDMX relative to p53/pDIQ-MDM2. Compared with p53-MDMX/MDM2, according to the analysis of the residue decomposition, the predicated total residue contributions are higher in pDIQ-MDMX/MDM2 than in p53-MDMX/MDM2, which explains why pDIQ has higher binding affinity than p53 with MDMX/MDM2. The current study provides theoretical guidance for understanding the binding mechanisms and designing a potent dual inhibitor that is targeted to MDMX/MDM2.

Drug-resistance mechanisms of three mutations in anaplastic lymphoma kinase against two inhibitors based on MM/PBSA combined with interaction entropy

As a promising drug target in the treatment of lung cancer, anaplastic lymphoma kinase (ALK) and its mutations have been studied widely. This work explored the origin of the resistance mechanism of the ALK mutants again two inhibitors.

Binding of ellagic acid and urolithin metabolites to the CK2 protein, based on the ONIOM method and molecular docking calculations

Using three-layer ONIOM and molecular docking calculations to investigate the binding of urolithins to the active site of the CK2 protein.

Effect of double mutations T790M/L858R on conformation and drug-resistant mechanism of epidermal growth factor receptor explored by molecular dynamics simulations

Epidermal growth factor receptor (EGFR) is one of the most promising targets for the treatment of cancers. Double mutations T790M/L858R lead to different degrees of drug resistance toward inhibitors. In this study, molecular dynamics (MD) simulations followed by principal component analysis are performed to study the conformational changes of EGFR induced by T790M/L858R. The results suggest that T790M/L858R cause obvious disturbance of the structural stability of EGFR. Molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) and residue-based free energy decomposition methods are integrated to explore the drug-resistant mechanism of T790M/L858R toward inhibitors. The results show that the decrease in van der Waals interactions of inhibitors with the mutated EFGR relative to the wild-type (WT) one is the main force inducing drug resistance of T790M/L858R toward inhibitors TAK-285, while drug resistance toward W2P and HKI-272 is dominated by the decrease in van der Waals interactions and the increase in polar interactions. We expect that the information obtained from this study can aid rational design of effective drugs to relieve drug resistance of EGFR induced by T790M/L858R.

The impact of interior dielectric constant and entropic change on HIV-1 complex binding free energy prediction

At present, the calculated binding free energy obtained using the molecular mechanics/Poisson-Boltzmann (Generalized-Born) surface area (MM/PB(GB)SA) method is overestimated due to the lack of knowledge of suitable interior dielectric constants in the simulation on the interaction of Human Immunodeficiency Virus (HIV-1) protease systems with inhibitors. Therefore, the impact of different values of the interior dielectric constant and the entropic contribution when using the MM/PB(GB)SA method to calculate the binding free energy was systemically evaluated. Our results show that the use of higher interior dielectric constants (1.4-2.0) can clearly improve the predictive accuracy of the MM/PBSA and MM/GBSA methods, and computational errors are significantly reduced by including the effects of electronic polarization and using a new highly efficient interaction entropy (IE) method to calculate the entropic contribution. The suitable range for the interior dielectric constant is 1.4-1.6 for the MM/PBSA method; within this range, the correlation coefficient fluctuates around 0.84, and the mean absolute error fluctuates around 2 kcal/mol. Similarly, an interior dielectric constant of 1.8-2.0 produces a correlation coefficient of approximately 0.76 when using the MM/GBSA method. In addition, the entropic contribution of each individual residue was further calculated using the IE method to predict hot-spot residues, and the detailed binding mechanisms underlying the interactions of the HIV-1 protease, its inhibitors, and bridging water molecules were investigated. In this study, the use of a higher interior dielectric constant and the IE method can improve the calculation accuracy of the HIV-1 system.

Exploring the Reasons for Decrease in Binding Affinity of HIV-2 Against HIV-1 Protease Complex Using Interaction Entropy Under Polarized Force Field

In this study, the differences of binding patterns between two type HIV (HIV-1 and HIV-2) protease and two inhibitors (darunavir and amprenavir) are analyzed and compared using the newly developed interaction entropy (IE) method for the entropy change calculation combined with the polarized force field. The functional role of protonation states in the two HIV-2 complexes is investigated and our study finds that the protonated OD1 atom of Asp25' in B chain is the optimal choice. Those calculated binding free energies obtained from the polarized force field combined with IE method are significantly consistent with the experimental observed. The bridging water W301 is favorable to the binding of HIV-1 complexes; however, it is unfavorable to the HIV-2 complexes in current study. The volume of pocket, B-factor of Cα atoms and the distance of flap tip in HIV-2 complexes are smaller than that of HIV-1 consistently. These changes may cause localized rearrangement of residues lining their surface and finally result in the different binding mode for the two types HIV. Predicated hot-spot residues (Ala28/Ala28', Ile50/Ile50', and Ile84/Ile84') are nearly same in the four systems. However, the contribution to the free energy of Asp30 residue is more favorable in HIV-1 system than in HIV-2 system. Current study, to some extent, reveals the origin for the decrease in binding affinity of inhibitors against HIV-2 compared with HIV-1 and will provides theoretical guidance for future design of potent dual inhibitors targeting two type HIV protease.

Toward amino acid typing for proteins in FFLUX

Continuing the development of the FFLUX, a multipolar polarizable force field driven by machine learning, we present a modern approach to atom-typing and building transferable models for predicting atomic properties in proteins. Amino acid atomic charges in a peptide chain respond to the substitution of a neighboring residue and this response can be categorized in a manner similar to atom-typing. Using a machine learning method called kriging, we are able to build predictive models for an atom that is defined, not only by its local environment, but also by its neighboring residues, for a minimal additional computational cost. We found that prediction errors were up to 11 times lower when using a model specific to the correct group of neighboring residues, with a mean prediction of ∼0.0015 au. This finding suggests that atoms in a force field should be defined by more than just their immediate atomic neighbors. When comparing an atom in a single alanine to an analogous atom in a deca-alanine helix, the mean difference in charge is 0.026 au. Meanwhile, the same difference between a trialanine and a deca-alanine helix is only 0.012 au. When compared to deca-alanine models, the transferable models are up to 20 times faster to train, and require significantly less ab initio calculation, providing a practical route to modeling large biological systems. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Effect of polarization on HIV-1protease and fluoro-substituted inhibitors binding energies by large scale molecular dynamics simulations

Molecular dynamics simulations in explicit water are carried out to study the binding of six inhibitors to HIV-1 protease (PR) for up to 700 ns using the standard AMBER force field and polarized protein-specific charge (PPC). PPC is derived from quantum mechanical calculation for protein in solution and therefore it includes electronic polarization effect. Our results show that in all six systems, the bridging water W301 drifts away from the binding pocket in AMBER simulation. However, it is very stable in all six complexes systems using PPC. Especially, intra-protease, protease-inhibitor hydrogen bonds are dynamic stabilized in MD simulation. The computed binding free energies of six complexes have a significantly linear correlation with those experiment values and the correlation coefficient is found to be 0.91 in PPC simulation. However, the result from AMBER simulation shows a weaker correlation with the correlation coefficient of −0.51 due to the lack of polarization effect. Detailed binding interactions of W301, inhibitors with PR are further analyzed and discussed. The present study provides important information to quantitative understanding the interaction mechanism of PR-inhibitor and PR-W301 and these data also emphasizes the importance of both the electronic polarization and the bridging water molecule in predicting precisely binding affinities.

Effect of electrostatic polarization and bridging water on CDK2–ligand binding affinities calculated using a highly efficient interaction entropy method

A new highly efficient interaction entropy (IE) method combined with the polarized protein-specific charge (PPC) force field is employed to investigate the interaction mechanism of CDK2–ligand binding and the effect of the bridging water.

FEREBUS: Highly parallelized engine for kriging training

FFLUX is a novel force field based on quantum topological atoms, combining multipolar electrostatics with IQA intraatomic and interatomic energy terms. The program FEREBUS calculates the hyperparameters of models produced by the machine learning method kriging. Calculation of kriging hyperparameters (θ and p) requires the optimization of the concentrated log-likelihood L̂(θ,p). FEREBUS uses Particle Swarm Optimization (PSO) and Differential Evolution (DE) algorithms to find the maximum of L̂(θ,p). PSO and DE are two heuristic algorithms that each use a set of particles or vectors to explore the space in which L̂(θ,p) is defined, searching for the maximum. The log-likelihood is a computationally expensive function, which needs to be calculated several times during each optimization iteration. The cost scales quickly with the problem dimension and speed becomes critical in model generation. We present the strategy used to parallelize FEREBUS, and the optimization of L̂(θ,p) through PSO and DE. The code is parallelized in two ways. MPI parallelization distributes the particles or vectors among the different processes, whereas the OpenMP implementation takes care of the calculation of L̂(θ,p), which involves the calculation and inversion of a particular matrix, whose size increases quickly with the dimension of the problem. The run time shows a speed-up of 61 times going from single core to 90 cores with a saving, in one case, of ∼98% of the single core time. In fact, the parallelization scheme presented reduces computational time from 2871 s for a single core calculation, to 41 s for 90 cores calculation. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Large-scale molecular dynamics simulation: Effect of polarization on thrombin-ligand binding energy

Molecular dynamics (MD) simulations lasting 500 ns were performed in explicit water to investigate the effect of polarization on the binding of ligands to human α-thrombin based on the standard nonpolarizable AMBER force field and the quantum-derived polarized protein-specific charge (PPC). The PPC includes the electronic polarization effect of the thrombin-ligand complex, which is absent in the standard force field. A detailed analysis and comparison of the results of the MD simulation with experimental data provided strong evidence that intra-protein, protein-ligand hydrogen bonds and the root-mean-square deviation of backbone atoms were significantly stabilized through electronic polarization. Specifically, two critical hydrogen bonds between thrombin and the ligand were broken at approximately 190 ns when AMBER force field was used and the number of intra-protein backbone hydrogen bonds was higher under PPC than under AMBER. The thrombin-ligand binding energy was computed using the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method, and the results were consistent with the experimental value obtained using PPC. Because hydrogen bonds were unstable, it was failed to predict the binding affinity under the AMBER force field. Furthermore, the results of the present study revealed that differences in the binding free energy between AMBER and PPC almost comes from the electrostatic interaction. Thus, this study provides evidence that protein polarization is critical to accurately describe protein-ligand binding.

A new algorithm for construction of coarse-grained sites of large biomolecules

The development of coarse-grained (CG) models for large biomolecules remains a challenge in multiscale simulations, including a rigorous definition of CG representations for them. In this work, we proposed a new stepwise optimization imposed with the boundary-constraint (SOBC) algorithm to construct the CG sites of large biomolecules, based on the s cheme of essential dynamics CG. By means of SOBC, we can rigorously derive the CG representations of biomolecules with less computational cost. The SOBC is particularly efficient for the CG definition of large systems with thousands of residues. The resulted CG sites can be parameterized as a CG model using the normal mode analysis based fluctuation matching method. Through normal mode analysis, the obtained modes of CG model can accurately reflect the functionally related slow motions of biomolecules. The SOBC algorithm can be used for the construction of CG sites of large biomolecules such as F-actin and for the study of mechanical properties of biomaterials.

Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation

Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.

Insights into the Nature of Anesthetic-Protein Interactions: An ONIOM Study

Anesthetics have been employed widely to relieve surgical suffering, but their mechanism of action is not yet clear. For over a century, the mechanism of anesthesia was previously thought to be via lipid bilayer interactions. In the present work, a rigorous three-layer ONIOM(M06-2X/6-31+G*:PM6:AMBER) method was utilized to investigate the nature of interactions between several anesthetics and actual protein binding sites. According to the calculated structural features, interaction energies, atomic charges, and electrostatic potential surfaces, the amphiphilic nature of anesthetic-protein interactions was demonstrated for both inhalational and injectable anesthetics. The existence of hydrogen and halogen bonding interactions between anesthetics and proteins was clearly identified, and these interactions served to assist ligand recognition and binding by the protein. Within all complexes of inhalational or injectable anesthetics, the polarization effects play a dominant role over the steric effects and induce a significant asymmetry in the otherwise symmetric atomic charge distributions of the free ligands in vacuo. This study provides new insight into the mechanism of action of general anesthetics in a more rigorous way than previously described. Future rational design of safer anesthetics for an aging and more physiologically vulnerable population will be predicated on this greater understanding of such specific interactions.

Understanding the basis of I50V-induced affinity decrease in HIV-1 protease via molecular dynamics simulations using polarized force field

Human immunodeficiency virus (HIV)-1 protease is one of the most promising drug target commonly utilized to combat Acquired Immune Deficiency Syndrome (AIDS). However, with the emergence of drug resistance arising from mutations, the efficiency of protease inhibitors (PIs) as a viable treatment for AIDS has been greatly reduced. I50V mutation as one of the most significant mutations occurring in HIV-1 protease will be investigated in this study. Molecular dynamics (MD) simulation was utilized to examine the effect of I50V mutation on the binding of two PIs namely indinavir and amprenavir to HIV-1 protease. Prior to the simulations conducted, the electron density distributions of the PI and each residue in HIV-1 protease are derived by combining quantum fragmentation approach molecular fractionation with conjugate caps and Poisson-Boltzmann solvation model based on polarized protein-specific charge scheme. The atomic charges of the binding complex are subsequently fitted using delta restrained electrostatic potential (delta-RESP) method to overcome the poor charge determination of buried atom. This way, both intraprotease polarization and the polarization between protease and the PI are incorporated into partial atomic charges. Through this study, the mutation-induced affinity variations were calculated and significant agreement between experiments and MD simulations conducted was observed for both HIV-1 protease-drug complexes. In addition, the mechanism governing the decrease in the binding affinity of PI in the presence of I50V mutation was also explored to provide insights pertaining to the design of the next generation of anti-HIV drugs.

Some practical approaches to treating electrostatic polarization of proteins

Conspectus Electrostatic interaction plays a significant role in determining many properties of biomolecules, which exist and function in aqueous solution, a highly polar environment. For example, proteins are composed of amino acids with charged, polar, and nonpolar side chains and their specific electrostatic properties are fundamental to the structure and function of proteins. An important issue that arises in computational study of biomolecular interaction and dynamics based on classical force field is lack of polarization. Polarization is a phenomenon in which the charge distribution of an isolated molecule will be distorted when interacting with another molecule or presented in an external electric field. The distortion of charge distribution is intended to lower the overall energy of the molecular system, which is counter balanced by the increased internal energy of individual molecules due to the distorted charge distributions. The amount of the charge redistribution, which characterizes the polarizability of a molecule, is determined by the level of the charge distortion. Polarization is inherently quantum mechanical, and therefore classical force fields with fixed atomic charges are incapable of capturing this important effect. As a result, simulation studies based on popular force fields, AMBER, CHARMM, etc., lack the polarization effect, which is a widely known deficiency in most computational studies of biomolecules today. Many efforts have been devoted to remedy this deficiency, such as adding additional movable charge on the atom, allowing atomic charges to fluctuate, or including induced multipoles. Although various successes have been achieved and progress at various levels has been reported over the past decades, the issue of lacking polarization in force field based simulations is far from over. For example, some of these methods do not always give converged results, and other methods require huge computational cost. This Account reviews recent work on developing polarized and polarizable force fields based on fragment quantum mechanical calculations for proteins. The methods described here are based on quantum mechanical calculations of proteins in solution, but with a different level of rigor and different computational efficiency for the molecular dynamics applications. In the general approach, a fragment quantum mechanical calculation for protein with implicit solvation is carried out to derive a polarized protein-specific charge (PPC) for any given protein structure. The PPC correctly reflects the polarization state of the protein in a given conformation, and it can also be dynamically changed as the protein changes conformation in dynamics simulations. Another approach that is computationally more efficient is the effective polarizable bond method in which only polar bonds or groups can be polarized and their polarizabilities are predetermined from quantum mechanical calculations of these groups in external electric fields. Both methods can be employed for applications in various situations by taking advantage of their unique features.

Protein's native structure is dynamically stabilized by electronic polarization

In this paper, molecular dynamics (MD) simulations were performed for a number of benchmark proteins using both the standard assisted model building with energy refinement (AMBER) charge and the dynamically adjusted polarized protein-specific charge (DPPC) from quantum fragment calculations to provide accurate electrostatic interactions. Our result shows that proteins' dynamic structures drifted away from the native structures in simulations under standard (nonpolarizable) AMBER force field. For comparison, proteins' native structures were dynamically stable after a long time simulation under DPPC. The free energy landscape reveals that the native structure is the lowest energy conformation under DPPC, while it is not under standard AMBER charge. To further investigate the polarization effect on the stability of native structures of proteins, we restarted from some decoy structures generated from simulations using standard AMBER charges and then carried out further MD simulation using DPPC to refine those structures. Our study shows that the native structures from these decoy structures can be mostly recovered using DPPC and that the dynamic structures with the highest population in cluster analysis are in close agreement with the corresponding native structures. The current study demonstrates the importance of electronic polarization of protein in stabilizing the native structure.

Importance of polarization effect in the study of metalloproteins: application of polarized protein specific charge scheme in predicting the reduction potential of azurin

Molecular dynamics (MD) simulation is commonly used in the study of protein dynamics, and in recent years, the extension of MD simulation to the study of metalloproteins is gaining much interest. Choice of force field is crucial in MD studies, and the inclusion of metal centers complicates the process of accurately describing the electrostatic environment that surrounds the redox centre. Herein, we would like to explore the importance of including electrostatic contribution from both protein and solvent in the study of metalloproteins. MD simulations with the implementation of thermodynamic integration will be conducted to model the reduction process of azurin from Pseudomonas aeruginosa. Three charge schemes will be used to derive the partial charges of azurin. These charge schemes differ in terms of the amount of immediate environment, respective to copper, considered during charge fitting, which ranges from the inclusion of copper and residues in the first coordination sphere during density functional theory charge fitting to the comprehensive inclusion of protein and solvent effect surrounding the metal centre using polarized protein-specific charge scheme. From the simulations conducted, the relative reduction potential of the mutated azurins respective to that of wild-type azurin (ΔEcal) were calculated and compared with experimental values. The ΔEcal approached experimental value with increasing consideration of environmental effect hence substantiating the importance of polarization effect in the study of metalloproteins. This study also attests the practicality of polarized protein-specific charge as a computational tool capable of incorporating both protein environment and solvent effect into MD simulations.

Electronic polarization stabilizes tertiary structure prediction of HP-36

Molecular dynamic (MD) simulations with both implicit and explicit solvent models have been carried out to study the folding dynamics of HP-36 protein. Starting from the extended conformation, the secondary structure of all three helices in HP-36 was formed in about 50 ns and remained stable in the remaining simulation. However, the formation of the tertiary structure was difficult. Although some intermediates were close to the native structure, the overall conformation was not stable. Further analysis revealed that the large structure fluctuation of loop and hydrophobic core regions was devoted mostly to the instability of the structure during MD simulation. The backbone root-mean-square deviation (RMSD) of the loop and hydrophobic core regions showed strong correlation with the backbone RMSD of the whole protein. The free energy landscape indicated that the distribution of main chain torsions in loop and turn regions was far away from the native state. Starting from an intermediate structure extracted from the initial AMBER simulation, HP-36 was found to generally fold to the native state under the dynamically adjusted polarized protein-specific charge (DPPC) simulation, while the peptide did not fold into the native structure when AMBER force filed was used. The two best folded structures were extracted and taken into further simulations in water employing AMBER03 charge and DPPC for 25 ns. Result showed that introducing polarization effect into interacting potential could stabilize the near-native protein structure.

Virtual screening of PRK1 inhibitors: ensemble docking, rescoring using binding free energy calculation and QSAR model development

Protein kinase C Related Kinase 1 (PRK1) has been shown to be involved in the regulation of androgen receptor signaling and has been identified as a novel potential drug target for prostate cancer therapy. Since there is no PRK1 crystal structure available to date, multiple PRK1 homology models were generated in order to address the protein flexibility. An in-house library of compounds tested on PRK1 was docked into the ATP binding site of the generated models. In most cases a correct pose of the inhibitors could be identified by ensemble docking, while there is still a challenge of finding a reasonable scoring function that is able to rank compounds according to their biological activity. We estimated the binding free energy for our data set of structurally diverse PRK1 inhibitors using the MM-PB(GB)SA and QM/MM-GBSA methods. The obtained results demonstrate that a correlation between calculated binding free energies and experimental IC50 values was found to be usually higher than using docking scores. Furthermore, the developed approach was tested on a set of diverse PRK1 inhibitors taken from literature, which resulted in a significant correlation. The developed method is computationally inexpensive and can be applied as a postdocking filter in virtual screening as well as for optimization of PRK1 inhibitors in order to prioritize compounds for further biological characterization.

Ab initio folding of extended α-helix: a theoretical study about the role of electrostatic polarization in the folding of helical structures

In this work, we report the ab initio folding of three different extended helical peptides namely 2khk, N36, and C34 through conventional molecular dynamics simulation at room temperature using implicit solvation model. Employing adaptive hydrogen bond specific charge (AHBC) scheme to account for the polarization effect of hydrogen bonds established during the simulation, the effective folding of the three extended helices were observed with best backbone RMSDs in comparison to the experimental structures over the helical region determined to be 1.30 Å for 2khk, 0.73 Å for N36 and 0.72 Å for C34. In this study, 2khk will be used as a benchmark case serving as a means to compare the ability of polarized (AHBC) and nonpolarized force field in the folding of an extended helix. Analyses conducted revealed the ability of the AHBC scheme in effectively folding the extended helix by promoting helix growth through the stabilization of backbone hydrogen bonds upon formation during the folding process. Similar observations were also noted when AHBC scheme was employed during the folding of C34 and N36. However, under Amber03 force field, helical structures formed during the folding of 2khk was not accompanied by stabilization thus highlighting the importance of electrostatic polarization in the folding of helical structures.

Predicting mutation-induced Stark shifts in the active site of a protein with a polarized force field

The electric field inside a protein has a significant effect on the protein structure, function, and dynamics. Recent experimental developments have offered a direct approach to measure the electric field by utilizing a nitrile-containing inhibitor as a probe that can deliver a unique vibration to the specific site of interest in the protein. The observed frequency shift of the nitrile stretching vibration exhibits a linear dependence on the electric field at the nitrile site, thus providing a direct measurement of the relative electric field. In the present work, molecular dynamics simulations were carried out to compute the electric field shift in human aldose reductase (hALR2) using a polarized protein-specific charge (PPC) model derived from fragment-based quantum-chemistry calculations in implicit solvent. Calculated changes of electric field in the active site of hALR2 between the wild type and mutants were directly compared with measured vibrational frequency shifts (Stark shifts). Our study demonstrates that the Stark shifts calculated using the PPC model are in much better agreement with the experimental data than widely used nonpolarizable force fields, indicating that the electronic polarization effect is important for the accurate prediction of changes in the electric field inside proteins.

SAR and Computer-Aided Drug Design Approaches in the Discovery of Peroxisome Proliferator-Activated Receptor γ Activators: A Perspective

Activators of PPARγ, Troglitazone (TGZ), Rosiglitazone (RGZ), and Pioglitazone (PGZ) were introduced for treatment of Type 2 diabetes, but TGZ and RGZ have been withdrawn from the market along with other promising leads due cardiovascular side effects and hepatotoxicity. However, the continuously improving understanding of the structure/function of PPARγ and its interactions with potential ligands maintain the importance of PPARγ as an antidiabetic target. Extensive structure activity relationship (SAR) studies have thus been performed on a variety of structural scaffolds by various research groups. Computer-aided drug discovery (CADD) approaches have also played a vital role in the search and optimization of potential lead compounds. This paper focuses on these approaches adopted for the discovery of PPARγ ligands for the treatment of Type 2 diabetes. Key concepts employed during the discovery phase, classification based on agonistic character, applications of various QSAR, pharmacophore mapping, virtual screening, molecular docking, and molecular dynamics studies are highlighted. Molecular level analysis of the dynamic nature of ligand-receptor interaction is presented for the future design of ligands with better potency and safety profiles. Recently identified mechanism of inhibition of phosphorylation of PPARγ at SER273 by ligands is reviewed as a new strategy to identify novel drug candidates.

Dynamical stability and assembly cooperativity of β-sheet amyloid oligomers--effect of polarization

The soluble intermediate oligomers of amyloidogenic proteins are suspected to be more cytotoxic than the mature fibrils in neurodegenerative disorders. Here, the dynamic stability and assembly cooperativity of a model oligomer of human islet amyloid polypeptide (hIAPP) segments were explored by means of all-atom molecular dynamics (MD) simulations under different force fields including AMBER99SB, OPLS, and polarized protein-specific charge (PPC) model. Simulation results show that the dynamic stability of β-sheet oligomers is seriously impacted by electrostatic polarization. Without inclusion of polarization (simulation under standard AMBER and OPLS force field), the β-sheet oligomers are dynamically unstable during MD simulation. For comparison, simulation results under PPC give significantly more stable dynamical structures of the oligomers. Furthermore, calculation of electrostatic interaction energy between the neighboring β strands with an approximate polarizable method produces energetic evidence for cooperative assembly of β-strand oligomers. This result supports a picture of downhill-like cooperative assembly of β strands during fibrillation process. The present study demonstrates the critical role of polarization in dynamic stability and assembly cooperativity of β-sheet-rich amyloid oligomers.

Design novel dual agonists for treating type-2 diabetes by targeting peroxisome proliferator-activated receptors with core hopping approach

Owing to their unique functions in regulating glucose, lipid and cholesterol metabolism, PPARs (peroxisome proliferator-activated receptors) have drawn special attention for developing drugs to treat type-2 diabetes. By combining the lipid benefit of PPAR-alpha agonists (such as fibrates) with the glycemic advantages of the PPAR-gamma agonists (such as thiazolidinediones), the dual PPAR agonists approach can both improve the metabolic effects and minimize the side effects caused by either agent alone, and hence has become a promising strategy for designing effective drugs against type-2 diabetes. In this study, by means of the powerful “core hopping” and “glide docking” techniques, a novel class of PPAR dual agonists was discovered based on the compound GW409544, a well-known dual agonist for both PPAR-alpha and PPAR-gamma modified from the farglitazar structure. It was observed by molecular dynamics simulations that these novel agonists not only possessed the same function as GW409544 did in activating PPAR-alpha and PPAR-gamma, but also had more favorable conformation for binding to the two receptors. It was further validated by the outcomes of their ADME (absorption, distribution, metabolism, and excretion) predictions that the new agonists hold high potential to become drug candidates. Or at the very least, the findings reported here may stimulate new strategy or provide useful insights for discovering more effective dual agonists for treating type-2 diabetes. Since the “core hopping” technique allows for rapidly screening novel cores to help overcome unwanted properties by generating new lead compounds with improved core properties, it has not escaped our notice that the current strategy along with the corresponding computational procedures can also be utilized to find novel and more effective drugs for treating other illnesses.

Effect of interprotein polarization on protein-protein binding energy

Molecular dynamics simulation in explicit water for the binding of the benchmark barnase-barstar complex was carried out to investigate the effect polarization of interprotein hydrogen bonds on its binding free energy. Our study is based on the AMBER force field but with polarized atomic charges derived from fragment quantum mechanical calculation for the protein complex. The quantum-derived atomic charges include the effect of polarization of interprotein hydrogen bonds, which was absent in the standard force fields that were used in previous theoretical calculations of barnase-barstar binding energy. This study shows that this polarization effect impacts both the static (electronic) and dynamic interprotein electrostatic interactions and significantly lowers the free energy of the barnase-barstar complex.

Electrostatic polarization is critical for the strong binding in streptavidin-biotin system

The origin of strong affinity of biotin and its analogs binding to (strept)avidin is still the subject of an ongoing controversy. In this work, thermodynamic integration is carried out to study of the difference of binding free energies between biotin and iminobiotin to streptavidin. Three atomic charge schemes are implemented and compared. One is the traditional AMBER charge, and the other two, termed the polarized protein-specific charge, are based on a linear scaling quantum mechanical method and a continuous solvation model and have polarization effect partially or fully included. The result indicates that when nonpolarized AMBER force field is applied, the result is much underestimated. When electronic polarization is gradually included, the difference of binding affinity increases along with it. Using the linear-response approximation to eliminate the error in self-charging process, the corrected binding affinity agrees well with the experimental observation. This study is direct evidence indicating that polarization effect is critical for the strong binding in streptavidin-biotin system.

Folding of a helix is critically stabilized by polarization of backbone hydrogen bonds: study in explicit water

Multiple single-trajectory molecular dynamics (MD) simulation at room temperature (300 K) in explicit water was carried out to study the folding dynamics of an α-helix (PDB 2I9M ) using a polarized charge scheme that includes electronic polarization of backbone hydrogen bonds. Starting from an extended conformation, the 17-residue peptide was successfully folded into the native structure (α-helix) between 80 and 130 ns with a root-mean-square deviation of ~1.0 Å. Analysis of the time-dependent trajectories revealed that helix formation of the peptide started at the terminals and progressed toward the center of the peptide. For comparison, MD trajectories generated under various versions of standard AMBER force fields failed to show any significant or stable helix formation in our simulation. Our result shows clear evidence that the electronic polarization of backbone hydrogen bonds energetically stabilizes the helix formation and is critical to the stable folding of the short helix structure.

Measurement and applications of long-range heteronuclear scalar couplings: recent experimental and theoretical developments

The use of long-range heteronuclear couplings, in association with (1)H-(1)H scalar couplings and NOE restraints, has acquired growing importance for the determination of the relative stereochemistry, and structural and conformational information of organic and biological molecules. However, the routine use of such couplings is hindered by the inherent difficulties in their measurement. Prior to the advancement in experimental techniques, both long-range homo- and heteronuclear scalar couplings were not easily accessible, especially for very large molecules. The development of a large number of multidimensional NMR experimental methodologies has alleviated the complications associated with the measurement of couplings of smaller strengths. Subsequent application of these methods and the utilization of determined J-couplings for structure calculations have revolutionized this area of research. Problems in organic, inorganic and biophysical chemistry have also been solved by utilizing the short- and long-range heteronuclear couplings. In this minireview, we discuss the advantages and limitations of a number of experimental techniques reported in recent times for the measurement of long-range heteronuclear couplings and a few selected applications of such couplings. This includes the study of medium- to larger-sized molecules in a variety of applications, especially in the study of hydrogen bonding in biological systems. The utilization of these couplings in conjunction with theoretical calculations to arrive at conclusions on the hyperconjugation, configurational analysis and the effect of the electronegativity of the substituents is also discussed.

Direct evaluation of individual hydrogen bond energy in situ in intra- and intermolecular multiple hydrogen bonds system

The results of evaluating the individual hydrogen bond (H-bond) strength are expected to be helpful for the rational design of new strategies for molecular recognition or supramolecular assemblies. Unfortunately, there is few obvious and unambiguous means of evaluating the energy of a single H-bond within a multiple H-bonds system. We present a local analytic model, ABEEMσπ H-bond energy (HBE) model based on ab initio calculations (MP2) as benchmark, to directly and rapidly evaluate the individual HBE in situ in inter- and intramolecular multiple H-bonds system. This model describes the HBE as the sum of electrostatic and van der Waals (vdW) interactions which all depend upon the geometry and environment, and the ambient environment of H-bond in the model is accounted fairly. Thus, it can fairly consider the cooperative effect and secondary effect. The application range of ABEEMσπ HBE model is rather wide. This work has discussed the individual H-bond in DNA base pair and protein peptide dimers. The results indicate that the interactions among donor H atom, acceptor atom as well as those atoms connected to them with 1,2 or 1,3 relationships are all important for evaluating the HBE, although the interaction between the donor H atom and the acceptor atom is large. Furthermore, our model quantitatively indicates the polarization ability of N, O, and S in a new style, and gives the percentage of the polarization effect in HBE, which can not be given by fixed partial charge force field.