In silico Assessment of Pharmacological Profile of Low Molecular Weight Oligo-Hydroxyalkanoates

Polyhydroxyalkanoates (PHAs) are a large class of polyesters that are biosynthesized by microorganisms at large molecular weights (Mw > 80 kDa) and have a great potential for medical applications because of their recognized biocompatibility. Among PHAs, poly(3-hydroxybutyrate), poly(4-hydroxybutyrate), poly(3-hydroxyvalerate), poly(4-hydroxyvalerate), and their copolymers are proposed to be used in biomedicine, but only poly(4-hydroxybutyrate) has been certified for medical application. Along with the hydrolysis of these polymers, low molecular weight oligomers are released typically. In this study, we have used a computational approach to assess the absorption, distribution, metabolism, and excretion (ADME)-Tox profiles of low molecular weight oligomers (≤32 units) consisting of 3-hydroxybutyrate, 4-hydroxybutyrate, 3-hydroxyvalerate, 4-hydroxyvalerate, 3-hydroxybutyrate-co-3-hydroxyvalerate, and the hypothetical PHA consisting of 4-hydroxybutyrate-co-4-hydroxyvalerate. According to our simulations, these oligomers do not show cardiotoxicity, hepatotoxicity, carcinogenicity or mutagenicity, and are neither substrates nor inhibitors of the cytochromes involved in the xenobiotic’s metabolism. They also do not affect the human organic cation transporter 2 (OCT2). However, they are considered to be inhibitors of the organic anion transporters OATP1B1, and OATP1B3. In addition, they may produce eye irritation, and corrosion, skin irritation and have a low antagonistic effect on the androgen receptor.


INTRODUCTION
Polyhydroxyalkanoates (PHAs) are polyesters produced by various microorganisms and serve them as intracellular carbon and energy storage compounds (Zinn et al., 2001). After chemical extraction and purification, PHAs have numerous industrial (e.g., packaging) but also very interesting high value applications in medicine due to their good biocompatibility and biodegradability. Thus, they have been proposed as drug delivery systems (Shrivastav et al., 2013;Radu et al., 2017;Butt et al., 2018;Zhang et al., 2018), implants for bone and cartilage regeneration (Puppi et al., 2010;Ray and Kalia, 2017), tissue engineering matrices (Ray and Kalia, 2017;Butt et al., 2018;Zhang et al., 2018), and anticancer agents (Ray and Kalia, 2017). In 2007, the United States Food and Drug Administration (FDA) approved the first PHA, poly(4-hydroxybutyrate) (P4HB) for use as a surgical suture (Williams et al., 2013).
It has been found that the chemical property of the monomeric unit(s) and their distribution in the polymer play an important role on the resulting material properties (crystallinity, tensile strength, melting endotherm, etc.). Consequently, degradation speed may vary and is also influenced by the dimensions of the PHA material (e.g., bone implant or nano-sized drug delivery beads). In contrast to poly(L-lactic acid), its degradation process is in addition to hydrolysis also enhanced by enzymatic degradation. As a first degradation product typically oligomeric hydroxyalkanoates (OHAs) and/or monomers are released and many of the latter ones have also been found as natural metabolites in animals (Butt et al., 2018;Zhang et al., 2018;Utsunomia et al., 2020). However, little is known about the interaction of OHAs with the human body and whether the different polymer sizes may result in a particular toxicity.
Thus, in the case of 3HB, it is known that it is a metabolic intermediate produced in the liver from the oxidation of fatty acids released from adipose tissue (Persson, 1970) and is readily consumed in normal mammalian metabolism being an essential carrier of energy and having some cellular signaling functions (Newman and Verdin, 2017). Similarly, 4HB is a natural human metabolite of γ-aminobutyric acid being a neurotransmitter or neuromodulator in the brain. It has a short half-lifetime of about 35 min in the human organism, being rapidly eliminated as exhaled CO 2 after its metabolic oxidation (Sendelbeck and Girdis, 1985;Maitre, 1997). The presence of 4HB in the human organism has been observed experimentally, trace amounts of 4HB have been found in blood, urine and tissues from healthy individuals that were not exposed to 4HB precursors and/or P4HB (Busardò and Jones, 2015). The chemical 4HB is known

Oligomer
Observed biological effects

3HB
The oral administration of the 3-hydroxybutyric acid salts increased ketonemia in Wistar rats and conducted to the improvement of body control by the reduction of fat mass and amelioration of serum lipid profile suggesting that ketone supplements could be helpful in the treatment of obesity and metabolic diseases (Caminhotto et al., 2017). Some gastrointestinal disturbances were observed in vivo in individuals who consumed a high dose of 3HB (Clarke et al., 2012). A high quantity of 3HB could be harmful because it may conduct to acidosis (Newman and Verdin, 2017). The 3HB may produce eye, skin and respiratory irritations (National Library of Medicine, 2020c).

4HB
The 4HB does not bind to a significant extent to plasma proteins. High doses of 4HB produced hypotension, bradycardia, tachycardia, hypothermia, unconsciousness, acute respiratory acidosis and gastrointestinal disturbances. After ingestion, 4HB is easily absorbed and it is able to cross the blood-brain barrier. It is rapidly metabolized and excreted through the lungs (Busardò and Jones, 2015). The 4HB may produce eye damage, drowsiness or dizziness (National Library of Medicine, 2020a).

3HV
The 3HV may produce eye, skin and respiratory irritations (National Library of Medicine, 2020b).

4HV
The 4HV is marketed as a dietary supplement and replacement for 4HB as at higher doses it shares some effects with 4HB such as sedation, catalepsy, and ataxia (Carter et al., 2005).

O3HB
Low molecular weight O3HB (smaller than 20 monomeric units) are constituents of cells and are covalently attached to proteins located within membranes and organelles (Reusch, 2013).
The in vitro cytotoxicity of the OHAs against mouse fibroblasts decreased with increasing OHAs side chain length (Sun et al., 2007 as a neurotransmitter and a psychoactive drug, a reason for why it is a controlled substance in the United States since 2000 and in the EU since 2001 (Riviello, 2010). Taking into account the lack of data concerning the biological effects of the PHA degradation products in general, this study aimed to use a computational approach to predict the Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) profiles and possible toxicological endpoints of the low molecular weight O3HB, O4HB, O3HV, O4HV, and their co-oligomers.
The Simplified Molecular-Input Line-Entry System (SMILES) formulas and structural mol files of the investigated oligomers have been obtained using ACD/ChemSketch (ACD ChemSketch 12.00 and Chemicalize from ChemAxon) utility (https:// chemicalize.com accessed in April 2020) and were used further to compute the molecular weights and partition coefficients of these oligomers and to predict their pharmacokinetics profiles as well as their biological effects. The SMILES and two-dimensional formulas of the monomeric units of O3HB, O3HV, O4HB, and O4HV are presented in Figure 1.
Specific literature is abundant and partially integrated in computational tools for predicting the pharmacokinetic profiles and the biological effects of chemicals. We have selected for this study computational tools that are freely accessible online or are open-source, are developed by recognized groups and benefit of a high number of citations, are based on models with a larger amount of training data, are continuously updated, are robust and their overall accuracy of predictions is usually higher than 70%. These methods have been successfully used for predicting TABLE 2 | In silico tools considered for assessing the ADMET profiles of investigated oligomers.

Computational tool
Predicted activity and accuracy of prediction admetSAR2.0 is a free available structure-activity relationship database that contains over 210,000 available properties data curated from literature for about 96,000 chemicals. It includes 22 qualitative classification models and 5 quantitative regression models allowing estimations quantitatively described by a probability output (Cheng et al., 2012;Yang et al., 2018).
As input data, this tool uses the SMILES formulas of investigated oligomers and we have extracted the probabilities of ADMET properties. ENDOCRINE DISRUPTOME tool uses the molecular docking approach to predict the interactions between the investigated chemicals and 12 human nuclear receptors (Kolsek et al., 2014).
As input data we have used the SMILES formulas of investigated oligomers and extracted data considering the color coded table that was displayed when the molecular docking was finished: green for compounds reflecting a low probability of interacting with the nuclear receptors, orange for compounds revealing a mean probability of interacting with the nuclear receptors and red for compounds revealing a high probability of affecting the nuclear receptors.
Pred-Skin 3.0 is a web-server allowing predictions concerning skin sensitization potential of chemicals based on QSAR models (Braga et al., 2017;Alves et al., 2018a). SMILES formulas have been used as input data and the probability of oligomers to illustrate skin sensitization or non-sensitization potentials, respectively, the probability maps illustrating the fragments contributions toward skin sensitization potential were retrieved.
It performs predictions based on five sensitization assay: in vivo (murine local lymph node assay, LLNA, accuracy 0.70-0.84), in chemico (Direct Peptide Reactivity Assay, DPRA, accuracy 0.73-0.76), in vitro (KeratinoSens and human Cell Line Activation Test, H-CLAT, accuracy 0.80-0.86), human repeated insult patch, HRIPT, test and human maximization test, HMT, accuracy 0.70-0.84. There also is a Bayesian consensus model that is generated by averaging the predictions of individual models. Pred-Skin 3.0 also outcomes a probability maps allowing the visualization of the contribution of predicted fragment toward skin sensitization.
Pred-hERG 4.2 is a web tool that builds predictive models of the ability of a chemical compound to inhibit the human Ether-à-go-go Related Gene (hERG) based on the QSAR approach (Braga et al., 2015). SMILES formulas have been used as input data and the probabilities of oligomers to illustrate hERG K + channel blocking or non-blocking potential, together with the probability maps illustrating the fragments contributions toward hERG blockage potential were retrieved.
The outcome is a binary prediction of hERG non-blocker or blocker (0.80). This tool also delivers a probability maps allowing the visualization of the contribution of predicted fragment toward hERG blockage.
CarcinoPred-EL (Carcinogenicity Prediction using Ensemble Learning methods) is a computational tool used for predictions concerning the carcinogenicity of chemicals ensemble classification models (Zhang et al., 2017). SMILES formulas have been used as input data and the software provides a binary prediction (Yes or No) for possible carcinogenicity of the investigated compounds.
Toxtree is an open-source application that performs predictions concerning carcinogenicity and mutagenicity by applying a decision tree approach (Patlewicz et al., 2008).
As input data we have used the SMILES formulas of investigated oligomers and we have retrieved the predictions (Yes or No) for the carcinogenic and mutagenic potential of investigated compounds.
The carcinogenic and mutagenic potential is predicted (0.70).
PASS (Prediction of Activity Spectra of Substances) is a computational tool that predicts biological activity spectra and toxic/side effects starting to the structural formulae of chemical compounds and using the QSAR approach (Poroikov et al., 2007). This tool also used the SMILES formulas as entry data and provides the probability that the investigated compound to be active for a given adverse biological activity.
PASS has been used to predict toxic and adverse effects (0.95).
Numbers in parantheses represent probabilities of predictions.
For every method, a short description, the predicted biological activity, the accuracy of predictions and references are provided. All these in silico tools use as input the SMILES formulas of investigated chemicals. The limits of these tools are mainly expressed by the fact that they do not allow to take into account the dose of the investigated compound. We have used these tools to obtain predictions concerning the ADMET profiles and possible toxicological endpoints of the considered oligomers. We must mention that there is no consideration concerning the chirality of the investigated compounds.
Usually, these predictions are expressed as probabilities for a given compound to have (positive values) or not to have (negative values) a tested biological activity. We have considered that a predicted biological effect has a chance to be detected experimentally if the computed probability is higher than 0.7. This threshold has been considered because in this case the investigated compound reveals a high similarity to compounds with known biological activity and used for building the statistical model (Filimonov et al., 2014). The outcomes of the used computational tools are compared with each other, the correlation/un-correlations are emphasized, and all the results are discussed taking into account the information available in scientific publications and the limitations of the in silico tools.

RESULTS
The molecular weights and the consensus values of partition coefficients for the studied oligomers, computed using SwissADME tool, are given in Table 3.
The oligomers considered in the present study have a molecular weight between 104.1 and 3221.72 g/mol (see Table 3). The values of partition coefficients are slightly distinct for oligomers containing the same number of monomers such as 3hydroxybutyrate and 4-hydroxybutyrate, and 3-hydroxyvalerate and 4-hydroxyvalerate, respectively. This is also true for co-oligomers of 3-hydroxybutyrate-co-3-hydroxyvalerate and 4-hydroxybutyrate-co-4-hydroxyvalerate containing the same number of monomeric units but with different sequence patterns. As the partition coefficient is an important physicochemical property for pharmacokinetics of xenobiotics, it underlines the significance of the pattern in the case of co-oligomers for their biological interactions.
Concerning the absorption and distribution profiles, admetSAR2.0 tool calculates for every investigated oligomer the probabilities to reveal gastrointestinal absorption (GI), to be able to penetrate the blood-brain barrier (BBB), to be a substrate or an inhibitor of P-glycoprotein (PgpS/PgpI) (Figure 2), and to be able to bind to plasma proteins (PPB) (Figure 3). The predicted probabilities take values between 0 and 1 when the investigated biological activity is considered to be present and between −1 and 0 when the biological activity is considered being absent. A probability value closer to 1 or −1 indicates a biological effect that is highly probable or highly improbable, respectively. Figure 2 suggests that there are high probabilities of gastrointestinal absorption of all OHAs, the capability of all investigated OHAs to penetrate the blood-brain barrier and that these compounds should be considered neither substrates nor inhibitors for the P-glycoprotein (P-gp). Figure 3 advocates a reduced potential of the investigated OHAs to bind to plasma proteins (PPB) as all values of the probabilities are lower than 0.7. This potential moderately increases with the molecular weight. Co-oligomers illustrate higher probabilities to bind to plasma proteins, and co-oligomers containing the same number of monomers but with different sequence patterns reveal distinct probabilities to bind to plasma proteins. It emphasizes the important role of the molecular weight and oligomer's pattern on the biological activity of investigated OHAs.
Regarding the effects of the investigated OHAs on the metabolism of endogenous and exogenous compounds, admetSAR2.0 tool suggests that none of the investigated FIGURE 2 | Absorption profiles of low molecular weight oligo-hydroxyalkanoates calculated with admetSAR2.0: u denotes the number of units in the oligomer, O3HB denotes the oligomer of 3HB, O3HV denotes the oligomer of 3HV, O4HB denotes the oligomer of 4HB, O4HV denotes the oligomer of 4HV. In the case of co-oligomers, BV, VB, BVB, VBV, BVV, VBB, BVBV, and VBVB, respectively, illustrate the succession of the 3-hydroxybutyrate/4-hydroxybutyrate (B) and 3-hydroxyvalerate/4-hydroxyvalerate (V) monomers in the co-oligomeric chain, respectively. The predicted probabilities for gastrointestinal absorption (GI), for penetration of the blood-brain barrier (BBB), or to be a substrate or an inhibitor pf P-glycoprotein (PgpS/PgpI) take values between 0 and 1 in the case of a biological activity that is present and between −1 and 0 when the activity is considered absent. For values closer to 1, the biological effects are highly probable and values closer to −1 correspond to highly improbable biological effects. oligomers has to be considered as a substrate for the human cytochromes involved in the metabolism of xenobiotics. Also, they are not able to inhibit these cytochromes (Supplementary Table 1).
Predictions obtained using admetSAR2.0 concerning toxicological endpoints of investigated OHAs suggest that they may not be considered neither as carcinogens not mutagens are not potential inhibitors of the hERG K + channel, and do not induce hepatotoxicity. However, they illustrate a moderate potential of producing eye corrosion and oligomers having low molecular weights may also produce eye irritations (Supplementary Table 2). We must mention that none of the hypothetical co-oligomers containing 4HB and 4HV that were part of this study have been predicted to cause eye irritation. The ability of the monomers 3HB, 3HV, 4HB, and 4HV to trigger eye, skin, and respiratory irritations has also been forecast using the PASS tool (see below) and few literature Abbreviations are used as explained in Figure 2.
Predictions obtained by using Pred-hERG 4.2 are illustrated in Figure 4 for the hypothetical co-oligomers O4HBVV and O4HVBB. This Figure 4 emphasizes the probability map illustrating the contribution of fragments (green regions) of the co-oligomers to produce hERG K + channel blockage.
The co-oligomers O4HBVV and O4HVBB contain fragments contributing to both increasing and decreasing hERG K + channel blockage. Synthesis of all results furnished by Pred-hERG 4.2 tool is presented in Table 5, suggesting that these compounds do not own a hERG blockage potential. Table 5 also contains the synthesis of the results regarding the skin sensitization potential (see below).
Pred-Skin 3.0 computational tool has been used to obtain prediction concerning the skin sensitization potential of the investigated oligomers. Figure 5 illustrates the probability maps reflecting the contribution of fragments (pink regions) of the hypothetical co-oligomers O4HBVV and O4HVBB on the skin sensitization potential based on the KeratinoSens model.  5 reveals that O4HBVV contains fragments that may contribute to decreasing the skin sensitization potential, but O4HVBB contains fragments suggesting a skin sensitization potential. This is in good agreement with the consensus model, as presented in Table 5. Information furnished by Pred-Skin3.0 online tool by using the five models based on skin sensitization assays is presented in Supplementary Table 3.
Predictions obtained using CarcinoPred-El and Toxtree (Supplementary Table 4) are in good correlation with each other and with those obtained using admetSAR2.0, advocating Non-cardiotoxic (0.6) Non-sensitizer The confidence of predictions obtained using Pred-hERG4.2 is also illustrated.
Abbreviations are used as explained in Figure 2.
that investigated OHAs should not be considered as reflecting carcinogenicity nor mutagenicity. ENDOCRINE DISRUPTOME computational tool suggested that investigated OHAs have only a weak potential to produce an antagonistic effect on the androgen receptor (AR) and oligomers containing 6-8 monomeric units may also weakly affect the glucocorticoid receptor (GR) (Supplementary Table 5). We also underline that oligomers with higher molecular weight (more than 8 monomeric units) were too big to accommodate in the binding cavities of considered nuclear receptors and therefore calculations have been aborted.
The predicted side effects of OHAs were obtained using the PASS tool and are listed in Table 6.
There are some side effects predicted for all investigated OHAs: toxicity by respiration, eye irritation, skin irritation, and conjunctivitis.

DISCUSSION
The high probabilities of gastrointestinal absorption of OHAs that are illustrated in Figure 2 correspond with known data indicating the rapid absorption of hydroxyalkanoates (Busardò and Jones, 2015;Caminhotto et al., 2017) and of 3HB (Cuenoud FIGURE 5 | Probability maps illustrating the contribution of fragments of co-oligomers O4HBVV (a) and O4HVBB (b) to skin sensitization potential based on KeratinoSens model. Fragments in green illustrate an increase in skin sensitization potential, the pink fragments contribute to decrease of skin sensitization. et al., 2020), the blood-brain permeability for ketone bodies, and the use of 4HB as neurotransmitter and psychoactive drug (National Library of Medicine, 2020a). Figure 2 also suggests that the analyzed compounds are neither substrates nor inhibitors for the P-glycoprotein (P-gp). This protein is expressed in various normal tissues and has a protective role against xenobiotic substances including toxins (Cascorbi, 2011). Our results suggest that the presence of OHAs in the human organism does not alter pharmacokinetics of drugs by being substrates and/or by inhibiting the P-gp.
There are some interesting differences between the predicted pharmacokinetics of oligomers and co-oligomers. Thus, gastrointestinal absorption and the potential of blood-brain barrier penetration are predicted to be higher for the cooligomers. Moreover, there are small differences between the The numbers in parentheses represent the probabilities that the mentioned side effects will occur. The predicted probabilities take values between 0 and 1. The values closer to 1 illustrate a high probability that these oligomers produce the mentioned side effects. Abbreviations are used as explained in Figure 2.
predicted probabilities of GI, BBB, PPB, and PgpS/PgpI for co-oligomers containing the same number of monomeric units but having distinct sequence patterns. The ability of OHAs to bind to plasma proteins has been predicted to be low. However, it was found that the probability to bind to plasma proteins gradually increases with the molecular weight. The co-oligomers containing the sequences 3-hydroxybutyrate-3-hydroxyvalerate or 4-hydroxybutyrate-4hydroxyvalerate revealed to have the highest probabilities to bind to plasma proteins. Low probabilities to bind to plasma proteins have been also predicted for low molecular weight chitosan based oligomers (Roman et al., 2019) and low molecular weight oligomers of lactic acid (Dascalu et al., 2020). This outcome is also in good agreement with the information that the binding of 4HB to plasma proteins is negligible and its distribution on the human organism depends on its distribution in the total body water (Busardò and Jones, 2015).
The human OATPs and OCTs handle a variety of endogenous and xenobiotic substrates playing a significant role in the absorption, distribution, and elimination of chemicals. OATP1B1, OATP1B3, OATP2B1, MATE1, and OCT2 are considered when investigating the interactions of chemicals with transporter polypeptides. This is of great importance because diverse classes of drugs are known as substrates of these transporters and knowledge about inhibition of these transporters is of great value in revealing actions in complex systems like hepatocytes (Karlgren et al., 2012;Belzer et al., 2013;Wittwer et al., 2013). Furthermore, inhibition of any of these transporters may result in potential drug-xenobiotics interactions, alterations in chemicals exposure and their accumulation in various tissues. The existence of high-quality data concerning the inhibition of these transporters allowed statistical modeling and, consequently, they are considered by numerous computational tools. OATP1B 1 and OATP1B3 are mainly expressed in the human liver (König et al., 2000) and are responsible for the uptake of bile acids and a broad range of other organic anions being involved in the hepatic clearance (Geier et al., 2007;Alam et al., 2018). OATP2B1 is a transporter with a narrow substrate specificity and is expressed in various human tissues including the liver, kidney, lung, and small intestine (Karlgren et al., 2012). MATE1 is widely expressed in the human body tissues and contributes to the tissue distribution and excretion of many drugs (Wittwer et al., 2013). OCT2 facilitates the renal secretion of organic cations being considered as a potential target for undesirable drug-drug interactions and numerous studies envisaged the understanding of the basis of its selectivity (Belzer et al., 2013). Data presented in Table 4 suggest that all investigated OHAs are potential inhibitors of OATP1B1 and OATP1B3. Decreased transport function of OATP1B1 and OATP1B3 may lead to drug-drug interactions and severe adverse events. Consequently, the US Food and Drug Administration recommended evaluating the inhibitory potential of new chemicals against OATP. It underlines the importance of predictions concerning the inhibition of OATPs by OHAs.
Data presented in Table 5 suggest the non-blocking potential of every of the investigated oligomers toward hERG K + channel and are in very good agreement with predictions obtained using admetSAR2.0 computational facility. Table 5 also advocates that, excepting the co-oligomers O3HBVBV, O4HVB, and O4HVBB, the other investigated oligomers do not show skin sensitization potential. This result underlines that co-oligomers have biological activities that are different compared with the ones found for homo-oligomers and that the biological effects also depend on the composition of the co-polymer with respect to monomeric units. This result is in good correlation with literature data presented by Bonartsev et al. (2019) and revealing that cooligomers have distinct properties and actions than homooligomers (Bonartsev et al., 2019).
Predictions presented in Table 6 are in good agreement with specific literature data (see Table 1). The side effects predicted for all investigated OHAs, i.e., toxic by respiration, eye irritation, skin irritation, may be important especially for professional exposure in case of powders. Furthermore, the results of this study highlight that co-oligomers show quite distinct biological effects compared with those of homo-oligomers and that the biological effects usually depend on the positions of HB and HV monomers in the chain.
Quite similar results with those emphasized in this study were obtained when assessing the pharmacological profiles of the short oligomers resulting from the degradation of other biopolymers that are used for medical applications: oligomers (up to 40 monomeric units) of lactic acid (OLAs) (Dascalu et al., 2020) and low molecular weight chito-oligosaccharides (Roman et al., 2019). OLAs were also predicted as being able to inhibit the organic anion transporting peptides OATP1B1 and/or OATP1B3.
In addition, they showed a minor probability of affecting the androgen and glucocorticoid receptors but their powders may produce eye injuries. Chito-oligosaccharides containing up to 8 monomeric units (glucosamine and/or N-acetyl glucosamine) were also suggested as potential inhibitors of the organic anion transporting polypeptides 1B1 and 1B3. They illustrated a minor probability of affecting the androgen receptor and they are considered to be able to produce cardiotoxicity and phospholipidosis.
The results of this study, and those of similar studies published in scientific literature, emphasize the importance of assessment of the physiological effects of PHAs and of their degradation products in the human organism and, in general, the necessity of designing novel polymers for medical applications with tailormade properties (e.g., P4HB4HV). The limitations of our study are the disregard of the chirality of the monomeric units of OHAs and of the concentration of investigated compounds when the computational predictions of ADMET profiles were realized.

CONCLUSION
Low molecular weight oligo-hydroxyalkanoates may result from the degradation of polyhydroxyalkanoate based implants or medical devices (e.g., drug delivery systems). Due to their presence in the human organism, it is necessary to investigate their pharmacokinetics profiles and to identify possible toxicological effects at a very early stage of research. The present study suggests favorable pharmacokinetic profiles for all types of OHAs. The investigated oligomers are not considered to reflect any hepatotoxicity, cardiotoxicity, mutagenicity nor carcinogenicity. There are a few predicted side effects for the investigated OHAs: the potential to produce eye, skin and respiratory irritation (in the case of exposure to their powders), antagonistic effects against the androgen receptor and inhibition of the OATP1B1 and OATP1B3. Few co-oligomers (O3HBVBV, O4HVB, O4HVBB) revealed high probabilities to produce a skin sensitization potential. The effects of eye, skin and respiratory irritations, and skin sensitization potential must be taken into consideration especially for people who are professionally exposed (e.g., during production or packaging of these polymers).
In our study, we could not find any significant differences between the predicted effects for distinct types of OHAs: O3HB, O4HB, O3HV, O4HV. There are dissimilarities in the predicted probabilities reflecting pharmacokinetics properties and side effects for oligomers by comparison to co-oligomers and for cooligomers containing the same number of monomeric units but having distinct sequence patterns.
There is a good agreement between the predictions obtained using various computational tools and concerning the biological/toxicological effects that have been investigated. Furthermore, there is a good correlation between the obtained predictions and literature experimental data, when available (which is not the case for the hypothetical P4HB4HV).
The results presented in this study suggest that the therapeutic use of PHAs offers new opportunities at moderate risks. Moreover, these results could be used to define or to guide in vitro and in vivo toxicity tests such as to contribute to the selection of best prototypes for medical applications. Implementation of defined co-polymeric structures, such as the co-polymer consisting of 4HB/4HV, could conduct to the further optimization of the properties of the biomaterial.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
DR performed the computations using Pred-Skin3.0, Pred-hERG4.2, CarcinoPredEL, Toxtree, and Endocrine Disruptome, and contributed to results analysis and manuscript editing. MF performed computations using admetSAR2.0 and PASS online and contributed to results analysis. VO furnished the SMILES and structural formulas of oligomers and contributed to the conception and design of the manuscript. AI and MZ conceived and planned the study and contributed to the analysis and conception of the results, design and editing of the manuscript. All authors contributed to the article and approved the submitted version.