Thermodynamic characterization of specific
interactions between the human Lon protease
and G-quartet DNA
Si-Han Chen
1,2, Carolyn K. Suzuki
3and Shih-Hsiung Wu
1,2,*
1
Institute of Biological Chemistry, Academia Sinica, Taipei 115,2Institute of Biochemical Science, National Taiwan University, Taipei 106, Taiwan and3Department of Biochemistry and Molecular Biology, New Jersey Medical School, University of Medicine and Dentistry of New Jersey, Newark, New Jersey 07103, USA
Received September 5, 2007; Revised November 24, 2007; Accepted December 7, 2007
ABSTRACT
Lon is an ATP-powered protease that binds DNA. However, the function of DNA binding by Lon remains elusive. Studies suggest that human Lon (hLon) binds preferentially to a G-rich single-stranded DNA (ssDNA) sequence overlapping the light strand promoter of mitochondrial DNA. This sequence is contained within a 24-base oligonu-cleotide referred to as LSPas. Here, we use biochemical and biophysical approaches to eluci-date the structural properties of ssDNAs bound by hLon, as well as the thermodynamics of DNA binding by hLon. Electrophoretic mobility shift assay and circular dichroism show that ssDNAs with a propensity for forming parallel G-quartets are specifically bound by hLon. Isothermal titration calorimetry demonstrates that hLon binding to LSPas is primarily driven by enthalpy change associated with a significant reduction in heat capacity. Differential scanning calorimetry pinpoints an excess heat capacity upon hLon binding to LSPas. By contrast, hLon binding to an 8-base G-rich core sequence is entropically driven with a relatively negligible change in heat capacity. A considerable enhancement of thermal stability accompanies hLon binding to LSPas as compared to the G-rich core. Taken together, these data support the notion that hLon binds G-quartets through rigid-body binding and that binding to LSPas is coupled with structural adaptation.
INTRODUCTION
The ATP-dependent Lon protease is a member of the AAA+ family of proteins (ATPases associated with various cellular activities) (1), which is present in archaea,
prokaryotes and eukaryotic mitochondria, peroxisomes and plastids (2–9). Lon has multiple cellular functions such as degrading abnormal polypeptides (2,4,10,11), as well as certain regulatory proteins (12–15) and metabolic enzymes (16–19), acting as a chaperone (20,21) and binding to nucleic acids (22–28). In addition, Lon may be important in maintaining genome integrity and expression, either by selectively degrading proteins involved in mtDNA metabolism, or by directly binding to DNA and/or RNA.
Recent work demonstrates that hLon binds to mtDNA in living cells and interacts preferentially with the control region for transcription initiation and mtDNA replication (28). It is possible that the sequence-specific binding of Lon to the mitochondrial genome may function to localize the protease at sites where it degrades protein substrates involved in mtDNA integrity. In bacteria, some of Lon’s protein substrates are involved in genome metabolism, which manifests the role of Lon in maintaining genome stability. For example, in Caulobacter crescentus, Lon mediates the ATP-dependent turnover of the CcrM DNA methyltransferase, thereby regulating DNA replication (15). In Escherichia coli, Lon selectively degrades the transcription regulator SoxS only when this factor is not bound to DNA or to RNA polymerase (29). In yeast, results demonstrate that the TCA enzyme mitochondrial aconitase is present in mitochondrial nucleoids and functions in maintaining mtDNA, independent of its enzymatic activity (30–32). In light of observations that mildly oxidized mammalian mitochondrial aconitase is subject to Lon-mediated proteolysis, it is possible that Lon participates in nucleoid remodeling by degrading mtDNA-bound aconitase in response to metabolic changes and environmental stress (19,33). However, no data as yet support the localization of mammalian aconitase to mito-chondrial nucleoids (34). The role of Lon in mitomito-chondrial homeostasis is indispensable in normal fibroblasts; depletion of the protease using an antisense morpholino oligonucleotide leads to apoptosis and collapse of
*To whom correspondence should be addressed. Tel: +886 2 27855696 7101; Fax: +886 2 27889759; Email: shwu@gate.sinica.edu.tw
ß 2008 The Author(s)
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/ by-nc/2.0/uk/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
mitochondrial structure (35). Nevertheless, understanding Lon’s physiological importance is enmeshed by its multi-ple functions as a protease and a DNA-binding protein. Thus, it remains unclear the extent to which mtDNA binding contributes to the essential role of Lon in mitochondrial biology.
Although the function and mechanism of Lon-mediated proteolysis is understood in detail, the physiological relevance and structural dynamics of nucleic acid binding by Lon remain unclear. In vitro studies demonstrate that both bacterial and mammalian Lon show a propensity for binding GT-rich oligonucleotides. Bacterial Lon binds to double-stranded DNA, whereas mouse and human Lon bind to ssDNA (22–25). Observations that the DNA-binding and enzymatic activities of bacterial and mam-malian Lon are coordinately regulated and that hLon binds to mtDNA in cells, support the notion that the association of Lon to DNA is physiologically relevant (25,36). In vitro experiments have shown that hLon specifically binds to a ssDNA sequence that is referred to as LSPas, which overlaps the light strand promoter anti-sense or G-rich strand mtDNA (23). Additional in vitro data also demonstrate that hLon binds G-rich oligonucleotides corresponding to regions present throughout the heavy strand of mtDNA but shows the highest affinity for LSPas (25,28). An hLon-binding consensus has been suggested by results obtained using a combination of in vitro and cell-based assays along with bioinformatics (28). Furthermore, in cultured mam-malian cells, hLon has been shown to associate with sites distributed primarily within one-half of the genome and preferentially with the control region for mtDNA replica-tion and transcripreplica-tion encompassing LSPas (28). In oxidatively stressed cells, hLon binding to mtDNA is reduced. However, it is unknown whether oxidative stress directly divests hLon of its DNA-binding property, or whether hLon binding is blocked by the oxidative modification of mtDNA-binding sites per se, or by other proteins recruited to repair these damaged sites.
Previous studies on the binding of hLon to ssDNA provided limited insight into the specificity, stability and dynamics of this interaction. It has been suggested that hLon may recognize G-rich sequences that form guanine tetraplexes, or quadruplexes, which are four-stranded structures based on the hydrogen-bonded guanine tetrad, or G-quartet (25). G-quartet DNA is a high-ordered structure that contains polyanionic groups, and thus might interact with basic proteins in a nonspecific manner. In this study, we demonstrate that mtDNA sequences bound by hLon form parallel G-quadruplexes. Results show that hLon selectively binds to DNA sequences that have a propensity for forming G-quartets. However, the presence of G-quartets does not confer specificity for hLon association. Human Lon complexed with LSPas exhibits a large negative heat capacity change, which is a thermodynamic signature of sequence-specific interactions. By contrast, the interaction between hLon and other G-rich sequences examined leads to little or no significant change in heat capacity, and is primarily entropy driven. Taken together, these energetic data show that the affinity of hLon for ssDNA is determined
by the propensity of the ssDNA sequence to form G-quartets, whereas the sequence specificity of hLon binding hinges on the extent of structural adaptation and release of interfacial solvent molecules between the protease and its DNA ligand.
MATERIALS AND METHODS
Expression and purification of human Lon protease
The expression plasmid of human lon was transformed into Rosetta 2 competent E. coli cells (Novagen), and recombinant hLon was expressed and purified as pre-viously described (25) with some modifications. Briefly, transformed cells were grown at 378C in medium contain-ing 3% tryptone, 2% yeast extract, 1% MOPS (pH 7.2), 100 mg/ml ampicillin and 34 mg/ml chloramphenicol. At mid-log phase (OD6000.6–0.8), hLon production was
induced by the addition of 1 mM IPTG for 1.5 h. The cells were harvested and resuspended in lysis buffer [50 mM Tris–HCl, pH 8.0, 300 mM NaCl, 10 mM imidazole, 10 mM b-Me, 20% glycerol (v/v) and 1% Triton X-100]. The lysate was sonicated on ice, and centrifuged at 48 000g at 48C for 1 h. The supernatant was filtered (0.22 mm) and then applied to a HiTrap HP column (Amersham) equilibrated with Buffer A [20 mM Tris–HCl, pH 8.0, 300 mM NaCl, 10 mM MgCl2, 10% glycerol (v/v) and
40 mM imidazole]. After serial washes with Buffer A and Wash Buffer B (the same as A but containing 80 mM imidazole), hLon was eluted in Elution Buffer (the same as Buffer A but containing 400 mM imidazole). Purified hLon was concentrated with Amicon-50000 (Millipore), and dialyzed extensively against a buffer containing 20 mM HEPES–NaOH (pH 7.5), 50 mM NaCl, 5 mM MgCl2, 0.1 mM DTT and 20% glycerol (v/v) unless
specified otherwise. The purification steps described above were conducted at 48C to avoid autodigestion of hLon. To exclude protein aggregates, the sample was centrifuged and filtered through a 0.22 mm-syringe filter. The protein concentration was determined by the Bradford method using bovine serum albumin (BSA) as a standard (Bio-Rad) (37).
Oligonucleotides
The oligonucleotides used in this study were purchased from ScinoPharm Biotech (Taiwan), and the sequences are listed in Table 1. DNA concentrations were determined spectrophotometrically using the extinction coefficients at 260 nm calculated from nearest-neighbor approxima-tion (38).
Polyacrylamide gel electrophoresis
The DNA oligonucleotides were analyzed by native and denaturing gel electrophoresis on a SE600 with a water-circulating system (Amersham, NJ, USA). The samples were prepared at a concentration of 50 mM by diluting the stock solutions with TB buffer [89 mM Tris–borate, pH 8.0, 100 mM NaCl and 5 mM MgCl2] in the presence
(denaturing) or absence (native) of 7 M urea. Prior to loading, the samples were heated at 958C for 10 min and
gently cooled down to room temperature. DNA electro-phoresis was performed on 20% polyacrylamide gels [1 : 19 (w/w) bisacrylamide: acrylamide] containing either native or denaturing TB buffer and run at constant power of 25 W. The nucleic acids were detected by UV shadowing by placing the gel on fluorescein-coated silica (Merck) (39). Electrophoresis mobility shift assay (EMSA)
50-Fluorescein-labeled LSPas or LSPasGA (1 mM) was
incubated with hLon (6.4 mM) in a 50 ml reaction contain-ing 20 mM HEPES–NaOH (pH 7.5), 50 mM NaCl, 5 mM MgCl2, 0.1 mM DTT and poly (dI-dC)(dI-dC)
(50 ng/reaction) (Amersham), in the presence and absence of 50-fold-excess unlabeled DNA competitor. After incubation for 20 min at 258C or specified otherwise, the reaction mixtures were resolved on 6% native PAGE containing TGE buffer [50 mM Tris, pH 8.5, 0.38 M glycine, 2 mM EDTA and 0.05 mM b-ME] run at a constant power of 8 W. The electrophoresis tempera-ture was controlled at 208C with water circulation. The gels were scanned with Typhoon 9200 Variable Mode Imager (Amersham) at blue-excited fluorescence mode (488 nm).
Circular dichroism (CD) spectroscopy
The CD experiments were carried out on a J-715 spectropolarimeter (Jasco, MD, USA) equipped with a
temperature controller. A water-jacketed cuvette of 0.1-cm path length was used. The cell-holder chamber was constantly purged with a stream of dry nitrogen gas. For CD measurements, the DNA samples (50 mM) were heated to 958C and cooled to room temperature prior to analysis. The data were recorded between 210 and 320 nm at a scanning rate of 20 nm/min. After background correction, the data are presented as molar ellipticity []M(degdmol1cm2) using
½ M¼ m8ð ÞMr
10cl 1
Mr is the sample molecular weight (g/mole), c is the
concentration in g/ml, and l is the path length in cm. The far-UV CD spectra of hLon (2.5 mM) were measured in 20-mM ammonium acetate (pH 7.5) and 20% glycerol. The data were collected between 200 and 250 nm at a rate of 20 nm/min. The CD measurements of hLon in the presence of different DNA oligonucleotides were carried out under identical condi-tions. After background subtraction (the buffer with and without the corresponding amount of DNA), the CD intensity was presented as mean residue ellipticity []MRE in degdmol1cm2, which was obtained by
averaging []Mout the number of amino acid residues in
hLon.
Table 1. The DNA sequences used in this study
Sequence Location of nucleotidesa Strand
LSPas 50-AATAATGTGTTAGTTGGGGGGTGA-30 455–432 Noncoding
LSPas18 50-GTGTTAGTTGGGGGGTGA-30 449–432 Noncoding
TG6T 50-TGGGGGGT-30 441–434 Noncoding
LSPs 50-TCACCCCCCAACTAACACATTATT-30 455–432 Coding
LSPs18 50-TCACCCCCCAACTAACAC-30 449–432 Coding
AC6A 50-ACCCCCCA-30 441–434 Coding
LS1as 50-GGCGTAGGTTTGGTCTAGGGT-30 7114–7134 Noncoding
LS2as 50-GTAGAGGGGGTGCTATAGGGT-30 8262–8282 Noncoding
LS3as 50-GTATGGGGGTAATTATGGTGG-30 8395–8415 Noncoding
LS4as 50-GGAGGGGGGTTGTTAGGGGGT-30 10 936–10 956 Noncoding
LS5as 50-CGAGGGTGGTAAGGATGGGGG-30 12 387–12 407 Noncoding
LS6as 50-GGGGAGGGGTGTTTAAGGGGGT-30 15 526–15 546 Noncoding
LSPasGA 50-AATAATGTGTTAGTTGGAAGGTGA-30 Mutant
18-mer ssDNA 50-CTCTTCTCCTCTCTTTCC-30 Control
Isothermal titration calorimetry (ITC)
The measurements of binding thermodynamics between hLon and DNA were performed isothermally on a VP-ITC MicroCalorimeter (MicroCal, MA, USA). To minimize the artifacts resulting from small differences in buffer composition, DNA (titrant) and hLon (titrand) were prepared with the identical lot of buffer consisting of 50 mM HEPES–NaOH (pH 7.5), 100 mM NaCl, 5 mM MgCl2, 10% glycerol (v/v) and 0.01% sodium azide. The
samples were filtered (0.22 mm) and degassed immediately before use. The calorimetric cell was loaded with hLon at 30 mM in a volume of 1.43 ml; LSPas or TG6T at a
concentration of 350 mM in the syringe was programmed to titrate into the reaction cell (12 ml/injection) with 6-min spacing intervals. The stirring rate of the injector was 300 r.p.m. The heat responses caused by the dilution of hLon were found negligible by titrating the buffer against hLon and vice versa. The ITC measurements of DNA against buffer were also carried out separately. After corrected for the dilution effects of DNA, the binding isotherms of hLon and DNA were normalized and fitted with Equations (2) and (3) by nonlinear least square minimiza-tion procedures (40).
qi¼n½titrandtotVcellH < 2
qi is the corrected heat response in the ith step of
titration, and n stands for the number of identical and independent binding sites for titrant to titrand. [titrand]tot
is the total concentration of titrand in the sample cell at a constant volume of Vcell. H represents the apparent
enthalpy change obtained from single experiment. The < in Equation (2) is given by
Y2i Yi 1 þ 1 nKa½titrandtot þ ½titrantitot n½titrandtot
þn½titrantitot½titrandtot¼0 3
where Yi is the degree of saturation defined by
Yi=[titranti]bound/[titrand]tot[titranti]totmeans the
accu-mulative concentration of titrant in the sample cell until injection i.
The free energy and entropy changes of binding were obtained using the relationship G = RTlnKa=
HTS. In addition, the heat capacity changes of association between hLon and DNA were given by Cp= @(H)/@T, where H is the binding enthalpy,
and T is the corresponding experiment temperature on ITC.
Differential scanning calorimetry (DSC)
DSC measures heat capacity (Cp) as a function of
temperature. The experiments were carried out on a VP-DSC MicroCalorimeter (MicroCal, MA, USA). The sample preparation steps were similar to those in ITC, including the buffer composition. The sample cell was loaded with a solution containing hLon at 25 mM in the presence or absence of DNA (37.5 mM). The reference cell was loaded with the buffer. The heat capacity was scanned from 108C to 1108C at a rate of 508C/h.
Scanning of the buffer against buffer as well as the DNA samples against buffer were performed in separate experiments. The heat capacity of the buffer effect was corrected from each thermogram according to the follow-ing equations (41). @Cp @T P ¼@Cp,Pð ÞT @T VPð ÞT VSð ÞT @Cp,Sð ÞT @T 4 @Cp @T L ¼@Cp,Lð ÞT @T VLð ÞT VSð ÞT @Cp,Sð ÞT @T 5 @Cp @T PL ¼@Cp,PLð ÞT @T VPLð ÞT VSð ÞT @Cp,Sð ÞT @T 6 @Cp/@T represents the slope of heat capacity in the
pretransition region of sample. Denotations of S, P, L and PL stand for the solvent, protein, DNA and mixture of protein/DNA, respectively. The V-bars are the partial-specific volumes (ml/mg) of each component.
The effect of DNA binding on the heat capacity of hLon may be substantial, which is reflected on the pretransition Cprather than the unfolding Cpbecause their interactions
are mostly noncovalent. The pretransition Cp can be
described by the following equation:
Cp¼Cp,0þ
@ PNi¼1XiHX,i
@T 7
where Cpis the partial heat capacity, and Cp,0is the partial
heat capacity of the native state; Xi represents the
population of solute in the ith state, and HX,i is the
enthalpy difference between the ith state and the native state (42). The heat capacities of proteins in the native state examined by DSC have always been found to increase linearly with temperature; they comprise not only an intrinsic variable of structural flexibility but also an instrumental artifact (42,43). To better assess the effect of DNA binding on hLon, we evaluated the difference in (@Cp/@T) between hLon and the DNA-bound hLon
[Equation (8)]. @Cp,PL @T ¼ @Cp,P,0 @T þ @Cp,L,0 @T þ @2 PN i¼1PLiHPL,i @T2 þ @2 PM j¼1PjHP,j @T2 þ @2 PQ k¼1LkHL,k @T2 8
Both @Cp,P,0/@T and @Cp,L,0/@T are nearly zero because
they represent the temperature dependence of the intrinsic heat capacities in the native state. Thus, the difference in @Cp,PL/@T is primarily accountable for the structural
flexibility of the protein, DNA and the complex.
RESULTS
The candidate-binding site of human Lon forms a parallel G-quartet structure in vitro
Recent work suggests that human Lon binds to mtDNA in living cells and preferentially associates with the con-trol region (CR) of the mitochondrial genome (28).
The putative mtDNA-binding target of hLon within CR overlaps the light strand promoter (LSP) and corresponds to the previously characterized oligonucleotide sequence referred to as LSPas. LSPas is a 24-base sequence located
on the heavy (or antisense) strand of mtDNA (50-AATA
ATGTGTTAGTTGGGGGGTGA-30). Notably, the
mtDNA heavy strand has a higher content of guanine and thymine residues, in contrast to the complementary light strand (44). LSPas contains six contiguous guanine bases, rendering a high-order structure demonstrated by electrophoretic analysis using both native and denaturing gels (Figure 1). Two forms of LSPas were observed (Figure 1A, lane 2): a monomeric form that migrated slightly slower than the complementary LSPs CA-contain-ing oligonucleotide (lane 1) and a high-order structure that migrated slower than the double-stranded LSP (compare lanes 2 and 3). Similar migration differences were also observed for an 18-mer counterpart (LSPas18) that lacks six bases from the 50end of LSPas but retains
the G-rich core sequence. The slower migrating forms of both LSPas and LSPas18 remained unchanged on denaturing polyacrylamide gels containing urea (Figure 1B, lanes 2 and 5, respectively), suggesting the presence of a structural motif, which was later identified by CD (Figure 2). Similarly, an 8-mer TG6T migrated
much slower than its complementary AC6A (Figure S1,
lanes 5 and 7, respectively). TG6T was unaffected by heat
treatment at 1208C for 2 h, whereas the slower migrating form of LSPas was disrupted after heat treatment (Figure S1, compare lane 5 with 6, and lane 1 with 2). The susceptibility to heat denaturation may be attrib-uted to the length of the ssDNA oligonucleotide, in addition to the number of bases flanking the G-rich core sequence.
Figure 1. Electrophoretic mobility of G-rich and control oligonu-cleotides. Oligonucleotides were analyzed under (A) native and (B) denaturing conditions using 20% polyacrylamide gels. G-rich oligonucleotides—LSPas, LSPas18 and TG6T; control oligonucleo-tides—LSPs, dsLSP, LSPs18, dsLSP18, an 18-mer ssDNA and 18-mer dsDNA. Arrows indicate the higher-ordered structures of LSPas and LSPas18. A 20-bp ladder DNA marker was loaded but not shown here under UV shadowing (stained with SYBR Gold).
Figure 2. Formation of G-quartets in LSPas and TG6T demonstrated by circular dichroism spectroscopy. CD spectra of LSPas (A and B) and TG6T (C and D) were monitored from 108C to 908C in the presence or absence of urea. The buffer used in A and C contains 10 mM sodium cacodylate (pH 7.5), 100 mM NaCl, 5 mM MgCl2and 0.1 mM EDTA. The buffer used in B and D consists of 10 mM sodium cacodylate (pH 7.5) and 1 M urea.
CD was employed to identify the possible structural differences in these G-rich ssDNA oligonucleotides. The CD spectra of LSPas in a native buffer containing 100 mM NaCl exhibited a strong positive peak at 265 nm and a negative one at 245 nm, reflecting the presence of a parallel G-quartet (Figure 2A). Typically, parallel G-quartets are characterized by a positive CD band at 260 nm and a negative one at 240 nm, whereas antiparallel G-quartets have a positive band at 295 nm and a negative one at 260 nm (45,46). The structure of LSPas remained mostly intact between 108C and 308C, the range of which was thus selected for monitoring its association with hLon in the ITC and DSC experiments (discussed later). LSPas became relatively unstructured at above 508C (Figure 2A). The CD profile of LSPas with respect to temperature in a Na+-containing denaturing buffer (the native buffer with 1-M urea; spectra not shown) was compared to that in the native buffer. However, the G-quartet signature of LSPas was dramatically weakened in the denaturing buffer without NaCl; the ellipticity not only decreased but also red shifted with the temperature increment (Figure 2B). This result suggested the importance of monovalent cations in stabilizing G-quartet LSPas DNA. It has been reported that K+, Na+, NH4+ and several other
monovalent cations stack between adjacent G-tetrads in G-quartet DNA to coordinate the eight carbonyl groups in every two tetrads, thereby making the G-quartet rigid (47–49). In particular, Na+and NH4+have been shown
to neutralize the peripheral phosphate groups of G-quartets (50). In contrast to LSPas, TG6T exhibited a
comparatively rigid G-quartet as revealed in the CD results; it was only mildly denatured by high temperature or by denaturing buffer (Figure 2, compare C with D; Figure S1, compare lane 5 with 6). TG6T showed a strong
positive peak at 262 nm and a negative one at 242 nm (Figure 2C). The sequences of both LSPas and TG6T
contain no contiguous guanines other than a 6-nt G cluster; the contexts hardly allow the formation of hairpins or intramolecular G-quartets. Based on the electrophoretic and CD data, LSPas and TG6T form
intermolecular parallel G-quartet structures in Na+ solutions. LSPas was more vulnerable than TG6T to
thermal and chemical denaturation because of the flanking sequence effect: the repulsion caused by the unpaired bases hinders the formation of LSPas quadruplex (see the schematic structures in Figure 2).
Human Lon selectively binds to G-quartets
The number of consecutive guanine bases in an inter-molecular G-quartet is the main determinant of its stability. When we replaced the third and fourth guanines within LSPas with adenines, the resulting ssDNA oligo-nucleotide referred to as LSPasGA did not form a G-quartet (Figure S1, lane 4). The EMSA was used to examine the binding of hLon to fluorescein-conjugated LSPas or LSPasGA (Figure 3A). The results showed that hLon selectively interacted with the G-quartet form of LSPas but not with the monomer. When 7-deaza-guanines (7-deaza-Gs) were substituted for the 6-nt G cluster of LSPas, the sequence failed to form a G-quartet because
the C-7 cannot be an electron donor for Hoogsteen hydrogen bonding. The 7-deaza-G modified LSPas was not bound by hLon (data not shown). The binding of hLon to fluoresceinated LSPas was competed by the addition of 50-fold excess of unlabeled LSPas but not LSPasGA. Moreover, hLon interacted only weakly with fluorescein-conjugated LSPasGA, which lacked G-quartet-forming ability.
It is noteworthy that the heavy strand of human mtDNA exhibits a high frequency of 4–12 contiguous guanine residues every 100 bases. The binding of hLon to G-rich sequences in vivo is supported by a combination of cell-based, in vitro and bioinformatics analyses (25,28). In vitro experiments examining hLon binding to ssDNA oligonucleotides containing several of these G-rich sites present within the mitochondrial genome, show that hLon binds with higher affinity to LSPas than to the other G-rich sequences (25,28). However, previous studies did not directly determine the extent to which any of these G-rich ssDNAs form G-quartets. For this purpose, we analyzed the G-quartet-forming capabilities of these G-rich oligonucleotides by native gel electrophoresis (Figure S2). Some of these G-rich ssDNAs formed extensive high order structures. Although spectroscopic analysis was not performed because of the structural complexities, the high order forms likely represented G-quartets resulting from 2–3 discrete G-clusters in their sequence contexts. Results showed that the sequences capable of forming high order structures were better able to compete with LSPas for hLon binding (Figure S3). Dose-dependent binding of LSPas as well as TG6T to
hLon is shown in Figure 3B.
The binding thermodynamics of hLon to DNA measured by ITC
Knowing that hLon binds to G-quartet DNA, we used ITC to investigate the energetics of hLon complexed with LSPas and TG6T. The binding reaction of hLon to LSPas
is exothermic over the temperature range inspected (Figure 4A). The amount of heat release was integrated and fitted nonlinearly (Figure 4B). The enthalpy of association of hLon with LSPas decreased from 1.48 to 13.72 kcal/mol over a temperature range from 108C to 308C (Table 2). In contrast to the downward trend in the enthalpy change, the entropy penalty (TS) increased from 5.49 to 6.21 kcal/mol with increasing temperature. The association constant (Ka) was 105M1 and was
essentially temperature independent. In addition, the free energy change corresponding to the Kawas not susceptible
to temperature changes, as a result of strong enthalpy– entropy compensation (Figure 4C).
The binding of hLon to LSPas showed a characteristic thermodynamic switch at Ts= 18.868C, where the
reac-tion entropy was zero. The zero-enthalpy temperature TH
was extrapolated to 6.988C. This binding reaction was enthalpy-driven and entropy-unfavorable when the tem-perature was above 18.868C, being entropy-driven and enthalpy-unfavorable when the temperature was below 6.988C. At temperatures between 6.988C and 18.868C, the reaction was both enthalpy- and entropy-driven.
This result demonstrates that the association of hLon with LSPas was dominated by a large negative enthalpy change and a reduction in entropy at physiologically relevant temperatures. The apparent heat capacity change evalu-ated by the temperature dependence of the enthalpy changes was 607.82 calK1mol1.
As demonstrated by EMSA, hLon formed a complex with TG6T (Figure 3B), which competed for hLon binding
to LSPas, albeit poorly (Figure S3A). We proceeded with a mechanistic study of hLon binding to TG6T under the
same conditions used to analyze LSPas binding. The thermodynamic data revealed that similar to hLon binding to LSPas, binding to TG6T was exothermic, and
exhibited moderate temperature dependence (Figure 4D and Table 2). However, in contrast to what was observed with LSPas, the heat capacity change for hLon binding to TG6T was signed positive and relatively small, giving a
value of 67.29 calK1mol1. The reaction entropy in the case of TG6T was positive over the temperature examined.
The two thermodynamic switch temperatures, THand TS,
were 38.658C and 30.868C, respectively (Figure 4F). Within the broad temperature range between 38.658C and 30.868C, the binding of hLon to TG6T was
enthalpy-and-entropy-driven, and was mostly dominated by the
entropy term. This reaction was mainly entropy-driven both near and above physiologically relevant tempera-tures. It was entropy-unfavorable at the temperature below 30.868C, and enthalpy-unfavorable at the tem-perature above 38.658C.
The magnitude of heat capacity change is generally regarded as a thermodynamic signature of hydrophobic interactions. It has been reported that ligand binding involving specific interactions is generally accompanied by a significant decrease in heat capacity (51–53). The heat capacity change for the binding of hLon to LSPas was large in comparison with that for the unfolding of globular proteins reported to be 12 calK1mol1 per residue, suggesting the possibility that the association of hLon with LSPas might be coupled with local folding of 51 amino acid residues (54,55). The large negative heat capacity change for their interactions reflects a significant decrease in the translational and rotational degrees of freedom, as a result of ensemble conformational restric-tion manifested by Equarestric-tion (9).
Cp¼ @H2 kT2 ¼ @S2 k 9
Figure 3. hLon binding to G-quartet DNA. (A) hLon binding to fluorescently labeled LSPas (lanes 1–4) and LSPasGA (lanes 5–8) probes was analyzed by electrophoresis mobility shift assay using 6% polyacrylamide gels. The unlabeled DNA competitors were added at a 50-fold excess the amount of the fluorescent probes. (B) A fixed amount of LSPas or TG6T (1 mM) was incubated with increasing concentration of hLon (0, 1.1, 5.5, 11, 16.5 and 22 mM). The reaction mixtures were analyzed on 6% polyacrylamide gels.
The scenario that reduction in thermal fluctuations as hLon binds to LSPas is consistent with the finding of this complex in temperature-dependent CD spectra (Figure S4A and B). The result indicated that the
association of hLon with LSPas provides a substantial enhancement of thermal stability, which was also observed for TG6T albeit to a lesser extent (Figure S4C). Moreover,
the CD spectral analysis of these hLon–DNA complexes
Figure 4. Isothermal titration calorimetric analysis for the binding of hLon to LSPas and TG6T. (A) LSPas and (D) TG6T were titrated (350 mM) into solutions containing hLon (30 mM); heat responses were registered isothermally at the temperatures as indicated. Integrals of the heat responses were normalized and fitted with a nonlinear least-square method to obtain thermodynamic parameters (B and E). The parameters including TS, H and G representing the association of hLon with LSPas (C) and TG6T (F) were plotted against temperature.
suggested that neither hLon nor the DNA molecules were subjected to global conformational transitions upon binding (Figure S5). The G-quartet structures recognized by hLon remain largely unperturbed as demonstrated by the CD spectra of DNA from 250 to 350 nm, within which DNA is sensitive to conformational changes whereas the protein has very weak signal (Figure S5). Although the CD intensity of DNA weakened in the presence of hLon, it underwent no structural transition (no chromatic shift was observed). The data suggest that the G-quartet structures bound by hLon are not disrupted; at least they are generally intact. However, it remains to be discovered whether the protein locally distorts or bends the G-quartet structures.
DSC analysis of the sequence-dependent heat capacity change
DSC was used to further investigate whether the significant decrease in heat capacity in the binding of hLon to LSPas was sequence-specific. The heat capacity of hLon as a function of temperature showed an accumula-tion of energy in the pretransiaccumula-tion region (Figure 5A). The transition temperature of hLon was 468C, and was followed by an abrupt decrease in heat capacity that indicated a step of dramatic exothermic reaction. The onset of this exothermic reaction occurred at 508C, which most likely resulted from the aggregation of hLon immediately after or partly during the thermal denatura-tion. When significant interactions between hLon and DNA occur, the complex can be populated with multiple microstates in the pretransition region [Equation (8)], which means an extra energy accumulation attributed to this binding. The effect of DNA binding on the overall heat capacity in the pretransition region from 108C to 308C was qualitatively assessed by Equation (8). Since the transition temperatures of the DNA oligonucleotides were all found beyond 808C, at which the heat capacity of buffer began to deviate from the baseline, we therefore undertook the evaluation by assuming that the differences caused by DNA were negligible in these cases.
The difference in the pretransition slope (@Cp/@T )
between hLon and the hLon/LSPas complex was +1.844 kcal/mol/K2 (Figure 5B and (@Cp/@T) in
Table 3). The positive value signified that more heat was absorbed to weaken the interactions between hLon and LSPas. In addition, the endothermic dissociation of the hLon-LSPas complex appeared to partly compensate the drastic decrease of heat capacity that arose from hLon aggregation. However, the gain of extra @Cp/@T for DNA
binding was not observed in other G-rich sequences. The (@Cp/@T) of the hLon–TG6T complex was 0.271 kcal/
mol/K2, signed negative and relatively small, and the hLon in complex with LS1as, LS2as and LS3as also showed little significant change in @Cp/@T (Figure 5C and
Table 3). The data of LS4as and LS5as are likely to be invalid because the hLon aggregated when either sequence was added in the solution prior to scan. In brief, these results revealed that hLon binding to LSPas was significantly involved in weak interactions as indicated by the relatively large positive change in the pretransition heat capacity relative to temperature, but it was not the case in other G-rich sequences although they formed G-quartets. The DSC data are fairly consistent with the ITC results: hLon binding to LSPas was accompanied by a highly negative heat capacity change, which was absent from the binding to TG6T. The combination of ITC and
DSC analyses sheds new light on the DNA-binding mechanism of hLon; whereby the specific binding of the protease to LSPas is coupled with structural adaptation as demonstrated by the large negative heat capacity change of the association. By contrast, the heat capacity measurements on hLon in complex with other G-rich sequences showed little or no significant change in heat capacity. Moreover, the dissociation constants for hLon binding to LSPas and TG6T at 208C were 3.30 and
7.35 mM, respectively. This result suggests that hLon has similar affinity for both sequences, but the binding to TG6T has a relatively negligible heat capacity change and
a lesser extent of structural stabilization, which can be referred to as rigid-body association.
DISCUSSION
Heat capacity change arising from a ligand-binding reaction can be attributed to polyelectrolyte and non-polyelectrolyte effects. Polyelectrolyte effects are a
Table 2. Thermodynamic parameters for the binding of hLon to LSPas and TG6T
Target DNA Temperature (8C) H (kcal/mol) TS (kcal/mol) Ka(10 5 M1) G (kcal/mol) Cp(cal/mol/K) LSPas 10 1.48 0.12 5.49 2.45 0.43 6.99 0.09 607.82 15 5.04 0.12 1.98 2.11 0.11 7.03 0.03 20 8.37 3.31 1.01 3.03 0.83 7.29 0.17 30 13.72 6.91 6.21 2.60 0.72 7.52 0.15 TG6T 5 2.04 0.75 4.14 1.24 0.58 6.49 0.21 67.29 10 1.84 0.88 4.71 1.35 0.89 6.58 0.36 20 1.30 0.14 5.69 1.36 0.22 6.80 0.09 30 0.63 0.02 7.05 3.48 0.36 7.69 0.06
Cp was obtained from linear regression of H against temperature. The binding thermodynamic parameters of both LSPas and TG6T were calculated on the monomeric concentrations of hLon and DNA. The LSPas constants corrected for the G-quartet fraction of LSPas are aH, T[S + (a/T + Rlna)], alnK, G-RTlna and aCp, respectively. a represents the reciprocal of the LSPas G-quartet fraction, which was estimated 1.15 as analyzed by ImageQuant ver. 5.2 (Molecular Dynamics) with Figure 3A lane 1 (G-quartet 0.87; monomer 0.13).
consequence of dehydration upon complex formation, and the relocation of water molecules from the solute hydration shell to the bulk solvent is entropy oriented. Nonpolyelectrolyte effects resulting from hydrophobic
and other weak interactions are directed by enthalpy changes. The influence of these effects on the binding of hLon to DNA is determined by thermodynamic-driving forces: the enthalpy and entropy changes. The polyelec-trolyte effect is noticeably temperature dependent com-pared to the nonpolyelectrolyte effect; apolar hydration heat capacity decreases with temperature increment, and thus the favorable gain of dehydration entropy for a complex formation is expected to increase with tempera-ture (56–58). Meanwhile, dehydration is likely to occur with a mild positive change in enthalpy. This inherent dependence of entropy and enthalpy on temperature gives rise to a partial positive heat capacity change. Nevertheless, hLon binding to LSPas was associated with a dramatic decrease in heat capacity; the ensemble enthalpy change was negative over the temperature examined, and the entropy change became quite negative at higher temperatures. The compensation temperature TC
(where the H and the TS converged) of hLon binding to LSPas was 12.798C, above which the enthalpy term predominated. By contrast, hLon binding to TG6T
was accompanied by a slight increase in heat capacity; its TCwas 5.158C, above which the reaction was primarily
driven by entropy. Taken together, these results suggested that hLon binding to LSPas revolves around nonpolyelec-trolyte effects, whereas the binding to TG6T is primarily
influenced by polyelectrolyte effects.
Despite striking differences in their heat capacity changes and major driving forces, hLon binding to LSPas and TG6T showed similar association constants
and free energy changes (Table 2). An underestimation of the LSPas parameters unavoidably occurs for hLon binding to the G-quartet form because because LSPas DNA was present in reactions as both monomeric and G-quartet forms. The thermodynamic parameters including the association constants and free energy changes are of the right order of magnitude; however, they should be corrected for the fraction of LSPas G-quartets present in order to compare more accurately the binding of LSPas
Figure 5. Heat capacity thermograms of hLon in the presence and absence of DNA. (A) The partial molar heat capacity of hLon is shown between 108C and 708C. The heat capacity of hLon in complex with LSPas (B) or TG6T (C) at a 1:1.5 molar ratio after buffer correction was normalized by hLon molarity.
Table 3. Determination of the pretransition slopes of heat capacity in hLon and hLon/DNA complexes
VP VL VPL h@Cp=@Ti (kcal/mol/K2) h@Cp=@Ti (kcal/mol/K2) hLon 0.745 0.871 hLon/LSPas 0.572 0.727 2.715 +1.844 hLon/TG6T 0.596 0.739 0.600 0.271 hLon/LS1as 0.581 0.730 0.711 0.160 hLon/LS2as 0.581 0.729 0.914 +0.043 hLon/LS3as 0.577 0.729 0.804 0.067 hLon/LS4as 0.589 0.730 0.199 0.672a hLon/LS5as 0.587 0.730 -0.556 1.427a Calculation of the partial-specific volume of hLon (VP) was performed with the program SEDNTERP (available at http://www.jphilo.mailway. com/download.htm). The V-bars of DNA molecules (VL) were estimated from the individual V-bars of bases, the values of which have been resolved by Dr Olwyn Byron (http://www.bbri.org/RASMB/ rasmb.html). The V-bars of the hLon/DNA mixtures (VPL) were calculated with the equation, (mpVPþmLVL)/(mp+ mL) (41). a
The of hLon/LS4as and hLon/LS5as are likely to be invalid because hLon aggregated upon binding.
and TG6T to hLon (see Table 2). In energetic terms, the
DNA-binding affinity of a protein is reflected by the free energy change, which is relatively insensitive to subtle changes in either binding component. It has been shown in recent studies that only a few crucial regions within a macromolecule complex determines the free energy of binding (55,59). The remaining regions are more relevant to binding specificity (stereochemistry), and contribute little to the overall G. The deletion or mutation of specificity-determining groups that do not damage the structural integrity of the complex would not likely affect the binding free energy; thus thermodynamic homeostasis could be maintained. Nonetheless, the H and the S would be susceptible to minor structural changes, even though these two terms may compensate in a similar fashion resulting in an invariant free energy change.
The thermodynamic ‘neutral’ was observed in the binding of hLon to TG6T, whose free energy change was
comparable to that of LSPas. This phenomenon revealed that TG6T, the core region of G-quartet formation
residing in LSPas, is essential for the stabilization of hLon–DNA complex by conferring most of the binding energy. As indicated by EMSA results, disrupting the G-quartet in LSPas by replacing two guanines with adenines interfered with binding to hLon. This result suggested that the DNA substrate determines the binding preference of hLon at free energy level. In other words, hLon binding to LSPas is significantly weakened if LSPas loses its ability to form G-quartet. However, the binding specificity should be addressed by its thermodynamic properties and heat capacity change. The weak depen-dence of entropy and enthalpy on temperature for hLon binding to TG6T parallels the association with LSPas at
low temperatures. But at high temperatures, weak bond formation apparently dominates the interactions between hLon and LSPas, rendering more negative the enthalpy as well as the entropy of reaction; while in the case of TG6T,
dehydration appears to dominate the thermodynamics of this complex formation. In addition, the drastic difference in heat capacity change between these two cases signifies their fundamental discrepancy in the degree of tempera-ture dependence for the individual thermodynamic para-meters including the enthalpy and entropy changes. For a ligand-binding reaction, heat capacity change is primarily concerned with hydrophobic interactions and dehydration of apolar groups. Summing up the thermodynamic results, we infer that the small change in heat capacity for the binding of hLon to TG6T mainly accounts for the extent
of dehydration, whereas the binding to LSPas associated with the large negative heat capacity change reflects a considerable involvement of hydrophobic interactions. Supported further by the DSC results, these findings suggest that hLon binds to most G-quartet-forming sequences present in mtDNA principally through electro-static interactions which are thermodynamically inade-quate for specific binding. By contrast, hLon recognizes LSPas as a specific target by enhancing the structural stability of the complex. In addition, these results raise the possibility that a sequence-specific recognition step might be preceded by nonspecific binding to G-quartet structures (60).
The human mitochondrial genome (mtDNA) consists of 16,569 bp. Although the base composition of mamma-lian mtDNA is biased with a high frequency of contiguous guanines on the heavy (H) strand as compared with the complementary light (L) strand, only a few studies have addressed the physiological implications of G-quartets in this closed circular genome (61,62). Analysis of human mtDNA with two G-quartet motif-searching programs, Quadfinder and QGRS Mapper, revealed 20 intramole-cular G-quartet-forming hits on the H strand, while none were found on the L strand (Table S1) (63,64). Fourteen out of these 20 sequences are located between LSP and OL
on the heavy strand. LSPas is contained within one of these hits positioned between 438 and 471.
One model of mtDNA replication as proposed by Clayton and the colleagues (65,66) holds that replication proceeds in an asymmetric manner. Synthesis is initiated at the origin of the H strand (OH) that is downstream from
the LSP. The single-stranded parental H strand is displaced during replication, and expands until approxi-mately two-thirds the genome is replicated reaching the origin of the L strand (OL), at which point L-strand
replication ensues in the opposite direction. Studies have found that replication of the leading H strand terminates prematurely at TAS (termination-associated sequences), resulting in a 7S DNA; this three-stranded structure referred to as D-loop (displacement-loop), consists of the parental L-strand template, the nascent prematurely terminated H strand and the displaced single-stranded parental H strand (Table 1) (67,68). It has been proposed that the rate of mtDNA replication is regulated by the turnover of the D-loop, which is likely mediated by the formation of secondary structures at certain conserved sequences within the D-loop (69,70) and some trans-acting factors binding to these regions (71–73). We suggest that the single-stranded parental H strand exposed to solution might allow the formation of G-quartets at various G-rich clusters between LSP and OLas identified by Quadfinder
(Table S1). Recent findings by Maizels and the colleagues (74) on the G-quartet-containing structures called G-loops may strengthen the validity of our suggestion. However, it is unclear whether hLon interacts with intramolecular LSPas quadruplex, and if the binding is accompanied by similar heat capacity effect: it remains to be clarified whether the specific recognition of the LSPas sequence requires the flanking bases to present quadruply. Although it has been debated on the biological relevance of intramolecular parallel G-quartets arising from topo-logical constraints, a FRET-based investigation reveals that parallel and antiparallel G-quartet conformations can coexist and interconvert under physiological relevant conditions (75). In our study, the combination of electrophoretic, spectroscopic and thermodynamic data indicates that hLon preferentially binds to intermolecular parallel G-quartet structures. Given that multiple copies of mtDNA molecules frequently concentrate to form heritable units of mtDNA or to facilitate intermitochon-drial recombination (76–78), formation of G-quartets between separate mtDNA molecules could be possible in vivo(see the schematic structures in Figure 2). We have found that the G-quartet-forming motifs on mtDNA are
thermodynamically favorable to form intermolecular G-quartets while kinetically prone to form intramolecular G-quartets in vitro (our unpublished data). However, whether mtDNA forms G-quartets intra- and/or inter-genomically remains to be determined.
This study reveals that hLon has a general preference for G-quartet DNA, and selectively binds to LSPas along with structural adaptations for the overall structural tightening. Our findings are consistent with recent work identifying a G-rich consensus sequence for hLon binding that matches well with a sequence contained within LSPas (28). Additional consensus sequences showing the best fit for hLon binding sites are primarily located downstream of OH and before OL. In vitro, purified hLon binds with
higher affinity to ssDNA oligonucleotides that conform to this G-rich consensus. Consistent with these findings, hLon binds preferentially in cultured human cells to the CR/D-loop region as well as to sites between OHand OL.
Our data support the notion that hLon binding to mtDNA sequences in vitro and in vivo is determined by the propensity of ssDNA to form G-quartets. We demonstrate that hLon is able to recognize its binding target with a notable reduction in thermal fluctuations. In addition, the high proportion of G-rich clusters in the H strand raises the possibility that these sites serve as hot spots of binding for hLon. We propose that hLon associates with the G-quartet-forming regions during the expansion of D-loop, and selectively binds to target sequence(s) on the basis of different binding heat capacities. Thus, the finding of specific target(s) by hLon to mtDNA might be facilitated by initial binding to these hot spots rather than random diffusion or one-dimensional sliding (79,80).
Potential G-quartet-forming regions that have been characterized in eukaryotic nuclear genome are principally found in telomeres, immunoglobulin heavy-chain switch regions, ribosomal DNA, G-rich microsatellites and certain promoter regions (81,82). Telomeres are guanine-rich tandem repeats at the ends of linear chromosomes, which tend to form G-quartet structures to prevent the 30
single-stranded overhang from nonhomologous end join-ing. Telomere integrity is regulated by some G-quartet destabilizing factors as well as telomerase. Telomere-binding proteins such as Pot1 and G-quartet-specific helicases like BLM and WRN syndrome proteins render the telomeric DNA accessible to telomerase for end extension (83). Deregulations of telomerase activity and the helicase-mediated unwinding of G-quartets at telo-meres are potentially associated with many types of cancer. On the other hand, many putative G-quadruplexes have been identified in nontelomeric regions, of which promoter sequences are intensely concerned with the correlation between the G-quartet stability and gene expression. Stabilization and destabilization of the G-quartet structures at the promoter regions of c-myc and Ki-ras proto-oncogenes have been revealed to be impor-tant to the gene regulation (84,85). It is of great interest to discover how gene expression is regulated by some G-quartet-binding proteins that are involved in the turnover of G-quartet structures at the control regions.
The specialized nucleic acid structures found in nuclear genome broaches an intriguing issue of whether they exist in organelle genomes, and if they are under similar or shared regulation. Unlike nuclear DNA, typical mito-chondrial genomes are circular and thus do not require telomere/telomerase system for an end-replication pro-blem. However, linear mtDNA molecules are found in certain yeast strains, which have mitochondrial telomere-binding proteins for maintaining their genome stability (86,87). Our data raise the possibility that some regions of mtDNA may form G-quadruplexes. Although this con-jecture awaits further experimental validation, recent findings on localization of potential G-quartet-binding proteins at mitochondria suggest they might interact with G-quartet structures on mtDNA. Telomere reverse transcriptase (TERT) that normally guards the end integrity of chromosomes in nuclei carries a mitochondrial targeting sequence and retains telomerase activity in mitochondria (88,89). Translocation of TERT from nuclei to mitochondria is triggered by oxidative stress and exacerbates mtDNA damage, in contrast to its protective role in nuclei. Topoisomerase I (Top1) has also been identified in mitochondria with yet unclear functions (90,91). The existence of mitochondrial Top1 might be relevant to the topologic problems encountered during mtDNA replication/transcription. Its nuclear counterpart interacts with and promotes the formation of G-quartet structures in vitro (92). However, whether the mitochondrial Top1 as well as TERT associates with G-quartet structures on mtDNA is unknown.
The dissociation constants of hLon binding to LSPas and TG6T were found at micromolar range by using ITC.
The binding strength seems to be weak, but reasonable because hLon most likely functions as a regulatory protein. It has been reported that hLon binding to DNA is inhibited by nucleotides whereas is stimulated by protein substrates (25). The physiological functions of mtDNA bound by hLon are still unclear. As hLon functions as an ATP-dependent protease and is a likely sensor of environmental stress in mitochondria, it is possible that hLon selectively degrades components of the mtDNA replication, repair and transcription machinery in response to changes in cellular metabolism. A variety of mitochondrial intermediate enzymes seemingly unrelated to mtDNA transactions such as aconitase, acetohy-droxyacid reductoisomerase and serine/threonine deami-nase have been identified in mitochondrial nucleoids (34,93–95), suggesting that they may be important in coordinating metabolic regulation with mtDNA main-tenance. Further investigation is required to determine whether these respective enzymatic activities are function-ally linked to their association with mtDNA binding, and whether they are selectively degraded and thereby regulated by hLon during cellular stress such as the increased generation of reactive oxygen species and/or hypoxia.
SUPPLEMENTARY DATA
ACKNOWLEDGEMENTS
The authors thank Dr Ming-Chya Wu for access to the differential scanning calorimeter. We are also grateful to Dr Lou-Sing Kan, Dr Chia-Ching Chang, Dr Jiunly Chir and Dr Kamal Singh for helpful comments, and to Drs Hui-Chuan Chang and Chris S.-C. Jao for technical support in circular dichronism and isothermal titration calorimeter. This study was funded by National Science Council, Taiwan (96-2311-B-001-010 to S.-H.W.). Funding to pay the Open Access publication charges for this article was provided by Academia Sinica, Taiwan. Conflict of interest statement. None declared.
REFERENCES
1. Hanson,P.I. and Whiteheart,S.W. (2005) AAA+ proteins: have engine, will work. Nat. Rev. Mol. Cell Biol., 6, 519–529. 2. Fukui,T., Eguchi,T., Atomi,H. and Imanaka,T. (2002) A membrane-bound archaeal Lon protease displays ATP-independent proteolytic activity towards unfolded proteins and ATP-dependent activity for folded proteins. J. Bacteriol., 184, 3689–3698.
3. Besche,H., Tamura,N., Tamura,T. and Zwickl,P. (2004)
Mutational analysis of conserved AAA+ residues in the archaeal Lon protease from Thermoplasma acidophilum. FEBS Lett., 574, 161–166.
4. Chung,C.H. and Goldberg,A.L. (1981) The product of the lon (capR) gene in Escherichia coli is the ATP-dependent protease, protease La. Proc. Natl Acad. Sci. USA, 78, 4931–4935. 5. Chin,D.T., Goff,S.A., Webster,T., Smith,T. and Goldberg,A.L.
(1988) Sequence of the lon gene in Escherichia coli. A heat-shock gene which encodes the ATP-dependent protease La. J. Biol. Chem., 263, 11718–11728.
6. Van Dyck,L., Pearce,D.A. and Sherman,F. (1994) PIM1 encodes a mitochondrial ATP-dependent protease that is required for mito-chondrial function in the yeast Saccharomyces cerevisiae. J. Biol. Chem., 269, 238–242.
7. Lee,A.Y., Tsay,S.S., Chen,M.Y. and Wu,S.H. (2004) Identification of a gene encoding Lon protease from Brevibacillus thermoruber WR-249 and biochemical characterization of its thermostable recombinant enzyme. Eur. J. Biochem., 271, 834–844.
8. Kikuchi,M., Hatano,N., Yokota,S., Shimozawa,N., Imanaka,T. and Taniguchi,H. (2004) Proteomic analysis of rat liver peroxisome: presence of peroxisome-specific isozyme of Lon protease. J. Biol. Chem., 279, 421–428.
9. Ostersetzer,O., Kato,Y., Adam,Z. and Sakamoto,W. (2007) Multiple intracellular locations of Lon protease in Arabidopsis: evidence for the localization of AtLon4 to chloroplasts. Plant Cell Physiol., 48, 881–885.
10. Charette,M.F., Henderson,G.W. and Markovitz,A. (1981) ATP hydrolysis-dependent protease activity of the lon (capR) protein of Escherichia coli K-12. Proc. Natl Acad. Sci. USA, 78, 4728–4732.
11. Waxman,L. and Goldberg,A.L. (1986) Selectivity of intracellular proteolysis: protein substrates activate the ATP-dependent protease (La). Science, 232, 500–503.
12. Mizusawa,S. and Gottesman,S. (1983) Protein degradation in Escherichia coli: the lon gene controls the stability of sulA protein. Proc. Natl Acad. Sci. USA, 80, 358–362.
13. Maurizi,M.R. (1987) Degradation in vitro of bacteriophage lambda N protein by Lon protease from Escherichia coli. J. Biol. Chem., 262, 2696–2703.
14. Bonnefoy,E., Almeida,A. and Rouviere-Yaniv,J. (1989) Lon-dependent regulation of the DNA binding protein HU in Escherichia coli. Proc. Natl Acad. Sci. USA, 86, 7691–7695. 15. Wright,R., Stephens,C., Zweiger,G., Shapiro,L. and Alley,M.R.
(1996) Caulobacter Lon protease has a critical role in cell-cycle control of DNA methylation. Genes Dev., 10, 1532–1542.
16. Torres-Cabassa,A.S. and Gottesman,S. (1987) Capsule synthesis in Escherichia coli K-12 is regulated by proteolysis. J. Bacteriol., 169, 981–989.
17. Stout,V., Torres-Cabassa,A., Maurizi,M.R., Gutnick,D. and Gottesman,S. (1991) RcsA, an unstable positive regulator of capsular polysaccharide synthesis. J. Bacteriol., 173, 1738–1747. 18. Whistler,C.A., Stockwell,V.O. and Loper,J.E. (2000) Lon
protease influences antibiotic production and UV tolerance of Pseudomonas fluorescens Pf-5. Appl. Environ. Microbiol., 66, 2718–2725.
19. Bota,D.A. and Davies,K.J. (2002) Lon protease preferentially degrades oxidized mitochondrial aconitase by an ATP-stimulated mechanism. Nat. Cell Biol., 4, 674–680.
20. Rep,M., van Dijl,J.M., Suda,K., Schatz,G., Grivell,L.A. and Suzuki,C.K. (1996) Promotion of mitochondrial membrane complex assembly by a proteolytically inactive yeast Lon. Science, 274, 103–106.
21. Hori,O., Ichinoda,F., Tamatani,T., Yamaguchi,A., Sato,N., Ozawa,K., Kitao,Y., Miyazaki,M., Harding,H.P. et al. (2002) Transmission of cell stress from endoplasmic reticulum to mito-chondria: enhanced expression of Lon protease. J. Cell Biol., 157, 1151–1160.
22. Fu,G.K., Smith,M.J. and Markovitz,D.M. (1997) Bacterial protease Lon is a site-specific DNA-binding protein. J. Biol. Chem., 272, 534–538.
23. Fu,G.K. and Markovitz,D.M. (1998) The human LON protease binds to mitochondrial promoters in a single-stranded, site-specific, strand-specific manner. Biochemistry, 37, 1905–1909.
24. Lu,B., Liu,T., Crosby,J.A., Thomas-Wohlever,J., Lee,I. and Suzuki,C.K. (2003) The ATP-dependent Lon protease of Mus musculus is a DNA-binding protein that is functionally conserved between yeast and mammals. Gene, 306, 45–55.
25. Liu,T., Lu,B., Lee,I., Ondrovicova,G., Kutejova,E. and Suzuki,C.K. (2004) DNA and RNA binding by the mitochondrial lon protease is regulated by nucleotide and protein substrate. J. Biol. Chem., 279, 13902–13910.
26. Nomura,K., Kato,J., Takiguchi,N., Ohtake,H. and Kuroda,A. (2004) Effects of inorganic polyphosphate on the proteolytic and DNA-binding activities of Lon in Escherichia coli. J. Biol. Chem., 279, 34406–34410.
27. Lee,A.Y., Hsu,C.H. and Wu,S.H. (2004) Functional domains of Brevibacillus thermoruber lon protease for oligomerization and DNA binding: role of N-terminal and sensor and substrate discrimination domains. J. Biol. Chem., 279, 34903–34912. 28. Lu,B., Yadav,S., Shah,P.G., Liu,T., Tian,B., Pukszta,S.,
Villaluna,N., Kutejova,E., Newlon,C.S. et al. (2007) Roles for the human ATP-dependent Lon protease in mitochondrial DNA maintenance. J. Biol. Chem., 282, 17363–17374.
29. Shah,I.M. and Wolf,R.E. Jr. (2006) Inhibition of Lon-dependent degradation of the Escherichia coli transcription activator SoxS by interaction with ‘soxbox’ DNA or RNA polymerase.
Mol. Microbiol., 60, 199–208.
30. Shadel,G.S. (2005) Mitochondrial DNA, aconitase ‘wraps’ it up. Trends Biochem. Sci., 30, 294–296.
31. Chen,X.J. and Butow,R.A. (2005) The organization and inheritance of the mitochondrial genome. Nat. Rev. Genet., 6, 815–825. 32. Chen,X.J., Wang,X., Kaufman,B.A. and Butow,R.A. (2005)
Aconitase couples metabolic regulation to mitochondrial DNA maintenance. Science, 307, 714–717.
33. Yan,L.J., Levine,R.L. and Sohal,R.S. (1997) Oxidative damage during aging targets mitochondrial aconitase. Proc. Natl Acad. Sci. USA, 94, 11168–11172.
34. Wang,Y. and Bogenhagen,D.F. (2006) Human mitochondrial DNA nucleoids are linked to protein folding machinery and metabolic enzymes at the mitochondrial inner membrane. J. Biol. Chem., 281, 25791–25802.
35. Bota,D.A., Ngo,J.K. and Davies,K.J. (2005) Downregulation of the human Lon protease impairs mitochondrial structure and function and causes cell death. Free Radic. Biol. Med., 38, 665–677.
36. Chung,C.H. and Goldberg,A.L. (1982) DNA stimulates ATP-dependent proteolysis and protein-ATP-dependent ATPase activity of protease La from Escherichia coli. Proc. Natl Acad. Sci. USA, 79, 795–799.
37. Bradford,M.M. (1976) A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem., 72, 248–254. 38. Fasman,G.D. and Chemical Rubber Company. (1975) Handbook of
Biochemistry and Molecular Biology: Nucleic Acids,3rd edn. CRC Press, Cleveland, p.589.
39. Hassur,S.M. and Whitlock,H.W. Jr. (1974) UV shadowing – a new and convenient method for the location of ultraviolet-absorbing species in polyacrylamide gels. Anal. Biochem., 59, 162–164. 40. Jelesarov,I. and Bosshard,H.R. (1999) Isothermal titration
calori-metry and differential scanning caloricalori-metry as complementary tools to investigate the energetics of biomolecular recognition. J. Mol. Recognit., 12, 3–18.
41. Hinz,H.J. (1986) Thermodynamic parameters for protein-protein and protein-ligand interaction by differential scanning microcalori-metry. Methods Enzymol., 130, 59–79.
42. Ladbury,J.E. and Chowdhry,B.Z. (1998) Biocalorimetry: Applications of Calorimetry in the Biological Sciences.J. Wiley, Chichester, New York, pp.183–205.
43. Privalov,P.L. and Potekhin,S.A. (1986) Scanning microcalorimetry in studying temperature-induced changes in proteins. Methods Enzymol., 131, 4–51.
44. Kasamatsu,H., Robberson,D.L. and Vinograd,J. (1971) A novel closed-circular mitochondrial DNA with properties of a replicating intermediate. Proc. Natl Acad. Sci. USA, 68, 2252–2257.
45. Williamson,J.R. (1994) G-quartet structures in telomeric DNA. Annu. Rev. Biophys. Biomol. Struct., 23, 703–730.
46. Keniry,M.A. (2000) Quadruplex structures in nucleic acids. Biopolymers, 56, 123–146.
47. Kang,C., Zhang,X., Ratliff,R., Moyzis,R. and Rich,A. (1992) Crystal structure of four-stranded Oxytricha telomeric DNA. Nature, 356, 126–131.
48. Laughlan,G., Murchie,A.I., Norman,D.G., Moore,M.H., Moody,P.C., Lilley,D.M. and Luisi,B. (1994) The high-resolution crystal structure of a parallel-stranded guanine tetraplex. Science, 265, 520–524.
49. Hud,N.V., Schultze,P., Sklenar,V. and Feigon,J. (1999) Binding sites and dynamics of ammonium ions in a telomere repeat DNA quadruplex. J. Mol. Biol., 285, 233–243.
50. Wong,A. and Wu,G. (2003) Selective binding of monovalent cations to the stacking G-quartet structure formed by guanosine 50-monophosphate: a solid-state NMR study. J. Am. Chem. Soc., 125, 13895–13905.
51. Ha,J.H., Spolar,R.S. and Record,M.T. Jr. (1989) Role of the hydrophobic effect in stability of site-specific protein-DNA com-plexes. J. Mol. Biol., 209, 801–816.
52. Spolar,R.S. and Record,M.T. Jr. (1994) Coupling of local folding to site-specific binding of proteins to DNA. Science, 263, 777–784. 53. Ladbury,J.E., Wright,J.G., Sturtevant,J.M. and Sigler,P.B. (1994)
A thermodynamic study of the trp repressor-operator interaction. J. Mol. Biol., 238, 669–681.
54. Cooper,A. (2000) Heat capacity of hydrogen-bonded networks: an alternative view of protein folding thermodynamics. Biophys. Chem., 85, 25–39.
55. Cooper,A., Johnson,C.M., Lakey,J.H. and Nollmann,M. (2001) Heat does not come in different colours: entropy-enthalpy compensation, free energy windows, quantum confinement, pressure perturbation calorimetry, solvation and the multiple causes of heat capacity effects in biomolecular interactions. Biophys. Chem., 93, 215–230.
56. Baldwin,R.L. (1986) Temperature dependence of the hydrophobic interaction in protein folding. Proc. Natl Acad. Sci. USA, 83, 8069–8072.
57. Privalov,P.L. and Makhatadze,G.I. (1990) Heat capacity of proteins. II. Partial molar heat capacity of the unfolded polypeptide chain of proteins: protein unfolding effects. J. Mol. Biol., 213, 385–391.
58. Prabhu,N.V. and Sharp,K.A. (2005) Heat capacity in proteins. Annu. Rev. Phys. Chem., 56, 521–548.
59. Bogan,A.A. and Thorn,K.S. (1998) Anatomy of hot spots in protein interfaces. J. Mol. Biol., 280, 1–9.
60. Kalodimos,C.G., Biris,N., Bonvin,A.M., Levandoski,M.M., Guennuegues,M., Boelens,R. and Kaptein,R. (2004) Structure and
flexibility adaptation in nonspecific and specific protein-DNA complexes. Science, 305, 386–389.
61. Anderson,S., Bankier,A.T., Barrell,B.G., de Bruijn,M.H., Coulson,A.R., Drouin,J., Eperon,I.C., Nierlich,D.P., Roe,B.A. et al. (1981) Sequence and organization of the human mitochondrial genome. Nature, 290, 457–465.
62. Nonin-Lecomte,S., Dardel,F. and Lestienne,P. (2005) Self-organisation of an oligodeoxynucleotide containing the G- and C-rich stretches of the direct repeats of the human mitochondrial DNA. Biochimie, 87, 725–735.
63. Scaria,V., Hariharan,M., Arora,A. and Maiti,S. (2006) Quadfinder: server for identification and analysis of quadruplex-forming motifs in nucleotide sequences. Nucleic Acids Res., 34, W683–W685. 64. Kikin,O., D’Antonio,L. and Bagga,P.S. (2006) QGRS Mapper: a
web-based server for predicting G-quadruplexes in nucleotide sequences. Nucleic Acids Res., 34, W676–W682.
65. Brown,T.A., Cecconi,C., Tkachuk,A.N., Bustamante,C. and Clayton,D.A. (2005) Replication of mitochondrial DNA occurs by strand displacement with alternative light-strand origins, not via a strand-coupled mechanism. Genes Dev., 19, 2466–2476.
66. Clayton,D.A. (2003) Mitochondrial DNA replication: what we know. IUBMB Life, 55, 213–217.
67. Doda,J.N., Wright,C.T. and Clayton,D.A. (1981) Elongation of displacement-loop strands in human and mouse mitochondrial DNA is arrested near specific template sequences. Proc. Natl Acad. Sci. USA, 78, 6116–6120.
68. Clayton,D.A. (1991) Replication and transcription of vertebrate mitochondrial DNA. Annu. Rev. Cell Biol., 7, 453–478.
69. Brown,G.G., Gadaleta,G., Pepe,G., Saccone,C. and Sbisa,E. (1986) Structural conservation and variation in the D-loop-containing region of vertebrate mitochondrial DNA. J. Mol. Biol., 192, 503–511.
70. Nass,M.M. (1995) Precise sequence assignment of replication origin in the control region of chick mitochondrial DNA relative to 50and 30 D-loop ends, secondary structure, DNA synthesis, and protein binding. Curr. Genet., 28, 401–409.
71. Madsen,C.S., Ghivizzani,S.C. and Hauswirth,W.W. (1993) Protein binding to a single termination-associated sequence in the mito-chondrial DNA D-loop region. Mol. Cell. Biol., 13, 2162–2171. 72. Shadel,G.S. and Clayton,D.A. (1997) Mitochondrial DNA
main-tenance in vertebrates. Annu. Rev. Biochem., 66, 409–435. 73. Roberti,M., Musicco,C., Polosa,P.L., Milella,F., Gadaleta,M.N.
and Cantatore,P. (1998) Multiple protein-binding sites in the TAS-region of human and rat mitochondrial DNA. Biochem. Biophys. Res. Commun., 243, 36–40.
74. Duquette,M.L., Handa,P., Vincent,J.A., Taylor,A.F. and Maizels,N. (2004) Intracellular transcription of G-rich DNAs induces formation of G-loops, novel structures containing G4 DNA. Genes Dev., 18, 1618–1629.
75. Ying,L., Green,J.J., Li,H., Klenerman,D. and Balasubramanian,S. (2003) Studies on the structure and dynamics of the human telomeric G quadruplex by single-molecule fluorescence resonance energy transfer. Proc. Natl Acad. Sci. USA, 100, 14629–14634.
76. Lockshon,D., Zweifel,S.G., Freeman-Cook,L.L., Lorimer,H.E., Brewer,B.J. and Fangman,W.L. (1995) A role for recombination junctions in the segregation of mitochondrial DNA in yeast. Cell, 81, 947–955.
77. MacAlpine,D.M., Perlman,P.S. and Butow,R.A. (2000) The num-bers of individual mitochondrial DNA molecules and mitochondrial DNA nucleoids in yeast are co-regulated by the general amino acid control pathway. EMBO J., 19, 767–775.
78. Gibson,T., Blok,V.C., Phillips,M.S., Hong,G., Kumarasinghe,D., Riley,I.T. and Dowton,M. (2007) The mitochondrial subgenomes of the nematode Globodera pallida are mosaics: evidence of recombi-nation in an animal mitochondrial genome. J. Mol. Evol., 64, 463–471.
79. Halford,S.E. and Marko,J.F. (2004) How do site-specific DNA-binding proteins find their targets? Nucleic Acids Res., 32, 3040–3052.
80. Flyvbjerg,H., Keatch,S.A. and Dryden,D.T. (2006) Strong physical constraints on sequence-specific target location by proteins on DNA molecules. Nucleic Acids Res., 34, 2550–2557.
81. Maizels,N. (2006) Dynamic roles for G4 DNA in the biology of eukaryotic cells. Nat. Struct. Mol. Biol., 13, 1055–1059. 82. Oganesian,L. and Bryan,T.M. (2007) Physiological relevance of
telomeric G-quadruplex formation: a potential drug target. Bioessays, 29, 155–165.
83. Opresko,P.L., Mason,P.A., Podell,E.R., Lei,M., Hickson,I.D., Cech,T.R. and Bohr,V.A. (2005) POT1 Stimulates recQ helicases WRN and BLM to unwind telomeric DNA substrates. J. Biol. Chem., 280, 32069–32080.
84. Siddiqui-Jain,A., Grand,C.L., Bearss,D.J. and Hurley,L.H. (2002) Direct evidence for a G-quadruplex in a promoter region and its targeting with a small molecule to repress c-MYC transcription. Proc. Natl Acad. Sci. USA, 99, 11593–11598.
85. Cogoi,S., Quadrifoglio,F. and Xodo,L.E. (2004) G-rich oligonu-cleotide inhibits the binding of a nuclear protein to the Ki-ras promoter and strongly reduces cell growth in human carcinoma pancreatic cells. Biochemistry, 43, 2512–2523.
86. Nosek,J., Tomaska,L., Pagacova,B. and Fukuhara,H. (1999) Mitochondrial telomere-binding protein from Candida
parapsilosis suggests an evolutionary adaptation of a nonspecific single-stranded DNA-binding protein. J. Biol. Chem., 274, 8850–8857.
87. Tomaska,L., Makhov,A.M., Nosek,J., Kucejova,B. and
Griffith,J.D. (2001) Electron microscopic analysis supports a dual role for the mitochondrial telomere-binding protein of Candida parapsilosis. J. Mol. Biol., 305, 61–69.
88. Santos,J.H., Meyer,J.N., Skorvaga,M., Annab,L.A. and Van Houten,B. (2004) Mitochondrial hTERT exacerbates free-radical-mediated mtDNA damage. Aging Cell, 3, 399–411.
89. Santos,J.H., Meyer,J.N. and Van Houten,B. (2006) Mitochondrial localization of telomerase as a determinant for hydrogen peroxide-induced mitochondrial DNA damage and apoptosis. Hum. Mol. Genet., 15, 1757–1768.
90. Tua,A., Wang,J., Kulpa,V. and Wernette,C.M. (1997)
Mitochondrial DNA topoisomerase I of Saccharomyces cerevisiae. Biochimie, 79, 341–350.
91. Zhang,H., Meng,L.H. and Pommier,Y. (2007) Mitochondrial topoisomerases and alternative splicing of the human TOP1mt gene. Biochimie, 89, 474–481.
92. Arimondo,P.B., Riou,J.F., Mergny,J.L., Tazi,J., Sun,J.S., Garestier,T. and Helene,C. (2000) Interaction of human DNA topoisomerase I with G-quartet structures. Nucleic Acids Res., 28, 4832–4838.
93. Kaufman,B.A., Newman,S.M., Hallberg,R.L., Slaughter,C.A., Perlman,P.S. and Butow,R.A. (2000) In organello formaldehyde crosslinking of proteins to mtDNA: identification of bifunctional proteins. Proc. Natl Acad. Sci. USA, 97, 7772–7777.
94. Bogenhagen,D.F., Wang,Y., Shen,E.L. and Kobayashi,R. (2003) Protein components of mitochondrial DNA nucleoids in higher eukaryotes. Mol. Cell Proteomics, 2, 1205–1216.
95. Malka,F., Lombes,A. and Rojo,M. (2006) Organization, dynamics and transmission of mitochondrial DNA: focus on vertebrate nucleoids. Biochim. Biophys. Acta, 1763, 463–472.
