New Insights into the Dynamics of Zwitterionic Micelles and Their Hydration Waters by Gigahertz-to-Terahertz Dielectric Spectroscopy

Specific Binding of Primary Ammoniums in Aqueous Media by Homooxacalixarenes Incorporated into Micelles

The development of artificial receptors for efficient recognition of analytes in water is a challenging task. Homooxacalix[3]arene-based receptor 1, which is selective toward primary ammoniums in organic solvents, was transferred into water following two different strategies: direct solubilization and micellar incorporation. Extensive ¹H NMR studies showed that recognition of ammoniums is only observed in the case of micellar incorporation, highlighting the beneficial effect of the microenvironment of the micellar core. The selectivity of the system for primary ammoniums over secondary and tertiary ones was also maintained. The hydrophobic effect plays an important role in the recognition properties, which are counterion-dependent due to the energy penalty for the dissociation of certain ammonium salts in the apolar micellar core. This study shows that the straightforward self-assembly process used for the encapsulation of artificial receptors in micelles is an efficient strategy for developing water-soluble nanosized supramolecular recognition systems.

Modeling Polyzwitterion-Based Drug Delivery Platforms: A Perspective of the Current State-of-the-Art and Beyond

Drug delivery platforms are anticipated to have biocompatible and bioinert surfaces. PEGylation of drug carriers is the most approved method since it improves water solubility and colloid stability and decreases the drug vehicles’ interactions with blood components. Although this approach extends their biocompatibility, biorecognition mechanisms prevent them from biodistribution and thus efficient drug transfer. Recent studies have shown (poly)zwitterions to be alternatives for PEG with superior biocompatibility. (Poly)zwitterions are super hydrophilic, mainly stimuli-responsive, easy to functionalize and they display an extremely low protein adsorption and long biodistribution time. These unique characteristics make them already promising candidates as drug delivery carriers. Furthermore, since they have highly dense charged groups with opposite signs, (poly)zwitterions are intensely hydrated under physiological conditions. This exceptional hydration potential makes them ideal for the design of therapeutic vehicles with antifouling capability, i.e., preventing undesired sorption of biologics from the human body in the drug delivery vehicle. Therefore, (poly)zwitterionic materials have been broadly applied in stimuli-responsive “intelligent” drug delivery systems as well as tumor-targeting carriers because of their excellent biocompatibility, low cytotoxicity, insignificant immunogenicity, high stability, and long circulation time. To tailor (poly)zwitterionic drug vehicles, an interpretation of the structural and stimuli-responsive behavior of this type of polymer is essential. To this end, a direct study of molecular-level interactions, orientations, configurations, and physicochemical properties of (poly)zwitterions is required, which can be achieved via molecular modeling, which has become an influential tool for discovering new materials and understanding diverse material phenomena. As the essential bridge between science and engineering, molecular simulations enable the fundamental understanding of the encapsulation and release behavior of intelligent drug-loaded (poly)zwitterion nanoparticles and can help us to systematically design their next generations. When combined with experiments, modeling can make quantitative predictions. This perspective article aims to illustrate key recent developments in (poly)zwitterion-based drug delivery systems. We summarize how to use predictive multiscale molecular modeling techniques to successfully boost the development of intelligent multifunctional (poly)zwitterions-based systems.

Probing terahertz dynamics of multidomain protein in cell-like confinement

The development of meaningful descriptions of multidomain proteins exhibiting complex inter-domain dynamics modes is a key challenge for understanding their roles in molecular recognition and signalling processes. Here we developed a generally applicable approach for probing the low frequency collective hydration dynamics of multidomain proteins that uses terahertz spectroscopy of a protein molecule confined in a phospholipid reverse micelles environment (named Droplet THz). With the combination of normal mode analysis, we demonstrated the binding of calcium ions modulates the local inter-domain motion of the human coagulant factor VIII protein in a concentration-dependent manner. These findings highlight the Droplet THz as a valuable tool for dissecting the ultrafast dynamics of domain motion in the multidomain proteins and suggest a modulating mechanism of calcium ions on the structural flexibility and function of human coagulant proteins.

Roles of interfacial water states on advanced biomedical material design

When biomedical materials come into contact with body fluids, the first reaction that occurs on the material surface is hydration; proteins are then adsorbed and denatured on the hydrated material surface. The amount and degree of denaturation of adsorbed proteins affect subsequent cell behavior, including cell adhesion, migration, proliferation, and differentiation. Biomolecules are important for understanding the interactions and biological reactions of biomedical materials to elucidate the role of hydration in biomedical materials and their interaction partners. Analysis of the water states of hydrated materials is complicated and remains controversial; however, knowledge about interfacial water is useful for the design and development of advanced biomaterials. Herein, we summarize recent findings on the hydration of synthetic polymers, supramolecular materials, inorganic materials, proteins, and lipid membranes. Furthermore, we present recent advances in our understanding of the classification of interfacial water and advanced polymer biomaterials, based on the intermediate water concept.

Probing Adaptation of Hydration and Protein Dynamics to Temperature

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

The origin and impact of bound water around intrinsically disordered proteins

Proteins and water couple dynamically over a wide range of time scales. Motivated by their central role in protein function, protein-water dynamics and thermodynamics have been extensively studied for structured proteins, where correspondence to structural features has been made. However, properties controlling intrinsically disordered protein (IDP)-water dynamics are not yet known. We report results of megahertz-to-terahertz dielectric spectroscopy and molecular dynamics simulations of a group of IDPs with varying charge content along with structured proteins of similar size. Hydration water around IDPs is found to exhibit more heterogeneous rotational and translational dynamics compared with water around structured proteins of similar size, yielding on average more restricted dynamics around individual residues of IDPs, charged or neutral, compared with structured proteins. The on-average slower water dynamics is found to arise from excess tightly bound water in the first hydration layer, which is related to greater exposure to charged groups. The more tightly bound water to IDPs correlates with the smaller hydration shell found experimentally, and affects entropy associated with protein-water interactions, the contribution of which we estimate based on the dielectric measurements and simulations. Water-IDP dynamic coupling at terahertz frequencies is characterized by the dielectric measurements and simulations.

Using Molecular Simulation to Guide Protein Engineering for Biocatalysis in Organic Solvents

Biocatalysis in organic solvents (OSs) is very appealing for the industry in producing bulk and/or fine chemicals, such as pharmaceuticals, biodiesel, and fragrances. The poor performance of enzymes in OSs (e.g., reduced activity, insufficient stability, and deactivation) negates OSs' excellent solvent properties. Molecular dynamics (MD) simulations provide a complementary method to study the relationship between enzymes dynamics and the stability in OSs. Here we describe computational procedure for MD simulation of enzymes in OSs with an example of Bacillus subtilis lipase A (BSLA) in dimethyl sulfoxide (DMSO) cosolvent with software GROMACS. We discuss main essential practical issues considered (such as choice of force field, parameterization, simulation setup, and trajectory analysis). The core part of this protocol (enzyme-OS system setup, analysis of structural-based and solvation-based observables) is transferable to other enzymes and any OS systems. Combining with experimental studies, the obtained molecular knowledge is most likely to guide researchers to access rational protein engineering approaches to tailor OS resistant enzymes and expand the scope of biocatalysis in OS media. Finally, we discuss potential solutions to overcome the remaining challenges of computational biocatalysis in OSs and briefly draw future directions for further improvement in this field.

Long-range DNA-water interactions

DNA functions only in aqueous environments and adopts different conformations depending on the hydration level. The dynamics of hydration water and hydrated DNA leads to rotating and oscillating dipoles that, in turn, give rise to a strong megahertz to terahertz absorption. Investigating the impact of hydration on DNA dynamics and the spectral features of water molecules influenced by DNA, however, is extremely challenging because of the strong absorption of water in the megahertz to terahertz frequency range. In response, we have employed a high-precision megahertz to terahertz dielectric spectrometer, assisted by molecular dynamics simulations, to investigate the dynamics of water molecules within the hydration shells of DNA as well as the collective vibrational motions of hydrated DNA, which are vital to DNA conformation and functionality. Our results reveal that the dynamics of water molecules in a DNA solution is heterogeneous, exhibiting a hierarchy of four distinct relaxation times ranging from ∼8 ps to 1 ns, and the hydration structure of a DNA chain can extend to as far as ∼18 Å from its surface. The low-frequency collective vibrational modes of hydrated DNA have been identified and found to be sensitive to environmental conditions including temperature and hydration level. The results reveal critical information on hydrated DNA dynamics and DNA-water interfaces, which impact the biochemical functions and reactivity of DNA.

Protein Solvent Shell Structure Provides Rapid Analysis of Hydration Dynamics

The solvation layer surrounding a protein is clearly an intrinsic part of protein structure-dynamics-function, and our understanding of how the hydration dynamics influences protein function is emerging. We have recently reported simulations indicating a correlation between regional hydration dynamics and the structure of the solvation layer around different regions of the enzyme Candida antarctica lipase B, wherein the radial distribution function (RDF) was used to calculate the pairwise entropy, providing a link between dynamics (diffusion) and thermodynamics (excess entropy) known as Rosenfeld scaling. Regions with higher RDF values/peaks in the hydration layer (the first peak, within 6 Å of the protein surface) have faster diffusion in the hydration layer. The finding thus hinted at a handle for rapid evaluation of hydration dynamics at different regions on the protein surface in molecular dynamics simulations. Such an approach may move the analysis of hydration dynamics from a specialized venture to routine analysis, enabling an informatics approach to evaluate the role of hydration dynamics in biomolecular function. This paper first confirms that the correlation between regional diffusive dynamics and hydration layer structure (via water center of mass around protein side-chain atom RDF) is observed as a general relationship across a set of proteins. Second, it seeks to devise an approach for rapid analysis of hydration dynamics, determining the minimum amount of information and computational effort required to get a reliable value of hydration dynamics from structural data in MD simulations based on the protein-water RDF. A linear regression model using the integral of the hydration layer in the water-protein RDF was found to provide statistically equivalent apparent diffusion coefficients at the 95% confidence level for a set of 92 regions within five different proteins. In summary, RDF analysis of 10 ns of data after simulation convergence is sufficient to accurately map regions of fast and slow hydration dynamics around a protein surface. Additionally, it is anticipated that a quick look at protein-water RDFs, comparing peak heights, will be useful to provide a qualitative ranking of regions of faster and slower hydration dynamics at the protein surface for rapid analysis when investigating the role of solvent dynamics in protein function.

Functional Hydration Behavior: Interrelation between Hydration and Molecular Properties at Lipid Membrane Interfaces

Water is an abundant commodity and has various important functions. It stabilizes the structure of biological macromolecules, controls biochemical activities, and regulates interfacial/intermolecular interactions. Common aspects of interfacial water can be obtained by overviewing fundamental functions and properties at different temporal and spatial scales. It is important to understand the hydrogen bonding and structural properties of water and to evaluate the individual molecular species having different hydration properties. Water molecules form hydrogen bonds with biomolecules and contribute to the adjustment of their properties, such as surface charge, hydrophilicity, and structural flexibility. In this review, the fundamental properties of water molecules and the methods used for the analyses of water dynamics are summarized. In particular, the interrelation between the hydration properties, determined by molecules, and the properties of molecules, determined by their hydration properties, are discussed using the lipid membrane as an example. Accordingly, interesting water functions are introduced that provide beneficial information in the fields of biochemistry, medicine, and food chemistry.

Probing relaxation dynamics of a cationic lipid based non-viral carrier: a time-resolved fluorescence study

Lipid vesicles composed of cationic dioctadecyldimethylammonium bromide (DODAB) and neutral 1-monooleoyl-rac-glycerol (MO) are promising non-viral carriers of nucleic acids for delivery into cells. Among them, higher cell transfection efficiency was displayed by DODAB-rich vesicles than those enriched with MO. Structural relaxation of these mixed lipid vesicles plays a key role in their cell transfection efficiency because structural organization of the DODAB-rich vesicles are different from that of the MO-rich vesicles. Polarization-gated anisotropy in conjunction with picosecond resolved emission transients of a novel fluorophore 6-acetyl-(2-((4-hydroxycyclohexyl)(methyl)amino)naphthalene) (ACYMAN) has been employed to probe relaxation dynamics in pure DODAB vesicles, and in mixed vesicles of DODAB with varying content of MO. Both orientational relaxation of ACYMAN and relaxation dynamics of its local environment are retarded significantly in mixed lipid vesicles with increasing MO content, from a mole fraction (χ_MO) of 0.2 to that of 0.8 which is consistent with increased rigidity of the MO-rich (χ_MO > 0.5) vesicles relative to the DODAB-rich (χ_MO < 0.5) vesicles. Therefore, higher structural rigidity of the MO-rich vesicles (χ_MO > 0.5) gives rise to their lower cell transfection efficiency than the more flexible DODAB-rich (χ_MO < 0.5) vesicles as observed in previous in vivo studies (Biochim. Biophys. Acta, Biomembr., 2014, 1838, 2555–2567).

Lipid vesicles composed of cationic dioctadecyldimethylammonium bromide (DODAB) and neutral 1-monooleoyl-rac-glycerol (MO) are promising non-viral carriers of nucleic acids for delivery into cells.

Specific Ion Effects on Zwitterionic Micelles Are Independent of Interfacial Hydration Changes

Zwitterionic micelles adsorb anions and several techniques were used to determine the specificity of this interaction. Although at a lower intensity, this adsorption can be compared to those observed in cationic micelles, which showed that interfacial dehydration is a fundamental property for the geometry and size of micelles. Because there is no information on the interfacial hydration of zwitterionic micelles, we used dielectric relaxation spectroscopy (DRS) together with molecular dynamics (MD) simulations to evaluate the importance of surface dehydration promoted by the binding of anions at the micellar interface (sodium bromide, sodium methanesulfonate, sodium trifluoroacetate, and sodium triflate) in N-dodecyl- N, N-dimethyl-3-ammonio-1-propanesulfonate (DPS) micelles. Our results, showing good agreement between DRS and MD simulations, strongly suggest that specific ion effects on zwitterionic micelles are unrelated to global changes in the interfacial hydration and depend on specific interactions of the headgroups with selected anions.

How Does Solvation Layer Mobility Affect Protein Structural Dynamics?

Solvation is critical for protein structural dynamics. Spectroscopic studies have indicated relationships between protein and solvent dynamics, and rates of gas binding to heme proteins in aqueous solution were previously observed to depend inversely on solution viscosity. In this work, the solvent-compatible enzyme Candida antarctica lipase B, which functions in aqueous and organic solvents, was modeled using molecular dynamics simulations. Data was obtained for the enzyme in acetonitrile, cyclohexane, n-butanol, and tert-butanol, in addition to water. Protein dynamics and solvation shell dynamics are characterized regionally: for each α-helix, β-sheet, and loop or connector region. Correlations are seen between solvent mobility and protein flexibility. So, does local viscosity explain the relationship between protein structural dynamics and solvation layer dynamics? Halle and Davidovic presented a cogent analysis of data describing the global hydrodynamics of a protein (tumbling in solution) that fits a model in which the protein's interfacial viscosity is higher than that of bulk water's, due to retarded water dynamics in the hydration layer (measured in NMR τ2 reorientation times). Numerous experiments have shown coupling between protein and solvation layer dynamics in site-specific measurements. Our data provides spatially-resolved characterization of solvent shell dynamics, showing correlations between regional solvation layer dynamics and protein dynamics in both aqueous and organic solvents. Correlations between protein flexibility and inverse solvent viscosity (1/η) are considered across several protein regions and for a rather disparate collection of solvents. It is seen that the correlation is consistently higher when local solvent shell dynamics are considered, rather than bulk viscosity. Protein flexibility is seen to correlate best with either the local interfacial viscosity or the ratio of the mobility of an organic solvent in a regional solvation layer relative to hydration dynamics around the same region. Results provide insight into the function of aqueous proteins, while also suggesting a framework for interpreting and predicting enzyme structural dynamics in non-aqueous solvents, based on the mobility of solvents within the solvation layer. We suggest that Kramers' theory may be used in future work to model protein conformational transitions in different solvents by incorporating local viscosity effects.

High-Precision Megahertz-to-Terahertz Dielectric Spectroscopy of Protein Collective Motions and Hydration Dynamics

The low-frequency collective vibrational modes in proteins as well as the protein-water interface have been suggested as dominant factors controlling the efficiency of biochemical reactions and biological energy transport. It is thus crucial to uncover the mystery of the hydration structure and dynamics as well as their coupling to collective motions of proteins in aqueous solutions. Here, we report dielectric properties of aqueous bovine serum albumin protein solutions as a model system using an extremely sensitive dielectric spectrometer with frequencies spanning from megahertz to terahertz. The dielectric relaxation spectra reveal several polarization mechanisms at the molecular level with different time constants and dielectric strengths, reflecting the complexity of protein-water interactions. Combining the effective-medium approximation and molecular dynamics simulations, we have determined collective vibrational modes at terahertz frequencies and the number of water molecules in the tightly bound and loosely bound hydration layers. High-precision measurements of the number of hydration water molecules indicate that the dynamical influence of proteins extends beyond the first solvation layer, to around 7 Å distance from the protein surface, with the largest slowdown arising from water molecules directly hydrogen-bonded to the protein. Our results reveal critical information of protein dynamics and protein-water interfaces, which determine biochemical functions and reactivity of proteins.

How do the interfacial properties of zwitterionic sulfobetaine micelles differ from those of cationic alkyl quaternary ammonium micelles? An excited state proton transfer study

The interface of a zwitterionic sulfobetaine micelle is more packed and less hydrated compared to a cationic alkyl-ammonium micelle with an identical alkyl tail.