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.

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.

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.

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.

From dilute to concentrated solutions of intrinsically disordered proteins: Interpretation and analysis of collected data.

Intrinsically disordered proteins (IDPs) have a broad energy landscape and consequently sample many different conformations in solution. The innate flexibility of IDPs is exploited in their biological function, and in many instances allows a single IDP to regulate a range of processes in vivo. Due to their highly flexible nature, characterizing the structural properties of IDPs is not straightforward. Often solution-based methods such as Nuclear Magnetic Resonance (NMR), Förster Resonance Energy Transfer (FRET), and Small-Angle X-ray Scattering (SAXS) are used. SAXS is indeed a powerful technique to study the structural and conformational properties of IDPs in solution, and from the obtained SAXS spectra, information about the average size, shape, and extent of oligomerization can be determined. In this chapter, we will introduce model-free methods that can be used to interpret SAXS data and introduce methods that can be used to interpret SAXS data beyond analytical models, for example, by using atomistic and different levels of coarse-grained models in combination with molecular dynamics (MD) and Monte Carlo simulations.

From dilute to concentrated solutions of intrinsically disordered proteins: Sample preparation and data collection.

It is well-known that an increasing proportion of proteins, protein regions, and partners of globular proteins are being recognized as having an intrinsic disorder, and therefore, not adopting a single three-dimensional structure in solution. For these proteins, small-angle X-ray scattering (SAXS) has become a premier method for examination, since it can provide information about the ensemble of the structural conformations as well as the intermolecular interactions. SAXS measurements can be performed from low to high protein concentrations under different physicochemical properties of the solution. The focus of this chapter is to introduce the basics of how to use SAXS for protein samples, for new and less experienced users, in a simple and concise manner, with emphasis on highly flexible proteins and regions. Methodological aspects in the sample preparation, experiment design, and data collection stages are raised that should be considered prior to attempting SAXS experiments. This is to ensure that high-quality SAXS data is obtained that enables accurate analysis. However, many of the points raised will also be worth considering for SAXS experiments of globular proteins.

Molecular dynamics simulations of the adsorption of an intrinsically disordered protein: Force field and water model evaluation in comparison with experiments.

This study investigates possible structural changes of an intrinsically disordered protein (IDP) when it adsorbs to a solid surface. Experiments on IDPs primarily result in ensemble averages due to their high dynamics. Therefore, molecular dynamics (MD) simulations are crucial for obtaining more detailed information on the atomistic and molecular levels. An evaluation of seven different force field and water model combinations have been applied: (A) CHARMM36IDPSFF + CHARMM-modified TIP3P, (B) CHARMM36IDPSFF + TIP4P-D, (C) CHARMM36m + CHARMM-modified TIP3P, (D) AMBER99SB-ILDN + TIP3P, (E) AMBER99SB-ILDN + TIP4P-D, (F) AMBERff03ws + TIP4P/2005, and (G) AMBER99SB-disp + disp-water. The results have been qualitatively compared with those of small-angle X-ray scattering, synchrotron radiation circular dichroism spectroscopy, and attenuated total reflectance Fourier transform infrared spectroscopy. The model IDP corresponds to the first 33 amino acids of the N-terminal of the magnesium transporter A (MgtA) and is denoted as KEIF. With a net charge of +3, KEIF is found to adsorb to the anionic synthetic clay mineral Laponite® due to the increase in entropy from the concomitant release of counterions from the surface. The experimental results show that the peptide is largely disordered with a random coil conformation, whereas the helical content (α- and/or 310-helices) increased upon adsorption. MD simulations corroborate these findings and further reveal an increase in polyproline II helices and an extension of the peptide conformation in the adsorbed state. In addition, the simulations provided atomistic resolution of the adsorbed ensemble of structures, where the arginine residues had a high propensity to form hydrogen bonds with the surface. Simulations B, E, and G showed significantly better agreement with experiments than the other simulations. Particularly noteworthy is the discovery that B and E with TIP4P-D water had superior performance to their corresponding simulations A and D with TIP3P-type water. Thus, this study shows the importance of the water model when simulating IDPs and has also provided an insight into the structural changes of surface-active IDPs induced by adsorption, which may play an important role in their function.

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.

Clustering and cross-linking of the wheat storage protein α-gliadin: A combined experimental and theoretical approach.

Our aim was to understand mechanisms for clustering and cross-linking of gliadins, a wheat seed storage protein type, monomeric in native state, but incorporated in network while processed. The mechanisms were studied utilizing spectroscopy and high-performance liquid chromatography on a gliadin-rich fraction, in vitro produced α-gliadins, and synthetic gliadin peptides, and by coarse-grained modelling, Monte Carlo simulations and prediction algorithms. In solution, gliadins with α-helix structures (dip at 205 nm in CD) were primarily present as monomeric molecules and clusters of gliadins (peaks at 650- and 700-s on SE-HPLC). At drying, large polymers (Rg 90.3 nm by DLS) were formed and ß-sheets increased (14% by FTIR). Trained algorithms predicted aggregation areas at amino acids 115-140, 150-179, and 250-268, and induction of liquid-liquid phase separation at P- and Poly-Q-sequences (Score = 1). Simulations showed that gliadins formed polymers by tail-to-tail or a hydrophobic core (Kratky plots and Ree = 35 and 60 for C- and N-terminal). Thus, the N-terminal formed clusters while the C-terminal formed aggregates by disulphide and lanthionine bonds, with favoured hydrophobic clustering of similar/exact peptide sections (synthetic peptide mixtures on SE-HPLC). Mechanisms of clustering and cross-linking of the gliadins presented here, contribute ability to tailor processing results, using these proteins.

Self-Diffusive Properties of the Intrinsically Disordered Protein Histatin 5 and the Impact of Crowding Thereon: A Combined Neutron Spectroscopy and Molecular Dynamics Simulation Study.

Intrinsically disordered proteins (IDPs) are proteins that, in comparison with globular/structured proteins, lack a distinct tertiary structure. Here, we use the model IDP, Histatin 5, for studying its dynamical properties under self-crowding conditions with quasi-elastic neutron scattering in combination with full atomistic molecular dynamics (MD) simulations. The aim is to determine the effects of crowding on the center-of-mass diffusion as well as the internal diffusive behavior. The diffusion was found to decrease significantly, which we hypothesize can be attributed to some degree of aggregation at higher protein concentrations, (≥100 mg/mL), as indicated by recent small-angle X-ray scattering studies. Temperature effects are also considered and found to, largely, follow Stokes-Einstein behavior. Simple geometric considerations fail to accurately predict the rates of diffusion, while simulations show semiquantitative agreement with experiments, dependent on assumptions of the ratio between translational and rotational diffusion. A scaling law that previously was found to successfully describe the behavior of globular proteins was found to be inadequate for the IDP, Histatin 5. Analysis of the MD simulations show that the width of the distribution with respect to diffusion is not a simplistic mirroring of the distribution of radius of gyration, hence, displaying the particular features of IDPs that need to be accounted for.

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.

Force Field Effects in Simulations of Flexible Peptides with Varying Polyproline II Propensity.

Five peptides previously suggested to possess polyproline II (PPII) structure have here been investigated by using atomistic molecular dynamics simulations to compare how well four different force fields known for simulating intrinsically disordered proteins relatively well (Amber ff99SB-disp, Amber ff99SB-ILDN, CHARM36IDPSFF, and CHARMM36m) can capture this secondary structure element. The results revealed that all force fields sample PPII structures but to different extents and with different propensities toward other secondary structure elements, in particular, the ß-sheet and "random coils". A cluster analysis of the simulations of histatin 5 also revealed that the conformational ensembles of the force fields are quite different. We compared the simulations to circular dichroism and nuclear magnetic resonance spectroscopy experiments and conclude that further experiments and methods for interpreting them are needed to assess the accuracy of force fields in determining PPII structure.

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.

Building polymer-like clusters from colloidal particles with isotropic interactions, in aqueous solution.

HYPOTHESIS: Colloidal particles that interact via a long-ranged repulsive barrier in combination with a very short-ranged attractive minimum can "polymerize" to form highly anisotropic structures. Motivated by previous experimental achievements in non-aqueous solvents, and recent theoretical predictions, we hypothesize that it is possible to construct clusters that resemble linear or branched polymers, in aqueous solution. If these clusters are not too large, they may even remain dispersed, but even if they grow large enough to sediment, they may be collected and used in future applications. EXPERIMENTS: In this work, we specifically synthesize poly (ethylene glycol) (PEG) chains, grafted onto poly (styrene) (PS) particles in aqueous solution, and adjust the conditions so that strongly anisotropic and isolated polymer-like clusters are formed. These conditions include a very low ionic strength (the particles are weakly charged), a relatively high temperature, and a low particle concentration. An important criterion is that the particle size is large enough to admit structural analyses via confocal laser scanning microscopy (CLSM). We have furthermore utilized Metropolis Monte Carlo (MC) simulation to generate theoretical predictions of these cluster formations. We have conducted such simulations of 3D as well as 2D systems, where the latter is also relevant, given that the clusters sometimes deposit onto the glass surfaces upon imaging. A simplistic particle-particle potential of mean force is adopted for the simulations, but we also invoke a more elaborate theoretical model, to demonstrate that similar interactions can be obtained when the grafted chains are treated explicitly. FINDINGS: According to our Zeta potential measurements, the particles indeed carry a weak negative charge, presumably due to ion specific adsorption. Furthermore, by ensuring that the ionic strength is very low, with a Debye length similar to the particle size, we could use temperature to control the hydrophobicity of the grafted PEG layer, and thus the strength of the short-ranged attraction. We were indeed able to establish highly anisotropic structures, that resemble linear or branched polymers, which we could image by CLSM. The average degree of polymerization could be adjusted by a variation of the particle concentration.

From Attraction to Repulsion to Attraction: Non-monotonic Temperature Dependence of Polymer-Mediated Interactions in Colloidal Dispersions.

In this work, we have synthesized polystyrene particles that carry short end-grafted polyethylene glycol (PEG) chains. We then added dissolved 100 kDa PEG polymers and monitored potential flocculation by confocal microscopy. Qualitative predictions, based on previous theoretical developments in this field (Xie, F.; et al. Soft Matter 2016, 12, 658), suggest a non-monotonic temperature response. These theories propose that the "free" (dissolved) polymers will mediate attractive depletion interactions at room temperature, with a concomitant clustering/flocculation at a sufficiently high polymer concentration. At high temperatures, where the solvent is poorer, this is predicted to be replaced by attractive bridging interactions, again resulting in particle condensation. Interestingly enough, our theoretical framework, based on classical density functional theory, predicts an intermediate temperature regime where the polymer-mediated interactions are repulsive! This obviously implies a homogeneous dispersion in this regime. These qualitative predictions have been experimentally tested and confirmed in this work, where flocs of particles start to form at room temperature for a high enough polymer dosage. At temperatures near 45 °C, the flocs redisperse, and we obtain a homogeneous sample. However, samples at about 75 °C will again display clusters and eventually phase separation. Using results from these studies, we have been able to fine-tune parameters of our coarse-grained theoretical model, resulting in predictions of temperature-dependent stability that display semiquantitative accuracy. A crucial aspect is that under "intermediate" conditions, where the polymers neither adsorb nor desorb at the particle surfaces, the polymer-mediated equilibrium interaction is repulsive.

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.

The Effects of Chain Length on the Structural Properties of Intrinsically Disordered Proteins in Concentrated Solutions.

Intrinsically disordered proteins (IDP) are proteins that sample a heterogeneous ensemble of conformers in solution. An estimated 25-30% of all eukaryotic proteins belong to this class. In vivo, IDPs function under conditions that are highly crowded by other biological macromolecules. Previous research has highlighted that the presence of crowding agents can influence the conformational ensemble sampled by IDPs, resulting in either compaction or expansion. The effects of self-crowding of the disordered protein Histatin 5 has, in an earlier study, been found to have limited influence on the conformational ensemble. In this study, it is examined whether the short chain length of Histatin 5 can explain the limited effects of crowding observed, by introducing (Histatin 5)2, a tandem repeat of Histatin 5. By utilizing small-angle X-ray scattering, it is shown that the conformational ensemble is conserved at high protein concentrations, in resemblance with Histatin 5, although with a lowered protein concentration at which aggregation arises. Under dilute conditions, atomistic molecular dynamics and coarse-grained Monte Carlo simulations, as well as an established scaling law, predicted more extended conformations than indicated by experimental data, hence implying that (Histatin 5)2 does not behave as a self-avoiding random walk.

The extracellular juncture domains in the intimin passenger adopt a constitutively extended conformation inducing restraints to its sphere of action.

Enterohemorrhagic and enteropathogenic Escherichia coli are among the most important food-borne pathogens, posing a global health threat. The virulence factor intimin is essential for the attachment of pathogenic E. coli to the intestinal host cell. Intimin consists of four extracellular bacterial immunoglobulin-like (Big) domains, D00-D2, extending into the fifth lectin subdomain (D3) that binds to the Tir-receptor on the host cell. Here, we present the crystal structures of the elusive D00-D0 domains at 1.5 Å and D0-D1 at 1.8 Å resolution, which confirms that the passenger of intimin has five distinct domains. We describe that D00-D0 exhibits a higher degree of rigidity and D00 likely functions as a juncture domain at the outer membrane-extracellular medium interface. We conclude that D00 is a unique Big domain with a specific topology likely found in a broad range of other inverse autotransporters. The accumulated data allows us to model the complete passenger of intimin and propose functionality to the Big domains, D00-D0-D1, extending directly from the membrane.