Nonadiabatic forces in ion-solid interactions:The initial stages of radiation damage

The Born-Oppenheimer approximation is the keystone for molecular dynamics simulations of radiation damage processes; however, actual materials response involves nonadiabatic energy exchange between nuclei and electrons. In this work, time dependent density functional theory is used to calculate the electronic excitations produced by energetic protons in Al. We study the influence of these electronic excitations on the interatomic forces and find that they differ substantially from the adiabatic case, revealing a nontrivial connection between electronic and nuclear stopping that is absent in the adiabatic case. These results unveil new effects in the early stages of radiation damage cascades.

Conformational Dynamics of the Human Propeller Telomeric DNA Quadruplex on a Microsecond Time Scale

The human telomeric DNA sequence with four repeats can fold into a parallel-stranded propeller-type topology. NMR structures solved under molecular crowding experiments correlate with the crystal structures found with crystal-packing interactions that are effectively equivalent to molecular crowding. This topology has been used for rationalization of ligand design and occurs experimentally in a number of complexes with a diversity of ligands, at least in the crystalline state. While G-quartet stems have been well characterised, the interactions of the TTA loop with the G-quartets are much less defined. To better understand the conformational variability and structural dynamics of the propeller-type topology, we performed molecular dynamics simulations in explicit solvent up to 1.5 µs. The analysis provides a detailed atomistic account of the dynamic nature of the TTA loops highlighting their interactions with the G-quartets including formation of an A:A base pair, triad, pentad and hexad. The results present a threshold in quadruplex simulations, with regards to understanding the flexible nature of the sugar-phosphate backbone in formation of unusual architecture within the topology. Furthermore, this study stresses the importance of simulation time in sampling conformational space for this topology.

Comparison of dynamics of wildtype and V94M human UDP-galactose 4-epimerase-A computational perspective on severe epimerase-deficiency galactosemia

UDP-galactose 4'-epimerase (GALE) catalyzes the interconversion of UDP-galactose and UDP-glucose, an important step in galactose catabolism. Type III galactosemia, an inherited metabolic disease, is associated with mutations in human GALE. The V94M mutation has been associated with a very severe form of type III galactosemia. While a variety of structural and biochemical studies have been reported that elucidate differences between the wildtype and this mutant form of human GALE, little is known about the dynamics of the protein and how mutations influence structure and function. We performed molecular dynamics simulations on the wildtype and V94M enzyme in different states of substrate and cofactor binding. In the mutant, the average distance between the substrate and both a key catalytic residue (Tyr157) and the enzyme-bound NAD(+) cofactor and the active site dynamics are altered making substrate binding slightly less stable. However, overall stability or dynamics of the protein is not altered. This is consistent with experimental findings that the impact is largely on the turnover number (kcat), with less substantial effects on Km. Active site fluctuations were found to be correlated in enzyme with substrate bound to just one of the subunits in the homodimer suggesting inter-subunit communication. Greater active site loop mobility in human GALE compared to the equivalent loop in Escherichia coli GALE explains why the former can catalyze the interconversion of UDP-N-acetylgalactosamine and UDP-N-acetylglucosamine while the bacterial enzyme cannot. This work illuminates molecular mechanisms of disease and may inform the design of small molecule therapies for type III galactosemia.

Ultrafast charge dynamics in an amino acid induced by attosecond pulses

In the past few years, attosecond techniques have been implemented for the investigation of ultrafast dynamics in molecules. The generation of isolated attosecond pulses characterized by a relatively high photon flux has opened up new possibilities in the study of molecular dynamics. In this paper, we report on experimental and theoretical results of ultrafast charge dynamics in a biochemically relevant molecule, namely, the amino acid phenylalanine. The data represent the first experimental demonstration of the generation and observation of a charge migration process in a complexmolecule, where electron dynamics precede nuclear motion. The application of attosecond technology to the investigation of electron dynamics in biologically relevant molecules represents a multidisciplinary work, which can open new research frontiers: those in which few-femtosecond and even subfemtosecond electron processes determine the fate of biomolecules. It can also open new perspectives for the development of new technologies, for example, in molecular electronics, where electron processes on an ultrafast temporal scale are essential to trigger and control the electron current on the scale of the molecule.

Nonconservative current-driven dynamics: beyond the nanoscale

Long metallic nanowires combine crucial factors for nonconservative current-driven atomic motion. These systems have degenerate vibrational frequencies, clustered about a Kohn anomaly in the dispersion relation, that can couple under current to form nonequilibrium modes of motion growing exponentially in time. Such motion is made possible by nonconservative current-induced forces on atoms, and we refer to it generically as the waterwheel effect. Here the connection between the waterwheel effect and the stimulated directional emission of phonons propagating along the electron flow is discussed in an intuitive manner. Nonadiabatic molecular dynamics show that waterwheel modes self-regulate by reducing the current and by populating modes in nearby frequency, leading to a dynamical steady state in which nonconservative forces are counter-balanced by the electronic friction. The waterwheel effect can be described by an appropriate effective nonequilibrium dynamical response matrix. We show that the current-induced parts of this matrix in metallic systems are long-ranged, especially at low bias. This nonlocality is essential for the characterisation of nonconservative atomic dynamics under current beyond the nanoscale.