Geometric somersaults of a polymer chain through cyclic twisting motions.

This study explores the significance of geometric angle shifts, which we call geometric somersaults, arising from cyclic twisting motions of a polymer chain. A five-bead polymer chain serves as a concise and minimal model of a molecular shaft throughout this study. We first show that this polymer chain can change its orientation about its longitudinal axis largely, e.g., 120^{∘}, under conditions of zero total angular momentum by changing the two dihedral angles in a cyclic manner. This phenomenon is an example of the so-called "falling cat" phenomenon, where a falling cat undergoes a geometric somersault by changing its body shape under conditions of zero total angular momentum. We then extend the geometric somersault of the polymer chain to a noisy and viscous environment, where the polymer chain is steered by external driving forces. This extension shows that the polymer chain can achieve an orientation change keeping its total angular momentum and total external torque fluctuating around zero in a noisy and viscous environment. As an application, we argue that the geometric somersault of the polymer chain by 120^{∘} may serve as a prototypical and coarse-grained model for the rotary motion of the central shaft of ATP synthase (F_{O}F_{1}-ATPase). This geometric somersault is in clear contrast to the standard picture for the rotary motion of the central shaft as a rigid body, which generally incurs nonzero total angular momentum and nonzero total external torque. The power profile of the geometric somersault implies a preliminary mechanism for elastic power transmission. The results of this study may be of fundamental interest in twisting and rotary motions of biomolecules.

Roles of dynamical symmetry breaking in driving oblate-prolate transitions of atomic clusters.

This paper explores the driving mechanisms for structural transitions of atomic clusters between oblate and prolate isomers. We employ the hyperspherical coordinates to investigate structural dynamics of a seven-atom cluster at a coarse-grained level in terms of the dynamics of three gyration radii and three principal axes, which characterize overall mass distributions of the cluster. Dynamics of gyration radii is governed by two kinds of forces. One is the potential force originating from the interactions between atoms. The other is the dynamical forces called the internal centrifugal forces, which originate from twisting and shearing motions of the system. The internal centrifugal force arising from twisting motions has an effect of breaking the symmetry between two gyration radii. As a result, in an oblate isomer, activation of the internal centrifugal force that has the effect of breaking the symmetry between the two largest gyration radii is crucial in triggering structural transitions into prolate isomers. In a prolate isomer, on the other hand, activation of the internal centrifugal force that has the effect of breaking the symmetry between the two smallest gyration radii is crucial in triggering structural transitions into oblate isomers. Activation of a twisting motion that switches the movement patterns of three principal axes is also important for the onset of structural transitions between oblate and prolate isomers. Based on these trigger mechanisms, we finally show that selective activations of specific gyration radii and twisting motions, depending on the isomer of the cluster, can effectively induce structural transitions of the cluster. The results presented here could provide further insights into the control of molecular reactions.

Chiral symmetry breaking of a double-stranded helical chain through bend-writhe coupling.

This paper explores asymmetric elasticity of a double-stranded helical chain, which serves as a minimal model of biopolymers. The model consists of two elastic chains that mutually intertwine in a right-handed manner, forming a double-stranded helix. A simple numerical experiment for structural relaxation, which reduces the total elastic energy of the model monotonically without thermal fluctuations, reveals possible asymmetric elasticity inherent in the helical chain. It is first shown that a short segment of the double-stranded helical chain has a tendency to unwind when it is bent. It is also shown that a short segment of the helical chain has a tendency to writhe in the left direction upon bending. This tendency gives rise to a propensity for a longer segment of the chain to form a left-handed superhelix spontaneously upon bending. Finally, this propensity of the helical chain to form a left-handed superhelix is proposed to be a possible origin of the uniform left-handed wrapping of DNA around nucleosome core particles in nature. The results presented here could provide deeper insights into the roles and significance of helical chirality of biopolymers.

Control of a model of DNA division via parametric resonance.

We study the internal resonance, energy transfer, activation mechanism, and control of a model of DNA division via parametric resonance. While the system is robust to noise, this study shows that it is sensitive to specific fine scale modes and frequencies that could be targeted by low intensity electro-magnetic fields for triggering and controlling the division. The DNA model is a chain of pendula in a Morse potential. While the (possibly parametrically excited) system has a large number of degrees of freedom and a large number of intrinsic time scales, global and slow variables can be identified by (1) first reducing its dynamic to two modes exchanging energy between each other and (2) averaging the dynamic of the reduced system with respect to the phase of the fastest mode. Surprisingly, the global and slow dynamic of the system remains Hamiltonian (despite the parametric excitation) and the study of its associated effective potential shows how parametric excitation can turn the unstable open state into a stable one. Numerical experiments support the accuracy of the time-averaged reduced Hamiltonian in capturing the global and slow dynamic of the full system.

Intramolecular energy transfer and the driving mechanisms for large-amplitude collective motions of clusters.

This paper uncovers novel and specific dynamical mechanisms that initiate large-amplitude collective motions in polyatomic molecules. These mechanisms are understood in terms of intramolecular energy transfer between modes and driving forces. Structural transition dynamics of a six-atom cluster between a symmetric and an elongated isomer is highlighted as an illustrative example of what is a general message. First, we introduce a general method of hyperspherical mode analysis to analyze the energy transfer among internal modes of polyatomic molecules. In this method, the (3n-6) internal modes of an n-atom molecule are classified generally into three coarse level gyration-radius modes, three fine level twisting modes, and (3n-12) fine level shearing modes. We show that a large amount of kinetic energy flows into the gyration-radius modes when the cluster undergoes structural transitions by changing its mass distribution. Based on this fact, we construct a reactive mode as a linear combination of the three gyration-radius modes. It is shown that before the reactive mode acquires a large amount of kinetic energy, activation or inactivation of the twisting modes, depending on the geometry of the isomer, plays crucial roles for the onset of a structural transition. Specifically, in a symmetric isomer with a spherical mass distribution, activation of specific twisting modes drives the structural transition into an elongated isomer by inducing a strong internal centrifugal force, which has the effect of elongating the mass distribution of the system. On the other hand, in an elongated isomer, inactivation of specific twisting modes initiates the structural transition into a symmetric isomer with lower potential energy by suppressing the elongation effect of the internal centrifugal force and making the effects of the potential force dominant. This driving mechanism for reactions as well as the present method of hyperspherical mode analysis should be widely applicable to molecular reactions in which a system changes its overall mass distribution in a significant way.

Elastic origin of chiral selection in DNA wrapping.

We investigated the mechanism that underlies the chiral selection on the direction of wrapping of DNA around a nucleosome core particle. A coarse-grained model for the formation of a nucleosome is introduced, in which DNA is treated as a semiflexible polymer and the histone core is modeled by a spherical particle. Asymmetric coupling between bending and twisting is incorporated into the model DNA, which is expected from the double-stranded helical structure of DNA. We show that the tendency of DNA to twist in a left-handed manner upon bending gives rise to the selective left-handed wrapping, provided that the size of the core particle is chosen appropriately. This result suggests the critical importance of the chiral asymmetry inherent in the bending-twisting elasticity of DNA as well as the size of the core in determining the handedness of wrapping in nucleosome formation.

Gyration-radius dynamics in structural transitions of atomic clusters.

This paper is concerned with the structural transition dynamics of the six-atom Morse cluster with zero total angular momentum, which serves as an illustrative example of the general reaction dynamics of isolated polyatomic molecules. It develops a methodology that highlights the interplay between the effects of the potential energy topography and those of the intrinsic geometry of the molecular internal space. The method focuses on the dynamics of three coarse variables, the molecular gyration radii. By using the framework of geometric mechanics and hyperspherical coordinates, the internal motions of a molecule are described in terms of these three gyration radii and hyperangular modes. The gyration radii serve as slow collective variables, while the remaining hyperangular modes serve as rapidly oscillating "bath" modes. Internal equations of motion reveal that the gyration radii are subject to two different kinds of forces: One is the ordinary force that originates from the potential energy function of the system, while the other is an internal centrifugal force. The latter originates from the dynamical coupling of the gyration radii with the hyperangular modes. The effects of these two forces often counteract each other: The potential force generally works to keep the internal mass distribution of the system compact and symmetric, while the internal centrifugal force works to inflate and elongate it. Averaged fields of these two forces are calculated numerically along a reaction path for the structural transition of the molecule in the three-dimensional space of gyration radii. By integrating the sum of these two force fields along the reaction path, an effective energy curve is deduced, which quantifies the gross work necessary for the system to change its mass distribution along the reaction path. This effective energy curve elucidates the energy-dependent switching of the structural preference between symmetric and asymmetric conformations. The present methodology should be of wide use for the systematic reduction of dimensionality as well as for the identification of kinematic barriers associated with the rearrangement of mass distribution in a variety of molecular reaction dynamics in vacuum.

Kinematic effects associated with molecular frames in structural isomerization dynamics of clusters.

Kinematic effects associated with movements of molecular frames, which specify instantaneous orientation of molecules, is investigated in structural isomerization dynamics of a triatomic cluster whose total angular momentum is zero. The principal-axis frame is employed to introduce the so-called principal-axis hyperspherical coordinates, with which the mechanism of structural isomerization dynamics of the cluster is systematically analyzed. A force called "democratic centrifugal force" is extracted from the associated kinematics. This force arises from an intrinsic non-Euclidean metric in the internal space and has an effect of distorting the triatomic cluster to a collapsed shape and of trapping the system around collinear transition states. The latter effect is particularly important in that the kinematics effectively makes a basin at the saddle (transition state) on the potential surface. Based on this framework, we study the effect of the gauge field associated with the Eckart frame in internal space, which has not been carefully examined in the conventional reaction rate theories. Numerical comparison between the dynamics with and without the gauge field has revealed that this field has an effect to suppress the rate of isomerization reaction to a considerable amount. Thus a theory neglecting this effect will significantly overestimate the rate of isomerization. We show the physical origin of this suppressing effect.