A unified model of protein dynamics

Decoupling of the onset of anharmonicity between a protein and its surface water around 200 K

The protein dynamical transition at ~200 K, where the biomolecule transforms from a harmonic, non-functional form to an anharmonic, functional state, has been thought to be slaved to the thermal activation of dynamics in its surface hydration water. Here, by selectively probing the dynamics of protein and hydration water using elastic neutron scattering and isotopic labeling, we found that the onset of anharmonicity in the two components around 200 K is decoupled. The one in protein is an intrinsic transition, whose characteristic temperature is independent of the instrumental resolution time, but varies with the biomolecular structure and the amount of hydration, while the one of water is merely a resolution effect.

Protein Retrieval via Integrative Molecular Ensembles (PRIME) through Extended Similarity Indices

Molecular dynamics (MD) simulations are ideally suited to describe conformational ensembles of biomolecules such as proteins and nucleic acids. Microsecond-long simulations are now routine, facilitated by the emergence of graphical processing units. Clustering, which groups objects based on structural similarity, is typically used to process ensembles, leading to different states, their populations, and the identification of representative structures. A popular pipeline combines hierarchical clustering for clustering and selecting the cluster centroid as representative of the cluster. Here, we propose to improve on this approach, by developing a module-Protein Retrieval via Integrative Molecular Ensembles (PRIME), that consists of tools to improve the prediction of the representative in the most populated cluster using extended continuous similarity. PRIME is integrated with our Molecular Dynamics Analysis with N-ary Clustering Ensembles (MDANCE) package and can be used as a postprocessing tool for arbitrary clustering algorithms, compatible with several MD suites. PRIME predictions produced structures that when aligned to the experimental structure were better superposed (lower RMSD). A further benefit of PRIME is its linear scaling─rather than the traditional O(N2) traditionally associated with comparisons of elements in a set.

Influence of the Water Model on the Structure and Interactions of the GPR40 Protein with the Lipid Membrane and the Solvent: Rigid versus Flexible Water Models

G protein-coupled receptors (GPCR) are responsible for modulating various physiological functions and are thus related to the pathophysiology of different diseases. Being potential therapeutic targets, multiple computational methodologies have been developed to analyze their behavior and interactions with other species. The solvent, on the other hand, has received much less attention. In this work, we analyzed the effect of four explicit water models on the structure and interactions of the GPR40 receptor in its apo form. We employed the rigid SPC/E and TIP4P models, and their flexible versions, the FBA/ϵ and TIP4P/ϵflex. We explored the structural changes and their correlation with some bulk dynamic properties of water. Our results showed an adverse effect on the conservation of the secondary structure of the receptor with all the models due to the breaking of the intramolecular hydrogen bond network, being more evident for the TIP4P models. Notably, all four models brought the receptor to states similar to the active one, modifying the intracellular part of the TM5 and TM6 domains in a "hinge" type movement, allowing the opening of the structure. Regarding the dynamic properties, the rigid models showed results comparable to those obtained in other studies on membrane systems. However, flexible models exhibit disparities in the molecular representation of systems. Surprisingly, the FBA/ϵ model improves the molecular picture of several properties, even though their agreement with bulk diffusion is poorer. These findings reinforce our idea that exploring other water models or improving the current ones, to better represent the membrane interface, can lead to a positive impact on the description of the signal transduction mechanisms and the search of new drugs by targeting these receptors.

Essential dynamics of ubiquitin in water and in a natural deep eutectic solvent

Natural deep eutectic solvents (NADESs) comprised of osmolytes are of interest as potential biomolecular (cryo)protectants. However, the way these solvents influence the structure and dynamics of biomolecules as well as the role of water remains poorly understood. We carried out principal component analysis of various secondary structure elements of ubiquitin in water and a betaine : glycerol : water (1 : 2 : ζ; ζ = 0, 1, 2, 5, 10, 20, 45) NADES, from molecular dynamics trajectories, to gain insight into the protein dynamics as it undergoes a transition from a highly viscous anhydrous to an aqueous environment. A crossover of the protein's essential dynamics at ζ ∼ 5, induced by solvent-shell coupled fluctuations, is observed, indicating that ubiquitin might (re)fold in the NADES upon water addition at ζ > ∼5. Further, in contrast to water, the anhydrous NADES preserves ubiquitin's essential modes at high temperatures explaining the protein's seemingly enhanced thermal stability.

Deeper Insights into Mixed Crowding through Enzyme Activity, Dynamics, and Crowder Diffusion

Given the fact that the cellular interior is crowded by many different kinds of macromolecules, it is important that in vitro studies be carried out in the presence of mixed crowder systems. In this regard, we have used binary crowders formed by the combination of some of the commonly used crowding agents, namely, Ficoll 70, Dextran 70, Dextran 40, and PEG 8000 (PEG 8), to study how these affect enzyme activity, dynamics, and crowder diffusion. The enzyme chosen is AK3L1, an isoform of adenylate kinase. To investigate its dynamics, we have carried out three single point mutations (A74C, A132C, and A209C) with the cysteine residues being labeled with a coumarin-based solvatochromic probe [CPM: (7-diethylamino-3-(4-maleimido-phenyl)-4-methylcoumarin)]. Both enzyme activity and dynamics decreased in the binary mixtures as compared with the sum of the individual crowders, suggesting a reduction in excluded volume (in the mixture). To gain deeper insights into the binary mixtures, fluorescence correlation spectroscopy studies were carried out using fluorescein isothiocyanate-labeled Dextran 70 and tetramethylrhodamine-labeled AK3L1 as the diffusion probes. Diffusion in binary mixtures was observed to be much more constrained (relative to the sum of the individual crowders) for the labeled enzyme as compared to the labeled crowder showing different environments being faced by the two species. This was further confirmed during imaging of the phase-separated droplets formed in the binary mixtures having PEG as one of the crowding agents. The interior of these droplets was found to be rich in crowders and densely packed, as shown by confocal and digital holographic microscopy images, with the enzymes predominantly residing outside these droplets, that is, in the relatively less crowded regions. Taken together, our data provide important insights into various aspects of the simplest form of mixed crowding, that is, composed of just two components, and also hint at the enhanced complexity that the cellular interior presents toward having a detailed and comprehensive understanding of the same.

Exploring the Synthesis of Self-Organization and Active Motion

Proteins, genetic material, and membranes are fundamental to all known organisms, yet these components alone do not constitute life. Life emerges from the dynamic processes of self-organization, assembly, and active motion, suggesting the existence of similar artificial systems. Against this backdrop, our Perspective explores a variety of chemical phenomena illustrating how nonequilibrium self-organization and micromotors contribute to life-like behavior and functionalities. After explaining key terms, we discuss specific examples including enzymatic motion, diffusiophoretic and bubble-driven self-propulsion, pattern-forming reaction-diffusion systems, self-assembling inorganic aggregates, and hierarchically emergent phenomena. We also provide a roadmap for combining self-organization and active motion and discuss possible outcomes through biological analogs. We suggest that this research direction, deeply rooted in physical chemistry, offers opportunities for further development with broad impacts on related sciences and technologies.

Comparison of Sucrose and Trehalose for Protein Stabilization Using Differential Scanning Calorimetry

The disaccharide trehalose is generally acknowledged as a superior stabilizer of proteins and other biomolecules in aqueous environments. Despite many theories aiming to explain this, the stabilization mechanism is still far from being fully understood. This study compares the stabilizing properties of trehalose with those of the structurally similar disaccharide sucrose. The stability has been evaluated for the two proteins, lysozyme and myoglobin, at both low and high temperatures by determining the glass transition temperature, Tg, and the denaturation temperature, Tden. The results show that the sucrose-containing samples exhibit higher Tden than the corresponding trehalose-containing samples, particularly at low water contents. The better stabilizing effect of sucrose at high temperatures may be explained by the fact that sucrose, to a greater extent, binds directly to the protein surface compared to trehalose. Both sugars show Tden elevation with an increasing sugar-to-protein ratio, which allows for a more complete sugar shell around the protein molecules. Finally, no synergistic effects were found by combining trehalose and sucrose. Conclusively, the exact mechanism of protein stabilization may vary with the temperature, as influenced by temperature-dependent interactions between the protein, sugar, and water. This variability can make trehalose to a superior stabilizer under some conditions and sucrose under others.

Membrane mimetic-dependence of GPCR energy landscapes

We leveraged variable-temperature 19F-NMR spectroscopy to compare the conformational equilibria of the human A2A adenosine receptor (A2AAR), a class A G protein-coupled receptor (GPCR), across a range of temperatures ranging from lower temperatures typically employed in 19F-NMR experiments to physiological temperature. A2AAR complexes with partial agonists and full agonists showed large increases in the population of a fully active conformation with increasing temperature. NMR data measured at physiological temperature were more in line with functional data. This was pronounced for complexes with partial agonists, where the population of active A2AAR was nearly undetectable at lower temperature but became evident at physiological temperature. Temperature-dependent behavior of complexes with either full or partial agonists exhibited a pronounced sensitivity to the specific membrane mimetic employed. Cellular signaling experiments correlated with the temperature-dependent conformational equilibria of A2AAR in lipid nanodiscs but not in some detergents, underscoring the importance of the membrane environment in studies of GPCR function.

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.

The inhibition of fibril formation of lysozyme by sucrose and trehalose

The two disaccharides, trehalose and sucrose, have been compared in many studies due to their structural similarity. Both possess the ability to stabilise and reduce aggregation of proteins. Trehalose has also been shown to inhibit the formation of highly structured protein aggregates called amyloid fibrils. This study aims to compare how the thermal stability of the protein lysozyme at low pH (2.0 and 3.5) is affected by the presence of the two disaccharides. We also address the anti-aggregating properties of the disaccharides and their inhibitory effects on fibril formation. Differential scanning calorimetry confirms that the thermal stability of lysozyme is increased by the presence of trehalose or sucrose. The effect is slightly larger for sucrose. The inhibiting effects on protein aggregation are investigated using small-angle X-ray scattering which shows that the two-component system consisting of lysozyme and water (Lys/H2O) at pH 2.0 contains larger aggregates than the corresponding system at pH 3.5 as well as the sugar containing systems. In addition, the results show that the particle-to-particle distance in the sugar containing systems (Lys/Tre/H2O and Lys/Suc/H2O) at pH 2.0 is longer than at pH 3.5, suggesting larger protein aggregates in the former. Finally, the characteristic distance separating β-strands in amyloid fibrils is observed for the Lys/H2O system at pH 2.0, using wide-angle X-ray scattering, while it is not clearly observed for the sugar containing systems. This study further shows that the two disaccharides stabilise the native fold of lysozyme by increasing the denaturation temperature. However, other factors, such as a weakening of hydrophobic interactions and hydrogen bonding between proteins, might also play a role in their inhibitory effect on amyloid fibril formation.

Identification of the Thermal Activation Network in Human 15-Lipoxygenase-2: Divergence from Plant Orthologs and Its Relationship to Hydrogen Tunneling Activation Barriers

The oxidation of polyunsaturated fatty acids by lipoxygenases (LOXs) is initiated by a C-H cleavage step in which the hydrogen atom is transferred quantum mechanically (i.e., via tunneling). In these reactions, protein thermal motions facilitate the conversion of ground-state enzyme-substrate complexes to tunneling-ready configurations and are thus important for transferring energy from the solvent to the active site for the activation of catalysis. In this report, we employed temperature-dependent hydrogen-deuterium exchange mass spectrometry (TDHDX-MS) to identify catalytically linked, thermally activated peptides in a representative animal LOX, human epithelial 15-LOX-2. TDHDX-MS of wild-type 15-LOX-2 was compared to two active site mutations that retain structural stability but have increased activation energies (Ea) of catalysis. The Ea value of one variant, V427L, is implicated to arise from suboptimal substrate positioning by increased active-site side chain rotamer dynamics, as determined by X-ray crystallography and ensemble refinement. The resolved thermal network from the comparative Eas of TDHDX-MS between wild-type and V426A is localized along the front face of the 15-LOX-2 catalytic domain. The network contains a clustering of isoleucine, leucine, and valine side chains within the helical peptides. This thermal network of 15-LOX-2 is different in location, area, and backbone structure compared to a model plant lipoxygenase from soybean that exhibits a low Ea value of catalysis compared to the human ortholog. The presented data provide insights into the divergence of thermally activated protein motions in plant and animal LOXs and their relationships to the enthalpic barriers for facilitating hydrogen tunneling.

Mechanism of Protein Aggregation Inhibition by Arginine: Blockage of Anionic Side Chains Favors Unproductive Encounter Complexes

Aggregation refers to the assembly of proteins into nonphysiological higher order structures. While amyloid has been studied extensively, much less is known about amorphous aggregation, a process that interferes with protein expression and storage. Free arginine (Arg+) is a widely used aggregation inhibitor, but its mechanism remains elusive. Focusing on myoglobin (Mb), we recently applied atomistic molecular dynamics (MD) simulations for gaining detailed insights into amorphous aggregation (Ng J. Phys. Chem. B 2021, 125, 13099). Building on that approach, the current work for the first time demonstrates that MD simulations can directly elucidate aggregation inhibition mechanisms. Comparative simulations with and without Arg+ reproduced the experimental finding that Arg+ significantly decreased the Mb aggregation propensity. Our data reveal that, without Arg+, protein-protein encounter complexes readily form salt bridges and hydrophobic contacts, culminating in firmly linked dimeric aggregation nuclei. Arg+ promotes the dissociation of encounter complexes. These "unproductive" encounter complexes are favored because Arg+ binding to D- and E- lowers the tendency of these anionic residues to form interprotein salt bridges. Side chain blockage is mediated largely by the guanidinium group of Arg+, which binds carboxylates through H-bond-reinforced ionic contacts. Our MD data revealed Arg+ self-association into a dynamic quasi-infinite network, but we found no evidence that this self-association is important for protein aggregation inhibition. Instead, aggregation inhibition by Arg+ is similar to that mediated by free guanidinium ions. The computational strategy used here should be suitable for the rational design of aggregation inhibitors with enhanced potency.

Structural and physical basis for the elasticity of elastin

Elastin function is to endow vertebrate tissues with elasticity so that they can adapt to local mechanical constraints. The hydrophobicity and insolubility of the mature elastin polymer have hampered studies of its molecular organisation and structure-elasticity relationships. Nevertheless, a growing number of studies from a broad range of disciplines have provided invaluable insights, and several structural models of elastin have been proposed. However, many questions remain regarding how the primary sequence of elastin (and the soluble precursor tropoelastin) governs the molecular structure, its organisation into a polymeric network, and the mechanical properties of the resulting material. The elasticity of elastin is known to be largely entropic in origin, a property that is understood to arise from both its disordered molecular structure and its hydrophobic character. Despite a high degree of hydrophobicity, elastin does not form compact, water-excluding domains and remains highly disordered. However, elastin contains both stable and labile secondary structure elements. Current models of elastin structure and function are drawn from data collected on tropoelastin and on elastin-like peptides (ELPs) but at the tissue level, elasticity is only achieved after polymerisation of the mature elastin. In tissues, the reticulation of tropoelastin chains in water defines the polymer elastin that bears elasticity. Similarly, ELPs require polymerisation to become elastic. There is considerable interest in elastin especially in the biomaterials and cosmetic fields where ELPs are widely used. This review aims to provide an up-to-date survey of/perspective on current knowledge about the interplay between elastin structure, solvation, and entropic elasticity.

Bridging Soft Interaction and Excluded Volume in Crowded Milieu through Subtle Protein Dynamics

The impact of macromolecular crowding on biological macromolecules has been elucidated through the excluded volume phenomenon and soft interactions. However, it has often been difficult to provide a clear demarcation between the two regions. Here, using temperature-dependent dynamics (local and global) of the multidomain protein human serum albumin (HSA) in the presence of commonly used synthetic crowders (Dextran 40, PEG 8, Ficoll 70, and Dextran 70), we have shown the presence of a transition that serves as a bridge between the soft and hard regimes. The bridging region is independent of the crowder identity and displays no apparent correlation with the critical overlap concentration of the polymeric crowding agents. Moreover, the dynamics of domains I and II and the protein gating motion respond differently, thereby bringing to the fore the asymmetry underlying the crowder influence on HSA. In addition, solvent-coupled and decoupled protein motions indicate the heterogeneity of the dynamic landscape in the crowded milieu. We also propose an intriguing correlation between protein stability and dynamics, with increased global stability being accompanied by eased local domain motion.

Osmotic stress studies of G-protein-coupled receptor rhodopsin activation

We summarize and critically review osmotic stress studies of the G-protein-coupled receptor rhodopsin. Although small amounts of structural water are present in these receptors, the effect of bulk water on their function remains uncertain. Studies of the influences of osmotic stress on the GPCR archetype rhodopsin have given insights into the functional role of water in receptor activation. Experimental work has discovered that osmolytes shift the metarhodopsin equilibrium after photoactivation, either to the active or inactive conformations according to their molar mass. At least 80 water molecules are found to enter rhodopsin in the transition to the photoreceptor active state. We infer that this movement of water is both necessary and sufficient for receptor activation. If the water influx is prevented, e.g., by large polymer osmolytes or by dehydration, then the receptor functional transition is back shifted. These findings imply a new paradigm in which rhodopsin becomes solvent swollen in the activation mechanism. Water thus acts as an allosteric modulator of function for rhodopsin-like receptors in lipid membranes.

Protein stability in a natural deep eutectic solvent: Preferential hydration or solvent slaving?

Deep eutectic solvents (DESs) emerged as potential alternative solvent media in multiple areas, including biomolecular (cryo)preservation. Herein, we studied the stability of a small protein (ubiquitin) in water and a betaine-glycerol-water (B:G:W) (1:2:ζ; ζ = 0, 1, 2, 5, 10) DES, through molecular dynamics. An AMBER-based model that accurately describes the density and shear viscosity of the DES is proposed. We find that water molecules are largely trapped in the solvent, precluding the formation of a full hydration layer, seemingly opposite to osmolytes' preferential exclusion/preferential hydration mechanism. Although the protein is stable in the DES, structural fluctuations are largely suppressed and only recovered upon sufficient hydration. This is explained by a solvent-slaving mechanism where β-fluctuations are key, with the non-monotonic hydration of some amino acids with the water content providing an explanation to the non-monotonic folding of some proteins in aqueous DESs. A major thermal stability enhancement in the DES is also observed, caused by a similar slowdown of the backbone torsional dynamics. Our results support a kinetic stabilization of the protein in the DES, whereas a possible thermodynamic stabilization does not follow a preferential hydration or water entrapment mechanism.

Temperature-dependent dielectric relaxation and hydrophobicity of aqueous alanine using time domain reflectometry

Physical, chemical and microbiological stability of the materials is affected by the rotational and translational mobility of free and hydrated water. The role of water in areas such as protein hydration and enzyme activity, food technology, lyophilization and polymers hydration is, therefore, important and can be well understood in terms of dielectric relaxation spectroscopy. Concentration and temperature-dependent hydrophobicity of amino acid is reflected in their tendencies to appear in appropriate positions in proteins. Therefore, to gain more insights on the temperature and concentration dependence of hydrophobicity and structural properties of amino acid, dielectric relaxation of aqueous alanine have been studied in the temperature region 303.15 K to 278.15 K. Time domain spectroscopy have been used in the frequency range of 10 MHz to 30 GHz and in the concentration range 0.18708 ≤ c/M ≤ 0.74831. Two relaxation processes namely the low-frequency relaxation (l) and the high-frequency relaxation (h) has been detected for the aqueous alanine. Dielectric parameters such as static dielectric constant (εj), relaxation time (τj) dipole moments (û) and correlation factor (g) have been studied to investigate molecular interaction between alanine and water. The number of water molecules irrotationally bond by the solute molecules (Zib) was also determined to examine the hydrophobicity of alanine which was found more hydrophobic towards low temperatures and concentrations. Thermodynamic parameters calculated are also supported well for the hydrophobic behaviour of alanine towards low temperatures and concentrations.Communicated by Ramaswamy H. Sarma.

Temperature-dependent hydration behavior of aqueous lysine: an approach towards protein binding through dielectric spectroscopy

Present work reports interaction between water and amino acid lysine for understanding the physicochemical properties that will be useful in the structure formation of protein. The dielectric relaxation of aqueous lysine was systematically investigated over a temperature range spanning from 298.15 K to 278.15 K, encompassing frequencies ranging from 10 MHz to 30 GHz, and across a concentration range of 0.152 M to 0.610 M. Within this study, aqueous lysine revealed the presence of two distinct relaxation modes. The low-frequency relaxation process (l-process) is primarily associated with the relaxation of lysine molecules, whereas the high-frequency relaxation process (h-process) is attributed to water molecules interacting with lysine. Several key dielectric parameters, including static dielectric constant (εj), relaxation time (τj), dipole moment (μj), correlation factor (gj), and the number of water molecules rotationally bonded by solute molecules (Zib), were meticulously determined. These parameters were interpreted in terms of molecular interactions, hydrogen bonding, hydrophobicity, and Lys-Lys binding. Additionally, various thermodynamic parameters such as molar enthalpy (ΔHj), molar entropy (ΔSj), and molar free energy (ΔFj) were calculated to provide further insights into the system's characteristics and behavior.Communicated by Ramaswamy H. Sarma.

Oligomeric and Fibrillar α-Synuclein Display Persistent Dynamics and Compressibility under Controlled Confinement

The roles of α-synuclein in neurotransmitter release in brain neurons and in the Parkinson's disease condition have challenged comprehensive description. To gain insight into molecular mechanistic properties that actuate α-synuclein function and dysfunction, the coupled protein and solvent dynamics of oligomer and fibril forms of human α-synuclein are examined in a low-temperature system that allows control of confinement and localization of a motionally sensitive electron paramagnetic resonance spin probe in the coupled solvent-protein regions. The rotational mobility of the spin probe resolves two distinct α-synuclein-associated solvent components for oligomers and fibrils, as for globular proteins, but with dramatically higher fluidities at each temperature, that are comparable to low-confinement, aqueous-cryosolvent mesophases. In contrast to the temperature-independent volumes of the solvent phases that surround globular and condensate-forming proteins, the higher-fluidity mesophase volume of α-synuclein oligomers and fibrils decreases with decreasing temperature, signaling a compression of this phase. This unique property and thermal hysteresis in the mobilities and component weights, together with previous high-resolution structural characterizations, suggest a model in which the dynamically disordered C-terminal domain of α-synuclein creates a compressible phase that maintains high fluidity under confinement. Robust dynamics and compressibility are fundamental molecular mechanical properties of α-synuclein oligomers and fibrils, which may contribute to dysfunction and inform about function.

Membrane Mimetic-Dependence of GPCR Energy Landscapes

Protein function strongly depends on temperature, which is related to temperature-dependent changes in the equilibria of protein conformational states. We leveraged variable-temperature 19F-NMR spectroscopy to interrogate the temperature dependence of the conformational landscape of the human A2A adenosine receptor (A2AAR), a class A GPCR. Temperature-induced changes in the conformational equilibria of A2AAR in lipid nanodiscs were markedly dependent on the efficacy of bound drugs. While antagonist complexes displayed only modest changes as the temperature rose, both full and partial agonist complexes exhibited substantial increases in the active state population. Importantly, the temperature-dependent response of complexes with both full and partial agonists exhibited a pronounced sensitivity to the specific membrane mimetic employed. In striking contrast to observations within lipid nanodiscs, in detergent micelles the active state population exhibited different behavior for A2AAR complexes with both full and partial agonists. This underscores the importance of the protein environment in understanding the thermodynamics of GPCR activation.

Native and nonnative reactions in ethanolamine ammonia-lyase are actuated by different dynamics

We address the contribution of select classes of solvent-coupled configurational fluctuations to the complex choreography involved in configurational and chemical steps in an enzyme by comparing native and nonnative reactions conducted at different protein internal sites. The low temperature, first-order kinetics of covalent bond rearrangement of the cryotrapped substrate radical in coenzyme B12-dependent ethanolamine ammonia-lyase (EAL) from Salmonella enterica display a kink, or increase in slope, of the Arrhenius plot with decreasing temperature. The event is associated with quenching of a select class of reaction-actuating collective fluctuations in the protein hydration layer. For comparison, a nonnative, radical reaction of the protein interior cysteine sulfhydryl group with hydrogen peroxide (H2O2) is introduced by cryotrapping EAL in an aqueous H2O2 eutectic system. The low-temperature aqueous H2O2 protein hydration and mesodomain solvent phases surrounding cryotrapped EAL are characterized by using TEMPOL spin probe electron paramagnetic resonance spectroscopy, including a freezing transition of the eutectic phase that orders the protein hydration layer. Kinetics of the cysteine-H2O2 reaction in the EAL protein interior are monitored by DEPMPO spin trapping of hydroxyl radical product. In contrast to the native reaction, the linear Arrhenius relation for the nonnative cysteine-H2O2 reaction is maintained through the solvent-protein ordering transition. The nonnative reaction is coupled to the generic local, incremental fluctuations that are intrinsic to globular proteins. The comparative approach supports the proposal that select coupled solvent-protein configurational fluctuations actuate the native reaction, and suggests that select dynamical coupling contributes to the degree of catalysis in enzymes.

Effect of solvent viscosity on the activation barrier of hydrogen tunneling in the lipoxygenase reaction

Hydrogen tunneling in enzyme reactions has played an important role in linking protein thermal motions to the chemical steps of catalysis. Lipoxygenases (LOXs) have served as model systems for such reactions, showcasing deep hydrogen tunneling mechanisms associated with enzymatic C-H bond cleavage from polyunsaturated fatty acids. Here, we examined the effect of solvent viscosity on the protein thermal motions associated with LOX catalysis using trehalose and glucose as viscogens. Kinetic analysis of the reaction of the paradigm plant orthologue, soybean lipoxygenase (SLO), with linoleic acid revealed no effect on the first-order rate constants, kcat, or activation energy, Ea. Further studies of SLO active site mutants displaying varying Eas, which have been used to probe catalytically relevant motions, likewise provided no evidence for viscogen-dependent motions. Kinetic analyses were extended to a representative fungal LOX from M. oryzae, MoLOX, and a human LOX, 15-LOX-2. While MoLOX behaved similarly to SLO, we show that viscogens inhibit 15-LOX-2 activity. The latter implicates viscogen sensitive, conformational motions in animal LOX reactions. The data provide insight into the role of water hydration layers in facilitating hydrogen (quantum) tunneling in LOX.

Local and Large-Scale Conformational Dynamics in Unfolded Proteins and IDPs. I. Effect of Solvent Viscosity and Macromolecular Crowding

Protein/solvent interactions largely influence protein dynamics, particularly motions in unfolded and intrinsically disordered proteins (IDPs). Here, we apply triplet-triplet energy transfer (TTET) to investigate the coupling of internal protein motions to solvent motions by determining the effect of solvent viscosity (η) and macromolecular crowding on the rate constants of loop formation (kc) in several unfolded polypeptide chains including IDPs. The results show that the viscosity dependence of loop formation depends on amino acid sequence, loop length, and co-solute size. Below a critical size (rc), co-solutes exert a maximum effect, indicating that under these conditions microviscosity experienced by chain motions matches macroviscosity of the solvent. rc depends on chain stiffness and reflects the length scale of the chain motions, i.e., it is related to the persistence length. Above rc, the effect of solvent viscosity decreases with increasing co-solute size. For co-solutes typically used to mimic cellular environments, a scaling of kc ∝ η-0.1 is observed, suggesting that dynamics in unfolded proteins are only marginally modulated in cells. The effect of solvent viscosity on kc in the small co-solute limit (below rc) increases with increasing chain length and chain flexibility. Formation of long and very flexible loops exhibits a kc ∝ η-1 viscosity dependence, indicating full solvent coupling. Shorter and less flexible loops show weaker solvent coupling with values as low as kc ∝ η-0.75 ± 0.02. Coupling of formation of short loops to solvent motions is very little affected by amino acid sequence, but solvent coupling of long-range loop formation is decreased by side chain sterics.

Long-Range Conductivity in Proteins Mediated by Aromatic Residues

Single-molecule measurements show that many proteins, lacking any redox cofactors, nonetheless exhibit electrical conductance on the order of a nanosiemen over 10 nm distances, implying that electrons can transit an entire protein in less than a nanosecond when subject to a potential difference of less than 1 V. This is puzzling because, for fast transport (i.e., a free energy barrier of zero), the hopping rate is determined by the reorganization energy of approximately 0.8 eV, and this sets the time scale of a single hop to at least 1 μs. Furthermore, the Fermi energies of typical metal electrodes are far removed from the energies required for sequential oxidation and reduction of the aromatic residues of the protein, which should further reduce the hopping current. Here, we combine all-atom molecular dynamics (MD) simulations of non-redox-active proteins (consensus tetratricopeptide repeats) with an electron transfer theory to demonstrate a molecular mechanism that can account for the unexpectedly fast electron transport. According to our MD simulations, the reorganization energy produced by the energy shift on charging (the Stokes shift) is close to the conventional value of 0.8 eV. However, the non-ergodic sampling of molecular configurations by the protein results in reaction-reorganization energies, extracted directly from the distribution of the electrostatic energy fluctuations, that are only ∼0.2 eV, which is small enough to enable long-range conductivity, without invoking quantum coherent transport. Using the MD values of the reorganization energies, we calculate a current decay with distance that is in agreement with experiment.

Reactivity Tracking of an Enzyme Progress Coordinate

The reactivity of individual solvent-coupled protein configurations is used to track and resolve the progress coordinate for the core reaction sequence of substrate radical rearrangement and hydrogen atom transfer in the ethanolamine ammonia-lyase (EAL) enzyme from Salmonella enterica. The first-order decay of the substrate radical intermediate is the monitored reaction. Heterogeneous confinement from sucrose hydrates in the mesophase solvent surrounding the cryotrapped protein introduces distributed kinetics in the non-native decay of the substrate radical pair capture substate, which arise from an ensemble of configurational microstates. Reaction rates increase by >103-fold across the distribution to approach that for the native enabled substate for radical rearrangement, which reacts with monotonic kinetics. The native progress coordinate thus involves a collapse of the configuration space to generate optimized reactivity. Reactivity tracking reveals fundamental features of solvent-protein-reaction configurational coupling and leads to a model that refines the ensemble paradigm of enzyme catalysis for strongly adiabatic chemical steps.

New insights into the protein stabilizing effects of trehalose by comparing with sucrose

Disaccharides are well known to be efficient stabilizers of proteins, for example in the case of lyophilization or cryopreservation. However, although all disaccharides seem to exhibit bioprotective and stabilizing properties, it is clear that trehalose is generally superior compared to other disaccharides. The aim of this study was to understand this by comparing how the structural and dynamical properties of aqueous trehalose and sucrose solutions influence the protein myoglobin (Mb). The structural studies were based on neutron and X-ray diffraction in combination with empirical potential structure refinement (EPSR) modeling, whereas the dynamical studies were based on quasielastic neutron scattering (QENS) and molecular dynamics (MD) simulations. The results show that the overall differences in the structure and dynamics of the two systems are small, but nevertheless there are some important differences which may explain the superior stabilizing effects of trehalose. It was found that in both systems the protein is preferentially hydrated by water, but that this effect is more pronounced for trehalose, i.e. trehalose forms less hydrogen bonds to the protein surface than sucrose. Furthermore, the rotational motion around dihedrals between the two glucose rings of trehalose is slower than in the case of the dihedrals between the glucose and fructose rings of sucrose. This leads to a less perturbed protein structure in the case of trehalose. The observations indicate that an aqueous environment closest to the protein molecules is beneficial for an efficient bioprotective solution.

Combining temperature perturbations with X-ray crystallography to study dynamic macromolecules: A thorough discussion of experimental methods

Temperature is an important state variable that governs the behavior of microscopic systems, yet crystallographers rarely exploit temperature changes to study the structure and dynamics of biological macromolecules. In fact, approximately 90% of crystal structures in the Protein Data Bank were determined under cryogenic conditions, because sample cryocooling makes crystals robust to X-ray radiation damage and facilitates data collection. On the other hand, cryocooling can introduce artifacts into macromolecular structures, and can suppress conformational dynamics that are critical for function. Fortunately, recent advances in X-ray detector technology, X-ray sources, and computational data processing algorithms make non-cryogenic X-ray crystallography easier and more broadly applicable than ever before. Without the reliance on cryocooling, high-resolution crystallography can be combined with various temperature perturbations to gain deep insight into the conformational landscapes of macromolecules. This Chapter reviews the historical reasons for the prevalence of cryocooling in macromolecular crystallography, and discusses its potential drawbacks. Next, the Chapter summarizes technological developments and methodologies that facilitate non-cryogenic crystallography experiments. Finally, the chapter discusses the theoretical underpinnings and practical aspects of multi-temperature and temperature-jump crystallography experiments, which are powerful tools for understanding the relationship between the structure, dynamics, and function of proteins and other biological macromolecules.

Dioxygen and glucose force motion of the electron-transfer switch in the iron(III) flavohemoglobin-type nitric oxide dioxygenase

Kinetic and structural investigations of the flavohemoglobin-type NO dioxygenase have suggested critical roles for transient Fe(III)O2 complex formation and O2-forced movements affecting hydride transfer to the FAD cofactor and electron-transfer to the Fe(III)O2 complex. Stark-effect theory together with structural models and dipole and internal electrostatic field determinations provided a semi-quantitative spectroscopic method for investigating the proposed Fe(III)O2 complex and O2-forced movements. Deoxygenation of the enzyme causes Stark effects on the ferric heme Soret and charge-transfer bands revealing the Fe(III)O2 complex. Deoxygenation also elicits Stark effects on the FAD that expose forces and motions that create a more restricted NADH access to FAD for hydride transfer and switch electron-transfer off. Glucose also forces the enzyme toward an off state. Amino acid substitutions at the B10, E7, E11, G8, D5, and F7 positions influence the Stark effects of O2 on resting heme spin states and FAD consistent with the proposed roles of the side chains in the enzyme mechanism. Deoxygenation of ferric myoglobin and hemoglobin A also induces Stark effects on the hemes suggesting a common 'oxy-met' state. The ferric myoglobin and hemoglobin heme spectra are also glucose-responsive. A conserved glucose or glucose-6-phosphate binding site is found bridging the BC-corner and G-helix in flavohemoglobin and myoglobin suggesting novel allosteric effector roles for glucose or glucose-6-phosphate in the NO dioxygenase and O2 storage functions. The results support the proposed roles of a ferric O2 intermediate and protein motions in regulating electron-transfer during NO dioxygenase turnover.

Nanocrystallites Modulate Intermolecular Interactions in Cryoprotected Protein Solutions

Studying protein interactions at low temperatures has important implications for optimizing cryostorage processes of biological tissue, food, and protein-based drugs. One of the major issues is related to the formation of ice nanocrystals, which can occur even in the presence of cryoprotectants and can lead to protein denaturation. The presence of ice nanocrystals in protein solutions poses several challenges since, contrary to microscopic ice crystals, they can be difficult to resolve and can complicate the interpretation of experimental data. Here, using a combination of small- and wide-angle X-ray scattering (SAXS and WAXS), we investigate the structural evolution of concentrated lysozyme solutions in a cryoprotected glycerol-water mixture from room temperature (T = 300 K) down to cryogenic temperatures (T = 195 K). Upon cooling, we observe a transition near the melting temperature of the solution (T ≈ 245 K), which manifests both in the temperature dependence of the scattering intensity peak position reflecting protein-protein length scales (SAXS) and the interatomic distances within the solvent (WAXS). Upon thermal cycling, a hysteresis is observed in the scattering intensity, which is attributed to the formation of nanocrystallites in the order of 10 nm. The experimental data are well described by the two-Yukawa model, which indicates temperature-dependent changes in the short-range attraction of the protein-protein interaction potential. Our results demonstrate that the nanocrystal growth yields effectively stronger protein-protein attraction and influences the protein pair distribution function beyond the first coordination shell.

Effects of Hydrogen/Deuterium Exchange on Protein Stability in Solution and in the Gas Phase

Mass spectrometry (MS)-based techniques are widely used for probing protein structure and dynamics in solution. H/D exchange (HDX)-MS is one of the most common approaches in this context. HDX is often considered to be a "benign" labeling method, in that it does not perturb protein behavior in solution. However, several studies have reported that D₂O pushes unfolding equilibria toward the native state. The origin, and even the existence of this protein stabilization remain controversial. Here we conducted thermal unfolding assays in solution to confirm that deuterated proteins in D₂O are more stable, with 2-4 K higher melting temperatures than unlabeled proteins in H₂O. Previous studies tentatively attributed this phenomenon to strengthened H-bonds after deuteration, an effect that may arise from the lower zero-point vibrational energy of the deuterated species. Specifically, it was proposed that strengthened water-water bonds (W···W) in D₂O lower the solubility of nonpolar side chains. The current work takes a broader view by noting that protein stability in solution also depends on water-protein (W···P) and protein-protein (P···P) H-bonds. To help unravel these contributions, we performed collision-induced unfolding (CIU) experiments on gaseous proteins generated by native electrospray ionization. CIU profiles of deuterated and unlabeled proteins were indistinguishable, implying that P···P contacts are insensitive to deuteration. Thus, protein stabilization in D₂O is attributable to solvent effects, rather than alterations of intraprotein H-bonds. Strengthening of W···W contacts represents one possible explanation, but the stabilizing effect of D₂O can also originate from weakened W···P bonds. Future work will be required to elucidate which of these two scenarios is correct, or if both contribute to protein stabilization in D₂O. In any case, the often-repeated adage that "D-bonds are more stable than H-bonds" does not apply to intramolecular contacts in native proteins.

Relaxation dynamics measure the aggregation propensity of amyloid-β and its mutants

Atomistic molecular dynamics simulations are employed to investigate the global and segmental relaxation dynamics of the amyloid-β protein and its causative and protective mutants. Amyloid-β exhibits significant global/local dynamics that span a broad range of length and time scales due to its intrinsically disordered nature. The relaxation dynamics of the amyloid-β protein and its mutants is quantitatively correlated with its experimentally measured aggregation propensity. The protective mutant has slower relaxation dynamics, whereas the causative mutants exhibit faster global dynamics compared with that of the wild-type amyloid-β. The local dynamics of the amyloid-β protein or its mutants is governed by a complex interplay of the charge, hydrophobicity, and change in the molecular mass of the mutated residue.

Bimodal 1/f Noise and Anticorrelation between DNA-Water and DNA-Ion Energy Fluctuations

The coupling between the conformational fluctuations of DNA and its surrounding environment, consisting of water and ions in solution, remains poorly understood and relatively less investigated as compared to proteins. Here, with the help of molecular dynamics simulations and statistical mechanical analyses, we explore the dynamical coupling among DNA, water, and counterions through correlations among respective energy fluctuations in both double- (ds-) and single-stranded (ss-) DNA solutions. Fluctuations in the collective DNA-water and DNA-ion interaction energies are found to be strongly anticorrelated across all the systems. The fluctuations of DNA self-energy, however, are weakly coupled to DNA-water and DNA-ion interactions in ds-DNA. An enhancement of the DNA-water coupling is observed in ss-DNA, where the system is less rigid. All the interaction energies exhibit 1/f noise in their energy power spectra with surprisingly prominent bimodality in the DNA-water and DNA-ion fluctuations. The nature of the energy spectra appears to be indifferent to the relative rigidity of the DNA. We discuss the role of the observed correlations in ion-water motions on a DNA duplex in the experimentally observed anomalous slow dielectric relaxation and solvation dynamics and in furthering our understanding of the DNA energy landscape.

Temporal and spatial resolution of distal protein motions that activate hydrogen tunneling in soybean lipoxygenase

The enzyme soybean lipoxygenase (SLO) provides a prototype for deep tunneling mechanisms in hydrogen transfer catalysis. This work combines room temperature X-ray studies with extended hydrogen-deuterium exchange experiments to define a catalytically-linked, radiating cone of aliphatic side chains that connects an active site iron center of SLO to the protein-solvent interface. Employing eight variants of SLO that have been appended with a fluorescent probe at the identified surface loop, nanosecond fluorescence Stokes shifts have been measured. We report a remarkable identity of the energies of activation (Ea) for the Stokes shifts decay rates and the millisecond C-H bond cleavage step that is restricted to side chain mutants within an identified thermal network. These findings implicate a direct coupling of distal protein motions surrounding the exposed fluorescent probe to active site motions controlling catalysis. While the role of dynamics in enzyme function has been predominantly attributed to a distributed protein conformational landscape, the presented data implicate a thermally initiated, cooperative protein reorganization that occurs on a timescale faster than nanosecond and represents the enthalpic barrier to the reaction of SLO.

Light-Induced Conformational Flexibility of the Orange Carotenoid Protein Studied by Quasielastic Neutron Scattering with In Situ Illumination

The orange carotenoid protein plays a vital role in the photoprotection of cyanobacteria and exhibits a significant structural change upon photoactivation. A rarely considered aspect is the importance of internal protein dynamics in facilitating the structural transition to the active state. In this study, we use quasielastic neutron scattering under (in situ) blue light illumination for the first time to directly probe the protein dynamics of the orange carotenoid protein in the dark-adapted and active states. This shows that the localized internal dynamics of amino acid residues is significantly enhanced upon photoactivation. This is attributed to the photoinduced structural changes exposing larger areas of the protein surface to the solvent, thus resulting in a higher degree of motional freedom. However, the flexibility of the W288A mutant assumed to mimic the active state structure is found to be different, thus highlighting the importance of in situ experiments.

From structure and dynamics to biomolecular functions: The ubiquitous role of solvent in biology

Biological activity requires a solvent that can provide a suitable environment, which satisfies the twin need for stability and the ability to change. Among all the solvents water plays the most important role. We review, analyze, and comment on recent works on the structure and dynamics of water around biomolecules and their role in specific biological functions. While studies in the past have focused on understanding the biomolecule-water interactions through a hydration layer; recently the attention has shifted towards understanding functions at a molecular level. Such a microscopic understanding clearly requires elucidation of detailed dynamical processes where solvent molecules play an important role. Finally, we comment on the advances made in understanding the role of water inside a biological cell.

Thermal tuning of protein hydration in a hyperthermophilic enzyme

Water at the protein surface is an active biological molecule that plays a critical role in many functional processes. Using NMR-restrained MD simulations, we here addressed how protein hydration is tuned at high biological temperatures by analysing homologous acylphosphatase enzymes (AcP) possessing similar structure and dynamics under very different thermal conditions. We found that the hyperthermophilic Sso AcP at 80°C interacts with a lower number of structured waters in the first hydration shell than its human homologous mt AcP at 37°C. Overall, the structural and dynamical properties of waters at the surface of the two enzymes resulted similar in the first hydration shell, including solvent molecules residing in the active site. By contrast the dynamical content of water molecules in the second hydration shell was found to diverge, with higher mobility observed in Sso AcP at 80°C. Taken together the results delineate the subtle differences in the hydration properties of mt AcP and Sso AcP, and indicate that the concept of corresponding states with equivalent dynamics in homologous mesophilic and hyperthermophylic proteins should be extended to the first hydration shell.

Conformational Heterogeneity and Frustration of the Tumor Suppressor p53 as Tuned by Punctual Mutations

The conformational heterogeneity of the p53 tumor suppressor, the wild-type (p53wt) and mutated forms, was investigated by a computational approach, including the modeling and all atoms of the molecular dynamics (MD) simulations. Four different punctual mutations (p53R175H, p53R248Q, p53R273H, and p53R282W) which are known to affect the DNA binding and belong to the most frequent hot-spot mutations in human cancers, were taken into consideration. The MD trajectories of the wild-type and mutated p53 forms were analyzed by essential dynamics to extract the relevant collective motions and by the frustration method to evaluate the degeneracy of the energy landscape. We found that p53 is characterized by wide collective motions and its energy landscape exhibits a rather high frustration level, especially in the regions involved in the binding to physiological ligands. Punctual mutations give rise to a modulation of both the collective motions and the frustration of p53, with different effects depending on the mutation. The regions of p53wt and of the mutated forms characterized by a high frustration level are also largely involved in the collective motions. Such a correlation is discussed also in connection with the intrinsic disordered character of p53 and with its central functional role.

Confinement dependence of protein-associated solvent dynamics around different classes of proteins, from the EPR spin probe perspective

Protein function is modulated by coupled solvent fluctuations, subject to the degree of confinement from the surroundings. To identify universal features of the external confinement effect, the temperature dependence of the dynamics of protein-associated solvent over 200-265 K for proteins representative of different classes and sizes is characterized by using the rotational correlation time (detection bandwidth, 10-10-10-7 s) of the electron paramagnetic resonance (EPR, X-band) spin probe, TEMPOL, which is restricted to regions vicinal to protein in frozen aqueous solution. Weak (protein surrounded by aqueous-dimethylsulfoxide cryosolvent mesodomain) and strong (no added crysolvent) conditions of ice boundary confinement are imposed. The panel of soluble proteins represents large and small oligomeric (ethanolamine ammonia-lyase, 488 kDa; streptavidin, 52.8 kDa) and monomeric (myoglobin, 16.7 kDa) globular proteins, an intrinsically disordered protein (IDP, β-casein, 24.0 kDa), an unstructured peptide (protamine, 4.38 kDa) and a small peptide with partial backbone order (amyloid-β residues 1-16, 1.96 kDa). Expanded and condensate structures of β-casein and protamine are resolved by the spin probe under weak and strong confinement, respectively. At each confinement condition, the soluble globular proteins display common T-dependences of rotational correlation times and normalized weights, for two mobility components, protein-associated domain, PAD, and surrounding mesodomain. Strong confinement induces a detectable PAD component and emulation of globular protein T-dependence by the amyloid-β peptide. Confinement uniformly impacts soluble globular protein PAD dynamics, and is therefore a generic control parameter for modulation of soluble globular protein function.

A protocol to measure slow protein dynamics of the cholera toxin B pentamers using broadband dielectric spectroscopy

The present protocol describes how to measure experimentally the slow protein dynamics that take place upon the thermal unfolding of the B subunit cholera toxin pentamers using broadband dielectric spectroscopy (BDS) in weakly hydrated and nanoconfined conditions. Transient unfolding intermediates, rarely identified otherwise, are revealed thanks to the B subunit's remarkable heat resistance up to 180°C and distinct molecular dynamics. The frequencies detected experimentally are consistent with the spatiotemporal scales of motions of molecular dynamics simulation.

For complete details on the use and execution of this protocol, please refer to Bourgeat et al. (2021, 2019).

•

Measure protein dynamics experimentally using BDS in nanoconfined conditions

•

Identify rare cholera toxin B subunit assembly and unfolding intermediates

•

Detect cholera toxin B subunits in temperatures up to 180°C

•

Match between protein molecular dynamics from experiments and simulations

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Structure and dynamics of supercooled water in the hydration layer of poly(ethylene glycol)

The statics and dynamics of supercooled water in the hydration layer of poly(ethylene glycol) (PEG) were studied by a combination of quasi-elastic neutron scattering (QENS) and molecular dynamics (MD) simulations. Two samples, that is, hydrogenated PEG/deuterated water (h-PEG/D₂O) and fully deuterated PEG/hydrogenated water (d-PEG/H₂O) with the same molar ratio of ethylene glycol (EG) monomer to water, 1:1, are compared. The QENS data of h-PEG/D₂O show the dynamics of PEG, and that of d-PEG/H₂O reveals the motion of water. The temperature-dependent elastic scattering intensity of both samples has shown transitions at supercooled temperature, and these transition temperatures depend on the energy resolution of the instruments. Therefore, neither one is a phase transition, but undergoes dynamic process. The dynamic of water can be described as an Arrhenius to super-Arrhenius transition, and it reveals the hydrogen bonding network relaxation of hydration water around PEG at supercooled temperature. Since the PEG-water hydrogen bond structural relaxation time from MD is in good agreement with the average relaxation time from QENS (d-PEG/H₂O), MD may further reveal the atomic pictures of the supercooled hydration water. It shows that hydration water molecules form a series of pools around the hydrophilic oxygen atom of PEG. At supercooled temperature, they have a more bond ordered structure than bulk water, proceed a trapping sites diffusion on the PEG surface, and facilitate the structural relaxation of PEG backbone.

Exploring the Limits of Biological Complexity Amenable to Studies by Incoherent Neutron Spectroscopy

The wavelengths of neutrons available at neutron scattering facilities are comparable with intra- and inter-molecular distances, while their energies are comparable with molecular vibrational energies, making such neutrons highly suitable for studies of molecular-level dynamics. The unmistakable trend in neutron spectroscopy has been towards measurements of systems of greater complexity. Several decades of studies of dynamics using neutron scattering have witnessed a progression from measurements of solids to liquids to protein complexes and biomembranes, which may exhibit properties characteristic of both solids and liquids. Over the last two decades, the frontier of complexity amenable to neutron spectroscopy studies has reached the level of cells. Considering this a baseline for neutron spectroscopy of systems of the utmost biological complexity, we briefly review what has been learned to date from neutron scattering studies at the cellular level and then discuss in more detail the recent strides into neutron spectroscopy of tissues and whole multicellular organisms.

Water Networks Repopulate Protein-Ligand Interfaces with Temperature

High-resolution crystal structures highlight the importance of water networks in protein-ligand interactions. However, as these are typically determined at cryogenic temperature, resulting insights may be structurally precise but not biologically accurate. By collecting 10 matched room-temperature and cryogenic datasets of the biomedical target Hsp90α, we identified changes in water networks that impact protein conformations at the ligand binding interface. Water repositioning with temperature repopulates protein ensembles and ligand interactions. We introduce Flipper conformational barcodes to identify temperature-sensitive regions in electron density maps. This revealed that temperature-responsive states coincide with ligand-responsive regions and capture unique binding signatures that disappear upon cryo-cooling. Our results have implications for discovering Hsp90 selective ligands, and, more generally, for the utility of hidden protein and water conformations in drug discovery.

Temperature-Dependent Dynamics at Protein-Solvent Interfaces

We perform differential scanning calorimetry, broadband dielectric spectroscopy (BDS), and nuclear magnetic resonance (NMR) studies to ascertain the molecular dynamics in mixtures of ethylene glycol with elastin or lysozyme over broad temperature ranges. To focus on the protein-solvent interface, we use mixtures with about equal numbers of amino acids and solvent molecules. The elastin and lysozyme mixtures show similar glass transition steps, which extend over a broad temperature range of 157-185K. The BDS and NMR studies yield fully consistent results for the fastest process P1, which is caused by the structural relaxation of ethylene glycol between the protein molecules and follows an Arrhenius law with an activation energy of Ea=0.63eV. It involves quasi-isotropic reorientation and is very similar in the elastin and lysozyme matrices but different from the alpha and beta relaxations of bulk ethylene glycol. Two slower BDS processes P2 and P3 have protein-dependent time scales, but exhibit a similar Arrhenius-like temperature dependence with an activation energy of Ea~0.81eV. However, P2 and P3 do not have a clear NMR signature. In particular, the NMR results for the lysozyme mixture reveal that the protein backbone does not show isotropic alpha-like motion on the P2 and P3 time scales but only restricted beta-like reorientation. The different activation energies of the P1 and P2/P3 processes do not support an intimate coupling of protein and ethylene glycol dynamics. The present results are compared with previous findings for mixtures of proteins with water or glycerol, implying qualitatively different dynamical couplings at various protein-solvent interfaces.

Probing Adaptation of Hydration and Protein Dynamics to Temperature

Protein dynamics is strongly influenced by the surrounding environment and physiological conditions. Here we employ broadband megahertz-to-terahertz spectroscopy to explore the dynamics of water and myoglobin protein on an extended time scale from femto- to nanosecond. The dielectric spectra reveal several relaxations corresponding to the orientational polarization mechanism, including the dynamics of loosely bound, tightly bound, and bulk water, as well as collective vibrational modes of protein in an aqueous environment. The dynamics of loosely bound and bulk water follow non-Arrhenius behavior; however, the dynamics of water molecules in the tightly bound layer obeys the Arrhenius-type relation. Combining molecular simulations and effective-medium approximation, we have determined the number of water molecules in the tightly bound hydration layer and studied the dynamics of protein as a function of temperature. The results provide the important impact of water on the biochemical functions of proteins.

Resolution and characterization of contributions of select protein and coupled solvent configurational fluctuations to radical rearrangement catalysis in coenzyme B12-dependent ethanolamine ammonia-lyase

Coenzyme B₁₂ (adenosylcobalamin) -dependent ethanolamine ammonia-lyase (EAL) is the signature enzyme in ethanolamine utilization metabolism associated with microbiome homeostasis and disease conditions in the human gut. The enzyme conducts a complex choreography of bond-making/bond-breaking steps that rearrange substrate to products through a radical mechanism, with themes common to other coenzyme B₁₂-dependent and radical enzymes. The methods presented are targeted to test the hypothesis that particular, select protein and coupled solvent configurational fluctuations contribute to enzyme function. The general approach is to correlate enzyme function with an introduced perturbation that alters the properties (for example, degree of concertedness, or collectiveness) of protein and coupled solvent dynamics. Methods for sample preparation and low-temperature kinetic measurements by using temperature-step reaction initiation and time-resolved, full-spectrum electron paramagnetic resonance spectroscopy are detailed. A framework for interpretation of results obtained in ensemble systems under conditions of statistical equilibrium within the reacting, globally unstable state is presented. The temperature-dependence of the first-order rate constants for decay of the cryotrapped paramagnetic substrate radical state in EAL, through the chemical step of radical rearrangement, displays a piecewise-continuous Arrhenius dependence from 203 to 295K, punctuated by a kinetic bifurcation over 219-220K. The results reveal the obligatory contribution of a class of select collective protein and coupled solvent fluctuations to the interconversion of two resolved, sequential configurational substates, on the decay time scale. The select class of collective fluctuations also contributes to the chemical step. The methods and analysis are generally applicable to other coenzyme B₁₂-dependent and related radical enzymes.

Hydration-mediated G-protein-coupled receptor activation

Although G-protein–coupled receptors (GPCRs) control vast physiological pathways, their activation remains chemically and physically enigmatic. Our osmotic stress studies of the visual receptor rhodopsin have redefined the standard model of GPCR signaling by revealing the essential role of bulk water. We show results consistent with a large number of water molecules flooding the rhodopsin interior during activation to stabilize the effector binding conformation. These results suggest a model of GPCR activation in which the receptor becomes solvent-swollen upon formation of the active state. We thus demonstrate the mechanism whereby water acts as a powerful allosteric modulator of a pharmacologically important membrane protein family.

The Rhodopsin family of G-protein–coupled receptors (GPCRs) comprises the targets of nearly a third of all pharmaceuticals. Despite structural water present in GPCR X-ray structures, the physiological relevance of these solvent molecules to rhodopsin signaling remains unknown. Here, we show experimental results consistent with the idea that rhodopsin activation in lipid membranes is coupled to bulk water movements into the protein. To quantify hydration changes, we measured reversible shifting of the metarhodopsin equilibrium due to osmotic stress using an extensive series of polyethylene glycol (PEG) osmolytes. We discovered clear evidence that light activation entails a large influx of bulk water (∼80–100 molecules) into the protein, giving insight into GPCR activation mechanisms. Various size polymer osmolytes directly control rhodopsin activation, in which large solutes are excluded from rhodopsin and dehydrate the protein, favoring the inactive state. In contrast, small osmolytes initially forward shift the activation equilibrium until a quantifiable saturation point is reached, similar to gain-of-function protein mutations. For the limit of increasing osmolyte size, a universal response of rhodopsin to osmotic stress is observed, suggesting it adopts a dynamic, hydrated sponge-like state upon photoactivation. Our results demand a rethinking of the role of water dynamics in modulating various intermediates in the GPCR energy landscape. We propose that besides bound water, an influx of bulk water plays a necessary role in establishing the active GPCR conformation that mediates signaling.

Structural, dynamic, and hydration properties of quercetin and its aggregates in solution

Quercetin is a flavonoid present in the human diet with multiple health benefits. Quercetin solutions are inhomogeneous even at very low concentrations due to quercetin's tendency to aggregate. We simulate, using molecular dynamics, three systems of quercetin solutions: infinite dilution, 0.22 M, and 0.46 M. The systems at the two highest concentrations represent regions of the quercetin aggregates, in which the concentration of this molecule is unusually high. We study the behavior of this molecule, its aggregates, and the modifications in the surrounding water. In the first three successive layers of quercetin hydration, the density of water and the hydrogen bonds formations between water molecules are smaller than that of bulk. Quercetin has a hydrophilic surface region that preferentially establishes donor hydrogen bonds with water molecules with relative frequencies from 0.12 to 0.46 at infinite dilution. Also, it has two hydrophobic regions above and below the planes of its rings, whose first hydration layers are further out from quercetin (≈0.3 Å) and their water molecules do not establish hydrogen bonds with it. Water density around the hydrophobic regions is smaller than that of the hydrophilic. Quercetin molecules aggregate inπ-stacking configurations, with a distance of ≈0.37 nm between the planes of their rings, and form bonds between their hydroxyl groups. The formation of quercetin aggregates decreases the hydrogen bonds between quercetin and the surrounding water and produces a subdiffusive behavior in water molecules. Quercetin has a subdiffusive behavior even at infinite dilution, which increases with the number of molecules within the aggregates and the time they remain within them.