Calculation of the total electrostatic energy of a macromolecular system: solvation energies, binding energies, and conformational analysis

In this report we describe an accurate numerical method for calculating the total electrostatic energy of molecules of arbitrary shape and charge distribution, accounting for both Coulombic and solvent polarization terms. In addition to the solvation energies of individual molecules, the method can be used to calculate the electrostatic energy associated with conformational changes in proteins as well as changes in solvation energy that accompany the binding of charged substrates. The validity of the method is examined by calculating the hydration energies of acetate, methyl ammonium, ammonium, and methanol. The method is then used to study the relationship between the depth of a charge within a protein and its interaction with the solvent. Calculations of the relative electrostatic energies of crystal and misfolded conformations of Themiste dyscritum hemerythrin and the VL domain of an antibody are also presented. The results indicate that electrostatic charge-solvent interactions strongly favor the crystal structures. More generally, it is found that charge-solvent interactions, which are frequently neglected in protein structure analysis, can make large contributions to the total energy of a macromolecular system.     

ATP binding to solubilized (Na+ + K+)-ATPase: The abolition of subunit-subunit interaction and the maximum weight of the nucleotide-binding unit

Membrane-bound (Na⁺ + K⁺)-ATPase from pig kidney outer medulla shows apparent heterogeneity in its ATP-binding site population when assays are carried out in the presence of K⁺. This finding has been interpreted as being due to interaction between (at least) two subunits, each containing an ATP-binding site. Treating the membrane-bound enzyme with the detergent, C₁₂E₈, has been shown to solubilize enzymatically active αβ-protomers. We show that in the dissolved enzyme all ATP-binding sites in the population are identical both in the absence and in the presence of K⁺, which would be consistent with an abolition of subunit-subunit interaction. This supports previous suggestions that enzyme solubilized by C₁₂E₈ is monomeric and that the membrane-bound enzyme is not. Differential extraction of enzyme-containing membranes with C₁₂E₈ yielded preparations with an ATP-binding capacity of up to 5.8 nmol per mg protein, measured by the method of Lowry et al. (Lowry, O.H., Rosebrough, N.J., Farr, A.L., and Randall, R.J. (1951) J. Biol. Chem. 193, 265–275), with bovine serum albumin as standard. Evidence is presented that makes it likely that preparations with an ATP-binding capacity of 7.5 nmol per mg protein (as determined by the above-mentioned assay) will be obtainable. This corresponds to an αβ-protomer molecular weight of 133 000 which approximates closely to the minimum value found in the literature for an αβ-protomer (i.e., 126 000).     

A comparison of stigmatellin conformations, free and bound to the photosynthetic reaction center and the cytochrome bc1 complex

We describe in detail the conformations of the inhibitor stigmatellin in its free form and bound to the ubiquinone-reducing (Q(B)) site of the reaction center and to the ubiquinol-oxidizing (Q(o)) site of the cytochrome bc(1) complex. We present here the first structures of a stereochemically correct stigmatellin in complexes with a bacterial reaction center and the yeast cytochrome bc1 complex. The conformations of the inhibitor bound to the two enzymes are not the same. We focus on the orientations of the stigmatellin side-chain relative to the chromone head group, and on the interaction of the stigmatellin side-chain with these membrane protein complexes. The different conformations of stigmatellin found illustrate the structural variability of the Q sites, which are affected by the same inhibitor. The free rotation about the chi1 dihedral angle is an essential factor for allowing stigmatellin to bind in both the reaction center and the cytochrome bc1 pocket.     

Structural insight into the binding mechanism of ATP to EGFR and L858R,and T790M and L858R/T790 mutants

Abstract

The L858R mutation in EGFR is particularly responsive to small tyrosine kinase inhibitors (TKIs) such as gefitinib and erlotinib. This efficacy decreases due to drug resistance conferred by a second mutation, T790M, which subsequently produces a double mutant, L858R/T790M. Although this resistance was initially attributed to steric blocking by the T790M mutation, experimental studies have demonstrated that differences in the binding affinities of TKIs to T790M and L858R/T790M mutants are more a result of the increased sensitivity of these mutants to ATP than to a decrease in the affinity to TKIs. Regrettably, detailed information at the atomic level on the origins of the increased binding affinity of mutants for ATP is lacking. In this study, we have combined structural data and molecular dynamics simulations with the MMGBSA approach to determine how the L858R, T790M and L858R/T790 mutations impact the binding mechanism of ATP with respect to wild-type EGFR. Structural and energetic analyses provided novel information that helps to explain the increased affinity of ATP to T790M and L858R/T790 mutants with respect to L858R and wild-type systems. In addition, it was observed that dimerization of the wild-type and mutant systems exerts dissimilar effects on the ATP binding affinity characteristic of negative cooperativity.

Communicated by Ramaswamy H. Sarma     

Application of molecular dynamics and free energy perturbation methods to metalloporphyrin-ligand systems II: CO and dioxygen binding to myoglobin.

The protein contribution to the relative binding affinity of the ligands CO and O2 toward myoglobin (Mb) has been simulated using free energy perturbation calculations. The tautomers of the His E7 residue are different for the oxymyoglobin (MbO2) and carboxymyoglobin (MbCO) systems. This was modeled by performing two-step calculations that mutate the ligand and mutate the His E7 tautomers in separate steps. Differences in hydrogen bonding to the O2 and CO ligands were incorporated into the model. The O2 complex was calculated to be 2-3 kcal/mol more stable than the corresponding CO complex when compared to the same difference in an isolated heme control. This value agrees well with the experimental value of 2.0 kcal/mol. In qualitative agreement with experiments, the Fe-C-O bond is found to be bent (theta = 159.8 degrees) with a small tilt (theta = 6.2 degrees). The contributions made by each of the 29 residues--within the 9.0-A radius of the iron atom--to the free energy difference are separated into van der Waals and electrostatic contributions; the latter contributions are dominant. Aside from the proximal histidine and the heme group, the residues having the largest difference in free energy in mutating MbO2-->MbCO are His E7, Phe CD1, Phe CD4, Val E11, and Thr E10.     

Multi-drug resistance profile of PR20 HIV-1 protease is attributed to distorted conformational and drug binding landscape: molecular dynamics insights

The PR20 HIV-1 protease, a variant with 20 mutations, exhibits high levels of multi-drug resistance; however, to date, there has been no report detailing the impact of these 20 mutations on the conformational and drug binding landscape at a molecular level. In this report, we demonstrate the first account of a comprehensive study designed to elaborate on the impact of these mutations on the dynamic features as well as drug binding and resistance profile, using extensive molecular dynamics analyses. Comparative MD simulations for the wild-type and PR20 HIV proteases, starting from bound and unbound conformations in each case, were performed. Results showed that the apo conformation of the PR20 variant of the HIV protease displayed a tendency to remain in the open conformation for a longer period of time when compared to the wild type. This led to a phenomena in which the inhibitor seated at the active site of PR20 tends to diffuse away from the binding site leading to a significant change in inhibitor–protein association. Calculating the per-residue fluctuation (RMSF) and radius of gyration, further validated these findings. MM/GBSA showed that the occurrence of 20 mutations led to a drop in the calculated binding free energies (ΔG_bind) by ~25.17 kcal/mol and ~5 kcal/mol for p2-NC, a natural peptide substrate, and darunavir, respectively, when compared to wild type. Furthermore, the residue interaction network showed a diminished inter-residue hydrogen bond network and changes in inter-residue connections as a result of these mutations. The increased conformational flexibility in PR20 as a result of loss of intra- and inter-molecular hydrogen bond interactions and other prominent binding forces led to a loss of protease grip on ligand. It is interesting to note that the difference in conformational flexibility between PR20 and WT conformations was much higher in the case of substrate-bound conformation as compared to DRV. Thus, developing analogues of DRV by retaining its key pharmacophore features will be the way forward in the search for novel protease inhibitors against multi-drug resistant strains.     

Evaluating the suitability of RNA intervention mechanism exerted by some flavonoid molecules against dengue virus MTase RNA capping site: a molecular docking,molecular dynamics simulation,and binding free energy study

Abstract

Dengue virus (DENV) is one of the most dangerous mosquito-borne human pathogens known to the mankind. Currently, no vaccines or standard therapy is avaliable to treate DENV infection. This makes the drug development against DENV more significant and challenging. The MTase domain of DENV RNA RdRp NS5 is a promising drug target, because this domain hosts the RNA capping process of DENV RNA to escape from human immune system. In the present study, we have analysed the RNA intervention mechanism exerted by flavoniod molecules against NS5 MTase RNA capping site by using molecular docking, molecular dynamics simulation and the binding free energy calculations. The results from the docking analysis confirmed that the RNA intervention mecanism is exerted by the quercetagetin (QGN) molecule with all necessary intermolecular interactions and high binding affinity. Notably, QGN forms strong hydrogen bonding interactions with Asn18, Leu20 and Ser150 residues and π???π stacking interaction with Phe25 residue. The apo and QGN bound NS5 MTase and QGN-NS5 MTase complex were used for MD simulation. The results of MD simulation reveal that the RMSD and RMSF values of QGN-MTase complex have increased on comparing the apo protein due to the effect of ligand binding. The binding free energy calulation includes prediction of total binding free energy of ligand-protein complex and per-residue free energy decomposition. The QGN binding to NS5 MTase affects it’s native motion, this result is found from Principal component analysis.

Communicated by Ramaswamy H. Sarma     

RNA targeting by DNA binding drugs: Structural,conformational and energetic aspects of the binding of quinacrine and DAPI to A-form and H-form of poly(rC)·poly(rG)

A key step in the rational design of new RNA binding small molecules necessitates a complete elucidation of the molecular aspects of the binding of existing molecules to RNA structures. This work focuses towards the understanding of the interaction of a DNA intercalator, quinacrine and a minor groove binder 4′,6-diamidino-2-phenylindole (DAPI) with the right handed Watson–Crick base paired A-form and the left-handed Hoogsteen base paired H^L-form of poly(rC)·poly(rG) evaluated by multifaceted spectroscopic and viscometric techniques. The energetics of their interaction has also been elucidated by isothermal titration calorimetry. Results of this study converge to suggest that (i) quinacrine intercalates to both A-form and H^L-form of poly(rC)·poly(rG); (ii) DAPI shows both intercalative and groove-binding modes to the A-form of the RNA but binds by intercalative mode to the H^L-form. Isothermal calorimetric patterns of quinacrine binding to both the forms of RNA and of DAPI binding to the H^L-form are indicative of single binding while the binding of DAPI to the A-form reveals two kinds of binding. The binding of both the drugs to both conformations of RNA is exothermic; while the binding of quinacrine to both conformations and DAPI to the A-form (first site) is entropy driven, the binding of DAPI to the second site of A-form and H^L-conformation is enthalpy driven. Temperature dependence of the binding enthalpy revealed that the RNA–ligand interaction reactions are accompanied by small heat capacity changes that are nonetheless significant. We conclude that the binding affinity characteristics and energetics of interaction of these DNA binding molecules to the RNA conformations are significantly different and may serve as data for the development of effective structure selective RNA-based antiviral drugs.