Structural, Vibrational, and Electronic Properties of the Glucoalkaloid Strictosidine: A Combined Experimental and Theoretical Study

Costa, Renyer Alves; Pinheiro, Maria Lucia Belem; Oliveira, Kelson Mota Teixeira de; Barison, Andersson; Salomé, Kahlil Schwanka; Iank, Júlio Rodolfo; da Silva, Noam Gadelha; Cabral, Tiara Souza; Costa, Emmanoel Vilaça

doi:https://doi.org/10.1155/2016/1752429

Journal of Chemistry

On this page

Abstract Introduction Results and Discussion Conclusion Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2016 | Article ID 1752429 | https://doi.org/10.1155/2016/1752429

Structural, Vibrational, and Electronic Properties of the Glucoalkaloid Strictosidine: A Combined Experimental and Theoretical Study

Renyer Alves Costa,¹Maria Lucia Belem Pinheiro,¹Kelson Mota Teixeira de Oliveira,¹Andersson Barison,²Kahlil Schwanka Salomé,²Júlio Rodolfo Iank,¹Noam Gadelha da Silva,¹Tiara Souza Cabral,³and Emmanoel Vilaça Costa¹

Academic Editor: Arturo Espinosa Ferao

Received21 Oct 2015

Accepted03 Jan 2016

Published08 Mar 2016

Abstract

A detailed structural analysis and spectral behavior of the glucoalkaloid strictosidine, a precursor of all monoterpene indole alkaloids, are discussed. The experimental NMR, FTIR, and UV results were compared to the theoretical DFT spectra calculated by Becke using the three-parameter Lee-Yang-Parr (B3LYP) function with 6-31G(d) and 6-311++G(2d,p) basis sets. The theoretical geometry optimization data were compared with the X-ray data for precursors and similar structures in the associated literature. The similarity between the theoretical and experimental coupling constants values made it possible to affirm the values of dihedral angles and their configuration, reinforcing findings from previous stereochemical studies. Theoretical UV analysis agreed well with the measured experimental data, with bands assigned. Calculated HOMO/LUMO gaps show low excitation energy for strictosidine, justifying its stability and reaction kinetics. The molecular electrostatic potential map shows opposite potentials regions that form hydrogen bonds that stabilize the dimeric form, which were confirmed by excellent agreement of the dimeric form theoretical wavenumbers with the experimental IR spectrum. ESI-MS/MS data revealed patterns for the fragmentation of the protonated strictosidine molecule outlined by an NBO study.

1. Introduction

Indole alkaloids play a very important role in the chemistry of natural products and are especially recognized for their use in clinical medicine as an adjunct to anesthetics. The finding of several clinic uses has driven intense study of this class of substances, and many antiplasmodial [1, 2], cytotoxic [3], antibacterial [4], antifungal [5], spasmodic [6], hypotensive [7], and anti-inflammatory [8] properties have been related to indole alkaloids. Strictosidine (Figure 1) is a key glucomonoterpene indole alkaloid precursor of all indole monoterpene alkaloids. This crucial molecule, which originates from a reaction between tryptamine and the monoterpene glycoside secologanin, is found in several plant species [9–13] and was first isolated from Rhazya stricta [14]. Strictosidine has also been obtained in cell suspension cultures and under biomimetic conditions [15, 16].

(a)

(b)

Studies discussing the structure and stereochemistry of strictosidine [16–18] have compared its spectral data with those of similar structures, confirming that its C3 atom has the S configuration or 3α [S]. This configuration is identical to that of known monoterpene indole alkaloids, in disagreement with the proposal that vincoside was the precursor of indole alkaloids with 3 [R] configuration, as previously thought. Patthy-Lukáts et al. [19] studied the stereochemistry of strictosidine based on experimental NMR analysis, determining its spatial configuration. However, there are no X-ray studies because this molecule has not been obtained yet in crystalline form [19] and no theoretical molecular modeling study that discusses the bond lengths and planar and dihedral angles was previously presented. The determination of the relationship between theoretical vibrational frequencies and experimental IR absorbance bands and between theoretical electronic transitions and experimental UV bands of strictosidine has not been investigated yet. Therefore, a detailed theoretical DFT and experimental investigation of the structure and spectral behavior of this molecule, providing a comprehensive description of strictosidine, have been reported. Initially, the alkaloid has been isolated from Strychnos amazonica and the molecule was characterized by NMR (¹H, HSQC, and HMBC), MS (ESI-MS/MS), UV, and FTIR. The theoretical data (optimized geometry, UV, IV, MEP, and NBO calculations) were compared with the experimental data to answer questions regarding structure, electronic transitions involved in the UV spectrum, vibrational assignments, and other physical properties of strictosidine.

2. Experimental Section

2.1. General Procedures

Nuclear magnetic resonance (NMR) spectra were recorded on a Bruker 600 AVANCE spectrometer with 12.4-Tesla magnetons at 295 K in deuteriomethanol containing 0.06% TMS as an internal standard. The ESI-M data were recorded on LQC fleet ion trap Thermo Scientific spectrometer equipped with an electrospray ionization (ESI) ion source in the positive mode with the following parameters: capillary voltage at 26 V; spray voltage at 5 kV; tube lens offset at 100 V; capillary temperature at 225°C; auxiliary gas at 5%; and sheet gas at 12%. The UV data were recorded in methanol using a PDA Detector plus Finnigan Surveyor (Thermo Scientific). The FTIR data were recorded in KBr pellet technique (solid phase) using an ABB FTLA200.

2.2. Plant Material

The leaves of Strychnos amazonica were collected in the Adolpho Ducke Forest Reserve located 25 km from the city of Manaus (2°56′01.0′′S, 59°57′45.8′′W), Amazonas, Brazil. The species was identified by the DNA barcoding technique, in which the distance matrices, which were calculated for rcbL and rpoC1 and for the concatenated gene sequences (rbcL + rpoC1), showed null values of genetic divergence between the collected specimen and a voucher specimen of Strychnos amazonica (INPA 216208) deposited in the herbarium of the National Research Institute of Amazonas (INPA). The rcbL and rpoC1 sequences were deposited in GenBank under accession numbers KJ764797 to KJ764819.

2.3. Isolation of Strictosidine

30 g of leaf extract in methanol was submitted to a silica-gel column using hexane/ethyl acetate and ethyl acetate/methanol as the mobile phases. Of the 31 fractions obtained, fractions 10 to 16 showed the presence of an alkaloid. These fractions were pooled and submitted to an ESI-TI-MS analysis, which showed an ion peak at m/z 531. CID analysis revealed fragments at m/z 514, m/z 369, and m/z 356, which are compatible with strictosidine [20]. This fraction was submitted to a silica-gel column using 100% ethyl acetate, 9 : 1 ethyl acetate/methanol, and 100% methanol in succession, resulting in an amorphous brown solid (20 mg) that appeared to be pure in CCD/Dragendorff and ESI-IT-MS analyses (with only one ion peak at m/z 531). Then, it was subjected to NMR (H¹, HMBC, and HSQC), UV, and IR analysis for structural confirmation (see Table S1 and Figures S1, S2, and S3 in the Supplementary Material data available online at http://dx.doi.org/10.1155/2016/1752429).

2.4. Computational Methods

The theoretical quantum chemical calculations were performed using the Gaussian 03 W Program (Revision E.01) on the Debian Linux (5.0 version) platform on an Intel Quadcore™ PC (8 GB RAM) [21]. The DFT approach was used to optimize the geometry using the 6-31G(d) and 6-311G++(2d,p) basis sets and the B3LYP functional. The molecular geometries were fully optimized by the force gradient method using Berny’s algorithm, and the potential energy surfaces were characterized using standard analytical harmonic vibrational analysis to confirm that the stationary points corresponded to the minima of the potential energy surfaces (no imaginary frequencies or negative eigenvalues were found). The theoretical coupling constants were calculated using the NMR protocols implemented using the DFT B3 LYP 6-31G(d) and 6-311++G(2d,p) basis sets in the Gaussian 03 software. These coupling constants were compared with the measured values that were experimentally obtained, showing RMSD values of 1.45 Hz for DFT B3LYP 6-31G(d) and 1.25 Hz for 6-311++G(2d,p). The UV spectra were calculated using the RTD-B3LYP-FC functional and 6-31G(d) and 6-311G++(2d,p) basis set [22, 23]. The NBO values were obtained with NBO 3.1, as implemented in the Gaussian 03 package using the 6-31G(d) basis set and the B3LYP functional. The harmonic frequencies were calculated at the B3LYP/6-31G(d) level using the optimized structural parameters. The assignments of the calculated wavenumbers were aided by the animation option of the GaussView program, which gives a visual presentation of the vibrational modes [24]. The potential energy distribution (PED) was calculated with the help of the VEDA4 software package [25].

3. Results and Discussion

3.1. Geometry Optimization

Because no crystallographic data for strictosidine are available, the geometry optimization data, which were calculated using the B3LYP/6-31G(d) and B3LYP/6-311++G(2d,p) basis set, were compared with the X-ray data for similar structures [26–29]. Due to the similarity of the data obtained from both basis sets (Table 1), the discussion that follows is based on the B3LYP 6-31G(d) values but should also apply to B3LYP 6-311++G(2d,p).

Initially, the bond lengths indicated a small distortion in the pentacyclic tryptophan ring, showing a greater distance for the C7-C8 (1.44 Å) and C8-C13 (1.42 Å) bonds compared to the N1-C2, C2-C7, and C13-N1 bonds. Significant distortions were observed in the second six-membered ring of the tryptophan portion, which was distinct for all connections. Similar distortions were observed in the dihydropyran ring: 1.52 Å (C15-C16), 1.35 Å (C16-C17), 1.34 Å (C7-O), 1.53 Å (C20-C21), and 1.56 Å (C20-C15). These distortions resulted from the presence of double bonds and heteroatoms, which make these bonds shorter than the C-C bonds.

Concerning the planar angles, the aromatic ring shows no large deformations, which was similar to the experimental data for previously analyzed indole alkaloids [26, 29]. The planar angles were 118.78° (C9-C8-C13), 119.16° (C8-C9-C10), 121.06° (C9-C10-C11), 121.17° (C10-C11-C12), 117.62° (C11-C12-C13), and 122.18° (C8-C12-C13). The pentacyclic ring showed a nonangular uniformity except for the angles between N1-C2-C7 (109.56°) and C13-C2-N1 (109.22°).

¹H-¹H NMR coupling constants are highly sensitive to the dihedral angles of hydrogen atoms and both the configuration and conformation of a structure can be validated if there is agreement between the experimental and theoretical values. The calculated initial dihedral angles for the hydrogen atoms in positions 5 and 6 of the third ring of the tryptophan portion were 162.78° (with theoretical Hz) and 46.76° (theoretical Hz). These values were consistent with the experimentally measured values for coupling constants, Hz and Hz, revealing a pseudoaxial-axial relationship between hydrogen atoms H5a-H6b and an axial-equatorial relationship between H5a and H6a. The obtained dihedral angles for hydrogen atoms of the C3-C14-C15 bridge also showed coupling constants values consistent with the experimental values: 172.28° for H3-H14proR with theoretical (experimental Hz), −72.41° for H3-H14proS with Hz ( Hz), and −87.53° for H14proR-H15 with Hz ( Hz). The conformation of the dihedral angles involving H15-H20 (−56.90°), H20-H21 (−178.69°), and H19-H20 (65.23°) is also plausible given the respective experimental hydrogen coupling constants Hz ( Hz), ( Hz), and Hz ( Hz). The similarity between the modeled DFT structure (Figure 2) (with 15S, 20S, and 21S configuration) and the experimental spectroscopic NMR data complements the conformational arrangement study of strictosidine. The 15S, 20S, and 21S configuration was reported in previous experimental X-ray studies of strictosidine precursors [19, 28–30]. Finally, due to the consistency between the geometry calculated using the DFT B3LYP/6-31G(d) and B3LYP/6-31++G(2d,p) methods and the coupling constant values measured from the experimental ¹H spectrum, the conformation of the C21-O-C1′ bridge (so far uncertain) was depicted with angles of 81.17° and lengths of 2.65 Å between H21 and H1′ (Figure 2). This depiction is consistent with the NOESY data provided by Patthy-Lukáts et al. [19]. Other interactions between segments found in the NOESY data are justified by the calculated distances between the hydrogen atoms as follows: H3-H5b (3.76 Å), H3-H14S (2.55 Å), H5a-H5b (1.76 Å), H6a-H6b (1.76 Å), H9-H10 (2.48 Å), H10-11 (2.477 Å), H11-H12 (2.488 Å), H14R-H14S (1.75 Å), H14R-H19 (2.83 Å), H14R-H15 (2.66 Å), H14S-H21 (2.59 Å), H15-H20 (2.41 Å), H18Z-H18E (1.85 Å), H18Z-H20 (3.57 Å), H18E-H19 (2.40 Å), and H20-H19 (2.52 Å) (for information about Mulliken charges and natural bond analysis, see Tables S2 and S3 in Supplementary Material data).

3.2. UV Analysis

The electronic spectrum of the molecule in a methanol solution was compared with the calculated spectrum (in the gas phase) at time dependent density functional using the B3LYP 6-31G(d) and B3LYP 6-311G++(2d,p) basis sets, as Figure 3 shows. The bands located at 222 and 272 nm could be experimentally observed and were in agreement with the presence of chromophores. Because strictosidine is an aromatic compound, transitions were involved but the presence of conjugations in the pentacyclic ring of the tryptophan portion and in the -alkoxyacrylate group suggested that transitions also occurred. The theoretical calculations predicted an intense electronic transition of 5.8964 eV, with an oscillator strength at 210.27 nm for B3LYP 6-31G(d). For B3LYP 6-311++G(2d,p), an electronic transition of 5.62 eV was predicted, with an oscillator strength at 220.62 nm. These results showed good agreement with the measured experimental data (222 nm) assigned to the and transitions of the indole portion. The calculations also predicted a weaker electronic transition at 265.72 nm (4.66 eV) for B3LYP 6-31G(d) and at 278.26 nm (4.45 eV) for B3LYP 6-311++G(2d,p), which were equivalent to the band at 272 nm in the experimental spectrum. These values were assigned to the sum of the and transitions of the indole and -alkoxyacrylate groups.

The highest occupied molecular orbitals (HOMO) and the lowest-lying unoccupied molecular orbitals (LUMO), also called frontier molecular orbitals (FMOs), play an important role in many properties of a compound as well as in its quantum chemistry and UV-Vis spectra. The energy gap between the HOMO and LUMO energies is the basis for the chemical stability and reactivity of a molecule. The quantum bonding features of strictosidine are depicted by a plot of the HOMO, HOMO−1, LUMO, and LUMO−1 in Figure 4. In the UV-Vis spectrum, the two maximum calculated absorption wavelengths corresponded to the contributions of the electronic transitions from HOMO LUMO+4 (50.7%) and HOMO−1 LUMO+1 (23.74%) for 210 nm and HOMO LUMO+1 (89.8%) for 265 nm in B3LYP 6-31G(d). In B3LYP 6-311G++(2d,p), the wavelengths corresponded to HOMO LUMO+16 (26%), HOMO LUMO+19 (14%), and HOMO LUMO+18 (5.7%) for 220.6 nm and HOMO LUMO+2 (19%) and HOMO LUMO+3 (70.81%) for 278.3 nm. It is clear from Figure 5 and Figures S5 and S6 (see Supplementary Information) that the major transitions were restricted to the indolic portion because the HOMO, HOMO−1, and LUMO+1 are located in this region and only the LUMO is located in the -alkoxyacrylate group. Transitions involving the LUMO contribute minimally to the absorption wavelengths of strictosidine, showing that the transitions in the -alkoxyacrylate group have no influence, in disagreement with the ancient proposal that -alkoxyacrylate group influences strictosidine UV spectrum [14]. This prediction was proven by the similarity of the experimental spectrum of strictosidine to the spectra of several other indole alkaloids that lack alkoxyacrylate groups [31].

3.3. Global and Local Reactive Descriptors

The energy gap between the HOMO and LUMO is very important for determining the electrical properties, kinetic stability, optical polarizability, and chemical reactivity descriptors, such as hardness and softness, of a molecule.

The concept of hardness and softness is related to a compound’s reactivity and is a property that measures the extent of chemical reactivity to which the addition of a charge stabilizes the system. The chemical potential provides a global reactivity index and is related to charge transfer from a system of higher chemical potential to one of lower chemical potential. Electronegativity is the power to attract electrons and is directly related to all the previously mentioned properties. All these properties are defined as follows [32, 33]:where and are the electronic energy and the external potential of an N-electron, respectively. Based on Koopmans theorem for closed-shell molecules, these global chemical reactivity descriptors can be simplified and defined as follows: where is the ionization potential and is the electron affinity of the molecule. The ionization energy and electron affinity are obtained from the HOMO and LUMO energies as and . In terms of chemical hardness, a large HOMO-LUMO gap indicates a hard molecule and is related to more stable molecules, whereas a small gap indicates a soft molecule and is related to a more reactive molecule.

Another important descriptor is the electrophilicity index (), a global maximum reactivity index that is similar to chemical hardness and chemical potential. The electrophilicity index measures the global electrophilic nature of a molecule and was proposed by Parr et al. [34, 35] as a measure of energy lowering due to charge transfer. The electrophilicity index is defined as follows:

This scale permits the classification of organic molecules as strong, eV, moderate, eV, and marginal, eV, electrophiles. On the other hand, there is a good correlation in the inverse of the electrophilicity (); thus, molecules located at the bottom of the electrophilicity scale are classified as marginal electrophiles, corresponding with good nucleophiles [36]. However, when the molecule bears more than one functional group with opposite electrical charge, its nucleophilic character cannot be straightforwardly associated with the inverse of the electrophilicity. Thus, the nucleophilicity index () appears as a different descriptor which gives more information about nucleophilicity and is defined as follows [37]:where tetracyanoethylene (TCE) is taken as reference. All these properties were calculated using these equations for strictosidine in methanol through B3LYP/6-31G(d) and B3LYP/6-311++G(2d,p) basis sets and their values are shown in Table 2. Both HOMO and LUMO are bonding orbitals, resulting in a low excitation energy for strictosidine. The excitation energies, which were calculated as 4.451 eV for 6-31G(d) and 4.408 eV for 6-311++G(2d,p), reflect the low hardness value (2.22) showing strictosidine as a soft molecule with high polarizability. Electronegativity and electrophilicity values of strictosidine indicate that this molecule has significative attractive electron power acting as an electrophile and in addition to its polarization becomes very reactive since the electrons are farther from the nucleus; however, the nucleophilicity index value indicates that strictosidine is a strong nucleophile too. Such characteristic is justified because strictosidine is a large molecule that has many reactive groups with different potentials which forms small polarized points over its surface (as discussed in Section 3.6), which makes strictosidine a versatile molecule in view of the variety of alkaloids which it forms through intramolecular reactions with groups acting as nucleophiles and others as electrophiles [10]. On the other hand, the calculated chemical potential values, −3.06 for B3LYP/6-31G(d) and −3.40 for B3LYP/6-311++G(2d,p), reveal certain stability, indicating that it does not decompose spontaneously; that is, strictosidine molecule is reactive but does not tend to degrade into the components that formed it.

3.4. Tandem Mass Identification

The tandem mass spectra of protonated strictosidine (Figure S4 in Supplementary Data) of m/z 531 showed major fragment ions of m/z 514, 369, and 352. The ion of m/z 514 [M + H—17 Da] originated by loss of the NH₃ group from the tryptophane portion in a manner that was similar to the fragmentation of aporphine alkaloids (Pathway A) [38]. The ESI-MS³ of this ion generated fragment ions of m/z 352 [M + H—162 Da] due to the hydrogen rearrangement in the glycoside followed by heterolytic cleavage of the C21-O-C1′ bridge (Pathway C). The fragment ion of m/z 369 [M + H—162 Da] possibly arose directly from the heterolytic breakage of the glycoside portion of strictosidine via the cleavage of the C21-O-C1′ bridge (Pathway B), which occurred through a similar mechanism to Pathway C. No fragment ions were observed in the MS³ spectrum of the ion of m/z 369. Guided by the ESI-MS² and MS³ data, a fragmentation mechanism is proposed for protonated strictosidine in Figure 5.

3.5. NBO Study

An NBO analysis describes the Lewis-like molecular bonding pattern of electron pairs (or of individual electrons in the open-shell case) in the optimally compact form of the molecule. More precisely, NBOs are orthonormal sets of localized “maximum occupancy” orbitals whose leading N/2 members (or N members in the open-shell case) give the most accurate possible Lewis-like description of the total N-electron density. The Lewis-type NBOs determine the localized natural Lewis structure (NLS) representation of the wave function, while the remaining “non-Lewis-type” NBOs complete the span of the basis and describe the residual “delocalization effects” (i.e., departures from a single localized Lewis structure). NBOs provide therefore a valence bond-type description of the wave function that is closely linked to classical Lewis structure concepts [39–42].

NBO analysis is a helpful tool for understanding the delocalization of the electron density from the occupied Lewis-type (donor) NBOs to properly unoccupied non-Lewis-type (acceptor) NBOs [43–47] within the molecule. This analysis uses the second-order perturbation energies [donor () acceptor ()] that involve the most important delocalization instances, which are given as follows:

Table 3 shows more-detailed breakdown of the Lewis and non-Lewis occupancies, confirming the quality of the natural Lewis structure description. The total Lewis occupancy was 98.395% and the non-Lewis occupancy was 1.605%. The NBO analysis revealed strong intramolecular interactions formed by the orbital overlaps between C-C bonding and C-C antibonding and by overlaps between the N and O lone pairs (LP) and C-C antibonding. These interactions led to intramolecular electron-density transfer that caused the stabilization of the molecular system (Table 4). The intramolecular hyperconjugative interactions between C2, C7, C8, C9, C10, C11, C12, C13, and N1 revealed strong stabilization of the indolic portion, principally by interactions between C11-C10 C8-C9 (20.07 kcal/mol), C9-C8 C11-C10 (20.77 kcal/mol), C9-C8 C13-C12 (20.01 kcal/mol), and C13-C12 C8-C9 (19 kcal/mol) and electron donation from LP(1) N1 to C2-C7 (36.44 kcal/mol) and C12-C13 (40.91 kcal/mol). Other hyperconjugative interactions, especially the donation of electron density from LP(2) O3 to C16-C17 (37.05 kcal/mol), LP(1) O2 to C22-O1 (32 kcal/mol), and LP(2) O1 to C22-O2, gave strong stabilization to the dihydropyran ring of strictosidine. In addition, the interactions of C8-C9 C7-C2 and C21-O2 C16-C17 provided enormous stabilization of 191.70 kcal/mol and 62.11 kcal/mol, respectively. Table 4 provides all significant values for the hyperconjugative interactions given by the second-order perturbation theory.

A relationship between the ESI-IT-MS study and the stabilization caused by orbital overlap between bonds in the second-order perturbation theory could be established. The strongest stabilization energies for strictosidine involved hyperconjugative interactions in the indolic portion (aromatic and pentacyclic rings) and in the dihydropyran portion, explaining the small number of cleavages and the absence of breakages in the indolic portion and in the dihydropyran ring. In the fragmentation mechanism proposed in Figure 2, the most stable fragments arose from the heterolytic breakage of the O4-C1′, N4-C3, and N4-C5 bonds, which have the following weak hyperconjugative interactions: N4-C3 N2-C3 (, 3.59 kcal/mol), C5-N4 C3 (, 1.5 kcal/mol), N4-H C5 (, 1.68 kcal/mol), N4-H C5-H (, 2.51 kcal/mol), O4-C1′ C20-C21 (, 1.46 kcal/mol), and O3-C21 C21-O4 (, 0.51 kcal/mol). These results showed that the NBO study complemented the mass fragmentation study.

3.6. Molecular Electrostatic Potential Surface

Figure 6 illustrates 3D plots of the molecular electrostatic potential (MEP) surface for the strictosidine molecule. An MEP is a plot of the electrostatic potential mapped onto the constant electron density surface and is used primarily for predicting sites and relative reactivities towards electrophilic and nucleophilic attacks. MEPs are used in studies of biological recognition and interactions between the same molecules (e.g., in forming clusters and crystal structures) or other molecules. MEPs also correlate and predict a wide range of macroscopic properties [48, 49]. The color code of these maps ranges from −0.08 a.u. (deepest red) to 0.08 a.u. (deepest blue), where blue indicates a minimal concentration of electrons and red indicates a high density of electrons.

The MEPs for strictosidine indicated regions with positive potentials over H1 (0.0627 a.u.) and over the hydrogen atom of the OH group in position 4′ (0.0862 a.u.). Regions with negative potentials were located over the aromatic ring (−0.0429 a.u.), over the carbonyl group O2 (−0.0576 a.u.), between O4 and the OH group of C2′ position (−0.0643 a.u.), over OH on position C6′ (−0.0471 a.u.), and over oxygen atom on C4′ position (−0.0201 a.u.). The predominance of light green color region indicates great charge dispersion. Strictosidine is a large molecule with polarized points scattered over its surface which promotes various possible forms of intramolecular and intermolecular interaction (between strictosidine molecules); in addition, the nonflat shape hinders the chain interactions between strictosidine molecules which facilitate formation of a packed crystal.

3.7. IR Analysis

Figure 7 shows the experimental and theoretical IV spectra. The differences can be attributed to the fact that the theoretical DFT calculations were made for the molecule in the gas phase, whereas intermolecular interactions occur in solution. The assignment of the experimental bands to the normal modes of vibration was made using the optimized structure with the lowest potential energy, considering the potential energy distribution (PED) by using the B3LYP/6-31G(d) level. A total of 210 normal vibration modes were obtained but were compared with the experimental spectrum only between 400 and 4000 cm⁻¹ (Table 5).

Modes between 3021 and 3640 cm⁻¹ were related to the following vibrations: the strong O-H stretching of the glycoside portion; the H-C stretching of the aromatic ring; the H-C stretching of dihydropyran ring; the N4-H1 and N1-H1 stretching; and the H3-C3, H5-C5, H6-C6, H14-C14, and H15-C15 stretching. The region from 2900 to 3000 cm⁻¹ showed a strong band at 2932 cm⁻¹ that was assigned to the symmetric CH₂ stretching in C6.

The region between 1060 and 1700 cm⁻¹ was related to the following vibrations: the C=O stretching of C22, C=C stretching modes of C16=C17, the H1-N1 bending mode, and the CH₂ scissoring modes (C6′ and C23) in 1400–1690 cm⁻¹; the H-C bending of the glycoside portion and the dihydropyran ring in 1250–1399 cm⁻¹; the O-C stretching of the glycoside portion between 1060 and 1211 cm⁻¹.

The region between 600 and 1030 cm⁻¹ was related to the C-C stretching of C10 and C11, the C-O and C-C stretching of the glycosidic moiety, the H-N bending of the N4 position (band at 891 cm⁻¹ was noteworthy), the out-of-plane mode of C22 (at 763 cm⁻¹), and the torsion modes of the entire structure.

The large differences from 3500 to 4000 cm⁻¹ are related to the H-N1 and O-H stretching, which are indicative that the interactions between strictosidine molecules occur between the tryptophanic and glycosidic regions. These interactions make sense based on the electrostatic potential map (Figure 6), which showed greater polarization of these two regions relative to the entire molecule. The optimized geometry of a strictosidine dimer (Figure 8) showed stabilization due to the existence of intermolecular hydrogen bonds (N1-H1---O-C4′ and C4′-OH---O2=C22) and the values assigned to the stretching of O-H (3556.00 cm⁻¹) and H-N1 (3493.18 cm⁻¹) groups are closer to the experimental ones (3430 and 3390 cm⁻¹). The value assigned to the stretching of carboline group in position C22 (1690 cm⁻¹) shows to be closer to the experimental one in the dimer too (1767 cm⁻¹ to the monomer and 1716 cm⁻¹ to the dimer form), implying that the interaction between the carbonyl and OH is plausible and decreasing the stretching frequency oscillator related to C=O bond. These interactions directly influenced the infrared spectrum by decreasing the stretching frequency oscillator related to these bonds, causing reduction in the wavelengths (see Figure 8) and in the RMSD values. For the monomer, the RMSD is 84.40 cm⁻¹; for the dimer, the value is 46.41. Applying the empirical scaling factor of 0.9613, the RMSD values feature a visible reduction, 66.33 cm⁻¹ for monomer and 30.66 cm⁻¹ for dimer.

4. Conclusion

The strictosidine alkaloid, which was isolated from Strychnos amazonica, was comprehensively characterized. The interatomic distances and angles proved to be plausible compared to the X-ray data for similar molecules. The similarity between the theoretical and experimental coupling constants values reveals that the theoretical hydrogen dihedral angles of the C14-C15-C20-C21-C19 positions are plausible, showing that the modeled structure justifies the experimental NMR data. The UV analysis was able to explain the similarity between the UV spectra of strictosidine and related indole alkaloids, showing that the transitions involving the indole moiety are energetically more significant and such characteristic can be used as a “fingerprint” for detecting indole alkaloids. The HOMO-LUMO gap is directly related to the reactivity of a compound reflecting the amount of important properties such as chemical hardness, electrophilicity, nucleophilicity index, and electronegativity. Strictosidine theoretically appears to be a nucleophile and electrophile that in addition to its polarizability behaves as a soft molecule. This indicates low charge states and faster reactions, which makes strictosidine a versatile molecule, justifying its high reactivity and its role as a precursor of indole alkaloids. The comparative IR studies revealed that interactions of strictosidine dimers (between the tryptophanic and glycosidic regions) influenced the infrared spectrum by decreasing the stretching frequency oscillator of the groups which forms hydrogen bonds and revealed several characteristic vibrations that may be used as a diagnostic tool for other indole alkaloids, simplifying their identification and structural characterization. The NBO calculations showed that the strongest stabilization energies for strictosidine involved hyperconjugative interactions in the indolic portion (aromatic and pentacyclic rings) and in the dihydropyran portion, justifying the few fragmentation modes for its protonated molecule observed in the M analysis, and these modes may contribute to the further characterization of strictosidine analogues and derivatives.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank Dr. Doriane Picanço of the Federal University of Amazonas for the support with the botanical identification, Professor Dr. Michel Rossman for the X-ray data for loganin, Dr. Elaine H. T. Oliveira for the spelling support, and Capes and FINEP for the financial support.

Supplementary Materials

For a better comprehension of the data obtained for strictosidine, additional information not presented in the manuscript (NMR signals of strictosidine, ESI-IT-MS spectra of strictosidine, UV transitions and Mulliken charges) are available in the supplementary material. Table S1 contains experimental 13C, 1H and HMBC signals of the molecule under study. Figure S1 shows the experimental 1H strictosidine spectrum with integrated peaks. Figure S2 Provides the HSQC spectrum for the molecule title. Figure S3 Provides a copy of the HMBC spectrum. Figure S4 shows ESI-MS spectra with MS² and MS³ fragmentations. Figures S5 and S6 features the most theoretically contributory electronic transitions.

Supplementary Material

References

M. Frederick, M.-P. Hayette, M. Tits, P. De Mol, and L. Angenot, “In vitro activities of Strychnos alkaloids and extracts against Plasmodium falciparum,” Antimicrobial Agents and Chemotherapy, vol. 43, no. 9, pp. 2328–2331, 1999.
View at: Google Scholar
P. Rasoanaivo, S. Ratsimamanga-Urverg, R. Milijaona et al., “In vitro and in vivo chloroquine-potentiating action of Strychnos myrtoides alkaloids against chloroquine-resistant strains of Plasmodium malaria,” Planta Medica, vol. 60, no. 1, pp. 13–16, 1994.
View at: Publisher Site | Google Scholar
R. Bassleer, M. C. Depauw-Gillet, B. Massart et al., “Effets de trois alcaloïdes extraits du Strychnos usambarensis sur des cellules cancéreuses en culture,” Planta Medica, vol. 45, pp. 123–126, 1982.
View at: Publisher Site | Google Scholar
C. Caron, M. J. Hoizey, L. Le Men-Olivier et al., “Antimicrobial and antifungal activities of quasi-dimeric and related alkaloids,” Planta Medica, vol. 54, no. 5, pp. 409–412, 1988.
View at: Publisher Site | Google Scholar
J. Quetin-Leclercq, A. Favel, G. Balansard, P. Regli, and L. Angenot, “Screening for in vitro antifungal activities of some indole alkaloids,” Planta Medica, vol. 61, no. 5, pp. 475–477, 1995.
View at: Publisher Site | Google Scholar
C. L. C. de Medeiros, G. Thomas, and R. Mukherjee, “The source of Ca²⁺ for the spasmolytic actions of longicaudatine, a bisindole alkaloid isolated from Strychnos trinervis (Vell.) Mart. (Loganiaceae),” Phytotherapy Research, vol. 5, no. 1, pp. 24–28, 1991.
View at: Publisher Site | Google Scholar
V. K. Kapoor, S. K. Sharma, K. K. Chagti, and M. Singh, “Synthesis and hypotensive activity of diaboline,” Indian Journal of Chemistry B: Organic Chemistry Included Medicinal. Chemistry, vol. 27, pp. 641–644, 1988.
View at: Google Scholar
M. Tits, J. Damas, J. Quetin-Leclercq, and L. Angenot, “From ethnobotanical uses of Strychnos henningsii to antiinflammatories, analgesics and antispasmodics,” Journal of Ethnopharmacology, vol. 34, no. 2-3, pp. 261–267, 1991.
View at: Publisher Site | Google Scholar
C. L. Cardoso, D. H. S. Silva, M. C. M. Young, I. Castro-Gamboa, and V. S. Bolzani, “Indole monoterpene alkaloids from Chimarrhis turbinata DC Prodr.: a contribution to the chemotaxonomic studies of the Rubiaceae family,” Revista Brasileira de Farmacognosia, vol. 18, no. 1, pp. 26–29, 2008.
View at: Publisher Site | Google Scholar
L. F. Szabó, “Rigorous biogenetic network for a group of indole alkaloids derived from strictosidine,” Molecules, vol. 13, no. 8, pp. 1875–1896, 2008.
View at: Publisher Site | Google Scholar
X. Huang, Y. Li, Y. Su, X. Chai, and S. Yan, “Monoterpene indole alkaloids and monoterpene diglycosides from the roots of Triosteum pinnatifidum,” Phytochemistry Letters, vol. 7, no. 1, pp. 30–34, 2014.
View at: Publisher Site | Google Scholar
Y.-J. Xu, K. Foubert, L. Dhooghe et al., “Chromatographic profiling and identification of two new iridoid-indole alkaloids by UPLC-MS and HPLC-SPE-NMR analysis of an antimalarial extract from Nauclea pobeguinii,” Phytochemistry Letters, vol. 5, no. 2, pp. 316–319, 2012.
View at: Publisher Site | Google Scholar
M. Tits, V. Brandt, J.-N. Wauters, C. Delaude, G. Llabres, and L. Angenot, “Glucoindole alkaloids from stem bark of Strychnos mellodora,” Planta Medica, vol. 62, no. 1, pp. 73–74, 1996.
View at: Publisher Site | Google Scholar
G. N. Smith, “Strictosidine: a key intermediate in the biogenesis of indole alkaloids,” Chemical Communications, no. 15, pp. 912–914, 1968.
View at: Publisher Site | Google Scholar
R. Verpoorte, R. van der Heijden, and P. R. H. Moreno, The Alkaloids, vol. 49, Academic Press, San Diego, Calif, USA, 1997.
T. M. Kutchan, “Strictosidine: from alkaloid to enzyme to gene,” Phytochemistry, vol. 32, no. 3, pp. 493–506, 1993.
View at: Publisher Site | Google Scholar
A. R. Battersby, A. R. Burnett, E. S. Hall, and P. G. Parsons, “The rearrangement process in indole alkaloid biosynthesis,” Chemical Communications, vol. 24, pp. 1582–1583, 1968.
View at: Publisher Site | Google Scholar
A. R. Battersby, A. R. Burnett, and P. G. Persons, “Alkaloid biosynthesis. Part XV. Partial synthesis and isolation of vincoside and isovincoside: biosynthesis of the three major classes of indole alkaloids from vincoside,” Journal of the Chemical Society C: Organic, vol. 20, no. 8, pp. 1193–1200, 1969.
View at: Google Scholar
Á. Patthy-Lukáts, L. Károlyházy, L. F. Szabó, and B. Podányi, “First direct and detailed stereochemical analysis of strictosidine,” Journal of Natural Products, vol. 60, no. 2, pp. 69–75, 1997.
View at: Publisher Site | Google Scholar
Y. Yamazaki, A. Urano, H. Sudo et al., “Metabolite profiling of alkaloids and strictosidine synthase activity in camptothecin producing plants,” Phytochemistry, vol. 62, no. 3, pp. 461–470, 2003.
View at: Publisher Site | Google Scholar
M. J. Frisch, G. W. Trucks, H. B. Schlegel, and G. E. Scuseriaet, Gaussian 03 Revision, B.05, Gaussian, Pittsburgh, Pa, USA, 2003.
C. Timmons and P. Wipf, “Density functional theory calculation of ¹³C NMR shifts of diazaphenanthrene alkaloids: reinvestigation of the structure of samoquasine A,” Journal of Organic Chemistry, vol. 73, no. 22, pp. 9168–9170, 2008.
View at: Publisher Site | Google Scholar
J. A. Dobado, N. Benkadour, S. Melchor, and D. Portal, “Structure and theoretical NMR chemical shifts of modified cyclodextrins,” Journal of Molecular Structure: THEOCHEM, vol. 672, no. 1–3, pp. 127–132, 2004.
View at: Publisher Site | Google Scholar
R. Dennington, T. Keith, and J. Millam, Gaussview Version 5, Semichem, Shawnee Missions, Kan, USA, 2009.
M. H. Jamroz, Vibrational Energy Distribution Analysis, VEDA 4 Computer Program, Warsaw, Poland, 2004.
L. Dupont and O. Dideberg, “Structure cristalline et moleulaire d'un nouvel alcaloide indolique, l'akagerine,” Acta Crystallographica B, vol. 31, pp. 2378–2383, 1975.
View at: Publisher Site | Google Scholar
P. Molina, P. M. Fresneda, M. A. Sanz, C. Foces-Foces, and M. C. R. de Arellano, “Investigative studies on the formation of the imidazo [4′,5′:3,4] pyrido [2,3-b]indole ring: formal synthesis of the alkaloids grossularines-1 and 2. X-Ray crystal structures of 5-indol-3-yl-imidazole and bisimidazocarbazole derivatives,” Tetrahedron, vol. 54, pp. 9623–9638, 1998.
View at: Google Scholar
S. C. Nyburg, P. Y. Siew, G. N. Saunders, J. R. Purdy, and S. McLean, “X-ray crystal structure analysis of a bisepoxide derived from secologanin: the structure and absolute configuration,” Canadian Journal of Chemistry, vol. 61, no. 2, pp. 282–283, 1983.
View at: Publisher Site | Google Scholar
P. J. Lentz and M. G. Rossmann, “The crystal structure of loganin penta-acetate monomethyl ether bromide,” Journal of the Chemical Society D: Chemical Communications, no. 21, p. 1269, 1969.
View at: Publisher Site | Google Scholar
C. R. Hutchinson, A. H. Heckendorf, J. L. Straugh, P. E. Daddona, and D. E. Cane, “Biosynthesis of camptothecin. 3. Definition of strictosamide as the penultimate biosynthetic precursor assisted by carbon-13 and deuterium NMR spectroscopy,” Journal of the American Chemical Society, vol. 101, pp. 3358–3369, 1979.
View at: Google Scholar
A. W. Sangster and K. L. Stuart, “Ultraviolet spectra of alkaloids,” Chemical Reviews, vol. 65, no. 1, pp. 69–130, 1965.
View at: Publisher Site | Google Scholar
N. Özdemir, S. Dayan, O. Dayan, M. Dinçer, and N. Kalaycıoglu, “Experimental and molecular modeling investigation of (E)-N- ${2 - [(2 -hydroxybenzylidene) amino] phenyl}$ benzenesulfonamide,” Journal Molecular Physics, vol. 111, no. 6, pp. 707–723, 2013.
View at: Publisher Site | Google Scholar
E. I. Paulraj and S. Muthu, “Molecular structure analysis and spectroscopic characterization of 5-ethyl-5-phenyl-1,3-diazinane-4,6-dione with experimental (FT-IR and FT-Raman) techniques and quantum chemical calculations,” Spectrochimica Acta—Part A: Molecular and Biomolecular Spectroscopy, vol. 106, pp. 310–320, 2013.
View at: Publisher Site | Google Scholar
R. G. Parr and W. Yang, Functional Theory of Atoms and Molecules, Oxford University Press, New York, NY, USA, 1989.
R. G. Parr, L. V. Szentpály, and S. Liu, “Electrophilicity index,” Journal of the American Chemical Society, vol. 121, no. 9, pp. 1922–1924, 1999.
View at: Publisher Site | Google Scholar
L. R. Domingo and P. Pérez, “The nucleophilicity N index in organic chemistry,” Organic and Biomolecular Chemistry, vol. 9, no. 20, pp. 7168–7175, 2011.
View at: Publisher Site | Google Scholar
L. R. Domingo, P. Perez, and E. Chamorro.
C. Stévigny, J.-L. H. Jiwan, R. Rozenberg, E. De Hoffmann, and J. Quetin-Leclercq, “Key fragmentation patterns of aporphine alkaloids by electrospray ionization with multistage mass spectrometry,” Rapid Communications in Mass Spectrometry, vol. 18, no. 5, pp. 523–528, 2004.
View at: Publisher Site | Google Scholar
E. R. Davidson, Reduced Density Matrices in Quantum Chemistry, Academic Press, New York, NY, USA, 1976.
J. E. Carpenter and F. Weinhold, “Analysis of the geometry of the hydroxymethyl radical by the ‘different hybrids for different spins’ natural bond orbital procedure,” Journal of Molecular Structure: THEOCHEM, vol. 169, pp. 41–62, 1988.
View at: Publisher Site | Google Scholar
B. C Carlson, J. M. Keller, and B. C Carlson, “Orthogonalization Procedures and the Localization of Wannier Functions,” Physical Reviews, vol. 105, p. 102, 1957.
View at: Google Scholar
A. E. Reed, R. B. Weinstock, and F. Weinhold, “Natural population analysis,” The Journal of Chemical Physics, vol. 83, no. 2, pp. 735–746, 1985.
View at: Publisher Site | Google Scholar
R. S. Mulliken, “Electronic population analysis on LCAO-MO molecular wave functions,” Journal of Chemical Physics, vol. 23, pp. 1833–1841, 1955.
View at: Google Scholar
R. S. Mulliken and W. C. Ermler, Diatomic Molecules: Results of Ab Initio Calculations, Academic Press, New York, NY, USA, 1977.
J. K. Badenhoop and F. Weinhold, “Natural bond orbital analysis of steric interactions,” Journal of Chemical Physics, vol. 107, no. 14, pp. 5406–5421, 1997.
View at: Publisher Site | Google Scholar
J. P. Foster and F. Weinhold, “Natural hybrid orbitals,” Journal of the American Chemical Society, vol. 102, no. 24, pp. 7211–7218, 1980.
View at: Publisher Site | Google Scholar
A. E. Reed and F. Weinhold, “Natural bond orbital analysis of near-Hartree-Fock water dimer,” The Journal of Chemical Physics, vol. 78, no. 6, pp. 4066–4073, 1983.
View at: Publisher Site | Google Scholar
V. Deval, A. Kumar, V. Gupta et al., “Molecular structure (monomeric and dimeric) and hydrogen bonds in 5-benzyl 2-thiohydantoin studied by FT-IR and FT-Raman spectroscopy and DFT calculations,” Spectrochimica Acta—Part A: Molecular and Biomolecular Spectroscopy, vol. 132, pp. 15–26, 2014.
View at: Publisher Site | Google Scholar
J. S. Murray and K. Sen, Molecular Electrostatic Potentials Concepts and Applications, Elsevier Science, Amsterdam, The Netherlands, 1996.

Copyright

Copyright © 2016 Renyer Alves Costa et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

4145

Downloads

1317

Citations