Quantitative determination of anthracycline antibiotic binding with DNA]

A new method of quantitative intravital assessment of anthracyclines accumulation and binding with DNA in alive cells has been developed. For this purpose DNA specific fluorescent dye Hoechst 33258 was used. Extent of fluorescence quenching of bound with DNA dye correlated with binding of daunomycin with cell DNA after mixing solutions of the drug and cells. It allowed us to determine quantity of drug molecules bound with DNA in the cells during the time.     

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.     

Electrostatic contributions to the binding free energy of the lambdacI repressor to DNA.

A model based on the nonlinear Poisson-Boltzmann (NLPB) equation is used to study the electrostatic contribution to the binding free energy of the lambdacI repressor to its operator DNA. In particular, we use the Poisson-Boltzmann model to calculate the pKa shift of individual ionizable amino acids upon binding. We find that three residues on each monomer, Glu34, Glu83, and the amino terminus, have significant changes in their pKa and titrate between pH 4 and 9. This information is then used to calculate the pH dependence of the binding free energy. We find that the calculated pH dependence of binding accurately reproduces the available experimental data over a range of physiological pH values. The NLPB equation is then used to develop an overall picture of the electrostatics of the lambdacI repressor-operator interaction. We find that long-range Coulombic forces associated with the highly charged nucleic acid provide a strong driving force for the interaction of the protein with the DNA. These favorable electrostatic interactions are opposed, however, by unfavorable changes in the solvation of both the protein and the DNA upon binding. Specifically, the formation of a protein-DNA complex removes both charged and polar groups at the binding interface from solvent while it displaces salt from around the nucleic acid. As a result, the electrostatic contribution to the lambdacI repressor-operator interaction opposes binding by approximately 73 kcal/mol at physiological salt concentrations and neutral pH. A variety of entropic terms also oppose binding. The major force driving the binding process appears to be release of interfacial water from the protein and DNA surfaces upon complexation and, possibly, enhanced packing interactions between the protein and DNA in the interface. When the various nonelectrostatic terms are described with simple models that have been applied previously to other binding processes, a general picture of protein/DNA association emerges in which binding is driven by the nonpolar interactions, whereas specificity results from electrostatic interactions that weaken binding but are necessary components of any protein/DNA complex.     

Molecular structure, stereochemistry and interactions of steffimycin B, and DNA binding anthracycline antibiotic

The crystal and molecular structure of anthracycline antibiotic steffimycin B(C29H32O13) has been determined by X-ray diffraction and the stereochemistry revealed. The orthorhombic crystals belong to space group P2(1)2(1)2(1), with the dimensions; a = 8.253 (2), b = 8.198 (2), c = 40.850 (8) A and Z = 4. Intensity data were collected for 2518 independent reflections. The structure was solved by direct methods and refined to an R value of 0.066 for 1410 reflections. The configuration in ring A is 7R,8S,9S. Ring A adopts half chair conformation, while the sugar ring has the regular chair conformation. The molecule most probably binds to double helical DNA through intercalation and hydrogen bonding.     

Kinetics of anthracycline antibiotic free radical formation and reductive glycosidase activity

Adriamycin free radical anion concentrations have been correlated with production of 7-deoxyadriamycin aglycone in a reaction catalyzed by NADPH-cytochrome c reductase. The free radical species is detected by electron spin resonance (ESR) spectroscopy and quantified by double integrations. The 7-deoxyaglycone product is isolated by thin-layer chromatography (TLC) and quantified by fluorometry. As the concentration of adriamycin increases, a concomitant increase in aglycone and free radical levels occurs. These results as well as those with inhibitors Vitamin K₃, Vitamin E, and 5,5-dimethyl-1-pyrroline-1-oxide (DMPO) point to an obligatory free radical intermediate in the metabolism of adriamycin. DMPO inhibits the reaction under aerobic conditions only, and shows no effect under anaerobiosis at the concentrations studied here. Vitamin E and aerobic DMPO act as free radical scavangers, while Vitamin K₃ competes for the reducing power of NADPH in the NADPH-cytochrome c reductase system.     

Comparison of the binding of anthracycline derivatives to purified DNA and to cell nuclei

Fluorescence method was used to study the interactions of anthracyclines with purified DNA and with cell nuclei at 37 degrees C, at pH ranging from 6.8 to 8. Four anthracyclines were used; adriamycin (ADR), 4'-o-tetrahydropyranyladriamycin (THP-ADR), daunorubicin (DNR) and aclacinomycin (ACM). The values of pKa of deprotonation of these four drugs in the pH range 6.5-8.5 are 8.4, 7.7, 8.4 and 7.0 for ADR, THP-ADR, DNR and ACM, respectively. The overall binding constants K* of these four drugs to purified DNA was determined at various pH values. The binding constants K0 and K+ of the respectively neutral form and once protonated form of the drugs to DNA were calculated. Using cell nuclei from K562 cells, the amount of drug intercalated (CN) within the nuclei of K562 cells and the amount of free drug (CE) in the solution were determined at various pH values: measuring at the same pH values, a linear correlation occurred between K* and CN/CE.     

Molecular dynamics simulations and binding free energy analysis of DNA minor groove complexes of curcumin

Curcumin is a natural phytochemical that exhibits a wide range of pharmacological properties, including antitumor and anticancer activities. The similarity in the shape of curcumin to DNA minor groove binding drugs is the motivation for exploring its binding affinity in the minor grooves of DNA sequences. Interactions of curcumin with DNA have not been extensively examined, while its pharmacological activities have been studied and documented in depth. Curcumin was docked with two DNA duplexes, d(GTATATAC)₂ and d(CGCGATATCGCG)₂, and molecular dynamics simulations of the complexes were performed in explicit solvent to determine the stability of the binding. In all systems, the curcumin is positioned in the minor groove in the A·T region, and was stably bound throughout the simulation, causing only minor modifications to the structural parameters of DNA. Water molecules were found to contribute to the stability of the binding of the ligand. Free energy analyses of the complexes were performed with MM-PBSA, and the binding affinities that were calculated are comparable to the values reported for other similar nucleic acid–ligand systems, indicating that curcumin is a suitable natural molecule for the development of minor groove binding drugs.     

Hydropathic interaction analyses of small organic activators binding to antithrombin

Recently we designed the first small organic ligands, sulfated flavanoids and flavonoids, that act as activators of antithrombin for accelerated inhibition of factor Xa, a key proteinase of the coagulation cascade [Gunnarsson and Desai, Bioorg. Med. Chem. Lett. (2003) 13:579]. To better understand the binding properties of these activators at a molecular level, we have utilized computerized hydropathic interaction (HINT) analyses of the sulfated molecules interacting in two plausible electropositive regions, the pentasaccharide- and extended heparin-binding sites, of antithrombin in its native and activated forms. HINT analyses indicate favorable multi-point interactions of the activators in both binding sites of the two forms of antithrombin. Yet, HINT predicts better interaction of most activators, except for (-)-catechin sulfate, with the activated form of antithrombin than with the native form supporting the observation in solution that these molecules function as activators of the inhibitor. Further, whereas (+)-catechin sulfate recognized the activated form of antithrombin better in both the pentasaccharide- and extended heparin- binding sites, the native form was better recognized by (-)-catechin sulfate, thus explaining its weaker binding and activation potential in solution. A reasonable linear correlation between the overall HINT score and the solution free energy of binding of the sulfated activators was evident. This investigation indicates that HINT is a useful tool in understanding interactions of antithrombin with small sulfated organic ligands at a molecular level, has some good predictive properties, and is likely to be useful for rational design purposes.     

Molecular dynamics simulations and free energy calculations of netropsin and distamycin binding to an AAAAA DNA binding site

Molecular dynamics simulations have been performed on netropsin in two different charge states and on distamycin binding to the minor groove of the DNA duplex d(CGCGAAAAACGCG)·d(CGCGTTTTTCGCG). The relative free energy of binding of the two non-covalently interacting ligands was calculated using the thermodynamic integration method and reflects the experimental result. From 2 ns simulations of the ligands free in solution and when bound to DNA, the mobility and the hydrogen-bonding patterns of the ligands were studied, as well as their hydration. It is shown that even though distamycin is less hydrated than netropsin, the loss of ligand–solvent interactions is very similar for both ligands. The relative mobilities of the ligands in their bound and free forms indicate a larger entropic penalty for distamycin when binding to the minor groove compared with netropsin, partially explaining the lower binding affinity of the distamycin molecule. The detailed structural and energetic insights obtained from the molecular dynamics simulations allow for a better understanding of the factors determining ligand–DNA binding.     

Use of Hoechst fluorescent probe 33258 for quantitative analysis of intracellular binding of anthracycline antibiotics to DNA]

The paper deals with development of a procedure for quantitative determination of the share of anthracycline antibiotics bound in cells directly to DNA. A DNA-specific Hoechst fluorescence dye 33258 was used for the purpose. The level of its quenching on DNA correlated with the quantity of the antibiotic bound to it. It was shown that the quenching of the Hoechst fluorescence dye bound to DNA was not due to the dye competition with the antibiotic for the site of bounding on DNA, as was suggested earlier. It was likely to be defined by reabsorption of the radiation by antibiotic molecules.     

Computation of entropy contribution to protein-ligand binding free energy

The entropy contribution ΔS to protein-ligand binding free energy is studied for nine protein-lipid complexes. The entropy effect from the loss of the translational/rotational degrees of freedom (ΔS _tr) is calculated using the ideal gas approach. The change in the vibrational entropy (ΔS _vib) is calculated using the effective quantum oscillator approach with frequencies derived from the coordinate covariance matrix, so the inharmonic effects are taken into account. The change in the entropy of solvation (ΔS _solv) is considered using the binomial cell model (developed by the authors) for the hydrophobic effect. The entropy contribution from loss of conformations that are available for the free ligand (ΔS _conf) is also estimated. It is revealed that the negative in view of binding term ΔS _tr is only partly compensated by increasing of ΔS _vib, so T(ΔS _tr + ΔS _vib + ΔS _conf) < 0 for all complexes under investigation, but taking into account ΔS _solv leads to significantly increased ΔS. For all complexes except biotin-streptavidin, the results are found to be in reasonable agreement with experimental data. Published in Russian in Biokhimiya, 2007, Vol. 72, No. 7, pp. 963–973.     

Contributions to the binding free energy of ligands to avidin and streptavidin

The free energy of binding of a ligand to a macromolecule is here formally decomposed into the (effective) energy of interaction, reorganization energy of the ligand and the macromolecule, conformational entropy change of the ligand and the macromolecule, and translational and rotational entropy loss of the ligand. Molecular dynamics simulations with implicit solvation are used to evaluate these contributions in the binding of biotin, biotin analogs, and two peptides to avidin and streptavidin. We find that the largest contribution opposing binding is the protein reorganization energy, which is calculated to be from 10 to 30 kcal/mol for the ligands considered here. The ligand reorganization energy is also significant for flexible ligands. The translational/rotational entropy is 4.5-6 kcal/mol at 1 M standard state and room temperature. The calculated binding free energies are in the correct range, but the large statistical uncertainty in the protein reorganization energy precludes precise predictions. For some complexes, the simulations show multiple binding modes, different from the one observed in the crystal structure. This finding is probably due to deficiencies in the force field but may also reflect considerable ligand flexibility.     

Study of the insulin dimerization: binding free energy calculations and per-residue free energy decomposition

A calculation of the binding free energy for the dimerization of insulin has been performed using the molecular mechanics-generalized Born surface area approach. The calculated absolute binding free energy is -11.9 kcal/mol, in approximate agreement with the experimental value of -7.2 kcal/mol. The results show that the dimerization is mainly due to nonpolar interactions. The role of the hydrogen bonds between the 2 monomers appears to give the direction of the interactions. A per-atom decomposition of the binding free energy has been performed to identify the residues contributing most to the self association free energy. Residues B24-B26 are found to make the largest favorable contributions to the dimerization. Other residues situated at the interface between the 2 monomers were found to make favorable but smaller contributions to the dimerization: Tyr B16, Val B12, and Pro B28, and to an even lesser extent, Gly B23. The energy decomposition on a per-residue basis is in agreement with experimental alanine scanning data. The results obtained from a single trajectory (i.e., the dimer trajectory is also used for the monomer analysis) and 2 trajectories (i.e., separate trajectories are used for the monomer and dimer) are similar.     

Thrombin inhibitors with novel P1 binding pocket functionality: free energy of binding analysis

The high incidence of thrombembolic diseases justifies the development of new antithrombotics. The search for a direct inhibitor has resulted in the synthesis of a considerable number of low molecular weight molecules that inhibit human α-thrombin potently. However, efforts to develop an orally active drug remain in progress as the most active inhibitors with a highly basic P1 moiety exhibit an unsatisfactory bioavailability profile. In our previous work we solved several X-ray structures of human α-thrombin in complexes with (1) novel bicyclic arginine mimetics attached to the glycylproline amide and pyridinone acetamide scaffold and (2) inhibitors with a novel aza scaffold and with charged or neutral P1 moieties. In the present contribution, we correlate the structures of the complex between these inhibitors and the protein with the calculated free energy of binding. The energy of solvation was calculated using the Poisson–Boltzmann approach. In particular, the requirements for successful recognition of an inhibitor at the protein’s active site pocket S1 are discussed. Figure We report here on free energy of binding analysis of thrombin inhibitors with novel aza scaffold and novel bicyclic arginine mimetics in S1 pocket of thrombin     

Modified solvent accessibility free energy prediction analysis of cyclic urea inhibitors binding to the HIV-1 protease

One of the most successful drug targets against AIDS in the last decade has been the HIV-1 protease (HIV-1 PR), an enzyme that processes the polyprotein gene products into active replicative viral proteins. In our quest for a wide-ranging, binding free energy function we have extended the solvent accessibility free energy predictor (SAFE_p) method, recently developed for peptidic HIV-1 PR inhibitors, to the study of the binding of cyclic urea (CU) HIV-1 PR inhibitors. Our results show that there is a need for a specific term depicting polar contacts to be added to the original SAFE_p analytical expression, an outcome not seen in our studies of HIV-1 PR peptidic inhibitors. Nevertheless, despite the higher profile of the electrostatic interactions in the binding of the CU inhibitors, our analysis indicates that CU inhibitor binding is still driven by the hydrophobic entropic contribution, as much as for the peptidic inhibitors.     

Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes

Absolute binding free energy calculations and free energy decompositions are presented for the protein-protein complexes H-Ras/C-Raf1 and H-Ras/RalGDS. Ras is a central switch in the regulation of cell proliferation and differentiation. In our study, we investigate the capability of the molecular mechanics (MM)-generalized Born surface area (GBSA) approach to estimate absolute binding free energies for the protein-protein complexes. Averaging gas-phase energies, solvation free energies, and entropic contributions over snapshots extracted from trajectories of the unbound proteins and the complexes, calculated binding free energies (Ras-Raf: -15.0(+/-6.3)kcal mol(-1); Ras-RalGDS: -19.5(+/-5.9)kcal mol(-1)) are in fair agreement with experimentally determined values (-9.6 kcal mol(-1); -8.4 kcal mol(-1)), if appropriate ionic strength is taken into account. Structural determinants of the binding affinity of Ras-Raf and Ras-RalGDS are identified by means of free energy decomposition. For the first time, computationally inexpensive generalized Born (GB) calculations are applied in this context to partition solvation free energies along with gas-phase energies between residues of both binding partners. For selected residues, in addition, entropic contributions are estimated by classical statistical mechanics. Comparison of the decomposition results with experimentally determined binding free energy differences for alanine mutants of interface residues yielded correlations with r(2)=0.55 and 0.46 for Ras-Raf and Ras-RalGDS, respectively. Extension of the decomposition reveals residues as far apart as 25A from the binding epitope that can contribute significantly to binding free energy. These "hotspots" are found to show large atomic fluctuations in the unbound proteins, indicating that they reside in structurally less stable regions. Furthermore, hotspot residues experience a significantly larger-than-average decrease in local fluctuations upon complex formation. Finally, by calculating a pair-wise decomposition of interactions, interaction pathways originating in the binding epitope of Raf are found that protrude through the protein structure towards the loop L1. This explains the finding of a conformational change in this region upon complex formation with Ras, and it may trigger a larger structural change in Raf, which is considered to be necessary for activation of the effector by Ras.     

Structure-activity relationships in engineered proteins: analysis of use of binding energy by linear free energy relationships

The activity of mutant enzymes can be analyzed quantitatively by structure-activity relationships in a manner analogous to Br?nsted or Hammett plots for simple organic reactions. The slopes of such plots, the beta values, indicate for the enzymatic reactions the fraction of the overall binding energy used in stabilizing particular complexes. In particular, information can be derived about the interactions between the enzyme and the transition state. The activities of many mutant tyrosyl-tRNA synthetases fit well simple linear free energy relationships. The formation of enzyme-bound tyrosyl adenylate (E.Tyr-AMP) from enzyme-bound tyrosine and ATP (E.Tyr-ATP) results in an increase in binding energy between the enzyme and the side chain of tyrosine and the ribose ring of ATP. Linear free energy plots of enzymes mutated in these positions give the fraction of the binding energy change that occurs on formation of the transition state for the chemical reaction and the various complexes. It is shown that groups that specifically stabilize the transition state of the reaction are characterized by beta values much greater than 1. This is found for residues that bind the gamma-phosphate of ATP (Thr-40 and His-45) and have previously been postulated to be involved in transition-state stabilization.(ABSTRACT TRUNCATED AT 250 WORDS)