Experimental and computational study of membrane affinity for selected energetic compounds

Bioconcentration potential of ionic liquids: New data on membrane partitioning and its comparison with predictions obtained by COSMOmic

Ionic liquids (ILs) have recently gained significant attention in both the scientific community and industry, but there is a limited understanding of the potential risks they might pose to the environment and human health, including their potential to accumulate in organisms. While membrane and storage lipids have been considered as primary sorption phases driving bioaccumulation, in this study we used an in vitro tool known as solid-supported lipid membranes (SSLMs) to investigate the affinity of ILs to membrane lipid - phosphatidylcholine and compare the results with an existing in silico model. Our findings indicate that ILs may have a strong affinity for the lipids that form cell membranes, with the key factor being the length of the cation's side chain. For quaternary ammonium cations, increase in membrane affinity (logMA) was observed from 3.45 ± 0.06 at 10 carbon atoms in chain to 4.79 ± 0.06 at 14 carbon atoms. We also found that the anion can significantly affect the membrane partitioning of the cation, even though the anions themselves tend to have weaker interactions with phospholipids than the cations of ILs. For 1-methyl-3-octylimidazolium cation the presence of tricyanomethanide anion caused increase in logMA to 4.23 ± 0.06. Although some of our data proved to be consistent with predictions made by the COSMOmic model, there are also significant discrepancies. These results suggest that further research is needed to improve our understanding of the mechanisms and structure-activity relationships involved in ILs bioconcentration and to develop more accurate predictive models.

Insight into the binding model of per- and polyfluoroalkyl substances to proteins and membranes

Legacy per- and polyfluoroalkyl substances (PFASs) have elicited much concern because of their ubiquitous distribution in the environment and the potential hazards they pose to wildlife and human health. Although an increasing number of effective PFAS alternatives are available in the market, these alternatives bring new challenges. This paper comprehensively reviews how PFASs bind to transport proteins (e.g., serum albumin, liver fatty acid transport proteins and organic acid transporters), nuclear receptors (e.g., peroxisome proliferator activated receptors, thyroid hormone receptors and reproductive hormone receptors) and membranes (e.g., cell membrane and mitochondrial membrane). Briefly, the hydrophobic fluorinated carbon chains of PFASs occupy the binding cavities of the target proteins, and the acid groups of PFASs form hydrogen bonds with amino acid residues. Various structural features of PFAS alternatives such as chlorine atom substitution, oxygen atom insertion and a branched structure, introduce variations in their chain length and hydrophobicity, which potentially change the affinity of PFAS alternatives for endogenous proteins. The toxic effects and mechanisms of action of legacy PFASs can be demonstrated and compared with their alternatives using binding models. In future studies, in vitro experiments and in silico quantitative structure-activity relationship modeling should be better integrated to allow more reliable toxicity predictions for both legacy and alternative PFASs.

The effects of molecular weight and orientation on the membrane permeation and partitioning of polycyclic aromatic hydrocarbons: a computational study

Membrane permeation and the partitioning of polycyclic aromatic hydrocarbons (PAHs) are crucial aspects affecting their carcinogenicity and mutagenicity. However, a clear understanding of these processes is still rare due to the difficulty of determining the details experimentally. Here, the interactions between PAHs and lipid bilayers were studied by molecular simulations, mainly to check the influence of molecular weight and orientation. The liposome-water partition coefficient (K_LW), transmembrane time (τ), and permeability coefficient (P) of the PAHs were calculated by integrating free energy profiles from umbrella sampling. For selected PAHs, the membrane adsorption is a spontaneous process. The preferred location is near the CC bond and the orientation is related to the molecular structure. The P values of all the PAHs are basically the same order of magnitude, which means that the molecular weight contributes little to the process. As for K_LW and τ, they show obvious increases with different molecular weights. Unconstrained simulations showed that a flat orientation on the membrane surface would prevent PAHs from being transported through the membrane. Highly hydrophobic driving forces are not always good for the absorption of PAHs, especially the formation of aggregates. In addition, the orientations and energetic barriers of PAHs near the midplane of the lipid bilayer explain the different transitions of high- and low-weight PAHs. This work provides molecular level details relating to the interactions of PAHs with lipid membranes, with significance for understanding the health effects of PAHs.

Bioconcentration potential and microbial toxicity of onium cations in photoacid generators

Despite the widespread utilization of onium salts as photoacid generators (PAGs) in semiconductor photolithography, their environmental, health, and safety (EHS) properties remain poorly understood. The present work reports the bioconcentration potential of five representative onium species (four sulfonium and one iodonium compound) by determining the octanol-water partition coefficient (P_OW) and lipid membrane affinity coefficient (K_MA); microbial toxicity was evaluated using the bioluminescent bacterium Aliivibrio fischeri (Microtox bioassay). Four of the oniums exhibited varying degrees of hydrophobic (lipophilic) partitioning (log P_OW: 0.08-4.12; K_MA: 1.70-5.62). A strong positive linear correlation was observed between log P_OW and K_MA (K_MA = log P_OW + 1.76, R² = 0.99). The bioconcentration factors (log BCF) estimated from P_OW and K_MA for the four oniums ranged from 0.13 to 3.67 L kg^-1. Bis-(4-tert-butyl phenyl)-iodonium and triphenylsulfonium had 50% inhibitory concentrations (IC₅₀) of 4.8 and 84.6 μM, whereas the IC₅₀ values of the other three oniums were not determined because these values were higher than their aqueous solubility. Given the increased regulatory scrutiny regarding the use and potential health impacts from onium PAGs, this study fulfills critical knowledge gaps concerning the EHS properties of PAG oniums, enabling more comprehensive evaluation of their environmental impacts and potential risk management strategies.

In vitro methods for predicting the bioconcentration of xenobiotics in aquatic organisms

The accumulation of anthropogenic chemical substances in aquatic organisms is an immensely important issue from the point of view of environmental protection. In the context of the increasing number and variety of compounds that may potentially enter the environment, there is a need for efficient and reliable solutions to assess the risks. However, the classic approach of testing with fish or other animals is not sufficient. Due to very high costs, significant time and labour intensity, as well as ethical concerns, in vivo methods need to be replaced by new laboratory-based tools. So far, many models have been developed to estimate the bioconcentration potential of chemicals. However, most of them are not sufficiently reliable and their predictions are based on limited input data, often obtained with doubtful quality. The octanol-water partition coefficient is still often used as the main laboratory tool for estimating bioconcentration. However, according to current knowledge, this method can lead to very unreliable results, both for neutral species and, above all, for ionic compounds. It is therefore essential to start using new, more advanced and credible solutions on a large scale. Over the last years, many in vitro methods have been newly developed or improved, allowing for a much more adequate estimation of the bioconcentration potential. Therefore, the aim of this work was to review the most recent laboratory methods for assessing the bioconcentration potential and to evaluate their applicability in further research.

A relationship between membrane permeation and partitioning of nitroaromatic explosives and their functional groups. A computational study

Nitroaromatic explosives, such as 2,4,6-trinitrotoluene, are representative aromatic compounds, which are generally highly toxic.

Collective absorption of 2,4,6-trinitrotoluene into lipid membranes and its effects on bilayer properties. A computational study

The widely used explosive, 2,4,6-trinitrotoluene (TNT), is a highly toxic chemical, which can cause hepatitis, cataracts, jaundice and so on, in humans. The interaction between TNT and biological membranes is crucial for understanding its toxic effects. Here, we mainly focused on molecular-level mechanisms for the collective adsorption of TNT into lipid membranes and the corresponding effects on bilayer properties by all-atom molecular dynamics simulations. We revealed that TNT can readily form an aggregate in the aqueous phase and quickly approach the surface of the membrane. At low concentrations of TNT (7 mol%), the aggregate is unstable and breaks up after several nanoseconds, and then the dispersed TNT molecules enter the membrane alone. At high concentrations (14 mol%), the aggregate is adsorbed as a whole and remains stable inside the membrane. After some of the TNT is absorbed by the membrane, the remaining TNT across the membrane would have greater permeability, i.e., the calculated permeability coefficient (P) is increased from 1.7 × 10^-2 to 18.3 cm s^-1. Correspondingly, a higher bioconcentration factor (BCF) was also observed. The increased level is more pronounced in the presence of TNT aggregates (i.e., high concentrations). This phenomenon is closely related to the strong interaction between TNT molecules. The results suggested that TNT molecules that have entered into the membrane can facilitate the membrane uptake, permeation and bioaccumulation of subsequent TNT molecules, exhibiting a synergistic effect. This work has a certain significance for understanding the toxicity of TNT.

Development of liposome/water partition coefficients predictive models for neutral and ionogenic organic chemicals

Membrane/water partition coefficient (K_m/w) is a vital parameter used to characterize the membrane permeability of compounds. Considering the K_m/w value is difficult to observe experimentally for real biological membranes, liposome/water partition coefficient (K_lip/w) is employed to approximate K_m/w. Here, quantitative structure property relationship (QSPR) models for logK_lip/w of the neutral organic chemicals and the neutral form of ionogenic organic chemicals (IOCs) (logK_{lip/w-neutral}), ionic form of IOCs (logK_lip/w-ionic), the speciation-corrected liposome-water distribution ratios at a pH = 7.40 (logD_{lip/w-(pH=7.40)}) were developed. In the modeling, two modeling methods (multiple linear regressions (MLR) and k-nearest neighbor (kNN)) were used. The predictive variables employed here could be calculated from the molecular structure directly. For logK_{lip/w-neutral} and logD_{lip/w-(pH=7.40)}, the logK_OW and logD_OW-based, non-logK_OW and non-logD_OW-based kNN-QSPR and MLR-QSPR models were developed, respectively. The evaluation results implied that the predictive performance of kNN-QSPR models is better than that of MLR-QSPR models. For logK_lip/w-ionic, only one acceptable MLR-QSPR model was developed for cation and anion, respectively. The model quality of the derived models was evaluated following the OECD QSPR models validation guideline. The determination coefficient (R²), leave-one-out cross validation Q² (Q²_LOO) and bootstrapping coefficient (Q²_BOOT), the external validation coefficient (Q²_EXT) of all the models met the acceptable criteria (Q² > 0.600, R² > 0.700); while the root-mean-square error (RMSE) range from 0.351 to 0.857. All the results implied that the models had good goodness-of-fit, robustness and predictive ability. Therefore, the developed models could be used to fill the data gap for substances within the applicability domain on their missing logK_{lip/w-neutral}, logK_lip/w-ionic, logD_{lip/w-(pH=7.40)} values.

Diffusion of energetic compounds through biological membrane: Application of classical MD and COSMOmic approximations

Computational studies of the potential biological impact of several energetic compounds were performed. The most commonly used explosives were considered in the present studies: trinitrotoluene (TNT), 2,4-dinitrotoluene (2,4-DNT), 2,4-dinitroanisole (DNAN), and 5-Nitro-2,4-dihydro-3H-1,2,4-triazol-3-one (NTO). The effect of such factors as ionic strength and presence of DMSO in the water solution on the structure of the membrane were considered using the POPC lipid bilayer as an example. Molecular dynamics (MD) simulations revealed that, even on a short-time scale, the influence of those additives is noticeable, and therefore those factors should always be taken into account. The MD and the COSMOmic approaches were used to elucidate the ability of the energetic compounds to penetrate the living cell. Calculated free energy profiles and partitioning coefficients revealed distributions of the compounds in the lipid bilayer as well as an overall ability to enter the cell. MD in this case provides a better representation of the free energy profile, while the COSMOmic approach works better to predict log(K_lipw) values. The effect of the functional group was observed for the profiles that were obtained using the MD method.