Molecular Dynamics Simulation of a PNA·DNA·PNA Triple Helix in Aqueous Solution

Molecular dynamics simulations have been used to explore the conformational flexibility of a PNA·DNA·PNA triple helix in aqueous solution. Three 1.05 ns trajectories starting from different but reasonable conformations have been generated and analyzed in detail. All three trajectories converge within about 300 ps to produce stable and very similar conformational ensembles, which resemble the crystal structure conformation in many details. However, in contrast to the crystal structure, there is a tendency for the direct hydrogen-bonds observed between the amide hydrogens of the Hoogsteen-binding PNA strand and the phosphate oxygens of the DNA strand to be replaced by water-mediated hydrogen bonds, which also involve pyrimidine O2 atoms. This structural transition does not appear to weaken the triplex structure but alters groove widths and so may relate to the potential for recognition of such structures by other ligands (small molecules or proteins). Energetic analysis leads us to conclude that the reason that the hybrid PNA/DNA triplex has quite different helical characteristics from the all-DNA triplex is not because the additional flexibility imparted by the replacement of sugar−phosphate by PNA backbones allows motions to improve base-stacking but rather that base-stacking interactions are very similar in both types of triplex and the driving force comes from weak but definate conformational preferences of the PNA strands.

Interplay Between Finite Resources and a Local Defect in an Asymmetric Simple Exclusion Process

When particle flux is regulated by multiple factors such as particle supply and varying transport rate, it is important to identify the respective dominant regimes. We extend the well-studied totally asymmetric simple exclusion model to investigate the interplay between a controlled entrance and a local defect site. The model mimics cellular transport phenomena where there is typically a finite particle pool and nonuniform moving rates due to biochemical kinetics. Our simulations reveal regions where, despite an increasing particle supply, the current remains constant while particles redistribute in the system. Exploiting a domain wall approach with mean-field approximation, we provide a theoretical ground for our findings. The results in steady-state current and density profiles provide quantitative insights into the regulation of the transcription and translation process in bacterial protein synthesis.