Sequence Modulates Polypeptoid Hydration Water Structure and Dynamics

Water-directed pinning is key to tau prion formation

Tau forms fibrillar aggregates that are pathological hallmarks of a family of neurodegenerative diseases known as tauopathies. The synthetic replication of disease-specific fibril structures is a critical gap for developing diagnostic and therapeutic tools. This study debuts a strategy of identifying a critical and minimal folding motif in fibrils characteristic of tauopathies and generating seeding-competent fibrils from the isolated tau peptides. The 19-residue jR2R3 peptide (295 to 313) which spans the R2/R3 splice junction of tau, and includes the P301L mutation, is one such peptide that forms prion-competent fibrils. This tau fragment contains the hydrophobic VQIVYK hexapeptide that is part of the core of all known pathological tau fibril structures and an intramolecular counterstrand that stabilizes the strand-loop-strand (SLS) motif observed in 4R tauopathy fibrils. This study shows that P301L exhibits a duality of effects: it lowers the barrier for the peptide to adopt aggregation-prone conformations and enhances the local structuring of water around the mutation site to facilitate site-directed pinning and dewetting around sites 300-301 to achieve in-register stacking of tau to cross β-sheets. We solved a 3 Å cryo-EM structure of jR2R3-P301L fibrils in which each protofilament layer contains two jR2R3-P301L copies, of which one adopts a SLS fold found in 4R tauopathies and the other wraps around the SLS fold to stabilize it, reminiscent of the three- and fourfold structures observed in 4R tauopathies. These jR2R3-P301L fibrils are competent to template full-length 4R tau in a prion-like manner.

Hydration-induced dynamical changes in lyophilised and weakly hydrated apoferritin: insights from molecular dynamics simulation

The dynamics and functionality of proteins are significantly influenced by their interaction with water. For lyophilised (i.e. h ≤ 0.05 where h = g of H2O per g of protein) and weakly hydrated systems (i.e. h ≤ 0.38) hydration generally enhances protein mobility above the so-called 'dynamical transition' temperature (Td > 220 K). However, water-induced mobility hindrance at low temperatures (T < 175 K) has been reported in various proteins of varying secondary structure; namely green fluorescent protein (GFP), pig liver esterase, lysozyme, ribonuclease A (RNAse A) and apoferritin. By focussing on the dynamic behaviour of the apoferritin molecule, this study proposes mechanisms driving these hydration-induced mobility changes, particularly the less understood hindrance at low temperatures. Using atomistic molecular dynamics (MD) simulations of horse spleen apoferritin in the lyophilised (h = 0.05) and weakly hydrated (h = 0.31) states, we report here the impact of water on protein dynamics as a function of temperature. Through residue-specific mean squared displacement (MSD), radial distribution function (RDF), solvent accessible surface area (SASA), local hydration degree and hydrogen bonding analyses, we demonstrate that while water proximity directly correlates with mobility enhancement at high temperatures, the hydration-induced mobility reduction observed at temperatures below 175 K is primarily propagated through the protein backbone.

Structure and Dynamics of Water in Polysaccharide (Alginate) Solutions and Gels Explained by the Core-Shell Model

In both biological and engineered systems, polysaccharides offer a means of establishing structural stiffness without altering the availability of water. Notable examples include the extracellular matrix of prokaryotes and eukaryotes, artificial skin grafts, drug delivery materials, and gels for water harvesting. Proper design and modeling of these systems require detailed understanding of the behavior of water confined in pores narrower than about 1 nm. We use molecular dynamics simulations to investigate the properties of water in solutions and gels of the polysaccharide alginate as a function of the water content and polymer cross-linking. We find that a detailed understanding of the nanoscale dynamics of water in alginate solutions and gels requires consideration of the discrete nature of water. However, we also find that the trends in tortuosity, permeability, dielectric constant, and shear viscosity can be adequately represented using the "core-shell" conceptual model that considers the confined fluid as a continuum.

Bridging the Gap in Cryopreservation Mechanism: Unraveling the Interplay between Structure, Dynamics, and Thermodynamics in Cryoprotectant Aqueous Solutions

Cryoprotectants play a crucial role in preserving biological material, ensuring their viability during storage and facilitating crucial applications such as the conservation of medical compounds, tissues, and organs for transplantation. However, the precise mechanism by which cryoprotectants modulate the thermodynamic properties of water to impede the formation and growth of ice crystals, thus preventing long-term damage, remains elusive. This is evident in the use of empirically optimized recipes for mixtures that typically contain DMSO, glycerol, and various sugar constituents. Here, we use terahertz calorimetry, Overhauser nuclear polarization, and molecular dynamics simulations to show that DMSO exhibits a robust structuring effect on water around its methyl groups, reaching a maximum at a DMSO mole fraction of XDMSO = 0.33. In contrast, glycerol exerts a smaller water-structuring effect, even at higher concentrations (Scheme 1). These results potentially suggest that the wrapped water around DMSO's methyl group, which can be evicted upon ligand binding, may render DMSO a more surface-active cryoprotectant than glycerol, while glycerol may participate more as a viscogen that acts on the entire sample. These findings shed light on the molecular intricacies of cryoprotectant solvation behavior and have potentially significant implications for optimizing cryopreservation protocols.

Inverse design of pore wall chemistry and topology through active learning of surface group interactions

Design of next-generation membranes requires a nanoscopic understanding of the effect of biologically inspired heterogeneous surface chemistries and topologies (roughness) on local water and solute behavior. In particular, the rejection of small, neutral solutes, such as boric acid, poses a heretofore unsolved challenge. In prior work, a computational inverse design technique using an evolutionary optimization successfully uncovered new surface design strategies for optimized transport of water over solutes in smooth, model pores consisting of two surface chemistries. However, extending such an approach to more complex (and realistic) scenarios involving many surface chemistries as well as surface roughness is challenging due to the expanded design space. In this work, we develop a new approach that uses active learning to optimize in a reduced feature space of surface group interactions, finding parameters that lead to their assembly into ordered, optimal patterns. This approach rapidly identifies novel surface functionalizations that maximize the difference in water and boric acid transport through the nanopore. Moreover, we find that the roughness of the nanopore wall, independent of its chemistry, can be leveraged to enhance transport selectivity: oscillations in the pore wall diameter optimally inhibit boric acid transport by creating energetic wells from which the solute must escape to transport down the pore. This proof-of-concept demonstrates the potential for active learning strategies, in concert with molecular simulations, to rapidly navigate complex design spaces of aqueous interfaces and is promising as a tool for engineering water-mediated surface interactions for a broad range of applications.

Nanoscale water-polymer interactions tune macroscopic diffusivity of water in aqueous poly(ethylene oxide) solutions

The separation and anti-fouling performance of water purification membranes is governed by both macroscopic and molecular-scale water properties near polymer surfaces. However, even for poly(ethylene oxide) (PEO) - ubiquitously used in membrane materials - there is little understanding of whether or how the molecular structure of water near PEO surfaces affects macroscopic water diffusion. Here, we probe both time-averaged bulk and local water dynamics in dilute and concentrated PEO solutions using a unique combination of experimental and simulation tools. Pulsed-Field Gradient NMR and Overhauser Dynamic Nuclear Polarization (ODNP) capture water dynamics across micrometer length scales in sub-seconds to sub-nanometers in tens of picoseconds, respectively. We find that classical models, such as the Stokes-Einstein and Mackie-Meares relations, cannot capture water diffusion across a wide range of PEO concentrations, but that free volume theory can. Our study shows that PEO concentration affects macroscopic water diffusion by enhancing the water structure and altering free volume. ODNP experiments reveal that water diffusivity near PEO is slower than in the bulk in dilute solutions, previously not recognized by macroscopic transport measurements, but the two populations converge above the polymer overlap concentration. Molecular dynamics simulations reveal that the reduction in water diffusivity occurs with enhanced tetrahedral structuring near PEO. Broadly, we find that PEO does not simply behave like a physical obstruction but directly modifies water's structural and dynamic properties. Thus, even in simple PEO solutions, molecular scale structuring and the impact of polymer interfaces is essential to capturing water diffusion, an observation with important implications for water transport through structurally complex membrane materials.

Relationships between Water's Structure and Solute Affinity at Polypeptoid Brush Surfaces

Excellent antifouling surfaces are generally thought to create a tightly bound layer of water that resists solute adsorption, and highly hydrophilic surfaces such as those with zwitterionic functionalities are of significant current interest as antifoulant strategies. However, despite significant proofs-of-concept, we still lack a fundamental understanding of how the nanoscopic structure of this hydration layer translates to reduced fouling, how surface chemistry can be tuned to achieve antifouling through hydration water, and why, in particular, zwitterionic surfaces seem so promising. Here, we use molecular dynamics simulations and free energy calculations to investigate the molecular relationships among surface chemistry, hydration water structure, and surface-solute affinity across a variety of surface-decorated chemistries. Specifically, we consider polypeptoid-decorated surfaces that display well-known experimental antifouling capabilities and that can be synthesized sequence specifically, with precise backbone positioning of, e.g., charged groups. Through simulations, we calculate the affinities of a range of small solutes to polypeptoid brush surfaces of varied side-chain chemistries. We then demonstrate that measures of the structure of surface hydration water in response to a particular surface chemistry signal solute-surface affinity; specifically, we find that zwitterionic chemistries produce solute-surface repulsion through highly coordinated hydration water while suppressing tetrahedral structuring around the solute, in contrast to uncharged surfaces that show solute-surface affinity. Based on the relationship of this structural perturbation to the affinity of small-molecule solutes, we propose a molecular mechanism by which zwitterionic surface chemistries enhance solute repulsion, with broader implications for the design of antifouling surfaces.

Tau P301L mutation promotes core 4R tauopathy fibril fold through near-surface water structuring and conformational rearrangement

Tau forms toxic fibrillar aggregates in a family of neurodegenerative diseases known as tauopathies. The faithful replication of tauopathy-specific fibril structures is a critical gap for developing diagnostic and therapeutic tools. This study debuts a strategy of identifying a critical segment of tau that forms a folding motif that is characteristic of a family of tauopathies and isolating it as a standalone peptide that form seeding-competent fibrils. The 19-residue jR2R3 peptide (295-313) spanning the R2/R3 splice junction of tau, in the presence of P301L, forms seeding-competent amyloid fibrils. This tau fragment contains the hydrophobic VQIVYK hexapeptide that is part of the core of every pathological tau fibril structure solved to-date and an intramolecular counter-strand that stabilizes the strand-loop-strand (SLS) motif observed in 4R tauopathy fibrils. This study shows that P301L exhibits a duality of effects: it lowers the barrier for the peptide to adopt aggregation-prone conformations and enhances the local structuring of water around the mutation site that facilitates site-specific dewetting and in-register stacking of tau to form cross β-sheets. We solve a 3 Å cryo-EM structure of jR2R3-P301L fibrils with a pseudo 2 1 screw symmetry in which each half of the fibril's cross-section contains two jR2R3-P301L peptides. One chain adopts a SLS fold found in 4R tauopathies that is stabilized by a second chain wrapping around the SLS fold, reminiscent of the 3-fold and 4-fold structures observed in 4R tauopathies. These jR2R3-P301L fibrils are able to template full length tau in a prion-like fashion.

SIGNIFICANCE STATEMENT

This study presents a first step towards designing a tauopathy specific aggregation pathway by engineering a minimal tau prion building block, jR2R3, that can template and propagate distinct disease folds. We present the discovery that P301L-among the widest used mutations in cell and animal models of Alzheimer's Disease-destabilizes an aggregation-prohibiting internal hairpin and enhances the local surface water structure that serves as an entropic hotspot to exert a hyper-localized effect in jR2R3. Our study suggests that P301L may be a more suitable mutation to include in modeling 4R tauopathies than for modelling Alzheimer's Disease, and that mutations are powerful tools for the purpose of designing of tau prion models as therapeutic tools.

Development of a Systematic and Extensible Force Field for Peptoids (STEPs)

Peptoids (N-substituted glycines) are a class of biomimetic polymers that have attracted significant attention due to their accessible synthesis and enzymatic and thermal stability relative to their naturally occurring counterparts (polypeptides). While these polymers provide the promise of more robust functional materials via hierarchical approaches, they present a new challenge for computational structure prediction for material design. The reliability of calculations hinges on the accuracy of interactions represented in the force field used to model peptoids. For proteins, structure prediction based on sequence and de novo design has made dramatic progress in recent years; however, these models are not readily transferable for peptoids. Current efforts to develop and implement peptoid-specific force fields are spread out, leading to replicated efforts and a fragmented collection of parameterized sidechains. Here, we developed a peptoid-specific force field containing 70 different side chains, using GAFF2 as starting point. The new model is validated based on the generation of Ramachandran-like plots from DFT optimization compared against force field reproduced potential energy and free energy surfaces as well as the reproduction of equilibrium cis/trans values for some residues experimentally known to form helical structures. Equilibrium cis/trans distributions (Kct) are estimated for all parameterized residues to identify which residues have an intrinsic propensity for cis or trans states in the monomeric state.

Local solvation structures govern the mixing thermodynamics of glycerol-water solutions

Glycerol is a major cryoprotective agent and is widely used to promote protein stabilization. By a combined experimental and theoretical study, we show that global thermodynamic mixing properties of glycerol and water are dictated by local solvation motifs. We identify three hydration water populations, i.e., bulk water, bound water (water hydrogen bonded to the hydrophilic groups of glycerol) and cavity wrap water (water hydrating the hydrophobic moieties). Here, we show that for glycerol experimental observables in the THz regime allow quantification of the abundance of bound water and its partial contribution to the mixing thermodynamics. Specifically, we uncover a 1 : 1 connection between the population of bound waters and the mixing enthalpy, which is further corroborated by the simulation results. Therefore, the changes in global thermodynamic quantity - mixing enthalpy - are rationalized at the molecular level in terms of changes in the local hydrophilic hydration population as a function of glycerol mole fraction in the full miscibility range. This offers opportunities to rationally design polyol water, as well as other aqueous mixtures to optimize technological applications by tuning mixing enthalpy and entropy based on spectroscopic screening.

Design of Soft Material Surfaces with Rationally Tuned Water Diffusivity

Water structure and dynamics can be key modulators of adsorption, separations, and reactions at soft material interfaces, but systematically tuning water environments in an aqueous, accessible, and functionalizable material platform has been elusive. This work leverages variations in excluded volume to control and measure water diffusivity as a function of position within polymeric micelles using Overhauser dynamic nuclear polarization spectroscopy. Specifically, a versatile materials platform consisting of sequence-defined polypeptoids simultaneously offers a route to controlling the functional group position and a unique opportunity to generate a water diffusivity gradient extending away from the polymer micelle core. These results demonstrate an avenue not only to rationally design the chemical and structural properties of polymer surfaces but also to design and tune the local water dynamics that, in turn, can adjust the local activity for solutes.

Relationships between Molecular Structural Order Parameters and Equilibrium Water Dynamics in Aqueous Mixtures

Water's unique thermophysical properties and how it mediates aqueous interactions between solutes have long been interpreted in terms of its collective molecular structure. The seminal work of Errington and Debenedetti [Nature 2001, 409, 318-321] revealed a striking hierarchy of relationships among the thermodynamic, dynamic, and structural properties of water, motivating many efforts to understand (1) what measures of water structure are connected to different experimentally accessible macroscopic responses and (2) how many such structural metrics are adequate to describe the collective structural behavior of water. Diffusivity constitutes a particularly interesting experimentally accessible equilibrium property to investigate such relationships because advanced NMR techniques allow the measurement of bulk and local water dynamics in nanometer proximity to molecules and interfaces, suggesting the enticing possibility of measuring local diffusivities that report on water structure. Here, we apply statistical learning methods to discover persistent structure-dynamic correlations across a variety of simulated aqueous mixtures, from alcohol-water to polypeptoid-water systems. We investigate a variety of molecular water structure metrics and find that an unsupervised statistical learning algorithm (namely, sequential feature selection) identifies only two or three independent structural metrics that are sufficient to predict water self-diffusivity accurately. Surprisingly, the translational diffusivity of water across all mixed systems studied here is strongly correlated with a measure of tetrahedral order given by water's triplet angle distribution. We also identify a separate small number of structural metrics that well predict an important thermodynamic property, the excess chemical potential of an idealized methane-sized hydrophobe in water. Ultimately, we offer a Bayesian method of inferring water structure by using only structure-dynamics linear regression models with experimental Overhauser dynamic nuclear polarization (ODNP) measurements of water self-diffusivity. This study thus quantifies the relationships among several distinct structural order parameters in water and, through statistical learning, reveals the potential to leverage molecular structure to predict fundamental thermophysical properties. In turn, these findings suggest a framework for solving the inverse problem of inferring water's molecular structure using experimental measurements such as ODNP studies that probe local water properties.

Lipidation Alters the Structure and Hydration of Myristoylated Intrinsically Disordered Proteins

Lipidated proteins are an emerging class of hybrid biomaterials that can integrate the functional capabilities of proteins into precisely engineered nano-biomaterials with potential applications in biotechnology, nanoscience, and biomedical engineering. For instance, fatty-acid-modified elastin-like polypeptides (FAMEs) combine the hierarchical assembly of lipids with the thermoresponsive character of elastin-like polypeptides (ELPs) to form nanocarriers with emergent temperature-dependent structural (shape or size) characteristics. Here, we report the biophysical underpinnings of thermoresponsive behavior of FAMEs using computational nanoscopy, spectroscopy, scattering, and microscopy. This integrated approach revealed that temperature and molecular syntax alter the structure, contact, and hydration of lipid, lipidation site, and protein, aligning with the changes in the nanomorphology of FAMEs. These findings enable a better understanding of the biophysical consequence of lipidation in biology and the rational design of the biomaterials and therapeutics that rival the exquisite hierarchy and capabilities of biological systems.