Molecular simulation of bovine beta-lactoglobulin adsorbed onto a positively charged solid surface

Molecular Modeling of the Adsorption of an Egg Yolk Protein on a Water-Oil Interface

Egg yolk contains several molecular species with emulsifying properties, such as proteins and phospholipids. In particular, these molecules have both polar and nonpolar parts and thus can act as surfactants. One of the most surface-active proteins from egg yolk low-density lipoproteins is the so-called Apovitellenin-1. Experimental studies have been hindered by difficulties in isolating individual species from egg yolk lipoproteins. The purpose of this work was to assess the emulsifying properties of Apovitellenin-1 and any potential cooperative or competitive behavior in the presence of phospholipids. To do so, molecular simulations were carried out in a liquid-liquid interfacial system consisting of water and soybean oil, with varying concentrations of phospholipids and for different spatial configurations. To evaluate the conformational stability of the protein at the water-oil interface, the Gibbs free energy was computed from Metadynamics simulations as a function of the distance from the interface and of the radius of gyration. Moreover, a detailed analysis was also performed to determine which peptide residues were responsible for the protein adsorption at the oil-water interface as well as the lowering of the interfacial tension. Lastly, we combined the simulation results with a thermodynamic model to predict the interfacial tension behavior at increasing protein bulk concentration, which cannot be measured experimentally.

Effect of substrate charge density on the adsorption of intrinsically disordered protein amyloid β40: a molecular dynamics study

The inhibitory effect of negatively charged gold nanoparticles (AuNPs) on amyloidogenic protein fibrillation has been established from experiments and computer simulations. Here, we investigate the effect of the charge density (σ) of gold (Au) surfaces on the adsorption of the intrinsically disordered amyloid β40 (Aβ40) monomer using molecular dynamics (MD) simulations. On the basis of the binding free energy, some key residues (ARG5, LYS16, LYS28, LEU17-ALA21, ILE31-VAL38) were found to be responsible for preventing the β-sheet formation, which is known to be a precursor for fibrillation. Until a critical charge density (σ_c) of -0.167 e nm^-2, the key residues remained adsorbed on the Au slab. A saturation in the number of condensed counterions (Na⁺) on Aβ40 was also observed at σ_c. Beyond σ_c, the condensation of Na⁺ occurs only on the Au slab, leading to competition between positively charged key residues and condensed ions. This competition was found to be responsible for the lack of adsorption of the key residues, leading to β-sheet formation for σ > -0.167 e nm^-2. This study suggests that if the key residues are not adsorbed, then β-sheet formation is observed, which can then lead to the development of proto-fibrils and subsequently fibrillation. Therefore the surface should have an optimal charge density to be an effective inhibitor of fibrillation.

Competitive protein adsorption on charge regulating silica-like surfaces: the role of protonation equilibrium

We develop a molecular thermodynamic theory to study the interaction of some proteins with a charge regulating silica-like surface under a wide range of conditions, including pH, salt concentration and protein concentration. Proteins are modeled using their three dimensional structure from crystallographic data and the average experimental pKa of amino acid residues. As model systems, we study single-protein and binary solutions of cytochrome c, green fluorescent protein, lysozyme and myoglobin. Our results show that protonation equilibrium plays a critical role in the interactions of proteins with these type of surfaces. The terminal hydroxyl groups on the surface display considerable extent of charge regulation; protein residues with titratable side chains increase protonation according to changes in the local environment and the drop in pH near the surface. This behavior defines protein-surface interactions and leads to the emergence of several phenomena: (i) a complex non-ideal surface charge behavior; (ii) a non-monotonic adsorption of proteins as a function of pH; and (iii) the presence of two spatial regions, a protein-rich and a protein-depleted layer, that occur simultaneously at different distances from the surface when pH is slightly above the isoelectric point of the protein. In binary mixtures, protein adsorption and surface-protein interactions cannot be predicted from single-protein solution considerations.

Adsorption of lysozyme into a charged confining pore

Several applications arise from the confinement of proteins on surfaces because their stability and biological activity are enhanced. It is also known that the way in which a protein adsorbs on the surface is important for its biological function since its active sites should not be obstructed. In this study, the adsorption properties of hen egg-white lysozyme, HEWL, into a negatively charged silica pore is examined by employing a coarse-grained model and constant-pH Monte Carlo simulations. The role of electrostatic interactions is taken into account via including the Debye-Hückel potentials into the Cα structure-based model. We evaluate the effects of pH, salt concentration, and pore radius on the protein preferential orientation and spatial distribution of its residues regarding the pore surface. By mapping the residues that stay closer to the pore surface, we find that the increase of pH leads to orientational changes of the adsorbed protein when the solution pH gets closer to the HEWL isoelectric point. Under these conditions, the pK_a shift of these important residues caused by the adsorption into the charged confining surface results in a HEWL charge distribution that stabilizes the adsorption in the observed protein orientation. We compare our observations to the results of the pK_a shift for HEWL available in the literature and to some experimental data.

Adsorption of the amyloid β40 monomer on charged gold nanoparticles and slabs: a molecular dynamics study

Negatively charged nanoparticles are known to inhibit the fibrillation of amyloidogenic protein amyloid β (Aβ40), though the overall charge on the protein is negative. In this work a molecular dynamics study is reported to investigate the interaction of Aβ40 on negatively charged gold nanoparticles (3-5 nm) and charged (positive and negative) and neutral gold slabs. The equilibrium structures of Aβ40 on gold surfaces are characterized using residue-specific contacts on the gold surface, secondary structure analysis and binding free energy calculations. The simulation results reveal that the Aβ40 protein in water interconverts into β-sheets, which are building blocks of the mature fibrils, whereas on gold nanoparticles Aβ40 unfolds and adsorbs. Both the negatively charged gold nanoparticles and gold slabs arrest the formation of β-sheets in Aβ40, whereas the positively charged gold slab does not inhibit the formation of β-sheets. The residue-specific interactions between Aβ40 and the gold surfaces are important in governing the adsorption of Aβ40 on charged surfaces.

Biomolecular Binding at Aqueous Interfaces of Langmuir Monolayers of Bioconjugated Amphiphilic Mesogenic Molecules: A Molecular Dynamics Study

We report a molecular dynamics (MD) simulation study of protein binding at the aqueous-liquid crystal (LC) interfaces of bioconjugated mesogenic molecules. As a simple model of these interfaces, we use monolayers composed of biotin-conjugated or biotin-free amphiphilic mesogenic molecules and streptavidin in water. The all-atom MD simulations reveal that the binding of streptavidin to the biotin mesogenic monolayer is significantly stronger than that to biotin-free mesogenic monolayers. Although specific protein binding marginally increases the overall orientational order and the tilt of the biotin-conjugated mesogenic molecules of the monolayer, significant changes in tilt were observed near the bound protein (in contrast to the protein interaction with the monolayer without biotin). We also observe that specific protein binding changes the dynamic properties of the mesogens within the monolayer (e.g., lateral diffusion coefficients) and associated water. Overall, these MD simulations advance our understanding of the molecular-level phenomena involved in the binding of biomolecules and subsequent dynamic changes at the aqueous-LC interfaces. These results provide guidance to future molecular-level designs of biofunctional LC interfaces.

Impact of Silane Monolayers on the Adsorption of Streptavidin on Silica and Its Subsequent Interactions with Biotin: Molecular Dynamics and Steered Molecular Dynamics Simulations

Protein adsorption on surfaces is used in analytical tools as an immobilization mean to trap the analyte to be detected. However, protein adsorption can lead to a conformational change in the protein structure, resulting in a loss of bioactivity. Here, we study adsorption of a streptavidin-biotin complex on amorphous SiO₂ surfaces functionalized with five different silane self-assembled monolayers by all-atom molecular dynamics simulations. We find that the streptavidin global conformational change, as well as the nature of residues with high mobility, depends on the alkyl chain length and head-group charge of silane molecules. Effects on interactions with biotin are further investigated by steered molecular dynamics (SMD) simulations, which mimics atomic force microscopy (AFM) with the biotin attached on the tip. We show the combined effects of adsorption-induced global conformational changes and of the position of residues with high mobility on the streptavidin-biotin rupture force. By comparing our results to experimental and SMD rupture forces obtained in water, without any surface, we conclude that silane with uncharged and short alkyl chains allows streptavidin immobilization, while keeping biotin interactions better than silanes with long alkyl chains or charged head groups.

Adsorption of Milk Proteins (β-Casein and β-Lactoglobulin) and BSA onto Hydrophobic Surfaces

Here, we study films of proteins over planar surfaces and protein-coated microspheres obtained from the adsorption of three different proteins ( β -casein, β -lactoglobulin and bovine serum albumin (BSA)). The investigation of protein films in planar surfaces is performed by combining quartz crystal microbalance (QCM) and atomic force microscopy (AFM) measurements with all-atomic molecular dynamics (MD) simulations. We found that BSA and β -lactoglobulin form compact monolayers, almost without interstices between the proteins. However, β -casein adsorbs forming multilayers. The study of the electrokinetic mobility of protein-coated latex microspheres shows substantial condensation of ions from the buffer over the complexes, as predicted from ion condensation theories. The electrokinetic behavior of the latex-protein complexes is dominated by the charge of the proteins and the phenomenon of ion condensation, whereas the charge of the latex colloids plays only a minor role.

Modeling and simulation of protein-surface interactions: achievements and challenges

Understanding protein-inorganic surface interactions is central to the rational design of new tools in biomaterial sciences, nanobiotechnology and nanomedicine. Although a significant amount of experimental research on protein adsorption onto solid substrates has been reported, many aspects of the recognition and interaction mechanisms of biomolecules and inorganic surfaces are still unclear. Theoretical modeling and simulations provide complementary approaches for experimental studies, and they have been applied for exploring protein-surface binding mechanisms, the determinants of binding specificity towards different surfaces, as well as the thermodynamics and kinetics of adsorption. Although the general computational approaches employed to study the dynamics of proteins and materials are similar, the models and force-fields (FFs) used for describing the physical properties and interactions of material surfaces and biological molecules differ. In particular, FF and water models designed for use in biomolecular simulations are often not directly transferable to surface simulations and vice versa. The adsorption events span a wide range of time- and length-scales that vary from nanoseconds to days, and from nanometers to micrometers, respectively, rendering the use of multi-scale approaches unavoidable. Further, changes in the atomic structure of material surfaces that can lead to surface reconstruction, and in the structure of proteins that can result in complete denaturation of the adsorbed molecules, can create many intermediate structural and energetic states that complicate sampling. In this review, we address the challenges posed to theoretical and computational methods in achieving accurate descriptions of the physical, chemical and mechanical properties of protein-surface systems. In this context, we discuss the applicability of different modeling and simulation techniques ranging from quantum mechanics through all-atom molecular mechanics to coarse-grained approaches. We examine uses of different sampling methods, as well as free energy calculations. Furthermore, we review computational studies of protein-surface interactions and discuss the successes and limitations of current approaches.

Deciphering β-Lactoglobulin Interactions at an Oil-Water Interface: A Molecular Dynamics Study

Protein adsorption at liquid-liquid interfaces is of immense relevance to many biological processes and dairy-based functional foods. Due to experimental limitations, however, there is still a remarkable lack of understanding of the adsorption mechanism, particularly at a molecular level. In this study, atomistic molecular dynamics simulations were used to elucidate the approach and adsorption mechanism of β-lactoglobulin (β-LG) at a decane-water interface. Through multiple independent simulations starting from three representative initial orientations of β-LG relative to the decane surface the rate at which β-LG approaches the oil/water interface is found to be independent of its initial orientation, and largely stochastic in nature. While the residues that first make contact with the decane and the final orientation of β-LG upon adsorption are similar in all cases, the adsorption process is driven predominantly by structural rearrangements that preserve the secondary structure but expose hydrophobic residues to the decane surface. This detailed characterization of the adsorption of β-LG at an oil/water interface should inform the design and development of novel encapsulation and delivery systems in the food and pharmaceutical sciences.

Interactions and diffusion in fine-stranded β-lactoglobulin gels determined via FRAP and binding

The effects of electrostatic interactions and obstruction by the microstructure on probe diffusion were determined in positively charged hydrogels. Probe diffusion in fine-stranded gels and solutions of β-lactoglobulin at pH 3.5 was determined using fluorescence recovery after photobleaching (FRAP) and binding, which is widely used in biophysics. The microstructures of the β-lactoglobulin gels were characterized using transmission electron microscopy. The effects of probe size and charge (negatively charged Na2-fluorescein (376Da) and weakly anionic 70kDa FITC-dextran), probe concentration (50 to 200 ppm), and β-lactoglobulin concentration (9% to 12% w/w) on the diffusion properties and the electrostatic interaction between the negatively charged probes and the positively charged gels or solutions were evaluated. The results show that the diffusion of negatively charged Na2-fluorescein is strongly influenced by electrostatic interactions in the positively charged β-lactoglobulin systems. A linear relationship between the pseudo-on binding rate constant and the β-lactoglobulin concentration for three different probe concentrations was found. This validates an important assumption of existing biophysical FRAP and binding models, namely that the pseudo-on binding rate constant equals the product of the molecular binding rate constant and the concentration of the free binding sites. Indicators were established to clarify whether FRAP data should be analyzed using a binding-diffusion model or an obstruction-diffusion model.

Multiscale simulations of protein G B1 adsorbed on charged self-assembled monolayers

The orientation of an antibody plays an important role in the development of immunosensors. Protein G is an antibody binding protein, which specifically targets the Fc fragment of an antibody. In this work, the orientation of prototypical and mutated protein G B1 adsorbed on positively and negatively charged self-assembled monolayers was studied by parallel tempering Monte Carlo and all-atom molecular dynamics simulations. Both methods present generally similar orientation distributions of protein G B1 for each kind of surface. The root-mean-square deviation, DSSP, gyration radius, eccentricity, dipole moment, and superimposed structures of protein G B1 were analyzed. Moreover, the orientation of binding antibody was also predicted in this work. Simulation results show that with the same orientation trends, the mutant exhibits narrower orientation distributions than does the prototype, which was mainly caused by the stronger dipole of the mutant. Both kinds of proteins adsorbed on charged surfaces were induced by the competition of electrostatic interaction and vdW interaction; the electrostatic interaction energy dominated the adsorption behavior. The protein adsorption was also largely affected by the distribution of charged residues within the proteins. Thus, the prototype could adsorb on a negatively charged surface, although it keeps a net charge of -4 e. The mutant has imperfect opposite orientation when it adsorbed on oppositely charged surfaces. For the mutant on a carboxyl-functionalized self-assembled monolayer (COOH-SAM), the orientation was the same as that inferred by experiments. While for the mutant on amine-functionalized self-assembled monolayer (NH2-SAM), the orientation was induced by the competition between attractive interactions (led by ASP40 and GLU56) and repulsive interactions (led by LYS10); thus, the perfect opposite orientation could not be obtained. On both surfaces, the adsorbed protein could retain its native conformation. The desired orientation of protein G B1, which would increase the efficiency of binding antibodies, could be obtained on a negatively charged surface adsorbed with the prototype. Further, we deduced that with the packing density of 12,076 protein G B1 domain per μm(2), the efficiency of the binding IgG would be maximized. The simulation results could be applied to control the orientation of protein G B1 in experiments and to provide a better understanding to maximize the efficiency of antibody binding.

Molecular simulation of fibronectin adsorption onto polyurethane surfaces

Poly(ethylene glycol)-based polyurethanes have been widely used in biomedical applications; however, they are prone to swelling. A natural polyol, castor oil, can be incorporated into these polyurethanes to control the degree of the swelling, which alters mechanical properties and protein adsorption characteristic of the polymers. In this work, we modeled poly(ethylene glycol) and castor oil copolymers of hexamethylene diisocyanate-based polyurethanes (PEG-HDI and CO-HDI, respectively) and compared their mechanisms for fibronectin adsorption using molecular mechanics and molecular dynamics simulations. Results showed that the interplay between the hydrophobic residues concentrated at the N-terminal end of the protein, the surface roughness, and the hydrophilicity of the polymer surface determine the overall protein adsorption affinity. Incorporating explicit water molecules in the simulations results in higher affinity for fibronectin adsorption to more hydrophobic surface of CO-HDI surfaces, emphasizing the role that water molecules play during adsorption. We also observed that the strain energies that are indicative of flexibility and consequently entropy are significantly affected by the changes in the patterns of β-sheet formation/breaking. Our study lends supports to the view that while castor oil controls the degree of swelling, it increases the adsorption of fibronectin to a limited extent due to the interplay between its hydrophobicity and its surface roughness, which needs to be taken into account during the design of polyurethane-based biomaterials.

Protein-decorated reduced oxide graphene composite and its application to SERS

A globular protein, β-lactoglobulin (BLG), was used to decorate reduced graphene oxide sheets (RGO) and the obtained BLG-RGO composite can be dispersed in aqueous solution with pH-sensitive solubility. The morphology of the BLG-RGO composite was studied by transmission electron microscopy (TEM) and atomic force microscopy (AFM). The results indicate that BLG-RGO is effectively exfoliated with an average thickness of 2.5 nm. UV-vis spectra were performed to examine the reduction degree and determine the optimum concentration of β-lactoglobulin and appropriate pH value. Furthermore, Raman spectra demonstrate that β-Lactoglobulin promotes the chemical reduction process of graphene oxide and benefits to repair the crystal defects. Due to the adsorption of β-Lactoglobulin on the surface of graphene sheets, the BLG-RGO composite was further used as template for Au nanoparticles assembly. These Au nanoparticles assembled on the BLG-RGO composite were shown to yield a large SERS enhancement for Rhodamine 6G.

Structural rearrangement of β-lactoglobulin at different oil-water interfaces and its effect on emulsion stability

Understanding the factors that control protein structure and stability at the oil-water interface continues to be a major focus to optimize the formulation of protein-stabilized emulsions. In this study, a combination of synchrotron radiation circular dichroism spectroscopy, front-face fluorescence spectroscopy, and dual polarization interferometry (DPI) was used to characterize the conformation and geometric structure of β-lactoglobulin (β-Lg) upon adsorption to two oil-water interfaces: a hexadecane-water interface and a tricaprylin-water interface. The results show that, upon adsorption to both oil-water interfaces, β-Lg went through a β-sheet to α-helix transition with a corresponding loss of its globular tertiary structure. The degree of conformational change was also a function of the oil phase polarity. The hexadecane oil induced a much higher degree of non-native α-helix compared to the tricaprylin oil. In contrast to the β-Lg conformation in solution, the non-native α-helical-rich conformation of β-Lg at the interface was resistant to further conformational change upon heating. DPI measurements suggest that β-Lg formed a thin dense layer at emulsion droplet surfaces. The effects of high temperature and the presence of salt on these β-Lg emulsions were then investigated by monitoring changes in the ζ-potential and particle size. In the absence of salt, high electrostatic repulsion meant β-Lg-stabilized emulsions were resistant to heating to 90 °C. Adding salt (120 mM NaCl) before or after heating led to emulsion flocculation due to the screening of the electrostatic repulsion between colloidal particles. This study has provided insight into the structural properties of proteins adsorbed at the oil-water interface and has implications in the formulation and production of emulsions stabilized by globular proteins.

Stability of two-dimensional growth of a packed body of proteins on a solid surface

Adsorption of proteins from the bulk is at times accompanied by a rearrangement which leads to the formation of closed packed bodies, that may or may not be crystalline. Mass transfer of protein molecules on a surface is modeled. Forced diffusion by van der Waals and electrostatic forces leads to segregation, which is eventually a different phase that is assumed to be thermodynamically favored. The net effective force in two-dimensions has been modeled approximately and shown to be much stronger and more long ranged than in the bulk: that is, under the same conditions, the protein molecules may not aggregate in the bulk they may aggregate on a surface. These forces have been used only indirectly but equivalently as an adsorption-desorption step at the interline. Eventually, a linear stability analysis of the growing body shows it to be unstable and would give rise to whiskers that are one molecule thick. This is what is observed experimentally. The conditions that give rise to the instability have been determined. The reverse case of rinsing of the protein molecules has also been studied experimentally and has been analyzed using the same mechanisms. Here it is seen that thicker inroads into the packed body cause the interline to take on a spongy appearance. It is conjectured that eventually islands will appear as seen in the experiments.

Understanding protein adsorption phenomena at solid surfaces

Protein adsorption at solid surfaces plays a key role in many natural processes and has therefore promoted a widespread interest in many research areas. Despite considerable progress in this field there are still widely differing and even contradictive opinions on how to explain the frequently observed phenomena such as structural rearrangements, cooperative adsorption, overshooting adsorption kinetics, or protein aggregation. In this review recent achievements and new perspectives on protein adsorption processes are comprehensively discussed. The main focus is put on commonly postulated mechanistic aspects and their translation into mathematical concepts and model descriptions. Relevant experimental and computational strategies to practically approach the field of protein adsorption mechanisms and their impact on current successes are outlined.

Adsorption mechanism of ribosomal protein L2 onto a silica surface: a molecular dynamics simulation study

A large-scale molecular dynamics simulation was carried out in order to investigate the adsorption mechanism of ribosomal protein L2 (RPL2) onto a silica surface at various pH values. RPL2 is a constituent protein of the 50S large ribosomal subunit, and a recent experimental report showed that it adsorbs strongly to silica surfaces and that it can be used to immobilize proteins on silica surfaces. The simulation results show that RPL2, especially domains 1 (residues 1-60) and 3 (residues 203-273), adsorbed more tightly to the silica surface above pH 7. We found that a major driving force for the adsorption of RPL2 onto the silica surface is the electrostatic interaction and that the structural flexibility of domains 1 and 3 may further contribute to the high affinity.