Stabilization of proteins by ligand binding: application to drug screening and determination of unfolding energetics

The observed stability of a protein is altered when ligands bind, which results in a shift in the melting temperature (T(m)). Binding to the native state in the absence of binding to the denatured state will necessarily lead to an increase in the T(m), while binding to the unfolded state in the absence of native state binding will decrease the T(m) relative to that of the protein in the absence of ligand. These effects are required by the thermodynamics of reversible folding. However, the relationship between binding affinity and the magnitude of the observed temperature shift is not a simple correlation (i.e., a larger shift in T(m) does not necessarily mean tighter binding) and is complicated by interaction with the denatured state. Using exact simulations, the range of behavior for the dependence of the observed T(m) shift on the energetics of ligand binding is investigated here. Specifically, differential scanning calorimetry (DSC) curves are simulated for protein unfolding in the presence of ligands binding to both the native and denatured states. The results have implications for drug screening and the determination of heat capacity changes for protein unfolding.     

Melting behavior and ligand binding of DNA intramolecular secondary structures

We use a variety of biophysical techniques to determine thermodynamic profiles, including hydration, for the unfolding of DNA stem-loop motifs (hairpin, a three-way junction and a pseudoknot) and their interaction with netropsin and random cationic copolymers. The unfolding thermodynamic data show that their helix-coil transition takes place according to their melting domains or sequences of their stems. All hairpins adopted the B-like conformation and their loop(s) contribute with an immobilization of structural water. The thermodynamic data of netropsin binding to the ^5′-AAATT-^3′/TTTAA site of each hairpin show affinities of ~ 10^6- 7 M^− 1, 1:1 stoichiometries, exothermic enthalpies of − 7 to − 12 kcal mol^− 1 (− 22 kcal mol^− 1 for the secondary site of the three-way junction), and water releases. Their interaction with random cationic copolymers yielded higher affinities of ~ 10⁶ M^− 1 with the more hydrophobic hairpins. This information should improve our current picture of how sequence and loops control the stability and melting behavior of nucleic acid molecules.     

Parallel-stranded duplex DNA containing dA · dU base pairs

DNA oligonucleotides with dA and dU residues can form duplexes with trans d(A · U) base pairing and the sugar-phosphate backbone in a parallel-stranded orientation, as previously established for oligonucleotides with d(A · T) base pairs. The properties of such parallel-stranded DNA (ps-DNA) 25-mer duplexes have been characterized by absorption (uv), CD, ir, and fluorescence spectroscopy, as well as by nuclease sensitivity. Comparisons were made with duplex molecules containing (a) dT in both strands, (b) dU in one strand and dT in the second, and (c) the same base combinations in reference antiparallel-stranded (aps) structures. Thermodynamic analysis revealed that total replacement of deoxythymine by deoxyuridine was accompanied by destabilization of the ps-helix (reduction in T_m by −13°C in 2 mM MgGl₂, 10 mM Na-cacodylate). The U-containing ps-helix (U1 · U2) also melted 14°C lower than the corresponding aps-helix under the same ionic conditions; this difference was very close to that observed between ps and aps duplexes with d(A · T) base pairs. Force field minimized structures of the various ps and aps duplexes with either d(A · T) or d(A · U) base pairs ps/aps and dT/dU combinations are presented. The energy-minimized helical parameters did not differ significantly between the DNAs containing dT and dU. © 1996 John Wiley & Sons, Inc.     

DNA oligonucleotide duplexes containing intramolecular platinated cross-links: energetics,hydration, sequence,and ionic effects

The anticancer activity of cisplatin arises from its ability to bind covalently to DNA, forming primarily intrastrand cross-links to adjacent purine residues; the most common adducts involve d(GpG) (65%) and d(ApG) (25%) intrastrand cross-links. The incorporation of these platinum adducts in a B-DNA helix induces local distortions, causing bending and unwinding of the DNA. In this work, we used temperature-dependent UV spectroscopy to investigate the unfolding thermodynamics, and associated ionic effects, of two sets of DNA decamer duplexes containing either cis-[Pt(NH(3))(2)[d(GpG]] or cis-[Pt(NH(3))(2) [d(ApG]] cross-links, and their corresponding unmodified duplexes. The platinated duplexes are less stable and unfold with lower T(M)s (and Delta G degrees s) in enthalpy-driven reactions, which indicates a loss of favorable base-pair stacking interactions. The folding thermodynamics and hydration effects for the first set of decamers containing the d(GpG) cross-link was investigated by a combination of titration calorimetry, density, and ultrasound techniques. The hydration parameters showed an uptake of structural water by the platinated duplex and a release of electrostricted water by the control duplex. Relative to the unmodified duplex, the folding of the platinated duplex at 20 degrees C yielded a positive Delta Delta G degrees term [and positive Delta Delta H-Delta(T Delta S) compensation] and a negative differential volume change. The opposite signs of the Delta Delta G degrees and Delta Delta V terms confirmed its uptake of structural water. Further, solvent-accessible surface areas calculations for a similar pair of dodecamer duplexes indicated that the modified duplex has a 503 oeA(2) higher polar and nonpolar surface area that is exposed to the solvent. Therefore, the incorporation of a platinum adduct in duplex DNA disrupts favorable base-pair stacking interactions, yielding a greater exposure of aromatic bases to the solvent, which in turn immobilizes structural water. The overall results correlate nicely with the results reported in the available structural data of nuclear magnetic resonance solution studies.     

DNA deformation energetics and protein binding

The formation of protein-DNA complexes often involves deformation of the DNA double helix. We have calculated the energy necessary to produce this deformation in 71 crystallographically determined complexes, using internal coordinate energy optimization with the JUMNA program and a generalized Born continuum solvent treatment. An analysis of the data allows deformation energy to be interpreted in terms of both local and global structural changes. We find that, in the majority of complexes, roughly 60% of the deformation energy corresponds to backbone distortion. It is also found that large changes in stacking and pairing energies are often compensated for by other, longer range, stabilizing factors. Some deformations, such as base opening, can be large, but only-produce local energetic effects. In terms of backbone distortions, the angle alpha, most often involved in alphagamma transitions, makes the most significant energetic contribution. This type of transition is twice as costly as those involving beta, or coupled epsilonzeta changes. Sugar amplitude changes are also energetically significant, in contrast to changes in phase angles.     

Individual-site binding data and the energetics of protein–DNA interactions

Individual-site isotherms for the binding of bacteriophage λ repressor to the left and right λ operators have been determined [D. F. Senear, M. Brenowitz, M. A. Shea, and G. K. Ackers (1986) Biochemistry, Vol. 25, pp. 7344–7354.] using the DNAse protection technique [ footprinting; D. J. Galas and A. Schmitz (1978) Nucleic Acids Research, Vol. 5, pp. 3157–3170]. These extensive data have been interpreted with a quantitative model that emphasized cooperative interactions between adjacently bound ligands [occupied ? occupied interactions; G. K. Ackers, A. D. Johnson, and M. A. Shea (1982) Proceedings of the National Academy of Science, USA, Vol. 79, pp. 1129–1133]. Overlooked in this model are the effects of cooperative interactions between a site containing a bound ligand and its neighboring unoccupied site (occupied ? unoccupied interactions). This paper reinterprets the existing data with a model that considers occupied ? unoccupied as well as occupied ? occupied interactions. The results yield parameters that differ substantially from those already reported. A discussion on the advisability of ignoring occupied ? unoccupied interactions is included. © 1993 John Wiley & Sons, Inc.     

Analysis of the two-state behavior of the thermal unfolding serum retinol binding protein containing a single retinol ligand.

Through the use of CD and DSC, the thermal unfolding of holo serum retinol binding protein containing a single, tightly bound retinol ligand was studied at pH 7.4. The DSC endotherm of the holoprotein ([retinol]/[protein] = 1) was asymmetric about the transition temperature of 78 degrees C. Using changes in ellipticity at 230 nm, the thermal unfolding curve was also asymmetric about the inflection point centered near 78 degrees C. van't Hoff enthalpies were determined by three means and compared to the calorimetric enthalpy (delta Hcal) of 200 kcal/mol. A van't Hoff enthalpy of 190 kcal/mol was determined from the dependence of transition temperature on the concentration of the ligand-bound protein. This value agreed well with the van't Hoff enthalpies found from fits of the DSC (delta HvH = 184 kcal/mol) and spectroscopic (delta HvH = 181 kcal/mol) curves to a two-state thermodynamic model that included ligand dissociation (NR in equilibrium with U+R, where NR is the native holoprotein, U is the unfolded apoprotein, and R is retinol). Poor agreement was obtained with a two-state model that ignored ligand dissociation (N in equilibrium with U). Furthermore, the NR in equilibrium with U+R model accounted for the asymmetry in both CD and DSC transitions and yielded a much improved fit of the data over the N in equilibrium with U model. From these considerations and simulations on other equilibrium models, it is suggested that the NR in equilibrium with U+R model is the simplest model that describes the thermal unfolding of this ligand-bound protein.(ABSTRACT TRUNCATED AT 250 WORDS)     

Interaction of Hoechst 33258 and echinomycin with nucleosomal DNA fragments containing isolated ligand binding sites

We have examined the interaction of Hoechst 33258 and echinomycin with nucleosomal DNA fragments which contain isolated ligand binding sites. A 145 base pair fragment was prepared on the basis of the sequence of tyrT DNA, which contained no CpG or (A/T)(4) binding sites for these ligands. Isolated binding sites were introduced into this fragment at discrete locations where the minor groove is known to face toward or away from the protein core when reconstituted onto nucleosome core particles. The interaction of ligands with target sites on these nucleosomal DNA fragments was assessed by DNase I footprinting. We find that Hoechst 33258 can bind to single nucleosomal sites which face both toward and away from the protein core, without affecting the nucleosome structure. Hoechst binding is also observed on nucleosomal fragments which contain two or more drug binding sites, though in these cases the footprints are accompanied by the presence of new cleavage products in positions which suggest that the ligand has caused a proportion of the DNA molecules to adopt a new rotational positioning on the protein surface. Hoechst 33258 does not affect nucleosome reconstitution with any of these fragments. In contrast, the bifunctional intercalating antibiotic echinomycin is not able to bind to single nucleosomal CpG sites. Echinomycin footprints are observed on nucleosomal fragments containing two or more CpG sites, but there are no changes in the cleavage patterns in the remainder of the fragment. Echinomycin abolishes nucleosome reconstitution when included in the reconstitution mixture.     

Synthesis, characterization and DNA binding of ruthenium(II) complexes containing the atatp ligand

Acenaphtheno[1,2-b]-1,4,8,9-tetraazatriphenylene (atatp) and its complexes [Ru(L)2atatp](ClO4)2 x nH2O (L = 2,2'-bipyridine (bpy), n=2 (1); 1,10-phenanthroline (phen), n=2 (2); and 2,9-dimethyl-1,10-phenanthroline (dmp), n=1 (3)) have been synthesized and characterized by elemental analyses and 1H NMR. The spectral and electrochemical properties of these complexes are also examined. Complexes 1 and 2 display bright luminescence in acetonitrile but very weak luminescence in water solution. However, complex 3 is not luminescent in either solvent. The interaction of the complexes with calf thymus DNA (CT-DNA) has been studied by absorption, emission and viscosity measurements. The intrinsic binding constants of complexes 1 and 2 are 7.6 x 10(4) and 8.8 x 10(4) M(-1) respectively. The relatively low affinities of complexes 1 and 2 with DNA may arise from the atatp ligand, indicating that the size and shape of the intercalated ligand have a marked effect on the strength of interaction. Complexes 1 and 2 bind with CT-DNA in an intercalative mode but complex 3 in a non-intercalative one, showing that changing the ancillary ligand affects not only the binding magnitude, but also the binding mode of the interaction.     

DNA-bending on ligand binding; change of DNA persistence length.

This paper simulates the helix-characteristic changes of apparent DNA persistence length caused by randomly distributed helix bends as induced, e.g., by DNA-bound ligand molecules. The parameters varied are the constant angle gamma of helix bending and the size alpha of the DNA drug binding site, but also the degree of DNA-ligand binding cooperativity and the helix-unwinding angle. If the size of the binding site is comparable with the helix pitch, the influence of phasing between helix bends and helix screw upon the apparent persistence length is obvious. In the accompanying paper experimental data are analyzed in terms of this theoretical background.     

Structure and energetics of ligand binding to proteins: Escherichia coli dihydrofolate reductase-trimethoprim, a drug-receptor system

A study of the binding of the antibacterial agent trimethoprim to Escherichia coli dihydrofolate reductase was carried out using energy minimization techniques with both a full, all-atom valence force field and a united atom force field. Convergence criteria ensured that no significant structural or energetic changes would occur with further minimization. Root-mean-square (RMS) deviations of both minimized structures with the experimental structure were calculated for selected regions of the protein. In the active site, the all-atom minimized structure fit the experimental structure much better than did the united atom structure. To ascertain what constitutes a good fit, the RMS deviations between crystal structures of the same enzyme either from different species or in different crystal environments were compared. The differences between the active site of the all-atom minimized structure and the experimental structure are similar to differences observed between crystal structures of the same protein. Finally, the energetics of ligand binding were analyzed for the all-atom minimized coordinates. Strain energy induced in the ligand, the corresponding entropy loss due to shifts in harmonic frequencies, and the role of specific residues in ligand binding were examined. Water molecules, even those not in direct contact with the ligand, were found to have significant interaction energies with the ligand. Thus, the inclusion of at least one shell of waters may be vital for accurate simulations of enzyme complexes.     

Chiral diruthenium(II) complexes containing a modified 4,4-bipyridine bridging ligand: syntheses, characterization and DNA binding

A pair of novel chiral dimeric ruthenium(II) complexes [ΔΔ-, ΛΛ-Ru(bpy)₂(btpb)Ru(bpy)₂]⁴⁺ (1; btpb=2,2^′-bis(1,2,4-triazino[5,6-f]phenanthren-3-yl)-4,4^′-bipyridine) have been synthesized and characterized by electrospray mass spectra, ¹H NMR, UV-Vis and circular dichroism spectra. Binding behaviors of the complexes with calf thymus DNA have been investigated by absorption spectra and viscosity measurements. The electronic absorption spectrum of ΔΔ-1 at 505.5 nm exhibits hypochromism of about 8.4% and bathochromism of 2.5 nm; ΛΛ-1 at 500.0 nm exhibits hypochromism of about 9.1% and bathochromism of 4.5 nm, respectively. The experiments suggest that ΔΔ-1 and ΛΛ-1 may be bound to DNA by non-intercalating binder.     

Iron binding effects on the kinetic stability and unfolding energetics of a thermophilic phenylalanine hydroxylase from Chloroflexus aurantiacus

The effects of non-heme iron binding on the function, structure, and stability of a monomeric phenylalanine hydroxylase from the thermophile Chloroflexus aurantiacus (caPAH) were investigated. Comparative studies on holo (iron-bound) and apo (iron-depleted) caPAH indicated that iron(II) binding does not significantly affect the overall structure of the enzyme. Thermal denaturation studies performed using differential scanning calorimetry showed that the unfolding reaction was kinetically controlled and that holo-caPAH displayed a large increase in thermal stability (approximately 15 °C upshift in the T _m value) compared with the apoenzyme. Analysis using a simple irreversible two-state model also showed a higher kinetic stability for holo-caPAH at optimal growth temperature (denaturing approximately 8 times more slowly than the apo form at 55 °C). Experiments performed in the presence of urea in combination with structure–energetics calculations suggest that iron binding reduces the change in accessible surface area exposed in the unfolding transition state (from approximately 36% to approximately 5% of the total change in accessible surface area) and also the surface involved in water-unsatisfied broken internal contacts (solvation barriers). Additional comparative analyses using phenylalanine hydroxylase from mesophilic and psychrophilic organisms suggest that, in addition to its catalytic role, the non-heme iron serves to enhance the kinetic stability of phenylalanine hydroxylase at the optimal growth temperature of the organism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Thermal unfolding of the archaeal DNA and RNA binding protein Ssh10

The reversible thermal unfolding of the archaeal histone-like protein Ssh10b from the extremophile Sulfolobus shibatae was studied using differential scanning calorimetry and circular dichroism spectroscopy. Analytical ultracentrifugation and gel filtration showed that Ssh10b is a stable dimer in the pH range 2.5–7.0. Thermal denaturation data fit into a two-state unfolding model, suggesting that the Ssh10 dimer unfolds as a single cooperative unit with a maximal melting temperature of 99.9 °C and an enthalpy change of 134 kcal/mol at pH 7.0. The heat capacity change upon unfolding determined from linear fits of the temperature dependence of ΔH^cal is 2.55 kcal/(mol K). The low specific heat capacity change of 13 cal/(mol K residue) leads to a considerable flattening of the protein stability curve (ΔG (T)) and results in a maximal ΔG of only 9.5 kcal/mol at 320 K and a ΔG of only 6.0 kcal/mol at the optimal growth temperature of Sulfolobus.     

Synthesis,characterization and binding studies of chromium(III) complex containing an intercalating ligand with DNA

A chromium(III) complex [Cr(DPPZ)(2)Cl(2)](+), where DPPZ is a planar bidentate ligand with an extended aromatic system, has been found to bind strongly to CT DNA with an apparent binding constant of (1.8+/-0.5)x10(7) M(-1). The effects of [Cr(DPPZ)(2)Cl(2)](+) on the melting temperature and the viscosity of DNA clearly show that the chromium(III) complex interacts with DNA intercalatively. Competitive binding study shows that the enhancement in emission intensity of ethidium bromide (EthBr) in the presence of DNA was quenched by [Cr(DPPZ)(2)Cl(2)](+) indicating that the Cr(III) complex displaces EthBr from its binding site in DNA. The binding of this complex has been found to bring about B to Z conformational transition in CT DNA as well as poly(dG-dC).poly(dG-dC). Molecular modeling study also shows that binding energy of the complex with d(GC)(12) is much higher than Dickerson model and d(AT)(12). Modeling studies show that [Cr(DPPZ)(2)Cl(2)](+) brings about twist in the DNA base pairs as well as phosphate ester backbone resulting in conformational transition in DNA.     

The interrelationship between ligand binding and thermal unfolding of the folate binding protein. The role of self-association and pH

The present study utilized a combination of DLS (dynamic light scattering) and DSC (differential scanning calorimetry) to address thermostability of high-affinity folate binding protein (FBP), a transport protein and cellular receptor for the vitamin folate. At pH 7.4 (pI = 7–8) ligand binding increased concentration-dependent self-association of FBP into stable multimers of holo-FBP. DSC of 3.3 μM holo-FBP showed Tm (76 °C) and molar enthalpy (146 kcal M^− 1) values increasing to 78 °C and 163 kcal M^− 1 at 10 μM holo-FBP, while those of apo-FBP were 55 °C and 105 kcal M^− 1. Besides ligand binding, intermolecular forces involved in concentration-dependent multimerization thus contribute to the thermostability of holo-FBP. Hence, thermal unfolding and dissociation of holo-FBP multimers occur simultaneously consistent with a gradual decrease from octameric to monomeric holo-FBP (10 μM) in DLS after a step-wise rise in temperature to 78 °C ≈ Tm. Stable holo-FBP multimers may protect naturally occurring labile folates against decomposition or bacterial utilization. DSC established an interrelationship between diminished folate binding at pH 5, especially in NaCl-free buffers, and low thermostability. Positively charged apo-FBP was almost completely unfolded and aggregated at pH 5 (Tm 38 °C) and holo-FBP, albeit more thermostable, was labile with aggregation tendency. Addition of 0.15 M NaCl increased thermostability of apo-FBP drastically, and even more so that of holo-FBP. Electrostatic forces thus seem to contribute to a diminished thermostability at low pH. Fluorescence spectroscopy after irreversible thermal unfolding of FBP revealed a weak-affinity folate binding.