<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Chem.</journal-id>
<journal-title>Frontiers in Chemistry</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Chem.</abbrev-journal-title>
<issn pub-type="epub">2296-2646</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1384832</article-id>
<article-id pub-id-type="doi">10.3389/fchem.2024.1384832</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Chemistry</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>3D computer modeling of inhibitors targeting the MCF-7 breast cancer cell line</article-title>
<alt-title alt-title-type="left-running-head">Zarougui et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fchem.2024.1384832">10.3389/fchem.2024.1384832</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Zarougui</surname>
<given-names>Sara</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2654439/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Er-Rajy</surname>
<given-names>Mohammed</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2654440/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Faris</surname>
<given-names>Abdelmoujoud</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2373562/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Imtara</surname>
<given-names>Hamada</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2371732/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>El fadili</surname>
<given-names>Mohamed</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Qurtam</surname>
<given-names>Ashraf Ahmed</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Nasr</surname>
<given-names>Fahd A.</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2700850/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Al-Zharani</surname>
<given-names>Mohammed</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1429601/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Elhallaoui</surname>
<given-names>Menana</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Laboratory of Engineering, Modelisation and Systems Analysis</institution>, <institution>Department of Chemical Sciences</institution>, <institution>Faculty of Sciences Dhar El Mahraz</institution>, <institution>Sidi Mohamed Ben Abdellah University</institution>, <addr-line>Fez</addr-line>, <country>Morocco</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Faculty of Medicine</institution>, <institution>Arab American University Palestine</institution>, <addr-line>Jenin</addr-line>, <country>Palestine</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Department of Biology</institution>, <institution>College of Science</institution>, <institution>Imam Mohammad Ibn Saud Islamic University (IMSIU)</institution>, <addr-line>Riyadh</addr-line>, <country>Saudi Arabia</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/308363/overview">Mohammad Ehtisham Khan</ext-link>, Jazan University, Saudi Arabia</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1217954/overview">Khaled Mohamed Darwish</ext-link>, Suez Canal University, Egypt</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/688234/overview">Vikram Dalal</ext-link>, Washington University at St. Louis, United States</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Hamada Imtara, <email>Hamada.imtara@aaup.edu</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>23</day>
<month>05</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>12</volume>
<elocation-id>1384832</elocation-id>
<history>
<date date-type="received">
<day>10</day>
<month>02</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>11</day>
<month>04</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Zarougui, Er-Rajy, Faris, Imtara, El fadili, Qurtam, Nasr, Al-Zharani and Elhallaoui.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Zarougui, Er-Rajy, Faris, Imtara, El fadili, Qurtam, Nasr, Al-Zharani and Elhallaoui</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>This study focused on developing new inhibitors for the MCF-7 cell line to contribute to our understanding of breast cancer biology and various experimental techniques. 3D QSAR modeling was used to design new tetrahydrobenzo[4, 5]thieno[2, 3-d]pyrimidine derivatives with good characteristics. Two robust 3D-QSAR models were developed, and their predictive capacities were confirmed through high correlations [CoMFA (Q<sup>2</sup> &#x3d; 0.62, <italic>R</italic>
<sup>2</sup> &#x3d; 0.90) and CoMSIA (Q<sup>2</sup> &#x3d; 0.71, <italic>R</italic>
<sup>2</sup> &#x3d; 0.88)] via external validations (R<sup>2</sup>
<sub>ext</sub> &#x3d; 0.90 and R<sup>2</sup>
<sub>ext</sub> &#x3d; 0.91, respectively). These successful evaluations confirm the potential of the models to provide reliable predictions. Six candidate inhibitors were discovered, and two new inhibitors were developed <italic>in silico</italic> using computational methods. The ADME-Tox properties and pharmacokinetic characteristics of the new derivatives were evaluated carefully. The interactions between the new tetrahydrobenzo[4, 5]thieno[2, 3-d]pyrimidine derivatives and the protein ER&#x3b1; (PDB code: 4XO6) were highlighted by molecular docking. Additionally, MM/GBSA calculations and molecular dynamics simulations provided interesting information on the binding stabilities between the complexes. The pharmaceutical characteristics, interactions with protein, and stabilities of the inhibitors were examined using various methods, including molecular docking and molecular dynamics simulations over 100&#xa0;ns, binding free energy calculations, and ADME-Tox predictions, and compared with the FDA-approved drug capivasertib. The findings indicate that the inhibitors exhibit significant binding affinities, robust stabilities, and desirable pharmaceutical characteristics. These newly developed compounds, which act as inhibitors to mitigate breast cancer, therefore possess considerable potential as prospective drug candidates.</p>
</abstract>
<kwd-group>
<kwd>MM/GBSA</kwd>
<kwd>3D-QSAR</kwd>
<kwd>molecular docking</kwd>
<kwd>molecular dynamics simulation</kwd>
<kwd>MCF-7 cell line</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Medicinal and Pharmaceutical Chemistry</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>All living beings, including humans, animals, and various organisms, can be affected by cancer; this disease does not exhibit favoritism and can affect individuals in different ways. Both men and women are vulnerable to cancer, but there are certain types of cancers that are more common in each gender. Prostate, lung, and colon cancers are frequently diagnosed in men, whereas breast, lung, and colon cancers are more prevalent in women. Young individuals are often affected by leukemia, brain tumors, and lymphomas as the most common types of cancer. Sadly, certain types of cancers consistently contribute toward a significant number of deaths, i.e., lung, stomach, liver, and colon cancers, which require urgent attention to reduce fatalities (<xref ref-type="bibr" rid="B43">Zaorsky et al., 2017</xref>). There are various cancer treatment options available that can be categorized as either local or systemic. Local treatments include surgery, hyperthermia therapy, and radiation, whereas systemic treatments comprise chemotherapy, hormone therapy, and targeted therapy. In most cases, individuals undergo a combination of treatments, such as surgery accompanied by chemotherapy and/or radiotherapy. Additionally, neoadjuvant chemotherapy, which involves using chemotherapy to reduce the size of a tumor before radiation or surgery, is another approach that is utilized frequently. Our focus here is on chemotherapy involving the use of drugs to address cancerous growths in various parts of the body (<xref ref-type="bibr" rid="B22">Falzone et al., 2018</xref>). The central fraction with anticancer properties is mainly composed of the pyrimidine group; additionally, the thieno[2,3-d]pyrimidine ring, which is a bioisostere of quinazoline, has been used as a central structure in the design and synthesis of various compounds with potential antitumor activities (<xref ref-type="bibr" rid="B12">Dai et al., 2005</xref>). The thienopyrimidine nucleus has also been utilized as a molecular framework in the development of numerous antimicrobial (<xref ref-type="bibr" rid="B17">El-Sherbeny et al., 1995</xref>; <xref ref-type="bibr" rid="B5">Aly et al., 2011</xref>), antiviral (<xref ref-type="bibr" rid="B17">El-Sherbeny et al., 1995</xref>; <xref ref-type="bibr" rid="B37">Rashad et al., 2010</xref>), anti-inflammatory (<xref ref-type="bibr" rid="B3">Alagarsamy et al., 2006</xref>), antidiabetic (<xref ref-type="bibr" rid="B15">Deng et al., 2011</xref>), antioxidant (<xref ref-type="bibr" rid="B29">Kotaiah et al., 2012</xref>), and anxiolytic (<xref ref-type="bibr" rid="B6">Amr et al., 2010</xref>) agents. Prior studies in literature have documented that the antitumor activity of the thieno[2,3-day]pyrimidine core can be maintained by fusing it with cycloalkyl lipophilic moieties with five, six, or seven members (<xref ref-type="bibr" rid="B45">Zhang et al., 2020</xref>; <xref ref-type="bibr" rid="B35">Ouassaf et al., 2021</xref>). This is evident in the case of the 4-substituted-5,6,7,8-tetrahydrobenzo[4,5]thieno[2,3-day]pyrimidine, which exhibits important anticancer activity (<xref ref-type="bibr" rid="B8">Aponte et al., 2010</xref>).</p>
<p>Breast cancer is now widely recognized as one of the leading types of cancer diagnosed in women across the globe. Malignant tumors in women present significant risks, often resulting in considerable illness and death. Moreover, statistics show that this type of cancer is relatively rare in men (<xref ref-type="bibr" rid="B16">Elbachiri et al., 2017</xref>). Breast cancer can be influenced by various risk factors, such as gender, diet, breast density, age, hormonal factors, radiation exposure, obesity, mutations, and family history (<xref ref-type="bibr" rid="B33">Maryam et al., 2017</xref>). The US Food and Drug Administration has approved the drug capivasertib (<xref ref-type="fig" rid="F1">Figure 1</xref>) as an AKT inhibitor that works by inhibiting the AKT1 enzyme. By targeting this key signaling pathway, it shuts down the machinery driving cancer cell growth and resilience. In prior studies, MCF-7 cell line cultures were subjected to varying doses of capivasertib to evaluate its effects; as the concentration of the tested agent increased, the cell viability decreased progressively, indicating a dose-dependent inhibitory effect. Higher doses of capivasertib resulted in greater inhibition. Subsequent experiments have revealed that capivasertib treatment reduces cell growth, causes cell cycle arrest, and facilitates cell death. Blocking the AKT1 pathway with capivasertib disrupts the key processes that the MCF-7 cells rely upon to grow and proliferate unchecked like cancer cells (<xref ref-type="bibr" rid="B27">Hopcroft et al., 2023</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Structures of the 2,3-disubstituted derivatives (3a-z and 4a-d), 2,4-disubstituted derivatives (6a-c), and FDA-approved drug capivasertib.</p>
</caption>
<graphic xlink:href="fchem-12-1384832-g001.tif"/>
</fig>
<p>In the present work, we use the three-dimensional quantitative structure&#x2013;activity relationship (3D-QSAR) with the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) descriptors to determine the key characteristics for the potential inhibition of MCF-7 cells by tetrahydrobenzo[4,5]thieno[2,3-d]pyrimidine derivatives (<xref ref-type="bibr" rid="B18">Er-rajy et al., 2022</xref>; <xref ref-type="bibr" rid="B19">Er-rajy et al., 2023a</xref>). In addition, the ADME-Tox (ADMET) studies and pharmacokinetic characteristics were used to predict the pharmaceutical features of the newly designed molecules for aspects crucial to drug effectiveness and pharmacokinetics, including solubility, permeability, and metabolic stability (<xref ref-type="bibr" rid="B4">Alqahtani, 2017</xref>). This research also involves investigating the interactions between the protein ER&#x3b1; and ligands using the molecular docking approach. The aim of this investigation is to discover promising compounds that can selectively inhibit the MCF-7 cell line. To achieve this, we conducted molecular dynamics (MD) and molecular mechanics (MM) simulations. Additionally, we utilized MM and generalized Born surface area (MM/GBSA) calculations (<xref ref-type="bibr" rid="B38">Sanachai et al., 2022</xref>). Specifically, we performed MD simulations for a duration of 100 ns to evaluate the stabilities of the derived complexes.</p>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>Materials and methods</title>
<sec id="s2-1">
<title>Dataset</title>
<p>Twenty-nine derivatives of the novel tetrahydrobenzo[4, 5]thieno[2, 3-d]pyrimidine molecule that showed breast cancer (MCF-7) inhibitory activity were collected from literature (<xref ref-type="bibr" rid="B1">Abbas et al., 2013</xref>). The database was randomly split into a training set of 24 derivatives and a test set of five derivatives. The half-maximal inhibitory concentration (IC<sub>50</sub>) values, which measure the abilities of the derivatives to inhibit breast cancer activity, were converted to negative logarithms (pIC<sub>50</sub> &#x3d; -log (IC<sub>50</sub>)) for use as the dependent variables during modeling. The investigated structures and their corresponding pIC<sub>50</sub> values are shown in <xref ref-type="table" rid="T1">Table 1</xref>. In addition, all molecular structures were analyzed with Powell&#x2019;s and Gasteiger&#x2013;Huckel partial atomic charge using the Tripos force field in SYBYL-X.2.1 software (<xref ref-type="bibr" rid="B11">Certara, 2017</xref>).</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Structures and anti-MCF-7 activities of the molecules studied.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">N&#xb0;</th>
<th align="center">Structure</th>
<th align="center">pIC<sub>50</sub>
</th>
<th align="center">N&#xb0;</th>
<th align="center">Structure</th>
<th align="center">pIC<sub>50</sub>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">3a</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx1.tif"/>
</td>
<td align="center">4.474</td>
<td align="center">3q</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx2.tif"/>
</td>
<td align="center">4.712</td>
</tr>
<tr>
<td align="center">3b</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx3.tif"/>
</td>
<td align="center">4.962</td>
<td align="center">3r</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx4.tif"/>
</td>
<td align="center">6.698</td>
</tr>
<tr>
<td align="center">3c</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx5.tif"/>
</td>
<td align="center">4.512</td>
<td align="center">3s</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx6.tif"/>
</td>
<td align="center">5.180</td>
</tr>
<tr>
<td align="center">3e</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx7.tif"/>
</td>
<td align="center">4.707</td>
<td align="center">3u</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx8.tif"/>
</td>
<td align="center">4.809</td>
</tr>
<tr>
<td align="center">3f</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx9.tif"/>
</td>
<td align="center">4.583</td>
<td align="center">3v</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx10.tif"/>
</td>
<td align="center">5.327</td>
</tr>
<tr>
<td align="center">3h</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx11.tif"/>
</td>
<td align="center">4.600</td>
<td align="center">3w</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx12.tif"/>
</td>
<td align="center">5.721</td>
</tr>
<tr>
<td align="center">3i</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx13.tif"/>
</td>
<td align="center">4.395</td>
<td align="center">3x</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx14.tif"/>
</td>
<td align="center">5.142</td>
</tr>
<tr>
<td align="center">3j</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx15.tif"/>
</td>
<td align="center">5.236</td>
<td align="center">3y</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx16.tif"/>
</td>
<td align="center">6.221</td>
</tr>
<tr>
<td align="center">3k</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx17.tif"/>
</td>
<td align="center">4.467</td>
<td align="center">3z</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx18.tif"/>
</td>
<td align="center">7</td>
</tr>
<tr>
<td align="center">3L</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx19.tif"/>
</td>
<td align="center">4.924</td>
<td align="center">4b</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx20.tif"/>
</td>
<td align="center">5.060</td>
</tr>
<tr>
<td align="center">3m</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx21.tif"/>
</td>
<td align="center">6.301</td>
<td align="center">4c</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx22.tif"/>
</td>
<td align="center">5.376</td>
</tr>
<tr>
<td align="center">3n</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx23.tif"/>
</td>
<td align="center">5.408</td>
<td align="center">4d</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx24.tif"/>
</td>
<td align="center">5.096</td>
</tr>
<tr>
<td align="center">3o</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx25.tif"/>
</td>
<td align="center">5.327</td>
<td align="center">6a</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx26.tif"/>
</td>
<td align="center">5.214</td>
</tr>
<tr>
<td align="center">3p</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx27.tif"/>
</td>
<td align="center">5.055</td>
<td align="center">6b</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx28.tif"/>
</td>
<td align="center">5.619</td>
</tr>
<tr>
<td align="center">6c</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx29.tif"/>
</td>
<td align="center">6.522</td>
<td align="left"/>
<td align="left"/>
<td align="left"/>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2-2">
<title>Molecular alignment</title>
<p>In the 3D-QSAR analysis, proper alignment of the molecules is of utmost importance. Molecular alignment plays a critical role in CoMFA and CoMSIA, and it has been widely recognized as a significant factor in these methods, demonstrating great sensitivity (<xref ref-type="bibr" rid="B44">Zarougui et al., 2023</xref>). As the first step in building the model, we selected the most potent compound (3z) with a pIC<sub>50</sub> value of 7 as the template structure. Subsequently, all other molecules were aligned to this template using the distill module in SYBYL 2.1 software by taking into consideration their structural similarities.</p>
</sec>
<sec id="s2-3">
<title>3D-QSAR modeling</title>
<p>To determine the key structural properties that affected the biological activities of the compounds under investigation, 3D-QSAR models were created using the CoMFA and CoMSIA techniques (<xref ref-type="bibr" rid="B7">Aparoy et al., 2011</xref>; <xref ref-type="bibr" rid="B39">Tong et al., 2021</xref>). These techniques enabled examination of various molecular characteristics, such as electrostatic, steric, hydrogen-bond donor, hydrophobic, and hydrogen-bond acceptor, and their potential effects on biological activity (<xref ref-type="bibr" rid="B35">Ouassaf et al., 2021</xref>). In this research, the quantitative relationships between the 3D field descriptors and MCF-7 breast cancer cell line activities (pIC<sub>50</sub>) were studied through the application of the partial least-squares (PLS) algorithm within the Tripos force field. Specifically, the Tripos force field utilized a reference input lattice defined at 2&#xa0;&#xc5; in all Cartesian directions, with the sp<sup>3</sup> hybrid carbon atoms serving as the source for the steric and electrostatic energy calculations, and the default cutoff energy set to 30&#xa0;kcal/mol (<xref ref-type="bibr" rid="B26">Goodhue et al., 2006</xref>). In addition, appropriate QSAR models representing the structure&#x2013;activity relationships of the tested compounds were chosen based on different statistical criteria, notably the coefficient of determination (<italic>R</italic>
<sup>2</sup>), coefficient of determination by cross-validation (LOO-CV Q<sup>2</sup>), standard error of estimate (SEE), and coefficient of determination calculated for all tests (R<sup>2</sup>
<sub>pred</sub>) (<xref ref-type="bibr" rid="B24">Golbraikh and Tropsha, 2000</xref>; <xref ref-type="bibr" rid="B25">Golbraikh and Tropsha, 2002</xref>).</p>
</sec>
<sec id="s2-4">
<title>Molecular docking</title>
<p>Our primary aim in the molecular docking study was to employ computational techniques to forecast the structures of the complexes formed between ligands and their receptors. This process entailed two main phases. First, we sampled various conformations of the ligand within the receptor&#x2019;s active site, followed by assessments of these conformations using a scoring function. For our investigation, we utilized AutoDock Tools 1.5.6 (<xref ref-type="bibr" rid="B34">My Biosoftware, 2019</xref>) and assigned Kollman charges to the protein (3.372) and Gasteiger charges to the 3M (&#x2212;1) and 5M (&#x2212;1) molecules (<xref ref-type="bibr" rid="B30">Kumari and Dalal, 2022</xref>), while incorporating polar hydrogen atoms. This enabled the removal of water molecules and facilitated visualization of the ligand&#x2013;protein interactions using Discovery studio (<xref ref-type="bibr" rid="B14">Dassault Syst&#xe8;mes, 2024</xref>). Molecular docking was implemented using the cocrystalline structure of the ER&#x3b1; protein (PDB code: 4XO6) as the model. Our focus was specifically on the newly designed compounds with high biological activities (pIC<sub>50</sub>), in addition to drug similarity assessments and ADMET studies, for the molecular docking simulations with the MCF-7 cell line. The 3D-grid defined using AUTOGRID algorithms was crucial for determining the strengths of attraction between the ligands and receptors. The grid center was located at (x &#x3d; 12.444, y &#x3d; &#x2212;0.171, z &#x3d; 2.943), with default dimensions of (x &#x3d; 60, y &#x3d; 60, z &#x3d; 60) and grid point spacing of 0.375&#xa0;&#xc5;. Finally, after preparing both the ligands and receptors, we conducted molecular docking on the two newly designed molecules, M3 and M5, using AutoDock Tools (<xref ref-type="bibr" rid="B34">My Biosoftware, 2019</xref>).</p>
</sec>
<sec id="s2-5">
<title>MD simulation</title>
<p>The newly created compounds that exhibited enhanced binding affinities with ER&#x3b1; were subjected to all-atom MD simulations using the Groningen Machine for Chemical Simulation (GROMACS 2021) software (<xref ref-type="bibr" rid="B41">Van Der Spoel et al., 2005</xref>; <xref ref-type="bibr" rid="B2">Abraham et al., 2015</xref>). Before implementing the MD simulations, the CHARMM-GUI web server (<xref ref-type="bibr" rid="B28">Jo et al., 2008</xref>) was used to generate the initial input parameters by applying the CHARMM36 force (<xref ref-type="bibr" rid="B47">Ziada et al., 2022</xref>). The simulations were carried out at pH 7. Before moving to the generating phase, each complex was solved in a rectangular grid encircled by TIP3P water molecules and complemented with the necessary counter ions (Na<sup>&#x2b;</sup>, Cl<sup>&#x2013;</sup>) to maintain a salt concentration of 0.15&#xa0;M, as determined by Monte Carlo ion shifting. Energy minimization was performed for each system using the steepest downslope algorithm with a top of 50,000 steps and top force of 10.0&#xa0;kJ/mol. The temperature and atmospheric pressure were set to 310&#xa0;K and 1.01325 bar, respectively. Canonical (NVT) and isothermal&#x2013;isobaric (NPT) assemblies were used to equilibrate each system, and the MD simulations were performed for a duration of 200 ns. To evaluate the structural stabilities of the newly designed molecules, several parameters like the root mean-squared deviation (RMSD), solvent-accessible surface area (SASA), radius of gyration (RoG), and root mean-squared flexibility (RMSF) were tested for the dynamic path results (<xref ref-type="bibr" rid="B31">Kumari et al., 2014</xref>; <xref ref-type="bibr" rid="B23">Faris et al., 2023</xref>).</p>
</sec>
<sec id="s2-6">
<title>Free binding energy (MM/GBSA)</title>
<p>The MM/GBSA technique is used to quantify the amount of energy released or absorbed when a ligand binds to a protein, which is known as the free binding energy of a complex (<xref ref-type="bibr" rid="B46">Zhang et al., 2017</xref>). We used AMBER 14 software (<xref ref-type="bibr" rid="B10">Case et al., 2023</xref>) to determine the binding free energy (<xref ref-type="bibr" rid="B9">Aqvist, 1990</xref>; <xref ref-type="bibr" rid="B46">Zhang et al., 2017</xref>) using the following equation: &#x394;G<sub>binding</sub> &#x3d; G<sub>complex</sub> &#x2013; (G<sub>receptor</sub> &#x2b; G<sub>ligand</sub>); here, G<sub>receptor</sub> is the free energy of the open receptor and G<sub>ligand</sub> is the ligand&#x2019;s boundaryless free energy that can be described using thermodynamics. The binding free energy (&#x394;G<sub>binding</sub>) is also given by the equation &#x394;G<sub>binding</sub> &#x3d; &#x394;H &#x2013; T&#x394;S, where &#x394;H represents the affinity enthalpy and &#x394;S is the change in entropy after binding (<xref ref-type="bibr" rid="B42">Wang et al., 2019</xref>; <xref ref-type="bibr" rid="B40">Vald&#xe9;s-Tresanco et al., 2021</xref>). The calculations were performed using the last 200 images by applying the following equations:<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mi mathvariant="bold">&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold">G</mml:mi>
<mml:mtext mathvariant="bold">bind</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold">G</mml:mi>
<mml:mtext mathvariant="bold">complex</mml:mtext>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi mathvariant="bold">G</mml:mi>
<mml:mtext mathvariant="bold">protein</mml:mtext>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi mathvariant="bold">G</mml:mi>
<mml:mtext mathvariant="bold">ligand</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:mi mathvariant="bold">&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold">G</mml:mi>
<mml:mtext mathvariant="bold">bind</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold">G</mml:mi>
<mml:mtext mathvariant="bold">gas</mml:mtext>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold">&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold">G</mml:mi>
<mml:mtext mathvariant="bold">sol</mml:mtext>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mo>&#x2013;</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="bold">T</mml:mi>
<mml:mi mathvariant="bold">&#x394;</mml:mi>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
<disp-formula id="e3">
<mml:math id="m3">
<mml:mrow>
<mml:mi mathvariant="bold">&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold">G</mml:mi>
<mml:mtext mathvariant="bold">gas</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mtext mathvariant="bold">Bond</mml:mtext>
<mml:mo>&#x2b;</mml:mo>
<mml:mtext mathvariant="bold">Angle</mml:mtext>
<mml:mo>&#x2b;</mml:mo>
<mml:mtext mathvariant="bold">Dihed</mml:mtext>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold">E</mml:mi>
<mml:mtext mathvariant="bold">EL</mml:mtext>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mtext mathvariant="bold">VdW</mml:mtext>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
<disp-formula id="e4">
<mml:math id="m4">
<mml:mrow>
<mml:mi mathvariant="bold">&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold">G</mml:mi>
<mml:mtext mathvariant="bold">sol</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold">E</mml:mi>
<mml:mtext mathvariant="bold">GB</mml:mtext>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold">&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold">E</mml:mi>
<mml:mtext mathvariant="bold">SURF</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>Here, &#x394;G<sub>bind</sub> represents the total binding energy change of the system, as given by Eq. <xref ref-type="disp-formula" rid="e1">(1),</xref> and is determined using Eq. <xref ref-type="disp-formula" rid="e2">(2)</xref>; T&#x394;S denotes the transformational entropy change resulting from ligand binding at a certain temperature; &#x394;G<sub>gas</sub> is the combined contributions from the bond, angle, dihedral, and van der Waals energies as well as the electrostatic component of the internal energy (E<sub>EL</sub>), as explained in Eq. <xref ref-type="disp-formula" rid="e3">(3)</xref>. The internal energy is linked to fluctuations and rotations in the torsion angles of the individual bonds. The solvation energy (&#x394;G<sub>sol</sub>) is composed of the polar property of solvation energy (&#x394;E<sub>GB</sub>) and non-polar property of solvation energy (&#x394;E<sub>SURF</sub>), as shown in Eq. <xref ref-type="disp-formula" rid="e4">(4)</xref>.</p>
</sec>
<sec id="s2-7">
<title>ADMET prediction</title>
<p>Following the stability studies conducted on the M3 and M5 molecules, their pharmacokinetic and pharmacodynamic properties were examined. Here, the first step involved confirming the feasibility of synthesizing these compounds, while the second step involved investigating the pharmaceutical properties of each compound. To evaluate the synthetic accessibility and pharmaceutical characteristics of the proposed compounds, the pkCSM (<xref ref-type="bibr" rid="B36">Pires et al., 2015</xref>) and SwissADME servers (<xref ref-type="bibr" rid="B13">Daina et al., 2017</xref>) were utilized.</p>
</sec>
</sec>
<sec sec-type="results|discussion" id="s3">
<title>Results and discussion</title>
<sec id="s3-1">
<title>Rigid molecular alignment method</title>
<p>Structural alignment of a molecule is key to the predictive precision of a 3D-QSAR model and accuracy of a contour model (<xref ref-type="bibr" rid="B21">Er-rajy et al., 2023c</xref>). The tetrahydrobenzo[4,5]thieno[2,4-day]pyrimidine and tetrahydrobenzo[4,5]thieno[2,3-day]pyrimidine molecules were first optimized and then aligned on the basis of their similarities with the template skeleton (3z) displayed in <xref ref-type="fig" rid="F2">Figure 2A</xref>. <xref ref-type="fig" rid="F2">Figure 2B</xref> demonstrates that all the examined molecular structures are appropriately aligned with the template structure. The shared core of all the overlaid molecules is depicted in <xref ref-type="fig" rid="F3">Figure 3</xref>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>
<bold>(A)</bold> Common structure; <bold>(B)</bold> two-dimensional template skeleton.</p>
</caption>
<graphic xlink:href="fchem-12-1384832-g002.tif"/>
</fig>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Molecular alignment.</p>
</caption>
<graphic xlink:href="fchem-12-1384832-g003.tif"/>
</fig>
</sec>
<sec id="s3-2">
<title>PLS analysis</title>
<p>PLS analysis was conducted to investigate the correlation between pIC<sub>50,</sub> and the 3D field descriptors obtained from the CoMFA and CoMSIA techniques. There are significant correlations between the inhibitory activities of the MCF-7 cells (measured as pIC<sub>50</sub>) and two descriptors: steric and electrostatic. The strong predictive abilities of our developed models are evident from the high R<sup>2</sup>
<sub>pred</sub> values [R<sup>2</sup>
<sub>pred</sub> (CoMFA) &#x3d; 0.90 and R<sup>2</sup>
<sub>pred</sub> (CoMSIA) &#x3d; 0.91] obtained from external validations (<xref ref-type="table" rid="T2">Table 2</xref>). This predictive power is further supported by the very low residuals between the observed pIC<sub>50</sub> values <italic>in vitro</italic> and predicted pIC<sub>50</sub> values <italic>in silico</italic> (<xref ref-type="table" rid="T3">Table 3</xref>).</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>CoMFA and CoMSIA statistics results based on PLS regression.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">&#x2003;</th>
<th align="center">Q<sup>2</sup>
</th>
<th align="center">
<italic>R</italic>
<sup>2</sup>
</th>
<th align="center">R<sup>2</sup>
<sub>pred</sub>
</th>
<th align="center">N</th>
<th align="center">SEE</th>
<th colspan="2" align="center">Fractions</th>
</tr>
<tr>
<th colspan="6" align="center">Metrics</th>
<th align="center">Steric</th>
<th align="center">Electrostatic</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">CoMFA</td>
<td align="center">0.62</td>
<td align="center">0.90</td>
<td align="center">0.90</td>
<td align="center">5</td>
<td align="center">0.24</td>
<td align="center">0.56</td>
<td align="center">0.43</td>
</tr>
<tr>
<td align="center">CoMSIA</td>
<td align="center">0.71</td>
<td align="center">0.88</td>
<td align="center">0.91</td>
<td align="center">9</td>
<td align="center">0.27</td>
<td align="center">0.55</td>
<td align="center">0.44</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>Q<sup>2</sup>, cross-validation correlation factor; N, optimum value of components determined by cross-validation (LOO-CV); SEE, standard error estimation; <italic>R</italic>
<sup>2</sup>, standard coefficient of determination; R<sup>2</sup>
<sub>pred</sub>, coefficient of determination based on the external test; fractions, participation of steric and electrostatic interactions.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>pIC<sub>50</sub> values observed and predicted by the 3D-QSAR (CoMFA and CoMSIA) models.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">Compound</th>
<th rowspan="2" align="center">MCF-7 pIC<sub>50</sub>
</th>
<th colspan="2" align="center">CoMFA</th>
<th colspan="2" align="center">CoMSIA</th>
</tr>
<tr>
<th align="center">Predicted (pIC<sub>50</sub>)</th>
<th align="center">Residual</th>
<th align="center">Predicted (pIC<sub>50</sub>)</th>
<th align="center">Residual</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">3a</td>
<td align="center">4.47</td>
<td align="center">4.42</td>
<td align="center">0.05</td>
<td align="center">4.44</td>
<td align="center">0.03</td>
</tr>
<tr>
<td align="center">3b<xref ref-type="table-fn" rid="Tfn1">
<sup>a</sup>
</xref>
</td>
<td align="center">4.92</td>
<td align="center">4.43</td>
<td align="center">0.49</td>
<td align="center">4.60</td>
<td align="center">0.32</td>
</tr>
<tr>
<td align="center">3c<xref ref-type="table-fn" rid="Tfn1">
<sup>a</sup>
</xref>
</td>
<td align="center">4.51</td>
<td align="center">4.41</td>
<td align="center">0.09</td>
<td align="center">4.69</td>
<td align="center">&#x2212;0.18</td>
</tr>
<tr>
<td align="center">3e</td>
<td align="center">4.70</td>
<td align="center">4.50</td>
<td align="center">0.19</td>
<td align="center">4.69</td>
<td align="center">0.016</td>
</tr>
<tr>
<td align="center">3f</td>
<td align="center">4.58</td>
<td align="center">4.41</td>
<td align="center">0.16</td>
<td align="center">4.54</td>
<td align="center">0.042</td>
</tr>
<tr>
<td align="center">3h<xref ref-type="table-fn" rid="Tfn1">
<sup>a</sup>
</xref>
</td>
<td align="center">4.60</td>
<td align="center">4.25</td>
<td align="center">0.34</td>
<td align="center">4.88</td>
<td align="center">&#x2212;0.28</td>
</tr>
<tr>
<td align="center">3i</td>
<td align="center">4.39</td>
<td align="center">4.60</td>
<td align="center">&#x2212;0.20</td>
<td align="center">4.47</td>
<td align="center">&#x2212;0.08</td>
</tr>
<tr>
<td align="center">3j</td>
<td align="center">5.23</td>
<td align="center">4.82</td>
<td align="center">0.41</td>
<td align="center">5.10</td>
<td align="center">0.13</td>
</tr>
<tr>
<td align="center">3k</td>
<td align="center">4.46</td>
<td align="center">4.75</td>
<td align="center">&#x2212;0.28</td>
<td align="center">4.45</td>
<td align="center">0.01</td>
</tr>
<tr>
<td align="center">3L</td>
<td align="center">4.92</td>
<td align="center">5.36</td>
<td align="center">&#x2212;0.43</td>
<td align="center">5.15</td>
<td align="center">&#x2212;0.23</td>
</tr>
<tr>
<td align="center">3m</td>
<td align="center">6.30</td>
<td align="center">6.38</td>
<td align="center">&#x2212;0.08</td>
<td align="center">6.30</td>
<td align="center">&#x2212;0.004</td>
</tr>
<tr>
<td align="center">3n<xref ref-type="table-fn" rid="Tfn1">
<sup>a</sup>
</xref>
</td>
<td align="center">5.40</td>
<td align="center">5.11</td>
<td align="center">0.29</td>
<td align="center">5.31</td>
<td align="center">0.096</td>
</tr>
<tr>
<td align="center">3o</td>
<td align="center">5.32</td>
<td align="center">5.13</td>
<td align="center">0.19</td>
<td align="center">5.38</td>
<td align="center">&#x2212;0.054</td>
</tr>
<tr>
<td align="center">3p</td>
<td align="center">5.05</td>
<td align="center">5.10</td>
<td align="center">&#x2212;0.05</td>
<td align="center">4.99</td>
<td align="center">0.058</td>
</tr>
<tr>
<td align="center">3q</td>
<td align="center">4.71</td>
<td align="center">5.14</td>
<td align="center">&#x2212;0.43</td>
<td align="center">4.73</td>
<td align="center">&#x2212;0.02</td>
</tr>
<tr>
<td align="center">3r<xref ref-type="table-fn" rid="Tfn1">
<sup>a</sup>
</xref>
</td>
<td align="center">6.69</td>
<td align="center">5.06</td>
<td align="center">1.63</td>
<td align="center">6.68</td>
<td align="center">0.01</td>
</tr>
<tr>
<td align="center">3s</td>
<td align="center">5.18</td>
<td align="center">5.09</td>
<td align="center">0.08</td>
<td align="center">5.17</td>
<td align="center">0.004</td>
</tr>
<tr>
<td align="center">3u</td>
<td align="center">4.80</td>
<td align="center">4.90</td>
<td align="center">&#x2212;0.09</td>
<td align="center">4.83</td>
<td align="center">&#x2212;0.02</td>
</tr>
<tr>
<td align="center">3v</td>
<td align="center">5.32</td>
<td align="center">5.32</td>
<td align="center">0.005</td>
<td align="center">5.26</td>
<td align="center">0.06</td>
</tr>
<tr>
<td align="center">3w</td>
<td align="center">5.72</td>
<td align="center">5.49</td>
<td align="center">0.22</td>
<td align="center">5.90</td>
<td align="center">&#x2212;0.18</td>
</tr>
<tr>
<td align="center">3x</td>
<td align="center">5.14</td>
<td align="center">5.27</td>
<td align="center">&#x2212;0.13</td>
<td align="center">5.17</td>
<td align="center">&#x2212;0.03</td>
</tr>
<tr>
<td align="center">3y</td>
<td align="center">6.22</td>
<td align="center">6.06</td>
<td align="center">0.15</td>
<td align="center">5.99</td>
<td align="center">0.22</td>
</tr>
<tr>
<td align="center">3z</td>
<td align="center">7</td>
<td align="center">6.78</td>
<td align="center">0.21</td>
<td align="center">6.97</td>
<td align="center">0.02</td>
</tr>
<tr>
<td align="center">4b</td>
<td align="center">5.06</td>
<td align="center">5.28</td>
<td align="center">&#x2212;0.22</td>
<td align="center">5.09</td>
<td align="center">&#x2212;0.03</td>
</tr>
<tr>
<td align="center">4c</td>
<td align="center">5.37</td>
<td align="center">5.21</td>
<td align="center">0.15</td>
<td align="center">5.35</td>
<td align="center">0.02</td>
</tr>
<tr>
<td align="center">4d</td>
<td align="center">5.09</td>
<td align="center">5.07</td>
<td align="center">0.02</td>
<td align="center">5.05</td>
<td align="center">0.04</td>
</tr>
<tr>
<td align="center">6a</td>
<td align="center">5.21</td>
<td align="center">5.06</td>
<td align="center">0.14</td>
<td align="center">5.12</td>
<td align="center">0.08</td>
</tr>
<tr>
<td align="center">6b</td>
<td align="center">5.61</td>
<td align="center">5.60</td>
<td align="center">0.01</td>
<td align="center">5.70</td>
<td align="center">&#x2212;0.08</td>
</tr>
<tr>
<td align="center">6c</td>
<td align="center">6.52</td>
<td align="center">6.62</td>
<td align="center">&#x2212;0.10</td>
<td align="center">6.53</td>
<td align="center">&#x2212;0.01</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn id="Tfn1">
<label>
<sup>a</sup>
</label>
<p>Test set.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>Based on the information in <xref ref-type="table" rid="T2">Table 2</xref>, it is observed that the electrostatic and steric forces that contribute 37% and 63%, respectively, to the overall variance represent the CoMFA descriptor. The Q<sup>2</sup> value obtained by adjusting the leave-one out (LOO) step in the PLS analysis is 0.62. Additionally, <italic>R</italic>
<sup>2</sup> is reported to be 0.90, with the number of components being five. This high <italic>R</italic>
<sup>2</sup> value also has strong external validation for the test set, with an R<sup>2</sup>
<sub>Test</sub> of 0.90. The CoMSIA descriptor is characterized by two fields, namely, the electrostatic and steric fields, which respectively contribute 44.6% and 55.4%. The Q<sup>2</sup> value adjusted through LOO-CV in the PLS method is 0.71. Additionally, <italic>R</italic>
<sup>2</sup> is evaluated to be 0.883, indicating a lower level of correlation. However, when considering the nine components, the high <italic>R</italic>
<sup>2</sup> value suggests robust external validation for the test set, with an R<sup>2</sup>
<sub>Test</sub> of 0.91.</p>
</sec>
<sec id="s3-3">
<title>Contour maps of CoMFA and CoMSIA</title>
<p>The contour maps in the QSAR models represent the relationships between the structures and activities of the series of chemicals. In this study, we analyzed the contour maps of the developed models to determine the factors influencing their anticancer activities (<xref ref-type="bibr" rid="B20">Er-rajy et al., 2023b</xref>). By conducting a three-dimensional analysis of these visualizations, we identified the regions in the template molecule that require modifications for designing new molecules with enhanced biological activities compared to the investigated derivatives. <xref ref-type="fig" rid="F4">Figure 4</xref> and <xref ref-type="fig" rid="F5">Figure 5</xref> depict the 3D contour maps of the template molecule structure generated using the CoMFA and CoMSIA models, respectively. The steric zone can be visualized using two colors; the presence of the green edges (which contribute to 80% of the zone) indicates regions where expansive clusters enhance the biological activity. On the other hand, the yellow areas (contribution amounting to 20%) represent regions where the voluminous clusters diminish the biological activity. The electrostatic zone can be visualized using two colors as well; the electrostatic field is deemed positive where the blue contours prevail (making up 80% of the zone), whereas the red edges (making up 20% of the zone) indicate the promotion of negative electrostatics. The CoMFA steric contours reveal a prominent green area surrounding the substituent in the Meta Position 2, suggesting that a bulky group is required for R substitution at the Meta position number 2 for the tetrahydrobenzo[4,5]thieno[2,3-day]pyrimidine derivative in this region to enhance the biological activity. The CoMFA electrostatic contours suggest that the introduction of strongly electropositive groups or atoms in this position may improve activity. In the CoMSIA electrostatic zone, we observed extensive yellow-colored contours surrounding the 9-methylacrinidin-3-amine substitution at the Meta position for the tetrahydrobenzo[4,5]thieno[2,3-day]pyrimidine derivatives. In addition, we noted red borders in the vicinities of the 9-methylacridin-3-amine substitutions in the Meta position number 3 of the compounds. This shows that the substitution of highly electronegative compounds or atoms can increase the activity of the tetrahydrobenzo[4,5]thieno[2,3-day]pyrimidine derivatives.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>
<bold>(A)</bold> CoMFA contour maps: electrostatic area interactions (blue &#x3d; promising/red &#x3d; unpromising) and steric area interactions (green &#x3d; promising/yellow &#x3d; unpromising). <bold>(B)</bold> CoMSIA contour maps: electrostatic area interactions (blue &#x3d; promising/red &#x3d; unpromising) and steric area interactions (green &#x3d; promising/yellow &#x3d; unpromising).</p>
</caption>
<graphic xlink:href="fchem-12-1384832-g004.tif"/>
</fig>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Mapping regions that are favorable for biological activity against the MCF-7 cell line.</p>
</caption>
<graphic xlink:href="fchem-12-1384832-g005.tif"/>
</fig>
</sec>
<sec id="s3-4">
<title>Design of new derivatives as inhibitors of MCF-7 cells</title>
<p>The ADMET analysis is crucial in the process of developing and discovering new drug compounds as it enables assessment of their pharmacological effects and potential risks more effectively. In our study, we specifically focused on investigating the effects of (3z), the most active compound, as well as the two newly designed compounds (M3 and M5), as presented in <xref ref-type="table" rid="T6">Table 6</xref>. Additionally, we examined the ADMET properties of capivasertib, an FDA-approved drug, to compare its results with those of the newly developed molecules. The ADMET findings are provided in <xref ref-type="table" rid="T4">Table 4</xref>. Here, 3z displayed a logS value of &#x2212;7.84, indicating significantly reduced solubility; the compound also exhibited an exceedingly high LogP value (7.088), suggesting its unfavorable hydrophobicity. In addition, 3z has low gastrointestinal (GI) absorption; this means that it is highly improbable for it to be successfully absorbed from the GI tract into the bloodstream, which can decrease the performance of a drug or limit the nutritional value of a food component. The AMES test (<italic>Salmonella</italic>/microsome) is a widely used assay to assess the potential of a chemical or substance to cause genetic mutations. Therefore, the AMES toxicity refers to the potential of a substance to be mutagenic based on its performance in the AMES test; it should be noted that this does not necessarily mean the substance is toxic in other ways, but it is highlighted for potential long-term health risks, particularly cancer 37. The compound 3z has been tested for its AMES toxicity, which means that it has been evaluated for its ability to induce genetic mutations in bacteria. From <xref ref-type="table" rid="T4">Table 4</xref>, although the molecule (3z) possesses favorable physicochemical properties, it is evident that it cannot be an orally active drug in humans according to the Lipinski rules. From these findings, it is concluded that molecule 3z does not meet the ADMET criteria. Based on the ADMET results, we modified the tetrahydrobenzo[4,5]thieno[2,3-day]pyrimidine structure by substituting new radicals at the Meta and Ortho positions while removing the ketone functional group to derive new compounds with good ADMET criteria.</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>ADMET analyses of the drug capivasertib, molecule object (3z), and newly designed molecules.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">ADMET/molecules</th>
<th align="center">Rule</th>
<th align="center">M3</th>
<th align="center">M5</th>
<th align="center">Capivasertib</th>
<th align="center">3z</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">
<bold>LogS</bold>
</td>
<td align="center">&#x2212;4 -0.5</td>
<td align="center">&#x2212;4.00</td>
<td align="center">&#x2212;4.50</td>
<td align="center">&#x2212;3.48</td>
<td align="center">&#x2212;7.84</td>
</tr>
<tr>
<td align="center">
<bold>LogP</bold>
</td>
<td align="center">0&#x2013;3</td>
<td align="center">2.691</td>
<td align="center">3.54</td>
<td align="center">2.14</td>
<td align="center">7.088</td>
</tr>
<tr>
<td align="center">
<bold>BBB permeability</bold>
</td>
<td align="left"/>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
</tr>
<tr>
<td align="center">
<bold>GI absorption</bold>
</td>
<td align="left"/>
<td align="center">High</td>
<td align="center">High</td>
<td align="center">High</td>
<td align="center">Low</td>
</tr>
<tr>
<td align="center">
<bold>D6-sub</bold>
</td>
<td align="left"/>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
</tr>
<tr>
<td align="center">
<bold>A4-sub</bold>
</td>
<td align="left"/>
<td align="center">Yes</td>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">Yes</td>
</tr>
<tr>
<td align="center">
<bold>A2-inh</bold>
</td>
<td align="left"/>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">Yes</td>
</tr>
<tr>
<td align="center">
<bold>C19-inh</bold>
</td>
<td align="left"/>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
</tr>
<tr>
<td align="center">
<bold>C9-inh</bold>
</td>
<td align="left"/>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
</tr>
<tr>
<td align="center">
<bold>D6-inh</bold>
</td>
<td align="left"/>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
</tr>
<tr>
<td align="center">
<bold>A4-inh</bold>
</td>
<td align="left"/>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
</tr>
<tr>
<td align="center">
<bold>HERG I-inh</bold>
</td>
<td align="left"/>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
</tr>
<tr>
<td align="center">
<bold>HEGR II-inh</bold>
</td>
<td align="left"/>
<td align="center">Yes</td>
<td align="center">No</td>
<td align="center">Yes</td>
<td align="center">No</td>
</tr>
<tr>
<td align="center">
<bold>Oral rat acute toxicity (LD</bold>
<sub>
<bold>50</bold>
</sub>
<bold>)</bold>
</td>
<td align="left"/>
<td align="center">2.061</td>
<td align="center">2.607</td>
<td align="center">2.305</td>
<td align="center">2.725</td>
</tr>
<tr>
<td align="center">
<bold>Oral rat chronic toxicity (LOAEL)</bold>
</td>
<td align="left"/>
<td align="center">1.258</td>
<td align="center">0.759</td>
<td align="center">0.643</td>
<td align="center">&#x2212;0.795</td>
</tr>
<tr>
<td align="center">
<bold>Skin sensitization</bold>
</td>
<td align="left"/>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
</tr>
<tr>
<td align="center">
<bold>AMES toxicity</bold>
</td>
<td align="left"/>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">No</td>
<td align="center">Yes</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>
<xref ref-type="table" rid="T5">Table 5</xref> represents the six new molecules that were designed to have high biological activity predictions (pIC<sub>50</sub>). Among these molecules, we selected M3 and M5 for further investigations based on the ADMET studies by considering their strongly predicted biological activities. We observed that compound M3 displayed a satisfactory logS value. However, neither of the newly designed molecules were able to penetrate the blood&#x2013;brain barrier (BBB) effectively. Additionally, both compounds exhibited high GI absorption. The M5 molecule and capivasertib served as substrates for Enzyme A4, while the compound M3 was not used. In addition to functioning as substrates, both M5 and capivasertib act as inhibitors of the HEGR II enzyme. According to the findings presented in <xref ref-type="table" rid="T5">Table 5</xref>, all compounds comply with the established guidelines for hydrogen-bond acceptor (nHA; ranging from 3 to 5), hydrogen-bond donor (nHD; ranging from 1 to 4), topological surface area (TPSA &#x3c; 140&#xa0;&#xc5;<sup>2</sup>), and molecular weight (ranging from 428.92 to 514.64&#xa0;g/mol). Furthermore, the two compounds adhere to Lipinski&#x2019;s criteria for potential oral bioavailability. In conclusion, the newly designed compounds satisfy the requirements for druglikeness and are promising as potential candidates for the development of novel oral medications.</p>
<table-wrap id="T5" position="float">
<label>TABLE 5</label>
<caption>
<p>Physicochemical characteristics and Lipinski&#x2019;s rule for the new compounds.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">P.P</th>
<th align="center">nHA</th>
<th align="center">nHD</th>
<th align="center">TPSA</th>
<th align="center">MW</th>
<th align="center">MR</th>
<th align="center">Synthetic accessibility</th>
<th align="center">Lipinski</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Rule</td>
<td align="center">0&#x2013;12</td>
<td align="center">0&#x2013;7</td>
<td align="center">0&#x2013;14</td>
<td align="center">100&#x2013;600</td>
<td align="left"/>
<td align="left"/>
<td align="left"/>
</tr>
<tr>
<td align="center">M3</td>
<td align="center">5</td>
<td align="center">4</td>
<td align="center">130.39</td>
<td align="center">446.52</td>
<td align="center">133.79</td>
<td align="center">4.66</td>
<td align="center">Yes</td>
</tr>
<tr>
<td align="center">M5</td>
<td align="center">5</td>
<td align="center">3</td>
<td align="center">121.60</td>
<td align="center">472.56</td>
<td align="center">143.02</td>
<td align="center">4.85</td>
<td align="center">Yes</td>
</tr>
<tr>
<td align="center">Capivasertib</td>
<td align="center">5</td>
<td align="center">4</td>
<td align="center">120.16</td>
<td align="center">428.92</td>
<td align="center">118.97</td>
<td align="center">-</td>
<td align="left"/>
</tr>
<tr>
<td align="center">3z</td>
<td align="center">3</td>
<td align="center">1</td>
<td align="center">102.04</td>
<td align="center">514.64</td>
<td align="center">160.20</td>
<td align="center">-</td>
<td align="center">No</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-5">
<title>Design of new molecules targeting the MCF-7 cancer cell line</title>
<p>We utilized 3D-QSAR modeling and conducted ADMET studies to design new compounds that demonstrate improved predicted biological activities. Furthermore, our analysis of the ADMET properties and physicochemical characteristics shows promising results, suggesting that the new compounds possess favorable pharmacokinetic properties (<xref ref-type="table" rid="T6">Table 6</xref>). Consequently, we believe that the designed molecules have the potential to serve as inhibitors of the MCF-7 cell line and may be considered as candidate drugs.</p>
<table-wrap id="T6" position="float">
<label>TABLE 6</label>
<caption>
<p>Design of new inhibitors against the MCF-7 cell line.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Name</th>
<th align="center">Structure</th>
<th align="center">pIC<sub>50</sub> predicted by 3D-QSAR</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">
<bold>M1</bold>
</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx30.tif"/>
</td>
<td align="center">6.21</td>
</tr>
<tr>
<td align="center">
<bold>M2</bold>
</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx31.tif"/>
</td>
<td align="center">6.14</td>
</tr>
<tr>
<td align="center">
<bold>M3</bold>
</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx32.tif"/>
</td>
<td align="center">6.21</td>
</tr>
<tr>
<td align="center">
<bold>M4</bold>
</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx33.tif"/>
</td>
<td align="center">6.18</td>
</tr>
<tr>
<td align="center">
<bold>M5</bold>
</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx34.tif"/>
</td>
<td align="center">6.18</td>
</tr>
<tr>
<td align="center">
<bold>M6</bold>
</td>
<td align="center">
<inline-graphic xlink:href="FCHEM_fchem-2024-1384832_wc_tfx35.tif"/>
</td>
<td align="center">6.18</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-6">
<title>Molecular docking studies</title>
<p>Molecular docking is based on algorithms that anticipate the most stable conformation of the ligand&#x2013;protein complex depending on the binding energy. We performed molecular docking studies for the most active compounds and capivasertib, which is an FDA-approved drug, to compare the docking interactions and performances between the new compounds with the receptor and drug capivasertib having the same receptor. We performed docking simulations to examine the interactions between our predicted molecules and the protein ER&#x3b1; (PDB ID: 4XO6) <xref ref-type="fig" rid="F6">Figure 6</xref>. The findings of the simulations for the two selected compounds (M3 and M5) and capivasertib are shown in <xref ref-type="table" rid="T7">Figure 7A&#x2013;C</xref>, respectively. The initial observations visually demonstrate that M3 has three hydrogen bonds involving residues Leu219, His222, and Ser217; this is attributed to the polarity of OH- and presence of a lone electron pair on the oxygen, which act as electron acceptors during formation of hydrogen bonds. However, M5 exhibits a distinct feature wherein it forms an extra hydrogen bond with residue Leu268; moreover, it engages in alkyl and &#x3c0;-alkyl interactions with Trp86, Leu308, Peh118, and Lys270. These interactions occur when the alkyl group of the compound interacts with a pi system within the receptor due to the hydrophobic nature of the cyclohexane and cyclopenta-1, 3-diene moieties. The presence of the penta-1, 3-diene aromatic group in the ligands (M3 and M5 molecules) leads to &#x3c0;-sigma interactions with a positively charged residue (Leu306), which affects the stability of the ligand&#x2013;molecule complex. These interactions are significant as they can greatly influence the binding affinities and specificities of the ligands.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Topology of the ER&#x3b1; protein (PDB ID: 4XO6, chain A).</p>
</caption>
<graphic xlink:href="fchem-12-1384832-g006.tif"/>
</fig>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>2D and 3D visualizations of the docking results for the newly designed compounds <bold>(A)</bold> M3 and <bold>(B)</bold> M5 against the MCF-7 breast cancer cell line. <bold>(C)</bold> Reference 2D and 3D visualizations of the docking results for the FDA-approved drug capivasertib.</p>
</caption>
<graphic xlink:href="fchem-12-1384832-g007.tif"/>
</fig>
<p>Other interactions are also noted, like the &#x3c0;-sulfur case between the thiophene and residues Asn167 and His117; these interactions occur when the delocalized &#x3c0;-orbitals of the thiophene ring interact with the lone electron pairs of the sulfur atom in the &#x3c0;-sulfur ligand. This may contribute to the affinity of the bond in the docking results by stabilizing the complex through orbital riding and electron delocalization.</p>
<p>The docking results of capivasertib shown in <xref ref-type="fig" rid="F7">Figure 7C</xref> and reveal four hydrogen bonds involving specific residues and functional groups. One hydrogen bond is formed between the residue Ser50 and OH functional group of capivasertib. Another hydrogen bond is formed between the residue Ser217 and Cl atom of capivasertib. Additionally, two hydrogen bonds are observed with the Tyr55 residue, one involving the OH functional group and the other involving the ketone functional group of capivasertib. These hydrogen bonds play crucial roles in the interactions between capivasertib and its target molecules. We observed that the interaction patterns between our newly designed drugs (M3/M5) and capivasertib are identical. Specifically, we note the presence of two interactions: a hydrogen bond with Ser217 and a &#x3c0;-&#x3c0; interaction involving Tyr24. By comparing the results obtained from M3 and M5 with the docking results of capivasertib, we deduce that both M5 and capivasertib interact with two similar residues, namely, Ser217 and Tyr24. On the other hand, M3 only interacts with one residue, which is Ser217. This suggests that the compounds have the potential to bind effectively with the target biomolecules. According to literature, the residues Leu219, His222, and Ser217 also play very important roles in studies on anticancer activities (<xref ref-type="bibr" rid="B48">Zouaghi et al., 2023</xref>). This is a positive sign for drug development because successful binding can influence the desired therapeutic effect.</p>
</sec>
<sec id="s3-7">
<title>MD simulation</title>
<p>MD simulations were conducted to comprehensively investigate the stabilities of the newly designed molecules in a complex with the protein. Such simulations have been used in many studies to explain atomic-level variations within the protein&#x2013;ligand complex, capturing the dynamic interactions between the molecules in a dynamic medium (<xref ref-type="bibr" rid="B32">Kumari et al., 2022</xref>). In the following study, we examined several parameters such as the RMSD, RMSF, RoG, SASA, and hydrogen-bond formation for the M3 and M5 complexes during 100 ns of molecular simulation.</p>
<p>The results achieved for the RMSD, RMSF, RoG, SASA, and H-bonds are presented in <xref ref-type="fig" rid="F8">Figure 8</xref>. The RMSD investigation for the D3M (M3&#x2013;ER&#x3b1;) and D5M (M5&#x2013;ER&#x3b1;) complexes show average distances between 1 and 2.3&#xa0;&#xc5;. Both complexes vary from 1 to 2.3&#xa0;&#xc5; and thereafter stabilize at 2&#xa0;&#xc5; during the first 20&#xa0;ns. The RMSD analysis demonstrates that D3M and D5M are more effective in promoting stability than capivasertib. However, some fluctuations are observed in the RMSD plot, such as a fluctuation in the D3M complex from 1.5 to 1.2&#xa0;&#xc5; at approximately 40&#xa0;ns. The D5M complex demonstrates strong stability over a period of 60&#xa0;ns, with low fluctuations from 1 to 2.3&#xa0;&#xc5; in the first 40&#xa0;ns, after which it remains stable at 2&#xa0;&#xc5;. In contrast, the reference molecule capivasertib experiences fluctuations from 1 to 2.8&#xa0;&#xc5; over a period of 40&#xa0;ns.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>RMSD, RMSF, SASA, RoG, and H-bonds of the D3M and D5M complexes, as well as capivasertib.</p>
</caption>
<graphic xlink:href="fchem-12-1384832-g008.tif"/>
</fig>
<p>Furthermore, the complex capivasertib&#x2013;ER&#x3b1; exhibits an increase in fluctuation from 2.7&#xa0;&#xc5; to 3.8&#xa0;&#xc5; between 45 ns and 75&#xa0;ns. Good stabilities of the D3M and D5M complexes are observed from 80 ns to 100 ns at 1.8&#xa0;&#xc5;, while the reference molecule capivasertib experiences a range of fluctuations from 1 to 3.8&#xa0;&#xc5; over a period of 100&#xa0;ns. The RMSD plots show that the compounds M3 and M5 exhibit good stabilities with the ER&#x3b1; protein compared to capivasertib.</p>
<p>The RMSF values were calculated to assess the flexibility of the protein structure after adjustment to the reference frame. Higher RMSF values suggest that the protein&#x2019;s loops and turns are poorly structured, whereas lower flexibility indicates the presence of more stable and organized regions. The RMSF analysis revealed similar stabilities for the residues of the complexes formed by the new compounds with ER&#x3b1;. <xref ref-type="table" rid="T7">Table 7</xref> presents the residues with RMSF values less than 3&#xa0;&#xc5; that contribute to the stability of the newly designed molecules. In particular, the Glu224 residue in capivasertib exhibits an RMSF value of 4.557291, followed by Glu227 with an RMSF value of 8.392049. These values explain the peaks observed in <xref ref-type="fig" rid="F8">Figure 8</xref>, which are located outside the active sites of ER&#x3b1; compared with D3M and D5M.</p>
<table-wrap id="T7" position="float">
<label>TABLE 7</label>
<caption>
<p>Residues with RMSF values less than 3 &#xc5;.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">N&#xb0;</th>
<th align="center">Residue</th>
<th align="center">D3M</th>
<th align="center">D5M</th>
<th align="center">Capivasertib</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">14</td>
<td align="left">ASP</td>
<td align="left">0.565390</td>
<td align="left">0.727055</td>
<td align="left">0.781049</td>
</tr>
<tr>
<td align="left">29</td>
<td align="left">SER</td>
<td align="left">1.496763</td>
<td align="left">1.682218</td>
<td align="left">1.531807</td>
</tr>
<tr>
<td align="left">44</td>
<td align="left">TYR</td>
<td align="left">0.700645</td>
<td align="left">0.682867</td>
<td align="left">0.678045</td>
</tr>
<tr>
<td align="left">59</td>
<td align="left">GLN</td>
<td align="left">0.932650</td>
<td align="left">0.799256</td>
<td align="left">0.681637</td>
</tr>
<tr>
<td align="left">74</td>
<td align="left">VAL</td>
<td align="left">0.801856</td>
<td align="left">0.731317</td>
<td align="left">0.909452</td>
</tr>
<tr>
<td align="left">98</td>
<td align="left">LEU</td>
<td align="left">0.738380373</td>
<td align="left">0.532187521</td>
<td align="left">0.5297966</td>
</tr>
<tr>
<td align="left">104</td>
<td align="left">ASN</td>
<td align="left">0.791040</td>
<td align="left">0.613014</td>
<td align="left">0.666166</td>
</tr>
<tr>
<td align="left">119</td>
<td align="left">HIS</td>
<td align="left">0.578250</td>
<td align="left">0.477338</td>
<td align="left">0.595151</td>
</tr>
<tr>
<td align="left">134</td>
<td align="left">PHE</td>
<td align="left">3.117864</td>
<td align="left">2.372437</td>
<td align="left">2.705154</td>
</tr>
<tr>
<td align="left">149</td>
<td align="left">PRO</td>
<td align="left">0.519812</td>
<td align="left">0.525314</td>
<td align="left">0.602962</td>
</tr>
<tr>
<td align="left">179</td>
<td align="left">PHE</td>
<td align="left">1.326364</td>
<td align="left">0.730985</td>
<td align="left">0.893380</td>
</tr>
<tr>
<td align="left">209</td>
<td align="left">ALA</td>
<td align="left">0.644678</td>
<td align="left">0.767425</td>
<td align="left">0.736610</td>
</tr>
<tr>
<td align="left">224</td>
<td align="left">GLU</td>
<td align="left">2.018936</td>
<td align="left">1.457529</td>
<td align="left">&#x3e;3</td>
</tr>
<tr>
<td align="left">239</td>
<td align="left">PRO</td>
<td align="left">1.502835</td>
<td align="left">1.847240</td>
<td align="left">1.512882</td>
</tr>
<tr>
<td align="left">254</td>
<td align="left">LYS</td>
<td align="left">0.813202</td>
<td align="left">0.625272</td>
<td align="left">0.764170</td>
</tr>
<tr>
<td align="left">269</td>
<td align="left">LYS</td>
<td align="left">0.524341</td>
<td align="left">0.672288</td>
<td align="left">0.597879</td>
</tr>
<tr>
<td align="left">284</td>
<td align="left">AGLU</td>
<td align="left">1.189149</td>
<td align="left">0.878250</td>
<td align="left">0.889363</td>
</tr>
<tr>
<td align="left">299</td>
<td align="left">ALA</td>
<td align="left">0.910628</td>
<td align="left">1.563716</td>
<td align="left">1.008637</td>
</tr>
<tr>
<td align="left">314</td>
<td align="left">LYS</td>
<td align="left">1.055927</td>
<td align="left">0.734300</td>
<td align="left">1.826706</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The plots of SASA and RoG exhibit that the compounds M3 and M5 maintain their structure stabilities during the 100 ns simulation; both molecules are observed to be stable on the SASA plot, with an average value of 16,800&#xa0;&#xc5;<sup>2</sup>. This pattern is also observed in the RoG calculation, which yields an average value of 19.8&#xa0;&#xc5;. However, in the case of capivasertib, the SASA calculation shows an average value of 18,300&#xa0;&#xc5;<sup>2</sup> with noticeable fluctuations; this remarkable variation significantly affects both the molecular surface of capivasertib and its interactions with the solvent. However, the RoG analysis reveals a fluctuation of 20&#xa0;&#xc5; initially, followed by stabilization at 19&#xa0;&#xc5;; this suggests that capivasertib is unable to maintain its compactness throughout the simulation. The data presented in <xref ref-type="fig" rid="F8">Figure 8</xref> demonstrate the favorable stabilities of the n-hydrogen bonds of the newly designed molecules as they do not exceed four hydrogen bonds over 100&#xa0;ns. Importantly, these bonds do not show any unnecessary or additional occupancies. However, capivasertib exhibits significant variation, particularly between 15&#xa0;ns and 25&#xa0;ns, where the count of hydrogen bonds reaches ten, allowing additional occupancies. Additional variations are observed between 30&#xa0;ns and 40&#xa0;ns and between 60&#xa0;ns and 70&#xa0;ns; these observations are further supported by the results of the RMSD, RMSF, SASA, and RoG analyses.</p>
<p>The MD simulation results of the RMSD, RMSF, SASA, and RoG calculations substantiate the strong affinities of the newly designed models towards ER&#x3b1;. Additionally, these simulations confirm the remarkable stabilities of the models, as evidenced by their consistent behaviors throughout the 100 ns of simulation. Essential dynamics analysis was performed to examine the significant movements that occur after a ligand is bound to the protein. By comparing the images at 0 ns, 50 ns, and 100&#xa0;ns, the dynamic differences among M3, M5, and capivasertib can be visualized easily.</p>
<p>Based on the 2D projection map in <xref ref-type="fig" rid="F9">Figure 9</xref>, it is clear that inhibitor binding leads to formation of stable complexes. Among the complexes, M3 is the most stable, as evidenced by its highly stable cluster and minimal occupation of the phase space compared with the other molecules. In contrast, capivasertib has the lowest stability out of all the compounds identified. This observation is further confirmed by the representative illustrations.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Mappings of the conformations of the M3, M5, and capivasertib complexes at 0&#xa0;ns, 50&#xa0;ns, and 100&#xa0;ns from the MD simulations.</p>
</caption>
<graphic xlink:href="fchem-12-1384832-g009.tif"/>
</fig>
</sec>
<sec id="s3-8">
<title>Free binding energy (MM/GBSA)</title>
<p>We input snapshots from the MD results obtained during a 100-ns trajectory to the MM-GBSA method. From <xref ref-type="table" rid="T8">Table 8</xref>, it is evident that &#x394;E<sub>EL</sub> for capivasertib is larger than the values of the new compounds. This leads us to the conclusion that capivasertib is relatively less stable compared to D3M and D5M. The newly designed compounds exhibit higher values of van der Waals interactions than capivasertib. This indicates that the strengths of the van der Waals interactions in D3M and D5M are larger. The negative energy surface interaction, &#x394;E<sub>SUR</sub>, typically indicates an attractive force. The D3M complex exhibits favorable energy surface interaction; the greater the negative energy of a compound in the gas phase, the better its gas-phase power. All compounds demonstrate favorable energies in the gas phase. Specifically, capivasertib exhibits a strong gas-phase interaction owing to its &#x394;E<sub>SURF</sub> value. Positive values of energy solvation, &#x394;E<sub>Sol</sub>, enhance the stabilities of the compounds; the new compounds exhibit positive energies when combined with the complexes D3M and D5M, unlike capivasertib that has a lower &#x394;E<sub>Sol</sub> value. This indicates that the new compounds possess lower solvation energies. The &#x394;E<sub>GB</sub> results show that the powers of all compounds are involved in cavitation/solvation; here, D5M exhibits a lower &#x394;E<sub>GB</sub> than D3M and capivasertib, which implies that D5M may have a decreased ability to solvate that in turn enhances its binding efficiency. Although capivasertib&#x2019;s total energy is positive, both D3M and D5M exhibit negative total energies. This suggests that D3M and D5M have more advantageous binding affinities than capivasertib, indicating their superior efficacies as ER&#x3b1; inhibitors in MCF-7 cells. The MM/GBSA results suggest that the newly designed compounds M3 and M5 show more beneficial interactions with the target protein than capivasertib. Molecule M3 exhibits highly promising outcomes, as evidenced by its favorable total energy and superior binding energy in comparison to the other compounds.</p>
<table-wrap id="T8" position="float">
<label>TABLE 8</label>
<caption>
<p>Free energy binding values of the complexes based on the MM/GBSA method.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">&#x394;E (kcal&#x2219;mol<sup>-1</sup>)</th>
<th align="center">DM3</th>
<th align="center">DM5</th>
<th align="center">Capivasertib</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">&#x394;E<sub>GB</sub>
</td>
<td align="center">291.88</td>
<td align="center">129.60</td>
<td align="center">393.17</td>
</tr>
<tr>
<td align="center">&#x394;E<sub>SURF</sub>
</td>
<td align="center">&#x2212;7.34</td>
<td align="center">&#x2212;5.93</td>
<td align="center">&#x2212;4.44</td>
</tr>
<tr>
<td align="center">&#x394;E<sub>EL</sub>
</td>
<td align="center">&#x2212;294.01</td>
<td align="center">&#x2212;98.85</td>
<td align="center">347.08</td>
</tr>
<tr>
<td align="center">&#x394;E<sub>Gas</sub>
</td>
<td align="center">&#x2212;337.84</td>
<td align="center">&#x2212;141.88</td>
<td align="center">&#x2212;367.25</td>
</tr>
<tr>
<td align="center">&#x394;E<sub>Sol</sub>
</td>
<td align="center">284.53</td>
<td align="center">123.66</td>
<td align="center">388.73</td>
</tr>
<tr>
<td align="center">&#x394;<sub>VWAALS</sub>
</td>
<td align="center">&#x2212;43.83</td>
<td align="center">&#x2212;43.02</td>
<td align="center">&#x2212;20.17</td>
</tr>
<tr>
<td align="center">&#x394;E<sub>Total</sub>
</td>
<td align="center">&#x2212;53.31</td>
<td align="center">&#x2212;18.21</td>
<td align="center">21.48</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec sec-type="conclusion" id="s4">
<title>Conclusion</title>
<p>In11 this study, new potential drugs were investigated using molecular modeling techniques against breast cancer, which is one of the most prevalent diseases worldwide. The primary aim of the study was to identify MCF-7 inhibitors, which are crucial cell lines in the treatment of breast cancer. The 3D-QSAR method was employed to predict new MCF-7 inhibitors, resulting in the generation of two strong models with good results in identifying groups that may interact with each other. The two established models were evaluated successfully using external validation methods. Analysis of the various contour maps of both models allowed us to identify the functional groups that could enhance the inhibitory abilities. We utilized ADMET studies and pharmacokinetics characteristics to design six inhibitors based on their CoMFA and CoMSIA predictions, with the addition of OH functional groups. The presence and position of the OH group can affect the binding affinity, selectivity, and efficacy of these molecules at the target sites, allowing strong intermolecular interactions with other polar molecules. Hydrogen bonding is critical for many biological processes such as protein folding. Out of the six inhibitors developed, two molecules exhibited good interactions with the target protein in molecular docking simulations. To validate these findings, the data were subjected to MD simulations and binding free energy calculations based on the MM/GBSA over a duration of 100&#xa0;ns. In conclusion, this research provides important insights into the development of novel drugs to combat breast cancer. The results of this study are expected to be of foundational value for further research and development of inhibitors as promising drugs for breast cancer treatment.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s5">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s6">
<title>Ethics statement</title>
<p>Ethical approval was not required for this study involving animals in accordance with the local legislation and institutional requirements because this study was conducted using computer programs to clarify some compounds and their activities on cancer cells.</p>
</sec>
<sec id="s7">
<title>Author contributions</title>
<p>SZ: conceptualization, data curation, formal analysis, methodology, software, writing&#x2013;original draft, and writing&#x2013;review and editing. ME-R: conceptualization, formal analysis, methodology, software, and writing&#x2013;original draft. AF: formal analysis, methodology, validation, and writing&#x2013;review and editing. HI: data curation and writing&#x2013;review and editing. ME: formal analysis, software, and writing&#x2013;review and editing. AQ: resources and writing&#x2013;review and editing. FN: resources, and writing&#x2013;review and editing. MA-Z: resources, validation, and writing&#x2013;review and editing. ME: supervision, validation, and writing&#x2013;review and editing.</p>
</sec>
<sec id="s8">
<title>Funding</title>
<p>The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported and funded by the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University (IMSIU) (grant number IMSIU-RG23122).</p>
</sec>
<ack>
<p>The authors are thankful for the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University (IMSIU) (grant number IMSIU-RG23122).</p>
</ack>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abbas</surname>
<given-names>S. E.</given-names>
</name>
<name>
<surname>Abdel Gawad</surname>
<given-names>N. M.</given-names>
</name>
<name>
<surname>George</surname>
<given-names>R. F.</given-names>
</name>
<name>
<surname>Akar</surname>
<given-names>Y. A.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Synthesis, antitumor and antibacterial activities of some novel tetrahydrobenzo[4,5]thieno[2,3-d]pyrimidine derivatives</article-title>. <source>Eur. J. Med. Chem.</source> <volume>65</volume>, <fpage>195</fpage>&#x2013;<lpage>204</lpage>. <pub-id pub-id-type="doi">10.1016/j.ejmech.2013.04.055</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abraham</surname>
<given-names>M. J.</given-names>
</name>
<name>
<surname>Murtola</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Schulz</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>P&#xe1;ll</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Smith</surname>
<given-names>J. C.</given-names>
</name>
<name>
<surname>Hess</surname>
<given-names>B.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers</article-title>. <source>SoftwareX</source> <volume>1</volume>, <fpage>19</fpage>&#x2013;<lpage>25</lpage>. <pub-id pub-id-type="doi">10.1016/j.softx.2015.06.001</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alagarsamy</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Meena</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Ramseshu</surname>
<given-names>K. V.</given-names>
</name>
<name>
<surname>Solomon</surname>
<given-names>V. R.</given-names>
</name>
<name>
<surname>Thirumurugan</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Dhanabal</surname>
<given-names>K.</given-names>
</name>
<etal/>
</person-group> (<year>2006</year>). <article-title>Synthesis, analgesic, anti-inflammatory, ulcerogenic index and antibacterial activities of novel 2-methylthio-3-substituted-5,6,7,8-tetrahydrobenzo (b) thieno[2,3-d]pyrimidin-4(3H)-ones</article-title>. <source>Eur. J. Med. Chem.</source> <volume>41</volume>, <fpage>1293</fpage>&#x2013;<lpage>1300</lpage>. <pub-id pub-id-type="doi">10.1016/j.ejmech.2006.06.005</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alqahtani</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>
<italic>In silico</italic> ADME-Tox modeling: progress and prospects</article-title>. <source>Expert Opin. Drug Metab. Toxicol.</source> <volume>13</volume>, <fpage>1147</fpage>&#x2013;<lpage>1158</lpage>. <pub-id pub-id-type="doi">10.1080/17425255.2017.1389897</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aly</surname>
<given-names>H. M.</given-names>
</name>
<name>
<surname>Saleh</surname>
<given-names>N. M.</given-names>
</name>
<name>
<surname>Elhady</surname>
<given-names>H. A.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Design and synthesis of some new thiophene, thienopyrimidine and thienothiadiazine derivatives of antipyrine as potential antimicrobial agents</article-title>. <source>Eur. J. Med. Chem.</source> <volume>46</volume>, <fpage>4566</fpage>&#x2013;<lpage>4572</lpage>. <pub-id pub-id-type="doi">10.1016/j.ejmech.2011.07.035</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Amr</surname>
<given-names>A. E.-G. E.</given-names>
</name>
<name>
<surname>Sherif</surname>
<given-names>M. H.</given-names>
</name>
<name>
<surname>Assy</surname>
<given-names>M. G.</given-names>
</name>
<name>
<surname>Al-Omar</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Ragab</surname>
<given-names>I.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Antiarrhythmic, serotonin antagonist and antianxiety activities of novel substituted thiophene derivatives synthesized from 2-amino-4,5,6,7-tetrahydro-N-phenylbenzo[b]thiophene-3-carboxamide</article-title>. <source>Eur. J. Med. Chem.</source> <volume>45</volume>, <fpage>5935</fpage>&#x2013;<lpage>5942</lpage>. <pub-id pub-id-type="doi">10.1016/j.ejmech.2010.09.059</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aparoy</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Suresh</surname>
<given-names>G. K.</given-names>
</name>
<name>
<surname>Kumar Reddy</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Reddanna</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>CoMFA and CoMSIA studies on 5-hydroxyindole-3-carboxylate derivatives as 5-lipoxygenase inhibitors: generation of homology model and docking studies</article-title>. <source>Bioorg. Med. Chem. Lett.</source> <volume>21</volume>, <fpage>456</fpage>&#x2013;<lpage>462</lpage>. <pub-id pub-id-type="doi">10.1016/j.bmcl.2010.10.119</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aponte</surname>
<given-names>J. C.</given-names>
</name>
<name>
<surname>Vaisberg</surname>
<given-names>A. J.</given-names>
</name>
<name>
<surname>Castillo</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Gonzalez</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Estevez</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Arevalo</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2010</year>). <article-title>Trypanoside, anti-tuberculosis, leishmanicidal, and cytotoxic activities of tetrahydrobenzothienopyrimidines</article-title>. <source>Bioorg. Med. Chem.</source> <volume>18</volume>, <fpage>2880</fpage>&#x2013;<lpage>2886</lpage>. <pub-id pub-id-type="doi">10.1016/j.bmc.2010.03.018</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aqvist</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>1990</year>). <article-title>Ion-water interaction potentials derived from free energy perturbation simulations</article-title>. <source>J. Phys. Chem.</source> <volume>94</volume>, <fpage>8021</fpage>&#x2013;<lpage>8024</lpage>. <pub-id pub-id-type="doi">10.1021/j100384a009</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Case</surname>
<given-names>D. A.</given-names>
</name>
<name>
<surname>Aktulga</surname>
<given-names>H. M.</given-names>
</name>
<name>
<surname>Belfon</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Cerutti</surname>
<given-names>D. S.</given-names>
</name>
<name>
<surname>Cisneros</surname>
<given-names>G. A.</given-names>
</name>
<name>
<surname>Cruzeiro</surname>
<given-names>V. W. D.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>AmberTools</article-title>. <source>J. Chem. Inf. Model.</source> <volume>63</volume>, <fpage>6183</fpage>&#x2013;<lpage>6191</lpage>. <pub-id pub-id-type="doi">10.1021/acs.jcim.3c01153</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="book">
<collab>Certara</collab> (<year>2017</year>). <source>SYBYL-X software</source>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.certara.com/sybyl-x-software/">https://www.certara.com/sybyl-x-software/</ext-link>(Accessed December 28, 2023)</comment>.</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dai</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Frey</surname>
<given-names>R. R.</given-names>
</name>
<name>
<surname>Ji</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Curtin</surname>
<given-names>M. L.</given-names>
</name>
<name>
<surname>Ahmed</surname>
<given-names>A. A.</given-names>
</name>
<etal/>
</person-group> (<year>2005</year>). <article-title>Thienopyrimidine ureas as novel and potent multitargeted receptor tyrosine kinase inhibitors</article-title>. <source>J. Med. Chem.</source> <volume>48</volume>, <fpage>6066</fpage>&#x2013;<lpage>6083</lpage>. <pub-id pub-id-type="doi">10.1021/jm050458h</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Daina</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Michielin</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Zoete</surname>
<given-names>V.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules</article-title>. <source>Sci. Rep.</source> <volume>7</volume>, <fpage>42717</fpage>. <pub-id pub-id-type="doi">10.1038/srep42717</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="book">
<collab>Dassault Syst&#xe8;mes</collab> (<year>2024</year>). <source>BIOVIA Discovery studio visualizer</source>. <publisher-name>Dassault Syst&#xe8;mes</publisher-name>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.3ds.com/products/biovia/discovery-studio/visualization">https://www.3ds.com/products/biovia/discovery-studio/visualization</ext-link> (Accessed April 9, 2024)</comment>.</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Deng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Lan</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>F.</given-names>
</name>
<etal/>
</person-group> (<year>2011</year>). <article-title>The highly potent and selective dipeptidyl peptidase IV inhibitors bearing a thienopyrimidine scaffold effectively treat type 2 diabetes</article-title>. <source>Eur. J. Med. Chem.</source> <volume>46</volume>, <fpage>71</fpage>&#x2013;<lpage>76</lpage>. <pub-id pub-id-type="doi">10.1016/j.ejmech.2010.10.016</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Elbachiri</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Fatima</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Bouchbika</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Benchekroun</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Jouhadi</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Tawfiq</surname>
<given-names>N.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>Breast cancer in men: about 40 cases and literature review</article-title>. <source>Pan Afr. Med. J.</source> <volume>28</volume>, <fpage>287</fpage>. <pub-id pub-id-type="doi">10.11604/pamj.2017.28.287.13527</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>El-Sherbeny</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>El-Ashmawy</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>El-Subbagh</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>El-Emam</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Badria</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>1995</year>). <article-title>Synthesis, antimicrobial and antiviral evaluation of certain thienopyrimidine derivatives</article-title>. <source>Eur. J. Med. Chem.</source> <volume>30</volume>, <fpage>445</fpage>&#x2013;<lpage>449</lpage>. <pub-id pub-id-type="doi">10.1016/0223-5234(96)88255-9</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Er-rajy</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>El Fadili</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Hadni</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Mrabti</surname>
<given-names>N. N.</given-names>
</name>
<name>
<surname>Zarougui</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Elhallaoui</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>2D-QSAR modeling, drug-likeness studies, ADMET prediction, and molecular docking for anti-lung cancer activity of 3-substituted-5-(phenylamino) indolone derivatives</article-title>. <source>Struct. Chem.</source> <volume>33</volume>, <fpage>973</fpage>&#x2013;<lpage>986</lpage>. <pub-id pub-id-type="doi">10.1007/s11224-022-01913-3</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Er-rajy</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>El fadili</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Imtara</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Saeed</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Ur Rehman</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Zarougui</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2023a</year>). <article-title>3D-QSAR studies, molecular docking, molecular dynamic simulation, and ADMET proprieties of novel pteridinone derivatives as PLK1 inhibitors for the treatment of prostate cancer</article-title>. <source>Life</source> <volume>13</volume>, <fpage>127</fpage>. <pub-id pub-id-type="doi">10.3390/life13010127</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Er-rajy</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>El fadili</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Mujwar</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Imtara</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Al kamaly</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Zuhair Alshawwa</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2023b</year>). <article-title>Design of novel anti-cancer agents targeting COX-2 inhibitors based on computational studies</article-title>. <source>Arabian J. Chem.</source> <volume>16</volume>, <fpage>105193</fpage>. <pub-id pub-id-type="doi">10.1016/j.arabjc.2023.105193</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Er-rajy</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>El fadili</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Mujwar</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zarougui</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Elhallaoui</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2023c</year>). <article-title>Design of novel anti-cancer drugs targeting TRKs inhibitors based 3D QSAR, molecular docking and molecular dynamics simulation</article-title>. <source>J. Biomol. Struct. Dyn.</source> <volume>41</volume>, <fpage>11657</fpage>&#x2013;<lpage>11670</lpage>. <pub-id pub-id-type="doi">10.1080/07391102.2023.2170471</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Falzone</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Salomone</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Libra</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Evolution of cancer pharmacological treatments at the turn of the third millennium</article-title>. <source>Front. Pharmacol.</source> <volume>9</volume>, <fpage>1300</fpage>. <pub-id pub-id-type="doi">10.3389/fphar.2018.01300</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Faris</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Hadni</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Ibrahim</surname>
<given-names>I. M.</given-names>
</name>
<name>
<surname>Elhallaoui</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>
<italic>In silico</italic> discovery of potent and selective Janus kinase 3 (JAK3) inhibitors through 3D-QSAR, covalent docking, ADMET analysis, molecular dynamics simulations, and binding free energy of pyrazolopyrimidine derivatives</article-title>. <source>J. Biomol. Struct. Dyn.</source> <volume>0</volume>, <fpage>1</fpage>&#x2013;<lpage>17</lpage>. <pub-id pub-id-type="doi">10.1080/07391102.2023.2222839</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Golbraikh</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Tropsha</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Predictive QSAR modeling based on diversity sampling of experimental datasets for the training and test set selection</article-title>. <source>Mol. Divers</source> <volume>5</volume>, <fpage>231</fpage>&#x2013;<lpage>243</lpage>. <pub-id pub-id-type="doi">10.1023/A:1021372108686</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Golbraikh</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Tropsha</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Beware of q2</article-title>. <source>J. Mol. Graph. Model.</source> <volume>20</volume>, <fpage>269</fpage>&#x2013;<lpage>276</lpage>. <pub-id pub-id-type="doi">10.1016/S1093-3263(01)00123-1</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Goodhue</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Lewis</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Thompson</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2006</year>). &#x201c;<article-title>PLS, small sample size, and statistical power in MIS research</article-title>,&#x201d; in <conf-name>Proceedings of the 39th Annual Hawaii International Conference on System Sciences (HICSS&#x2019;06)</conf-name>, <conf-loc>Kauai, HI, USA</conf-loc>, <conf-date>04-07 January 2006</conf-date> (<publisher-name>IEEE</publisher-name>), <fpage>202b</fpage>.</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hopcroft</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wigmore</surname>
<given-names>E. M.</given-names>
</name>
<name>
<surname>Williamson</surname>
<given-names>S. C.</given-names>
</name>
<name>
<surname>Ros</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Eberlein</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Moss</surname>
<given-names>J. I.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Combining the AKT inhibitor capivasertib and SERD fulvestrant is effective in palbociclib-resistant ER&#x2b; breast cancer preclinical models</article-title>. <source>NPJ Breast Cancer</source> <volume>9</volume>, <fpage>64</fpage>. <pub-id pub-id-type="doi">10.1038/s41523-023-00571-w</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jo</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Iyer</surname>
<given-names>V. G.</given-names>
</name>
<name>
<surname>Im</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>CHARMM-GUI: a web-based graphical user interface for CHARMM</article-title>. <source>J. Comput. Chem.</source> <volume>29</volume>, <fpage>1859</fpage>&#x2013;<lpage>1865</lpage>. <pub-id pub-id-type="doi">10.1002/jcc.20945</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kotaiah</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Harikrishna</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Nagaraju</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Venkata Rao</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Synthesis and antioxidant activity of 1,3,4-oxadiazole tagged thieno[2,3-d]pyrimidine derivatives</article-title>. <source>Eur. J. Med. Chem.</source> <volume>58</volume>, <fpage>340</fpage>&#x2013;<lpage>345</lpage>. <pub-id pub-id-type="doi">10.1016/j.ejmech.2012.10.007</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kumari</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Dalal</surname>
<given-names>V.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Identification of potential inhibitors for LLM of <italic>Staphylococcus aureus</italic>: structure-based pharmacophore modeling, molecular dynamics, and binding free energy studies</article-title>. <source>J. Biomol. Struct. Dyn.</source> <volume>40</volume>, <fpage>9833</fpage>&#x2013;<lpage>9847</lpage>. <pub-id pub-id-type="doi">10.1080/07391102.2021.1936179</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kumari</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Kumar</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Lynn</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>g_mmpbsa&#x2014;a GROMACS tool for high-throughput MM-PBSA calculations</article-title>. <source>J. Chem. Inf. Model.</source> <volume>54</volume>, <fpage>1951</fpage>&#x2013;<lpage>1962</lpage>. <pub-id pub-id-type="doi">10.1021/ci500020m</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kumari</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Rathi</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Pathak</surname>
<given-names>S. R.</given-names>
</name>
<name>
<surname>Dalal</surname>
<given-names>V.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Structural-based virtual screening and identification of novel potent antimicrobial compounds against YsxC of <italic>Staphylococcus aureus</italic>
</article-title>. <source>J. Mol. Struct.</source> <volume>1255</volume>, <fpage>132476</fpage>. <pub-id pub-id-type="doi">10.1016/j.molstruc.2022.132476</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Maryam</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Faegheh</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Majid</surname>
<given-names>A.-S.</given-names>
</name>
<name>
<surname>Kazem</surname>
<given-names>N.-K.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Effect of quercetin on secretion and gene expression of leptin in breast cancer</article-title>. <source>J. Traditional Chin. Med.</source> <volume>37</volume>, <fpage>321</fpage>&#x2013;<lpage>325</lpage>. <pub-id pub-id-type="doi">10.1016/S0254-6272(17)30067-5</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="web">
<collab>My Biosoftware</collab> (<year>2019</year>). <article-title>AutoDock 4.2.6/AutoDockTools 1.5.6 &#x2013; suite of automated docking tools &#x2013; my biosoftware &#x2013; bioinformatics softwares blog</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://mybiosoftware.com/autodock-4-2-3-autodocktools-1-5-6-suite-automated-docking-tools.html">https://mybiosoftware.com/autodock-4-2-3-autodocktools-1-5-6-suite-automated-docking-tools.html</ext-link> (Accessed December 29, 2023)</comment>.</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ouassaf</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Belaidi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Khamouli</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Belaidi</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Chtita</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Combined 3D-QSAR and molecular docking analysis of thienopyrimidine derivatives as <italic>Staphylococcus aureus</italic> inhibitors</article-title>. <source>Acta Chim. Slov.</source> <volume>68</volume>, <fpage>289</fpage>&#x2013;<lpage>303</lpage>. <pub-id pub-id-type="doi">10.17344/acsi.2020.5985</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pires</surname>
<given-names>D. E. V.</given-names>
</name>
<name>
<surname>Blundell</surname>
<given-names>T. L.</given-names>
</name>
<name>
<surname>Ascher</surname>
<given-names>D. B.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>pkCSM: predicting small-molecule pharmacokinetic and toxicity properties using graph-based signatures</article-title>. <source>J. Med. Chem.</source> <volume>58</volume>, <fpage>4066</fpage>&#x2013;<lpage>4072</lpage>. <pub-id pub-id-type="doi">10.1021/acs.jmedchem.5b00104</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rashad</surname>
<given-names>A. E.</given-names>
</name>
<name>
<surname>Shamroukh</surname>
<given-names>A. H.</given-names>
</name>
<name>
<surname>Abdel-Megeid</surname>
<given-names>R. E.</given-names>
</name>
<name>
<surname>Mostafa</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>El-Shesheny</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Kandeil</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2010</year>). <article-title>Synthesis and screening of some novel fused thiophene and thienopyrimidine derivatives for anti-avian influenza virus (H5N1) activity</article-title>. <source>Eur. J. Med. Chem.</source> <volume>45</volume>, <fpage>5251</fpage>&#x2013;<lpage>5257</lpage>. <pub-id pub-id-type="doi">10.1016/j.ejmech.2010.08.044</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sanachai</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Mahalapbutr</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Hengphasatporn</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Shigeta</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Seetaha</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Tabtimmai</surname>
<given-names>L.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Pharmacophore-based virtual screening and experimental validation of pyrazolone-derived inhibitors toward janus kinases</article-title>. <source>ACS Omega</source> <volume>7</volume>, <fpage>33548</fpage>&#x2013;<lpage>33559</lpage>. <pub-id pub-id-type="doi">10.1021/acsomega.2c04535</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tong</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>T.-H.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>6-amide-2-aryl benzoxazole/benzimidazole derivatives as VEFGR-2 inhibitors in two-and three-dimensional QSAR studies: topomer CoMFA and HQSAR</article-title>. <source>Chem. Pap.</source> <volume>75</volume>, <fpage>3551</fpage>&#x2013;<lpage>3562</lpage>. <pub-id pub-id-type="doi">10.1007/s11696-021-01588-w</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vald&#xe9;s-Tresanco</surname>
<given-names>M. S.</given-names>
</name>
<name>
<surname>Vald&#xe9;s-Tresanco</surname>
<given-names>M. E.</given-names>
</name>
<name>
<surname>Valiente</surname>
<given-names>P. A.</given-names>
</name>
<name>
<surname>Moreno</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>gmx_MMPBSA: a new tool to perform end-state free energy calculations with GROMACS</article-title>. <source>J. Chem. Theory Comput.</source> <volume>17</volume>, <fpage>6281</fpage>&#x2013;<lpage>6291</lpage>. <pub-id pub-id-type="doi">10.1021/acs.jctc.1c00645</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Van Der Spoel</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Lindahl</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Hess</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Groenhof</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Mark</surname>
<given-names>A. E.</given-names>
</name>
<name>
<surname>Berendsen</surname>
<given-names>H. J. C.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>GROMACS: fast, flexible, and free</article-title>. <source>J. Comput. Chem.</source> <volume>26</volume>, <fpage>1701</fpage>&#x2013;<lpage>1718</lpage>. <pub-id pub-id-type="doi">10.1002/jcc.20291</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J. Z. H.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>End-point binding free energy calculation with MM/PBSA and MM/GBSA: strategies and applications in drug design</article-title>. <source>Chem. Rev.</source> <volume>119</volume>, <fpage>9478</fpage>&#x2013;<lpage>9508</lpage>. <pub-id pub-id-type="doi">10.1021/acs.chemrev.9b00055</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zaorsky</surname>
<given-names>N. G.</given-names>
</name>
<name>
<surname>Churilla</surname>
<given-names>T. M.</given-names>
</name>
<name>
<surname>Egleston</surname>
<given-names>B. L.</given-names>
</name>
<name>
<surname>Fisher</surname>
<given-names>S. G.</given-names>
</name>
<name>
<surname>Ridge</surname>
<given-names>J. A.</given-names>
</name>
<name>
<surname>Horwitz</surname>
<given-names>E. M.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>Causes of death among cancer patients</article-title>. <source>Ann. Oncol.</source> <volume>28</volume>, <fpage>400</fpage>&#x2013;<lpage>407</lpage>. <pub-id pub-id-type="doi">10.1093/annonc/mdw604</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zarougui</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Er-rajy</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Faris</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Imtara</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>El fadili</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Al kamaly</surname>
<given-names>O.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>QSAR, DFT studies, docking molecular and simulation dynamic molecular of 2-styrylquinoline derivatives through their anticancer activity</article-title>. <source>J. Saudi Chem. Soc.</source> <volume>27</volume>, <fpage>101728</fpage>. <pub-id pub-id-type="doi">10.1016/j.jscs.2023.101728</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>L.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Recent advances in iodine-promoted C&#x2212;S/N&#x2212;S bonds formation</article-title>. <source>Chem. &#x2013; A Eur. J.</source> <volume>26</volume>, <fpage>17289</fpage>&#x2013;<lpage>17317</lpage>. <pub-id pub-id-type="doi">10.1002/chem.202001414</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Perez-Sanchez</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Lightstone</surname>
<given-names>F. C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>A comprehensive docking and MM/GBSA rescoring study of ligand recognition upon binding antithrombin</article-title>. <source>Curr. Top. Med. Chem.</source> <volume>17</volume>, <fpage>1631</fpage>&#x2013;<lpage>1639</lpage>. <pub-id pub-id-type="doi">10.2174/1568026616666161117112604</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ziada</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Diharce</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Raimbaud</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Aci-S&#xe8;che</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Ducrot</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Bonnet</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Estimation of drug-target residence time by targeted molecular dynamics simulations</article-title>. <source>J. Chem. Inf. Model.</source> <volume>62</volume>, <fpage>5536</fpage>&#x2013;<lpage>5549</lpage>. <pub-id pub-id-type="doi">10.1021/acs.jcim.2c00852</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zouaghi</surname>
<given-names>M. O.</given-names>
</name>
<name>
<surname>Amri</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Hassen</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Arfaoui</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>&#xd6;zdemir</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>&#xd6;zdemir</surname>
<given-names>I.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Biological determination, molecular docking and Hirshfeld surface analysis of rhoduim(I)-<italic>N</italic>-heterocyclic carbene complex: synthesis, crystal structure, DFT calculations, optical and non linear optical properties</article-title>. <source>Inorganica Chim. Acta</source> <volume>551</volume>, <fpage>121459</fpage>. <pub-id pub-id-type="doi">10.1016/j.ica.2023.121459</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>