Water Determines the Structure and Dynamics of Proteins

Exploring the effect of high-pressure processing conditions on the deaggregation of natural major royal jelly proteins (MRJPs) fibrillar aggregates

High pressure processing is a safe and green novel non-thermal processing technique for modulating food protein aggregation behavior. However, the systematic relationship between high pressure processing conditions and protein deaggregation has not been sufficiently investigated. Major royal jelly proteins, which are naturally highly fibrillar aggregates, and it was found that the pressure level and exposure time could significantly promote protein deaggregation. The 100-200 MPa treatment favoured the deaggregation of proteins with a significant decrease in the sulfhydryl group content. Contrarily, at higher pressure levels (>400 MPa), the exposure time promoted the formation of disordered agglomerates. Notably, the inter-conversion of α-helix and β-strands in major royal jelly proteins after high pressure processing eliminates the solvent-free cavities inside the aggregates, which exerts a 'collapsing' effect on the fibrillar aggregates. Furthermore, the first machine learning model of the high pressure processing conditions and the protein deaggregation behaviour was developed, which provided digital guidance for protein aggregation regulation.

Freezing of Water Solvation Dynamics in Nanoconfinement by Reverse Micelles at Room Temperature

Much attention has recently been paid to anomalously low dielectric constants of nanoconfined water between two slabs at room temperature (Fumagalli et al. Science, 2018, 360, 1339). These low values imply that the dipole rotation of the interfacial water on the slab is completely suppressed. Such freezing has so far been observed for water confined between solids. In contrast, it remains unclear whether this holds for water in soft confinement, which is omnipresent naturally and artificially. Here, we address this question using encapsulated reverse micelles with a dye molecule, allowing us to study water sandwiched between the surfactant and dye molecules in solution. Moreover, we examine the solvation related to the dielectric property of water, which is reorientational motion in the hydration layer of the dye molecule, by persistent hole-burning spectroscopy. We first show that the dye molecule is surrounded by water without contact with the surfactant and that the dye molecule has two or three hydration layers on average. We next demonstrate that the solvation dynamics is frozen below the water droplet size of ∼4 nm, whereas they become liquid-like when the RM size is further increased. The average gap distance (∼1.5 nm) for freezing the solvation agrees with the gap distance with no rotational water motions between slabs. Our findings may have biological relevance, providing a new aspect for understanding biological function in cells.

MANUSCRIPT Local Crowd, Local Probe: Strengths and Drawbacks of Azidohomoalanine as a Site-Specific Crowding Probe

Every residue on a protein can be characterized by its interaction with water, in lack or in excess, as water is the matrix of biological systems. Infrared spectroscopy and the implementation of local azidohomoalanine (AHA) probes allow us to move beyond an ensemble or surface-driven conceptualization of water behavior and toward a granular, site-specific picture. In this paper, we examined the role of crowding in modulating both global and local behavior on the β-hairpin, TrpZip2 using a combination of Fourier-transform infrared spectroscopy (FTIR) spectroscopy, two-dimensional infrared (2D IR) spectroscopy, and molecular dynamics simulations. We found that, at the amino acid level, crowding drove dehydration of both sheet and turn peptide sites as well as free AHA. However, the subpicosecond dynamics showed highly individualized responses based on the local environment. Interestingly, while steady-state FTIR measurements revealed similar responses at the amino-acid level to hard versus soft crowding (dehydration), we found that PEG and glucose had opposite stabilizing and destabilizing effects on the protein secondary structure, emphasizing an important distinction in understanding the impact of crowding on protein structure as well as the role of crowding across length scales.

Exploration of the Existence Forms and Patterns of Dissolved Oxygen Molecules in Water

The structure of liquid water is primarily composed of three-dimensional networks of water clusters formed by hydrogen bonds, and dissolved oxygen is one of the most important indicators for assessing water quality. In this work, distilled water with different concentration of dissolved oxygen were prepared, and a clear negative correlation between the size of water clusters and dissolved oxygen concentration was observed. Besides, a phenomenon of rapid absorption and release of oxygen at the water interfaces was unveiled, suggesting that oxygen molecules predominantly exist at the interfaces of water clusters. Oxygen molecules can move rapidly through the interfaces among water clusters, allowing dissolved oxygen to quickly reach a saturation level at certain partial pressure of oxygen and temperature. Further exploration into the mechanism by molecular dynamics simulations of oxygen and water clusters found that oxygen molecules can only exist stably at the interfaces among water clusters. A semi-empirical formula relating the average number of water molecules in a cluster (n) to 17O NMR half-peak width (W) was summarized: n = 0.1 W + 0.85. These findings provide a foundation for exploring the structure and properties of water.

Nanoporous Gold Thin Films as Substrates to Analyze Liquids by Cryo-atom Probe Tomography

Cryogenic atom probe tomography (cryo-APT) is being developed to enable nanoscale compositional analyses of frozen liquids. Yet, the availability of readily available substrates that allow for the fixation of liquids while providing sufficient strength to their interface is still an issue. Here, we propose the use of 1-2-µm-thick binary alloy film of gold-silver sputtered onto flat silicon, with sufficient adhesion without an additional layer. Through chemical dealloying, we successfully fabricate a nanoporous substrate, with an open-pore structure, which is mounted on a microarray of Si posts by lift-out in the focused-ion beam system, allowing for cryogenic fixation of liquids. We present cryo-APT results obtained after cryogenic sharpening, vacuum cryo-transfer, and analysis of pure water on the top and inside the nanoporous film. We demonstrate that this new substrate has the requisite characteristics for facilitating cryo-APT of frozen liquids, with a relatively lower volume of precious metals. This complete workflow represents an improved approach for frozen liquid analysis, from preparation of the films to the successful fixation of the liquid in the porous network, to cryo-APT.

Cross-Linking Methods of the Silk Protein Hydrogel in Oral and Craniomaxillofacial Tissue Regeneration

BACKGROUND

Craniomaxillofacial tissue defects are clinical defects involving craniomaxillofacial and oral soft and hard tissues. They are characterized by defect-shaped irregularities, bacterial and inflammatory environments, and the need for functional recovery. Conventional clinical treatments are currently unable to achieve regeneration of high-quality oral craniomaxillofacial tissue. As a natural biomaterial, silk fibroin (SF) has been widely studied in biomedicine and has broad prospects for use in tissue regeneration. Hydrogels made of SF showed excellent water retention, biocompatibility, safety and the ability to combine with other materials.

METHODS

To gain an in-depth understanding of the current development of SF, this article reviews the structure, preparation and application prospects in oral and craniomaxillofacial tissue regenerative medicine. It first briefly introduces the structure of SF and then summarizes the principles, advantages and disadvantages of the different cross-linking methods (physical cross-linking, chemical cross-linking and double network structure) of SF. Finally, the existing research on the use of SF in tissue engineering and the prospects of using SF with different cross-linking methods in oral and craniomaxillofacial tissue regeneration are also discussed.

CONCLUSIONS

This review is intended to show the advantages of SF hydrogels in tissue engineering and provides theoretical support for establishing novel and viable silk protein hydrogels for regeneration.

Water for agriculture: more crop per drop

Global crop production in agriculture depends on water availability. Future scenarios predict increasing occurrence of flash floods and rapidly developing droughts accompanied by heatwaves in humid regions that rely on rain-fed agriculture. It is challenging to maintain high crop yields, even in arid and drought-prone regions that depend on irrigation. The average water demand of crops varies significantly, depending on plant species, development stage, and climate. Most crops, such as maize and wheat, require relatively more water during the vegetative phase compared to the ripening phase. In this review, we explain WUE and options to improve water use and thus crop yield. Nutrient management might represent another possibility to manipulate water uptake and use by plants. An emerging topic involves agroforest co-cultivation, where trees in the system facilitate water transfer through hydraulic lift, benefiting neighbouring crops. Other options to enhance crop yield per water use are discussed.

Helical reorganization in the context of membrane protein folding: Insights from simulations with bacteriorhodopsin (BR) fragments

Membrane protein folding is distinct from folding of soluble proteins. Conformational acquisition in major membrane protein subclasses can be delineated into insertion and folding processes. An exception to the "two stage" folding, later developed to "three stage" folding, is observed within the last two helices in bacteriorhodopsin (BR), a system that serves as a model membrane protein. We employ a reductionist approach to understand interplay of molecular factors underlying the apparent defiance. Leveraging available solution NMR structures, we construct, sample in silico, and analyze partially (PIn) and fully inserted (FIn) BR membrane states. The membrane lateral C-terminal helix (CH) in PIn is markedly prone to transient structural distortions over microsecond timescales; a disorder prone region (DPR) is thereby identified. While clear transmembrane propensities are not acquired, the distortions induce alterations in local membrane curvature and area per lipid. Importantly, energetic decompositions reveal that overall, the N-terminal helix (NH) is thermodynamically more stable in the PIn. Higher overall stability of the FIn arises from favorable interactions between the NH and the CH. Our results establish lack of spontaneous transition of the PIn to the FIn, and attributes their partitioning to barriers that exceed those accessible with thermal fluctuations. This work paves the way for further detailed studies aimed at determining the thermo-kinetic roles of the initial five helices, or complementary external factors, in complete helical folding and insertion in BR. We comment that complementing such efforts with the growing field of machine learning assisted energy landscape searches may offer unprecedented insights.

Rotational Dynamics of Water near Osmolytes by Molecular Dynamics Simulations

The behavior of water molecules around organic molecules has attracted considerable attention as a crucial factor influencing the properties and functions of soft matter and biomolecules. Recently, it has been suggested that the change in protein stability upon the addition of small organic molecules (osmolytes) is dominated by the change in the water dynamics caused by the osmolyte, where the dynamics of not only the directly interacting water molecules but also the long-range hydration layer affect the protein stability. However, the relation between the long-range structure of hydration water in various solutions and the water dynamics remains unclear at the molecular level. We performed density-functional tight-binding molecular dynamics simulations to elucidate the varying rotational dynamics of water molecules in 15 osmolyte solutions. A positive correlation was observed between the rotational relaxation time and our proposed normalized parameter obtained by dividing the number of hydrogen bonds between water molecules by the number of nearest-neighbor water molecules. For the 15 osmolyte solutions, an increase or a decrease in the value of the normalized parameter for the second hydration shell tended to result in water molecules with slow and fast rotational dynamics, respectively, thus illustrating the importance of the second hydration shell for the rotational dynamics of water molecules. Our simulation results are anticipated to advance the current understanding of water dynamics around organic molecules and the long-range structure of water molecules.

The Surface of Colloidal Metal-Organic Framework Nanoparticles Revealed by Vibrational Sum Frequency Scattering Spectroscopy

Solvation shells strongly influence the interfacial chemistry of colloidal systems, from the activity of proteins to the colloidal stability and catalysis of nanoparticles. Despite their fundamental and practical importance, solvation shells have remained largely undetected by spectroscopy. Furthermore, their ability to assemble at complex but realistic interfaces with heterogeneous and rough surfaces remains an open question. Here, we apply vibrational sum frequency scattering spectroscopy (VSFSS), an interface-specific technique, to colloidal nanocrystals with porous metal-organic frameworks (MOFs) as a case study. Due to the porous nature of the solvent-particle boundary, MOF particles challenge conventional models of colloidal and interfacial chemistry. Their multiweek colloidal stability in the absence of conventional surface ligands suggests that stability may arise in part from solvation forces. Spectra of colloidally stable Zn(2-methylimidazolate)2 (ZIF-8) in polar solvents indicate the presence of ordered solvation shells, solvent-metal binding, and spontaneous ordering of organic bridging linkers within the MOF. These findings help explain the unexpected colloidal stability of MOF colloids, while providing a roadmap for applying VSFSS to wide-ranging colloidal nanocrystals in general.

Fast collective motions of backbone in transmembrane α helices are critical to water transfer of aquaporin

Fast collective motions are widely present in biomolecules, but their functional relevance remains unclear. Herein, we reveal that fast collective motions of backbone are critical to the water transfer of aquaporin Z (AqpZ) by using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and molecular dynamics (MD) simulations. A total of 212 residue site-specific dipolar order parameters and 158 15N spin relaxation rates of the backbone are measured by combining the 13C- and 1H-detected multidimensional ssNMR spectra. Analysis of these experimental data by theoretic models suggests that the small-amplitude (~10°) collective motions of the transmembrane α helices on the nanosecond-to-microsecond timescales are dominant for the dynamics of AqpZ. The MD simulations demonstrate that these collective motions are critical to the water transfer efficiency of AqpZ by facilitating the opening of the channel and accelerating the water-residue hydrogen bonds renewing in the selectivity filter region.

Suppressed terahertz dynamics of water confined in nanometer gaps

Nanoconfined waters exhibit low static permittivity mainly due to interfacial effects that span about one nanometer. The characteristic length scale may be much longer in the terahertz (THz) regime where long-range collective dynamics occur; however, the THz dynamics have been largely unexplored because of the lack of a robust platform. Here, we use metallic loop nanogaps to sharply enhance light-matter interactions and precisely measure real and imaginary THz refractive indices of nanoconfined water at gap widths ranging from 2 to 20 nanometers, spanning mostly interfacial waters all the way to quasi-bulk waters. We find that, in addition to the well-known interfacial effect, the confinement effect also contributes substantially to the decrease in the complex refractive indices of the nanoconfined water by cutting off low-energy vibrational modes, even at gap widths as large as 10 nanometers. Our findings provide valuable insights into the collective dynamics of water molecules which is crucial to understanding water-mediated processes.

Water will Find Its Way: Transport through Narrow Tunnels in Hydrolases

An aqueous environment is vital for life as we know it, and water is essential for nearly all biochemical processes at the molecular level. Proteins utilize water molecules in various ways. Consequently, proteins must transport water molecules across their internal network of tunnels to reach the desired action sites, either within them or by functioning as molecular pipes to control cellular osmotic pressure. Despite water playing a crucial role in enzymatic activity and stability, its transport has been largely overlooked, with studies primarily focusing on water transport across membrane proteins. The transport of molecules through a protein's tunnel network is challenging to study experimentally, making molecular dynamics simulations the most popular approach for investigating such events. In this study, we focused on the transport of water molecules across three different α/β-hydrolases: haloalkane dehalogenase, epoxide hydrolase, and lipase. Using a 5 μs adaptive simulation per system, we observed that only a few tunnels were responsible for the majority of water transport in dehalogenase, in contrast to a higher diversity of tunnels in other enzymes. Interestingly, water molecules could traverse narrow tunnels with subangstrom bottlenecks, which is surprising given the commonly accepted water molecule radius of 1.4 Å. Our analysis of the transport events in such narrow tunnels revealed a markedly increased number of hydrogen bonds formed between the water molecules and protein, likely compensating for the steric penalty of the process. Overall, these commonly disregarded narrow tunnels accounted for ∼20% of the total water transport observed, emphasizing the need to surpass the standard geometrical limits on the functional tunnels to properly account for the relevant transport processes. Finally, we demonstrated how the obtained insights could be applied to explain the differences in a mutant of the human soluble epoxide hydrolase associated with a higher incidence of ischemic stroke.

Anomalously Large Heat Generation of Hydration Water under Microwave Irradiation

Much attention has been paid to the biological effects of microwave irradiation. The hydration water surrounding a biomolecule is crucial in its biological reactions and functions. Therefore, it is important to know the response of hydration water to microwaves to understand their biological effects; however, the scarcity of studies about it often leads to speculations and debates about that effect. In this study, we have made real-time temperature measurements of reverse micellar solutions with their water droplet size from ∼2.3 to ∼9.5 nm using a waveguide system combined with a microwave generator at 2.45 GHz. The heat generated by water in reverse micelles has been observed to depend on their size. It is about 10 times larger than that of liquid water at their small sizes (<∼3.5 nm) and diminishes with further enlarging the size, approaching the water's value at their large sizes (∼10 nm). These results indicate that the heat generation behavior has an interfacial effect; specifically, the hydration water on the surfactant layer produces heat 10 times larger than bulk water. Moreover, the hydration number per surfactant molecule decreases in a core-shell model with increasing the reverse micelle size. These features are also reflected in the heat generation rate. Our findings may offer a new and fundamental perspective for studies on the biological effects of microwave irradiation.

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the λ-exchange of free energy perturbation (FEP). Two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/λ-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.

Entropy Tug-of-War Determines Solvent Effects in the Liquid-Liquid Phase Separation of a Globular Protein

Liquid-liquid phase separation (LLPS) plays a key role in the compartmentalization of cells via the formation of biomolecular condensates. Here, we combined atomistic molecular dynamics (MD) simulations and terahertz (THz) spectroscopy to determine the solvent entropy contribution to the formation of condensates of the human eye lens protein γD-Crystallin. The MD simulations reveal an entropy tug-of-war between water molecules that are released from the protein droplets and those that are retained within the condensates, two categories of water molecules that were also assigned spectroscopically. A recently developed THz-calorimetry method enables quantitative comparison of the experimental and computational entropy changes of the released water molecules. The strong correlation mutually validates the two approaches and opens the way to a detailed atomic-level understanding of the different driving forces underlying the LLPS.

HydraProt: A New Deep Learning Tool for Fast and Accurate Prediction of Water Molecule Positions for Protein Structures

Water molecules are integral to the structural stability of proteins and vital for facilitating molecular interactions. However, accurately predicting their precise position around protein structures remains a significant challenge, making it a vibrant research area. In this paper, we introduce HydraProt (deep Hydration of Proteins), a novel methodology for predicting precise positions of water molecule oxygen atoms around protein structures, leveraging two interconnected deep learning architectures: a 3D U-net and a Multi-Layer Perceptron (MLP). Our approach starts by introducing a coarse voxel-based representation of the protein, which allows for rapid sampling of candidate water positions via the 3D U-net. These water positions are then assessed by embedding the water-protein relationship in the Euclidean space by means of an MLP. Finally, a postprocessing step is applied to further refine the MLP predictions. HydraProt surpasses existing state-of-the-art approaches in terms of precision and recall and has been validated on large data sets of protein structures. Notably, our method offers rapid inference runtime and should constitute the method of choice for protein structure studies and drug discovery applications. Our pretrained models, data, and the source code required to reproduce these results are accessible at https://github.com/azamanos/HydraProt.

Molecular Structure of Salicylic Acid and Its Hydrates: A Rotational Spectroscopy Study

We present a study of salicylic acid and its hydrates, with up to four water molecules, done by employing chirped-pulse Fourier transform microwave spectroscopy. We employed the spectral data set of the parent, 13C, and 2H isotopologues to determine the molecular structure and characterize the intra- and intermolecular interactions of salicylic acid and its monohydrate. Complementary theoretical calculations were done to support the analysis of the experimental results. For the monomer, we analyzed structural properties, such as the angular-group-induced bond alternation (AGIBA) effect. In the microsolvates, we analyzed their main structural features dominated by the interaction of water with the carboxylic acid group. This work contributes to seeding information on how water molecules accumulate around this group. Moreover, we discussed the role of cooperative effects further stabilizing the observed inter- and intramolecular hydrogen bond interactions.

Engineering diaryl alcohol dehydrogenase KpADH reveals importance of retaining hydration shell in organic solvent tolerance

Alcohol dehydrogenases (ADHs) are synthetically important biocatalysts for the asymmetric synthesis of chiral alcohols. The catalytic performance of ADHs in the presence of organic solvents is often important since most prochiral ketones are highly hydrophobic. Here, the organic solvent tolerance of KpADH from Kluyveromyces polyspora was semi-rationally evolved. Using tolerant variants obtained, meticulous experiments and computational studies were conducted to explore properties including stability, activity and kinetics in the presence of various organic solvents. Compared with WT, variant V231D exhibited 1.9-fold improvement in ethanol tolerance, while S237G showed a 6-fold increase in catalytic efficiency, a higherT 50 15 $$ {\mathrm{T}}_{50}^{15} $$ , as well as 15% higher tolerance in 7.5% (v/v) ethanol. Based on 3 × 100 ns MD simulations, the increased tolerance of V231D and S237G against ethanol may be ascribed to their enhanced ability in retaining water molecules and repelling ethanol molecules. Moreover, 6.3-fold decreased KM value of V231D toward hydrophilic ketone substrate confirmed its capability of retaining hydration shell. Our results suggest that retaining hydration shell surrounding KpADH is critical for its tolerance to organic solvents, as well as catalytic performance. This study provides useful guidance for engineering organic solvent tolerance of KpADH and other ADHs.

Sulfation of hyaluronic acid reconfigures the mechanistic pathway of bone morphogenetic protein-2 aggregation

Bone morphogenetic protein-2 (BMP-2) is a critical growth factor of bone extracellular matrix (ECM), pivotal for osteogenesis. Glycosaminoglycans (GAGs), another vital ECM biomolecules, interact with growth factors, affecting signal transduction. Our study primarily focused on hyaluronic acid (HA), a prevalent GAG, and its sulfated derivative (SHA). We explored their impact on BMP-2's conformation, aggregation, and mechanistic pathways of aggregation using diverse optical and rheological methods. In the presence of HA and SHA, the secondary structure of BMP-2 underwent a structured transformation, characterized by a substantial increase in beta sheet content, and a detrimental alteration, manifesting as a shift towards unstructured content, respectively. Although both HA and SHA induced BMP-2 aggregation, their mechanisms differed. SHA led to rapid amorphous aggregates, while HA promoted amyloid fibrils with a lag phase and sigmoidal kinetics. Aggregate size and shape varied; HA produced larger structures, SHA smaller. Each aggregation type followed distinct pathways influenced by viscosity and excluded volume. Higher viscosity, low diffusivity of protein and higher excluded volume In the presence of HA promotes fibrillation having size in micrometer range. Low viscosity, high diffusivity of protein and lesser excluded volume leads to amorphous aggregate of size in nanometer range.

Unraveling Water Solvation Effects with Quantum Mechanics/Molecular Mechanics Semiclassical Vibrational Spectroscopy: The Case of Thymidine

We introduce a quantum mechanics/molecular mechanics semiclassical method for studying the solvation process of molecules in water at the nuclear quantum mechanical level with atomistic detail. We employ it in vibrational spectroscopy calculations because this is a tool that is very sensitive to the molecular environment. Specifically, we look at the vibrational spectroscopy of thymidine in liquid water. We find that the C═O frequency red shift and the C═C frequency blue shift, experienced by thymidyne upon solvation, are mainly due to reciprocal polarization effects, that the molecule and the water solvent exert on each other, and nuclear zero-point energy effects. In general, this work provides an accurate and practical tool to study quantum vibrational spectroscopy in solution and condensed phase, incorporating high-level and computationally affordable descriptions of both electronic and nuclear problems.

Computational insights into dynamics and conformational stability of N-acetylmannosamine kinase mutations

The activity of UDP-GlcNAc 2-epimerase/ManNAc kinase (GNE) is essential for the biosynthesis of sialic acid, which is involved in cellular processes in health and diseases. GNE contains an N-terminal epimerase domain and a C-terminal kinase domain (N-acetylmannosamine kinase, MNK). Mutations of the GNE protein led to hypoactivity of the enzyme and cause sialurea or autosomal recessive inclusion body myopathy/Nonaka myopathy. Here, we used all-atom molecular dynamics (MD) simulations to comprehend the folding, dynamics and conformational stability of MNK variants, including the wild type (WT) and three mutants (H677R, V696M and H677R/V696M). The deleterious and destabilizing nature of MNK mutants were predicted using different prediction tools. Results predicted that mutations modulate the stability, flexibility and function of MNK. The effect of mutations on the conformational stability and dynamics of MNK was next studied through the free-energy landscape (FEL), hydrogen-bonds and secondary structure changes. The FEL results show that the mutations interfere with various conformational transitions in both WT and mutants, exposing the structural underpinnings of protein destabilization and unfolding brought on by mutation. We discover that, when compared to the other two mutations, V696M and H677R/V696M, H677R has the most harmful effects. These findings have a strong correlation with published experimental studies that demonstrate how these mutations disrupt MNK activity. Hence, this computational study describes the structural details to unravel the mutant effects at the atomistic resolution and has implications for understanding the GNE's physiological and pathological role.Communicated by Ramaswamy H. Sarma.

Conformational energies of biomolecules in solution: Extending the MPCONF196 benchmark with explicit water molecules

A prerequisite for the computational prediction of molecular properties like conformational energies of biomolecules is a reliable, robust, and computationally affordable method usually selected according to its performance for relevant benchmark sets. However, most of these sets comprise molecules in the gas phase and do not cover interactions with a solvent, even though biomolecules typically occur in aqueous solution. To address this issue, we introduce a with explicit water molecules solvated version of a gas-phase benchmark set containing 196 conformers of 13 peptides and other relevant macrocycles, namely MPCONF196 [J. Řezáč et al., JCTC 2018, 14, 1254-1266], and provide very accurate PNO-LCCSD(T)-F12b/AVQZ' reference values. The novel solvMPCONF196 benchmark set features two additional challenges beyond the description of conformers in the gas phase: conformer-water and water-water interactions. The overall best performing method for this set is the double hybrid revDSDPBEP86-D4/def2-QZVPP yielding conformational energies of almost coupled cluster quality. Furthermore, some (meta-)GGAs and hybrid functionals like B97M-V and ω B97M-D with a large basis set reproduce the coupled cluster reference with an MAD below 1 kcal mol- 1  . If more efficient methods are required, the composite DFT-method r2 SCAN-3c (MAD of 1.2 kcal mol- 1 ) is a good alternative, and when conformational energies of polypeptides or macrocycles with more than 500-1000 atoms are in the focus, the semi-empirical GFN2-xTB or the MMFF94 force field (for very large systems) are recommended.

Deciphering Density Fluctuations in the Hydration Water of Brownian Nanoparticles via Upconversion Thermometry

We investigate the intricate relationship among temperature, pH, and Brownian velocity in a range of differently sized upconversion nanoparticles (UCNPs) dispersed in water. These UCNPs, acting as nanorulers, offer insights into assessing the relative proportion of high-density and low-density liquid in the surrounding hydration water. The study reveals a size-dependent reduction in the onset temperature of liquid-water fluctuations, indicating an augmented presence of high-density liquid domains at the nanoparticle surfaces. The observed upper-temperature threshold is consistent with a hypothetical phase diagram of water, validating the two-state model. Moreover, an increase in pH disrupts the organization of water molecules, similar to external pressure effects, allowing simulation of the effects of temperature and pressure on hydrogen bonding networks. The findings underscore the significance of the surface of suspended nanoparticles for understanding high- to low-density liquid fluctuations and water behavior at charged interfaces.

Cholesterol modulates the structural dynamics of the paddle motif loop of KvAP voltage sensor

KvAP is a prokaryotic Kv channel, which has been widely used as a model system to understand voltage- and lipid-dependent gating mechanisms. In phospholipid membranes, the KvAP-VSD adopts the activated/'Up' conformation, whereas the presence of non-phospholipids in membranes favours the structural transition to resting/'Down' state. The S3b-S4 paddle motif loop of KvAP-VSD is functionally important as this participates in protein-protein interactions and is the target for animal toxins. In this study, we have monitored the modulatory role of cholesterol - the physiologically-relevant non-phospholipid - on the organization and dynamics of the S3b-S4 loop of the isolated KvAP-VSD in membranes by site-directed fluorescence approaches using the environmental sensitivity of 7-nitrobenz-2-oxa-1,3-diazol-4-yl-ethylenediamine (NBD) fluorescence. Our results show that cholesterol alters the dynamic nature (rotational and hydration dynamics) of S3b-S4 loop in a segmental fashion, i.e., the residues 110 to 114 and 115 to 117 behave differently in the presence of cholesterol, which is accompanied by considerable change in conformational heterogeneity. Further, quantitative depth measurements using the parallax quenching method reveal that the sensor loop is located at the shallow interfacial region of cholesterol-containing membranes, suggesting that the sensor loop organization is not directly correlated with S4 helix movement. Our results clearly show that cholesterol-induced changes in bilayer properties may not be the predominant factor for the sensor loop's altered structural dynamics, but can be attributed to the conformational change of the KvAP-VSD in cholesterol-containing membranes. Overall, these results are relevant for gating mechanisms, particularly the lipid-dependent gating, of Kv channels in membranes.

Solvent Suppression in Pure Shift NMR

Intense solvent signals in 1H solution-state NMR experiments typically cause severe distortion of spectra and mask nearby solute signals. It is often infeasible or undesirable to replace a solvent with its perdeuterated form, for example, when analyzing formulations in situ, when exchangeable protons are present, or for practical reasons. Solvent signal suppression techniques are therefore required. WATERGATE methods are well-known to provide good solvent suppression while enabling retention of signals undergoing chemical exchange with the solvent signal. Spectra of mixtures, such as pharmaceutical formulations, are often complicated by signal overlap, high dynamic range, the narrow spectral width of 1H NMR, and signal multiplicity. Here, we show that by combining WATERGATE solvent suppression with pure shift NMR, ultrahigh-resolution 1H NMR spectra can be acquired while suppressing intense solvent signals and retaining exchangeable 1H signals. The new method is demonstrated in the analysis of cyanocobalamin, a vitamin B12 supplement, and of an eye-drop formulation of atropine.

Diffusive dynamics of a model protein chain in solution

A Markov state model is a powerful tool that can be used to track the evolution of populations of configurations in an atomistic representation of a protein. For a coarse-grained linear chain model with discontinuous interactions, the transition rates among states that appear in the Markov model when the monomer dynamics is diffusive can be determined by computing the relative entropy of states and their mean first passage times, quantities that are unchanged by the specification of the energies of the relevant states. In this paper, we verify the folding dynamics described by a diffusive linear chain model of the crambin protein in three distinct solvent systems, each differing in complexity: a hard-sphere solvent, a solvent undergoing multi-particle collision dynamics, and an implicit solvent model. The predicted transition rates among configurations agree quantitatively with those observed in explicit molecular dynamics simulations for all three solvent models. These results suggest that the local monomer-monomer interactions provide sufficient friction for the monomer dynamics to be diffusive on timescales relevant to changes in conformation. Factors such as structural ordering and dynamic hydrodynamic effects appear to have minimal influence on transition rates within the studied solvent densities.

Hydrogen bond energy estimation (H-BEE) in large molecular clusters: A Python program for quantum chemical investigations

A procedure, derived from the fragmentation-based molecular tailoring approach (MTA), has been proposed and extensively applied by Deshmukh and Gadre for directly estimating the individual hydrogen bond (HB) energies and cooperativity contributions in molecular clusters. However, the manual fragmentation and high computational cost of correlated quantum chemical methods make the application of this method to large molecular clusters quite formidable. In this article, we report an in-house developed software for automated hydrogen bond energy estimation (H-BEE) in large molecular clusters. This user-friendly software is essentially written in Python and executed on a Linux platform with the Gaussian package at the backend. Two approximations to the MTA-based procedure, viz. the first spherical shell (SS1) and the Fragments-in-Fragments (Frags-in-Frags), enabling cost-effective, automated evaluation of HB energies and cooperativity contributions, are also implemented in this software. The software has been extensively tested on a variety of molecular clusters and is expected to be of immense use, especially in conjunction with correlated methods such as MP2, CCSD(T), and so forth.

Olive Oil-Based Reverse Microemulsion for Stability and Topical Delivery of Methotrexate: In Vitro

Hydrolysis of pharmaceutically active molecules can be in control under a confined environment of water-in-oil microemulsion. Stability of model drug methotrexate (MTX) in a sodium bis(2-ethylhexyl) sulfosuccinate (AOT) and olive oil microemulsion system has been evaluated. The physicochemical properties of AOT-MTX-water-olive oil reverse microemulsion (MTX-RM) were examined by UV-vis, Fourier transform infrared, and X-ray diffraction techniques, and the hydrodynamic size was determined by dynamic light scattering techniques and morphologies were characterized by a transmission electron microscope and atomic force microscope. In vitro permeation of MTX-RM through treated skin and its mechanism are evaluated by a UV-visible spectrophotometer, confocal laser scanning microscope, differential scanning calorimeter, and attenuated total reflecting infrared spectroscopy (ATR). The interaction of MTX with the AOT headgroup in confined environment RM enhanced the stability of MTX without affecting the molecular integrity at room temperature. Chemical stability of MTX in MTX-RM (W0 = 5) is significantly higher (∼97%) at room temperature for the study period of 1 year than in MTX-RM (W0 = 15) (∼72%). Interaction of MTX with the AOT headgroup is also visualized by a high-resolution transmission electron microscope and is in correlation with FT-IR data of MTX-RM. The skin fluxes of MTX are 15.1, 19.75, and 22.75 times higher at water content (W0) of 5, 10, and 15, respectively, in MTX-RM in comparison to aqueous solution of MTX. The enhanced amounts of the MTX were detected using CLSM in hair follicles, sweat glands, and epidermis layer of the skin. Merging of T2, T3, and T4 thermal peaks in one broad peak in treated skin endothermograph shows that carrier MTX-RM affects the lipid as well protein structure of the treated skin. ATR data of treated skin showed an increase in the intensity of the carbonyl peak at 1750 cm-1 (lipid), shifting of the amide II peaks, and separation of peaks in the range of 1060 to 1000 cm-1 (vibration mode of -CH2OH, C-O stretching, and C-OH bending peak of the carbohydrate) in comparison to control skin, which indicates that MTX-RM interacts with glycolipid and glycoprotein through carbohydrate hydroxy groups.

Improved Protein Dynamics and Hydration in the Martini3 Coarse-Grain Model

The Martini coarse-grain force-field has emerged as an important framework to probe cellular processes at experimentally relevant time- and length-scales. However, the recently developed version, the Martini3 force-field with the implemented Go̅ model (Martini3Go̅), as well as previous variants of the Martini model have not been benchmarked and rigorously tested for globular proteins. In this study, we consider three globular proteins, ubiquitin, lysozyme, and cofilin, and compare protein dynamics and hydration with observables from experiments and all-atom simulations. We show that the Martini3Go̅ model is able to accurately model the structural and dynamic features of small globular proteins. Overall, the structural integrity of the proteins is maintained, as validated by contact maps, radii of gyration (Rg), and SAXS profiles. The chemical shifts predicted from the ensemble sampled in the simulations are consistent with the experimental data. Further, a good match is observed in the protein-water interaction energetics, and the hydration levels of the residues are similar to atomistic simulations. However, the protein-water interaction dynamics is not accurately represented and appears to depend on the protein structural complexity, residue specificity, and water dynamics. Our work is a step toward testing and assessing the Martini3Go̅ model and provides insights into future efforts to refine Martini models with improved solvation effects and better correspondence to the underlying all-atom systems.

Atomic-scale structure of interfacial water on gel and liquid phase lipid membranes

Hydration of biological membranes is essential to a wide range of biological processes. In particular, it is intrinsically linked to lipid thermodynamic properties, which in turn influence key cell functions such as ion permeation and protein mobility. Experimental and theoretical studies of the surface of biomembranes have revealed the presence of an interfacial repulsive force, which has been linked to hydration or steric effects. Here, we directly characterise the atomic-scale structure of water near supported lipid membranes of 1,2-dimyristoyl-sn-glycero-3-phosphocholine in their gel and liquid phase through three-dimensional atomic force microscopy (3D AFM). First, we demonstrate the ability to probe the morphology of interfacial water of lipid bilayers in both phases with sub-molecular resolution by using ultrasharp tips. We then visualise the molecular arrangement of water at the lipid surface at different temperatures. Our experiments reveal that water is organised in multiple hydration layers on both the solid-ordered and liquid-disordered lipid phases. Furthermore, we observe a monotonic repulsive force, which becomes relevant only in the liquid phase. These results offer new insights into the water structuring near soft biological surfaces, and demonstrate the importance of investigating it with vertical and lateral sub-molecular resolution.

Subdomain dynamics enable chemical chain reactions in non-ribosomal peptide synthetases

Many peptide-derived natural products are produced by non-ribosomal peptide synthetases (NRPSs) in an assembly-line fashion. Each amino acid is coupled to a designated peptidyl carrier protein (PCP) through two distinct reactions catalysed sequentially by the single active site of the adenylation domain (A-domain). Accumulating evidence suggests that large-amplitude structural changes occur in different NRPS states; yet how these molecular machines orchestrate such biochemical sequences has remained elusive. Here, using single-molecule Förster resonance energy transfer, we show that the A-domain of gramicidin S synthetase I adopts structurally extended and functionally obligatory conformations for alternating between adenylation and thioester-formation structures during enzymatic cycles. Complementary biochemical, computational and small-angle X-ray scattering studies reveal interconversion among these three conformations as intrinsic and hierarchical where intra-A-domain organizations propagate to remodel inter-A-PCP didomain configurations during catalysis. The tight kinetic coupling between structural transitions and enzymatic transformations is quantified, and how the gramicidin S synthetase I A-domain utilizes its inherent conformational dynamics to drive directional biosynthesis with a flexibly linked PCP domain is revealed.

Martini on the Rocks: Can a Coarse-Grained Force Field Model Crystals?

Computational chemistry is an important tool in numerous scientific disciplines, including drug discovery and structural biology. Coarse-grained models offer simple representations of molecular systems that enable simulations of large-scale systems. Because there has been an increase in the adoption of such models for simulations of biomolecular systems, critical evaluation is warranted. Here, the stability of the amyloid peptide and organic crystals is evaluated using the Martini 3 coarse-grained force field. The crystals change shape drastically during the simulations. Radial distribution functions show that the distance between backbone beads in β-sheets increases by ∼1 Å, breaking the crystals. The melting points of organic compounds are much too low in the Martini force field. This suggests that Martini 3 lacks the specific interactions needed to accurately simulate peptides or organic crystals without imposing artificial restraints. The problems may be exacerbated by the use of the 12-6 potential, suggesting that a softer potential could improve this model for crystal simulations.

Energy-Scaled Debye-Hückel Theory for the Electrostatic Solvation Free Energy in Size-Asymmetric Electrolyte Solutions

In this report, an energy-scaled Debye-Hückel theory is developed for fast and accurate evaluation of the electrostatic solvation free energy in size-asymmetric electrolyte solutions. A size-asymmetric electrolyte solution is mapped to a dielectric continuum medium with Debye-Hückel-like response. Based on the scaling relation of the electrostatic energy of a spherical ion in the small and large size limits, a Padé polynomial is used to interpolate the electrostatic energy at finite size. The Padé polynomial is further interpreted as the electrostatic energy of an effective Debye-Hückel mean field model, depicted by a modified Debye parameter and a surface charge density due to the size asymmetry of the solvent ions. This theory can distinguish the electrostatic energies and the electrostatic solvation free energies of solutes with the same size but opposite charges. Application to charged hard and charged soft spheres demonstrates the accuracy of our approach.

Lysozyme aggregation and unfolding in ionic liquid solvents: Insights from small angle X-ray scattering and high throughput screening

Understanding protein behaviour is crucial for developing functional solvent systems. Ionic liquids (ILs) are designer salts with versatile ion combinations, where some suppress unfavourable protein behaviour. This work utilizes small angle X-ray scattering (SAXS) to investigate the size and shape changes of model protein hen egg white lysozyme (HEWL) in 137 IL and salt solutions. Guinier, Kratky, and pair distance distribution analysis were used to evaluate the protein size, shape, and aggregation changes in these solvents. At low IL and salt concentration (1 mol%), HEWL remained monodispersed and globular. Most ILs increased HEWL size compared to buffer, while the nitrate and mesylate anions induced the most significant size increases. IL cation branching, hydroxyl groups, and longer alkyl chains counteracted this size increase. Common salts exhibited specific ion effects, while the IL effect varied with concentration due to complex ion-pairing. Protein aggregation and unfolding occurred at 10 mol% IL, altering the protein shape, especially for ILs with multiple alkyl chains on the cation, or with a mesylate/nitrate anion. This study highlights the usefulness of adopting a high-throughput SAXS strategy for understanding IL effects on protein behaviour and provides insights on controlling protein aggregation and unfolding with ILs.

Hydrodynamic Control of Alzheimer Aβ Fibrillation with Glucosaminic Acid Containing Click-Cyclized β-Bodies

It is well established that the dynamic hydration shell plays a vital role in macromolecular functions such as protein-ligand, protein-protein, protein-DNA, and protein-lipid interactions. Here we investigate how the water modality affects conformational changes, solubility, and motion of fibrillar proteins. The hypothesis is that the introduction of a poly hydroxyl amino acid would increase solvation of the fibril forming peptides, preventing their misfolding and aggregation. For the amyloid β (Aβ) peptide, which is considered to be connected with nervous system diseases, including dementia and cognitive decline in Alzheimer's disease, the formation of β-sheet fibrils always occurs with a conformational change and a decrease in the dynamic hydration shell around Aβ(1-42). We present novel cyclic d-amino acid peptides that effectively inhibit fibrillation through affecting the dynamic hydration shell of Aβ(1-42) in vitro. Using de novo design within the software Molecular Operating Environment (MOE), five different peptides that recognize Alzheimer's fibrils were designed and synthesized. Three of them were cyclic all-d-amino acid peptides incorporating the same polyhydroxy building block derived from d-glucosaminic acid (GA). One peptide was the parent cyclic all d-amino acid inhibitor with no GA incorporated, and another was an all l-amino acid linear fibrillation inhibitor. The GA-containing peptides were found to show significantly improved inhibition of Aβ(1-42) aggregation. The inhibition was dramatically improved by the synergistic application of two GA peptides targeting each end of the growing fibril. The present study may facilitate future developments of intervention strategies for Alzheimer's disease and similar neurodegenerative diseases.

Chiral Spinodal-like Ordering of Homoimmiscible Water at Interface between Water and Chiral Ice III

Diversity in structures of water endowed by a hydrogen-bonding network plays crucial roles in wide varieties of phenomena in nature. Chiral ordering of water molecules is an intriguing phenomenon from the viewpoint of bimolecular functions. However, experimental reports on chiral ordering have been limited to the water molecules interacting with biomolecules on the molecular scale. It remains unclear whether pure liquid water forms long-range chiral ordering without any interaction with biomolecules. Here, we show that chiral anisotropy can be observed in the macro/mesoscopic network pattern of an unknown water layer formed via spinodal phase separation-like dynamics at the interface between water and ice III with a chiral crystal structure. We named this unknown water homoimmiscible water. Our observations infer that the unknown water is a chiral liquid crystal. This possibility opens new avenues for a wide variety of research fields such as liquid polymorphism, biology, earth and planetary science, and so forth from the perspective of chirality.

Water inside the Selectivity Filter of a K+ Ion Channel: Structural Heterogeneity, Picosecond Dynamics, and Hydrogen Bonding

Water inside biological ion channels regulates the key properties of these proteins, such as selectivity, ion conductance, and gating. In this article, we measure the picosecond spectral diffusion of amide I vibrations of an isotope-labeled KcsA potassium channel using two-dimensional infrared (2D IR) spectroscopy. By combining waiting time (100-2000 fs) 2D IR measurements of the KcsA channel including 13C18O isotope-labeled Val76 and Gly77 residues with molecular dynamics simulations, we elucidated the site-specific dynamics of water and K+ ions inside the selectivity filter of KcsA. We observe inhomogeneous 2D line shapes with extremely slow spectral diffusion. Our simulations quantitatively reproduce the experiments and show that water is the only component with any appreciable dynamics, whereas K+ ions and the protein are essentially static on a picosecond timescale. By analyzing simulated and experimental vibrational frequencies, we find that water in the selectivity filter can be oriented to form hydrogen bonds with adjacent or nonadjacent carbonyl groups with the reorientation timescales being three times slower and comparable to that of water molecules in liquid, respectively. Water molecules can reside in the cavity sufficiently far from carbonyls and behave essentially like "free" gas-phase-like water with fast reorientation times. Remarkably, no interconversion between these configurations was observed on a picosecond timescale. These dynamics are in stark contrast with liquid water, which remains highly dynamic even in the presence of ions at high concentrations.

Unveiling the Role of Water on π-π Stacking Through Microwave Spectroscopy of (Thiophene)2-(Water)1-2 Clusters

There has been an ongoing debate about whether water enhances or hinders π-π stacking, a phenomenon crucial in various biological and chemical systems. In this study, the influence of water on π-π stacking is investigated by microwave spectroscopic observation of gas-phase hydrated clusters of thiophene dimers. Two isomers of (C4H4S)2-H2O and two isomers of (C4H4S)2-(H2O)2 have been unambiguously identified. These identifications are supported by quantum chemistry calculations and isotopic measurements. In each of these conformations, water molecules are situated between aromatic pairs, forming distinctive interactions. Water molecules engage with thiophene molecules either as hydrogen bond donors through OH···π interactions or as hydrogen bond acceptors through CH···O interactions. The energy decomposition analysis indicates that the bonding pattern of water molecules significantly affects the π···π interactions between aromatic rings. These findings offer valuable structural insights into the role of water in shaping π-π stacking.

Quantifying the energy landscape in weakly and strongly disordered frictional media

We investigate the "roughness" of the energy landscape of a system that diffuses in a heterogeneous medium with a random position-dependent friction coefficient α(x). This random friction acting on the system stems from spatial inhomogeneity in the surrounding medium and is modeled using the generalized Caldira-Leggett model. For a weakly disordered medium exhibiting a Gaussian random diffusivity D(x) = kBT/α(x) characterized by its average value ⟨D(x)⟩ and a pair-correlation function ⟨D(x1)D(x2)⟩, we find that the renormalized intrinsic diffusion coefficient is lower than the average one due to the fluctuations in diffusivity. The induced weak internal friction leads to increased roughness in the energy landscape. When applying this idea to diffusive motion in liquid water, the dissociation energy for a hydrogen bond gradually approaches experimental findings as fluctuation parameters increase. Conversely, for a strongly disordered medium (i.e., ultrafast-folding proteins), the energy landscape ranges from a few to a few kcal/mol, depending on the strength of the disorder. By fitting protein folding dynamics to the escape process from a metastable potential, the decreased escape rate conceptualizes the role of strong internal friction. Studying the energy landscape in complex systems is helpful because it has implications for the dynamics of biological, soft, and active matter systems.

Thermoresponsive Polymers under Pressure with a Focus on Poly(N-isopropylacrylamide) (PNIPAM)

Pressure is a key variable in the phase behavior of responsive polymers, both for applications and from a fundamental point of view. In this feature article, we review recent developments, particularly applications of neutron techniques such as small-angle neutron scattering (SANS) and quasi-elastic neutron scattering (QENS), across the temperature-pressure phase diagram. These are complemented by kinetic SANS experiments following pressure jumps. In the prototype system poly(N-isopropylacrylamide) (PNIPAM), QENS revealed the pressure-dependent characteristics of hydration water around the lower critical solution temperature transition. The size, water content, and inner structure of the mesoglobules formed in the two-phase region depend strongly on pressure, as shown by SANS. Beside these changes at the phase transition, the mesoglobule formation at low pressure is determined by kinetic factors, namely the formation of a polymer-rich, rigid shell, which hampers further growth by coalescence. At high pressure, in contrast, the growth proceeds by diffusion-limited coalescence without any kinetic hindrance. The disintegration of the mesoglobules evolves either via chain release from their surface or via swelling, depending on the osmotic pressure of the water. Moreover, we report on the profound influence of pressure on the cononsolvency effect. In the temperature-pressure frame, the one-phase region is hugely expanded upon the addition of the cosolvent methanol. SANS experiments unveil the enthalpic and entropic contributions to the effective Flory-Huggins interaction parameter between the segments and the solvent mixture. QENS experiments demonstrate an increase in polymer associated water with pressure, whereas methanol is released. Correspondingly, the solvent phase becomes enriched in methanol, providing a mechanism for the breakdown of cononsolvency at a high pressure. Finally, we outline future opportunities for high-pressure studies of thermoresponsive polymers, with a focus on neutron methods.

Alterations in Lens Free Water Distribution Are Associated with Shape Deformation in Accommodation

Objective

To investigate whether a redistribution of water within the crystalline lens is associated with the shape deformation that occurs during accommodation.

Design

Observational, cross sectional study.

Subjects

Eleven young adults without presbyopia (aged 18-39 years) and 9 middle-aged adults with presbyopia (aged 40-55 years).

Methods

Magnetic resonance imaging (MRI) scans of the lens were acquired on a 3 Tesla clinical MRI scanner, without and with the presentation of a 3 Diopter accommodative stimulus. The MRIs were postprocessed using established methods to extract the geometric dimensions and spatial maps of water distribution of the lens.

Main Outcome Measures

Accommodative changes in the full 3-dimensional description of lens shape, the lens total-water distribution profile, and the lens free-water distribution profile.

Results

Viewing of an accommodative stimulus by young subjects elicited an elastic shape deformation of the lens consistent with accommodation that was associated with an elevated, smoother free-water distribution, primarily in the anterior region of the lens. In contrast, viewing of an accommodative stimulus by presbyopic subjects produced an atypical shape deformation of the lens that was instead associated with a lowered free-water distribution, primarily in the anterior region of the lens. No discernible changes to the lens total-water distribution were observed in response to the accommodative stimulus in either subject cohort.

Conclusions

The present study suggests that protein-mediated alterations in the free-water distribution of the anterior region of the lens influence the shape deformation in accommodation, presenting pharmacological modulation of free-water distribution as an attractive novel approach for treating presbyopia.

Financial Disclosures

The authors have no proprietary or commercial interest in any materials discussed in this article.

Water network chemistry to exploit the nature of catalytic water molecules in Mtb DNA gyrase: a computational study to understand the binding mechanism of fluoroquinolones

The dynamics of DNA gyrase and mutants of DNA gyrA such as G88A, A90V, S91P, D94A, D94G, D94N, D94Y; and double-point mutant (S91P-D94G), are meticulously investigated using computational approaches. Molecular dynamics (MD) and hydration thermodynamics have shed light on the fundamental, mechanistic basis of mutations on the conformational stability of Quinolone Binding Pocket (QBP) of DNA gyrase. Analysis of MD results revealed the displacement of a single crystal water molecule (HOH201) from the catalytic site of wild-type (WT) and mutants of DNA gyrA. This prompted our research group to probe the five crystal water molecules present in the QBP of the enzyme using water thermodynamics. Hydration thermodynamics analysis revealed the displacement of HOH201 due to unstable thermodynamic signatures. Further, the analysis highlighted significant changes in thermodynamic signatures and locations of five crystal water hydration sites upon mutation. Integrated MD simulations and water thermodynamics provided promising insights into the conformational changes and inaccessibility of the catalytic water molecule that can influence the design of DNA gyrase inhibitors.Communicated by Ramaswamy H. Sarma.

Single-molecule analysis of osmolyte-mediated nanomechanical unfolding behavior of a protein domain

The small organic molecules, known as osmolytes being ubiquitously present in different cell types, affect protein folding, stability and aggregation. However, it is unknown how the osmolytes affect the nanomechanical unfolding behavior of protein domain. Here, we show the osmolyte-dependent mechanical unfolding properties of protein titin immunoglobulin-27 (I27) domain using an atomic force microscopy (AFM)-based single-molecule force spectroscopy. We found that amines and methylamines improved the mechanical stability of I27 domain, whereas polyols had no effect. Interestingly, glycine betaine (GB) or trimethylamine-N-oxide (TMAO) increased the average unfolding force of the protein domain. The kinetic parameters analyzed at single-molecule level reveal that stabilizing effect of osmolytes is due to a decrease in the unfolding rate constant of I27, which was confirmed by molecular dynamics simulations. Our study reveals different effects that diverse osmolytes have on the mechanical properties of the protein, and suggests the potential use of osmolytes in modulating the mechanical stability of proteins required for various nano-biotechnological applications.

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.

Changes in the structure and hydration properties of high-temperature peanut protein induced by cold plasma oxidation

To improve the hydration properties of high-temperature pressed peanut protein isolate (HPPI), we investigated the effect of cold plasma (CP) oxidation on functional and structural properties. Compared to HPPI, the hydrated molecules number and the surface contact angle were significantly decreased at 70 W, from 77.2 × 109 to 17.7 × 109 and from 85.74° to 57.81°, respectively. The reduction of the sulfhydryl content and the increase of the disulfide bond and di-tyrosine content indicated that the structural transformation was affected by the oxidation effect. In terms of structural changes, a stretched tertiary structure, ordered secondary structure, and rough apparent structure were observed after CP treatment. Additionally, the enhancement of surface free energy and group content such as -COOH, -CO and -OH on the surface of HPPI contributed to the formation of hydrated crystal structures. In general, the oxidation effect of CP effectively improved the hydration properties of HPPI and broaden its application field.

Markov State Model of Solvent Features Reveals Water Dynamics in Protein-Peptide Binding

In this work, we investigate the role of solvent in the binding reaction of the p53 transactivation domain (TAD) peptide to its receptor MDM2. Previously, our group generated 831 μs of explicit-solvent aggregate molecular simulation trajectory data for the MDM2-p53 peptide binding reaction using large-scale distributed computing and subsequently built a Markov State Model (MSM) of the binding reaction (Zhou et al. 2017). Here, we perform a tICA analysis and construct an MSM with similar hyperparameters while using only solvent-based structural features. We find a remarkably similar landscape but accelerated implied timescales for the slowest motions. The solvent shells contributing most to the first tICA eigenvector are those centered on Lys24 and Thr18 of the p53 TAD peptide in the range of 3-6 Å. Important solvent shells were visualized to reveal solvation and desolvation transitions along the peptide-protein binding trajectories. Our results provide a solvent-centric view of the hydrophobic effect in action for a realistic peptide-protein binding scenario.

Ultrafast Spectroscopy Reveals Slow Water Dynamics in Biocondensates

Cells achieve high spatiotemporal control over biochemical processes through compartmentalization to membrane-bound as well as membraneless organelles that assemble by liquid-liquid phase separation. Characterizing the balance of forces within these environments is essential to understanding their stability and function, and water is an integral part of the condensate, playing an important role in mediating electrostatic and hydrogen-bonding interactions. Here, we investigate the ultrafast, picosecond hydrogen-bond dynamics of a model biocondensate consisting of a peptide poly-l-arginine (Poly-R) and the nucleic acid adenosine monophosphate (AMP) using coherent two-dimensional infrared (2D IR) spectroscopy. We investigated three vibrational modes: the arginine side-chain C═N stretches, an AMP ring mode, and the amide backbone carbonyl stretching modes. Dynamics slow considerably between the dilute phase and the condensate phase for each vibrational probe. For example, the arginine side-chain C═N modes slow from 0.38 to 2.26 ps due to strong electrostatic interactions. All-atom molecular dynamics simulations provide an atomistic interpretation of the H-bond network disruption resulting from electrostatic contributions as well as collapse within the condensate. Simulations predict that a fraction of water molecules are highly constrained within the condensate, explaining the observed slowdown in the H-bond dynamics.

Determination of the Thickness of Interfacial Water by Time-Resolved Sum-Frequency Generation Vibrational Spectroscopy

The physics and chemistry of a charged interface are governed by the structure of the electrical double layer (EDL). Determination of the interfacial water thickness (diw) of the charged interface is crucial to quantitatively describe the EDL structure, but it can be utilized with very scarce experimental methods. Here, we propose and verify that the vibrational relaxation time (T1) of the OH stretching mode at 3200 cm-1, obtained by time-resolved sum frequency generation vibrational spectroscopy with ssp polarizations, provides an effective tool to determine diw. By investigating the T1 values at the SiO2/NaCl solution interface, we established a time-space (T1-diw) relationship. We find that water has a T1 lifetime of ≥0.5 ps for diw ≤ 3 Å, while it displays bulk-like dynamics with T1 ≤ 0.2 ps for diw ≥ 9 Å. T1 decreases as diw increases from ∼3 Å to 9 Å. The hydration water at the DPPG lipid bilayer and LK15β protein interfaces has a thickness of ≥9 Å and shows a bulk-like feature. The time-space relationship will provide a novel tool to pattern the interfacial topography and heterogeneity in Ångstrom-depth resolution by imaging the T1 values.

Scrutinizing the protein hydration shell from molecular dynamics simulations against consensus small-angle scattering data

Biological macromolecules in solution are surrounded by a hydration shell, whose structure differs from the structure of bulk solvent. While the importance of the hydration shell for numerous biological functions is widely acknowledged, it remains unknown how the hydration shell is regulated by macromolecular shape and surface composition, mainly because a quantitative probe of the hydration shell structure has been missing. We show that small-angle scattering in solution using X-rays (SAXS) or neutrons (SANS) provide a protein-specific probe of the protein hydration shell that enables quantitative comparison with molecular simulations. Using explicit-solvent SAXS/SANS predictions, we derived the effect of the hydration shell on the radii of gyration Rg of five proteins using 18 combinations of protein force field and water model. By comparing computed Rg values from SAXS relative to SANS in D2O with consensus SAXS/SANS data from a recent worldwide community effort, we found that several but not all force fields yield a hydration shell contrast in remarkable agreement with experiments. The hydration shell contrast captured by Rg values depends strongly on protein charge and geometric shape, thus providing a protein-specific footprint of protein-water interactions and a novel observable for scrutinizing atomistic hydration shell models against experimental data.