Introduction

The use of nanoparticles (NPs) in biological and industrial applications such as bio-imaging, medication administration, plasmonics, photonics, chemical sensors, anti-microbial activities, cell electrodes, optical devices, anti-microbial coating and hyperthermia therapy is being investigated progressively1,2,3,4,5. Thus, the experimental and theoretical studies concerning the interaction of nanoparticles with biological systems including nucleic acids, neurotransmitters, biomolecules, peptides, and proteins are very crucial in order to have a better understanding of the adsorption mechanism of aforementioned biological systems over the NPs surface. This can help in safe and proficient biological and industrial application of nanoparticles based bio-molecular conjugates. Among various NPs, the silver nano-particles (AgNPs) have drawn great attention from the scientific community and gaining popularity as suitable nanomaterials for biological applications owing to their remarkable properties including good biocompatibility with biological molecules (i.e. amino acid, peptides and nucleic acids), non-toxicity as well as shape-dependent electronic and optical features6,7,8,9,10. Furthermore, the AgNPs have unique optical, electrical, and thermal properties and a high surface-to-volume ratio that make them ideal for cutting-edge applications in the fields of medicine, biology, mechanics, packaging, food science, electronics, and information technology, among others11.

Tyramine is an essential trace biogenic amine and acts as a neurotransmitter12. Trace biogenic amines play an important role in many human disorders that involve mood, emotion, cognition, and attention. Tyramine is being synthesized by the mammalian brain and peripheral nerve tissues13,14. It is found in excessive extent in assorted variety of raw and fermented foods together with fish, meat, fruits, cheese, soybean products, and wine15. A large amount of consumption of tyramine rich food shows toxic effect on human health by producing the symptoms like flush, rash, hypertonia, vomiting, palpitation and tachycardia. Besides, it also leads to dangerous blood pressure increase and chance of strong migraine headaches16. The remarkable set of properties of AgNPs and their application in biology, medicine, chemical sensors and bio-imaging and the scientific usefulness of tyramine motivated us to perform a detailed in silico investigation on the modeling of AgNPs@Tyramine bio-molecular conjugate to have a deep insight into the molecular interaction of neurotransmitter with AgNPs cluster and its electronic properties. The interaction of AgNPs cluster and tyramine can be a crucial perquisite in drug delivery applications of tyramine related diseases and many other future industrial and biological investigations.

Literature survey revealed that the interaction of silver nano-particles (AgNPs) with biological systems have been studied by a number of authors in many different aspects6,17,18,19,20,21,22,23,24. Ando and coworkers used silver, gold, and silver-gold alloy as imaging probe to trail motions of isolated biomolecules17. In order to track the motion of isolated biomolecules, they designed the multicolor high-speed single particle tracking system of silver, gold, and silver-gold alloy. The interaction of silver nano-particles with monosaccharide was performed with the help of density functional theory by Gallegos and coworkers in order to explore a better reducing agent among α-Dglucose, α-d-ribose, D-erythrose, and glyceraldehyde in green synthesis approach18. Besides, they also performed theoretical computations on electron localization function analysis which helped in the explaining the reduction procedure in the formation of AgNPs. A comparative investigation of the interaction mechanisms between various sized, shaped, and surface functionalized silver and gold nanomaterials to commercially available human transferrin, a glycosylated protein, and to its nonglycosylated recombinant form was reported by Barber et al.19. Fahy and his co-workers used silver-engineered nanomaterials in evaluation of surface charge role in protein interactions and cellular cytotoxicity20. The report on kinetic study of AgNPs under acidic conditions and impact of protein interactions i.e., bovine serum albumin (BSA) to the colloidal stability of silver nano-particles was made by Tai et al.21. A detailed study on the interaction of ionic liquids with AgNPs was reported by Banjare and his group22. It was reported that AgNPs could be effective drug carriers for anti-malarial drugs for treatment of malaria and Covid-19 by performing theoretical calculations at DFTB3LYP/6-311++ g(d,p) level6. Venkatesh and coworkers reported the geometrical, vibrational and physicochemical properties of 2-deoxy-D-glucose by performing the DFT calculations on its interaction with AgNPs23. Moreover, Jadoon with the help of theoretical calculations showed that silver (Ag6) decorated Coronene quantum dot is an effective non-enzymatic sensor for glucose and H2O2 detection and in another study they proposed silver cluster decorated graphene nanoflakes can be an effective detector for detection of nitroaniline isomers24,25. Adsorption of p-aminophenol over silver-graphene composite was explored by Ayub and co-workers26. DFT study on the sensitivity of silver-graphene quantum dots for vital and harmful analytes were also performed27. Moreover, silver cluster based sensors has been proposed to expose various other molecules as well28,29. We have reported possible low lying energy conformers, vibrational dynamics and NBO analysis of isolated tyramine and its interaction with HCl in our earlier reported investigation12. Now we are extending our work on tyramine to investigate its adsorption and charge transfer mechanism over the surface of Ag6 nanoparticle cluster.

Computational details

All the theoretical calculations have been performed at Gaussian09 suite30 by employing density functional theory (DFT) through B3LYP functional. The Lanl2DZ and 6-31G(d,p) basis sets are applied for Ag atoms and tyramine-molecule. The adsorption energies of the most favorable Ag6@tyramine complex have been computed as follows:

$$\user2{ E}_{{{\varvec{ads}}}} = {\varvec{E}}\left( {{\mathbf{Ag}}_{{\mathbf{6}}} - {\mathbf{Tyramine}}} \right) \, {-}{\varvec{E}}\left( {{\mathbf{Ag}}_{{\mathbf{6}}} } \right) \, {-}{\varvec{E}}\left( {{\mathbf{Tyramine}}} \right) \, + \left( {{\mathbf{BSSE}}} \right)$$
(1)

where, E(Ag6-Tyramine) is the energy of Ag6@tyramine complex. The E(Ag6) and E(Tyramine) terms denote the energies of relaxed Ag6 cluster and tyramine molecule, respectively. The BSSE which is basis set superposition error (BSSE) for the adsorption energy was computed by applying counterpoise method. Besides, the dispersion calculations which account interactions including short and long range contributions were also computed in case of the most favorable configuration of Ag6@tyramine bio-molecular conjugate through D3 version of Grimme's dispersion with the original D3 damping function. Accordingly, a negative adsorption energy implies to the thermodynamic stability of the obtained complex. Furthermore, natural bond orbital (NBO) analysis has been performed for the most favorable configuration of Ag6@tyramine in the gas phase at the same level using NBO3.1 program implemented in Gaussian09 suite. NBO method gives the information about the intra and intermolecular binding interaction, donor–acceptor interactions, delocalization of charge density and hyper-conjugative interaction energies of any molecular system. The second-order hyper-conjugative interaction energy (E(2)) has been computed using the following equation31,32,33,34:

$$\user2{ E}^{\left( 2 \right)} = {{\varvec{\Delta}}}{\varvec{E}}_{{{\varvec{ij}}}} = \user2{ q}_{{\varvec{i}}} \frac{{{\varvec{F}}_{{{\varvec{ij}}}}^{2} }}{{{{\varvec{\upvarepsilon}}}_{{\varvec{j}}} \user2{ } - {{\varvec{\upvarepsilon}}}_{{\varvec{i}}} }}$$
(2)

The qi is the ith donor orbital occupancy, εi and εj are the diagonal elements (orbital energies), and Fij is the off-diagonal NBO Fock Matrix element.

Results and discussion

Tyramine adsorption over Ag6 cluster

The optimization of Ag6@tyramine in various possible configurations is carried out to find out the most favorable adsorption configurations of Ag6@tyramine bio-molecular conjugate and these optimized configurations have been used for further analysis. Consequently, five possible configurations for Ag6@tyramine conjugates are obtained and depicted in Fig. 1. The adsorption energy, which is very crucial for any nanomaterials based bio-molecular conjugate to utilize it practically in various biological and industrial applications is also computed for each configuration and listed in Table 1. All the five possible configurations have been simulated in gas and aqueous phase as well. Moreover, Ag6@tyramine bio-molecular conjugate contains even number electrons, and thereby it offers singlet multiplicity and neutral charge. Figure 1 reveals that Ag6 cluster interacts with amine group of side chain (i.e. in case of configuration-A) then the bond length of N-Ag and N–H…Ag are 2.37 and 2.80 Å, respectively with the adsorption energy of − 14.36 kcal/mol (without BSSE term) and the BSSE is computed to be 3.12 kcal/mol. Thus including BSSE term the adsorption energy is − 11.24 kcal/mol for this configuration. Moreover, the adsorption energy for this configuration after applying D3 version of Grimme’s dispersion with the original D3 damping function is computed to be − 21.65 kcal/mol. In configuration B, the interaction of oxygen atom (O) of hydroxyl group with Ag6 cluster results − 6.01 kcal/mol adsorption energy, which demonstrates that tyramine molecule, interacts more effectively through amine group with Ag6 cluster rather than OH group. The interaction distance between oxygen atom (O) and Ag is computed to be 2.48 Å in configuration B. In case of configurations-C, D and E the adsorption energies are computed to be − 1.02, − 5.35, and − 3.21 kcal/mol, respectively. The prominent interaction distances between different atoms of tyramine and Ag atoms have been shown in Fig. 1. A literature survey35,36,37,38,39,40 reveals that the magnitude of adsorption energy neither should be very high nor very low as desorption of drugs from nanomaterials can occur easily which is very necessary in drug delivery and sensing applications. The magnitude of adsorption energy in case of the most favorable configuration of Ag6@tyramine bio-molecular conjugate is good to utilize it in various biological and industrial applications. Moreover, the charge transfer between Ag6 cluster and tyramine is also computed using NBO approach for the most favorable adsorption configuration (i.e. A). NBO calculations indicate that there is a charge transfer of 0.46e from tyramine to Ag6 cluster. Thus, the direction of charge transfer is from drug to cluster and the obtained complex is stabilized through charge transfer process.

Figure 1
figure 1

Various optimized configurations of AgNPs@tyramine.

Table 1 Value of adsorption energies for various configurations of AgNPs@tyramine.

Some important electronic parameters viz HOMO energy (EH), LUMO energy (EL), energy gap (Eg) and Fermi level (EF) for Ag6 cluster before and after tyramine adsorption are summarized in Table 2. The study of highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) are very crucial in quantum chemistry as these provides details about electronic characteristics of molecular species41. Furthermore, graphical representation of HOMO and LUMO iso-surfaces is depicted in Fig. 2.

Table 2 Calculated HOMO energy (EH), LUMO energy (EL), HOMO–LUMO energy gap (Eg), Fermi level energy (EF) and work function (ϕ). All data are in eV.
Figure 2
figure 2

HOMO and LUMO surfaces of AgNPs cluster, Tyramine and AgNPs@tyramine.

This figure infers that the HOMO is mainly located over the surface of Ag6 cluster and nitrogen atom (N) of tyramine rather than whole tyramine molecule. However, the LUMO is located only on the Ag6 cluster. Table 2 shows that the value of HOMO energy (EH), LUMO energy (EL), and energy gap (Eg) in the most favorable configuration is computed to be − 4.59, − 2.06 and 2.53 eV respectively. The corresponding values are − 5.23 and − 2.54 eV with an energy gap of 2.69 eV for isolated Ag6 cluster. The magnitude of energy gap in the most favorable adsorption configuration decrease by 0.16 eV (i.e. 6%.) after tyramine adsorption. Although, this decrement in the energy gap is not remarkable but it is still noteworthy. The decrement in the magnitude of energy gap will enhance the electrical conductivity in exponential manner as energy gap is directly proportional to the electrical conductivity by the following equation42:

$${\upsigma } \propto {\mathbf{exp}} \, \left( {\frac{{ - {\varvec{E}}_{{\varvec{g}}} }}{{{\varvec{2KT}}}}} \right)$$
(3)

This electrical conductivity can be converted into electrical signals to expose tyramine. The energy gap in rest configurations is tabulated in Table 2. Nhat et al.42 performed the interaction of thioguanine with gold cluster (Au6) and reported the change in energy gap by 5.9%. Moreover, the reduction in the energy gap in the present study in case of the most favorable Ag6@tyramine is mainly due to the interaction on nitrogen (N) atom of amine group present in the side chain of tyramine with silver atoms of Ag6 cluster. Total density of states (DOS) is studied in inspecting the electronic structure of tyramine on the surface silver cluster using GaussSum program43. Figure 3 illustrates a basic representation of the parameters of the molecular orbitals in a certain energy range and exhibits population analysis per orbital. The effect of tyramine interaction with Ag6 cluster on Fermi level and work function. The Fermi level is explored accordingly:

$$\user2{ E}_{{\varvec{F}}} = \user2{ E}_{{{\varvec{HOMO}}}} + \frac{{({\varvec{E}}_{{{\varvec{LUMO}}}} - {\varvec{E}}_{{{\varvec{HOMO}}}} )}}{2}\user2{ }$$
(4)
Figure 3
figure 3

The density of states (DOS) for (a) Tyramine, (b) AgNPs and (c) AgNPs@tyramine.

The calculations infer that energy of Fermi level in Ag6 cluster is − 3.88 eV which gets shift towards higher value in all the obtained Ag6@tyramine configurations except configuration-E when tyramine interacts Ag6 cluster. Besides, the work function (Φ = − \({E}_{F}\)) of Ag6 cluster before tyramine adsorption is computed to be 3.88 eV. However, the adsorption of tyramine lowers the magnitude of work function which shows that the field emission properties of Ag6 cluster may alter upon adsorption of tyramine. Moreover, the zero point vibration energy in case of the most favorable configuration is computed to be 115.28 kcal/mol.

Solvent effect

The effect of aqueous media on adsorption energy and electronic properties of the most favorable Ag6@tyraminebio-molecular conjugate is also explored using polarizable continuum model (PCM)44 in order to understand the impact of the real biological situation on adsorbate-adsorbent configuration. Water has been considered as solvent as it acts key role in human body. In case of aqueous phase, the adsorption energy of Ag6@tyramine conjugate in its most favorable configuration has been computed to be − 11.06 kcal/mol which shows strong effect on the adsorption energy of Ag6@tyramine conjugate in aqueous medium. Furthermore, the negative sign of adsorption energy indicates that the adsorption of tyramine over the surface of Ag6 cluster is of spontaneous nature. The values of dipole moment, the energy of HOMO–LUMO orbitals, energy gap and Fermi energy in the aqueous phase are tabulated in Table 3. This table reveals an increment in the energy gap between HOMO and LUMO by 0.18 eV in the presence of aqueous medium which shows that reactivity increase slightly in aqueous phase compared to the gas phase. The increased energy gap is attributed to the fact that interaction of water molecules with tyramine stabilizes HOMO and destabilizes LUMO. Dipole moment which is crucial parameter to explore the symmetricity of the molecular system and privileged magnitude of dipole moment indicates a strong interaction between the adsorbate and adsorbent. The magnitude of dipole moment in gas phase for the most favorable configuration is computed to be 5.51 Debye while it is 6.81 Debye in aqueous phase. This increment in dipole moment infers that the solubility of Ag6@tyramine conjugate gets enhanced by varying the medium from gas to aqueous. Additionally, the Gibbs (free) energy has been calculated using the formula below30:

$${\mathbf{DG}}^{{\mathbf{0}}} \left( {{\mathbf{298}}} \right) \, = \sum ({\varvec{\varepsilon}}_{0} + {\varvec{G}}_{{{\varvec{corr}}}} )_{{{\mathbf{final}}}} - \sum ({\varvec{\varepsilon}}_{0} + {\varvec{G}}_{{{\varvec{corr}}}} )_{{{\mathbf{Inial}}}}$$
(5)

Here, \({(\varepsilon }_{0}+{G}_{corr})\) is the total sum of electronic and thermal free energies. In gas phase Gibbs (free) energy is computed to be − 3.98 and − 0.87 kcal/mol in aqueous phase corresponding to the most favorable configuration. The negative sign indicates that adsorption of tyramine over Ag6 cluster is of exothermic nature.

Table 3 Calculated HOMO energy (EH), LUMO energy (EL), HOMO–LUMO energy gap (Eg), Fermi level energy (EF) and work function (ϕ) in aqueous phase.

Recovery time

The recovery time could provide a theoretical estimation for light controlled desorption mechanism. Based on the conventional transition state theory, the recovery time can be computed using the Arrhenius-type equation40:

$$\tau = \nu_{0}^{ - 1} \exp \left( {\frac{{{-}E_{bin} }}{KT}} \right)$$
(6)

The K, υ denotes the Boltzmann’s constant and attempt frequency, respectively. The Ebin could be assumed as the potential barrier for the desorption process. Herein, T has been inserted at room temperature (298.15 K). The recovery times for desorption of tyramine from the studied Ag6 cluster in the presence of different light (UV, Vis, IR) are tabulated in Table 4 for the most stable configuration with adsorption energy of − 14.36 kcal/mol. The obtained recovery times in microsecond in the exposure of different light, which can provide light release control for tyramine desorption. Moreover the biosensor can be recovered rapidly.

Table 4 Recovery time for tyramine from Ag6 cluster.

NBO analysis

The second order Fock matrix was performed to determine the donor–acceptor interactions in the NBO for which we used the formula of second order perturbation approach to deduce the hyper conjugative energy (E (2))31,32,33,34. In present study for the tyramine interaction with Ag6 particles, only the interaction energies of core and lone pair orbitals to anti-bonding orbitals are considered and collected in Table S1. The intramolecular interactions are formed through overlapping of bonding and anti bonding orbitals, which impel the intramolecular charge transfer (ICT) causing stabilization of the system. The required energy for the stabilization of the assorted allowed transitions within bonding and anti-bonding levels is called hyper conjugative energy. The most probable transitions have higher magnitude of hyper conjugative energy (E(2)) as compared to other transitions. The transitions LP1(O20)→σ*(C4-C5), LP1(O20)→π*(C4-H5), LP1(N17)→σ*(C14-H15), LP1(N17)→σ*(C14-H16) having interaction energies 4.60 kcal/mol, 25.17 kcal/mol, 5.32 kcal/mol, 1.63 kcal/mol, respectively, indicating the transition LP1(O20) →π*(C4-H5) is highly intense causing charge transfer among them. There is another intensive transition CR1(C4) → σ*(C4-O20) having energy 1.43 kcal/mol. The E(2) energy for the interaction between tyramine and Ag atoms have also been collected in Table S1. There are some transitions which are very intense like LP1(N17) → LP*6(Ag25), LP1(N17) → LP*8(Ag25), σ (C14-N17) → LP*6(Ag25) with transition energy 18.66 kcal/mol, 4.92 kcal/mol and 0.89 kcal/mol, respectively. There are some other intensive interactions with Ag22, Ag23 and Ag24 with tyramine.

The occupancies of different bonding, anti-bonding, core and lone pair orbitals of carbon, nitrogen and oxygen atoms are collected in supplementary materials (see Table S2). Results indicate maximum occupancy in bonding, core and lone pair orbitals of carbon, nitrogen and oxygen atoms whereas it is very less in anti-bonding orbitals. In present investigation, the occupancy of bonding orbital of π type i.e. π(C1–C2) shows that the charge has been delocalized from this bonding orbital to its anti-bonding orbital. Moreover, the lone pair orbitals of nitrogen and oxygen LP(1) (N17) and LP(1) (O23) are calculated with 1.89e and ~ 1.98e. The occupancies of all the Ag atoms are very high close to 1e which indicates that the electrons are highly localized in Ag atoms; but when it interacts with tyramine molecule except for LP5(Ag22) which occupancy is 1.97e, i.e., this electron is resonating more.

Conclusion

In conclusion, this research utilized the density functional theory (DFT) method to investigate the binding of tyramine to clusters of silver nanoparticle. Results could be summarized as follows:

  • Tyramine binds strongly to clusters of silver nanoparticles, with adsorption energies of − 14.36 kcal/mol in the gas phase and − 11.06 kcal/mol in the aqueous phase.

  • The amine group of the tyramine side chain interacts more effectively with the silver cluster compared to other sites.

  • The most favorable adsorption configuration results in a 6% decrease in the energy gap, leading to increased electrical conductivity.

  • The Gibbs (free) energy is negative in both the gas and aqueous phases, indicating an exothermic adsorption process.

  • NBO calculations reveal a net charge transfer of 0.46e from tyramine to the silver cluster.

Consequently, these findings provide new insights into the interactions between tyramine and silver nanoparticles, laying the ground work for further research. The study's results have implications for the development of drug delivery systems, biosensors, and other formulations with enhanced efficiency. The investigation encourages future experimental studies on the functionalization of silver nanoparticle clusters with bio-molecular species.