Slow nucleic acid unzipping kinetics from sequence-defined barriers

Recent experiments on unzipping of RNA helix-loop structures by force have shown that ≈40-base molecules can undergo kinetic transitions between two well-defined “open” and “closed” states, on a timescale ≈1 sec [Liphardt et al., Science 297, 733-737 (2001)]. Using a simple dynamical model, we show that these phenomena result from the slow kinetics of crossing large free energy barriers which separate the open and closed conformations. The dependence of barriers on sequence along the helix, and on the size of the loop(s) is analyzed. Some DNA and RNA sequences that could show dynamics on different time scales, or three(or more)-state unzipping, are proposed. Our dynamical model is also applied to the unzipping of long (kilo-basepair) DNA molecules at constant force. Received 29 July 2002 / Received in final form 5 February 2003 Published online: 16 April 2003 RID="a" ID="a"e-mail: cocco@ldfc.u-strasbg.fr RID="b" ID="b"e-mail: jmarko@uic.edu RID="c" ID="c"e-mail: monasson@lpt.ens.fr     

Structure-property relationships in major ampullate spider silk as deduced from polarized FTIR spectroscopy

Polarized Fourier Transform Infrared (FTIR) spectroscopy is employed to study structure-property relationships in major ampullate spider silk being exposed to an external mechanical strain. From the measured infrared dichroism of aminoacid-residue - specific bands the molecular order parameter, the frequency width at half-maximum (FWHM) and the spectral position of the absorption maximum are determined in dependence on the external strain. For the highly ordered alanine-rich β sheets a change in the vibrational potential is found for macroscopic strains as low as a few percent. It can be quantitatively described by a quantum-mechanical approach in which the mechanical strain is treated as a weak external perturbation. The immediate microscopic response to the external field proves that β -sheeted crystals are tightly interconnected by pre-stretched chains as suggested recently (Y. Liu et al., Nat. Mater. 4, 901 (2005)).     

Probing DNA and RNA single molecules with a double optical tweezer

A double-tweezer setup is used to induce mechanical stress in systems of molecular biology. A double strand of DNA is first stretched and the data is compared to precedent experiments to check the experimental setup. Then a short foldable fragment of RNA is probed; the typical unfolding/refolding hysteresis behaviour of this kind of construction is shown and followed by a study of its elasticity and a comparison to a worm-like chain model. Eventually, we describe the unfolding of a larger RNA structure, which unfolds by multiple steps. We show that this unfolding is not reversible and that it presents numerous unfolding pathways.     

Localization properties of electronic states in a polaron model of poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA polymers

We numerically investigate localization properties of electronic states in a static model of poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA polymers with realistic parameters obtained by quantum-chemical calculation. The randomness in the on-site energies caused by the electron-phonon coupling is completely correlated to the off-diagonal parts. In the single electron model, the effect of the hydrogen-bond stretchings, the twist angles between the base pairs and the finite system size effects on the energy dependence of the localization length and on the Lyapunov exponent are given. The localization length is reduced by the influence of the fluctuations in the hydrogen bond stretchings. It is also shown that the helical twist angle affects the localization length in the poly(dG)-poly(dC) DNA polymer more strongly than in the poly(dA)-poly(dT) one. Furthermore, we show resonance structures in the energy dependence of the localization length when the system size is relatively small.     

Quantum diffusion in polaron model of poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA polymers

We numerically investigate quantum diffusion of an electron in a model of poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA polymers with fluctuation of the parameters due to the impact of colored noise. The randomness is introduced by fluctuations of distance between two consecutive bases along the stacked base pairs. We demonstrate that in the model the decay time of the correlation can control the spread of the electronic wavepacket. Furthermore it is shown that in a motional narrowing regime the averaging over fluctuation causes ballistic propagation of the wavepacket, and in the adiabatic regime the electronic states are affected by localization.     

A DFT study of electron or hole localization in a peptide containing asparagin

The mechanisms of protein degradation induced by ionisation are of great interest for radiobiology, improvement of mass spectroscopy and industrial processes such as radio sterilisation. Sequences containing asparagin are very sensitive especially if surrounded by glycine. Very few techniques allow a satisfying understanding of the processes induced by creation of an anionic or cationic site in a peptide. We used the methods of quantum chemistry (DFT/B3LYP with 6-31G^* basis set) to characterise the geometry modifications induced in the cations or in the anions derived from peptide Gly Asn Gly. The cationic sites are localised mostly close to the first peptidic bond and induce a lengthening of the Ca-C(O) bond. Conversely the anionic sites are localised on a carbonyl function. Implications are discussed considering the radiolytic products and the proposed mechanisms. Received 21 January 2002 Published online 13 September 2002     

Counterion vibrations in the DNA low-frequency spectra

The vibrations of univalent metal cations with respect to phosphate groups of the DNA backbone are described using the four-mass model approach (S.N. Volkov, S.N. Kosevich, J. Biomol. Struct. Dyn. 8, 1069 (1991)) extended in this paper. The force constant of the counterion-phosphate interaction is determined by considering the DNA with counterions as a lattice of ion crystal. For such ion-phosphate lattice the Madelung constant and the dielectric constant are estimated. The obtained value of the Madelung constant is lower than for the NaCl crystal, and its value is about 1.3. The dielectric constant is within 2.3-2.7 depending on the counterion type and form of the double helix. The calculations of the low-frequency spectra show that for the DNA with metal cations Na⁺ , K⁺ , Rb⁺ and Cs⁺ the frequency of ion-phosphate vibrations decreases from 174 to 96cm^-1 as the counterion mass increases. The obtained frequencies agree well with the vibrational spectra of polynucleotides in a dry state which prove our suggestion about the existence of the ion-phosphate lattice around the DNA double helix. The amplitudes of conformational vibrations for DNA in B -form are calculated as well. The results demonstrate that light counterions ( Na⁺ do not disturb the internal dynamics of the DNA. However, heavy counterions ( Cs⁺ have effect on the internal vibrations of the DNA structural elements.     

Effects of NH3, O2, and N2 co-implantation on Cu out-diffusion and antimicrobial properties of copper plasma-implanted polyethylene

Metal antibacterial reagents are effective in the enhancement of the antimicrobial properties of medical polymers. However, incorporation of metal antibacterial reagents into polymers using conventional methods usually results in unstable antimicrobial effects. Our previous research demonstrates that plasma immersion ion implantation (PIII) can be used to effectively incorporate metal antibacterial reagents such as Cu into polyethylene (PE) in the near surface region up to several hundred nanometers without causing noticeable damage to the polymer matrix. In this work, various gases including NH₃, O₂, and N₂ were plasma-implanted in concert with Cu plasma immersion ion implantation to study the effects of these gas species on the release rate of Cu from the substrate. Our experimental results reveal that the copper depth profiles are not affected significantly by NH₃, O₂, or N₂ co-implantation and these gas elements have similar depth profiles as Cu. Chemical analyses demonstrate that polar functional CO, CO, CN, CN, and CN bonds formed in the substrate play an important role in regulating Cu out-diffusion. Among the three gas species, N₂ shows the best effects in regulating Cu out-diffusion and produces the best long term antibacterial properties. The Cu retention and out-diffusion mechanism in the ion-implanted polyethylene is described.     

Sphere packing, helices and the polytope {3, 3, 5}

The packing of tetrahedra in face contact is well-known to be relevant to atomic clustering in many complex alloys. We briefly review some of the structures that can arise in this way, and introduce methods of dealing with the geometry of the polytope {3, 3, 5}, which is highly relevant to an understanding of these structures. Finally, we present a method of projection from S₃ to E₃ that enables coordinates for the key vertices of the collagen model of Sadoc and Rivier to be calculated. Received 27 March 2001     

Activated dynamics of semiflexible polymers on structured substrates

We study the thermally activated motion of semiflexible polymers in double-well potentials using field-theoretic methods. Shape, energy, and effective diffusion constant of kink excitations are calculated, and their dependence on the bending rigidity of the semiflexible polymer is determined. For symmetric potentials, the kink motion is purely diffusive whereas kink motion becomes directed in the presence of a driving force. We determine the average velocity of the semiflexible polymer based on the kink dynamics. The Kramers escape over the potential barriers proceeds by nucleation and diffusive motion of kink-antikink pairs, the relaxation to the straight configuration by annihilation of kink-antikink pairs. We consider both uniform and point-like driving forces. For the case of point-like forces the polymer crosses the potential barrier only if the force exceeds a critical value. Our results apply to the activated motion of biopolymers such as DNA and actin filaments or of synthetic polyelectrolytes on structured substrates.     

Ground state and glass transition of the RNA secondary structure

RNA molecules form a sequence-specific self-pairing pattern at low temperatures. We analyze this problem using a random pairing energy model as well as a random sequence model that includes a base stacking energy in favor of helix propagation. The free energy cost for separating a chain into two equal halves offers a quantitative measure of sequence specific pairing. In the low temperature glass phase, this quantity grows quadratically with the logarithm of the chain length, but it switches to a linear behavior of entropic origin in the high temperature molten phase. Transition between the two phases is continuous, with characteristics that resemble those of a disordered elastic manifold in two dimensions. For designed sequences, however, a power-law distribution of pairing energies on a coarse-grained level may be more appropriate. Extreme value statistics arguments then predict a power-law growth of the free energy cost to break a chain, in agreement with numerical simulations. Interestingly, the distribution of pairing distances in the ground state secondary structure follows a remarkable power-law with an exponent -4/3, independent of the specific assumptions for the base pairing energies.     

Electrostatic image effects for counterions between charged planar walls

We study the effect of dielectric inhomogeneities on the interaction between two planparallel charged surfaces with oppositely charged mobile charges in between. The dielectric constant between the surfaces is assumed to be different from the dielectric constant of the two semiinfinite regions bounded by the surfaces, giving rise to electrostatic image interactions. We show that on the weak-coupling level the image charge effects are generally small, making their mark only in the second-order fluctuation term. However, in the strong-coupling limit, the image effects are large and fundamental. They modify the interactions between the two surfaces in an essential way. Our calculations are particularly useful in the regime of parameters where computer simulations would be difficult and extremely time consuming due to the complicated nature of the long-range image potentials.     

Unfolding designable structures

Among an infinite number of possible folds, nature has chosen only about 1000 distinct folds to form protein structures. Theoretical studies suggest that selected folds are intrinsically more designable than others; these selected folds are unusually stable, a property called the designability principle. In this paper we use the 2D hydrophobic-polar lattice model to classify structures according to their designability, and Langevin dynamics to account for their time evolution in the presence of shear flow. We demonstrate that, among all possible folds, the more designable ones are easier to unfold due to their large number of surface-core bonds.     

Log-periodic oscillations for biased diffusion of a polymer chain in a porous medium

A coarse-grained off-lattice bead-spring model is used to reveal the complex dynamics of a polymer chain in a quenched porous medium in the presence of an external field B. The behavior of the mean square displacement (MSD) of the center chain bead and that of the center of mass of the chain as a function of time is studied at different values of the barrier concentration C, the field strength B and the chain length N. In a field, important information on the way in which chains move between obstacles and overcome them is gained from the MSD vs. time analysis in the directions parallel and perpendicular to the flow. Instead of a steady approach to uniform drift-like motion at low C, for sufficiently strong field B we observe logarithmic oscillations in the effective exponents describing the time dependence of the MSD along and perpendicular to field. A common nature of this phenomenon with oscillatory behavior, observed earlier for biased diffusion of tracers on random lattices, is suggested. Received 7 August 1998     

The photophysics of isolated protein chromophores

Gas-phase absorption properties of chromophores of several photoactive proteins have been studied experimentally at the electrostatic heavy-ion storage ring ELISA in Aarhus. The absorption wavelength has been calculated using an augmented effective Hamiltonian technique based on the multiconfigurational quasi-degenerate perturbation theory. The results have been compared to those of widely used state-specific second-order perturbation theory formalisms and their multistate extensions and also to ground-state linear response methods. It would appear that ab initio theory is now at a stage where the intrinsic properties of the chromophore molecules may be predicted with reasonable precision. There is evidence that in terms of absorption there is almost vacuum-like conditions in the hydrophobic interior of some proteins like the green fluorescent protein (GFP). In others, like for example the visual opsins, some significant perturbations are responsible for colour tuning.     

Is there a universality of the helix-coil transition in protein models?

The similarity in the thermodynamic properties of two completely different theoretical models for the helix-coil transition is examined critically. The first model is an all-atomic representation for a poly-alanine chain, while the second model is a minimal helix-forming model that contains no system specifics. Key characteristics of the helix-coil transition, in particular, the effective critical exponents of these two models agree with each other, within a finite-size scaling analysis. Received 8 December 1999     

Interactions of like-charged rods at low temperatures: Analytical theory vs. simulations

We investigate a system consisting of two like-charged infinitely long rods and neutralizing counterions at low temperatures, using both analytic theory and simulations. With some reasonable approximations we can analytically solve for several ground-state structures of the model, starting with states where all counterions are lined up in the gap between the rods, over planar configurations, where the counterions are divided up into a fraction which resides between the rods, and counterions which are located on the outer surfaces, up to configurations which cover the full rod surfaces. Using parallel tempering simulations, we are able to study the system over a wide range of temperatures. At low temperatures we find good agreement with our T = 0 results. At higher temperatures, the strong coupling (SC) theory delivers qualitatively better results. We furthermore demonstrate that for the SC theory and our ground-state approximations to yield quantitative agreement, three parameters are required to be large, the strong-coupling parameter Ξ, the Rouzina-Bloomfield parameter, and the ratio of the average distance of the counterions to the radius of the rods. In the case of the latter ratio being small, our T = 0 results show better agreement with the simulation data at very low temperatures.     

Direct evidence of a magnetization transfer between laser-polarized xenon and protons of a cage-molecule in water

In a dedicated experimental setup, we directly prepare liquid-state NMR samples containing laser-polarized xenon with nuclear polarization larger than 5% at pressures up to 4 bars. Coating of the NMR tube surface allows us to increase the self-relaxation time of xenon in the gaseous phase to approximately 4.5 hours. Using a modified SPINOE pulse sequence, we present the first direct detection of a regioselective proton signal enhancement of a molecule -cyclodextrin) dissolved in water resulting from cross-polarization between laser-polarized xenon and protons. Received 16 March 2000 and Received in final form 22 May 2000