Protein dynamical transition does not require protein structure

Anharmonic transitions in nearly dry L-cysteine I

Two special dynamical transitions of universal character have recently been observed in macromolecules (lysozyme, myoglobin, bacteriorhodopsin, DNA and RNA) at T* ~100-150 K and T(D) ~180-220 K. The underlying mechanisms governing these transitions have been the subject of debate. In the present work, a survey is reported on the temperature dependence of structural, vibrational and thermodynamical properties of a nearly anhydrous amino acid (orthorhombic polymorph of the amino acid l-cysteine at a hydration level of 3.5%). The temperature dependence of x-ray powder diffraction patterns, Raman spectra and specific heat revealed these two transitions at T* = 70 K and T(D) = 230 K for this sample. The data were analyzed considering amino acid-amino acid, amino acid-water, water-water phonon-phonon interactions and molecular rotor activation. Our results indicated that the two referred temperatures define the triggering of very simple and particular events that govern all the interactions of the biomolecular: activation of CH(2) rigid rotors (T < T* ), phonon-phonon interactions between specific amino acid and water dimer vibrational modes (T* < T < T(D)), and water rotational barriers surpassing (T > T(D)).

Protein dynamical transition at 110 K

Proteins are known to undergo a dynamical transition at around 200 K but the underlying mechanism, physical origin, and relationship to water are controversial. Here we report an observation of a protein dynamical transition as low as 110 K. This unexpected protein dynamical transition precisely correlated with the cryogenic phase transition of water from a high-density amorphous to a low-density amorphous state. The results suggest that the cryogenic protein dynamical transition might be directly related to the two liquid forms of water proposed at cryogenic temperatures.

Electrostatics of the protein-water interface and the dynamical transition in proteins

Atomic displacements of hydrated proteins are dominated by phonon vibrations at low temperatures and by dissipative large-amplitude motions at high temperatures. A crossover between the two regimes is known as a dynamical transition. Recent experiments indicate a connection between the dynamical transition and the dielectric response of the hydrated protein. We analyze two mechanisms of the coupling between the protein atomic motions and the protein-water interface. The first mechanism considers viscoelastic changes in the global shape of the protein plasticized by its coupling to the hydration shell. The second mechanism involves modulations of the local motions of partial charges inside the protein by electrostatic fluctuations. The model is used to analyze mean-square displacements of iron of metmyoglobin reported by Mössbauer spectroscopy. We show that high displacement of heme iron at physiological temperatures is dominated by electrostatic fluctuations. Two onsets, one arising from the viscoelastic response and the second from electrostatic fluctuations, are seen in the temperature dependence of the mean-square displacements when the corresponding relaxation times enter the instrumental resolution window.

Terahertz and far infrared spectroscopy of alanine-rich peptides having variable ellipticity

Terahertz spectra of four alanine-rich peptides with known secondary structures were studied by terahertz time domain spectroscopy (THz-TDS) and by Fourier transform infrared spectroscopy (FTIR) using a synchrotron light source and a liquid-helium cooled bolometer. At ambient temperatures the usable bandwidth was restricted to 0.2-1.5 THz by the absorbance of water. The existence of a solvation shell around the peptide in solution was observed and its size estimated to be between 11 and 17 Å. By cooling the peptide solution to 80 K in order to reduce the water absorbance the bandwidth was increased to 0.1-3.0 THz for both THz-TDS and FTIR. Spectra were consistent with monotonic absorbance of the peptide and the existence of a solid amorphous low density solvation shell.

The role of structured water in mediating general anesthetic action on alpha4beta2 nAChR

Water is an essential component for many biological processes. Pauling proposed that water might play a critical role in general anesthesia by forming water clathrates around anesthetic molecules. To examine potential involvement of water in general anesthesia, we analyzed water within alpha4beta2 nAChR, a putative protein target hypersensitive to volatile anesthetics. Experimental structure-derived closed- and open-channel nAChR systems in a fully hydrated lipid bilayer were examined using all-atom molecular dynamics simulations. At the majority of binding sites in alpha4beta2 nAChR, halothane replaced the slow-exchanging water molecules and caused a regional water population decrease. Only two binding sites had an increased quantity of water in the presence of halothane, where water arrangements resemble clathrate-like structures. The small number of such clathrate-like water clusters suggests that the formation of water clathrates is unlikely to be a primary cause for anesthesia. Despite the decrease in water population at most of the halothane binding sites, the number of sites that can be occupied transiently by water is actually increased in the presence of halothane. Many of these water sites were located between two subunits or in regions containing agonist binding sites or critical structural elements for transducing agonist binding to channel gating. Changes in water sites in the presence of halothane affected water-mediated protein-protein interactions and the protein dynamics, which can have direct impact on protein function. Collectively, water contributes to the action of anesthetics in proteins by mediating interactions between protein subunits and altering protein dynamics, instead of forming water clathrates around anesthetics.

Protein functional landscapes, dynamics, allostery: a tortuous path towards a universal theoretical framework

Energy landscape theories have provided a common ground for understanding the protein folding problem, which once seemed to be overwhelmingly complicated. At the same time, the native state was found to be an ensemble of interconverting states with frustration playing a more important role compared to the folding problem. The landscape of the folded protein – the native landscape – is glassier than the folding landscape; hence, a general description analogous to the folding theories is difficult to achieve. On the other hand, the native basin phase volume is much smaller, allowing a protein to fully sample its native energy landscape on the biological timescales. Current computational resources may also be used to perform this sampling for smaller proteins, to build a ‘topographical map’ of the native landscape that can be used for subsequent analysis. Several major approaches to representing this topographical map are highlighted in this review, including the construction of kinetic networks, hierarchical trees and free energy surfaces with subsequent structural and kinetic analyses. In this review, we extensively discuss the important question of choosing proper collective coordinates characterizing functional motions. In many cases, the substates on the native energy landscape, which represent different functional states, can be used to obtain variables that are well suited for building free energy surfaces and analyzing the protein's functional dynamics. Normal mode analysis can provide such variables in cases where functional motions are dictated by the molecule's architecture. Principal component analysis is a more expensive way of inferring the essential variables from the protein's motions, one that requires a long molecular dynamics simulation. Finally, the two popular models for the allosteric switching mechanism, ‘preexisting equilibrium’ and ‘induced fit’, are interpreted within the energy landscape paradigm as extreme points of a continuum of transition mechanisms. Some experimental evidence illustrating each of these two models, as well as intermediate mechanisms, is presented and discussed.

Effects of water on the structure and low/high temperature stability of confined proteins

In this study well-characterized model proteins were confined in silica nanoporous matrices. Confinement of the proteins in silica matrices allowed us to explore the role of water hydrogen bonding on the structures of the proteins in a broad range of temperatures (-120 degrees C to 95 degrees C). At low temperatures confinement suppressed freezing of water, which remained in the liquid state. We obtained direct evidence that the changes in the hydrogen bonding of water induced changes in the structure of confined proteins. At high temperatures, a reduction of hydrogen bonding of water facilitated protein-silica interactions and the confined proteins underwent denaturation. However, the incorporation of the osmolyte, trehalose, reduced protein-silica interactions, and altered the hydrogen bonding of water. As a result, the high temperature thermal stability of the confined proteins was greatly improved.

Dielectric relaxation dynamics of water in model membranes probed by terahertz spectroscopy

We study hydrated model membranes, consisting of stacked bilayers of 1,2-dioleoyl-sn-glycero-3-phosphocholine lipids, using terahertz time-domain spectroscopy and infrared spectroscopy. Terahertz spectroscopy enables the investigation of water dynamics, owing to its sensitivity to dielectric relaxation processes associated with water reorientation. By controlling the number of water molecules per lipid molecule in the system, we elucidate how the interplay between the model membrane and water molecules results in different water dynamics. For decreasing hydration levels, we observe the appearance of new types of water dynamics: the collective bulklike dynamics become less pronounced, whereas an increased amount of both very slowly reorienting (i.e., irrotational) and very rapidly reorienting (i.e., fast) water molecules appear. Temperature-dependent measurements reveal the interconversion between the three distinct types of water present in the system.

Mean-squared atomic displacements in hydrated lysozyme, native and denatured

We use elastic neutron scattering to demonstrate that a sharp increase in the mean-squared atomic displacements, commonly observed in hydrated proteins above 200 K and often referred to as the dynamical transition, is present in the hydrated state of both native and denatured lysozyme. A direct comparison of the native and denatured protein thus confirms that the presence of the transition in the mean-squared atomic displacements is not specific to biologically functional molecules.

Nanosecond motions in proteins impose bounds on the timescale distributions of local dynamics

We elucidate the physics of protein dynamical transition via 10-100-ns molecular dynamics simulations at temperatures spanning 160-300 K. By tracking the energy fluctuations, we show that the protein dynamical transition is marked by a crossover from nonstationary to stationary processes that underlie the dynamics of protein motions. A two-timescale function captures the nonexponential character of backbone structural relaxations. One timescale is attributed to the collective segmental motions and the other to local relaxations. The former is well defined by a single-exponential, nanosecond decay, operative at all temperatures. The latter is described by a set of processes that display a distribution of timescales. Although their average remains on the picosecond timescale, the distribution is markedly contracted at the onset of the transition. It is shown that the collective motions impose bounds on timescales spanned by local dynamical processes. The nonstationary character below the transition implicates the presence of a collection of substates whose interactions are restricted. At these temperatures, a wide distribution of local-motion timescales, extending beyond that of nanoseconds, is observed. At physiological temperatures, local motions are confined to timescales faster than nanoseconds. This relatively narrow window makes possible the appearance of multiple channels for the backbone dynamics to operate.

Dynamical transition in a small helical peptide and its implication for vibrational energy transport

The two-dimensional infrared spectrum of an octameric helical peptide in chloroform was measured as a function of temperature. Isotope labeling of the carbonyl group of one of the amino acids was used to obtain information for an isolated vibration. The antidiagonal width of the 2D-IR signal, which is a measure of the homogeneous dephasing time T(2), is constant from 220 to 260 K (within experimental error), and increases steeply above. The homogeneous dephasing time of the carbonyl vibration is attributed to the flexibility of the system and/or its immediate surrounding. The system undergoes a dynamical transition at about 270 K, with similarities to the protein dynamical transition. Furthermore, the temperature dependence of the antidiagonal width strongly resembles that of the efficiency of vibrational energy transport along the helix, which has been studied in a recent paper (J. Phys. Chem. B 2008, 112, 15487). The connection between the two processes, structural flexibility and energy transport mechanism, is discussed.

Deconstructing the native state: energy landscapes, function, and dynamics of globular proteins

Proteins are highly complex molecules with features exquisitely selected by nature to carry out essential biological functions. Physical chemistry and polymer physics provide us with the tools needed to make sense of this complexity. Upon translation, many proteins fold to a thermodynamically stable form known as the native state. The native state is not static, but consists of a hierarchy of conformations, that are continuously explored through dynamics. In this review we provide a brief introduction to some of the core concepts required in the discussion of the protein native dynamics using energy landscapes ideas. We first discuss recent works which have challenged the structure-function paradigm by demonstrating function in disordered proteins. Next we examine the hierarchical organization in the energy landscapes using atomistic molecular dynamics simulations and principal component analysis. In particular, the role of direct and water-mediated contacts in sculpting the landscape is elaborated. Another approach to studying the native state ensemble is based on choosing high-resolution order parameters for computing one- or two-dimensional free energy surfaces. We demonstrate that 2D free energy surfaces provide rich thermodynamic and kinetic information about the native state ensemble. Brownian dynamics simulations on such a surface indicate that protein conformational dynamics is weakly activated. Finally, we briefly discuss implicit and coarse-grained protein models and emphasize the solvent role in determining native state structure and dynamics.

Bendable, low-loss Topas fibers for the terahertz frequency range

We report on a new class of polymer photonic crystal fibers for low-loss guidance of THz radiation. The use of the cyclic olefin copolymer Topas, in combination with advanced fabrication technology, results in bendable THz fibers with unprecedented low loss and low material dispersion in the THz regime.We demonstrate experimentally how the dispersion may be engineered by fabricating both high- and low-dispersion fibers with zero-dispersion frequency in the regime 0.5-0.6 THz. Near-field, frequency-resolved characterization with high spatial resolution of the amplitude and phase of the modal structure proves that the fiber is single-moded over a wide frequency range, and we see the onset of higher-order modes at high frequencies as well as indication of microporous guiding at low frequencies and high porosity of the fiber. Transmission spectroscopy demonstrates low-loss propagation (< 0.1 dB/cm loss at 0.6 THz) over a wide frequency range.