Multiprotein interactions during surface adsorption: a molecular dynamics study of lysozyme aggregation at a charged solid surface

The Adsorption of P2X2 Receptors Interacting with IgG Antibodies Revealed by Combined AFM Imaging and Mechanical Simulation

The adsorption of proteins onto surfaces significantly impacts biomaterials, medical devices, and biological processes. This study aims to provide insights into the irreversible adsorption process of multiprotein complexes, particularly focusing on the interaction between anti-His6 IgG antibodies and the His6-tagged P2X2 receptor. Traditional approaches to understanding protein adsorption have centered around kinetic and thermodynamic models, often examining individual proteins and surface coverage, typically through Molecular Dynamics (MD) simulations. In this research, we introduce a computational approach employing Autodesk Maya 3D software for the investigation of multiprotein complexes' adsorption behavior. Utilizing Atomic Force Microscopy (AFM) imaging and Maya 3D-based mechanical simulations, our study yields real-time structural and kinetic observations. Our combined experimental and computational findings reveal that the P2X2 receptor-IgG antibody complex likely undergoes absorption in an 'extended' configuration. Whereas the P2X2 receptor is less adsorbed once is complexed to the IgG antibody compared to its individual state, the opposite is observed for the antibody. This insight enhances our understanding of the role of protein-protein interactions in the process of protein adsorption.

Interactions between glucosides of the tip of the S1 subunit of SARS-CoV-2 spike protein and dry and wet surfaces of CuO and Cu-A model for the surfaces of coinage metals

Despite their importance there is little knowledge at the atomic scale on the interactions between fragments of SARS-CoV-2 and inorganic materials. Such knowledge is important to understand the survival of the virus at surfaces and for the development of antiviral materials. Here is reported a study of the interactions between glucoside monomers of the tip of the S1 subunit of SARS-CoV-2 spike protein with dry and wet surfaces of CuO and Cu, performed with dispersion corrected density functional theory—DFT. The three glucoside monomers that constitute the tip of S1: 6VSB, 6VXX and 6X6P, were adsorbed onto dry and wet CuO(111) and Cu(110) with different orientations and surface alignments.

There are large differences—of up to 1.3 eV—in binding energies between these monomers and the surfaces. These differences depend on: the type of surface; if the surface is wet or dry; if the glucosidic O-atom points towards or away from the surfaces; and to a smaller extent on the surface alignment of the monomers. All monomers bind strongly to the surfaces via molecular adsorption that does not involve bond breaking in the monomers at this stage. 6VSB has the larger adsorption energies—that reach 2.2 eV—due to its larger dipole moment. Both materials bind the monomers more strongly when their surfaces are dry. At Cu(110) the bonds are on average 1 eV stronger when the surface is dry when compared to wet. The difference between dry and wet CuO(111) is smaller, in the order of 0.2 eV. Overall, it is here shown that the stability of the monomers of the tip of the spike protein of the virus is very different at different surfaces. For a given surface the larger binding energies in dry conditions could explain the differences in the surface stability of the spike protein depending on the presence of moisture.

Simulated and experimental force spectroscopy of lysozyme on silica

The force spectra of proteins detaching from oxide surfaces measured by atomic force microscopy (AFM) often present complex patterns of peaks, which are difficult to correlate with individual bond-breaking events at the atomic scale. In this work we rationalize experimental AFM force spectra of hen-egg-white lysozyme detaching from silica by means of all-atom steered molecular dynamics (SMD) simulations. In particular, we demonstrate that the native tertiary structure of lysozyme is preserved if, and only if, its four intramolecular disulfide bridges are intact. Otherwise, the protein pulled off the surface undergoes severe unfolding, which is well captured by SMD simulations in explicit solvent. Implicit solvent simulations, on the contrary, wrongly predict protein unfolding even in the presence of S-S bridges, due to the lack of additional structural stabilization provided by the water's hydrogen-bond network within and surrounding the protein. On the basis of our combined experimental and theoretical findings, we infer that the rugged force spectra characteristic of lysozyme/silica interfaces are not due to the successive breaking of internal disulfide bonds leading to partial unfolding events. Rather, they reflect the detachment of several molecules bound to the same AFM tip, each anchored to the surface via multiple hydrogen and ionic bonds.

Molecular dynamics simulations of lysozyme-lipid systems: probing the early steps of protein aggregation

Using the molecular dynamics simulation, the role of lipids in the lysozyme transition into the aggregation-competent conformation has been clarified. Analysis of the changes of lysozyme secondary structure upon its interactions with the model bilayer membranes composed of phosphatidylcholine and its mixtures with phosphatidylglycerol (10, 40, and 80 mol%) within the time interval of 100 ns showed that lipid-bound protein is characterized by the increased content of β-structures. Along with this, the formation of protein-lipid complexes was accompanied by the increase in the gyration radius and the decrease in RMSD of polypeptide chain. The results obtained were interpreted in terms of the partial unfolding of lysozyme molecule on the lipid matrix, with the magnitude of this effect being increased with increasing the fraction of anionic lipids. Based on the results of molecular dynamics simulation, a hypothetical model of the nucleation of lysozyme amyloid fibrils in a membrane environment was suggested.

Molecular Dynamics Simulation of a RNA Aptasensor

Single-stranded RNA aptamers have emerged as novel biosensor tools. However, the immobilization procedure of the aptamer onto a surface generally induces a loss of affinity. To understand this molecular process, we conducted a complete simulation study for the Flavin mononucleotide aptamer for which experimental data are available. Several molecular dynamics simulations (MD) of the Flavin in complex with its RNA aptamer were conducted in solution, linked with six thymidines (T6) and, finally, immobilized on an hexanol-thiol-functionalized gold surface. First, we demonstrated that our MD computations were able to reproduce the experimental solution structure and to provide a meaningful estimation of the Flavin free energy of binding. We also demonstrated that the T6 linkage, by itself, does not generate a perturbation of the Flavin recognition process. From the simulation of the complete biosensor system, we observed that the aptamer stays oriented parallel to the surface at a distance around 36 Å avoiding, this way, interaction with the surface. We evidenced a structural reorganization of the Flavin aptamer binding mode related to the loss of affinity and induced by an anisotropic distribution of sodium cationic densities. This means that ionic diffusion is different between the surface and the aptamer than above this last one. We suggest that these findings might be extrapolated to other nucleic acids systems for the future design of biosensors with higher efficiency and selectivity.

How Negatively Charged Proteins Adsorb to Negatively Charged Surfaces: A Molecular Dynamics Study of BSA Adsorption on Silica

How proteins adsorb to inorganic material surfaces is critically important for the development of new biotechnologies, since the orientation and structure of the adsorbed proteins impacts their functionality. While it is known that many negatively charged proteins readily adsorb to negatively charged oxide surfaces, a detailed understanding of how this process occurs is lacking. In this work we study the adsorption of BSA, an important transport protein that is negatively charged at physiological conditions, to a model silica surface that is also negatively charged. We use fully atomistic molecular dynamics to provide detailed understanding of the noncovalent interactions that bind the BSA to the silica surface. Our results provide new insight into the competing roles of long-range electrostatics and short-range forces, and the consequences this has for the orientation and structure of the adsorbed proteins.

Steering protein adsorption at charged surfaces: electric fields and ionic screening

Protein adsorption at charged surfaces is a common process in the development of functional technological devices.

When the Label Matters: Adsorption of Labeled and Unlabeled Proteins on Charged Surfaces

Fluorescent labels are often attached to proteins to monitor binding and adsorption processes. Docking simulations for native hen egg white lysozyme (HEWL) and HEWL labeled with fluorescein isothiocyanate show that these adsorb differently on charged surfaces. Attachment of even a small label can significantly change the interaction properties of a protein. Thus, the results of experiments with fluorescently labeled proteins should be interpreted by modeling the structures and computing the interaction properties of both labeled and unlabeled species.

Lysozyme adsorption at a silica surface using simulation and experiment: effects of pH on protein layer structure

Hen Egg White Lysozyme (HEWL) is a widely used exemplar to study protein adsorption on surfaces and interfaces. Here we use fully atomistic Molecular Dynamics (MD) simulations, Multi-Parametric Surface Plasmon Resonance (MP-SPR), contact angle and zeta potential measurements to study HEWL adsorption at a silica surface. The simulations provide a detailed description of the adsorption mechanism and indicate that at pH7 the main adsorption driving force is electrostatics, supplemented by weaker hydrophobic forces. Moreover, they reveal the preferred orientation of the adsorbed protein and show that its structure is only slightly altered at the interface with the surface. This provides the basis for interpreting the experimental results, which indicate the surface adsorbs a close-packed monolayer at about pH10 where the surface has a large negative zeta potential and the HEWL is positively charged. At higher pH, the adsorption amount of the protein layer is greatly reduced due to the loss of charge on the protein. At lower pH, the smaller zeta potential of the surface leads to lower HEWL adsorption. These interpretations are complemented by the contact angle measurements that show how the hydrophobicity of the surface is greatest when the surface coverage is highest. The simulations provide details of the hydrophobic residues exposed to solution by the adsorbed HEWL, completing the picture of the protein layer structure.

Role of arginine in mediating protein-carbon nanotube interactions

Arginine-rich proteins (e.g., lysozyme) or poly-L-arginine peptides have been suggested as solvating and dispersing agents for single-wall carbon nanotubes (CNTs) in water. In addition, protein structure-function in porous and hydrophobic materials is of broad interest. The amino acid residue, arginine (Arg(+)), has been implicated as an important mediator of protein/peptide-CNT interactions. To understand the structural and thermodynamic aspects of this interaction at the molecular level, we employ molecular dynamics (MD) simulations of the protein lysozyme in the interior of a CNT, as well as of free solutions of Arg(+) in the presence of a CNT. To dissect the Arg(+)-CNT interaction further, we also perform simulations of aqueous solutions of the guanidinium ion (Gdm(+)) and the norvaline (Nva) residue in the presence of a CNT. We show that the interactions of lysozyme with the CNT are mediated by the surface Arg(+) residues. The strong interaction of Arg(+) residue with the CNT is primarily driven by the favorable interactions of the Gdm(+) group with the CNT wall. The Gdm(+) group is not as well-hydrated on its flat sides, which binds to the CNT wall. This is consistent with a similar binding of Gdm(+) ions to a hydrophobic polymer. In contrast, the Nva residue, which lacks the Gdm(+) group, binds to the CNT weakly. We present details of the free energy of binding, molecular structure, and dynamics of these solutes on the CNT surface. Our results highlight the important role of Arg(+) residues in protein-CNT or protein-carbon-based material interactions. Such interactions could be manipulated precisely through protein engineering, thereby offering control over protein orientation and structure on CNTs, graphene, or other hydrophobic interfaces.

A bottom-up approach to understanding protein layer formation at solid-liquid interfaces

A common goal across different fields (e.g. separations, biosensors, biomaterials, pharmaceuticals) is to understand how protein behavior at solid-liquid interfaces is affected by environmental conditions. Temperature, pH, ionic strength, and the chemical and physical properties of the solid surface, among many factors, can control microscopic protein dynamics (e.g. adsorption, desorption, diffusion, aggregation) that contribute to macroscopic properties like time-dependent total protein surface coverage and protein structure. These relationships are typically studied through a top-down approach in which macroscopic observations are explained using analytical models that are based upon reasonable, but not universally true, simplifying assumptions about microscopic protein dynamics. Conclusions connecting microscopic dynamics to environmental factors can be heavily biased by potentially incorrect assumptions. In contrast, more complicated models avoid several of the common assumptions but require many parameters that have overlapping effects on predictions of macroscopic, average protein properties. Consequently, these models are poorly suited for the top-down approach. Because the sophistication incorporated into these models may ultimately prove essential to understanding interfacial protein behavior, this article proposes a bottom-up approach in which direct observations of microscopic protein dynamics specify parameters in complicated models, which then generate macroscopic predictions to compare with experiment. In this framework, single-molecule tracking has proven capable of making direct measurements of microscopic protein dynamics, but must be complemented by modeling to combine and extrapolate many independent microscopic observations to the macro-scale. The bottom-up approach is expected to better connect environmental factors to macroscopic protein behavior, thereby guiding rational choices that promote desirable protein behaviors.

Mesoscopic coarse-grained simulations of lysozyme adsorption

Coarse-grained simulations are adopted to study the adsorption behavior of lysozyme on different (hydrophobic, neutral hydrophilic, zwitterionic, negatively charged, and positively charged) surfaces at the mesoscopic microsecond time scale (1.2 μs). Simulation results indicate the following: (i) the conformation change of lysozyme on the hydrophobic surface is bigger than any other studied surfaces; (ii) the active sites of lysozyme are faced to the hydrophobic surface with a "top end-on" orientation, while they are exposed to the liquid phase on the hydrophilic surface with a "back-on" orientation; (iii) the neutral hydrophilic surface can induce the adsorption of lysozyme, while the nonspecific protein adsorption can be resisted by the zwitterionic surface; (iv) when the solution ionic strength is low, lysozyme can anchor on the negatively charged surface easily but cannot adsorb on the positively charged surface; (v) when the solution ionic strength is high, the positively charged lysozyme can also adsorb on the like-charged surface; (vi) the major positive potential center of lysozyme, especially the residue ARG128, plays a vital role in leading the adsorption of lysozyme on charged surfaces; (vii) when the ionic strength is high, a counterion layer is formed above the positively charged surface, which is the key factor why lysozyme can adsorb on a like-charged surface. The coarse-grained method based on the MARTINI force field for proteins and the BMW water model could provide an efficient way to understand protein interfacial adsorption behavior at a greater length scale and time scale.

Spontaneous membrane-translocating peptide adsorption at silica surfaces: a molecular dynamics study

Spontaneous membrane-translocating peptides (SMTPs) have recently been shown to directly penetrate cell membranes. Adsorption of a SMTP, and some engineered extensions, at model silica surfaces is studied herein using fully atomistic molecular dynamics simulations in order to assess their potential to construct novel drug delivery systems. The simulations are designed to reproduce the electric fields above single, siloxide-rich charged surfaces, and the trajectories indicate that the main driving force for adsorption is electrostatic. An increase in the salt concentration slows down but does not prevent adsorption of the SMTP to the surface; it also does not result in peptide desorption, suggesting additional binding via hydrophobic forces. The results are used to design extensions to the peptide sequence which we find enhance adsorption but do not affect the adsorbed conformation. We also investigate the effect of surface hydroxylation on the peptide adsorption. In all cases, the final adsorbed conformations are with the peptide flattened to the surface with arginine residues, which are key to the peptide’s function, anchoring it to the surface so that they are not exposed to solution. This conformation could impact their role in membrane translocation and thus has important implications for the design of future drug delivery vehicles.

Protein diffusion and long-term adsorption states at charged solid surfaces

The diffusion pathways of lysozyme adsorbed to a model charged ionic surface are studied using fully atomistic steered molecular dynamics simulation. The simulations start from existing protein adsorption trajectories, where it has been found that one particular residue, Arg128 at the N,C-terminal face, plays a crucial role in anchoring the lysozyme to the surface [Langmuir 2010 , 26 , 15954 - 15965]. We first investigate the desorption pathway for the protein by pulling the Arg128 side chain away from the surface in the normal direction, and its subsequent readsorption, before studying diffusion pathways by pulling the Arg128 side chain parallel to the surface. We find that the orientation of this side chain plays a decisive role in the diffusion process. Initially, it is oriented normal to the surface, aligning in the electrostatic field of the surface during the adsorption process, but after resorption it lies parallel to the surface, being unable to return to its original orientation due to geometric constraints arising from structured water layers at the surface. Diffusion from this alternative adsorption state has a lower energy barrier of ∼0.9 eV, associated with breaking hydrogen bonds along the pathway, in reasonable agreement with the barrier inferred from previous experimental observation of lysozyme surface clustering. These results show the importance of studying protein diffusion alongside adsorption to gain full insight into the formation of protein clusters and films, essential steps in the future development of functionalized surfaces.

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.

Interaction of biological molecules with clay minerals: a combined spectroscopic and sorption study of lysozyme on saponite

The interaction of hen egg white lysozyme (HEWL) with Na- and Cs-exchanged saponite was investigated using sorption, structural, and spectroscopic methods as a model system to study clay-protein interactions. HEWL sorption to Na- and Cs-saponite was determined using the bicinchoninic acid (BCA) assay, thermogravimetric analysis, and C and N analysis. For Na-saponite, the TGA and elemental analysis-derived sorption maximum was 600 mg/g corresponding to a surface coverage of 0.85 ng/mm(2) with HEWL occupying 526 m(2)/g based on a cross-sectional area of 13.5 nm(2)/molecule. HEWL sorption on Na-saponite was accompanied by the release of 9.5 Na(+) ions for every molecule of HEWL sorbed consistent with an ion exchange mechanism between the positively charged HEWL (IEP 11) and the negatively charged saponite surface. The d-spacing of the HEWL-Na-saponite complex increased to a value of 4.4 nm consistent with the crystallographic dimensions of HEWL of 3 × 3 × 4.5 nm. In the case of Cs-saponite, there was no evidence of interlayer sorption; however, sorption of HEWL to the "external" surface of Cs-saponite showed a high affinity isotherm. FTIR and Raman analysis of the amide I region of the HEWL-saponite films prepared from water and D(2)O showed little perturbation to the secondary structure of the protein. The overall hydrophilic nature of the HEWL-Na-saponite complex was determined by water vapor sorption measurements. The clay retained its hydrophilic character with a water content of 18% at high humidity corresponding to 240 H(2)O molecules per molecule of HEWL.