Radiation damage to a DNA-binding protein. Combined circular dichroism and molecular dynamics simulation analysis.

The E. coli lactose operon, the paradigm of gene expression regulation systems, is the best model for studying the effect of radiation on such systems. The operon function requires the binding of a protein, the repressor, to a specific DNA sequence, the operator. We have previously shown that upon irradiation the repressor loses its operator binding ability. The main radiation-induced lesions of the headpiece have been identified by mass spectrometry. All tyrosine residues are oxidized into 3,4-dihydroxyphenylalanine (DOPA). In the present study we report a detailed characterization of the headpiece radiation-induced modification. An original approach combining circular dichroism measurements and the analysis of molecular dynamics simulation of headpieces bearing DOPA-s instead of tyrosines has been applied. The CD measurements reveal an irreversible modification of the headpiece structure and stability. The molecular dynamics simulation shows a loss of stability shown by an increase in internal dynamics and allows the estimation of the modifications due to tyrosine oxidation for each structural element of the protein. The changes in headpiece structure and stability can explain at least in part the radiation-induced loss of binding ability of the repressor to the operator. This conclusion should hold for all proteins containing radiosensitive amino acids in their DNA-binding site.

Molecular dynamics simulation for probing the flexibility of the 35 nucleotide SL1 sequence kissing complex from HIV-1Lai genomic RNA.

The SL1 stem-loop located in the encapsidation domain is responsible for initiating the dimerisation of HIV-1 genomic RNA by means of a loop-loop interaction known as Kissing Complex (KC). The SL1 secondary structure has been predicted as a 35 nucleotides [K. G. Murti, M. Bondurant, and A. Tereba. J Virol 37, 411-419 (1981)] stem-loop composed of a 4 base pairs (bp) terminal duplex, a 4 nt asymmetrical internal loop, a 7 bp internal duplex, and a 9 nt apical loop. Several high resolution structures of the monomer and of KC of a 23 nt sequence containing only the internal duplex and the apical loop of SL1 are available in the literature. No experimental high resolution structure of the complete native SL1 sequence has been reported so far, either for the monomer or for KC. The asymmetrical internal loop has been described from NMR studies of different monomeric hairpin sequences, leading to divergent results, which suggests its high flexibility. In this work, we built a SL1(35) KC model which was submitted to a 31 ns molecular dynamics simulation (MD). Our results allows to describe the internal dynamics of SL1(35) KC and the differences of behavior of the different parts of the dimer. Thus, we could show the stability of the interactions between the two apical loops and of the terminal duplexes, the destabilization of the internal duplexes and the high flexibility of the asymmetrical internal loops.

Conformational pathway for the kissing complex-->extended dimer transition of the SL1 stem-loop from genomic HIV-1 RNA as monitored by targeted molecular dynamics techniques.

HIV-1 retroviral genomic RNA dimerization is initiated by loop-loop interactions between the SL1 stem-loops of two identical RNA molecules. The SL1-SL1 unstable resulting kissing complex (KC) then refolds irreversibly into a more stable complex called extended dimer (ED). Although the structures of both types of complex have been determined, very little is known about the conformational pathway corresponding to the transition, owing to the difficulty of observing experimentally intermediate conformations. In this study, we applied targeted molecular dynamics simulation techniques (TMD) to the phosphorus atoms for monitoring this pathway for the backbone, and a two-step strategy was adopted. In a first step, called TMD(-1), the dimer structure was constrained to progressively move away from KC without indicating the direction, until the RMSD from KC reaches 36A. A total of 20 TMD(-1) simulations were performed under different initial conditions and different simulation parameters. For RMSD ranging between 0 and 22A, the whole set of TMD(-1) simulations follows a similar pathway, then divergences are observed. None of the simulations leads to the ED structure. At RMSD=22A, the dimers look like two parallel Us, still linked by the initial loop-loop interaction, but the strands of the stems (the arms of the Us) are positioned in such a manner that they can form intramolecular as well as intermolecular Watson-Crick base-pairs. This family of structure is called UU. In a second step (TMD simulations), 18 structures were picked up along the pathways generated with TMD(-1) and were constrained to move toward ED by decreasing progressively their RMSD from ED. We found that only structures from the UU family are able to easily reach ED-like conformations of the backbones without exhibiting a large constraint energy.