**1. Introduction**

HSA is a major transporter and the most abundant protein in the plasma. It is responsible for balancing the osmotic pressure and is a major regulator of fluid distribution in the body. HSA exhibits astonishing ligand-binding capabilities, acting as a warehouse and transporter of many endo- and exogenous compounds [1]. This promiscuous, nonspecific affinity can lead to sudden changes in concentrations caused by displacement when two or more compounds compete for binding to the same molecular site [2]. HSA, a major carrier in the human circulatory system performs a key function of circulating various compounds such as drugs [3]. HSA interactions with drugs alter proteins' pharmacokinetics and pharmacodynamics as well. Free drugs diffuse inactively across the membranes via specific

**Citation:** Shamsi, A.; Shahwan, M.; Khan, M.S.; Alhumaydhi, F.A.; Alsagaby, S.A.; Al Abdulmonem, W.; Abdullaev, B.; Yadav, D.K. Mechanistic Insight into Binding of Huperzine A with Human Serum Albumin: Computational and Spectroscopic Approaches. *Molecules* **2022**, *27*, 797. https://doi.org/ 10.3390/molecules27030797

Academic Editor: Anna Maria Almerico

Received: 1 December 2021 Accepted: 14 January 2022 Published: 25 January 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

transporters to interact with their respective protein targets [4]. A pivotal step in the domain of drug discovery is the investigation of the pharmacokinetics and pharmacodynamics of drugs. Pharmacological profiling of drugs offers understanding of the interactions of vital therapeutic drugs or derivatives with either plasma or target tissue proteins [5]. The binding of a drug to HSA is a critical factor determining its pharmacological profiling and distribution [6]. In medicinal chemistry, studies pertaining to plasma proteins and drugs binding are attracting researchers across the globe, because these studies provide a platform to study drugs' behavior and action, thereby delineating their transport and distribution characteristics in the circulatory system. Furthermore, it is also vital to study the protein–protein or protein–drug interactions to make progressive inroads in the advancements made in pharmaceutical industry [7]. HSA is the main plasma protein, and it is imperative to study binding of drugs with HSA. Thus, this study aimed at investigating the binding mechanism of HpzA to HSA. HSA is a one-chained polypeptide weighing 66.5 kDa. The 585 amino acid polypeptide assembles in a heart-shaped structure [1]. The protein displays two specific sites for drug binding as specified using Sudlow's fluorescent probe displacement method. According to the Sudlow's classification, there are two main binding sites for drug ligands of HSA, namely, subdomain IIA (site IIA) and IIIA (site IIIA) [8]. Site IIA, positioned in subdomain IIA, concerns with the selective binding of heterocyclic anions, while aromatic carboxylates bind to the site IIIA that is located in subdomain IIIA,. For example, ibuprofen binds to site IIIA and warfarin to site IIA [9]. The binding of all the ligands cannot be considered under Sudlow's model. Site IB, a third binding pocket, located in subdomain IB, is a hydrophobic D-shaped cavity, is also used for interaction by some compounds such as bilirubin [10]. A study [11] suggested site IB as the third major site having the potential to bind drug ligands of HSA.

AD, the sixth leading cause of claiming lives worldwide, is a progressive, neurodegenerative disorder. The condition is characterized by impairments in mental functions, amnesia, and dementia [12]. The major hallmarks of the condition are accumulated amyloid β-peptide (Aβ) in plaques, the formation of neurofibrillary tangles (NFTs), and degeneration of neurons. The mechanism by which neurons degenerate is still unclear and needs extensive research. Aβ accumulation and deposition can cause neurodegeneration by many mechanisms, such as inflammation, oxidative stress, and apoptosis. In addition to the genetic factors, environmental factors such as neurotoxins, stress, etc. play a major role in AD progression [13]. Reactive oxygen species (ROS) creates a condition of oxidative stress that leads to cell death, oxidative bursts, accumulation of spare free metals, etc. These have been hypothesized as major mediators in AD progression and neurodegeneration. Postmortem reports have indicated apoptosis as a major event in neurodegeneration and AD [14].

HpzA, a Lycopodium alkaloid, is a Chinese medicine isolated from *Qiang Ceng Ta*, the whole plant of *Huperzia serrata*. The plant belongs to the Huperziaceae family, and it has gained wide popularity due to its anticholinesterase (AChE) and anti-AD properties. The alkaloid also acts as a remedy for various ailments such as strains and swellings and neurological disorders such as schizophrenia and neurodegenerative disorders [15]. HpzA, a well-known AChE inhibitor, has properties superior to those of FDA-approved drugs such as donepezil, galantamine, etc. [16]. Figure S1 describes the 2D and 3D structural features of HpzA. Extensive clinical trials for the HpzA conducted in China showed enhanced memory in elderly individuals, subjects with issues related to amnesia, and AD patients. The naturally extracted ligand has also been through various extensive clinical trials on vascular dementia, showing improved cognitive function in subjects with AD [17]. In this introductory study, we report the plausible binding of HpzA with HSA. A molecular docking study was performed primarily to explore the possible interactions and binding affinity of HpzA towards HSA and was followed by MD simulation studies for 100 ns. Computational studies suggested plausible binding of HpzA with HSA, which was further validated by fluorescence-based binding assays. Thus, this study delineates the mechanism of binding of HpzA with HSA by employing spectroscopic and in silico approaches.

### **2. Results 2. Results**  *2.1. Molecular Docking*

approaches.

### *2.1. Molecular Docking* The functional activity of a protein changes broadly upon conformational alterations

*Molecules* **2021**, *26*, x 3 of 13

The functional activity of a protein changes broadly upon conformational alterations induced by ligand binding. The molecular docking approach can systematically consider potential binding modes for a protein–ligand system [18]. The blind docking approach was used to find all possible interactions between HpzA and HSA. HpzA has multiple binding sites on has, with varying affinities. All the docking sites and their docking energies are depicted in Figures S2 and S3, respectively. Out of nine docked conformations of HpzA in the splitting process, the preferable docked conformation was taken and explored in detail. The selected conformation of the docked HpzA showed a considerable binding score of −7.2 kcal/mol. It had LE and pKi values 0.4 (kcal/mol/non-H atom) and 5.28, respectively. The docking score affirmed HpzA as a plausible binding partner of HSA. Figure 1 shows the binding mode of HpzA with HSA. HpzA was shown to have a decent structural complementarity with good binding affinity. It showed various interactions with the key residues of the HSA binding pocket, such as one hydrogen bond with His266 and a few other noncovalent interactions with important HSA residues (Figure 1). The binding site of HpzA is located in domain II of HSA, but Tyr172, Tyr174, and Glu177 of domain I also participate in the interaction. HpzA binds at the Sudlow site I of HSA. The surface representation of HSA showed HpzA fitted inside the deep cavity of the binding site, which might cause significant changes to the conformational activity of HSA (Figure 1). induced by ligand binding. The molecular docking approach can systematically consider potential binding modes for a protein–ligand system [18]. The blind docking approach was used to find all possible interactions between HpzA and HSA. HpzA has multiple binding sites on has, with varying affinities. All the docking sites and their docking energies are depicted in Figures S2 and S3, respectively. Out of nine docked conformations of HpzA in the splitting process, the preferable docked conformation was taken and explored in detail. The selected conformation of the docked HpzA showed a considerable binding score of −7.2 kcal/mol. It had LE and pKi values 0.4 (kcal/mol/non-H atom) and 5.28, respectively. The docking score affirmed HpzA as a plausible binding partner of HSA. Figure 1 shows the binding mode of HpzA with HSA. HpzA was shown to have a decent structural complementarity with good binding affinity. It showed various interactions with the key residues of the HSA binding pocket, such as one hydrogen bond with His266 and a few other noncovalent interactions with important HSA residues (Figure 1). The binding site of HpzA is located in domain II of HSA, but Tyr172, Tyr174, and Glu177 of domain I also participate in the interaction. HpzA binds at the Sudlow site I of HSA. The surface representation of HSA showed HpzA fitted inside the deep cavity of the binding site, which might cause significant changes to the conformational activity of HSA (Figure 1).

was further validated by fluorescence-based binding assays. Thus, this study delineates the mechanism of binding of HpzA with HSA by employing spectroscopic and in silico

**Figure 1.** The binding of HpzA with HSA. Ribbon and surface representation of HSA protein showing various interactions and docking fit of HpzA (UniProt ID: P02768). **Figure 1.** The binding of HpzA with HSA. Ribbon and surface representation of HSA protein showing various interactions and docking fit of HpzA (UniProt ID: P02768).

The selected docking pose of HpzA with HSA was further subjected to a detailed analysis of possible interactions. The generated 2D plot showed that HpzA occupied the hydrophobic pocket of HSA bordered by residues Tyr174, Glu177, Lys219, Gln220, Lys223, Leu262, His266, Cys269, Arg281, His312, and Ala315 (Figure 2A). The plot suggested that HpzA was tightly bound with HSA with the help of one hydrogen bond with His266 and several hydrophobic contacts. HpzA showed multiple van der Waals interactions with Glu177, Gln220, Leu262, and His312 of HSA. The phenolic moiety at the end of HpzA occupied a crucial binding cleft of HSA formed by different residues (Figure 2B). The selected docking pose of HpzA with HSA was further subjected to a detailed analysis of possible interactions. The generated 2D plot showed that HpzA occupied the hydrophobic pocket of HSA bordered by residues Tyr174, Glu177, Lys219, Gln220, Lys223, Leu262, His266, Cys269, Arg281, His312, and Ala315 (Figure 2A). The plot suggested that HpzA was tightly bound with HSA with the help of one hydrogen bond with His266 and several hydrophobic contacts. HpzA showed multiple van der Waals interactions with Glu177, Gln220, Leu262, and His312 of HSA. The phenolic moiety at the end of HpzA occupied a crucial binding cleft of HSA formed by different residues (Figure 2B). As a whole, the HpzA binding with HSA suggested that it could disrupt the structural conformation of the protein, which may further result in altering HSA's binding capability.

**Figure 2.** Detailed interactions of HSA with HpzA as (**A**) 2D and (**B**) 3D plots. **Figure 2.** Detailed interactions of HSA with HpzA as (**A**) 2D and (**B**) 3D plots.

#### *2.2. MD Simulations 2.2. MD Simulations*

ity.

Many studies have used the MD simulation approach to investigate ligand binding with proteins [7,19–21]. All-atom MD simulations were carried out for the docked complex of HSA–HpzA and compared with those for the free state of HSA. Plots for root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (*R*g), and solvent accessible surface area (SASA) fluctuations, along with hydrogen bond analysis, were generated and analyzed in detail. The analysis of RMSD assists in understanding the stability of a protein and protein–ligand complex. To explore structural deviations in HSA, we used RMSD examination of the simulated protein and protein–ligand complex. As examined, both HSA and the HSA–HpzA complex reached the equilibrium phase without any major shift in the RMSD pattern (Figure 3A, left panel). The protein–ligand complex was stable throughout 100 ns simulation. The analysis of the generated plots suggested a small amount of fluctuation at the beginning in the RMSD values of the HSA– HpzA that was possibly due to the initial adjustment of the system. As a whole, HSA showed a small number of fluctuations in its backbone in the simulation, but no significant shift was observed (Figure 3A, right panel). Exploring RMSFs is helpful in understanding the fluctuations in the protein residues during the course of the simulation [22,23]. The RMSFs were recorded for each residue in Many studies have used the MD simulation approach to investigate ligand binding with proteins [7,19–21]. All-atom MD simulations were carried out for the docked complex of HSA–HpzA and compared with those for the free state of HSA. Plots for root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (*R*g), and solvent accessible surface area (SASA) fluctuations, along with hydrogen bond analysis, were generated and analyzed in detail. The analysis of RMSD assists in understanding the stability of a protein and protein–ligand complex. To explore structural deviations in HSA, we used RMSD examination of the simulated protein and protein–ligand complex. As examined, both HSA and the HSA–HpzA complex reached the equilibrium phase without any major shift in the RMSD pattern (Figure 3A, left panel). The protein–ligand complex was stable throughout 100 ns simulation. The analysis of the generated plots suggested a small amount of fluctuation at the beginning in the RMSD values of the HSA–HpzA that was possibly due to the initial adjustment of the system. As a whole, HSA showed a small number of fluctuations in its backbone in the simulation, but no significant shift was observed (Figure 3A, right panel). *Molecules* **2021**, *26*, x 5 of 13

As a whole, the HpzA binding with HSA suggested that it could disrupt the structural conformation of the protein, which may further result in altering HSA's binding capabil-

**Figure 3.** Dynamics performance of HSA. (**A**) RMSD plot of HSA and its docked complex with HpzA. (**B**) Residual fluctuation plot of HSA and its docked complex with HpzA. **Figure 3.** Dynamics performance of HSA. (**A**) RMSD plot of HSA and its docked complex with HpzA. (**B**) Residual fluctuation plot of HSA and its docked complex with HpzA.

that of the *Rg* values of HSA and the HSA–HpzA docked complex, but with a minor increase (Figure 4A, left panel). This minor increase could have resulted from the occupancy of intramolecular space of HSA by HpzA. The distribution plot of *Rg* values also indicated a slight increase in the *Rg* of the HSA–HpzA docked complex compared to that of the free

SASA delineates the folding mechanism of a protein structure by analyzing the protein's solvent accessibility [26]. SASA values were calculated and plotted to investigate the surface area of HSA and the HSA–HpzA docked complex during the simulation, the (Figure 4B). The results indicated that the SASA values were minorly increased for HSA in ligand-bound state. This slight increase in the values of SASA indicated that a number of internal residues might be exposed to the surface post the binding of HpzA (Figure 4B, left panel). The distribution plot of SASA values also indicated a slight increase in SASA

Calculating and examining the *Rg* of a protein is useful to measure the shape of the protein based on its hydrodynamic radius [18,23–25]. The *Rg*s of HSA and the HSA–HpzA docked complex were calculated to assess their structural compactness during the simu-

of the HSA–HpzA docked complex (Figure 4B, right panel).

state of HSA (Figure 4A, right panel).

Exploring RMSFs is helpful in understanding the fluctuations in the protein residues during the course of the simulation [22,23]. The RMSFs were recorded for each residue in HSA and plotted for analysis purposes. The RMSF analysis emphasized possible residual movements in HSA in free and ligand-bound states (Figure 3B). The compared RMSF values between the free and ligand-bound HSA suggested that the residual fluctuations were decreased and compacted the ligand-bound state of HSA during the simulation (Figure 3B, left panel). The results suggested that the fluctuations in the residues were minimized upon binding of HpzA. However, a few increased fluctuations were observed in the HSA residues corresponding to the loop region far from the ligand site. RMSF values distribution plot revealed that the residual fluctuations were minimized after HpzA binding, which further confirmed a robust constancy of the protein–ligand complex (Figure 3B, right panel).

Calculating and examining the *R<sup>g</sup>* of a protein is useful to measure the shape of the protein based on its hydrodynamic radius [18,23–25]. The *Rg*s of HSA and the HSA– HpzA docked complex were calculated to assess their structural compactness during the simulation (Figure 4A). The average *R<sup>g</sup>* values of free HSA and the HSA–HpzA docked complex were 2.76 nm and 2.80 nm, respectively. The pattern of *R<sup>g</sup>* values of HSA was like that of the *R<sup>g</sup>* values of HSA and the HSA–HpzA docked complex, but with a minor increase (Figure 4A, left panel). This minor increase could have resulted from the occupancy of intramolecular space of HSA by HpzA. The distribution plot of *R<sup>g</sup>* values also indicated a slight increase in the *R<sup>g</sup>* of the HSA–HpzA docked complex compared to that of the free state of HSA (Figure 4A, right panel). *Molecules* **2021**, *26*, x 6 of 13

**Figure 4.** Structural compactness of HSA as a function of time. (**A**) Time evolution of the radius of gyration. (**B**) SASA plot of HSA as a function of time. **Figure 4.** Structural compactness of HSA as a function of time. (**A**) Time evolution of the radius of gyration. (**B**) SASA plot of HSA as a function of time.

Dynamics of Intra-/Intermolecular Hydrogen Bonding Hydrogen bonds are necessary for maintaining the structural conformation of a protein [27]. Analysis of hydrogen bonding is useful to assess the stability of a protein and protein–ligand complex [27]. The stability of HSA and the HSA–HpzA docked complex was examined by exploring the formation of intramolecular hydrogen bonds (Figure 5). The analysis of the intramolecular hydrogen bonds formed within HSA indicated the formation of several hydrogen bonds, which maintained the integrity of HSA's three-dimensional structure. The average hydrogen bonds estimated within HSA before and after the HpzA binding were 480 and 467, respectively. The plot indicated a significant decline in SASA delineates the folding mechanism of a protein structure by analyzing the protein's solvent accessibility [26]. SASA values were calculated and plotted to investigate the surface area of HSA and the HSA–HpzA docked complex during the simulation, the (Figure 4B). The results indicated that the SASA values were minorly increased for HSA in ligand-bound state. This slight increase in the values of SASA indicated that a number of internal residues might be exposed to the surface post the binding of HpzA (Figure 4B, left panel). The distribution plot of SASA values also indicated a slight increase in SASA of the HSA–HpzA docked complex (Figure 4B, right panel).

### hydrogen bonds in HSA when it formed a complex with HpzA, which was possibly due to the ligand's occupancy of intramolecular space in the binding pocket of HSA. The bind-Dynamics of Intra-/Intermolecular Hydrogen Bonding

ing of HpzA might cause a hindrance in the formation of intramolecular hydrogen bonding within HSA. However, the plot suggested that this decrement in the formation did not lead to any major shift, which resulted in maintaining the geometry of the HSA structure Hydrogen bonds are necessary for maintaining the structural conformation of a protein [27]. Analysis of hydrogen bonding is useful to assess the stability of a protein and protein–ligand complex [27]. The stability of HSA and the HSA–HpzA docked complex

during the simulation. The PDF distribution also indicated a decrease in intramolecular

**Figure 5.** Intramolecular hydrogen bonds during the simulation. (**A**) Time evolution and stability of hydrogen bonds formed intramolecular within HSA. (**B**) The probability of distribution of hydrogen

The formation of intermolecular hydrogen bonds plays a crucial role in the stability and directionality of a protein–ligand complex [28]. The intermolecular hydrogen bonds can be explored to measure the strength of a ligand towards the binding pocket of a protein. The hydrogen bonds formed between HpzA and HSA were recorded and plotted to

bonding as PDF.

hydrogen bonding within HSA when bound to HpzA (Figure 5B).

was examined by exploring the formation of intramolecular hydrogen bonds (Figure 5). The analysis of the intramolecular hydrogen bonds formed within HSA indicated the formation of several hydrogen bonds, which maintained the integrity of HSA's three-dimensional structure. The average hydrogen bonds estimated within HSA before and after the HpzA binding were 480 and 467, respectively. The plot indicated a significant decline in hydrogen bonds in HSA when it formed a complex with HpzA, which was possibly due to the ligand's occupancy of intramolecular space in the binding pocket of HSA. The binding of HpzA might cause a hindrance in the formation of intramolecular hydrogen bonding within HSA. However, the plot suggested that this decrement in the formation did not lead to any major shift, which resulted in maintaining the geometry of the HSA structure during the simulation. The PDF distribution also indicated a decrease in intramolecular hydrogen bonding within HSA when bound to HpzA (Figure 5B). was examined by exploring the formation of intramolecular hydrogen bonds (Figure 5). The analysis of the intramolecular hydrogen bonds formed within HSA indicated the formation of several hydrogen bonds, which maintained the integrity of HSA's three-dimensional structure. The average hydrogen bonds estimated within HSA before and after the HpzA binding were 480 and 467, respectively. The plot indicated a significant decline in hydrogen bonds in HSA when it formed a complex with HpzA, which was possibly due to the ligand's occupancy of intramolecular space in the binding pocket of HSA. The binding of HpzA might cause a hindrance in the formation of intramolecular hydrogen bonding within HSA. However, the plot suggested that this decrement in the formation did not lead to any major shift, which resulted in maintaining the geometry of the HSA structure during the simulation. The PDF distribution also indicated a decrease in intramolecular hydrogen bonding within HSA when bound to HpzA (Figure 5B).

**Figure 4.** Structural compactness of HSA as a function of time. (**A**) Time evolution of the radius of

Hydrogen bonds are necessary for maintaining the structural conformation of a protein [27]. Analysis of hydrogen bonding is useful to assess the stability of a protein and protein–ligand complex [27]. The stability of HSA and the HSA–HpzA docked complex

*Molecules* **2021**, *26*, x 6 of 13

gyration. (**B**) SASA plot of HSA as a function of time.

Dynamics of Intra-/Intermolecular Hydrogen Bonding

**Figure 5.** Intramolecular hydrogen bonds during the simulation. (**A**) Time evolution and stability of hydrogen bonds formed intramolecular within HSA. (**B**) The probability of distribution of hydrogen bonding as PDF. **Figure 5.** Intramolecular hydrogen bonds during the simulation. (**A**) Time evolution and stability of hydrogen bonds formed intramolecular within HSA. (**B**) The probability of distribution of hydrogen bonding as PDF.

The formation of intermolecular hydrogen bonds plays a crucial role in the stability and directionality of a protein–ligand complex [28]. The intermolecular hydrogen bonds can be explored to measure the strength of a ligand towards the binding pocket of a protein. The hydrogen bonds formed between HpzA and HSA were recorded and plotted to The formation of intermolecular hydrogen bonds plays a crucial role in the stability and directionality of a protein–ligand complex [28]. The intermolecular hydrogen bonds can be explored to measure the strength of a ligand towards the binding pocket of a protein. The hydrogen bonds formed between HpzA and HSA were recorded and plotted to explore their formation and breakdown during the simulation (Figure 6A). Analysis of the generated plot indicated the formation of up to three hydrogen bonds during the simulation time. The analysis clearly indicated that at least one hydrogen bond was formed during the simulation time. This strong hydrogen bond could be correlated with the molecular docking observation that His266 was involved in the conventional hydrogen bond formation. The analysis suggested the importance of intermolecular hydrogen bonding in the ligand binding to form a stable complex. The PDF plot indicated that at least one hydrogen bond was formed with higher probability and stability in the HSA–HpzA docked complex (Figure 6B). *Molecules* **2021**, *26*, x 7 of 13 explore their formation and breakdown during the simulation (Figure 6A). Analysis of the generated plot indicated the formation of up to three hydrogen bonds during the simulation time. The analysis clearly indicated that at least one hydrogen bond was formed during the simulation time. This strong hydrogen bond could be correlated with the molecular docking observation that His266 was involved in the conventional hydrogen bond formation. The analysis suggested the importance of intermolecular hydrogen bonding in the ligand binding to form a stable complex. The PDF plot indicated that at least one hydrogen bond was formed with higher probability and stability in the HSA–HpzA docked complex (Figure 6B).

**Figure 6.** Intermolecular hydrogen bonds during the simulation. (**A**) Time evolution of hydrogen bonds formed between HpzA and HSA. (**B**) The probability of distribution of hydrogen bonds. **Figure 6.** Intermolecular hydrogen bonds during the simulation. (**A**) Time evolution of hydrogen bonds formed between HpzA and HSA. (**B**) The probability of distribution of hydrogen bonds.

#### *2.3. Principal Component and Free Energy Landscape Analyses 2.3. Principal Component and Free Energy Landscape Analyses*

Principal component analysis (PCA) is an arithmetic technique to reduce the complexity of MD trajectories by obtaining the collective motion of Cα atoms. It is vital to evaluate the stability of proteins and protein–ligand complexes [29,30]. PCA was per-Principal component analysis (PCA) is an arithmetic technique to reduce the complexity of MD trajectories by obtaining the collective motion of Cα atoms. It is vital to evaluate the stability of proteins and protein–ligand complexes [29,30]. PCA was performed to

formed to evaluate the conformational projection and transition dynamics of HSA with HpzA. The 2D projections of the first two principal components (PCs), PC1 and PC2, for

first two eigenvectors (EVs), EV1 and EV2, showed the time evolution of the HSA projection, which indicated a similar pattern of trajectories (Figure 7B). The PCA showed that HSA and HSA–HpzA had similar correlated motions without significant change. However, the ligand-bound HSA explored a wider phase space on both EVs, which indicated its higher dynamics. The PCA confirmed that the HSA–HpzA complex was quite stable

**Figure 7.** Principal component analysis. **(A)** Two-dimensional projections of HSA trajectories on eigenvector 1 and eigenvector 2. (**B**) Time evaluation of conformational projections of trajectories on

during the simulation and behaved like the free state of HSA.

eigenvector 1 and eigenvector 2.

evaluate the conformational projection and transition dynamics of HSA with HpzA. The 2D projections of the first two principal components (PCs), PC1 and PC2, for HSA and HSA–HpzA is presented in Figure 7A. The analysis revealed that the 2D projection of HSA showed smaller phase space than that of HSA–HpzA. The projection of the first two eigenvectors (EVs), EV1 and EV2, showed the time evolution of the HSA projection, which indicated a similar pattern of trajectories (Figure 7B). The PCA showed that HSA and HSA–HpzA had similar correlated motions without significant change. However, the ligand-bound HSA explored a wider phase space on both EVs, which indicated its higher dynamics. The PCA confirmed that the HSA–HpzA complex was quite stable during the simulation and behaved like the free state of HSA. plexity of MD trajectories by obtaining the collective motion of Cα atoms. It is vital to evaluate the stability of proteins and protein–ligand complexes [29,30]. PCA was performed to evaluate the conformational projection and transition dynamics of HSA with HpzA. The 2D projections of the first two principal components (PCs), PC1 and PC2, for HSA and HSA–HpzA is presented in Figure 7A. The analysis revealed that the 2D projection of HSA showed smaller phase space than that of HSA–HpzA. The projection of the first two eigenvectors (EVs), EV1 and EV2, showed the time evolution of the HSA projection, which indicated a similar pattern of trajectories (Figure 7B). The PCA showed that HSA and HSA–HpzA had similar correlated motions without significant change. However, the ligand-bound HSA explored a wider phase space on both EVs, which indicated its higher dynamics. The PCA confirmed that the HSA–HpzA complex was quite stable during the simulation and behaved like the free state of HSA.

*Molecules* **2021**, *26*, x 7 of 13

*2.3. Principal Component and Free Energy Landscape Analyses* 

complex (Figure 6B).

explore their formation and breakdown during the simulation (Figure 6A). Analysis of the generated plot indicated the formation of up to three hydrogen bonds during the simulation time. The analysis clearly indicated that at least one hydrogen bond was formed during the simulation time. This strong hydrogen bond could be correlated with the molecular docking observation that His266 was involved in the conventional hydrogen bond formation. The analysis suggested the importance of intermolecular hydrogen bonding in the ligand binding to form a stable complex. The PDF plot indicated that at least one hydrogen bond was formed with higher probability and stability in the HSA–HpzA docked

**Figure 6.** Intermolecular hydrogen bonds during the simulation. (**A**) Time evolution of hydrogen bonds formed between HpzA and HSA. (**B**) The probability of distribution of hydrogen bonds.

Principal component analysis (PCA) is an arithmetic technique to reduce the com-

**Figure 7.** Principal component analysis. **(A)** Two-dimensional projections of HSA trajectories on eigenvector 1 and eigenvector 2. (**B**) Time evaluation of conformational projections of trajectories on eigenvector 1 and eigenvector 2. **Figure 7.** Principal component analysis. **(A)** Two-dimensional projections of HSA trajectories on eigenvector 1 and eigenvector 2. (**B**) Time evaluation of conformational projections of trajectories on eigenvector 1 and eigenvector 2. *Molecules* **2021**, *26*, x 8 of 13

Free energy landscape (FEL) analysis based on PCA has been utilized in representing protein conformations in a time evolution manner [31]. The FEL distinguishes the kinetic and thermodynamic estates of proteins and protein–ligand complexes. The FELs depend on the prospect of an arrangement of specific data points and converting them to free-energy values. FELs analysis was performed to investigate the conformational stability underlying HSA–HpzA binding. Figure 8 shows the FEL plots projected onto PC1 and PC2 of HSA and the HSA–HpzA complex for Cα atoms. The centralized blue zones suggested that the subsequent system was stable at that moment. The dimension and appearance of the minimum energy space, called the global minimum (in blue) in the FELs, indicated that the HSA–HpzA system was more stable than HSA. HSA and HSA–HpzA displayed a single global minimum confined within a large basin. The FEL analysis further validated our previous observations that HSA, when bound with HpzA, was stable and reached a more thermodynamically stable conformational state (Figure 8B). Free energy landscape (FEL) analysis based on PCA has been utilized in representing protein conformations in a time evolution manner [31]. The FEL distinguishes the kinetic and thermodynamic estates of proteins and protein–ligand complexes. The FELs depend on the prospect of an arrangement of specific data points and converting them to freeenergy values. FELs analysis was performed to investigate the conformational stability underlying HSA–HpzA binding. Figure 8 shows the FEL plots projected onto PC1 and PC2 of HSA and the HSA–HpzA complex for Cα atoms. The centralized blue zones suggested that the subsequent system was stable at that moment. The dimension and appearance of the minimum energy space, called the global minimum (in blue) in the FELs, indicated that the HSA–HpzA system was more stable than HSA. HSA and HSA–HpzA displayed a single global minimum confined within a large basin. The FEL analysis further validated our previous observations that HSA, when bound with HpzA, was stable and reached a more thermodynamically stable conformational state (Figure 8B).

**Figure 8.** The Gibbs free energy landscapes of (**A**) free HSA and (**B**) HSA–HpzA obtained during 100 ns MD simulation. **Figure 8.** The Gibbs free energy landscapes of (**A**) free HSA and (**B**) HSA–HpzA obtained during 100 ns MD simulation.

Fluorescence spectroscopy serves as an important technique to provide an insight

idea about the strength of interaction between a protein and ligand [32]. Intrinsic fluorescence depicts changes in the local microenvironment of aromatic amino acid residues and this serves as an important technique in determining the protein-ligand complex formation [33,34]. When a decrease in the fluorescence of native protein is observed with increasing concentration of the ligands, it is referred to as fluorescence quenching [35]. Fluorescence quenching was performed at three different temperatures (15, 20, and 25 °C). Fluorescence emission spectra of free HSA and HSA with different HpzA concentrations (0–11 µM) at various temperatures are shown in Figure 9A–C. Native HSA showed fluorescence emission maxima at 346 nm. There was a noticeable decrease HSA's fluorescence intensity with increasing HpzA concentration, with no peak shift. Fluorescence quenching can be either static, dynamic, or a combination of both [36]. The deviation of the binding parameters with temperature delineates the operative quenching for a particular interaction. The quenching data obtained were fitted into the Stern–Volmer, modified Stern– Volmer, and van 't Hoff equations to obtain various binding and thermodynamic parameters for HSA–HpzA interaction as per previously published reports [37,38]. Table 1 shows the obtained values of the Stern–Volmer constant (*K*sv) at various temperatures, which were found to increase with increasing temperature, revealing the existence of the dynamic mode. Moreover, to confirm the quenching mode, a biomolecular quenching rate constant (*K*q) was calculated using *K*sv = *K*q/τ0 (τ0 = 2.7 × 10−9 s), and the value was found to be higher than that of the maximum dynamic quenching constant (nearly 1010 M−1 s−1) [39],

*2.4. Fluorescence-Based Binding* 
