The pH of liquid water is determined by the infrequent process in which water molecules split into short-lived hydroxide and hydronium ions. This reaction is difficult to probe experimentally and challenging to simulate. One of the open questions is whether the local water structure around a slightly stretched OH bond is actually initiating the eventual breakage of this bond or whether this event is driven by a global ordering that involves many water molecules far away from the reaction center. Here, we investigated the self-ionization of water at room temperature by rare-event ab initio molecular dynamics and obtained autoionization rates and activation energies in good agreement with experiments. Based on the analysis of thousands of molecular trajectories, we identified a couple of local order parameters and show that if a bond stretch occurs when all these parameters are around their ideal range, the chance for the first dissociation step (double-proton jump) increases from [Formula: see text] to 0.4. Understanding these initiation triggers might ultimately allow the steering of chemical reactions.

A molecular dynamics modeling and simulation approach is presented and employed to construct porous dextran polymer ion-exchange adsorbent media. Both the activation step of the surface of the pores of the dextran polymer layer grafted on an agarose surface and the immobilization of charged ligands on the activated surface of the porous dextran polymer layer are considered. For the systems studied in this work, the activation step modifies slightly the pore structure of the base, nonactivated porous dextran polymer, while the immobilization of the ligands on the activated pore surface of the dextran layer changes significantly the pore structure of the activated dextran layer. The density distributions of the counterions and immobilized charged ligands along the direction of net transport in the adsorbent media constructed in this study are found to be nonuniform. The variables that affect the shape and magnitude of the density distributions of the counterions and immobilized charged ligands as well as the total number of charged ligands that can be immobilized on the activated porous dextran layer are identified and presented in this work. Furthermore, the data clearly show that there is local nonelectroneutrality in the porous dextran polymer ion-exchange adsorbent media, and this result has very important practical implications for the operation and performance of separation systems involving ion-exchange adsorbent media (e.g., ion-exchange chromatography systems). Also, the results of this work suggest approaches for (1) controlling the immobilization process of charged ligands and (2) constructing and studying the behavior of chromatographic polymeric monoliths and packed bed columns having a gradient of density of functionalities along the axis of the chromatographic polymeric monolith or packed bed column.

The construction and use of nonflat agarose surfaces in a simulation box, together with the employment of criteria for the immobilization of a set of dextran polymer chains on the nonflat agarose surfaces whose mathematical physics is compatible with that of the criteria used for the immobilization of the same set of dextran polymer chains on flat agarose surfaces, are shown to generate, through the use of molecular dynamics simulations whose simulation box has linear dimensions along the lateral directions that are the same when flat and nonflat agarose surfaces are used, dextran porous polymer structures whose pore sizes at the outermost surface and in the vicinity of the outermost surface of the porous medium can be controlled by an indirect manner through the variation of the parameters that characterize the nonflat surface. The use of a nonflat surface for the generation of desired large pores requires only a small or modest increase in the number of solvent molecules in the simulation box, while the use of a flat surface for the construction of the same desired large pores requires significant increases in the size of the linear dimensions of the flat surface. This increases so substantially the number of solvent molecules that the computational loads become intractable. The results in this work show that through the use of nonflat surfaces porous dextran polymer layers having pores of desired sizes can be effectively constructed, and this approach could be used for the design and construction of polymer-based porous adsorbent media that could effectively facilitate the transport and adsorption of an adsorbate biomolecule of interest that must be separated from a mixture of components. A useful definition about the properties that a porous polymer structure must have in order to become, for an adsorbate biomolecule of interest of known molecular size, a useful adsorbent medium, is presented and is used to (1) evaluate the porous polymer structures generated through the employment of different nonflat surface models and (2) determine and select the nonflat surface model from a set of nonflat surface models that is effective in producing promising porous structures. Then a procedure is presented by which a set of porous polymer media is generated through the use of the selected nonflat surface model, and the desired porous structure from this set is determined and could be considered to be used for the transport and immobilization of the selected affinity groups/ligands and the subsequent transport and adsorption of the desired to be separated adsorbate.

Ion-exchange porous adsorbent media having intermediate and low surface densities of dextran polymer grafted on the surface of the throughpores of polymeric monoliths are constructed and characterized by a molecular dynamics modeling and simulation approach that has also been shown to be effective in the construction and characterization of porous ion-exchange adsorbent media whose number of immobilized dextran polymer chains per unit surface area is high. The activation step that prepares the surface of the pores of the dextran polymer layer for the immobilization of the charged ligands insignificantly affected the pore structure of the dextran polymer layer, while this was found to not be the case for previously studied systems that involved high dextran polymer surface densities. Compared to the high dextran polymer density system studied previously, the intermediate dextran polymer density system can generate significantly larger pores but still possesses relatively high interconnection and mutual steric support between dextran chains to exhibit similar structural characteristics and responses to charged ligand immobilization, including dextran layer thickness, stability, monomer distribution, ligand-induced compact chain structures, dextran layer shrinkage, distributions of ligands and counterions, and local nonelectroneutrality. The low dextran polymer density system having relatively isolated dextran chains and insufficient mutual steric support can result in even larger pores than those obtained in the intermediate dextran polymer density system, but a significantly thinner porous dextran polymer layer and different dextran monomer distributions are obtained in the low dextran polymer density system. More importantly, the gradient of the local nonelectroneutrality after the immobilization of the charged ligands is significantly smaller in magnitude in the low dextran polymer density system than that obtained in the system having intermediate dextran polymer density, and, despite a lack of porous layer depth to accommodate adsorbate biomolecules in large amounts, it could still be useful for the effective transport and adsorption of very large biomolecules. Compared with the polymeric monoliths without a porous dextran polymer layer grafted on the surface of their throughpores, the intermediate and low dextran polymer density systems explored and studied in this work provide pore structures with desirable characteristics for the effective transport of adsorbate biomolecules and substantially larger effective surface areas and throughput capacities for the adsorption of the adsorbate biomolecules.

Adaptive immune receptor repertoires (AIRR) are key targets for biomedical research as they record past and ongoing adaptive immune responses. The capacity of machine learning (ML) to identify complex discriminative sequence patterns renders it an ideal approach for AIRR-based diagnostic and therapeutic discovery. To date, widespread adoption of AIRR ML has been inhibited by a lack of reproducibility, transparency, and interoperability. immuneML ( immuneml.uio.no ) addresses these concerns by implementing each step of the AIRR ML process in an extensible, open-source software ecosystem that is based on fully specified and shareable workflows. To facilitate widespread user adoption, immuneML is available as a command-line tool and through an intuitive Galaxy web interface, and extensive documentation of workflows is provided. We demonstrate the broad applicability of immuneML by (i) reproducing a large-scale study on immune state prediction, (ii) developing, integrating, and applying a novel method for antigen specificity prediction, and (iii) showcasing streamlined interpretability-focused benchmarking of AIRR ML. 1.