Abstract

Pyrazine is a kind of compound containing pyrazine ring structure, which has attracted the attention of researchers because of its special aroma and healthcare function. It is currently used in spices, pharmaceuticals, food additives, and other industries. Human serum albumin (HSA) is a transport and storage protein in the human circulatory system. Previous research reported that pyrazine can interact with albumin, but the interaction between pyrazine and HSA has not been reported. Consequently, this study is based on the multispectral method and molecular dynamics (MD) simulation, exploring the interaction between different pyrazines and HSA. The results show that the four pyrazines mainly quench the endogenous fluorescence of HSA through the static quenching mechanism. Molecular docking indicates that hydrophobic forces play a major role in the binding of pyrazines with HSA. At the same time, synchronous fluorescence spectroscopy, circular dichroism (CD), three-dimensional fluorescence spectroscopy, and resonance light scattering (RLS) showed that each pyrazine compound can change the conformation of HSA. MD simulation results show that the combination of four pyrazines and HSA can enhance the stability of HSA. Among them, 2,3,5,6-tetramethylpyrazine (TTMP) and 2,3,5-trimethylpyrazine (TMTP) binding is the most stable and the most unstable combination with HSA, respectively. In addition, the four pyrazines could affect the physiological activity of HSA. This study is of great significance to explore the mechanism of action of pyrazine compounds in vivo and can provide useful information for the application of pyrazine compounds in food.

1. Introduction

Human serum albumin (HSA) is a highly water-soluble spherical monomeric plasma protein, which has a variety of important functions and roles in the body. HSA can delay the metabolism or excretion of lipophilic substances in the human body and increase their solubility and half-life in plasma [1]. And it has been used as a biopharmaceutical in the clinical treatment of low blood volume, bleeding, shock, trauma, burns, and other diseases [2]. After some small drug molecules enter the human blood system, they reversibly bind to HSA before being stored and transported to exert their pharmacological effects [3, 4]. In addition, HSA also has the functions of osmotic adjustment, immune regulation, and maintaining acid-base balance [5, 6]. Therefore, HSA is widely used as a model protein for the study of protein structure.

Pyrazine is a kind of compound containing a pyrazine ring structure, which is widely found in nature, and can be utilized by microorganisms, plants, animals, and humans [7]. Currently, pyrazine and its derivatives are widely used in the spice, pharmaceutical, and food additive industries [8, 9]. 2,3-Dimethylpyrazine (2,3-DMTP), 2,5-dimethylpyrazine (2,5-DMTP), 2,3,5-trimethylpyrazine (TMTP), and 2,3,5,6-tetramethylpyrazine (TTMP) are pyrazine compounds with different position substituents (Figure 1(b)). These pyrazines have a pleasant nutty and barbecue flavour that produces a unique aroma [10], which is widely used in the food industry [1113]. At the same time, these pyrazines are also important flavour substances in liquor, which can be detected in fermented grains by GC-MS [14, 15] (Figure 1(a)). Except for aroma contribution, these pyrazine compounds also play an important role in the pharmaceutical industry. 2,5-DMTP is an important intermediate for the synthesis of new drugs such as amoxicillin and can be used to treat various diseases [16]. TTMP and its derivatives can capture superoxide anion and reduce the production of nitrogen oxide in human granulocytes. It is a clinically active ingredient used in the treatment of cardiocerebral vascular disease [17]. Studies have shown that there is an interaction between pyrazine drugs and albumin [18]. Some pyrazine drugs can bind to albumin to form drug-protein binding complexes, thereby changing the distribution of drugs and tissue targeting [19, 20].

However, there is no report on the interaction between pyrazines and HSA. This study investigated the interaction between pyrazine-flavoured compounds containing different substituents and HSA through multispectral analysis, thermodynamic analysis, binding sites, and molecular dynamics simulation and compared them. In addition, the effects of different pyrazines on the physiological activity of HSA were also investigated. Not only it is of great significance for revealing in vivo pharmacokinetics, guiding clinical rational drug use, and developing new drugs but it can also provide useful information for the application of pyrazines in food.

2. Materials and Methods

2.1. Materials

HSA (≥96%) was provided by Sigma-Aldrich Company (Shanghai, China). All the pyrazine flavour substances (≥98%, in pure substance form) were purchased from Macklin Biochemical Co., Ltd. (Shanghai, China). Tris-HCl (pH 7.4, 1.0 mol·L−1) was purchased from Solarbio Bioscience & Technology Co., Ltd. (Shanghai, China). Dimethyl sulfoxide (DMSO, 99%) was purchased from Shanghai Aladdin Biochemical Technology Co., Ltd. (Shanghai, China). HSA was dissolved in Tris-HCl (0.05 mol·L−1, pH = 7.40) buffer to obtain 1.00 × 10−4 mol·L−1 of HSA. 2,3-DMTP and 2,5-DMTP were dissolved in pure water to obtain solutions with concentrations of 1.875 × 10−3 mol·L−1. TMTP and TTMP were dissolved in ethanol and water, respectively, to obtain the solution with concentrations of 1.875 × 10−3 mol·L−1.

2.2. Fluorescent Quenching

A certain volume of 2,3-DMTP, 2,5-DMTP, TMTP, and TTMP (1.875 × 10−3 mol·L−1) was added to a cuvette containing 2.5 ml HSA (1.00 × 10−4 mol·L−1). The 2,3-DMTP, 2,5-DMTP, TMTP, and TTMP concentrations were 0, 7.50, 15.00, 22.50, 30.00, 37.50, 45.00, 52.50, 60.00, and 67.50 (×10−6 mol·L−1). The excitation (Ex) wavelength was 285 nm, and the emission (Em) wavelength was recorded from 300 to 420 nm. The slit widths of the Ex and Em wavelengths were set at 15 nm. The fluorescence emission spectra of HSA in different reaction systems were measured at 3 different temperatures (291, 301, and 310 K).

2.3. Three-Dimensional Fluorescence Spectrum

50 μL of 2,3-DMTP, 2,5-DMTP, TMTP, and TTMP was added to a cuvette containing 2.5 ml of HSA. The Ex was 200 nm, and the Em was set from 200 to 500 nm, and scanned 30 times. The three-dimensional fluorescence spectrum was measured at 301 K.

2.4. Circular Dichroism Spectrum

A certain volume of pyrazine (1.875 × 10−3 mol·L−1) was added to a quartz cell (a quartz cell of 2 mm path length) containing 350 μL of HSA (2.00 × 10−6 mol·L−1) with an optical path of 2 mm. The 2,3-DMTP, 2,5-DMTP, TMTP, and TTMP concentrations were 0, 0.46, 0.94, and 1.41 (×10−4 mol·L−1). The spectra of HSA in different systems in the range of 200–300 nm were collected at 301 K. Each spectrum presented was the average of three consecutive scans at a scan speed of 100 nm·min−1 and a bandwidth of 1 nm.

2.5. Molecular Dynamics Simulation

Molecular dynamics (MD) simulation was performed for the complex using the GROMACS 2022 software and AMBER 19 force field. Initially, the periodic boundary condition was applied and a cubic box was placed around all atoms with a distance of 10 Å. TIP3P water molecules with a density of 0.10 g·mL−1 were used to solvate the system. The simulation temperature and pressure were controlled at 298 K and 1 bar, respectively. The long-range electrostatic interactions were calculated using the particle mesh Ewald method with an 8 Å cutoff. The time step was 2 fs, and simulation snapshots were saved every 25 ps during the MD simulation. The MD simulation of each system was performed for 200 ns. Finally, the analysis of complex RMSD, SASA, , and protein RMSF was obtained from the MD trajectories, and the analysis of binding energy was conducted using gmx_MMPBSA.

2.6. The Effect of Different Pyrazines on the Physiological Activity of HSA

HSA was mixed with different concentrations of pyrazines for 30 min, and pNPA was added to measure the absorbance of the reaction product at 400 nm. At the same time, the esterase activity of pure HSA was determined under the same conditions [21, 22]. The effects of different pyrazines on HSA esterase activity were obtained. Ellman’s reagent was used instead of pNPA to determine the absorbance of the reaction product at 412 nm. The free radical scavenging ability of pure HSA was determined under the same conditions, and the effects of different pyrazine compounds on the free radical scavenging ability of HSA were obtained.

2.7. Data Processing

The internal filter effect is an important reason for reducing fluorescence intensity. When the absorbance is less than 0.3, the internal filter effect is corrected by the following formula:where and are the corrected and observed fluorescence intensities, respectively, and and are the absorption of the system at the Ex and Em wavelengths, respectively.

3. Results and Discussion

3.1. The Fluorescence Quenching and Quenching Mechanism

Fluorescence quenching refers to the decrease in fluorescence intensity caused by the interaction between fluorescent material molecules and solvent molecules [23]. Figure 2 shows the fluorescence quenching spectra of various pyrazines on HSA at 301 K. With the increase of pyrazine, the fluorescence peak shape of HSA did not change significantly, but the fluorescence intensity decreased gradually. After interaction with four pyrazines, the maximum absorption wavelength of HSA showed different degrees of blueshift. The results showed that 2,3-DMTP, 2,5-DMTP, TMTP, and TTMP could interact with HSA and induce its fluorescence quenching [24]. Fluorescence quenching can be divided into static quenching and dynamic quenching [25]. The fluorescence quenching process is analyzed by the Stern–Volmer equation (2) [26] as follows:where and F are the fluorescence intensities of HSA without quencher and with quencher, respectively is the dynamic quenching constant, [Q] is the concentration of the quencher, is the bimolecular quenching rate constant, and is the average life of HSA without quencher (a value of 10−8 s for HSA).

The Stern–Volmer plots for HSA fluorescence quenching by different concentrations of 2,3-DMTP, 2,5-DMTP, TMTP, and TTMP are displayed in Figure 3. The values of and could be obtained from the slopes of a linear adjustment, and the results are summarised in Table 1. The value decreases with the increase in temperature. Meanwhile, the values after the interaction of pyrazines and HSA were greater than the maximum scatter collision quenching constant (2.00 × 1010 L·mol−1·s−1), which indicated that fluorescence quenching of HSA is caused by static quenching [27]. In addition, at the same temperature, the interaction between TTMP and HSA is the strongest.

3.2. UV-Vis Absorption Spectrum

The UV-visible absorption spectrum is caused by the transition of valence electrons after molecules/ions absorb the ultraviolet light or visible light [28]. In the protein molecule, the benzene ring of Tyr, Phe, and Trp residues contains conjugated double bonds, which makes the protein have ultraviolet absorption at 280 nm [29]. The exact mechanism of the quenching process of different pyrazines and HSA was further revealed by the UV-vis absorption spectra (Figure 4). Upon the addition of pyrazines, there were distinct changes in the absorption intensity and peak position of HSA, and the maximum absorption wavelength was blue-shifted. It is indicated that pyrazine compounds can change the conformation of HSA and increase the hydrophobicity of aromatic amino acid hydrophobic groups in HSA molecules [30]. Since the dynamic quenching only affects the excited fluorescent molecules, it does not change the absorption intensity of the absorption spectrum. However, static quenching can interact with fluorescent molecules to form complexes, causing changes in the absorption intensity of fluorescent molecules [31]. Therefore, the quenching process of different pyrazines and HSA is static quenching.

3.3. Combining Properties and Thermodynamic Parameters

The binding constant () is a parameter for quantifying the interaction between small molecules and proteins. The binding site (n) refers to the specific site in the molecule that can form a stable interaction with the ligand. The and n can be calculated by equation (3) [32].

There are four fundamental forces that affect the interaction between proteins and ligands, namely, hydrogen bond, electrostatic force, hydrophobic interaction, and van der Waals force, which can be evaluated by thermodynamic parameters. The thermodynamic parameters were calculated by Van’t Hoff equations (4)–(6) [33].where is the binding constant that reflects the strength of the interaction of proteins acting with small molecules at the corresponding temperatures, n is the number of binding sites, [Q] is the concentration of FA or FAM, T is the temperature, and R is the gas constant. ∆H, ∆S, and ∆G are the enthalpy change, entropy change, and Gibbs free energy, respectively.

Figure 5 shows the double logarithmic curve of the interaction between pyrazine compounds and HSA. Based on the slope and constant, the binding site and the binding constant can be calculated (Table 2). increases with the temperature, indicating that the binding ability of pyrazine compounds to HSA is strong. At the same temperature, the binding force of pyrazine compounds with HSA can be described as follows: TMTP < 2,5-DMTP < 2,3-DMTP < TTMP. This result indicated that the more substituents, the easier it is to bind to HSA. The position of the substituent also affects the binding effect with HSA. Compared with the parasubstituent, the orthosubstituent is easier to bind to HSA. In addition, n is close to 1, and there is only one binding site between HSA and pyrazine compounds. The thermodynamic parameters are shown in Table 2. ∆G < 0, indicating that the reaction can proceed spontaneously. ΔH > 0, ΔS > 0, indicating that hydrophobic force is the main force between pyrazine compounds and HSA [34], which once more proves that pyrazine compounds can interact with HSA.

3.4. Energy Transfer in the Interaction between HSA and Pyrazines

Energy transfer can be used to evaluate the distance between the ligand and Trp residue in the protein [35]. Based on Förster’s nonradiative energy transfer theory, the distance between pyrazines and HSA luminescent groups can be calculated by using the following equations [36, 37]:where is the dipole spatial orientation factor, N is the refractive index of the medium, Φ is the fluorescence quantum yield of the donor in the presence of no acceptor, F (λ) is the fluorescence intensity of the donor at wavelength λ, and ε (λ) is the molar absorption of the acceptor at wavelength λ. Among them,  = 2/3, N = 1.336, and Φ = 0.118.

When the concentration of pyrazines and HSA is 1 : 1, the overlap of HSA fluorescence spectra and different pyrazine absorption spectra is shown in Figure 6. The relationship between energy overlap integral (J), energy transfer efficiency (E), critical distance (), and binding distance (r) is shown in Table 3. The binding distances of 2,3-DMTP, 2,5-DMTP, TTMP, and TMTP with HSA were calculated to be 2.99 nm, 2.82 nm, 3.49 nm, and 4.22 nm, respectively. r < 7 nm and 0.5 < r < 1.5 indicate that the chromophore in HSA can transfer energy to pyrazine compounds by nonradiation. of the interaction between 2,3-DMTP, 2,5-DMTP, TMTP, and HSA indicated that the contribution of nonradiative energy transfer to fluorescence quenching was small in these three systems. In summary, the endogenous fluorescence quenching mechanism of pyrazine compounds on HSA is mainly static quenching, which is consistent with the conclusion of the fluorescence quenching mechanism [38].

3.5. Synchronous Fluorescence Spectrum

Synchronous fluorescence measurement can provide microenvironment information on molecules near fluorophore functional groups [39]. When Δλ = 15 nm and 50 nm, the synchronous fluorescence reflects the fluorescence characteristics of Tyr and Trp residues, respectively. As shown in Figure 7, the fluorescence intensity of HSA decreased with the increase in pyrazine concentration, that is, the change in amino acid environment and hydrophobicity of HSA. At ∆λ = 15 nm, the maximum emission wavelength of HSA in the TMTP-HSA, 2,3-DMTP-HSA, 2,5-DMTP-HSA, and TTMP-HSA systems showed a significant redshift, indicating that 2,3-DMTP, 2,5-DMTP, TMTP, and TTMP increase the polarity around Tyr residues in HSA [40, 41]. At ∆λ = 50 nm, the maximum emission wavelength of HSA in the 2,3-DMTP-HSA and 2,5-DMTP-HSA systems showed a significant redshift, and there is no significant change in the maximum emission wavelength of HSA interacting with TTMP and TMTP, indicating that 2,3-DMTP and 2,5-DMTP increase the polarity around Tyr residues in HSA. In summary, pyrazine compounds can enhance the polarity of the microenvironment around amino acid residues in HSA, change the secondary structure of HSA protein, and lead to changes in the conformation of HSA. Moreover, the redshift of HSA was more significant when ∆λ = 15 nm, indicating that the fluorescence quenching intensity of pyrazine compounds to Tyr was greater than that to Trp.

3.6. Circular Dichroism (CD)

In order to further explore the effect of pyrazine compounds on the conformation of HSA, the circular dichroism (CD) of different pyrazine compounds on HSA was determined. As shown in Figure 8, two negative peaks appear in the spectrum, and it is the characteristic peak of the α-helix structure in HSA [42]. The intensity of the two peaks gradually decreases as the concentration of pyrazines in the HSA solution increases. The SELCON3 algorithm was used to calculate the proportion of different types of secondary structures in HSA molecules with the addition of pyrazine compounds. The results are shown in Table 4. With the addition of pyrazines, the α-helix structure in HSA molecules decreased, while the content of β-Sheet, β-Turn, and random coil structure in protein molecules increased. The change of secondary structure indicated that pyrazine compounds changed the microenvironment around HSA and changed the secondary structure of HSA. At the same time, this further indicates that the interaction between pyrazine and amino acid residues on the peptide chain of HSA destroys the network structure of hydrogen bonds and loosens the conformation of HSA.

3.7. Three-Dimensional Fluorescence Spectrum

The three-dimensional fluorescence spectroscopy can simultaneously obtain fluorescence intensity information from two different dimensions of excitation and emission, and the changes in HSA structure can be observed more comprehensively and intuitively [43]. Figure 9 shows the three-dimensional fluorescence spectrum of 2,3-DMTP, 2,5-DMTP, TMTP, and TTMP interacting with HSA. Among them, peak 1 and peak 2 are Rayleigh scattering peaks. The peak a is the fluorescence peak of Trp and Tyr residues, which mainly reflects the changes in Trp and Tyr spectra caused by transitions. The peak b is the characteristic peak of the peptide backbone, which reflects the spectral changes of the HSA peptide backbone caused by transitions [44]. It can be seen from the figure that the addition of pyrazine significantly reduced the fluorescence intensity of peak a and peak b. At the same time, the emission wavelength corresponding to the maximum absorption peak of peak b shifted significantly. It is fully demonstrated that pyrazine compounds lead to changes in the spatial conformation of HSA, which changes the internal microenvironment of the molecule and reduces the luminous efficiency of Trp and Tyr residues.

3.8. Resonance Light Scattering (RLS)

When the excitation light is located in or close to the molecular absorption band, the electron strongly absorbs the light energy due to resonance, resulting in resonance light [45]. RLS is very sensitive to electrostatic interactions, hydrogen bonding, and hydrophobic interactions [46], which can be used to characterize changes in the concentration and aggregation state of the species. With the addition of pyrazine compounds, the RLS intensity of HSA showed a downward trend (Figure 10). After the combination of pyrazine compounds with HSA, the protective water layer on the surface of HSA was destroyed, which made the originally more dispersed protein become more dispersed, and the protein size in HSA decreased along with the degree of aggregation.

3.9. Molecular Docking

Molecular docking is the process of identifying two or more molecules by geometric matching and energy matching to find the best matching mode, which is helpful in further understanding the interaction between molecules [47]. The docking energy (Table 5) and interaction (Figure 11) of different substituent pyrazines binding to HSA can be obtained by molecular docking. The docking energy represents the binding capacity of each pyrazine compound and HSA, and the lower the docking energy, the more stable the binding between the ligand and the receptor [48]. It is generally believed that the docking energy value less than −4.25 kcal/mol−1 indicates that there is a certain binding activity between the two, less than −5.0 kcal/mol−1 indicates a good binding activity, and less than −7.0 kcal/mol−1 indicates a strong binding activity [49]. It can be seen from the table that the binding energy of different substituent pyrazine compounds and HSA is between −7.0 kcal·mol−1 and −5.0 kcal·mol−1, and the binding activity of pyrazines to HSA is TTMP > TMTP > 2,3-DMTP > 2,5-DMTP. This also proves that the more the substituents of pyrazine compounds, the easier it is to bind to HSA. Compared with the parasubstituent, the orthosubstituent is easier to bind to HSA.

During the interaction, four pyrazines formed hydrophobic interactions with Ala158, Tyr138, Ile142, Tyr161, and Phe165 residues in HSA (Figure 11). In addition, 2,3-DMTP, TMTP, and TTMP also formed hydrophobic interactions with Phe134 residues (Figures 11(a), 11(c), and 11(d)). TMTP formed a hydrophobic interaction with Met123 (Figure 11(c)). It can be seen from the above that 2,3-DMTP, 2,5-DMTP, TMTP, TTMP, and HSA are mainly combined by hydrophobic force. This further verifies the thermodynamic results. In the part of three-dimensional fluorescence spectroscopy, we found that the fluorescence peak of Trp and Tyr residues decreased significantly. Furthermore, molecular docking proved that the contribution of Tyr was greater than that of Trp in the interaction with pyrazine compounds.

3.10. Molecular Dynamics (MD) Simulation

The trend of root-mean-square displacement (RMSD) of proteins and ligands is an important characterization to determine whether the simulation is stable [50]. The height of RMSD reflects the overall deviation of the trajectory from the reference molecule [51]. The RMSD calculation results are shown in Figure 12(a). The overall structure of the complex was balanced in an MD simulation time of 200 ns. The RMSD of free HSA fluctuates more than the complex. At the same time, HSA after binding to different pyrazines is more stable than protein (Figure 12(a-1)). It is shown that pyrazines made the structure of HSA more stable. In addition, the RMSD values of HSA and TTMP complexes were lower and fluctuated less than those of other pyrazines and HSA complexes, indicating that TTMP and HSA bound stronger. TMTP binds weakest to HSA, and 2,5-DMTP and 2,3-DMTP are not much different. These are consistent with the previous analysis of the associative constant.

To gain insight into the conformational changes of HSA, root-mean-square fluctuations (RMSFs) were used to measure the intensity of change in the secondary structure of HSA [52]. Compared with the RMSF of different complexes after pyrazines binding to HSA, the RMSF fluctuations simulated by free HSA molecular dynamics were greater (Figure 12(b)), indicating that HSA was more stable after binding to pyrazines. At the same time, the RMSF value decreased to a certain extent after HSA was combined with different pyrazines. Among them, the RMSF value of the complex was lower than that of other pyrazine compounds after HSA was combined with TTMP, indicating that TTMP greatly increased the stability of HSA. The reason for this result may be related to the number and position of substituents.

The radius of rotation () is a physical quantity that describes the structural compactness of lactoglobulin, and smaller values indicate greater lactoglobulin density [53]. reacts to the volume and shape of proteins, helping to study protein compactness [54]. simulated by free HSA dynamics was greater than that of the complex, indicating that the complex was more stable and the system was tighter after binding to different pyrazines (Figure 12(c)). Among them, TTMP was the most tightly bound to HSA, followed by 2,3-DMTP and 2,5-DMTP, and TMTP was the least tightly bound to HSA. This structure was consistent with RMSD and RMSF analyses.

Solvent-accessible surface area (SASA) is the surface area of biomolecules accessible to solvents [55]. The SASA of the complex was analyzed (Figure 12(d)). After different pyrazines and HSA binding, the SASA of the complex decreased compared with the SASA of free HSA, indicating that the protein structure stability of HSA increased. At the same time, the conformations before and after MD are superimposed to generate porcupine diagrams to represent the changes in protein structure after small molecule binding (Figure 13) and the arrow indicates the movement direction of the protein amino acid residue, and the length of the arrow represents the size of the travel distance. After small molecules bind to the active site of HSA, HSA has strong movement, indicating that different pyrazines have a greater effect on the structure of HSA.

3.11. The Effect of Different Pyrazines on the Physiological Activity of HSA

HSA has esterase activity, can hydrolyze pNPA to produce p-nitrophenol, and has the maximum absorption at 400 nm [56]. The esterase activity of HSA decreases with the increase of pyrazine concentration (Figure 14(a)). The reason for this result is that the addition of pyrazines increases the folding degree of HSA, exposing some hydrolyzed pNPA groups in its structure, resulting in reduced HSA esterase activity [57]. Among them, TTMP decreased the fastest, indicating that TTMP had the greatest effect on HSA esterase activity, which was related to the stability of TTMP and HSA binding.

Ellman’s reagent can react with the thiol group of the free amino acid cysteine (Cys34) in HSA to form a yellow complex with a maximum absorption at 412 nm [58]. With the increase in the concentrations of 2,3-DMTP, 2,5-DMTP, and TMTP, the absorbance showed a downward trend (Figure 14(b)). It indicates that the interaction of 2,3-DMTP, 2,5-DMTP, and TMTP can reduce the radical scavenging capacity of HSA. The speculated reason for this may be that during the interaction between pyrazines and HSA, part of the active sites of HSA are occupied by small molecules of pyrazines, which makes it unable to combine with free radicals, thus reducing its antioxidant capacity. At the same time, the steric hindrance after combining with pyrazines is greater, which is not conducive to the free binding of HSA with free radicals [59, 60]. Contrary to the results of esterase activity, the free radical scavenging ability of HSA decreased only after TTMP reached a certain concentration, and TTMP had the least effect on the free radical scavenging ability of HSA.

4. Conclusions

This research employed the multispectral method and molecular dynamics simulation to investigate the interaction between different substituents of pyrazine and HSA. Binding constants, thermodynamic parameters, and molecular docking analysis showed that pyrazine mainly quenched the endogenous fluorescence of HSA through a static quenching mechanism. The hydrophobic force plays a major role in the binding of pyrazine compounds to HSA. According to Förster’s nonradiative energy transfer theory, the binding distances of 2,3-DMTP, 2,5-DMTP, TTMP, and TMTP with HSA were calculated to be 2.99 nm, 2.82 nm, 3.49 nm, and 4.22 nm, respectively. The nonradiative energy transfer occurs between HSA and different pyrazines. The docking energy of different pyrazines with HSA is demonstrated as follows: TTMP > TMTP > 2,3-DMTP > 2,5-DMTP. It was between −7.0 kcal·mol−1 and −5.0 kcal·mol−1, indicating that HSA had good binding activity. Meanwhile, MD simulation results show that the combination of four pyrazines and HSA can enhance the stability of HSA. Among them, TTMP and HSA binding is the most stable, and TMTP is the most unstable combination with HSA. 2,3-DMTP and 2,5-DMTP have roughly the same effect on HSA stability. In addition, different pyrazines and HSA can reduce the esterase activity and free radical scavenging ability of HSA. It indicated that pyrazine could affect the physiological activity of HSA. Pyrazine compounds are widely used in food, spices, pharmaceuticals, and other industries. This paper deeply studies the interaction between pyrazines and HSA, which is helpful in understanding the reasons for the differences between pyrazines and protein functions in blood transport and in in vivo and in vitro biological activities, and provides a theoretical basis for the rational use of pyrazines and additives in food.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Wenhua Tong and Shuqin Wang designed the study and wrote the manuscript. Zongwei Qiao and Ying Yang supervised the study and revised the manuscript. Wenhua Tong, Shuqin Wang, and Haoran Xu collected the data and visualized the study. Haoran Xu and Liming Zhao analyzed the data. All authors reviewed the manuscript and have read and approved the final version of the manuscript.

Acknowledgments

This research was funded by the Key Research and Development Project of Sichuan Provincial Science and Technology Department (no. 2022YFS0549), the Special Support of Sichuan Province 2022 Postdoctoral Research Project (no. 294074), the Talent Introduction Project of Sichuan University of Science and Engineering (no. 2019RC31), and the Graduate Innovation Fund of Sichuan University of Science and Engineering (nos. Y2022090 and Y2023246).