Next Article in Journal
Highly Effective Removal of Ofloxacin from Water with Copper-Doped ZIF-8
Next Article in Special Issue
Multispectral and Molecular Docking Studies Reveal Potential Effectiveness of Antidepressant Fluoxetine by Forming π-Acceptor Complexes
Previous Article in Journal
Effect of Yakae-Prajamduen-Jamod Traditional Thai Remedy on Cognitive Impairment in an Ovariectomized Mouse Model and Its Mechanism of Action
Previous Article in Special Issue
Screening of Specific and Common Pathways in Breast Cancer Cell Lines MCF-7 and MDA-MB-231 Treated with Chlorophyllides Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spectroscopic and Molecular Docking Studies of Cu(II), Ni(II), Co(II), and Mn(II) Complexes with Anticonvulsant Therapeutic Agent Gabapentin

1
Department of Chemistry, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
2
Department of Biology, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
3
Centre of Biomedical Sciences Research (CBSR), Deanship of Scientific Research, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
4
Department of Pharmacology and Toxicology, College of Pharmacy, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
5
Addiction and Neuroscience Research Unit, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
6
College of Medicine, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
7
Department of Clinical Laboratories Sciences, the Faculty of Applied Medical Sciences, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
8
Department of Chemistry, Faculty of Science, Aligarh Muslim University, Aligarh 202002, India
9
Department of Chemistry, Faculty of Science, Zagazig University, Zagazig 44519, Egypt
*
Authors to whom correspondence should be addressed.
Molecules 2022, 27(13), 4311; https://doi.org/10.3390/molecules27134311
Submission received: 9 May 2022 / Revised: 30 June 2022 / Accepted: 1 July 2022 / Published: 5 July 2022
(This article belongs to the Special Issue Metal-Based Drugs: Past, Present and Future)

Abstract

:
New Cu(II), Ni(II), Co(II), and Mn(II) complexes of the gabapentin (Gpn) bidentate drug ligand were synthesized and studied using elemental analyses, melting temperatures, molar conductivity, UV–Vis, magnetic measurements, FTIR, and surface morphology (scanning (SEM) and transmission (TEM) electron microscopes).The gabapentin ligand was shown to form monobasic metal:ligand (1:1) stoichiometry complexes with the metal ions Cu(II), Ni(II), Co(II), and Mn(II). Molar conductance measurements in dimethyl-sulfoxide solvent with a concentration of 10−3 M correlated to a non-electrolytic character for all of the produced complexes. A deformed octahedral environment was proposed for all metal complexes. Through the nitrogen atom of the –NH2 group and the oxygen atom of the carboxylate group, the Gpn drug chelated as a bidentate ligand toward the Mn2+, Co2+, Ni2+, and Cu2+ metal ions. This coordination behavior was validated by spectroscopic, magnetic, and electronic spectra using the formulas of the [M(Gpn)(H2O)3(Cl)]·nH2O complexes (where n = 2–6).Transmission electron microscopy was used to examine the nanostructure of the produced gabapentin complexes. Molecular docking was utilized to investigate the comparative interaction between the Gpn drug and its four metal [Cu(II), Ni(II), Co(II), and Mn(II)] complexes as ligands using serotonin (6BQH) and dopamine (6CM4) receptors. AutoDock Vina results were further refined through molecular dynamics simulation, and molecular processes for receptor–ligand interactions were also studied. The B3LYP level of theory and LanL2DZ basis set was used for DFT (density functional theory) studies. The optimized geometries, along with the MEP map and HOMO → LUMO of the metal complexes, were studied.

1. Introduction

GABAPENTIN (Gpn; Figure 1), is sold under the brand name neurontin, and is commonly named as 2-[1-(aminomethyl)cyclohexyl]acetic acid. Gpn is a structural analogue of the neurotransmitter gamma aminobutyric acid, which has mostly been investigated for its inhibitory effect on the central nervous system [1]. Since Gpn has both acidic groups (COOH) and a basic group (NH2), it is an artificial amino acid [1,2]. It is a antiepileptic drug that is utilized as both a supplement and a stand-alone treatment for people suffering from partial seizures [2]. Gpn has also been used in the treatment of neuropathic pain. For a growing number of individuals with epilepsy, a safe and effective seizure is a major concern. Epilepsy has a significant economic impact [3]. The effectiveness, tolerability, and safety of an antiepileptic medicine are all factors to consider when choosing one. Gpn can also be regarded as an emergent solution to the “pain puzzle”. Double-blind and more randomized studies comparing analgesic medications to Gpn may be useful in determining the first-line treatment for chronic and acute pain relief [4]. As a result, finding compounds to treat epilepsy is necessary [5,6]. Metal-binding or metal-recognition sites are found in a variety of drugs and potential pharmacological agents. These sites can bind or interact with metal ions, influencing their bioactivities and possibly causing damage to their target biomolecules. The literature [7,8,9,10] has numerous examples of these “metallodrugs” and “metallopharmaceuticals” and their actions. Metals and metal complexes have played an important part in the evolution of contemporary chemotherapy. Anticancer platinum drugs, for example, are used in more chemotherapy regimens than any other class of anticancer drugs, and they have played a significant role in the success of cancer treatment over the last three decades [11].
After coordinating to a metal, ligands can play an essential role in changing the pharmacological characteristics of existing medications. Because the resulting prodrugs have varied physical and pharmacological properties, they can be delivered in a regulated manner or at a specific site [7]. This method can be used to save medications that have failed due to poor pharmacology or high toxicity. The complexation of nonsteroidal anti-inflammatory medications to copper, for example, eliminates some of the drugs’ gastrointestinal adverse effects. In the hypoxic zones of solid tumors, the release of cytotoxins such as nitrogen mustards from redox-active metals such as cobalt has the potential to improve medication effectiveness and reduce toxicity [8].
Metal complexes have received less attention in medical chemistry than organic molecules [9]. Many organic chemicals employed in medicine, in reality, do not have a wholly organic mechanism of action and require the residues of metal ions for activation or biotransformation, either directly or indirectly. The status of metal ions and their complexes with biomolecules in the body affects our health, aging, physiological problems, and diseases. Metals make up about 0.03% of the weight of the human body. The concentrations of Cd, Cr, Ti, V, Cu, Se, and Zn in malignant sections of the kidney were shown to be lower than in noncancerous parts [10]. Coordination bonds can be formed by ligands with electron donor atoms such as N, O, S, and P. Chelation alters the biological properties of both ligands and metal moieties, and in many circumstances, it produces a synergistic interaction between the metal ion and the ligand [10].
Mechanosynthesis has been used to create a variety of Gpn coordination networks with Mn2+, Y3+, La3+, Er3+, Nd3+, and Ce3+ [11,12]. In the literature survey, there appears to be a lack of study on the chelation behavior of gabapentin and its derivatives toward metal ions [13,14]. Therefore, in this paper, the experimental plan aimed to synthesize some new Cu(II), Ni(II), Co(II), and Mn(II) complexes of gabapentin by refluxing some transition metal(II) chlorides with a basic solution of gabapentin in a 1:1 molar ratio in a methanol solvent. The structure of the obtained complexes was studied using spectral, magnetic, and morphological techniques. AutoDock Vina software was used to study the interactions between the receptors (serotonin and dopamine) and ligands [gabapentin (Gpn) and its four synthesized metal complexes ([Cu(II)–(Gpn)], [Ni(II)–(Gpn)], [Co(II)–(Gpn)], and [Mn(II)–(Gpn)]) were studied theoretically by the molecular docking method. For more clarity, the binding energy along with the aromatic, hydrogen bond, hydrophobic, SAS, ionizability, and interpolated charge surfaces at the interaction site were also generated. The molecular docking outputs were further investigated for receptor–ligand interactions through MDS (molecular dynamic simulation), which used a 100 ns run at 300 K. Structural stability, hydrogen bond interactions, SASA, compactness of structure, and residue flexibility of the complexes were studied to compare their dynamic features. The optimized geometries of ([Cu(II)–(Gpn)], [Ni(II)–(Gpn)], [Co(II)–(Gpn)], and [Mn(II)–zx(Gpn)] were obtained through the B3LYP level of theory and LanL2DZ basis set of the DFT calculations. Important parameters such as chemical, structural, and spectroscopic properties of the [Cu(II)–(Gpn)] complex were also obtained.

2. Materials and Methods

2.1. Chemicals

Metal chloride salts (MnCl2.4H2O (≥99%), CoCl2.6H2O (98%), CuCl2.2H2O (≥99%), and NiCl2.6H2O (≥98%)) and gabapentin (≥99%) were obtained from Sigma-Aldrich, USA and utilized in the preparation as received.

2.2. Instruments

The following is a list of the different types of analysis and their respective models:
Type of AnalysisModels
SEMQuanta FEG 250 equipment
TEMJEOL 100s microscopy
Magnetic momentMagnetic Susceptibility Balance
Electronic spectraUV2 Unicam UV/Vis Spectrophotometer
FTIR spectraBruker FTIR Spectrophotometer
ConductanceJenway 4010 conductivity meter
Elemental analyses (C,H,N)Perkin Elmer CHN 2400
Metal ions (Mn, Co, Cu, Ni)An atomic absorption spectrometer model PYE-UNICAM SP 1900

2.3. Synthesis of Gpn Metal Complexes

2.3.1. Synthesis of the [Mn(Gpn)(H2O)3(Cl)]·4H2O Complex

The synthesis of the manganese(II)complex was carried out by adding MnCl2·4H2O (1.98 g, 1.0 mmol) to the basic solution of the gabapentin ligand (Gpn) (1.72 g, 1.0 mmol) by neutralization (pH = ~7.0) using NH3 solution. The mixture was stirred at 80 °C for 3 h. The dark brown precipitate of the manganese(II) complex was filtered and washed with distilled water three times, dried, and kept in desiccators. Yield: 77%. C9H30ClMnNO9, Mwt. 386.73 g/mol (Found: C 27.87; H 7.76; N 3.59; Mn 14.08; Mwt. 386.73, calcd. C 27.95; H 7.82; N 3.62; Mn 14.21; m.p. > 234 °C).

2.3.2. Synthesis of [Co(Gpn)(H2O)3(Cl)]·6H2O Complex

The synthesis of the cobalt(II) complex was carried out by adding CoCl2·6H2O (2.38 g, 1.0 mmol) to the basic solution of the gabapentin ligand (Gpn) (1.72 g, 1.0 mmol) by neutralization (pH = ~7.0) using the NH3 solution. The mixture was stirred at 80 °C for 2 h. The dark blue precipitate of the cobalt(II) complex was filtered and washed with distilled water three times, dried, and kept in desiccators. Yield: 70%. C9H34ClCoNO11, Mwt. 426.75 g/mol (Found: C 25.23; H 7.98; N 3.24; Co13.76; Mwt. 426.75, calcd. C 25.33; H 8.03; N 3.28; Co13.81; m.p. > 237 °C).

2.3.3. Synthesis of [Cu(Gpn)(H2O)3(Cl)]·2H2O Complex

The synthesis of copper(II) complex was carried out by adding CuCl2·2H2O (1.71 g, 1.0 mmol) to the basic solution of gabapentin ligand (Gpn) (1.72 g, 1.0 mmol) by neutralization (pH = ~7.0) using the NH3 solution. The mixture was stirred at 80 °C for 2 h. The dark brown precipitate of the copper(II) complex was filtered and washed with distilled water three times, dried, and kept in desiccators. Yield: 72%. C9H26ClCuNO7, Mwt. 359.30 g/mol (Found: C 30.03; H 7.16; N 3.87; Cu17.54; Mwt. 359.30, calcd. C 30.08; H 7.29; N 3.90; Cu17.69; m.p. > 232 °C).

2.3.4. Synthesis of [Ni(Gpn)(H2O)3(Cl)]·3H2O Complex

The synthesis of the nickel(II) complex was carried out by adding NiCl2·6H2O (2.38 g, 1.0 mmol) to the basic solution of the gabapentin ligand (Gpn) (1.72 g, 1.0 mmol) by neutralization (pH = ~7.0) using the NH3 solution. The mixture was stirred at 80 °C for 2 h. The green precipitate of the nickel(II) complex was filtered and washed with distilled water three times, dried, and kept in desiccators. Yield: 75%. C9H28ClNNiO8, Mwt. 372.47 (Found: C 29.00; H 7.46; N 3.71; Ni 17.65; Mwt. 372.47, calcd. C 29.02; H 7.58; N 3.76; Ni 15.76; m.p. > 238 °C).

2.4. Computational Studies

2.4.1. Molecular Docking

All of the docking was conducted through an Intel(R) Core(TM) i5-4200U CPU-2.10 GHz, 64 bit processor. AutoDock Vina [15] program was used to perform the docking. The receptors serotonin and dopamine with PDB ID: 6BQH and 6CM4, respectively, were obtained through an online PDB bank [16]. Serotonin and dopamine were set up by cleaning the ligand and water molecules through Discovery Studio (DS) software (https://www.3ds.com/products-services/biovia/, accessed on 1 Febrauary 2022). DTool [17] was used to add polar hydrogen atoms and Kollman charges, while the Geistenger method was used to assign the partial charges. The ligands [gabapentin (Gpn) and its four complexes—([Mn(II)–(Gpn)], [Co(II)–(Gpn)], [Ni(II)–(Gpn)], and [Cu(II)–(Gpn)])]—in the PDBQT format were obtained using the OpenBabelIGUI program [18]. The energy of all of the ligands was minimized through PyRx-Python at 200 steps using the conjugate gradient optimization algorithm and MMFF94 force field [19]. Discovery Studio software was used to analyze the docking results.

2.4.2. Molecular Dynamics Simulation (MDS) Studies

To conduct the simulation, an “Intel(R) Xeon(R) CPU E5-2680 v4-2.40GHz, 64 bit” processor was used. MDS was performed using GROMACS version 2019.2 with a GROMOS96-43a1 force field. To initiate the MDS analysis, results obtained from molecular docking for Gpn and [Cu(II)–(Gpn)] with a high docking score were used. Topology and parameter files were created through CHARMM-GUI with the latest CGenFF [20,21]. To solve receptor–ligand structures in a triclinic box, SPC water models were used [22]. Twenty-seven Cl- and 28 Na+ ions were added for neutralization at 0.15 M salt (Figure S1). A Leap-frog MD integrator for a 100 ns simulation time in the NPT/NVT equilibration run systems was exposed to the periodic boundary conditions (300 K and 1.0 bar) [23,24]. The RMSD (root mean square deviation), hydrogen bonding, SASA, and gyration radius were examined using gmx rms, gmxhbond, gmxsasa, and gmx gyrate tools, respectively [25]. PyMol/VMD and Grace Software [12,13] were used to prepare the different plots.

2.4.3. Density Functional Theory (DFT)

The Gaussian 09RevD.01 package [26] was used for the DFT/TD-DFT calculations. The molecules were optimized in their ground spin state B3LYP density functional to know the exact position of the hydrogen atoms, as B3LYP, a hybrid method containing elements from DFT and the Hartree–Fock theory, was the first DFT exchange correlation functional to convince computational chemists that DFT could predict molecule physicochemical properties and reaction barriers with an accuracy comparable to some wave function-based methods but with much improved computational efficiency. The Los Alamos Effective Core Potentials lanL2DZ basis set was employed for the Mn, Co, Ni, and Cu atom while the split-valence 6-31G(d) basis set was applied for the other atoms [27]. A molecular electrostatic potential map (MEP) of the synthesized metal complexes were studied [28]. The IR frequencies were computed with full accuracy using animated modes of vibrations. To establish a good comparison between the experimental and theoretical wavenumbers, the scaling factor of 0.8522 was used. Moreover, structure-based molecular properties such as bond lengths, bond angles, atomic charges, total energy, electronic properties, and frontier molecular orbitals energy were calculated by this theory in the gas phase. ChemCraft 1.5 software [29] was used for visualization.

3. Results and Discussion

Herein, our paper aimed to synthesize a complexation between the gabapentin pure drug and four transition metal ions with a 1:1 molar ratio in alkaline media and check the suggested structures of the resulting complexes by spectroscopic characterizations. The difference between our synthesized complexes and others mentioned in the literature is the isolated 1:1 molar ratio complexation and the isolation at pH = 7. The articles mentioned in the literature dealing with the mixed gabapentin ligands and the isolation of the 1:2 molar ratio were different from our synthesized complexes. We looked at a lot of crystallites but could not find one that was adequate for full data capture. As a result, a complete dataset of nickel(II), cobalt(II), copper(II), and manganese(II) gabapentin complexes was not collected, and the discussion was based on the elemental analysis, melting points, molar conductivity, UV–Vis, magnetic measurements, FTIR, and surface morphology (scanning (SEM) and transmission (TEM) electron microscopes) as well as the theoretical study.

3.1. Elemental and Conductance Measurements

The gabapentin metal complexes of nickel(II), cobalt(II) copper(II), and manganese(II) are commonly soluble in DMSO (dimethyl sulfoxide) and DMF (dimethylformamide). The synthesized complexes had a stoichiometry ratio of 1:1 (metal:ligand), according to the elemental analytical data; the postulated structure of the complexes is well-supported by analytical data and other spectrum studies (Figure 2). The color, yield, and melting point data of all the compounds were introduced in experimental Section 2.3. In DMSO, the molar conductance of the solutions of all the complexes were in the range of 11–14 Ω−1cm2 mole−1. These experimental data indicate that complexes are non-electrolytes [30] in DMSO (10−3 M) at room temperature. Metal ions and Gpn ligand reactions with a 1:1 molar ratio yielded bidentate complexes. The Gpn ligand on reaction with Cu(II), Ni(II), Co(II), and Mn(II) salt yielded complexes that corresponded to the formulas [Cu(Gpn)(H2O)3(Cl)]·2H2O, [Ni(Gpn)(H2O)3(Cl)]·3H2O, [Co(Gpn)(H2O)3(Cl)]·6H2O, and [Mn(Gpn)(H2O)3(Cl)]·4H2O. The general compositions for the complexes are [M(Gpn)(H2O)3(Cl)]·nH2O complexes (where n = 2–6). Elemental analyses of the prepared compounds indicate the formation of six coordination complexes containing one molecule of gabapentin as a bidentate ligand, one chlorine atom, and three water molecules. The absence of an absorption band at ca. 3400 cm−1 in the IR spectrum of the Gpn ligand shows that free amino groups were absent. A broad band was observed at ca. 1650–1680 cm−1, which indicates the presence of water molecules from the complexes.
The FTIR spectra of the pure gabapentin and each of the Cu(II), Ni(II), Co(II), and Mn(II)–Gpn complexes were also obtained (Figure 3 and Table 1).
Regarding the infrared spectrum of the free gabapentin ligand, it did not show a peak of the stretching –NH within the region (3500–3300 cm−1), since it is a zwitterion in the solid state [16,17]. The two main characteristic bands at 2928 and 2861 cm−1 were assigned to the –NH3+ stretching vibrations [31,32]. The peak at 2150 cm−1 was for the CN group side chain stretching vibration [32]. At 1614 and 1545 cm−1, the peaks represented the stretching asymmetric vibration and vibrations of the NH3+ bending deformation of the carboxylate group, respectively [32]. In the 400–4000 cm−1 region, the FTIR spectra of the gabapentin complexes exhibited several bands. The δ(NH2) bond vibration appeared at 1552–1568 cm−1 for the copper complex, which showed a shift to lower frequencies to free Gpn [33]. To identify the type of carboxylate binding to a transition metal ion, the value Δν = [νas(COO) − νs(COO)] was employed. Generally, unidentate carboxylation is indicated by the difference in Δν between the symmetric and asymmetric (COO) absorption frequencies of >200 cm−1. In the nickel(II), cobalt(II), copper(II), and manganese(II) complexes, the Δν values were 283, 274, 270, and 240, cm−1, indicating monodendate coordination of the carboxylate group [34]. In these synthesized complexes, it can be attributed to the Gpn ligand, which is linked to the metal ion by carboxylic acid atoms and the nitrogen of amino groups. Furthermore, the existence of ν(M–O) and ν(M–N) medium intensity bands at 612–620 cm−1 and 468–519, respectively, suggested the generation of Gpn complexes [35]. In the FTIR spectra, the characteristic bands of bending to water hydrated, δ(H2O) at ca. 1650 cm−1 and broad absorption bands of (OH) with maximum at 3404–3507 cm−1 in synthesized metal complexes confirmed the presence of coordinated water molecules.

3.2. Electronic Spectra, and Magnetic Susceptibility Studies

3.2.1. Manganese(II) Complex

The molar conductance measurements showed that the Mn(II) complex is a non-electrolyte. Thus, it is formulated as [Mn(Gpn)(H2O)3(Cl)] 4H2O. It showed a magnetic moment corresponding to five unpaired electrons (5.87 B.M.) at room temperature, which is close to the spin only value (5.92 B.M.). The electronic spectrum of this complex displayed weak absorption bands at 568 nm (ν1), 388 nm (ν2), 378 nm (ν3), and 369 nm (ν4), characteristic of octahedral geometry [36]. These bands may be assigned as 6A1g4T1g, 6A1g4Eg, 4A1g, 6A1g4Eg and 6A1g4T1g transitions, respectively.

3.2.2. Cobalt(II) Complex

The molar conductance measurements showed that the complex is also a non-electrolyte. Thus, it may be formulated as [Co(Gpn)(H2O)3(Cl)] 6H2O. At room temperature, it showed a magnetic moment of 4.90 B.M., a value in tune with a high-spin configuration showing the presence of a distorted octahedral environment around the cobalt(II) ion in the complex. The electronic spectrum of the complex displayed absorption bands at 610 and 542 nm. A distorted octahedral geometry, this complex is indicated through the investigation of this electronic spectral data [37]. The assignment of the spectral bands may be given as 4T1(F) → 4A2(F) and 4T1(F) → 4T1(P) transitions, respectively.

3.2.3. Copper(II) Complex

The molar conductance measurements showed that the complex is a non-electrolyte and that it may be formulated as [Cu(Gpn)(H2O)3(Cl)] 2H2O. The magnetic moment of the complex was 1.96 B.M. The six-coordinate copper(II) complex had either D4h or C4v symmetry, and the Eg and T2g levels of the 2D free ion split into B1g, A1g, B2g, and Eg levels, respectively. Thus, three spin-allowed transitions are expected in the visible and near-IR region, but only a few complexes are known in which such bands are resolved [38]. These bands were assigned to the following transitions, in order of increasing energy, 2B1g2A1g, 2B1g2B2g, and 2B1g2Eg. The energy level sequence will depend on the amount of tetragonal distortion due to ligand-field and Jahn–Teller effects [39]. The electronic spectrum of the complex showed a band at 722 nm and a well-defined shoulder at 620 nm, which may be assigned to the 2B1g2A1g and 2B1g2Eg transitions, respectively.

3.2.4. Nickel(II) Complex

Nickel(II) has a (d8) configuration with the following orgel diagram found at 780 nm, 398 nm, and 384 nm, respectively, for the nickel complex. Based on the orgel diagram, the first peak at 697 nm was assigned to 3A2g3T2g, the second at 346 nm for 3A2g3T1g (F), and the third at 315 nm was due to the 3A2g3T1g (P) transition, respectively. This confirms the presence of an octahedral geometry for the nickel complex. The magnetic moment of the complex was 3.28 B.M. at room temperature, a value in tune with a high-spin configuration, showing the presence of a distorted octahedral environment around the nickel(II) ion in the complex. On the basis of magnetic, molar conductivity, molar ratio, UV/Vis., and the IR results, the proposed molecular structure of the synthesized complexes is octahedral (Figure 3) [37].

3.3. Scanning and Transmission Electron Microscopes

The SEM analysis was carried out to check the surface morphology of the selected complexes, and the micrographs obtained are given in Figure 4. The micrographs of the manganese(II), copper(II), cobalt(II), and nickel(II)–gabapentin complexes are given in Figure 5I–IV; it can be seen that these gave an appearance of pressed chips (a wooden board has several layers), a coral reef-like, a block of rock, and rectangular panels of different sizes, respectively.
The TEM image (Figure 5I–IV) shows that the Mn(II)–Gpn, Cu(II)–Gpn, Co(II)–Gpn, and Ni(II)–Gpn complexes were aggregates with an irregular shape. The particle sizes of the [Mn(Gpn)(H2O)3(Cl)] 4H2O, [Co(Gpn)(H2O)3(Cl)]·6H2O, [Cu(Gpn)(H2O)3(Cl)] 2H2O, and [Ni(Gpn)(H2O)3(Cl)] 3H2O complexes were calculated from the TEM images with the ranges of 13–27 nm, 12–18 nm, 10–18 nm, and 10–20 nm, respectively (Figure 5).

3.4. Molecular Docking Studies

The interactions between the synthesized metal complexes—([Cu(II)–(Gpn), [Ni(II)–(Gpn)], [Co(II)–(Gpn)], and [Mn(II)–(Gpn)]]—with the prepared serotonin and dopamine receptors were studied, and the best docking poses were analyzed. For comparative purposes, the Gpn drug was used as the control. The results revealed that the potential binding energy of all of the metal complexes was higher than Gpn in both the receptors (Table 2).
Among the four metal complexes screened, [Cu(II)–(Gpn)] showed the highest docking energy. Docking of [Cu(II)–(Gpn)] with serotonin and dopamine gave the potential binding energy of −7.2 and −6.5 kcal/mol, respectively. The greater value of the binding energy in the case of [Cu(II)–(Gpn)]-serotonin (CuGS) suggests a stronger interaction than with dopamine. The docking data are given in Table 3, and the best docking pose (CuGS) is shown in Figure 6.
An illustration of the molecular docking for the interactions of the ligands and receptor is depicted in Figure 7a,b.
As shown in Figure 8a, CuGS showed that the residues Arg173, Thr109, His182, and Asn187 had established hydrogen bonds. Additionally, Asp172 (Attractive charge) and Ala176 and Arg173 (Alkyl) interactions could also be seen [40,41]. Similarly, docking of the Gpn drug with serotonin and dopamine gave the potential binding energy of −5.1 and −4.8 kcal/mol, respectively. The greater value of the binding energy in the case of Gpn-serotonin (GpnS) suggests a stronger interaction than that with dopamine. Figure 8b shows the interaction between Gpn and serotonin, revealing that the residues Thr109 and Ala108 had established hydrogen bonds. Additionally, Ala321 showed alkyl interactions. These findings suggest that the metal [Cu(II)–(Gpn)] binds with the receptors more efficiently than the Gpn drug. The 2D representations of the interactions between the metal complex/drug and receptor are shown in Figure 8. The other binding details of the interactions are tabularized in Table 4 and Table 5.

3.5. Solvent Accessible Surface, Aromatic, Ionizability, Hydrophobicity, and Hydrogen Bond Surfaces

Aromatic, hydrogen bond, hydrophobic, solvent accessible surface (SAS), ionizability, and interpolated charge surfaces at the interaction site (Figure 9 and Figure S2) were studied through DS software [42]. The aromatic shown in Figure S2a shows the face surface as an orange color and the edge surface as a blue color. The green color in the hydrogen bond surface represents the acceptor area and the pink color is the donor area of the amino acid residues (Figure 9). The presence of hydrophilicity features of the receptor around the ligand was confirmed by the hydrophobicity surface (Figure S2b). The solvent accessibility surface (SAS) is the surface area of the receptor (Figure S2c) that is reachable to a solvent [43], the green color is the poorly accessible area, and the blue color is the highly accessible area [43]. The ionization surface reflects the acidic and basic propensity (Figure S2d, blue color = basic and red color = acidic).

3.6. MDS Analysis

The top docking results (CuGS and GpnS) were used for the MDS analysis of the 100 ns run. After studying the RMSD, it was observed that both CuGS and GpnS established a stable conformation after ~30 and ~75 ns, respectively, having appropriate values of RMSD of 2.15 and 2.75 Å, respectively (Figure 10).
The average distance and standard deviation for all amino acid pairs between the two conformations were investigated using RR distance maps, which are two-dimensional representations of a protein’s 3D structure [46]. Figure 12 represents the RR distance maps, which plot patterns of spatial interactions [47,48].
The zero distance between two residues is shown by the white diagonal on the map, whereas the red and blue components reflect residue pairs with the highest distance differences in the two conformations. The average radius of gyration (Rg) value of 25.75 and 26.75 Å was observed for GpnS and CuGS, respectively. It can be seen that over the simulation time, the Rg for CuGS decreased, indicating that the structures became more compact (Figure 13).
The grid search at 25 × 11 × 14 = grid and rcut = 0.35 revealed the H-bond between the receptor and ligand (GpnS and CuGS), plotted against time (Figure 14). On calculating the hydrogen bonds between the ligand (Gpn and [Cu(II)–(Gpn)]) and protein (3706 atoms), 508 donors and 987 acceptors were observed. The number of H-bonds per time frame on average were found to be 2.107 and 1.104 out of a possible 250,698 for CuGS and GpnS, respectively. Overall, it was discovered that the protein–receptor interaction increased the number of hydrogen bonds substantially, which was more in CuGS. Figure 15 shows that the SASA (solvent accessibility surface area) values changed due to the binding of the ligand (Gpn and [Cu(II)–(Gpn)]) to the receptor (serotonin). The decreased SASA value indicates the reduced pocket size and alteration in the protein structure conformation with increased hydrophobicity.

3.7. DFT Calculations

The structures of all four synthesized metal complexes ([Mn(II)–(Gpn)], [Co(II)–(Gpn)], [Ni(II)–(Gpn)], and [Cu(II)–(Gpn)]) were optimized using the B3LYP level of theory and LanL2DZ basis set. The minimum SCF energy after 33, 22, 22, and 40 optimization steps was found to be −906.0136, −937.0486, −961.2607, and −988.0661 a.u., respectively, for ([Mn(II)–(Gpn)], [Co(II)–(Gpn)], [Ni(II)–(Gpn)], and [Cu(II)–(Gpn)], respectively. The optimized geometry of the four metal complexes is given in Figure 16.
The strength of the electrostatic potentials for ([Mn(II)–(Gpn)], [Co(II)–(Gpn)], [Ni(II)–(Gpn)], and [Cu(II)–(Gpn)] was represented through the MEP map (Figure 17), in order to investigate the most electron rich and poor regions and to rationalize the noncovalent interactions. The electropositive regions are displayed in a blue color and the electronegative in a red color [49]. It was found that the area around the Cl and O atoms had strong negative electrostatic potential, and strong positive electrostatic potential could be seen around H2O, which shows the preferential binding sites over the molecule. The MEP surface map is represented in the color scale from deep red to deep blue [50]. From the optimized structure of [Cu(II)–(Gpn)], the bond lengths and bond angles were obtained (Figure 18, Table 6). The Mulliken charges for [Cu(II)–(Gpn)] were also calculated (Table 7). The shift in the Mulliken charges were observed in atoms of the synthesized metal complexes from the reactant moieties, which suggests the formation of different complexes. The result of the MEP is in agreement with the Mulliken charges.
The IR frequencies were investigated through B3LYP/LanL2DZ (data not shown). The slight difference between the experimental and theoretical frequencies was due to the experimental values being obtained in the solid phase and theoretical values being acquired in the gas phase. The calculated vibrational frequencies varied to a smaller extent from the experimental results due to neglecting the incompleteness and anharmonicity of the basis set [51]. TD-DFT method was used to explore the nature of the electronic transitions in [Cu(II)–(Gpn)] in the gas phase. One broad electronic absorption band was obtained from TD-DFT at 569 nm. The HOMO (−4.2308 eV) to LUMO (−1.6620 eV) energy gap ( E) was calculated to be 2.5687 eV. The spatial arrangements of the HOMO–LUMO, associated energies, and gap are represented in Figure 19 [52,53,54]. Some molecular parameters in the gas phase, based on HOMO–LUMO and optimized geometry, are presented in Table 8.

4. Conclusions

The gabapentin drug coordinated in a 1:1 ligand to metal ratio as a monobasic bidentate (ON) donor in all of the complexes. The analytical, magnetic, molar conductance, infrared vibrational motions, and electronic spectral study suggest that the structures shown in Figure 2 have an octahedral arrangement. The TEM morphology showed that the complexes were in nanosize ranges. The molecular docking results showed that the [Cu(II)–(Gpn)] metal complex interacted with both receptors more efficiently than the gabapentin drug and among all of them, [Cu(II)–(Gpn)]–serotonin (CuGS) had the highest binding energy value. A molecular dynamic (MD) simulation with a 100 ns run revealed that the CuGS complex possesses a more stable complex with the serotonin receptor than GpnS. The theoretical data obtained by DFT calculations agreed with the experimental data.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27134311/s1, Figure S1: Receptor-ligand complex (a) CuGS and (b) GpnS in triclinic box solvated with water molecules and neutralized with 28 Na+ and 27 Cl- ions (0.15 M salt). Figure S2: Representation of (a) aromatic surface, (b) hydrophobic surface, (c) solvent accessible surface, (d) ionizability surface, and (e) Interpolated charge; between serotonin and metal complex [Cu(II)-(Gpn)].

Author Contributions

Conceptualization, formal analysis, software, and data curation. M.S.R., S.S. and A.G.; writing—original draft preparation, M.S.R. and S.S.; Validation, resources and investigation, Y.S.A., H.A. and W.F.A.; visualization, A.M.A.A., M.I.K. and K.A.A.; Funding acquisition, M.S.R. and W.F.A.; Writing—review and editing, W.F.A. and A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia, project number 1-442-58.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All of the data supporting the reported results are available in the manuscript.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia for funding this work through project number 1-442-58.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Loscher, W.C. Current status and future in the pharmacology of epilepsy. Trends Pharmacol. Sci. 2002, 23, 113–118. [Google Scholar] [CrossRef]
  2. Ananda, K.; Aravinda, S.; Vasudev, P.G.; Raja, K.M.P.; Sivaramakrishnan, H.; Nagarajan, K.; Shamala, N.; Balaram, P. Stereochemistry of gabapentin and several derivatives: Solid state conformations and solution equilibria. Curr. Sci. 2003, 85, 1002–1011. [Google Scholar]
  3. McCabe, P.H. Role of levetiracetam in the treatment of epilepsy. Expert Opin. Pharmacother. 2000, 1, 633–674. [Google Scholar] [CrossRef]
  4. Elmastas, M.; Gulcin, I.; Beydemir, S.; Kufreviogluo Aboul-Enein, H.Y. A study on the in vitro antioxidant activity of juniper seeds extracts. Anal. Lett. 2006, 39, 47–65. [Google Scholar] [CrossRef]
  5. Regesta, G.; Tanganelli, P. Clinical aspects and biological bases of drug-resistant epilepsies. Epilepsy Res. 1999, 34, 109–122. [Google Scholar] [CrossRef]
  6. Moretti, R.; Antonello, R.M.; Torre, P.; Cazzato, G. Headache and neck pain: Gabapentin as a possible treatment. J. Headache Pain 2000, 1, 155–161. [Google Scholar] [CrossRef] [Green Version]
  7. Gielen, M.; Tiekink, E.R.T. (Eds.) Metallotherapeutic Drugs and Metal-Based Diagnostic Agents, the Use of Metals in Medicine; Wiley: Chichester, UK, 2005. [Google Scholar]
  8. Weder, J.E.; Dillon, C.T.; Hambley, T.W.; Kennedy, B.J.; Lay, P.A.; Biffin, J.R.; Regtop, H.L.; Daview, N.M. Copper complexes of non-steroidal anti-inflammatory drugs: An opportunity yet to be realized. Coord. Chem. Rev. 2002, 232, 95–126. [Google Scholar] [CrossRef]
  9. Nakai, M.; Sekiguchi, F.; Obata, M.; Ohtsuki, C.; Adachi, Y.; Sakurai, H.; Orvig, C.; Rehder, D.; Yano, S.J. Synthesis and insulin-mimetic activities of metal complexes with 3-hydroxypyridine-2-carboxylic acid. Inorg. Biochem. 2005, 99, 1275–1782. [Google Scholar] [CrossRef]
  10. Sadler, P.J.; Guo, Z. Metal complexes in medicine: Design and mechanism of action. Pure Appl. Chem. 1998, 70, 863–871. [Google Scholar] [CrossRef] [Green Version]
  11. Nadgi, A.; Mhlong, N.; Esoliman, M. Metal complexes in cancer therapy update from drug design perspective. Drug Des. Dev. Ther. 2017, 11, 599–618. [Google Scholar]
  12. Hu, P.; Snyder, S.A. Enantiospecific Total Synthesis of the Highly Strained (−)-Presilphiperfolan-8-ol via a Pd-Catalyzed Tandem Cyclization. J. Am. Chem. Soc. 2017, 14, 5007–5017. [Google Scholar] [CrossRef] [PubMed]
  13. Al-Mudhaffar, D.M.H.; Al-Mowali, A.H.; Al-Badran, U.H. Synthesis, Characterization, Biological Activities and Antioxidant Study of Gabapentin Derivatives and Their Complexes with [Cu (II), Zn (II]. Int. J. Appl. 2018, 8, 105–116. [Google Scholar] [CrossRef]
  14. Abd-Alrassol, K.S.; Qasim, Q.A.; AL-Salman HN, K.; Hussein, H.H. The Spectrophotometric Determination of Antiepileptic Drug in Standard and Pharmaceutical Formulations by Diazotization Coupling Reaction and Some Metals Complexes. Syst. Rev. Pharm. 2020, 11, 247–260. [Google Scholar]
  15. Eberhardt, J.; Santos-Martins, D.; Tillack, A.F.; Forli, S. AutoDock Vina 1.2.0: New Docking Methods, Expanded Force Field, and Python Bindings. J. Chem. Inf. Model. 2021, 61, 3891–3898. [Google Scholar] [CrossRef] [PubMed]
  16. Chu, C.-H.; Li, K.-M.; Lin, S.-W.; Chang, M.D.-T.; Jiang, T.-Y.; Sun, Y.-J. Crystal structures of starch binding domain from Rhizopus oryzae glucoamylase in complex with isomaltooligosaccharide: Insights into polysaccharide binding mechanism of CBM21 family. Proteins Struct. Funct. Bioinform. 2014, 82, 1079–1085. [Google Scholar] [CrossRef]
  17. Rehman, M.T.; AlAjmi, M.F.; Hussain, A. Natural compounds as inhibitors of SARS-CoV-2 main protease (3CLpro): A molecular docking and simulation approach to combat COVID-19. Curr. Pharm. Des. 2021, 27, 3577–3589. [Google Scholar] [CrossRef]
  18. O’Boyle, N.M.; Banck, M.; James, C.A.; Morley, C.; Vandermeersch, T.; Hutchison, G.R. Open Babel: An open chemical toolbox. J. Cheminform. 2011, 3, 33. [Google Scholar] [CrossRef] [Green Version]
  19. Dallakyan, S. PyRx-Python Prescription v. 0.8; Scripps Research Institute: La Jolla, CA, USA, 2008; Volume 2010. [Google Scholar]
  20. Bekker, H.; Berendsen, H.J.C.; Dijkstra, E.J.; Achterop, S.; Vondrumen, R.; Vanderspoel, D.; Sijbers, A.; Keegstra, H.; Renardus, M. Gromacs-a parallel computer for molecular-dynamics simulations. In Proceedings of the 4th International Conference on Computational Physics (PC 92); World Scientific Publishing: Singapore, 1993; pp. 252–256. [Google Scholar]
  21. Vanommeslaeghe, K.; Hatcher, E.; Acharya, C.; Kundu, S.; Zhong, S.; Shim, J.; Darian, E.; Guvench, O.; Lopes, P.; Vorobyov, I.; et al. CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additivebiological force fields. J. Comput. Chem. 2010, 31, 671–690. [Google Scholar] [CrossRef] [Green Version]
  22. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1, 19–25. [Google Scholar] [CrossRef] [Green Version]
  23. Allen, M.P.; Tildesley, D.J. Computer Simulations of Liquids; Clarendon Press: Oxford, UK, 1987. [Google Scholar]
  24. Oostenbrink, C.; Villa, A.; Mark, A.E.; Van Gunsteren, W.F. A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force-field parameter sets 53A5 and 53A6. J. Comput. Chem. 2004, 25, 1656–1676. [Google Scholar] [CrossRef]
  25. Steinbach, P.J.; Brooks, B.R. New Spherical-Cutoff Methods for Long-Range Forces in Macromolecular Simulation. J. Comput. Chem. 1994, 15, 667–683. [Google Scholar] [CrossRef]
  26. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 28–33. [Google Scholar] [CrossRef]
  27. DeLano, W.L. The PyMOL Molecular Graphics System. PyMOL. DeLano Scientific. 2002. Available online: http://www.pymol.org (accessed on 15 March 2022).
  28. Frisch, M.J.; Clemente, F.R. Gaussian 09, Revision a. 01, Mjfrisch, gw Trucks, Hbschlegel, Gescuseria, ma Robb, Jrcheeseman; Scalmani, G., Barone, V., Mennucci, B., Petersson, G.A., Nakatsuji, H., Caricato, M., Li, X., Hratchian, H.P., Izmaylov, A.F., Bloino, J., et al., Eds.; Gaussian Inc.: Wallingford, CT, USA, 2009; pp. 20–44. [Google Scholar]
  29. Zhurko, G.A.; Zhurko, D.A. Chemcraft—Graphical Program for Visualization of Quantum Chemistry Computations; Academic Version 1.5: Ivanovo, Russia, 2004. [Google Scholar]
  30. Geary, W.J. The use of conductivitymeasurements in organic solvents for the characterisation of coordination compounds. Coord. Chem. Rev. 1971, 7, 81–122. [Google Scholar] [CrossRef]
  31. Hariharan, P.C.; Pople, J.A. The effect of d-functions on molecular orbital energies for hydrocarbons. Chem. Phys. Lett. 1972, 16, 217–219. [Google Scholar] [CrossRef]
  32. Chimatadara, S.A.; Basavaraja, T.; Thabaja, K.A.; Sharanappa, T. Ruthenium (III) catalysed oxidation of gabapentin (neurontin) by diperiodatocuprate (III) in aqueous alkaline medium: A kinetic and mechanistic study. J. Mol. Catal. A Chem. 2007, 267, 65–71. [Google Scholar] [CrossRef]
  33. Hsu, C.H.; Ke, W.T.; Lin, S.Y. Progressive steps of polymorphic transformation of gabapentin polymorphs studied by hot-stage FTIR microspectroscopy. J. Pharm. Pharmaceut. Sci. 2010, 13, 67–77. [Google Scholar] [CrossRef] [Green Version]
  34. Kakanejadifard, A.; Esna-Ashari, F.; Hashemi, P.; Zabardasti, A. Spectrochim. Acta A 2013, 106, 80–85. [Google Scholar] [CrossRef]
  35. Nakamoto, K. Infrared Spectra of Inorganic and Coordination Compounds, 6th ed.; Wiley: Hoboken, NJ, USA, 2009; Volume 978, p. 232. [Google Scholar]
  36. Smith, H.E.; Neergaard, J.R.; Burrows, E.P.; Chen, F.M. Optically active amines. XVI. Exciton chirality method applied to the salicylidenimino chromophore. Salicylidenimino chirality rule. J. Am. Chem. Soc. 1974, 96, 2908–2916. [Google Scholar] [CrossRef]
  37. Lever, A.B.P. Inorganic Electronic Spectroscopy; Elsevier: Amsterdam, The Netherlands, 1968; Volume 249. [Google Scholar]
  38. Rana, V.B.; Singh, D.B.; Singh, P.; Teotia, M.P. Divalent nickel, cobalt and copper complexes of a tetradentate N6 macrocyclic ligand. Transit. Met. Chem. 1981, 6, 36–39. [Google Scholar] [CrossRef]
  39. Procter, I.M.; Hathaway, B.J.; Nicholls, P.J. The electronic properties and stereochemistry of the copper(II) ion. Part I. Bis(ethylenediamine)copper(II) complexes. Chem. Soc. A 1968, 1678–1684. [Google Scholar] [CrossRef]
  40. Akram, M.; Lal, H.; Shakya, S.; Kabir-ud-Di. Multispectroscopic and Computational Analysis Insight into the Interaction of Cationic Diester-Bonded Gemini Surfactants with Serine Protease α-Chymotrypsin. ACS Omega 2020, 5, 3624–3637. [Google Scholar] [CrossRef] [PubMed]
  41. Shakya, B.; Shakya, S.; Siddique, Y.H. Effect of geraniol against arecoline induced toxicity in the third instar larvae of transgenic Drosophila melanogaster (hsp70-lacZ) Bg9. Toxicol. Mech. Methods 2019, 29, 187–202. [Google Scholar] [CrossRef] [PubMed]
  42. Khan, I.M.; Shakya, S.; Islam, M.; Khan, S.; Najnin, H. Synthesis and spectrophotometric studies of CT complex between 1, 2-dimethylimidazole and picric acid in different polar solvents: Exploring antimicrobial activities and molecular (DNA) docking. Phys. Chem. Liq. 2021, 59, 753–769. [Google Scholar] [CrossRef]
  43. Krivák, R.; Jendele, L.; Hoksza, D. Peptide-Binding Site Prediction From Protein Structure via points on the Solvent Accessible Surface. In Proceedings of the 2018 ACM International Conference on Bioinformatics, Computational Biology, and Health Informatics, Washington, DC, USA, 29 August–1 September 2018. [Google Scholar]
  44. Ranjbar, A.; Jamshidi, M.; Torabi, S. Molecular modelling of the antiviral action of Resveratrol derivatives against the activity of two novel SARS-CoV-2 and 2019-nCoV receptors. Eur. Rev. Med. Pharmacol. Sci. 2020, 24, 7834–7844. [Google Scholar] [PubMed]
  45. Kufareva, I.; Abagyan, R. Methods of protein structure comparison. Methods Mol. Biol. 2012, 857, 231–257. [Google Scholar] [PubMed] [Green Version]
  46. Khan, M.D.; Shakya, S.; Vu HH, T.; Habte, L.; Ahn, J.W. Low concentrated phosphorus sorption in aqueous medium on aragonite synthesized by carbonation of seashells: Optimization, kinetics, and mechanism study. J. Environ. Manag. 2021, 280, 111652. [Google Scholar] [CrossRef]
  47. Dahham, S.S.; Ali, M.N.; Tabassum, H.; Khan, M. Studies on antibacterial and antifungal activity of pomegranate (Punica granatum L.). Am. Eurasian J. Agric. Environ. Sci. 2010, 9, 273–281. [Google Scholar]
  48. Sohn, H.-Y.; Kwon, C.-S.; Son, K.-H. Fungicidal effect of prenylated flavonol, papyriflavonol A, isolated from Broussonetiapapyrifera (L.) vent. against Candida albicans. J. Microbiol. Biotechnol. 2010, 20, 1397–1402. [Google Scholar] [CrossRef] [Green Version]
  49. Edziri, H.; Mastouri, M.; Mahjoub, M.A.; Mighri, Z.; Mahjoub, A.; Verschaeve, L. Antibacterial, Antifungal and Cytotoxic Activities of Two Flavonoids from Retama raetam Flowers. Molecules 2012, 17, 7284–7293. [Google Scholar] [CrossRef] [Green Version]
  50. Shakya, S.; Khan, I.M.; Ahmad, M. Charge transfer complex based real-time colorimetric chemosensor for rapid recognition of dinitrobenzene and discriminative detection of Fe2+ ions in aqueous media and human hemoglobin. J. Photochem. Photobiol. A Chem. 2020, 392, 112402. [Google Scholar] [CrossRef]
  51. Foresman, J.B. Exploring Chemistry with Electronic Structure Methods: A Guide to Using Gaussian; Frisch, E., Ed.; Gaussian, Incorporated: Pittsburg, PA, USA, 1996. [Google Scholar]
  52. Islam, M.R.; Shakya, S.; Selim, A.; Alam, M.S.; Ali, M. Solvatochromic Absorbance and Fluorescence Probe Behavior within Ionic Liquid+ γ-Butyrolactone Mixture. J. Chem. Eng. Data 2019, 64, 4169–4180. [Google Scholar] [CrossRef]
  53. Murugavel, S.; Ravikumar, C.; Jaabil, G.; Alagusundaram, P. Synthesis, crystal structure analysis, spectral investigations (NMR, FT-IR, UV), DFT calculations, ADMET studies, molecular docking and anticancer activity of 2-(1-benzyl-5-methyl-1H-1, 2, 3-triazol-4-yl)-4-(2-chlorophenyl)-6-methoxypyridine–a novel potent human topoisomerase IIα inhibitor. J. Mol. Struct. 2019, 1176, 729–742. [Google Scholar]
  54. Khan, I.M.; Shakya, S. Exploring Colorimetric Real-Time Sensing Behavior of a Newly Designed CT Complex toward Nitrobenzene and Co2+: Spectrophotometric, DFT/TD-DFT, and Mechanistic Insights. ACS Omega 2019, 4, 9983–9995. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The zwitterion structure of the gabapentin (Gpn) drug.
Figure 1. The zwitterion structure of the gabapentin (Gpn) drug.
Molecules 27 04311 g001
Figure 2. The structures of the synthesized Mn(II), Co(II), Ni(II), and Cu(II)–Gpn complexes, where = 2 for Cu(II), 3 for Ni(II), 4 for Mn(II), and 6 for Co(II).
Figure 2. The structures of the synthesized Mn(II), Co(II), Ni(II), and Cu(II)–Gpn complexes, where = 2 for Cu(II), 3 for Ni(II), 4 for Mn(II), and 6 for Co(II).
Molecules 27 04311 g002
Figure 3. The infrared spectra of (A) Gpn free drug and (B) Mn(II), Co(II), Ni(II), and Cu(II)–Gpn complexes.
Figure 3. The infrared spectra of (A) Gpn free drug and (B) Mn(II), Co(II), Ni(II), and Cu(II)–Gpn complexes.
Molecules 27 04311 g003
Figure 4. The SEM images of (I) Mn(II), (II) Cu(II), (III) Co(II), and (IV) Ni(II)–Gpn.
Figure 4. The SEM images of (I) Mn(II), (II) Cu(II), (III) Co(II), and (IV) Ni(II)–Gpn.
Molecules 27 04311 g004
Figure 5. The TEM images of the (I) Mn(II), (II) Cu(II), (III) Co(II), and (IV) Ni(II)–Gpn complexes.
Figure 5. The TEM images of the (I) Mn(II), (II) Cu(II), (III) Co(II), and (IV) Ni(II)–Gpn complexes.
Molecules 27 04311 g005
Figure 6. The best docked pose showing a helical model of serotonin (PDB ID: 6BQH) docked with (a) the Gpn drug and (b) metal complex [Cu(II)–(Gpn)].
Figure 6. The best docked pose showing a helical model of serotonin (PDB ID: 6BQH) docked with (a) the Gpn drug and (b) metal complex [Cu(II)–(Gpn)].
Molecules 27 04311 g006
Figure 7. The 3D representation of the interactions for serotonin (PDB ID: 6BQH) docked with (a) the Gpn drug and (b) metal complex [Cu(II)–(Gpn)].
Figure 7. The 3D representation of the interactions for serotonin (PDB ID: 6BQH) docked with (a) the Gpn drug and (b) metal complex [Cu(II)–(Gpn)].
Molecules 27 04311 g007
Figure 8. The 2D representation of interactions for serotonin (PDB ID: 6BQH) docked with (a) the Gpn drug and (b) metal complex [Cu(II)–(Gpn)].
Figure 8. The 2D representation of interactions for serotonin (PDB ID: 6BQH) docked with (a) the Gpn drug and (b) metal complex [Cu(II)–(Gpn)].
Molecules 27 04311 g008
Figure 9. Representation of the hydrogen binding surface between the serotonin and metal complex [Cu(II)–(Gpn)].
Figure 9. Representation of the hydrogen binding surface between the serotonin and metal complex [Cu(II)–(Gpn)].
Molecules 27 04311 g009
Figure 10. The root mean square deviation (RMSD) of the solvated receptor backbone and ligand complex during the 100 ns MD simulation [unbound serotonin receptor (yellow), GpnS complex (black), and CuGS complex (brown)]. According to the literature, a RMSD value under <3.0 Å is the most acceptable [44]. The drop in the RMSD value for CuGS showed an alteration in the secondary structure conformation of the protein due to the [Cu(II)–(Gpn)] interaction. This finding shows that CuGS developed a more stable combination. The results confirm that this interaction brings protein chains closer and reduces the gap between them, as shown in Figure 11 [45].
Figure 10. The root mean square deviation (RMSD) of the solvated receptor backbone and ligand complex during the 100 ns MD simulation [unbound serotonin receptor (yellow), GpnS complex (black), and CuGS complex (brown)]. According to the literature, a RMSD value under <3.0 Å is the most acceptable [44]. The drop in the RMSD value for CuGS showed an alteration in the secondary structure conformation of the protein due to the [Cu(II)–(Gpn)] interaction. This finding shows that CuGS developed a more stable combination. The results confirm that this interaction brings protein chains closer and reduces the gap between them, as shown in Figure 11 [45].
Molecules 27 04311 g010
Figure 11. The superimposed structure of (a) the unbound serotonin receptor (yellow) and serotonin receptor after simulation (brown) for GpnS and (b) the unbound serotonin receptor (yellow) and serotonin receptor after simulation (green) for CuGS.
Figure 11. The superimposed structure of (a) the unbound serotonin receptor (yellow) and serotonin receptor after simulation (brown) for GpnS and (b) the unbound serotonin receptor (yellow) and serotonin receptor after simulation (green) for CuGS.
Molecules 27 04311 g011
Figure 12. The RR distance map displaying patterns of spatial interactions between (a) the unbound serotonin receptor and the serotonin receptor after the simulation for GpnS; and (b) the unbound serotonin receptor and serotonin receptor after simulation for CuGS, showing the average distance and standard deviation for all amino acid pairs. (MDS = Molecular Dynamics simulation).
Figure 12. The RR distance map displaying patterns of spatial interactions between (a) the unbound serotonin receptor and the serotonin receptor after the simulation for GpnS; and (b) the unbound serotonin receptor and serotonin receptor after simulation for CuGS, showing the average distance and standard deviation for all amino acid pairs. (MDS = Molecular Dynamics simulation).
Molecules 27 04311 g012
Figure 13. The radius of gyration (Rg) for the unbound serotonin receptor (yellow), GpnS complex (black), and CuGS complex (brown) during the 100 ns simulation time.
Figure 13. The radius of gyration (Rg) for the unbound serotonin receptor (yellow), GpnS complex (black), and CuGS complex (brown) during the 100 ns simulation time.
Molecules 27 04311 g013
Figure 14. The number of average hydrogen bonding interactions between (a) the CuGS complex and (b) GpnS complex during the 100 ns simulation time.
Figure 14. The number of average hydrogen bonding interactions between (a) the CuGS complex and (b) GpnS complex during the 100 ns simulation time.
Molecules 27 04311 g014
Figure 15. The solvent accessible surface area analysis for the unbound serotonin receptor (yellow), GpnS complex (black), and CuGS complex (brown) during the 100 ns simulation time.
Figure 15. The solvent accessible surface area analysis for the unbound serotonin receptor (yellow), GpnS complex (black), and CuGS complex (brown) during the 100 ns simulation time.
Molecules 27 04311 g015
Figure 16. The optimized structure of the synthesized metal complexes (a) [Mn(II)–(Gpn)], (b) [Co(II)–(Gpn)], (c) [Ni(II)–(Gpn)], and (d) [Cu(II)–(Gpn)] with the Mulliken atom numbering scheme.
Figure 16. The optimized structure of the synthesized metal complexes (a) [Mn(II)–(Gpn)], (b) [Co(II)–(Gpn)], (c) [Ni(II)–(Gpn)], and (d) [Cu(II)–(Gpn)] with the Mulliken atom numbering scheme.
Molecules 27 04311 g016
Figure 17. The MEP surface map of the synthesized metal complexes (a) [Mn(II)–(Gpn)], (b) [Co(II)–(Gpn)], (c) [Ni(II)–(Gpn)], and (d) [Cu(II)–(Gpn)] with the respective color scales.
Figure 17. The MEP surface map of the synthesized metal complexes (a) [Mn(II)–(Gpn)], (b) [Co(II)–(Gpn)], (c) [Ni(II)–(Gpn)], and (d) [Cu(II)–(Gpn)] with the respective color scales.
Molecules 27 04311 g017
Figure 18. The optimized structure of the synthesized metal complex [Cu(II)–(Gpn)] showing bond lengths.
Figure 18. The optimized structure of the synthesized metal complex [Cu(II)–(Gpn)] showing bond lengths.
Molecules 27 04311 g018
Figure 19. The spatial plot of HOMO and LUMO with their energy gap for the synthesized metal complex [Cu(II)–(Gpn)].
Figure 19. The spatial plot of HOMO and LUMO with their energy gap for the synthesized metal complex [Cu(II)–(Gpn)].
Molecules 27 04311 g019
Table 1. Infrared spectral data (cm−1) of Gpn and its complexes.
Table 1. Infrared spectral data (cm−1) of Gpn and its complexes.
CompoundsFrequencies, cm−1
νas(COO)νs(COO)δ(NH2)ν(M-O)ν(M-N)
Mn(II)165014101568617502
Cu(II)166013901558612516
Co(II)167213981552620519
Ni(II)167813951553615468
Table 2. The docking score of the four metal complexes docked with two receptors [serotonin (6BQH) and dopamine (6CM4)].
Table 2. The docking score of the four metal complexes docked with two receptors [serotonin (6BQH) and dopamine (6CM4)].
S. No.Receptor ComplexBinding Free Energy (kcal/mol)
PDB: 6BQHPDB: 6CM4
1[Mn(II)–(Gpn)]−6.9−6.6
2[Co(II)–(Gpn)]−7.0−6.7
3[Ni(II)–(Gpn)]−6.9−6.4
4[Cu(II)–(Gpn)]−7.2−6.5
5Gpn−5.1−4.8
Table 3. The interactions of the [Cu(II)–(Gpn)] complex and Gpn with serotonin (6BQH).
Table 3. The interactions of the [Cu(II)–(Gpn)] complex and Gpn with serotonin (6BQH).
S. No.ReceptorBinding Free Energy (kcal/mol)Interactions
H-BondOthers
1[Cu(II)–(Gpn)]−7.2Arg173, Thr109, His182, and Asn187Asp172 (Attractive charge) and Ala176, Arg173 (Alkyl)
2Gpn−5.1Thr109 and Ala108Ala321 (Alkyl)
Table 4. The [Cu(II)–(Gpn)]–serotonin interaction results by Discovery Studio.
Table 4. The [Cu(II)–(Gpn)]–serotonin interaction results by Discovery Studio.
NameDistanceCategory
THR109:HG1–[Cu(II)–Gpn]:O2.50Hydrogen Bond
ARG173:HH11–[Cu(II)–Gpn]:O2.63Hydrogen Bond
ARG173:HH12–[Cu(II)–Gpn]:O2.78Hydrogen Bond
ASN187:HD22–[Cu(II)–Gpn]:O2.43Hydrogen Bond
[Cu(II)-Gpn]:H26–HIS182:O2.59Hydrogen Bond
[Cu(II)-Gpn]:H26–HIS182:O2.49Hydrogen Bond
[Cu(II)-Gpn]:H25–ASN187:OD12.35Hydrogen Bond
[Cu(II)-Gpn]:H25–ASN187:OD12.50Hydrogen Bond
ALA321–[Cu(II)–Gpn]5.31Hydrophobic
Table 5. The Gpn–serotonin interaction results by Discovery Studio.
Table 5. The Gpn–serotonin interaction results by Discovery Studio.
NameDistanceCategory
Gpn:N–ASP172:OD25.34Electrostatic
ALA108:HN–Gpn:O2.30Hydrogen Bond
THR109:HN–Gpn:O1.85Hydrogen Bond
THR109:HG1–Gpn:O1.97Hydrogen Bond
Cu(II)-Gpn:H–Gpn:O2.91Hydrogen Bond
ARG173–Gpn4.93Hydrophobic
ALA176–Gpn5.04Hydrophobic
Table 6. The bond lengths of the metal complex [Cu(II)–(Gpn)] obtained through DFT.
Table 6. The bond lengths of the metal complex [Cu(II)–(Gpn)] obtained through DFT.
S. No.[Cu(II)–(Gpn)] (B3LYP/LanL2DZ)
Atom No.Bond Length (Å)Atom No.Bond Length (Å)
1R(1–28)1.258R(16–17)1.104
2R(2–28)1.400R(16–18)1.103
3R(2–39)1.944R(16–19)1.56
4R(3–4)1.053R(19–20)1.105
5R(3–5)1.064R(19–21)1.103
6R(3–22)1.529R(22–23)1.109
7R(3–39)2.104R(22–24)1.104
8R(6–7)1.578R(25–26)1.103
9R(6–19)1.577R(25–27)1.109
10R(6–22)1.574R(25–28)1.574
11R(6–25)1.578R(30–31)1.024
12R(7–8)1.103R(30–32)1.01
13R(7–9)1.104R(30–39)2.222
14R(7–10)1.560R(33–34)1.008
15R(10–11)1.104R(33–35)1.009
16R(10–12)1.103R(33–39)2.088
17R(10–13)1.559R(36–37)1.025
18R(13–14)1.104R(36–38)1.008
19R(13–15)1.104R(36–39)2.102
20R(13–16)1.559R(29–39)3.798
Table 7. The Mulliken atomic charges of the metal complex [Cu(II)–(Gpn)] atoms.
Table 7. The Mulliken atomic charges of the metal complex [Cu(II)–(Gpn)] atoms.
S. No.Synthesized [Cu(II)-(Gpn)] Complex
Mulliken Atomic NumbersMulliken Atomic ChargesMullikenAtomic NumbersMulliken Atomic Charges
11O−0.1939121H0.06922
22O−0.1170322C−0.06068
33N−0.3464623H0.08103
44H0.2060624H0.09072
55H0.2381725C−0.16982
66C0.0292626H0.08199
77C−0.1347127H0.11175
88H0.0709728C0.21118
99H0.0721229Cl−0.7276
1010C−0.1335130O−0.28504
1111H0.0716631H0.24647
1212H0.0644632H0.20389
1313C−0.1324733O−0.25834
1414H0.0671734H0.23457
1515H0.0699235H0.23367
1616C−0.132636O−0.28583
1717H0.0764337H0.25414
1818H0.0657438H0.21341
1919C−0.13937C39u−0.03369
2020H0.08705
Table 8. The various other theoretical molecular parameters of the metal complex [Cu(II)–(Gpn)].
Table 8. The various other theoretical molecular parameters of the metal complex [Cu(II)–(Gpn)].
ParametersRB3LYP/LanL2DZ
Minimum SCF energy (a.u.)−982.219370
Polarizability (a) (a.u.)84.074494
Dipole moment (Debye)7.662013
Zero point vibrational energy (kcal/mol)241.62872
Total thermal energy (kcal/mol)253.148
Electronic spatial extent (a.u.)5378.5506
Frontier MO energies (eV)
LUMO−1.6620
HOMO−4.2308
Gap (HOMO–LUMO)2.5687
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Refat, M.S.; Gaber, A.; Althobaiti, Y.S.; Alyami, H.; Alsanie, W.F.; Shakya, S.; Adam, A.M.A.; Kobeasy, M.I.; Asla, K.A. Spectroscopic and Molecular Docking Studies of Cu(II), Ni(II), Co(II), and Mn(II) Complexes with Anticonvulsant Therapeutic Agent Gabapentin. Molecules 2022, 27, 4311. https://doi.org/10.3390/molecules27134311

AMA Style

Refat MS, Gaber A, Althobaiti YS, Alyami H, Alsanie WF, Shakya S, Adam AMA, Kobeasy MI, Asla KA. Spectroscopic and Molecular Docking Studies of Cu(II), Ni(II), Co(II), and Mn(II) Complexes with Anticonvulsant Therapeutic Agent Gabapentin. Molecules. 2022; 27(13):4311. https://doi.org/10.3390/molecules27134311

Chicago/Turabian Style

Refat, Moamen S., Ahmed Gaber, Yusuf S. Althobaiti, Hussain Alyami, Walaa F. Alsanie, Sonam Shakya, Abdel Majid A. Adam, Mohamed I. Kobeasy, and Kareem A. Asla. 2022. "Spectroscopic and Molecular Docking Studies of Cu(II), Ni(II), Co(II), and Mn(II) Complexes with Anticonvulsant Therapeutic Agent Gabapentin" Molecules 27, no. 13: 4311. https://doi.org/10.3390/molecules27134311

APA Style

Refat, M. S., Gaber, A., Althobaiti, Y. S., Alyami, H., Alsanie, W. F., Shakya, S., Adam, A. M. A., Kobeasy, M. I., & Asla, K. A. (2022). Spectroscopic and Molecular Docking Studies of Cu(II), Ni(II), Co(II), and Mn(II) Complexes with Anticonvulsant Therapeutic Agent Gabapentin. Molecules, 27(13), 4311. https://doi.org/10.3390/molecules27134311

Article Metrics

Back to TopTop