Search PubMed⌕ Search

Biomedical subjects

T Schlick

Publications and source records attributed to T Schlick.

25 records · Page 2Linked to original sources
PubMed1995

Modeling superhelical DNA: recent analytical and dynamic approaches.

During the past year, a variety of diverse and complementary approaches have been presented for modeling superhelical DNA, offering new physical and biological insights into fundamental functional processes of DNA. Analytical approaches have probed deeper into the effects of entropy and thermal fluctuations on DNA structure and on various topological constraints induced by DNA-binding proteins. In tandem, new kinetic approaches--by molecular, Langevin and Brownian dynamics, as well as extensions of elastic-rod theory--have begun to offer dynamic information associated with supercoiling. Such dynamic approaches, along with other equilibrium studies, are refining the basic elastic-rod and polymer framework and incorporating more realistic treatments of salt and sequence-specific features. These collective advances in modeling large DNA molecules, in concert with technological innovations, are pointing to an exciting interplay between theory and experiment on the horizon.

Computer Graphics↗
PubMed1994

On higher buckling transitions in supercoiled DNA.

A combination of detailed energy minimization and molecular dynamics studies of closed circular DNA offers here new information that may be relevant to the dynamics of short DNA chains and/or low superhelical densities. We find a complex dependence of supercoiled DNA energies and geometries on the linking number difference delta Lk as physiological superhelical densities (magnitude of sigma approximately 0.06) are approached. The energy minimization results confirm and extend predictions of classical elasticity theory for the equilibria of elastic rods. The molecular dynamics results suggest how these findings may affect the dynamics of supercoiled DNA. The minimization reveals sudden higher order configurational transitions in addition to the well-known catastrophic buckling from the circle to the figure-8. The competition among the bending, twisting, and self-contact forces leads to different families of supercoiled forms. Some of those families begin with configurations of near-zero twist. This offers the intriguing possibility that nicked DNA may relax to low-twist forms other than the circle, as generally assumed. Furthermore, for certain values of delta Lk, more than one interwound DNA minimum exists. The writhing number as a function of delta Lk is discontinuous in some ranges; it exhibits pronounced jumps as delta Lk is increased from zero, and it appears to level off to a characteristic slope only at higher values of delta Lk. These findings suggest that supercoiled DNA may undergo systematic rapid interconversions between different minima that are both close in energy and geometry. Our molecular dynamics simulations reveal such transitional behavior. We observe the macroscopic bending and twisting fluctuations of interwound forms about the global helix axis as well as the end-over-end tumbling of the DNA as a rigid body. The overall mobility can be related to magnitude of sigma and to the bending, twisting, and van der Waals energy fluctuations. The general character of molecular motions is thus determined by the types of energy minima found at a given delta Lk. Different time scales may be attributed to each type of motion: The overall chain folding occurs on a time scale almost an order of magnitude faster than the end-over-end tumbling. The local bending and twisting of individual chain residues occur at an even faster rate, which in turn correspond to several cycles of local variations for each large-scale bending and straightening motion of the DNA.

DNA, Superhelical↗
PubMed1994

The influence of salt on the structure and energetics of supercoiled DNA.

We present a detailed computational study of the influence of salt on the configurations, energies, and dynamics of supercoiled DNA. A potential function that includes both elastic and electrostatic energy components is employed. Specifically, the electrostatic term, with salt-dependent coefficients, is modeled after Stigter's pioneering work on the effective diameter of DNA as a function of salt concentration. Because an effective charge per unit length is used, the electrostatic formulation does not require explicit modeling of phosphates and can be used to study long DNAs at any desired resolution of charge. With explicit consideration of the electrostatic energy, an elastic bending constant corresponding to the nonelectrostatic part of the bending contribution to the persistence length is used. We show, for a series of salt concentrations ranging from 0.005 to 1.0 M sodium, how configurations and energies of supercoiled DNA (1000 and 3000 base pairs) change dramatically with the simulated salt environment. At high salt, the DNA adopts highly compact and bent interwound states, with the bending energy dominating over the other components, and the electrostatic energy playing a minor role in comparison to the bending and twisting terms. At low salt, the DNA supercoils are much more open and loosely interwound, and the electrostatic components are dominant. Over the range of three decades of salt examined, the electrostatic energy changes by a factor of 10. The buckling transition between the circle and figure-8 is highly sensitive to salt concentration: this transition is delayed as salt concentration decreases, with a particularly sharp increase below 0.1 M. For example, for a bending-to-twisting force constant ratio of A/C = 1.5, the linking number difference (delta LK) corresponding to equal energies for the circle and figure-8 increases from 2.1 to 3.25 as salt decreases from 1.0 to 0.005 M. We also present in detail a family of three-lobed supercoiled DNA configurations that are predicted by elasticity theory to be stable at low delta Lk. To our knowledge, such three-dimensional structures have not been previously presented in connection with DNA supercoiling. These branched forms have a higher bending energy than the corresponding interwound configurations at the same delta Lk but, especially at low salt, this bending energy difference is relatively small in comparison with the total energy, which is dominated by the electrostatic contributions. Significantly, the electrostatic energies of the three-lobed and (straight) interwound forms are comparable at each salt environment. We show how the three-lobed configurations change slowly with ALk, resulting in branched interwound forms at higher salt. In longer chains, the branched forms are highly interwound, with bent arms. At low salt, the branched supercoils are asymmetric, with a longer interwound stem and two shorter arms. From molecular dynamics simulations we observe differences in the motions of the DNA as a function of salt. At high salt, the supercoiled chain is quite compact but fairly rigid, whereas at low salt the DNA is loosely coiled but more dynamic. Especially notable at low salt are the large-scale opening and closing of the chain as a whole and the rapid "slithering"of individual residues past one another. Toroidal forms are not detected under these conditions. However, the overall features of the open, loose supercoils found at low salt are more similar to those of toroidal than interwound configurations. Indeed,simulated x-ray scattering profiles reveal the same trends observed experimentally and are consistent with a change from closed to open forms as salt is decreased. Like the minimization studies, the dynamics reveal a critical point near 0.1 M associated with the collapse of loose to tight supercoils. Near this physiological concentration, enhanced flexibility of the DNA is noted. The collective observations suggest a potential regulatory role for salt on supercoiled DNA function, not only for closed circular DNA,but also for linear DNA with small looped regions.

Bacteriophage phi X 174↗
PubMed1992

Trefoil knotting revealed by molecular dynamics simulations of supercoiled DNA.

Computer simulations of the supercoiling of DNA, largely limited to stochastic search techniques, can offer important information to complement analytical models and experimental data. Through association of an energy function, minimum-energy supercoiled conformations, fluctuations about these states, and interconversions among forms may be sought. In theory, the observation of such large-scale conformational changes is possible, but modeling and numerical considerations limit the picture obtained in practice. A new computational approach is reported that combines an idealized elastic energy model, a compact B-spline representation of circular duplex DNA, and deterministic minimization and molecular dynamics algorithms. A trefoil knotting result, made possible by a large time-step dynamics scheme, is described. The simulated strand passage supports and details a supercoiled-directed knotting mechanism. This process may be associated with collective bending and twisting motions involved in supercoiling propagation and interwound branching. The results also demonstrate the potential effectiveness of the Langevin/implicit-Euler dynamics scheme for studying biomolecular folding and reactions over biologically interesting time scales.

Chemical Phenomena↗
PubMed1992

Supercoiled DNA energetics and dynamics by computer simulation.

A new formulation is presented for investigating supercoiled DNA configurations by deterministic techniques. Thus far, the computational difficulties involved in applying deterministic methods to supercoiled DNA studies have generally limited computer simulations to stochastic approaches. While stochastic methods, such as simulated annealing and Metropolis-Monte Carlo sampling, are successful at generating a large number of configurations and estimating thermodynamic properties of topoisomer ensembles, deterministic methods offer an accurate characterization of the minima and a systematic following of their dynamics. To make this feasible, we model circular duplex DNA compactly by a B-spline ribbon-like model in terms of a small number of control vertices. We associate an elastic deformation energy composed of bending and twisting integrals and represent intrachain contact by a 6-12 Lennard Jones potential. The latter is parameterized to yield an energy minimum at the observed DNA-helix diameter inclusive of a hydration shell. A penalty term to ensure fixed contour length is also included. First and second partial derivatives of the energy function have been derived by using various mathematical simplifications. First derivatives are essential for Newton-type minimization as well as molecular dynamics, and partial second-derivative information can significantly accelerate minimization convergence through preconditioning. Here we apply a new large-scale truncated-Newton algorithm for minimization and a Langevin/implicit-Euler scheme for molecular dynamics. Our truncated-Newton method exploits the separability of potential energy functions into terms of differing complexity. It relies on a preconditioned conjugate gradient method that is efficient for large-scale problems to solve approximately for the search direction at every step. Our dynamics algorithm is numerically stable over large time steps. It also introduces a frequency-discriminating mechanism so that vibrational modes with frequencies greater than a chosen cutoff frequency are essentially frozen by the method. With these tools, we rapidly identify corresponding circular and interwound energy minima for small DNA rings for a series of imposed linking-number differences. These structures are consistent with available electron microscopy data. The energetic exchange of stability between the circle and the figure-8, in very good agreement with analytical results, is also detailed. Molecular dynamics trajectories at 100 femtosecond time steps then reveal the rapid folding of the unstable circular state into supercoiled forms. Significant bending and twisting motions of the interwound structures are also observed. Such information may be useful for understanding transition states along the folding pathway and the role of enzymes that regulate supercoiling.(ABSTRACT TRUNCATED AT 400 WORDS)

Algorithms↗
PubMedDate not supplied

An efficient projection protocol for chemical databases: singular value decomposition combined with truncated-newton minimization.

A rapid algorithm for visualizing large chemical databases in a low-dimensional space (2D or 3D) is presented as a first step in database analysis and design applications. The projection mapping of the compound database (described as vectors in the high-dimensional space of chemical descriptors) is based on the singular value decomposition (SVD) combined with a minimization procedure implemented with the efficient truncated-Newton program package (TNPACK). Numerical experiments on four chemical datasets with real-valued descriptors (ranging from 58 to 27 255 compounds) show that the SVD/TNPACK projection duo achieves a reasonable accuracy in 2D, varying from 30% to about 100% of pairwise distance segments that lie within 10% of the original distances. The lowest percentages, corresponding to scaled datasets, can be made close to 100% with projections onto a 10-dimensional space. We also show that the SVD/TNPACK duo is efficient for minimizing the distance error objective function (especially for scaled datasets), and that TNPACK is much more efficient than a current popular approach of steepest descent minimization in this application context. Applications of our projection technique to similarity and diversity sampling in drug design can be envisioned.

Database Management Systems↗