Molecular modeling of closed circular DNA thermodynamic ensembles

Many modeling studies of supercoiled DNA are based on equilibrium structures from theoretical calculations or energy minimization. Since closed circular DNAs are flexible, it is possible that errors are introduced by calculating properties from a single minimum energy structure, rather than from a complete thermodynamic ensemble. We have investigated this question using molecular dynamics simulations on a low resolution molecular mechanics model in which each base pair is represented by three points (a plane). This allows the inclusion of sequence-dependent variations of tip, inclination, and twist. Three kinds of sequences were tested: (1) homogeneous DNA, in which all base pairs have the helicoidal parameters of an ideal, average B-DNA; (2) random sequence DNA; and (3) curved DNA. We examined the rate of convergence of various structural parameters. Convergence for most of these is slowest for homogeneous sequences, more rapid for random sequences, and most rapid for curved sequences. The most slowly converging parameter is the antipodes profile. In a plasmid with N base pairs (bp), the antipodes distance is the distance d _ij from base pair i to base pair j halfway around the plasmid, j = i + N/2. The antipodes profile at time t is a plot of d _ij over the range i = 1, N/2. In a homogeneous plasmid, convergence requires that the antipodes profile averaged over time must be flat. Even in the small plasmids examined here, the average properties of the ensembles were found to differ from those of static equilibrium structures. These effects will be even more dramatic for larger plasmids. Further, average and dynamic properties are affected by both plasmid size and sequence. © 1996 John Wiley & Sons, Inc.     

A multiscale dynamic model of DNA supercoil relaxation by topoisomerase IB

In this study, we report what we believe to be the first multiscale simulation of the dynamic relaxation of DNA supercoils by human topoisomerase IB (topo IB). We leverage our previous molecular dynamics calculations of the free energy landscape describing the interaction between a short DNA fragment and topo IB. Herein, this landscape is used to prescribe boundary conditions for a computational, elastodynamic continuum rod model of a long length of supercoiled DNA. The rod model, which accounts for the nonlinear bending, twisting, and electrostatic interaction of the (negatively charged) DNA backbone, is extended to include the hydrodynamic drag induced by the surrounding physiological buffer. Simulations for a 200-bp-long DNA supercoil in complex with topo IB reveal a relaxation timescale of ∼0.1–1.0 μs. The relaxation follows a sequence of cascading reductions in the supercoil linking number (Lk), twist (Tw), and writhe (Wr) that follow companion cascading reductions in the supercoil elastic and electrostatic energies. The novel (to our knowledge) multiscale modeling method may enable simulations of the entire experimental setup that measures DNA supercoiling and relaxation via single molecule magnetic trapping.     

The dynamics of camphor in the cytochrome P450 CYP101D2

The recent crystal structures of CYP101D2, a cytochrome P450 protein from the oligotrophic bacterium Novosphingobium aromaticivorans DSM12444 revealed that both the native (substrate‐free) and camphor‐soaked forms have open conformations. Furthermore, two other potential camphor‐binding sites were also identified from electron densities in the camphor‐soaked structure, one being located in the access channel and the other in a cavity on the surface near the F‐helix side of the F‐G loop termed the substrate recognition site. These latter sites may be key intermediate positions on the pathway for substrate access to or product egress from the active site. Here, we show via the use of unbiased atomistic molecular dynamics simulations that despite the open conformation of the native and camphor‐bound crystal structures, the underlying dynamics of CYP101D2 appear to be very similar to other CYP proteins. Simulations of the native structure demonstrated that the protein is capable of sampling many different conformational substates. At the same time, simulations with the camphor positioned at various locations within the access channel or recognition site show that movement towards the active site or towards bulk solvent can readily occur on a short timescale, thus confirming many previously reported in silico studies using steered molecular dynamics. The simulations also demonstrate how the fluctuations of an aromatic gate appear to control access to the active site. Finally, comparison of camphor‐bound simulations with the native simulations suggests that the fluctuations can be of similar level and thus are more representative of the conformational selection model rather than induced fit.     

Curcumin stably interacts with DNA hairpin through minor groove binding and demonstrates enhanced cytotoxicity in combination with FdU nucleotides

We report, based on biophysical studies and molecular mechanical calculations that curcumin binds DNA hairpin in the minor groove adjacent to the loop region forming a stable complex. UV–Vis and fluorescence spectroscopy indicated interaction of curcumin with DNA hairpin. In this novel binding motif, two ? H of curcumin heptadiene chain are closely positioned to the A₁₆-H8 and A₁₇-H8, while G₁₂-H8 is located in the close proximity of curcumin α H. Molecular dynamics (MD) simulations suggest, the complex is stabilized by noncovalent forces including; π-π stacking, H-bonding and hydrophobic interactions. Nuclear magnetic resonance (NMR) spectroscopy in combination with molecular dynamics simulations indicated curcumin is bound in the minor groove, while circular dichroism (CD) spectra suggested minute enhancement in base stacking and a little change in DNA helicity, without significant conformational change of DNA hairpin structure. The DNA:curcumin complex formed with FdU nucleotides rather than Thymidine, demonstrated enhanced cytotoxicity towards oral cancer cells relative to the only FdU substituted hairpin. Fluorescence co-localization demonstrated stability of the complex in biologically relevant conditions, including its cellular uptake. Acridine orange/EtBr staining further confirmed the enhanced cytotoxic effects of the complex, suggesting apoptosis as mode of cell death. Thus, curcumin can be noncovalently complexed to small DNA hairpin for cellular delivery and the complex showed increased cytotoxicity in combination with FdU nucleotides, demonstrating its potential for advanced cancer therapy.     

Structural Dynamics of Human Telomeric G-Quadruplex Loops Studied by Molecular Dynamics Simulations

Loops which are linkers connecting G-strands and supporting the G-tetrad core in G-quadruplex are important for biological roles of G-quadruplexes. TTA loop is a common sequence which mainly resides in human telomeric DNA (hTel) G-quadruplex. A series of molecular dynamics (MD) simulations were carried out to investigate the structural dynamics of TTA loops. We found that (1) the TA base pair formed in TTA loops are very stable, the occupied of all hydrogen bonds are more than 0.95. (2) The TA base pair makes the adjacent G-quartet more stable than others. (3) For the edgewise loop and the diagonal loop, most loop bases are stacking with others, only few bases have considerable freedom. (4) The stabilities of these stacking structures are distinct. Part of the loops, especially TA base pairs, and bases stacking with the G-quartet, maintain certain stable conformations in the simulation, but other parts, like TT and TA stacking structures, are not stable enough. For the first time, spontaneous conformational switches of TTA edgewise loops were observed in our long time MD simulations. (5) For double chain reversal loop, it is really hard to maintain a stable conformation in the long time simulation under present force fields (parm99 and parmbsc0), as it has multiple conformations with similar free energies.     

Identifying molecular dynamics in single-molecule FRET experiments with burst variance analysis

Histograms of single-molecule Förster resonance energy transfer (FRET) efficiency are often used to study the structures of biomolecules and relate these structures to function. Methods like probability distribution analysis analyze FRET histograms to detect heterogeneities in molecular structure, but they cannot determine whether this heterogeneity arises from dynamic processes or from the coexistence of several static structures. To this end, we introduce burst variance analysis (BVA), a method that detects dynamics by comparing the standard deviation of FRET from individual molecules over time to that expected from theory. Both simulations and experiments on DNA hairpins show that BVA can distinguish between static and dynamic sources of heterogeneity in single-molecule FRET histograms and can test models of dynamics against the observed standard deviation information. Using BVA, we analyzed the fingers-closing transition in the Klenow fragment of Escherichia coli DNA polymerase I and identified substantial dynamics in polymerase complexes formed prior to nucleotide incorporation; these dynamics may be important for the fidelity of DNA synthesis. We expect BVA to be broadly applicable to single-molecule FRET studies of molecular structure and to complement approaches such as probability distribution analysis and fluorescence correlation spectroscopy in studying molecular dynamics.     

Influence of base-pair changes and cooperativity parameters on the melting curves of short DNAs

Theoretical melting curves were calculated for four DNA restriction fragments, 157–257 base pairs (bp), and a series of hypothetical block DNAs with sequences d(C_2xA_xC_2x). d(C_2xT_xG_2x), 5 ? x ? 40. These DNAs provided a mixture of A·T/G·C sequence distributions with which to investigate the effects of parameters and base-pair changes on the melting of short DNAs. The sensitivity of DNA melting curves to changes in internal loop melting parameters σ and κ was examined. As Expected, theoretical melting curves of short DNAs with a quasirandom base-pair sequence vary little with changes in internal loop parameters. End melting dominates the transition behaviour of these moleucles. This was also observed for the block DNAs up to x = 22. Beyond this length, melting curves are highly sensitive to the internal loop parameters. Sensitivity is also predicted for a 157-bp fragment with a block distribution of A·T and G·C pairs. These results indicate that accurate evaluation of internal loop parameters is possible with short DNAs (100–200 bp) containing a G·C/A·T/G·C block distribution with at least 22 bp in each block. Duplex-to-single-strands dissociation parameters were reevaluated form experimental melting curve data of eight DNA fragments using a least squares fit approach. This analysis confirmed parameter values previously found with a simplified dissociation model. A Priori predictions are made on the effects of base-pair changes on the melting curves of three characterized DNA restriction fragments. Single base-pair changes are predicted to induce small but measurable changes in the melting curves. The characteristics of the altered melting curves depend on the location of the base-pair change.     

Solution structure of duplex DNA containing an extrahelical abasic site analog determined by NMR spectroscopy and molecular dynamics.

Translesional DNA synthesis past abasic sites proceeds with the preferential incorporation of dAMP opposite the lesion and, depending on the sequence context, one or two base deletions. High-resolution NMR spectroscopy and molecular dynamics simulations were used to determine the three-dimensional structure of a DNA heteroduplex containing a synthetic abasic site (tetrahydrofuran) residue positioned in a sequence that promotes one base deletions. Analysis of NMR spectra indicates that the stem region of the duplex adopts a right-handed helical structure and the glycosidic torsion angle is in anti orientation for all residues. NOE interactions establish Watson-Crick alignments for all canonical base pairs of the duplex. Measurement of distance interactions at the lesion site shows the abasic residue excluded from the helix. Restrained molecular dynamics simulations generated three-dimensional models in excellent agreement with the spectroscopic data. These structures show a regular duplex region and a slight bend at the lesion site. The tetrahydrofuran residue extrudes from the helix and is highly flexible. The model reported here, in conjunction with a previous study performed on abasic sites, explains the structural bias of one-base deletion mutations.     

Binding of Net-Fla,a Netropsin-Flavin Hybrid Molecule,to DNA: Molecular Mechanics and Dynamics Studies in vacuo and in Water Solution

Abstract

We have studied the binding of the hybrid netropsin-flavin (Net-Fla) molecule onto four sequences containing four A.T base pairs. Molecular mechanics minimizations in vacuo show numerous minimal conformations separated by one base pair. 400 ps molecular dynamics simulations in vacuo have been performed using the lowest minima as the starting conformations. During these simulations, the flavin moiety of the drug makes two hydrogen bonds with an amino group of a neighboring guanine. A 200 ps molecular dynamics simulation in explicit water solution suggests that the binding of Net-Fla upon the DNA substrate is enhanced by water bridges. A water molecule bridging the amidinium of Net-Fla to the N3 atom of an adenine seems to be stuck in the dmg-DNA complex during the whole simulation. The fluctuations of the DNA helical parameters and of the torsion angles of the sugar-phosphate backbone are very similar in the simulations in vacuo and in water. The time auto-correlation functions for the DNA helical parameters decrease rapidly in the picosecond range in vacuo. The same functions computed from the water solution molecular dynamics simulations seem to have two modes: the rapid mode is similar to the behavior in vacuo, and is followed by a slower mode in the 10 ps range.     

Dynamical allosterism in the mechanism of action of DNA mismatch repair protein MutS

The multidomain protein Thermus aquaticus MutS and its prokaryotic and eukaryotic homologs recognize DNA replication errors and initiate mismatch repair. MutS actions are fueled by ATP binding and hydrolysis, which modulate its interactions with DNA and other proteins in the mismatch-repair pathway. The DNA binding and ATPase activities are allosterically coupled over a distance of ∼70 Å, and the molecular mechanism of coupling has not been clarified. To address this problem, all-atom molecular dynamics simulations of ∼150 ns including explicit solvent were performed on two key complexes—ATP-bound and ATP-free MutS⋅DNA(+T bulge). We used principal component analysis in fluctuation space to assess ATP ligand-induced changes in MutS structure and dynamics. The molecular dynamics-calculated ensembles of thermally accessible structures showed markedly small differences between the two complexes. However, analysis of the covariance of dynamical fluctuations revealed a number of potentially significant interresidue and interdomain couplings. Moreover, principal component analysis revealed clusters of correlated atomic fluctuations linking the DNA and nucleotide binding sites, especially in the ATP-bound MutS⋅DNA(+T) complex. These results support the idea that allosterism between the nucleotide and DNA binding sites in MutS can occur via ligand-induced changes in motion, i.e., dynamical allosterism.     

Differential homogenization and amplification of two satellite DNAs in the genus Cucurbita (Cucurbitaceae)

Two different satellite DNAs exist in the genus Cucurbita which are different with respect to repeat length (350 by and 170 bp), array size, and sequence homogenization. Whereas the 350-bp satellite DNA is prominent and very homogeneous in all species investigated except for C. maxima and C. lundelliana, the 170-bp satellite is rather evenly distributed in all species. In C. maxima and C. lundelliana the 350-bp satellite is present only in small amounts, but detectable by the sensitive PCR method. These repeats are also very homogeneous, reflecting a silent stage of satellite DNA. In contrast, the 170-bp satellite DNA is intra- and interspecifically heterogeneous. It is striking that the species with no detectable amount of 350-bp satellite contain 170-bp satellite DNA clusters with the highest degree of homogeneity. The evolution of satellite DNA repeats within cultivated and wild species in the genus Cucurbita is elucidated using the sequence data of both satellite DNAs from all species investigated. The value of satellite DNA for phylogenetic analysis between closely related species is discussed. Correspondence to: V. Hemleben     

Thermodynamics of site-specific small molecular ion interactions with DNA duplex: a molecular dynamics study

The stability and dynamics of a double-stranded DNA (dsDNA) is affected by the preferential occupancy of small monovalent molecular ions. Small metal and molecular ions such as sodium and alkyl ammonium have crucial biological functions in human body, affect the thermodynamic stability of the duplex DNA and exhibit preferential binding. Here, using atomistic molecular dynamics simulations, we investigate the preferential binding of metal ion such as Na⁺ and molecular ions such as tetramethyl ammonium (TMA⁺) and 2-hydroxy-N,N,N-trimethylethanaminium (CHO⁺) to double-stranded DNA. The thermodynamic driving force for a particular molecular ion-DNA interaction is determined by decomposing the free energy of binding into its entropic and enthalpic contributions. Our simulations show that each of these molecular ions preferentially binds to the minor groove of the DNA and the extent of binding is highest for CHO⁺. The ion binding processes are found to be entropically favourable. In addition, the contribution of hydrophobic effects towards the entropic stabilisation (in case of TMA⁺) and the effect of hydrogen bonding contributing to enthalpic stabilisation (in case of CHO⁺) have also been investigated.     

Molecular structure of the lampbrush loops nooses of the Y chromosome of Drosophila hydei

The molecular structure of the lampbrush loopforming fertility gene nooses from the short arm of the Y chromosome of Drosophila hydei is described on the basis of cloned DNA sequences which are characteristic for the sequence organization in the lampbrush loop. Y chromosomal lampbrush loops are organized into tandem repeat clusters of loop-specific repetitive DNA sequences and in interspersed repetitive DNA sequences with homologies elsewhere in the genome. In this paper, the basic properties of a repeat unit of the tandemly repeated sequence family ay1 are described. Moreover, it is shown that a loop contains several different domains carrying repeat clusters of the same repeated DNA family but with divergent sequence character. One of these clusters is characterized by an internal duplication of the basic repeat unit. We propose that the tandem repeat DNA family ay1 forms a frame of the lampbrush loop which is required for structural and functional reasons.     

Fine-tuning of Voltage Sensitivity of the Kv1.2 Potassium Channel by Interhelix Loop Dynamics

Many proteins function by changing conformation in response to ligand binding or changes in other factors in their environment. Any change in the sequence of a protein, for example during evolution, which alters the relative free energies of the different functional conformations changes the conditions under which the protein will function. Voltage-gated ion channels are membrane proteins that open and close an ion-selective pore in response to changes in transmembrane voltage. The charged S4 transmembrane helix transduces changes in transmembrane voltage into a change in protein internal energy by interacting with the rest of the channel protein through a combination of non-covalent interactions between adjacent helices and covalent interactions along the peptide backbone. However, the structural basis for the wide variation in the V₅₀ value between different voltage-gated potassium channels is not well defined. To test the role of the loop linking the S3 helix and the S4 helix in voltage sensitivity, we have constructed a set of mutants of the rat K_v1.2 channel that vary solely in the length and composition of the extracellular loop that connects S4 to S3. We evaluated the effect of these different loop substitutions on the voltage sensitivity of the channel and compared these experimental results with molecular dynamics simulations of the loop structures. Here, we show that this loop has a significant role in setting the precise V₅₀ of activation in K_v1 family channels.     

Theory of DNA melting curves

Exact algorithms for the calculation of melting curves of heterogeneous DNA with N base pairs apparently require computer time proportional to N². However, it is shown that a decomposition of the loop entropy factor into a sum of I exponential functions (1) gives an extremely accurate approximation to the loop entropy factor for small values of I and (2) makes the computer time for the exact algorithms proportional to I·N. In effect, exact results for melting curves and lengths of helix or coil stretches are obtained with computer time comparable to that required for the Frank-Kamenetskii approximation. The remarkable accuracy of the latter for the fraction of helical content (errors of 0.01–0.05) is confirmed, but appreciably larger errors are found for the lengths of helix or coil stretches (typical errors of 30–100%).     

Molecular docking and dynamics simulations on the interaction of cationic porphyrin–anthraquinone hybrids with DNA G-quadruplexes

A series of cationic porphyrin–anthraquinone hybrids bearing either pyridine, imidazole, or pyrazole rings at the meso-positions have been investigated for their interaction with DNA G-quadruplexes by employing molecular docking and molecular dynamics simulations. Three types of DNA G-quadruplexes were utilized, which comprise parallel, antiparallel, and mixed hybrid topologies. The porphyrin hybrids have a preference to bind with parallel and mixed hybrid structures compared to the antiparallel structure. This preference arises from the end stacking of porphyrin moiety following G-stem and loop binding of anthraquinone tail, which is not found in the antiparallel due to the presence of diagonal and lateral loops that crowd the G-quartet. The binding to the antiparallel, instead, occurred with poorer affinity through both the loop and wide groove. All sites of porphyrin binding were confirmed by 6 ns molecular dynamics simulation, as well as by the negative value of the total binding free energies that were calculated using the MMPBSA method. Free energy analysis shows that the favorable contribution came from the electrostatic term, which supposedly originated from the interaction of either cationic pyridinium, pyrazole, or imidazole groups and the anionic phosphate backbone, and also from the van der Waals energy, which primarily contributed through end stacking interaction.     

Homology modeling and molecular dynamics simulations of HgiDII methyltransferase in complex with DNA and S-adenosyl-methionine: Catalytic mechanism and interactions with DNA

M.HgiDII is a methyltransferase (MTase) from Herpetosiphon giganteus that recognizes the sequence GTCGAC. This enzyme belongs to a group of MTases that share a high degree of amino acid similarity, albeit none of them has been thoroughly characterized. To study the catalytic mechanism of M.HgiDII and its interactions with DNA, we performed molecular dynamics simulations with a homology model of M.HgiDII complexed with DNA and S-adenosyl-methionine. Our results indicate that M.HgiDII may not rely only on Glu119 to activate the cytosine ring, which is an early step in the catalysis of cytosine methylation; apparently, Arg160 and Arg162 may also participate in the activation by interacting with cytosine O2. Another residue from the catalytic site, Val118, also played a relevant role in the catalysis of M.HgiDII. Val118 interacted with the target cytosine and kept water molecules from accessing the region of the catalytic pocket where Cys79 interacts with cytosine, thus preventing water-mediated disruption of interactions in the catalytic site. Specific recognition of DNA was mediated mainly by amino acids of the target recognition domain, although some amino acids (loop 80–88) of the catalytic domain may also contribute to DNA recognition. These interactions involved direct contacts between M.HgiDII and DNA, as well as indirect contacts through water bridges. Additionally, analysis of sequence alignments with closely related MTases helped us to identify a motif in the TRD of M.HgiDII that may be relevant to specific DNA recognition.     

Molecular dynamics simulations of a helicase

Helicases are ubiquitous enzymes involved in nucleic acid metabolism. The PcrA DNA helicase is an essential bacterial protein involved in rolling circle plasmid replication and DNA repair. Recent crystal structures of PcrA bound to DNA indicate that a flexible loop mediates a functionally important rigid-body-domain rotation. In this study, we report stochastic boundary molecular dynamics simulations focused on this region for wild-type and mutants designed to increase the rigidity of the region. Residues in the region that were helix-disfavoring, such as glycine, threonine, and others, were mutated to alanine. The simulated dynamics, analyzed with a variety of measures of structure and mobility, indicate that a few point mutations will substantially increase helix formation in this region. Subnanosecond stochastic boundary molecular dynamics simulations at several temperatures offer a rapid protocol for assessing large numbers of mutants and provides a novel strategy for the design of experiments to test the role of this flexible loop region in the function of PcrA.