PvuII-endonuclease induces structural alterations at the scissile phosphate group of its cognate DNA

Role of Mg2+ Ions in DNA Hydrolysis by EcoRV, Studied by the 3D-Reference Interaction Site Model and Molecular Dynamics

The role of Mg²⁺ ions during precursor formation in DNA hydrolysis by the homodimeric restriction enzyme EcoRV was elucidated based on the 3D-reference interaction site model (RISM) theory and the molecular dynamics (MD) simulation. From an analysis of the spatial distribution of Mg²⁺ in an active site using 3D-RISM, we identified a new position for Mg²⁺ in the X-ray EcoRV-DNA complex structure ( 1rvb ). We refer to the position as site IV^†. Site IV^† is almost at the same position as that of a Ca²⁺ ion in the superimposed X-ray crystallographic active-site structure of the PvuII-DNA complex ( 1f0o ). 3D-RISM was also used to locate the position of water molecules, including the water nucleophile at the active site. MD simulations were carried out with the initial structure having two Mg²⁺ ions at site IV^† and at site I*, experimentally identified by Horton et al., to find a stable complex structure in which the DNA fragment was rearranged to orient the scissile bond direction toward the water nucleophile. The equilibrium active-site structure of the EcoRV-DNA complex obtained from the MD simulation was similar to the superimposed X-ray crystallographic structure of the BamHI-DNA complex ( 2bam ). In the active-site structure, two metal ions have almost the same position (≤1.0 Å) as that of 2bam , and the scissile phosphate is twisted to orient the scissile bond toward the water nucleophile, as is the case in 2bam . We propose the equilibrium active-site structure obtained in this study as a precursor for the hydrolysis reaction of EcoRV.

Structure and dynamics of phosphate ion in aqueous solution: an ab initio QMCF MD study

A simulation of phosphate in aqueous solution was carried out employing the new QMCF MD approach which offers the possibility to investigate composite systems with the accuracy of a QMMM method but without the time consuming creation of solute-solvent potential functions. The data of the simulations give a clear picture of the hydration shells of the phosphate anion. The first shell consists of 13 water molecules and each oxygen of the phosphate forms in average three hydrogens bonds to different solvent molecules. Several structural parameters such as radial distribution functions and coordination number distributions allow to fully characterize the embedding of the highly charged phosphate ion in the solvent water. The dynamics of the hydration structure of phosphate are described by mean residence times of the solvent molecules in the first hydration shell and the water exchange rate.

Molecular dynamics studies on free and bound targets of the bovine papillomavirus type I e2 protein: the protein binding effect on DNA and the recognition mechanism

Molecular dynamics simulations of a total duration of 30 ns in explicit solvent were carried out on the BPV-1-E2 protein complexed to a high-affinity DNA target containing the two hydrogen-bonded ACCG.CGGT half-sites separated by the noncontacted ACGT sequence. The analysis of the trajectories focuses on the DNA structure and on the dynamics. The data are compared to those issued from recent simulations made on three free targets that recognize E2 with different affinities. E2 does not drastically perturb the mechanic properties of the free DNA: the structural relationships between the BI/BII backbone substates and some helical parameters are preserved in the complex despite a severe slowing down of the phosphate group motions. The structures of both free and bound half-sites are very close to each other although the conformational space explored by these regions is narrowed when they are contacted by the protein. The enhanced plasticity found in the best free target spacers, mainly manifested through the backbone motions, allows a clear overlap between several free and bound global DNA features such as the base displacement. Furthermore, this flexibility is preserved in the complex. Our results support the hypothesis that E2 takes advantage of free predistorted structures that may minimize the DNA deformation cost. In addition, we observe that E2 is far from totally stiffening the DNA, suggesting that the entropic penalty inherent in the complex formation could be limited.

Towards an understanding of DNA recognition by the methyl-CpG binding domain 1

CpG methylation determines a variety of biological functions of DNA. The methylation signal is interpreted by proteins containing a methyl-CpG binding domain (MBDs). Based on the NMR structure of MBD1 complexed with methylated DNA we analysed the recognition mode by means of molecular dynamics simulations. As the protein is monomeric and recognizes a symmetrically methylated CpG step, the recognition mode is an asymmetric one. We find that the two methyl groups do not contribute equally to the binding energy. One methyl group is associated with the major part of the binding energy and the other one nearly does not contribute at all. The contribution of the two cytosine methyl groups to binding energy is calculated to be -3.6 kcal/mol. This implies a contribution of greater than two orders of magnitude to the binding constant. The conserved amino acid Asp32 is known to be essential for DNA binding by MBD1, but so far no direct contact with DNA has been observed. We detected a direct DNA base contact to Asp32. This could be the main reason for the importance of this amino acid. MBD contacts DNA exclusively in the major groove, the minor groove is reserved for histone contacts. We found a deformation of the minor groove shape due to complexation by MBD1, which indicates an information transfer between the major and the minor groove.

Pre-organized structure of viral DNA at the binding-processing site of HIV-1 integrase

The integration of the human immunodeficiency virus type 1 DNA into the host cell genome is catalysed by the viral integrase (IN). The reaction consists of a 3'-processing [dinucleotide released from each 3' end of the viral long terminal repeat (LTR)] followed by a strand transfer (insertion of the viral genome into the human chromosome). A 17 base pair oligonucleotide d(GGAAAATCTCTAGCAGT), d(ACTGCTAGAGATTTTCC) reproducing the U5-LTR extremity of viral DNA that contains the IN attachment site was analysed by NMR using the classical NOEs and scalar coupling constants in conjunction with a small set of residual dipolar coupling constants (RDCs) measured at the 13C/15N natural abundance. The combination of these two types of parameters in calculations significantly improved the DNA structure determination. The well-known features of A-tracts were clearly identified by RDCs in the first part of the molecule. The binding/cleavage site at the viral DNA end is distinguishable by a loss of regular base stacking and a distorted minor groove that can aid its specific recognition by IN.

Single-molecule kinetics of lambda exonuclease reveal base dependence and dynamic disorder

We used a multiplexed approach based on flow-stretched DNA to monitor the enzymatic digestion of lambda-phage DNA by individual bacteriophage lambda exonuclease molecules. Statistical analyses of multiple single-molecule trajectories observed simultaneously reveal that the catalytic rate is dependent on the local base content of the substrate DNA. By relating single-molecule kinetics to the free energies of hydrogen bonding and base stacking, we establish that the melting of a base from the DNA is the rate-limiting step in the catalytic cycle. The catalytic rate also exhibits large fluctuations independent of the sequence, which we attribute to conformational changes of the enzyme-DNA complex.