Ion-Specific Effects on Ion and Polyelectrolyte Solvation.

Ion-specific effects on aqueous solvation of monovalent counter ions, Na + ${^+ }$ , K + ${^+ }$ , Cl - ${^- }$ , and Br - ${^- }$ , and two model polyelectrolytes (PEs), poly(styrene sulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA) were here studied with ab initio molecular dynamics (AIMD) and classical molecular dynamics (MD) simulations based on the OPLS-aa force-field which is an empirical fixed point-charge force-field. Ion-specific binding to the PE charge groups was also characterized. Both computational methods predict similar response for the solvation of the PEs but differ notably in description of ion solvation. Notably, AIMD captures the experimentally observed differences in Cl - ${^- }$ and Br - ${^- }$ anion solvation and binding with the PEs, while the classical MD simulations fail to differentiate the ion species response. Furthermore, the findings show that combining AIMD with the computationally less costly classical MD simulations allows benefiting from both the increased accuracy and statistics reach.

Exploration of Chemical Diversity in Intercellular Quorum Sensing Signalling Systems in Prokaryotes.

Quorum sensing (QS) serves as a vital means of intercellular signalling in a variety of prokaryotes, which enables single cells to act in multicellular configurations. The potential to control community-wide responses has also sparked numerous recent biotechnological innovations. However, our capacity to utilize intercellular communication is hindered due to a scarcity of complementary signalling systems and a restricted comprehension of interconnections between these systems caused by variations in their dynamic range. In this study, we utilize uniform manifold approximation and projection and extended-connectivity fingerprints to explore the available chemical space of QS signalling molecules. We investigate and experimentally characterize a set of closely related QS signalling ligands, consisting of N-acyl homoserine lactones and the aryl homoserine lactone p-coumaroyl, as well as a set of more widely diverging QS ligands, consisting of photopyrones, dialkylresorcinols, 3,5-dimethylpyrazin-2-ol and autoinducer-2, and define their performance. We report on a set of six signal- and promoter-orthogonal intercellular QS signalling systems, significantly expanding the toolkit for engineering community-wide behaviour. Furthermore, we demonstrate that ligand diversity can serve as a statistically significant tool to predict much more complicated ligand-receptor interactions. This approach highlights the potential of dimensionality reduction to explore chemical diversity in microbial dynamics.

A Perspective on the Glass Transition and the Dynamics of Polyelectrolyte Multilayers and Complexes.

Polyelectrolyte multilayers (PEMs) or polyelectrolyte complexes (PECs), formed by layer-by-layer assembly or the mixing of oppositely charged polyelectrolytes (PEs) in aqueous solution, respectively, have potential applications in health, energy, and the environment. PEMs and PECs are very tunable because their structure and properties are influenced by factors such as pH, ionic strength, salt type, humidity, and temperature. Therefore, it is increasingly important to understand how these factors affect PECs and PEMs on a molecular level. In this Feature Article, we summarize our contributions to the field in the development of approaches to quantify the swelling, thermal properties, and dynamic mechanical properties of PEMs and PECs. First, the role of water as a plasticizer and in the glass-transition temperature (Tg) in both strong poly(diallyldimethylammonium)/poly(sodium 4-styrenesulfonate) (PDADMA/PSS) and weak poly(allylamine hydrochloride)/poly(acrylic acid) (PAH/PAA) systems is presented. Then, factors influencing the dynamics of PECs and PEMs are discussed. We also reflect on the swelling of PEMs in response to different salts and solvent additives. Last, the nature of water's microenvironment in PEMs/PECs is discussed. A special emphasis is placed on experimental techniques, along with molecular simulations. Taken together, this review presents an outlook and offers recommendations for future research directions, such as studying the additional effects of hydrogen-bonding hydrophobic interactions.

pH dependence of the assembly mechanism and properties of poly(L-lysine) and poly(L-glutamic acid) complexes.

We show by extensive experimental characterization combined with molecular simulations that pH has a major impact on the assembly mechanism and properties of poly(L-lysine) (PLL) and poly(L-glutamic acid) (PGA) complexes. A combination of dynamic light scattering (DLS) and laser Doppler velocimetry (LDV) is used to assess the complexation, charge state, and other physical characteristics of the complexes, isothermal titration calorimetry (ITC) is used to examine the complexation thermodynamics, and circular dichroism (CD) is used to extract the polypeptides' secondary structure. For enhanced analysis and interpretation of the data, analytical ultracentrifugation (AUC) is used to define the precise molecular weights and solution association of the peptides. Molecular dynamics simulations reveal the associated intra- and intermolecular binding changes in terms of intrinsic vs. extrinsic charge compensation, the role of hydrogen bonding, and secondary structure changes, aiding in the interpretation of the experimental data. We combine the data to reveal the pH dependency of PLL/PGA complexation and the associated molecular level mechanisms. This work shows that not only pH provides a means to control complex formation but also that the associated changes in the secondary structure and binding conformation can be systematically used to control materials assembly. This gives access to rational design of peptide materials via pH control.

SARS-CoV-2 Spike Protein (RBD) Subunit Adsorption at Abiotic Surfaces and Corona Formation at Polymer Particles.

The adsorption kinetics of the SARS-CoV-2 spike protein subunit with the receptor binding domain at abiotic surfaces was investigated. A combination of sensitive methods was used such as atomic force microscopy yielding a molecular resolution, a quartz microbalance, and optical waveguide lightmode spectroscopy. The two latter methods yielded in situ information about the protein adsorption kinetics under flow conditions. It was established that at pH 3.5-4 the protein adsorbed on mica and silica surfaces in the form of compact quasi-spherical aggregates with an average size of 14 nm. The maximum coverage of the layers was equal to 3 and 1 mg m-2 at pH 4 and 7.4, respectively. The experimental data were successfully interpreted in terms of theoretical results derived from modeling. The experiments performed for flat substrates were complemented by investigations of the protein corona formation at polymer particles carried out using in situ laser Doppler velocimetry technique. In this way, the zeta potential of the protein layers was acquired as a function of the coverage. Applying the electrokinetic model, these primary data were converted to the dependence of the subunit zeta potential on pH. It was shown that a complete acid-base characteristic of the layer can be acquired only using nanomolar quantities of the protein.

SARS-CoV-2 virion physicochemical characteristics pertinent to abiotic substrate attachment.

The structure, size, and main physicochemical characteristics of the SARS-CoV-2 virion with the spike transmembrane protein corona were discussed. Using these data, diffusion coefficients of the virion in aqueous media and in air were calculated. The structure and dimensions of the spike protein derived from molecular dynamic modeling and thorough cryo-electron microscopy measurements were also analyzed. The charge distribution over the molecule was calculated and shown to be largely heterogeneous. Although the stalk part is negatively charged, the top part of the spike molecule, especially the receptor binding domain, remains positively charged for a broad range of pH. It is underlined that such a charge distribution promotes the spike corona stability and enhances the virion attachment to receptors and surfaces, mostly negatively charged. The review is completed by the analysis of experimental data pertinent to the spike protein adsorption at abiotic surfaces comprising nanoparticle carrier particles. It is argued that these theoretical and experimental data can be used for developing quantitative models of virus attachment to surfaces, facilitating adequate analysis of future experimental results.

Self-Assembly of Silk-like Protein into Nanoscale Bicontinuous Networks under Phase-Separation Conditions.

Liquid-liquid phase separation of biomacromolecules is crucial in various inter- and extracellular biological functions. This includes formation of condensates to control, e.g., biochemical reactions and structural assembly. The same phenomenon is also found to be critically important in protein-based high-performance biological materials. Here, we use a well-characterized model triblock protein system to demonstrate the molecular level formation mechanism and structure of its condensate. Large-scale molecular modeling supported by analytical ultracentrifuge characterization combined with our earlier high magnification precision cryo-SEM microscopy imaging leads to deducing that the condensate has a bicontinuous network structure. The bicontinuous network rises from the proteins having a combination of sites with stronger mutual attraction and multiple weakly attractive regions connected by flexible, multiconfigurational linker regions. These attractive sites and regions behave as stickers of varying adhesion strength. For the examined model triblock protein construct, the ß-sheet-rich end units are the stronger stickers, while additional weaker stickers, contributing to the condensation affinity, rise from spring-like connections in the flexible middle region of the protein. The combination of stronger and weaker sticker-like connections and the flexible regions between the stickers result in a versatile, liquid-like, self-healing structure. This structure also explains the high flexibility, easy deformability, and diffusion of the proteins, decreasing only 10-100 times in the bicontinuous network formed in the condensate phase in comparison to dilute protein solution. The here demonstrated structure and condensation mechanism of a model triblock protein construct via a combination of the stronger binding regions and the weaker, flexible sacrificial-bond-like network as well as its generalizability via polymer sticker models provide means to not only understand intracellular organization, regulation, and cellular function but also to identify direct control factors for and to enable engineering improved protein and polymer constructs to enhance control of advanced fiber materials, smart liquid biointerfaces, or self-healing matrices for pharmaceutics or bioengineering materials.

Fourier transform infrared spectroscopy investigation of water microenvironments in polyelectrolyte multilayers at varying temperatures.

Polyelectrolyte multilayers (PEMs) are thin films formed by the alternating deposition of oppositely charged polyelectrolytes. Water plays an important role in influencing the physical properties of PEMs, as it can act both as a plasticizer and swelling agent. However, the way in which water molecules distribute around and hydrate ion pairs has not been fully quantified with respect to both temperature and ionic strength. Here, we examine the effects of temperature and ionic strength on the hydration microenvironments of fully immersed poly(diallyldimethylammonium)/polystyrene sulfonate (PDADMA/PSS) PEMs. This is accomplished by tracking the OD stretch peak using attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectroscopy at 0.25-1.5 M NaCl and 35-70 °C. The OD stretch peak is deconvoluted into three peaks: (1) high frequency water, which represents a tightly bound microenvironment, (2) low frequency water, which represents a loosely bound microenvironment, and (3) bulk water. In general, the majority of water absorbed into the PEM exists in a bound state, with little-to-no bulk water observed. Increasing temperature slightly reduces the amount of absorbed water, while addition of salt increases the amount of absorbed water. Finally, a van't Hoff analysis is applied to estimate the enthalpy (11-22 kJ mol-1) and entropy (48-79 kJ mol-1 K-1) of water exchanging from low to high frequency states.

Myoglobin molecule charging in electrolyte solutions.

All atom molecular dynamic modeling was applied in order to determine water molecule and electrolyte ion concentration profiles around and inside the myoglobin molecule at various pH values. Significant penetration of counter ions into the molecule was confirmed. The electric potential distribution within and outside the molecule was quantitatively described using the non-linear Poisson-Boltzmann (PB) approach. Using this model, calculations were performed, yielding the surface and zeta potential for various physicochemical parameters, comprising pH, the electric permittivity, the ion penetration depth and the protein volume fraction (crowding effect). The theoretical results were used for the interpretation of experimental data acquired under different ionic strengths and temperatures by electrophoretic mobility measurements. It is confirmed that the experimental data are adequately reflected for acidic pH values by the non-linear PB model where the nominal molecule charge was calculated from the H++ model. The deviations occurring for larger pH values were accounted for by considering additional non-electrostatic interactions stemming from the van der Waals and ion-induced dipole forces. In this way, it is both experimentally and theoretically confirmed that the effective charge of the myoglobin molecule in electrolyte solutions is considerably smaller than the nominal, structure-based, predicted charge. As a result, under physiological conditions prevailing, e.g. in skeletal muscles, the effective charge of the myoglobin molecule should practically vanish. One can expect that the approach developed in this work can be applied for predicting charging mechanisms of other protein molecules characterized by an analogous charge vs. pH characteristic, e.g., the SARS-CoV-2 virus spike proteins, and for soft particles with pH responsive characteristics.

Comparing water-mediated hydrogen-bonding in different polyelectrolyte complexes.

All-atom molecular dynamics simulations are used to investigate the polyelectrolyte-specific influence of hydration and temperature on water diffusion in hydrated polyelectrolyte complexes (PECs). Two model PECs were compared: poly(allylamine hydrochloride) (PAH)-poly(sodium 4-styrenesulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA)-poly(acrylic acid) (PAA). The findings show that the strength of the hydrogen bonding i.e. polyelectrolyte water interaction has enormous influence on the water mobility, which has implications for PEC structure and properties. A 10-fold difference in the average water diffusion coefficient between PAH-PSS and PDADMA-PAA PECs at the same hydration level is observed. The vast majority of the water molecules hydrating the PDADMA-PAA PECs, for hydrations in the range of 26-38 wt%, are effectively immobilized, whereas for PAH-PSS PECs the amount of immobilized water decreases with hydration. This points to the polyelectrolyte-specific character of the PE-water hydrogen bonding relationship with temperature. PAA-water hydrogen bonds are found to be significantly less sensitive to temperature than for PSS-water. The polyelectrolyte-water interactions, investigated via radial distribution function, hydrogen bond distance and angle distributions, are connected with resulting structure of the PECs. The PDADMA-PAA and PAH-PSS PECs are prepared experimentally and the states of water at different hydration levels is determined using differential scanning calorimetry (DSC). Experiments confirm the differences between PDADMA-PAA and PAH-PSS PECs observed in the theoretical modelling. The results suggest that the initial predictions of the PEC's bonding with water can be based on simple molecular-level considerations.

Particulate Coatings via Evaporation-Induced Self-Assembly of Polydisperse Colloidal Lignin on Solid Interfaces.

Polydisperse smooth and spherical biocolloidal particles were suspended in aqueous media and allowed to consolidate via evaporation-induced self-assembly. The stratification of the particles at the solid-air interface was markedly influenced, but not monotonically, by the drying rate. Cross-sectional imaging via electron microscopy indicated a structured coating morphology that was distinctive from that obtained by using particles with a mono- or bimodal distribution. Segregation patterns were found to derive from the interplay of particle diffusion, interparticle forces, and settling dynamics. Supporting our experimental findings, computer simulations showed an optimal drying rate for achieving maximum segregation. Overall, stratified coatings comprising nano- and microparticles derived from lignin are expected to open opportunities for multifunctional structures that can be designed and predicted on the basis of experimental Péclet numbers and computational order.

QCM-D Investigation of Swelling Behavior of Layer-by-Layer Thin Films upon Exposure to Monovalent Ions.

Polyelectrolyte multilayers and layer-by-layer assemblies are susceptible to structural changes in response to ionic environment. By altering the salt type and ionic strength, structural changes can be induced by disruption of intrinsically bound ion pairs within the multilayer network via electrostatic screening. Notably, high salt concentrations have been used for the purposes of salt-annealing and self-healing of LbL assemblies with KBr, in particular, yielding a remarkably rapid response. However, to date, the structural and swelling effects of various monovalent ion species on the behavior of LbL assemblies remain unclear, including a quantitative view of ion content in the LbL assembly and thickness changes over a wide concentration window. Here, we investigate the effects of various concentrations of KBr (0 to 1.6 M) on the swelling and de-swelling of LbL assemblies formed from poly(diallyldimethylammonium) polycation (PDADMA) and poly(styrene sulfonate) polyanion (PSS) in 0.5 M NaCl using quartz-crystal microbalance with dissipation (QCM-D) monitoring as compared to KCl, NaBr, and NaCl. The ion content after salt exchange is quantified using neutron activation analysis (NAA). Our results demonstrate that Br- ions have a much greater effect on the structure of as-prepared thin films than Cl- at ionic strengths above assembly conditions, which is possibly caused by the more chaotropic nature of Br-. It is also found that the anion in general dominates the swelling response as compared to the cation because of the excess PDADMA in the multilayer. Four response regimes are identified that delineate swelling due to electrostatic repulsion, slight contraction, swelling due to doping, and film destruction as ionic strength increases. This understanding is critical if such materials are to be used in applications requiring submersion in chemically dynamic environments such as sensors, coatings on biomedical implants, and filtration membranes.

Ability of the Poisson-Boltzmann equation to capture molecular dynamics predicted ion distribution around polyelectrolytes.

Here, we examine polyelectrolyte (PE) and ion chemistry specificity in ion condensation via all-atom molecular dynamics (MD) simulations and assess the ability of the Poisson-Boltzmann (PB) equation to describe the ion distribution predicted by the MD simulations. The PB model enables the extraction of parameters characterizing ion condensation. We find that the modified PB equation which contains the effective PE radius and the energy of the ion-specific interaction as empirical fitting parameters describes ion distribution accurately at large distances but close to the PE, especially when strongly localized charge or specific ion binding sites are present, the mean field description of PB fails. However, the PB model captures the MD predicted ion condensation in terms of the Manning radius and fraction of condensed counterions for all the examined PEs and ion species. We show that the condensed ion layer thickness in our MD simulations collapses on a single master curve for all the examined simple, monovalent ions (Na+, Br+, Cs+, Cl-, and Br-) and PEs when plotted against the Manning parameter (and consequently the PE line charge density). The significance of this finding is that, contrary to the Manning radius extracted from the mean field PB model, the condensed layer thickness in the all atom detail MD modelling does not depend on the PE chemistry or counterion type. Furthermore, the fraction of condensed counterions in the MD simulations exceeds the PB theory prediction. The findings contribute toward understanding and modelling ion distribution around PEs and other charged macromolecules in aqueous solutions, such as DNA, functionalized nanotubes, and viruses.

Counter Anion Type Influences the Glass Transition Temperature of Polyelectrolyte Complexes.

Salt acts as a plasticizer in polyelectrolyte complexes (PECs), which impacts the physical, thermal, and mechanical properties, thus having implications in applications, such as drug delivery, energy storage, and smart coatings. Added salt disrupts polycation-polyanion intrinsic ion pairs, lowering a hydrated PEC's glass transition temperature (Tg). However, the relative influence of counterion type on the PEC's Tg is not well understood. Here, the effect of anion type (NaCl, NaBr, NaNO3, and NaI) on the Tg of solid-like, hydrated PECs composed of poly(diallydimethylammonium) (PDADMA)-poly(styrenesulfonate) (PSS) is investigated. With increasing the chaotropic nature of the salt anion, the Tg decreases. The relative differences are attributed to the doping level, the amount of bound water, the mobility of water molecules within the PECs, and the strength of interactions between the PEs. For all studied salt concentrations and salt types, the Tg followed the scaling of -1/Tg ≈ ln([IP]/[H2O]), in which [IP]/[H2O] is the ratio of intrinsic pairs to water. The scaling estimates that about 7 to 17% of the intrinsic ion pairs should be weakened for the PEC to partake in a glass transition. Put together, this study highlights that the Tg in PECs is impacted by the salt anion, but the mechanism of the glass transition remains unchanged.

Solid-Liquid-Solution Phases in Poly(diallyldimethylammonium)/Poly(acrylic acid) Polyelectrolyte Complexes at Varying Temperatures.

The coacervation and complexation of oppositely charged polyelectrolytes are dependent on numerous environmental and preparatory factors, but temperature is often overlooked. Temperature effects remain unclear because the temperature dependence of both the dielectric constant and polymer-solvent interaction parameter can yield lower and/or upper critical solution phase behaviors for PECs. Further, secondary interactions, such as hydrogen bonding, can affect the temperature response of a PEC. That is, mixtures of oppositely charged polyelectrolytes can exhibit phase separation upon lowering and/or increasing the mixture's temperature. Here, the phase behavior of poly(diallylmethylammonium)/poly(acrylic acid) (PDADMA/PAA) complexes under varying KBr ionic strengths, mixing ratios, and temperatures at a fixed pH (in which PAA hydrogen bonding can occur) is examined. At room temperature, the PDADMA/PAA PECs exhibit four different phase states: precipitate, coexisting precipitate and coacervate, solid-like gel, and coacervate. Variable-temperature optical microscopy reveals the upper critical solution temperature (UCST) at which each phase transitioned to a solution state. Interestingly, the UCST value is highly dependent on the original phase of the PEC, in which solid-like precipitates exhibit higher UCST values. Large-scale all-atom molecular dynamics (MD) simulations support that precipitates exhibit kinetic trapping, which may contribute to the higher UCST values observed in the experiment. Taken together, this study highlights the significance of temperature on the phase behavior of PECs, which may play a larger role in stimuli-responsive materials, membraneless organelles, and separations applications.

Accelerated Engineering of ELP-Based Materials through Hybrid Biomimetic-De Novo Predictive Molecular Design.

Efforts to engineer high-performance protein-based materials inspired by nature have mostly focused on altering naturally occurring sequences to confer the desired functionalities, whereas de novo design lags significantly behind and calls for unconventional innovative approaches. Here, using partially disordered elastin-like polypeptides (ELPs) as initial building blocks this work shows that de novo engineering of protein materials can be accelerated through hybrid biomimetic design, which this work achieves by integrating computational modeling, deep neural network, and recombinant DNA technology. This generalizable approach involves incorporating a series of de novo-designed sequences with α-helical conformation and genetically encoding them into biologically inspired intrinsically disordered repeating motifs. The new ELP variants maintain structural conformation and showed tunable supramolecular self-assembly out of thermal equilibrium with phase behavior in vitro. This work illustrates the effective translation of the predicted molecular designs in structural and functional materials. The proposed methodology can be applied to a broad range of partially disordered biomacromolecules and potentially pave the way toward the discovery of novel structural proteins.

Mechanisms of Fibroblast Growth Factor 21 Adsorption on Macroion Layers: Molecular Dynamics Modeling and Kinetic Measurements.

Molecular dynamic modeling and various experimental techniques, including multi-angle dynamic light scattering (MADLS), streaming potential, optical waveguide light spectroscopy (OWLS), quartz crystal microbalance with dissipation (QCM), and atomic force microscopy (AFM), were applied to determine the basic physicochemical parameters of fibroblast growth factor 21 in electrolyte solutions. The protein size and shape, cross-section area, dependence of the nominal charge on pH, and isoelectric point of 5.3 were acquired. These data enabled the interpretation of the adsorption kinetics of FGF 21 on bare and macrocation-covered silica investigated by OWLS and QCM. It was confirmed that the protein molecules irreversibly adsorbed on the latter substrate, forming layers with controlled coverage up to 0.8 mg m-2, while their adsorption on bare silica was much smaller. The viability of two cell lines, CHO-K1 and L-929, on both bare and macrocation/FGF 21-covered substrates was also determined. It is postulated that the acquired results can serve as useful reference systems for designing complexes that can extend the half-life of FGF 21 in its active state.

Nanoscale insights into the local structural rearrangements of amyloid-ß induced by bexarotene.

A better understanding of the abnormal protein aggregation and the effect of anti-aggregation agents on the fibrillation pathways and the secondary structure of aggregates can determine strategies for the early treatment of dementia. Herein, we present a combination of experimental and theoretical studies providing new insights into the influence of the anti-aggregation drug bexarotene on the secondary structure of individual amyloid-ß aggregates and its primary aggregation. The molecular rearrangements and the spatial distribution of ß-sheets within individual aggregates were monitored at the nanoscale with infrared nanospectroscopy. We observed that bexarotene limits the parallel ß-sheets formation, known to be highly abundant in fibrils at later phases of the amyloid-ß aggregation composed of in-register cross-ß structure. Moreover, we applied molecular dynamics to provide molecular-level insights into the investigated system. Both theoretical and experimental results revealed that bexarotene slows down the protein aggregation process via steric effects, largely prohibiting the antiparallel to parallel ß-sheet rearrangement. We also found that bexarotene interacts not only via the single hydrogen bond formation with the peptide backbone but also with the amino acid side residue via a hydrophobic effect. The studied model of the drug-amyloid-ß interaction contributes to a better understanding of the inhibition mechanism of the amyloid-ß aggregation by the small molecule drugs. However, our nanoscale findings need to meet in vivo research requiring different analytical approaches.

Modeling of LbL multilayers with controlled thickness, roughness, and specific surface area.

We present computer simulation results of the layer by layer self-assembling process of colloidal particles. We have generated five multilayer structures of monodisperse spherical particles according to a generalized model of random sequential adsorption of hard spheres. The multilayers, each created at a different single-layer surface coverage, are of similar thickness. We have compared the transparency of the five multilayers and the structure of their outer layers in terms of the two-dimensional pair-correlation function. We have analyzed the variation of multilayer thickness with the number of adsorbed layers. We have also calculated the root-mean-square roughness of the multilayers as a function of the number of adsorption cycles. Finally, we have determined the specific surface area of the porous films as a function of the distance from the solid substrate. Our results suggest that in the limit of low porosity the multilayer transparency decreases exponentially with its porosity. The multilayer thickness is directly proportional to the number of adsorption cycles. The average single-layer thickness grows asymptotically with the single-layer coverage. We have also found that with the number of adsorbed layers the multilayer roughness increases to an asymptotic value. We have observed oscillatory variations of the multilayer specific surface area, decaying exponentially with the distance from the substrate. The decay length of the oscillation increases exponentially with the surface coverage. We have also determined the particle layer interpenetration for each multilayer and we have found that it decreases exponentially with the increase of the coverage. Our results suggest that all the film characteristics strongly depend on the method of its preparation and can be controlled by manipulating the single-layer surface coverage or deposition time. The results can be useful for efficient designing multilayers with desired properties.

Random sequential adsorption: An efficient tool for investigating the deposition of macromolecules and colloidal particles.

Random Sequential Adsorption (RSA) is one of the most efficient theoretical models used to investigate adsorption of macromolecules and particles, with a long-standing tradition in the field of colloid and interface science. In the first part of this paper, we demonstrate how the RSA model can be applied to interpret the experimental data and extract information about the density of the adsorption monolayer, the kinetics of its growth, and microstructural properties such as pair-correlation function and monolayer roughness. We briefly summarized the most important generalizations of the RSA model for monolayers and reviewed its extensions considering, e.g., various particle shapes, the introduction of electrostatic interaction, or adsorption on non-uniform substrates. We thoroughly scrutinized the extended RSA model developed for bilayer and multilayer formation. We collected the mean saturated packing fractions of various two- and three-dimensional objects and provided the most accurate result for two-dimensional disk packing. In the second part of this paper, we summarize various numerical algorithms and techniques that allow one to effectively implement RSA algorithms. We describe efficient methods for detecting intersections of various shapes and techniques enabling generation of strictly saturated RSA packings built of a wide range of different shapes. We hinted at how an inherently sequential RSA scheme can be parallelized. Finally, we critically discuss the limitations of the model and possible directions for future studies.