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

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.

Is It Possible to Find an Antimicrobial Peptide That Passes the Membrane Bilayer with Minimal Force Resistance? An Attempt at a Predictive Approach by Molecular Dynamics Simulation

There is still no answer to the mechanism of penetration of AMP peptides through the membrane bilayer. Several mechanisms for such a process have been proposed. It is necessary to understand whether it is possible, using the molecular dynamics method, to determine the ability of peptides of different compositions and lengths to pass through a membrane bilayer. To explain the passage of a peptide through a membrane bilayer, a method for preparing a membrane phospholipid bilayer was proposed, and 656 steered molecular dynamics calculations were carried out for pulling 7 amyloidogenic peptides with antimicrobial potential, and monopeptides (homo-repeats consisting of 10 residues of the same amino acid: Poly (Ala), Poly (Leu), Poly (Met), Poly (Arg), and Poly (Glu)) with various sequences through the membrane. Among the 15 studied peptides, the peptides exhibiting the least force resistance when passing through the bilayer were found, and the maximum reaction occurred at the boundary of the membrane bilayer entry. We found that the best correlation between the maximum membrane reaction force and the calculated parameters corresponds to the instability index (the correlation coefficient is above 0.9). One of the interesting results of this study is that the 10 residue amyloidogenic peptides and their extended peptides, with nine added residue cell-penetrating peptides and four residue linkers, both with established antimicrobial activity, have the same bilayer resistance force. All calculated data are summarized and posted on the server.

Functionalisation of Inorganic Material Surfaces with Staphylococcus Protein A: A Molecular Dynamics Study

Staphylococcus protein A (SpA) is found in the cell wall of Staphylococcus aureus bacteria. Its ability to bind to the constant Fc regions of antibodies means it is useful for antibody extraction, and further integration with inorganic materials can lead to the development of diagnostics and therapeutics. We have investigated the adsorption of SpA on inorganic surface models such as experimentally relevant negatively charged silica, as well as positively charged and neutral surfaces, by use of fully atomistic molecular dynamics simulations. We have found that SpA, which is itself negatively charged at pH7, is able to adsorb on all our surface models. However, adsorption on charged surfaces is more specific in terms of protein orientation compared to a neutral Au (111) surface, while the protein structure is generally well maintained in all cases. The results indicate that SpA adsorption is optimal on the siloxide-rich silica surface, which is negative at pH7 since this keeps the Fc binding regions free to interact with other species in solution. Due to the dominant role of electrostatics, the results are transferable to other inorganic materials and pave the way for new diagnostic and therapeutic designs where SpA might be used to conjugate antibodies to nanoparticles.

Simulating Peptide Monolayer Formation: GnRH-I on Silica

Molecular dynamics (MD) simulations can provide a detailed view of molecule behaviour at an atomic level, which can be useful when attempting to interpret experiments or design new systems. The decapeptide gonadotrophin-releasing hormone I (GnRH-I) is known to control fertility in mammals for both sexes. It was previously shown that inoculation with silica nanoparticles (SiNPs) coated with GnRH-I makes an effective anti-fertility vaccine due to how the peptide adsorbs to the nanoparticle and is presented to the immune system. In this paper, we develop and employ a protocol to simulate the development of a GnRH-I peptide adlayer by allowing peptides to diffuse and adsorb in a staged series of trajectories. The peptides start the simulation in an immobile state in solution above the model silica surface, and are then released sequentially. This facile approach allows the adlayer to develop in a natural manner and appears to be quite versatile. We find that the GnRH-I adlayer tends to be sparse, with electrostatics dominating the interactions. The peptides are collapsed to the surface and are seemingly free to interact with additional solutes, supporting the interpretations of the GNRH-I/SiNP vaccine system.

Water Bridges in Clay Nanopores: Mechanisms of Formation and Impact on Hydrocarbon Transport

Clays are prevalent in the earth's crust and usually deposited in the presence of water. An unusual finding in clays is that under certain conditions, water molecules can collectively form a bridge across a clay-hosted pore. However, there are relatively few studies focused on the formation mechanism of the water bridge in clay nanopores. In this work, we use molecular dynamics simulations to investigate the formation of the water bridge and its influence on fluid transport in slit-shaped illite nanopores. Two different basal illite surface chemistries are constructed: potassium-hydroxyl (P-H) and hydroxyl-hydroxyl (H-H) structures. Because pore size and water concentration are expected to control the formation of the water bridge, our simulations span a wide range of pore sizes and water concentrations. Generally, positive potassium layers and negative hydroxyl groups in P-H nanopore can induce partial charges which in return produce instant and local electric fields, favoring the formation of the water bridge. In P-H nanopores, the water bridge happens at a relatively low water concentration. However, in H-H nanopores, the water bridge only forms at high water concentrations. Additionally, smaller pore sizes favor the formation of water bridges. However, the presence of an electric field promotes the formation of a water bridge even in larger pore sizes in P-H pores. The results also indicate that in both P-H and H-H nanopores, water adsorption films initially create a smooth surface to promote the hydrocarbon flow. In P-H nanopores, further increases in the water concentration causes a sharp decline in the self-diffusion coefficients of the hydrocarbon and water due to the formation of the water bridge. The presence of electric fields in P-H pores can however weaken the confinement effect of illite and promote the hydrocarbon flow. In contrast, in H-H nanopores, the self-diffusion coefficients decline slowly with the increase of water concentration. This is because no water bridge is formed at low water concentrations in H-H nanopores.

Rationalising drug delivery using nanoparticles: a combined simulation and immunology study of GnRH adsorbed to silica nanoparticles

Silica nanoparticles (SiNPs) have been shown to have significant potential for drug delivery and as adjuvants for vaccines. We have simulated the adsorption of GnRH-I (gonadotrophin releasing hormone I) and a cysteine-tagged modification (cys-GnRH-I) to model silica surfaces, as well as its conjugation to the widely-used carrier protein bovine serum albumin (BSA). Our subsequent immunological studies revealed no significant antibody production was caused by the peptide-SiNP systems, indicating that the treatment was not effective. However, the testosterone response with the native peptide-SiNPs indicated a drug effect not found with cys-GnRH-I-SiNPs; this behaviour is explained by the specific orientation of the peptides at the silica surface found in the simulations. With the BSA systems, we found significant testosterone reduction, particularly for the BSA-native conjugates, and an antibody response that was notably higher with the SiNPs acting as an adjuvant; this behaviour again correlates well with the epitope presentation predicted by the simulations. The range of immunological and hormone response can therefore be interpreted and understood by the simulation results and the presentation of the peptides to solution, paving the way for the future rational design of drug delivery and vaccine systems guided by biomolecular simulation.

Adsorption of Fibronectin Fragment on Surfaces Using Fully Atomistic Molecular Dynamics Simulations

The effect of surface chemistry on the adsorption characteristics of a fibronectin fragment (FNIII^8⁻10) was investigated using fully atomistic molecular dynamics simulations. Model surfaces were constructed to replicate self-assembled monolayers terminated with methyl, hydroxyl, amine, and carboxyl moieties. It was found that adsorption of FNIII^8⁻10 on charged surfaces is rapid, specific, and driven by electrostatic interactions, and that the anchoring residues are either polar uncharged or of opposing charge to that of the targeted surfaces. On charged surfaces the presence of a strongly bound layer of water molecules and ions hinders FNIII^8⁻10 adsorption. In contrast, adsorption kinetics on uncharged surfaces are slow and non-specific, as they are driven by van der Waals interactions, and the anchoring residues are polar uncharged. Due to existence of a positively charged area around its cell-binding region, FNIII^8⁻10 is available for subsequent cell binding when adsorbed on a positively charged surface, but not when adsorbed on a negatively charged surface. On uncharged surfaces, the availability of the fibronectin fragment's cell-binding region is not clearly distinguished because adsorption is much less specific.

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.

Bovine Serum Albumin Adsorption at a Silica Surface Explored by Simulation and Experiment

Molecular details of BSA adsorption on a silica surface are revealed by fully atomistic molecular dynamics (MD) simulations (with a 0.5 μs trajectory), supported by dynamic light scattering (DLS), zeta potential, multiparametric surface plasmon resonance (MP-SPR), and contact angle experiments. The experimental and theoretical methods complement one another and lead to a wider understanding of the mechanism of BSA adsorption across a range of pH 3-9. The MD results show how the negatively charged BSA at pH7 adsorbs to the negatively charged silica surface, and reveal a unique orientation with preserved secondary and tertiary structure. The experiments then show that the protein forms complete monolayers at ∼ pH6, just above the protein's isoelectric point (pH5.1). The surface contact angle is maximum when it is completely coated with protein, and the hydrophobicity of the surface is understood in terms of the simulated protein conformation. The adsorption behavior at higher pH > 6 is also consistently interpreted using the MD picture; both the contact angle and the adsorbed protein mass density decrease with increasing pH, in line with the increasing magnitude of negative charge on both the protein and the surface. At lower pH < 5 the protein starts to unfold, and the adsorbed mass dramatically decreases. The comprehensive picture that emerges for the formation of oriented protein films with preserved native conformation will help guide efforts to create functional films for new technologies.

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.

Self-Assembled Monolayers of an Azobenzene Derivative on Silica and Their Interactions with Lysozyme

The capability of the photoresponsive isomerization of azobenzene derivatives in self-assembled monolayer (SAM) surfaces to control protein adsorption behavior has very promising applications in antifouling materials and biotechnology. In this study, we performed an atomistic molecular dynamics (MD) simulation in combination with free-energy calculations to study the morphology of azobenzene-terminated SAMs (Azo-SAMs) grafted on a silica substrate and their interactions with lysozyme. Results show that the Azo-SAM surface morphology and the terminal benzene rings' packing are highly correlated with the surface density and the isomer state. Higher surface coverage and the trans-isomer state lead to a more ordered polycrystalline backbone as well as more ordered local packing of benzene rings. On the Azo-SAM surface, water retains a high interfacial diffusivity, whereas the adsorbed lysozyme is found to have extremely low mobility but a relative stable secondary structure. The moderate desorption free energy (∼60 kT) from the trans-Azo-SAM surface was estimated by using both the nonequilibrium-theorem-based Jarzynski's equality and equilibrium umbrella sampling.

Molecular Dynamics Simulations of the Initial Adsorption Stages of Fibrinogen on Mica and Graphite Surfaces

Fibrinogen, a blood glycoprotein of vertebrates, plays an essential role in blood clotting by polymerizing into fibrin when activated. Upon adsorption on material surfaces, it also contributes to determine their biocompatibility and has been implicated in the onset of thrombosis and inflammation at medical implants. Here we present the first fully atomistic simulations of the initial stages of the adsorption process of fibrinogen on mica and graphite surfaces. The simulations reveal a weak adsorption on mica that allows frequent desorption and reorientation events. This adsorption is driven by electrostatic interactions between the protein and the silicate surface as well as the counterion layer. Preferred adsorption orientations for the globular regions of the protein are identified. The adsorption on graphite is found to be stronger with fewer reorientation and desorption events and shows the onset of denaturation of the protein.

Cell Penetrating Peptide Adsorption on Magnetite and Silica Surfaces: A Computational Investigation

Magnetic nanoparticles (MNPs) represent one of the most promising materials as they can act as a versatile platform in the field of bionanotechnology for enhanced imaging, diagnosis, and treatment of various diseases. Silica is the most common compound for preparing coated iron oxide NPs since it improves colloidal stability and the binding affinity for various organic molecules. Biomolecules such as cell penetrating peptides (CPPs) might be employed to decorate MNPs, combining their promising physicochemical properties with a cell penetrating ability. In this work, a computational investigation on adsorption of Antennapedia homeodomain-derived penetrating peptide (pAntp) on silica and magnetite (MAG) surfaces is presented. By employing umbrella sampling molecular dynamics, we provided a quantitative estimation of the pAntp-surface adsorption free energy to highlight the influence of surface hydroxylation state on the adsorption mechanism. The interaction between peptide and surface has shown to be mainly driven by electrostatics. In case of MAG surface, also an important contribution of van der Waals (VdW) attraction was observed. Our data suggest that a competitive mechanism between MNPs and cell membrane might partially inhibit the CPP to carry out its membrane penetrating function.