Counting charges on membrane-bound peptides

Disentangling the molecular polarizability and first hyperpolarizability of methanol-air interfaces

Liquid-air interfaces have extensive implications in different areas of interest because the dynamical processes at the interface can be different from those in bulk. Thus, its characterization, understanding, and control may be pivotal in advancing discoveries. However, characterizing the interface requires special and selective tools to avoid signals from the bulk region. This surface specificity and versatility is achieved by using the second harmonic generation (SHG) responses. This study adopts multiscale simulation methods to evaluate the surface SHG responses of methanol-air interfaces with submonolayer resolution tackled by sequentially using classical molecular dynamics simulations under different temperatures and then employing quantum chemistry methods to compute the molecular first hyperpolarizabilities (β). This approach ensures the configurational diversity required to evaluate the average β values. The main achievements are (i) a quasi-absence of surface sensitivity of the mean polarizability 〈α〉 with values about 2% larger than those obtained in bulk, (ii) conversely, smooth variations on the polarizability anisotropy Δα are observed up to the fourth molecular layer at around 20 Å from the interface, and (iii) narrow interfacial effects on the SHG responses, β(-2ω;ω,ω), which are limited to the first molecular layer (∼3.0 Å) and characterized by a high contrast in the βZZZ(-2ω;ω,ω) tensor component between the first and the subsequent layers. Similar trends are obtained at different temperatures or when increasing the number of methanol molecules treated at the quantum chemistry level, indicating the robustness of the approach for describing the dipolar molecular responses of air-liquid interfaces.

Why local and non-local terms are essential for second harmonic generation simulation?

Second Harmonic Generation (SHG) today represents one of the most powerful techniques to selectively probe all types of interfaces. However, the origin of the SHG signal at a molecular level is still debated since the local dipole contribution, which is strongly correlated to the molecular orientation can be counterbalanced by non-local quadrupole contributions. Here, we propose a method to simulate the SHG signal of a model water/air interface from the molecular response of each contribution. This method includes both local and non-local terms, which are represented, respectively, by the dependency of the polarisability and hyperpolarisability upon the chemical environment of the molecule and by the bulk quadrupole response. The importance of both terms for the sound simulation of the SHG signals and their interpretation is assessed. We demonstrate that the sole dipole term is unable to simulate a SHG signal, even if the dependency of the hyperpolarisability on the local environment is considered. The inclusion of the bulk quadrupole contribution, which largely dominates the dipole contribution, is essential to predict the SHG response, although the accuracy of the prediction is increased when the dependency upon the local environment is considered.

Deciphering second harmonic generation signals

Second harmonic generation (SHG) has emerged as one of the most powerful techniques used to selectively monitor surface dynamics and reactions for all types of interfaces as well as for imaging non-centrosymmetric structures, although the molecular origin of the SHG signal is still poorly understood. Here, we present a breakthrough approach to predict and interpret the SHG signal at the atomic level, which is freed from the hyperpolarisability concept and self-consistently considers the non-locality and the coupling with the environment. The direct ab initio method developed here shows that a bulk quadrupole contribution significantly overwhelms the interface dipole term in the purely interfacial induced second-order polarisation for water/air interfaces. The obtained simulated SHG responses are in unprecedented agreement with the experimental signal. This work not only paves the road for the prediction of SHG response from more complex interfaces of all types, but also suggests new insights in the interpretation of the SHG signal at a molecular level. In particular, it highlights the modest influence of the molecular orientation and the high significance of the bulk quadrupole contribution, which does not depend on the interface, in the total experimental response.

Second harmonic generation is one of the most powerful techniques used to selectively probe interfaces of all types. The direct ab initio method developed here allows predicting the signal and highlights the importance of local and non-local effects.

First hyperpolarizability of water at the air-vapor interface: a QM/MM study questions standard experimental approximations

Surface Second-Harmonic Generation (S-SHG) experiments provide a unique approach to probe interfaces. One important issue for S-SHG is how to interpret the S-SHG intensities at the molecular level. Established frameworks commonly assume that each molecule emits light according to an average molecular hyperpolarizability tensor β(-2ω,ω,ω). However, for water molecules, this first hyperpolarizability is known to be extremely sensitive to their environment. We have investigated the molecular first hyperpolarizability of water molecules within the liquid-vapor interface, using a quantum description with explicit, inhomogeneous electrostatic embedding. The resulting average molecular first hyperpolarizability tensor depends on the distance relative to the interface, and it practically respects the Kleinman symmetry everywhere in the liquid. Within this numerical approach, based on the dipolar approximation, the water layer contributing to the Surface Second Harmonic Generation (S-SHG) intensity is less than a nanometer. The results reported here question standard interpretations based on a single, averaged hyperpolarizability for all molecules at the interface. Not only the molecular first hyperpolarizability tensor significantly depends on the distance relative to the interface, but it is also correlated to the molecular orientation. Such hyperpolarizability fluctuations may impact the S-SHG intensity emitted by an aqueous interface.

Heterogeneity in Lateral Distribution of Polycations at the Surface of Lipid Membrane: From the Experimental Data to the Theoretical Model

Natural and synthetic polycations of different kinds attract substantial attention due to an increasing number of their applications in the biomedical industry and in pharmacology. The key characteristic determining the effectiveness of the majority of these applications is the number of macromolecules adsorbed on the surface of biological cells or their lipid models. Their study is complicated by a possible heterogeneity of polymer layer adsorbed on the membrane. Experimental methods reflecting the structure of the layer include the electrokinetic measurements in liposome suspension and the boundary potential of planar bilayer lipid membranes (BLM) and lipid monolayers with a mixed composition of lipids and the ionic media. In the review, we systematically analyze the methods of experimental registration and theoretical description of the laterally heterogeneous structures in the polymer layer published in the literature and in our previous studies. In particular, we consider a model based on classical theory of the electrical double layer, used to analyze the available data of the electrokinetic measurements in liposome suspension with polylysines of varying molecular mass. This model suggests a few parameters related to the heterogeneity of the polymer layer and allows determining the conditions for its appearance at the membrane surface. A further development of this theoretical approach is discussed.

Water at charged interfaces

The ubiquity of aqueous solutions in contact with charged surfaces and the realization that the molecular-level details of water-surface interactions often determine interfacial functions and properties relevant in many natural processes have led to intensive research. Even so, many open questions remain regarding the molecular picture of the interfacial organization and preferential alignment of water molecules, as well as the structure of water molecules and ion distributions at different charged interfaces. While water, solutes and charge are present in each of these systems, the substrate can range from living tissues to metals. This diversity in substrates has led to different communities considering each of these types of aqueous interface. In this Review, by considering water in contact with metals, oxides and biomembranes, we show the essential similarity of these disparate systems. While in each case the classical mean-field theories can explain many macroscopic and mesoscopic observations, it soon becomes apparent that such theories fail to explain phenomena for which molecular properties are relevant, such as interfacial chemical conversion. We highlight the current knowledge and limitations in our understanding and end with a view towards future opportunities in the field.

Inhomogeneity of polylysine adsorption layers on lipid membranes revealed by theoretical analysis of electrokinetic data and molecular dynamics simulations

The adsorption of large polycations on a charged lipid membrane is qualitatively different from the small inorganic cations, which almost uniformly populate the membrane surface. We assume that the polycationic adsorption layer might be laterally inhomogeneous starting from a certain polymer length, and this effect can be more visible for membranes with low anionic lipid content. To study systems with inhomogeneous adsorption layers, we carried out electrokinetic measurements of mobility of liposomes containing anionic and neutral phospholipids in the presence of polylysine molecules. Some of these systems were simulated by all-atom molecular dynamics. Here we proposed a theoretical approach accounting for the formation of separated regions at the membrane surface, which differ in charge density and surface potential. Our model allowed us to determine the adsorption layer's geometric parameters such as surface coverage and surface-bound monomer fraction of polymer, which correlate with the molecular dynamics (MD) simulations. We demonstrated that the configuration polylysine adopts on the membrane surface (tall or planar) depends on the polymer/membrane charge ratio. Both theory and MD indicate a decrease in the anionic lipid content, alongside with a decrease in the bound monomer fraction and corresponding increase in the extension length of the adsorbed polymers.

Influence of Sensor Coating and Topography on Protein and Nanoparticle Interaction with Supported Lipid Bilayers

Supported lipid bilayers (SLBs) have proven to be valuable model systems for studying the interactions of proteins, peptides, and nanoparticles with biological membranes. The physicochemical properties (e.g., topography, coating) of the solid substrate may affect the formation and properties of supported phospholipid bilayers, and thus, subsequent interactions with biomolecules or nanoparticles. Here, we examine the influence of support coating (SiO₂ vs Si₃N₄) and topography [sensors with embedded vs protruding gold nanodisks for nanoplasmonic sensing (NPS)] on the formation and subsequent interactions of supported phospholipid bilayers with the model protein cytochrome c and with cationic polymer-wrapped quantum dots using quartz crystal microbalance with dissipation monitoring and NPS techniques. The specific protein and nanoparticle were chosen because they differ in the degree to which they penetrate the bilayer. We find that bilayer formation and subsequent non-penetrative association with cytochrome c were not significantly influenced by substrate composition or topography. In contrast, the interactions of nanoparticles with SLBs depended on the substrate composition. The substrate-dependence of nanoparticle adsorption is attributed to the more negative zeta-potential of the bilayers supported by the silica vs the silicon nitride substrate and to the penetration of the cationic polymer wrapping the nanoparticles into the bilayer. Our results indicate that the degree to which nanoscale analytes interact with SLBs may be influenced by the underlying substrate material.

Interfacial water and ion distribution determine ζ potential and binding affinity of nanoparticles to biomolecules

The molecular features that dictate interactions between functionalized nanoparticles and biomolecules are not well understood. This is in part because for highly charged nanoparticles in solution, establishing a clear connection between the molecular features of surface ligands and common experimental observables such as ζ potential requires going beyond the classical models based on continuum and mean field models. Motivated by these considerations, molecular dynamics simulations are used to probe the electrostatic properties of functionalized gold nanoparticles and their interaction with a charged peptide in salt solutions. Counterions are observed to screen the bare ligand charge to a significant degree even at the moderate salt concentration of 50 mM. As a result, the apparent charge density and ζ potential are largely insensitive to the bare ligand charge densities, which fall in the range of ligand densities typically measured experimentally for gold nanoparticles. While this screening effect was predicted by classical models such as the Manning condensation theory, the magnitudes of the apparent surface charge from microscopic simulations and mean-field models are significantly different. Moreover, our simulations found that the chemical features of the surface ligand (e.g., primary vs. quaternary amines, heterogeneous ligand lengths) modulate the interfacial ion and water distributions and therefore the interfacial potential. The importance of interfacial water is further highlighted by the observation that introducing a fraction of hydrophobic ligands enhances the strength of electrostatic binding of the charged peptide. Finally, the simulations highlight that the electric double layer is perturbed upon binding interactions. As a result, it is the bare charge density rather than the apparent charge density or ζ potential that better correlates with binding affinity of the nanoparticle to a charged peptide. Overall, our study highlights the importance of molecular features of the nanoparticle/water interface and underscores a set of design rules for the modulation of electrostatic driven interactions at nano/bio interfaces.

Combined tetraphenylethylene fluorogens with positive charge for imaging capsule-covered pathogens

Capsule-covered pathogens can cause serious infectious diseases, and are highly pathogenic to humans. Herein, we developed four positively charged tetraphenylethylene derivatives (PC-TPEgens) that in certain combinations were applied to identify capsule-bearing pathogens using fluorescence imaging. The dual-charged probes were used to visualize the entire process of phagocytosis of pathogens into macrophages.

Electrostatics, Hydrogen Bonding, and Molecular Structure at Polycation and Peptide:Lipid Membrane Interfaces

Polycation and peptide-modified surfaces represent opportunities for developing potentially novel biocidal materials in a growing effort to combat bacterial resistance to traditional bactericides. It is well-known that the positive charge of these compounds is crucial to their function in biofouling prevention and as antimicrobials; however, methods for quantifying the number of positive charges on surface-bound polycations and peptides are necessary to predict, control, and optimize the design and therefore the utility of these compounds. This Spotlight on Applications reports on such an approach that combines second harmonic generation (SHG) spectroscopy, quartz crystal microbalance with dissipation monitoring (QCM-D), and atomistic simulations to obtain mechanistic insight into polycation-membrane interactions using supported lipid bilayers (SLBs) as our model system. We find that at high surface coverage, the large polycations we surveyed feature a considerably smaller percentage of ionization when compared to the smaller polycations and peptides. At these high charge densities, we suspect a pK_a shift of the charged groups to lower charge-charge repulsion as well as the formation of a looplike conformation such that less monomeric units form contact-ion pairs with the bilayer. Our sum frequency generation (SFG) spectroscopy results complement our understanding of the polycation-membrane interaction. At a high density of the polycation poly(allylamine hydrochloride) (PAH), second-order spectral line shapes are consistent with the expulsion of interfacial water molecules possibly due to contact-ion pair formation between PAH and the lipid bilayer. This finding could be essential for understanding the underlying first steps of cell lysis and penetration by polycations and should be explored further.

All-Factor Analysis and Correlations on the Transmembrane Process for Arginine-Rich Cell-Penetrating Peptides

Currently, arginine-rich cell-penetrating peptides (CPPs), due to their little cytotoxicity and high transmembrane efficiency, are considered as one of the important intracellular carriers. Although the mechanism of the transmembrane process for arginine-rich CPPs was proposed, the quantitative correlations and the key factors involved in this process still deserve further investigation. In this study, all-atom molecular dynamics and the umbrella sampling technique were employed to study the arginine-rich CPPs transmembrane process. In the adsorption process of CPPs from solution to the surface of the lipid bilayer, the adsorption free energy (ΔG_A) is found to be linearly related to the interaction energy change (ΔE_A): ΔG_A = 0.0426ΔE_A + 36.7, R² = 0.92. In the CPPs transmembrane process, the transmembrane free energy barrier (ΔG_B) is roughly correlated with the corresponding interaction energy change (ΔE_B): ΔG_B = 0.108ΔE_B +135, R² = 0.73. The multiple salt bridges of guanidinium-PO₄ account for 65% of the overall interaction energy, so the increased negative charges of the lipid bilayer or more salt bridges would facilitate CPPs adsorption and transmembrane processes. Also, the increased negative charges of the lipid bilayer would reduce the amount of water to be carried into the pore and further reduce the ΔG_B. The peptide backbone would not have a direct impact on transmembrane efficiency. The ΔG_B is also found to be related to the length of the pore (L): ΔG_B = 46.2L - 31.3, R² = 0.92, which makes the transmembrane efficiency estimable. This work is expected to deliver an in-depth understanding and help the optimization of CPPs.

Rational design of novel amphipathic antimicrobial peptides focused on the distribution of cationic amino acid residues

Antimicrobial peptides (AMPs) have garnered much attention as novel therapeutic agents against infectious diseases. They exhibit antimicrobial activity through microbial membrane disruption based on their amphipathic properties. In this study, we rationally designed and synthesized a series of novel AMPs Block, Stripe, and Random, and revealed that Stripe exhibits potent antimicrobial activity against Gram-positive and Gram-negative microbes. Moreover, we also demonstrated that Stripe disrupts both Gram-positive and Gram-negative mimetic bacterial membranes. Finally, we investigated the hemolytic activity and cytotoxicity in human blood cells and human cell lines, and found that Stripe exhibited neither. These data indicated that Stripe is a promising antimicrobial reagent that does not display significant cytotoxicity.

Electrostatics and Interactions of an Ionizable Silica Nanoparticle Approaching a Plasma Membrane

The surface charge of the plasma membrane (PM) and the large salt content of the extracellular space ensure a significant role of the electrostatic effect dictating the interaction between the PM and an approaching nanoparticle (NP). In this article, we theoretically study the case of an ionizable silica NP approaching the PM. We witness that the surface charge of the silica NP, dictated by the surface ionization of the silica in the electrostatic environment created by the PM surface charge and the extracellular ion concentration, decreases as it approaches the PM. In other words, a silica NP is more negative away from the PM than in close proximity to the PM. Accordingly, we witness a significantly lower repulsion between the PM and NP favoring the approach and the interactions of the silica NP with the PM. Additionally, the presence of the silica NP in the vicinity of the PM induces a large nonisopotentiality, even across a fully permeable PM. We anticipate that these findings will be critically important in the better design of the widely used silica NPs for targeted drug and gene deliveries.

Hydrogen-Bond Networks near Supported Lipid Bilayers from Vibrational Sum Frequency Generation Experiments and Atomistic Simulations

We report vibrational sum frequency generation (SFG) spectra in which the C-H stretches of lipid alkyl tails in fully hydrogenated single- and dual-component supported lipid bilayers are detected along with the O-H stretching continuum above the bilayer. As the salt concentration is increased from ∼10 μM to 0.1 M, the SFG intensities in the O-H stretching region decrease by a factor of 2, consistent with significant absorptive-dispersive mixing between χ⁽²⁾ and χ⁽³⁾ contributions to the SFG signal generation process from charged interfaces. A method for estimating the surface potential from the second-order spectral lineshapes (in the OH stretching region) is presented and discussed in the context of choosing truly zero-potential reference states. Aided by atomistic simulations, we find that the strength and orientation distribution of the hydrogen bonds over the purely zwitterionic bilayers are largely invariant between submicromolar and hundreds of millimolar concentrations. However, specific interactions between water molecules and lipid headgroups are observed upon replacing phosphocholine (PC) lipids with negatively charged phosphoglycerol (PG) lipids, which coincides with SFG signal intensity reductions in the 3100-3200 cm^-1 frequency region. The atomistic simulations show that this outcome is consistent with a small, albeit statistically significant, decrease in the number of water molecules adjacent to both the lipid phosphate and choline moieties per unit area, supporting the SFG observations. Ultimately, the ability to probe hydrogen-bond networks over lipid bilayers holds the promise of opening paths for understanding, controlling, and predicting specific and nonspecific interactions between membranes and ions, small molecules, peptides, polycations, proteins, and coated and uncoated nanomaterials.