Browsing by Author "Kumar, Rajendra"

The do_x3dna package has been developed to analyze the structural fluctuations of DNA or RNA during molecular dynamics simulations. It extends the capability of the 3DNA package to GROMACS MD trajectories and includes new methods to calculate the global-helical axis of DNA and bending fluctuations during simulations. The package also includes a Python module dnaMD to perform and visualize statistical analyses of complex data obtained from the trajectories.

The DNA packaging motor of the bacteriophage ϕ29, comprising head-tail connector, ATPase, and pRNA, transports the viral DNA inside the procapsid against pressure differences of up to ∼60 atm during replication. Several models for the DNA packaging mechanism have been proposed, which attribute different roles to the connector, and require specific mechanical properties of the connector. To characterize these properties at the atomic level, and to understand how the connector withstands this large pressure, we have carried out molecular dynamics simulations of the whole connector both in equilibrium and under mechanical stress. The simulations revealed a quite heterogeneous distribution of stiff and soft regions, resembling that of typical composite materials that are also optimized to resist mechanical stress. In particular, the conserved middle α-helical region is found to be remarkably stiff, similar only to structural proteins forming viral shell, silk, or collagen. In contrast, large parts of the peripheral interface to the ϕ29 procapsid turned out to be rather soft. Force probe and umbrella sampling simulations showed that large connector deformations are remarkably reversible, and served to calculate the free energies required for these deformations. In particular, for an untwisting deformation by 12°, as postulated by the untwist-twist model, more than four times’ larger energy is required than is available from hydrolysis of one ATP molecule. Combined with previous experiments, this result is incompatible with the untwist-twist model. In contrast, our simulations support the recently proposed one-way revolution model and suggest in structural terms how the connector blocks DNA leakage. In particular, conserved loops at the rim of the central channel, which are in direct contact with the DNA, are found to be rather flexible and tightly anchored to the rigid central region. These findings suggest a check-valve mechanism, with the flexible loops obstructing the channel by interacting with the viral DNA.

During replication of the ϕ29 bacteriophage inside a bacterial host cell, a DNA packaging motor transports the viral DNA into the procapsid against a pressure difference of up to 40 ± 20 atm. Several models have been proposed for the underlying molecular mechanism. Here we have used molecular dynamics simulations to examine the role of the connector part of the motor, and specifically the one-way revolution and the push-roll model. We have focused at the structure and intermolecular interactions between the DNA and the connector, for which a near-complete structure is available. The connector is found to induce considerable DNA deformations with respect to its canonical B-form. We further assessed by force-probe simulations to which extent the connector is able to prevent DNA leakage and found that the connector can act as a partial one-way valve by a check-valve mechanism via its mobile loops. Analysis of the geometry, flexibility, and energetics of channel lysine residues suggested that this arrangement of residues is incompatible with the observed DNA packaging step-size of ∼2.5 bp, such that the step-size is probably determined by the other components of the motor. Previously proposed DNA revolution and rolling motions inside the connector channel are both found implausible due to structural entanglement between the DNA and connector loops that have not been resolved in the crystal structure. Rather, in the simulations, the connector facilitates minor DNA rotation during the packaging process compatible with recent optical-tweezers experiments. Combined with the available experimental data, our simulation results suggest that the connector acts as a check-valve that prevents DNA leakage and induces DNA compression and rotation during DNA packaging.

Experimental and computational dynamic force spectroscopy is widely used to determine the mechanical properties of single biomolecules. Whereas so far the focus has mainly been on rupture or unfolding forces, recent force-probe molecular dynamics simulations have revealed a strong loading rate dependence of biomolecular elasticities, which cannot be explained by the established one-dimensional transition-state treatments. We show that this nonequilibrium behavior can be explained by a theory that includes relaxation effects. For three structurally and mechanically quite diverse systems, a single relaxation mode suffices to quantitatively describe their loading-rate-dependent elastic behavior. Atomistic simulations of these systems revealed the microscopic nature of the respective relaxation modes. This result suggests a new type of "elasticity spectroscopy" experiment, which should render nonequilibrium properties of structured macromolecules accessible to single-molecule force spectroscopy.