Biological water at the protein surface: dynamical solvation probed directly with femtosecond resolution

Terahertz calorimetry spotlights the role of water in biological processes

Terahertz (THz) calorimetry is a framework that allows for the deduction and quantification of changes in solvation entropy and enthalpy associated with biological processes in real-time. Fundamental biological processes are inherently non-equilibrium, and a small imbalance in free energy can trigger protein condensation or folding. Although biophysical techniques typically focus mainly on structural characterization, water is often ignored. Being a generic solvent, the intermolecular protein-water interactions act as a strong competitor for intramolecular protein-protein interactions, leading to a delicate balance between functional structure formation and complete solvation. Characteristics for biological processes are large, but competing enthalpic and entropic solvation contributions to the total Gibbs free energy lead to subtle energy differences of only a few kJ mol-1 that are capable of dictating biological functions. THz calorimetry spotlights these intermolecular coupled protein-water interactions. With experimental advances in THz technology, a new frequency window has opened, which is ideally suited to probe these low-frequency intermolecular interactions. The future impact of these studies is based on the belief that the observed changes in solvation entropy and enthalpy are not secondary effects but dictate biological function.

Photophysics of a Nucleoside Analogue: 4-Cyanoindole-2'-deoxyribonucleoside

It has been shown that 4-cyanoindole-2'-deoxyribonucleoside (4CNI-NS) is a versatile spectroscopic probe of DNA structure and dynamics, as it can pair with all four natural DNA bases. However, its photophysics have not been examined in detail. Herein, we employed multiple techniques, including static fluorescence spectroscopy, time-resolved fluorescence spectroscopy, transient infrared spectroscopy, and theoretical calculations, to assess the photophysical properties of this nucleoside analogue in a series of solvents. We found that (1) its fluorescence has a large quantum yield (0.85 ± 0.10) and a relatively long decay lifetime (10.3 ± 1.0 ns) in H2O, both of which become smaller in other solvents examined, (2) its maximum fluorescence emission frequency exhibits a linear dependence on the relative polarity index of the solvent, (3) the oscillator strength of its C≡N stretching vibration is increased by a factor of >10 upon transition to its excited electronic state (S1), (4) besides solvent relaxation, a rotational motion around the single bond connecting the pentose group and the indole ring is also present in the S1 state, and (5) in aprotic solvents both processes lead the C≡N stretching frequency (νESA) to shift toward higher frequencies, whereas in protic solvents the effects of these processes are more complex. Taken together, these findings provide a molecular basis for interpreting the spectroscopic signals of this nucleoside analogue in practical applications.

Evaluating long-range orientational ordering of water around proteins: signature of a tug-of-war scenario

Long-range perturbations of water structure and dynamics by biomolecules are of great interest owing to their potential role in biomolecular recognition. In this article, we examined the local and long-range orientational structure of water molecules surrounding proteins with different total charges (+8, 0 and -8), both with and without the presence of a physiological salt environment. A prominent population of in-oriented water molecules was observed in the first hydration shell of the proteins, irrespective of their total charges. Starting from the third hydration layer, water molecules primarily reflected the total charge of the respective protein. This long-range ordering persisted up to the ninth hydration layer without a physiological salt environment and vanished beyond the fifth hydration shell in the presence of a physiological salt environment. Long-range orientational ordering around different types of surface atoms of a protein showed a particularly rich and heterogeneous behaviour. When the surface atom's charge and the protein's total charge were opposite, a clear signature of a tug-of-war was demonstrated in the long-range orientational ordering of water molecules. While water molecules reported the surface atom's charge at shorter distances, at longer distances, water molecules reported the total charge of the protein, with a crossover occurring around 10 Å. This phenomenon persisted even in the presence of a physiological salt environment. Evidence of destructive/constructive superposition of water-mediated orientation waves originating from two individual proteins with similar/opposite total charges was also demonstrated. These results are important for understanding long-range water-mediated recognition phenomena among biomolecules (e.g., protein-protein, protein-ligand, and protein-DNA interactions).

Multi-Site Red-Edge Excitation Shift Reveals the Residue-Specific Solvation Dynamics during the Native to Amyloid-like Transition of an Amyloidogenic Protein

Changes in water-protein interactions are crucial for proteins to achieve functional and nonfunctional conformations during structural transitions by modulating local stability. Amyloid-like protein aggregates in deteriorating neurons are hallmarks of neurodegenerative disorders. These aggregates form through significant structural changes, transitioning from functional native conformations to supramolecular cross-β-sheet structures via misfolded and oligomeric intermediates in a multistep process. However, the site-specific dynamics of water molecules from the native to misfolded conformations and further to oligomeric and compact amyloid structures remain poorly understood. In this study, we used the fluorescence method known as red-edge excitation shift (REES) to investigate the solvation dynamics at specific sites in various equilibrium conformations en route to the misfolding and aggregation of the functional domain of the TDP-43 protein (TDP-43tRRM). We generated three single tryptophan-single cysteine mutants of TDP-43tRRM, with the cysteines at different positions and tryptophan at a fixed position. Each sole cysteine was fluorescently labeled and used as a site-specific fluorophore along with the single tryptophan, creating four monitorable sites for REES studies. By investigating the site-specific extent of REES, we developed a residue-specific solvation dynamics map of TDP-43tRRM during its misfolding and aggregation. Our observations revealed that solvation dynamics progressively became more rigid and heterogeneous to varying extents at different sites during the transition from native to amyloid-like conformations.

Critical Role of Water beyond the Media to Maintain Protein Stability and Activity in Hydrated Deep Eutectic Solvent

Hydrated deep eutectic solvents (DESs) are recognized for their potential in biocatalysis due to their tunability, biocompatibility, greenness, and ability to keep protein stable and active. However, the mechanisms governing enzyme stability and activity in DES remain poorly understood. Herein, using bromelain as the model enzyme and acetamide (0.5)/urea(0.3)/sorbitol(0.2) as the model DES, we provide experimental evidence that modulation of associated water plays a key role in dictating protein stability and activity in hydrated DES. Specifically, rigid associated water at higher DES concentrations (beyond 40% v/v) stabilizes bromelain through entropy but destabilizes it through enthalpy. On the other hand, flexible associated water dynamics at lower DES concentrations result in an opposite thermodynamic outcome. Importantly, the bulk water dynamics cannot explain the stability trend, which emphasizes the critical role of water near the protein surface. Strikingly, associated water dynamics also correlates strongly with bromelain's proteolytic activity. An increasing flexibility of the associated water dynamics leads to the enhancement of the activity. This is the first study to experimentally link associated water dynamics to enzyme behavior in hydrated DES, offering insights that could guide future developments in solvent engineering for enzyme catalysis.

Peculiarities of the Dynamical Hydration Shell of Native Conformation Protein Using a Bovine Serum Albumin Example

This paper describes an approach based on the method of terahertz time-domain spectroscopy, which allows the analysis of dynamical hydration shells of proteins with a thickness of 1-2 nm. Using the example of bovine serum albumin in three conformations, it is shown that the hydration shells of the protein are characterized by increased binding of water molecules in the primary hydration layers, and in more distant areas of hydration, on the contrary, the water structure is somewhat destroyed. The fraction of free or weakly bound molecules, usually observed in the structure of liquid water in hydration shells, become more numerous but its average binding is greater than in undisturbed water. The energy distribution of hydrogen bonds in hydration shells is narrowed compared to undisturbed water. All these manifestations of hydration are most pronounced for the native conformation of the protein. Also, the hydration shells of the native protein are characterized by a smaller number of hydrogen bonds and a tendency to decrease their average energy compared to non-native conformations. The fact of a pronounced peculiarity of the hydration shells of the protein in the native conformation has been noted for different proteins before. However, the methodological approach used in this work for the first time allowed this peculiarity to be described by specific parameters of the intermolecular structure and dynamics of water.

Water is a radiation protection agent for ionised pyrrole

Radiation-induced damage of biological matter is an ubiquitous problem in nature. The influence of the hydration environment is widely discussed, but its exact role remains elusive. Utilising well defined solvated-molecule aggregates, we experimentally observed a hydrogen-bonded water molecule acting as a radiation protection agent for ionised pyrrole, a prototypical aromatic biomolecule. Pure samples of pyrrole and pyrrole(H2O) were outer-valence ionised and the subsequent damage and relaxation processes were studied. Bare pyrrole ions fragmented through the breaking of C-C or N-C covalent bonds. However, for pyrrole(H2O)+, we observed a strong protection of the pyrrole ring through the dissociative release of neutral water or by transferring an electron or proton across the hydrogen bond. Overall, a single water molecule strongly reduces the fragmentation probability and thus the persistent radiation damage of singly-ionised pyrrole.

Chasing weakly-bound biological water in aqueous environment near the peptide backbone by ultrafast 2D infrared spectroscopy

There has been a long-standing debate as to how many hydrogen bonds a peptide backbone amide can form in aqueous solution. Hydrogen-bonding structural dynamics of N-ethylpropionamide (a β-peptide model) in water was examined using infrared (IR) spectroscopy. Two amide-I sub bands arise mainly from amide C=O group that forms strong H-bonds with solvent water molecules (SHB state), and minorly from that involving one weak H-bond with water (WHB state). This picture is supported by molecular dynamics simulations and ab-initio calculations. Further, thermodynamics and kinetics of the SHB and WHB species were examined mainly by chemical-exchange two-dimensional IR spectroscopy, yielding an activation energy for the SHB-to-WHB exchange of 13.25 ± 0.52 kJ mol‒1, which occurs in half picosecond at room temperature. Our results provided experimental evidence of an unstable water molecule near peptide backbone, allowing us to gain more insights into the dynamics of the protein backbone hydration.

Directed ultrafast conformational changes accompany electron transfer in a photolyase as resolved by serial crystallography

Charge-transfer reactions in proteins are important for life, such as in photolyases which repair DNA, but the role of structural dynamics remains unclear. Here, using femtosecond X-ray crystallography, we report the structural changes that take place while electrons transfer along a chain of four conserved tryptophans in the Drosophila melanogaster (6-4) photolyase. At femto- and picosecond delays, photoreduction of the flavin by the first tryptophan causes directed structural responses at a key asparagine, at a conserved salt bridge, and by rearrangements of nearby water molecules. We detect charge-induced structural changes close to the second tryptophan from 1 ps to 20 ps, identifying a nearby methionine as an active participant in the redox chain, and from 20 ps around the fourth tryptophan. The photolyase undergoes highly directed and carefully timed adaptations of its structure. This questions the validity of the linear solvent response approximation in Marcus theory and indicates that evolution has optimized fast protein fluctuations for optimal charge transfer.

Tunable Surface Wettability via Terahertz Electrowave Controlled Vicinal Subnanoscale Water Layer

Achieving timely, reversible, and long-range remote tunability over surface wettability is highly demanded across diverse fields, including nanofluidic systems, drug delivery, and heterogeneous catalysis. Herein, using molecular dynamic simulations, we show, for the first time, a theoretical design of electrowetting to achieve remotely controllable surface wettability via using a terahertz wave. The key idea driving the design is the unique terahertz collective vibration identified in the vicinal subnanoscale water layer, which is absent in bulk water, enabling efficient energy transfer from the terahertz wave to the rotational motion of the vicinal subnanoscale water layer. Consequently, a frequency-specific alternating terahertz electric field near the critical strength can significantly affect the local hydrogen-bonding network of the contact water layer on the solid surface, thereby achieving tunable surface wettability.

Role of Exchange Cations and Layer Charge on the Dynamics of Confined Water

Describing the dynamic behavior of water confined in clay minerals is a fascinating challenge and crucial in many research areas, ranging from materials science and geotechnical engineering to environmental sustainability. Water is the most abundant resource on Earth, and the high reactivity of naturally occurring hydrous clay minerals used since prehistoric times for a variety of applications means that water-clay interaction is a ubiquitous phenomenon in nature. We have attempted to experimentally distinguish the rotational dynamics and translational diffusion of two distinct populations of interlayer water, confined and ultraconfined, in the sodium (Na) forms of two smectite clay minerals, montmorillonite (Mt) and hectorite (Ht). Samples hydrated at a pseudo one-layer hydration (1LH) state under ambient conditions were studied with quasi-elastic neutron scattering (QENS) between 150 and 300 K. Using a simplified revised jump-diffusion and rotation-diffusion model (srJRM), we observed that while interlayer water near the ditrigonal cavity in Ht forms strong H-bonds to both adjacent surface O and structural OH, H-bonding of other more prevalent interlayer water with the surface O is weaker compared to Mt, inducing a higher temperature for dynamical changes of confined water. Given the lower layer charge and faster dynamics observed for Ht compared to Mt, we consider this strong evidence confirming the influence of the interlayer cation and surfaces on confined water dynamics.

Understanding the intricacy of protein in hydrated deep eutectic solvent: Solvation dynamics, conformational fluctuation dynamics, and stability

Deep eutectic solvents (DESs) are potential biocatalytic media due to their easy preparation, fine-tuneability, biocompatibility, and most importantly, due to their ability to keep protein stable and active. However, there are many unanswered questions and gaps in our knowledge about how proteins behave in these alternate media. Herein, we investigated solvation dynamics, conformational fluctuation dynamics, and stability of human serum albumin (HSA) in 0.5 Acetamide/0.3 Urea/0.2 Sorbitol (0.5Ac/0.3Ur/0.2Sor) DES of varying concentrations to understand the intricacy of protein behaviour in DES. Our result revealed a gradual decrease in the side-chain flexibility and thermal stability of HSA beyond 30 % DES. On the other hand, the associated water dynamics around domain-I of HSA decelerate only marginally with increasing DES content, although viscosity rises considerably. We propose that even though macroscopic solvent properties are altered, a protein feels only an aqueous type of environment in the presence of DES. This is probably the first experimental study to delineate the role of the associated water structure of the enzyme for maintaining its stability inside DES. Although considerable effort is necessary to generalize such claims, it might serve as the basis for understanding why proteins remain stable and active in DES.

Osmolyte induced protein stabilization: modulation of associated water dynamics might be a key factor

The mechanism of protein stabilization by osmolytes remains one of the most important and long-standing puzzles. The traditional explanation of osmolyte-induced stability through the preferential exclusion of osmolytes from the protein surface has been seriously challenged by the observations like the concentration-dependent reversal of osmolyte-induced stabilization/destabilization. The more modern explanation of protein stabilization/destabilization by osmolytes considers an indirect effect due to osmolyte-induced distortion of the water structure. It provides a general mechanism, but there are numerous examples of protein-specific effects, i.e., a particular osmolyte might stabilize one protein, but destabilize the other, that could not be rationalized through such an explanation. Herein, we hypothesized that osmolyte-induced modulation of associated water might be a critical factor in controlling protein stability in such a medium. Taking different osmolytes and papain as a protein, we proved that our proposal could explain protein stability in osmolyte media. Stabilizing osmolytes rigidify associated water structures around the protein, whereas destabilizing osmolytes make them flexible. The strong correlation between the stability and the associated water dynamics, and the fact that such dynamics are very much protein specific, established the importance of considering the modulation of associated water structures in explaining the osmolyte-induced stabilization/destabilization of proteins. More interestingly, we took another protein, bromelain, for which a traditionally stabilizing osmolyte, sucrose, acts as a stabilizer at higher concentrations but as a destabilizer at lower concentrations. Our proposal successfully explains such observations, which is probably impossible by any known mechanisms. We believe this report will trigger much research in this area.

Water molecule ordering on the surface of an intrinsically disordered protein

The dynamics and local structure of the hydration water on surfaces of folded proteins have been extensively investigated. However, our knowledge of the hydration of intrinsically disordered proteins (IDPs) is more limited. Here, we compare the local structure of water molecules hydrating a globular protein, lysozyme, and the intrinsically disordered N-terminal of c-Src kinase (SH4UD) using molecular dynamics simulation. The radial distributions from the protein surface of the first and the second hydration shells are similar for the folded protein and the IDP. However, water molecules in the first hydration shell of both the folded protein and the IDP are perturbed from the bulk. This perturbation involves a loss of tetrahedrality, which is, however, significantly more marked for the folded protein than the IDP. This difference arises from an increase in the first hydration shell of the IDP of the fraction of hydration water molecules interacting with oxygen. The water ordering is independent of the compactness of the IDP. In contrast, the lifetimes of water molecules in the first hydration shell increase with IDP compactness, indicating a significant impact of IDP configuration on water surface pocket kinetics, which here is linked to differential pocket volumes and polarities.

Primary Photophysics of Nicotinamide Chromophores in Their Oxidized and Reduced Forms

Nicotinamide adenine dinucleotide (NADH) is an important enzyme cofactor with emissive properties that allow it to be used in fluorescence microscopies to study cell metabolism. Its oxidized form NAD+, on the other hand, is considered to produce negligible fluorescence. In this contribution, we describe the photophysics of the isolated nicotinamidic system in both its reduced and oxidized states. This was achieved through the study of model molecules that do not carry the adenine nucleotide since its absorbance would overlap with the absorption spectrum of the nicotinamidic chromophores. We studied three model molecules: nicotinamide (niacinamide, an oxidized form without nitrogen substitution), the oxidized chromophore 1-benzyl-3-carbamoyl-pyridinium bromide (NBzOx), and its reduced form 1-benzyl-1,4-dihydronicotinamide (NBz). For a full understanding of the dynamics, we performed both femtosecond-resolved emission and transient absorption experiments. The oxidized systems, nicotinamide and NBzOx, have similar photophysics, where the originally excited bright state decays on an ultrafast timescale of less than 400 fs. The depopulation of this state is followed by excited-state positive absorption signals, which evolve in two timescales: the first one is from 1 to a few picoseconds and is followed by a second decaying component of 480 ps for nicotinamide in water and of 80-90 ps for nicotinamide in methanol and NBzOx in aqueous solution. The long decay times are assigned as the S1 lifetimes populated from the original higher-lying bright singlet, where this state is nonemissive but can be detected by transient absorption. While for NBzOx in aqueous solution and for nicotinamide in methanol, the S1 signal decays to the solvent-only level, for the aqueous solutions of nicotinamide, a small transient absorption signal remains after the 480 ps decay. This residual signal was assigned to a small population of triplet states formed during the slower S1 decay for nicotinamide in water. The experimental results were complemented by XMS-CASPT2 calculations, which reveal that in the oxidized forms, the rapid evolution of the initial π-π* state is due to a direct crossing with lower-energy dark n-π* singlet states. This coincides with the experimental observation of long-lived nonemissive states (80 to 480 ps depending on the system). On the other hand, the reduced model compound NBz has a long-lived emissive π-π* S1 state, which decays with a 510 ps time constant, similarly to the parent compound NADH. This is consistent with the XMS-CASPT2 calculations, which show that for the reduced chromophore, the dark states lie at higher energies than the bright π-π* S1 state.

Terahertz spectroscopy as a method for investigation of hydration shells of biomolecules

The hydration of biomolecules is one of the fundamental processes underlying the construction of living matter. The formation of the native conformation of most biomolecules is possible only in an aqueous environment. At the same time, not only water affects the structure of biomolecules, but also biomolecules affect the structure of water, forming hydration shells. However, the study of the structure of biomolecules is given much more attention than their hydration shells. A real breakthrough in the study of hydration occurred with the development of the THz spectroscopy method, which showed that the hydration shell of biomolecules is not limited to 1-2 layers of strongly bound water, but also includes more distant areas of hydration with altered molecular dynamics. This review examines the fundamental features of the THz frequency range as a source of information about the structural and dynamic characteristics of water that change during hydration. The applied approaches to the study of hydration shells of biomolecules based on THz spectroscopy are described. The data on the hydration of biomolecules of all main types obtained from the beginning of the application of THz spectroscopy to the present are summarized. The emphasis is placed on the possible participation of extended hydration shells in the realization of the biological functions of biomolecules and at the same time on the insufficient knowledge of their structural and dynamic characteristics.

The relayed nuclear Overhauser effect in magnetization transfer and chemical exchange saturation transfer MRI

Magnetic resonance (MR) is a powerful technique for noninvasively probing molecular species in vivo but suffers from low signal sensitivity. Saturation transfer (ST) MRI approaches, including chemical exchange saturation transfer (CEST) and conventional magnetization transfer contrast (MTC), allow imaging of low-concentration molecular components with enhanced sensitivity using indirect detection via the abundant water proton pool. Several recent studies have shown the utility of chemical exchange relayed nuclear Overhauser effect (rNOE) contrast originating from nonexchangeable carbon-bound protons in mobile macromolecules in solution. In this review, we describe the mechanisms leading to the occurrence of rNOE-based signals in the water saturation spectrum (Z-spectrum), including those from mobile and immobile molecular sources and from molecular binding. While it is becoming clear that MTC is mainly an rNOE-based signal, we continue to use the classical MTC nomenclature to separate it from the rNOE signals of mobile macromolecules, which we will refer to as rNOEs. Some emerging applications of the use of rNOEs for probing macromolecular solution components such as proteins and carbohydrates in vivo or studying the binding of small substrates are discussed.

Ultrafast Energy Transfer from Photoexcited Tryptophan to the Haem in Cytochrome c

We report femtosecond Fe K-edge absorption (XAS) and nonresonant X-ray emission (XES) spectra of ferric cytochrome C (Cyt c) upon excitation of the haem (>300 nm) or mixed excitation of the haem and tryptophan (<300 nm). The XAS and XES transients obtained in both excitation energy ranges show no evidence for electron transfer processes between photoexcited tryptophan (Trp) and the haem, but rather an ultrafast energy transfer, in agreement with previous ultrafast optical fluorescence and transient absorption studies. The reported (J. Phys. Chem. B 2011, 115 (46), 13723-13730) decay times of Trp fluorescence in ferrous (∼350 fs) and ferric (∼700 fs) Cyt c are among the shortest ever reported for Trp in a protein. The observed time scales cannot be rationalized in terms of Förster or Dexter energy transfer mechanisms and call for a more thorough theoretical investigation.

Ultrafast Förster resonance energy transfer from tyrosine to tryptophan in monellin: potential intrinsic spectroscopic ruler

Ultrafast Förster Resonance Energy Transfer (FRET) between tyrosine (Tyr) and tryptophan (Trp) residues in the protein monellin has been investigated using picosecond and femtosecond time-resolved fluorescence spectroscopy. Decay associated spectra (DAS) and time-resolved emission spectra (TRES) taken with the different excitation wavelengths of 275, 290 and 295 nm were constructed via global analysis. At two of those three excitation loci (275 and 290 nm), earmarks of energy transfer from Tyr to Trp in monellin are seen, and particularly when the excitation is 275 nm, the energy transfer between Tyr and Trp clearly changes the signature emission DAS shape to that indicating excited state reaction (especially on the red side of fluorescence emission, near 380 nm). Those FRET signatures may overlap with the conventional signatory DAS in heterogeneous systems. When overlap and addition occur between FRET type DAS and "full positive" QSSQ (quasi-static self-quenching), mixed DAS shapes will emerge that still show "positive blue side and negative red sides", just with zero crossing shifted. In addition, excitation decay associated spectra (EDAS) taken with the different emission wavelengths of 330, 350 and 370 nm were constructed. In the study of protein dynamics, ultrafast FRET between Tyr and Trp could provide a basis for an intrinsic (non-perturbing) "spectroscopic ruler", potentially a powerful tool to detect even slight changes in protein structures.

Polymorphic protein phase transitions driven by surface anisotropy

Phase transitions of proteins are strongly influenced by surface chemical modifications or mutations. Human γD-crystallin (HGD) single-mutants have been extensively studied because they are associated with the onset of juvenile cataract. However, they have also provided a rich library of molecules to examine how specific inter-protein interactions direct protein assembly, providing new insights and valuable experimental data for coarse-grained patchy-particle models. Here, we demonstrate that the addition of new inter-protein interactions by mutagenesis is additive and increases the number and variety of condensed phases formed by proteins. When double mutations incorporating two specific single point mutations are made, the properties of both single mutations are retained in addition to the formation of a new condensed phase. We find that the HGD double-mutant P23VC110M self-assembles into spherical particles with retrograde solubility, orthorhombic crystals, and needle/plate shape crystals, while retaining the ability to undergo liquid-liquid phase separation. This rich polymorphism is only partially predicted by the experimental data on the constituent single mutants. We also report a previously un-characterized amorphous protein particle, with unique properties that differ from those of protein spherulites, protein particulates previously described. The particles we observe are amorphous, reversible with temperature, tens of microns in size, and perfectly spherical. When they are grown on pristine surfaces, they appear to form by homogeneous nucleation, making them unique, and we believe a new form of protein condensate. This work highlights the challenges in predicting protein behavior, which has frustrated rational assembly and crystallization but also provides rich data to develop new coarse-grained models to explain the observed polymorphism.

Molecular Dynamics Simulations of Rhodamine B Zwitterion Diffusion in Polyelectrolyte Solutions

Polyelectrolytes continue to find wide interest and application in science and engineering, including areas such as water purification, drug delivery, and multilayer thin films. We have been interested in the dynamics of small molecules in a variety of polyelectrolyte (PE) environments; in this paper, we report simulations and analysis of the small dye molecule rhodamine B (RB) in several very simple polyelectrolyte solutions. Translational diffusion of the RB zwitterion has been measured in fully atomistic, 2 μs long molecular dynamics simulations in four different polyelectrolyte solutions. Two solutions contain the common polyanion sodium poly(styrene sulfonate) (PSS), one with a 30-mer chain and the other with 10 trimers. The other two solutions contain the common polycation poly(allyldimethylammonium) chloride (PDDA), one with two 15-mers and the other with 10 trimers. RB diffusion was also simulated in several polymer-free solutions to verify its known experimental value for the translational diffusion coefficient, D_RB, of 4.7 × 10^-6 cm²/s at 300 K. RB diffusion was slowed in all four simulated PE solutions, but to varying degrees. D_RB values of 3.07 × 10^-6 and 3.22 × 10^-6 cm²/s were found in PSS 30-mer and PSS trimer solutions, respectively, whereas PDDA 15-mer and trimer solutions yielded values of 2.19 × 10^-6 and 3.34 × 10^-6 cm²/s. Significant associations between RB and the PEs were analyzed and interpreted via a two-state diffusion model (bound and free diffusion) that describes the data well. Crowder size effects and anomalous diffusion were also analyzed. Finally, RB translation along the polyelectrolytes during association was characterized.

Ultrafast light-induced dynamics in the microsolvated biomolecular indole chromophore with water

Interactions between proteins and their solvent environment can be studied in a bottom-up approach using hydrogen-bonded chromophore-solvent clusters. The ultrafast dynamics following UV-light-induced electronic excitation of the chromophores, potential radiation damage, and their dependence on solvation are important open questions. The microsolvation effect is challenging to study due to the inherent mix of the produced gas-phase aggregates. We use the electrostatic deflector to spatially separate different molecular species in combination with pump-probe velocity-map-imaging experiments. We demonstrate that this powerful experimental approach reveals intimate details of the UV-induced dynamics in the near-UV-absorbing prototypical biomolecular indole-water system. We determine the time-dependent appearance of the different reaction products and disentangle the occurring ultrafast processes. This approach ensures that the reactants are well-known and that detailed characteristics of the specific reaction products are accessible - paving the way for the complete chemical-reactivity experiment.

Non-Aromatic Fluorescence in Biological Matter: The Exception or the Rule?

While in the vast majority of cases fluorescence in biological matter has been attributed to aromatic or conjugated groups, peptides associated with neurodegenerative diseases, such as Alzheimer's, Parkinson's, or Huntington's, have been recently shown to display an intrinsic visible fluorescence even in the absence of aromatic residues. This has called the attention of researchers from many different fields, trying to understand the origin of this peculiar behavior and, at the same time, motivating the search for novel strategies to control the optical properties of new biophotonic materials. Today, after nearly 15 years of its discovery, there is a growing consensus about the mechanism underlying this phenomenon, namely, that electronic interactions between non-optically active molecules can result in supramolecular assemblies that are fluorescent. Despite this progress, many aspects of this phenomenon remain uncharted territory. In this Perspective, we lay down the state-of-the-art in the field highlighting the open questions from both experimental and theoretical fronts in this fascinating emerging area of non-aromatic fluorescence.

Ultrafast Förster resonance energy transfer between tyrosine and tryptophan: potential contributions to protein-water dynamics measurements

Ultrafast Förster Resonance Energy Transfer (FRET) between Tyrosine (Tyr, Y) and Tryptophan (Trp, W) in the model peptides Trp-(Pro)_n-Tyr (WPnY) has been investigated using a femtosecond up-conversion spectrophotofluorometer. The ultrafast energy transfer process (<100 ps) in short peptides (WY, WPY and WP2Y) has been resolved. In fact, this FRET rate is found to be mixed with the rates of solvent relaxation (SR), ultrafast population decay (QSSQ) and other lifetime components. To further dissect and analyze the FRET, a spectral working model is constructed, and the contribution of a FRET lifetime is separated by reconciling the shapes of decay associated spectra (DAS). Surprisingly, FRET efficiency did not decrease monotonically with the growth of the peptide chain (as expected) but increased first and then decreased. The highest FRET efficiency occurred in peptide WPY. The kinetic results have been accompanied with molecular dynamics simulations that reconcile and explain this strange phenomenon: due to the strong interaction between amino acids, the distance between the donor and receptor in peptide WPY is actually closest, resulting in the fastest FRET. In addition, the FRET lifetimes (τ_cal) were estimated within the molecular dynamics simulations, and they were consistent with the lifetimes (τ_exp) separated out by the experimental measurements and the DAS working model. This benchmark study has implications for both previous and future studies of protein ultrafast dynamics. The approach taken can be generalized for the study of proximate tyrosine and tryptophan in proteins and it suggests spectral strategies for extracting mixed rates in other complex FRET problems.

Stochastic Analysis of Molecular Dynamics Reveals the Rotation Dynamics Distribution of Water around Lysozyme

Water dynamics is essential to biochemical processes by mediating all such reactions, including biomolecular degeneration in solutions. To disentangle the molecular-scale distribution of water dynamics around a solute biomolecule, we investigated here the rotational dynamics of water around lysozyme by combining molecular dynamics (MD) simulations and broadband dielectric spectroscopy (BDS). A statistical analysis using the relaxation times and trajectories of every single water molecule was proposed, and the two-dimensional probability distribution of water at a distance from the lysozyme surface with a rotational relaxation time was given. For the observed lysozyme solutions of 34-284 mg/mL, we discovered that the dielectric relaxation time obtained from this distribution agrees well with the measured γ relaxation time, which suggests that rotational self-correlation of water molecules underlies the gigahertz domain of the dielectric spectra. Regardless of protein concentration, water rotational relaxation time versus the distance from the lysozyme surface revealed that the water rotation is severely retarded within 3 Å from the lysozyme surface and is nearly comparable to pure water when farther than 10 Å. The dimension of the first hydration layer was subsequently identified in terms of the relationship between the acceleration of water rotation and the distance from the protein surface.

Ultraslow Biological Water-Like Dynamics in Waterless Liquid Protein

Fluorescence correlation spectroscopy and time-dependent fluorescence Stokes shift have been employed to elucidate dynamics in different time scales, ranging from picoseconds to nanoseconds, for human serum albumin, in its native and cationized forms as well as in the self-assembled complex of the cationized protein with the polymer surfactant (PS) glycolic acid ethoxylate lauryl ether. The effect of crowding in this complex, especially in the waterless condition, is of prime importance in this context. Excellent correlation of the dynamics with the structures, obtained by circular dichroism and Fourier transform infrared spectroscopy, has been observed. Slow solvation, associated classically with biological water, has been observed in these systems, even in the waterless condition. This apparently intriguing observation has been rationalized by the relaxation of segments of the protein and the PS in the microenvironment of the fluorescent probe.

Ultrafast Excited-State Proton Transfer of a Cationic Superphotoacid in a Nanoscopic Water Pool

The excited-state proton transfer (ESPT) of a cationic superphotoacid, N-methyl-7-hydroxyquinolium, was studied within the water pool of an anionic aerosol-OT (AOT), bis(2-ethylhexyl) sulfosuccinate, reverse micelle (RM). Previously, we had found that the cationic photoacid residing at the anionic AOT interface was conducive to ESPT to the bound water having concentric heterogeneity on the time scale of hundreds of picoseconds to nanoseconds. In our present study, on the time scale of hundreds of femtoseconds to a few tens of picoseconds, the photoacid underwent an ultrafast ESPT influenced by mobile water constituting the core of the RM. The two subpopulations of the core water molecules that determine the ultrafast biphasic deprotonation of the photoacid on time scales differing by an order of magnitude were identified. The core water molecules solvating the counteranion of the photoacid showed a higher basicity than typical water clusters in bulk resulting in ESPT on a subpicosecond time scale. Bare water clusters sensed by the photoacid showed a slower ESPT, over several picoseconds, as typically limited by the rotational motion of water molecules for similar types of the photoacid.

Conformational and Solvation Dynamics of an Amyloidogenic Intrinsically Disordered Domain of a Melanosomal Protein

The conformational plasticity of intrinsically disordered proteins (IDPs) allows them to adopt a range of conformational states that can be important for their biological functions. The driving force for the conformational preference of an IDP emanates from an intricate interplay between chain-chain and chain-solvent interactions. Using ultrafast femtosecond and picosecond time-resolved fluorescence measurements, we characterized the conformational and solvation dynamics around the N- and C-terminal segments of a disordered repeat domain of a melanosomal protein Pmel17 that forms functional amyloid responsible for melanin biosynthesis. Our time-resolved fluorescence anisotropy results revealed slight compaction and slower rotational dynamics around the amyloidogenic C-terminal segment when compared to the proline-rich N-terminal segment of the repeat domain. The compaction of the C-terminal region was also associated with the restrained mobility of hydration water as indicated by our solvation dynamics measurements. Our findings indicate that sequence-dependent chain-solvent interactions govern both the conformational and solvation dynamics that are crucial in directing the conversion of a highly dynamic IDP into an ordered amyloid assembly. Such an interplay of amino acid composition-dependent conformational and solvation dynamics might have important physicochemical consequences in specific water-protein, ion-protein, and protein-protein interactions involved in amyloid formation and phase transitions.

Learning the Hydrophobic, Hydrophilic, and Aromatic Character of Amino Acids from Thermal Relaxation and Interfacial Thermal Conductance

In this study, the thermal relaxation of the 20 naturally occurring amino acids in water and in the protein lysozyme is investigated using transient nonequilibrium molecular dynamics simulations. By modeling the thermal relaxation process, the relaxation times of the amino acids in water occurs over a time scale covering 2-5 ps. For the hydrophobic amino acids, the relaxation time is controlled by the size of the hydrocarbon side chain, while for hydrophilic amino acids, the number of hydrogen bonds does not significantly affect the time scales of the heat dissipation. Our results show that the interfacial thermal conductance at the amino acid-water interface is in the range of 40-80 MW m^-2 K^-1. Hydrophobic and aromatic amino acids tend to have a lower interfacial thermal conductance. Notably, we show that amino acids can be correlated with their thermal relaxation times and molar masses, into simply connected phases with the same hydrophilicity, hydrophobicity, and aromaticity. The thermal relaxation slows down by a factor of up to five in the protein relative to that in water. In the case of the hydrophobic amino acids in the protein lysozyme, the slow down in the thermal relaxation relative to that in water appears to be controlled primarily by the size of the side chain.

The Role of Hydrogen Bonding in the Preferential Solvation of 5-Aminoquinoline in Binary Solvent Mixtures

5-Aminoquinoline (5AQ) has been used as a fluorescent probe of preferential solvation (PS) in binary solvent mixtures in which the nonpolar component is diethyl ether and the polar component is protic (methanol) or aprotic (acetonitrile). Hence, the roles of solvent polarity and solute-solvent hydrogen bonding have been delineated. Positive deviations of spectral shifts from a linear dependence on the concentration of the polar component, signifying PS, are markedly more pronounced in case of the protic solvent. Solvation dynamics on a nanosecond time scale mark the formation of the solvation shell around the fluorescent probe. Time-resolved area-normalized emission spectra indicate the occurrence of the continuous solvation of the excited state when the polar component is acetonitrile. In contrast, two distinct states were observed when the polar component was methanol, the second state being the hydrogen bonded one. Translational diffusion is the rate-determining step for formation of the solvation shell. The time constant associated with it has been estimated from rise times observed in fluorescence transients monitored at the red end of the fluorescence spectra and also from the time evolution of the spectral width of time-resolved emission spectra.

What Does Time-Dependent Fluorescence Shift (TDFS) in Biomembranes (and Proteins) Report on?

The organization of biomolecules and bioassemblies is highly governed by the nature and extent of their interactions with water. These interactions are of high intricacy and a broad range of methods based on various principles have been introduced to characterize them. As these methods view the hydration phenomena differently (e.g., in terms of time and length scales), a detailed insight in each particular technique is to promote the overall understanding of the stunning “hydration world.” In this prospective mini-review we therefore critically examine time-dependent fluorescence shift (TDFS)—an experimental method with a high potential for studying the hydration in the biological systems. We demonstrate that TDFS is very useful especially for phospholipid bilayers for mapping the interfacial region formed by the hydrated lipid headgroups. TDFS, when properly applied, reports on the degree of hydration and mobility of the hydrated phospholipid segments in the close vicinity of the fluorophore embedded in the bilayer. Here, the interpretation of the recorded TDFS parameters are thoroughly discussed, also in the context of the findings obtained by other experimental techniques addressing the hydration phenomena (e.g., molecular dynamics simulations, NMR spectroscopy, scattering techniques, etc.). The differences in the interpretations of TDFS outputs between phospholipid biomembranes and proteins are also addressed. Additionally, prerequisites for the successful TDFS application are presented (i.e., the proper choice of fluorescence dye for TDFS studies, and TDFS instrumentation). Finally, the effects of ions and oxidized phospholipids on the bilayer organization and headgroup packing viewed from TDFS perspective are presented as application examples.

Structural Stability of Insulin Oligomers and Protein Association-Dissociation Processes: Free Energy Landscape and Universal Role of Water

Association and dissociation of proteins are important biochemical events. In this Feature Article, we analyze the available studies of these processes for insulin oligomers in aqueous solution. We focus on the solvation of the insulin monomer in water, stability and dissociation of its dimer, and structural integrity of the hexamer. The intricate role of water in solvation of the dimer- and hexamer-forming surfaces, in long-range interactions between the monomers and the stability of the oligomers, is discussed. Ten water molecules inside the central cavity stabilize the structure of the insulin hexamer. We discuss how different order parameters can be used to understand the dissociation of the insulin dimer. The calculation of the rate using a recently computed multidimensional free energy provides considerable insight into the interplay between protein and water dynamics.

The functional role of the hemoglobin-water interface

The interface between hemoglobin (Hb) and its environment, in particular water, is of great physiological relevance. Here, results from in vitro, in vivo, and computational experiments (molecular dynamics simulations) are summarized and put into perspective. One of the main findings from the computations is that the stability of the deoxy, ligand-free T-state (T₀) can be stabilized relative to the deoxy R-state (R₀) only in sufficiently large simulation boxes for the hydrophobic effect to manifest itself. This effect directly influences protein stability and is operative also under physiological conditions. Furthermore, molecular simulations provide a dynamical interpretation of the Perutz model for Hb function. Results from experiments using higher protein concentrations and realistic cellular environments are also discussed. One of the next great challenges for computational studies, which as we show is likely to be taken up in the near future, is to provide a molecular-level understanding of the dynamics of proteins in such crowded environments.

Machine-learning model selection and parameter estimation from kinetic data of complex first-order reaction systems

Dealing with a system of first-order reactions is a recurrent issue in chemometrics, especially in the analysis of data obtained by spectroscopic methods applied on complex biological systems. We argue that global multiexponential fitting, the still common way to solve such problems, has serious weaknesses compared to contemporary methods of sparse modeling. Combining the advantages of group lasso and elastic net-the statistical methods proven to be very powerful in other areas-we created an optimization problem tunable from very sparse to very dense distribution over a large pre-defined grid of time constants, fitting both simulated and experimental multiwavelength spectroscopic data with high computational efficiency. We found that the optimal values of the tuning hyperparameters can be selected by a machine-learning algorithm based on a Bayesian optimization procedure, utilizing widely used or novel versions of cross-validation. The derived algorithm accurately recovered the true sparse kinetic parameters of an extremely complex simulated model of the bacteriorhodopsin photocycle, as well as the wide peak of hypothetical distributed kinetics in the presence of different noise levels. It also performed well in the analysis of the ultrafast experimental fluorescence kinetics data detected on the coenzyme FAD in a very wide logarithmic time window. We conclude that the primary application of the presented algorithms-implemented in available software-covers a wide area of studies on light-induced physical, chemical, and biological processes carried out with different spectroscopic methods. The demand for this kind of analysis is expected to soar due to the emerging ultrafast multidimensional infrared and electronic spectroscopic techniques that provide very large and complex datasets. In addition, simulations based on our methods could help in designing the technical parameters of future experiments for the verification of particular hypothetical models.

Cooperativity and Anticooperativity in Ion-Water Interactions: Implications for the Aqueous Solvation of Ions

Non-additive effects in hydrogen bonds (HB) take place as a consequence of electronic charge transfers. Therefore, it is natural to expect cooperativity and anticooperativity in ion-water interactions. Nevertheless, investigations on this matter are scarce. This paper addresses the interactions of (i) the cations Li⁺ , Na⁺ , K⁺ , Be²⁺ , Mg²⁺ , and Ca²⁺ together with (ii) the anions F^- , Cl^- , Br^- , NO₃ ^- and SO₄ ^2- with water clusters (H₂ O)_n , n=1-8, and the effects of these ions on the HBs within the complete molecular adducts. We used quantum chemical topology tools, specifically the quantum theory of atoms in molecules and the interacting quantum atoms energy partition to investigate non-additive effects among the interactions studied herein. Our results show a decrease on the interaction energy between ions and the first neighbouring water molecules with an increment of the coordination number. We also found strong cooperative effects in the interplay between HBs and ion-dipole interactions within the studied systems. Such cooperativity affects considerably the interactions among ions with their first and second solvation shells in aqueous environments. Overall, we believe this article provides valuable information about how ion-dipole contacts interact with each other and how they relate to other interactions, such as HBs, in the framework of non-additive effects in aqueous media.

THE IMS PARADOX: A PERSPECTIVE ON STRUCTURAL ION MOBILITY-MASS SPECTROMETRY

Studies of large proteins, protein complexes, and membrane protein complexes pose new challenges, most notably the need for increased ion mobility (IM) and mass spectrometry (MS) resolution. This review covers evolutionary developments in IM-MS in the authors' and key collaborators' laboratories with specific focus on developments that enhance the utility of IM-MS for structural analysis. IM-MS measurements are performed on gas phase ions, thus "structural IM-MS" appears paradoxical-do gas phase ions retain their solution phase structure? There is growing evidence to support the notion that solution phase structure(s) can be retained by the gas phase ions. It should not go unnoticed that we use "structures" in this statement because an important feature of IM-MS is the ability to deal with conformationally heterogeneous systems, thus providing a direct measure of conformational entropy. The extension of this work to large proteins and protein complexes has motivated our development of Fourier-transform IM-MS instruments, a strategy first described by Hill and coworkers in 1985 (Anal Chem, 1985, 57, pp. 402-406) that has proved to be a game-changer in our quest to merge drift tube (DT) and ion mobility and the high mass resolution orbitrap MS instruments. DT-IMS is the only method that allows first-principles determinations of rotationally averaged collision cross sections (CSS), which is essential for studies of biomolecules where the conformational diversities of the molecule precludes the use of CCS calibration approaches. The Fourier transform-IM-orbitrap instrument described here also incorporates the full suite of native MS/IM-MS capabilities that are currently employed in the most advanced native MS/IM-MS instruments. © 2020 John Wiley & Sons Ltd. Mass Spec Rev.

A Comparative Study on DMSO-Induced Modulation of the Structural and Dynamical Properties of Model Bilayer Membranes

Modulating the structures and properties of biomembranes via permeation of small amphiphilic molecules is immensely important, having diverse applications in cell biology, biotechnology, and pharmaceuticals, because their physiochemical and biological interactions lead to new pathways for transdermal drug delivery and administration. In this work, we have elucidated the role of dimethyl sulfoxide (DMSO), broadly used as a penetration-enhancing agent and cryoprotective agent on model lipid membranes, using a combination of fluorescence microscopy and time-resolved fluorescence spectroscopy. Spatially resolved fluorescence lifetime imaging microscopy (FLIM) has been employed to unravel how the fluidity of the DMSO-induced bilayer regulates the structural alteration of the vesicles. Moreover, we have also shown that the dehydration effect of DMSO leads to weakening of the hydrogen bond between lipid headgroups and water molecules and results in faster solvation dynamics as demonstrated by femtosecond time-resolved fluorescence spectroscopy. It has been gleaned that the water dynamics becomes faster because bilayer rigidity decreases in the presence of DMSO, which is also supported by time-resolved rotational anisotropy measurements. The enhanced diffusivity and increased membrane fluidity in the presence of DMSO are further ratified at the single-molecule level through fluorescence correlation spectroscopy (FCS) measurements. Our results indicate that while the presence of DMSO significantly affects the 1,2-dimyristoyl-rac-glycero-3-phosphocholine (DMPC) and 1,2-dipalmitoyl-rac-glycero-3-phosphatidylcholine (DPPC) bilayers, it has a weak effect on 1,2-dimyristoyl-sn-glycero-3-phospho-rac-glycerol (DMPG) vesicles, which might explain the preferential interaction of DMSO with the positively charged choline group present in DMPC and DPPC vesicles. The experimental findings have also been further verified with molecular dynamics simulation studies. Moreover, it has been observed that DMSO is likely to have a differential effect on heterogeneous bilayer membranes depending on the structure and composition of their headgroups. Our results illuminate the importance of probing the lipid structure and composition of cellular membranes in determining the effects of cryoprotective agents.

Active-Site Environmental Factors Customize the Photophysics of Photoenzymatic Old Yellow Enzymes

The development of non-natural photoenzymatic systems has reinvigorated the study of photoinduced electron transfer (ET) within protein active sites, providing new and unique platforms for understanding how biological environments affect photochemical processes. In this work, we use ultrafast spectroscopy to compare the photoinduced electron transfer in known photoenzymes. 12-Oxophytodienoate reductase 1 (OPR1) is compared to Old Yellow Enzyme 1 (OYE1) and morphinone reductase (MR). The latter enzymes are structurally homologous to OPR1. We find that slight differences in the amino acid composition of the active sites of these proteins determine their distinct electron-transfer dynamics. Our work suggests that the inside of a protein active site is a complex/heterogeneous dielectric network where genetically programmed heterogeneity near the site of biological ET can significantly affect the presence and lifetime of various intermediate states. Our work motivates additional tunability of Old Yellow Enzyme active-site reorganization energy and electron-transfer energetics that could be leveraged for photoenzymatic redox approaches.

Protein-Water and Water-Water Long-Time Relaxations in Protein Hydration Water upon Cooling-A Close Look through Density Correlation Functions

We report results on the translational dynamics of the hydration water of the lysozyme protein upon cooling obtained by means of molecular dynamics simulations. The self van Hove functions and the mean square displacements of hydration water show two different temperature activated relaxation mechanisms, determining two dynamic regimes where transient trapping of the molecules is followed by hopping phenomena to allow to the structural relaxations. The two caging and hopping regimes are different in their nature. The low-temperature hopping regime has a time scale of tenths of nanoseconds and a length scale on the order of 2–3 water shells. This is connected to the nearest-neighbours cage effect and restricted to the supercooling, it is absent at high temperature and it is the mechanism to escape from the cage also present in bulk water. The second hopping regime is active at high temperatures, on the nanoseconds time scale and over distances of nanometers. This regime is connected to water displacements driven by the protein motion and it is observed very clearly at high temperatures and for temperatures higher than the protein dynamical transition. Below this temperature, the suppression of protein fluctuations largely increases the time-scale of the protein-related hopping phenomena at least over 100 ns. These protein-related hopping phenomena permit the detection of translational motions of hydration water molecules longly persistent in the hydration shell of the protein.

Effect of Relaxation Times in a Preservative Solution on Protein Deterioration Rate

In this study, the relationship between the deterioration of the enzymatic activity of the protein in a preservative solution and the molecular mobility of the preservative solution was experimentally investigated by dielectric spectroscopy and a protein deterioration test. Dielectric spectroscopy was used to observe the molecular mobilities in the preservative solutions including various ratios of two protective agents, trehalose and ε-poly-l-lysine, at 40 °C. We also examined the enzymatic activity of l-lactate dehydrogenase (LDH) at 40 °C for 40 days to obtain the deterioration rate of LDH in the same preservative solutions. Our experimental results reveal that (1) three relaxation times of water molecular motion are detected by dielectric spectroscopy that we categorize into bulky water (relaxation time of 10 ps), weakly hydration water (relaxation time of tens to hundreds of picoseconds), and strongly hydration water (relaxation time of hundreds of picoseconds to tens of microseconds) and (2) the deterioration rate of LDH has a power-law relationship with the relaxation times of bulky and hydration water with specific power indices. The results also support the protein stabilization theory of high viscosity and the practical advantage of predicting the shelf life of proteins in the preservative solution by the relaxation time of water measured by dielectric spectroscopy.

Solvation Dynamics in Water. 4. On the Initial Regime of Solvation Relaxation

It is shown, by means of numerical and analytic work, that initial molecular momenta play little significant role in the initial fast solvation relaxation that follows electronic excitation of, and charge creation for, a standard model system of a solute in water. Instead, the nonequilibrium dynamics are predominantly described by noninertial "steering" by the torques directly generated by the newly created charge distribution. It is this process that largely overcomes inertia and drives the relaxation dynamics on a time scale of a few tens of femtoseconds in the key initial regime of the dynamics. These results are discussed in the context of commonly employed descriptions such as inertial, Gaussian, and underdamped dynamical behavior.

Length-scale dependence of protein hydration-shell density

Here we present a computational approach based on molecular dynamics (MD) simulation to study the dependence of the protein hydration-shell density on the size of the protein molecule. The hydration-shell density of eighteen different proteins, differing in size, shape and function (eight of them are antifreeze proteins), is calculated. The results obtained show that an increase in the hydration-shell density, relative to that of the bulk, is observed (in the range of 4-14%) for all studied proteins and that this increment strongly correlates with the protein size. In particular, a decrease in the density increment is observed for decreasing protein size. A simple model is proposed in which the basic idea is to approximate the protein molecule as an effective ellipsoid and to partition the relevant parameters, i.e. the solvent-accessible volume and the corresponding solvent density, into two regions: inside and outside the effective protein ellipsoid. It is found that, within the model developed here, almost all of the hydration-density increase is located inside the protein ellipsoid, basically corresponding to pockets within, or at the surface of the protein molecule. The observed decrease in the density increment is caused by the protein size only and no difference is found between antifreeze and non-antifreeze proteins.

Exploring the 2D-IR repertoire of the -SCN label to study site-resolved dynamics and solvation in the calcium sensor protein calmodulin

The calcium sensor protein calmodulin is ubiquitous among eukaryotes. It translates intracellular Ca2+ influx (by a decrease of conformational flexibility) into increased target recognition affinity. Here we demonstrate that by using the IR reporter -SCN in combination with 2D-IR spectroscopy, global structure changes and local dynamics, degree of solvent exposure and protein-ligand interaction can be characterised in great detail. The long vibrational lifetime of the -SCN label allows for centerline slope analysis of the 2D-IR line shape up to 120 ps to deduce the frequency-frequency correlation function (FFCF) of the -SCN label in various states and label positions in the protein. Based on that we show clear differences between a solvent exposed site, the environment close to the Ca2+ binding motif and three highly conserved positions for ligand binding. Furthermore, we demonstrate how these dynamics are affected by conformational change induced by the addition of Ca2+ ions and by interaction with a short helical peptide mimicking protein binding. We show that the binding mode is strongly heterogeneous among the probed key binding methionine residues. SCN's vibrational relaxation is dominated by intermolecular contributions. Changes in the vibrational lifetime upon changing between H2O and D2O buffer therefore provide a robust measure for water accessibility of the label. Characterising -SCN's extinction coefficient, vibrational lifetime in light and heavy water and its FFCF we demonstrate the vast potential it has as a label especially for nonlinear spectroscopies, such as 2D-IR spectroscopy.

Exploring long-range proton conduction, the conduction mechanism and inner hydration state of protein biopolymers

Proteins are the main proton mediators in various biological proton circuits. Using proteins for the formation of long-range proton conductors is offering a bioinspired approach for proton conductive polymers. One of the main challenges in the field of proton conductors is to explore the local environment within the polymers, along with deciphering the conduction mechanism. Here, we show that the protonic conductivity across a protein-based biopolymer can be hindered using straightforward chemical modifications, targeting carboxylate- or amine-terminated residues of the protein, as well as exploring the effect of surface hydrophobicity on proton conduction. We further use the natural tryptophan residue as a local fluorescent probe for the inner local hydration state of the protein surface and its tendency to form hydrogen bonds with nearby water molecules, along with the dynamicity of the process. Our electrical and spectroscopic measurements of the different chemically-modified protein materials as well as the material at different water-aprotic solvent mixtures result in our fundamental understanding of the proton mediators within the material and gaining important insights on the proton conduction mechanism. Our biopolymer can be used as an attractive platform for the study of bio-related protonic circuits as well as a proton conducting biopolymer for various applications, such as protonic transistors, ionic transducers and fuel cells.

In situ molecular imaging of adsorbed protein films in water indicating hydrophobicity and hydrophilicity

In situ molecular imaging of protein films adsorbed on a solid surface in water was realized by using a vacuum compatible microfluidic interface and time-of-flight secondary ion mass spectrometry (ToF-SIMS). Amino acid fragments from such hydrated protein films are observed and identified in the positive ion mode and the results are in agreement with reported works on dry protein films. Moreover, water clusters from the hydrated protein films have been observed and identified in both the positive and negative ion mode for a series protein films. Thus, the detailed composition of amino acids and water molecules in the hydrated protein films can be characterized, and the protein water microstructures can be revealed by the distinct three-dimensional spatial distribution reconstructed from in situ liquid ToF-SIMS molecular imaging. Furthermore, spectral principal component analysis of amino acid fragment peaks and water cluster peaks provides unique insights into the water cluster distribution, hydrophilicity, and hydrophobicity of hydrated adsorbed protein films in water.