Mosaic energy landscapes of liquids and the control of protein conformational dynamics by glass-forming solvents

Glycerol-slaved 1H-1H NMR cross-relaxation in quasi-native lysozyme

1H-1H nuclear cross-relaxation experiments have been carried out with lysozyme in variable glycerol viscosity to study intramolecular motion, self-diffusion, and isotropic rigid-body rotational tumbling at 298 K, pH 3.8. Dynamics of intramolecular 1H-1H cross-relaxation rates, the increase in internuclear spatial distances, and lateral and rotational diffusion coefficients all show fractional viscosity dependence with a power law exponent κ in the 0.17-0.83 range. The diffusion coefficient of glycerol Ds with the bulk viscosity itself is non-Stokesian, having a fractional viscosity dependence on the medium viscosity (Ds ∼ η-κ, κ ≈ 0.71). The concurrence and close similarity of the fractional viscosity dependence of glycerol diffusion on the one hand, and diffusion and intramolecular cross-relaxation rates of the protein on the other lead to infer that relaxation of glycerol slaves protein relaxations. Glycerol-transformed native lysozyme to a quasi-native state does not affect the conclusion that both global and internal fluctuations are slaved to glycerol relaxation.

Metastability of Protein Solution Structures in the Absence of a Solvent: Rugged Energy Landscape and Glass-like Behavior

Native ion mobility/mass spectrometry is well-poised to structurally screen proteomes but characterizes protein structures in the absence of a solvent. This raises long-standing unanswered questions about the biological significance of protein structures identified through ion mobility/mass spectrometry. Using newly developed computational and experimental ion mobility/ion mobility/mass spectrometry methods, we investigate the unfolding of the protein ubiquitin in a solvent-free environment. Our data suggest that the folded, solvent-free ubiquitin observed by ion mobility/mass spectrometry exists in a largely native fold with an intact β-grasp motif and α-helix. The ensemble of folded, solvent-free ubiquitin ions can be partitioned into kinetically stable subpopulations that appear to correspond to the structural heterogeneity of ubiquitin in solution. Time-resolved ion mobility/ion mobility/mass spectrometry measurements show that folded, solvent-free ubiquitin exhibits a strongly stretched-exponential time dependence, which simulations trace to a rugged energy landscape with kinetic traps. Unfolding rate constants are estimated to be approximately 800 to 20,000 times smaller than in the presence of water, effectively quenching the unfolding process on the time scale of typical ion mobility/mass spectrometry measurements. Our proposed unfolding pathway of solvent-free ubiquitin shares substantial characteristics with that established for the presence of solvent, including a polarized transition state with significant native content in the N-terminal β-hairpin and α-helix. Our experimental and computational data suggest that (1) the energy landscape governing the motions of folded, solvent-free proteins is rugged in analogy to that of glassy systems; (2) large-scale protein motions may at least partially be determined by the amino acid sequence of a polypeptide chain; and (3) solvent facilitates, rather than controls, protein motions.

Decarbonization in Biobanking: A Potential New Scientific Area

Calls to reduce or entirely remove the carbon footprint of ongoing activities, collectively termed as decarbonization, have become increasingly more vocal in health care with a number of recent, high profile consensus statements. These calls encourage the biobanking field, as one of the foundational health care research infrastructures, to consider decarbonization as a potential novel research area both in terms of the molecules and the equipment used in research. The current article provides a summary of the roundtable discussion during the 2022 ISBER Annual Meeting and Exhibits, highlighting the current knowledge gaps, challenges, and opportunities in this field. In particular, technological innovation, a greater awareness of the current situation, and behavioral change are important pieces of the puzzle to improving the future of decarbonization in biobanking, even if the eventually implemented routes between resource-abundant and resource-restricted settings might be distinctly different. This article sets the foundation for raising awareness of the subject and of subsequent steps that need to be undertaken.

Glass Dynamics Deep in the Energy Landscape

When a liquid is cooled, progress down the energy landscape is arrested near the glass transition temperature T_g. In principle, lower energy states can be accessed by waiting for further equilibration, but the rough energy landscape of glasses quickly leads to kinetics on geologically slow time scales below T_g. Over the past decade, progress has been made probing deeper into the energy landscape via several techniques. By looking at bulk and surface diffusion, using layered deposition that promotes equilibration, imaging glass surfaces with faster dynamics below T_g, and optically exciting glasses, experiments have moved into a regime of ultrastable, low energy glasses that was difficult to access in the past. At the same time, both simulations and energy landscape theory based on a random first order transition (RFOT) have tackled systems that include surfaces, optical excitation, and interfacial dynamics. Here we review some of the recent experimental work, and how energy landscape theory illuminates glassy dynamics well below the glass transition temperature by making direct connections between configurational entropy, energy landscape barriers, and the resulting dynamics.

Implicit water model within the Zimm-Bragg approach to analyze experimental data for heat and cold denaturation of proteins

Studies of biopolymer conformations essentially rely on theoretical models that are routinely used to process and analyze experimental data. While modern experiments allow study of single molecules in vivo, corresponding theories date back to the early 1950s and require an essential update to include the recent significant progress in the description of water. The Hamiltonian formulation of the Zimm-Bragg model we propose includes a simplified, yet explicit model of water-polypeptide interactions that transforms into the equivalent implicit description after performing the summation of solvent degrees of freedom in the partition function. Here we show that our model fits very well to the circular dichroism experimental data for both heat and cold denaturation and provides the energies of inter- and intra-molecular H-bonds, unavailable with other processing methods. The revealed delicate balance between these energies determines the conditions for the existence of cold denaturation and thus clarifies its absence in some proteins.

Moving beyond static snapshots: Protein dynamics and the Protein Data Bank

Proteins are the molecular machines of living systems. Their dynamics are an intrinsic part of their evolutionary selection in carrying out their biological functions. Although the dynamics are more difficult to observe than a static, average structure, we are beginning to observe these dynamics and form sound mechanistic connections between structure, dynamics, and function. This progress is highlighted in case studies from myoglobin and adenylate kinase to the ribosome and molecular motors where these molecules are being probed with a multitude of techniques across many timescales. New approaches to time-resolved crystallography are allowing simple “movies” to be taken of proteins in action, and new methods of mapping the variations in cryo-electron microscopy are emerging to reveal a more complete description of life’s machines. The results of these new methods are aided in their dissemination by continual improvements in curation and distribution by the Protein Data Bank and their partners around the world.

Kramers’ Theory and the Dependence of Enzyme Dynamics on Trehalose-Mediated Viscosity

The disaccharide trehalose is accumulated in the cytoplasm of some organisms in response to harsh environmental conditions. Trehalose biosynthesis and accumulation are important for the survival of such organisms by protecting the structure and function of proteins and membranes. Trehalose affects the dynamics of proteins and water molecules in the bulk and the protein hydration shell. Enzyme catalysis and other processes dependent on protein dynamics are affected by the viscosity generated by trehalose, as described by the Kramers’ theory of rate reactions. Enzyme/protein stabilization by trehalose against thermal inactivation/unfolding is also explained by the viscosity mediated hindering of the thermally generated structural dynamics, as described by Kramers’ theory. The analysis of the relationship of viscosity–protein dynamics, and its effects on enzyme/protein function and other processes (thermal inactivation and unfolding/folding), is the focus of the present work regarding the disaccharide trehalose as the viscosity generating solute. Finally, trehalose is widely used (alone or in combination with other compounds) in the stabilization of enzymes in the laboratory and in biotechnological applications; hence, considering the effect of viscosity on catalysis and stability of enzymes may help to improve the results of trehalose in its diverse uses/applications.

Increasing the Time Resolution of Single-Molecule Experiments with Bayesian Inference

Many time-resolved single-molecule biophysics experiments seek to characterize the kinetics of biomolecular systems exhibiting dynamics that challenge the time resolution of the given technique. Here, we present a general, computational approach to this problem that employs Bayesian inference to learn the underlying dynamics of such systems, even when they are much faster than the time resolution of the experimental technique being used. By accurately and precisely inferring rate constants, our Bayesian inference for the analysis of subtemporal resolution dynamics approach effectively enables the experimenter to super-resolve the poorly resolved dynamics that are present in their data.

Mechanism of Solvent Control of Protein Dynamics

We find that the coupled interactions between protein and water polarization fluctuations play a dominant role in driving the configuration space random walk of solvated proteins. We perform atomistic molecular dynamics simulations on five proteins. Owing to a very low dielectric constant of protein, its dipolar groups experience forces from water along with local forces due to protein atoms. Energy fluctuations reveal a pronounced anticorrelation between protein and water contributions. The protein energy spectrum shows bimodal 1/f noise, which can be attributed to the influence of water on the dynamics of protein.

Small-Angle Neutron Scattering Reveals Energy Landscape for Rhodopsin Photoactivation

Knowledge of the activation principles for G-protein-coupled receptors (GPCRs) is critical to development of new pharmaceuticals. Rhodopsin is the archetype for the largest GPCR family, yet the changes in protein dynamics that trigger signaling are not fully understood. Here we show that rhodopsin can be investigated by small-angle neutron scattering (SANS) in fully protiated detergent micelles under contrast matching to resolve light-induced changes in the protein structure. In SANS studies of membrane proteins, the zwitterionic detergent [(cholamidopropyl)dimethylammonio]-propanesulfonate (CHAPS) is advantageous because of the low contrast difference between the hydrophobic core and hydrophilic head groups as compared with alkyl glycoside detergents. Combining SANS results with quasielastic neutron scattering reveals how changes in volumetric protein shape are coupled (slaved) to the aqueous solvent. Upon light exposure, rhodopsin is swollen by the penetration of water into the protein core, allowing interactions with effector proteins in the visual signaling mechanism.

The role of momentum transfer during incoherent neutron scattering is explained by the energy landscape model

We recently introduced a model of incoherent quasielastic neutron scattering (QENS) that treats the neutrons as wave packets of finite length and the protein as a random walker in the free energy landscape. We call the model ELM for "energy landscape model." In ELM, the interaction of the wave packet with a proton in a protein provides the dynamic information. During the scattering event, the momentum [Formula: see text] is transferred by the wave packet to the struck proton and its moiety, exerting the force [Formula: see text] The resultant energy [Formula: see text] is stored elastically and returned to the neutron as it exits. The energy is given by [Formula: see text], where [Formula: see text] is the ambient temperature and [Formula: see text] ([Formula: see text] 91 K Å) is a new elastobaric coefficient. Experiments yield the scattering intensity (dynamic structure factor) [Formula: see text] as a function of [Formula: see text] and [Formula: see text] To test our model, we use published data on proteins where only thermal vibrations are active. ELM competes with the currently accepted theory, here called the spatial motion model (SMM), which explains [Formula: see text] by motions in real space. ELM is superior to SMM: It can explain the experimental angular and temperature dependence, whereas SMM cannot do so.

Reaction mechanism and reaction coordinates from the viewpoint of energy flow

Reaction coordinates are of central importance for correct understanding of reaction dynamics in complex systems, but their counter-intuitive nature made it a daunting challenge to identify them. Starting from an energetic view of a reaction process as stochastic energy flows biased towards preferred channels, which we deemed the reaction coordinates, we developed a rigorous scheme for decomposing energy changes of a system, both potential and kinetic, into pairwise components. The pairwise energy flows between different coordinates provide a concrete statistical mechanical language for depicting reaction mechanisms. Application of this scheme to the C7eq → C7ax transition of the alanine dipeptide in vacuum revealed novel and intriguing mechanisms that eluded previous investigations of this well studied prototype system for biomolecular conformational dynamics. Using a cost function developed from the energy decomposition components by proper averaging over the transition path ensemble, we were able to identify signatures of the reaction coordinates of this system without requiring any input from human intuition.

A benchmark for reaction coordinates in the transition path ensemble

The molecular mechanism of a reaction is embedded in its transition path ensemble, the complete collection of reactive trajectories. Utilizing the information in the transition path ensemble alone, we developed a novel metric, which we termed the emergent potential energy, for distinguishing reaction coordinates from the bath modes. The emergent potential energy can be understood as the average energy cost for making a displacement of a coordinate in the transition path ensemble. Where displacing a bath mode invokes essentially no cost, it costs significantly to move the reaction coordinate. Based on some general assumptions of the behaviors of reaction and bath coordinates in the transition path ensemble, we proved theoretically with statistical mechanics that the emergent potential energy could serve as a benchmark of reaction coordinates and demonstrated its effectiveness by applying it to a prototypical system of biomolecular dynamics. Using the emergent potential energy as guidance, we developed a committor-free and intuition-independent method for identifying reaction coordinates in complex systems. We expect this method to be applicable to a wide range of reaction processes in complex biomolecular systems.

Dynamical and orientational structural crossovers in low-temperature glycerol

Mean-square displacements of hydrogen atoms in glass-forming materials and proteins, as reported by incoherent elastic neutron scattering, show kinks in their temperature dependence. This crossover, known as the dynamical transition, connects two approximately linear regimes. It is often assigned to the dynamical freezing of subsets of molecular modes at the point of equality between their corresponding relaxation times and the instrumental observation window. The origin of the dynamical transition in glass-forming glycerol is studied here by extensive molecular dynamics simulations. We find the dynamical transition to occur for both the center-of-mass translations and the molecular rotations at the same temperature, insensitive to changes of the observation window. Both the translational and rotational dynamics of glycerol show a dynamic crossover from the structural to a secondary relaxation at the temperature of the dynamical transition. A significant and discontinuous increase in the orientational Kirkwood factor and in the dielectric constant is observed in the same range of temperatures. No indication is found of a true thermodynamic transition to an ordered low-temperature phase. We therefore suggest that all observed crossovers are dynamic in character. The increase in the dielectric constant is related to the dynamic freezing of dipolar domains on the time scale of simulations.

Protein Structural Memory Influences Ligand Binding Mode(s) and Unbinding Rates

The binding of small molecules (e.g., natural ligands, metabolites, and drugs) to proteins governs most biochemical pathways and physiological processes. Here, we use molecular dynamics to investigate the unbinding of dimethyl sulfoxide (DMSO) from two distinct states of a small rotamase enzyme, the FK506-binding protein (FKBP). These states correspond to the FKBP protein relaxed with and without DMSO in the active site. Since the time scale of ligand unbinding (2-20 ns) is faster than protein relaxation (100 ns), a novel methodology is introduced to relax the protein without having to introduce an artificial constraint. The simulation results show that the unbinding time is an order of magnitude longer for dissociation from the DMSO-bound state (holo-relaxed). That is, the actual rate of unbinding depends on the state of the protein, with the protein having a long-lived memory. The rate thus depends on the concentration of the ligand as the apo and holo states reflect low and high concentrations of DMSO, respectively. Moreover, there are multiple binding modes in the apo-relaxed state, while a single binding mode dominates the holo-relaxed state in which DMSO acts as hydrogen bond acceptor from the backbone NH of Ile56, as in the crystal structure of the DMSO/FKBP complex. The solvent relaxes very fast (∼1 ns) close to the NH of Ile56 and with the same time scale of the protein far away from the active site. These results have implications for high-throughput docking, which makes use of a rigid structure of the protein target.

Comprehensive analysis of the green-to-blue photoconversion of full-length Cyanobacteriochrome Tlr0924

Cyanobacteriochromes are members of the phytochrome superfamily of photoreceptors and are of central importance in biological light-activated signaling mechanisms. These photoreceptors are known to reversibly convert between two states in a photoinitiated process that involves a basic E/Z isomerization of the bilin chromophore and, in certain cases, the breakage of a thioether linkage to a conserved cysteine residue in the bulk protein structure. The exact details and timescales of the reactions involved in these photoconversions have not been conclusively shown. The cyanobacteriochrome Tlr0924 contains phycocyanobilin and phycoviolobilin chromophores, both of which photoconvert between two species: blue-absorbing and green-absorbing, and blue-absorbing and red-absorbing, respectively. Here, we followed the complete green-to-blue photoconversion process of the phycoviolobilin chromophore in the full-length form of Tlr0924 over timescales ranging from femtoseconds to seconds. Using a combination of time-resolved visible and mid-infrared transient absorption spectroscopy and cryotrapping techniques, we showed that after photoisomerization, which occurs with a lifetime of 3.6 ps, the phycoviolobilin twists or distorts slightly with a lifetime of 5.3 μs. The final step, the formation of the thioether linkage with the protein, occurs with a lifetime of 23.6 ms.

A wave-mechanical model of incoherent quasielastic scattering in complex systems

Quasielastic incoherent neutron scattering (QENS) is an important tool for the exploration of the dynamics of complex systems such as biomolecules, liquids, and glasses. The dynamics is reflected in the energy spectra of the scattered neutrons. Conventionally these spectra are decomposed into a narrow elastic line and a broad quasielastic band. The band is interpreted as being caused by Doppler broadening due to spatial motion of the target molecules. We propose a quantum-mechanical model in which there is no separate elastic line. The quasielastic band is composed of sharp lines with twice the natural line width, shifted from the center by a random walk of the protein in the free-energy landscape of the target molecule. The walk is driven by vibrations and by external fluctuations. We first explore the model with the Mössbauer effect. In the subsequent application to QENS we treat the incoming neutron as a de Broglie wave packet. While the wave packet passes the protons in the protein and the hydration shell it exchanges energy with the protein during the passage time of about 100 ns. The energy exchange broadens the ensemble spectrum. Because the exchange involves the free-energy landscape of the protein, the QENS not only provides insight into the protein dynamics, but it may also illuminate the free-energy landscape of the protein-solvent system.

Crowding induced collective hydration of biological macromolecules over extended distances

Ultrafast two-dimensional infrared (2D-IR) spectroscopy reveals picosecond protein and hydration dynamics of crowded hen egg white lysozyme (HEWL) labeled with a metal-carbonyl vibrational probe covalently attached to a solvent accessible His residue. HEWL is systematically crowded alternatively with polyethylene glycol (PEG) or excess lysozyme in order to distinguish the chemically inert polymer from the complex electrostatic profile of the protein crowder. The results are threefold: (1) A sharp dynamical jamming-like transition is observed in the picosecond protein and hydration dynamics that is attributed to an independent-to-collective hydration transition induced by macromolecular crowding that slows the hydration dynamics up to an order of magnitude relative to bulk water. (2) The interprotein distance at which the transition occurs suggests collective hydration of proteins over distances of 30-40 Å. (3) Comparing the crowding effects of PEG400 to our previously reported experiments using glycerol exposes fundamental differences between small and macromolecular crowding agents.

Dielectric spectroscopy study of myoglobin in glycerol-water mixtures

Due to the interest in protein dynamics, there are numerous dielectric relaxation studies of proteins in water and in glass-forming aqueous solvents such as glycerol-water mixtures. In the regime of low frequencies, the inevitable dc-conductivity of such systems limits the resolution of dynamics that are slow compared with the solvent relaxation. Solutions of myoglobin in glycerol/water mixtures of various compositions are measured by dielectric spectroscopy in the frequency range from 10mHz to 10MHz. The resolution of low frequency modes is improved by two approaches: electrical 'cleaning' and the analysis of the derivative of the real component of permittivity, which shows no direct signature of dc-conductivity. Effects of internal interfacial polarization are also addressed by measuring the same solvents in confinement as well as mixed with glass beads. We find two processes, the structural relaxation of the solvent and the slower rotational mode of the protein, with no indication at even lower frequencies of a dielectric signature of fluctuations associated with protein dynamics.

Reverse micelles as a tool for probing solvent modulation of protein dynamics: Reverse micelle encapsulated hemoglobin

Hydration waters impact protein dynamics. Dissecting the interplay between hydration waters and dynamics requires a protein that manifests a broad range of dynamics. Proteins in reverse micelles (RMs) have promise as tools to achieve this objective because the water content can be manipulated. Hemoglobin is an appropriate tool with which to probe hydration effects. We describe both a protocol for hemoglobin encapsulation in reverse micelles and a facile method using PEG and cosolvents to manipulate water content. Hydration properties are probed using the water-sensitive fluorescence from Hb bound pyranine and covalently attached Badan. Protein dynamics are probed through ligand recombination traces derived from photodissociated carbonmonoxy hemoglobin on a log scale that exposes the potential role of both α and β solvent fluctuations in modulating protein dynamics. The results open the possibility of probing hydration level phenomena in this system using a combination of NMR and optical probes.

Biomolecular electrostatics and solvation: a computational perspective

An understanding of molecular interactions is essential for insight into biological systems at the molecular scale. Among the various components of molecular interactions, electrostatics are of special importance because of their long-range nature and their influence on polar or charged molecules, including water, aqueous ions, proteins, nucleic acids, carbohydrates, and membrane lipids. In particular, robust models of electrostatic interactions are essential for understanding the solvation properties of biomolecules and the effects of solvation upon biomolecular folding, binding, enzyme catalysis, and dynamics. Electrostatics, therefore, are of central importance to understanding biomolecular structure and modeling interactions within and among biological molecules. This review discusses the solvation of biomolecules with a computational biophysics view toward describing the phenomenon. While our main focus lies on the computational aspect of the models, we provide an overview of the basic elements of biomolecular solvation (e.g. solvent structure, polarization, ion binding, and non-polar behavior) in order to provide a background to understand the different types of solvation models.

Coupling of protein and environment fluctuations

We review the concepts of protein dynamics developed over the last 35years and extend applications of the unified model of protein dynamics to heat flow and spatial fluctuations in hydrated myoglobin (Mb) powders. Differential scanning calorimetry (DSC) and incoherent neutron scattering (INS) data on hydration Mb powders are explained by the temperature-dependence of the hydration-shell β(h) process measured by dielectric relaxation spectroscopy (DRS). The unified model explains the temperature dependence of DSC and INS data as a kinetic effect due to a fixed experimental time window and a broad distribution of hydration-shell β(h) fluctuation rates. We review the slaving of large scale protein motions to the bulk solvent α process, and the metastability of Mb molecules in glass forming bulk solvent at low temperatures. This article is part of a Special Issue entitled: "Protein Dynamics: Experimental and Computational Approaches".

Effect of solvation-related interaction on the low-temperature dynamics of proteins

The effect of solvation-related interaction on the low-temperature dynamics of proteins is studied by taking into account the desolvation barriers in the interactions of native contacts. It is found out that about the folding transition temperature, the protein folds in a cooperative manner, and the water molecules are expelled from the hydrophobic core at the final stage in the folding process. At low temperature, however, the protein would generally be trapped in many metastable conformations with some water molecules frozen inside the protein. The desolvation takes an important role in these processes. The number of frozen water molecules and that of frozen states of proteins are further analyzed with the methods based on principal component analysis (PCA) and the clustering of conformations. It is found out that both the numbers of frozen water molecules and the frozen states of the protein increase quickly below a certain temperature. Especially, the number of frozen states of the protein increases exponentially following the decrease in the temperature, which resembles the basic features of glassy dynamics. Interestingly, it is observed that the freezing of water molecules and that of protein conformations happen at almost the same temperature. This suggests that the solvation-related interaction performs an important role for the low-temperature dynamics of the model protein.

Simultaneous measurements of solvent dynamics and functional kinetics in a light-activated enzyme

Solvent fluctuations play a key role in controlling protein motions and biological function. Here, we have studied how individual steps of the reaction catalyzed by the light-activated enzyme protochlorophyllide oxidoreductase (POR) couple with solvent dynamics. To simultaneously monitor the catalytic cycle of the enzyme and the dynamical behavior of the solvent, we designed temperature-dependent UV-visible microspectrophotometry experiments, using flash-cooled nanodroplets of POR to which an exogenous soluble fluorophore was added. The formation and decay of the first two intermediates in the POR-catalyzed reaction were measured, together with the solvent glass transition and the buildup of crystalline ice at cryogenic temperatures. We find that formation of the first intermediate occurs below the glass transition temperature (T(g)), and is not affected by changes in solvent dynamics induced by modifying the glycerol content. In contrast, formation of the second intermediate occurs above T(g) and is influenced by changes in glycerol concentration in a manner remarkably similar to the buildup of crystalline ice. These results suggest that internal, nonslaved protein motions drive the first step of the POR-catalyzed reaction whereas solvent-slaved motions control the second step. We propose that the concept of solvent slaving applies to complex enzymes such as POR.

A unified model of protein dynamics

Protein functions require conformational motions. We show here that the dominant conformational motions are slaved by the hydration shell and the bulk solvent. The protein contributes the structure necessary for function. We formulate a model that is based on experiments, insights from the physics of glass-forming liquids, and the concepts of a hierarchically organized energy landscape. To explore the effect of external fluctuations on protein dynamics, we measure the fluctuations in the bulk solvent and the hydration shell with broadband dielectric spectroscopy and compare them with internal fluctuations measured with the Mössbauer effect and neutron scattering. The result is clear. Large-scale protein motions are slaved to the fluctuations in the bulk solvent. They are controlled by the solvent viscosity, and are absent in a solid environment. Internal protein motions are slaved to the beta fluctuations of the hydration shell, are controlled by hydration, and are absent in a dehydrated protein. The model quantitatively predicts the rapid increase of the mean-square displacement above approximately 200 K, shows that the external beta fluctuations determine the temperature- and time-dependence of the passage of carbon monoxide through myoglobin, and explains the nonexponential time dependence of the protein relaxation after photodissociation.

On the surface of glasses

Dynamics near the surface of glasses is generally much faster than in the bulk. Neglecting static perturbations of structure at the surface, we use random first order transition (RFOT) theory to show the free energy barrier for activated motion near a free surface should be half that of the bulk at the same temperature. The increased mobility allows the surface layers to descend much further on the energy landscape than the bulk ordinarily does. The simplified RFOT calculation, however, predicts a limiting value for the configurational entropy a vapor deposited glass may reach as a function of deposition rate. We sketch how mode coupling effects extend the excess free surface mobility into the bulk so that the glass transition temperature is measurably perturbed at depths greater than the naive length scale of dynamic cooperativity.

Quantum efficiency distributions of photo-induced side-pathway donor oxidation at cryogenic temperature in photosystem II

We monitored illuminated-minus-dark absorption difference spectra in the range of 450-1100 nm induced by continuous illumination at 8 K of photosystem II (PSII) core complexes from Thermosynechococcus elongatus. The photo-induced oxidation of the side-path donors Cytb(559), beta-carotene and chlorophyll Z, as well as the concomitant stable (t(1/2) > 1 s) reduction of the first plastoquinone electron acceptor, Q(A) (monitored by the well-known 'C550' shift), were quantified as a function of the absorbed photons per PSII. The Q(A) photo-induced reduction data can be described by three distinct quantum efficiency distributions: (i) a very high efficiency of approximately 0.5-1, (ii) a middle efficiency with a very large range of approximately 0.014-0.2, and (iii) a low efficiency of approximately 0.002. Each of the observed side-path donors exhibited similar quantum efficiency distributions, which supports a branched pathway model for side-path oxidation where beta-carotene is the immediate electron donor to the photo-oxidized chlorophyll (P680(+)). The yields of the observed side-path donors account quantitatively for the wide middle efficiency range of photo-induced Q(A) reduction, but not for the PSII fractions that exhibit the highest and lowest efficiencies. The high-efficiency component may be due to Tyr(Z) oxidation. A donor that does not exhibit an identified absorption in the visible-near-IR region is mainly responsible for the lowest efficiency component.

Temperature and pressure dependence of protein stability: the engineered fluorescein-binding lipocalin FluA shows an elliptic phase diagram

We have measured the equilibrium constant for the denaturation transition of the engineered fluorescein-binding lipocalin FluA as a function of pressure and temperature, taking advantage of the fact that the ligand's fluorescence is almost fully quenched when complexed with the folded protein, but reversibly reappears on denaturation. From the equilibrium constant as a function of pressure and temperature all of the involved thermodynamic parameters of protein folding, in particular the changes in entropy and volume, compressibility, thermal expansion, and specific heat, were deduced in a global fitting procedure. Assuming that these parameters are independent of temperature and pressure, we can demonstrate from the ratio of Deltabeta, Deltaalpha(2), DeltaC(p) that the phase diagram of protein folding assumes an elliptic shape. Furthermore, we can show that the thermodynamic condition for such an elliptic phase diagram is related to the degree of correlation between the fluctuations of the changes in volume and enthalpy at the phase boundary. For the protein investigated this correlation is low, as generally expected for highly degenerate systems. Our study suggests that the elliptic phase diagram is a consequence of the inherent conformational disorder of proteins and that it may be viewed as the thermodynamic manifestation of the high degeneracy of conformational energies that is characteristic for this class of macromolecules.

Quaternary relaxations in sol-gel encapsulated hemoglobin studied via NIR and UV spectroscopy

In this work, we study the kinetics of the R --> T transition in hemoglobin using a combination of near-infrared and near-ultraviolet spectroscopy. We use a sol-gel encapsulation protocol to decelerate the conformational transitions and to avoid spectral perturbations arising from ligand migration and recombination. We monitor two spectroscopic markers: band III in the near-IR, which is a fine probe of the heme pocket conformation, and the tryptophan band in the near-UV, which probes the formation of the Trpbeta37-Aspalpha94 hydrogen bond, characteristic of the T structure, at the critical alpha1beta2 subunit interface. The time evolution of these two bands is monitored after deoxygenation of encapsulated oxyhemoglobin, obtained by diffusion of a reducing agent into the porous silica matrix. Characteristic spectral shifts are observed: comparison with myoglobin enables us to assign them to quaternary structure relaxations. Band III spectral relaxation is clearly nonexponential, and analysis with the Maximum Entropy Method enables us to identify three processes. On the other hand, near-UV spectral relaxation follows an exponential decay with a time constant closely corresponding to the second process observed in the near IR. Very interestingly, the rates of all processes markedly depend on the viscosity of the co-encapsulated solvent, following a power law. Our results reveal correlations between heme pocket relaxations, induced by the R --> T transition, and structural event(s) occurring at the alpha1beta2 interface and highlight their solvent dependence. The power law viscosity dependence of relaxation rates suggests that the observed protein relaxations are "slaved" to the co-encapsulated solvent. The stepwise character of the quaternary transition is also evidenced.

Ligand recombination and a hierarchy of solvent slaved dynamics: the origin of kinetic phases in hemeproteins

Ligand recombination studies play a central role both for characterizing different hemeproteins and their conformational states but also for probing fundamental biophysical processes. Consequently, there is great importance to providing a foundation from which one can understand the physical processes that give rise to and modulate the large range of kinetic patterns associated with ligand recombination in myoglobins and hemoglobins. In this work, an overview of cryogenic and solution phase recombination phenomena for COMb is first reviewed and then a new paradigm is presented for analyzing the temperature and viscosity dependent features of kinetic traces in terms of multiple phases that reflect which tier(s) of solvent slaved protein dynamics is (are) operative on the photoproduct population during the time course of the measurement. This approach allows for facile inclusion of both ligand diffusion among accessible cavities and conformational relaxation effects. The concepts are illustrated using kinetic traces and MEM populations derived from the CO recombination process for wild type and mutant myoglobins either in sol-gel matrices bathed in glycerol or in trehalose-derived glassy matrices.

Minimizing frustration by folding in an aqueous environment

Although life as we know it evolved in an aqueous medium, the properties of water are not completely understood. In this review, we focus on the role of water in guiding protein folding and stability. Specifically, we discuss the mechanisms of protein folding in an aqueous environment, the effects of water on the folding energy landscape as well as the transition state ensemble, and interactions of water with the folded state. We show that water cannot be viewed as a passive solvent, but rather, plays a very active role in the life of a protein.

Charge and momentum transfer in supercooled melts: Why should their relaxation times differ?

The steady-state values of the viscosity and the intrinsic ionic conductivity of quenched melts are computed, in terms of independently measurable quantities. The frequency dependence of the ac dielectric response is estimated. The discrepancy between the corresponding characteristic relaxation times is only apparent; it does not imply distinct mechanisms, but stems from the intrinsic barrier distribution for alpha-relaxation in supercooled fluids and glasses. This type of intrinsic "decoupling" is argued not to exceed four orders in magnitude for known glassformers. The origin of the discrepancy between the stretching exponent beta, as extracted from epsilon(omega) and the dielectric modulus data, is explained. The actual width of the barrier distribution always grows with lowering the temperature. The contrary is an artifact of the large contribution of the dc-conductivity component to the modulus data. The methodology allows one to single out other contributions to the conductivity, as in "superionic" liquids or when charge carriers are delocalized, implying that in those systems, charge transfer does not require structural reconfiguration.

Nanosecond relaxation dynamics of protein GB1 identified by the time-dependent red shift in the fluorescence of tryptophan and 5-fluorotryptophan

The B1 domain of Streptococcal protein G (GB1) is a small, thermostable protein containing a single tryptophan residue. We recorded time-resolved fluorescence of the wild-type GB1 and its 5-fluorotryptophan (5FTrp) variant at more than 30 emission wavelengths between 300 and 470 nm. The time-resolved emission spectra reveal no signs of heterogeneity, but show a time-dependent red shift characteristic of microscopic dielectric relaxation. This is true for both 5FTrp and unmodified Trp in GB1. The time-dependent red shifts in the fluorescence of 5FTrp and unmodified Trp are essentially identical, confirming that the shift is caused by the relaxation of the protein matrix rather than by the fluorophore itself. The total amplitude (but not the rate) of the time-dependent red shift depends on the fluorophore, specifically, on the magnitude of the vector difference between its excited state and ground state electric dipole moments; for 5FTrp this is estimated to be about 88% of that for the unmodified Trp. The decay of the excited state fluorophore population is not monoexponential for either fluorophore; however, the deviation from the monoexponential decay law is larger in the case of unmodified Trp. The relaxation dynamics of GB1 was found to be considerably faster than that of other proteins studied previously, consistent with the small size, tightly packed core, and high thermodynamic stability of GB1.

Protein folding is slaved to solvent motions

Proteins, the workhorses of living systems, are constructed from chains of amino acids, which are synthesized in the cell based on the instructions of the genetic code and then folded into working proteins. The time for folding varies from microseconds to hours. What controls the folding rate is hotly debated. We postulate here that folding has the same temperature dependence as the alpha-fluctuations in the bulk solvent but is much slower. We call this behavior slaving. Slaving has been observed in folded proteins: Large-scale protein motions follow the solvent fluctuations with rate coefficient k(alpha) but can be slower by a large factor. Slowing occurs because large-scale motions proceed in many small steps, each determined by k(alpha). If conformational motions of folded proteins are slaved, so a fortiori must be the motions during folding. The unfolded protein makes a Brownian walk in the conformational space to the folded structure, with each step controlled by k(alpha). Because the number of conformational substates in the unfolded protein is extremely large, the folding rate coefficient, k(f), is much smaller than k(alpha). The slaving model implies that the activation enthalpy of folding is dominated by the solvent, whereas the number of steps n(f) = k(alpha)/k(f) is controlled by the number of accessible substates in the unfolded protein and the solvent. Proteins, however, undergo not only alpha- but also beta-fluctuations. These additional fluctuations are local protein motions that are essentially independent of the bulk solvent fluctuations and may be relevant at late stages of folding.