Comparison of solution structure of free and complexed lac operator by molecular modelling with NMR constraints

The experimental folding landscape of monomeric lactose repressor, a large two-domain protein, involves two kinetic intermediates

To probe the experimental folding behavior of a large protein with complex topology, we created a monomeric variant of the lactose repressor protein (MLAc), a well characterized tetrameric protein that regulates transcription of the lac operon. Purified MLAc is folded, fully functional, and binds the inducer isopropyl beta-d-thiogalactoside with the same affinity as wild-type LacI. Equilibrium unfolding of MLAc induced by the chemical denaturant urea is a reversible, apparent two-state process (pH 7.5, 20 degrees C). However, time-resolved experiments demonstrate that unfolding is single-exponential, whereas refolding data indicate two transient intermediates. The data reveal the initial formation of a burst-phase (tau < ms) intermediate that corresponds to approximately 50% of the total secondary-structure content. This step is followed by a rearrangement reaction that is rate-limited by an unfolding process (tau approximately 3 s; pH 7.5, 20 degrees C) and results in a second intermediate. This MLAc intermediate converts to the native structure (tau approximately 30 s; pH 7.5, 20 degrees C). Remarkably, the experimental folding-energy landscape for MLAc is in excellent agreement with theoretical predictions using a simple topology-based C(alpha)-model as presented in a companion article in this issue.

Distortions of the DNA double helix induced by 1,3-trans-diamminedichloroplatinum(II)-intrastrand cross-link: an internal coordinate molecular modeling study

A trans-diamminedichloroplatinum(II) (trans-DDP) intrastrand adduct within the sequence d(TCTG*TG*TC).d(GACACAGA) (where G* represents a platinated guanine) is modeled on the basis of qualitative experimental data concerning global unwinding and curvature as well as information on base pairing. Modeling is performed using the internal coordinate JUMNA program, specific to nucleic acids, and modified to include the possibility of covalently bound ligands. Calibration of the energy functions representing the Pt-N7 bond with guanine is described. The platinum atom and the platinum-nitrogen bonds are parameterized for use in the Hückel Del Re method to calculate monopoles at each atom. These monopoles are consistent with the Flex force field included in Jumna. By developing an appropriate minimization protocol we are able to generate stable, distorted three-dimensional structures compatible with the experimental data and including an unusually high global unwinding. No a priori geometric assumptions are made in generating these structures.

Reorganization of terminator DNA upon binding replication terminator protein: implications for the functional replication fork arrest complex

Termination of DNA replication in Bacillus subtilis involves the polar arrest of replication forks by a specific complex formed between the replication terminator protein (RTP) and DNA terminator sites. While determination of the crystal structure of RTP has facilitated our understanding of how a single RTP dimer interacts with terminator DNA, additional information is required in order to understand the assembly of a functional fork arrest complex, which requires an interaction between two RTP dimers and the terminator site. In this study, we show that the conformation of the major B.subtilis DNA terminator,TerI, becomes considerably distorted upon binding RTP. Binding of the first dimer of RTP to the B site of TerI causes the DNA to become slightly unwound and bent by approximately 40 degrees. Binding of a second dimer of RTP to the A site causes the bend angle to increase to approximately 60 degrees . We have used this new data to construct two plausible models that might explain how the ternary terminator complex can block DNA replication in a polar manner. In the first model, polarity of action is a consequence of the two RTP-DNA half-sites having different conformations. These different conformations result from different RTP-DNA contacts at each half-site (due to the intrinsic asymmetry of the terminator DNA), as well as interactions (direct or indirect) between the RTP dimers on the DNA. In the second model, polar fork arrest activity is a consequence of the different affinities of RTP for the A and B sites of the terminator DNA, modulated significantly by direct or indirect interactions between the RTP dimers.

Selectively 13C-enriched DNA: evidence from 13C1' relaxation rate measurements of an internal dynamics sequence effect in the lac operator

In order to study some internal dynamic processes of the lac operator sequence, the 13C-labeled duplex 5'd(C0G1C2T3C4A5C6A7A8T9T10).d(A10A9T8T7G6T5G4A3G2C1G0)3' was used. The spreading of both the 1H1' and 13C1' resonances brought about an excellent dispersion of the 1H1'-13C1' correlations. The spin-lattice relaxation parameters R(Cz), R(Cx,y) and R(Hz --> Cz) were measured for each residue of the two complementary strands, except for the 3'-terminal residues which were not labeled. Variation of the relaxation rates was found along the sequence. These data were analyzed in the context of the model-free formalism proposed by Lipari and Szabo [(1982) J. Am. Chem. Soc., 104, 4546-4570] and extended to three parameters by Clore et al. [(1990) Biochemistry, 29, 7387-7401; and (1990) J. Am. Chem. Soc., 112, 4989-4991]. A careful analysis using a least-squares program showed that our data must be interpreted in terms of a three-parameter spectral density function. With this approach, the global correlation time was found to be the same for each residue. All the C1'-H1' fragments exhibited both slow (tau s = 1.5 ns) and fast (tau f = 20 ps) restricted libration motions (Ss2 = 0.74 to 1.0 and Sf2 = 0.52 to 0.96). Relaxation processes were described as governed by the motion of the sugar relative to the base and in terms of bending of the whole duplex. The possible role played by the special structure of the AATT sequence is discussed. No evident correlation was found between the amplitude motions of the complementary residues. The 5'-terminal residues showed large internal motions (S2 = 0.5), which describe the fraying of the double helix. Global examination of the microdynamical parameters Sf2 and Ss2 along the nucleotide sequence showed that the adenine residues exhibit more restricted fast internal motions (Sf2 = 0.88 to 0.96) than the others, whereas the measured relaxation rates of the four nucleosides in solution were mainly of dipolar origin. Moreover, the fit of both R(Cz) and R(Hz --> Cz) experimental relaxation rates using an only global correlation time for all the residues, gave evidence of a supplementary relaxation pathway affecting R(Cx,y) for the purine residues in the (5' --> 3') G4A3 and A10A9T8T7 sequences. This relaxation process was analyzed in terms of exchange stemming from motions of the sugar around the glycosidic bond on the millisecond time scale. It should be pointed out that these residues gave evidence of close contacts with the protein in the complex with the lac operator [Boelens et al. (1987) J. Mol. Biol., 193, 213-216] and that these motions could be implied in the lac-operator-lac-repressor recognition process.

Solution structure of oligonucleotides covalently linked to a psoralen derivative

Psoralen (pso) was attached via its C-5 position to the 5'-phosphate group of an oligodeoxynucleotide d(TAAGCCG) by a hexamethylene linker (m6). Complex formation between pso-m6-d(TAAGCCG) and the complementary strands d(CGGCTTA)[7-7mer] or d(CGGCTTAT)[7-8mer] was investigated by nuclear magnetic resonance in aqueous solution. Structural informations derived from DQF-COSY and NOESY maps, revealed that the mini double helix adopts a B-form conformation and that the deoxyriboses preferentially adopt a C2'-endo conformation. The nOe connectivities observed between the protons of the bases or the sugars in each duplex, and the protons of the psoralen and the hexamethylene chain, led us to propose a model involving an equilibrium between two conformations due to different locations of the psoralen. Upon UV-irradiation, the psoralen moiety cross-linked the two DNA strands at the level of 5'TpA3' sequences. NMR studies of the single major photo-cross-linked duplex pso-m6-d(TAAGCCG) and d(CGGCTTA) were performed. The stereochemistry of the diadduct is indeed cis-syn at both cyclobutane rings. In addition, the effects of this diadduct on the helical structure are analyzed in detail.