Dynamics of tRNA at different levels of hydration

Inelastic neutron scattering and spectroscopy methods to characterize dynamics in colloidal and soft matter systems

Colloidal systems and soft materials are well suited to neutron scattering, and the community has readily adopted elastic scattering techniques to investigate their structure. Due to their unique properties, neutrons may also be used to characterize the dynamics of soft materials over a wide range of length and time scales in situ. Both static structures and an understanding of how molecules move about their equilibrium positions is essential if we are to deliver on the promise of rationally designing soft materials. In this review we introduce the basics of neutron spectroscopy and explore the ways in which inelastic neutron scattering can be used to study colloidal and soft materials. Illustrative examples are chosen that highlight the phenomena suitable for investigation using this suite of techniques.

Influence of ion and hydration atmospheres on RNA structure and dynamics: insights from advanced theoretical and computational methods

RNA, a highly charged biopolymer composed of negatively charged phosphate groups, defies electrostatic repulsion to adopt well-defined, compact structures. Hence, the presence of positively charged metal ions is crucial not only for RNA's charge neutralization, but they also coherently decorate the ion atmosphere of RNA to stabilize its compact fold. This feature article elucidates various modes of close RNA-ion interactions, with a special emphasis on Mg2+ as an outer-sphere and inner-sphere ion. Through examples, we highlight how inner-sphere chelated Mg2+ stabilizes RNA pseudoknots, while outer-sphere ions can also exert long-range electrostatic interactions, inducing groove narrowing, coaxial helical stacking, and RNA ring formation. In addition to investigating the RNA's ion environment, we note that the RNA's hydration environment is relatively underexplored. Our study delves into its profound interplay with the structural dynamics of RNA, employing state-of-the-art atomistic simulation techniques. Through examples, we illustrate how specific ions and water molecules are associated with RNA functions, leveraging atomistic simulations to identify preferential ion binding and hydration sites. However, understanding their impact(s) on the RNA structure remains challenging due to the involvement of large length and long time scales associated with RNA's dynamic nature. Nevertheless, our contributions and recent advances in coarse-grained simulation techniques offer insights into large-scale structural changes dynamically linked to the RNA ion atmosphere. In this connection, we also review how different cutting-edge computational simulation methods provide a microscopic lens into the influence of ions and hydration on RNA structure and dynamics, elucidating distinct ion atmospheric components and specific hydration layers and their individual and collective impacts.

Probing the structure of human tRNA3Lys in the presence of ligands using docking, MD simulations and MSM analysis

The tRNA3Lys, which acts as a primer for human immunodeficiency virus type 1 (HIV-1) reverse transcription, undergoes structural changes required for the formation of a primer-template complex. Small molecules have been targeted against tRNA3Lys to inhibit the primer-template complex formation. The present study aims to understand the kinetics of the conformational landscape spanned by tRNA3Lys in apo form using molecular dynamics simulations and Markov state modeling. The study is taken further to investigate the effect of small molecules like 1,4T and 1,5T on structural conformations and kinetics of tRNA3Lys, and comparative analysis is presented. Markov state modeling of tRNA3Lys apo resulted in three metastable states where the conformations have shown the non-canonical structures of the anticodon loop. Based on analyses of ligand-tRNA3Lys interactions, crucial ion and water mediated H-bonds and free energy calculations, it was observed that the 1,4-triazole more strongly binds to the tRNA3Lys compared to 1,5-triazole. However, the MSM analysis suggest that the 1,5-triazole binding to tRNA3Lys has brought rigidity not only in the binding pocket (TΨC arm, D-TΨC loop) but also in the whole structure of tRNA3Lys. This may affect the easy opening of primer tRNA3Lys required for HIV-1 reverse transcription.

Hierarchical Hydration Dynamics of RNA with Nano-Water-Pool at Its Core

Many functional RNAs fold into a compact, roughly globular shape by minimizing the electrostatic repulsion between their negatively charged phosphodiester backbone. The fold of such close, compact RNA architecture is often so designed that its outer surface and complex core both are predominately populated by phosphate groups loosely sequestering bases in the intermediate layers. A number of helical junctions maintain the RNA core and its nano-water-pool. While the folding of RNA is manifested by its counterion environment composed of mixed mono- and divalent salts, the concerted role of ion and water in maintaining an RNA fold is yet to be explored. In this work, detailed atomistic simulations of SAM-I and Add Adenine riboswitch aptamers, and subgenomic flavivirus RNA (sfRNA) have been performed in a physiological mixed mono- and divalent salt environment. All three RNA systems have compact folds with a core diameter of range 1-1.7 nm. The spatiotemporal heterogeneity of RNA hydration was probed in a layer-wise manner by distinguishing the core, the intermediate, and the outer layers. The layer-wise decomposition of hydrogen bonds and collective single-particle reorientational dynamics reveal a nonmonotonic relaxation pattern with the slowest relaxation observed at the intermediate layers that involves functionally important tertiary motifs. The slowness of this intermediate layer is attributed to two types of long-resident water molecules: (i) water from ion-hydration layers and (ii) structurally trapped water (distant from ions). The relaxation kinetics of the core and the surface water essentially exposed to the phosphate groups show well-separated time scales from the intermediate layers. In the slow intermediate layers, site-specific ions and water control the functional dynamics of important RNA motifs like kink-turn, observed in different structure-probing experiments. Most interestingly, we find that as the size of the RNA core increases (SAM1 core < sfRNAcore < Add adenine core), its hydration tends to show faster relaxation. The hierarchical hydration and the layer-wise base-phosphate composition uniquely portray the globular RNA to act like a soft vesicle with a quasi-dynamic nano-water-pool at its core.

Coherent X-ray Scattering Reveals Nanoscale Fluctuations in Hydrated Proteins

Hydrated proteins undergo a transition in the deeply supercooled regime, which is attributed to rapid changes in hydration water and protein structural dynamics. Here, we investigate the nanoscale stress-relaxation in hydrated lysozyme proteins stimulated and probed by X-ray Photon Correlation Spectroscopy (XPCS). This approach allows us to access the nanoscale dynamics in the deeply supercooled regime (T = 180 K), which is typically not accessible through equilibrium methods. The observed stimulated dynamic response is attributed to collective stress-relaxation as the system transitions from a jammed granular state to an elastically driven regime. The relaxation time constants exhibit Arrhenius temperature dependence upon cooling with a minimum in the Kohlrausch-Williams-Watts exponent at T = 227 K. The observed minimum is attributed to an increase in dynamical heterogeneity, which coincides with enhanced fluctuations observed in the two-time correlation functions and a maximum in the dynamic susceptibility quantified by the normalized variance χ_T. The amplification of fluctuations is consistent with previous studies of hydrated proteins, which indicate the key role of density and enthalpy fluctuations in hydration water. Our study provides new insights into X-ray stimulated stress-relaxation and the underlying mechanisms behind spatiotemporal fluctuations in biological granular materials.

Low-temperature librations and dynamical transition in proteins at differing hydration levels

Hydration of water affects the dynamics and in turn the activity of biomacromolecules. We investigated the dependence of the librational oscillations and the dynamical transition on the hydrating conditions of two globular proteins with different structure and size, namely β-lactoglobulin (βLG) and human serum albumin (HSA), by spin-label electron paramagnetic resonance (EPR) in the temperature range of 120–270 K. The proteins were spin-labeled with 5-maleimide spin-label on free cysteins and prepared in the lyophilized state, at low (h = 0.12) and full (h = 2) hydration levels in buffer. The angular amplitudes of librations are small and almost temperature independent for both lyophilized proteins. Therefore, in these samples, the librational dynamics is restricted and the dynamical transition is absent. In the small and compact beta-structured βLG, the angular librational amplitudes increase with temperature and hydrating conditions, whereas hydration-independent librational oscillations whose amplitudes rise with temperature are recorded in the large and flexible alpha-structured HSA. Both βLG and HSA at low and fully hydration levels undergo the dynamical transition at about 230 K. The overall results indicate that protein librational dynamics is activated at the low hydration level h = 0.12 and highlight biophysical properties that are common to other biosamples at cryogenic temperatures.

Influence of hydration on segmental chain librations and dynamical transition in lipid bilayers

Continuous wave electron paramagnetic resonance spectroscopy of chain-labeled phospholipids is used to investigate the effects of hydration on the librational oscillations and the dynamical transition of phospholipid membranes in the low-temperature range 120-270 K. Bilayers of dipalmitoylphostatidiycholine (DPPC) spin-labeled at the first acyl chain segments and at the methyl ends and prepared at full, low, and very low hydration are considered. The segmental mean-square angular amplitudes of librations, 〈α²〉, are larger in the bilayer interior than at the polar/apolar interface and larger in the fully and low hydrated than in the very low hydrated membranes. For chain segments at the beginning of the hydrocarbon region, 〈α²〉-values are markedly restricted and temperature independent in DPPC with the lowest water content, whereas they increase with temperature in the low and fully hydrated bilayers, particularly at the highest temperatures. For chain segments at the chain termini, the librational amplitudes increase progressively, first slowly and then more rapidly with temperature in bilayers at any level of hydration. From the temperature dependence of the mean-square librational amplitude, the dynamical transition is detected around 240 K at the polar/apolar interface in fully and low hydrated DPPC and at around 225 K at the inner hydrocarbon region for bilayers at any hydration condition. At the dynamical transition the bilayers cross low energy barriers of activation energy in the range 10-20 kJ/mol. The results highlight biophysical properties of DPPC bilayers at low-temperature and provide evidence of the effects of the hydration on the dynamical transition in bilayers.

Probing Small-Angle Molecular Motions with EPR Spectroscopy: Dynamical Transition and Molecular Packing in Disordered Solids

Disordered molecular solids present a rather broad class of substances of different origin—amorphous polymers, materials for photonics and optoelectronics, amorphous pharmaceutics, simple molecular glass formers, and others. Frozen biological media in many respects also may be referred to this class. Theoretical description of dynamics and structure of disordered solids still does not exist, and only some phenomenological models can be developed to explain results of particular experiments. Among different experimental approaches, electron paramagnetic resonance (EPR) applied to spin probes and labels also can deliver useful information. EPR allows probing small-angle orientational molecular motions (molecular librations), which intrinsically are inherent to all molecular solids. EPR is employed in its conventional continuous wave (CW) and pulsed—electron spin echo (ESE)—versions. CW EPR spectra are sensitive to dynamical librations of molecules while ESE probes stochastic molecular librations. In this review, different manifestations of small-angle motions in EPR of spin probes and labels are discussed. It is shown that CW-EPR-detected dynamical librations provide information on dynamical transition in these media, similar to that explored with neutron scattering, and ESE-detected stochastic librations allow elucidating some features of nanoscale molecular packing. The possible EPR applications are analyzed for gel-phase lipid bilayers, for biological membranes interacting with proteins, peptides and cryoprotectants, for supercooled ionic liquids (ILs) and supercooled deep eutectic solvents (DESs), for globular proteins and intrinsically disordered proteins (IDPs), and for some other molecular solids.

Universal dynamical onset in water at distinct material interfaces

Interfacial water remains liquid and mobile much below 0 °C, imparting flexibility to the encapsulated materials to ensure their diverse functions at subzero temperatures. However, a united picture that can describe the dynamical differences of interfacial water on different materials and its role in imparting system-specific flexibility to distinct materials is lacking. By combining neutron spectroscopy and isotope labeling, we explored the dynamics of water and the underlying substrates independently below 0 °C across a broad range of materials. Surprisingly, while the function-related anharmonic dynamical onset in the materials exhibits diverse activation temperatures, the surface water presents a universal onset at a common temperature. Further analysis of the neutron experiment and simulation results revealed that the universal onset of water results from an intrinsic surface-independent relaxation: switching of hydrogen bonds between neighboring water molecules with a common energy barrier of ∼35 kJ mol⁻¹.

We demonstrated that the dynamical onset of interfacial water is an intrinsic property of water itself, resulting from a surface independent relaxation process in water with an approximately universal energy barrier of ∼35 kJ mol⁻¹.

Functional-State Dependence of Picosecond Protein Dynamics

We examine temperature-dependent picosecond dynamics of two benchmarking proteins lysozyme and cytochrome c using temperature-dependent terahertz permittivity measurements. We find that a double Arrhenius temperature dependence with activation energies E₁ ∼ 0.1 kJ/mol and E₂ ∼ 10 kJ/mol fits the folded and ligand-free state response. The higher activation energy is consistent with the so-called protein dynamical transition associated with beta relaxations at the solvent-protein interface. The lower activation energy is consistent with correlated structural motions. When the structure is removed by denaturing, the lower-activation-energy process is no longer present. Additionally, the lower-activation-energy process is diminished with ligand binding but not for changes in the internal oxidation state. We suggest that the lower-energy activation process is associated with collective structural motions that are no longer accessible with denaturing or binding.

Structural alphabets for conformational analysis of nucleic acids available at dnatco.datmos.org

A detailed description of the dnatco.datmos.org web server implementing the universal structural alphabet of nucleic acids is presented. It is capable of processing any mmCIF- or PDB-formatted files containing DNA or RNA molecules; these can either be uploaded by the user or supplied as the wwPDB or PDB-REDO structural database access code. The web server performs an assignment of the nucleic acid conformations and presents the results for the intuitive annotation, validation, modeling and refinement of nucleic acids.

Temperature Induced Dynamical Transition of Biomolecules in Polarizable and Nonpolarizable TIP3P Water

Temperature induced dynamical transition (DT), associated with a sharp rise in molecular flexibility, is well-known to be exhibited between 270 and 280 K in glycerol to 200-230 K in hydrated biomolecules and is controlled by diffusivity (viscosity) of the solvation layer. In the molecular dynamics (MD) community, especially for water as a solvent, this has been an intense area of research despite decades of investigations. However, in general, water in these studies is described by empirical nonpolarizable force fields in which electronic polarizability is treated implicitly with effective charges and related parameters. This might have led to the present trait of discovery that DTs of biomolecules, irrespective of the potential functions for water models used, occur within a narrow band of temperature variation (30-40 K). Whereas a water molecule in a biomolecular surface and one in bulk are polarized differently, therefore explicit treatment of water polarizability would be a powerful approach toward the treatment of hydration water, believed to cause the DT manifestation. Using MD simulations, we investigated the effects of polarizable water on the DT of biomolecules and the dynamic properties of hydration water. We chose two types of solutes: globular protein (lysozyme) and more open and flexible RNAs (a hairpin and a riboswitch) with different natures of hydrophilic sites than proteins in general. We found that the characteristic temperature of DT ( T_DT) for the solutes in polarizable water is always higher than that in its nonpolarizable counterpart. In particular, for RNAs, the variations are found to be ∼45 K between the two water models, whereas for the more compact lysozyme, it is only ∼4 K. The results are discussed in light of the enormous increase in relaxation times of a liquid upon cooling in the paradigm of dynamic switchover in hydration water with liquid-liquid phase transition, derived from the existence of the second critical point. Our result supports the idea that structures of biomolecules and their interactions with the hydration water determines T_DT and provides evidence for the decisive role of polarizable water on the onset of DT, which has been hitherto ignored.

Dynamical transition in molecular glasses and proteins observed by spin relaxation of nitroxide spin probes and labels

In glassy substances and biological media, dynamical transitions are observed in neutron scattering that manifests itself as deviations of the translational mean-squared displacement, 〈x²〉, of hydrogen atoms from harmonic dynamics. In biological media, the deviation occurs at two temperature intervals, at ∼100-150 K and at ∼170-230 K, and it is attributed to the motion of methyl groups in the former case and to the transition from harmonic to anharmonic or diffusive motions in the latter case. In this work, electron spin echo (ESE) spectroscopy-a pulsed version of electron paramagnetic resonance-is applied to study the spin relaxation of nitroxide spin probes and labels introduced in molecular glass former o-terphenyl and in protein lysozyme. The anisotropic contribution to the rate of the two-pulse ESE decay, ΔW, is induced by spin relaxation appearing because of restricted orientational stochastic molecular motion; it is proportional to 〈α²〉τ_c, where 〈α²〉 is the mean-squared angle of reorientation of the nitroxide molecule around the equilibrium position and τ_c is the correlation time of reorientation. The ESE time window allows us to study motions with τ_c < 10^-7 s. For glassy o-terphenyl, the 〈α²〉τ_c temperature dependence shows a transition near 240 K, which is in agreement with the literature data on 〈x²〉. For spin probes of essentially different size, the obtained data were found to be close, which evidences that motion is cooperative, involving a nanocluster of several neighboring molecules. For the dry lysozyme, the 〈α²〉τ_c values below 260 K were found to linearly depend on the temperature in the same way as it was observed in neutron scattering for 〈x²〉. As spin relaxation is influenced only by stochastic motion, the harmonic motions seen in ESE must be overdamped. In the hydrated lysozyme, ESE data show transitions near 130 K for all nitroxides, near 160 K for the probe located in the hydration layer, and near 180 K for the label in the protein interior. For this system, the two latter transitions are not observed in neutron scattering. The ESE-detected transitions are suggested to be related with water dynamics in the nearest hydration shell: with water glass transition near 130 K and with the onset of overall water molecular reorientations near 180 K; the disagreement with neutron scattering is ascribed to the larger time window for ESE-detected motions.

Determination of Dynamical Heterogeneity from Dynamic Neutron Scattering of Proteins

Motional displacements of hydrogen (H) in proteins can be measured using incoherent neutron-scattering methods. These displacements can also be calculated numerically using data from molecular dynamics simulations. An enormous amount of data on the average mean-square motional displacement (MSD) of H as a function of protein temperature, hydration, and other conditions has been collected. H resides in a wide spectrum of sites in a protein. Some H are tightly bound to molecular chains, and the H motion is dictated by that of the chain. Other H are quite independent. As a result, there is a distribution of motions and MSDs of H within a protein that is denoted dynamical heterogeneity. The goal of this paper is to incorporate a distribution of MSDs into models of the H incoherent intermediate scattering function, I(Q,t), that is calculated and observed. The aim is to contribute information on the distribution as well as on the average MSD from comparison of the models with simulations and experiment. For example, we find that simulations of I(Q,t) in lysozyme are well reproduced if the distribution of MSDs is bimodal with two broad peaks rather than a single broad peak.

Dynamic Neutron Scattering by Biological Systems

Dynamic neutron scattering directly probes motions in biological systems on femtosecond to microsecond timescales. When combined with molecular dynamics simulation and normal mode analysis, detailed descriptions of the forms and frequencies of motions can be derived. We examine vibrations in proteins, the temperature dependence of protein motions, and concepts describing the rich variety of motions detectable using neutrons in biological systems at physiological temperatures. New techniques for deriving information on collective motions using coherent scattering are also reviewed.

Low-Temperature Dynamics of Chain-Labeled Lipids in Ester- and Ether-Linked Phosphatidylcholine Membranes

Continuous wave electron paramagnetic resonance spectroscopy and two-pulse echo detected spectra of chain-labeled lipids are used to study the dynamics of frozen lipid membranes over the temperature range 77-260 K. Bilayers of ester-linked dihexadecanoylphosphatidylcholine (DPPC) with noninterdigitated chains and ether-linked dihexadecyl phosphatidylcholine (DHPC) with interdigitated chains are considered. Rapid stochastic librations of small angular amplitude are found in both lipid matrices. In noninterdigitated DPPC bilayers, the mean-square angular amplitude, [Formula: see text], of the motion increases with temperature and it is larger close to the chain termini than close to the polar/apolar interface. In contrast, in interdigitated DHPC lamellae, [Formula: see text] is small and temperature and label-position independent at low temperature and increases steeply at high temperature. The rotational correlation time, τ_c, of librations lies in the subnanosecond range for DPPC and in the nanosecond range for DHPC. In all membrane samples, the temperature dependence of [Formula: see text] resembles that of the mean-square atomic displacement revealed by neutron scattering and a dynamical transition is detected in the range 210-240 K. The results highlight the librational oscillations and the glass-like behavior in bilayer and interdigitated lipid membranes.

Microscopic diffusion in hydrated encysted eggs of brine shrimp

BACKGROUND

We have studied microscopic diffusion of water in fully hydrated encysted eggs of brine shrimp (Artemia).

METHODS

We have utilized quasielastic neutron scattering.

RESULTS

Dry eggs of brine shrimp were rehydrated using (1) water without additives, (2) eutectic mixture of water and dimethyl sulfoxide, and (3) a concentrated aqueous solution of lithium chloride. Despite the complexity of the hydrated multicellular organism, measurable microscopic diffusivity of water is rather well defined. Pure hydration water in eggs exhibits freezing temperature depression, whereas hydration water in eggs mixed with dimethyl sulfoxide or lithium chloride does not crystallize at all.

CONCLUSIONS

The characteristic size of the voids occupied by water or aqueous solvents in hydrated brine shrimp eggs is between 2 and 10nm. Those voids are accessible to co-solvents such as dimethyl sulfoxide and lithium chloride. There is no evidence of intracellular water in the hydrated eggs.

GENERAL SIGNIFICANCE

The lack of intracellular water in the fully hydrated (but still under arrested development) state must be linked to the unique resilience against adverse environmental factors documented not only for the anhydrous, but also hydrated encysted eggs of brine shrimp.

Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

BACKGROUND

Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics.

SCOPE OF REVIEW

This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties.

MAJOR CONCLUSIONS

In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration.

GENERAL SIGNIFICANCE

An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.

Comparative Structural Dynamics of tRNA(Phe) with Respect to Hinge Region Methylated Guanosine: A Computational Approach

Transfer RNAs (tRNAs) contain various uniquely modified nucleosides thought to be useful for maintaining the structural stability of tRNAs. However, their significance for upholding the tRNA structure has not been investigated in detail at the atomic level. In this study, molecular dynamic simulations have been performed to assess the effects of methylated nucleic acid bases, N (2)-methylguanosine (m(2)G) and N (2)-N (2)-dimethylguanosine (m 2 (2) G) at position 26, i.e., the hinge region of E. coli tRNA(Phe) on its structure and dynamics. The results revealed that tRNA(Phe) having unmodified guanosine in the hinge region (G26) shows structural rearrangement in the core of the molecule, resulting in lack of base stacking interactions, U-turn feature of the anticodon loop, and TΨC loop. We show that in the presence of the unmodified guanosine, the overall fold of tRNA(Phe) is essentially not the same as that of m(2)G26 and m 2 (2) G26 containing tRNA(Phe). This structural rearrangement arises due to intrinsic factors associated with the weak hydrogen-bonding patterns observed in the base triples of the tRNA(Phe) molecule. The m(2)G26 and m 2 (2) G26 containing tRNA(Phe) retain proper three-dimensional fold through tertiary interactions. Single-point energy and molecular electrostatics potential calculation studies confirmed the structural significance of tRNAs containing m(2)G26 and m 2 (2) G26 compared to tRNA with normal G26, showing that the mono-methylated (m(2)G26) and dimethylated (m 2 (2) G26) modifications are required to provide structural stability not only in the hinge region but also in the other parts of tRNA(Phe). Thus, the present study allows us to better understand the effects of modified nucleosides and ionic environment on tRNA folding.

Atomistic details of protein dynamics and the role of hydration water

BACKGROUND

The importance of protein dynamics for their biological activity is now well recognized. Different experimental and computational techniques have been employed to study protein dynamics, hierarchy of different processes and the coupling between protein and hydration water dynamics. Yet, understanding the atomistic details of protein dynamics and the role of hydration water remains rather limited.

SCOOP OF REVIEW

Based on overview of neutron scattering, molecular dynamic simulations, NMR and dielectric spectroscopy results we present a general picture of protein dynamics covering time scales from faster than ps to microseconds and the influence of hydration water on different relaxation processes.

MAJOR CONCLUSIONS

Internal protein dynamics spread over a wide time range from faster than picosecond to longer than microseconds. We suggest that the structural relaxation in hydrated proteins appears on the microsecond time scale, while faster processes present mostly motion of side groups and some domains. Hydration water plays a crucial role in protein dynamics on all time scales. It controls the coupled protein-hydration water relaxation on 10-100ps time scale. This process defines the friction for slower protein dynamics. Analysis suggests that changes in amount of hydration water affect not only general friction, but also influence significantly the protein's energy landscape.

GENERAL SIGNIFICANCE

The proposed atomistic picture of protein dynamics provides deeper understanding of various relaxation processes and their hierarchy, similarity and differences between various biological macromolecules, including proteins, DNA and RNA. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo".

Charge screening in RNA: an integral route for dynamical enhancements

Electrostatic interactions of RNA are at the center of determining the dynamical flexibility and structural stability. By analysing neutron scattering spectroscopy, we show that fast dynamics of hydrated tRNA on ps to ns timescales increases with stronger charge screening, while its structural stability either increases or remains largely unchanged. An unprecedented electrostatic threshold for the onset of additional flexibility is induced from the correlation between the charge-screening density of counterions and the promoted dynamical properties. The results demonstrate that the enhanced dynamical flexibility of tRNA originates from local conformational relaxation coupled with stabilized charge screening rather than governed by fluctuation of hydrated counterions. The present study casts light on the specificity of electrostatic interactions in the thermodynamic balance between the dynamical flexibility and structural stability of RNA.

Elasticity and Inverse Temperature Transition in Elastin

Elastin is a structural protein and biomaterial that provides elasticity and resilience to a range of tissues. This work provides insights into the elastic properties of elastin and its peculiar inverse temperature transition (ITT). These features are dependent on hydration of elastin and are driven by a similar mechanism of hydrophobic collapse to an entropically favorable state. Using neutron scattering, we quantify the changes in the geometry of molecular motions above and below the transition temperature, showing a reduction in the displacement of water-induced motions upon hydrophobic collapse at the ITT. We also measured the collective vibrations of elastin gels as a function of elongation, revealing no changes in the spectral features associated with local rigidity and secondary structure, in agreement with the entropic origin of elasticity.

MD Simulations of tRNA and Aminoacyl-tRNA Synthetases: Dynamics, Folding, Binding, and Allostery

While tRNA and aminoacyl-tRNA synthetases are classes of biomolecules that have been extensively studied for decades, the finer details of how they carry out their fundamental biological functions in protein synthesis remain a challenge. Recent molecular dynamics (MD) simulations are verifying experimental observations and providing new insight that cannot be addressed from experiments alone. Throughout the review, we briefly discuss important historical events to provide a context for how far the field has progressed over the past few decades. We then review the background of tRNA molecules, aminoacyl-tRNA synthetases, and current state of the art MD simulation techniques for those who may be unfamiliar with any of those fields. Recent MD simulations of tRNA dynamics and folding and of aminoacyl-tRNA synthetase dynamics and mechanistic characterizations are discussed. We highlight the recent successes and discuss how important questions can be addressed using current MD simulations techniques. We also outline several natural next steps for computational studies of AARS:tRNA complexes.

Protein dynamics: from rattling in a cage to structural relaxation

We present an overview of protein dynamics based mostly on results of neutron scattering, dielectric relaxation spectroscopy and molecular dynamics simulations. We identify several major classes of protein motions on the time scale from faster than picoseconds to several microseconds, and discuss the coupling of these processes to solvent dynamics. Our analysis suggests that the microsecond backbone relaxation process might be the main structural relaxation of the protein that defines its glass transition temperature, while faster processes present some localized secondary relaxations. Based on the overview, we formulate a general picture of protein dynamics and discuss the challenges in this field.

Slowdown of Interhelical Motions Induces a Glass Transition in RNA

RNA function depends crucially on the details of its dynamics. The simplest RNA dynamical unit is a two-way interhelical junction. Here, for such a unit--the transactivation response RNA element--we present evidence from molecular dynamics simulations, supported by nuclear magnetic resonance relaxation experiments, for a dynamical transition near 230 K. This glass transition arises from the freezing out of collective interhelical motional modes. The motions, resolved with site-specificity, are dynamically heterogeneous and exhibit non-Arrhenius relaxation. The microscopic origin of the glass transition is a low-dimensional, slow manifold consisting largely of the Euler angles describing interhelical reorientation. Principal component analysis over a range of temperatures covering the glass transition shows that the abrupt slowdown of motion finds its explanation in a localization transition that traps probability density into several disconnected conformational pools over the low-dimensional energy landscape. Upon temperature increase, the probability density pools then flood a larger basin, akin to a lakes-to-sea transition. Simulations on transactivation response RNA are also used to backcalculate inelastic neutron scattering data that match previous inelastic neutron scattering measurements on larger and more complex RNA structures and which, upon normalization, give temperature-dependent fluctuation profiles that overlap onto a glass transition curve that is quasi-universal over a range of systems and techniques.

Motional displacements in proteins: The origin of wave-vector-dependent values

The average mean-square displacement, 〈r(2)〉, of H atoms in a protein is frequently determined using incoherent neutron-scattering experiments. 〈r(2)〉 is obtained from the observed elastic incoherent dynamic structure factor, S(i)(Q,ω=0), assuming the form S(i)(Q,ω=0) =exp(-Q(2)〈r(2)〉/3). This is often referred to as the Gaussian approximation (GA) to S(i)(Q,ω=0). 〈r(2)〉 obtained in this way depends on the value of the wave vector, Q considered. Equivalently, the observed S(i)(Q,ω=0) deviates from the GA. We investigate the origin of the Q dependence of 〈r(2)〉 by evaluating the scattering functions in different approximations using molecular dynamics (MD) simulation of the protein lysozyme. We find that keeping only the Gaussian term in a cumulant expansion of S(Q,ω) is an accurate approximation and is not the origin of the Q dependence of 〈r(2)〉. This is demonstrated by showing that the term beyond the Gaussian is negligible and that the GA is valid for an individual atom in the protein. Rather, the Q dependence (deviation from the GA) arises from the dynamical heterogeneity of the H in the protein. Specifically it arises from representing, in the analysis of data, this diverse dynamics by a single average scattering center that has a single, average 〈r(2)〉. The observed Q dependence of 〈r(2)〉 can be used to provide information on the dynamical heterogeneity in proteins.

Dynamics and rigidity in an intrinsically disordered protein, β-casein

The emergence of intrinsically disordered proteins (IDPs) as a recognized structural class has forced the community to confront a new paradigm of structure, dynamics, and mechanical properties for proteins. We present novel data on the similarities and differences in the dynamics and nanomechanical properties of IDPs and other biomacromolecules on the picosecond time scale. An IDP, β-casein (CAS), has been studied in a calcium bound and unbound state using neutron and light scattering techniques. We show that CAS partially folds and stiffens upon calcium binding, but in the unfolded state, it is softer than folded proteins such as green fluorescence protein (GFP). We also see that some localized diffusive motions in CAS have a larger amplitude than in GFP at this time scale but are still smaller than those observed in tRNA. In spite of these differences, CAS dynamics are consistent with the classes of motions seen in folded protein on this time scale.

NMR Water Self-Diffusion and Relaxation Studies on Sodium Polyacrylate Solutions and Gels in Physiologic Ionic Solutions

Water self-diffusion coefficients and longitudinal relaxation rates in sodium polyacrylate solutions and gels were measured by NMR, as a function of polymer content and structure in a physiological concentration range of monovalent and divalent cations, Ca2+ and Na+. Several physical models describing the self-diffusion of the solvent were applied and compared. A free-volume model was found to be in good agreement with the experimental results over a wide range of polymer concentrations. The longitudinal relaxation rate exhibited linear dependence on polymer concentration below a critical concentration and showed non-linear behavior at higher concentrations. Both the water self-diffusion and relaxation were less influenced by the polymer in the gel state than in the uncrosslinked polymer solutions. The effect of Na+ on the mobility of water molecules was practically undetectable. By contrast, addition of Ca2+ strongly increased the longitudinal relaxation rate while its effect on the self-diffusion coefficient was much less pronounced.

Structural intermediates and folding events in the early assembly of the ribosomal small subunit

Using all-atom explicit solvent molecular dynamics (MD) simulations, we investigated the early structural intermediates of the 5' domain of the 16S rRNA in Escherichia coli upon the removal of the primary binding r-proteins S4, S17, and S20 and the secondary binding r-protein S16. Removal of each r-protein corresponded to the disappearance of subdomains with correlated dynamics. Correlation-based network analysis of the MD trajectories of the naked rRNA showed that the different subdomains are connected via multiple pathways with high betweenness. These pathways cross at the internal loop of helix 17 (h17) in the five-way junction (5WJ). The structure of the internal loop is disrupted by the binding of S17 and rescued by the addition of S16, suggesting an important function of the secondary binding protein in biasing the rRNA folding landscape toward the native basin. Using structure-based Gō simulations, we investigated the folding barriers of the lower four-way junction (4WJ) with h6, which is the primary binding site of S20 and the first to be transcribed. The folding of the 4WJ is consistent with the protection patterns observed in hydroxyl radical footprinting. Results from the all-atom simulations show that the fluctuations in the 5WJ are independent of the fluctuations in the 4WJ, suggesting that the subdomains fold independently and are stabilized by primary r-proteins.

Temperature Dependence of Logarithmic-like Relaxational Dynamics of Hydrated tRNA

The dynamics of RNA within the β-relaxation region of 10 ps to 1 ns is crucial to its biological function. Because of its simpler chemical building blocks and the lack of the side methyl groups, faster relaxational dynamics of RNA compared to proteins can be expected. However, the situation is actually opposite. In this work, the relaxational dynamics of tRNA is measured by quasielastic neutron scattering and analyzed using the mode coupling theory, originally developed for glass-forming liquids. Our results reveal that the dynamics of tRNA follows a log-decay within the β-relaxation region, which is an important trait demonstrated by the dynamics of proteins. The dynamics of hydrated tRNA and lysozyme compared in the time domain further demonstrate that the slower dynamics of tRNA relative to proteins originates from the difference in the folded states of tRNA and proteins, as well as the influence of their hydration water.

Dynamics of protein and its hydration water: neutron scattering studies on fully deuterated GFP

We present a detailed analysis of the picosecond-to-nanosecond motions of green fluorescent protein (GFP) and its hydration water using neutron scattering spectroscopy and hydrogen/deuterium contrast. The analysis reveals that hydration water suppresses protein motions at lower temperatures (<~ 200 K), and facilitates protein dynamics at high temperatures. Experimental data demonstrate that the hydration water is harmonic at temperatures <~ 180-190 K and is not affected by the proteins' methyl group rotations. The dynamics of the hydration water exhibits changes at ~ 180-190 K that we ascribe to the glass transition in the hydrated protein. Our results confirm significant differences in the dynamics of protein and its hydration water at high temperatures: on the picosecond-to-nanosecond timescale, the hydration water exhibits diffusive dynamics, while the protein motions are localized to <~3 Å. The diffusion of the GFP hydration water is similar to the behavior of hydration water previously observed for other proteins. Comparison with other globular proteins (e.g., lysozyme) reveals that on the timescale of 1 ns and at equivalent hydration level, GFP dynamics (mean-square displacements and quasielastic intensity) are of much smaller amplitude. Moreover, the suppression of the protein dynamics by the hydration water at low temperatures appears to be stronger in GFP than in other globular proteins. We ascribe this observation to the barrellike structure of GFP.

Physical origin of anharmonic dynamics in proteins: new insights from resolution-dependent neutron scattering on homomeric polypeptides

Neutron scattering reveals a complex dynamics in polypeptide chains, with two main onsets of anharmonicity whose physical origin and biological role are still debated. In this study the dynamics of strategically selected homomeric polypeptides is investigated with elastic neutron scattering using different energy resolutions and compared with that of a real protein. Our data spotlight the dependence of anharmonic transition temperatures and fluctuation amplitudes on energy resolution, which we quantitatively explain in terms of a two-site model for the protein-hydration water energy landscape. Experimental data strongly suggest that the protein dynamical transition is not a mere resolution effect but is due to a real physical effect. Activation barriers and free energy values obtained for the protein dynamical transition allow us to make a connection with the two-well interaction potential of supercooled-confined water proposed to explain a low-density→high-density liquid-liquid transition.

Dynamic behavior of oligomeric inorganic pyrophosphatase explored by quasielastic neutron scattering

The purpose of this investigation is to determine whether a large oligomeric protein, inorganic pyrophosphatase (IPPase) from Thermococcus thioreducens with quaternary structural complexity, would have distinguishable dynamic characteristics compared to those of the small simple monomeric model protein, lysozyme. In this study, the β-relaxational dynamics of the two proteins, IPPase and lysozyme, are compared in the 10 ps to 0.5 ns time interval using quasi-elastic neutron scattering (QENS). Both of the protein dynamics show a characteristic logarithmic-like decay in the intermediate scattering function (ISF) of the hydrogen atoms. Distinguishable dynamical behavior found between two proteins reveals local flexibility and conformational substates unique to oligomeric structures. Moreover, the temperature dependence of the mean square displacement (MSD) of the hydrogen atoms in protein molecules, which is a traditional way to determine the "softness" of the protein molecule, is measured and shows no difference for the two proteins.

Hydrostatic and osmotic pressure study of the RNA hydration

The tertiary structure of nucleic acids results from an equilibrium between electrostatic interactions of phosphates, stacking interactions of bases, hydrogen bonds between polar atoms and water molecules. Water interactions with ribonucleic acid play a key role in its structure formation, stabilization and dynamics. We used high hydrostatic pressure and osmotic pressure to analyze changes in RNA hydration. We analyzed the lead catalyzed hydrolysis of tRNA^Phe from S. cerevisiae as well as hydrolytic activity of leadzyme. Pb(II) induced hydrolysis of the single phosphodiester bond in tRNA^Phe is accompanied by release of 98 water molecules, while other molecule, leadzyme releases 86.

Controllable self-assembly of organic-inorganic amphiphiles containing Dawson polyoxometalate clusters

An organic-inorganic molecular hybrid containing the Dawson polyoxometalate, ((C(4)H(9))(4)N)(5)H[P(2)V(3)W(15)O(59)(OCH(2))(3)CNHCOC(15)H(31)], was synthesized and its surfactant-like amphiphilic properties, represented by the formation of bilayer vesicles, were studied in polar solvents. The vesicle size decreases with both decreasing hybrid concentration and with increasing polarity of the solvent, independently. The self-assembly behavior of this hybrid can be controlled by introducing different counterions into the acetonitrile solutions. The addition of ZnCl(2) and NaI can cause a gradual decrease and increase of vesicular sizes, respectively. Tetraalkylammonium bromide is found to disassemble the vesicle assemblies. Moreover, the original counterions of the hybrid can be replaced with protons, resulting in pH-dependent formation of vesicles in aqueous solutions. The hybrid surfactant can further form micro-needle structures in aqueous solutions upon addition of Ca(2+) ions.