Plasmon Interpretation of 25 cm−1 Mode in DNA

Abstract

An extension of the effective field approach for the normal mode dynamics of dissolved DNA polymers has been applied to study the vibrational modes of DNA-hydration sheath-counterion system, to include the effect of site bound counterions on the system dynamics. An alternative interpretation has been suggested for a 25 cm^?1 mode recently observed in DNA samples and interpreted earlier as an interhelical mode. Analysing the eigenvectors this mode is found to possess a large electric dipole moment with longitudinal collective oscillations of the system. These characteristics identify this mode as a collective plasmon mode. Possible physical reasons for the existence of this character have been presented.     

Long-range DNA-water interactions

DNA functions only in aqueous environments and adopts different conformations depending on the hydration level. The dynamics of hydration water and hydrated DNA leads to rotating and oscillating dipoles that, in turn, give rise to a strong megahertz to terahertz absorption. Investigating the impact of hydration on DNA dynamics and the spectral features of water molecules influenced by DNA, however, is extremely challenging because of the strong absorption of water in the megahertz to terahertz frequency range. In response, we have employed a high-precision megahertz to terahertz dielectric spectrometer, assisted by molecular dynamics simulations, to investigate the dynamics of water molecules within the hydration shells of DNA as well as the collective vibrational motions of hydrated DNA, which are vital to DNA conformation and functionality. Our results reveal that the dynamics of water molecules in a DNA solution is heterogeneous, exhibiting a hierarchy of four distinct relaxation times ranging from ∼8 ps to 1 ns, and the hydration structure of a DNA chain can extend to as far as ∼18 Å from its surface. The low-frequency collective vibrational modes of hydrated DNA have been identified and found to be sensitive to environmental conditions including temperature and hydration level. The results reveal critical information on hydrated DNA dynamics and DNA-water interfaces, which impact the biochemical functions and reactivity of DNA.     

Collective motion in DNA and its role in drug intercalation

The effects of collective motion in DNA as reflected by resonance coupling among its intact segments have been discussed for both linear and circular DNA molecules. The results indicate that due to the effects of this kind of internal collective motion, the energy will be at times highly concentrated at some spots. As a result of the overfocus of energy, the stress built up along the direction of hydrogen bonds between complementary base pairs will be dramatically increased, rupturing a series of consecutive hydrogen bonds simultaneously and resulting in a suddenly free jerk, such that the DNA molecule will undergo a local “quake.” The “hole” formed by this kind of quake-like motion will be large enough for bulky drugs to gain entrance and intercalate into DNA. Even for smaller drugs, this local quake-like motion can also provide a significant mode of entry for intercalation. Energy minimizations carried out for DNA–drug complexes indicate that, for most drugs, a distortion or disruption of 2 to 4 base pairs occurs at the intercalation site in DNA molecules.     

Filtering molecular dynamics trajectories to reveal low-frequency collective motions: phospholipase A2

A novel method for analysing molecular dynamics trajectories has been developed, which filters out high frequencies using digital signal processing techniques and facilitates focusing on the low-frequency collective motions of proteins. These motions involve low energy slow motions, which lead to important biological phenomena such as domain closure and allosteric effects in enzymes. The filtering method treats each of the atomic trajectories obtained from the molecular dynamics simulation as a "signal". The trajectories of each of the atoms in the system (or any subset of interest) are Fourier transformed to the frequency domain, a filtering function is applied and then an inverse transformation back to the time domain yields the filtered trajectory. The filtering method has been used to study the dynamics of the enzyme phospholipase A2. In the filtered trajectory, all the high frequency bond and valence angle vibrations were eliminated, leaving only low-frequency motion, mainly fluctuations in torsions and conformational transitions. Analysis of this trajectory revealed interesting motions of the protein, including concerted movements of helices, and changes in shape of the active site cavity. Unlike normal mode analysis, which has been used to study the motion of proteins, this method does not require converged minimizations or diagonalization of a matrix of second derivatives. In addition, anharmonicity, multiple minima and conformational transitions are treated explicitly. Thus, the filtering method avoids most of the approximations implicit in other investigations of the dynamic behaviour of large systems.     

Cloning of Environmental Genomic Fragments as Physical Markers for Monitoring Microbial Populations in Coking Wastewater Treatment System

The association between community functional shift and dynamics of genomic DNA composition can be used to identify functionally relevant populations as indicator organisms for systems monitoring. In this work, fingerprinting-based community DNA hybridization was used to monitor community structural dynamics and identify genomic fragments whose abundance shifts were concomitant to changes in COD removal capacity in a reactor. A laboratory-scale anaerobic–anoxic–oxic fixed biofilm system treating coking wastewater was operated with (LR mode) or without effluent recirculation (LNR mode). The contribution to total chemical oxygen demand (COD) removal by the anoxic reactor increased from 4% in LNR mode to 26% in LR mode. Long primer RAPD (randomly amplified polymorphic DNA) community fingerprints of the anoxic reactor also changed most significantly from the one similar to the anaerobic reactor to one similar to the oxic reactor. DNA hybridization revealed one signature band of 2.1 kb shared by the anoxic and oxic reactors in LR, but not LNR mode. Clone library profiling of this band resulted in one predominant 2.1-kb genomic fragment (B3) with no homologous sequences in GenBank. Real-time polymerase chain reaction indicated that copy numbers of B3 in the anoxic reactor under LR mode were 69 times higher than that under LNR mode, concomitant to a significant increase in COD removal capacity in this reactor. The different patterns of distribution of B3 in the laboratory system and a comparable malfunctioning industrial system demonstrated the potential of this genomic fragment as physical markers in systems monitoring. In addition, this genomic fragment may allow sequence-guided isolation of the host microbe.     

Coarse-grained normal mode analysis in structural biology

The realization that experimentally observed functional motions of proteins can be predicted by coarse-grained normal mode analysis has renewed interest in applications to structural biology. Notable applications include the prediction of biologically relevant motions of proteins and supramolecular structures driven by their structure-encoded collective dynamics; the refinement of low-resolution structures, including those determined by cryo-electron microscopy; and the identification of conserved dynamic patterns and mechanically key regions within protein families. Additionally, hybrid methods that couple atomic simulations with deformations derived from coarse-grained normal mode analysis are able to sample collective motions beyond the range of conventional molecular dynamics simulations. Such applications have provided great insight into the underlying principles linking protein structures to their dynamics and their dynamics to their functions.     

Selective excitation of native fluctuations during thermal unfolding simulations: horse heart cytochrome c as a case study

The effect of temperature on the activation of native fluctuation motions during molecular dynamics unfolding simulations of horse heart cytochrome c has been studied. Essential dynamics analysis has been used to analyze the preferred directions of motion along the unfolding trajectories obtained by high temperature simulations. The results of this study have evidenced a clear correlation between the directions of the deformation motions that occur in the first stage of the unfolding process and few specific essential motions characterizing the 300 K dynamics of the protein. In particular, one of those collective motions, involved in the fluctuation of a loop region, is specifically excited in the thermal denaturation process, becoming progressively dominant during the first 500 ps of the unfolding simulations. As further evidence, the essential dynamics sampling performed along this collective motion has shown a tendency of the protein to promptly unfold. According to these results, the mechanism of thermal induced denaturation process involves the selective excitation of one or few specific equilibrium collective motions.     

On the role of arginine-glutamic acid ion pair in the ATP hydrolysis

The complex between adenosine triphosphate (ATP) and 4-guanidinobutyric acid (GBA) has been studied by infrared spectroscopy dry and hydrated (60% relative humidity). Partial nonenzymic hydrolysis has been detected, as deduced from characteristic bands of adenosine diphosphate (ADP) and inorganic orthophosphate formation. An infrared continuum, which increases upon hydration, demonstrates that the hydrogen bonded system in this complex has a large proton polarizability due to collective proton fluctuation. On this basis, a mechanism for splitting of lytic water molecules is also discussed.     

Dimerization affects collective dynamics of triosephosphate isomerase

Molecular dynamics simulations (30-60 ns runs) are performed on free/apo triosephosphate isomerase (TIM) to determine any correlation between collective motions and loop 6 dynamics. Native TIM is reported to be active only as a homodimer even though cooperativity has not been observed between the two identical subunits. Both dimeric and monomeric (isolated from dimer) forms of TIM are simulated in explicit water at 300 K and 1 bar to inspect any differences between the structures in terms of fluctuation dynamics and functionally important loop 6 dynamics/closure. Significant cross-correlations between residue fluctuations are observed in the dimer, which result from the global counter-rotations of the two identical subunits in the essential modes of the dimer. Specifically, the first essential mode contributing to 34% of overall motion of the dimer is strongly coupled to the loop 6's closure over the active site. In contrast, such significant correlations cannot be observed in the monomeric structure, which maintains relatively localized motions of the loops in the essential modes. Thus, the onset of collective motions at ns time scale due to dimerization has functional implications as to the coordination of loop 6 closure.     

Dynamics of charge carrier transfer in a DNA molecule

An equivalent electric circuit has been developed which describes the charge transfer in DNA molecule. A computer simulation of the charge carrier transfer dynamics in the molecule has been performed based on this circuit. It was found that the switching time of a molecular junction lies in the femtosecond range and depends on the frequency of the input electric signal. An increase in the frequency of the input signal in the range from 1 GHz to 4 THz and a reduction of temperature lead to a decrease in the current passing through the DNA molecule. It has been shown that the sequence of the DNA base pairs defines the rate of localization and delocalization of holes and controls the signal propagation rate in the DNA molecule.     

Dynamics of double stranded DNA reptation from bacteriophage.

The dynamics of dsDNA release process from a phage head has been analyzed theoretically. The process was considered as dsDNA reptation through the phage tail. The driving force is assumed to be the decrease of the DNA globule free energy on its releasing from the head in the surrounding medium. The results of the equilibrium theory on an intraphage DNA globule were applied. Three possible sources of friction were examined. The first one is in the inner channel of the tail. The second is the friction of DNA segments in the whole globule volume, which is essential when the globule decondensation is a collective process, at simultaneous moving of all the turns (mechanism 1). The third is the globule friction with the capsid inner surface, that is most important when decondensation proceeds via the globule rotation as a whole spool (mechanism 2). Mechanism 1 would require a lot of time for ejection. Mechanism 2 would lead to different ejection dynamics of short- and long-tailed phages. Comparison of the theoretical results with the published experimental data argues in favor of mechanism 2.     

Inherent Dynamics of the Acid-Sensing Ion Channel 1 Correlates with the Gating Mechanism

The acid-sensing ion channel 1 (ASIC1) is a key receptor for extracellular protons. Although numerous structural and functional studies have been performed on this channel, the structural dynamics underlying the gating mechanism remains unknown. We used normal mode analysis, mutagenesis, and electrophysiological methods to explore the relationship between the inherent dynamics of ASIC1 and its gating mechanism. Here we show that a series of collective motions among the domains and subdomains of ASIC1 correlate with its acid-sensing function. The normal mode analysis result reveals that the intrinsic rotation of the extracellular domain and the collective motions between the thumb and finger induced by proton binding drive the receptor to experience a deformation from the extracellular domain to the transmembrane domain, triggering the channel pore to undergo “twist-to-open” motions. The movements in the transmembrane domain indicate that the likely position of the channel gate is around Leu440. These motion modes are compatible with a wide body of our complementary mutations and electrophysiological data. This study provides the dynamic fundamentals of ASIC1 gating.     

Domain motions in bacteriophage T4 lysozyme: A comparison between molecular dynamics and crystallographic data

A comparison of a series of extended molecular dynamics (MD) simulations of bacteriophage T4 lysozyme in solvent with X-ray data is presented. Essential dynamics analyses were used to derive collective fluctuations from both the simulated trajectories and a distribution of crystallographic conformations. In both cases the main collective fluctuations describe domain motions. The protein consists of an N- and C-terminal domain connected by a long helix. The analysis of the distribution of crystallographic conformations reveals that the N-terminal helix rotates together with either of these two domains. The main domain fluctuation describes a closure mode of the two domains in which the N-terminal helix rotates concertedly with the C-terminal domain, while the domain fluctuation with second largest amplitude corresponds to a twisting mode of the two domains, with the N-terminal helix rotating concertedly with the N-terminal domain. For the closure mode, the difference in hinge-bending angle between the most open and most closed X-ray structure along this mode is 49 degrees. In the MD simulation that shows the largest fluctuation along this mode, a rotation of 45 degrees was observed. Although the twisting mode has much less freedom than the closure mode in the distribution of crystallographic conformations, experimental results suggest that it might be functionally important. Interestingly, the twisting mode is sampled more extensively in all MD simulations than it is in the distribution of X-ray conformations. Proteins 31:116–127, 1998. © 1998 Wiley-Liss, Inc.     

Laser-Driven Hybridization of a Surface Plasmon Resonance Collective Mode in a Monolayer of Silver Nanoparticles

Plasmonics - Laser-driven hybridization of a collective surface plasmon mode of a monolayer of silver nanoparticles has been studied as a function of irradiation power density. Two collective...     

Dynamical modeling of collective behavior from pigeon flight data: flock cohesion and dispersion

Several models of flocking have been promoted based on simulations with qualitatively naturalistic behavior. In this paper we provide the first direct application of computational modeling methods to infer flocking behavior from experimental field data. We show that this approach is able to infer general rules for interaction, or lack of interaction, among members of a flock or, more generally, any community. Using experimental field measurements of homing pigeons in flight we demonstrate the existence of a basic distance dependent attraction/repulsion relationship and show that this rule is sufficient to explain collective behavior observed in nature. Positional data of individuals over time are used as input data to a computational algorithm capable of building complex nonlinear functions that can represent the system behavior. Topological nearest neighbor interactions are considered to characterize the components within this model. The efficacy of this method is demonstrated with simulated noisy data generated from the classical (two dimensional) Vicsek model. When applied to experimental data from homing pigeon flights we show that the more complex three dimensional models are capable of simulating trajectories, as well as exhibiting realistic collective dynamics. The simulations of the reconstructed models are used to extract properties of the collective behavior in pigeons, and how it is affected by changing the initial conditions of the system. Our results demonstrate that this approach may be applied to construct models capable of simulating trajectories and collective dynamics using experimental field measurements of herd movement. From these models, the behavior of the individual agents (animals) may be inferred.