Inter-residue Interactions Research Articles

Weighted families of contact maps to characterize conformational ensembles of (highly-)flexible proteins.

Characterizing the structure of flexible proteins, particularly within the realm of intrinsic disorder, presents a formidable challenge due to their high conformational variability. Currently, their structural representation relies on (possibly large) conformational ensembles derived from a combination of experimental and computational methods. The detailed structural analysis of these ensembles is a difficult task, for which existing tools have limited effectiveness. This study proposes an innovative extension of the concept of contact maps to the ensemble framework, incorporating the intrinsic probabilistic nature of disordered proteins. Within this framework, a conformational ensemble is characterized through a weighted family of contact maps. To achieve this, conformations are first described using a refined definition of contact that appropriately accounts for the geometry of the inter-residue interactions and the sequence context. Representative structural features of the ensemble naturally emerge from the subsequent clustering of the resulting contact-based descriptors. Importantly, transiently populated structural features are readily identified within large ensembles. The performance of the method is illustrated by several use cases and compared with other existing approaches, highlighting its superiority in capturing relevant structural features of highly flexible proteins. An open-source implementation of the method is provided together with an easy-to-use Jupyter notebook, available at https://gitlab.laas.fr/moma/WARIO.

Molecular modelling of Pro-Gly glyproline and its complexes

In the present work, Pro-Gly glyproline and its biologically active complexes with heparin and with the antithrombin-heparin protein complex have been studied using molecular modelling methods. This molecule is used as a drug and exhibits antithrombotic, anticoagulant immunomodulatory, antiulcer, and antidiabetic effects. The conformational profiles of the dipeptide were studied within the molecular mechanics framework, the geometry and energy of intra- and interresidual interactions were estimated for the most stable states of the investigated dipeptide. A comparative analysis of the extended and folded structures of this molecule optimized by the DFT/B3LYP/6-31+G(d,p) method was carried out, the frontier molecular orbitals were calculated, the surface of the molecular electrostatic potential was studied, the effective charges on the atoms were determined, and the dipole moment and polarizability values were calculated. Molecular docking revealed that the folded structure of Pro-Gly glyproline showed a higher affinity value for both heparin and the antithrombin-heparin protein complex. A model of the pharmacophore of the Pro-Gly molecule for its interaction with a specific receptor was proposed. The obtained data form the basis for the development of effective analogs of glyprolines in the complex with heparin, which are more specific concerning their putative target — antithrombin.

SPATIAL STRUCTURE OF THE HEPTAPEPTIDE ANALOGUE OF NOCICEPTIN MOLECULE

<p>It is known that nociceptins are a new type of regulatory peptide. Knowledge of these peptide molecules' structural and functional properties is of great practical importance for medicine and pharmacology. Their mechanisms of action are considered anti-opioid. This scientific work is devoted to studying the spatial structure of the heptapeptide H-Phe1-Gly2-Gly3-Phe4-Val5-Gly6-Pro7-OH. It examines the conformational capabilities of this heptapeptide molecule. This neuropeptide molecule is a stable analog of the nociceptin. The biologically active conformation of the peptide molecule, which is realized upon interaction with the receptor, is included in the set of low-energy structures. Therefore, studying the spatial structure of peptide molecules is of great interest. Theoretical conformational analysis about nonvalent, electrostatic, and torsional interactions, the energy of the hydrogen bonds, and a special computer program carried out the calculations. The 10 low-energy conformations of this molecule and the values of the dihedral angles of the main chain and side chains are found, and the energy of the intra- and inter-residue interactions is estimated. It is revealed that low energy conformations of this molecule have the half-folded and folded type of backbone. The side chains of the Phe1 and Phe4 amino acids in low-energy conformations carry out effective interactions and are conformationally labile amino acids; they bring together the regions of the main chain and the side chains of the amino acids included in the heptapeptide. These folded forms bring parts of the backbone and the amino acids' side chains together, resulting in important interactions.</p>

СТРУКТУРНАЯ ОРГАНИЗАЦИЯ ГЛИПРОЛИНОВОЙ ГЕКСАПЕПТИДНОЙ МОЛЕКУЛЫ

Computer modeling based on the use of the method of theoretical conformational analysis and programs that allow obtaining a graphic image of the spatial structures of biomolecules was performed for a hexapeptide molecule – glyproline H-Phe-Gly-Gly-Phe-Gly-Pro-OH. Glyprolines are short peptides whose amino acid sequences contain glycine and proline residues. Their mechanisms of action are currently poorly understood. Using the method of molecular mechanics, the spatial structure and conformational possibilities of this hexapeptide molecule were determined. Its potential energy was estimated as the sum of non-valente, electrostatic, torsion interactions and the energy of hydrogen bonds. 9 low-energy structures were found for the glyproline hexapeptide, the values of the dihedral angles of the main and side chains of the amino acids included in the molecule. The energy of intra- and interresidual interactions was estimated. The calculation showed that the folded forms of the main chain are low-energy for the hexapeptide. In them, the side chains of amino acids phenylalanine and proline, which are conformationally rigid, carry out effective interactions with all parts of the hexapeptide molecule.

СТРУКТУРНАЯ ОРГАНИЗАЦИЯ МОЛЕКУЛ СОЙМОРФИНОВ

The conformational capabilities of soymorphine-5 (Thr1-Pro2-Phe3-Val4-Val5-NH2), soymorphine-6 (Tyr1-Pro2-Phe3-Val4-Val5-Asn6-NH2) and soymorphine-7 (Tyr1- Pro2-Tyr3-Val4-Val5-Asn6-Ala7-NH2) molecules have been studied by the method of theoretical conformational analysis. The potential function of the system is chosen as the sum of non-valence, electrostatic and torsion interactions and the energy of hydrogen bonds. The low-energy conformations of soymorphine-5, soymorphine-6 and soymorphine-7 molecules were found, the dihedral angles of the main and side chains of amino acid residues that make up the molecule were found, and the energy of intra- and interresidual interactions was estimated. Thus, the spatial structure of soymorphine-5, soymorphine-6 and soymorphine-7 molecules can be represented by eight structural types. It can be assumed that the molecules perform their physiological functions in these structures. Comparison of the low-energy structures of soymorphins shows that in all molecules the first four low-energy conformations are representatives of the structural types efef, efee, efff, effe for soymorphine-5, effff, efeff, efffe, effee for soymorphine-6, efffff, efeffe, efffef, effeee for soymorphine-6. soymorphine-7. On the basis of these structures, it is possible to propose their artificial analogues for synthesis. It was shown that the spatial structure of soymorphine-5, soymorphine-6 and soymorphine-7 molecules is represented by the conformations of eight shapes of the peptide skeleton. The results obtained can be used to elucidate the structural and structural-functional organization of soymorphine molecules.

КОНФОРМАЦИОННЫЕ СВОЙСТВА И ЭЛЕКТРОННАЯ СТРУКТУРА ПРОТИВООПУХОЛЕВОГО АГЕНТА TYR-SER-LEU

Molecular modeling methods were used to study the features of the spatial and electronic structure of the antitumor tripeptide YSL (Tyr-Ser-Leu), developed by Chinese scientists. The conformational analysis of the molecule revealed a limited set of its energetically preferable conformational states in a certain range of relative energy. The nature of the forces stabilizing the low-energy conformations of the tripeptide molecule was determined. As a result of the study, the energetically preferable ranges of dihedral angles, the energy contributions of interresidual interactions and hydrogen bonds, as well as the mutual arrangement of residues and their side chains in low-energy conformations of the tripeptide were also determined. Using the methods of molecular mechanics, the energy contributions of intramolecular interactions in low-energy conformational states of the molecule were obtained. Based on quantum-chemical calculations, the distribution of electron density and the values of dipole moments of the most optimal spatial structures of the YSL tripeptide molecule were determined. The reactivity of the tripeptide was also studied by quantum chemical calculations based on the obtained electronic characteristics of each low-energy conformation of the molecule. Using the calculated coordinates of the atoms of the energetically preferable structures of the molecule, their molecular models were built, the comparison of which makes it possible to identify the structural criteria necessary to create a drug suitable for clinical use.

Microbial Biofilm Inhibition Using Magnetic Cross-Linked Polyphenol Oxidase Aggregates.

Microbial biofilm accumulation poses a serious threat to the environment, presents significant challenges to different industries, and exhibits a large impact on public health. Since there has not been a conclusive answer found despite various efforts, the potential green and economical methods are being focused on, particularly the innovative approaches that employ biochemical agents. In the present study, we propose a bio-nanotechnological method using magnetic cross-linked polyphenol oxidase aggregates (PPO m-CLEA) for inhibition of microbial biofilm including multidrug resistant bacteria. Free PPO solution showed only 55-60% biofilm inhibition, whereas m-CLEA showed 70-75% inhibition, as confirmed through microscopic techniques. The carbohydrate and protein contents in biofilm extracellular polymeric substances (EPSs) were reduced significantly. The m-CLEA demonstrated reusability up to 5 cycles with consistent efficiency in biofilm inhibition. Computational work was also done where molecular docking of PPO with microbial proteins associated with biofilm formation was conducted, resulting in favorable binding scores and inter-residual interactions. Overall, both in vitro and in silico results suggest that PPO interferes with microbial cell attachment and EPS formation, thereby preventing biofilm colonization.

Spatial Structure of the Casoxin D Molecule

ABSTRACT The conformational possibilities of the casoxin D (Tyr1-Val2-Pro3-Phe4-Pro5-Pro6-Phe7) molecule were studied by theoretical conformational analysis. The potential function of the system is chosen as the sum of non-valence, electrostatic and torsion interactions and the energy of hydrogen bonds. The low-energy conformations of the casoxin D molecule, the dihedral angles of the main and side chains of the amino acid residues that make up the molecule were found, and the energy of intra- and interresidual interactions was estimated. It has been shown that the spatial structure of the casoxin D molecule is represented by conformations of eight shapes of the peptide skeleton. The results obtained can be used to elucidate the structural and structural-functional organisation of casoxin molecules.

Importance of Inter-residue Contacts for Understanding Protein Folding and Unfolding Rates, Remote Homology, and Drug Design.

Inter-residue interactions in protein structures provide valuable insights into protein folding and stability. Understanding these interactions can be helpful in many crucial applications, including rational design of therapeutic small molecules and biologics, locating functional protein sites, and predicting protein-protein and protein-ligand interactions. The process of developing machine learning models incorporating inter-residue interactions has been improved recently. This review highlights the theoretical models incorporating inter-residue interactions in predicting folding and unfolding rates of proteins. Utilizing contact maps to depict inter-residue interactions aids researchers in developing computer models for detecting remote homologs and interface residues within protein-protein complexes which, in turn, enhances our knowledge of the relationship between sequence and structure of proteins. Further, the application of contact maps derived from inter-residue interactions is highlighted in the field of drug discovery. Overall, this review presents an extensive assessment of the significant models that use inter-residue interactions to investigate folding rates, unfolding rates, remote homology, and drug development, providing potential future advancements in constructing efficient computational models in structural biology.

Elucidating the molecular mechanism of noncompetitive inhibition of acetylcholinesterase by an antidiabetic drug chlorpropamide: identification of new allosteric sites.

Acetylcholinesterase (AChE) has emerged as an important drug target for the treatment of neurodegenerative disorders such as Alzheimer's disease (AD). Recent experimental studies indicate that certain antidiabetic drugs can be repurposed as potent AChE inhibitors. Enzymatic kinetic assays suggest that the antidiabetic drug chlorpropamide (CPM) acts as a noncompetitive inhibitor, but the mechanism of action and the binding site(s) of interaction with AChE are not known. In this work, we have carried out molecular dynamics (MD) simulations to discover a new allosteric site in addition to the known peripheral anionic site (PAS) as a potential binding site of this noncompetitive inhibitor. We show that the conformational ensemble of the catalytic triad, particularly the HIS447, undergoes a significant population shift on ligand binding that is responsible for deactivation of the enzyme. We also elucidate the pathway of the allosteric signaling in terms of locally correlated domains of the inter-residue interaction network. Thus, our work identifies a new allosteric site for AChE inhibition and eludiates the underlying mechanistic principles. These results would be useful for the rational design of new noncompetitive inhibitors for AChE.

Spatial structure of soymorphin-6 molecule

By the method of the theoretic conformational analysis the conformational capabilities of the soymorphin-6 (Tyr1-Pro2-Phe3-Val4-Val5-Asn6-NH2) molecule were studied. The potential energy of this molecule is given as the sum of nonvalent, electrostatic and torsional interactions and the energy of hydrogen bonds. Low-energy conformations of the soymorphin-6 molecule, the values of the dihedral angles of the main and side chains of amino acid residues that make up the molecules are found, the energy of intra- and inter-residual interactions is estimated. It is shown that the spatial structure of the soymorphin-6 molecule can be represented by eight conformations. The results obtained can be used to study the structure and structural-functional organization of the soymorphin molecules.

Comparison of structural networks across homologous proteins.

Protein sequence determines its structure and function. The indirect relationship between protein function and structure lies deep-rooted in the structural topology that has evolved into performing optimal function. The evolution of structure and its interconnectivity has been conventionally studied by comparing the root means square deviation between protein structures at the backbone level. Two factors that are necessary for the quantitative comparison of non-covalent interactions are (a) explicit inclusion of the coordinates of side-chain atoms and (b) consideration of multiple structures from the conformational landscape to account for structural variability. We have recently addressed these fundamental issues by investigating the alteration of inter-residue interactions across an ensemble of protein structure networks through a graph spectral approach. In this study, we have developed a rigorous method to compare the structure networks of homologous proteins, with a wide range of sequence identity percentages. A range of dissimilarity measures that show the extent of change in the network across homologous structures are generated, which also includes the comparison of the protein structure variability. We discuss in detail, scenarios where the variation of structure is not accompanied by loss or gain of the overall network and its vice versa. The sequence-based phylogeny among the homologs is also compared with the lineage obtained from information from such a robust structure comparison. In summary, we can obtain a quantitative comparison score for the structure networks of homologous proteins, which also enables us to study the evolution of protein function based on the variation of their topologies.

Allosteric Determinants in High Temperature Requirement A Enzymes Are Conserved and Regulate the Population of Active Conformations.

High temperature requirement A (HtrA) are allosterically regulated enzymes wherein effector binding to the PDZ domain triggers proteolytic activity. Yet, it remains unclear if the inter-residue network governing allostery is conserved across HtrA enzymes. Here, we investigated and identified the inter-residue interaction networks by molecular dynamics simulations on representative HtrA proteases, Escherichia coli DegS and Mycobacterium tuberculosis PepD, in effector-bound and free forms. This information was used to engineer mutations that could potentially perturb allostery and conformational sampling in a different homologue, M. tuberculosis HtrA. Mutations in HtrA perturbed allosteric regulation─a finding consistent with the hypothesis that the inter-residue interaction network is conserved across HtrA enzymes. Electron density from data collected on cryo-protected HtrA crystals revealed that mutations altered the topology of the active site. Ensemble models fitted into electron density calculated from room-temperature diffraction data showed that only a fraction of these models had a catalytically competent active site conformation alongside a functional oxyanion hole thus providing experimental evidence that these mutations influenced conformational sampling. Mutations at analogous positions in the catalytic domain of DegS perturbed the coupling between effector binding and proteolytic activity, thus confirming the role of these residues in the allosteric response. The finding that a perturbation in the conserved inter-residue network alters conformational sampling and the allosteric response suggests that an ensemble allosteric model best describes regulated proteolysis in HtrA enzymes.

Designing the 5HT2BR structure and its modulation as a therapeutic target for repurposing approach in drug-resistant epilepsy

The study intends to repurpose FDA drugs and investigate the mechanism of (5HT2BR) activation by comprehending inter-residue interactions. The 5HT2BR is a novel thread, and its role in reducing seizures in Dravet syndrome is emerging. The crystal structure (5HT2BR) is a chimera with mutations; hence 3D-structure is modeled (4IB4: 5HT2BRM). The structure is cross-validated to simulate the human receptor using enrichment analysis (ROC: 0.79) and SAVESv6.0. Virtual screening of 2456 approved drugs yielded the best hits that are subjected to MM/GBSA and molecular dynamic (MD) simulations. The 2 top drugs Cabergoline (−53.44 kcal/mol) and Methylergonovine (−40.42 kcal/mol), display strong binding affinity, and ADMET/SAR analysis also suggests their non-mutagenic or non-carcinogenic nature. Methylergonovine has a weaker binding affinity and lower potency than standards [Ergotamine (agonist) and Methysergide (antagonist)] due to its higher Ki (1.32 M) and Kd (6.44 ×10−8 M) values. Compared to standards, Cabergoline has moderate binding affinity and potency [Ki = 0.85 M and Kd = 5.53 × 10–8 M]. The top 2 drugs primarily interact with conserved residues (ASP135, LEU209, GLY221, ALA225, and THR140) as in agonists, unlike the antagonist. The top 2 drugs, upon binding to the 5HT2BRM, modify the helices VI, V, and III and shift the RMSD 2.48 Å and 3.07 Å. LEU209 forms a latch with residues 207–214 (forms a lid) in the 5HT2BRM receptor, which enhances agonist binding and prevents drug escape. Methylergonovine and Cabergoline interact more stongly with ALA225 than the antagonist. The post-MD analysis of Cabergoline suggests a better MM/GBSA value (−89.21 kcal/mol) than Methylergonovine (−63.54 kcal/mol). In this study, Cabergoline and Methylergonovine's agonistic mechanism and solid binding properties suggest their strong role in regulating the 5HT2BR and might target drug-resistant epilepsy.

Creatinase: Using Increased Entropy to Improve the Activity and Thermostability.

Improving protein thermostability in mutagenesis-based enzyme engineering was often achieved by enhancing interresidue interactions via mutation to increase the enthalpy penalty of unfolding. However, this approach may trade off the functional activity due to the loss of structural flexibility of the biomolecule. Here, by performing X-ray crystallography, enzymatic kinetic experiments, neutron scattering, and thermodynamical measurements, we compared the structures, catalytic behaviors, dynamics, and thermostability between a wild-type creatinase and its four-point mutant. We found that the mutant is an entropy-driven thermostable protein with higher structural flexibility, i.e., higher conformational entropy, in the folded state compared to the wild type. The increased conformational entropy of the mutant in the folded state can reduce the entropy gain during unfolding and thus renders it greater thermostability. Moreover, the increased structural flexibility, particularly around the catalytic site, can broaden the mutant's working temperature range and considerably improve its activity at ambient conditions, which is crucial for its application in diagnosing kidney diseases. Complementary all-atom molecular dynamics simulations indicated that the four mutations replaced several of the strong interresidue interactions (electrostatic interactions and hydrogen bonds) with weak hydrophobic interactions. These substitutions not only release the structural flexibility to promote the thermostability and enzymatic activity of the protein but they also preserve the protein structure from collapsing. Our findings may pave a route for the entropy-driven strategy to design proteins with high thermostability and activity.

Comparision of structural network between homologous proteins with variability derived from multiple conformers.

Protein sequence determines its structure and hence its function. This indirect relationship between protein function and its structure lies deep rooted in the structural topology that has evolved into performance of optimal function. The evolution of structure and its interconnectivity have been studied predominantly using protein structure comparison metrics that quantify deviation from one another. A fundamental issue in computational structural biology is that the comparison of proteins by recruiting representative states of their structures is insufficient, as it does not account for their structural variability. We know how dynamic proteins are and it is necessary to include the variation in their topology while comparing them with their homologs. This fundamental issue is addressed here by looking into the alteration of inter-residue interactions across an ensemble of protein structure networks, instead of a single representative structure. A graph spectral comparison of the protein structural networks is employed to quantify the extent of dissimilarity in the structural topologies. We discuss in detail, scenarios where the variation of structure is not accompanied by loss or gain of overall network and its vice versa. Also, a newly obtained metric of topological change collated from the network dissimilarity score is used to obtain the phylogeny between the homologs and is compared with the sequence-based phylogeny. The method can also help to study the evolution of protein function based on the variation of their topologies.

Introducing electrostatic stimulation to normal mode analysis.

Normal Mode Analysis (NMA) has proven to be a powerful tool to understand the macromolecular movement of protein complexes. While less precise than an all-atom molecular dynamics simulation, it allows to predict movements on timescales that are not yet accessible to MD simulations on a routine basis. In this work, we build on the previous successful implementation EnCoM to include electrostatic interactions in the NMA analysis. The motivation behind this step is to study the conformational rearrangements of voltage-gated ion channels that are sensitive to the electric field over the membrane. We show here preliminary results on the extension of the inter-residue interactions to include electrostatics and first results on the effect of an electric field on the normal modes of a voltage-gated potassium (Kv) channel.

Revisiting the Sweet Taste Receptor T1R2-T1R3 through Molecular Dynamics Simulations Coupled with a Noncovalent Interactions Analysis.

It is nowadays widely accepted that sweet taste perception is elicited by the activation of the heterodimeric complex T1R2-T1R3, customarily known as sweet taste receptor (STR). However, the interplay between STR and sweeteners has not yet been fully clarified. Here through a methodology coupling molecular dynamics and the independent gradient model (igm) approach we determine the main interacting signatures of the closed (active) conformation of the T1R2 Venus flytrap domain (VFD) toward aspartame. The igm methodology provides a rapid and reliable quantification of noncovalent interactions through a score (Δginter) based on the attenuation of the electronic density gradient when two molecular fragments approach each other. Herein, this approach is coupled to a 100 ns molecular dynamics simulation (MD-igm) to explore the ligand-cavity contacts on a per-residue basis as well as a series of key inter-residue interactions that stabilize the closed form of VFD. We also apply an atomic decomposition scheme of noncovalent interactions to quantify the contribution of the ligand segments to the noncovalent interplay. Finally, a series of structural modification on aspartame are conducted in order to obtain guidelines for the rational design of novel sweeteners. Given that innovative methodologies to reliably quantify the extent of ligand-protein coupling are strongly demanded, this approach combining a noncovalent analysis and MD simulations represents a valuable contribution, that can be easily applied to other relevant biomolecular systems.

Identification of an Evolutionarily Conserved Allosteric Network in Steroid Receptors.

Allosteric pathways in proteins describe networks comprising amino acid residues which may facilitate the propagation of signals between distant sites. Through inter-residue interactions, dynamic and conformational changes can be transmitted from the site of perturbation to an allosteric site. While sophisticated computational methods have been developed to characterize such allosteric pathways linking specific sites on proteins, few attempts have been made to apply these approaches toward identifying new allosteric sites. Here, we use molecular dynamics simulations and suboptimal path analysis to discover new allosteric networks in steroid receptors with a focus on evolutionarily conserved pathways. Using modern receptors and a reconstructed ancestral receptor, we identify networks connecting several sites to the activation function surface 2 (AF-2), the site of coregulator recruitment. One of these networks is conserved across the entire family, connecting a predicted allosteric site located between helices 9 and 10 of the ligand-binding domain. We investigate the basis of this conserved network as well as the importance of this site, discovering that the site lies in a region of the ligand-binding domain characterized by conserved inter-residue contacts. This study suggests an evolutionarily importance of the helix 9-helix 10 site in steroid receptors and identifies an approach that may be applied to discover previously unknown allosteric sites in proteins.

A molecular dynamics and quantum mechanical investigation of intermolecular interaction and electron-transfer mechanism between copper-containing nitrite reductase and redox partner pseudoazurin.

Much of biological electron transfer occurs between proteins. These molecular processes usually involve molecular recognition and intermolecular electron transfer (inter-ET). The inter-ET reaction between copper-containing nitrite reductase (CuNiR) and partner protein pseudoazurin (PAz) is the first step in denitrification, which is affected by intermolecular association. However, the transient interaction between CuNiR and PAz and the indistinct inter-ET pathway pose challenges for people to understand the biological functions of the CuNiR-PAz complex. Thus, molecular dynamics simulation and quantum mechanical calculation were used to investigate the question in this study. The interaction of the interface residues was determined through hydrogen bonds, root-mean-square deviation, root-mean-square fluctuation, the dynamics cross-correlation matrix, and molecular mechanics Poisson-Boltzmann surface area of molecular dynamics simulations. The interactions among the residues Glu89, Gly200, Asp205, Asn91, Glu204, Thr92, and Met141 on CuNiR and the residues Lys109, Ala15, Lys10, Asn9, Ile110, Met84, and Met16 on PAz are responsible for the stabilization of the complex. The binding free energy is up to -25.33 kcal mol-1. We compared the wild-type and mutant (M84A) interfacial optimized complex models at the CAM-B3LYP level with Grimme dispersion corrections (GD3) to confirm Met84 as a relay station for promoting the inter-ET. Additionally, to test whether Met84 may combine with the adjacent Met141 to form a special two-center, three-electron (S∴S)+ structure to promote the inter-ET, QM/MM was further performed to discuss the possibility of generating an electron stepping stone. Our study will promote a deep understanding of the stable protein-protein interaction, and the identified inter-residue interaction will be theoretical guidance for enhancing the catalytic activity of CuNiR in denitrification.