Determination of small crystal structures from a minimum set of diffraction intensities by homotopy continuation

No deterministic approach to obtaining a crystal structure from a set of diffraction intensities exists, despite significant progress in traditional probabilistic direct methods. One of the biggest hurdles in determining a crystal structure algebraically is solving a system of many polynomial equations of high power on intensities in terms of atomic coordinates. In this study, homotopy continuation is used for exhaustive investigation of such systems and an optimized homotopy continuation method is developed with random restarts to determine small (N< 5) crystal structures from a minimum set of error-free intensities.

Structural insight into MtmC, a bifunctional ketoreductase-methyltransferase involved in the assembly of the mithramycin trisaccharide chain

More and more post-PKS tailoring enzymes are recognized as being multifunctional and codependent on other tailoring enzymes. One of the recently discovered intriguing examples is MtmC, a bifunctional TDP-4-keto-d-olivose ketoreductase-methyltransferase, which-in codependence with glycosyltransferase MtmGIV-is a key contributor to the biosynthesis of the critical trisaccharide chain of the antitumor antibiotic mithramycin (MTM), produced by Streptomyces argillaceus. We report crystal structures of three binary complexes of MtmC with its methylation cosubstrate SAM, its coproduct SAH, and a nucleotide TDP as well as crystal structures of two ternary complexes, MtmC-SAH-TDP-4-keto-d-olivose and MtmC-SAM-TDP, in the range of 2.2-2.7 Å resolution. The structures reveal general and sugar-specific recognition and catalytic structural features of MtmC. Depending on the catalytic function that is conducted by MtmC, it must bind either NADPH or SAM in the same cofactor binding pocket. A tyrosine residue (Tyr79) appears as a lid covering the sugar moiety of the substrate during the methyl transfer reaction. This residue swings out of the active site by ~180° in the absence of the substrate. This unique conformational change likely serves to release the methylated product and, possibly, to open the active site for binding the bulkier cosubstrate NADPH prior to the reduction reaction.

Structure of inorganic pyrophosphatase from Staphylococcus aureus reveals conformational flexibility of the active site

Cytoplasmic inorganic pyrophosphatase (PPiase) is an enzyme essential for survival of organisms, from bacteria to human. PPiases are divided into two structurally distinct families: family I PPiases are Mg(2+)-dependent and present in most archaea, eukaryotes and prokaryotes, whereas the relatively less understood family II PPiases are Mn(2+)-dependent and present only in some archaea, bacteria and primitive eukaryotes. Staphylococcus aureus (SA), a dangerous pathogen and a frequent cause of hospital infections, contains a family II PPiase (PpaC), which is an attractive potential target for development of novel antibacterial agents. We determined a crystal structure of SA PpaC in complex with catalytic Mn(2+) at 2.1Å resolution. The active site contains two catalytic Mn(2+) binding sites, each half-occupied, reconciling the previously observed 1:1 Mn(2+):enzyme stoichiometry with the presence of two divalent metal ion sites in the apo-enzyme. Unexpectedly, despite the absence of the substrate or products in the active site, the two domains of SA PpaC form a closed active site, a conformation observed in structures of other family II PPiases only in complex with substrate or product mimics. A region spanning residues 295-298, which contains a conserved substrate binding RKK motif, is flipped out of the active site, an unprecedented conformation for a PPiase. Because the mutant of Arg295 to an alanine is devoid of activity, this loop likely undergoes an induced-fit conformational change upon substrate binding and product dissociation. This closed conformation of SA PPiase may serve as an attractive target for rational design of inhibitors of this enzyme.

Antimycobacterial activity of DNA intercalator inhibitors of Mycobacterium tuberculosis primase DnaG

Owing to the rise in drug resistance in tuberculosis combined with the global spread of its causative pathogen, Mycobacterium tuberculosis (Mtb), innovative anti mycobacterial agents are urgently needed. Recently, we developed a novel primase-pyrophosphatase assay and used it to discover inhibitors of an essential Mtb enzyme, primase DnaG (Mtb DnaG), a promising and unexplored potential target for novel antituberculosis chemotherapeutics. Doxorubicin, an anthracycline antibiotic used as an anticancer drug, was found to be a potent inhibitor of Mtb DnaG. In this study, we investigated both inhibition of Mtb DnaG and the inhibitory activity against in vitro growth of Mtb and M. smegmatis (Msm) by other anthracyclines, daunorubicin and idarubicin, as well as by less cytotoxic DNA intercalators: aloe-emodin, rhein and a mitoxantrone derivative. Generally, low-μM inhibition of Mtb DnaG by the anthracyclines was correlated with their low-μM minimum inhibitory concentrations. Aloe-emodin displayed threefold weaker potency than doxorubicin against Mtb DnaG and similar inhibition of Msm (but not Mtb) in the mid-μM range, whereas rhein (a close analog of aloe-emodin) and a di-glucosylated mitoxantrone derivative did not show significant inhibition of Mtb DnaG or antimycobacterial activity. Taken together, these observations strongly suggest that several clinically used anthracyclines and aloe-emodin target mycobacterial primase, setting the stage for a more extensive exploration of this enzyme as an antibacterial target.

2',6'-Dihalostyrylanilines, pyridines, and pyrimidines for the inhibition of the catalytic subunit of methionine S-adenosyltransferase-2

Inhibition of the catalytic subunit of the heterodimeric methionine S-adenosyl transferase-2 (MAT2A) with fluorinated N,N-dialkylaminostilbenes (FIDAS agents) offers a potential avenue for the treatment of liver and colorectal cancers where upregulation of this enzyme occurs. A study of structure–activity relationships led to the identification of the most active compounds as those with (1) either a 2,6-difluorostyryl or 2-chloro-6-fluorostyryl subunit, (2) either an N-methylamino or N,N-dimethylamino group attached in a para orientation relative to the 2,6-dihalostyryl subunit, and (3) either an N-methylaniline or a 2-(N,N-dimethylamino)pyridine ring. These modifications led to FIDAS agents that were active in the low nanomolar range, that formed water-soluble hydrochloride salts, and that possessed the desired property of not inhibiting the human hERG potassium ion channel at concentrations at which the FIDAS agents inhibit MAT2A. The active FIDAS agents may inhibit cancer cells through alterations of methylation reactions essential for cancer cell survival and growth.

Finding of widespread viral and bacterial revolution dsDNA translocation motors distinct from rotation motors by channel chirality and size

Background

Double-stranded DNA translocation is ubiquitous in living systems. Cell mitosis, bacterial binary fission, DNA replication or repair, homologous recombination, Holliday junction resolution, viral genome packaging and cell entry all involve biomotor-driven dsDNA translocation. Previously, biomotors have been primarily classified into linear and rotational motors. We recently discovered a third class of dsDNA translocation motors in Phi29 utilizing revolution mechanism without rotation. Analogically, the Earth rotates around its own axis every 24 hours, but revolves around the Sun every 365 days.

Results

Single-channel DNA translocation conductance assay combined with structure inspections of motor channels on bacteriophages P22, SPP1, HK97, T7, T4, Phi29, and other dsDNA translocation motors such as bacterial FtsK and eukaryotic mimiviruses or vaccinia viruses showed that revolution motor is widespread. The force generation mechanism for revolution motors is elucidated. Revolution motors can be differentiated from rotation motors by their channel size and chirality. Crystal structure inspection revealed that revolution motors commonly exhibit channel diameters larger than 3 nm, while rotation motors that rotate around one of the two separated DNA strands feature a diameter smaller than 2 nm. Phi29 revolution motor translocated double- and tetra-stranded DNA that occupied 32% and 64% of the narrowest channel cross-section, respectively, evidencing that revolution motors exhibit channel diameters significantly wider than the dsDNA. Left-handed oriented channels found in revolution motors drive the right-handed dsDNA via anti-chiral interaction, while right-handed channels observed in rotation motors drive the right-handed dsDNA via parallel threads. Tethering both the motor and the dsDNA distal-end of the revolution motor does not block DNA packaging, indicating that no rotation is required for motors of dsDNA phages, while a small-angle left-handed twist of dsDNA that is aligned with the channel could occur due to the conformational change of the phage motor channels from a left-handed configuration for DNA entry to a right-handed configuration for DNA ejection for host cell infection.

Conclusions

The revolution motor is widespread among biological systems, and can be distinguished from rotation motors by channel size and chirality. The revolution mechanism renders dsDNA void of coiling and torque during translocation of the lengthy helical chromosome, thus resulting in more efficient motor energy conversion.

A random sequential mechanism of aminoglycoside acetylation by Mycobacterium tuberculosis Eis protein

An important cause of bacterial resistance to aminoglycoside antibiotics is the enzymatic acetylation of their amino groups by acetyltransferases, which abolishes their binding to and inhibition of the bacterial ribosome. Enhanced intracellular survival (Eis) protein from Mycobacterium tuberculosis (Mt) is one of such acetyltransferases, whose upregulation was recently established as a cause of resistance to aminoglycosides in clinical cases of drug-resistant tuberculosis. The mechanism of aminoglycoside acetylation by MtEis is not completely understood. A systematic analysis of steady-state kinetics of acetylation of kanamycin A and neomycin B by Eis as a function of concentrations of these aminoglycosides and the acetyl donor, acetyl coenzyme A, reveals that MtEis employs a random-sequential bisubstrate mechanism of acetylation and yields the values of the kinetic parameters of this mechanism. The implications of these mechanistic properties for the design of inhibitors of Eis and other aminoglycoside acetyltransferases are discussed.

Chemical and structural insights into the regioversatility of the aminoglycoside acetyltransferase Eis

A recently discovered cause of tuberculosis resistance to a drug of last resort, the aminoglycoside kanamycin, results from modification of this drug by the enhanced intracellular survival (Eis) protein. Eis is a structurally and functionally unique acetyltransferase with an unusual capability of acetylating aminoglycosides at multiple positions. The extent of this regioversatility and its defining protein features are unclear. Herein, we determined the positions and order of acetylation of five aminoglycosides by NMR spectroscopy. This analysis revealed unprecedented acetylation of the 3''-amine of kanamycin, amikacin, and tobramycin, and the γ-amine of the 4-amino-2-hydroxybutyryl group of amikacin. A crystal structure of Eis in complex with coenzyme A and tobramycin revealed how tobramycin can be accommodated in the Eis active site in two binding modes, consistent with its diacetylation. These studies, describing chemical and structural details of acetylation, will guide future efforts towards designing aminoglycosides and Eis inhibitors to overcome resistance in tuberculosis.

Discovery of inhibitors of Bacillus anthracis primase DnaG

Primase DnaG is an essential bacterial enzyme that synthesizes short ribonucleotide primers required for chromosomal DNA replication. Inhibitors of DnaG can serve as leads for development of new antibacterials and biochemical probes. We recently developed a nonradioactive in vitro primase-pyrophosphatase assay to identify and analyze DnaG inhibitors. Application of this assay to DnaG from Bacillus anthracis (Ba DnaG), a dangerous pathogen, yielded several inhibitors, which include agents with DNA intercalating properties (doxorubicin and tilorone) as well as those that do not intercalate into DNA (suramin). A polyanionic agent and inhibitor of eukaryotic primases, suramin, identified by this assay as a low-micromolar Ba DnaG inhibitor, was recently shown to be also a low-micromolar inhibitor of Mycobacterium tuberculosis DnaG (Mtb DnaG). In contrast, another low-micromolar Ba DnaG inhibitor, tilorone, is much more potent against Ba DnaG than against Mtb DnaG, despite homology between these enzymes, suggesting that DnaG can be targeted selectively.

Quantifying additive interactions of the osmolyte proline with individual functional groups of proteins: comparisons with urea and glycine betaine, interpretation of m-values

To quantify interactions of the osmolyte l-proline with protein functional groups and predict their effects on protein processes, we use vapor pressure osmometry to determine chemical potential derivatives dμ2/dm3 = μ23, quantifying the preferential interactions of proline (component 3) with 21 solutes (component 2) selected to display different combinations of aliphatic or aromatic C, amide, carboxylate, phosphate or hydroxyl O, and amide or cationic N surface. Solubility data yield μ23 values for four less-soluble solutes. Values of μ23 are dissected using an ASA-based analysis to test the hypothesis of additivity and obtain α-values (proline interaction potentials) for these eight surface types and three inorganic ions. Values of μ23 predicted from these α-values agree with the experiment, demonstrating additivity. Molecular interpretation of α-values using the solute partitioning model yields partition coefficients (Kp) quantifying the local accumulation or exclusion of proline in the hydration water of each functional group. Interactions of proline with native protein surfaces and effects of proline on protein unfolding are predicted from α-values and ASA information and compared with experimental data, with results for glycine betaine and urea, and with predictions from transfer free energy analysis. We conclude that proline stabilizes proteins because of its unfavorable interactions with (exclusion from) amide oxygens and aliphatic hydrocarbon surfaces exposed in unfolding and that proline is an effective in vivo osmolyte because of the osmolality increase resulting from its unfavorable interactions with anionic (carboxylate and phosphate) and amide oxygens and aliphatic hydrocarbon groups on the surface of cytoplasmic proteins and nucleic acids.

Substrate specificity of MarP, a periplasmic protease required for resistance to acid and oxidative stress in Mycobacterium tuberculosis

The transmembrane serine protease MarP is important for pH homeostasis in Mycobacterium tuberculosis (Mtb). Previous structural studies revealed that MarP contains a chymotrypsin fold and a disulfide bond that stabilizes the protease active site in the substrate-bound conformation. Here, we determined that MarP is located in the Mtb periplasm and showed that this localization is essential for function. Using the recombinant protease domain of MarP, we identified its substrate specificity using two independent assays: positional-scanning synthetic combinatorial library profiling and multiplex substrate profiling by mass spectrometry. These methods revealed that MarP prefers bulky residues at P4, tryptophan or leucine at P2, arginine or hydrophobic residues at P1, and alanine or asparagine at P1'. Guided by these data, we designed fluorogenic peptide substrates and characterized the kinetic properties of MarP. Finally, we tested the impact of mutating MarP cysteine residues on the peptidolytic activity of recombinant MarP and its ability to complement phenotypes of Mtb ΔMarP. Taken together, our studies provide insight into the enzymatic properties of MarP, its substrate preference, and the importance of its transmembrane helices and disulfide bond.

DksA2, a zinc-independent structural analog of the transcription factor DksA

Transcription factor DksA contains a four-Cys Zn(2 +)-finger motif thought to be responsible for structural integrity and the relative disposition of its domains. Pseudomonas aeruginosa encodes an additional DksA paralog (DksA2) that is expressed selectively under Zn(2+) limitation. Although DksA2 does not bind Zn(2+), it complements the Escherichia coli dksA deletion and has similar effects on transcription in vitro. In this study, structural and biochemical analyses reveal that DksA2 has a similar fold, domain structure and RNA polymerase binding properties to those of the E. coli DksA despite the lack of the stabilizing metal ion.

A novel non-radioactive primase-pyrophosphatase activity assay and its application to the discovery of inhibitors of Mycobacterium tuberculosis primase DnaG

Bacterial DNA primase DnaG synthesizes RNA primers required for chromosomal DNA replication. Biochemical assays measuring primase activity have been limited to monitoring formation of radioactively labelled primers because of the intrinsically low catalytic efficiency of DnaG. Furthermore, DnaG is prone to aggregation and proteolytic degradation. These factors have impeded discovery of DnaG inhibitors by high-throughput screening (HTS). In this study, we expressed and purified the previously uncharacterized primase DnaG from Mycobacterium tuberculosis (Mtb DnaG). By coupling the activity of Mtb DnaG to that of another essential enzyme, inorganic pyrophosphatase from M. tuberculosis (Mtb PPiase), we developed the first non-radioactive primase–pyrophosphatase assay. An extensive optimization of the assay enabled its efficient use in HTS (Z′ = 0.7 in the 384-well format). HTS of 2560 small molecules to search for inhibitory compounds yielded several hits, including suramin, doxorubicin and ellagic acid. We demonstrate that these three compounds inhibit Mtb DnaG. Both suramin and doxorubicin are potent (low-µM) DNA- and nucleotide triphosphate-competitive priming inhibitors that interact with more than one site on Mtb DnaG. This novel assay should be applicable to other primases and inefficient DNA/RNA polymerases, facilitating their characterization and inhibitor discovery.

New insights into fluoroquinolone resistance in Mycobacterium tuberculosis: functional genetic analysis of gyrA and gyrB mutations

Fluoroquinolone antibiotics are among the most potent second-line drugs used for treatment of multidrug-resistant tuberculosis (MDR TB), and resistance to this class of antibiotics is one criterion for defining extensively drug resistant tuberculosis (XDR TB). Fluoroquinolone resistance in Mycobacterium tuberculosis has been associated with modification of the quinolone resistance determining region (QRDR) of gyrA. Recent studies suggest that amino acid substitutions in gyrB may also play a crucial role in resistance, but functional genetic studies of these mutations in M. tuberculosis are lacking. In this study, we examined twenty six mutations in gyrase genes gyrA (seven) and gyrB (nineteen) to determine the clinical relevance and role of these mutations in fluoroquinolone resistance. Transductants or clinical isolates harboring T80A, T80A+A90G, A90G, G247S and A384V gyrA mutations were susceptible to all fluoroquinolones tested. The A74S mutation conferred low-level resistance to moxifloxacin but susceptibility to ciprofloxacin, levofloxacin and ofloxacin, and the A74S+D94G double mutation conferred cross resistance to all the fluoroquinolones tested. Functional genetic analysis and structural modeling of gyrB suggest that M330I, V340L, R485C, D500A, D533A, A543T, A543V and T546M mutations are not sufficient to confer resistance as determined by agar proportion. Only three mutations, N538D, E540V and R485C+T539N, conferred resistance to all four fluoroquinolones in at least one genetic background. The D500H and D500N mutations conferred resistance only to levofloxacin and ofloxacin while N538K and E540D consistently conferred resistance to moxifloxacin only. Transductants and clinical isolates harboring T539N, T539P or N538T+T546M mutations exhibited low-level resistance to moxifloxacin only but not consistently. These findings indicate that certain mutations in gyrB confer fluoroquinolone resistance, but the level and pattern of resistance varies among the different mutations. The results from this study provide support for the inclusion of the QRDR of gyrB in molecular assays used to detect fluoroquinolone resistance in M. tuberculosis.

Aminoglycoside multiacetylating activity of the enhanced intracellular survival protein from Mycobacterium smegmatis and its inhibition

The enhanced intracellular survival (Eis) protein improves the survival of Mycobacterium smegmatis (Msm) in macrophages and functions as the acetyltransferase responsible for kanamycin A resistance, a hallmark of extensively drug-resistant (XDR) tuberculosis, in a large number of Mycobacterium tuberculosis (Mtb) clinical isolates. We recently demonstrated that Eis from Mtb (Eis_Mtb) efficiently multiacetylates a variety of aminoglycoside (AG) antibiotics. Here, to gain insight into the origin of substrate selectivity of AG multiacetylation by Eis, we analyzed AG acetylation by Eis_Msm, investigated its inhibition, and compared these functions to those of Eis_Mtb. Even though for several AGs the multiacetylation properties of Eis_Msm and Eis_Mtb are similar, there are three major differences. (i) Eis_Msm diacetylates apramycin, a conformationally constrained AG, which Eis_Mtb cannot modify. (ii) Eis_Msm triacetylates paromomycin, which can be only diacetylated by Eis_Mtb. (iii) Eis_Msm only monoacetylates hygromycin, a structurally unique AG that is diacetylated by Eis_Mtb. Several nonconserved amino acid residues lining the AG-binding pocket of Eis are likely responsible for these differences between the two Eis homologues. Specifically, we propose that because the AG-binding pocket of Eis_Msm is more open than that of Eis_Mtb, it accommodates apramycin for acetylation in Eis_Msm, but not in Eis_Mtb. We also demonstrate that inhibitors of Eis_Mtb that we recently discovered can inhibit Eis_Msm activity. These observations help define the structural origins of substrate preference among Eis homologues and suggest that Eis_Mtb inhibitors may be applied against all pathogenic mycobacteria to overcome AG resistance caused by Eis upregulation.

Enumeration of one-dimensional crystal structures obtained from a minimum of diffraction intensities

A central problem in crystallography is crystal structure determination directly from diffraction intensities. For structures of small molecules, this problem has been addressed by probabilistic direct methods that allow one to obtain the structure coordinates with a high degree of certainty given a sufficiently large set of intensities. In contrast, deterministic algebraic methods that could guarantee a solution and may be applicable to macromolecules have not yet emerged. In this study a basic algebraic question is posed: how many crystal structures can be obtained from a given set of intensities? Recently, by using a new origin definition and the method of elementary symmetrical polynomials, all small (N ≤ 4 atoms) one-dimensional crystal structures that could be obtained from the minimum set of N - 1 lowest-resolution intensities were enumerated. Here, by using methods of modern algebraic geometry the maximum number of one-dimensional crystal structures that can be determined from the minimum set of intensities for N > 4 is obtained. It is demonstrated that this ambiguity increases exponentially with the increasing number of atoms in the structure N (~4(N)/N(3/2) for N >> 1) and includes non-homometric structures. Therefore, a minimum set of intensities, even in principle, is insufficient for structure determination for all but very small structures.

The structural basis for substrate versatility of chloramphenicol acetyltransferase CATI

Novel antibiotics are needed to overcome the challenge of continually evolving bacterial resistance. This has led to a renewed interest in mechanistic studies of once popular antibiotics like chloramphenicol (CAM). Chloramphenicol acetyltransferases (CATs) are enzymes that covalently modify CAM, rendering it inactive against its target, the ribosome, and thereby causing resistance to CAM. Of the three major types of CAT (CAT(I-III)), the CAM-specific CAT(III) has been studied extensively. Much less is known about another clinically important type, CAT(I). In addition to inactivating CAM and unlike CAT(III), CAT(I) confers resistance to a structurally distinct antibiotic, fusidic acid. The origin of the broader substrate specificity of CAT(I) has not been fully elucidated. To understand the substrate binding features of CAT(I), its crystal structures in the unbound (apo) and CAM-bound forms were determined. The analysis of these and previously determined CAT(I)-FA and CAT(III)-CAM structures revealed interactions responsible for CAT(I) binding to its substrates and clarified the broader substrate preference of CAT(I) compared to that of CAT(III).

Structural and thermodynamic signatures of DNA recognition by Mycobacterium tuberculosis DnaA

An essential protein, DnaA, binds to 9-bp DNA sites within the origin of replication oriC. These binding events are prerequisite to forming an enigmatic nucleoprotein scaffold that initiates replication. The number, sequences, positions, and orientations of these short DNA sites, or DnaA boxes, within the oriCs of different bacteria vary considerably. To investigate features of DnaA boxes that are important for binding Mycobacterium tuberculosis DnaA (MtDnaA), we have determined the crystal structures of the DNA binding domain (DBD) of MtDnaA bound to a cognate MtDnaA-box (at 2.0 Å resolution) and to a consensus Escherichia coli DnaA-box (at 2.3 Å). These structures, complemented by calorimetric equilibrium binding studies of MtDnaA DBD in a series of DnaA-box variants, reveal the main determinants of DNA recognition and establish the [T/C][T/A][G/A]TCCACA sequence as a high-affinity MtDnaA-box. Bioinformatic and calorimetric analyses indicate that DnaA-box sequences in mycobacterial oriCs generally differ from the optimal binding sequence. This sequence variation occurs commonly at the first 2 bp, making an in vivo mycobacterial DnaA-box effectively a 7-mer and not a 9-mer. We demonstrate that the decrease in the affinity of these MtDnaA-box variants for MtDnaA DBD relative to that of the highest-affinity box TTGTCCACA is less than 10-fold. The understanding of DnaA-box recognition by MtDnaA and E. coli DnaA enables one to map DnaA-box sequences in the genomes of M. tuberculosis and other eubacteria.

Conditional ambiguity of one-dimensional crystal structures determined from a minimum of diffraction intensity data

When the number of intensities greatly exceeds the number of unknown atomic coordinates, the problem of obtaining a crystal structure from the intensities is overdetermined and, for a sufficiently small structure, a chemically meaningful solution can be found by direct methods. A difficulty in determining a structure has been historically attributed to the non-uniqueness of such a structure owing to multiple, or homometric, structures that yield the same set of intensities. The number of homometric structures has not been rigorously analyzed owing to the complexity of this problem. By using the method of elementary symmetric polynomials with a new origin definition, one-dimensional crystal structures of a small number of identical atoms (N < 5), determined from a minimum (N - 1) of the lowest-resolution intensities, are enumerated. It is demonstrated that such a structure is unique for N ≤ 3. Interestingly, for N = 4, the structure can be determined either uniquely or twofold ambiguously, depending on the intensity values. These results suggest that, even for larger structures, a minimum set of (or not many more) accurately measured intensities can yield a unique structure.

Structural insight into serine protease Rv3671c that Protects M. tuberculosis from oxidative and acidic stress

Rv3671c, a putative serine protease, is crucial for persistence of Mycobacterium tuberculosis in the hostile environment of the phagosome. We show that Rv3671c is required for M. tuberculosis resistance to oxidative stress in addition to its role in protection from acidification. Structural and biochemical analyses demonstrate that the periplasmic domain of Rv3671c is a functional serine protease of the chymotrypsin family and, remarkably, that its activity increases on oxidation. High-resolution crystal structures of this protease in an active strained state and in an inactive relaxed state reveal that a solvent-exposed disulfide bond controls the protease activity by constraining two distant regions of Rv3671c and stabilizing it in the catalytically active conformation. In vitro biochemical studies confirm that activation of the protease in an oxidative environment is dependent on this reversible disulfide bond. These results suggest that the disulfide bond modulates activity of Rv3671c depending on the oxidative environment in vivo.

A new scaffold of an old protein fold ensures binding to the bisintercalator thiocoraline

Thiocoraline is a thiodepsipeptide with potent antitumor activity. TioX, a protein with an unidentified function, is encoded by a gene of the thiocoraline biosynthetic gene cluster. The crystal structure of the full-length TioX protein at 2.15 A resolution reveals that TioX protomer shares an ancient betaalphabetabetabeta fold motif with glyoxalase I and bleomycin resistance protein families, despite a very low sequence homology. Intriguingly, four TioX monomers form a unique 2-fold symmetric tetrameric assembly that is stabilized by four intermolecular disulfide bonds formed cyclically between Cys60 and Cys66 of adjacent monomers. The arrangement of two of the four monomers in the TioX tetramer is analogous to that in dimeric bleomycin resistance proteins. This analogy indicates that this novel higher-order structural scaffold of TioX may have evolved to bind thiocoraline. Our equilibrium titration studies demonstrate the binding of a thiocoraline chromophore analog, quinaldic acid, to TioX, thereby substantiating this model. Furthermore, a strain of Streptomyces albus containing an exogenous thiocoraline gene cluster devoid of functional tioX maintains thiocoraline production, albeit with a lower yield. Taken together, these observations rule out a direct enzymatic function of TioX and suggest that TioX is involved in thiocoraline resistance or secretion.

The XPA-binding domain of ERCC1 is required for nucleotide excision repair but not other DNA repair pathways

The endonuclease ERCC1-XPF incises the damaged strand of DNA 5' to a lesion during nucleotide excision repair (NER) and has additional, poorly characterized functions in interstrand cross-link repair, double-strand break repair, and homologous recombination. XPA, another key factor in NER, interacts with ERCC1 and recruits it to sites of damage. We identified ERCC1 residues that are critical for the interaction with XPA and assessed their importance for NER in vitro and in vivo. Mutation of two conserved residues (Asn-110 and Tyr-145) located in the XPA-binding site of ERCC1 dramatically affected NER but not nuclease activity on model DNA substrates. In ERCC1-deficient cells expressing ERCC1(N110A/Y145A), the nuclease was not recruited to sites of UV damage. The repair of UV-induced (6-4)photoproducts was severely impaired in these cells, and they were hypersensitive to UV irradiation. Remarkably, the ERCC1(N110A/Y145A) protein rescues the sensitivity of ERCC1-deficient cells to cross-linking agents. Our studies suggest that ERCC1-XPF engages in different repair pathways through specific protein-protein interactions and that these functions can be separated through the selective disruption of these interactions. We discuss the impact of these findings for understanding how ERCC1 contributes to resistance of tumor cells to therapeutic agents such as cisplatin.

The tail of KdsC: conformational changes control the activity of a haloacid dehalogenase superfamily phosphatase

The phosphatase KdsC cleaves 3-deoxy-D-manno-octulosonate 8-phosphate to generate a molecule of inorganic phosphate and Kdo. Kdo is an essential component of the lipopolysaccharide envelope in Gram-negative bacteria. Because lipopolysaccharide is an important determinant of bacterial resistance and toxicity, KdsC is a potential target for novel antibacterial agents. KdsC belongs to the broad haloacid dehalogenase superfamily. In haloacid dehalogenase superfamily enzymes, substrate specificity and catalytic efficiency are generally dictated by a fold feature called the cap domain. It is therefore not clear why KdsC, which lacks a cap domain, is catalytically efficient and highly specific to 3-deoxy-D-manno-octulosonate 8-phosphate. Here, we present a set of seven structures of tetrameric Escherichia coli KdsC (ranging from 1.4 to 3.06 A in resolution) that model different intermediate states in its catalytic mechanism. A crystal structure of product-bound E. coli KdsC shows how the interface between adjacent monomers defines the active site pocket. Kdo is engaged in a network of polar and nonpolar interactions with residues at this interface, which explains substrate specificity. Furthermore, this structural and kinetic analysis strongly suggests that the binding of the flexible C-terminal region (tail) to the active site makes KdsC catalytically efficient by facilitating product release.

A nonradioactive high-throughput assay for screening and characterization of adenylation domains for nonribosomal peptide combinatorial biosynthesis

Adenylation domains are critical enzymes that dictate the identity of the amino acid building blocks to be incorporated during nonribosomal peptide (NRP) biosynthesis. NRPs display a wide range of biological activities and are some of the most important drugs currently used in clinics. Traditionally, activity of adenylation domains has been measured by radioactive ATP-[32P]pyrophosphate (PP(i)) exchange assays. To identify adenylation domains for future combinatorial production of novel NRPs as potential drugs, we report a convenient high-throughput nonradioactive method to measure activity of these enzymes. In our assay, malachite green is used to measure orthophosphate (P(i)) concentrations after degradation by inorganic pyrophosphatase of the PP(i) released during aminoacyl-AMP formation by action of the adenylation domains. The assay is quantitative, accurate, and robust, and it can be performed in 96- and 384-well plate formats. The performance of our assay was tested by using NcpB-A(4), one of the seven adenylation domains involved in nostocyclopeptide biosynthesis. The kinetics of pyrophosphate release monitored by this method are much slower than those measured by a traditional ATP-[32P]PP(i) exchange assay. This observation indicates that the formation of the adenylated amino acid and its release are the rate-limiting steps during the catalytic turnover.

Hexameric ring structure of the N-terminal domain of Mycobacterium tuberculosis DnaB helicase

Hexameric DnaB helicase unwinds the DNA double helix during replication of genetic material in bacteria. DnaB is an essential bacterial protein; therefore, it is an important potential target for antibacterial drug discovery. We report a crystal structure of the N-terminal region of DnaB from the pathogen Mycobacterium tuberculosis (MtDnaBn), determined at 2.0 A resolution. This structure provides atomic resolution details of formation of the hexameric ring of DnaB by two distinct interfaces. An extensive hydrophobic interface stabilizes a dimer of MtDnaBn by forming a four-helix bundle. The other, less extensive, interface is formed between the dimers, connecting three of them into a hexameric ring. On the basis of crystal packing interactions between MtDnaBn rings, we suggest a model of a helicase-primase complex that explains previously observed effects of DnaB mutations on DNA priming.

Domain-swapped dimerization of the HIV-1 capsid C-terminal domain

Assembly of the HIV and other retroviruses is primarily driven by the oligomerization of the Gag polyprotein, the major viral structural protein capable of forming virus-like particles even in the absence of all other virally encoded components. Several critical determinants of Gag oligomerization are located in the C-terminal domain of the capsid protein (CA-CTD), which encompasses the most conserved segment in the highly variable Gag protein called the major homology region (MHR). The CA-CTD is thought to function as a dimerization module, although the existing model of CA-CTD dimerization does not readily explain why the conserved residues of the MHR are essential for retroviral assembly. Here we describe an x-ray structure of a distinct domain-swapped variant of the HIV-1 CA-CTD dimer stabilized by a single amino acid deletion. In the domain-swapped structure, the MHR-containing segment forms a major part of the dimerization interface, providing a structural mechanism for the enigmatic function of the MHR in HIV assembly. Our observations suggest that swapping of the MHR segments of adjacent Gag molecules may be a critical intermediate in retroviral assembly.

A flexible interface between DNA ligase and PCNA supports conformational switching and efficient ligation of DNA

DNA sliding clamps encircle DNA and provide binding sites for many DNA-processing enzymes. However, it is largely unknown how sliding clamps like proliferating cell nuclear antigen (PCNA) coordinate multistep DNA transactions. We have determined structures of Sulfolobus solfataricus DNA ligase and heterotrimeric PCNA separately by X-ray diffraction and in complex by small-angle X-ray scattering (SAXS). Three distinct PCNA subunits assemble into a protein ring resembling the homotrimeric PCNA of humans but with three unique protein-binding sites. In the absence of nicked DNA, the Sulfolobus solfataricus DNA ligase has an open, extended conformation. When complexed with heterotrimeric PCNA, the DNA ligase binds to the PCNA3 subunit and ligase retains an open, extended conformation. A closed, ring-shaped conformation of ligase catalyzes a DNA end-joining reaction that is strongly stimulated by PCNA. This open-to-closed switch in the conformation of DNA ligase is accommodated by a malleable interface with PCNA that serves as an efficient platform for DNA ligation.

Crystal structure and DNA binding functions of ERCC1, a subunit of the DNA structure-specific endonuclease XPF-ERCC1

Human XPF-ERCC1 is a DNA endonuclease that incises a damaged DNA strand on the 5' side of a lesion during nucleotide excision repair and has additional role(s) in homologous recombination and DNA interstrand crosslink repair. We show that a truncated form of XPF lacking the N-terminal helicase-like domain in complex with ERCC1 exhibits a structure-specific endonuclease activity with similar specificity to that of full-length XPF-ERCC1. Two domains of ERCC1, a central domain and a C-terminal tandem helix-hairpin-helix (HhH2) dimerization domain, bind to ssDNA. The central domain of ERCC1 binds ssDNA/dsDNA junctions with a defined polarity, preferring a 5' single-stranded overhang. The XPF-ERCC1 HhH2 domain heterodimer contains two independent ssDNA-binding surfaces, which are revealed by a crystal structure of the protein complex. A crystal structure of the central domain of ERCC1 shows its fold is strikingly similar to that of the nuclease domains of the archaeal Mus81/XPF homologs, despite very low sequence homology. A groove lined with basic and aromatic residues on the surface of ERCC1 has apparently been adapted to interact with ssDNA. On the basis of these crystallographic and biochemical studies, we propose a model in which XPF-ERCC1 recognizes a branched DNA substrate by binding the two ssDNA arms with the two HhH2 domains of XPF and ERCC1 and by binding the 5'-ssDNA arm with the central domain of ERCC1.

Organized unidirectional waves of ATP hydrolysis within a RecA filament

The RecA protein forms nucleoprotein filaments on DNA, and individual monomers within the filaments hydrolyze ATP. Assembly and disassembly of filaments are both unidirectional, occurring on opposite filament ends, with disassembly requiring ATP hydrolysis. When filaments form on duplex DNA, RecA protein exhibits a functional state comparable to the state observed during active DNA strand exchange. RecA filament state was monitored with a coupled spectrophotometric assay for ATP hydrolysis, with changes fit to a mathematical model for filament disassembly. At 37 degrees C, monomers within the RecA-double-stranded DNA (dsDNA) filaments hydrolyze ATP with an observed k(cat) of 20.8 +/- 1.5 min(-1). Under the same conditions, the rate of end-dependent filament disassembly (k(off)) is 123 +/- 16 monomers per minute per filament end. This rate of disassembly requires a tight coupling of the ATP hydrolytic cycles of adjacent RecA monomers. The relationship of k(cat) to k(off) infers a filament state in which waves of ATP hydrolysis move unidirectionally through RecA filaments on dsDNA, with successive waves occurring at intervals of approximately six monomers. The waves move nearly synchronously, each one transiting from one monomer to the next every 0.5 s. The results reflect an organization of the ATPase activity that is unique in filamentous systems, and could be linked to a RecA motor function.

Kinetic studies and structural models of the association of E. coli sigma(70) RNA polymerase with the lambdaP(R) promoter: large scale conformational changes in forming the kinetically significant intermediates

The kinetics of interaction of Esigma(70) RNA polymerase (R) with the lambdaP(R) promoter (P) were investigated by filter binding over a broad range of temperatures (7.3-42 degrees C) and concentrations of RNA polymerase (1-123 nM) in large excess over promoter DNA. Under all conditions examined, the kinetics of formation of competitor-resistant complexes (I(2), RP(o)) are single-exponential with first order rate constant beta(CR). Interpretation of the polymerase concentration dependence of beta(CR) in terms of the three step mechanism of open complex formation yields the equilibrium constant K(1) for formation of the first kinetically significant intermediate (I(1)) and the forward rate constant (k(2)) for the conformational change converting I(1) to the second kinetically significant intermediate I(2): R + P-->(K(1))<--I(1)(k(2))-->I(2). Use of rapid quench mixing allows K(1) and k(2) to be individually determined over the entire temperature range investigated, previously not possible at this promoter using manual mixing. Given the large (>60 bp) interface formed in I(1), its relatively small binding constant K(1) at 37 degrees C at this [salt] (approximately 6 x 10(6) M(-1)) strongly argues that binding free energy is used to drive large-scale structural changes in polymerase and/or promoter DNA or other coupled processes. Evidence for coupling of protein folding is provided by the large and negative activation heat capacity of k(a)[DeltaC(o,++)(a)= -1.5(+/-0.2)kcal K(-1)], now shown to originate directly from formation of I(1) [DeltaC(o)(1)= -1.4(+/-0.3)kcal K(-1)] rather than from the formation of I(2) as previously proposed. The isomerization I(1)-->I(2) exhibits relatively slow kinetics and has a very large temperature-independent Arrhenius activation energy [E(act)(2)= 34(+/-2)kcal]. This kinetic signature suggests that formation of the transition state (I(1)-I(2)++ involves large conformational changes dominated by changes in the exposure of polar and/or charged surface to water. Structural and biochemical data lead to the following hypotheses to interpret these results. We propose that formation of I(1) involves coupled folding of unstructured regions of polymerase (beta, beta' and sigma(70)) and bending of promoter DNA (in the -10 region). We propose that interactions with region 2 of sigma(70) and possibly domain 1 of beta induce a kink at the -11/-12 base pairs of the lambdaP(R) promoter which places the downstream DNA (-5 to +20) in the jaws of the beta and beta' subunits of polymerase in I(1). These early interactions of beta and beta' with the DNA downstream of position -5 trigger jaw closing (with coupled folding) and subsequent steps of DNA opening.

Novel computer program for fast exact calculation of accessible and molecular surface areas and average surface curvature

New computer programs, SurfRace and FastSurf, perform fast calculations of the solvent accessible and molecular (solvent excluded) surface areas of macromolecules. Program SurfRace also calculates the areas of cavities inaccessible from the outside. We introduce the definition of average curvature of molecular surface and calculate average molecular surface curvatures for each atom in a structure. All surface area and curvature calculations are analytic and therefore yield exact values of these quantities. High calculation speed of this software is achieved primarily by avoiding computationally expensive mathematical procedures wherever possible and by efficient handling of surface data structures. The programs are written initially in the language C for PCs running Windows 2000/98/NT, but their code is portable to other platforms with only minor changes in input-output procedures. The algorithm is robust and does not ignore either multiplicity or degeneracy of atomic overlaps. Fast, memory-efficient and robust execution make this software attractive for applications both in computationally expensive energy minimization algorithms, such as docking or molecular dynamics simulations, and in stand-alone surface area and curvature calculations.

Interaction of RNA polymerase with forked DNA: evidence for two kinetically significant intermediates on the pathway to the final complex

RNA polymerase forms competitor-resistant complexes with "forked DNA" templates that are double-stranded from the -35 promoter region through the first base pair of the -10 region, with an additional unpaired A at the 3' end of the nontemplate strand. These types of substrates were introduced recently as model templates for the study of DNA-protein interactions occurring in the early stages of the formation of RNA polymerase-promoter open complexes. We have performed kinetic and equilibrium measurements of interactions of wild-type and mutant RNA polymerases bearing substitutions in the sigma(70) initiation factor, with forked DNA of wild-type and mutant sequence. Our observations reveal that formation of a competitor-resistant complex between RNA polymerase and forked DNA, similar to the formation of open complexes at promoters, is a multistep process, and some of the sequentially formed intermediates along the two pathways share common properties. This work establishes, for the forked template, progression through these intermediates in the absence of downstream DNA and validates the use of forked DNA to determine the effects of changes in promoter or RNA polymerase sequence on the process of open complex formation.

Asymptotic solution of the cylindrical nonlinear Poisson-Boltzmann equation at low salt concentration: analytic expressions for surface potential and preferential interaction coefficient

The analytic solution to the nonlinear Poisson-Boltzmann equation describing the ion distributions surrounding a nucleic acid or other cylindrical polyions as a function of polyion structural quantities and salt concentration ([salt]) has been sought for more than 80 years to predict the effect of these quantities on the thermodynamics of polyion processes. Here we report an accurate asymptotic solution of the cylindrical nonlinear Poisson-Boltzmann equation at low to moderate concentration of a symmetrical electrolyte (< or = 0.1 M 1:1 salt). The approximate solution for the potential is derived as an asymptotic series in the small parameter var epsilon(-1), where var epsilon identical with kappa(-1)/a, the ratio of the Debye length (kappa(-1)) to the polyion radius (a). From the potential at the polyion surface, we obtain the coulombic contribution to the salt-polyelectrolyte preferential interaction (Donnan) coefficient (Gamma(u)coul) per polyion charge at any reduced axial charge density xi. Gamma(u)coul is the sum of the previously recognized low-salt limiting value and a salt-dependent contribution, analytically derived here in the range of low-salt concentrations. As an example of the application of this solution, we obtain an analytic expression for the derivative of the midpoint temperature of a nucleic acid conformational transition with respect to the logarithm of salt concentration (dT(m)/d ln[salt]) in terms of [salt] and nucleic acid structural quantities. This expression explains the experimental observation that this derivative is relatively independent of salt concentration but deviates significantly from its low-salt limiting value in the range 0.01-0.1 M.

Analytic binding isotherms describing competitive interactions of a protein ligand with specific and nonspecific sites on the same DNA oligomer

Many studies of specific protein-nucleic acid binding use short oligonucleotides or restriction fragments, in part to minimize the potential for nonspecific binding of the protein. However, when the specificity ratio is low, multiple nonspecifically bound proteins may occupy the region of DNA corresponding to one specific site; this situation was encountered in our recent calorimetric study of binding of integration host factor (IHF) protein to its specific 34-bp H' DNA site. Here, beginning from the analytical McGhee and von Hippel infinite-lattice nonspecific binding isotherm, we derive a novel analytic isotherm for nonspecific binding of a ligand to a finite lattice. This isotherm is an excellent approximation to the exact factorial-based Epstein finite lattice isotherm even for short lattices and therefore is of great practical significance for analysis of experimental data and for analytic theory. Using this isotherm, we develop an analytic treatment of the competition between specific and nonspecific binding of a large ligand to the same finite lattice (i.e., DNA oligomer) containing one specific and multiple overlapping nonspecific binding sites. Analysis of calorimetric data for IHF-H' DNA binding using this treatment yields enthalpies and binding constants for both specific and nonspecific binding and the nonspecific site size. This novel analysis demonstrates the potential contribution of nonspecific binding to the observed thermodynamics of specific binding, even with very short DNA oligomers, and the need for reverse (constant protein) titrations or titrations with nonspecific DNA to resolve specific and nonspecific contributions. The competition treatment is useful in analyzing low-specificity systems, including those where specificity is weakened by mutations or the absence of specificity factors.

Specific and non-specific interactions of integration host factor with DNA: thermodynamic evidence for disruption of multiple IHF surface salt-bridges coupled to DNA binding

Site-specific DNA binding of architectural protein integration host factor (IHF) is involved in formation of functional multiprotein-DNA assemblies in Escherichia coli, while non-specific binding of IHF and other histone-like proteins serves to structure the nucleoid. Here, we report an isothermal titration calorimetry study of the thermodynamics of binding IHF to a 34 bp fragment composed entirely of the specific H' site from lambda-phage DNA. At low to moderate [K(+)] (60-100 mM), strong competition is observed between specific and non-specific binding as a result of a low specificity ratio (approximately 10(2)) and a very small non-specific site size. In this [K(+)] range, both specific and non-specific binding are enthalpy-driven, with large negative enthalpy, entropy and heat capacity changes and binding constants that are insensitive to [K(+)]. Above 100 mM K(+), only specific binding is observed, and both the binding constant and the magnitudes of enthalpy, entropy and heat capacity changes all decrease strongly with increasing [K(+)]. When interpreted in the context of the structure of the specific complex, the thermodynamics provide compelling evidence for a previously unrecognized design principle by which proteins that form extensive binding interfaces with nucleic acids control binding constants, binding site sizes and effects of temperature and ion concentrations on stability and specificity. We propose that up to 22 of the 23 IHF cationic side-chains that are located within 6 A of DNA phosphate oxygen atoms in the complex, are masked in the absence of DNA by pairing with anionic carboxylate groups in intramolecular salt-bridges (dehydrated ion-pairs). These salt-bridges increase in stability with increasing temperature and decreasing [K(+)]. To explain the unusual thermodynamics of IHF-DNA interactions, we propose that both specific and non-specific binding at low [K(+)] require disruption of salt-bridges (as many as 18 for specific binding) whereupon many of the unmasked charged groups hydrate and the cationic groups interact with DNA. From structural or thermodynamic parallels with IHF, we propose that large-scale coupling of disruption of protein salt-bridges to DNA binding is significant for other large-interface DNA wrapping proteins including the nucleosome, lac repressor core tetramer, RNA polymerase core protein, HU and SSB.

Wrapping of flanking non-operator DNA in lac repressor-operator complexes: implications for DNA looping

In our studies of lac repressor tetramer (T)-lac operator (O) interactions, we observed that the presence of extended regions of non-operator DNA flanking a single lac operator sequence embedded in plasmid DNA produced large and unusual cooperative and anticooperative effects on binding constants (Kobs) and their salt concentration dependences for the formation of 1:1 (TO) and especially 1:2 (TO2) complexes. To explore the origin of this striking behavior we report and analyze binding data on 1:1 (TO) and 1:2 (TO2) complexes between repressor and a single O(sym) operator embedded in 40 bp, 101 bp, and 2514 bp DNA, over very wide ranges of [salt]. We find large interrelated effects of flanking DNA length and [salt] on binding constants (K(TO)obs, K(TO2)obs) and on their [salt]-derivatives, and quantify these effects in terms of the free energy contributions of two wrapping modes, designated local and global. Both local and global wrapping of flanking DNA occur to an increasing extent as [salt] decreases. Global wrapping of plasmid-length DNA is extraordinarily dependent on [salt]. We propose that global wrapping is driven at low salt concentration by the polyelectrolyte effect, and involves a very large number (>/similar 20) of coulombic interactions between DNA phosphates and positively charged groups on lac repressor. Coulombic interactions in the global wrap must involve both the core and the second DNA-binding domain of lac repressor, and result in a complex which is looped by DNA wrapping. The non-coulombic contribution to the free energy of global wrapping is highly unfavorable ( approximately +30-50 kcal mol(-1)), which presumably results from a significant extent of DNA distortion and/or entropic constraints. We propose a structural model for global wrapping, and consider its implications for looping of intervening non-operator DNA in forming a complex between a tetrameric repressor (LacI) and one multi-operator DNA molecule in vivo and in vitro. The existence of DNA wrapping in LacI-DNA interactions motivates the proposal that most if not all DNA binding proteins may have evolved the capability to wrap and thereby organize flanking regions of DNA.

Quantitative analysis of the kinetics of end-dependent disassembly of RecA filaments from ssDNA

On linear single-stranded DNA, RecA filaments assemble and disassemble in the 5' to 3' direction. Monomers (or other units) associate at one end and dissociate from the other. ATP hydrolysis occurs throughout the filament. Dissociation can result when ATP is hydrolyzed by the monomer at the disassembly end. We have developed a comprehensive model for the end-dependent filament disassembly process. The model accounts not only for disassembly, but also for the limited reassembly that occurs as DNA is vacated by disassembling filaments. The overall process can be monitored quantitatively by following the resulting decline in DNA-dependent ATP hydrolysis. The rate of disassembly is highly pH dependent, being negligible at pH 6 and reaching a maximum at pH values above 7. 5. The rate of disassembly is not significantly affected by the concentration of free RecA protein within the experimental uncertainty. For filaments on single-stranded DNA, the monomer kcat for ATP hydrolysis is 30 min-1, and disassembly proceeds at a maximum rate of 60-70 monomers per minute per filament end. The latter rate is that predicted if the ATP hydrolytic cycles of adjacent monomers are not coupled in any way.

General method of analysis of kinetic equations for multistep reversible mechanisms in the single-exponential regime: application to kinetics of open complex formation between Esigma70 RNA polymerase and lambdaP(R) promoter DNA

A novel analytical method based on the exact solution of equations of kinetics of unbranched first- and pseudofirst-order mechanisms is developed for application to the process of Esigma70 RNA polymerase (R)-lambdaPR promoter (P) open complex formation, which is described by the minimal three-step mechanism with two kinetically significant intermediates (I1, I2), [equation: see text], where the final product is an open complex RPo. The kinetics of reversible and irreversible association (pseudofirst order, [R] >> [P]) to form long-lived complexes (RPo and I2) and the kinetics of dissociation of long-lived complexes both exhibit single exponential behavior. In this situation, the analytical method provides explicit expressions relating observed rate constants to the microscopic rate constants of mechanism steps without use of rapid equilibrium or steady-state approximations, and thereby provides a basis for interpreting the composite rate constants of association (ka), isomerization (ki), and dissociation (kd) obtained from experiment for this or any other sequential mechanism of any number of steps. In subsequent papers, we apply this formalism to analyze kinetic data obtained in the reversible and irreversible binding regimes of Esigma70 RNA polymerase (R)-lambdaP(R) promoter (P) open complex formation.

DNA footprints of the two kinetically significant intermediates in formation of an RNA polymerase-promoter open complex: evidence that interactions with start site and downstream DNA induce sequential conformational changes in polymerase and DNA

Kinetic studies of formation and dissociation of open-promoter complexes (RPo) involving Esigma70 RNA polymerase (R) and the lambdaPR promoter (P) demonstrate the existence of two kinetically significant intermediates, designated I1 and I2, and facilitate the choice of conditions under which each accumulates. For such conditions, we report the results of equilibrium and transient DNase I and KMnO4 footprinting studies which characterize I1 and I2. At 0 degreesC, where extrapolation of equilibrium data indicates I1 is the dominant complex, DNA bases in the vicinity of the transcription start site (+1) do not react with KMnO4, indicating that this region is closed in I1. However, the DNA backbone in I1 is extensively protected from DNase I cleavage; the DNase I footprint extends approximately 30 bases downstream and at least approximately 40 bases upstream from the start site. I1 has a short lifetime (</=15 seconds), based on its sensitivity to competition with heparin. Shortly after a temperature downshift from 37 degreesC to 0 degreesC, in the time-range where we conclude that the dominant, transiently accumulated complex is I2, DNase I and KMnO4 footprinting reveal a complex with a closed-start site and an extended DNase I footprint like that of I1. However, unlike I1, I2 is insensitive to heparin competition and has a much longer dissociation lifetime at 0 degreesC. Based on footprinting, kinetic and thermodynamic studies, we conclude that in the short-lived intermediate I1 the promoter start site and downstream region are bound in a cleft defined by the open clamp-like jaws of Esigma70. We propose that binding of the start site and downstream DNA in this cleft triggers massive, relatively slow conformational changes which likely include RNA polymerase jaw closing with coupled folding. These proposed conformational changes occur prior to opening of the promoter start site region, and are responsible for the much longer lifetime of I2. Closing of the jaws of polymerase around the downstream region of promoter DNA appears to trigger opening of the start site region. From a quantitative analysis of the biphasic decay of KMnO4 reactivity of RPo at 0 degreesC, we obtain the equilibrium constant K3 for the conversion of I2 to RPo and the rate constant k-2 for the conversion of I2 to I1 (i.e. jaw opening). These quantitative results were previously unavailable at any temperature, and are necessary for the dissection of dissociation kinetic data at higher temperatures.

Quantitative analysis of multiple-hit footprinting studies to characterize DNA conformational changes in protein-DNA complexes: application to DNA opening by Esigma70 RNA polymerase

Formation of many site-specific protein-nucleic acid complexes involves sequential conformational changes subsequent to initial binding which create functionally active assemblies. Characterization of population distributions and structural characteristics of intermediate and product conformations is necessary to understand both the mechanisms and the thermodynamics of these processes. For these purposes, here we develop the quantitative method of multiple hit footprinting (MHF), where chemical or enzymatic probing is performed as a function of either concentrations of the footprinting agent and/or time of exposure to it, in the multiple hit regime where many of the population or subpopulation of reactive DNA molecules are modified at more than one site. Properly controlled MHF experiments yield both the population distribution of different conformers and reactivity rate constants of the footprinting agent at all reactive positions in each conformer, which may be interpreted in terms of the accessibility of the site or the local concentration of the reagent. MHF experiments are particularly well-suited for dissecting effects at sites where unbound DNA is non-reactive and bound DNA is reactive with base-specific probes (e.g. KMnO4, DMS). We suggest that this method will also be applicable to analysis of enhancements in reactivity of other footprinting agents (e.g. DNase I, HO.). To demonstrate the utility of the MHF analysis, we quantify fragment distributions and individual site reactivities from multiple-hit KMnO4 footprinting of the non-template strand of Esigma70 RNA polymerase-lambdaPR promoter DNA complexes populated at binding equilibrium at 37 degreesC and transiently populated at a fixed time after a temperature downshift from 37 degreesC to 0 degreesC. For this system, a MHF analysis directly addresses the following questions: (i) what fraction of the population of promoter DNA molecules is open in the vicinity of the transcription start site (RPo) both at 37 degreesC and (transiently) after a downshift to 0 degreesC; (ii) does opening of the start site region in RPo occur entirely in one mechanistic step at the lambdaPR promoter and (iii) does the structure of RPo vary with temperature? In addition, we use the MHF-determined population distribution of KMnO4-reactive (RPo) and non-reactive promoter DNA to normalize the biphasic kinetics of decay of RPo to free promoter DNA after a 37 degrees to 0 degreesC temperature downshift, and thereby characterize the kinetics of the conformational changes involved in forming RPo.

Thermodynamics of the interactions of lac repressor with variants of the symmetric lac operator: effects of converting a consensus site to a non-specific site

What are the thermodynamic consequences of the stepwise conversion of a highly specific (consensus) protein-DNA interface to one that is nonspecific? How do the magnitudes of key favorable contributions to complex stability (burial of hydrophobic surfaces and reduction of DNA phosphate charge density) change as the DNA sequence of the specific site is detuned? To address these questions we investigated the binding of lac repressor (LacI) to a series of 40 bp fragments carrying symmetric (consensus) and variant operator sequences over a range of temperatures and salt concentrations. Variant DNA sites contained symmetrical single and double base-pair substitutions at positions 4 and/or 5 [sequence: see text] in each 10 bp half site of the symmetric lac operator (Osym). Non-specific interactions were examined using a 40 bp non-operator DNA fragment. Disruption of the consensus interface by a single symmetrical substitution reduces the observed equilibrium association constant (K(obs)) for Osym by three to four orders of magnitude; double symmetrical substitutions approach the six orders in magnitude difference between specific and non-specific binding to a 40 bp fragment. At these adjacent positions in the consensus site, the free energy effects of multiple substitutions are non-additive: the first reduces /deltaG(obs)o/ by 3 to 5 kcal mol(-1), approximately halfway to the non-specific level, whereas the second is less deleterious, reducing /deltaG(obs)o/ by less than 3 kcal mol(-1). Variant-specific dependences of K(obs) on temperature and salt concentration characterize these LacI-operator interactions. In general, binding constants and standard free energies of binding both exhibit characteristic extrema near 290 K. As a consequence, both the enthalpic and entropic contributions to stability of Osym and variant complexes change from positive (i.e. entropy driven) at lower temperatures to negative (i.e. enthalpy driven) at higher temperatures, indicating that the heat capacity change upon binding, deltaC(obs)o, is large and negative. In general, /deltaC(obs)o/ decreases as the specificity and stability of the variant complex decreases. Stabilities of complexes of LacI with Osym and all variant operators are strongly [salt]-dependent. Binding constants for the variant complexes exhibit a power-dependence on [salt] that is larger in magnitude (i.e. more negative) than for Osym, but no obvious trend relates changes in contributions from the polyelectrolyte effect and the observed reductions in stability (delta deltaG(obs)o). These variant-specific thermodynamic signatures provide novel insights into the consequences of converting a consensus interface to a less specific one; such insights are not obtained from comparisons at the level of delta deltaG(obs)o. We propose that this variant-specific behavior arises from a strong effect of operator sequence on the extent of induced conformational changes in the protein (and possibly also in the DNA site) which accompany binding.

Cooperative and anticooperative effects in binding of the first and second plasmid Osym operators to a LacI tetramer: evidence for contributions of non-operator DNA binding by wrapping and looping

The interaction of lac operator DNA with lac repressor (LacI) is a classic example of a genetic regulatory switch. To dissect the role of stoichiometry, subunit association, and effects of DNA length in positioning this switch, we have determined binding isotherms for the interaction of LacI with a high affinity (Osym) operator on linearized plasmid (2500 bp) DNA over a wide range of macromolecular concentrations (10(-14) to 10(-8) M). Binding data were analyzed using a thermodynamic model involving four equilibria: dissociation of tetramers (T) into dimers (D), and binding of operator-containing plasmid DNA (O) to dimers and tetramers to form three distinct complexes, DO, TO, and TO2. Over the range of concentrations of repressor, operator, and salt (0.075 M K+ to 0.40 M K+) investigated, we find no evidence for any significant thermodynamic effect of LacI dimers. Instead, all isotherms can be interpreted in terms of just two equilibria, involving only T and the TO and TO2 complexes. As a reference binding equilibrium, which we propose must approximate the DO binding interaction, we compare the plasmid Osym results with our extensive studies of the binding of a 40 bp Osym DNA fragment to LacI. On this basis, we obtain a lower bound on the LacI dimer-tetramer equilibrium constant and values of the equilibrium constants for formation of TO and TO2 complexes. At a salt concentration of 0.40 M, the Osym plasmid binding data are consistent with a model with two independent and identical binding sites for operator per LacI tetramer, in which the binding to a site on the tetramer is only slightly more favorable than the reference binding interaction. Increasingly large deviations from the independent-site model are observed as the salt concentration is reduced; binding of a second operator to from TO2 becomes strongly disfavored relative to formation of TO at low salt concentrations (0.075 to 0.125 M). In addition, binding of both the first and second plasmid operator DNA molecules to the tetramer becomes increasingly more favorable than the reference binding interaction as [K+] is reduced from 0.40 M to 0.125 M. At 0.075 M K+, however, the strength of binding of the second plasmid operator DNA to the LacI tetramer is dramatically reduced; this interaction is much less favorable than binding the first plasmid operator DNA, and becomes much less favorable than the reference binding interaction. We propose that these differences arise from changes in the nature of the TO and TO2 complexes with decreasing salt concentration. At low salt concentration, we suggest the hypothesis that flanking non-operator sequences bind non-specifically (coulombically) by local wrapping, and that distant regions of non-operator DNA occupy the second operator-binding site by looping. We propose that wrapping stabilizes both 1:1 and 2:1 complexes at low salt concentration, and that looping stabilizes the 1:1 complex but competitively destabilizes the 2:1 TO2 complex at low salt concentration. These effects must play a role in adjusting the stability and structure of the LacI-lac operator repression complex as the cytoplasmic [K+] varies in response to changes in extracellular osmolarity.