Diblock polyampholytes grafted onto spherical particles: Monte Carlo simulation and lattice mean-field theory

Spherical brushes composed of diblock polyampholytes (diblock copolymers with oppositely charged blocks) grafted onto solid spherical particles in aqueous solution are investigated by using the primitive model solved with Monte Carlo simulations and by lattice mean-field theory. Polyampholyte chains of two compositions are considered: a copolymer with a long and a short block, A100B10, and a copolymer with two blocks of equal length, A50B50. The B block is end-grafted onto the surface, and its charge is varied, whereas the charge of the A block is fixed. Single-chain properties, radial and lateral spatial distributions of different types, and structure factors are analyzed. The brush structure strongly depends on the charge of the B block. In the limit of an uncharged B block, the chains are stretched and form an extended polyelectrolyte brush. In the other limit with the charges of the blocks compensating each other, the chains are collapsed and form a polyelectrolyte complex surrounding the particles. At intermediate charge conditions, a polyelectrolyte brush and a polyelectrolyte complex coexist and constitute two substructures of the spherical brush. The differences of the brush structures formed by the A100B10 and A50B50 polyampholytes are also analyzed. Finally, a comparison of the predictions of the two theoretical approaches is made.     

Swelling of cross-linked polyelectrolyte gels

A model of a cross-linked polyelectrolyte gel has been examined using Monte Carlo simulations. The simple model contained a charged defect-free network represented by linked charged beads and explicit counterions. Pressure-density relations for the polyelectrolyte gel, a corresponding non-ionic polymer gel, and several partly or fully degraded gels have been determined. The polyelectrolyte gel displayed a very large swelling capacity, in agreement with experiments. The swelling mechanism and chain properties are discussed and foundations of current theories on gel swelling are examined.     

Stepwise disproportionation in polyelectrolyte complexes

Structural properties and the topology of polyelectrolyte complexes (PECs) formed in solution have been investigated under different conditions by Monte Carlo simulations using a coarse-grained model. The extension of individual polyions has been characterized by their radius of gyration, whereas the composition of the complexes has been investigated by their net charge and their internal topological structure by a novel analysis describing how the shorter polycations link to monomers of the longer polyanion. Conditions have been found at which the polyanion and a given number of polycations form distinguishable complexes differing in (i) the polyanion conformation and (ii) the fraction of polycations being in extended and collapsed states. Thus, at equilibrium, these PECs display a stepwise variation of the degree of intrachain disproportionation within the polyanion (also referred to as intrachain segregation), concomitant with the interchain disproportionation of the polycations, which is in agreement with previous theoretical predictions. The coexistence of the different polyelectrolyte complex structures appears, generally, at mixing ratios close to but different from charge equivalence and, as a consequence, broad polyelectrolyte size distributions are commonly obtained.     

Structure of microemulsion-ABA triblock copolymer networks

Structural equilibrium properties of transient networks formed by microemulsion droplets and ABA triblock copolymers in solution have been studied by Monte Carlo simulation. The droplets were represented by soft spheres, and the polymers were represented by junctions connected by harmonic bonds with an angular potential regulating the intrinsic chain stiffness. The interaction parameters were selected such that the end A-blocks were localized inside the droplets and the middle B-block in the continuous phase. The influence of (i) the polymer concentration, (ii) the polymer stiffness, and (iii) the contour length of the middle B-block on the formation and the structure of the microemulsion-polymer network were investigated using polymer end-to-end separation probability distribution functions, droplet radial distribution functions, droplet-droplet nearest-neighbor probability distribution functions, and network connectivity indicators. An increase of the polymer-droplet number ratio had a strong impact on the network formation. Under typical conditions and at an intermediate polymer-droplet number ratio, (i) the fraction of polymers forming bridges between droplets increased from essentially zero to unity and (ii) the fraction of polymers that were forming loops decreased as the ratio of the polymer end-to-end separation and the surface-to-surface separation between neighboring droplets for a hypothetical homogeneous droplet distribution was increased from 0.5 to 2. For long and flexible polymers, a mesoscopic segregation triggered by a depletion attraction between droplets appeared, and, furthermore, for sufficiently stiff chains, only bridge conformations occurred. The percolation probability could be represented as a function of the average droplet cluster size only, across all systems.     

Molecular design of specific metal-binding peptide sequences from protein fragments: theory and experiment

A novel strategy is presented for designing peptides with specific metal-ion chelation sites, based on linking computationally predicted ion-specific combinations of amino acid side chains coordinated at the vertices of the desired coordination polyhedron into a single polypeptide chain. With this aim, a series of computer programs have been written that 1) creates a structural combinatorial library containing Z(i)-(X)(n)-Z(j) sequences (n=0-14; Z: amino acid that binds the metal through the side chain; X: any amino acid) from the existing protein structures in the non-redundant Protein Data Bank; 2) merges these fragments into a single Z(1)-(X)(n(1) )-Z(2)-(X)(n(2) )-Z(3)-(X)(n(3) )--Z(j) polypeptide chain; and 3) automatically performs two simple molecular mechanics calculations that make it possible to estimate the internal strain in the newly designed peptide. The application of this procedure for the most M(2+)-specific combinations of amino acid side chains (M: metal; see L. Rulísek, Z. Havlas J. Phys. Chem. B 2003, 107, 2376-2385) yielded several peptide sequences (with lengths of 6-20 amino acids) with the potential for specific binding with six metal ions (Co(2+), Ni(2+), Cu(2+), Zn(2+), Cd(2+) and Hg(2+)). The gas-phase association constants of the studied metal ions with these de novo designed peptides were experimentally determined by MALDI mass spectrometry by using 3,4,5-trihydroxyacetophenone as a matrix, whereas the thermodynamic parameters of the metal-ion coordination in the condensed phase were measured by isothermal titration calorimetry (ITC), chelatometry and NMR spectroscopy methods. The data indicate that some of the computationally predicted peptides are potential M(2+)-specific metal-ion chelators.     

Energy-driven asymmetric partitioning of a semiflexible polymer between interconnected cavities

The distribution of a semiflexible chain in the volume of two interconnected spherical cavities of equal size has been investigated by using Monte Carlo simulations. The chain possessed an extension exceeding that of the cavity, leading to large probabilities of translocated states despite the entropic penalty of passing the narrow passage. Furthermore, an asymmetric state with unequal subchain lengths in the two cavities was more favorable than the symmetric state. The preference for the asymmetric state is driven by the bending energy. Basically, in the symmetric state both subchains are forced to be bent, whereas in the asymmetric case only one of the subchains must bend, leading to an overall smaller bending penalty and overall smaller free energy of the asymmetric state. These results are in contrast to the entropy-controlled partitioning of polymers into confinement and the symmetric translocation state appearing for flexible polymers.     

The effect of nanoparticles on amyloid aggregation depends on the protein stability and intrinsic aggregation rate

Nanoparticles interfere with protein amyloid formation. Catalysis of the process may occur due to increased local protein concentration and nucleation on the nanoparticle surface, whereas tight binding or a large particle/protein surface area may lead to inhibition of protein aggregation. Here we show a clear correlation between the intrinsic protein stability and the nanoparticle effect on the aggregation rate. The results were reached for a series of five mutants of single-chain monellin differing in intrinsic stability toward denaturation, for which a correlation between protein stability and aggregation propensity has been previously documented by Szczepankiewicz et al. [Mol. Biosyst.20107 (2), 521-532]. The aggregation process was monitored by thioflavin T fluorescence in the absence and presence of copolymeric nanoparticles with different hydrophobic characters. For mutants with a high intrinsic stability and low intrinsic aggregation rate, we find that amyloid fibril formation is accelerated by nanoparticles. For mutants with a low intrinsic stability and high intrinsic aggregation rate, we find the opposite--a retardation of amyloid fibril formation by nanoparticles. Moreover, both catalytic and inhibitory effects are most pronounced with the least hydrophobic nanoparticles, which have a larger surface accessibility of hydrogen-bonding groups in the polymer backbone.     

Coarse-grained modeling of proline rich protein 1 (PRP-1) in bulk solution and adsorbed to a negatively charged surface

Structural properties of the acidic proline rich protein PRP-1 of salivary origin in bulk solution and adsorbed onto a negatively charged surface have been studied by Monte Carlo simulations. A simple model system with focus on electrostatic interactions and short-ranged attractions among the uncharged amino acids has been used. In addition to PRP-1, some mutants were considered to assess the role of the interactions in the systems. Contrary to polyelectrolytes, the protein has a compact structure in salt-free bulk solutions, whereas at high salt concentration the protein becomes more extended. The protein adsorbs to a negatively charged surface, although its net charge is negative. The adsorbed protein displays an extended structure, which becomes more compact upon addition of salt. Hence, the conformational response upon salt addition in the adsorbed state is the opposite as compared to that in bulk solution. The conformational behavior of PRP-1 in bulk solution and at charged surfaces as well as its propensity to adsorb to surfaces with the same net charge are rationalized by the block polyampholytic character of the protein. The presence of a triad of positively charged amino acids in the C-terminal was found to be important for the adsorption of the protein.