A QM protein–ligand investigation of antipsychotic drugs with the dopamine D2 Receptor (D2R)

The dopamine D2 Receptor (D2R) is a member of the G-Protein-Coupled Receptor family and plays a critical role in neurotransmission activities in the human brain. Dysfunction in dopamine receptor signaling may lead to mental health illnesses such as schizophrenia and Parkinson’s disease. D2R is the target protein of the commonly used antipsychotic drugs such as risperidone, clozapine, aripiprazole, olanzapine, ziprasidone, and quetiapine. Due to their significant side effects and non-selective profiles, the discovery of novel drugs has become a challenge for researchers working in this field. Recently, our group has focused on the interactions of these drug molecules in the active site of the D2R using different in silico approaches. We here compare the performances of different approaches in estimating the drug binding affinities using quantum chemical approaches. Conformations of drug molecules (ligands) at the binding site of the D2R taken from the preliminary docking studies and molecular dynamics simulations were used to generate protein–ligand interaction models. In a first approach, the BSSE-corrected interaction energies of the ligands with the most critical amino acid Asp114 and with the other amino acids closest to ligands in the binding cavity were calculated separately by density functional theory method in implicit water environment at the M06-2X/6-31 g(d,p) level of the theory. In a second approach, ligand binding affinities were calculated by taking into consideration not only the interaction energies but also deformation and desolvation energies of ligands with surrounding amino acid residues, in a radius of 5 Å of the protein-bound ligand. The quantum mechanically obtained results were compared with the experimentally obtained binding affinity values. We concluded that although H-bond interactions of ligands with Asp114 are the most dominant interaction in the binding site, if van der Waals and steric interactions of ligands which have cumulative effect on the ligand binding are not included in the calculations, the interaction energies are overestimated.     

Theoretical and Molecular Docking Study of Ketoconazole on Heptakis(2,3,6‐tri‐O‐methyl)‐β‐cyclodextrin as Chiral Selector

A molecular docking study, using molecular mechanics calculations with AutoDock and semi‐empirical PM3 calculations, was used to predict the enantiodiscrimination of heptakis(2,3,6‐tri‐O‐methyl)‐β‐cyclodextrin (TMβCD) and ketoconazole (KTZ) enantiomers. A Density Functional Theory (DFT) single‐point calculation at the level of B3LYP/6‐311G (d,p) was performed for the PM3‐optimized complexes to obtain more accurate binding energy and the electronic structures of the complexes. The difference in energies of the inclusion complexes between the KTZ enantiomers and TMβCD is probably a measure of chiral discrimination, which results in the separation of the enantiomers as observed in the experimental studies. Chirality 28:209–214, 2016. © 2015 Wiley Periodicals, Inc.     

Time dependent – density functional theory characterization of organic dyes for dye-sensitized solar cells

We aim at providing better insight into the parameters that govern the intramolecular charge transfer (ICT) and photo-injection processes in dyes for dye-sensitised solar cells (DSSC). Density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations are utilized to study the geometry, electronic structure, electrostatic potential (ESP) and absorption spectrum, for a representative donor-π bridge-acceptor (D–π–A) dye for DSSC. The coplanar geometry of the dye (D1) facilitates strong conjugation and considerable delocalization originating the π CT interaction from donor to acceptor orbitals and the hyper-conjugative interactions involving Rydberg states. A model simulating the adsorption of the dye on the TiO₂ surface is utilized to estimate binding energies. The effect of fluorine substituents in the π-spacer on the quantum efficiency of DSSCs was investigated. Gibb’s free energy values, redox potentials, excited state lifetime, non-linear optical properties (NLO) and driving forces for D1 and its fluorinated derivatives were computed.     

Density functional theory calculations on entire proteins for free energies of binding: Application to a model polar binding site

In drug optimization calculations, the molecular mechanics Poisson‐Boltzmann surface area (MM‐PBSA) method can be used to compute free energies of binding of ligands to proteins. The method involves the evaluation of the energy of configurations in an implicit solvent model. One source of errors is the force field used, which can potentially lead to large errors due to the restrictions in accuracy imposed by its empirical nature. To assess the effect of the force field on the calculation of binding energies, in this article we use large‐scale density functional theory (DFT) calculations as an alternative method to evaluate the energies of the configurations in a “QM‐PBSA” approach. Our DFT calculations are performed with a near‐complete basis set and a minimal parameter implicit solvent model, within the self‐consistent calculation, using the ONETEP program on protein–ligand complexes containing more than 2600 atoms. We apply this approach to the T4‐lysozyme double mutant L99A/M102Q protein, which is a well‐studied model of a polar binding site, using a set of eight small aromatic ligands. We observe that there is very good correlation between the MM and QM binding energies in vacuum but less so in the solvent. The relative binding free energies from DFT are more accurate than the ones from the MM calculations, and give markedly better agreement with experiment for six of the eight ligands. Furthermore, in contrast to MM‐PBSA, QM‐PBSA is able to correctly predict a nonbinder. Proteins 2014; 82:3335–3346. © 2014 Wiley Periodicals, Inc.     

Computational design of a lipase for catalysis of the Diels-Alder reaction

Combined molecular docking, molecular dynamics (MD) and density functional theory (DFT) studies have been employed to study catalysis of the Diels-Alder reaction by a modified lipase. Six variants of the versatile enzyme Candida Antarctica lipase B (CALB) have been rationally engineered in silico based on the specific characteristics of the pericyclic addition. A kinetic analysis reveals that hydrogen bond stabilization of the transition state and substrate binding are key components of the catalytic process. In the case of substrate binding, which has the greater potential for optimization, both binding strength and positioning of the substrates are important for catalytic efficiency. The binding strength is determined by hydrophobic interactions and can be tuned by careful selection of solvent and substrates. The MD simulations show that substrate positioning is sensitive to cavity shape and size, and can be controlled by a few rational mutations. The well-documented S105A mutation is essential to enable sufficient space in the vicinity of the oxyanion hole. Moreover, bulky residues on the edge of the active site hinders the formation of a sandwich-like nearattack conformer (NAC), and the I189A mutation is needed to obtain enough space above the face of the α,β-double bond on the dienophile. The double mutant S105A/I189A performs quite well for two of three dienophiles. Based on binding constants and NAC energies obtained from MD simulations combined with activation energies from DFT computations, relative catalytic rates (v _cat /v _uncat ) of up to 103 are predicted.     

Free energies of binding of <Emphasis Type="Italic">R</Emphasis>- and <Emphasis Type="Italic">S</Emphasis>-propranolol to wild-type and F483A mutant cytochrome P450 2D6 from molecular dynamics simulations

Detailed molecular dynamics (MD) simulations have been performed to reproduce and rationalize the experimental finding that the F483A mutant of CYP2D6 has lower affinity for R-propranolol than for S-propranolol. Wild-type (WT) CYP2D6 does not show this stereospecificity. Four different approaches to calculate the free energy differences have been investigated and were compared to the experimental binding data. From the differences between calculations based on forward and backward processes and the closure of thermodynamic cycles, it was clear that not all simulations converged sufficiently. The approach that calculates the free energies of exchanging R-propranolol with S-propranolol in the F483A mutant relative to the exchange free energy in WT CYP2D6 accurately reproduced the experimental binding data. Careful inspection of the end-points of the MD simulations involved in this approach, allowed for a molecular interpretation of the observed differences. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Analysis of cyclosporin A and a set of analogs as inhibitors of a T. cruzi cyclophilin by docking and molecular dynamics

Cyclophilins (CyPs) are enzymes involved in protein folding. In Trypanosoma cruzi (T. cruzi), the most abundantly expressed CyP is the isoform TcCyP19. It has been shown that TcCyP19 is inhibited by the immunosuppressive drug cyclosporin A (CsA) and analogs, which also proved to have potent trypanosomicidal activity in vitro. In this work, we continue and expand a previous study on the molecular interactions of CsA, and a set of analogs modeled in complexes with TcCyP19. The modeled complexes were used to evaluate binding free energies by molecular dynamics (MD), applying the Linear Interaction Energy (LIE) method. In addition, putative binding sites were identified by molecular docking. In our analysis, the binding free energy calculations did not correlate with experimental data. The heterogeneity of the non-bonded energies and the variation in the pattern of hydrogen bonds suggest that the systems may not be suitable for the application of the LIE method. Further, the docking calculations identified two other putative binding sites with comparable scoring energies to the active site, a fact that may also explain the lack of correlation found. Kinetic experiments are needed to confirm or reject the multiple binding sites hypothesis. In the meantime, MD simulations at the alternative sites, employing other methods to compute binding free energies, might be successful at finding good correlations with the experimental data.     

Study of a ligand complexed with Cdk2/Cdk4 by computer simulation

Cyclin-dependent kinases (Cdks) play important roles in the regulation of the cell cycle. Their inhibitors have entered clinical trials to treat cancer. Very recently, Davis et al. (Nat Struct Biol 9:745–749, 2002) have found a ligand NU6102, which has a high affinity with cyclin-dependent kinase 2 (K _i =6 nM) but a low affinity with cyclin-dependent kinase 4 (K _i =1,600 nM). To understand the selectivity, we use homology modeling, molecular docking, molecular dynamics and free-energy calculations to analyze the interactions. A rational 3D model of the Cdk4–NU6102 complex is built. Asp86 is a key residue that recognizes NU6102 more effectively with Cdk2 rather than Cdk4. Good binding free energies are obtained. Energetic analysis reveals that van der Waals interaction and nonpolar contributions to solvent are favorable in the formation of complexes and the sulfonamide group of the ligand plays a crucial role for binding selectivity between Cdk2 and Cdk4. Figure Two-dimensional representative for the interacting model of NU6102 complexed with the Cdk4 from a predicted structure by LIGPLOT.      

Theoretical study of the ligand-CYP2B4 complexes: effect of structure on binding free energies and heme spin state

The molecular origins of temperature-dependent ligand-binding affinities and ligand-induced heme spin state conversion have been investigated using free energy analysis and DFT calculations for substrates and inhibitors of cytochrome P450 2B4 (CYP2B4), employing models of CYP2B4 based on CYP2C5(3LVdH)/CYP2C9 crystal structures, and the results compared with experiment. DFT calculations indicate that large heme-ligand interactions (ca. -15 kcal/mol) are required for inducing a high to low spin heme transition, which is correlated with large molecular electrostatic potentials (approximately -45 kcal/mol) at the ligand heteroatom. While type II ligands often contain oxygen and nitrogen heteroatoms that ligate heme iron, DFT results indicate that BP and MF heme complexes, with weak substrate-heme interactions (ca. -2 kcal/mol), and modest MEPS minima (>-35 kcal/mol) are high spin. In contrast, heme complexes of the CYP2B4 inhibitor, 4PI, the product of benzphetamine metabolism, DMBP, and water are low spin, have substantial heme-ligand interaction energies (<-15 kcal/mol) and deep MEPS minima (<-45 kcal/mol) near their heteroatoms. MMPBSA analysis of MD trajectories were made to estimate binding free energies of these ligands at the heme binding site of CYP2B4. In order to initially assess the realism of this approach, the binding free energy of 4PI inhibitor was computed and found to be a reasonable agreement with experiment: -7.7 kcal/mol [-7.2 kcal/mol (experiment)]. BP was determined to be a good substrate [-6.3 kcal/mol (with heme-ligand water), -7.3 kcal/mol (without ligand water)/-5.8 kcal/mol (experiment)], whereas the binding of MF was negligible, with only marginal binding binding free energy of -1.7 kcal/mol with 2-MF bound [-3.8 kcal/mol (experiment)], both with and without retained heme-ligand water. Analysis of the free energy components reveal that hydrophobic/nonpolar contributions account for approximately 90% of the total binding free energy of these substrates and are the source of their differential and temperature-dependent CYP2B4 binding. The results indicate the underlying origins of the experimentally observed differential binding affinities of BP and MF, and indicate the plausibility of the use of models derived from moderate sequence identity templates in conjunction with approximate free energy methods in the estimation of ligand-P450 binding affinities.     

Density functional theory study on the interaction between keto-9H guanine and aspartic acid

A theoretical study was performed using density functional theory (DFT) to investigate hydrogen bonding interactions in signature complexes formed between keto-9H guanine (Gua) and aspartic acid (Asp) at neutral pH. Optimized geometries, binding energies and the theoretical IR spectra of guanine, aspartic acid and their corresponding complexes (Gua-Asp) were calculated using the B3LYP method and the 6-31+G(d) basis set. Stationary points found to be at local minima on the potential energy surface were verified by second derivative harmonic vibrational frequency calculations at the same level of theory. AIM theory was used to analyze the hydrogen bonding characteristics of these DNA base complex systems. Our results show that the binding motif for the most stable complex is strikingly similar to a Watson-Crick motif observed in the guanine-cytosine base pair. We have found a range of hydrogen bonding interactions between guanine and aspartic acid in the six complexes. This was further verified by theoretical IR spectra of ω(C-H---O-H) cm⁻¹ stretches for the Gua-Asp complexes. The electron density plot indicates strong hydrogen bonding as shown by the 2p _z dominant HOMO orbital character.     

Do electrostatic interactions destabilize protein–nucleic acid binding?

The negatively charged phosphates of nucleic acids are often paired with positively charged residues upon binding proteins. It was thus counter-intuitive when previous Poisson-Boltzmann (PB) calculations gave positive energies from electrostatic interactions, meaning that they destabilize protein-nucleic acid binding. Our own PB calculations on protein-protein binding have shown that the sign and the magnitude of the electrostatic component are sensitive to the specification of the dielectric boundary in PB calculations. A popular choice for the boundary between the solute low dielectric and the solvent high dielectric is the molecular surface; an alternative is the van der Waals (vdW) surface. In line with results for protein-protein binding, in this article, we found that PB calculations with the molecular surface gave positive electrostatic interaction energies for two protein-RNA complexes, but the signs are reversed when the vdW surface was used. Therefore, whether destabilizing or stabilizing effects are predicted depends on the choice of the dielectric boundary. The two calculation protocols, however, yielded similar salt effects on the binding affinity. Effects of charge mutations differentiated the two calculation protocols; PB calculations with the vdW surface had smaller deviations overall from experimental data.     

MP2, DFT and DFT-D study of the dimers of diazanaphthalenes: a comparative study of their structures,stabilisation and binding energies

Diazanaphthalenes (DAPs) are a broad class of N-heteroaromatic compounds with several technological and biological applications. Some of these applications are attributed to the ability of DAP molecules to form associated dimers through non-covalent interactions. A study of the types and strength of the interactions involved is crucial for understanding the preferred geometries and energetics of the dimers. In this study, the dimers of five DAPs are investigated by means of Møller–Plesset second order perturbation theory, hybrid meta-GGA [density functional theory methods (DFT): DFT/MPWB1K, DFT/M05-2X and DFT/M06-2X] and DFT dispersion-corrected (DFT-D/ωB97XD) methods to elucidate their dimers' preferred geometries, relative energies and nature of the interactions between monomer units. The results indicate that the monomer units of the dimers are held by either intermolecular hydrogen bonds or π…π stacking interactions, and that the preferred dimers are those in which the monomer units interact through π…π stacking interactions. A comparison across structures suggests that the position of the N atom in the ring has significant role in determining the relative energy and binding strength of the dimers. A comparison among the different methods utilised for the study indicates that DFT/M06-2X method provides binding energies that are close to those of DFT-CCSD(T) correction scheme and could therefore be considered as the best method for describing the binding properties of DAP dimers.     

Molecular dynamics,free energy,and SPR analyses of the interactions between the SH2 domain of Grb2 and ErbB phosphotyrosyl peptides

We studied the interactions between the SH2 domain of growth factor receptor binding protein 2 (Grb2) and ErbB receptor-derived phosphotyrosyl peptides using molecular dynamics, free energy calculations, and surface plasmon resonance (SPR) analysis. Binding free energies for nine phosphotyrosyl peptides were calculated using the MM-PBSA continuum solvent method, and excellent qualitative agreement with the SPR experimental data, with a correlation coefficient of 0.92, was obtained. Consistent with previous experimental findings, phosphotyrosyl peptides with the consensus sequence pYXNX showed favorable binding affinity for the Grb2. Unexpectedly, phosphotyrosyl peptides with the consensus sequence pYQQD, which had not shown any specific binding affinity for the Grb2 in earlier studies, also showed favorable binding affinity for the Grb2 in our experimental and computational analyses. Component analysis of the calculated binding free energies revealed that van der Waals interaction between the Grb2 and the phosphotyrosyl peptide was the dominant factor for specificity and binding affinity. These results indicate that current methods of estimating binding free energies are efficient for obtaining important information about protein-protein interactions, which are essential for the transmission of signals in cellular signaling pathways.     

Binding free energy calculations of adenosine deaminase inhibitors

The interactions between four inhibitors and adenosine deaminase (ADA) were examined by calculating their binding free energies after molecular dynamics simulations. A bonded model was used to represent the electrostatic potentials of the zinc coordination site. The charge distribution of the model was derived by using a two-stage electrostatic potential fitting calculations. The calculated binding free energies are in good agreement with the experimental data and the ranking of binding affinities is well reproduced. Notably, our findings suggest that non-polar contributions play an important role for ADA-inhibitor interactions.     

Computational Evaluation of Intermolecular Interactions of a Universal Base 3-Nitropyrrole in Stacked Dimers and DNA Duplexes

Abstract

The stacking interactions between a universal base of 3-nitropyrrole (3NP) and four canonical nucleobases were studied by means of ab initio molecular orbital calculations. The stabilities of the complexes are comparable to those of the stacked dimers of canonical bases reported previously. The detailed analysis of the interaction energies revealed the importance of the dipole-dipole interaction included in the Hartree-Fock terms to determine the geometry dependence of the stacking energies. It was also clarified that the dispersion energies included in the electron-correlation terms were essential to obtain adequate stabilities. The contribution of the nitro group was evaluated by the comparative studies of pyrrole and 3NP. The increased molecular dipole moment and surface are expected to account for the enhancement of the stability of the stacked dimers containing 3NP. The force field parameters required for calculation of the molecular mechanics of 3NP were obtained for 3NP on the basis of these molecular orbital calculations. The energy-minimized structures obtained by the molecular mechanics calculations of 3NP accorded with those obtained by the molecular orbital calculations described above. A DNA duplex structure containing 3NP-A, 3NP-T, or 3NP-C was calculated by use of these force field parameters. In the case of 3NP-A, the computationally calculated structure was in good agreement with that previously determined by use of ¹H-NMR except for the orientation of the nitro group.     

Understanding microscopic binding of macrophage migration inhibitory factor with phenolic hydrazones by molecular docking, molecular dynamics simulations and free energy calculations

Macrophage migration inhibitory factor (MIF), an immunoregulatory protein, is a potential target for a number of inflammatory diseases. In the current work, the interactions between MIF and a series of phenolic hydrazones were studied by molecular docking, molecular dynamics (MD) simulations, binding free energy calculations, and binding energy decomposition analysis to determine the structural requirement for achieving favorable biological activity of phenolic hydrazones. First, molecular docking was used to predict the binding modes of inhibitors in the binding site of MIF. The good correlation between the predicted docking scores and the experimental activities shows that the binding conformations of the inhibitors in the active site of MIF are well predicted. Moreover, our results suggest that the flexibility of MIF is essential in ligand binding process. Then, MD simulations and MM/GBSA free energy calculations were employed to determine the dynamic binding process and compare the binding modes of the inhibitors with different activities. The predicted binding free energies given by MM/GBSA are not well correlated with the experimental activities for the two subsets of the inhibitors; however, for each subset, a good correlation between the predicted binding free energies and the experimental activities is achieved. The MM/GBSA free energy decomposition analysis highlights the importance of hydrophobic residues for the MIF binding of the studied inhibitors. Based on the essential factors for MIF-inhibitor interactions derived from the theoretical predictions, some derivatives were designed and the higher inhibitory activities of several candidates were confirmed by molecular docking studies. The structural insights obtained from our study are useful for designing potent inhibitors of MIF.     

A simple model for predicting the free energy of binding between anthracycline antibiotics and DNA

A theoretical model for predicting the free energy of binding between anthracycline antibiotics and DNA was developed using the electron density functional (DFT) and molecular mechanics (MM) methods. Partial DFT-ESP charges were used in calculating the MM binding energies for complexes formed between anthracycline antibiotics and oligodeoxynucleotides. These energies were then compared with experimental binding free energies. The good correlation between the experimental and theoretical energies allowed us to propose a model for predicting the binding free energy for derivatives of anthracycline antibiotics and for quickly screening new anthracycline derivatives.     

Relating Frontier Orbital Energies from Voltammetry and Photoelectron Spectroscopy to the Open‐Circuit Voltage of Organic Solar Cells

For 19 diketopyrrolopyrrole polymers, the highest occupied molecular orbital (HOMO) energies are determined from i) the oxidation potential with square‐wave voltammetry (SWV), ii) the ionization potential using ultraviolet photoelectron spectroscopy (UPS), and iii) density functional theory (DFT) calculations. The SWV HOMO energies show an excellent linear correlation with the open‐circuit voltage (V_oc) of optimized solar cells in which the polymers form blends with a fullerene acceptor ([6,6]‐phenyl‐C₆₁‐butyl acid methyl ester or [6,6]‐phenyl‐C₇₁‐butyl acid methyl ester). Remarkably, the slope of the best linear fit is 0.75 ± 0.04, i.e., significantly less than unity. A weaker correlation with V_oc is found for the HOMO energies obtained from UPS and DFT. Within the experimental error, the SWV and UPS data are correlated with a slope close to unity. The results show that electrochemically determined oxidation potentials provide an excellent method for predicting the V_oc of bulk heterojunction solar cells, with absolute deviations less than 0.1 V.