Molecular dynamics simulation of trp-repressor/operator complex: analysis of hydrogen bond patterns of protein–DNA interaction

Long-time instability in the Runge-Kutta algorithm for a Nosé-Hoover heat bath of a harmonic chain and its stabilization

In this paper, we investigate the Runge-Kutta algorithm for the Nosé-Hoover heat bath of a harmonic chain. The Runge-Kutta algorithm is found to be unstable in long-time calculations, with the system temperature growing exponentially. The growth rate increases if time step size is chosen larger. By analyzing the Fourier spectra in both space (wave number) and time (frequency), we discover that the growth is caused by spurious energy accumulation, particularly at the largest wave number. Such accumulation may be explained by von Neumann analysis for an infinite chain, with the nonlinear heat bath being ignored. Furthermore, we propose to add a filter to remove excessive energy, which effectively stabilizes the algorithm.

Microscopic understanding of the conformational features of a protein–DNA complex

Protein–DNA interactions play crucial roles in different stages of genetic activities, such as replication of genome, initiation of transcription,etc.

Recovery of the poisoned topoisomerase II for DNA religation: coordinated motion of the cleavage core revealed with the microsecond atomistic simulation

Type II topoisomerases resolve topological problems of DNA double helices by passing one duplex through the reversible double-stranded break they generated on another duplex. Despite the wealth of information in the cleaving operation, molecular understanding of the enzymatic DNA ligation remains elusive. Topoisomerase poisons are widely used in anti-cancer and anti-bacterial therapy and have been employed to entrap the intermediates of topoisomerase IIβ with religatable DNA substrate. We removed drug molecules from the structure and conducted molecular dynamics simulations to investigate the enzyme-mediated DNA religation. The drug-unbound intermediate displayed transitions toward the resealing-compliant configuration: closing distance between the cleaved DNA termini, B-to-A transformation of the double helix, and restoration of the metal-binding motif. By mapping the contact configurations and the correlated motions between enzyme and DNA, we identified the indispensable role of the linker preceding winged helix domain (WHD) in coordinating the movements of TOPRIM, the nucleotide-binding motifs, and the bound DNA substrate during gate closure. We observed a nearly vectorial transition in the recovery of the enzyme and identified the previously uncharacterized roles of Asn508 and Arg677 in DNA rejoining. Our findings delineate the dynamic mechanism of the DNA religation conducted by type II topoisomerases.

Conformational stability and DNA binding energetics of the rat thyroid transcription factor 1 homeodomain

The conformational stability of the rat thyroid transcription factor 1 homeodomain, TTF-1HD, has been investigated by means of circular dichroism (CD) and differential scanning calorimetry (DSC) measurements at pH 5.0 as a function of KCl concentration. Thermal unfolding of TTF-1HD is a reversible two-state transition. The protein is not stable against temperature, showing a denaturation temperature of 32 degrees C in the absence of salt and 50 degrees C at 75 mM KCl. The binding energetics of TTF-1HD to its target DNA sequence has been characterized by means of isothermal titration calorimetry (ITC) measurements, complemented with CD data. At 25 degrees C, pH 5.0 and 75 mM KCl, the binding constant amounts to 1.5 x 10(8)M(-1) and the binding enthalpy change amounts to -41 kJ mol(-1). The process is enthalpy driven, but also the entropy change is favorable to complex formation. To gain a molecular level understanding of the interactions determining the association of TTF-1HD to the target DNA sequence structural information would be requested, but it is not yet available. Therefore, structural models of two complexes, TTF-1HD with the target DNA sequence and TTF-1HD with a modified DNA sequence, have been constructed by using as a template the NMR structure of the complex between NK-2 HD and its target DNA, and by performing molecular dynamics simulations 3.5 ns long. Analysis of these models allows one to shed light on the origin of the DNA binding specificity characteristic of TTF-1HD.

A study of collective atomic fluctuations and cooperativity in the U1A-RNA complex based on molecular dynamics simulations

Cooperative interactions play an important role in recognition and binding in macromolecular systems. In this study, we find that cross-correlated atomic fluctuations can be used to identify cooperative networks in a protein-RNA system. The dynamics of the RRM-containing protein U1A-stem loop 2 RNA complex have been calculated theoretically from a 10 ns molecular dynamics (MD) simulation. The simulation was analyzed by calculating the covariance matrix of all atomic fluctuations. These matrix elements are then presented in the form of a two-dimensional grid, which displays fluctuations on a per residue basis. The results indicate the presence of strong, selective cross-correlated fluctuations throughout the RRM in U1A-RNA. The atomic fluctuations correspond well with previous biophysical studies in which a multiplicity of cooperative networks have been reported and indicate that the various networks identified in separate individual experiments are fluctuationally correlated into a hyper-network encompassing most of the RRM. The calculated results also correspond well with independent results from a statistical covariance analysis of 330 aligned RRM sequences. This method has significant implications as a predictive tool regarding cooperativity in the protein-nucleic acid recognition process.

The role of water in protein-DNA recognition

Is it by design or by default that water molecules are observed at the interfaces of some protein-DNA complexes? Both experimental and theoretical studies on the thermodynamics of protein-DNA binding overwhelmingly support the extended hydrophobic view that water release from interfaces favors binding. Structural and energy analyses indicate that the waters that remain at the interfaces of protein-DNA complexes ensure liquid-state packing densities, screen the electrostatic repulsions between like charges (which seems to be by design), and in a few cases act as linkers between complementary charges on the biomolecules (which may well be by default). This review presents a survey of the current literature on water in protein-DNA complexes and a critique of various interpretations of the data in the context of the role of water in protein-DNA binding and principles of protein-DNA recognition in general.

Water-mediated contacts in thetrp-repressor operator complex recognition process

Water-mediated contacts are known as an important recognition tool in trp-repressor operator systems. One of these contacts involves two conserved base pairs (G(6).C(-6) and A(5). T(-5)) and three amino acids (Lys 72, Ile 79, and Ala 80). To investigate the nature of these contacts, we analyzed the X-ray structure (PDB code: 1TRO) of the trp-repressor operator complex by means of molecular dynamics simulations. This X-ray structure contains two dimers that exhibit structural differences. From these two different starting structures, two 10 ns molecular dynamics simulations have been performed. Both of our simulations show an increase of water molecules in the major groove at one side of the dimer, while the other side remains unchanged compared to the X-ray structure. Though the maximum residence time of the concerned water molecules decreases with an increase of solvent at the interface, these water molecules continue to play an important role in mediating DNA-protein contacts. This is shown by new stable amino acids-DNA distances and a long water residence time compared to free DNA simulation. To maintain stability of the new contacts, the preferential water binding site on O6(G6) is extended. This extension agrees with mutation experiment data on A5 and G6, which shows different relative affinity due to mutation on these bases [A. Joachimiak, T. E. Haran, P. B. Sigler, EMBO Journal 1994, Vol. 13, No. (2) pp. 367-372]. Due to the rearrangements in the system, the phosphate of the base G6 is able to interconvert to the B(II) substate, which is not observed on the other half side of the complex. The decrease of the number of hydrogen bonds between protein and DNA backbone could be the initial step of the dissociation process of the complex, or in other words an intermediate complex conformation of the association process. Thus, we surmise that these features show the importance of water-mediated contacts in the trp-repressor operator recognition process.

Simulations of nucleic acids and their complexes

Recent years have seen considerable progress in simulations of nucleic acids. Improvements in force fields, simulation techniques and protocols, and increasing computer power have all contributed to making nanosecond-scale simulations of both DNA and RNA commonplace. The results are already helping to explain how nucleic acids respond to their environment and to their base sequence and to reveal the factors underlying recognition processes by probing biologically important nucleic acid-protein interactions and medically important nucleic acid-drug complexation. This Account summarizes methodological progress and applications of molecular dynamics to nucleic acids over the past few years and tries to identify remaining challenges.

Experimental and theoretical high pressure strategies for investigating protein-nucleic acid assemblies

A method was developed to investigate the stability of protein-nucleic acid complexes using hydrostatic pressure during electrophoretic gel mobility shift analysis. The initial system probed by this technique was the well-characterized cognate BamHI-DNA complex. Band shift analysis at several elevated pressures found the equilibrium dissociation (K(d)) constant to be dependent on pressure, which allowed the volume change of dissociation (deltaV) to be calculated. In order to describe the effects of pressure on the specific BamHI-DNA complex at the molecular level, molecular dynamics simulations at both ambient and elevated pressure was performed. Comparison of the simulation trajectories identified several individual BamHI-DNA contacts that are disrupted due to pressure. The disruption of these contacts can be attributed to an observed pressure-induced increase in hydration at the protein-DNA interface during the elevated pressure simulation.