Hydration dependence of myoglobin dynamics studied with elastic neutron scattering, differential scanning calorimetry and broadband dielectric spectroscopy

A Long Journey into the Investigation of the Structure–Dynamics–Function Paradigm in Proteins through the Activities of the Palermo Biophysics Group

An overview of the biophysics activity at the Department of Physics and Chemistry Emilio Segrè of the University of Palermo is given. For forty years, the focus of the research has been on the protein structure–dynamics–function paradigm, with the aim of understanding the molecular basis of the relevant mechanisms and the key role of solvent. At least three research lines are identified; the main results obtained in collaboration with other groups in Italy and abroad are presented. This review is dedicated to the memory of Professors Massimo Ugo Palma, Maria Beatrice Palma Vittorelli, and Lorenzo Cordone, which were the founders of the Palermo School of Biophysics. We all have been, directly or indirectly, their pupils; we miss their enthusiasm for scientific research, their deep physical insights, their suggestions, their strict but always constructive criticisms, and, most of all, their friendship. This paper is dedicated also to the memory of Prof. Hans Frauenfelder, whose pioneering works on nonexponential rebinding kinetics, protein substates, and energy landscape have inspired a large part of our work in the field of protein dynamics.

Unravelling the Adaptation Mechanisms to High Pressure in Proteins

Life is thought to have appeared in the depth of the sea under high hydrostatic pressure. Nowadays, it is known that the deep biosphere hosts a myriad of life forms thriving under high-pressure conditions. However, the evolutionary mechanisms leading to their adaptation are still not known. Here, we show the molecular bases of these mechanisms through a joint structural and dynamical study of two orthologous proteins. We observed that pressure adaptation involves the decoupling of protein–water dynamics and the elimination of cavities in the protein core. This is achieved by rearranging the charged residues on the protein surface and using bulkier hydrophobic residues in the core. These findings will be the starting point in the search for a complete genomic model explaining high-pressure adaptation.

Dynamics of Apolipoprotein B-100 in Interaction with Detergent Probed by Incoherent Neutron Scattering

Apolipoprotein B-100 (apo B-100) is the protein moiety of both low- and very-low-density lipoproteins, whose role is crucial to cholesterol and triglyceride transport. Aiming at the molecular dynamics' details of apo B-100, scarcely studied, we performed elastic and quasi-elastic incoherent neutron scattering (EINS, QENS) experiments combining different instruments and time scales. Similar to classical membrane proteins, the solubilization results in remaining detergent, here Nonidet P-40 (NP40). Therefore, we propose a framework for QENS studies of protein-detergent complexes, with the introduction of a combined model, including the experimental apo B-100/NP40 ratio. Relying on the simultaneous analysis of all QENS amplitudes, this approach is sensitive enough to separate both contributions. Its application identified two points: (i) apo B-100 slow dynamics and (ii) the acceleration of NP40 dynamics in the presence of apo B-100. Direct translation of the exposed methodology now makes the investigation of more membrane proteins by neutron spectroscopy achievable.

Effect of Polymer Hydrophilicity and Molar Mass on the Properties of the Protein in Protein-Polymer Conjugates: The Case of PPEylated Myoglobin

Polyphosphoesters (PPEs), a versatile class of biodegradable and biocompatible polymers, have been proposed as alternatives to poly(ethylene glycol) (PEG), which is suspected to be responsible for anaphylactic reactions in some patients after the administration of PEGylated compounds, e.g., in the current Covid-19 vaccines. We present the synthesis and characterization of a novel set of protein-polymer conjugates using the model protein myoglobin and a set of PPEs with different hydrophilicity and molar mass. We report an extensive evaluation of the (bio)physical properties of the protein within the conjugates, studying its conformation, residual activity, and thermal stability by complementary techniques (UV-vis spectroscopy, nano-differential scanning calorimetry, and fluorometry). The data underline the systematic influence of polymer hydrophilicity on protein properties. The more hydrophobic polymers destabilize the protein, the more hydrophilic PPEs protect against thermally induced aggregation and proteolytic degradation. This basic study aims at guiding the design of future PPEylated drugs and protein conjugates.

Quasielastic neutron scattering studies on couplings of protein and water dynamics in hydrated elastin

Performing quasielastic neutron scattering measurements and analyzing both elastic and quasielasic contributions, we study protein and water dynamics of hydrated elastin. At low temperatures, hydration-independent methyl group rotation dominates the findings. It is characterized by a Gaussian distribution of activation energies centered at about E_m = 0.17 eV. At ∼195 K, coupled protein-water motion sets in. The hydration water shows diffusive motion, which is described by a Gaussian distribution of activation energies with E_m = 0.57 eV. This Arrhenius behavior of water diffusion is consistent with previous results for water reorientation, but at variance with a fragile-to-strong crossover at ∼225 K. The hydration-related elastin backbone motion is localized and can be attributed to the cage rattling motion. We speculate that its onset at ∼195 K is related to a secondary glass transition, which occurs when a β relaxation of the protein has a correlation time of τ_β ∼ 100 s. Moreover, we show that its temperature-dependent amplitude has a crossover at the regular glass transition T_g = 320 K of hydrated elastin, where the α relaxation of the protein obeys τ_α ∼ 100 s. By contrast, we do not observe a protein dynamical transition when water dynamics enters the experimental time window at ∼240 K.

Insight into Protein-Polymer Conjugate Relaxation Dynamics: The Importance of Polymer Grafting

The bio and chemical physics of protein-polymer conjugates are related to parameters that characterize each component. With this work, it is intended to feature the dynamical properties of the protein-polymer conjugate myoglobin (Mb)-poly(ethyl ethylene phosphate), in the ps and ns time scales, in order to understand the respective roles of the protein and of the polymer size in the dynamics of the conjugate. Elastic and quasi-elastic neutron scattering is performed on completely hydrogenated samples with variable number of polymer chains covalently attached to the protein. The role of the polymer length in the protein solvation and internal dynamics is investigated using two conjugates formed by polymers of different molecular weight. It is confirmed that the flexibility of the complex increases with the number of grafted polymer chains and that a sharp dynamical transition appears when either grafting density or polymer molecular weight are high. It is shown that protein size is crucial for the polymer structural organization and interaction on the protein surface and it is established that the glass properties of the polymer change upon conjugation. The results give a better insight of the equivalence of the polymer coating and the role of water on the surface of proteins.

Segmental Motions of Proteins under Non-native States Evaluated Using Quasielastic Neutron Scattering

Characterization of the dynamics of disordered polypeptide chains is required to elucidate the behavior of intrinsically disordered proteins and proteins under non-native states related to the folding process. Here we develop a method using quasielastic neutron scattering, combined with small-angle X-ray scattering and dynamic light scattering, to evaluate segmental motions of proteins as well as diffusion of the entire molecules and local side-chain motions. We apply this method to RNase A under the unfolded and molten-globule (MG) states. The diffusion coefficients arising from the segmental motions are evaluated and found to be different between the unfolded and MG states. The values obtained here are consistent with those obtained using the fluorescence-based techniques. These results demonstrate not only feasibility of this method but also usefulness to characterize the behavior of proteins under various disordered states.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

The Role of Hydrogen Bonding in the Folding/Unfolding Process of Hydrated Lysozyme: A Review of Recent NMR and FTIR Results

The biological activity of proteins depends on their three-dimensional structure, known as the native state. The main force driving the correct folding mechanism is the hydrophobic effect and when this folding kinetics is altered, aggregation phenomena intervene causing the occurrence of illnesses such as Alzheimer and Parkinson’s diseases. The other important effect is performed by water molecules and by their ability to form a complex network of hydrogen bonds whose dynamics influence the mobility of protein amino acids. In this work, we review the recent results obtained by means of spectroscopic techniques, such as Fourier Transform Infrared (FTIR) and Nuclear Magnetic Resonance (NMR) spectroscopies, on hydrated lysozyme. In particular, we explore the Energy Landscape from the thermal region of configurational stability up to that of the irreversible denaturation. The importance of the coupling between the solute and the solvent will be highlighted as well as the different behaviors of hydrophilic and hydrophobic moieties of protein amino acid residues.

Direct Experimental Characterization of Contributions from Self-Motion of Hydrogen and from Interatomic Motion of Heavy Atoms to Protein Anharmonicity

One fundamental challenge in biophysics is to understand the connection between protein dynamics and its function. Part of the difficulty arises from the fact that proteins often present local atomic motions and collective dynamics on the same time scales, and challenge the experimental identification and quantification of different dynamic modes. Here, by taking lyophilized proteins as the example, we combined deuteration technique and neutron scattering to separate and characterize the self-motion of hydrogen and the collective interatomic motion of heavy atoms (C, O, N) in proteins on the pico-to-nanosecond time scales. We found that hydrogen atoms present an instrument-resolution-dependent onset for anharmonic motions, which can be ascribed to the thermal activation of local side-group motions. However, the protein heavy atoms exhibit an instrument-resolution-independent anharmonicity around 200 K, which results from unfreezing of the relaxation of the protein structures on the laboratory equilibrium time (100-1000 s), softening of the entire bio-macromolecules.

Spectroscopic methods for assessing the molecular origins of macroscopic solution properties of highly concentrated liquid protein solutions

In cases of subcutaneous injection of therapeutic monoclonal antibodies, high protein concentrations (>50 mg/ml) are often required. During the development of these high concentration liquid formulations (HCLF), challenges such as aggregation, gelation, opalescence, phase separation, and high solution viscosities are more prone compared to low concentrated protein formulations. These properties can impair manufacturing processes, as well as protein stability and shelf life. To avoid such unfavourable solution properties, a detailed understanding about the nature of these properties and their driving forces are required. However, the fundamental mechanisms that lead to macroscopic solution properties, as above mentioned, are complex and not fully understood, yet. Established analytical methods for assessing the colloidal stability, i.e. the ability of a native protein to remain dispersed in solution, are restricted to dilute conditions and provide parameters such as the second osmotic virial coefficient, B₂₂, and the diffusion interaction coefficient, k_D. These parameters are routinely applied for qualitative estimations and identifications of proteins with challenging solution behaviours, such as high viscosities and aggregation, although the assays are prepared for low protein concentration conditions, typically between 0.1 and 20 mg/ml ("ideal" solution conditions). Quantitative analysis of samples of high protein concentration is difficult and it is hard to obtain information about the driving forces of such solution properties and corresponding protein-protein self-interactions. An advantage of using specific spectroscopic methods is the potential of directly analysing highly concentrated protein solutions at different solution conditions. This allows for collecting/gaining valuable information about the fundamental mechanisms of solution properties of the high protein concentration regime. In addition, the derived parameters might be more predictive as compared to the parameters originating from assays which are optimized for the low protein concentration range. The provided information includes structural data, molecular dynamics at various timescales and protein-solvent interactions, which can be obtained at molecular resolution. Herein, we provide an overview about spectroscopic techniques for analysing the origins of macroscopic solution behaviours in general, with a specific focus on pharmaceutically relevant high protein concentration and formulation conditions.

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.

Molecular Dynamics Simulations of a Powder Model of the Intrinsically Disordered Protein Tau

The tau protein, whose aggregates are involved in Alzheimer's disease, is an intrinsically disordered protein (IDP) that regulates microtubule activity in neurons. An IDP lacks a single, well-defined structure and, rather, constantly exchanges among multiple conformations. In order to study IDP dynamics, the combination of experimental techniques, such as neutron scattering, and computational techniques, such as molecular dynamics (MD) simulations, is a powerful approach. Amorphous hydrated powder samples have been very useful for studying protein internal dynamics experimentally, e.g., using neutron scattering. Thus, there is demand for realistic in silico models of hydrated protein powders. Here we present an MD simulation analysis of a powder hydrated at 0.4 g water/g protein of the IDP tau in the temperature range 20-300 K. By comparing with neutron scattering data, we identify the protein-water interface as the predominant feature determining IDP dynamics. The so-called protein dynamical transition is shown to be attenuated, but not suppressed, in the parts of the protein that are not exposed to the solvent. In addition, we find similarities in the mean-squared displacements of the core of a globular protein and "dry" clusters formed by the IDP in hydrated powders. Thus, the ps to ns dynamics of proteins in hydrated powders originate mainly from those residues in contact with solvent. We propose that by measuring the dynamics of protein assemblies, such as aggregates, one might assess qualitatively their state of hydration.

Experimental evidence for a liquid-liquid crossover in deeply cooled confined water

In this work we investigate, by means of elastic neutron scattering, the pressure dependence of mean square displacements (MSD) of hydrogen atoms of deeply cooled water confined in the pores of a three-dimensional disordered SiO2 xerogel; experiments have been performed at 250 and 210 K from atmospheric pressure to 1200 bar. The "pressure anomaly" of supercooled water (i.e., a mean square displacement increase with increasing pressure) is observed in our sample at both temperatures; however, contrary to previous simulation results and to the experimental trend observed in bulk water, the pressure effect is smaller at lower (210 K) than at higher (250 K) temperature. Elastic neutron scattering results are complemented by differential scanning calorimetry data that put in evidence, besides the glass transition at about 170 K, a first-order-like endothermic transition occurring at about 230 K that, in view of the neutron scattering results, can be attributed to a liquid-liquid crossover. Our results give experimental evidence for the presence, in deeply cooled confined water, of a crossover occurring at about 230 K (at ambient pressure) from a liquid phase predominant at 210 K to another liquid phase predominant at 250 K; therefore, they are fully consistent with the liquid-liquid transition hypothesis.

The importance of individual protein molecule dynamics in developing and assessing solid state protein preparations

Processing protein solutions into the solid state is a common approach for generating stable amorphous protein mixtures that are suitable for long-term storage. Great care is typically given to protecting the protein native structure during the various drying steps that render it into the amorphous solid state. However, many studies illustrate that chemical and physical degradations still occur in spite of this amorphous material having good glassy properties and it being stored at temperatures below its glass transition temperature (Tg). Because of these persistent issues and recent biophysical studies that have refined the debate ascribing meaning to the molecular dynamical transition temperature and Tg of protein molecules, we provide an updated discussion on the impact of assessing and managing localized, individual protein molecule nondiffusive motions in the context of proteins being prepared into bulk amorphous mixtures. Our aim is to bridge the pharmaceutical studies addressing bulk amorphous preparations and their glassy behavior, with the biophysical studies historically focused on the nondiffusive internal protein dynamics and a protein's activity, along with their combined efforts in assessing the impact of solvent hydrogen-bonding networks on local stability. We also provide recommendations for future research efforts in solid-state formulation approaches.