Abstract

According to the World Health Organization, breast cancer represents 16% of all cancer cases in women and is the second most common cancer. In the past decades, the mortality among patients with metastasis breast cancer has been reduced significantly via drug delivery by means of nanodiamond therapies, which are both biocompatible and scalable. In this study, we determined a theoretical pathway for the construction of a nanodiamond–tamoxifen complex that will act as a drug delivery vehicle for targeting tumor tissues of breast cancer. The tamoxifen pharmacophore was defined and the binding zone was identified for the electrostatic interaction between tamoxifen and a functionalized site of a nanodiamond particle allowing for attachment of the payload (this drug) to the surface of the nanodiamond particle. In addition, an analysis of the intermolecular interaction between the nanodiamond and tamoxifen was conducted, showing three hydrogen bonds complying fully with Lipinski’s rule of five, which states that a compound should have five or fewer hydrogen bonds to be permeating and easily absorbed by the body (qualitative prediction). All calculations were performed using the conceptual Density Functional Theory with the M06 functional and the basis set 6-31G(d). The solvent effect has been taken into account by an implicit model, the conductor like polarizable continuum model.

1. Introduction

Breast cancer is an important and increasingly common neoplastic disease in women and is actually treated with different drugs, with tamoxifen (TAM) being one of them. It is an antiestrogen agent [1, 2] belonging to the family of selective estrogen receptor modulators (SERMs); TAM is most often used for the treatment of estrogen receptor-positive (ER) breast cancer in pre- and postmenopausal women and has helped to reduce breast cancer death rates [3]. TAM is also used for prevention of breast cancer in women at high risk of this disease [4]. Although aromatase inhibitors are currently available for breast cancer treatment in postmenopausal women [5], TAM is still the drug most frequently used for breast cancer treatment [4]. Nonetheless, prolonged use of TAM causes adverse effects such as endometrial cancer, thromboembolism, and menopause symptoms [6, 7].

Conventional cancer therapies face great challenges such as poor bioavailability and intrinsic toxicity. Therapeutic efficacy of many useful drugs is compromised by their toxicity. Nanomaterials with pharmacological and therapeutic alterations have overcome some of these conventional limitations [8]. One of the nanomaterials used to deliver chemotherapeutic drugs is nanodiamond (ND), a carbon-based nanomaterial that has caught the attention of investigators for its biological applications, physical and chemical properties, stability, scalability, small size, good adsorption, biocompatibility [911], and easy functionalization [12], allowing researchers to generate groups (ether –COC, peroxide –COO–, carbonyl –C=O, or carboxyl –COOH) on the ND surface [13, 14]. This approach gives investigators the opportunity to use these derivatives for specific or nonspecific binding of nucleic acids and proteins [13, 14]: a possible application to biomedicine. In addition, this surface can be modified with biologically active molecules by adsorption or covalent or noncovalent chemical immobilization [15, 16]. Furthermore, the functional groups on the surface of an ND particle can interact electrostatically or chemically with a bimolecular appurtenant in a sample under study for cell-specific interaction and targeting [16]. Within functionalization platforms of ND, some investigators have also reported the use of polylactides in general and poly(L-lactide) (PLLA) in particular, which may be viewed as a second-generation biomaterial approved by the US Food and Drug Administration (FDA). PLLA has been widely used in biodegradable polymer microspheres because of its great advantages of biodegradability and biocompatibility [17, 18], due to the controlled polymer chains that can be broken down to lactic acid monomers (LA). The latter can be metabolized by the body [19] to CO2 and water in the Krebs cycle [20]; the abovementioned data are suggestive of good safety, that is, limited toxicity or tissue reactions [17].

According to studies by Ho’s group at the University of California at Los Angeles (UCLA) [2124], NDs have several unique properties that make them a promising nanomaterial for biomedical applications; these properties include unique electrostatic properties, a chemically inert core, and a tunable surface [21]. Some studies indicated that drug molecules can be adsorbed by ND thus forming the complex drug–ND, and with these approaches, water-insoluble drugs may be delivered to cells to effectively exert biological action. The main drugs involved in these studies on drug delivery vehicles, including those for cancer treatment, are doxorubicin hydrochloride (DOX) [22, 25], purvalanol A, and 4-hydroxytamoxifen [24] as well as insulin therapeutics for diabetes [23].

It has been suggested that the molecules of chemotherapeutic drugs bond with the ND surface mainly through electrostatic interactions or hydrogen bonds (HBs) between the drug to be studied and a carboxyl or hydroxyl group of functionalized ND [2225]. Theoretical studies on chemical properties of the ND–TAM complex are important for understanding the relation between this complex and interaction with an estrogen receptor. Recent Density Functional Theory (DFT) studies have shown that the electronic orbitals in nanoscale diamonds are not homogeneously distributed within individual particles or grains [26]. In addition, theoretical studies based on ab initio calculations within the DFT framework were performed to determine optical and electronic properties of ND [27]. The loading of daunorubicin (DNR) on ND was optimized by adjusting reaction parameters such as acidity and concentration through molecular dynamics simulations [28].

The main purpose of this study was to load a drug on the surface of ND by molecular modeling techniques, without altering the drug’s pharmacodynamics. This loading can be achieved through electrostatic interactions (HBs) between the surface of ND and TAM to deliver the drug specifically to estrogen receptors.

2. Materials and Methods

2.1. Computational Details

Characterization of the TAM molecule was performed to determine the core unit (of the drug) that is responsible for the biological activity. In addition, analysis of the esterification of lactic acid with the functionalized ND enabled us to generate a five-atom linear chain. Furthermore, characterization of the ND–TAM complex was performed.

The theoretical study was carried out using DFT [29, 30] in the Gaussian 09 software [31]. The hybrid meta-GGA functional considered in this study is M06 [32, 33] (developed by the Truhlar group at the University of Minnesota) combined with basis set 6-31G(d) (proposed by People [34]) and the conductor like polarizable continuum model [35] using water as a solvent. This methodology was previously validated by us. Optimization of the ND–TAM complex was performed by semiempirical method PM6 [36] (to achieve convergence of the multiatom complex) and by taking into account the good results on structural parameters of the method [3739] with water as a solvent in the conductor like polarizable continuum model. The rest of the electronic properties were calculated by means of M06/6-31G(d). Density functional methodology has been widely used to elucidate the chemical reactivity defined as a set of chemical concepts such as the ionization potential and electron affinity (EA), which are obtained by the energy difference and , respectively, where indicates the parent molecule and N − 1 and N + 1 correspond to the cation and anion radicals generated after electron transfer. Once the ionization potential and EA are obtained, other reactivity parameters are calculated: chemical hardness () [40], electrophilicity () [41], chemical potential () [42], electronegativity () [43, 44], and chemical softness [45].

DFT has been used to understand the selectivity of molecular systems [4648]. The reactivity of a particular site of a molecular species can be explained by local quantities such as electron density, Fukui function (FF), and local softness (), which have been used successfully in studies on site selectivity in a molecule [49].

FF are mathematical expressions that define the sensitivity that a molecular system has to have to undergo changes in its electronic density, at different sites in its structure [50]. FF is defined asFor practical purposes, quantification of FF is possible through a condensation scheme in an atomic region of the molecule. This reaction can proceed under the influence of an electrophilic, nucleophilic, or radical attack at a particular reaction site. The use of procedures of population analysis results in the following equations [51]:In this study, we used exclusively the equation that describes the electrophilic attack (3) because of the nature of the interaction within the ND–TAM complex. Similar to condensed FF for an electrophilic attack and global softness to an electrophilic attack, local softness () is given bywhere S is the global softness and is FF as defined in (3).

3. Results and Discussion

Modeling of the ND–TAM complex requires several steps. First, the geometry optimization of the molecules involved was defined, and the TAM density distribution was evaluated to determine the binding site within the nanostructure. The esterification of the nanocarrier with lactic acid was selected. Also, the calculation of reactivity parameters for both systems was performed. Then, the ND–TAM equilibrium distance calculation was carried out. These phases are explained below.

3.1. Geometry Optimization

The geometry was optimized for TAM, functionalized ND, and the complex of functionalized ND with lactic acid in the aqueous phase; the frequency calculations were performed to guarantee that the structures are at their lowest energy level. Optimized geometries were obtained by the M06/6-31G(d) computational methodology. The structures are shown in Figure 1.

3.2. Frontier Orbital Analysis in TAM

Local softness as a reactivity descriptor was used successfully for predicting the selectivity site in TAM. The local softness for electrophilic attack () of TAM is presented in Figure 2. Because the 4C atom belongs to the alkene functional group and has the highest value (0.0163), this site is the most susceptible to an electrophilic attack.

According to Thomas [52], knowledge of the electron density makes it possible to determine the pharmacophore of a drug, explaining how structurally diverse ligands can bind to a common receptor site. Thus, using the minimal energy TAM structure, an analysis of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) was conducted to identify the pharmacophore of TAM. The results of mapping of the frontier molecular orbitals showed that the electron-saturated zone of TAM is in the phenyl, ethyl, and alkene functional groups, and according to local softness, the 4C atom (which belongs to the alkene functional group that is inside the electron-saturated zone) is the possible binding site of the drug. This means that the electrostatic binding site for ND’s functionalization is the chain localized outside the electron density area (see Figure 3).

3.3. Nanodiamond Esterification

The stability of various phases of nanoscale carbon has been the focus of numerous theoretical and computational studies. It was discovered that the morphology plays an important role in the stability of NDs by influencing the surface reconstruction and formation of sp2 carbon [53]. Mochalin et al. agree that the NDs that have a range of 1.5–2.0 nm (or <1 nm) are stable [5355]. In the present study, the ND was analyzed via DFT and was found to have an octahedral surface showing a transition from carbon to carbon (when it is not functionalized) at 1.5 nm diameter with (111) facets and with 351 atoms (Figure 1).

The functionalization of the nanostructure begins with acid treatment that creates carboxyl functional groups on its surface [56, 57]. In addition, esterification is done using a lactic acidic monomer of PLLA, which has proven to be biocompatible, and the carboxyl functional group was found on the surface in the carrier vehicle [58]. Next, those authors generated a linear chain of five atoms to ensure better interaction between the ND and TAM. The esterification enabled stronger electrostatic interaction (HBs) between the linear chain of five atoms generated in the carrier vehicle with the linear chain of TAM defined by the electron density analysis composed of ether and tertiary amine located in the drug, thereby forming the ND–TAM complex as shown in Figure 4. This proposed drug delivery allows for better biodistribution and controlled release of the drug at the specific targeted site. This feature makes this vehicle an ideal carrier for drug delivery and targeting applications.

3.4. Reactivity Parameters

Reactivity parameters that we analyzed in TAM and ND are ionization potential, which is defined as the energy needed to remove an electron from a molecule; EA, which measures the ability of a molecule to accept electrons or form anions; electronegativity (), which represents the tendency of atoms or molecules to attract electrons; electrophilicity (), which predicts stabilization in energy when the system acquires electrons from the environment up to saturation; chemical hardness (), that is, resistance to changes in the electronic configuration; and chemical potential (), which measures the tendency of electrons to be released from a system, in such a way that a large chemical hardness difference between two systems favors the electron transference; also, the difference between HOMO and LUMO energy levels (), which reflects biological activity of the molecule in question [59]. The results are presented in Table 1.

Chemical hardness approximates reactivity according to the principle of maximum hardness [60, 61], which tells us that the most reactive systems have low hardness, whereas less reactive systems have high hardness. According to the results of this study, the system with lower chemical hardness is ND with 1.35 eV; therefore, this is the most reactive system according to the chemical hardness concept mentioned above.

Another important concept that can characterize a bond is the electronegativity difference. When it is greater than 1.7 eV, electronic attraction is so great that electrons are transferred from one atom to another; this situation is referred to as an electrostatic interaction [62] (ionic bonds, HBs). Results in Table 1 show a difference in the electronegativity of 1.91 eV between TAM and ND; this indicate an electrostatic attraction.

Electron density of the HOMO and LUMO plays an important role in several chemical and pharmacological processes. According to Galindo-Murillo et al. [63], the greater , the greater the ability to give electrons. The smaller , the lower the resistance to accept electrons. TAM is the compound with the highest electronegativity and lowest value in this case; therefore, TAM will bring electrons to ND. Also, Helal et al. [64] in their molecular modeling study (on a series of pyridine derivatives used against cancer) concluded that the molecules with low do not donate their electrons easily; on the other hand, a higher energy implies that the molecule is a good electron donor. In their work, they found that the lower the gap, the easier the charge transfer interaction taking place within the molecule. The statements above indicate that the biological activity is related to the energy difference between the molecular orbitals HOMO and LUMO of the molecule; those authors concluded that their pyridine derivatives with an energy gap between 7.54 and 7.69 eV have a good ability to interact with cancerous tissues. In the present study, TAM showed an HOMO–LUMO gap energy of 4.98 eV, which implies better capacity for interaction with cancerous tissues.

3.5. Calculation of the ND–TAM Equilibrium Distance

A rigid potential-energy surface scan was performed on the optimized geometry of ND and TAM by the method mentioned above in Computational Details. A rigid potential-energy surface that consists of single-point energy assessments across a rectangular grid involving a scan of selected internal coordinates was realized in 15 steps with approximations of 1 Å. The energy curve for ND and TAM shown in Figure 5 indicates the minimal energy state between these structures, with −8204.01479 Hartrees defining the equilibrium distance between TAM and ND of 4 Å. Nonetheless, given that the potential-energy surface scan calculations do not include an optimized geometry [65], the ND–TAM complex was optimized by the semiempirical PM6 method using water as a solvent, followed by calculation of frequencies to verify that the ND–TAM complex is located at the lowest energy level. After optimization, the distance between TAM and ND remained 3.8 Å.

3.6. Hydrogen Bonds (HBs)

Once the complex is completely constructed, its electronic configuration is analyzed, particularly the HBs. These HBs are responsible for many properties of biological molecules. They provide the flexibility and specificity required in biological processes and are responsible for the regulation of rates and directions of biological reactions allowing for implementation of a specific biological objective [66]. Furthermore, HBs are used to build supermolecular capsules or reversibly associated polymers because they enable self-assembly and molecular recognition [67].

A comparison of the theoretical IR spectra of a free A–H group (no-HB group) and an A–H⋯B HB group was conducted. The A–H group corresponds to the TAM drug whose stoichiometry is C26H29NO, and the A–H⋯B HB group corresponds to the ND–TAM complex with stoichiometry C201H213NO5. According to Desiraju, when comparing the IR spectrum of a free A–H group and another one with an HB, in the second one, the A–H stretching frequency is reduced and indicates hydrogen bonding [68].

In the analysis of the IR spectra, it was found that one of the C–H stretching frequencies in TAM and in the ND–TAM complex is 2915 and 2721 cm−1, respectively. The other C–H stretching is at 2935 cm−1 in the TAM and 2769 cm−1 in the ND–TAM complex. The R–O–R vibration is present at 1098 cm−1 in TAM and at 1065 cm−1 in the ND–TAM complex. The reduction of all the frequencies in the ND–TAM complex indicates the presence of HBs. The IR spectra are presented in Figure 6.

Once the presence of HBs was corroborated, a method for analysis of the different types of HBs was developed. The theoretical bond lengths and angles of HBs C=O⋯H–C, H–O⋯H, and C–H⋯O (all of them present in the ND–TAM complex) suggested that these HBs can be considered weak according to Jeffrey’s classification [69]. Results are shown in Table 2.

The HBs formed in the ND–TAM complex include the C–H⋯O interaction, where groups C–H and O–H act as donors. This is possible, in agreement with the study by Stainer [70], who states that the atom connected with the hydrogen to form an HB does not need high electronegativity; it simply has to be slightly polar. This requirement includes such groups as C–H, aryl–H, and some metal hydrides (Figure 7).

The number of HBs is also relatively important because, according to Lipinski’s rule of five, this number can determine qualitatively whether a compound is more permeating and easily absorbed by the body [71]; thus, a compound should have five or fewer HBs. The ND–TAM complex showed three HBs here, complying fully with this rule.

4. Conclusion

The TAM pharmacophore was determined and a theoretical functionalization method was applied to the ND–TAM complex. Stabilization of the ND–TAM complex was achieved through electrostatic interactions between the surface of the carrier vehicle and the drug.

Functionalization was implemented by means of TAM and the ND esterified beforehand with lactic acid. This procedure creates a linear five-atom chain that coincides with the linear five-atom chain in the drug; this situation generates a sufficient electrostatic force for ND to act as a carrier of TAM, and the complex will be absorbed by the cell for eventual release and pharmacological action on cancer cells.

The complex formed due to electrostatic interactions between the ND and TAM generates three HBs: interactions between C=O and O–H groups and the oxygen of ND; interaction between the C–H group and oxygen of TAM. This situation ensures stability in the complex. Furthermore, analysis of the calculated IR spectrum between ND and the ND–TAM complex shows the presence of HBs that are weak according to Jeffrey’s classification.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by Consejo Nacional de Ciencia y Tecnología (CONACYT) and Centro de Investigación de Materiales Avanzados, SC (CIMAV). Linda-Lucila Landeros-Martínez gratefully acknowledge a fellowship from CONACYT. Erasmo Orrantia-Borunda is a researcher at CIMAV and CONACYT, and Norma Flores-Holguin is a researcher at CIMAV and CONACYT. David Chavez-Flores is a researcher at UACH and CONACYT.