Atomistic model of DNA: phonons and base-pair opening

Sub-THz acoustic excitation of protein motion

The application of terahertz radiation has been shown to affect both protein structure and cellular function. As the key to such structural changes lies in the dynamic response of a protein, we focus on the susceptibility of the protein's internal dynamics to mechanical stress induced by acoustic pressure waves. We use the open-boundary molecular dynamics method, which allows us to simulate the protein exposed to acoustic waves. By analyzing the dynamic fluctuations of the protein subunits, we demonstrate that the protein is highly susceptible to acoustic waves with frequencies corresponding to those of the internal protein vibrations. This is confirmed by changes in the compactness of the protein structure. As the amplitude of the pressure wave increases, even larger deviations from average positions and larger changes in protein compactness are observed. Furthermore, performing the mode-projection analysis, we show that the breathing-like character of collective modes is enhanced at frequencies corresponding to those used to excite the protein.

Collective Domain Motion Facilitates Water Transport in SGLT1

The human sodium-glucose cotransporter protein (SGLT1) is an important representative of the sodium solute symporters belonging to the secondary active transporters that are critical to the homeostasis of sugar, sodium, and water in the cell. The underlying transport mechanism of SGLT1 is based on switching between inward- and outward-facing conformations, known as the alternating access model, which is crucial for substrate transport, and has also been postulated for water permeation. However, the nature of water transport remains unclear and is disputed along the passive and active transport, with the latter postulating the presence of the pumping effect. To better examine the water transport in SGLT1, we performed a series of equilibrium all-atom molecular dynamics simulations, totaling over 6 μs of sample representative conformational states of SGLT1 and its complexes, with the natural substrates, ions, and inhibitors. In addition to elucidating the basic physical factors influencing water permeation, such as channel openings and energetics, we focus on dynamic flexibility and its relationship with domain motion. Our results clearly demonstrate a dependence of instantaneous water flux on the channel opening and local water diffusion in the channel, strongly supporting the existence of a passive water transport in SGLT1. In addition, a strong correlation found between the local water diffusion and protein domain motion, resembling the "rocking-bundle" motion, reveals its facilitating role in the water transport.

Low-frequency vibrational modes in G-quadruplexes reveal the mechanical properties of nucleic acids

Low-frequency vibrations play an essential role in biomolecular processes involving DNA such as gene expression, charge transfer, drug intercalation, and DNA–protein recognition. However, understanding the vibrational basis of these mechanisms relies on theoretical models due to the lack of experimental evidence. Here we present the low-frequency vibrational spectra of G-quadruplexes (structures formed by four strands of DNA) and B-DNA characterized using femtosecond optical Kerr-effect spectroscopy. Contrary to expectation, we found that G-quadruplexes show several strongly underdamped delocalized phonon-like modes that have the potential to contribute to the biology of the DNA at the atomic level. In addition, G-quadruplexes present modes at a higher frequency than B-DNA demonstrating that changes in the stiffness of the molecule alter its gigahertz to terahertz vibrational profile.

Low-frequency vibrations play an essential role in biomolecular processes involving DNA such as gene expression, charge transfer, drug intercalation, and DNA–protein recognition.

In silico prediction of the in vitro behavior of polymeric gene delivery vectors

Non-viral gene delivery vectors have increasingly come under the spotlight, but their performaces are still far from being satisfactory. Therefore, there is an urgent need for forecasting tools and screening methods to enable the development of ever more effective transfectants. Here, coarse-grained (CG) models of gold standard transfectant poly(ethylene imine)s (PEIs) have been profitably used to investigate and highlight the effect of experimentally-relevant parameters, namely molecular weight (2 vs. 10 kDa) and topologies (linear vs. branched), protonation state, and ammine-to-phosphate ratios (N/Ps), on the complexation and the gene silencing efficiency of siRNA molecules. The results from the in vitro screening of cationic polymers and conditions were used to validate the in silico platform that we developed, such that the hits which came out of the CG models were of high practical relevance. We show that our in silico platform enables to foresee the most suitable conditions for the complexation of relevant siRNA-polycation assemblies, thereby providing a reliable predictive tool to test bench transfectants in silico, and foster the design and development of gene delivery vectors.

Kinetic theory for DNA melting with vibrational entropy

By treating DNA as a vibrating nonlinear lattice, an activated kinetic theory for DNA melting is developed to capture the breakage of the hydrogen bonds and subsequent softening of torsional and bending vibration modes. With a coarse-grained lattice model, we identify a key bending mode with GHz frequency that replaces the hydrogen vibration modes as the dominant out-of-phase phonon vibration at the transition state. By associating its bending modulus to a universal in-phase bending vibration modulus at equilibrium, we can hence estimate the entropic change in the out-of-phase vibration from near-equilibrium all-atom simulations. This and estimates of torsional and bending entropy changes lead to the first predictive and sequence-dependent theory with good quantitative agreement with experimental data for the activation energy of melting of short DNA molecules without intermediate hairpin structures.

Observation of coherent delocalized phonon-like modes in DNA under physiological conditions

Underdamped terahertz-frequency delocalized phonon-like modes have long been suggested to play a role in the biological function of DNA. Such phonon modes involve the collective motion of many atoms and are prerequisite to understanding the molecular nature of macroscopic conformational changes and related biochemical phenomena. Initial predictions were based on simple theoretical models of DNA. However, such models do not take into account strong interactions with the surrounding water, which is likely to cause phonon modes to be heavily damped and localized. Here we apply state-of-the-art femtosecond optical Kerr effect spectroscopy, which is currently the only technique capable of taking low-frequency (GHz to THz) vibrational spectra in solution. We are able to demonstrate that phonon modes involving the hydrogen bond network between the strands exist in DNA at physiologically relevant conditions. In addition, the dynamics of the solvating water molecules is slowed down by about a factor of 20 compared with the bulk.

First-principles investigation of the lattice vibrational properties of inorganic double helical XY (X = Li, Na, K, Rb, Cs; Y = P, As, Sb)

The vibrational frequencies of the newly discovered inorganic double helical compounds XY (X = Li, Na, K, Rb, Cs; Y = P, As, Sb) are sensitive to either cation or anion or both of them.

Terahertz Dynamics in Human Cells and Their Chromatin

The terahertz dynamics of human cells of the U937 line and their chromatin has been investigated by high-resolution inelastic X-ray scattering. To highlight its dynamical features in situ, nuclear DNA has been stained by uranyl-acetate salt. The general behavior of the collective dynamics of the whole cell is quite similar to that of bulk water, with a nearly wavevector-independent branch located at about 5 meV and a propagating mode with a linear trend corresponding to a speed of sound of 2900 ± 100 m/s. We provide the first experimental evidence for the existence of two branches also in the dispersion curves of chromatin. The high-energy mode displays an acoustic-like behavior with a sound velocity similar to unstained cells, but in this case the branch likely originates from the superposition of intramolecular DNA optic modes. A low-energy optic-like branch, distinctive of the chromatin moiety, is found at about 2.5 meV.

Competition between Na⁺ and Rb⁺ in the minor groove of DNA

The competition between Na⁺ and Rb⁺ ions in the minor groove of a synthetic B-DNA dodecamer d (CGCGAATTCGCG) is studied using molecular dynamics simulations as the ratio of these two ions changing from 9:1 to 1:9 with the DNA merged into the solvent of water molecule at 298 K. When the concentration of Rb⁺ ions increases, from minority to majority, Na⁺ ions are gradually released from the A tract, and the binding sites in the minor groove are occupied by Rb⁺ ions, extending from the A tract to the whole minor groove. Comparing Na⁺ with Rb⁺ ions, the conformation of the minor groove is influenced strongly by Na⁺ ions.

Direct determination of the base-pair force constant of DNA from the acoustic phonon dispersion of the double helix

Quantifying the molecular elasticity of DNA is fundamental to our understanding of its biological functions. Recently different groups, through experiments on tailored DNA samples and numerical models, have reported a range of stretching force constants (0.3 to 3 N/m). However, the most direct, microscopic measurement of DNA stiffness is obtained from the dispersion of its vibrations. A new neutron scattering spectrometer and aligned, wet spun samples have enabled such measurements, which provide the first data of collective excitations of DNA and yield a force constant of 83 N/m. Structural and dynamic order persists unchanged to within 15 K of the melting point of the sample, precluding the formation of bubbles. These findings are supported by large scale phonon and molecular dynamics calculations, which reconcile hard and soft force constants.

Vibrational softening of a protein on ligand binding

Neutron scattering experiments have demonstrated that binding of the cancer drug methotrexate softens the low-frequency vibrations of its target protein, dihydrofolate reductase (DHFR). Here, this softening is fully reproduced using atomic detail normal-mode analysis. Decomposition of the vibrational density of states demonstrates that the largest contributions arise from structural elements of DHFR critical to stability and function. Mode-projection analysis reveals an increase of the breathing-like character of the affected vibrational modes consistent with the experimentally observed increased adiabatic compressibility of the protein on complexation.

Acoustic spectroscopy of DNA in the gigahertz range

We find a parametric resonance in the gigahertz range of DNA dynamics, generated by pumping hypersound. The resonance may be accompanied by the formation of localized phonon modes due to the random structure of elastic modulii of DNA.

Prediction of RNA multiloop and pseudoknot conformations from a lattice-based, coarse-grain tertiary structure model

We present a semiquantitative lattice model of RNA folding, which is able to reproduce complex folded structures such as multiloops and pseudoknots without relying on the frequently employed ad hoc generalization of the Jacobson-Stockmayer loop entropy. We derive the model parameters from the Turner description of simple secondary structural elements and pay particular attention to the unification of mismatch and coaxial stacking parameters as well as of border and nonlocal loop parameters, resulting in a reduced, unified parameter set for simple loops of arbitrary type and size. For elementary structures, the predictive power of the model is comparable to the standard secondary structure approaches, from which its parameters are derived. For complex structures, our approach offers a systematic treatment of generic effects of chain connectivity as well as of excluded volume or attractive interactions between and within all elements of the secondary structure. We reproduce the native structures of tRNA multiloops and of viral frameshift signal pseudoknots.

Overstretching of a 30 bp DNA duplex studied with steered molecular dynamics simulation: effects of structural defects on structure and force-extension relation

Single-molecule experiments on polymeric DNA show that the molecule can be overstretched at nearly constant force by about 70% beyond its relaxed contour length. In this publication we use steered molecular dynamics (MD) simulation to study the effect of structural defects on force-extension curves and structures at high elongation in a 30 base pair duplex pulled by its torsionally unconstrained 5' -5' ends. The defect-free duplex shows a plateau in the force-extension curve at 120 pN in which large segments with inclined and paired bases ("S-DNA") near both ends of the duplex coexist with a central B-type segment separated from the former by small denaturation bubbles. In the presence of a base mismatch or a nick, force-extension curves are very similar to the ones of the defect-free duplex. For the duplex with a base mismatch, S-type segments with highly inclined base pairs are not observed; rather, the overstretched duplex consists of B-type segments separated by denaturation bubbles. The nicked duplex evolves, via a two-step transition, into a two-domain structure characterized by a large S-type segment coexisting with several short S-type segments which are separated by short denaturation bubbles. Our results suggest that in the presence of nicks the force-extension curve of highly elongated duplex DNA might reflect locally highly inhomogeneous stretching.