Polymer translocation driven by longitudinal and transversal time-dependent end-pulling forces.

In this article, we simulate the translocation of a semiflexible homopolymer through an extended pore, driven by both a constant and a time-dependent end-pulled force, employing a model introduced in previous studies. The time dependence is simplistically modeled as a cosine function, and we distinguish between two scenarios for the driving--longitudinal force and transversal force-depending on the relative orientation of the force, parallel or perpendicular, respectively, with respect to the pore axis. Besides some key differences between the two drivings, the mean translocation times present a large minimum region as a function of the frequency of the force that is typical of the resonant activation effect. The presence of the minimum is independent on the elastic characteristics of the polymeric chains and reveals a linear relation between the optimum mean translocation time and the corresponding period of the driving. The mean translocation times show different scaling exponents with the polymer length for different flexibilities. Lastly, we derive an analytical expression of the mean translocation time for low driving frequency, which clearly agrees with the simulations.

Mesoscopic model for DNA G-quadruplex unfolding.

Genomes contain rare guanine-rich sequences capable of assembling into four-stranded helical structures, termed G-quadruplexes, with potential roles in gene regulation and chromosome stability. Their mechanical unfolding has only been reported to date by all-atom simulations, which cannot dissect the major physical interactions responsible for their cohesion. Here, we propose a mesoscopic model to describe both the mechanical and thermal stability of DNA G-quadruplexes, where each nucleotide of the structure, as well as each central cation located at the inner channel, is mapped onto a single bead. In this framework we are able to simulate loading rates similar to the experimental ones, which are not reachable in simulations with atomistic resolution. In this regard, we present single-molecule force-induced unfolding experiments by a high-resolution optical tweezers on a DNA telomeric sequence capable of adopting a G-quadruplex conformation. Fitting the parameters of the model to the experiments we find a correct prediction of the rupture-force kinetics and a good agreement with previous near equilibrium measurements. Since G-quadruplex unfolding dynamics is halfway in complexity between secondary nucleic acids and tertiary protein structures, our model entails a nanoscale paradigm for non-equilibrium processes in the cell.

Active translocation of a semiflexible polymer assisted by an ATP-based molecular motor.

In this work we study the assisted translocation of a polymer across a membrane nanopore, inside which a molecular motor exerts a force fuelled by the hydrolysis of ATP molecules. In our model the motor switches to its active state for a fixed amount of time, while it waits for an ATP molecule which triggers the motor, during an exponentially distributed time lapse. The polymer is modelled as a beads-springs chain with both excluded volume and bending contributions, and moves in a stochastic three dimensional environment modelled with a Langevin dynamics at a fixed temperature. The resulting dynamics shows a Michaelis-Menten translocation velocity that depends on the chain flexibility. The scaling behavior of the mean translocation time with the polymer length for different bending values is also investigated.

Thermal and inertial resonances in DNA unzipping.

Single-molecule experiments combined with alternate forces are able to provide useful information not present in standard constant-force and -velocity pulling protocols. Here, we study the effects of such forces in the DNA mechanical unzipping by using an extension of the Peyrard-Bishop-Dauxois model. By changing the damping regime in the dynamical equations, we obtained two resonant mechanisms in both the mean time and the mean force of unzipping. One is thermally assisted and it is characterized by a matching between the period of the external force and the mean unzipping time of the DNA chain, while the other depends on the inertial properties of the system. Both mechanisms are studied systematically under different opening protocols and different parameters of the system. The main results here presented contribute in characterizing and finding optimized conditions in DNA unzipping experiments.

Active polymer translocation in the three-dimensional domain.

In this work we study the translocation process of a polymer through a nanochannel where a time dependent force is acting. Two conceptually different types of driving are used: a deterministic sinusoidal one and a random telegraph noise force. The mean translocation time presents interesting resonant minima as a function of the frequency of the external driving. For the computed sizes, the translocation time scales with the polymer length according to a power law with the same exponent for almost all the frequencies of the two driving forces. The dependence of the translocation time with the polymer rigidity, which accounts for the persistence length of the molecule, shows a different low frequency dependence for the two drivings.

Mechanical unfolding of a simple model protein goes beyond the reach of one-dimensional descriptions.

We study the mechanical unfolding of a simple model protein. The Langevin dynamics results are analyzed using Markov-model methods which allow to describe completely the configurational space of the system. Using transition-path theory we also provide a quantitative description of the unfolding pathways followed by the system. Our study shows a complex dynamical scenario. In particular, we see that the usual one-dimensional picture: free-energy vs end-to-end distance representation, gives a misleading description of the process. Unfolding can occur following different pathways and configurations which seem to play a central role in one-dimensional pictures are not the intermediate states of the unfolding dynamics.

Modeling the interaction of DNA with alternating fields.

We study the influence of a terahertz field on thermal properties of DNA molecules. A Peyrard-Bishop-Dauxois model with the inclusion of a solvent interaction term is considered. The terahertz field is included as a sinusoidal driven force in the equation of motion. We show how under certain field and system parameters, the melting transition and bubble formation are modified.

Directional motion of forced polymer chains with hydrodynamic interaction.

We study the propulsion of a one-dimensional (1D) polymer chain under sinusoidal external forces in the overdamped (low Reynolds number) regime. We show that, when hydrodynamical interactions are included, the polymer presents directional motion which depends on the phase differences of the external force applied along the chain. Moreover, the velocity shows a maximum as a function of the frequency. We discuss the relevance of all these results in light of recent nanotechnology experiments.

Mesoscopic model for free-energy-landscape analysis of DNA sequences.

A mesoscopic model which allows us to identify and quantify the strength of binding sites in DNA sequences is proposed. The model is based on the Peyrard-Bishop-Dauxois model for the DNA chain coupled to a Brownian particle which explores the sequence interacting more importantly with open base pairs of the DNA chain. We apply the model to promoter sequences of different organisms. The free energy landscape obtained for these promoters shows a complex structure that is strongly connected to their biological behavior. The analysis method used is able to quantify free energy differences of sites within genome sequences.

Thermal and mechanical properties of a DNA model with solvation barrier.

We study the thermal and mechanical behaviors of DNA denaturation in the frame of the mesoscopic Peyrard-Bishop-Dauxois model with the inclusion of solvent interaction. By analyzing the melting transition of a homogeneous A-T sequence, we are able to set suitable values of the parameters of the model and study the formation and stability of bubbles in the system. Then, we focus on the case of the P5 promoter sequence and use the principal component analysis of the trajectories to extract the main information on the dynamical behavior of the system. We find that this analysis method gives an excellent agreement with previous biological results.

Model for hand-over-hand motion of molecular motors.

A simple flashing ratchet model in two dimensions is proposed to simulate the hand-over-hand motion of two head molecular motors such as kinesin. Extensive Langevin simulations of the model are performed. Good qualitative agreement with the expected behavior is observed. We discuss different regimes of motion and efficiency depending on model parameters.

Computational diagnosis of protein conformational diseases: short molecular dynamics simulations reveal a fast unfolding of r-LDL mutants that cause familial hypercholesterolemia.

The molecular basis of conformational diseases frequently resides in mutant proteins constituting a subset of the vast mutational space. While the subtleties of protein structure point to molecular dynamics (MD) techniques as promising tools for an efficient exploration of such a space, the average size of proteins and the time scale of unfolding events make this goal difficult with present computational capabilities. We show here, nevertheless, that an efficient approach is already feasible for modular proteins. Familial hypercholesterolemia (FH) is a conformational disease linked to mutations in the gene encoding the low density lipoprotein receptor. A high percentage of these mutations has been found in the seven small modular binding repeats of the receptor. Taking advantage of its small size, we have performed an in depth MD study of the fifth binding repeat. Fast unfolding dynamics have been observed in the absence of a structural bound calcium ion, which agrees with its reported essential role in the stability of the module. In addition, several mutations detected in FH patients have been analyzed, starting from the native conformation. Our results indicate that in contrast with the wild type protein and an innocuous control mutant, disease-related mutants experience, in short simulation times (2-8 ns), gross departures from the native state that lead to unfolded conformations and, in some cases, to binding site desorganization deriving in calcium release. Computational diagnosis of mutations leading to conformational diseases seems thus feasible, at least for small or modular pathogenic proteins.

Mode-locking of mobile discrete breathers.

We study numerically synchronization phenomena of mobile discrete breathers in dissipative nonlinear lattices periodically forced. When varying the driving intensity, the breather velocity generically locks at rational multiples of the driving frequency. In most cases, the locking plateau coincides with the linear stability domain of the resonant mobile breather and desynchronization occurs by the regular appearance of type-I intermittencies. However, some plateaus also show chaotic mobile breathers with locked velocity in the locking region. The addition of a small subharmonic driving tames the locked chaotic solution and enhances the stability of resonant mobile breathers.

Can we model DNA at the mesoscale?

Modelling DNA is useful for understanding its properties better but it is also challenging because many of these properties involve hundreds of base pairs or more, or time scales which are much longer than the time scales accessible to molecular dynamics. It is therefore necessary to develop models at a mesoscale, which include enough details to describe the properties of interest, for instance the biological sequence, while staying sufficiently simple and realistic.We discuss here two examples: a dynamical model to study the mechanical denaturation of DNA, which probes the sequence on various scales, and a model for the self assembly of DNA which describes the formation of hairpins and allows us to study its kinetics.

Dissipative discrete breathers: periodic, quasiperiodic, chaotic, and mobile.

The properties of discrete breathers in dissipative one-dimensional lattices of nonlinear oscillators subject to periodic driving forces are reviewed. We focus on oscillobreathers in the Frenkel-Kontorova chain and rotobreathers in a ladder of Josephson junctions. Both types of exponentially localized solutions are easily obtained numerically using adiabatic continuation from the anticontinuous limit. Linear stability (Floquet) analysis allows the characterization of different types of bifurcations experienced by periodic discrete breathers. Some of these bifurcations produce nonperiodic localized solutions, namely, quasiperiodic and chaotic discrete breathers, which are generally impossible as exact solutions in Hamiltonian systems. Within a certain range of parameters, propagating breathers occur as attractors of the dissipative dynamics. General features of these excitations are discussed and the Peierls-Nabarro barrier is addressed. Numerical scattering experiments with mobile breathers reveal the existence of two-breather bound states and allow a first glimpse at the intricate phenomenology of these special multibreather configurations.

Discrete breathers in dissipative lattices.

We study the properties of discrete breathers, also known as intrinsic localized modes, in the one-dimensional Frenkel-Kontorova lattice of oscillators subject to damping and external force. The system is studied in the whole range of values of the coupling parameter, from C=0 (uncoupled limit) up to values close to the continuum limit (forced and damped sine-Gordon model). As this parameter is varied, the existence of different bifurcations is investigated numerically. Using Floquet spectral analysis, we give a complete characterization of the most relevant bifurcations, and we find (spatial) symmetry-breaking bifurcations that are linked to breather mobility, just as it was found in Hamiltonian systems by other authors. In this way moving breathers are shown to exist even at remarkably high levels of discreteness. We study mobile breathers and characterize them in terms of the phonon radiation they emit, which explains successfully the way in which they interact. For instance, it is possible to form "bound states" of moving breathers, through the interaction of their phonon tails. Over all, both stationary and moving breathers are found to be generic localized states over large values of C, and they are shown to be robust against low temperature fluctuations.

Mirror symmetry breaking through an internal degree of freedom leading to directional motion.

We analyze here the minimal conditions for directional motion (net flow in phase space) of a molecular motor placed on a mirror-symmetric environment and driven by a center-symmetric and time-periodic force field. The complete characterization of the deterministic limit of the dissipative dynamics of several realizations of this minimal model reveals a complex structure in the phase diagram in parameter space, with intertwined regions of pinning (closed orbits) and directional motion. This demonstrates that the mirror symmetry breaking, which is needed for directional motion to occur, can operate through an internal degree of freedom coupled to the translational one.