Interaction of an ionic complementary peptide with a hydrophobic graphite surface

The impact of chemical properties of the solid-liquid-adsorbate interfaces on the entropy-enthalpy compensation involved in adsorption

Despite extensive studies on the thermodynamic mechanism governing molecular adsorption at the solid-water interface, a comprehensive understanding of the crucial role of interface properties in mediating the entropy-enthalpy compensation during adsorption is lacking, particularly at a quantitative level. Herein, we employed two types of surface models (hydroxyapatite and graphene) along with a series of amino acids to successfully elucidate how distinct interfacial features dictate the delicate balance between entropy and enthalpy variations. The adsorption of all amino acids on the hydroxyapatite surface is an enthalpy-dominated process, where the water-induced enthalpic component of the free energy and the surface-adsorbate electrostatic interaction term alternatively act as the driving force for adsorption in different regions of the surface. Although favorable interactions are observed between amino acids and the graphene surface, the entropy-enthalpy compensation exhibits dependence on the molecular size of the adsorbates. For small amino acids, favorable enthalpy changes predominantly determine their adsorption behavior; however, larger amino acids tend to bind more tightly with the graphene surface, which is thermodynamically dominated by the entropy variations despite the structural characteristics of amino acids. This study reveals specific entropy-enthalpy mechanisms underlying amino acid adsorption at the solid-liquid interface, providing guidance for surface design and synthesis of new biomolecules.

Recent Advances of Studies on Cell-Penetrating Peptides Based on Molecular Dynamics Simulations

With the ability to transport cargo molecules across cell membranes with low toxicity, cell-penetrating peptides (CPPs) have become promising candidates for next generation peptide-based drug delivery vectors. Over the past three decades since the first CPP was discovered, a great deal of work has been done on the cellular uptake mechanisms and the applications for the delivery of therapeutic molecules, and significant advances have been made. But so far, we still do not have a precise and unified understanding of the structure-activity relationship of the CPPs. Molecular dynamics (MD) simulations provide a method to reveal peptide-membrane interactions at the atomistic level and have become an effective complement to experiments. In this paper, we review the progress of the MD simulations on CPP-membrane interactions, including the computational methods and technical improvements in the MD simulations, the research achievements in the CPP internalization mechanism, CPP decoration and coupling, and the peptide-induced membrane reactions during the penetration process, as well as the comparison of simulated and experimental results.

Hierarchical Nanomaterials Assembled from Peptoids and Other Sequence-Defined Synthetic Polymers

In nature, the self-assembly of sequence-specific biopolymers into hierarchical structures plays an essential role in the construction of functional biomaterials. To develop synthetic materials that can mimic and surpass the function of these natural counterparts, various sequence-defined bio- and biomimetic polymers have been developed and exploited as building blocks for hierarchical self-assembly. This review summarizes the recent advances in the molecular self-assembly of hierarchical nanomaterials based on peptoids (or poly-N-substituted glycines) and other sequence-defined synthetic polymers. Modern techniques to monitor the assembly mechanisms and characterize the physicochemical properties of these self-assembly systems are highlighted. In addition, discussions about their potential applications in biomedical sciences and renewable energy are also included. This review aims to highlight essential features of sequence-defined synthetic polymers (e.g., high stability and protein-like high-information content) and how these unique features enable the construction of robust biomimetic functional materials with high programmability and predictability, with an emphasis on peptoids and their self-assembled nanomaterials.

Hierarchical Assembly of Peptoid-Based Cylindrical Micelles Exhibiting Efficient Resonance Energy Transfer in Aqueous Solution

Herein we show that by appending bulky β-cyclodextrin (CD) groups onto sheet-forming peptoids, we obtain cylindrical micelles that further assembly into membranes and intertwined ribbons on substrates in aqueous solution, depending on the choice of solution and substrate conditions. In situ atomic force microscopy (AFM) shows that micelle assembly occurs in two steps, starting with "precursor" particles that transform into worm-like micelles, which extend and coalesce to form the higher order structures with a rate and a degree of cooperativity dependent on pH and Ca²⁺ concentration. After co-assembly with hydrophobic 4-(2-hydroxyethylamino)-7-nitro-2,1,3-benzoxadiazole (NBD) donors that occupy the hydrophobic core, followed by exposure to hydrophilic Rhodamine B as acceptors that insert into cyclodextrin, the micelles exhibit highly efficient Förster resonance energy transfer efficiency in aqueous solution, thereby mimicking natural light harvesting systems.

Potential disruption of protein-protein interactions by graphene oxide

Graphene oxide (GO) is a promising novel nanomaterial with a wide range of potential biomedical applications due to its many intriguing properties. However, very little research has been conducted to study its possible adverse effects on protein-protein interactions (and thus subsequent toxicity to human). Here, the potential cytotoxicity of GO is investigated at molecular level using large-scale, all-atom molecular dynamics simulations to explore the interaction mechanism between a protein dimer and a GO nanosheet oxidized at different levels. Our theoretical results reveal that GO nanosheet could intercalate between the two monomers of HIV-1 integrase dimer, disrupting the protein-protein interactions and eventually lead to dimer disassociation as graphene does [B. Luan et al., ACS Nano 9(1), 663 (2015)], albeit its insertion process is slower when compared with graphene due to the additional steric and attractive interactions. This study helps to better understand the toxicity of GO to cell functions which could shed light on how to improve its biocompatibility and biosafety for its wide potential biomedical applications.

Rubredoxin refolding on nanostructured hydrophobic surfaces: evidence for a new type of biomimetic chaperones

Rubredoxins (Rds) are small proteins containing a tetrahedral Fe(SCys)4 site. Folded forms of metal free Rds (apoRds) show greatly impaired ability to incorporate iron compared with chaotropically unfolded apoRds. In this study, formation of the Rd holoprotein (holoRd) on addition of iron to a structured, but iron-uptake incompetent apoRd was investigated in the presence of polystyrene nanoparticles (NP). In our rationale, hydrophobic contacts between apoRd and the NP surface would expose protein regions (including ligand cysteines) buried in the structured apoRd, allowing iron incorporation and folding to the native holoRd. Burial of the hydrophobic regions in the folded holoRd would allow its detachment from the NP surface. We found that both rate and yield of holoRd formation increased significantly in the presence of NP and were influenced by the NP concentration and size. Rates and yields had an optimum at "catalytic" NP concentrations (0.2 g/L NP) when using relatively small NP (46 nm diameter). At these optimal conditions, only a fraction of the apoRd was bound to the NP, consistent with the occurrence of turnover events on the NP surface. Lower rates and yields at higher NP concentrations or when using larger NP (200 nm) suggest that steric effects and molecular crowding on the NP surface favor specific "iron-uptake-competent" conformations of apoRd on the NP surface. This bio-mimetic chaperone system may be applicable to other proteins requiring an unfolding step before cofactor-triggered refolding, particularly when over-expressed under limited cofactor accessibility.

pH-dependent self-assembly of EAK16 peptides in the presence of a hydrophobic surface: coarse-grained molecular dynamics simulation

Self-assembly behavior of the three types of ionic peptide, EAK16, is studied in the presence of a hydrophobic surface using coarse-grained molecular dynamics simulations at three pH ranges of the solution. It is found that the peptide chains of all the three types assemble on the hydrophobic surface. EAK16-I and EAK16-II peptides assemble into ribbon-like structures, regardless of the value of pH. EAK16-IV peptide chains, however, assemble into ribbon-like structures at low and high pH ranges and form disc-shaped assemblies on the hydrophobic surface at the isoelectric point, pH = 7. Strong intra-chain electrostatic interactions in the case of EAK16-IV peptide play the main role in dependence of its self-assembly behavior on pH and the different morphology of its assembly relative to those of the two other types. Kinetics of growth of the assemblies on the hydrophobic surface is also studied.

How do proteins unfold upon adsorption on nanoparticle surfaces?

Owing to their many outstanding features, such as small size, large surface area, and cell penetration ability, nanoparticles have been increasingly used in medicine and biomaterials as drug carriers and diagnostic or therapeutic agents. However, our understanding of the interactions of biological entities, especially proteins, with nanoparticles is far behind the explosive development of nanotechnology. In typical protein-nanoparticle interactions, two processes (i.e., surface binding and conformational changes in proteins) are intermingled with each other and have not yet been quantitatively described. Here, by using a stopped-flow fast mixing technique, we were able to shed light on the kinetics of the adsorption-induced protein unfolding on nanoparticle surfaces in detail. We observed a biphasic denaturation behavior of protein GB1 on latex nanoparticle surfaces. Such kinetics can be adequately described by a fast equilibrium adsorption followed by a slow reversible unfolding of GB1. On the basis of this model, we quantitatively measured all rate constants that are involved in this process, from which the free-energy profile is constructed. This allows us to evaluate the effects of environmental factors, such as pH and ionic strength, on both the adsorption and the conformational change in GB1 on the latex nanoparticle surface. These studies provide a general physical picture of the adsorption-induced unfolding of proteins on nanoparticle surfaces and a quantitative description of the energetics of each transition. We anticipate that it will greatly advance our current understanding of protein-nanoparticle interactions and will be helpful for the rational control of such interactions in biomedical applications.

Elastic properties of a protein-polymer-grafted surface

Surfaces grafted with poly(methyl methacrylate) (PMMA) and streptavidin were synthesized through click chemistry to investigate the role of surface stiffness on protein adsorption as the hydrophilic and hydrophobic surface coverage of the substituents vary. Surface topographies coupled with the nanoindentation results indicated that, with the appropriate selections of polymer coverage and chain length, the extent of non-specific protein adhesion could be controlled by the hydrophobic interactions between PMMA, biotin, and streptavidin. It was shown that, when the molecular weight and stiffness of PMMA was close to that of streptavidin, patchy PMMA morphologies were obtained, which help inhibit the non-specific adsorption of streptavidin.

Efficient and Unbiased Sampling of Biomolecular Systems in the Canonical Ensemble: A Review of Self-Guided Langevin Dynamics

This review provides a comprehensive description of the self-guided Langevin dynamics (SGLD) and the self-guided molecular dynamics (SGMD) methods and their applications. Example systems are included to provide guidance on optimal application of these methods in simulation studies. SGMD/SGLD has enhanced ability to overcome energy barriers and accelerate rare events to affordable time scales. It has been demonstrated that with moderate parameters, SGLD can routinely cross energy barriers of 20 kT at a rate that molecular dynamics (MD) or Langevin dynamics (LD) crosses 10 kT barriers. The core of these methods is the use of local averages of forces and momenta in a direct manner that can preserve the canonical ensemble. The use of such local averages results in methods where low frequency motion "borrows" energy from high frequency degrees of freedom when a barrier is approached and then returns that excess energy after a barrier is crossed. This self-guiding effect also results in an accelerated diffusion to enhance conformational sampling efficiency. The resulting ensemble with SGLD deviates in a small way from the canonical ensemble, and that deviation can be corrected with either an on-the-fly or a post processing reweighting procedure that provides an excellent canonical ensemble for systems with a limited number of accelerated degrees of freedom. Since reweighting procedures are generally not size extensive, a newer method, SGLDfp, uses local averages of both momenta and forces to preserve the ensemble without reweighting. The SGLDfp approach is size extensive and can be used to accelerate low frequency motion in large systems, or in systems with explicit solvent where solvent diffusion is also to be enhanced. Since these methods are direct and straightforward, they can be used in conjunction with many other sampling methods or free energy methods by simply replacing the integration of degrees of freedom that are normally sampled by MD or LD.

Force-momentum-based self-guided Langevin dynamics: a rapid sampling method that approaches the canonical ensemble

The self-guided Langevin dynamics (SGLD) is a method to accelerate conformational searching. This method is unique in the way that it selectively enhances and suppresses molecular motions based on their frequency to accelerate conformational searching without modifying energy surfaces or raising temperatures. It has been applied to studies of many long time scale events, such as protein folding. Recent progress in the understanding of the conformational distribution in SGLD simulations makes SGLD also an accurate method for quantitative studies. The SGLD partition function provides a way to convert the SGLD conformational distribution to the canonical ensemble distribution and to calculate ensemble average properties through reweighting. Based on the SGLD partition function, this work presents a force-momentum-based self-guided Langevin dynamics (SGLDfp) simulation method to directly sample the canonical ensemble. This method includes interaction forces in its guiding force to compensate the perturbation caused by the momentum-based guiding force so that it can approximately sample the canonical ensemble. Using several example systems, we demonstrate that SGLDfp simulations can approximately maintain the canonical ensemble distribution and significantly accelerate conformational searching. With optimal parameters, SGLDfp and SGLD simulations can cross energy barriers of more than 15 kT and 20 kT, respectively, at similar rates for LD simulations to cross energy barriers of 10 kT. The SGLDfp method is size extensive and works well for large systems. For studies where preserving accessible conformational space is critical, such as free energy calculations and protein folding studies, SGLDfp is an efficient approach to search and sample the conformational space.

Molecular dynamics simulation of free and forced BSA adsorption on a hydrophobic graphite surface

The adsorption of bovine serum albumin (BSA) onto a hydrophobic graphite surface is studied using molecular-dynamics simulation. In addition to the free, that is, unsteered, adsorption, we also investigate forced adsorption, in which the action of an AFM tip pushing the protein with constant force to the surface is modeled. Using an implicit inviscid water model, the adsorption dynamics and energetics are monitored for two different initial protein orientations toward the surface. In all cases, we find that the protein partially unfolds and spreads on the surface. The spreading is in agreement with the well-known high biocompatibility of graphite-based implants. The denaturation is, however, greatly enhanced in the case of forced adsorption. We follow the position of the so-called lipid-binding pocket found in subdomain IIIA (Sudlow site II) during adsorption and find that it is tilted and moved toward the graphite surface in all cases, in agreement with its hydrophobic character. The relevance of our findings for the common measurement procedure of studying protein adhesion using AFM experiments is discussed.

Diffusion of hydrophobin proteins in solution and interactions with a graphite surface

Background

Hydrophobins are small proteins produced by filamentous fungi that have a variety of biological functions including coating of spores and surface adhesion. To accomplish these functions, they rely on unique interface-binding properties. Using atomic-detail implicit solvent rigid-body Brownian dynamics simulations, we studied the diffusion of HFBI, a class II hydrophobin from Trichoderma reesei, in aqueous solution in the presence and absence of a graphite surface.

Results

In the simulations, HFBI exists in solution as a mixture of monomers in equilibrium with different types of oligomers. The oligomerization state depends on the conformation of HFBI. When a Highly Ordered Pyrolytic Graphite (HOPG) layer is present in the simulated system, HFBI tends to interact with the HOPG layer through a hydrophobic patch on the protein.

Conclusions

From the simulations of HFBI solutions, we identify a tetrameric encounter complex stabilized by non-polar interactions between the aliphatic residues in the hydrophobic patch on HFBI. After the formation of the encounter complex, a local structural rearrangement at the protein interfaces is required to obtain the tetrameric arrangement seen in HFBI crystals. Simulations performed with the graphite surface show that, due to a combination of a geometric hindrance and the interaction of the aliphatic sidechains with the graphite layer, HFBI proteins tend to accumulate close to the hydrophobic surface.