Design of Two New Sulfur Derivatives of Perezone: In Silico Study Simulation Targeting PARP-1 and In Vitro Study Validation Using Cancer Cell Lines

Poly-ADP-Ribose Polymerase (PARP-1) is an overexpressed enzyme in several carcinomas; consequently, the design of PARP-1 inhibitors has acquired special attention. Hence, in the present study, three compounds (8–10) were produced through a Michael addition protocol, using phenylmethanethiol, 5-fluoro-2-mercaptobenzyl alcohol, and 4-mercaptophenylacetic acid, respectively, as nucleophiles and perezone as the substrate, expecting them to be convenient candidates that inhibit PARP-1. It is convenient to note that in the first stage of the whole study, the molecular dynamics (MD) simulations and the quantum chemistry studies of four secondary metabolites, i.e., perezone (1), perezone angelate (2), hydroxyperezone (3), and hydroxyperezone monoangelate (4), were performed, to investigate their interactions in the active site of PARP-1. Complementarily, a docking study of a set of eleven sulfur derivatives of perezone (5–15) was projected to explore novel compounds, with remarkable affinity to PARP-1. The molecules 8–10 provided the most adequate results; therefore, they were evaluated in vitro to determine their activity towards PARP-1, with 9 having the best IC50 (0.317 µM) value. Additionally, theoretical calculations were carried out using the density functional theory (DFT) with the hybrid method B3LYP with a set of base functions 6-311++G(d,p), and the reactivity properties were compared between the natural derivatives of perezone and the three synthesized compounds, and the obtained results exhibited that 9 has the best properties to bind with PARP-1. Finally, it is important to mention that 9 displays significant inhibitory activity against MDA-MB-231 and MCF-7 cells, i.e., 145.01 and 83.17 µM, respectively.


Introduction
Cancer is the principal cause of death worldwide in countries of all income levels [1].Cells turn cancerous after the accumulation of several mutations in their genes that control cell division and survival [2].A wide variety of proteins have been implicated in the establishment and progression of cancer.Recently, particular attention has been paid to Poly-ADP-Ribose Polymerase 1 (PARP-1), due to its implication in the development and progress of cancer [3].PARP-1 is a nuclear enzyme, which interacts with DNA and regulate several cell functions, such as DNA repair, recombination, proliferation, and genomic stability [4].It has been well documented that PARP-1 is overexpressed in various carcinomas.Hence, the development of PARP-1 inhibitors has acquired importance since it has been demonstrated that this kind of molecules could induce cell death in cancer cells [5].PARP-1 is an enzyme that facilitates the post-translational modification of proteins in response to DNA damage and genotoxic stress [6].PARP-1 is formed by three domains: (1) a DNA-binding domain, (2) an auto-modification domain (AMD) with a role in the interaction between PARP-1 and its partner proteins, and (3) a catalytic domain [7].Clinical PARP-1 inhibitors bind to a catalytic pocket, where they directly interfere with ADPribosylation [8].The active site of PARP-1 is formed by a sequence of 50 conserved amino acids, and this region is identical for all PARP-1 sequences of vertebrates [6].The catalytic pocket is formed by a conserved catalytic triad consisting of Hys201, Tyr235, and Glu327, which are required for the catalytic activity of PARP-1 [4].Thus, classically, PARP-1 inhibitors are designed to target these amino acid residues.
Perezone represents a phyto compound occurring in the roots of plant specimens of the genus Acourtia [9].It has been demonstrated that perezone (1) and natural compounds chemically related to it, such as perezone angelate (2), hydroxyperezone (3), and hydroxyperezone monoangelate (4) (Figure 1), induce apoptosis in several cancer cell lines in a range of micromolar concentrations [10][11][12][13][14]. Through in vitro and in silico studies, it has also been demonstrated that two mechanisms explain the pro-apoptotic activity of perezone: PARP-1 inhibition and the induction of an oxidative stress state [12].To develop novel PARP-1 inhibitors, in the present study, molecular dynamics (MD) and quantum chemistry studies of 1-4 in the active site of PARP-1 were performed to investigate the interactions that drive the protein-ligand recognition, and this knowledge is crucial to design novel PARP-1 inhibitors.Then, a set of perezone derivatives were designed, highlighting that only three compounds, with the highest affinity to PARP-1 based on docking studies, were selected.Finally, the PARP-1 inhibition and its effects on several cancer cell lines were evaluated for selected compounds.
Int. J. Mol.Sci.2024, 25, x FOR PEER REVIEW 2 of 21 establishment and progression of cancer.Recently, particular attention has been paid to Poly-ADP-Ribose Polymerase 1 (PARP-1), due to its implication in the development and progress of cancer [3].PARP-1 is a nuclear enzyme, which interacts with DNA and regulate several cell functions, such as DNA repair, recombination, proliferation, and genomic stability [4].It has been well documented that PARP-1 is overexpressed in various carcinomas.Hence, the development of PARP-1 inhibitors has acquired importance since it has been demonstrated that this kind of molecules could induce cell death in cancer cells [5].PARP-1 is an enzyme that facilitates the post-translational modification of proteins in response to DNA damage and genotoxic stress [6].PARP-1 is formed by three domains: (1) a DNA-binding domain, (2) an auto-modification domain (AMD) with a role in the interaction between PARP-1 and its partner proteins, and (3) a catalytic domain [7].Clinical PARP-1 inhibitors bind to a catalytic pocket, where they directly interfere with ADP-ribosylation [8].The active site of PARP-1 is formed by a sequence of 50 conserved amino acids, and this region is identical for all PARP-1 sequences of vertebrates [6].The catalytic pocket is formed by a conserved catalytic triad consisting of Hys201, Tyr235, and Glu327, which are required for the catalytic activity of PARP-1 [4].Thus, classically, PARP-1 inhibitors are designed to target these amino acid residues.Perezone represents a phyto compound occurring in the roots of plant specimens of the genus Acourtia [9].It has been demonstrated that perezone (1) and natural compounds chemically related to it, such as perezone angelate (2), hydroxyperezone (3), and hydroxyperezone monoangelate (4) (Figure 1), induce apoptosis in several cancer cell lines in a range of micromolar concentrations [10][11][12][13][14]. Through in vitro and in silico studies, it has also been demonstrated that two mechanisms explain the pro-apoptotic activity of perezone: PARP-1 inhibition and the induction of an oxidative stress state [12].To develop novel PARP-1 inhibitors, in the present study, molecular dynamics (MD) and quantum chemistry studies of 1-4 in the active site of PARP-1 were performed to investigate the interactions that drive the protein-ligand recognition, and this knowledge is crucial to design novel PARP-1 inhibitors.Then, a set of perezone derivatives were designed, highlighting that only three compounds, with the highest affinity to PARP-1 based on docking studies, were selected.Finally, the PARP-1 inhibition and its effects on several cancer cell lines were evaluated for selected compounds.

MD Simulations
MD allows the monitoring of binding interactions for a specific ligand with a particular protein during a defined time.As it has been previously demonstrated [13,14], 1-4 display IC 50 values of 181.5 µM, 5.25 µM, 320.3 µM, and 28.7 µM, respectively.The analysis of IC 50 of selected compounds allows the observation that the hydroxyl substituent in 3 diminishes the potency as PARP-1 inhibitor in comparison with 1.However, the addition of an α,β-unsaturated ester moiety as an angelate, 2 and 4, increases the PARP-1 inhibitory activity in comparison with 1, highlighting the higher PARP-1 inhibitory activity of 2. Consequently, the four protein-ligand complexes were subjected to 50 ns of MD simulations to investigate their interactions in the active site of PARP-1 at the molecular level.
Figure 2a shows the RMSD values of PARP-1 in the absence and the presence of compounds 1-4.The determination of the RMSD of the protein alone and its complex with the studied compounds allows the quantitative evaluation of the changes in the protein structure, throughout the simulation time [15], thus determining if the simulation reaches the convergence [16].As can be seen in Figure 2a (green), the RMSD of PARP-1 fluctuated from 2.9 to 4.1 Å, with a mean value of 3.3 Å, reaching an equilibrium at 10 ns.In comparison, the RMSD of the protein in complex with 1 fluctuated from 2.7 to 4.1 Å with a mean value of 3.1 Å, reaching an equilibrium at 10 ns (Figure 2a, pink), while the RMSD of the protein in complex with 3 ranged from 2.7 to 3.5 Å, with a mean value of 2.9 Å, reaching an equilibrium at 8 ns (Figure 2a, red).In contrast, the RMSD of the protein in complex with 2 fluctuated from 2.9 to 3.7 Å, without reaching an equilibrium (Figure 2a, gray), while the RMSD of the protein in complex with 4 ranged from 3.2 to 5.0 Å, without reaching an equilibrium (Figure 2a, blue).RMSD values showed that the presence of compounds with the highest PARP-1 inhibitory activity (2 and 4) prevent PARP-1 from reaching the equilibrium, and this is evidenced by higher RMSD values.

MD Simulations
MD allows the monitoring of binding interactions for a specific ligand with a particular protein during a defined time.As it has been previously demonstrated [13,14], 1-4 display IC50 values of 181.5 μM, 5.25 μM, 320.3 μM, and 28.7 μM, respectively.The analysis of IC50 of selected compounds allows the observation that the hydroxyl substituent in 3 diminishes the potency as PARP-1 inhibitor in comparison with 1.However, the addition of an α,β-unsaturated ester moiety as an angelate, 2 and 4, increases the PARP-1 inhibitory activity in comparison with 1, highlighting the higher PARP-1 inhibitory activity of 2. Consequently, the four protein-ligand complexes were subjected to 50 ns of MD simulations to investigate their interactions in the active site of PARP-1 at the molecular level.
Figure 2a shows the RMSD values of PARP-1 in the absence and the presence of compounds 1-4.The determination of the RMSD of the protein alone and its complex with the studied compounds allows the quantitative evaluation of the changes in the protein structure, throughout the simulation time [15], thus determining if the simulation reaches the convergence [16].As can be seen in Figure 2a (green), the RMSD of PARP-1 fluctuated from 2.9 to 4.1 Å, with a mean value of 3.3 Å, reaching an equilibrium at 10 ns.In comparison, the RMSD of the protein in complex with 1 fluctuated from 2.7 to 4.1 Å with a mean value of 3.1 Å, reaching an equilibrium at 10 ns (Figure 2a, pink), while the RMSD of the protein in complex with 3 ranged from 2.7 to 3.5 Å, with a mean value of 2.9 Å, reaching an equilibrium at 8 ns (Figure 2a, red).In contrast, the RMSD of the protein in complex with 2 fluctuated from 2.9 to 3.7 Å, without reaching an equilibrium (Figure 2a, gray), while the RMSD of the protein in complex with 4 ranged from 3.2 to 5.0 Å, without reaching an equilibrium (Figure 2a, blue).RMSD values showed that the presence of compounds with the highest PARP-1 inhibitory activity (2 and 4) prevent PARP-1 from reaching the equilibrium, and this is evidenced by higher RMSD values.Based on the RMSF results, it is possible to demonstrate that the studied compounds induce great conformational changes in PARP-1, mainly at the N-terminal region (Figure Based on the RMSF results, it is possible to demonstrate that the studied compounds induce great conformational changes in PARP-1, mainly at the N-terminal region (Figure 2b), and this is evidenced by higher RMSF values.Interestingly, residues from the catalytic triad (Hys201, Tyr235, and Glu327) show low RMSF values (Figure 2b, red arrows).
Figure 3 shows the binding mode of the studied compounds to the catalytic site of PARP-1, showing that larger compounds 2 and 4 fit deep in the catalytic site in contrast to compounds 1 and 3. Additionally, the binding mode of 2 is maintained throughout the 50 ns of the simulation, thus explaining its good PARP-1 inhibitory activity.
Figure 3 shows the binding mode of the studied compounds to the catalytic site of PARP-1, showing that larger compounds 2 and 4 fit deep in the catalytic site in contrast to compounds 1 and 3. Additionally, the binding mode of 2 is maintained throughout the 50 ns of the simulation, thus explaining its good PARP-1 inhibitory activity.As shown in Table S1 (Supplementary Material), the molecule 2 interacts with Hys201 and Tyr235 located at the catalytic site, and interestingly, these interactions are maintained during the entire simulation, supporting the higher PARP-1 inhibitory activity that was previously mentioned.
The analysis of the 50 ns snapshot of MD simulations shows that 1 exhibits hydrophobic interactions with the following amino acid residues: Ala167, Asn347, Phe348, and Asn347.In addition, both a hydrogen bond interaction with Asn347 (2.7 Å) and one π-π type interaction, between the quinonic ring of perezone and the amino acid residue ring of Tyr168, are observed.Compound 2 displays four interactions of hydrogen bonds with the following amino acid residues: Asn207 (1.7 Å), Gly202 (2.6 Å), Gly233 (2.3 Å), and Arg217 (2.2 Å).Also, it presents hydrophobic interactions with the amino acid residues Ser203, Ile234, Tyr246, Gly233, Tyr235, Hys201, and Asn106.Compound 3 has two interactions of hydrogen bonds, with the following amino acid residues: Gly227 (2.4 Å) and Asp105 (1.8 Å); in addition, the quinonic ring of hydroxyperezone displays π-π type As shown in Table S1 (Supplementary Material), the molecule 2 interacts with Hys201 and Tyr235 located at the catalytic site, and interestingly, these interactions are maintained during the entire simulation, supporting the higher PARP-1 inhibitory activity that was previously mentioned.

Quantum Chemistry 2.2.1. Molecular Electrostatic Potential (MPE)
The MEP shows the electron density of a molecule, and it is used to identify the sites of positive and negative electrostatic potentials for nucleophilic and electrophilic attacks.It is also a powerful tool to identify biological recognition processes and hydrogen bonding interactions [17].Figure 4 shows the MEP of the molecules 1-4, and red areas represent the regions of high electron density, and this is principally visualized in both carbonyl groups of the quinonic ring of all molecules, and for 3 and 4, the density is also displayed in the carbonyl of the ester moiety.The blue areas represent the regions of low electron density; this can be seen in the hydrogen atoms of the hydroxyl groups.MEP is indicative of sites, both negative (red) or positive (blue), on molecules (ligands) where they could favor the formation of hydrophobic interactions, π-π type interactions, or hydrogen bonds towards a protein.Thus, due to a high electron density and a higher MEP value of 2, these features contributed to the development of interactions with amino acids, highlighting the four hydrogen bonds described above and resulting in good PARP-1 inhibitory activity.In contrast, 1, which shows a low electron density zone, displays poor interactions with amino acids, highlighting the four hydrophobic interactions and diminishing the PARP-1 inhibitory activity.

Molecular Electrostatic Potential (MPE)
The MEP shows the electron density of a molecule, and it is used to identify the sites of positive and negative electrostatic potentials for nucleophilic and electrophilic attacks.It is also a powerful tool to identify biological recognition processes and hydrogen bonding interactions [17].Figure 4 shows the MEP of the molecules 1-4, and red areas represent the regions of high electron density, and this is principally visualized in both carbonyl groups of the quinonic ring of all molecules, and for 3 and 4, the density is also displayed in the carbonyl of the ester moiety.The blue areas represent the regions of low electron density; this can be seen in the hydrogen atoms of the hydroxyl groups.MEP is indicative of sites, both negative (red) or positive (blue), on molecules (ligands) where they could favor the formation of hydrophobic interactions, π-π type interactions, or hydrogen bonds towards a protein.Thus, due to a high electron density and a higher MEP value of 2, these features contributed to the development of interactions with amino acids, highlighting the four hydrogen bonds described above and resulting in good PARP-1 inhibitory activity.In contrast, 1, which shows a low electron density zone, displays poor interactions with amino acids, highlighting the four hydrophobic interactions and diminishing the PARP-1 inhibitory activity.

Charge
Charge analysis is a powerful tool to investigate the intramolecular charge transfer that occurs between atoms in a molecule, to identify atoms or groups of atoms that play the role of electron donor or acceptor [18].In the electrostatic theory, atoms with more negative charges donate electrons to atoms with higher positive charges [19].A charge delocalization is observed on the carbon atoms of the quinone ring; therefore, the oxygen atoms O1, O2, and O3 (Figure 1) show the most negative values for the natural products derived from perezone, Table S2 (Supplementary Material).The hydrogen atom of the hydroxyl group of 1 is the most deficient in charge (0.508 e − ).Additionally, the oxygens of the carbonyl groups of the quinone ring and the hydroxyl group display higher charges (−0.494 to −0.654 e − ).Consequently, the hydrogen atom of the hydroxyl group can form hydrogen bonds, confirming the interaction by hydrogen bond with the amino acid residues of PARP-1.In addition, oxygen from the carbonyl of the ester group of compounds 2 and 4 interacts with amino acid residues Arg217 and Tyr235 to form hydrogen bonds, thereby confirming the hydrogen bond interaction with the amino acid residues of PARP-1 and supporting the good PARP-1 inhibitory activity.
Thus, taking into account the molecular dynamics, molecular electrostatic potentials, and charge results, in the next section, the design of a set of sulfur derivatives of perezone is described, to improve the PARP-1 inhibitory activity and to compare the activity with the natural derivatives of perezone (2-4).

Design of Compounds
According to the obtained results, it was possible to make several observations that allow the explanation of the affinity of the studied compounds towards the PARP-1 protein.
In the first instance, the presence of a large and deep cavity in the protein allowed compounds to fit deeply; this fact is promoted by the establishment of non-bonding interactions such as π-π interaction, hydrogen bonds, and hydrophobic contacts.The presence of the angelate moiety in 2 and 4 contributes to an increase in the affinity to PARP-1.Thus, the angelate substituent establishes interactions with a non-polar aliphatic moiety of amino acid residues existent in the binding cavity, indicating that the size of the protein cavity is important for compound selectivity as has been reported for other enzymes [13,14].Also, compounds 2 and 4 establish additional hydrogen bonding, π-π type, and hydrophobic interactions with Tyr235, and Hys201, which could be correlated with their high potency in comparison to 1 and 3.It is also worth noting that these interactions are important because they favor an affinity between the protein and the ligand [13,14].The analysis of the computational studies allows us to propose that hydrophobic moiety substituents could allow the orientation of the quinonic ring of perezone with Hys201 to generate a better π-π interaction, thus developing novel PARP-1 inhibitors, and aliphatic and aromatic sulfur substituents were selected to design a family of perezone derivatives, Figure 5.

Docking Studies
From this study, it is possible to predict and computationally calculate the most favorable position of the interaction of the sulfur perezone (ligand) molecules with the enzyme PARP-1 from their three-dimensional representations [20,21].Figure 6 shows the free binding energies (ΔG) between the PARP-1 protein with each of the derivatives of

Docking Studies
From this study, it is possible to predict and computationally calculate the most favorable position of the interaction of the sulfur perezone (ligand) molecules with the enzyme PARP-1 from their three-dimensional representations [20,21].Figure 6 shows the free binding energies (∆G) between the PARP-1 protein with each of the derivatives of sulfur perezone (−11.78 to −5.71 kcal/mol).This parameter represents the energy that is released when intermolecular interactions between a protein and a ligand occur.In addition, it is observed that the compounds (8-10) show a lower free binding energy, and consequently, they display a higher affinity in the interaction (ligand-protein).In contrast, compounds 5-7 and 11-15 display higher free binding energies; thus, these compounds exhibit a lower interaction with the protein.Therefore, 8-10 were selected to continue with the next stage of study.
. Sulfur radical of designed compounds.

Docking Studies
From this study, it is possible to predict and computationally calculate the most favorable position of the interaction of the sulfur perezone (ligand) molecules with the enzyme PARP-1 from their three-dimensional representations [20,21].Figure 6 shows the free binding energies (ΔG) between the PARP-1 protein with each of the derivatives of sulfur perezone (−11.78 to −5.71 kcal/mol).This parameter represents the energy that is released when intermolecular interactions between a protein and a ligand occur.In addition, it is observed that the compounds (8-10) show a lower free binding energy, and consequently, they display a higher affinity in the interaction (ligand-protein).In contrast, compounds 5-7 and 11-15 display higher free binding energies; thus, these compounds exhibit a lower interaction with the protein.Therefore, 8-10 were selected to continue with the next stage of study.The ligands that have a greater stability of interaction with the protein are those that have a lower energy requirement.In one sense, the sulfur perezone derivative 8 displays the best ∆G value, i.e., −11.78 kcal/mol, followed by 9 and 10, with −9.79 and −9.53 kcal/mol, respectively.Table S3 (Supplementary Material) shows the interaction among 8-10 sulfur perezone molecules and the amino acid residues of PARP-1.It is important to highlight that these molecules interact with two amino acid residues, Hys201 and Tyr235, as was previously described for natural derivatives (2-4).
According to the structures of the sulfur derivatives of perezone, 8-10, all of them exhibit a lower free binding energy, and it is important to highlight that their main feature is that the three compounds bear an aromatic ring; consequently, the π-π interactions could be improved, see above.Moreover, for 8 and 9, there are hydrogen bond interaction types that do not exist in 2-4.In general, it is important to remark that hydrophobic alkyl-π interactions are generated between the double bond of the side chain of perezone and the alkyl moiety of the amino acid residue.Hydrogen bonds are formed, for 8 and 9, between the hydroxyl and carbonyl groups and the hydrogen of the amino groups of the amino acid residues, and for 10, the hydrogen interactions are developed between the side hydrogen of the aliphatic chain of perezone and the carbonyl group of Gly202 and between the carbonyl groups and the hydrogen of the amine group and the alkyl chain of the amino acid residue.Compound 10, in Figure 9, shows three hydrogen bonding interactions with the following amino acid residues: Gly202 (2.3 Å), Gly233 (1.9 Å), and Lys242 (2.0 Å); moreover, one π-π interaction with Tyr246 is detected.Finally, five hydrophobic interactions are observed between the side chain of the ligand with the following amino acid residues: Hys201, Tyr235, Ile234, Gly233, and Ser203.According to the structures of the sulfur derivatives of perezone, 8-10, all of them exhibit a lower free binding energy, and it is important to highlight that their main feature is that the three compounds bear an aromatic ring; consequently, the π-π interactions could be improved, see above.Moreover, for 8 and 9, there are hydrogen bond interaction types that do not exist in 2-4.In general, it is important to remark that hydrophobic alkylπ interactions are generated between the double bond of the side chain of perezone and the alkyl moiety of the amino acid residue.Hydrogen bonds are formed, for 8 and 9, between the hydroxyl and carbonyl groups and the hydrogen of the amino groups of the amino acid residues, and for 10, the hydrogen interactions are developed between the side Additionally, the amino acid residue Tyr235 presents three types of interactions with 8: hydrogen bond, π-π, and hydrophobic.The Hys201 amino acid residue displays two types of interactions with 9: hydrogen bond and hydrophobic.Compound 10 only displays hydrophobic interactions with both amino acids.Thus, it is logical to think that these compounds could offer a better inhibitory activity.In comparison with Olaparib, the compounds 8-10 interact with two amino acid residues from the catalytic pocket, Tyr235 and Hys201, and Olaparib only interacts with Tyr235; however, it is important to highlight that this residue is bonded by a peptidic bond with some amino acids (Ile234, Phe236, Gly233, and Ala237), favoring all types of interactions (hydrogen bond, π-π, and hydrophobic) with almost all Olaparib structures, and probably this fact favors their lower inhibitory concentration, see below.

Synthesis and Spectroscopic Attribution of 9 and 10
To offer a better understanding of both the synthesis and the spectroscopical characterization of the new molecules 9 and 10, Scheme 1 is provided.
The corresponding experimental data of 8 has been previously described by our research group [22].Consequently, the structural attribution for the new molecules 9 and 10 is presented.The compound 9 is obtained as a red oil with 45% yield of pure compound.The infrared spectrophotometric data show a broad medium band centered at 3303 cm −1 , which is assigned to OH groups, a weak band at 2050 cm −1 is assigned to the C=C of the aromatic moiety, two medium bands are observed at 1724 and 1600 cm −1 and are assigned to C=O group and C=C system, and at 739 cm −1 , a medium band is assigned to the C-S bond.The corresponding 1 H NMR shows a double doublet signal at 7.38-7.33ppm and is assigned to H21, and this signal displays ortho coupling with H20, with a J = 6 Hz.At 7.28 ppm, a singlet signal was observed for the OH groups.In the range of 7.25-7.24ppm, a doublet signal is observed for H18, showing a meta coupling with H20 (J = 3 Hz), and at 6.96-6.89ppm, a triplet of doublet signal displaying a meta coupling with H18 (J = 3 Hz) is assigned to H20.A triplet signal at 5.28 ppm is assigned to H12.Singlet signals for H22, H7, H14, and H15 are located at 4.74, 2.19, 1.64, and 1.50 ppm, respectively.Two multiplets are assigned for H8 and H10-H11 at 2.95 and 1.80-1.64ppm, respectively.A doublet signal is observed at 1.09 ppm and is assigned to H9.In the 13 C NMR spectrum, the main signals that agree with the structure of 9 are assigned to the carbonyl groups at 165.4 and 162.1 ppm and the carbonyl groups of the quinonic ring.Also, the signal for aromatic carbons is from 145.7 to 115.0 ppm.The signal at 62.6 ppm is assigned to C22 for methylene carbon.Additionally, the carbon-sulfur base, C6, is observed at 145.8 ppm.Finally, the corresponding mass spectrum (EIMS) shows the ion fragment m/z 404 related to the molecular ion, and this ion is validated by HRMS-DART + data, offering an elemental composition of C 22 H 26 F 1 O 4 S 1 , in agreement with a fragment ion [M + 1] + , and this fragment is correlated with an exact value of 405.15313Da, a precise value of 405.15358Da, and an error of −1.12 ppm, and complementarily, the provided unsaturation data, 10.5, is in agreement with the structure.
Additionally, the amino acid residue Tyr235 presents three types of interactions with 8: hydrogen bond, π-π, and hydrophobic.The Hys201 amino acid residue displays two types of interactions with 9: hydrogen bond and hydrophobic.Compound 10 only displays hydrophobic interactions with both amino acids.Thus, it is logical to think that these compounds could offer a better inhibitory activity.In comparison with Olaparib, the compounds 8-10 interact with two amino acid residues from the catalytic pocket, Tyr235 and Hys201, and Olaparib only interacts with Tyr235; however, it is important to highlight that this residue is bonded by a peptidic bond with some amino acids (Ile234, Phe236, Gly233, and Ala237), favoring all types of interactions (hydrogen bond, π-π, and hydrophobic) with almost all Olaparib structures, and probably this fact favors their lower inhibitory concentration, see below.

Synthesis and Spectroscopic Attribution of 9 and 10
To offer a better understanding of both the synthesis and the spectroscopical characterization of the new molecules 9 and 10, Scheme 1 is provided.The corresponding experimental data of 8 has been previously described by our research group [22].Consequently, the structural attribution for the new molecules 9 and 10 is presented.The compound 9 is obtained as a red oil with 45% yield of pure compound.The infrared spectrophotometric data show a broad medium band centered at 3303 cm −1 , which is assigned to OH groups, a weak band at 2050 cm −1 is assigned to the C=C of the aromatic moiety, two medium bands are observed at 1724 and 1600 cm −1 and are assigned to C=O group and C=C system, and at 739 cm −1 , a medium band is assigned to the C-S bond.The corresponding 1 H NMR shows a double doublet signal at 7.38-7.33ppm and is assigned to H21, and this signal displays ortho coupling with H20, with a J = 6 Hz.At 7.28 ppm, a singlet signal was observed for the OH groups.In the range of 7.25-7.24ppm, a doublet signal is observed for H18, showing a meta coupling with H20 (J = 3 Hz), and at 6.96-6.89ppm, a triplet of doublet signal displaying a meta coupling with H18 (J = 3 Hz) is assigned to H20.A triplet signal at 5.28 ppm is assigned to H12.Singlet signals for H22, H7, H14, and H15 are located at 4.74, 2.19, 1.64, and 1.50 ppm, respectively.Two multiplets Scheme 1. Production and numeration assigned to sulfur derivatives of perezone, 9-10.
Compound 10 is produced as a red with a 65% yield of pure compound.The infrared spectrophotometric data show a broad medium band centered at 3372 cm −1 , which is assigned to the CO 2 H group, a weak band at 2110 cm −1 that is assigned to the C=C of the aromatic moiety, two intense bands are observed at 1707 and 1637 cm −1 and are assigned to the C=O group and the C=C system, and at 738 cm −1 , a band is assigned to the C-S bond.The corresponding 1 H NMR spectrum displays a broad singlet signal at 13.83 ppm, which is unequivocally assigned to the proton of carboxylic acid, and in a range of 7.22-7.19and 7.14-7.12ppm, two doublet signals are presented and assigned to aromatic protons with the coupling data of J = 6 Hz, representative for ortho coupling.At 7.06 ppm, the singlet signal for the OH group is observed.Triplet, singlet, and multiplet signals are presented at 4.95, 3.54, and 2.94 ppm, and are assigned to H12, H22, and H8, respectively.The singlets for H7, H14, and H15 are localized at 2.02, 1.56, and 1.44 ppm, respectively.The signal for H10 and H11 is observed as a multiplet at 1.83-1.63ppm; finally, a doublet signal for H9 is located at 1.07 ppm.In the 13 C NMR spectrum, the main signals that agree with the structure of 10 are assigned to the carbonyl groups at 182.5, 181.7, and 176.7 ppm and are assigned to the carbonyl groups of the quinonic ring and the carboxylic moiety, respectively.Also, the signal at 40.5 ppm is assigned to C22, a methylene carbon type.Additionally, the signal at 150.9 ppm is observed for C6, a carbon-sulfur base.Finally, the EIMS shows the ion fragment m/z 414 in agreement with the molecular ion, and this ion is unequivocally corroborated by the HRMS-DART + study, providing an elemental composition of C 23 H 27 O 5 S 1 , in agreement with a fragment ion [M + 1] + , and this fragment is correlated with an exact value of 415.15669Da and a precise value 415.15792Da (error: −2.96 ppm); complementarily, the provided unsaturation data, 11.5, is in agreement with the structure.of 20

In Vitro Cytotoxic Assay
In terms of the cytotoxic activity, compounds 8 and 9 display an inhibitory activity against MCF-7 and Hela cancer cell lines, Table 1.The inhibitory concentration 50 (IC 50 ) of compound 8 for Hela cells is 10.18 µM, which in comparison with 9 and 10 offers a better response.But compound 9, showing an IC 50 = 83.17µM against MCF-7 cells, when compared with 8 and 10, displays a good performance.In addition, the concentrationresponse charts are provided in Supplementary Material (Figures S1 and S2).

PARP-1 Inhibition Assay
It is well known that the major inhibitor of the PARP-1 enzyme is Olaparib, with an IC 50 of 5.5 nM (±0.5 µM) [13].In one sense, 8-10 exhibit excellent PARP-1 inhibitory activity with IC 50 of 5.205, 0.317, and 0.320 µM, respectively, Table 2. Thus, they show a higher potency than 1 (181.5 µM) and their natural perezone derivatives 2-4 (5.25, 320.3, and 28.7 µM) [13,14].With the aforementioned commentaries, it is important to mention that with the addition of an aromatic substituent to 1, in particular with groups that can form hydrogen bonds or π interactions, the concentration needed to inhibit PARP-1 could diminish.

Quantum Properties and HOMO-LUMO Surface of 2-4 and 8-10
Table 3 shows the results obtained from the energies of the frontier orbitals, which are mainly involved in the chemical stability of the compound.The highest occupied molecular orbital (HOMO) represents the region that donates electrons, and the lowest unoccupied molecular orbital (LUMO) represents the region that accept electrons [23,24].In addition, the bandgap energy (Gap) is calculated.Band gap energy is an indicator of chemical reactivity, hardness, and biological activity of molecular systems.Thus, if the bandgap energy is relatively low, it implies that a molecule can be easily polarized and attains a greater biological activity [25].For a higher band gap energy, the molecule shows a high hardness.According to Table 3, The Gap energies are located between 2.554 and 2.853 eV.Compound 9 has the lowest band gap energy of 2.554 eV; consequently, this molecule is the most reactive of all molecules, according to the order from the highest to the lowest Gap energies the studied molecules: 8 > 2 > 3 > 10 > 4 > 1 > 9.In addition, according to the Principle of Maximum Hardness, it is established that the more reactive systems present low hardness values, while the less reactive systems show high hardness values.The most reactive molecules show a hardness ranging from 2.80 to 2.85 (8)(9)(10).According to the Gap and the hardness data of the seven molecules, they could show chemical and possibly biological activity since they can be reactive to change [25].In general, for all studied compounds 1-4 and 8-10, the HOMO orbital is located on the aliphatic side chain, in addition to the C2 and C3 atoms of the quinonic ring (Figure 10).Moreover, for 9 and 10, the HOMO is present in the quinonic ring and also the sulfur atom.In contrast, the LUMO orbital of all molecules is distributed in the quinone system (Figure 10).These results confirm that 9 is the molecule with a better inhibitory activity towards the PARP-1 enzyme.2.9.Reactivity Properties of 2-4 and 8-10 Table 3 shows the energies of chemical systems in a neutral state and with negative and positive ions.From these data, the reactivity parameters are determined in the DFT scheme.In addition, the electronic affinity (EA) values were calculated for the isolated molecules in the gas phase.The method used for calculating the electron affinity utilizes vertical electron affinities, which are obtained as the difference between the energies of the appropriate neutral form (E Neutral ) and the anion form (E Anion ) at optimized neutral geometries.EA is positive for all molecules explored in this study (1.79-2.14eV).Therefore, the molecules release energy when an electron is added to a neutral molecule.It is known that the ionization potential (IP) is associated with the antioxidant capacity [26] of some compounds since this quantity is directly related to one of the mechanisms that drives the oxidation process, i.e., electron transfer.Thus, as seen in Table 3, all compounds exhibit values for this property, in the range of 7.60-8.13kcal/mol.Due to the lowest IP of 8, it is the easiest to oxidize.The average value of the HOMO and LUMO energies is related to the electronegativity (χ) defined by Mulliken [27] with χ = (IP + EA)/2.Electronegativity is associated with the free movement of electrons in compounds, i.e., electronegativity shows how the electrons will flow from high electronic density regions in a molecule to other sites of lower electronic density.Low global electronegativity values or high chemical potential values mean that an electron is delocalized on the molecule and the molecule can easily supply electrons to interact intra-or inter-molecularly [28].Chemical reactivity increases as electronegativity increases [25].The range of electronegativity is from 4.75 to 4.98.This means that 9 has a high tendency to attract electrons during interaction with another chemical compound.The global electrophilicity index, ω, measures the stabilization in energy when the molecule acquires an additional electronic charge from the environment.The electrophilicity index encompasses both the propensity of the electrophile to acquire an additional electronic charge driven by µ 2 and the resistance of the system to exchange charge with the environment described by η, simultaneously [29].In a way, a good electrophile, such as 9, is characterized by a high value of µ (−4.98 eV) and a low value of η (2.84 eV).The structures of this study show the following order of reactivity: Finally, considering all these reactivity properties, it is appropriate to comment that the two new compounds, 9 and 10, display a better PARP-1 inhibitory activity, and these properties support the previous results, highlighting that 9 is a good PARP-1 inhibitor candidate.

Molecular Dynamics Simulations of Four PARP-1 Inhibitors Related to Perezone Derivatives
The nature of ligand-receptor interaction, the binding affinity of the ligand towards the receptor, and the ligand stability in the active site of the receptor are basic requirements for rational drug design.For this purpose, the three-dimensional (3D) structure of PARP-1 was retrieved from the Protein Data Bank (PDB ID: 1UK0).The water molecules, co-crystalized ligands, and peptides were eliminated.The structures of perezone, perezone angelate, hydroxyperezone, and hydroxyperezone monoangelate were optimized using quantum chemistry.The starting structure of the complexes was obtained by docking studies employing Autodock 4.2 [30].A GRID-based procedure was utilized to prepare the structural inputs and define all the binding sites.A rectangular lattice (70 × 70 × 70 Å) with points separated by 0.375 Å was centered on the active site of PARP-1.All docking simulations were conducted using the hybrid Lamarckian genetic algorithm with an initial population of 100 randomly placed individuals and a maximum of 1.0 × 10 7 energy evaluations.All other parameters were maintained at their default settings.The lowest energy cluster for each protein-ligand was selected as the starting structure for MD simulations using the Nanoscale Molecular Dynamics (NAMD) Simulation 2.6 program (http://www.ks.uiuc.edu/Research/namd/accessed on 17 October 2023) (Champaign-Urbana, IL, USA) [31], employing the CHARMM27 force field [32].
For the MD simulations of PARP-1 in absence of the ligands, hydrogen atoms were added to the initial coordinates of the protein using the psfgen program included in the Visual Molecular Dynamic (VMD) 1.8 program [33].Afterwards, the system was solvated using TIP3P water molecules, and these structures were neutralized using X Na + .The equilibration protocol consisted of 1500 minimization steps followed by 30 ps of MD simulations at 10 K with fixed protein atoms.Subsequently, the entire system was minimized over 1500 steps (at 0 K), followed by gradual heating from 10 to 310 K using temperature reassignment during the initial 60 ps of the 100 ps equilibration molecular dynamics without restraints.
The RMSD was analyzed to determine whether the protein had undergone conformational changes because this value reflects the distance between the pairs of the same atoms, represented by δ with time, Equation (1).Furthermore, RMSD values depicted the equilibration phases that determine the MD simulation quality.The RMSF indicated the displacement of each particle (Cα) to the structure, Equation (2).Therefore, the RMSF average was taken over time to provide a displacement for each particle (Cα per residue), whereas the RMSD average was calculated over the whole particle to provide time-specific values.
All computational work was performed using CUDA on an Intel Core i7-980x 3.33 GHz.Linux workstation with 12 Gb of RAM, (2×) NVIDIA GeForce GTX580 video cards, and a (1×) NVIDIA GeForce GTX680 video card.The parameters for perezone and related compounds to perform the MD simulations of the ligand-protein complex were obtained by employing Swiss Param Server [35].

Design of Compounds
According to MD simulations results to evaluate the interactions that drive the recognition of PARP-1 to compound 1 to 4, a family of the perezone derivatives was designed.For this purpose, the sulfur substituent, Figure 5, was incorporated instead of the hydrogen atom on C6 of 1. Perezone derivatives were submitted to docking studies with PARP-1, according to the previously mentioned procedure.The lowest energy cluster for each ligand was subjected to further free energy and binding geometry analyses, as was previously reported [53].Conformations with the lowest free energy binding (∆G) and the highest frequency were selected employing Autodock tools [30].The images were created employing PyMol [54].Starting materials phenylmethanethiol, 4-mercaptophenylacetic acid, and 5-fluoro-2mercaptobenzyl alcohol were purchased from Sigma Aldrich Chemistry (St. Louis, MO, USA) and used as received.The solvents n-hexane, ethyl acetate, and absolute methanol were technical grade and were used without further purification from Materiales y Abastos Especializados S.A. de C.V. (Zapopan, Jalisco, Mexico).Perezone was isolated from the Acourtia platyphilla specimen and identified by their corresponding spectroscopical data correlated with the literature [55,56].The reactions were monitored using thin-layer chromatography (TLC) performed on pre-coated Merck silica gel 60F254 aluminum sheets, and the visualization was achieved using a 254 nm UV lamp (UVLS-24, Upland, CA, USA), and purified with column chromatography using silica gel (0.063-0.200 mm, 70-230 mesh ASTM, acquired from Merck-Millipore, Germany). 1 H and 13 C (CDCl 3 ) spectra were recorded using a Bruker Advance III spectrometer (San Francisco, CA, USA) at 300 MHz and 75 MHz, respectively.The multiplicities are reported as singlet (s), broad singlet (bs), doublet (d), triplet (t), double of double (dd), triplet of double (td), and multiplet (m).The corresponding δ-chemical shifts are given in ppm, employing TMS as internal reference.The EIMS (70 eV) were determined using a JEOL JMS-700 MStation mass spectrometer (JEOL, Tokyo, Japan).The HRMS-DART + (19 eV) data were obtained using a JEOL Ac-cuTOF JMS-T100LC (Direct Analysis in Real Time) spectrometer.The measurements were performed using the DART + experiment with PEG (polyethylene glycol) 400 as an internal reference at 6000 resolutions and triplet helium as a carrier gas at 350 • C. In the first orifice, the temperature and voltage were 120 • C and 15 V, respectively, and the voltage in the second orifice was 5 V. Elemental composition was calculated within a mass range ±10 ppm from the measured mass.The IR spectra were recorded on a Bruker TENSOR 27 spectrophotometer (Bruker, Billerica, MA, USA).

Synthesis of 9 and 10
In an appropriate glass vessel, 1.0080 mmol (250 mg) of perezone, 0.98101 mmol (155 mg) of 5-fluoro-2-mercaptobenzyl alcohol or 1.0119 mmol (170 mg) of 4-mercaptophenylacetic acid, and 10 mL of absolute methanol were mixed, and then, the mixtures were treated under reflux by 8 h.The reaction progress was monitored using TLC (n-hexane/ethyl acetate, The corresponding products were purified with column chromatography using the same mobile phase as the TLC and were obtained as red oils, Scheme 1.

Cytotoxic Evaluation on Cancer Cell Lines
To determine the inhibitory concentration 50 (IC 50 ) of the compounds, the MTT assay (Creative Biolabs, Shirley, NY, USA, life technologies M6494) was performed on the A-549, MDA-MB-231, MCF-7, and Hela cell lines, which were cultured on 96-well microtiter plates (Glendale, AZ, USA, Corning 009017) with DMEM F-12 medium (Thermo Fisher Scientific, Waltham, MA, USA, gibco 11320-033) with 1% Antibiotic-Antimycotic (100-079) enriched with 10% FBS (fetal bovine serum, gibco 16141-079) and were incubated at 37 • C and 5% CO 2 until reaching 90% confluency.Eight concentrations of the compounds 8-10 dissolved in culture medium were added in 1:2 dilutions, starting from 400 µM (400, 200, 100, 50, 25, 12.5, 6.25, and 3.125 µM) and keeping the remaining wells as controls (DMEM F-12 without FBS) and incubated for 48 h.At the end of this time, 20 µL of the MTT solution (5 was added to each well and incubated for 4 h in the dark at 37 • C and 5% CO 2 .The supernatant was decanted, and the resulting formazan was dissolved in isopropyl alcohol with 4 mM HCl solution and was read in a UV/Vis reader at 590 nm.The results obtained were plotted, and the IC 50 was obtained through linear regression [13].

Conclusions
In the present study, in the first stage, the perezone (1), perezone angelate (2), hydroxyperezone (3), and hydroxyperezone monoangelate (4) were studied with molecular dynamic to verify their inhibitory activity on PARP-1 protein and analyze their potential as antineoplastic agents.In addition, the stability of protein-ligand complex structures during the simulation period of 50 ns was confirmed.Also, the compounds were subjected to theoretical studies: molecular electrostatic potential and charges that were obtained with the Density Functional Theory using the B3LYP hybrid method and with the set of base functions 6-311++G(d,p), to know the sites that could interact with the atoms of the amino acid residues of the PARP-1 protein.In terms of the obtained results from the molecular dynamics and quantum chemistry studies of compounds 1-4, a set of eleven sulfur derivatives of perezone was proposed.Thus, the docking studies of the strategic compound set indicated that molecules 8-10 showed a higher ∆G values.Consequently, they were selected to be synthesized, characterized, evaluated in vitro, and analyzed at the quantum level (molecular orbitals and reactivity properties).In vitro testing with PARP-1 showed that the IC 50 of PARP-1 of compounds 9 (0.317 µM) was lower compared to that of compounds 1-4, 8, and 10.In addition, 9 exhibited inhibitory activities against MDA-MB-231 and MCF-7 cells.
According to the studies of the reactivity properties and the frontier orbitals of systems 1-4 and 8-10, it was observed that compounds 8-10 were more reactive, and this improved their affinity for PARP-1, and therefore, the IC 50 showed by the sulfur derivates of perezone was lower.With the development of this study, compounds 8-10 were determined to be prospective antineoplastic candidates.

Data Availability Statement:
The data reported in this study are available upon request to nicovain@yahoo.com.mx.

Figure 2 .
Figure 2. Molecular dynamics simulations of PARP-1 in the absence or the presence of perezone and their related compounds: (a) RMSD and (b) RMSF.

Figure 2 .
Figure 2. Molecular dynamics simulations of PARP-1 in the absence or the presence of perezone and their related compounds: (a) RMSD and (b) RMSF.

Figure 6 .
Figure 6.Docking results, and ∆G values were obtained through the docking studies of the sulfur perezone derivatives' ligands with the protein PARP-1.

Figure 7 .
Figure 7. (a) It shows the catalytic site of the interaction of the protein PARP-1 with ligand 8. (b) It shows amino acid residues that interact with ligand 8.

Figure 7 . 21 Figure 8 .
Figure 7. (a) It shows the catalytic site of the interaction of the protein PARP-1 with ligand 8. (b) It shows amino acid residues that interact with ligand 8. Int.J. Mol.Sci.2024, 25, x FOR PEER REVIEW 9 of 21

Figure 8 .
Figure 8.(a) It shows the catalytic site of the interaction of the protein PARP-1 with ligand 9. (b) It shows amino acid residues that interact with ligand 9.

Figure 8 .
Figure 8.(a) It shows the catalytic site of the interaction of the protein PARP-1 with ligand 9. (b) It shows amino acid residues that interact with ligand 9.

Figure 9 .
Figure 9. (a) It shows the catalytic site of the interaction of the protein PARP-1 with ligand 10.(b) It shows amino acid residues that interact with ligand 10.

Figure 9 .
Figure 9. (a) It shows the catalytic site of the interaction of the protein PARP-1 with ligand 10.(b) It shows amino acid residues that interact with ligand 10.
): elemental composition of C 22 H 26 F 1 O 4 S 1 [M + 1] + for protonated molecule, exact value of 405.15313Da, and precise value of 405.15358Da, with an error of −1.12 ppm, and complementarily, the provided unsaturation data of 10.5.