Next Article in Journal
Understanding the Sensitive Skin Subject to Achieve a More Holistic Diagnosis
Previous Article in Journal
The Potential of Sky Fruit as an Anti-Aging and Wound Healing Cosmeceutical Agent
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Empirical Analysis Revealing Privileged Chemical Space of Cosmetic Preservatives

Ashland Specialties Deutschland GmbH, 22851 Norderstedt, Germany
*
Author to whom correspondence should be addressed.
Cosmetics 2021, 8(3), 80; https://doi.org/10.3390/cosmetics8030080
Submission received: 5 August 2021 / Revised: 24 August 2021 / Accepted: 26 August 2021 / Published: 31 August 2021

Abstract

:
Most cosmetic products require preservation to prevent microbial contamination and to ensure consumer safety. Due to regulatory restrictions and rejection by consumers, preservative options have become limited and the development of novel solutions is needed. This search can be guided by knowledge about favorable chemical space for cosmetic preservatives. Therefore, we used preservatives allowed in the EU as training set and calculated various molecular properties. Empirical analysis revealed two separated areas of privileged chemical space with the net charge as distinctive property. The first area comprises the group of neutral and anionic preservatives and is characterized by low molecular size as well as limited hydrogen-bonding capacity, polarity, and flexibility. The second area includes cationic preservatives, which are rather diffusely distributed regarding molecular weight and hydrogen-bonding, however, all members share high flexibility. Both groups significantly differ from antibiotics, reflecting the specific requirement of cosmetic preservation. The molecular properties defining the privileged chemical space are easy to calculate, and thus, can provide guidance for the development of novel preservatives.

Graphical Abstract

1. Introduction

Microbial contamination of cosmetic products is limiting shelf life and can pose a significant health risk for consumers [1]. To prevent spoilage by contamination with microorganisms, preservatives are typically added to ensure that the cosmetic composition remains free of microbial contamination for an adequately long period. Thereby, antimicrobial activity should be effective over a broad-spectrum of microorganisms and should persist longer than expected shelf-life plus the usage time. However, according to the Rapid Alert System (RAPEX) of the European Commission (EC), 61 cosmetic products were recalled during 2016 and 2020, as they present a serious risk due to microbiological contamination [2].
In the EU, Annex V of the Regulation (EC) No. 1223/2009 on cosmetic products provides a list of substances, which can be added to cosmetic products with the intention to preserve the composition from microbial contamination [3]. However, use concentrations are restricted, in many cases resulting in insufficient antimicrobial protection. Many ingredients negatively impact odor or color of the cosmetic formulation, and disturbance of the emulsifying system can lead to significant change of viscosity or even ‘breaking’ of emulsions. Furthermore, some preservatives are, despite of their documented safety, not well accepted by consumers. Formaldehyde releasers are a major class of Annex V preservatives, however, reports about carcinogenicity at high levels [4] and skin sensitizing potential [5], which were regulatory considered by maximum use concentrations, led to the disappearance in cosmetic products. Both the antimicrobial effect, as well as toxicological effects on humans, are based on alkylation reactions reflecting the high electrophilicity of formaldehyde [6,7]. Similarly, isothiazolinones, sulfur dioxide released by sulfites, organomercuric agents, bronidox, bronopol, and silver react with biopolymers, in particular with thiol groups of proteins [8,9,10,11,12,13].
Combining limited applicability of many preservatives and rejection by consumers, there is a depletion of usable preservatives, and thus, modern options that are not based on chemical reactivity are needed to ensure adequate preservation of cosmetics. The search for novel antimicrobials can be guided by analysis of preservatives with accepted mechanisms-of-action.
Molecular properties of substances are determinant factors for interaction with their immediate surrounding. This becomes impressively clear by the pioneering work of Lipinski and coworkers, who analyzed the influence of various chemical and physical properties on oral bioavailability of drugs in humans [14]. Thereby, they formulated the “rule of five” for molecular properties of drugs, stating that poor absorption or permeation after oral intake in humans is more likely when (1) there are more than 5 H-bond donors; (2) the molecular weight is over 500 g/mol; (3) the logP is over 5; (4) there are more than 10 H-bond acceptors (expressed as the sum of nitrogen and oxygen atoms). Various adaptations and extensions were proposed in the field of drug discovery. Veber et al. noticed the importance of the number of rotatable bonds (nrotb) and topical polar surface area (tPSA) as determining factors for good oral bioavailability in rats [15].
Regarding preservation of cosmetic formulations, it seems intuitive to us that molecular properties of substances have direct influence on antimicrobial activity on two levels:
(1) Availability of antimicrobial substance at their target site: most cosmetic formulations are complex systems composed of multiple phases. For example, many creams and lotions are o/w emulsions containing a liquid oily phase dispersed in a continuous aqueous phase stabilized by emulsifiers. Shampoos and body washes are also inhomogeneous systems as they contain surfactants, which can form aggregates in water, such as micelles. To be effective against microbes, preservatives need to get in close contact to microbial cells to allow adequate interaction, i.e., they must be present in sufficient concentration in the bulk water. Various studies revealed lower antimicrobial activity in cosmetic formulations compared to aqueous solution [16]. For example, Tschierske et al. has shown that individual interactions between an antimicrobial agent and specific emulsifiers determine the degree of binding, which reduces the availability for target interaction and leads to a loss in activity [17].
(2) Preservative/microbe interaction: the interaction between the antimicrobial substance and the microbe is expected to strongly depend on molecular properties. Even if preservatives in cosmetics differ in their specific mechanism-of-action, they all have in common that there is not a specific microbial binding site [18]. Many antimicrobial agents rather bind promiscuously to microbial cell components as proteins or membranes. The interaction of non-covalent binders is rather short-term and multiple molecules can be involved. This is in strong contrast to most drugs, which intentionally attack a specific binding site of a single molecular target as enzymes or receptors. Fine-tuned shape and specific interaction lead to molecular recognition. Taken together, this indicates that the interaction between microbes and preservatives examined in this article is less driven by the presence of specific functional groups in a specific spatial constellation (as symbolized by pharmacophore model), but rather by more general molecular properties.
The objective of this study was to identify beneficial molecular characteristics of antimicrobial substances, which are successfully used in cosmetic products for preservations purposes. Therefore, we studied molecular characteristics of preservatives allowed in the EU. Our results highlight, that there are two distinctive privileged areas in chemical space, which can be defined by few simple parameters. These findings provide a guide for identification and optimization of new antimicrobial substances applicable for cosmetic preservation.

2. Materials and Methods

The list of preservatives allowed in cosmetic products (Annex V, Regulation 1223/2009/EC on Cosmetic Products [3], version of March 2018) was used as starting point for the analysis. When entries contain various salt forms, only one form was considered as no significant influence on activity is expected. Some entries contain derivatives of different alkyl chain lengths. In these cases, representative chain lengths (e.g., C12, C16, and C22 for Alkyl (C12–C22) trimethyl ammonium) were included in the analysis.
Substances were transformed into two relevant species: (1) the predominant species in cosmetic formulations at pH 6. Deprotonated species are used if pKa < 6. Dissociation of salts is expected, if no opposite information is available; (2) the active species interacting with microbials. Exemplary, organic acids as benzoic acid are thought to permeate as neutral species, even if they predominantly exist as anion at pH 6. Substances used for analysis as well as structure of predominant and active species can be found in Supporting Information Table S1.
Substances with a reactivity-based mechanism-of-action were excluded from closer empirical analysis. Information about the acceptance of underlying mechanisms are specified in the Supporting Information Table S1.
Various molecular properties were considered in the empirical analysis. Structural descriptors were derived from Lewis structure of an object. The formal net charge (Q) and molecular weight (MW) were calculated using ChemDraw® Version 19.1.1.21 (PerkinElmer Informatics, Inc., Waltham, MA, USA). The number of rotatable bonds (nrotb) was calculated using Chem3D® Version 19.1.1.21 (PerkinElmer Informatics, Inc., Waltham, MA, USA). The number of hydrogen-bond donors (HBD) was determined by counting all H-atoms bound to nitrogen or oxygen atoms in the Lewis structure. The number of hydrogen-bond acceptors (HBA) was determined by counting all nitrogen or oxygen atoms bearing lone electron pair in the Lewis structure to participate in hydrogen bonds. Hydrogen bonding capacity (HBDnA) was defined as sum of HBD and HBA. Descriptors of surface and shape as Connolly molecular surface area (MSA) and topological polar surface area (tPSA) were calculated using Chem3D® Version 19.1.1.21 (PerkinElmer Informatics, Inc., Waltham, MA, USA). Before calculating MSA, structures were optimized using MM2 energy minimization using Minimum RMS Gradient of 0.010. MSA was calculated using a probe radius of 1.4. Physico-chemical parameters including distribution coefficient (logD) and aqueous solubility (logS) were calculated by ChemAxon logD calculator [19] or Chem3D® Version 19.1.1.21 (PerkinElmer Informatics, Inc., Waltham, MA, USA), respectively. For logS calculation, the predominant species in water at pH 6 was used to calculate solubility, as the carboxylate mainly contributes to solubility and availability in water phase and fast conversion to the active species is expected. For visualization, molecular characteristics were plotted using Microsoft 365® Excel® (Microsoft Corporation).

3. Results

3.1. Substance Set Description

All entries listed as preservatives in Annex V of the Regulation (EC) No. 1223/2009 [3] were analyzed in this study (Supporting Information Table S1). Many substances are hardly accepted by consumers and several of these substances have a chemical reactivity-based mechanism-of-action. Thus, aldehydes and formaldehyde releaser, isothiazolinones, sulfur dioxide releasing sulfites, and substances based on heavy metals (such as mercury and silver) were excluded from analysis, since their mechanisms-of-action are irrelevant for future developments.

3.2. Evaluation of Chemical Space Reveals Two Distinctive Areas

Forty substances with accepted mechanism-of-action were analyzed regarding molecular properties (Table 1). A plot of total polar surface area (tPSA) versus Connolly molecular surface area (MSA) shows two distinctive groups (Figure 1). The first condensed group A (blue) is characterized by MSA between 49 and 263 Å2 and tPSA between 20 and 92 Å2. The second group B (red) is rather diffuse revealing MSA ≥ 320 Å2 and no limitations for tPSA. No upper limit is given due to the polymer polyaminopropyl biguanide with a wide scattering of MSA and tPSA depending on varying chain length.
The distinctive property separating both groups is apparently the net charge (Q), as members of group B are cationic (Q ≥ 1). Members of group A are either neutral or anionic in their predominant form at pH 6. As all anionics of group A have pKa values (of the conjugated acid) of at least 3.9, a significant percentage exist in their un-ionized form.

3.3. Characterization of Molecular Chemical Space: Structural Descriptors

3.3.1. Molecular Size

MSA is a measure of the size of a molecule. Alternatively, molecular weight (MW) can be used as indicator. Thereby, the distinction of both groups is less pronounced compared to MSA (Figure 2a). MW spans 45 to 316 g/mol for group A, excluding bromochlorophene wearing two heavy bromine atoms, and MW ≥ 257 g/mol for cationics of group B. However, molecular weight is very easy to calculate and might be valuable as size indicator when searching for novel preservatives in both privileged areas.

3.3.2. Flexibility

The number of rotatable bond (nrotb) is a measure for a molecule’s flexibility. Conformational adaption might facilitate target binding, however, restriction of rotatable bonds upon target binding leads to reduced affinity due to an entropic penalty [20]. Furthermore, flexibility allows conformations with different polar surface areas, having influence on aqueous solubility of a substance and solubilization of a substance-target adduct. Group A preservatives reveal nrotb ≤ 6 with the exemption of undec-10-enoic acid, which shows higher flexibility (nrotb = 9). In contrast, cationic substances of group B contain at least 11 rotatable bonds (Figure 2b).
Table 1. Calculated molecular descriptors of preservatives, which were used for empirical analysis of privileged chemical space.
Table 1. Calculated molecular descriptors of preservatives, which were used for empirical analysis of privileged chemical space.
Annex V
Ref. No.
Preservative 1Q 2,3MSA 1
2)
tPSA 1
2)
MW 1
(g/mol)
nrotb 1HBD 1HBA 1HBDnA 1logDlogS 1MoA 3Group
1Benzoic acid and its salts−1 (0)119401211022−0.29−1.39Protonophores; causing membrane depolariz. [21] A
2Propionic acid and its salts−1 (0)8940731022−0.790.12
3Salicylic acid and its salts−1 (0)108601371133−1.06−1.55
4Sorbic acid and its salts−1 (0)1294011120220.20−1.07
8Pyrithione zinc (here: pyrithione)−1 (0)1153516200220.14−0.64
13Dehydroacetic acid and its salts−1 (0)155661671044−0.16−0.73
14Formic acid and its salts−1 (0)4940450022−2.000.54
18Undec-10-enoic acid and its salts−1 42324018390222.71−3.01- 4A
12 bisButyl 4-hydroxybenzoate02094719251343.00−2.64Mechanism
is unclear 5
A
12 bisPropyl 4-hydroxybenzoate01904718041342.55−2.22
12Methylparaben01514715221341.67−1.60
12Ethylparaben01714716631342.03−1.94
7Biphenyl-2-ol01692017011123.32−3.34Solubilization of lipids and
denaturing proteins [21]
A
11Chlorobutanol01252017711121.75−1.95
222,4-Dichlorobenzyl alcohol01532017711122.41−2.44
24Chlorocresol01372014201122.79−2.21
26Chloroxylenol01542015701123.30−2.44
292-Phenoxyethanol01532913831231.13−1.10
34Benzyl alcohol01252010811121.21−1.02
37Bromochlorophen02624042722245.48−5.91
384-Isopropyl-m-cresol01712015011123.43−2.79
40Chlorophene02052021921124.37−3.92
431-Phenoxypropan-2-ol01692915231231.54−1.49
50Chlorphenesin01935020242351.10−1.50
563-Iodo-2-propynylbutylcarbamate02183828161342.54−3.14-A
23Triclocarban02314131622354.93−5.48Intracell. enzyme inhibitors [22,23,24,25]A
25Triclosan02243029021234.97−5.19
32Climbazol02634229350443.07−3.96
35Piroctone and its monoethanolamine salt−12464323640223.13−3.62
19Hexetidine145637341123252.75−4.67B
15Dibromohexamidine and its salts2412122514118412−0.53−8.08Interaction with membrane component
destroying membrane integrity [21]
B
28Polyaminopropyl biguanide (here: n = 22)2249321881400017613288220-−2.82
42Chlorhexidin and its salts24881815071312820−3.90−7.52
44Alkyl (C12-C22) trimethyl ammonium salts 6 (here: C12)13200257130000.91−3.25
44Alkyl (C12-C22) trimethyl ammonium salts 6 (here: C16)13940285150002.69−4.61
44Alkyl (C12-C22) trimethyl ammonium salts 6 (here: C22)15070369210005.35−6.64
47Hexamidine and its salts2371122356118412−2.06−6.50
53Benzethonium chloride139619448120221.78−5.32
54Benzalkonium salts 7 (here: C12)13790304130003.63−5.83
58Ethyl Lauroyl Arginate Hydrochloride14761213861966121.12−5.48
1 Predominant species in water at pH 6 was used for calculation. Active species only differs significantly for protonophores. 2 Net charge Q of active form is given in bracket, if differing from net charge of predominant form. 3 Since preservative act unspecific, mechanism-of-action and cellular changes resulting from exposure to the substance were roughly categorized. 4 Undec-10-enoic acid is often described to act as protonophore in analogy to other organic acids [21], but alternative mechanism are under discussion [26]. Since we did not observe pH-dependency as expected for protonophores (data not shown), we consider the carboxylate as active species. 5 For parabens, various mechanism-of-actions are discusses as disrupting of membrane transport processes, inhibition of DNA and RNA synthesis, and inhibition of ATPases and phosphotransferases has been discussed [27] leading to membrane depolarization [21]. 6 Chain lengths C12-C22 are authorized for use as preservatives. Three derivatives representing the full range of chain length were used here. 7 Benzalkonium salts (alkyl dimethyl benzyl ammonium salts with alkyl residues of ranging from C8 to C18 are authorized for use as preservatives. However, greatest biocide activity is associated with C12-C14 derivatives, which are the main components [28].

3.3.3. Hydrogen Bonding

Hydrogen bonding strongly contributes to the interaction capacity of a substance with a biomolecular target and surrounding solvent. This has impact on affinity to microbial targets, but also on solubility and distribution within a cosmetic formulation. Therefore, we analyzed the overall hydrogen-bonding capacity (HBDnA) as calculated by the sum of hydrogen-bond donors (HBD) and hydrogen-bond acceptors (HBA) (Figure 2c). Group A revealed HBDnA between 2 and 6. Cationic group B did not show limitations as HBDnA can be very high for the polymer polyaminopropyl biguanide, whereas quaternary ammonium compounds cannot form hydrogen bonds. Individual numbers of HBD or HBA revealed same pattern of distribution (Supporting Information Figure S2), and thus, do not provide additional value.

3.4. Physicochemical Properties: Lipophilicity and Solubility

3.4.1. Lipophilicity

Lipophilicity is a crucial property as it represents the affinity of a molecule for a lipophilic environment. Keeping in mind that cosmetic formulations usually consist of multiple phases, a substance’s lipophilicity has direct impact on the distribution and availability at the target site. Furthermore, lipophilicity affects the affinity of a molecule for microbial binding sites.
Many preservatives analyzed in this study are ionizable. Since neutral and ionic forms exhibit different polarities, a substance’s lipophilicity depends on pH. Thus, we calculated the distribution coefficient logD for evaluation of lipophilicity, as it considers all species of a molecule at a given pH, which was representatively set to pH 6 in this study. Both groups A and B revealed similar distribution of logD values in the range between −4 and 5.5 (Figure 2d).

3.4.2. Solubility

Aqueous solubility is an upper limit for the concentration of a substance in aqueous phase. Since antimicrobials are expected to be located in aqueous phase to be efficient [17], solubility is crucial for antimicrobial activity. Calculated logS values lie between −6 and 0.5 for group A, and between −2.8 and 8.1 for cationic substances of group B (Figure 2d). However, some cationics have amphiphilic properties and might form aggregates in bulk water. Furthermore, it is important to note that calculated solubility might be inaccurate in some cases compared to experimental values. Nevertheless, the calculated logS range of preservatives can be used in the search for novel preservatives when compared to calculated logS of candidates.

4. Discussion

4.1. Calculation and Analysis of Molecular Desciptors

In analogy to the concept of drug-likeness, we expected that there is a privileged chemical space for cosmetic preservatives. Therefore, various molecular descriptors and properties were analyzed.
Preservatives in cosmetic products might exist in various forms as different protonation states or even might undergo chemical reactions. Many of the molecular descriptors examined depend on the molecular species used for calculation. There are two forms pivotal for antimicrobial activity in cosmetics: (1) the predominant species in cosmetic formulations, which determines the substance’s location within a cosmetic formulation, and thus, availability at the target site; (2) the active species interacting with microbes, which is the crucial form for antimicrobial activity. The original form as added to the cosmetic formulation is considered as irrelevant as fast conversion to the predominant species is expected.
Conspicuously, for all substances with accepted mechanism-of-action, the predominant and active species differed in protonation state only. As there are only minor differences in molecular descriptors, we have focused on the predominant form if not stated otherwise. For calculation of aqueous solubility, which might be considered as an upper limit for the concentration in the aqueous phase, use of the predominant species is reasonable as this species is, per definition, the thermodynamically favored one contributing most to the overall solubility.
To measure lipophilicity of drug candidates, Lipinski and others used the partition coefficient logP, which refers to the partitioning of un-ionized compound between octan-1-ol and water. For drugs this appears reasonable due to the hypothesis, that absorption of ionizable drugs takes place where the local pH provides the maximum concentration of non-ionized form [29]. However, the effectiveness of preservatives in cosmetics depends on distribution into the aqueous phase [17]. As fast conversion between ionized and un-ionized form is expected, all species are relevant. Thus, the distribution coefficient logD, which considers all forms of a compounds (ionized plus un-ionized) at the given pH (here: pH 6) was used as descriptor for lipophilicity in this study. One might argue that for interaction with the microbial target, only the lipophilicity of the active form of an antimicrobial is relevant. This applies to organic acids only, as they penetrate as neutral molecules into microbials cells even if the predominant species is ionized. Indeed, clogP values of organic acids were increased compared to ionized carboxylates but remain in a typical range for group A members.

4.2. Dependencies of Molecular Descriptors

Among various molecular properties analyzed in this study, close correlations were found. As stated above, molecular weight and molecular surface area are both measures of size, and thus, intrinsically correlate to some degree. To keep the rules easily applicable for future searches, definition of privileged chemical space should be limited to independent properties providing additional value. From this perspective, consideration of size can be limited to molecular weight, which is easily accessible.
In a modification of Lipinski’s rules, number of rotatable bonds rather than molecular weight was proposed as crucial measure for oral drug bioavailability [15]. As the chance for high number of rotatable bonds increases with number of bonds, some degree of dependency between nrotb and MW can be expected. However, we found a wide distribution for nrotb/MW, demonstrating independence of flexibility and size (Supporting Information Figure S2). In analogy to nrotb, we found low nrotb/MW for neutral and anionic molecules, and significantly higher nrotb/MW for cationics, legitimating the use of nrotb to define the areas of preservatives in chemical space.
A correlation between total polar surface area and interaction capacity, more precisely the sum of hydrogen-bonding donors and acceptors, was identified. tPSA as typically calculated arises from oxygen and nitrogen atoms including hydrogen attached to oxygen and nitrogen [30], which are the moieties contributing to hydrogen bonding. Thus, there is a causal association and consideration can be limited to tPSA, as hydrogen-bonding alone does not give additional value [20]. Taking into account further interactions, e.g., ionic or van der Waals forces, leading to a hypothetic overall interaction capacity is expected to reduce dependency on tPSA.

4.3. Mechanism-of-Action and Rules for Privileged Chemical Space

We suppose that a set of compounds sharing the same mechanism-of-action also share beneficial molecular characteristics. However, in contrast to antibiotics acting on specific sites of biosynthetic processes, preservatives have multiple targets [10]. At microbicidal concentrations, direct and indirect cellular effects can hardly be distinguished. Thus, we categorized all preservatives in rough classes (Table 1) and analyzed their position in chemical space.
As stated above, the cationic group interacting with membrane components occupy a distinctive area in chemical space. Hexetidine, which almost entirely exist in cationic form (pKa 8.3), shares key characteristics of this group as charge, high molecular size, and high flexibility, and thus, strongly differs from the members of intracellular inhibitor groups (see Supporting Information Figure S3).
All mechanism-of-action categories bearing neutral and anionic substances are in a common area of chemical space. There are only trends as protonophores tend to be small and share low logD and high logS values due to predominant existence as carboxylate. Small size and pKa values, which allow a minor part to exist as neutral molecules, enable penetration into microbial cells. However, there is no clear distinction but rather strong overlap with other categories.
Interestingly, undec-10-enoic acid reveals significant higher flexibility compared to other group members (nrotb = 9 vs. ≤6). Furthermore, unlike for protonophores we did not observe tremendous activity drop when pH is increased above 6, which indicates that the ionized carboxylate is the active species (data not shown). Both observations might support alternative mechanism-of-actions as cytoplasmic membrane disruption [26]. In any case, anionic substances bearing a high flexibility might represent an under-exploited category.
Taking our findings together, we postulate a privileged area for neutral or anionic cosmetic preservation narrowed down to MW = 40–320 g/mol, nrotb ≤ 6, tPSA = 20–100 Å2, logD = −4–6, logS = −6–6. For cationics, privileged area is characterized by MW ≥ 250 g/mol, nrotb ≥ 11, no limitations for tPSA, logD = −4–6, logS = −3–9. This does not mean that substances fulfilling these criteria are necessarily strong antimicrobials; however, substances in privileged chemical space have been shown to comply with requirements of cosmetic preservatives, which includes availability at target site. Thus, the criteria can serve as guide for future developments.

4.4. Comparison of Cosmetic Preservatives with Marketed Drugs and Antibiotics

Lipinski’s work marked the starting point for the analysis of physicochemical properties of various substance sets including marketed antibiotics [5,31,32]. Keeping the fundamental differences of application between antibiotics and cosmetic preservatives in mind, we compared molecular properties of both substance sets (Table 2).
Obviously, compared to antibiotic substance sets, neutral or anionic cosmetic preservatives are simple molecules characterized by small size and low hydrogen-bonding capabilities. Antibiotics, in particular if active against Gram-positive bacteria, have higher molecular weight, and higher numbers of hydrogen-bond donors and acceptors. An overall increased polarity is also reflected by increased tPSA and decreased logD.
Cationic preservatives have similar size as antibiotics for Gram-negative bacteria. However, they are characterized by higher proportion of rotatable bonds, and lower polarity as revealed by total polar surface area, hydrogen-bonding capacity, and distribution coefficient logD.
Thus, cosmetic preservatives occupy a completely different area in chemical space than antibiotics do. This might reflect the differences in mechanism-of-action of antibiotics and preservatives, as the latter do not need multiple interactions for specific target recognition. Molecular properties of the preservative’s privileged chemical spaces rather combine two crucial requirements, availability at the site of action and favorable interactions with microbes. This can guide future efforts to identify new ingredients for cosmetic preservation.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cosmetics8030080/s1, Table S1: List of substances including acceptance of mechanism-of-action and structures; Figure S2: Additional plots; Figure S3: Plots of MoA categories.

Author Contributions

Conceptualization, M.P.S.; Methodology, M.P.S.; Writing—original draft preparation, M.P.S.; Writing—review and editing, M.P.S. and L.H.; Visualization, M.P.S.; Project administration, M.P.S. and L.H.; Funding Acquisition, M.P.S. and L.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Union—Regional Development Fund, grant number 095 18 001, project number LPW-E/1.2.3/895.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Becks, V.E.; Lorenzoni, N.M. Pseudomonas aeruginosa Outbreak in a Neonatal Intensive Care Unit: A Possible Link to Contaminated Hand Lotion. Am. J. Infect. Control 1995, 23, 396–398. [Google Scholar] [CrossRef]
  2. Safety Gate: The EU Rapid Alert System for Dangerous Non-Food Products. Available online: https://ec.europa.eu/safety-gate-alerts/screen/webReportusingsearchforalertoption (accessed on 22 April 2021).
  3. The European Parliament and the Council of the European Union. Regulation (EC) No. 1223/2009 of the European Parliament and of the Council of 30 November 2009 on Cosmetic Products–Consolidated Version of 25/12/2017. Off. J. Eur. Union Lex 2009, 342, 59.
  4. McLaughlin, J.K. Formaldehyde and cancer: A critical review. Int. Arch. Occup. Environ. Heath 1994, 66, 295–301. [Google Scholar] [CrossRef]
  5. De Groot, A.C.; Flyvholm, M.-A.; Lensen, G.; Menné, T.; Coenraads, P.-J. Formaldehyde-releasers: Relationship to formal-dehyde contact allergy. Contact allergy to formaldehyde and inventory of formaldehyde-releasers. Contact Dermat. 2009, 61, 63–85. [Google Scholar] [CrossRef] [Green Version]
  6. Denyer, S.P. Mechanisms of Action of Antibacterial Biocides. Int. Biodeterior. 1995, 36, 221–245. [Google Scholar] [CrossRef]
  7. Swenberg, J.A.; Moeller, B.C.; Lu, K.; Rager, J.E.; Fry, R.C.; Starr, T.B. Formaldehyde Carcinogenicity Research: 30 Years and Counting for Mode of Action, Epidemiology, and Cancer Risk Assessment. Toxicol. Pathol. 2013, 41, 181–189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Williams, T.M. The Mechanism of Action of Isothiazolone Biocides. Power Plant Chem. 2007, 9, 14–22. [Google Scholar]
  9. Irwin, S.V.; Fisher, P.; Graham, E.; Malek, A.; Robidoux, A. Sulfites Inhibit the Growth of Four Species of Beneficial Gut Bacteria at Concentrations Regarded as Safe for Food. PLoS ONE 2017, 12, e0186629. [Google Scholar] [CrossRef] [Green Version]
  10. Halla, N.; Fernandes, I.P.; Heleno, S.A.; Costa, P.; Boucherit-Otmani, Z.; Boucherit, K.; Rodrigues, A.E.; Ferreira, I.C.F.R.; Barreiro, M.F. Cosmetics Preservation: A Review on Present Strategies. Molecules 2018, 23, 1571. [Google Scholar] [CrossRef] [Green Version]
  11. Shepherd, J.A.; Waigh, R.D.; Gilbert, P. Antibacterial Action of 2-Bromo-2-Nitropropane-1,3-Diol (Bronopol). Antimicrob. Agents Chemother. 1988, 32, 1693–1698. [Google Scholar] [CrossRef] [Green Version]
  12. Ghannoum, M.; Thomson, M.; Bowman, W.; Al-Khalil, S. Mode of Action of the Antimicrobial Compound 5-Bromo-5-nitro-1,3-dioxane (Bronidox). Folia Microbiol. 1986, 31, 19–31. [Google Scholar] [CrossRef] [PubMed]
  13. Yin, E.X.; Zhang, J.; Zhao, E.S.; Mei, M.L.; Li, Q.; Chu, C.H. The Antibacterial Mechanism of Silver Nanoparticles and Its Application in Dentistry. Int. J. Nanomed. 2020, 15, 2555–2562. [Google Scholar] [CrossRef] [Green Version]
  14. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and Computational Approaches to Estimate Solubility and Permeability in Drug Discovery and Development Settings. Adv. Drug Deliv. Rev. 1997, 23, 3–25. [Google Scholar] [CrossRef]
  15. Veber, D.F.; Johnson, S.R.; Cheng, H.-Y.; Smith, B.R.; Ward, K.W.; Kopple, K.D. Molecular Properties That Influence the Oral Bioavailability of Drug Candidates. J. Med. Chem. 2002, 45, 2615–2623. [Google Scholar] [CrossRef] [PubMed]
  16. Maillard, J.-Y. Factors Affecting the Activities of Microbicides. In Hugo and Ayliffe’s Principles and Practice of Disinfection, Preservation and Sterilization, 5th ed.; Fraise, A.P., Maillard, J.-Y., Sattar, S.A., Eds.; John Wiley & Sons: Hoboken, NJ, USA, 2013; pp. 71–86. [Google Scholar]
  17. Tschierske, N.; Emeis, D.; Biel, S.S. Efficacy Loss of Antimicrobial Actives in Emulsions. SOFW-J. 2012, 4, 1–4. [Google Scholar]
  18. Denyer, S.P. Mechanism of Action of Biocides. Int. Biodeterior. 1990, 26, 89–100. [Google Scholar] [CrossRef]
  19. Available online: https://disco.chemaxon.com/calculators/demo/plugins/logd/ (accessed on 15 April 2021).
  20. Bradley, C.D.; Kihlberg, J. Drug Discovery Beyond the Rule of 5-Opportunities and Challenges. Expert Opin. Drug Discov. 2017, 12, 115–119. [Google Scholar] [CrossRef] [Green Version]
  21. Geis, P.A. Appendix: Common Cosmetic Preservatives. In Cosmetic Microbiology A Practical Guide, 2nd ed.; Geis, P.A., Ed.; Taylor & Francis Group: New York, NY, USA, 2006; pp. 228–281. [Google Scholar]
  22. Macsics, R.; Hackl, M.W.; Fetzer, C.; Mostert, D.; Bender, J.; Layer, F.; Sieber, S.A. Comparative Target Analysis of Chlorinated Biphenyl Antimicrobials Highlights MenG as a Molecular Target of Triclocarban. Appl. Environ. Microbiol. 2020, 86, e00933-20. [Google Scholar] [CrossRef]
  23. Russell, D. Whither triclosan? J. Antimicrob. Chemother. 2004, 53, 693–695. [Google Scholar] [CrossRef] [Green Version]
  24. Richter, E.; Wick, A.; Ternes, T.A.; Coors, A. Ecotoxicity of Climbazole, a Fungicide Contained in Antidandruff Shampoo. Environ. Toxicol. Chem. 2013, 32, 2816–2825. [Google Scholar] [CrossRef]
  25. do Couto, F.M.; do Nascimento, S.C.; Júnior, S.F.; da Silva, V.K.; Leal, A.F.; Neves, R.P. Antifungal Activity of the Piroctone Olamine in Experimental Intra-Abdominal Candidiasis. SpringerPlus 2016, 5, 468. [Google Scholar] [CrossRef] [Green Version]
  26. Casillas-Vargas, G.; Ocasio-Malavé, C.; Medina, S.; Morales-Guzmán, C.; Del Valle, R.G.; Carballeira, N.M.; Sanabria-Ríos, D.J. Antibacterial Fatty Acids: An update of Possible Mechanisms of Action and Implications in the Development of the Next-generation of Antibacterial Agents. Prog. Lipid Res. 2021, 82, 101093. [Google Scholar] [CrossRef]
  27. Valkova, N.; Lépine, F.; Valeanu, L.; Dupont, M.; Labrie, L.; Bisaillon, J.G.; Beaudet, R.; Shareck, F.; Villemur, R. Hydrolysis of 4-Hydroxybenzoic Acid Esters (Parabens) and Their Aerobic Transformation into Phenol by the Resistant Enterobacter cloacae Strain EM. Appl. Environ. Microbiol. 2001, 67, 2404–2409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Pereira, B.M.P.; Tagkopoulos, I. Benzalkonium Chlorides: Uses, Regulatory Status, and Microbial Resistance. Appl. Environ. Microbiol. 2019, 85, e00377-19. [Google Scholar]
  29. Constantinescu, T.; Lungu, C.N.; Lung, I. Lipophilicity as a Central Component of Drug-Like Properties of Chalchones and Flavonoid Derivatives. Molecules 2019, 24, 1505. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Caron, G.; Ermondi, G. Molecular Descriptors for Polarity: The Need for Going Beyond Polar Surface Area. Future Med. Chem. 2016, 8, 2013–2016. [Google Scholar] [CrossRef] [PubMed]
  31. O’Shea, R.; Moser, H.E. Physicochemical Properties of Antibacterial Compounds: Implications for Drug Discovery. J. Med. Chem. 2008, 51, 2871–2878. [Google Scholar] [CrossRef] [PubMed]
  32. Ebejer, J.-P.; Charlton, M.H.; Finn, P.W. Are the Physicochemical Properties of Antibacterial Compounds Really Different From Other Drugs? J. Cheminform. 2016, 8, 30. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Two-dimensional plots of total polar surface area (tPSA) vs. Connolly molecular surface area (MSA): position of preservatives reveal two distinctive areas in chemical space for neutral or anionic (blue) and cationic (red) substances.
Figure 1. Two-dimensional plots of total polar surface area (tPSA) vs. Connolly molecular surface area (MSA): position of preservatives reveal two distinctive areas in chemical space for neutral or anionic (blue) and cationic (red) substances.
Cosmetics 08 00080 g001
Figure 2. Two-dimensional plots including position of neutral or anionic (blue) and cationic (red) preservatives in chemical space: (a) plot of Connolly molecular surface area (MSA) vs. molecular weight (MW) as easily available alternative measure of size; (b) plot of number of rotatable bonds (nrotb) vs. MW; (c) plot of hydrogen bonding capacity (HBDnA) vs. MW; (d) plot of distribution coefficient (logD) vs. aqueous solubility (logS).
Figure 2. Two-dimensional plots including position of neutral or anionic (blue) and cationic (red) preservatives in chemical space: (a) plot of Connolly molecular surface area (MSA) vs. molecular weight (MW) as easily available alternative measure of size; (b) plot of number of rotatable bonds (nrotb) vs. MW; (c) plot of hydrogen bonding capacity (HBDnA) vs. MW; (d) plot of distribution coefficient (logD) vs. aqueous solubility (logS).
Cosmetics 08 00080 g002
Table 2. Average values of molecular characteristics. Comparison between cosmetic preservatives and antimicrobial drugs.
Table 2. Average values of molecular characteristics. Comparison between cosmetic preservatives and antimicrobial drugs.
Molecular
Property
Annex V 1
(Q = −1; 0)
Annex V 1,2
(Q ≥ 1)
Antibiotics 3
(Gram-Positive)
Antibiotics 3
(Gram-Negative)
Non-Antibiotics 3
MW184377813414338
nrotb2.314.013.86.85.9
tPSA366024316570
HBDnA2.96.323.414.56.5
logD1.981.17−0.20−2.791.59
logS−2.33−5.79---
1 Cosmetic preservatives with accepted mechanism-of-action as analyzed in this study and listed in Table 1. 2 For calculation of average values, the polymer polyaminopropyl biguanide was not considered. 3 Average values as calculated by O’Shea et al. [31]. Minor differences between the calculation process might result in slightly different numbers.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Storz, M.P.; Holsten, L. Empirical Analysis Revealing Privileged Chemical Space of Cosmetic Preservatives. Cosmetics 2021, 8, 80. https://doi.org/10.3390/cosmetics8030080

AMA Style

Storz MP, Holsten L. Empirical Analysis Revealing Privileged Chemical Space of Cosmetic Preservatives. Cosmetics. 2021; 8(3):80. https://doi.org/10.3390/cosmetics8030080

Chicago/Turabian Style

Storz, Michael P., and Lea Holsten. 2021. "Empirical Analysis Revealing Privileged Chemical Space of Cosmetic Preservatives" Cosmetics 8, no. 3: 80. https://doi.org/10.3390/cosmetics8030080

APA Style

Storz, M. P., & Holsten, L. (2021). Empirical Analysis Revealing Privileged Chemical Space of Cosmetic Preservatives. Cosmetics, 8(3), 80. https://doi.org/10.3390/cosmetics8030080

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop