Urea-Induced Drying of Hydrophobic Nanotubes: Comparison of Different Urea Models

Using Metadynamics To Explore the Free Energy of Dewetting in Biologically Relevant Nanopores

Water confined within hydrophobic spaces can undergo cooperative dewetting transitions due to slight changes in water density and pressure that push water toward the vapor phase. Many transmembrane protein ion channels contain nanoscale hydrophobic pores that could undergo dewetting transitions, sometimes blocking the flow of ions without physical blockages. Standard molecular dynamics simulations have been extensively applied to study the behavior of water in nanoscale pores, but the large free energy barriers of dewetting often prevent direct sampling of both wet and dry states and quantitative studies of the hydration thermodynamics of biologically relevant pores. Here, we describe a metadynamics protocol that uses the number of waters within the pore as the collective variable to drive many reversible transitions between relevant hydration states and calculate well-converged free energy profiles of pore hydration. By creating model nanopore systems and changing their radius and morphology and including various cosolvents, we quantify how these pore properties and cosolvents affect the dewetting transition. The results reveal that the dewetting free energy of nanoscale pores is determined by two key thermodynamic parameters, namely, the effective surface tension coefficients of water-air and water-pore interfaces. Importantly, while the effect of salt can be fully captured in the water activity dependence, amphipathic cosolvents such as alcohols modify both dry and wet states of the pore and dramatically shift the wet-dry equilibrium. The metadynamics approach could be applied to studies of dewetting transitions within nanoscale pores of proteins and provide new insights into why different pore properties evolved in biological systems.

Electric Field Induced Dewetting of Hydrophobic Nanocavities at Ambient Temperature

The understanding of water dewetting in nanoporous materials is of great importance in various fields of science and technology. Herein, we report molecular dynamics simulation results of dewetting of water droplet in hydrophobic nanocavities between graphene walls under the influence of electric field. At ambient temperature, the rate of dewetting induced by electric field is significantly large. Whereas, it is a very low rate of dewetting induced by high temperature (423 K) due to the strong interaction of the hydrogen-bonding networks of water droplets in nanocavities. In addition, the electric filed induced formation of a water column has been found in a vacuum chamber. When the electric field is turned off, the water column will transform into a water droplet. Importantly, the results demonstrate that the rate of electric field-induced dewetting increases with growth of the electric field. Overall, our results suggest that electric field may have a great potential application for nanomaterial dewetting.

An investigation into non-covalent functionalization of a single-walled carbon nanotube and a graphene sheet with protein G:A combined experimental and molecular dynamics study

Investigation of non-covalent interaction of hydrophobic surfaces with the protein G (PrG) is necessary due to their frequent utilization in immunosensors and ELISA. It has been confirmed that surfaces, including carbonous-nanostructures (CNS) could orient proteins for a better activation. Herein, PrG interaction with single-walled carbon nanotube (SWCNT) and graphene (Gra) nanostructures was studied by employing experimental and MD simulation techniques. It is confirmed that the PrG could adequately interact with both SWCNT and Gra and therefore fine dispersion for them was achieved in the media. Results indicated that even though SWCNT was loaded with more content of PrG in comparison with the Gra, the adsorption of the PrG on Gra did not induce significant changes in the IgG tendency. Several orientations of the PrG were adopted in the presence of SWCNT or Gra; however, SWCNT could block the PrG-FcR. Moreover, it was confirmed that SWCNT reduced the α-helical structure content in the PrG. Reduction of α-helical structure of the PrG and improper orientation of the PrG-SWCNT could remarkably decrease the PrG tendency to the Fc of the IgG. Importantly, the Gra could appropriately orient the PrG by both exposing the PrG-FcR and also by blocking the fragment of the PrG that had tendency to interact with Fab in IgG.

Exploring the binding mechanism between human profilin (PFN1) and polyproline-10 through binding mode screening

The large magnitude of protein-protein interaction (PPI) pairs within the human interactome necessitates the development of predictive models and screening tools to better understand this fundamental molecular communication. However, despite enormous efforts from various groups to develop predictive techniques in the last decade, PPI complex structures are in general still very challenging to predict due to the large number of degrees of freedom. In this study, we use the binding complex of human profilin (PFN1) and polyproline-10 (P10) as a model system to examine various approaches, with the aim of going beyond normal protein docking for PPI prediction and evaluation. The potential of mean force (PMF) was first obtained from the time-consuming umbrella sampling, which confirmed that the most stable binding structure identified by the maximal PMF difference is indeed the crystallographic binding structure. Moreover, crucial residues previously identified in experimental studies, W3, H133, and S137 of PFN1, were found to form favorable hydrogen bonds with P10, suggesting a zipping process during the binding between PFN1 and P10. We then explored both regular molecular dynamics (MD) and steered molecular dynamics (SMD) simulations, seeking for better criteria of ranking the PPI prediction. Despite valuable information obtained from conventional MD simulations, neither the commonly used interaction energy between the two binding parties nor the long-term root mean square displacement correlates well with the PMF results. On the other hand, with a sizable collection of trajectories, we demonstrated that the average and minimal rupture works calculated from SMD simulations correlate fairly well with the PMFs (R ² = 0.67), making this a promising PPI screening method.

Contrasting Effects of Guanidinium Chloride and Urea on the Activity and Unfolding of Lysozyme

Cosolvents play an important role in regulating the stability and function of proteins present in the cell. We studied the role of cosolvents, urea and guanidinium chloride (GdmCl), which act as protein denaturants, in the catalytic activity and structural stability of the protein lysozyme using activity measurements, spectroscopy, and molecular dynamics simulations. We find that the activity of lysozyme increases on the addition of urea, whereas it decreases sharply on the addition of GdmCl. At low GdmCl concentrations ([GdmCl] < 4 M), the activity of lysozyme decreases, even though there is no significant perturbation in the structure of the lysozyme folded state. We find that this is due to the strong interaction of the Gdm⁺ ion with the residues Asp52 and Glu35, which are present in the lysozyme catalytic site. In contrast, urea interacts with Trp63 present in the loop region present near the active site of lysozyme, inducing minor conformational changes in lysozyme, which can increase the activity of lysozyme. At higher denaturant concentrations, experiments show that GdmCl completely denatures the protein, whereas the folded state is stable in the presence of urea. We further show that GdmCl denatures lysozyme with the disulfide bonds intact in the protein, whereas urea denatures the protein only when the disulfide bonds are broken using reducing agents.

Mechanism of Urea Crystal Dissolution in Water from Molecular Dynamics Simulation

Molecular dynamics simulations are used to determine the mechanism of urea crystal dissolution in water under sink conditions. Crystals of cubic and tablet shapes are considered, and results are reported for four commonly used water models. The dissolution rates for different water models can differ considerably, but the overall dissolution mechanism remains the same. Urea dissolution occurs in three stages: a relatively fast initial stage, a slower intermediate stage, and a final stage. We show that the long intermediate stage is well described by classical rate laws, which assume that the dissolution rate is proportional to the active surface area. By carrying out simulations at different temperatures, we show that urea dissolution is an activated process, with an activation energy of ∼32 kJ mol^-1. Our simulations give no indication of a significant diffusion layer, and we conclude that the detachment of molecules from the crystal is the rate-determining step for dissolution. The results we report for urea are consistent with earlier observations for the dissolution of NaCl crystals. This suggests that the three-stage mechanism and classical rate laws might apply to the dissolution of other ionic and molecular crystals.

Overcoming the crystallization and designability issues in the ultrastable zirconium phosphonate framework system

Metal-organic frameworks (MOFs) based on zirconium phosphonates exhibit superior chemical stability suitable for applications under harsh conditions. These compounds mostly exist as poorly crystallized precipitates, and precise structural information has therefore remained elusive. Furthermore, a zero-dimensional zirconium phosphonate cluster acting as secondary building unit has been lacking, leading to poor designability in this system. Herein, we overcome these challenges and obtain single crystals of three zirconium phosphonates that are suitable for structural analysis. These compounds are built by previously unknown isolated zirconium phosphonate clusters and exhibit combined high porosity and ultrastability even in fuming acids. SZ-2 possesses the largest void volume recorded in zirconium phosphonates and SZ-3 represents the most porous crystalline zirconium phosphonate and the only porous MOF material reported to survive in aqua regia. SZ-2 and SZ-3 can effectively remove uranyl ions from aqueous solutions over a wide pH range, and we have elucidated the removal mechanism.

Zirconium phosphonate based metal-organic frameworks often exhibit superior chemical stabilities, but typically exist as poorly crystalline or amorphous materials. Here the authors exploit an ionothermal method to obtain highly porous and remarkably stable single crystalline zirconium phosphonate frameworks that can efficiently remove uranyl ions from aqueous solutions.

Graphene-Induced Pore Formation on Cell Membranes

Examining interactions between nanomaterials and cell membranes can expose underlying mechanisms of nanomaterial cytotoxicity and guide the design of safer nanomedical technologies. Recently, graphene has been shown to exhibit potential toxicity to cells; however, the molecular processes driving its lethal properties have yet to be fully characterized. We here demonstrate that graphene nanosheets (both pristine and oxidized) can produce holes (pores) in the membranes of A549 and Raw264.7 cells, substantially reducing cell viability. Electron micrographs offer clear evidence of pores created on cell membranes. Our molecular dynamics simulations reveal that multiple graphene nanosheets can cooperate to extract large numbers of phospholipids from the membrane bilayer. Strong dispersion interactions between graphene and lipid-tail carbons result in greatly depleted lipid density within confined regions of the membrane, ultimately leading to the formation of water-permeable pores. This cooperative lipid extraction mechanism for membrane perforation represents another distinct process that contributes to the molecular basis of graphene cytotoxicity.

Robust Denaturation of Villin Headpiece by MoS2 Nanosheet: Potential Molecular Origin of the Nanotoxicity

MoS₂ nanosheet, a new two-dimensional transition metal dichalcogenides nanomaterial, has attracted significant attentions lately due to many potential promising biomedical applications. Meanwhile, there is also a growing concern on its biocompatibility, with little known on its interactions with various biomolecules such as proteins. In this study, we use all-atom molecular dynamics simulations to investigate the interaction of a MoS₂ nanosheet with Villin Headpiece (HP35), a model protein widely used in protein folding studies. We find that MoS₂ exhibits robust denaturing capability to HP35, with its secondary structures severely destroyed within hundreds of nanosecond simulations. Both aromatic and basic residues are critical for the protein anchoring onto MoS₂ surface, which then triggers the successive protein unfolding process. The main driving force behind the adsorption process is the dispersion interaction between protein and MoS₂ monolayer. Moreover, water molecules at the interface between some key hydrophobic residues (e.g. Trp-64) and MoS₂ surface also help to accelerate the process driven by nanoscale drying, which provides a strong hydrophobic force. These findings might have shed new light on the potential nanotoxicity of MoS₂ to proteins with atomic details, which should be helpful in guiding future biomedical applications of MoS₂ with its nanotoxicity mitigated.

DNA translocation through single-layer boron nitride nanopores

Ultra-thin nanopores have become promising biological sensors because of their outstanding signal-to-noise ratio and spatial resolution. Here, we show that boron nitride (BN), which is a new two-dimensional (2D) material similar to graphene, could be utilized for making a nanopore with an atomic thickness. Using an all-atom molecular dynamics simulation, we investigated the dynamics of DNA translocation through the BN nanopore. The results of our simulations demonstrated that it is possible to detect different double-stranded DNA (dsDNA) sequences from the recording of ionic currents through the pore during the DNA translocation. Surprisingly, opposite to results for a graphene nanopore, we found the calculated blockage current for poly(A-T)40 in a BN nanopore to be less than that for poly(G-C)40. Also in contrast with the case of graphene nanopores, dsDNA models moved smoothly and in an unimpeded manner through the BN nanopores in the simulations, suggesting a potential advantage for using BN nanopores to design stall-free sequencing devices. BN nanopores, which display several properties (such as being hydrophilic and non-metallic) that are superior to those of graphene, are thus expected to find applications in the next generation of high-speed and low-cost biological sensors.

Protein corona mitigates the cytotoxicity of graphene oxide by reducing its physical interaction with cell membrane

Many recent studies have shown that the way nanoparticles interact with cells and biological molecules can vary greatly in the serum-containing or serum-free culture medium. However, the underlying molecular mechanisms of how the so-called "protein corona" formed in serum medium affects nanoparticles' biological responses are still largely unresolved. Thus, it is critical to understand how absorbed proteins on the surfaces of nanoparticles alter their biological effects. In this work, we have demonstrated with both experimental and theoretical approaches that protein BSA coating can mitigate the cytotoxicity of graphene oxide (GO) by reducing its cell membrane penetration. Our cell viability and cellular uptake experiments showed that protein corona decreased cellular uptake of GO, thus significantly mitigating the potential cytotoxicity of GO. The electron microscopy images also confirmed that protein corona reduced the cellular morphological damage by limiting GO penetration into the cell membrane. Further molecular dynamics (MD) simulations validated the experimental results and revealed that the adsorbed BSA in effect weakened the interaction between the phospholipids and graphene surface due to a reduction of the available surface area plus an unfavorable steric effect, thus significantly reducing the graphene penetration and lipid bilayer damaging. These findings provide new insights into the underlying molecular mechanism of this important graphene protein corona interaction with cell membranes, and should have implications in future development of graphene-based biomedical applications.

The role of basic residues in the adsorption of blood proteins onto the graphene surface

With its many unique properties, graphene has shown great potential in various biomedical applications, while its biocompatibility has also attracted growing concerns. Previous studies have shown that the formation of protein-graphene corona could effectively reduce its cytotoxicity; however, the underlying molecular mechanism remains not well-understood. Herein, we use extensive molecular dynamics simulations to demonstrate that blood proteins such as bovine fibrinogen (BFG) can absorb onto the graphene surface quickly and tightly to form a corona complex. Aromatic residues contributed significantly during this adsorption process due to the strong π−π stacking interactions between their aromatic rings and the graphene sp²-carbons. Somewhat surprisingly, basic residues like arginine, also played an equally or even stronger role during this process. The strong dispersion interactions between the sidechains of these solvent-exposed basic residues and the graphene surface provide the driving force for a tight binding of these basic residues. To the best of our knowledge, this is the first study with blood proteins to show that, in addition to the aromatic residues, the basic residues also play an important role in the formation of protein-graphene corona complexes.

Rotation motion of designed nano-turbine

Construction of nano-devices that can generate controllable unidirectional rotation is an important part of nanotechnology. Here, we design a nano-turbine composed of carbon nanotube and graphene nanoblades, which can be driven by fluid flow. Rotation motion of nano-turbine is quantitatively studied by molecular dynamics simulations on this model system. A robust linear relationship is achieved with this nano-turbine between its rotation rate and the fluid flow velocity spanning two orders of magnitude, and this linear relationship remains intact at various temperatures. More interestingly, a striking difference from its macroscopic counterpart is identified: the rotation rate is much smaller (by a factor of ~15) than that of the macroscopic turbine with the same driving flow. This discrepancy is shown to be related to the disruption of water flow at nanoscale, together with the water slippage at graphene surface and the so-called “dragging effect”. Moreover, counterintuitively, the ratio of “effective” driving flow velocity increases as the flow velocity increases, suggesting that the linear dependence on the flow velocity can be more complicated in nature. These findings may serve as a foundation for the further development of rotary nano-devices and should also be helpful for a better understanding of the biological molecular motors.

Effect of urea concentration on aggregation of amyloidogenic hexapeptides (NFGAIL)

We have performed large-scale all-atom molecular dynamics (MD) simulations to study the aggregation behavior of four NFGAIL hexapeptides in the aqueous urea solution, with a urea concentration ranging from 0 to 5 M. We find that urea in general suppresses the peptide aggregation, but suppression slows down in the intermediation concentration regime around 3 M. Two competing mechanisms of urea are determined: urea molecules accumulated near the first solvation shell (FSS) tend to unfold the hexapeptide, which favors aggregation; on the other hand, the tight hydrogen bonds formed between urea and peptide mainchains hinder the association of peptides which disfavors the formation of the β-sheet. Furthermore, the different nonlinear urea concentration dependences of the urea-peptide and peptide-peptide hydrogen bonds lead to a nonmonotonic behavior, with a weak enhancement in the peptide aggregation around 3 M.

Structure and dynamics of urea/water mixtures investigated by vibrational spectroscopy and molecular dynamics simulation

Urea/water is an archetypical "biological" mixture and is especially well-known for its relevance to protein thermodynamics as urea acts as a protein denaturant at high concentration. This behavior has given rise to an extended debate concerning urea's influence on water structure. On the basis of a variety of methods and of definitions of the water structure, urea has been variously described as a structure-breaker, a structure-maker, or as remarkably neutral toward water. Because of its sensitivity to microscopic structure and dynamics, vibrational spectroscopy can help resolve these debates. We report experimental and theoretical spectroscopic results for the OD stretch of HOD/H2O/urea mixtures (linear IR, 2DIR, and pump-probe anisotropy decay) and for the CO stretch of urea-D4/D2O mixtures (linear IR only). Theoretical results are obtained using existing approaches for water and a modification of a frequency map developed for acetamide. All absorption spectra are remarkably insensitive to urea concentration, consistent with the idea that urea only very weakly perturbs the water structure. Both this work and experiments by Rezus and Bakker, however, show that water's rotational dynamics are slowed down by urea. Analysis of the simulations casts doubt on the suggestion that urea immobilizes particular doubly hydrogen bonded water molecules.

Toward an atomistic description of the urea-denatured state of proteins

We present here the characterization of the structural, dynamics, and energetics of properties of the urea-denatured state of ubiquitin, a small prototypical soluble protein. By combining state-of-the-art molecular dynamics simulations with NMR and small-angle X-ray scattering data, we were able to: (i) define the unfolded state ensemble, (ii) understand the energetics stabilizing unfolded structures in urea, (iii) describe the dedifferential nature of the interactions of the fully unfolded proteins with urea and water, and (iv) characterize the early stages of protein refolding when chemically denatured proteins are transferred to native conditions. The results presented herein are unique in providing a complete picture of the chemically unfolded state of proteins and contribute to deciphering the mechanisms that stabilize the native state of proteins, as well as those that maintain them unfolded in the presence of urea.

Collapse of unfolded proteins in a mixture of denaturants

Both urea and guanidinium chloride (GdmCl) are frequently used as protein denaturants. Given that proteins generally adopt extended or unfolded conformations in either aqueous urea or GdmCl, one might expect that the unfolded protein chains will remain or become further extended due to the addition of another denaturant. However, a collapse of denatured proteins is revealed using atomistic molecular dynamics simulations when a mixture of denaturants is used. Both hen egg-white lysozyme and protein L are found to undergo collapse in the denaturant mixture. The collapse of the protein conformational ensembles is accompanied by a decreased solubility and increased non-native self-interactions of hydrophobic residues in the urea/GdmCl mixture. The increase of non-native interactions rather than the native contacts indicates that the proteins experience a simple collapse transition from the fully denatured states. During the protein collapse, the relatively stronger denaturant GdmCl displays a higher tendency to be absorbed onto the protein surface due to their stronger electrostatic interactions with proteins. At the same time, urea molecules also accumulate near the protein surface, resulting in an enhanced "local crowding" for the protein near its first solvation shell. This rearrangement of denaturants near the protein surface and crowded local environment induce the protein collapse, mainly by burying their hydrophobic residues. These findings from molecular simulations are then further explained by a simple analytical model based on statistical mechanics.

Coherent microscopic picture for urea-induced denaturation of proteins

In a previous study, we explored the mechanism of urea-induced denaturation of proteins by performing molecular dynamics (MD) simulations of hen lysozyme in 8 M urea and supported the "direct interaction mechanism" whereby urea denatures protein via dispersion interaction (Hua, L.; Zhou, R. H.; Thirumalai, D.; Berne, B. J. Proc. Natl. Acad. Sci. U.S.A. 2008, 105, 16928). Here we perform large scale MD simulations of five representative protein/peptide systems in aqueous urea to investigate if the above mechanism is common to other proteins. In all cases, accumulations of urea around proteins/peptide are observed, suggesting that urea denatures proteins by directly attacking protein backbones and side chains rather than indirectly disrupting water structure as a "water breaker". Consistent with our previous case study of lysozyme, the current energetic analyses with five protein/peptide systems reveal that urea's preferential binding to proteins mainly comes from urea's stronger dispersion interactions with proteins than with bulk solution, whereas the electrostatic (hydrogen-bonded) interactions only play a relatively minor (even negative) role during this denaturation process. Furthermore, the simulations of the peptide system at different urea concentrations (8 and 4.5 M), and with different force fields (CHARMM and OPLSAA) suggest that the above mechanism is robust, independent of the urea concentration and force field used. Last, we emphasize the importance of periodic boundary conditions in pairwise energetic analyses. This article provides a comprehensive study on the physical mechanism of urea-induced protein denaturation and suggests that the "dispersion-interaction-driven" mechanism should be general.

Water-induced ethanol dewetting transition

The dewetting transitions of two hydrophobic plates immersed in pure water, aqueous ethanol solutions with concentrations from 25% to 90%, and pure ethanol were investigated by molecular dynamics simulations, where the dewetting transition was analogous to a first-order phase transition from liquid to vapor. It was found that the dewetting transitions occurred except that in the pure ethanol system. Although the ethanol molecules prefer to locate in the vicinity of the two plates, the inter-plate region is unfavorable for water molecules, due to losing more than one hydrogen bond. Moreover, each inter-plate water molecule forms hydrogen bonds on average with about two ethanol molecules. These intermolecular hydrogen bonds cause water and ethanol to cooperatively fill or exit the inter-plate region. Thus, water molecules play a more important role in the inter-plate filling/empty process, and induce the ethanol dewetting transition. Our results provide insight into the effect of water on the ethanol dewetting phenomena.

Effect of cosolvents on nano-confined water: a molecular dynamics study

We present results from atomistic molecular dynamics simulations to characterize the effects of cosolvents, such as urea and guanidinium (Gdm) salts, on the water confined in hydrophobic carbon nanotubes. We observed complete drying of the nanotube interiors of diameter ranging from 8 to 17 Å in urea. In contrast, the water population within nanotube cores smaller than 12 Å remains unaffected in GdmCl solution, whereas larger nanotube interiors become partially dehydrated with prevailing presence of stable Gdm(+)-Gdm(+) dimers. The molecular arrangement and the lifetime inside the nanotube were found to be characteristics of a particular cosolvent. In both urea and GdmCl solutions, preferential cosolvent intrusion resulting in nanotube dehydration is driven by the stronger dispersion interaction of cosolvent than water with the nanotube. The partial drying of the hydrophobic core is attributed to guanidinium's better hydration and weaker self-association propensity compared to urea, as well as to its moderate ion-pairing with strongly hydrated chloride ions. The Gdm(+) induced dehydration varies with the charge density of counter-ions, as the presence of high charge-density sulfate ions impedes penetration of guanidinium, and consequent dehydration of the nanotube. These findings provide important insights into the effect of cosolvents on the nano-confined water in a hydrophobic environment.

Molecular wire of urea in carbon nanotube: a molecular dynamics study

We perform molecular dynamics simulations of narrow single-walled carbon nanotubes (SWNTs) in aqueous urea to investigate the structure and dynamical behavior of urea molecules inside the SWNT. Even at low urea concentrations (e.g., 0.5 M), we have observed spontaneous and continuous filling of SWNT with a one-dimensional urea wire (leaving very few water molecules inside the SWNT). The urea wire is structurally ordered, both translationally and orientationally, with a contiguous hydrogen-bonded network and concerted urea's dipole orientations. Interestingly, despite the symmetric nature of the whole system, the potential energy profile of urea along the SWNT is asymmetric, arising from the ordering of asymmetric urea partial charge distribution (or dipole moment) in confined environment. Furthermore, we study the kinetics of confined urea and find that the permeation of urea molecules through the SWNT decreases significantly (by a factor of ∼20) compared to that of water molecules, due to the stronger dispersion interaction of urea with SWNT than water, and a maximum in urea permeation happens around a concentration of 5 M. These findings might shed some light on the better understanding of unique properties of molecular wires (particularly the wires formed by polar organic small molecules) confined within both artificial and biological nanochannels, and are expected to have practical applications such as the electronic devices for signal transduction and multiplication at the nanoscale.