DNA curvature influences the internal motions of supercoiled DNA.

We present evidence that short curved DNA segments can act as mediators for the ordering of large domains in superhelical DNA. Using a non-invasive solution method (dynamic light scattering), we investigated the effect of permanently curved inserts on the solution structure and on the internal motions of superhelical plasmid DNA. We find that the dynamics of superhelical DNA are strongly influenced by sequence- or protein-induced bending: in superhelical plasmids containing curved inserts the amplitude of the internal motion is lower than that of non-curved controls. Furthermore, the relative arrangement of curved sequences in the plasmids can influence the overall shape of the superhelical DNA. On linearized forms of the plasmids, these effects are not observed.     

Origin of low-frequency motions in biological macromolecules. A view of recent progress in the quasi-continuity model

The recent progress in the quasi-continuity model and its applications in studying the low-frequency internal motions of biological macromolecules have been surveyed. Emphasis is placed on revealing the origin of this kind of internal collective motion, which involves many atoms and has significant biological functions. In light of such a line, the low-frequency motions in alpha-helix structure, beta-structure (including beta-sheet and beta-barrel), and DNA double-helix structure, the three most fundamental component elements in biological macromolecules, are discussed, and the corresponding physical pictures described. It turns out that the low-frequency motion in biological macromolecules originates from their two common intrinsic characteristics, i.e., they possess a series of weak bonds, such as hydrogen bonds and salt bridges, and a substantial mass distributed over the region containing those weak bonds.     

Anisotropic overall and internal motions of short DNA fragments.

Anisotropic motions of DNA fragments in the size range 6-118 base pairs are studied by the steady-state fluorescence polarization of different excitation transitions in the intercalated ethidium cation. Calculated effective tumbling and twisting times are found to be shorter than predicted for overall motions of rigid DNA, indicating that internal motions and/or dye wobbling contribute to the depolarization. The data are consistent with a model where the DNA fragments are considered to be rigid against bending but torsionally flexible, and where the dye can wobble within the intercalated site. We also discuss the possibility of correlated out-of-plane motions of the dye and the DNA bases.     

Simplified force field calculations of the low-frequency motions of the alpha-helix

Using valence force field constants and a highly simplified model of the α-helix consisting of only two masses per repeat unit we have been able to very nearly reproduce the low frequency phonon dispersion curves which we had obtained earlier with the use of a seven mass model and a Urey-Bradley force field. Since the frequencies obtained from this latter model agree well with the observed Raman frequencies it appears that the low frequency chain motions of polypeptide chains may be calculated from a very simple model involving the amino acid residues as essentially a single mass. Such techniques appear to offer the possibility of making normal mode calculations for the chain dynamics of polypeptide chains of low symmetry such as proteins.     

Theory of low-frequency vibrations in DNA macromolecules

To describe low-frequency dynamics of DNA macromolecules a model is developed taking into account the hydrogen bond stretching in base pairs, the backbone flexibility and intranucleoside mobility. For double-stranded DNA the normal vibrations are found and the structure of low-frequency spectrum is determined. The agreement between theory and Raman spectroscopy data for B-DNA is demonstrated. Conformational dependences of vibration spectrum during the B----A and helix----coil DNA transitions are studied. The contribution coming from low-frequency mobility to the nucleic-protein recognition processes is discussed.     

Phosphorus-31 nuclear magnetic resonance of highly oriented DNA fibers. 2. Molecular motions in hydrated DNA

31P nuclear magnetic resonances (NMR) of salmon sperm DNA, poly(rA).poly(rU), and poly(rA).poly(dT) fibers were measured as a function of relative humidity. The results indicated that the spectra were strongly perturbed by the molecular motions occurring in the hydrated fibers. The humidity dependence of the spectra at a number of orientations of the fibers relative to the magnetic field was reasonably explained by taking into account at least three motional modes, namely, conformational fluctuations, restricted rotation about a tilted axis, and rotational diffusion about the helical axis. The rotational diffusion about the helical axis was found to perturb the spectral line shapes most strongly, and its constants were 1.5 X 10(4) and 5.0 X 10(4) S-1 for DNA fibers at 92% and 98% relative humidities, respectively. A DNA-RNA hybrid, poly(rA).poly(dT), has been shown to adopt different conformations on two strands at high relative humidity [Zimmerman, S. B., & Pheiffer, B. H. (1981) Proc. Natl. Acad. Sci. U.S.A. 78, 78-82], which was unquestionably confirmed in the present study: that is, the 31P NMR spectra from the hydrated form of this polymer were clearly explained by assuming that one strand had an A-like conformation and the other a B-like conformation.     

Dynamic light scattering studies of internal motions in DNA. II. Clean viral DNAs

Internal Brownian motions of clean ?29 and λ-DNAs have been studied using photon-correlation techniques at both visible (λ₀ = 632.8 nm) and uv (λ₀ = 363.8 nm) wavelengths. The present dynamic light scattering data, which extend to K² = 19 × 10¹⁰ cm^?2, can in every case be satisfactorily simulated by a Rouse-Zimm model polymer with an appropriate choice of the three model parameters. The effects of pH, salt concentration, single-strand breaks, and molecular weight on those model parameters have also been investigated. Intact clean DNAs exhibit surprisingly little variation with pH from 7.85 to 10.25, with salt concentration from 0.01 NaCl to 5.4M NH₄Cl, or with molecular weight or GC content. The single-strand breaks have no effect at pH 9.46, but produce dramatic changes in the model parameters at pH 10.0 and 10.25, indicating the introduction of titratable joints at those pHs. The failure of either single-strand breaks or a large change in GC content to alter the model parameters in the neutral pH range is a strong indication that local denaturation is not required for those flexions and torsions that dominate the relaxation of fluctuations in the scattered light. The Langevin relaxation time for the slowest internal mode of a particular Rouse-Zimm model derived from the dynamic light scattering data is compared with pertinent literature data extrapolated to the same molecular weight. The present algorithm for determining model parameters from the light-scattering D_app vs K² curve actually yields a Langevin time in fairly good agreement with the literature value. For unknown reasons the light-scattering D₀ values generally exceed those obtained from the molecular weight and sedimentation coefficient by about 20%.     

Filtering molecular dynamics trajectories to reveal low-frequency collective motions: phospholipase A2

A novel method for analysing molecular dynamics trajectories has been developed, which filters out high frequencies using digital signal processing techniques and facilitates focusing on the low-frequency collective motions of proteins. These motions involve low energy slow motions, which lead to important biological phenomena such as domain closure and allosteric effects in enzymes. The filtering method treats each of the atomic trajectories obtained from the molecular dynamics simulation as a "signal". The trajectories of each of the atoms in the system (or any subset of interest) are Fourier transformed to the frequency domain, a filtering function is applied and then an inverse transformation back to the time domain yields the filtered trajectory. The filtering method has been used to study the dynamics of the enzyme phospholipase A2. In the filtered trajectory, all the high frequency bond and valence angle vibrations were eliminated, leaving only low-frequency motion, mainly fluctuations in torsions and conformational transitions. Analysis of this trajectory revealed interesting motions of the protein, including concerted movements of helices, and changes in shape of the active site cavity. Unlike normal mode analysis, which has been used to study the motion of proteins, this method does not require converged minimizations or diagonalization of a matrix of second derivatives. In addition, anharmonicity, multiple minima and conformational transitions are treated explicitly. Thus, the filtering method avoids most of the approximations implicit in other investigations of the dynamic behaviour of large systems.     

A fluorescence photobleaching study of the microsecond reorientational motions of DNA.

We have conducted a polarized fluorescence photobleaching recovery (FPR) study of the rotational dynamics of ethidium azide labeled DNA. Polarized photobleaching experiments provide data on microsecond and millisecond molecular reorientation that complement the information available from nanosecond fluorescence depolarization studies. In polarized FPR experiments an anisotropic angular concentration of fluorophore is created by bleaching dye molecules in a preferred orientation with a short, intense pulse of polarized light. The sample is then weakly illuminated, and the temporal variation in the emitted fluorescence is monitored. The fluorescence signal will systematically change as molecules undergo post-bleach reorientation and the angular distribution of dye tends toward isotropy. We have observed that the time dependence of our microsecond FPR curves is also determined in part by nonrotational phenomena. To isolate the reorientational recovery we conduct our FPR experiments in two modes (called parallel and perpendicular) that differ only in the polarization of the bleaching light. A quotient function, R(t), is constructed from the data obtained in these two modes; the variation with time of this new quantity is governed solely by processes that are sensitive to the polarization of the incident light (e.g., molecular rotation). It is found experimentally that R(t) remains constant, as expected, for rotationally restricted DNA systems despite a temporal recovery in the parallel and perpendicular FPR curves. We also follow the dynamics of solutions of phage lambda DNA as revealed in the temporal dependence of R(t). This DNA system rotationally relaxes after approximately 100 microseconds and the dye/DNA complex reorients substantially during the 10-microseconds bleach period. Our FPR data are interpreted in terms of dynamic models of DNA motion.     

Rapid segmental and subdomain motions of DNA polymerase beta

DNA polymerase (pol) beta is a two-domain DNA repair enzyme that undergoes structural transitions upon binding substrates. Crystallographic structures indicate that these transitions include movement of the amino-terminal 8-kDa lyase domain relative to the 31-kDa polymerase domain. Additionally, a polymerase subdomain moves toward the nucleotide-binding pocket after nucleotide binding, resulting in critical contacts between alpha-helix N and the nascent base pair. Kinetic and structural characterization of pol beta has suggested that these conformational changes participate in stabilizing the ternary enzyme-substrate complex facilitating chemistry. To probe the microenvironment and dynamics of both the lyase domain and alpha-helix N in the polymerase domain, the single native tryptophan (Trp-325) of wild-type enzyme was replaced with alanine, and tryptophan was strategically substituted for residues in the lyase domain (F25W/W325A) or near the end of alpha-helix N (L287W/W325A). Influences of substrate on the fluorescence anisotropy decay of these single tryptophan forms of pol beta were determined. The results revealed that the segmental motion of alpha-helix N was rapid ( approximately 1 ns) and far more rapid than the step that limits chemistry. Binding of Mg(2+) and/or gapped DNA did not cause a noticeable change in the rotational correlation time or angular amplitude of tryptophan in alpha-helix N. More important, binding of a correct nucleotide significantly limited the angular range of the nanosecond motion within alpha-helix N. In contrast, the segmental motion of the 8-kDa domain was "frozen" upon DNA binding alone, and this restriction did not increase further upon nucleotide binding. The dynamics of alpha-helix N are discussed from the perspective of the "open" to "closed" conformational change of pol beta deduced from crystallography, and the results are more generally discussed in the context of reaction cycle-regulated flexibility for proteins acting as molecular motors.     

Dynamic light-scattering studies of internal motions in DNA. I. Applicability of the rouse-zimm model

The intermediate scattering function G(K,t) for any polymer model obeying a linear separable Langevin equation can be expressed in terms of the eigenvalues and eigenvectors of its normal coordinate transformation. An algorithm for the extract numerical evaluation of G(K,t) for linear Rouse-Zimm chains in the presence of hydrodynamic interaction has been developed. The computed G(K,t)² were fit to C(t) = A exp(?t/τA) + B, and apparent diffusion coefficients calculated according to D_app ≡ 1/(2τ_AK²). G(K,t)² was surprisingly well-fit by single-exponential decays, especially at both small and large values of Kb, where K is the scattering vector and b the root-mean-squared subunit extension. Plots of D_app vs K² in-variably showed a sigmoidal rise from D₀ at K² = O up to a constant plateau value at large K²b². Analytical expression for G(K,t), exact in the limit of short times, were obtained for circular Rouse-Zimm chains with and without hydrodynamic interaction, and also for free-draining linear chains, and in addition for the independent-segment-mean-force (ISMF) model. The predicted behaviors for G(K,t) at large Kb (or KR_G) was found in all cases to be single-exponential with 1/τ ∝ K² at large Kb, in agreement with the computational results. A simple procedure for estamating all parameter of the Rouse-Zimm model from a plot of D_app vs K² is proposed. Experimental data for both native and pH-denatured calf-thymus DNA in 1.0M Nacl with and without EDTA clearly plateau behavior of D_app at large values of K, in harmony with the present Rouse-Zimm and ISMF theories, and in sharp contrast to previous predictions based on the Rouse-Zimm model.     

Structural models for non-helical DNA.

Structural modelling techniques are employed to explore the energetic requirements for the transformation of classical B DNA into unwound yet double-stranded DNA structures. Structural idealization using CORELS computer program of Sussman et al. followed by energy minimization using the EREF program of Levitt, leads to two regular non-helical models. In both models, the bases are conventionally paired and stacked, yet there is no net rotation between successive base pairs. One model, N1, has a 1-bp repeating unit; the second, N2, has a 2-bp repeating unit. The dihedral angles of the backbone all have values found either in the B or the Z form of DNA, except for the P-O5'-C5'-C4' angle, which is in the unprecedented g+ or g- domains. The energy difference found between the two N form models and B form DNA are 6.6 and 3.4 kcal/mol/nucleotide for N1 and N2 respectively. These relatively low energy differences encourage the idea that non-helical forms of DNA may contribute to the alternate DNA structures found in S1 nuclease sensitive and other regulatory regions of active genes.     

DNA polymerase alpha and models for proofreading.

Using a modified system to measure fidelity at an amber site in phi X174, we have employed DNA polymerase alpha to test different mechanisms for proofreading. DNA polymerase alpha does not exhibit the characteristics of "kinetic proofreading" seen with procaryotic polymerases. Polymerase alpha shows no evidence for a "next nucleotide" effect, and added deoxynucleoside monophosphates do not alter fidelity. Pyrophosphate, which increases error rates with a procaryotic polymerase, appears to weakly improve polymerase alpha fidelity. DNA polymerase alpha does exhibit a dramatic increase in error rate in the presence of a deoxycytidine thiotriphosphate (dCTP alpha S), but this enhanced mutagenesis also occurs under conditions where kinetic proofreading should be otherwise defeated. This particular effect with dCTP alpha S appears specific for DNA polymerase alpha and is not seen with the other polymerases tested.     

Charge transport in DNA model with vibrational and rotational coupling motions

The dynamics of the Peyrard-Bishop model for vibrational motion of DNA dynamics, which has been extended by taking into account the rotational motion for the nucleotides (Silva et al., J. Biol. Phys. 34, 511–519, 2018) is studied. We report on the presence of the modulational instability (MI) of a plane wave for charge migration in DNA and the generation of soliton-like excitations in DNA nucleotides. We show that the original differential-difference equation for the DNA dynamics can be reduced in the continuum approximation to a set of three coupled nonlinear equations. The linear stability analysis of continuous wave solutions of the coupled systems is performed and the growth rate of instability is found numerically. Numerical simulations show the validity of the analytical approach with the generation of wave packets provided that the wave numbers fall in the instability domain.     

Nucleotide-induced DNA polymerase active site motions accommodating a mutagenic DNA intermediate

DNA polymerases occasionally insert the wrong nucleotide. For this error to become a mutation, the mispair must be extended. We report a structure of DNA polymerase beta (pol beta) with a DNA mismatch at the boundary of the polymerase active site. The structure of this complex indicates that the templating adenine of the mispair stacks with the primer terminus adenine while the templating (coding) cytosine is flipped out of the DNA helix. Soaking the crystals of the binary complex with dGTP resulted in crystals of a ternary substrate complex. In this case, the templating cytosine is observed within the DNA helix and forms Watson-Crick hydrogen bonds with the incoming dGTP. The adenine at the primer terminus has rotated into a syn-conformation to interact with the opposite adenine in a planar configuration. Yet, the 3'-hydroxyl on the primer terminus is out of position for efficient nucleotide insertion.     

Internal motions in myosin.

High-resolution proton nuclear magnetic resonance (1H NMR) measurements were made on myosin, heavy meromyosin (HMM), myosin subfragment 1 (S1), light meromyosin (LMM), and actin. A strong signal from amino acid side chains undergoing motions too fast to be accounted for by simple rotations of groups on a rigid backbone was obtained from myosin. Comparison of myosin, HMM, S1, and LMM showed that the mobile region is located almost entirely in S1 and accounts for approximately 22% of its structure. Adenosine triphosphate (ATP) and ATP analogues had no measurable effect on the S1 spectrum. Actin, on the other hand, quenched the internal motions of S1. When S1 was titrated with actin, an association was obtained which was in agreement with other measured values. The actin effect was reversed by adding magnesium pyrophosphate (MgPPi) or adenyl-5'-yl imidophosphate (MgAMPPNP). Quantitative treatment of the broad signals from myosin and its subfragments substantiated the existence of two flexible regions in myosin. The highly mobile portion of myosin may be located in the "swivel" between S1 and the rest of myosin or in the actin binding site or in both. These possibilites are discussed, and a new possible mechanism for muscle cross bridge elasticity is proposed.     

Structure determination of a DNA octamer in solution by NMR spectroscopy. Effect of fast local motions.

NMR structures of biomolecules are primarily based on nuclear Overhauser effects (NOEs) between protons. For the interpretation of NOEs in terms of distances, usually the assumption of a single rotational correlation time corresponding to a rigid molecule approximation is made. Here we investigate the effect of fast internal motions of the interproton vectors in the context of the relaxation matrix approach for structure determination of biomolecules. From molecular dynamics simulations generalized order parameters were calculated for the DNA octamer d(GCGTTCGC).d(CGCAACGC), and these were used in the calculation of NOE intensities. The magnitudes of the order parameters showed some variation for the different types of interproton vectors. The lowest values were observed for the interresidue base H6/H8-H2" proton vectors (S2 = 0.60), while the cytosine H5-H6 interproton vectors were among the most motionally restricted (S2 = 0.92). Inclusion of the motion of the interproton vectors resulted in a much better agreement between theoretically calculated NOE spectra and the experimental spectra measured by 2D NOE spectroscopy. The interproton distances changed only slightly, with a maximum of 10%; nevertheless, the changes were significant and resulted in constraints that were better satisfied. The structure of the DNA octamer was determined by using restrained molecular dynamics simulations with H2O as a solvent, with and without the inclusion of local internal motions. Starting from A- or B-DNA, the structures showed a high local convergence (0.86 A), while the global convergence for the octamer was ca. 2.6 A.