Translocation of Antimicrobial Peptides across Model Membranes: The Role of Peptide Chain Length.

Cushioned lipid bilayers are structures consisting of a lipid bilayer supported on a solid substrate with an intervening layer of soft material. They offer possibilities for studying the behavior and interactions of biological membranes more accurately under physiological conditions. In this work, we continue our studies of cushion formation induced by histatin 5 (24Hst5), focusing on the effect of the length of the peptide chain. 24Hst5 is a short, positively charged, intrinsically disordered saliva peptide, and here, both a shorter (14Hst5) and a longer (48Hst5) peptide variant were evaluated. Experimental surface active techniques were combined with coarse-grained Monte Carlo simulations to obtain information about these peptides. Results show that at 10 mM NaCl, both the shorter and the longer peptide variants behave like 24Hst5 and a cushion below the bilayer is formed. At 150 mM NaCl, however, no interaction is observed for 24Hst5. On the contrary, a cushion is formed both in the case of 14Hst5 and 48Hst5, and in the latter, an additional thick, diffuse, and highly hydrated layer of peptide and lipid molecules is formed, on top of the bilayer. Similar trends were observed from the simulations, which allowed us to hypothesize that positively charged patches of the amino acids lysine and arginine in all three peptides are essential for them to interact with and translocate over the bilayer. We therefore hypothesize that electrostatic interactions are important for the interaction between the solid-supported lipid bilayers and the peptide depending on the linear charge density through the primary sequence and the positively charged patches in the sequence. The understanding of how, why, and when the cushion is formed opens up the possibility for this system to be used in the research and development of new drugs and pharmaceuticals.

Deeper Insight of the Conformational Ensemble of Intrinsically Disordered Proteins.

It is generally known that, unlike structured proteins, intrinsically disordered proteins, IDPs, exhibit various structures and conformers, the so-called conformational ensemble, CoE. This study aims to better understand the conformers that make up the IDP ensemble by decomposing the CoE into groups separated by their radius of gyration, Rg. A common approach to studying CoE for IDPs is to use low-resolution techniques, such as small-angle scattering, and combine those with computer simulations on different length scales. Herein, the well-studied antimicrobial saliva protein histatin 5 was utilized as a model peptide for an IDP; the average intensity curves were obtained from small-angle X-ray scattering; and compared with fully atomistic, explicit water, molecular dynamics simulations; then, the intensity curve was decomposed with respect to the different Rg values; and their secondary structure propensities were investigated. We foresee that this approach can provide important information on the CoE and the individual conformers within; in that case, it will serve as an additional tool for understanding the IDP structure-function relationship on a more detailed level.

Anion-Specific Adsorption of Carboxymethyl Cellulose on Cellulose.

Integration of fiber modification step with a modern pulp mill is a resource efficient way to produce functional fibers. Motivated by the need to integrate polymer adsorption with the current pulping system, anion-specific effects in carboxymethylcellulose (CMC) adsorption have been studied. The QCM-D adsorption experiments revealed that CMC adsorption to the cellulose model surface is prone to anion-specific effects. A correlation was observed between the adsorbed CMC and the degree of hydration of the co-ions present in the magnesium salts. The presence of a chaotropic co-ion such as nitrate increased the adsorption of CMC on cellulose compared to the presence of the kosmotropic sulfate co-ion. However, anion-specificity was not significant in the case of salts containing zinc cations. The hydration of anions determines the distribution of the ions at the interface. Chaotropic ions, such as nitrates, are likely to be distributed near the chaotropic cellulose surface, causing changes in the ordering of water molecules and resulting in greater entropy gain once released from the surface, thus increasing CMC adsorption.

Interaction of a Histidine-Rich Antimicrobial Saliva Peptide with Model Cell Membranes: The Role of Histidines.

Histatin 5 is a histidine-rich, intrinsically disordered, multifunctional saliva protein known to act as a first line of defense against oral candidiasis caused by Candida albicans. An earlier study showed that, upon interaction with a common model bilayer, a protein cushion spontaneously forms underneath the bilayer. Our hypothesis is that this effect is of electrostatic origin and that the observed behavior is due to proton charge fluctuations of the histidines, promoting attractive electrostatic interactions between the positively charged proteins and the anionic surfaces, with concomitant counterion release. Here we are investigating the role of the histidines in more detail by defining a library of variants of the peptide, where the former have been replaced by the pH-insensitive amino acid glutamine. By using experimental techniques such as circular dichroism, small angle X-ray scattering, quartz crystal microbalance with dissipation monitoring, and neutron reflectometry, it was determined that changing the number of histidines in the peptide sequence did not affect the structure of the peptide dissolved in solution. However, it was shown to affect the penetration depth of the peptide into the bilayer, where all variants except the one with zero histidines were found below the bilayer. A decrease in the number of histidine from the original seven to zero decreases the ability of the peptide to penetrate the bilayer, and the peptide is then also found residing within the bilayer. We hypothesize that this is due to the ability of the histidines to charge titrate, which charges up the peptide, and enables it to penetrate and translocate through the lipid bilayer.

Comparative Performance of Computer Simulation Models of Intrinsically Disordered Proteins at Different Levels of Coarse-Graining.

Coarse-graining is commonly used to decrease the computational cost of simulations. However, coarse-grained models are also considered to have lower transferability, with lower accuracy for systems outside the original scope of parametrization. Here, we benchmark a bead-necklace model and a modified Martini 2 model, both coarse-grained models, for a set of intrinsically disordered proteins, with the different models having different degrees of coarse-graining. The SOP-IDP model has earlier been used for this set of proteins; thus, those results are included in this study to compare how models with different levels of coarse-graining compare. The sometimes naive expectation of the least coarse-grained model performing best does not hold true for the experimental pool of proteins used here. Instead, it showed the least good agreement, indicating that one should not necessarily trust the otherwise intuitive notion of a more advanced model inherently being better in model choice.

PED in 2021: a major update of the protein ensemble database for intrinsically disordered proteins.

The Protein Ensemble Database (PED) (https://proteinensemble.org), which holds structural ensembles of intrinsically disordered proteins (IDPs), has been significantly updated and upgraded since its last release in 2016. The new version, PED 4.0, has been completely redesigned and reimplemented with cutting-edge technology and now holds about six times more data (162 versus 24 entries and 242 versus 60 structural ensembles) and a broader representation of state of the art ensemble generation methods than the previous version. The database has a completely renewed graphical interface with an interactive feature viewer for region-based annotations, and provides a series of descriptors of the qualitative and quantitative properties of the ensembles. High quality of the data is guaranteed by a new submission process, which combines both automatic and manual evaluation steps. A team of biocurators integrate structured metadata describing the ensemble generation methodology, experimental constraints and conditions. A new search engine allows the user to build advanced queries and search all entry fields including cross-references to IDP-related resources such as DisProt, MobiDB, BMRB and SASBDB. We expect that the renewed PED will be useful for researchers interested in the atomic-level understanding of IDP function, and promote the rational, structure-based design of IDP-targeting drugs.

Molecular dynamics simulations and solid-state nuclear magnetic resonance spectroscopy measurements of C-H bond order parameters and effective correlation times in a POPC-GM3 bilayer.

Glycolipids such as gangliosides affect the properties of lipid membranes and in extension the interactions between membranes and other biomolecules like proteins. To better understand how the properties of individual lipid molecules can contribute to shape the functional aspects of a membrane, the spatial restriction and dynamics of C-H bond segments can be measured using nuclear magnetic resonance (NMR) spectroscopy. We combine solid-state NMR spectroscopy with all-atom molecular dynamics (MD) simulations to investigate how ganglioside GM3 affects the bilayer structure and dynamics of C-H bond segments. These two methods yield reorientational correlation functions, molecular profiles of C-H bond order parameters |SCH| and effective correlation times τe, which we compare for lipids in POPC bilayers with and without 30 mol% GM3. Our results revealed that all C-H segments of POPC reorient slower in the presence of GM3 and that the defining features of the GM3-POPC bilayer lie in the GM3 headgroup; it gives the bilayer an extended headgroup layer with high order (|SCH| up to 0.3-0.4) and slow dynamics (τe up to 100 ns), a character that may be mechanistically important in ganglioside interactions with other biomolecules.

Impact of Arginine-Phosphate Interactions on the Reentrant Condensation of Disordered Proteins.

Re-entrant condensation results in the formation of a condensed protein regime between two critical ion concentrations. The process is driven by neutralization and inversion of the protein charge by oppositely charged ions. Re-entrant condensation of cationic proteins by the polyvalent anions, pyrophosphate and tripolyphosphate, has previously been observed, but not for citrate, which has similar charge and size compared to the polyphosphates. Therefore, besides electrostatic interactions, other specific interactions between the polyphosphate ions and proteins must contribute. Here, we show that additional attractive interactions between arginine and tripolyphosphate determine the re-entrant condensation and decondensation boundaries of the cationic, intrinsically disordered saliva protein, histatin 5. Furthermore, we show by small-angle X-ray scattering (SAXS) that polyvalent anions cause compaction of histatin 5, as would be expected based solely on electrostatic interactions. Hence, we conclude that arginine-phosphate-specific interactions not only regulate solution properties but also influence the conformational ensemble of histatin 5, which is shown to vary with the number of arginine residues. Together, the results presented here provide further insight into an organizational mechanism that can be used to tune protein interactions in solution of both naturally occurring and synthetic proteins.

The Effect of Multisite Phosphorylation on the Conformational Properties of Intrinsically Disordered Proteins.

Intrinsically disordered proteins are involved in many biological processes such as signaling, regulation, and recognition. A common strategy to regulate their function is through phosphorylation, as it can induce changes in conformation, dynamics, and interactions with binding partners. Although phosphorylated intrinsically disordered proteins have received increased attention in recent years, a full understanding of the conformational and structural implications of phosphorylation has not yet been achieved. Here, we present all-atom molecular dynamics simulations of five disordered peptides originated from tau, statherin, and ß-casein, in both phosphorylated and non-phosphorylated state, to compare changes in global dimensions and structural elements, in an attempt to gain more insight into the controlling factors. The changes are in qualitative agreement with experimental data, and we observe that the net charge is not enough to predict the impact of phosphorylation on the global dimensions. Instead, the distribution of phosphorylated and positively charged residues throughout the sequence has great impact due to the formation of salt bridges. In statherin, a preference for arginine-phosphoserine interaction over arginine-tyrosine accounts for a global expansion, despite a local contraction of the phosphorylated region, which implies that also non-charged residues can influence the effect of phosphorylation.

Molecular Dynamics Simulations of Phosphorylated Intrinsically Disordered Proteins: A Force Field Comparison.

Phosphorylation is a common post-translational modification among intrinsically disordered proteins and regions, which helps regulate function by changing the protein conformations, dynamics, and interactions with binding partners. To fully comprehend the effects of phosphorylation, computer simulations are a helpful tool, although they are dependent on the accuracy of the force field used. Here, we compared the conformational ensembles produced by Amber ff99SB-ILDN+TIP4P-D and CHARMM36m, for four phosphorylated disordered peptides ranging in length from 14-43 residues. CHARMM36m consistently produced more compact conformations with a higher content of bends, mainly due to more stable salt bridges. Based on comparisons with experimental size estimates for the shortest and longest peptide, CHARMM36m appeared to overestimate the compactness. The difference between the force fields was largest for the peptide showing the greatest separation between positively charged and phosphorylated residues, in line with the importance of charge distribution. For this peptide, the conformational ensemble did not change significantly upon increasing the ionic strength from 0 mM to 150 mM, despite a reduction of the salt-bridging probability in the CHARMM36m simulations, implying that salt concentration has negligible effects in this study.

Spontaneous Formation of Cushioned Model Membranes Promoted by an Intrinsically Disordered Protein.

In this article, it is shown that by exposing commonly used lipids for biomembrane mimicking studies, to a solution containing the histidine-rich intrinsically disordered protein histatin 5, a protein cushion spontaneously forms underneath the bilayer. The underlying mechanism is attributed to have an electrostatic origin, and it is hypothesized that the observed behavior is due to proton charge fluctuations promoting attractive electrostatic interactions between the positively charged proteins and the anionic surfaces, with concomitant counterion release. Hence, we anticipate that this novel "green" approach of forming cushioned bilayers can be an important tool to mimic the cell membrane without the disturbance of the solid substrate, thereby achieving a further understanding of protein-cell interactions.

Self-association of a highly charged arginine-rich cell-penetrating peptide.

Small-angle X-ray scattering (SAXS) measurements reveal a striking difference in intermolecular interactions between two short highly charged peptides-deca-arginine (R10) and deca-lysine (K10). Comparison of SAXS curves at high and low salt concentration shows that R10 self-associates, while interactions between K10 chains are purely repulsive. The self-association of R10 is stronger at lower ionic strengths, indicating that the attraction between R10 molecules has an important electrostatic component. SAXS data are complemented by NMR measurements and potentials of mean force between the peptides, calculated by means of umbrella-sampling molecular dynamics (MD) simulations. All-atom MD simulations elucidate the origin of the R10-R10 attraction by providing structural information on the dimeric state. The last two C-terminal residues of R10 constitute an adhesive patch formed by stacking of the side chains of two arginine residues and by salt bridges formed between the like-charge ion pair and the C-terminal carboxyl groups. A statistical analysis of the Protein Data Bank reveals that this mode of interaction is a common feature in proteins.

Polyelectrolyte-Nanoplatelet Complexation: Is It Possible to Predict the State Diagram?

The addition of polyelectrolytes (PEs) to suspensions of charged colloids, such as nanoplatelets (NPs), is of great interest due to their specific feature of being either a stabilizing or a destabilizing agent. Here, the complexation between a PE and oppositely charged NPs is studied utilizing coarse-grained molecular dynamics simulations based on the continuum model. The complex formation is evaluated with respect to the stoichiometric charge-ratio within the system, as well as by the alternation of the chain properties. It is found that the formed complexes can possess either an extended or a compact shape. Moreover, it is observed that the chain can become overcharged by the oppositely charged NPs. With an increase in chain length, or a decrease in chain flexibility, the complex obtains a more extended shape, where the NPs are less tightly bound to the PE. The latter is also true when reducing the total charge of the chain by varying the linear charge density, whereas in this case, the chain contracts. With our coarse-grained model and molecular dynamics simulations, we are able to predict the composition and the shape of the formed complex and how it is affected by the characteristics of the chain. The take-home message is that the complexation between PEs and NPs results in a versatile and rich state diagram, which indeed is difficult to predict, and dependent on the properties of the chain and the model used. Thus, we propose that the present model can be a useful tool to achieve an understanding of the PE-NPs complexation, a system commonly used in industrial and in technological processes.

Structural Characterization of Histatin 5-Spermidine Conjugates: A Combined Experimental and Theoretical Study.

Histatin 5 (Hst5) is a naturally occurring antimicrobial peptide that acts as the first line of defense against oral candidiasis. It has been shown that conjugation of the active Hst5 fragment, Hst54-15, and the polyamine spermidine (Spd) improves the candidacidal effect. Knowledge about the structure of these conjugates is, however, very limited. Thus, the aim of this study was to characterize the structural properties of the Hst54-15-Spd conjugates by performing atomistic molecular dynamics simulations in combination with small-angle X-ray scattering. It was shown that the Hst54-15-Spd conjugates adopt extended and slightly rigid random coil conformations without any secondary structure in aqueous solution. It is hypothesized that the increased fungal killing potential of Hst54-15-Spd, in comparison with the Spd-Hst54-15 conjugate, is due to the more extended conformations of the former, which cause the bonded Spd molecule to be more accessible for recognition by polyamine transporters in the cell.

Coarse-grained modeling of the intrinsically disordered protein Histatin 5 in solution: Monte Carlo simulations in combination with SAXS.

Monte Carlo simulations and coarse-grained modeling have been used to analyze Histatin 5, an unstructured short cationic salivary peptide known to have anticandidical properties. The calculated scattering functions have been compared with intensity curves and the distance distribution function P(r) obtained from small angle X-ray scattering (SAXS), at both high and low salt concentrations. The aim was to achieve a molecular understanding and a physico-chemical insight of the obtained SAXS results and to gain information of the conformational changes of Histatin 5 due to altering salt content, charge distribution, and net charge. From a modeling perspective, the accuracy of the electrostatic interactions are of special interest. The used coarse-grained model was based on the primitive model in which charged hard spheres differing in charge and in size represent the ionic particles, and the solvent only enters the model through its relative permittivity. The Hamiltonian of the model comprises three different contributions: (i) excluded volumes, (ii) electrostatic, and (iii) van der Waals interactions. Even though the model can be considered as gross omitting all atomistic details, a great correspondence is obtained with the experimental results. Proteins 2016; 84:777-791. © 2016 Wiley Periodicals, Inc.

Confined polyelectrolytes: The complexity of a simple system.

The interaction between polyelectrolytes and counterions in confined situations and the mutual relationship between chain conformation and ion condensation is an important issue in several areas. In the biological field, it assumes particular relevance in the understanding of the packaging of nucleic acids, which is crucial in the design of gene delivery systems. In this work, a simple coarse-grained model is used to assess the cooperativity between conformational change and ion condensation in spherically confined backbones, with capsides permeable to the counterions. It is seen that the variation on the degree of condensation depends on counterion valence. For monovalent counterions, the degree of condensation passes through a minimum before increasing as the confining space diminishes. In contrast, for trivalent ions, the overall tendency is to decrease the degree of condensation as the confinement space also decreases. Most of the particles reside close to the spherical wall, even for systems in which the density is higher closer to the cavity center. This effect is more pronounced, when monovalent counterions are present. Additionally, there are clear variations in the charge along the concentric layers that cannot be totally ascribed to polyelectrolyte behavior, as shown by decoupling the chain into monomers. If both chain and counterions are confined, the formation of a counterion rich region immediately before the wall is observed. Spool and doughnut-like structures are formed for stiff chains, within a nontrivial evolution with increasing confinement.

Role of histidine for charge regulation of unstructured peptides at interfaces and in bulk.

Histidine-rich, unstructured peptides adsorb to charged interfaces such as mineral surfaces and microbial cell membranes. At a molecular level, we investigate the adsorption mechanism as a function of pH, salt, and multivalent ions showing that (1) proton charge fluctuations are-in contrast to the majority of proteins-optimal at neutral pH, promoting electrostatic interactions with anionic surfaces through charge regulation and (2) specific zinc(II)-histidine binding competes with protons and ensures an unusually constant charge distribution over a broad pH interval. In turn, this further enhances surface adsorption. Our analysis is based on atomistic molecular dynamics simulations, coarse grained Metropolis Monte Carlo, and classical polymer density functional theory. This multiscale modeling provides a consistent picture in good agreement with experimental data on Histatin 5, an antimicrobial salivary peptide. Biological function is discussed and we suggest that charge regulation is a significant driving force for the remarkably robust activity of histidine-rich antimicrobial peptides.

Exploring the Role of Globular Domain Locations on an Intrinsically Disordered Region of p53: A Molecular Dynamics Investigation.

The pre-tetramerization loop (PTL) of the human tumor suppressor protein p53 is an intrinsically disordered region (IDR) necessary for the tetramerization process, and its flexibility contributes to the essential conformational changes needed. Although the IDR can be accurately simulated in the traditional manner of molecular dynamics (MD) with the end-to-end distance (EEdist) unhindered, we sought to explore the effects of restraining the EEdist to the values predicted by electron microscopy (EM) and other distances. Simulating the PTL trajectory with a restrained EEdist , we found an increased agreement of nuclear magnetic resonance (NMR) chemical shifts with experiments. Additionally, we observed a plethora of secondary structures and contacts that only appear when the trajectory is restrained. Our findings expand the understanding of the tetramerization of p53 and provide insight into how mutations could make the protein impotent. In particular, our findings demonstrate the importance of restraining the EEdist in studying IDRs and how their conformations change under different conditions. Our results provide a better understanding of the PTL and the conformational dynamics of IDRs in general, which are useful for further studies regarding mutations and their effects on the activity of p53.

Exploring the Functional Landscape of the p53 Regulatory Domain: The Stabilizing Role of Post-Translational Modifications.

This study focuses on the intrinsically disordered regulatory domain of p53 and the impact of post-translational modifications. Through fully atomistic explicit water molecular dynamics simulations, we show the wealth of information and detailed understanding that can be obtained by varying the number of phosphorylated amino acids and implementing a restriction in the conformational entropy of the N-termini of that intrinsically disordered region. The take-home message for the reader is to achieve a detailed understanding of the impact of phosphorylation with respect to (1) the conformational dynamics and flexibility, (2) structural effects, (3) protein interactivity, and (4) energy landscapes and conformational ensembles. Although our model system is the regulatory domain p53 of the tumor suppressor protein p53, this study contributes to understanding the general effects of intrinsically disordered phosphorylated proteins and the impact of phosphorylated groups, more specifically, how minor changes in the primary sequence can affect the properties mentioned above.